Alternative splicing of beta-tropomyosin pre-mRNA: multiple cis-elements can contribute to the use of the 5'- and 3'-splice sites of the nonmuscle/smooth muscle exon 6.

We previously found that the splicing of exon 5 to exon 6 in the rat beta-TM gene required that exon 6 first be joined to the downstream common exon 8 (Helfman et al., Genes and Dev. 2, 1627-1638, 1988). Pre-mRNAs containing exon 5, intron 5 and exon 6 are not normally spliced in vitro. We have carried out a mutational analysis to determine which sequences in the pre-mRNA contribute to the inability of this precursor to be spliced in vitro. We found that mutations in two regions of the pre-mRNA led to activation of the 3'-splice site of exon 6, without first joining exon 6 to exon 8. First, introduction of a nine nucleotide poly U tract upstream of the 3'-splice site of exon 6 results in the splicing of exon 5 to exon 6 with as little as 35 nucleotides of exon 6. Second, introduction of a consensus 5'-splice site in exon 6 led to splicing of exon 5 to exon 6. Thus, three distinct elements can act independently to activate the use of the 3'-splice site of exon 6: (1) the sequences contained within exon 8 when joined to exon 6, (2) a poly U tract in intron 5, and (3) a consensus 5'-splice site in exon 6. Using biochemical assays, we have determined that these sequence elements interact with distinct cellular factors for 3'-splice site utilization. Although HeLa cell nuclear extracts were able to splice all three types of pre-mRNAs mentioned above, a cytoplasmic S100 fraction supplemented with SR proteins was unable to efficiently splice exon 5 to exon 6 using precursors in which exon 6 was joined to exon 8. We also studied how these elements contribute to alternative splice site selection using precursors containing the mutually exclusive, alternatively spliced cassette comprised of exons 5 through 8. Introduction of the poly U tract upstream of exon 6, and changing the 5'-splice site of exon 6 to a consensus sequence, either alone or in combination, facilitated the use of exon 6 in vitro, such that exon 6 was spliced more efficiently to exon 8. These data show that intron sequences upstream of an exon can contribute to the use of the downstream 5'-splice, and that sequences surrounding exon 6 can contribute to tissue-specific alternative splice site selection.     

Updating of transposable element annotations from large wheat genomic sequences reveals diverse activities and gene associations

Triticeae species (including wheat, barley and rye) have huge and complex genomes due to polyploidization and a high content of transposable elements (TEs). TEs are known to play a major role in the structure and evolutionary dynamics of Triticeae genomes. During the last 5 years, substantial stretches of contiguous genomic sequence from various species of Triticeae have been generated, making it necessary to update and standardize TE annotations and nomenclature. In this study we propose standard procedures for these tasks, based on structure, nucleic acid and protein sequence homologies. We report statistical analyses of TE composition and distribution in large blocks of genomic sequences from wheat and barley. Altogether, 3.8 Mb of wheat sequence available in the databases was analyzed or re-analyzed, and compared with 1.3 Mb of re-annotated genomic sequences from barley. The wheat sequences were relatively gene-rich (one gene per 23.9 kb), although wheat gene-derived sequences represented only 7.8% (159 elements) of the total, while the remainder mainly comprised coding sequences found in TEs (54.7%, 751 elements). Class I elements [mainly long terminal repeat (LTR) retrotransposons] accounted for the major proportion of TEs, in terms of sequence length as well as element number (83.6% and 498, respectively). In addition, we show that the gene-rich sequences of wheat genome A seem to have a higher TE content than those of genomes B and D, or of barley gene-rich sequences. Moreover, among the various TE groups, MITEs were most often associated with genes: 43.1% of MITEs fell into this category. Finally, the TRIM and copia elements were shown to be the most active TEs in the wheat genome. The implications of these results for the evolution of diploid and polyploid wheat species are discussed. Electronic Supplementary Material Supplementary material is available for this article at     

snRNA interactions at 5' and 3' splice sites monitored by photoactivated crosslinking in yeast spliceosomes.

Splice site recognition and catalysis of the transesterification reactions in the spliceosome are accompanied by a dynamic series of interactions involving conserved or invariant sequences in the spliceosomal snRNAs. We have used site-specific photoactivated crosslinking in yeast spliceosomes to monitor interactions between snRNAs and exon sequences near the 5' and 3' splice sites. The last nucleotide of the 5' exon can be crosslinked to an invariant loop sequence in U5 SnRNA before and after 5' splice site cleavage. The first nucleotide of the 3' exon can also be crosslinked to the same U5 loop sequence, but this contact is only detectable after the first transesterification. These results are in close agreement with earlier data from mammalian splicing extracts, and they are consistent with a model in which U5 snRNA aligns the 5' and 3' exons for the second transesterification. After the first catalytic step of splicing, the first nucleotide of the 3' exon can also crosslink to nt U23 in U2 snRNA. This is one of a cluster of residues in U2-U6 helix I implicated by mutational analysis in the second catalytic step of splicing. The crosslinking data suggest that these residues in U2-U6 helix I are in close proximity to the scissile phosphodiester bond at the 3' splice site prior to the second transesterification. These results constitute the first biochemical evidence for a direct interaction between the 3' splice site and U2 snRNA.     

Comparative analysis of the PCOLCE region in Fugu rubripes using a new automated annotation tool

The Japanese pufferfish Fugu rubripes with a genome of about 400 Mb is becoming increasingly recognized as a vertebrate model organism for comparative gene analysis (see Elgar 1996 for review). We have isolated and sequenced two Fugu cosmids spanning a genomic region of 66 kb containing the Fugu homolog to the human PCOLCE-I (Gl?ckner et al. 1998). We then examined if RUMMAGE-DP, a newly developed analysis tool for gene discovery which was designed for human and mouse genomic DNA, can be used for automatic annotation of Fugu genomic sequence. The exon prediction programs contained in RUMMAGE-DP performed better overall for the human sequence than for the Fugu contig. The GENSCAN program was the only exon prediction programme that performed equally well for both organisms. We show that RUMMAGE-DP is very useful in automatic analysis of Fugu sequences. Comparative analysis of the genomic structure of the PCOLCE-I genes in Fugu and human reveals that the exon/intron structure throughout the protein coding region is almost identical. We defined an additional domain based on the high degree of similarity of 26 aa between mammals and Fugu. The PCOLCE-I protein in both organisms contains two highly conserved CUB domains. Exons 6 and 7 are the only coding exons that differ in length between the two species. We assume that these exons do not code for any catalytic domain of the protein. Analysis of the remaining five Fugu genes within the 66 kb interval revealed no conserved synteny with the corresponding human 7q22 region. Received: 13 October 1998 / Accepted: 25 July 1999     

Report of a chimeric origin of transposable elements in a bovine-coding gene

Despite the wide distribution of transposable elements (TEs) in mammalian genomes, part of their evolutionary significance remains to be discovered. Today there is a substantial amount of evidence showing that TEs are involved in the generation of new exons in different species. In the present study, we searched 22,805 genes and reported the occurrence of TE-cassettes in coding sequences of 542 cow genes using the RepeatMasker program. Despite the significant number (542) of genes with TE insertions in exons only 14 (2.6%) of them were translated into protein, which we characterized as chimeric genes. From these chimeric genes, only the FAST kinase domains 3 (FASTKD3) gene, present on chromosome BTA 20, is a functional gene and showed evidence of the exaptation event. The genome sequence analysis showed that the last exon coding sequence of bovine FASTKD3 is approximately 85% similar to the ART2A retrotransposon sequence. In addition, comparison among FASTKD3 proteins shows that the last exon is very divergent from those of Homo sapiens, Pan troglodytes and Canis familiares. We suggest that the gene structure of bovine FASTKD3 gene could have originated by several ectopic recombinations between TE copies. Additionally, the absence of TE sequences in all other species analyzed suggests that the TE insertion is clade-specific, mainly in the ruminant lineage.     

Genome-Wide Comparative Analyses of Microsatellites in Papaya

Microsatellites, or simple sequence repeats (SSRs), are highly polymorphic and universally distributed in eukaryotes. SSRs have been used extensively as sequence tagged markers in genetic studies. Recently, the functional and evolutionary importance of SSRs has received considerable attention. Here we report the mining and characterization of the SSRs in papaya genome. We analyzed SSRs from 277.4 Mb of whole genome shotgun (WGS) sequences, 51.2 Mb bacterial artificial chromosome (BAC) end sequences (BES), and 13.4 Mb expressed sequence tag (EST) sequences. The papaya SSR density was one SSR per 0.7 kb of DNA sequence in the WGS, which was higher than that in BES and EST sequences. SSR abundance was dramatically reduced as the repeat length increased. According to SSR motif length, dinucleotide repeats were the most common motif in class I, whereas hexanucleotides were the most copious in class II SSRs. The tri- and hexanucleotide repeats of both classes were greater in EST sequences compared to genomic sequences. In class I SSR, AT and AAT were the most frequent motifs in BES and WGS sequences. By contrast, AG and AAG were the most abundant in EST sequences. For SSR marker development, 9,860 primer pairs were surveyed for amplification and polymorphism. Successful amplification and polymorphic rates were 66.6% and 17.6%, respectively. The highest polymorphic rates were achieved by AT, AG, and ATG motifs. The genome wide analysis of microsatellites revealed their frequency and distribution in papaya genome, which varies among plant genomes. This complete set of SSRs markers throughout the genome will assist diverse genetic studies in papaya and related species.     

Design principles for bifunctional targeted oligonucleotide enhancers of splicing

Controlling the patterns of splicing of specific genes is an important goal in the development of new therapies. We have shown that the splicing of a refractory exon, SMN2 exon 7, could be increased in fibroblasts derived from patients with spinal muscular atrophy by using bifunctional targeted oligonucleotide enhancers of splicing (TOES) oligonucleotides that anneal to the exon and contain a 'tail' of enhancer sequences that recruit activating proteins. We show here that there are striking agreements between the effects of oligonucleotides on splicing in vitro and on both splicing and SMN2 protein expression in patient-derived fibroblasts, indicating that the effects on splicing are the major determinant of success. Increased exon inclusion depends on the number, sequence and chemistry of the motifs that bind the activator protein SRSF1, but it is not improved by increasing the strength of annealing to the target site. The optimal oligonucleotide increases protein levels in transfected fibroblasts by a mean value of 2.6-fold (maximum 4.6-fold), and after two rounds of transfection the effect lasted for a month. Oligonucleotides targeted to the upstream exon (exon 6 in SMN) are also effective. We conclude that TOES oligonucleotides are highly effective reagents for restoring the splicing of refractory exons and can act across long introns.     

The gene for human acidic fibroblast growth factor encodes two upstream exons alternatively spliced to the first coding exon

We have isolated two cDNA clones encoding human acidic fibroblast growth factor (aFGF) which represent the utilization of alternative upstream exons in aFGF mRNA. Isolation and sequence analysis of genomic clones spanning the first coding exon and each of the upstream sequences confirms that the divergent 5' sequences are separate exons, spliced alternatively to the first coding exon 34 nucleotides upstream of the initiator AUG codon. Restriction mapping of the genomic clones provides a minimum size estimate of 45 kilobase pairs for the aFGF locus.     

Genome-Wide Comparative Analysis of Microsatellites in Pineapple

Pineapple (Ananas comosus (L.) Merrill) is the second most important tropical fruit in term of international trade. The availability of whole genomic sequences and expressed sequence tags (ESTs) offers an opportunity to identify and characterize microsatellite or simple sequence repeat (SSR) markers in pineapple. A total of 278,245 SSRs and 41,962 SSRs with an overall density of 728.57 SSRs/Mb and 619.37 SSRs/Mb were mined from genomic and ESTs sequences, respectively. 5′-untranslated regions (5′-UTRs) had the greatest amount of SSRs, 3.6–5.2 fold higher SSR density than other regions. For repeat length, 12 bp was the predominant repeat length in both assembled genome and ESTs. Class I SSRs were underrepresented compared with class II SSRs. For motif length, dinucleotide repeats were the most abundant in genomic sequences, whereas trinucleotides were the most common motif in ESTs. Tri- and hexanucleotides of total SSRs were more prevalent in ESTs than in the whole genome. The SSR frequency decreased dramatically as repeat times increased. AT was the most frequent single motif across the entire genome while AG was the most abundant motif in ESTs. Across six examined plant species, the pineapple genome displayed the highest density, substantially more than the second-place cucumber. Annotation and expression analyses were also conducted for genes containing SSRs. This thorough analysis of SSR markers in pineapple provided valuable information on the frequency and distribution of SSRs in the pineapple genome. This genomic resource will expedite genomic research and pineapple improvement.     

Massively parallel resequencing of the isogenic Drosophila melanogaster strain w1118; iso-2; iso-3 identifies hotspots for mutations in sensory perception genes

We used the Illumina reversible-short sequencing technology to obtain 17-fold average depth (s.d.~8) of ~94% of the euchromatic genome and ~1-5% of the heterochromatin sequence of the Drosophila melogaster isogenic strain w1118; iso-2; iso-3. We show that this strain has a ~9 kb deletion that uncovers the first exon of the white (w) gene, ~4 kb of downstream promoter sequences, and most of the first intron, thus demonstrating that whole-genome sequencing can be used for mutation characterization. We chose this strain because there are thousands of transposon insertion lines and hundreds of isogenic deficiency lines available with this genetic background, such as the Exelixis, Inc., and the DrosDEL collections. We compared our sequence to Release 5 of the finished reference genome sequence which was made from the isogenic strain y1; cn1 bw1 sp1 and identified ~356,614 candidate SNPs in the ~117 Mb unique sequence genome, which represents a substitution rate of ~1/305 nucleotides (~0.30%). The distribution of SNPs is not uniform, but rather there is a ~2-fold increase in SNPs on the autosome arms compared with the X chromosome and a ~7-fold increase when compared to the small 4th chromosome. This is consistent with previous analyses that demonstrated a correlation between recombination frequency and SNP frequency. An unexpected finding was a SNP hotpot in a ~20Mb central region of the 4th chromosome, which might indicate higher than expected recombination frequency in this region of this chromosome. Interestingly, genes involved in sensory perception are enriched in SNP hotspots and genes encoding developmental genes are enriched in SNP coldspots, which suggests that recombination frequencies might be proportional to the evolutionary selection coefficient. There are currently 12 Drosophila species sequenced, and this represents one of many isogenic Drosophila melanogaster genome sequences that are in progress. Because of the dramatic increase in power in using isogenic lines rather than outbred individuals, the SNP information should be valuable as a test bed for understanding genotype-by-environment interactions in human population studies.