High intron sequence conservation across three mammalian orders suggests functional constraints

Several studies have demonstrated high levels of sequence conservation in noncoding DNA compared between two species (e.g., human and mouse), and interpreted this conservation as evidence for functional constraints. If this interpretation is correct, it suggests the existence of a hidden class of abundant regulatory elements. However, much of the noncoding sequence conserved between two species may result from chance or from small-scale heterogeneity in mutation rates. Stronger inferences are expected from sequence comparisons using more than two taxa, and by testing for spatial patterns of conservation in addition to primary sequence similarity. We used a Bayesian local alignment method to compare approximately 10 kb of intron sequence from nine genes in a pairwise manner between human, whale, and seal to test whether the degree and pattern of conservation is consistent with neutral divergence. Comparison of the three sets of conserved gapless pairwise blocks revealed the following patterns: The proportion of identical intron nucleotides averaged 47% in pairwise comparisons and 28% across the three taxa. Proportions of conserved sequence were similar in unique sequence and general mammalian repetitive elements. We simulated sequence evolution under a neutral model using published estimates of substitution rate heterogeneity for noncoding DNA and found pairwise identity at 33% and three-taxon identity at 16% of nucleotide sites. Spatial patterns of primary sequence conservation were also nonrandomly distributed within introns. Overall, segments of intron sequence closer to flanking exons were significantly more conserved than interior intron sequence. This level of intron sequence conservation is above that expected by chance and strongly suggests that intron sequences are playing a larger functional role in gene regulation than previously realized.     

Cooperative assembly of an hnRNP complex induced by a tissue-specific homolog of polypyrimidine tract binding protein

Splicing of the c-src N1 exon in neuronal cells depends in part on an intronic cluster of RNA regulatory elements called the downstream control sequence (DCS). Using site-specific cross-linking, RNA gel shift, and DCS RNA affinity chromatography assays, we characterized the binding of several proteins to specific sites along the DCS RNA. Heterogeneous nuclear ribonucleoprotein (hnRNP) H, polypyrimidine tract binding protein (PTB), and KH-type splicing-regulatory protein (KSRP) each bind to distinct elements within this sequence. We also identified a new 60-kDa tissue-specific protein that binds to the CUCUCU splicing repressor element of the DCS RNA. This protein was purified, partially sequenced, and cloned. The new protein (neurally enriched homolog of PTB [nPTB]) is highly homologous to PTB. Unlike PTB, nPTB is enriched in the brain and in some neural cell lines. Although similar in sequence, nPTB and PTB show significant differences in their properties. nPTB binds more stably to the DCS RNA than PTB does but is a weaker repressor of splicing in vitro. nPTB also greatly enhances the binding of two other proteins, hnRNP H and KSRP, to the DCS RNA. These experiments identify specific cooperative interactions between the proteins that assemble onto an intricate splicing-regulatory sequence and show how this hnRNP assembly is altered in different cell types by incorporating different but highly related proteins.     

The murine gene encoding the highly conserved Sm B protein contains a nonfunctional alternative 3' splice site.

The cDNA and partial genomic nucleotide (nt) sequences were derived for the mouse Sm B polypeptide and compared to the cDNA and genomic sequences encoding human Sm B. The deduced amino acid (aa) sequences from the mouse and human genes are identical with the exception of a single conserved aa substitution, accounting for the ability of anti-Sm antibodies to recognize the Sm polypeptides from a broad range of species. The genomic sequence of mouse B gene is similar to the human B genomic locus that extends from exon 6 to exon 7. These loci include conservation of both 3' alternative splice sites and putative branch points required to process B and B' mRNAs in human cells. However, the nt sequence downstream from the putative distal 3' splice junction and single nt flanking the 3' splice site consensus sequence, differ between mouse and human B. This results in a murine mRNA with a different predicted secondary structure around the distal 3' splice site when compared to humans. Thus, secondary structural constraints in the mRNA or changes in the exon sequence might prevent recognition of this alternative splice site to form B' mRNA in murine tissues.     

Genomic checkpoints for exon 10 usage in the luteinizing hormone receptor type 1 and type 2

Alternative splicing is a hallmark of glycoprotein hormone receptor gene regulation, but its molecular mechanism is unknown. The LH receptor (LHR) gene possesses 11 exons, but exon 10 is constitutively skipped in the New World monkey lineage (LHR type 2), whereas it is constitutively spliced in the human (LHR type 1). This study identifies the regulatory elements of exon 10 usage. Sequencing of genomic marmoset DNA revealed that the cryptic LHR exon 10 is highly homologous to exon 10 from other species and displays intact splice sites. Functional studies using a minigene approach excluded the contribution of intronic, marmoset-specific long interspersed nucleotide-1 elements to exon 10 skipping. Sequencing of the genomic regions surrounding exon 10 from several primate lineages, sequence comparisons including the human and mouse LHR gene, revealed the presence of unique nucleotides at 3'-intronic position -19 and -10 and at position +26 within exon 10 of the marmoset LHR. Exon trap experiments and in vitro mutagenesis of these nucleotides resulted in the identification of a composite regulatory element of splicing consisting of cis-acting elements represented by two polypyrimidine tracts and a trans-acting element within exon 10, which affect the secondary RNA structure. Changes within this complex resulted either in constitutive exon inclusion, constitutive skipping, or alternative splicing of exon 10. This work delineates the molecular pathway leading to intronization of exon 10 in the LHR type 2 and reveals, for the first time, the essential function of regulatory and structural elements involved in glycoprotein hormone receptor splicing.     

Isolation and characterization of the human apolipoprotein A-II gene. Electron microscopic analysis of RNA:DNA hybrids, nucleotide sequence, identification of a polymorphic MspI site, and general structural organization of apolipoprotein genes

Three cloned apolipoprotein A-II genes were isolated from a human genomic cosmid library constructed in our laboratory. An approximately 3-kilobase HindIII insert containing the entire gene was analyzed by RNA:DNA hybridization and electron microscopy. The apo-A-II gene was found to consist of 4 exons and 3 intervening sequences (IVS), and the lengths of each exon and IVS were estimated by direct observation of the hybrids. The entire approximately 3-kilobase HindIII insert was sequenced. The 5' end of the gene was determined by primer extension. The DNA sequence confirms the presence of 4 exons and 3 IVS: exon 1, 34 nucleotides; exon 2, 76 nucleotides; exon 3, 133 nucleotides; exon 4, 230 nucleotides; IVS-I, 169 nucleotides; IVS-II, 299 nucleotides; and IVS-III, 396 nucleotides. A "TATA box" is located at position -29 from the CAP site. A "CAT box" is present at position -78. A "TG" element consisting of (TG)19 is identified at the 3' end of IVS-III. Furthermore, an enhancer core sequence, CTTTCCA, is identified at position -355 in the 5' flanking sequence. At positions -497 to -471 upstream from the CAP site is a stretch of 27 nucleotides that show high homology to stretches of 5' flanking sequences in the apo-C-II, apo-A-I, apo-E, and apo-C-III genes. An Alu dimer sequence is located approximately 300 nucleotides from the 3' end of the gene. Within this Alu sequence, we have identified a polymorphic MspI site. Restriction fragment length polymorphism involving this site has been previously shown to correlate with apo-A-II levels and high density lipoprotein structure. Analysis of conformation by Chou-Fasman analysis and by the helical hydrophobic moment of Eisenberg et al. (Eisenberg, D., Weiss, R. M., and Tergwillager, T. C. (1982). Nature (Lond.) 299, 371-374) indicates that in all of the 5 apolipoproteins characterized at the nucleotide level to date, i.e. apo-C-II, apo-A-II, apo-E, apo-A-I, and apo-C-III, the 2 IVS within the peptide coding regions of the gene tend to occur at regions corresponding to the surface of the polypeptide chain and divide the protein into distinct functional domains.     

Identification of a novel intron and 4 polymorphisms in the gene encoding the gamma subunit of the epithelial sodium channel.

The amiloride-sensitive epithelial sodium channel is a highly selective sodium channel that constitutes the rate-limiting step of sodium reabsorption in distal nephrons. It consists of at least 3 subunits (alpha, beta, and gamma) of similar structure and plays an important role in sodium and fluid homeostasis. Defects of this channel have been critically implicated in Liddle syndrome (pseudoaldosteronism) and pseudohypo-aldosteronism type 1. A sample of 48 individuals from 23 nuclear families was selected from Anhui, China. We sequenced 12 exons and flanking intronic sequences and discovered a new 207-bp intron located in the previously described exon X of Thomas et al. (1996). In addition, 4 novel single nucleotide polymorphisms were identified; 3 were in exon 3 and 1 was in exon 13. Furthermore, 2 base substitutions in exon 13 were present in all the Chinese subjects compared with the published European SCNN1G DNA sequence.     

Human and mouse introns are linked to the same processes and functions through each genome's most frequent non-conserved motifs

We identified the most frequent, variable-length DNA sequence motifs in the human and mouse genomes and sub-selected those with multiple recurrences in the intergenic and intronic regions and at least one additional exonic instance in the corresponding genome. We discovered that these motifs have virtually no overlap with intronic sequences that are conserved between human and mouse, and thus are genome-specific. Moreover, we found that these motifs span a substantial fraction of previously uncharacterized human and mouse intronic space. Surprisingly, we found that these genome-specific motifs are over-represented in the introns of genes belonging to the same biological processes and molecular functions in both the human and mouse genomes even though the underlying sequences are not conserved between the two genomes. In fact, the processes and functions that are linked to these genome-specific sequence-motifs are distinct from the processes and functions which are associated with intronic regions that are conserved between human and mouse. The findings show that intronic regions from different genomes are linked to the same processes and functions in the absence of underlying sequence conservation. We highlight the ramifications of this observation with a concrete example that involves the microsatellite instability gene MLH1.     

Identification of a conserved sequence in the non-coding regions of many human genes.

We have analyzed a sequence of approximately 70 base pairs (bp) that shows a high degree of similarity to sequences present in the non-coding regions of a number of human and other mammalian genes. The sequence was discovered in a fragment of human genomic DNA adjacent to an integrated hepatitis B virus genome in cells derived from human hepatocellular carcinoma tissue. When one of the viral flanking sequences was compared to nucleotide sequences in GenBank, more than thirty human genes were identified that contained a similar sequence in their non-coding regions. The sequence element was usually found once or twice in a gene, either in an intron or in the 5' or 3' flanking regions. It did not share any similarities with known short interspersed nucleotide elements (SINEs) or presently known gene regulatory elements. This element was highly conserved at the same position within the corresponding human and mouse genes for myoglobin and N-myc, indicating evolutionary conservation and possible functional importance. Preliminary DNase I footprinting data suggested that the element or its adjacent sequences may bind nuclear factors to generate specific DNase I hypersensitive sites. The size, structure, and evolutionary conservation of this sequence indicates that it is distinct from other types of short interspersed repetitive elements. It is possible that the element may have a cis-acting functional role in the genome.     

The emergence of alternative 3' and 5' splice site exons from constitutive exons

Alternative 3' and 5' splice site (ss) events constitute a significant part of all alternative splicing events. These events were also found to be related to several aberrant splicing diseases. However, only few of the characteristics that distinguish these events from alternative cassette exons are known currently. In this study, we compared the characteristics of constitutive exons, alternative cassette exons, and alternative 3'ss and 5'ss exons. The results revealed that alternative 3'ss and 5'ss exons are an intermediate state between constitutive and alternative cassette exons, where the constitutive side resembles constitutive exons, and the alternative side resembles alternative cassette exons. The results also show that alternative 3'ss and 5'ss exons exhibit low levels of symmetry (frame-preserving), similar to constitutive exons, whereas the sequence between the two alternative splice sites shows high symmetry levels, similar to alternative cassette exons. In addition, flanking intronic conservation analysis revealed that exons whose alternative splice sites are at least nine nucleotides apart show a high conservation level, indicating intronic participation in the regulation of their splicing, whereas exons whose alternative splice sites are fewer than nine nucleotides apart show a low conservation level. Further examination of these exons, spanning seven vertebrate species, suggests an evolutionary model in which the alternative state is a derivative of an ancestral constitutive exon, where a mutation inside the exon or along the flanking intron resulted in the creation of a new splice site that competes with the original one, leading to alternative splice site selection. This model was validated experimentally on four exons, showing that they indeed originated from constitutive exons that acquired a new competing splice site during evolution.     

RNA helicase p68 (DDX5) regulates tau exon 10 splicing by modulating a stem-loop structure at the 5' splice site

Regulation of tau exon 10 splicing plays an important role in tauopathy. One of the cis elements regulating tau alternative splicing is a stem-loop structure at the 5' splice site of tau exon 10. The RNA helicase(s) modulating this stem-loop structure was unknown. We searched for splicing regulators interacting with this stem-loop region using an RNA affinity pulldown-coupled mass spectrometry approach and identified DDX5/RNA helicase p68 as an activator of tau exon 10 splicing. The activity of p68 in stimulating tau exon 10 inclusion is dependent on RBM4, an intronic splicing activator. RNase H cleavage and U1 protection assays suggest that p68 promotes conformational change of the stem-loop structure, thereby increasing the access of U1snRNP to the 5' splice site of tau exon 10. This study reports the first RNA helicase interacting with a stem-loop structure at the splice site and regulating alternative splicing in a helicase-dependent manner. Our work uncovers a previously unknown function of p68 in regulating tau exon 10 splicing. Furthermore, our experiments reveal functional interaction between two splicing activators for tau exon 10, p68 binding at the stem-loop region and RBM4 interacting with the intronic splicing enhancer region.     

Structure of the human pancreatic cholesterol esterase gene.

The gene for human pancreatic cholesterol esterase consists of 11 exons and 10 introns and is 9.2 kb in length. The last and longest exon (841 nucleotides) is unique to the human gene. Functional amino acids are encoded on separate exons. The leader sequence is encoded by a single exon which carries two additional N-terminal amino acids of the mature functional protein. A positive TATA element is identified 43 nucleotides from the start codon. Pulse-field gel electrophoresis and hybridization with various cDNA probes and direct sequence data revealed the existence of a CEase-like gene. Partial sequence analysis of this gene from a human cosmid library and human genomic DNA showed a premature stop signal in exon 10, shortly after the codon for the active-site histidine. Both the functional gene and the CEase-like gene have a polyadenylation signal in the 3'-untranslated region. Thus, the complex gene structure for this intestinally active enzyme may provide in part a potential molecular explanation for the well-known heterogeneity of the intestinal absorption of cholesterol.     

Splicing factors induce cystic fibrosis transmembrane regulator exon 9 skipping through a nonevolutionary conserved intronic element

In monosymptomatic forms of cystic fibrosis such as congenital bilateral absence of vas deferens, variations in the TG(m) and T(n) polymorphic repeats at the 3' end of intron 8 of the cystic fibrosis transmembrane regulator (CFTR) gene are associated with the alternative splicing of exon 9, which results in a nonfunctional CFTR protein. Using a minigene model system, we have previously shown a direct relationship between the TG(m)T(n) polymorphism and exon 9 splicing. We have now evaluated the role of splicing factors in the regulation of the alternative splicing of this exon. Serine-arginine-rich proteins and the heterogeneous nuclear ribonucleoprotein A1 induced exon skipping in the human gene but not in its mouse counterpart. The effect of these proteins on exon 9 exclusion was strictly dependent on the composition of the TG(m) and T(n) polymorphic repeats. The comparative and functional analysis of the human and mouse CFTR genes showed that a region of about 150 nucleotides, present only in the human intron 9, mediates the exon 9 splicing inhibition in association with exonic regulatory elements. This region, defined as the CFTR exon 9 intronic splicing silencer, is a target for serine-arginine-rich protein interactions. Thus, the nonevolutionary conserved CFTR exon 9 alternative splicing is modulated by the TG(m) and T(n) polymorphism at the 3' splice region, enhancer and silencer exonic elements, and the intronic splicing silencer in the proximal 5' intronic region. Tissue levels and individual variability of splicing factors would determine the penetrance of the TG(m)T(n) locus in monosymptomatic forms of cystic fibrosis.