DOI: 10.2337/db06-0402 © 2007 by the American Diabetes Association
Evaluation of Polymorphic Splicing in the Mechanism of the Association of the Insulin Gene With DiabetesFrom the Endocrine Genetics Laboratory, McGill University Health Center Research Institute, Montréal, Québec, Canada Address correspondence and reprint requests to Constantin Polychronakos, MD, McGill University Health Center, Montreal Childrens Hospital, 2300 Tupper St., Ste. C244, Montreal, Quebec, Canada H3H 1P3. E-mail: constantin.polychronakos{at}mcgill.ca
Abbreviations:
LD, linkage disequilibrium; SNP, single nucleotide polymorphism; UTR, untranslated region; VNTR, variable number of tandem repeats
The association of type 1 diabetes with the insulin gene (IDDM2 locus) has been mapped to a short haplotype encompassing two single nucleotide polymorphisms (SNPs) in perfect linkage disequilibrium (r2 = 1) with each other and with the two allele classes at the variable number of tandem repeats (VNTR) polymorphism upstream of the transcription site. Although it is believed that the genetic effect is mediated through transcriptional effects of the VNTR, an alternative mechanism has been recently proposed: In transfected cells, the common A allele at one of the SNPs (–23A T, in relation to the translation-initiation codon) weakens the splicing of intron 1, resulting in a minor ( 15% of total RNA) transcript with a longer 5' untranslated region and sixfold enhanced translational efficiency. The purpose of our study was to confirm these findings in RNA from normal human pancreas and thymus. We report that pancreas does contain the alternative transcript in an allele-dependent manner but at a very low proportion (<5% of total INS mRNA). We believe that this level would have a minor, if any, biological effect involved in the mechanism of the IDDM2 locus.
A haplotype within a 2-kb linkage disequilibrium (LD) block encompassing the human insulin gene (INS) is strongly associated with type 1 diabetes (1–3) and, possibly, also with plasma insulin levels (4), juvenile obesity (5), and type 2 diabetes (6). The association maps to a haplotype defined by a repeat polymorphism upstream of INS that consists of variable number of tandem repeats (VNTR). The short class I alleles predispose to type 1 diabetes, whereas the long class III alleles have a dominantly protective effect (2,3). However, because of tight LD, the effect cannot be genetically dissected from two single nucleotide polymorphisms (SNPs), rs689 (–23A
In the absence of coding polymorphisms, the association must be due to allelic effects on transcription, splicing, and/or RNA stability. In pancreas, steady-state mRNA derived from class III chromosomes is indeed
Exon 1 of INS is untranslated and the entire coding sequence is contained in exons 2 and 3 (Fig. 1). SNP –23A
However, splicing is tissue specific, and isoforms found in unrelated cells, or even tumors or cell lines derived from the tissue of interest, may be irrelevant. As an example, insulinomas preferentially use a cryptic splicing site that results in an isoform retaining 22 nucleotides of intron 1 (Fig. 1), only traces of which are found in normal ß-cells (13). Interestingly, this cryptic splicing site is activated in the presence of a 4-bp insertion (12) and is frequent in Africans ( 25%) but absent in the European-descent populations (14) in which the IDDM2 effect was observed. The aims of this study were to confirm the allelic effect in RNA from human pancreas and thymus and to assess its quantitative importance.
We first examined the available human expressed sequence tag (EST) database (Online Table 1, available at http://dx.doi.org/10.2337/db06-0402). Besides the two libraries analyzed in ref. 12, there were four more containing at least five clones informative for splicing and having at least one A allele at –23A Transcripts splicing out part of coding sequences, observed in vitro (transcripts 1, 2, 3, and 5 of ref. 12) were rare—9 of 2,009 ESTs—underscoring the pitfall in extrapolating from in vitro studies. In our examination of RNA from pancreas and thymus, we first confirmed the allele specificity of intron 1 retention by use of RT-PCR primers that specifically amplified mRNA retaining intron 1 (transcript 6 in Fig. 1, abbreviated T6 to maintain the numbering of ref. 12). All samples showed amplification product (Fig. 2, top panel), and genotype had no apparent effect on band intensity. However, because PCR was carried beyond the logarithmic phase, quantitative comparison between samples is not meaningful; alleles can still be compared within heterozygous samples (competitive RT-PCR). Thus, in AT heterozygotes, relative allele quantification by single nucleotide primer extension clearly showed that most of the T6 transcripts carry the A allele. After correction for labeling and fluorochrome efficiency using heterozygous DNA as standard for 1:1 stoichiometry, the A allele was five- to sixfold more abundant than the T allele (mean ± SD = 5.7 ± 1.7–fold in four samples). Therefore, we confirm the observation (12) that intron 1 retention is driven mostly by A alleles.
Given the allele dependence of splicing, we then asked whether our previous finding of allelic imbalance in the thymus might be due to the presence of different, allele-dependent isoforms in the 3' UTR RT-PCR used to demonstrate it. To answer this question, we amplified transcripts containing all three exons, using primers S1 and AS1 from a thymus sample heterozygous for –23A T (Fig. 3). The band corresponding to the 451-bp canonical transcript T4 was gel-purified and examined for relative expression of the two alleles. There was no band corresponding to T6 (Fig. 3). Because –23A T is spliced out in these transcripts, we used two 3' UTR SNPs as proxies. In +1140A C, allele C is in perfect (r2 = 1) LD with VNTR class I (A at –23) and A with VNTR class III (T at –23) (3). In +1127C T (rs3842752), allele T invariably predicts a class III VNTR (T at –23) (3,8) and, in individuals heterozygous at both –23A T and +1,127 C T, C predicts class I-A with virtual certainty (99.93% of cases, see METHODS in ref. 8). By both SNPs, transcripts from the class III chromosome are approximately twice as abundant as those from the class I chromosome. Therefore, the allelic imbalance is present in gel-purified T4 transcripts and cannot be due to differential presence of allele-dependent splicing isoforms. A VNTR transcriptional effect remains the most compelling explanation. To confirm the absence of the T6 transcript in thymus, we amplified an additional five RNA samples using primers S1 and AS1. To maximize visualization of faint bands, RT-PCR products were subjected to Southern blot (Fig. 4). T4, the canonical transcript, was seen in all samples, but no band corresponding to the intron-containing T6 isoform was present. Therefore, intron 1 retention in thymus, if it occurs at all, is quantitatively insignificant.
Finally, to address the question of whether the enhanced translational efficiency of T6 may be contributing to the IDDM2 effect, we defined its relative representation within total pancreatic mRNA, which should determine the proportion of insulin peptide synthesis attributable to this isoform. Fig. 5 shows the results of competitive RT-PCR, which we used as the most direct way of defining the relative quantitative contribution of different isoforms amplified with the same primers in the same reaction. Full-length INS cDNA was amplified using primers S1 and AS1 (Fig. 1) from eight pancreas RNA samples (Fig. 5A). Transcript identity as inferred from size was confirmed in all cases by extracting and sequencing PCR bands. The canonical transcript (T4) is the only clearly visible band. No band corresponding to the size of T6 can be seen. The extra band seen in the seventh sample is due to retention of 22 bp at the 5' end of intron 1, resulting from activation of the cryptic splice site (12) in an individual carrying the TTGC insertion (12) (confirmed by sequencing). Thus, we confirm the effect of +TTGC on splicing as previously reported (12). This insertion, common in Africans ( 25%), is never seen in Europeans, in whom the IDDM2 effect was observed. Association of this SNP with type 1 diabetes is worth studying in Africans, but it cannot explain the IDDM2 effect observed in Europeans.
In a further attempt to co-amplify T6 with T4, RNA from six human pancreas samples was RT-PCR amplified with primers S1 and AS3, going from exon 1 to exon 2 (Fig. 5B). Because the sequence of RT-PCR products retaining intron 1 is identical with genomic DNA, we are also showing PCR products from the same RNA samples, processed identically except for the omission of reverse transcriptase. A distinct band corresponding to T6 is seen only in the second sample, which is clearly contaminated with genomic DNA. Of the remaining five samples, two AA homozygotes show no T6 bands, while an AT heterozygote has a very faint band. Thus, the relative contribution of T6 to the total amount of insulin mRNA is very low, even if multiplied by a sixfold greater translational efficiency. However, because T6 is longer than T4 and may be amplified with lower efficiency, we had to evaluate the sensitivity of our technique to detect T6 at low proportions in competitive PCR. Figure 5C shows PCR amplification of templates generated by spiking T4 DNA with T6 at the proportions indicated. We easily detected T6 at 5% of total transcripts. This allows the conclusion that T6, when it is present at all, is well below 5% of total transcript, much lower than might be suggested by splicing in cultured cells (12).
Therefore, we confirm the allelic effect on splicing described in ref. 12. We also show that the in vitro system used in that report seriously overestimated the abundance of the allele-dependent alternative transcript and, hence, its biological significance. If in a few individuals T6 represents as much as 3–4% of INS mRNA (the most generous estimate from our data), the enhanced translation from T6 would result in a It would also be interesting to speculate on the immune consequences, in view of its genetic behavior, of increased insulin synthesis associated with the A allele. Enhancing expression of the ß-cell autoantigen (insulin) might promote the type 1 diabetes autoimmune process, as is compatible with the predisposing effect of the haplotype. At the level of the thymus, however, where allelic differences may be more important, higher levels would result in better insulin-specific self-tolerance. This is difficult to reconcile with the increased diabetes risk conferred by the A allele haplotype. We believe that our findings place this allelic alternative splicing in quantitative perspective and suggest minimal effect, if any at all.
Analysis of INS libraries. INS EST libraries were accessed at the National Center for Biotechnology Information (NCBI) Web site (http://www.ncbi.nlm.nih.gov/UniGene/lbrowse2.cgi?TAXID=9606).
Tissue samples and RNA preparation.
Primers and PCR conditions. The primers were as follows: For all isoforms: S1 5'ATCAGAAGAGGCCATCAAGC3'; AS1 5'TTCCATCTCTCTCGGTGCAG3'. For isoform 6 specific: S2 5'GAAGCATGTGGGGGTGAG3'; AS2 5'CACAATGCCACGCTTCTG3'. For the shorter PCR product of T4 and T6, S1 with this antisense primer: AS3 5'CCCCGCACACTAGGTAGAGA3'. PCR cycling conditions were 35 cycles at 94°C for 20 s, 54°C for 20 s, and 72°C for 30 s.
Allelic imbalance measurements.
This work was supported by a grant to C.P. from Genome Québec and Genome Canada. We thank Dr. Cindy G. Goodyer for providing the tissue samples.
Additional information for this article can be found in an online appendix at http://dx.doi.org/10.2337/db06-0402. The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact. Received for publication March 26, 2006 and accepted in revised form December 18, 2006
This article has been cited by other articles:
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||