Variants in the human insulin gene that affect pre-mRNA splicing: is -23HphI a functional single nucleotide polymorphism at IDDM2?

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

Shifts in isoform usage underlie transcriptional differences in regulatory T cells in type 1 diabetes

Genome-wide association studies have identified numerous loci with allelic associations to Type 1 Diabetes (T1D) risk. Most disease-associated variants are enriched in regulatory sequences active in lymphoid cell types, suggesting that lymphocyte gene expression is altered in T1D. Here we assay gene expression between T1D cases and healthy controls in two autoimmunity-relevant lymphocyte cell types, memory CD4+/CD25+ regulatory T cells (Treg) and memory CD4+/CD25- T cells, using a splicing event-based approach to characterize tissue-specific transcriptomes. Limited differences in isoform usage between T1D cases and controls are observed in memory CD4+/CD25- T-cells. In Tregs, 402 genes demonstrate differences in isoform usage between cases and controls, particularly RNA recognition and splicing factor genes. Many of these genes are regulated by the variable inclusion of exons that can trigger nonsense mediated decay. Our results suggest that dysregulation of gene expression, through shifts in alternative splicing in Tregs, contributes to T1D pathophysiology.

Presence of immunogenic alternatively spliced insulin gene product in human pancreatic delta cells

AIMS/HYPOTHESIS

Transcriptome analyses revealed insulin-gene-derived transcripts in non-beta endocrine islet cells. We studied alternative splicing of human INS mRNA in pancreatic islets.

METHODS

Alternative splicing of insulin pre-mRNA was determined by PCR analysis performed on human islet RNA and single-cell RNA-seq analysis. Antisera were generated to detect insulin variants in human pancreatic tissue using immunohistochemistry, electron microscopy and single-cell western blot to confirm the expression of insulin variants. Cytotoxic T lymphocyte (CTL) activation was determined by MIP-1β release.

RESULTS

We identified an alternatively spliced INS product. This variant encodes the complete insulin signal peptide and B chain and an alternative C-terminus that largely overlaps with a previously identified defective ribosomal product of INS. Immunohistochemical analysis revealed that the translation product of this INS-derived splice transcript was detectable in somatostatin-producing delta cells but not in beta cells; this was confirmed by light and electron microscopy. Expression of this alternatively spliced INS product activated preproinsulin-specific CTLs in vitro. The exclusive presence of this alternatively spliced INS product in delta cells may be explained by its clearance from beta cells by insulin-degrading enzyme capturing its insulin B chain fragment and a lack of insulin-degrading enzyme expression in delta cells.

CONCLUSIONS/INTERPRETATION

Our data demonstrate that delta cells can express an INS product derived from alternative splicing, containing both the diabetogenic insulin signal peptide and B chain, in their secretory granules. We propose that this alternative INS product may play a role in islet autoimmunity and pathology, as well as endocrine or paracrine function or islet development and endocrine destiny, and transdifferentiation between endocrine cells. INS promoter activity is not confined to beta cells and should be used with care when assigning beta cell identity and selectivity.

DATA AVAILABILITY

The full EM dataset is available via www.nanotomy.org (for review: http://www.nanotomy.org/OA/Tienhoven2021SUB/6126-368/ ). Single-cell RNA-seq data was made available by Segerstolpe et al [13] and can be found at https://sandberglab.se/pancreas . The RNA and protein sequence of INS-splice was uploaded to GenBank (BankIt2546444 INS-splice OM489474).

Effects of 14 F9 synonymous codon variants on hemophilia B expression: Alteration of splicing along with protein expression

There is growing evidence that synonymous codon variants (SCVs) can cause disease through the disruption of different processes of protein production. The aim of the study is to investigate whether the 14 SCVs reported in the F9 variant database were the pathogenic causes of hemophilia B. The impacts of SCVs on splicing and protein expression were detected using a combination of in silico prediction, in vitro minigene splicing assay and cell expression detection. The splicing transcripts were identified and quantified by co-amplification fluorescent PCR. The mechanism of splicing was verified by a modified pU1snRNA and pU7snRNA approach. Aberrant splicing patterns were found in eight SCVs. Five of the 8 SCVs produced almost all aberrant splicing isoforms, which were expected to truncate protein, three of them presented a partial defect on both splicing and protein secretion, the overall effects were consistent with the residual Factor IX activity of the affected cases. Neither the pre-messenger RNA (mRNA) splicing process nor the protein function was impaired in the rest six SCVs. In conclusion, our study firstly revealed the pathogenic mechanism of the 14 F9 SCVs and highlighted the importance of performing mRNA splicing analysis and protein expression studies of SCVs in inherited disorders.

Regulatory roles and mechanisms of alternative RNA splicing in adipogenesis and human metabolic health

Alternative splicing (AS) regulates gene expression patterns at the post-transcriptional level and generates a striking expansion of coding capacities of genomes and cellular protein diversity. RNA splicing could undergo modulation and close interaction with genetic and epigenetic machinery. Notably, during the adipogenesis processes of white, brown and beige adipocytes, AS tightly interplays with the differentiation gene program networks. Here, we integrate the available findings on specific splicing events and distinct functions of different splicing regulators as examples to highlight the directive biological contribution of AS mechanism in adipogenesis and adipocyte biology. Furthermore, accumulating evidence has suggested that mutations and/or altered expression in splicing regulators and aberrant splicing alterations in the obesity-associated genes are often linked to humans’ diet-induced obesity and metabolic dysregulation phenotypes. Therefore, significant attempts have been finally made to overview novel detailed discussion on the prospects of splicing machinery with obesity and metabolic disorders to supply featured potential management mechanisms in clinical applicability for obesity treatment strategies.

DNA methylation near the INS gene is associated with INS genetic variation (rs689) and type 1 diabetes in the Diabetes Autoimmunity Study in the Young

OBJECTIVE

Mechanisms underlying the role of non-human leukocyte antigen (HLA) genetic risk variants in type 1 diabetes (T1D) are poorly understood. We aimed to test the association between methylation and non-HLA genetic risk.

METHODS

We conducted a methylation quantitative trait loci (mQTL) analysis in a nested case-control study from the Dietary Autoimmunity Study in the Young. Controls (n = 83) were frequency-matched to T1D cases (n = 83) based on age, race/ethnicity, and sample availability. We evaluated 13 non-HLA genetic markers known be associated with T1D. Genome-wide methylation profiling was performed on peripheral blood samples collected prior to T1D using the Illumina 450 K (discovery set) and infinium methylation EPIC beadchip (EPIC validation) platforms. Linear regression models, adjusting for age and sex, were used to test to each single nucleotide polymorphism (SNP) -probe combination. Logistic regression models were used to test the association between T1D and methylation levels among probes with a significant mQTL. A meta-analysis was used to combine odds ratios from the two platforms.

RESULTS

We identified 10 SNP-methylation probe pairs (false discovery rate (FDR) adjusted P < .05 and validation P < .05). Probes were associated with the GSDMB, C1QTNF6, IL27, and INS genes. The cg03366382 (OR: 1.9, meta-P = .0495), cg21574853 (OR: 2.5, meta-P = .0232), and cg25336198 (odds ratio: 6.6, meta-P = .0081) probes were significantly associated with T1D. The three probes were located upstream from the INS transcription start site.

CONCLUSIONS

We confirmed an association between DNA methylation and rs689 that has been identified in related studies. Measurements in our study preceded the onset of T1D suggesting methylation may have a role in the relationship between INS variation and T1D development.

PUF60-activated exons uncover altered 3' splice-site selection by germline missense mutations in a single RRM

PUF60 is a splicing factor that binds uridine (U)-rich tracts and facilitates association of the U2 small nuclear ribonucleoprotein with primary transcripts. PUF60 deficiency (PD) causes a developmental delay coupled with intellectual disability and spinal, cardiac, ocular and renal defects, but PD pathogenesis is not understood. Using RNA-Seq, we identify human PUF60-regulated exons and show that PUF60 preferentially acts as their activator. PUF60-activated internal exons are enriched for Us upstream of their 3′ splice sites (3′ss), are preceded by longer AG dinucleotide exclusion zones and more distant branch sites, with a higher probability of unpaired interactions across a typical branch site location as compared to control exons. In contrast, PUF60-repressed exons show U-depletion with lower estimates of RNA single-strandedness. We also describe PUF60-regulated, alternatively spliced isoforms encoding other U-bound splicing factors, including PUF60 partners, suggesting that they are co-regulated in the cell, and identify PUF60-regulated exons derived from transposed elements. PD-associated amino-acid substitutions, even within a single RNA recognition motif (RRM), altered selection of competing 3′ss and branch points of a PUF60-dependent exon and the 3′ss choice was also influenced by alternative splicing of PUF60. Finally, we propose that differential distribution of RNA processing steps detected in cells lacking PUF60 and the PUF60-paralog RBM39 is due to the RBM39 RS domain interactions. Together, these results provide new insights into regulation of exon usage by the 3′ss organization and reveal that germline mutation heterogeneity in RRMs can enhance phenotypic variability at the level of splice-site and branch-site selection.

Thioflavin T Monitoring of Guanine Quadruplex Formation in the rs689-Dependent INS Intron 1

The human proinsulin gene (INS) contains a thymine-to-adenine variant (rs689) located in the 3′ splice site (3′ ss) recognition motif of the first intron. The adenine at rs689 is strongly associated with type 1 diabetes. By weakening the polypyrimidine tract, the adenine allele reduces the efficiency of intron 1 splicing, which can be ameliorated by antisense oligonucleotides blocking a splicing silencer located upstream of the 3′ ss. The silencer is surrounded by guanine-rich tracts that may form guanine quadruplexes (G4s) and modulate the accessibility of the silencer. Here, we employed thioflavin T (ThT) to monitor G4 formation in synthetic DNAs and RNAs derived from INS intron 1. We show that the antisense target is surrounded by ThT-positive segments in each direction, with oligoribonucleotides exhibiting consistently higher fluorescence than their DNA counterparts. The signal was reduced for ThT-positive oligonucleotides that were extended into the silencer, indicating that flanking G4s have a potential to mask target accessibility. Real-time monitoring of ThT fluorescence during INS transcription in vitro revealed a negative correlation with ex vivo splicing activities of corresponding INS constructs. Together, these results provide a better characterization of antisense targets in INS primary transcripts for restorative strategies designed to improve the INS splicing defect associated with type 1 diabetes.

Homogentisate 1,2-dioxygenase (HGD) gene variants, their analysis and genotype-phenotype correlations in the largest cohort of patients with AKU

Alkaptonuria (AKU) is a rare metabolic disorder caused by a deficient enzyme in the tyrosine degradation pathway, homogentisate 1,2-dioxygenase (HGD). In 172 AKU patients from 39 countries, we identified 28 novel variants of the HGD gene, which include three larger genomic deletions within this gene discovered via self-designed multiplex ligation-dependent probe amplification (MLPA) probes. In addition, using a reporter minigene assay, we provide evidence that three of eight tested variants potentially affecting splicing cause exon skipping or cryptic splice-site activation. Extensive bioinformatics analysis of novel missense variants, and of the entire HGD monomer, confirmed mCSM as an effective computational tool for evaluating possible enzyme inactivation mechanisms. For the first time for AKU, a genotype-phenotype correlation study was performed for the three most frequent HGD variants identified in the Suitability Of Nitisinone in Alkaptonuria 2 (SONIA2) study. We found a small but statistically significant difference in urinary homogentisic acid (HGA) excretion, corrected for dietary protein intake, between variants leading to 1% or >30% residual HGD activity. There was, interestingly, no difference in serum levels or absolute urinary excretion of HGA, or clinical symptoms, indicating that protein intake is more important than differences in HGD variants for the amounts of HGA that accumulate in the body of AKU patients.

RNASEH1 gene variants are associated with autoimmune type 1 diabetes in Colombia

BACKGROUND

In a previous work, we found linkage and association of type 1 diabetes (T1D) to a 12 known gene region at chromosome 2p25 in Colombian families. Here, we present further work on this candidate region.

MATERIALS AND METHODS

Seventeen SNPs located on the 12 candidate genes, in 100 familial trios set, were tested by ARMS-tetraprimer-PCR or PCR-RFLP. Five extra SNPs in the vicinity of rs10186193 were typed. A replica phase included 97 novel familial trios, in whom diabetes-related auto-antibodies (AABs) were tested in sera of the patients. In addition to transmission disequilibrium tests, haplotype analyses were carried out using the unphased software.

RESULTS

SNP rs10186193 (at RNASEH1 gene) showed association with T1D (P = 0.005). The additional five SNPs revealed that rs7607888 (P = 2.03 × 10^-7), rs55981318 (P = 0.018), and rs1136545 (P = 1.93 × 10^-9) were also associated with T1D. Haplotype analysis showed association for rs55981318-rs10186193 (P = 0.0005), rs7563960-rs7607888 (P = 0.0007), rs7607888-rs1136545 (P = 9.21 × 10^-10), and rs1136545-rs11538545 (P = 6.67 × 10^-8). In contrast, the new set of 97 familial trios tested for SNPs rs55981318, rs10186193, and rs7607888 did not support the previous finding; however, by combining the sample (197 trios), evidence of association of T1D with rs55981318 and rs7607888 was conclusive. In addition, a two-loci haplotype analysis of the combined sample showed significant association of RNASEH1 with T1D (P = 3.1 × 10^-5).

CONCLUSION

In conclusion, our analyses suggest that RNASEH1 gene variants associate with susceptibility/protection to T1D in Colombia.

Disease-specific biases in alternative splicing and tissue-specific dysregulation revealed by multitissue profiling of lymphocyte gene expression in type 1 diabetes

Genome-wide association studies (GWAS) have identified multiple, shared allelic associations with many autoimmune diseases. However, the pathogenic contributions of variants residing in risk loci remain unresolved. The location of the majority of shared disease-associated variants in noncoding regions suggests they contribute to risk of autoimmunity through effects on gene expression in the immune system. In the current study, we test this hypothesis by applying RNA sequencing to CD4+, CD8+, and CD19+ lymphocyte populations isolated from 81 subjects with type 1 diabetes (T1D). We characterize and compare the expression patterns across these cell types for three gene sets: all genes, the set of genes implicated in autoimmune disease risk by GWAS, and the subset of these genes specifically implicated in T1D. We performed RNA sequencing and aligned the reads to both the human reference genome and a catalog of all possible splicing events developed from the genome, thereby providing a comprehensive evaluation of the roles of gene expression and alternative splicing (AS) in autoimmunity. Autoimmune candidate genes displayed greater expression specificity in the three lymphocyte populations relative to other genes, with significantly increased levels of splicing events, particularly those predicted to have substantial effects on protein isoform structure and function (e.g., intron retention, exon skipping). The majority of single-nucleotide polymorphisms within T1D-associated loci were also associated with one or more cis-expression quantitative trait loci (cis-eQTLs) and/or splicing eQTLs. Our findings highlight a substantial, and previously underrecognized, role for AS in the pathogenesis of autoimmune disorders and particularly for T1D.

Evidence of selection on splicing-associated loci in human populations and relevance to disease loci mapping

We performed a whole-genome scan of genetic variants in splicing regulatory elements (SREs) and evaluated the extent to which natural selection has shaped extant patterns of variation in SREs. We investigated the degree of differentiation of single nucleotide polymorphisms (SNPs) in SREs among human populations and applied long-range haplotype- and multilocus allelic differentiation-based methods to detect selection signatures. We describe an approach, sampling a large number of loci across the genome from functional classes and using the consensus from multiple tests, for identifying candidates for selection signals. SRE SNPs in various SNP functional classes show different patterns of population differentiation compared with their non-SRE counterparts. Intronic regions display a greater enrichment for extreme population differentiation among the potentially tissue-dependent transcript ratio quantitative trait loci (trQTLs) than SRE SNPs in general and includ outlier trQTLs for cross-population composite likelihood ratio, suggesting that incorporation of context annotation for regulatory variation may lead to improved detection of signature of selection on these loci. The proportion of extremely rare SNPs disrupting SREs is significantly higher in European than in African samples. The approach developed here will be broadly useful for studies of function and disease-associated variation in the human genome.

Single-nucleotide polymorphism of INS, INSR, IRS1, IRS2, PPAR-G and CAPN10 genes in the pathogenesis of polycystic ovary syndrome

Polycystic ovary syndrome (PCOS) is the most common and a complex female endocrine disorder, and is one of the leading cause of female infertility. Here, we aimed to investigate the association of single-nucleotide polymorphism of INS, INSR, IRS1, IRS2, PPAR-G and CAPN10 gene in the pathogenesis of PCOS. A hospital-based, observational case-control study was carried on 169 PCOS and 169 control women in the southern region of India. Genotype was carried out by real-time polymerase chain reaction. A chi-square (χ²) test was performed and the genotypes were verified to comply with the Hardy- Weinberg equilibrium. Odds ratio and 95% confidence interval were calculated to assess the relative risk. Comparison of clinical characteristics of women with PCOS and controls reveal an increase in body mass index (BMI), luteinizing hormone / follicle stimulating hormone (LH/FSH) ratio, glucose levels, insulin, testosterone, hirsutism and antral follicular count in PCOS women. The variant rs1801278 (P = 0.002; OR = 2.88; 95% CI = 1.43, 5.80) show an association with PCOS. In the genotypic (P = 0.0002) and allelic models (P = 0.000), significance persisted even after Bonferroni correction. The genotypes of SNPs strongly influence BMI, LH, LH/FSH ratio, ovarian volume and antral follicular count in PCOS women. The study results were suggestive of a positive association between Gly972Arg of IRS1 and PCOS in the south Indian population, while INS, IRS2, PPAR-G and CAPN10 failed to show any association with PCOS in our studied population. Further studies focussing the role of IRS1 are warranted to delineate its implication towards PCOS.

Alternative splicing of U2AF1 reveals a shared repression mechanism for duplicated exons

The auxiliary factor of U2 small nuclear ribonucleoprotein (U2AF) facilitates branch point (BP) recognition and formation of lariat introns. The gene for the 35-kD subunit of U2AF gives rise to two protein isoforms (termed U2AF35a and U2AF35b) that are encoded by alternatively spliced exons 3 and Ab, respectively. The splicing recognition sequences of exon 3 are less favorable than exon Ab, yet U2AF35a expression is higher than U2AF35b across tissues. We show that U2AF35b repression is facilitated by weak, closely spaced BPs next to a long polypyrimidine tract of exon Ab. Each BP lacked canonical uridines at position -2 relative to the BP adenines, with efficient U2 base-pairing interactions predicted only for shifted registers reminiscent of programmed ribosomal frameshifting. The BP cluster was compensated by interactions involving unpaired cytosines in an upstream, EvoFold-predicted stem loop (termed ESL) that binds FUBP1/2. Exon Ab inclusion correlated with predicted free energies of mutant ESLs, suggesting that the ESL operates as a conserved rheostat between long inverted repeats upstream of each exon. The isoform-specific U2AF35 expression was U2AF65-dependent, required interactions between the U2AF-homology motif (UHM) and the α6 helix of U2AF35, and was fine-tuned by exon Ab/3 variants. Finally, we identify tandem homologous exons regulated by U2AF and show that their preferential responses to U2AF65-related proteins and SRSF3 are associated with unpaired pre-mRNA segments upstream of U2AF-repressed 3′ss. These results provide new insights into tissue-specific subfunctionalization of duplicated exons in vertebrate evolution and expand the repertoire of exon repression mechanisms that control alternative splicing.

Identification of insulin gene variants in neonatal diabetes

OBJECTIVES

Permanent neonatal diabetes (PNDM) is caused by mutations in the genes responsible for the synthesis of different proteins that are important for the normal behavior of beta cells in the pancreas. Mutations in the insulin gene (INS) are considered as one of the causes of diabetes in neonates. This study aimed to investigate the genetic variations in the INS gene in a group of Egyptian infants diagnosed with PNDM.

METHODS

We screened exons 2 and 3 with intronic boundaries of the INS gene by direct gene sequencing in 30 PNDM patients and 20 healthy controls. A detailed clinical phenotyping of the patients was carried out to specify the diabetes features in those found to carry an INS variant.

RESULTS

We identified five variants (four SNPs and one synonymous variant), c(0).187 + 11T > C, c.-17-6T > A, c.*22A > C, c.*9C > T, and c.36G > A (p.A12A), with allelic frequencies of 96.7%, 80%, 75%, 5%, and 1.7%, respectively. All showed no statistically significance difference compared with the controls, with the exception of c.*22A > C.

CONCLUSION

Genetic screening for the INS gene did not reveal an evident role in the diagnosis of PNDM.

Identification of U2AF(35)-dependent exons by RNA-Seq reveals a link between 3' splice-site organization and activity of U2AF-related proteins

The auxiliary factor of U2 small nuclear RNA (U2AF) is a heterodimer consisting of 65- and 35-kD proteins that bind the polypyrimidine tract (PPT) and AG dinucleotides at the 3′ splice site (3′ss). The gene encoding U2AF35 (U2AF1) is alternatively spliced, giving rise to two isoforms U2AF35a and U2AF35b. Here, we knocked down U2AF35 and each isoform and characterized transcriptomes of HEK293 cells with varying U2AF35/U2AF65 and U2AF35a/b ratios. Depletion of both isoforms preferentially modified alternative RNA processing events without widespread failure to recognize 3′ss or constitutive exons. Over a third of differentially used exons were terminal, resulting largely from the use of known alternative polyadenylation (APA) sites. Intronic APA sites activated in depleted cultures were mostly proximal whereas tandem 3′UTR APA was biased toward distal sites. Exons upregulated in depleted cells were preceded by longer AG exclusion zones and PPTs than downregulated or control exons and were largely activated by PUF60 and repressed by CAPERα. The U2AF(35) repression and activation was associated with a significant interchange in the average probabilities to form single-stranded RNA in the optimal PPT and branch site locations and sequences further upstream. Although most differentially used exons were responsive to both U2AF subunits and their inclusion correlated with U2AF levels, a small number of transcripts exhibited distinct responses to U2AF35a and U2AF35b, supporting the existence of isoform-specific interactions. These results provide new insights into function of U2AF and U2AF35 in alternative RNA processing.

Optimal antisense target reducing INS intron 1 retention is adjacent to a parallel G quadruplex

Splice-switching oligonucleotides (SSOs) have been widely used to inhibit exon usage but antisense strategies that promote removal of entire introns to increase splicing-mediated gene expression have not been developed. Here we show reduction of INS intron 1 retention by SSOs that bind transcripts derived from a human haplotype expressing low levels of proinsulin. This haplotype is tagged by a polypyrimidine tract variant rs689 that decreases the efficiency of intron 1 splicing and increases the relative abundance of mRNAs with extended 5' untranslated region (5' UTR), which curtails translation. Co-expression of haplotype-specific reporter constructs with SSOs bound to splicing regulatory motifs and decoy splice sites in primary transcripts revealed a motif that significantly reduced intron 1-containing mRNAs. Using an antisense microwalk at a single nucleotide resolution, the optimal target was mapped to a splicing silencer containing two pseudoacceptor sites sandwiched between predicted RNA guanine (G) quadruplex structures. Circular dichroism spectroscopy and nuclear magnetic resonance of synthetic G-rich oligoribonucleotide tracts derived from this region showed formation of a stable parallel 2-quartet G-quadruplex on the 3' side of the antisense retention target and an equilibrium between quadruplexes and stable hairpin-loop structures bound by optimal SSOs. This region interacts with heterogeneous nuclear ribonucleoproteins F and H that may interfere with conformational transitions involving the antisense target. The SSO-assisted promotion of weak intron removal from the 5' UTR through competing noncanonical and canonical RNA structures may facilitate development of novel strategies to enhance gene expression.

Glial cell line-derived neurotrophic factor (GDNF) as a novel candidate gene of anxiety

Glial cell line-derived neurotrophic factor (GDNF) is a neurotrophic factor for dopaminergic neurons with promising therapeutic potential in Parkinson's disease. A few association analyses between GDNF gene polymorphisms and psychiatric disorders such as schizophrenia, attention deficit hyperactivity disorder and drug abuse have also been published but little is known about any effects of these polymorphisms on mood characteristics such as anxiety and depression. Here we present an association study between eight (rs1981844, rs3812047, rs3096140, rs2973041, rs2910702, rs1549250, rs2973050 and rs11111) GDNF single nucleotide polymorphisms (SNPs) and anxiety and depression scores measured by the Hospital Anxiety and Depression Scale (HADS) on 708 Caucasian young adults with no psychiatric history. Results of the allele-wise single marker association analyses provided significant effects of two single nucleotide polymorphisms on anxiety scores following the Bonferroni correction for multiple testing (p = 0.00070 and p = 0.00138 for rs3812047 and rs3096140, respectively), while no such result was obtained on depression scores. Haplotype analysis confirmed the role of these SNPs; mean anxiety scores raised according to the number of risk alleles present in the haplotypes (p = 0.00029). A significant sex-gene interaction was also observed since the effect of the rs3812047 A allele as a risk factor of anxiety was more pronounced in males. In conclusion, this is the first demonstration of a significant association between the GDNF gene and mood characteristics demonstrated by the association of two SNPs of the GDNF gene (rs3812047 and rs3096140) and individual variability of anxiety using self-report data from a non-clinical sample.

Posttranscriptional regulation of insulin family ligands and receptors

Insulin system including ligands (insulin and IGFs) and their shared receptors (IR and IGFR) are critical regulators of insulin signaling and glucose homeostasis. Altered insulin system is associated with major pathological conditions like diabetes and cancer. The mRNAs encoding for these ligands and their receptors are posttranscriptionally controlled by three major groups of regulators; (i) alternative splicing regulatory factors; (ii) turnover and translation regulator RNA-binding proteins (TTR-RBPs); and (iii) non-coding RNAs including miRNAs and long non-coding RNAs (lncRNAs). In this review, we discuss the influence of these regulators on alternative splicing, mRNA stability and translation. Due to the pathological impacts of insulin system, we also discussed the possibilities of discovering new potential regulators which will improve understanding of insulin system and associated diseases.

Very low PSA concentrations and deletions of the KLK3 gene

BACKGROUND

Prostate-specific antigen (PSA), a widely used biomarker for prostate cancer (PCa), is encoded by a kallikrein gene (KLK3, kallikrein-related peptidase 3). Serum PSA concentrations vary in the population, with PCa patients generally showing higher PSA concentrations than control individuals, although a small proportion of individuals in the population display very low PSA concentrations. We hypothesized that very low PSA concentrations might reflect gene-inactivating mutations in KLK3 that lead to abnormally reduced gene expression.

METHODS

We have sequenced all KLK3 exons and the promoter and searched for gross deletions or duplications in KLK3 in the 30 individuals with the lowest observed PSA concentrations in a sample of approximately 85 000 men from the Prostate Testing for Cancer and Treatment (ProtecT) study. The ProtecT study examines a community-based population of men from across the UK with little prior PSA testing.

RESULTS

We observed no stop codons or frameshift mutations, but we did find 30 single-base genetic variants, including 3 variants not described previously. These variants included missense variants that could be functionally inactivating and splicing variants. At this stage, however, we cannot confidently conclude whether these variants markedly lower PSA concentration or activity. More importantly, we identified 3 individuals with different large heterozygous deletions that encompass all KLK3 exons. The absence of a functional copy of KLK3 in these individuals is consistent with their reduced serum PSA concentrations.

CONCLUSIONS

The clinical interpretation of the PSA test for individuals with KLK3 gene inactivation could lead to false-negative PSA findings used for screening, diagnosis, or monitoring of PCa.

Both polymorphic variable number of tandem repeats and autoimmune regulator modulate differential expression of insulin in human thymic epithelial cells

OBJECTIVE

Polymorphic INS-VNTR plays an important role in regulating insulin transcript expression in the human thymus that leads to either insulin autoimmunity or tolerance. The molecular mechanisms underlying the INS-VNTR haplotype-dependent insulin expression are still unclear. In this study, we determined the mechanistic components underlying the differential insulin gene expression in human thymic epithelial cells, which should have profound effects on the insulin autoimmune tolerance induction.

RESEARCH DESIGN AND METHODS

A repetitive DNA region designated as a variable number of tandem repeats (VNTR) is located upstream of the human insulin gene and correlates with the incidence of type 1 diabetes. We generated six class I and two class III VNTR constructs linked to the human insulin basal promoter or SV40 heterologous promoter/enhancer and demonstrated that AIRE protein modulates the insulin promoter activities differentially through binding to the VNTR region.

RESULTS

Here we show that in the presence of the autoimmune regulator (AIRE), the class III VNTR haplotype is responsible for an average of three-fold higher insulin expression than class I VNTR in thymic epithelial cells. In a protein-DNA pull-down experiment, AIRE protein is capable of binding to VNTR class I and III probes. Further, the transcriptional activation of the INS-VNTR by AIRE requires the insulin basal promoter. The VNTR sequence loses its activation activity when linked to a heterologous promoter and/or enhancer.

CONCLUSIONS

These findings demonstrate a type 1 diabetes predisposition encoded by the INS-VNTR locus and a critical function played by AIRE, which constitute a dual control mechanisms regulating quantitative expression of insulin in human thymic epithelial cells.

The insulin polymorphism -23Hph increases the risk for type 1 diabetes mellitus in the Romanian population

The insulin -23Hph and IGF2 Apa polymorphisms were genotyped in Romanian patients with T1DM (n = 204), T2DM (n = 215) or obesity (n = 200) and normoponderal healthy subjects (n = 750). The genotypes of both polymorphisms were distributed in concordance with Hardy-Weinberg equilibrium in all groups. The -23Hph AA genotype increased the risk for T1DM (OR: 3.22, 95%CI: 2.09-4.98, p < 0,0001), especially in patients without macroalbuminuria (OR: 4.32, 95%CI: 2.54-7.45, p < 0,0001). No other significant association between the alleles or genotypes of insulin -23Hph and IGF2 Apa and diabetes or obesity was identified.

Genome-wide data-mining of candidate human splice translational efficiency polymorphisms (STEPs) and an online database

BACKGROUND

Variation in pre-mRNA splicing is common and in some cases caused by genetic variants in intronic splicing motifs. Recent studies into the insulin gene (INS) discovered a polymorphism in a 5' non-coding intron that influences the likelihood of intron retention in the final mRNA, extending the 5' untranslated region and maintaining protein quality. Retention was also associated with increased insulin levels, suggesting that such variants--splice translational efficiency polymorphisms (STEPs)--may relate to disease phenotypes through differential protein expression. We set out to explore the prevalence of STEPs in the human genome and validate this new category of protein quantitative trait loci (pQTL) using publicly available data.

METHODOLOGY/PRINCIPAL FINDINGS

Gene transcript and variant data were collected and mined for candidate STEPs in motif regions. Sequences from transcripts containing potential STEPs were analysed for evidence of splice site recognition and an effect in expressed sequence tags (ESTs). 16 publicly released genome-wide association data sets of common diseases were searched for association to candidate polymorphisms with HapMap frequency data. Our study found 3324 candidate STEPs lying in motif sequences of 5' non-coding introns and further mining revealed 170 with transcript evidence of intron retention. 21 potential STEPs had EST evidence of intron retention or exon extension, as well as population frequency data for comparison.

CONCLUSIONS/SIGNIFICANCE

Results suggest that the insulin STEP was not a unique example and that many STEPs may occur genome-wide with potentially causal effects in complex disease. An online database of STEPs is freely accessible at http://dbstep.genes.org.uk/.

DBASS3 and DBASS5: databases of aberrant 3'- and 5'-splice sites

DBASS3 and DBASS5 provide comprehensive repositories of new exon boundaries that were induced by pathogenic mutations in human disease genes. Aberrant 5'- and 3'-splice sites were activated either by mutations in the consensus sequences of natural exon-intron junctions (cryptic sites) or elsewhere ('de novo' sites). DBASS3 and DBASS5 currently contain approximately 900 records of cryptic and de novo 3'- and 5'-splice sites that were produced by over a thousand different mutations in approximately 360 genes. DBASS3 and DBASS5 data can be searched by disease phenotype, gene, mutation, location of aberrant splice sites in introns and exons and their distance from authentic counterparts, by bibliographic references and by the splice-site strength estimated with several prediction algorithms. The user can also retrieve reference sequences of both aberrant and authentic splice sites with the underlying mutation. These data will facilitate identification of introns or exons frequently involved in aberrant splicing, mutation analysis of human disease genes and study of germline or somatic mutations that impair RNA processing. Finally, this resource will be useful for fine-tuning splice-site prediction algorithms, better definition of auxiliary splicing signals and design of new reporter assays. DBASS3 and DBASS5 are freely available at http://www.dbass.org.uk/.

Allele-specific recognition of the 3' splice site of INS intron 1

Genetic predisposition to type 1 diabetes (T1D) has been associated with a chromosome 11 locus centered on the proinsulin gene (INS) and with differential steady-state levels of INS RNA from T1D-predisposing and -protective haplotypes. Here, we show that the haplotype-specific expression is determined by INS variants that control the splicing efficiency of intron 1. The adenine allele at IVS1-6 (rs689), which rapidly expanded in modern humans, renders the 3' splice site of this intron more dependent on the auxiliary factor of U2 small nuclear ribonucleoprotein (U2AF). This interaction required both zinc fingers of the 35-kD U2AF subunit (U2AF35) and was associated with repression of a competing 3' splice site in INS exon 2. Systematic mutagenesis of reporter constructs showed that intron 1 removal was facilitated by conserved guanosine-rich enhancers and identified additional splicing regulatory motifs in exon 2. Sequencing of intron 1 in primates revealed that relaxation of its 3' splice site in Hominidae coevolved with the introduction of a short upstream open reading frame, providing a more efficient coupled splicing and translation control. Depletion of SR proteins 9G8 and transformer-2 by RNA interference was associated with exon 2 skipping whereas depletion of SRp20 with increased representation of transcripts containing a cryptic 3' splice site in the last exon. Together, these findings reveal critical interactions underlying the allele-dependent INS expression and INS-mediated risk of T1D and suggest that the increased requirement for U2AF35 in higher primates may hinder thymic presentation of autoantigens encoded by transcripts with weak 3' splice sites.

Insulin: a critical autoantigen and potential therapeutic agent in Type 1 diabetes

Insulin is a polypeptide hormone secreted by pancreatic beta-cells and is critical for glucose homeostasis. Abnormalities in insulin secretion result in various forms of diabetes. Type 1A diabetes is an autoimmune form in which insulin has been identified as a critical autoantigen. Recent studies have identified genetic determinants of insulin-specific autoimmune responses and insulin epitopes targeted by autoreactive T lymphocytes. The study of insulin as an autoantigen has also led to discoveries about basic mechanisms of immunological tolerance and autoimmunity. Experimental and clinical evidence suggests that insulin and insulin-derived peptides may delay and perhaps prevent the development of diabetes. Further clinical trials may identify effective treatment modalities for inhibiting diabetogenic autoimmunity and preventing disease development.

Dual effect of a single nucleotide polymorphism in the first intron of the porcine secreted phosphoprotein 1 gene: allele-specific binding of C/EBP beta and activation of aberrant splicing

Background

Secreted phosphoprotein 1 (SPP1 or Osteopontin, OPN) is a multifunctional matricellular glycoprotein involved in development and regeneration of skeletal muscle. Previously, we have demonstrated that porcine SPP1 shows breed-related differential mRNA expression during myogenesis. With the aim to identify putative contributing cis-regulatory DNA variation we resequenced the 5' upstream region of the gene in the respective breeds Pietrain and Duroc. We found two single nucleotide polymorphisms (SNP; [GenBank:M84121]: g.1804C>T and g.3836A>G). We focused our investigation on the SNP g.3836A>G, because in silico analysis and knowledge about the regulation of SPP1 suggested an effect of this SNP on a CCAAT/enhancer binding protein beta (C/EBPβ) responsive transcriptional enhancer.

Results

Using electrophoretic mobility shift assay we demonstrated that, similar to human SPP1, the 3' terminal end of the first intron of porcine SPP1 harbors a C/EBPβ binding site and showed that this binding site is negatively affected by the mutant G allele. Genotyping of 48 fetuses per breed revealed that the G allele segregated exclusively in Duroc fetuses with a frequency of 57 percent. Using real-time quantitative PCR we showed that, consistent with its negative effect on a transcriptional enhancer element, the G allele tends to decrease mRNA abundance of SPP1 in the fetal musculus longissimus dorsi (~1.3 fold; P ≥ 0.1).

Moreover, we showed that the SNP g.3836A>G leads to ubiquitous aberrant splicing of the first intron by generating a de novo and activating a cryptic splice acceptor site. Aberrantly spliced transcripts comprise about half of the SPP1 messages expressed by the G allele. Both aberrant splice variants differ from the native transcript by insertions in the leader sequences which do not change the reading frame of SPP1.

Conclusion

At the 3' terminal end of the first intron of the porcine SPP1 we identified a unique, dually functional SNP g.3836A>G. This SNP affects the function of the SPP1 gene at the DNA level by affecting a C/EBPβ binding site and at the RNA level by activating aberrant splicing of the first intron, and thus represents an interesting DNA-marker to study phenotypic effects of SPP1 DNA-variation.

Common polymorphic variation in the genetically diverse African insulin gene and its association with size at birth

The insulin variable number of tandem repeats (INS VNTR) has been variably associated with size at birth in non-African populations. Small size at birth is a major determinant of neonatal mortality, so the INS VNTR may influence survival. We tested the hypothesis, therefore, that genetic variation around the INS VNTR in a rural Gambian population, who experience seasonal variation in nutrition and subsequently birth weight, may be associated with foetal and early growth. Six polymorphisms flanking the INS VNTR were genotyped in over 2,500 people. Significant associations were detected between the maternally inherited SNP 27 (rs689) allele and birth length [effect size 17.5 (5.2-29.8) mm; P = 0.004; n = 361]. Significant associations were also found between the maternally inherited African-specific SNP 28 (rs5506) allele and post-natal weight gain [effect size 0.19 (0.05-0.32) z score points/year; P = 0.005; n = 728). These results suggest that in the Gambian population studied there are associations between polymorphic variation in the genetically diverse INS gene and foetal and early growth characteristics, which contribute to overall polygenic associations with these traits.

Genome-wide survey of allele-specific splicing in humans

Background

Accurate mRNA splicing depends on multiple regulatory signals encoded in the transcribed RNA sequence. Many examples of mutations within human splice regulatory regions that alter splicing qualitatively or quantitatively have been reported and allelic differences in mRNA splicing are likely to be a common and important source of phenotypic diversity at the molecular level, in addition to their contribution to genetic disease susceptibility. However, because the effect of a mutation on the efficiency of mRNA splicing is often difficult to predict, many mutations that cause disease through an effect on splicing are likely to remain undiscovered.

Results

We have combined a genome-wide scan for sequence polymorphisms likely to affect mRNA splicing with analysis of publicly available Expressed Sequence Tag (EST) and exon array data. The genome-wide scan uses published tools and identified 30,977 SNPs located within donor and acceptor splice sites, branch points and exonic splicing enhancer elements. For 1,185 candidate splicing polymorphisms the difference in splicing between alternative alleles was corroborated by publicly available exon array data from 166 lymphoblastoid cell lines. We developed a novel probabilistic method to infer allele-specific splicing from EST data. The method uses SNPs and alternative mRNA isoforms mapped to EST sequences and models both regulated alternative splicing as well as allele-specific splicing. We have also estimated heritability of splicing and report that a greater proportion of genes show evidence of splicing heritability than show heritability of overall gene expression level. Our results provide an extensive resource that can be used to assess the possible effect on splicing of human polymorphisms in putative splice-regulatory sites.

Conclusion

We report a set of genes showing evidence of allele-specific splicing from an integrated analysis of genomic polymorphisms, EST data and exon array data, including several examples for which there is experimental evidence of polymorphisms affecting splicing in the literature. We also present a set of novel allele-specific splicing candidates and discuss the strengths and weaknesses of alternative technologies for inferring the effect of sequence variants on mRNA splicing.

Association between paternally inherited haplotypes upstream of the insulin gene and umbilical cord IGF-II levels

The insulin (INS) and IGF 2 (IGF2) genes are in close proximity to each other and undergo maternal imprinting during fetal growth. We investigated the association between maternal and umbilical cord IGF 2 protein (IGF-II) levels and single nucleotide polymorphisms (SNPs) in the INS and IGF2 genes in 207 healthy African-American mother-newborn pairs. No association was found between maternal IGF-II levels and polymorphism in the INS-IGF2 locus. A significant association was found between newborn IGF-II levels and two SNPs (rs3842738 and rs689) at the 5' end of the INS-IGF2 locus. Analyses of haplotypes inferred from these two SNPs demonstrate a significant relationship between paternally transmitted haplotypes and newborn IGF-II levels, but no association with maternally transmitted haplotypes.

Global control of aberrant splice-site activation by auxiliary splicing sequences: evidence for a gradient in exon and intron definition

Auxiliary splicing signals play a major role in the regulation of constitutive and alternative pre-mRNA splicing, but their relative importance in selection of mutation-induced cryptic or de novo splice sites is poorly understood. Here, we show that exonic sequences between authentic and aberrant splice sites that were activated by splice-site mutations in human disease genes have lower frequencies of splicing enhancers and higher frequencies of splicing silencers than average exons. Conversely, sequences between authentic and intronic aberrant splice sites have more enhancers and less silencers than average introns. Exons that were skipped as a result of splice-site mutations were smaller, had lower SF2/ASF motif scores, a decreased availability of decoy splice sites and a higher density of silencers than exons in which splice-site mutation activated cryptic splice sites. These four variables were the strongest predictors of the two aberrant splicing events in a logistic regression model. Elimination or weakening of predicted silencers in two reporters consistently promoted use of intron-proximal splice sites if these elements were maintained at their original positions, with their modular combinations producing expected modification of splicing. Together, these results show the existence of a gradient in exon and intron definition at the level of pre-mRNA splicing and provide a basis for the development of computational tools that predict aberrant splicing outcomes.

Association between small for gestational age and paternally inherited 5' insulin haplotypes

OBJECTIVE

To test the association between small for gestational age and polymorphisms in the insulin gene in newborns and their mothers, as well as the effect of the parental transmission of haplotypes.

SUBJECTS

Pairs of healthy African-American full-term newborns (N=207) and mothers were recruited from Memphis TN and Jackson MS with birth weights ranging from 2210 to 4735 g.

METHODS

Six single nucleotide polymorphisms (SNPs) located in the insulin (INS) and insulin-like growth factor 2 (IGF2) genes were genotyped in mothers and newborns. Haplotypes composed of three SNPs in the 5' region of the INS-IGF2 locus were computationally inferred. Odds ratios for risk of small for gestational age (SGA) birth were calculated for individual SNPs and inferred haplotypes in the newborns and in the mothers using logistic regression. For 162 mother--newborn pairs the parental transmission of the haplotypes could be inferred, and the risks for SGA birth were calculated for the three common haplotypes in this sample.

RESULTS

Three INS SNPs exhibited significant association with risk for SGA birth. The SNP alleles associated with increased risk for SGA were opposite in the maternal and newborn genomes, implying opposing influences on the rate of fetal growth. Consistent with these results, haplotypes composed of complementary nucleotide sequences (CAC at rs3842738, rs689 and rs3842748, respectively, in the newborn versus GTG in the mother) were significantly associated with risk for SGA birth. In analyses of haplotypes according to parental transmission, the same trend in risk for SGA was observed for both maternally and paternally transmitted haplotypes, although a significant difference in risk was observed only for paternally transmitted haplotypes.

CONCLUSION

Polymorphisms near the 5' end of the INS-IGF2 locus are significantly associated with risk for SGA birth with a major effect due to the paternally transmitted haplotype, which is preferentially expressed due to imprinting.

Aberrant 5' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization

Despite a growing number of splicing mutations found in hereditary diseases, utilization of aberrant splice sites and their effects on gene expression remain challenging to predict. We compiled sequences of 346 aberrant 5′splice sites (5′ss) that were activated by mutations in 166 human disease genes. Mutations within the 5′ss consensus accounted for 254 cryptic 5′ss and mutations elsewhere activated 92 de novo 5′ss. Point mutations leading to cryptic 5′ss activation were most common in the first intron nucleotide, followed by the fifth nucleotide. Substitutions at position +5 were exclusively G>A transitions, which was largely attributable to high mutability rates of C/G>T/A. However, the frequency of point mutations at position +5 was significantly higher than that observed in the Human Gene Mutation Database, suggesting that alterations of this position are particularly prone to aberrant splicing, possibly due to a requirement for sequential interactions with U1 and U6 snRNAs. Cryptic 5′ss were best predicted by computational algorithms that accommodate nucleotide dependencies and not by weight-matrix models. Discrimination of intronic 5′ss from their authentic counterparts was less effective than for exonic sites, as the former were intrinsically stronger than the latter. Computational prediction of exonic de novo 5′ss was poor, suggesting that their activation critically depends on exonic splicing enhancers or silencers. The authentic counterparts of aberrant 5′ss were significantly weaker than the average human 5′ss. The development of an online database of aberrant 5′ss will be useful for studying basic mechanisms of splice-site selection, identifying splicing mutations and optimizing splice-site prediction algorithms.

Quantitated transcript haplotypes (QTH) of AGTR1, reduced abundance of mRNA haplotypes containing 1166C (rs5186:A>C), and relevance to metabolic syndrome traits

The angiotensin II type 1 receptor (AGTR1) is the main target through which angiotensin II influences cardiovascular tone, cell growth, and fluid and electrolyte balance. AGTR1 polymorphism has been reported to associate with hypertension, myocardial infarction (MI), and metabolic traits. Here we describe a novel approach to quantitation of transcript haplotypes (QTH) of AGTR1. To determine relative allelic expression from haplotypes, within-individual-between-allele ratiometric analyses in placental cDNA were developed for the transcribed SNPs rs5182:C>T (encoding p.L191) and rs5186:A>C (3'-noncoding "A1166C"). Additionally, between-individual comparisons were made using TaqMan assays applied to both homozygous and heterozygous genotypes and haplotypes. In conjunction, linkage disequilibrium (LD) and genomic haplotype associations with metabolic syndrome were examined. There was no significant difference of mRNA level for alleles of rs5182:C>T, but allele and mRNA haplotypes carrying 1166C exhibited reduced abundance. The effect was much greater in CC homozygotes than in heterozygotes. The promoter region was confirmed to be in a separate haplotype block from the AGTR1 3' region containing rs5182:C>T and rs5186:A>C. Metabolic syndrome trait associations were strongest for the 3' block generally and for the C allele of rs5186:A>C specifically. All effects were much more prominent in homozygotes, possibly reflecting interallelic interaction through feedback loops of mRNA regulation. Differential abundance of AGTR1 mRNA haplotypes may mediate clinical phenotypic observations of the AGTR1 genotype.

Evaluation of polymorphic splicing in the mechanism of the association of the insulin gene with diabetes

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 (r(2) = 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 ( approximately 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 study of association between common variation in the growth hormone-chorionic somatomammotropin hormone gene cluster and adult fasting insulin in a UK Caucasian population

Background

Reduced growth during infancy is associated with adult insulin resistance. In a UK Caucasian cohort, the CSH1.01 microsatellite polymorphism in the growth hormone-chorionic somatomammotropin hormone gene cluster was recently associated with increases in adult fasting insulin of approximately 23 pmol/l for TT homozygote males compared to D1D1 or D2D2 homozygotes (P = 0.001 and 0.009; n = 206 and 92, respectively), but not for females. TT males additionally had a 547-g lower weight at 1 year (n = 270; P = 0.008) than D2D2 males. We sought to replicate these data in healthy UK Caucasian subjects. We genotyped 1396 subjects (fathers, mothers and children) from a consecutive birth study for the CSH1.01 marker and analysed genotypes for association with 1-year weight in boys and fasting insulin in fathers.

Results

We found no evidence for association of CSH1.01 genotype with adult male fasting insulin concentrations (TT/D1D1 P = 0.38; TT/D2D2 P = 0.18) or weight at 1 year in boys (TT/D1D1 P = 0.76; TT/D2D2 P = 0.85). For fasting insulin, our data can exclude the previously observed effect sizes as the 95 % confidence intervals for the differences observed in our study exclude increases in fasting insulin of 9.0 and 12.6 pmol/l for TT relative to D1D1 and D2D2 homozygotes, respectively. Whilst we have fewer data on boys' 1-year weight than the original study, our data can exclude a reduction in 1-year weight greater than 557 g for TT relative to D2D2 homozygotes.

Conclusion

We have not found association of the CSH1.01 genotype with fasting insulin or weight at 1 year. We conclude that the original study is likely to have over-estimated the effect size for fasting insulin, or that the difference in results reflects the younger age of subjects in this study relative to those in the previous study.

Refined Association Mapping for a Quantitative Trait: Weight in the H19-IGF2-INS-TH Region

Previous analyses have provided evidence for one or more loci affecting body weight in the H19-IGF2-INS-TH region on chromosome 11p15. To identify the location of a possible causal locus or loci we applied association analysis by composite likelihood to a large cohort under the Malecot model for body weight. A random sample of 2731 men in the UK were typed for eleven single nucleotide polymorphisms (SNPs) in IGF2, two SNPs in H19, one SNP in INS and one microsatellite marker in the TH genes. Using F tests appropriate to small marker sets, the superiority of regression over correlation was confirmed. All the evidence for association came from IGF2, with P= 0.007 for height-adjusted weight and P= 0.019 for weight additionally adjusted for smoking and alcohol drinking. Although the estimated point location for the suspected causal variant was close to IGF2 ApaI, the 95% confidence and support intervals covered most of IGF2 but none of the other loci. Identification of the causal SNP or SNPs within IGF2 will require typing of more variants in this region.

Aberrant 3' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization

The frequency distribution of mutation-induced aberrant 3' splice sites (3'ss) in exons and introns is more complex than for 5' splice sites, largely owing to sequence constraints upstream of intron/exon boundaries. As a result, prediction of their localization remains a challenging task. Here, nucleotide sequences of previously reported 218 aberrant 3'ss activated by disease-causing mutations in 131 human genes were compared with their authentic counterparts using currently available splice site prediction tools. Each tested algorithm distinguished authentic 3'ss from cryptic sites more effectively than from de novo sites. The best discrimination between aberrant and authentic 3'ss was achieved by the maximum entropy model. Almost one half of aberrant 3'ss was activated by AG-creating mutations and approximately 95% of the newly created AGs were selected in vivo. The overall nucleotide structure upstream of aberrant 3'ss was characterized by higher purine content than for authentic sites, particularly in position -3, that may be compensated by more stringent requirements for positive and negative nucleotide signatures centred around position -11. A newly developed online database of aberrant 3'ss will facilitate identification of splicing mutations in a gene or phenotype of interest and future optimization of splice site prediction tools.

Model of stimulation-responsive splicing and strategies in identification of immunogenic isoforms of tumor antigens and autoantigens

We recently proposed a novel model of stimulation-responsive splicing for the selection of autoantigens and self-tumor antigens. Our model theorizes that the significantly higher rates of alternative splicing of autoantigen and self-tumor antigen transcripts that occur in response to stimuli could induce extra-thymic expression of untolerized antigen epitopes for elicitation of autoimmune and anti-tumor responses. To facilitate the identification of immunogenic isoforms of antigens, we have developed strategies using improved SEREX in conjunction with database-mining and immunogenic isoform mapping. Identification of immunogenic isoforms of autoantigens and self-tumor antigens is very important for the development of novel therapeutics and diagnostic tools for autoimmune diseases and tumors, such as: (1) autoantigen isoform microarrays for disease diagnosis and prognosis; (2) autoantigen isoform-specific tolerizing vaccines and splicing-redirection therapies, as well as (3) immunogenic antigen isoform-specific immunotherapy for tumors.

The variable number of tandem repeats upstream of the insulin gene is a susceptibility locus for latent autoimmune diabetes in adults

The etiopathological relationship between latent autoimmune diabetes in adults (LADA) and classical type 1 (insulin dependent) diabetes remains unclear. Variation at the insulin gene variable number of tandem repeats (VNTR) minisatellite influences susceptibility to type 1 diabetes, but studies in LADA have been small and inconsistent. We examined the role of insulin gene variation (using flanking variants as surrogates for VNTR subtypes) in the largest case-control study of LADA to date (400 case and 332 control subjects). Highly significant associations were identified with disease, with dominant protective effects of the T allele at -23HphI (odds ratio [OR] 0.42 [95% CI 0.31-0.58], P=2.4 x 10(-8)), A allele at +1,404Fnu4HI (0.50 [0.36-0.70], P=3.2 x 10(-5)), and C allele at +3,580MspI (0.55 [0.35-0.85], P=0.0046). As with type 1 diabetes, the -23HphI variant (a surrogate for the subdivision of VNTR into class I and III alleles) most clearly defined susceptibility in LADA. However, there was no association with age at diagnosis or requirement for insulin therapy 6 years postdiagnosis. This study establishes that variation within the insulin gene region does influence susceptibility to LADA, with the direction and magnitude of effect indistinguishable from that previously reported for type 1 diabetes. In conclusion, differences in VNTR-encoded susceptibility do not explain the differences in clinical presentation that distinguish classical type 1 diabetes and LADA.

Phenotypic consequences of branch point substitutions

The branch point sequence (BPS) is a conserved splicing signal important for spliceosome assembly and lariat intron formation. BPS mutations may result in aberrant pre-mRNA splicing and genetic disorders, but their phenotypic consequences have been difficult to predict, largely due to a highly degenerate nature of the BPS consensus. Here, we have examined the splicing pattern of nine reporter pre-mRNAs that have previously been shown to give rise to human hereditary diseases as a result of single-nucleotide substitutions in the predicted BPS. Increased exon skipping and intron retention observed in vivo were recapitulated for each mutated pre-mRNA, but the reproducibility of cryptic splice site activation was lower. BP mutations in reporter pre-mRNAs frequently induced aberrant 3' splice sites and also activated a cryptic 5' splice site. Systematic mutagenesis of BP adenosines showed that in most pre-mRNAs, the expression of canonical transcripts was lower for BP transitions than BP transversions. Differential splicing outcome for transitions vs. transversions was abrogated or reduced if introns were truncated to 200 nt or less, suggesting that the nature of the BP residue is less critical for interactions across very short introns. Together, these results improve prediction of phenotypic consequences of point mutations upstream of splice acceptor sites and suggest that the overrepresentation of disease-causing adenosine-to-guanosine BP substitutions observed in Mendelian disorders is due to more profound defects of gene expression at the level of pre-mRNA splicing.