ESEfinder: A web resource to identify exonic splicing enhancers

Mutation spectrum leading to an attenuated phenotype in dystrophinopathies

Although Becker muscular dystrophy (BMD; MIM 300376) is mainly caused by gross deletions of the dystrophin gene, the nature of the mutations involved in the remaining cases is of importance because of the milder clinical course of Becker. We have extensively characterized the mRNA changes associated with five novel point mutations giving rise to a Becker phenotype, which confirm that Becker arises largely due to alterations in splicing. In two cases the milder phenotype arises because of exon skipping, leading to an in-frame deletion (c.1603-2A>C and c.4250T>A). In further two cases intronic mutations (c.4519-5C>G and c.961-5925A>C) result in complex splicing changes, but with some residual normal transcripts. The last case, c.10412T>A (p.Leu3471X), results in a truncated transcript missing only part of the COOH terminal of the protein, suggesting that this region is not crucial for dystrophin function. The detection of a low amount of dystrophin in this patient could be attributable to a reduced efficiency of nonsense-mediated decay. The results emphasize that mRNA analysis is important in defining Becker mutations and will be of value in assessing various gene therapy strategies.

Hereditary angioedema: The mutation spectrum ofSERPING1/C1NHin a large Spanish cohort

Hereditary angioedema (HAE) is a disease caused by defects in the C1 inhibitor gene (SERPING1/C1NH). We screened the entire C1NH gene for mutations in a large series of 87 Spanish families (77 with type I, and 10 with type II HAE) by SSCP, sequencing, Southern blotting, and quantitative multiplex PCR of short fluorescent fragments (QMPSF), and we characterized several defects at the mRNA level. We found large rearrangements in 13 families, and point mutations or microdeletions/insertions in 74 families. The 13 large rearrangements included nine exon deletions, of which at least eight were distinct, two were distinct exon duplications, and two were rearrangements whose precise nature could not be determined. We confirmed that exon 4 is particularly prone to rearrangements. Thirty-six mutations were unreported, and included 10 microdeletions/insertions, 10 missense, five nonsense, eight splicing, and three splicing or missense mutations. Moreover, we detected six novel uncharacterized sequence variants (USV). RT-PCR studies showed that in addition to several intronic splice site mutations tested, the exonic mutations c.882C>G and c.884T>G, located near the 3' end of exon 5, also produced exon skipping. This is the first evidence of SERPING1/C1NH mutations in coding regions that differ from the canonical splice sites that affect splicing, which suggests the presence of an exonic splicing enhancer (ESE) in exon 5.

Familial adenomatous polyposis: aberrant splicing due to missense or silent mutations in the APC gene

Familial adenomatous polyposis (FAP) is caused by germline mutations in the tumor suppressor gene APC. To date, the relevance of rare exonic single-base substitutions at nucleotide positions close to splice sites that are predicted to result in missense or silent (SNP) variants or substitutions in introns at splice-site positions that are not highly conserved has not been systematically examined in FAP patients. In 34 index patients, we identified 26 different heterozygous single-base substitutions at or close to the splice sites. We characterized five exonic mutations in exon 4 (c.423G>T), exon 14 (c.1956C>T, c.1957A>G, and c.1957A>C), and exon 15 (c.1959G>A) by transcript analysis and by splice-prediction programs (BDGP, SpliceSiteFinder, and ESEfinder). The splicing patterns of these variants were compared to those of 16 different substitutions at highly or less-conserved intronic splice-site positions, and to normal controls. In addition, we analyzed cosegregation of the variants with affected family members and examined the genotype-phenotype correlation. We could demonstrate that the four unclear variants in exon 4 and 14 that are predicted to result in missense or silent mutations in fact lead to complete exon skipping due to aberrant splicing; one possible explanation for this observed effect might be the disruption of exonic splicing enhancer (ESE) motifs. In contrast, the substitution at the first position of exon 15 seems to actually be a silent variant. We present the first systematic evaluation of different single-base substitutions in APC at or close to splice sites at transcript level. We show that the consequence of exonic mutations cannot be evaluated only by the predicted change in amino acid sequence but rather by the change at DNA level. The functional analysis of variants with unknown pathogenic effect plays an important role in increasing the mutation detection rate and achieving validation of predictive testing.

FAST DB: a website resource for the study of the expression regulation of human gene products

Human genes use various mechanisms to generate different transcripts having different exon content, which in turn generate multiple protein isoforms having differential and even opposite biological activities. To understand the biological consequences of gene transcriptional activity modulation, it is necessary to integrate the capability of genes to generate distinct functional products, particularly because transcriptional stimuli also affect the exon content of their target gene products. For this purpose, we have developed a bioinformatics suite, FAST DB, which defines easily and accurately the exon content of all known transcripts produced by human genes. In addition, several tools have been developed, including a graphical presentation of all gene products, a sequence multi-alignment of all gene transcripts and an in silico PCR computer program. The FAST DB interface also offers extensive links to website resources for promoter analysis and transcription factor binding site prediction, splicing regulatory sequence prediction, as well as 5′- and 3′-untranslated region analysis. FAST DB has been designed to facilitate studies that integrate transcriptional and post-transcriptional events to investigate the expression regulation of human gene products.

The Alu insertion in the CLCN5 gene of a patient with Dent's disease leads to exon 11 skipping

Alu sequences are short, interspersed elements that have generated more than one million copies in the human genome. They propagate by transcription followed by reverse transcription and integration, causing mutations, recombination, and changes in pre-mRNA splicing. We have recently identified a 345-bp long Alu Ya5 element inserted in codon 650 within exon 11 of the chloride channel ClC-5 gene (CLCN5) of a patient with Dent's disease. A microsatellite pedigree analysis indicated that the insertion occurred in the germline of the maternal grandfather. Dent's disease is an X-linked renal tubular disorder characterized by low-molecular-weight proteinuria, hypercalciuria, nephrolithiasis, and nephrocalcinosis. Here, we found, by RT-PCR amplification of RNA extracted from the patient's blood and subsequent DNA sequencing, that the Alu insertion led to an aberrant splicing of the CLCN5 pre-mRNA that skipped exon 11. Using the ESE finder and RESCUE-ESE Web interfaces, we identified two high-score exonic splicing enhancer (ESE) sequences in the site of insertion. The functional significance of these ESE motifs is suggested by our observation that these sequences are highly conserved among mammal CLCN5 genes. Therefore, we suggest that the Alu insertion causes exon skipping by interfering with splicing regulatory elements. The altered splicing would predict a truncated ClC-5 protein that lacks critical domains for sorting and chloride channel function.

Exonization of AluYa5 in the human ACE gene requires mutations in both 3' and 5' splice sites and is facilitated by a conserved splicing enhancer

Ancient Alu elements have been shown to be included in mature transcripts by point mutations that improve their 5′ or 3′ splice sites. We have examined requirements for exonization of a younger, disease-associated AluYa5 in intron 16 of the human ACE gene. A single G>C transversion in position −3 of the new Alu exon was insufficient for Alu exonization and a significant inclusion in mRNA was only observed when improving several potential splice donor sites in the presence of 3′ CAG. Since complete Alu exonization was not achieved by optimizing traditional splicing signals, including the branch site, we tested whether auxiliary elements in AluYa5 were required for constitutive inclusion. Exonization was promoted by a SELEX-predicted heptamer in Alu consensus sequence 222–228 and point mutations in highly conserved nucleotides of this heptamer decreased Alu inclusion. In addition, we show that Alu exonization was facilitated by a subset of serine/arginine-rich (SR) proteins through activation of the optimized 3′ splice site. Finally, the haplotype- and allele-specific ACE minigenes generated similar splicing patterns in both ACE-expressing and non-expressing cells, suggesting that previously reported allelic association with plasma ACE activity and cardiovascular disease is not attributable to differential splicing of introns 16 and 17.

Understanding alternative splicing: towards a cellular code

In violation of the 'one gene, one polypeptide' rule, alternative splicing allows individual genes to produce multiple protein isoforms - thereby playing a central part in generating complex proteomes. Alternative splicing also has a largely hidden function in quantitative gene control, by targeting RNAs for nonsense-mediated decay. Traditional gene-by-gene investigations of alternative splicing mechanisms are now being complemented by global approaches. These promise to reveal details of the nature and operation of cellular codes that are constituted by combinations of regulatory elements in pre-mRNA substrates and by cellular complements of splicing regulators, which together determine regulated splicing pathways.

A mutation-created novel intra-exonic pre-mRNA splice site causes constitutive activation of KIT in human gastrointestinal stromal tumors

We report a new mechanism of aberrant pre-mRNA splicing resulting in constitutive activation of a mis-spliced oncoprotein (KIT) leading to malignancy (gastrointestinal stromal tumor) in contrast to loss of function of mis-spliced proteins resulting in diverse human diseases in the literature. The mechanisms of three consecutive molecular events, deletion of noncoding and coding regions encompassing the 3' authentic splice site, creation of a novel intra-exonic pre-mRNA 3' splice acceptor site leading to in-frame loss of 27 nucleotides (nine amino acids; Lys550-Lys558), and the mechanism of constitutive activation of the mis-spliced KIT are elucidated. Loss of a peptide in a critical location unleashed the protein from autoinhibition (as evidenced by three-dimensional structural analysis), causing KIT to become constitutively activated and resulting in the GIST phenotype. We also demonstrated that only one of the following two exonic splicing enhancers is sufficient for inclusion of the KIT exon 11 in vivo: AACCCATGT (nucleotides 2-10 from the 5' end, which are recognized by SC35, SRp55, and SF2/ASF) or GGTTGTTGAGG (nucleotides 27-37 from the 5' end, which are recognized by SC35 and SRp55), suggestive of exonic enhancer redundancy.

Potential impact of a single nucleotide polymorphism in the hyaluronan synthase 1 gene in Waldenstrom's macroglobulinemia

The hyaluronan synthase 1 (HAS1) gene encodes a plasma membrane protein that synthesizes hyaluronan, an extracellular matrix molecule. Previously, in patients with Waldenstrom's macroglobulinemia (WM), we detected upregulation of HAS1 transcripts and identified aberrant splice variants of this gene. Aberrant splicing of HAS1 results from activation of cryptic splice sites. In turn, activation of cryptic donor and acceptor splice sites can be promoted by mutations occurring upstream of these sites and/or at the branch point of slicing. We measured the frequency of the HAS1 833A/G polymorphism (ie, single-nucleotide polymorphism; SNP) in patients with WM and healthy donors. Additionally, HAS1 gene expression was evaluated in the same group of patients. Our observations so far suggest that HAS1 833A/G SNPs contribute to aberrant splicing of this gene; this idea is supported by the fact that 833A/G SNP is located on an exonic splicing enhancer motif. Based on the results obtained thus far, we speculate that individuals with HAS1 833G/G genotype are predisposed toward aberrant HAS1 splicing and expression of HAS1 variants, resulting in an enhanced risk of developing WM. Study of a larger group of patients and healthy donors is needed to confirm these speculations and to evaluate the prognostic significance of these findings.

A duplication in the canine beta-galactosidase gene GLB1 causes exon skipping and GM1-gangliosidosis in Alaskan huskies

GM(1)-gangliosidosis is a lysosomal storage disease that is inherited as an autosomal recessive disorder, predominantly caused by structural defects in the beta-galactosidase gene (GLB1). The molecular cause of GM(1)-gangliosidosis in Alaskan huskies was investigated and a novel 19-bp duplication in exon 15 of the GLB1 gene was identified. The duplication comprised positions +1688-+1706 of the GLB1 cDNA. It partially disrupted a potential exon splicing enhancer (ESE), leading to exon skipping in a fraction of the transcripts. Thus, the mutation caused the expression of two different mRNAs from the mutant allele. One transcript contained the complete exon 15 with the 19-bp duplication, while the other transcript lacked exon 15. In the transcript containing exon 15 with the 19-bp duplication a premature termination codon (PTC) appeared, but due to its localization in the last exon of canine GLB1, nonsense-mediated RNA decay (NMD) did not occur. As a consequence of these molecular events two different truncated GLB1 proteins are predicted to be expressed from the mutant GLB1 allele. In heterozygous carrier animals the wild-type allele produces sufficient amounts of the active enzyme to prevent clinical signs of disease. In affected homozygous dogs no functional GLB1 is synthesized and G(M1)-gangliosidosis occurs.

Inter-individual variation of several cytochrome P450 2D6 splice variants in human liver

To examine the possibility that inter-individual differences in splicing partially explain the observed differences in CYP2D6 activity, we amplified its full-length cDNA in 96 human liver RNA samples and discovered five splice variants: intron 5 retention, intron 6 retention, intron 5 and intron 6 double retention, exon 3 skipping, and partial intron 1 retention. All of the CYP2D6 splice variants we identified are probably nonfunctional transcripts. Substantial inter-individual variation in the proportions of the CYP2D6 transcript represented by splice variants, measured by real-time PCR, suggests that the presence of these splice variants contributes to the population variation in CYP2D6 activity. Relatively high levels of intron 6 retention were not correlated with the newly discovered single nucleotide polymorphism 2988G > A in intron 6 (CYP2D6*41) but did correlate with the more common CYP2D6*34 allele. Our study prompts further investigations to explore the effect of these splice variants on drug metabolism.

Selective modification of alternative splicing by indole derivatives that target serine-arginine-rich protein splicing factors

The prevalence of alternative splicing as a target for alterations leading to human genetic disorders makes it highly relevant for therapy. Here we have used in vitro splicing reactions with different splicing reporter constructs to screen 4,000 chemical compounds for their ability to selectively inhibit spliceosome assembly and splicing. We discovered indole derivatives as potent inhibitors of the splicing reaction. Importantly, compounds of this family specifically inhibit exonic splicing enhancer (ESE)-dependent splicing, because they interact directly and selectively with members of the serine-arginine-rich protein family. Treatment of cells expressing reporter constructs with ESE sequences demonstrated that selected indole derivatives mediate inhibition of ESE usage in vivo and prevent early splicing events required for HIV replication. This discovery opens the exciting possibility of a causal pharmacological treatment of aberrant splicing in human genetic disorders and development of new antiviral therapeutic approaches.

Mutations in theNDUFS4gene of mitochondrial complex I alter stability of the splice variants

The effect on the stability of alternative transcripts of different mutations of the NDUFS4 gene in patients with Leigh syndrome with complex I deficiency is presented. Normally, two NDUFS4 splice variants are degraded by nonsense mediated mRNA decay (NMD) while a third form does not trigger NMD degradation. In a patient with a premature termination codon in exon 1, all the three splice variants are up-regulated. The present is the first case of a nonsense mutation leading to the abrogation of NMD, which can represent an additional event to be considered in the evaluation of clinically relevant mutations.

Candidate-gene screening and association analysis at the autism-susceptibility locus on chromosome 16p: evidence of association at GRIN2A and ABAT

Autism is a highly heritable neurodevelopmental disorder whose underlying genetic causes have yet to be identified. To date, there have been eight genome screens for autism, two of which identified a putative susceptibility locus on chromosome 16p. In the present study, 10 positional candidate genes that map to 16p11-13 were examined for coding variants: A2BP1, ABAT, BFAR, CREBBP, EMP2, GRIN2A, MRTF-B, SSTR5, TBX6, and UBN1. Screening of all coding and regulatory regions by denaturing high-performance liquid chromatography identified seven nonsynonymous changes. Five of these mutations were found to cosegregate with autism, but the mutations are not predicted to have deleterious effects on protein structure and are unlikely to represent significant etiological variants. Selected variants from candidate genes were genotyped in the entire International Molecular Genetics Study of Autism Consortium collection of 239 multiplex families and were tested for association with autism by use of the pedigree disequilibrium test. Additionally, genotype frequencies were compared between 239 unrelated affected individuals and 192 controls. Patterns of linkage disequilibrium were investigated, and the transmission of haplotypes across candidate genes was tested for association. Evidence of single-marker association was found for variants in ABAT, CREBBP, and GRIN2A. Within these genes, 12 single-nucleotide polymorphisms (SNPs) were subsequently genotyped in 91 autism trios (one affected individual and two unaffected parents), and the association was replicated within GRIN2A (Fisher's exact test, P<.0001). Logistic regression analysis of SNP data across GRIN2A and ABAT showed a trend toward haplotypic differences between cases and controls.

Dichotomous splicing signals in exon flanks

Intronic elements flanking the splice-site consensus sequences are thought to play a role in pre-mRNA splicing. However, the generality of this role, the catalog of effective sequences, and the mechanisms involved are still lacking. Using molecular genetic tests, we first showed that the approximately 50-nt intronic flanking sequences of exons beyond the splice-site consensus are generally important for splicing. We then went on to characterize exon flank sequences on a genomic scale. The G+C content of flanks displayed a bimodal distribution reflecting an exaggeration of this base composition in flanks relative to the gene as a whole. We divided all exons into two classes according to their flank G+C content and used computational and statistical methods to define pentamers of high relative abundance and phylogenetic conservation in exon flanks. Upstream pentamers were often common to the two classes, whereas downstream pentamers were totally different. Upstream and downstream pentamers were often identical around low G+C exons, and in contrast, were often complementary around high G+C exons. In agreement with this complementarity, predicted base pairing was more frequent between the flanks of high G+C exons. Pseudo exons did not exhibit this behavior, but rather tended to form base pairs between flanks and exon bodies. We conclude that most exons require signals in their immediate flanks for efficient splicing. G+C content is a sequence feature correlated with many genetic and genomic attributes. We speculate that there may be different mechanisms for splice site recognition depending on G+C content.

Disruption of an exon splicing enhancer in exon 3 of MLH1 is the cause of HNPCC in a Quebec family

BACKGROUND

A 3 bp deletion located at the 5' end of exon 3 of MLH1, resulting in deletion of exon 3 from RNA, was recently identified.

HYPOTHESIS

That this mutation disrupts an exon splicing enhancer (ESE) because it occurs in a purine-rich sequence previously identified as an ESE in other genes, and ESEs are often found in exons with splice signals that deviate from the consensus signals, as does the 3' splice signal in exon 3 of MLH1.

DESIGN

The 3 bp deletion and several other mutations were created by polymerase chain reaction mutagenesis and tested using an in vitro splicing assay. Both mutant and wild type exon 3 sequences were cloned into an exon trapping vector and transiently expressed in Cos-1 cells.

RESULTS

Analysis of the RNA indicates that the 3 bp deletion c.213_215delAGA (gi:28559089, NM_000249.2), a silent mutation c.216T-->C, a missense mutation c.214G-->C, and a nonsense mutation c.214G-->T all cause varying degrees of exon skipping, suggesting the presence of an ESE at the 5' end of exon 3. These mutations are situated in a GAAGAT sequence 3 bp downstream from the start of exon 3.

CONCLUSIONS

The results of the splicing assay suggest that inclusion of exon 3 in the mRNA is ESE dependent. The exon 3 ESE is not recognised by all available motif scoring matrices, highlighting the importance of RNA analysis in the detection of ESE disrupting mutations.

The SR protein SC35 is responsible for aberrant splicing of the E1alpha pyruvate dehydrogenase mRNA in a case of mental retardation with lactic acidosis

Pyruvate dehydrogenase (PDH) complex deficiency is a major cause of lactic acidosis and Leigh's encephalomyelopathies in infancy and childhood, resulting in early death in the majority of patients. Most of the molecular defects have been localized in the coding regions of the E1alpha PDH gene. Recently, we identified a novel mutation of the E1alpha PDH gene in a patient with an encephalopathy and lactic acidosis. This mutation, located downstream of exon 7, activates a cryptic splice donor and leads to the retention of intronic sequences. Here, we demonstrate that the mutation results in an increased binding of the SR protein SC35. Consistently, ectopic overexpression of this splicing factor enhanced the use of the cryptic splice site, whereas small interfering RNA-mediated reduction of the SC35 protein levels in primary fibroblasts from the patient resulted in the almost complete disappearance of the aberrantly spliced E1alpha PDH mRNA. Our findings open the exciting prospect for a novel therapy of an inherited disease by altering the level of a specific splicing factor.

The homozygous state for Hb Crete [β129 (H7) Ala→Pro] is associated with a complex phenotype including erythrocytosis and functional anemia

Hb Crete, an electrophoretically neutral, unstable, high oxygen affinity variant, was characterized by protein and DNA analyses in the homozygous state in a 32-year-old woman from Crete, with erythrocytosis and microcytosis. The proband and members of her family over 3 generations, including 5 carriers of Hb Crete, were subject to clinical, hematological and biochemical investigations, and DNA, RNA and protein studies were carried out. The proband demonstrated features associated with disturbed hemoglobin (Hb) structure and function, including erythrocytosis and additionally a state of functional anemia, the latter reflected by increased erythropoetin levels and cardiac output. In addition, all the carriers surprisingly had hematological and biosynthetic findings more usually associated with thalassemia trait. The structural change in Hb Crete only partly explains all the pathological manifestations of this variant, and other mechanisms are discussed.

A combinatorial code for splicing silencing: UAGG and GGGG motifs

Alternative pre-mRNA splicing is widely used to regulate gene expression by tuning the levels of tissue-specific mRNA isoforms. Few regulatory mechanisms are understood at the level of combinatorial control despite numerous sequences, distinct from splice sites, that have been shown to play roles in splicing enhancement or silencing. Here we use molecular approaches to identify a ternary combination of exonic UAGG and 5'-splice-site-proximal GGGG motifs that functions cooperatively to silence the brain-region-specific CI cassette exon (exon 19) of the glutamate NMDA R1 receptor (GRIN1) transcript. Disruption of three components of the motif pattern converted the CI cassette into a constitutive exon, while predominant skipping was conferred when the same components were introduced, de novo, into a heterologous constitutive exon. Predominant exon silencing was directed by the motif pattern in the presence of six competing exonic splicing enhancers, and this effect was retained after systematically repositioning the two exonic UAGGs within the CI cassette. In this system, hnRNP A1 was shown to mediate silencing while hnRNP H antagonized silencing. Genome-wide computational analysis combined with RT-PCR testing showed that a class of skipped human and mouse exons can be identified by searches that preserve the sequence and spatial configuration of the UAGG and GGGG motifs. This analysis suggests that the multi-component silencing code may play an important role in the tissue-specific regulation of the CI cassette exon, and that it may serve more generally as a molecular language to allow for intricate adjustments and the coordination of splicing patterns from different genes.

Rat nicastrin gene: cDNA isolation, mRNA variants and expression pattern analysis

Nicastrin is a type 1 transmembrane glycoprotein that interacts with presenilin, Aph-1, and Pen-2 proteins to form a high molecular complex with gamma secretase activity. Then, nicastrin has a central role in presenilin-mediated processing of beta-amyloid precursor protein and in some aspects of Notch/glp-1 signaling in vivo. Here, we isolated a rat nicastrin cDNA and investigated gene expression in embryonic and adult rat tissues. The predicted amino acid sequence is comprised of 708 residues and showed a high degree of identity with other vertebrate orthologs. Besides full-length nicastrin mRNA, we identified an alternative spliced variant lacking the whole exon 3 and predicted to encode a 62-residue-long truncated protein. Full-length nicastrin mRNA was observed to be ubiquitously expressed, while the spliced variant was preferentially transcribed in the nervous system, whether in embryonic or adult neural tissues. Studies performed on primary cell cultures demonstrated that the short isoform was expressed in neurons, but not in astrocyte and microglial cells. Further experiments performed to verify the presence of the variant in neuroblastoma culture failed to show any truncated protein. Treatments by cyclohexamide showed the involvement of a quality control-based surveillance mechanism, which selectively degrades the exon 3-skipped isoform. In summary, this is the first report describing a novel skipped isoform of nicastrin which may suggest a new possible control mechanism based on the alternative splicing and nonsense-mediated mRNA decay to regulate brain protein expression and provide newer insights into potential implication in Alzheimer's disease.

A combinatorial code for splicing silencing: UAGG and GGGG motifs

Alternative pre-mRNA splicing is widely used to regulate gene expression by tuning the levels of tissue-specific mRNA isoforms. Few regulatory mechanisms are understood at the level of combinatorial control despite numerous sequences, distinct from splice sites, that have been shown to play roles in splicing enhancement or silencing. Here we use molecular approaches to identify a ternary combination of exonic UAGG and 5′-splice-site-proximal GGGG motifs that functions cooperatively to silence the brain-region-specific CI cassette exon (exon 19) of the glutamate NMDA R1 receptor (GRIN1) transcript. Disruption of three components of the motif pattern converted the CI cassette into a constitutive exon, while predominant skipping was conferred when the same components were introduced, de novo, into a heterologous constitutive exon. Predominant exon silencing was directed by the motif pattern in the presence of six competing exonic splicing enhancers, and this effect was retained after systematically repositioning the two exonic UAGGs within the CI cassette. In this system, hnRNP A1 was shown to mediate silencing while hnRNP H antagonized silencing. Genome-wide computational analysis combined with RT-PCR testing showed that a class of skipped human and mouse exons can be identified by searches that preserve the sequence and spatial configuration of the UAGG and GGGG motifs. This analysis suggests that the multi-component silencing code may play an important role in the tissue-specific regulation of the CI cassette exon, and that it may serve more generally as a molecular language to allow for intricate adjustments and the coordination of splicing patterns from different genes.

Many genes are alternatively spliced, but the signals that regulate the process are unclear. These authors have found a sequence motif that appears to function at many alternatively spliced genes

A molecular code for splicing silencing: configurations of guanosine-rich motifs

Alternative pre-mRNA splicing is frequently used to expand the protein-coding capacity of genomes, and to regulate gene expression at the post-transcriptional level. It is a significant challenge to decipher the molecular language of tissue-specific splicing because the inherent flexibility of these mechanisms is specified by numerous short sequence motifs distributed in introns and exons. In the present study, we employ the glutamate NMDA (N-methyl-D-aspartate) R1 receptor (GRIN1) transcript as a model system to identify the molecular determinants for a brain region-specific exon silencing mechanism. We identify a set of guanosine-rich motifs that function co-operatively to regulate the CI cassette exon in a manner consistent with its in vivo splicing pattern. Whereas hnRNP (heterogeneous nuclear ribonucleoprotein) A1 mediates silencing of the CI cassette exon in conjunction with the guanosine-rich motifs, hnRNP H functions as an antagonist to silencing. Genome-wide analysis shows that, while this motif pattern is rarely present in human and mouse exons, those exons for which the pattern is conserved are generally found to be skipped exons. The identification of a similar arrangement of guanosine-rich motifs in transcripts of the hnRNP H family of splicing factors has implications for their co-ordinate regulation at the level of splicing.

SpliceInfo: an information repository for mRNA alternative splicing in human genome

We have developed an information repository named SpliceInfo to collect the occurrences of the four major alternative-splicing (AS) modes in human genome; these include exon skipping, 5′-alternative splicing, 3′-alternative splicing and intron retention. The dataset is derived by comparing the nucleotide and protein sequences available for a given gene for evidence of AS. Additional features such as the tissue specificity of the mRNA, the protein domain contained by exons, the GC-ratio of exons, the repeats contained within the exons, and the Gene Ontology are annotated computationally for each exonic region that is alternatively spliced. Motivated by a previous investigation of AS-related motifs such as exonic splicing enhancer and exonic splicing silencer, this resource also provides a means of identifying motifs candidates and this should help to identify potential regulatory mechanisms within a particular exonic sequence set and its two flanking intronic sequence sets. This is carried out using motif discovery tools to identify motif candidates related to alternative splicing regulation and together with a secondary structure prediction tool, will help in the identification of the structural properties of such regulatory motifs. The integrated resource is now available on http://SpliceInfo.mbc.NCTU.edu.tw/.

DG-CST (Disease Gene Conserved Sequence Tags), a database of human-mouse conserved elements associated to disease genes

The identification and study of evolutionarily conserved genomic sequences that surround disease-related genes is a valuable tool to gain insight into the functional role of these genes and to better elucidate the pathogenetic mechanisms of disease. We created the DG-CST (Disease Gene Conserved Sequence Tags) database for the identification and detailed annotation of human–mouse conserved genomic sequences that are localized within or in the vicinity of human disease-related genes. CSTs are defined as sequences that show at least 70% identity between human and mouse over a length of at least 100 bp. The database contains CST data relative to over 1088 genes responsible for monogenetic human genetic diseases or involved in the susceptibility to multifactorial/polygenic diseases. DG-CST is accessible via the internet at http://dgcst.ceinge.unina.it/ and may be searched using both simple and complex queries. A graphic browser allows direct visualization of the CSTs and related annotations within the context of the relative gene and its transcripts.

A germline DNA polymorphism enhances alternative splicing of the KLF6 tumor suppressor gene and is associated with increased prostate cancer risk

Prostate cancer is a leading and increasingly prevalent cause of cancer death in men. Whereas family history of disease is one of the strongest prostate cancer risk factors and suggests a hereditary component, the predisposing genetic factors remain unknown. We first showed that KLF6 is a tumor suppressor somatically inactivated in prostate cancer and since then, its functional loss has been further established in prostate cancer cell lines and other human cancers. Wild-type KLF6, but not patient-derived mutants, suppresses cell growth through p53-independent transactivation of p21. Here we show that a germline KLF6 single nucleotide polymorphism, confirmed in a tri-institutional study of 3,411 men, is significantly associated with an increased relative risk of prostate cancer in men, regardless of family history of disease. This prostate cancer-associated allele generates a novel functional SRp40 DNA binding site and increases transcription of three alternatively spliced KLF6 isoforms. The KLF6 variant proteins KLF6-SV1 and KLF6-SV2 are mislocalized to the cytoplasm, antagonize wtKLF6 function, leading to decreased p21 expression and increased cell growth, and are up-regulated in tumor versus normal prostatic tissue. Thus, these results are the first to identify a novel mechanism of self-encoded tumor suppressor gene inactivation and link a relatively common single nucleotide polymorphism to both regulation of alternative splicing and an increased risk in a major human cancer.

AHSG gene variant is associated with leanness among Swedish men

Alpha(2) Heremans-Schmid glycoprotein (AHSG) is a plasma protein inhibiting the activity of the insulin receptor tyrosine kinase. Ahsg knock-out mice have increased insulin sensitivity and are resistant to diet-induced obesity. We hypothesized that functional variants of the AHSG gene segregating in the human population would reflect variation in body mass index (BMI). We genotyped 356 overweight or obese (BMI: 37.2 [25.0-66.5] kg/m(2)) and 148 lean (BMI: 23.7 [23.4-24.9] kg/m(2)) otherwise healthy Swedish men for three non-synonymous single-nucleotide polymorphisms (SNPs) within exon 6 (rs4917) and exon 7 (rs4918 and Arg299Cys) and one SNP in intron 1 (rs2593813) of the AHSG gene. The G/G genotype for rs2593813 was more common among lean than among obese and overweight individuals (odds ratio = 2.01, P = 0.009), whereas rs2593813 was in strong linkage disequilibrium (|D'| > or = 0.97) with rs4917 and rs4918. Homozygosity for the rs2593813:G-rs4917:Met-rs4918:Ser haplotype conferred an increased risk for leanness (odds ratio=1.90, P = 0.027). rs4917:Met and rs4918:Ser have previously been associated with lower AHSG protein level. A common variant of AHSG, previously associated with a lower AHSG protein level, is thus more common among lean than obese and overweight men, supporting the results from Ahsg knock-out mice, namely, that AHSG modulates body mass.

SIMPLEmutations in Charcot-Marie-Tooth disease and the potential role of its protein product in protein degradation

Charcot-Marie-Tooth (CMT) disease is a clinically and genetically heterogeneous group of inherited peripheral neuropathies characterized by progressive weakness and atrophy of distal limb muscles. Recently, SIMPLE/LITAF was shown to be responsible for an autosomal dominant demyelinating form of CMT linked to 16p (CMT1C). Although two transcripts encoding different proteins (SIMPLE and LITAF) have been reported from the same gene, we could not confirm the existence of LITAF. Here we show that the LITAF transcript appears to result from a DNA sequencing error. We screened the SIMPLE gene for mutations in a cohort of 192 patients with CMT or related neuropathies, each of whom tested negative for other known genetic causes of CMT. In 16 unrelated CMT families we identified nine different nucleotide variations in SIMPLE that were not detected in control chromosomes. SIMPLE mutations can occur de novo, associated with sporadic CMT1 and may convey both demyelinating and axonal forms. Bioinformatics analyses and other observations of SIMPLE suggest that 1) it could be a member of the RING finger motif-containing subfamily of E3 ubiquitin ligases that are associated with the ubiquitin-mediated proteasome processing pathway, 2) it could interact through its PPXY motifs with a WW domain containing protein, for instance with NEDD4, an E3 ubiquitin ligase, and 3) it could interact through the PSAP motif with TSG10, a protein associated with endosomal multivesicular protein sorting. Since both SIMPLE and Hrs are endosomal proteins and have both PPXY and P(S/T)AP motifs, we hypothesize that SIMPLE, like Hrs, is potentially a clathrin adaptor aiding in the retention of ubiquitinated proteins on to the endosomes. Thus the potential E3 ubiquitin ligase activity of SIMPLE, alteration in its interactions with NEDD4 or TSG101, or changes in its properties as a clathrin coat adaptor may underlie the pathogenesis of Charcot-Marie-Tooth disease.

Proteomic analysis of reaper 5' untranslated region-interacting factors isolated by tobramycin affinity-selection reveals a role for La antigen in reaper mRNA translation

Translational control is a key step in gene expression regulation during apoptosis. To understand the mechanisms of mRNA translation of a pro-apoptotic gene, reaper (rpr), we adapted the tobramycin-aptamer technique described by Hartmuth et al. (Proc. Natl. Acad. Sci. USA 2002, 99, 16719-16724) for the analysis of proteins interacting with rpr 5' untranslated region (UTR). We assembled ribonucleoprotein complexes in vitro using translation extracts derived from Drosophila embryos and purified the RNA-protein complexes for mas spectrometry analysis. We identified the proteins bound to the 5' UTR of rpr. One of them, the La antigen, was validated by RNA-crosslinking experiments using recombinant protein and by the translation efficiency of reporter mRNAs in Drosophila cells after RNAinterference experiments. Our data provide evidence of the involvement of La antigen in the translation of rpr and set a protocol for purification of tagged-RNA-protein complexes from cytoplasmic extracts.

Nova autoregulation reveals dual functions in neuronal splicing

The Nova family of neuron-specific RNA-binding proteins were originally identified as targets in an autoimmune neurologic disease characterized by failure of motor inhibition. Nova-1 regulates alternative splicing of pre-mRNAs encoding the inhibitory neurotransmitter receptor subunits GABA(A)Rgamma2 and GlyRalpha2 by directly binding intronic elements, resulting in enhancement of exon inclusion. Here we identify exon E4 in the Nova-1 pre-mRNA itself, encoding a phosphorylated protein domain, as an additional target of Nova-dependent splicing regulation in the mouse spinal cord. Nova binding to E4 is necessary and sufficient for Nova-dependent exon exclusion. E4 harbors five repeats of the known Nova-binding tetranucleotide YCAY and mutation of these elements destroys Nova-dependent regulation. Furthermore, swapping of the sites from Nova-1 and GABA(A)Rgamma2 indicates that the ability of Nova to enhance or repress alternative exon inclusion is dependent on the position of the Nova-binding element within the pre-mRNA. These studies demonstrate that in addition to its previously described role as a splicing activator, Nova autoregulates its own expression by acting as a splicing repressor.

Hrp59, an hnRNP M protein in Chironomus and Drosophila, binds to exonic splicing enhancers and is required for expression of a subset of mRNAs

Here, we study an insect hnRNP M protein, referred to as Hrp59. Hrp59 is relatively abundant, has a modular domain organization containing three RNA-binding domains, is dynamically recruited to transcribed genes, and binds to premRNA cotranscriptionally. Using the Balbiani ring system of Chironomus, we show that Hrp59 accompanies the mRNA from the gene to the nuclear envelope, and is released from the mRNA at the nuclear pore. The association of Hrp59 with transcribed genes is not proportional to the amount of synthesized RNA, and in vivo Hrp59 binds preferentially to a subset of mRNAs, including its own mRNA. By coimmunoprecipitation of Hrp59–RNA complexes and microarray hybridization against Drosophila whole-genome arrays, we identify the preferred mRNA targets of Hrp59 in vivo and show that Hrp59 is required for the expression of these target mRNAs. We also show that Hrp59 binds preferentially to exonic splicing enhancers and our results provide new insights into the role of hnRNP M in splicing regulation.

Screening of the ARX gene in 682 retarded males

The newly identified gene, ARX, when mutated has been shown to cause both syndromic and nonsyndromic forms of mental retardation. It seems that the less severe forms are due to polyalanine expansions and missense mutations in the gene. We screened 682 developmentally retarded males for polyalanine expansions in ARX in order to examine the contribution of ARX mutations to the causes of developmental retardation. We also reinvestigated 11 putative MRX and three MR families where no cause of mental retardation had been found, by mutational analysis of ARX. Mutational analysis was also performed in 11 probands with autism from families with two or more affected males. We find that previously described polyalanine expansions of ARX are not a common cause of mental retardation.

Viable adenovirus vaccine prototypes: high-level production of a papillomavirus capsid antigen from the major late transcriptional unit

Safe, effective, orally delivered, live adenovirus vaccines have been in use for three decades. Recombinant derivatives of the live adenovirus vaccines may prove an economical alternative to current vaccines for a variety of diseases. To explore that possibility, we constructed a series of recombinants that express the major capsid protein (L1) of canine oral papillomavirus (COPV), a model for mucosal human papillomavirus (HPV) infection. Vaccination with virus-like particles (VLPs) composed of recombinant HPV L1 completely prevents persistent HPV infection [Koutsky, L. A., Ault, K. A., Wheeler, C. M., Brown, D. R., Barr, E., Alvarez, F. B., Chiacchierini, L. M. & Jansen, K. U. (2002) N. Engl. J. Med. 347, 1645-1651], suggesting that L1 expressed from recombinant adenoviruses might provide protective immunity. In our recombinants, COPV L1 is incorporated into adenovirus late region 5 (Ad L5) and is expressed as a member of the adenoviral major late transcriptional unit (MLTU). COPV L1 production by the most prolific recombinant is comparable to that of the most abundant adenoviral protein, hexon. COPV L1 production by recombinants is influenced by Ad L5 gene order, the specific mRNA processing signals associated with COPV L1, and the state of a putative splicing inhibitor in the COPV L1 gene. Recombinant COPV L1 protein assembles into VLPs that react with an antibody specific for conformational epitopes on native COPV L1 protein that correlate with protection in vivo. The designs of these recombinants can be applied directly to the production of recombinants appropriate for assessing immunogenicity and protective efficacy in animal models and in human trials.

Functional consequences of PRODH missense mutations

PRODH maps to 22q11 in the region deleted in the velocardiofacial syndrome/DiGeorge syndrome (VCFS/DGS) and encodes proline oxidase (POX), a mitochondrial inner-membrane enzyme that catalyzes the first step in the proline degradation pathway. At least 16 PRODH missense mutations have been identified in studies of type I hyperprolinemia (HPI) and schizophrenia, 10 of which are present at polymorphic frequencies. The functional consequences of these missense mutations have been inferred by evolutionary conservation, but none have been tested directly. Here, we report the effects of these mutations on POX activity. We find that four alleles (R185Q, L289M, A455S, and A472T) result in mild (<30%), six (Q19P, A167V, R185W, D426N, V427M, and R431H) in moderate (30%-70%), and five (P406L, L441P, R453C, T466M, and Q521E) in severe (>70%) reduction in POX activity, whereas one (Q521R) increases POX activity. The POX encoded by one severe allele (T466M) shows in vitro responsiveness to high cofactor (flavin adenine dinucleotide) concentrations. Although there is limited information on plasma proline levels in individuals of known PRODH genotype, extant data suggest that severe hyperprolinemia (>800 microM) occurs in individuals with large deletions and/or PRODH missense mutations with the most-severe effect on function (L441P and R453C), whereas modest hyperprolinemia (300-500 microM) is associated with PRODH alleles with a moderate reduction in activity. Interestingly, three of the four alleles associated with or found in schizophrenia (V427M, L441P, and R453C) resulted in severe reduction of POX activity and hyperprolinemia. These observations plus the high degree of polymorphism at the PRODH locus are consistent with the hypothesis that reduction in POX function is a risk factor for schizophrenia.

Enhancer-dependent splicing of FGFR1 alpha-exon is repressed by RNA interference-mediated down-regulation of SRp55

The FGFR1 gene transcript is alternatively processed to produce functionally different receptor forms. Previously, we identified a 69-nucleotide exonic splicing enhancer (ESE) required for alpha-exon inclusion in JEG3 cells. In the present study, we found that this sequence is composed of three independent elements, two smaller ESE sequences flanking an exonic splicing silencer sequence. Ultraviolet cross-linking and immunoprecipitation identified ESE-specific binding of the splicing regulator SRp55. A RNA interference-mediated decrease in SRp55 confirmed the significance of this interaction. There was a 6- to 14-fold decrease in exon inclusion on ablation of SRp55. In SNB19 glioblastoma cells, which normally skip this exon, SRp55 was also demonstrated to play a role in exon inclusion after the removal of intronic splicing silencer sequences. These observations indicate that SRp55 plays a major role in maintaining normal FGFR1 alpha-exon inclusion, which is subject to dominant intronic splicing silencer-mediated and exonic splicing silencer-mediated inhibition in SNB19 cells.

Tau exons 2 and 10, which are misregulated in neurodegenerative diseases, are partly regulated by silencers which bind a SRp30c.SRp55 complex that either recruits or antagonizes htra2beta1

Tau is a microtubule-associated protein whose transcript undergoes complex regulated splicing in the mammalian nervous system. Exon 2 modulates the tau N-terminal domain, which interacts with the axonal membrane. Exon 10 codes for a microtubule binding domain, increasing the affinity of tau for microtubules. Both exons are excluded from fetal brain, but their default behavior is inclusion, suggesting that silencers are involved in their regulation. Exon 2 is significantly reduced in myotonic dystrophy type 1, whose symptoms include dementia. Mutations that affect exon 10 splicing cause frontotemporal dementia (FTDP). In this study, we investigated three regulators of exon 2 and 10 splicing: serine/arginine-rich (SR) proteins SRp55, SRp30c, and htra2beta1. The first two inhibit both exons; htra2beta1 inhibits exon 2 but activates exon 10. By deletion analysis, we identified splicing silencers located at the 5' end of each exon. Furthermore, we demonstrated that SRp30c and SRp55 bind to both silencers and to each other. In exon 2, htra2beta1 binds to the inhibitory heterodimer through its RS1 domain but not to exon 2, whereas in exon 10 the heterodimer may sterically interfere with htra2beta1 binding to a purine-rich enhancer (defined by FTDP mutation E10-Delta5 = Delta280K) directly downstream of the silencer. Increased exon 10 inclusion in FTDP mutant ENH (N279K) may arise from abolishing SRp30c binding. Also, htra2beta3, a naturally occurring variant of htra2beta1, no longer inhibits exon 2 splicing but can partially rescue splicing of exon 10 in FTDP mutation E10-Delta5. This work provides interesting insights into the splicing regulation of the tau gene.

The splicing regulatory element, UGCAUG, is phylogenetically and spatially conserved in introns that flank tissue-specific alternative exons

Previous studies have identified UGCAUG as an intron splicing enhancer that is frequently located adjacent to tissue-specific alternative exons in the human genome. Here, we show that UGCAUG is phylogenetically and spatially conserved in introns that flank brain-enriched alternative exons from fish to man. Analysis of sequence from the mouse, rat, dog, chicken and pufferfish genomes revealed a strongly statistically significant association of UGCAUG with the proximal intron region downstream of brain-enriched alternative exons. The number, position and sequence context of intronic UGCAUG elements were highly conserved among mammals and in chicken, but more divergent in fish. Control datasets, including constitutive exons and non-tissue-specific alternative exons, exhibited a much lower incidence of closely linked UGCAUG elements. We propose that the high sequence specificity of the UGCAUG element, and its unique association with tissue-specific alternative exons, mark it as a critical component of splicing switch mechanism(s) designed to activate a limited repertoire of splicing events in cell type-specific patterns. We further speculate that highly conserved UGCAUG-binding protein(s) related to the recently described Fox-1 splicing factor play a critical role in mediating this specificity.

Alternative splicing induced by nonsense mutations in the immunoglobulin mu VDJ exon is independent of truncation of the open reading frame

In addition to triggering nonsense-mediated mRNA decay (NMD), premature translation-termination codons (PTCs) frequently induce alternative splicing, an observation referred to as nonsense-associated alternative splicing (NAS). In many cases, NAS is induced because the nonsense mutation alters a splicing signal, such as inactivating an exonic splicing enhancer. However, for a few genes, NAS was reported to be PTC specific, implying that a translation signal could influence splicing. Here, we investigated whether production of a previously undetected alternatively spliced transcript from immunoglobulin mu (Ig-mu) depends on premature termination of the open reading frame. We show that PTCs at different positions in the VDJ exon of an Ig-mu minigene activate usage of an alternative 3' splice site, generating an alternative transcript that lacks the initial PTC and a previously identified NMD-promoting element (NPE), but contains new PTCs because of a frame shift. Corroborating the importance of the NPE for maximal NMD response, the alternative transcript is only moderately down-regulated by NMD. We further demonstrate that NAS of Ig-mu minigene transcripts is not PTC specific. This finding, together with our results that contradict the previously reported frame dependence of TCR-beta NAS, challenges the idea that cells might possess mechanisms that would allow regulation of splice site selection in response to premature termination of the ORF.

Evidence for the regulation of alternative splicing via complementary DNA sequence repeats

MOTIVATION

While the mechanism for regulating alternative splicing is poorly understood, secondary structure has been shown to be integral to this process. Due to their propensity for forming complementary hairpin loops and their elevated mutation rates, tandem repeated sequences have the potential to influence splicing regulation.

RESULTS

An analysis of human intronic sequences reveals a strong correlation between alternative splicing and the prevalence of mono- through hexanucleotide tandem repeats that may engage in complementary pairing in introns that flank alternatively spliced exons. While only 44% of the 18 173 genes in the Human Alternative Splicing Database are known to be alternatively spliced, they contain 84% of the 694 237 intronic complementary repeat pairs. Significantly, the normalized frequency and distribution of repeat sequences, independent of their potential for pairing, are indistinguishable between alternatively spliced and non-alternatively spliced genes. Thus, the increased prevalence of repeats with pairing potential in alternatively spliced genes is not merely a consequence of more repeats or repeat composition bias. These results suggest that complementary repeats may play a role in the regulation of alternative splicing.

CONTACT

harold.garner@utsouthwestern.edu.

Functional significance of a deep intronic mutation in the ATM gene and evidence for an alternative exon 28a

Screening for ATM mutations is usually performed using genomic DNA as a template for PCR amplification across exonic regions, with the consequence that deep intronic sequences are not analyzed. Here we report a novel pseudoexon-retaining deep intronic mutation (IVS28-159A>G; g.75117A>G based on GenBank U82828.1) in a patient with ataxia-telangiectasia (A-T), as well as the identification of a previously unrecognized alternative exon in the ATM gene (exon 28a) expressed in lymphoblastoid cell lines (LCL) derived from normal individuals. cDNA analysis using the A-T patient's LCL showed the retention of two aberrant intronic segments of 112 and 190 nt between exons 28 and 29. Minigenes were constructed to determine the functional significance of two genomic changes in the region of aberrant splicing: IVS28-193C>T (g.75083C>T) and IVS28-159A>G, revealing that: 1) the first is a polymorphism; 2) IVS28-159A>G weakens the 5' splice site of the alternative exon 28a and activates a cryptic 5' splice site (ss) 83 nt downstream; and 3) wild-type constructs also retain a 29-nt segment (exon 28a) as part of both the 112- and 190-nt segments. Maximum entropy estimates of ss strengths corroborate the cDNA and minigene findings. Such mutations may prove relevant in planning therapy that targets specific splicing aberrations.

Cytochrome P450 oxidoreductase gene mutations and Antley-Bixler syndrome with abnormal genitalia and/or impaired steroidogenesis: molecular and clinical studies in 10 patients

We report on molecular and clinical findings in 10 Japanese patients (four males and six females) from eight families (two pairs of siblings and six isolated cases) with Antley-Bixler syndrome accompanied by abnormal genitalia and/or impaired steroidogenesis. Direct sequencing was performed for all the 15 exons of cytochrome P450 oxidoreductase gene (POR), showing two missense mutations (R457H and Y578C), a 24-bp deletion mutation resulting in loss of nine amino acids and creation of one amino acid (L612_W620delinsR), a single bp insertion mutation leading to frameshift (I444fsX449), and a silent mutation (G5G). R457H has previously been shown to be a pathologic mutation, and computerized modeling analyses indicated that the 15A>G for G5G could disturb an exonic splicing enhancer motif, and the remaining three mutations should affect protein conformations. Six patients were compound heterozygotes, and three patients were R457H homozygotes; no mutation was identified on one allele of the remaining one patient. Clinical findings included various degrees of skeletal features, such as brachycephaly, radiohumeral synostosis, and digital joint contractures in patients of both sexes, normal-to-poor masculinization during fetal and pubertal periods in male patients, virilization during fetal life and poor pubertal development without worsening of virilization in female patients, and relatively large height gain and delayed bone age from the pubertal period in patients of both sexes, together with maternal virilization during pregnancy. Blood cholesterol was grossly normal, and endocrine studies revealed defective CYP17A1 and CYP21A2 activities. The results suggest that Antley-Bixler syndrome with abnormal genitalia and/or impaired steroidogenesis is caused by POR mutations, and that clinical features are variable and primarily explained by impaired activities of POR-dependent CYP51A1, CYP17A1, CYP21A2, and CYP19A1.

Identification of two further splice variants of GABABR1 characterizes the conserved micro-exon 4 as a hot spot for regulated splicing in the rat brain

Inhibitory neurotransmission in the mammalian brain is principally mediated by gamma-aminobutyric acid (GABA) acting through different subtypes of cell membrane GABA receptor (GABAR). The expression of one GABAR gene, GABABR1, is distinguished by the expression of multiple splice variants that encode different isoforms of the receptor. In the present study, we have identified two novel GABABR1 variants, GABABR1h (R1h) and GABABR1i (R1i), which appear to arise from alternative splicing of the GABABR1 gene. The expression of R1h and R1i is differentially regulated in brain and peripheral tissues, but expression is not altered in the brain of a genetic model of absence epilepsy (GAERS rat [genetic absence epilepsy rat from Strasbourg]). Both the R1h and R1i variants exhibit a novel 80-bp insert downstream of exon 4 that is flanked by consensus splice sites, and both encode C-terminal-truncated proteins. The new insight into the family of GABABR1 variants gained from this study identifies exon 4 as a preferred locus, or hot spot for regulated splicing in the GABABR1 gene. This finding correlates with the micro-exonic nature of exon 4 (21 bp). Bioinformatic analysis of micro-exon 4 and its flanking pre-mRNA sequences has revealed multiple, potentially competitive, exonic splicing enhancers that provide a mechanistic basis for the preponderance of alternative splicing events at this locus. Conservation of GABABR1 micro-exon 4 across species suggests a conserved functional role, facilitating either N-terminal protein production or post-transcriptional gene regulation through regulated splicing coupled to transcript decay.

Alpha-synuclein missense and multiplication mutations in autosomal dominant Parkinson's disease

Missense mutations and genomic multiplications of the alpha-synuclein gene (SNCA) have been linked to autosomal dominant familial Parkinson's disease. We screened 50 probands of families with autosomal dominant parkinsonism for alpha-synuclein mutations by exon sequencing. No known or novel mutations were found. We also analyzed the genomic DNA for multiplications of the SNCA locus using multiplex panels of microsatellite markers. All samples were diploid with two normal copies of the SNCA locus. Hence, alpha-synuclein missense mutations and SNCA genomic multiplications remain a rare cause of disease.

Genetic mutations in surfactant protein C are a rare cause of sporadic cases of IPF

BACKGROUND

While idiopathic pulmonary fibrosis (IPF) is one of the most common forms of interstitial lung disease, the aetiology of IPF is poorly understood. Familial cases of pulmonary fibrosis suggest a genetic basis for some forms of the disease. Recent reports have linked genetic mutations in surfactant protein C (SFTPC) with familial forms of pulmonary fibrosis, including one large family in which a number of family members were diagnosed with usual interstitial pneumonitis (UIP), the pathological correlate to IPF. Because of this finding in familial cases of pulmonary fibrosis, we searched for SFTPC mutations in a cohort of sporadic cases of UIP and non-specific interstitial pneumonitis (NSIP).

METHODS

The gene for SFTPC was sequenced in 89 patients diagnosed with UIP, 46 patients with NSIP, and 104 normal controls.

RESULTS

Ten single nucleotide polymorphisms in the SFTPC sequence were found in IPF patients and not in controls. Only one of these created an exonic change resulting in a change in amino acid sequence. In this case, a T to C substitution resulted in a change in amino acid 73 of the precursor protein from isoleucine to threonine. Of the remaining polymorphisms, one was in the 5' UTR, two were exonic without predicted amino acid sequence changes, and six were intronic. One intronic mutation suggested a potential enhancement of a splicing site.

CONCLUSIONS

Mutations in SFTPC are identified infrequently in this patient population. These findings indicate that SFTPC mutations do not contribute to the pathogenesis of IPF in the majority of sporadic cases.

An erythroid differentiation-specific splicing switch in protein 4.1R mediated by the interaction of SF2/ASF with an exonic splicing enhancer

Protein 4.1R is a vital component of the red blood cell membrane cytoskeleton. Promotion of cytoskeletal junctional complex stability requires an erythroid differentiation stage-specific splicing switch promoting inclusion of exon 16 within the spectrin/actin binding domain. We showed earlier that an intricate combination of positive and negative RNA elements controls exon 16 splicing. In this report, we further identified 3 putative exonic splicing enhancers within exon 16 and investigated the function of the sequence CAGACAT in the regulation of exon 16 splicing. Mutation of these sequences leads to increased exclusion of exon 16 in both in vivo and in vitro splicing assays, indicating that CAGACAT is a functional exonic splicing enhancer. UV cross-linking further detects an approximately 33-kDa protein that specifically binds to the CAGACAT-containing transcript. An anti-SF2/ASF antibody specifically immunoprecipitates the approximately 33-kDa protein. Furthermore, SF2/ASF stimulates exon 16 inclusion in both in vitro complementation assays and minigene-transfected mouse erythroleukemia cells (MELCs). Finally, SF2/ASF expression is up-regulated and correlates with exon 16 inclusion in differentiated MELCs. These results suggest that increased splicing factor 2/alternative splicing factor (SF2/ASF) expression in differentiated mouse erythroleukemia mediates a differentiation stage-specific exon 16 splicing switch through its interaction with the exonic splicing enhancer.

A novel type 2A von Willebrand factor mutation located at the last nucleotide of exon 26 (3538G>A) causes skipping of 2 nonadjacent exons

In this manuscript, we describe a case of type 2A von Willebrand disease (VWD) caused by the novel heterozygous G&gt;A transition at nucleotide 3538, which should result in the putative, nonconservative substitution of G1180R. This mutation was reproduced by site-directed mutagenesis; however, the recombinant mutant protein was efficiently secreted from cells and assembled correctly into multimers. Because the substitution is located at the last nucleotide of exon 26, the patient's platelet von Willebrand factor (VWF) mRNA was analyzed and 3 transcripts were observed: the normal transcript without the 3538G&gt;A transition, a transcript with the in-frame deletion of exon 26, and a transcript with the in-frame deletions of exons 23 and 26. These deletion VWF cDNA constructs were created and the resulting recombinant proteins were analyzed following transfection into COS-7 cells. Cotransfection results demonstrate that the exon-skipped transcripts led to intracellular retention, and the levels of VWF antigen (VWF:Ag) produced by these constructs were as follows: del23/26&lt;del26&lt; G1180R≤wild type. The homozygous exon-skipped transcripts show the presence of only the lowest molecular weight multimers. The G&gt;A transition at nucleotide (nt) 3538 does not result in the expression of the G1180R missense mutation, but rather leads to exon skipping, which is the pathogenic basis of the patient's phenotype. This is the first report of a coding region mutation resulting in the skipping of 2 nonadjacent exons. (Blood. 2004;104:2739-2745)

Over-representation of exonic splicing enhancers in human intronless genes suggests multiple functions in mRNA processing

The human transcriptome is constituted of a great majority of intron-containing and a minority of intron-lacking mRNAs; given the different processing these transcripts undergo, they are expected to carry, intermingled with coding properties, very different editing information. Here we applied a computational approach to compare intronless and intron-containing coding sequences. Hexamer composition comparison allowed the definition of over- and under-represented motifs in intronless genes; surprisingly, experimental testing revealed that intron-lacking coding sequences are enriched rather than depleted in elements with splicing enhancement ability. Similarly, we show evidence that intronless transcripts display a significantly higher frequency of both shuttling and non-shuttling SR protein binding sites compared to intron-containing sequences. These observations suggest that SR proteins (and possibly other splicing factors) play a role in cellular processes distinct from splicing.

The role of common single-nucleotide polymorphisms on exon 9 and exon 12 skipping in nonmutated CFTR alleles

Classic cystic fibrosis (CF) is caused by two loss-of-function mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene, whereas patients with nonclassic CF have at least one copy of a mutant gene that retains partial function of the CFTR protein. In addition, there are several other phenotypes associated with CFTR gene mutations, such as idiopathic chronic pancreatitis. In CFTR-associated disorders and in nonclassic CF, often only one CFTR mutation or no CFTR mutations can be detected. In this study, we screened 23 patients with CFTR-associated disorders for CFTR mutations by complete gene testing and quantitative transcript analysis. Mutations were found in 10 patients. In cells from respiratory epithelium, we detected aberrant splicing of CFTR mRNA in all investigated individuals. We observed a highly significant association between the presence of coding single-nucleotide polymorphisms (coding SNPs, or cSNPs) and increased skipping of exon 9 and 12. This association was found both in patients and in normal individuals carrying the same cSNPs. The cSNPs c.1540A>G, c.2694T>G, and c.4521G>A may have affected pre-mRNA splicing by changing regulatory sequence motifs of exonic splice enhancers, leading to lower amounts of normal transcripts. The analysis of CFTR exons indicated that less frequent and weak exonic splicing enhancer (ESE) motifs make exon 12 vulnerable to skipping. The number of splice variants in individuals with cSNPs was similar to previously reported values for the T5 allele, suggesting that cSNPs may enhance susceptibility to CFTR related diseases. In addition, cSNPs may be responsible for variation in the phenotypic expression of CFTR mutations. Quantitative approaches rather than conventional genomic analysis are required to interpret the role of cSNPs.

Branch site haplotypes that control alternative splicing

We show that the allele-dependent expression of transcripts encoding soluble HLA-DQbeta chains is determined by branchpoint sequence (BPS) haplotypes in DQB1 intron 3. BPS RNAs associated with low inclusion of the transmembrane exon in mature transcripts showed impaired binding to splicing factor 1 (SF1), indicating that alternative splicing of DQB1 is controlled by differential BPS recognition early during spliceosome assembly. We also demonstrate that naturally occurring human BPS point mutations that alter splicing and lead to recognizable phenotypes cluster in BP and in position -2 relative to BP, implicating impaired SF1-BPS interactions in disease-associated BPS substitutions. Coding DNA variants produced smaller fluctuations of exon inclusion levels than random exonic substitutions, consistent with a selection against coding mutations that alter their own exonization. Finally, proximal splicing in this multi-allelic reporter system was promoted by at least seven SR proteins and repressed by hnRNPs F, H and I, supporting an extensive antagonism of factors balancing the splice site selection. These results provide the molecular basis for the haplotype-specific expression of soluble DQbeta, improve prediction of intronic point mutations and indicate how extraordinary, selection-driven DNA variability in HLA affects pre-mRNA splicing.

Mutations in PITX2 may contribute to cases of omphalocele and VATER-like syndromes

Omphalocele is a congenital anomaly with substantial morbidity. Rieger syndrome, an autosomal dominant disorder, is characterized by craniofacial abnormalities and abdominal wall defects. PITX2 mutations are etiologic in >40% of cases of Rieger syndrome. We demonstrate that the birth prevalence of omphalocele is significantly higher in Rieger syndrome than in the general population, with omphaloceles found in 0.03% in the Iowa newborn population and 4.3% of patients with Rieger syndrome. Our objective was to screen coding and conserved non-coding regions of PITX2 for mutations in 209 patients with omphalocele. We identified remarkable evolutionarily conserved regions by comparing the 3'UTR of Pitx2 in 13 vertebrate and 3 invertebrate species. No mutations changing the amino acid sequence were found within the omphalocele population. In one case of omphalocele with VATER-like additional anomalies, a three nucleotide deletion was found in the 3'UTR. This deletion was not seen in 1,186 controls. Also in the 3'UTR, we identified a single nucleotide polymorphism at a highly conserved residue. Our findings suggest additional studies of PITX2 conserved regions will be valuable. We also screened the omphalocele cases for mutations in exon 5 of the gene FLNA. Mutations in FLNA have been shown to cause a broad range of congenital malformations, including otopalatodigital syndrome type 2 in which a missense mutation occurring in exon 5 of FLNA results in omphalocele as part of the phenotype. We did not find any mutations in exon 5 of FLNA in 179 omphalocele cases studied.

Identification of a human mutation of DMT1 in a patient with microcytic anemia and iron overload

Divalent metal transporter 1 (DMT1) is a transmembrane protein crucial for duodenal iron absorption and erythroid iron transport. DMT1 function has been elucidated largely in studies of the mk mouse and the Belgrade rat, which have an identical single nucleotide mutation of this gene that affects protein processing, stability, and function. These animals exhibit hypochromic microcytic anemia due to impaired intestinal iron absorption, and defective iron utilization in red cell precursors. We report here the first human mutation of DMT1 identified in a female with severe hypochromic microcytic anemia and iron overload. This homozygous mutation in the ultimate nucleotide of exon 12 codes for a conservative E399D amino acid substitution; however, its pre-dominant effect is preferential skipping of exon 12 during processing of pre-messenger RNA (mRNA). The lack of full-length mRNA would predict deficient iron absorption in the intestine and deficient iron utilization in erythroid precursors; however, unlike the animal models of DMT1 mutation, the patient is iron overloaded. This does not appear to be due to up-regulation of total DMT1 mRNA. DMT1 protein is easily detectable by immunoblotting in the patient's duodenum, but it is unclear whether the protein is properly processed or targeted.