Mutations in the hepatocyte nuclear factor-1α (HNF-1a) gene cause maturity-onset diabetes of the young (MODY). Approximately 30% of these mutations generate mRNA transcripts harboring premature termination codons (PTCs). Degradation of such transcripts by the nonsense-mediated decay (NMD) pathway has been reported for many genes. To determine whether PTC mutant transcripts of the HNF-1α gene elicit NMD, we have developed a novel quantitative RT-PCR assay. We performed quantification of ectopically expressed mutant transcripts relative to normal transcripts in lymphoblastoid cell lines using a coding single nucleotide polymorphism (cSNP) as a marker. The nonsense mutations R171X, I414G415ATCG→CCA, and P291fsinsC showed reduced mutant mRNA expression to 40% (P = 0.009), <0.01% (P ≤ 0.0001), and 6% (P = 0.001), respectively, of the normal allele. Transcript levels were restored using the translation inhibitor cycloheximide, indicating that the instability arises from NMD. The missense mutations G207D and R229P did not show NMD although R229P exhibited moderate RNA instability. This study provides the first evidence that HNF-1α PTC mutations may be subject to NMD. Mutations that result in significant reduction of protein levels due to NMD will not have dominant-negative activity in vivo. Haploinsufficiency is therefore likely to be the most important mutational mechanism of HNF-1α mutations causing MODY.
Maturity-onset diabetes of the young (MODY) is a subtype of diabetes characterized by early age of onset (usually <25 years), autosomal dominant inheritance, and a progressive defect in β-cell function (1). In U.K. populations, mutations in the hepatocyte nuclear factor-1α (HNF-1α) gene account for ∼65% of cases (2). Over 120 different HNF-1α mutations have been reported (3), of which ∼30% are nonsense or frameshift and lead to the production of premature termination codons (PTCs) (4). The most common HNF-1α mutation results from the insertion of a C nucleotide in a polyC tract in exon 4 (2,5–7). This mutation, P291fsinsC, accounts for ∼20% of families with HNF-1α mutations (4,8,9).
HNF-1α mutations might produce the MODY phenotype by haploinsufficiency or a dominant-negative mutational mechanism. There is considerable support for haploinsufficiency; a mutation in the HNF-1α promoter that disrupts the HNF-4α binding site results in a phenotype indistinguishable from mutations in the coding region (10). Many mutations when tested in vitro have little transactivation activity and do not act as a dominant negative (9). It has been suggested that some mutations may act through a dominant-negative mutational mechanism. The HNF-1α protein acts as a homodimer (or heterodimer with HNF-1β); therefore, mutant proteins lacking transactivation capability but retaining DNA binding and dimerization might exert a dominant-negative effect through the sequestration of wild-type proteins in inactive dimers. The commonest mutation, P291fsinsC, shows dominant-negative activity when overexpressed in MIN6 mouse insulinoma, C33 human epithelial cervical carcinoma, and INS-1 rat cells in vitro (8,9,11), but not in HeLa cells (12).
It has not been possible to analyze mRNA in patients with diabetes caused by HNF-1α mutations because tissues with high levels of HNF-1α expression (liver, pancreas, kidney, and gut) are inaccessible (13). Highly complex mechanisms occurring in living cells may be poorly represented by in vitro systems. For example, many transcripts harboring PTCs do not produce truncated proteins since they are degraded by the nonsense-mediated decay (NMD) pathway (14). This process is dependent on translation and requires the presence of at least one intron 3′ of the PTC for its action (15). Reductions in the levels of mutant mRNA transcripts in fibroblasts or lymphoblastoid cell lines have been shown to be associated with PTC mutations in the fibrillin (16), β-globin (17), Col1A1 (18), and BRCA1 genes (19). Although only a limited number of human genes have been studied, NMD is likely to be a global process that affects the majority of genes (20). To date, there have been no studies to determine whether HNF-1α mutant mRNA transcripts are subject to NMD.
In a recent publication we used a novel RT-PCR-based method to amplify HNF-1α transcripts from lymphoblastoid cells to show that three splice site mutations resulted in PTCs (21). Interestingly, the mutant transcripts were identified by sequence analysis in only one of three replicate RT-PCRs. This suggested that the mutant mRNA transcripts were less abundant than their wild-type counterparts. We hypothesized that mutant mRNAs harboring PTCs may be subject to degradation by the NMD pathway and sought to develop a quantitative method to investigate this further.
We now report a quantitative real-time PCR method to determine the relative amount of mutant and wild-type mRNA transcripts in an allele-specific manner using a coding single nucleotide polymorphism (cSNP) as a marker. Since amplification of HNF-1α mRNA from peripheral lymphocytes was unsuccessful, lymphoblastoid cell lines were established from diabetic HNF-1α mutation carriers and unaffected nonmutation carriers. All subjects selected were heterozygous for the L459L C>T polymorphism in exon 7, and it was possible to determine from family cosegregation studies which allele was inherited with the mutation. The assay is based on the fluorescent resonance energy transfer principle (Fig. 1A) (22) and relies upon the differential melting points of matched and mismatched probe-to-template duplexes for allele specificity (Fig. 1B). The use of two probes with different emission spectra allows the quantitation of both transcripts in an individual heterozygous for the cSNP. The assay was validated and tested for specificity and was shown to be sensitive and accurate over a wide range of template concentrations (data not shown). No significant differences in the relative abundance of alleles were seen between RNA extractions, reverse transcriptions, or PCR runs in control lines (P = 0.912, 0.796, and 0.222, respectively, as tested pairwise by the Mann-Whitney statistic).
We quantified the levels of mutant mRNA transcripts relative to wild-type transcripts for three mutations predicted to lead to the formation of PTCs at codons 171, 315, and 466 (mutations R171X, P291fsinsC, and I414G415ATCG→CCA), two missense mutations (G207D and R229P), and two unaffected controls. P values were determined by comparison of the median ΔCt value of mutated cell lines with the median ΔCt value obtained from the control cell lines by Mann-Whitney U test. Using control cell lines with the L459L polymorphism, we demonstrated that both alleles were expressed at similar levels (P = 0.156). All three nonsense mutant transcripts showed reduced expression (to 40%, P = 0.009; <0.01%, P ≤ 0.0001; and 6%, P = 0.001 for R171X, I414G415ATCG→CCA, and P291fsinsC, respectively) compared with the normal transcript (Table 1). The mutant transcript levels were increased by treatment with the NMD inhibitor cycloheximide by 1.8-fold (P = 0.05), 3.7 × 105-fold (P ≤ 0.0001), and 81-fold (P = 0.003) for R171X, I414G415ATCG→CCA, and P291fsinsC, respectively, when compared with the solvent control (Table 2). These results indicate that the observed reduction in mutant allele frequency is due to NMD.
The variable levels of NMD seen in HNF-1α (mutant allele present at <0.01–40% of wild type) are comparable with those reported for PTC-containing transcripts in the fibrillin and BRCA1 genes where reductions from <10 to 50% of the normal allele have been described (16,19). The reduced fibrillin mutant mRNA levels were classified as severely depleted (<10% of wild type), moderately depleted (10–29% of wild type), or mildly depleted (30–50% of wild type), as determined by pulse-chase analysis (16). By this classification scheme, the R171X nonsense mutation would be described as mildly depleted and the I414G415ATCG→CCA and P291fsinsC frameshift mutations as severely depleted. The efficiency of NMD may be affected by both the type of stop codon, the surrounding sequence, and differential read-through (23).
The missense mutation G207D showed no reduction in mutant transcripts that were present at ∼93% of the normal allele (P = 0.546). The transcript carrying the missense mutation R229P showed an unexpected reduction to 52% of the normal allele (P = 0.03). This reduction was not restored by cycloheximide, indicating that the instability is independent of NMD. Reduced protein levels were also observed for the missense mutants Y122C, R159Q, K205Q, and P447L in in vitro transient transfection experiments (9). The mechanisms resulting in RNA instability are not yet clear, but missense mutations can sometimes cause alterations in mRNA folding that leads to mRNA instability (recently reported for the dopamine receptor (24) or may disrupt sequences involved in stabilization or polyadenylation (25).
In the absence of pancreatic tissue samples from patients with HNF-1α gene mutations, we have exploited the illegitimate transcription of the HNF-1α gene in lymphoblastoid cell lines derived from patients. This approach has previously been used to study other genes that are also expressed in inaccessible tissues, e.g., BRCA1 (19). Confirmation of our results using pancreatic mRNA from mutation carriers would be desirable but is clearly not possible, so we cannot be certain that the same degree of NMD occurs in the pancreas.
The advantage of using lymphoblastoid cells over engineered cell lines, transient transfections, or transgenic mice is that the mutations are present in the heterozygous state and the genomic structure of the gene is intact. Previous reports showed that the overexpression of P291fsinsC inhibited the endogenous activity of HNF-1α in C33 human epithelial cervical carcinoma, INS-1 rat, and MIN6 mouse insulinoma cells in vitro (8,9,11). However, this dominant-negative effect was not observed at 50:50 ratios of mutant-to-wild-type protein, and the absence of introns in the cDNA constructs may preclude the degradation of PTC mutations by NMD. Transgenic mice expressing P291fsinsC in pancreatic β-cells develop progressive hyperglycemia (11) and provide a closer model for HNF-1α mutations in humans than the homozygous knockout mouse that exhibits liver and kidney dysfunction (13) in addition to diabetes (26). However, the expression levels of HNF-1α protein in the transgenic mice cannot be directly compared with humans since the cDNA constructs were present in between 5 and 20 copies per genome and are intronless (27). In keeping with there not being a significant dominant-negative effect in vivo, there is no evidence to suggest a more severe phenotype in patients with diabetes caused by mutations that have a dominant-negative effect in vitro (8,12).
The previous reports of mutations in the HNF-1α gene promoter and dimerization domain provide evidence that haploinsufficiency can cause MODY (10). PTC mutations are also likely to result in haploinsufficiency as a consequence of a reduction in mRNA level due to NMD. We would suggest that even the common HNF-1α mutation P291fsinsC, which has a dominant-negative effect in vitro (8,9,11), may result in diabetes as the result of haploinsufficiency because insufficient truncated protein would be translated in order to have a dominant-negative effect in vivo.
In conclusion, our study utilized ectopic transcription in lymphoblastoid cell lines to show that premature truncating codon mutations are subject to NMD in vivo. This is likely to be an important disease mechanism in the 30% of nonsense, frameshift, and splicing mutations that result in a PTC. Haploinsufficiency may be the common mutational mechanism for the majority of HNF-1α mutations. However, final proof will require the analysis of mRNA expression in legitimately expressing tissues from patients with HNF-1α mutations.
RESEARCH DESIGN AND METHODS
Informed consent was obtained from patients and unaffected family members. Mutations were detected by direct sequencing of genomic DNA (2). Cell lines were established from families DUK05 (R171X) (2), DUK226 (I414G415ATCG→CCA) (28), and DUK582 (P291fsinsC). The novel missense mutations were identified in families DUK16 (R229P-CGA>CCA) (29) and DUK18 (G207D-GGC>GAC). The latter family included nine affected members in four generations who were diagnosed between the ages of 11 and 50 years. All patients and control subjects selected for this study were heterozygous for the exon 7 cSNP 1375C>T (L459L), as demonstrated by sequence analysis and PstI restriction digest. In all cases, the mutation was in phase with the C allele
Transformation of cell lines and tissue culture.
The cell lines MY0104 (R171X), ABO154 (I414G415ATCG→CCA), MYO414 (P291fsinsC), ABO169 (G207D), ABO269 (R229P), ABO179 (control 1), and ABO133 (control 2) were established from Epstein Barr virus transformation of peripheral blood lymphocytes from mutation carriers and unaffected family members by the European Collection of Cell Cultures (Porton Down, Salisbury, U.K.). Cell lines were maintained in 1× RPMI-1640 (Gibco Life Technologies, Paisley, U.K.), supplemented with 10% FCS (Gibco Life Technologies).
RNA extraction and reverse transcription.
Total RNA was extracted from ∼1 × 106 Epstein Barr virus-transformed lymphoblastoid cells using the Perfect RNA Mini RNA kit (Eppendorf, Hamburg, Germany). cDNA was synthesized from 4.5 μg of the total RNA using the Thermoscript RT-PCR system (Gibco Life Technologies) with an incubation temperature of 50°C.
Allele-specific real-time RT-PCR.
Ectopic HNF-1α mRNA transcripts were amplified from lymphoblastoid cell lines in a nested PCR. The first-round reaction was carried out using 2 μl of cDNA in a conventional thermocycler using primers EX5FP and EX10RP, as described previously (21). Ten replicate second-round reactions were carried out for each cell line using 1 μl of a 1:10 dilution of the first-round products with HNF1AFS and HNF1ARS (Table 3) on the Roche Lightcycler generating a cDNA-specific amplicon of 307 bp. PCR products were detected by the use of allele-specific probes that were identical except for the base at the site of the 1375 C>T polymorphism (Table 3). The sensor probes also contained 3′ mismatches (shown underlined) to increase the destabilization of heteroduplexes and hence increase specificity. The reaction included 0.25 μmol/l each of fluorescein-labeled anchor probe Anchor F, LC red-640-labeled sensor C, and LC red-705–labeled sensor T probes (Tib Molbiol, Berlin, Germany) in a total volume of 20 μl (2.5 units) of Taq polymerase (Abgene, Epsom, U.K.) were used per reaction. Lightcycler PCR conditions were initial denaturation 30 s at 90°C followed by 40 cycles of 95°C denaturation for 5 s, 55°C annealing for 10 s, and 72°C elongation for 16 s. The emission spectra of LC red-640 and LC red-705 have a slight overlap, so the two signals were corrected by the Roche Lightcycler color compensation kit (Roche Diagnostics, Lewes, U.K.).
Relative quantitation of mutant cell lines.
Relative quantitation of the two transcripts was carried out using the equation 2−ΔΔCt described by Applied Biosystems (Foster City, CA) (30), which refers to the efficiency of the PCR and the difference between the PCR crossing points (the point at which the signal becomes visible over the background) for each allele normalized to the value for a control cell line. This is an accurate measurement of relative abundance that is independent of all other factors.
Inhibition of NMD.
To determine whether any observed reductions in the amount of the mutant transcript were due to NMD, cells were split into two cultures. Each subculture was treated with either solvent alone (1% DMSO) or cycloheximide dissolved in solvent (Sigma, Poole, U.K.) at a concentration of 100 μg/ml for 4 h. After the incubation, the cells were washed once with PBS (Gibco Life Technologies) and harvested by centrifugation.
Statistical analysis.
Nonparametric methods were used, as the data were not normally distributed about the mean. Summary data are therefore presented as median values (interquartile range). For the validation studies, we examined the data using the Kruskal-Wallis test. This tested for any consistent sources of variation within the validation study as a whole. Pairwise analysis between groups compared the median ΔCT values of mutant and control cell lines by the nonparametric Mann-Whitney test. Statistical analysis was performed using the SSPS package (SSPS, Chicago, IL).
A: Mechanism of fluorescent resonance energy transfer. The head-to-tail alignment of an anchor probe labeled 3′ with fluorescein (FITC) and a sensor probe labeled 5′ with an acceptor dye (LC red-640 or LC red-705), allows the transfer of energy from the donor dye to the acceptor dye (in this example, LC red-640), which then fluoresces as shown by the star. The anchor probe has a higher melting point than the sensor probe, meaning that any variations in fluorescence are due to the binding and melting of the sensor probe. The incoming light from the diode source excites the donor dye, but the acceptor dye does not fluoresce unless both probes are bound. B: Differential melting point of homoduplexes and heteroduplexes. The use of a derivative melting curve allows the measurement of the melting temperature of perfectly matched and mismatched duplexes. Using a probe to allele 1, the perfectly matched allele 1/allele 1 homoduplexes (solid line) will be more stable and thus have a higher melting point than the mismatched allele 1/allele 2 heteroduplexes (dashed line). Therefore, the characteristic melting profile of a mismatched heteroduplex will be shifted to the left on the melting plot.
A: Mechanism of fluorescent resonance energy transfer. The head-to-tail alignment of an anchor probe labeled 3′ with fluorescein (FITC) and a sensor probe labeled 5′ with an acceptor dye (LC red-640 or LC red-705), allows the transfer of energy from the donor dye to the acceptor dye (in this example, LC red-640), which then fluoresces as shown by the star. The anchor probe has a higher melting point than the sensor probe, meaning that any variations in fluorescence are due to the binding and melting of the sensor probe. The incoming light from the diode source excites the donor dye, but the acceptor dye does not fluoresce unless both probes are bound. B: Differential melting point of homoduplexes and heteroduplexes. The use of a derivative melting curve allows the measurement of the melting temperature of perfectly matched and mismatched duplexes. Using a probe to allele 1, the perfectly matched allele 1/allele 1 homoduplexes (solid line) will be more stable and thus have a higher melting point than the mismatched allele 1/allele 2 heteroduplexes (dashed line). Therefore, the characteristic melting profile of a mismatched heteroduplex will be shifted to the left on the melting plot.
Relative quantitation of mutant and wild-type mRNA transcript levels.
. | Control 2 . | R171X . | I414G415ATCG→CCA . | P291fsinsC . | G207D . | R229P . |
---|---|---|---|---|---|---|
Median ΔCt | 0.145 | 1.095 | 18.01 | 3.845 | −0.105 | 0.745 |
Interquartile range | 0–0.46 | 0.37–1.81 | 6.47–20.13 | 0.19–12.08 | 0.81–1.09 | −0.22 to 26.62 |
ΔΔCt | 0.36 | 1.305 | 18.22 | 4.055 | 0.105 | 0.955 |
% Wild type | 78 | 40 | <0.01 | 6 | 93 | 52 |
P value | 0.156 | 0.009 | <0.0001 | 0.001 | 0.546 | 0.03 |
. | Control 2 . | R171X . | I414G415ATCG→CCA . | P291fsinsC . | G207D . | R229P . |
---|---|---|---|---|---|---|
Median ΔCt | 0.145 | 1.095 | 18.01 | 3.845 | −0.105 | 0.745 |
Interquartile range | 0–0.46 | 0.37–1.81 | 6.47–20.13 | 0.19–12.08 | 0.81–1.09 | −0.22 to 26.62 |
ΔΔCt | 0.36 | 1.305 | 18.22 | 4.055 | 0.105 | 0.955 |
% Wild type | 78 | 40 | <0.01 | 6 | 93 | 52 |
P value | 0.156 | 0.009 | <0.0001 | 0.001 | 0.546 | 0.03 |
ΔCt refers to the difference between the crossing point for the C allele and the T allele. ΔΔCt refers to this figure normalized to the controls. The median ΔCt values and interquartile range are shown for control 2, R171X, 1414G415ATCG→CCA, P291fsinsC, R229P, and G207D. The calculated levels of the mutant transcripts are shown with their P values as determined by Mann-Whitney test.
Effect of cycloheximide treatment on mRNA transcript levels.
Treatment . | ΔCt . | ΔΔCt . | % Normal transcript . | P value . |
---|---|---|---|---|
R171X (CHX −) | 1.225 | 1.435 | 37 | 0.002 |
R171X (CHX +) | 0.375 | 0.585 | 67 | 0.185 |
I414G415ATCG→CCA (CHX −) | 18.68 | 18.89 | <0.01 | 0.001 |
I414G415ATCG→CCA (CHX +) | 0.18 | 0.39 | 76 | 0.603 |
P291fsinsC (CHX −) | 5.98 | 6.19 | 2 | 0.001 |
P291fsinsC (CHX +) | −0.38 | −0.17 | 112 | 0.776 |
R229P (CHX −) | 0.775 | 0.985 | 51 | 0.03 |
R229P (CHX +) | 1.245 | 1.455 | 36 | 0.003 |
G207D (CHX −) | 0.02 | 0.23 | 86 | 0.580 |
G207D (CHX +) | −0.58 | −0.37 | 129 | 0.656 |
Treatment . | ΔCt . | ΔΔCt . | % Normal transcript . | P value . |
---|---|---|---|---|
R171X (CHX −) | 1.225 | 1.435 | 37 | 0.002 |
R171X (CHX +) | 0.375 | 0.585 | 67 | 0.185 |
I414G415ATCG→CCA (CHX −) | 18.68 | 18.89 | <0.01 | 0.001 |
I414G415ATCG→CCA (CHX +) | 0.18 | 0.39 | 76 | 0.603 |
P291fsinsC (CHX −) | 5.98 | 6.19 | 2 | 0.001 |
P291fsinsC (CHX +) | −0.38 | −0.17 | 112 | 0.776 |
R229P (CHX −) | 0.775 | 0.985 | 51 | 0.03 |
R229P (CHX +) | 1.245 | 1.455 | 36 | 0.003 |
G207D (CHX −) | 0.02 | 0.23 | 86 | 0.580 |
G207D (CHX +) | −0.58 | −0.37 | 129 | 0.656 |
ΔCt refers to the difference between the crossing point for the C allele and the T allele. ΔΔCt refers to this figure normalized to the controls. Median ΔCt values and SDs are shown for all mutations with and without cycloheximide treatment. The calculated levels of the mutant transcript are given together with their P values as determined by Mann-Whitney test.
Primer and probe sequences used in this study.
Primer/probe name . | Sequence . | Label . |
---|---|---|
HNF1AFS | 5′ GTC CTA CGT TCA CCA ACA CAG 3′ | none |
HNF1ARS | 5′ GTG ATG AGC ATA GTC TGC GG 3′ | none |
Anchor F | 5′ TGT GCC GGT CAT CAA CAG CAT GGG CAG 3′ | 3′ FITC |
Sensor C | 5′ AGC CTG ACC ACA̅ CTG CAG C 3′ | 5′ LC red-640 3′ phosphate |
Sensor T | 5′ AGC CTG ACC ACA̅ TTG CAG C 3′ | 5′ LC red-705 3′ phosphate |
Primer/probe name . | Sequence . | Label . |
---|---|---|
HNF1AFS | 5′ GTC CTA CGT TCA CCA ACA CAG 3′ | none |
HNF1ARS | 5′ GTG ATG AGC ATA GTC TGC GG 3′ | none |
Anchor F | 5′ TGT GCC GGT CAT CAA CAG CAT GGG CAG 3′ | 3′ FITC |
Sensor C | 5′ AGC CTG ACC ACA̅ CTG CAG C 3′ | 5′ LC red-640 3′ phosphate |
Sensor T | 5′ AGC CTG ACC ACA̅ TTG CAG C 3′ | 5′ LC red-705 3′ phosphate |
The inserted mismatches in the probes Sensor C and Sensor T are underlined, and the site of the cSNP is marked in bold type.
Article Information
We are grateful to Diabetes U.K., the European Union (contract number QLG-CT-1999-00546 [GIFT]), and the Wellcome Trust for financial support. A.T.H is a Wellcome Trust Research Leave Fellow.
We would also like to acknowledge Dr. Michael Greco from the Research and Development Support Unit at the University of Exeter and Michael Weedon for statistical advice. We thank Dr. Michael Bulman for help in the initial stages of the study, The European Collection of Cell Cultures for providing the transformed cell lines, and all the patients who contributed to the study.