A diversity of U1 small nuclear RNAs in the silk moth Bombyx mori

Developmental analysis of spliceosomal snRNA isoform expression

Pre-mRNA splicing is a critical step in eukaryotic gene expression that contributes to proteomic, cellular, and developmental complexity. Small nuclear (sn)RNAs are core spliceosomal components; however, the extent to which differential expression of snRNA isoforms regulates splicing is completely unknown. This is partly due to difficulties in the accurate analysis of the spatial and temporal expression patterns of snRNAs. Here, we use high-throughput RNA-sequencing (RNA-seq) data to profile expression of four major snRNAs throughout Drosophila development. This analysis shows that individual isoforms of each snRNA have distinct expression patterns in the embryo, larva, and pharate adult stages. Expression of these isoforms is more heterogeneous during embryogenesis; as development progresses, a single isoform from each snRNA subtype gradually dominates expression. Despite the lack of stable snRNA orthologous groups during evolution, this developmental switching of snRNA isoforms also occurs in distantly related vertebrate species, such as Xenopus, mouse, and human. Our results indicate that expression of snRNA isoforms is regulated and lays the foundation for functional studies of individual snRNA isoforms.

U1 snDNA clusters in grasshoppers: chromosomal dynamics and genomic organization

The spliceosome, constituted by a protein set associated with small nuclear RNA (snRNA), is responsible for mRNA maturation through intron removal. Among snRNA genes, U1 is generally a conserved repetitive sequence. To unveil the chromosomal/genomic dynamics of this multigene family in grasshoppers, we mapped U1 genes by fluorescence in situ hybridization in 70 species belonging to the families Proscopiidae, Pyrgomorphidae, Ommexechidae, Romaleidae and Acrididae. Evident clusters were observed in all species, indicating that, at least, some U1 repeats are tandemly arrayed. High conservation was observed in the first four families, with most species carrying a single U1 cluster, frequently located in the third or fourth longest autosome. By contrast, extensive variation was observed among Acrididae, from a single chromosome pair carrying U1 to all chromosome pairs carrying it, with occasional occurrence of two or more clusters in the same chromosome. DNA sequence analysis in Eyprepocnemis plorans (species carrying U1 clusters on seven different chromosome pairs) and Locusta migratoria (carrying U1 in a single chromosome pair) supported the coexistence of functional and pseudogenic lineages. One of these pseudogenic lineages was truncated in the same nucleotide position in both species, suggesting that it was present in a common ancestor to both species. At least in E. plorans, this U1 snDNA pseudogenic lineage was associated with 5S rDNA and short interspersed elements (SINE)-like mobile elements. Given that we conclude in grasshoppers that the U1 snDNA had evolved under the birth-and-death model and that its intragenomic spread might be related with mobile elements.

U1 small nuclear RNA variants differentially form ribonucleoprotein particles in vitro

The U1 small nuclear (sn)RNA participates in splicing of pre-mRNAs by recognizing and binding to 5' splice sites at exon/intron boundaries. U1 snRNAs associate with 5' splice sites in the form of ribonucleoprotein particles (snRNPs) that are comprised of the U1 snRNA and 10 core components, including U1A, U1-70K, U1C and the 'Smith antigen', or Sm, heptamer. The U1 snRNA is highly conserved across a wide range of taxa; however, a number of reports have identified the presence of expressed U1-like snRNAs in multiple species, including humans. While numerous U1-like molecules have been shown to be expressed, it is unclear whether these variant snRNAs have the capacity to form snRNPs and participate in splicing. The purpose of the present study was to further characterize biochemically the ability of previously identified human U1-like variants to form snRNPs and bind to U1 snRNP proteins. A bioinformatics analysis provided support for the existence of multiple expressed variants. In vitro gel shift assays, competition assays, and immunoprecipitations (IPs) revealed that the variants formed high molecular weight assemblies to varying degrees and associated with core U1 snRNP proteins to a lesser extent than the canonical U1 snRNA. Together, these data suggest that the human U1 snRNA variants analyzed here are unable to efficiently bind U1 snRNP proteins. The current work provides additional biochemical insights into the ability of the variants to assemble into snRNPs.

Evolution of spliceosomal snRNA genes in metazoan animals

While studies of the evolutionary histories of protein families are commonplace, little is known on noncoding RNAs beyond microRNAs and some snoRNAs. Here we investigate in detail the evolutionary history of the nine spliceosomal snRNA families (U1, U2, U4, U5, U6, U11, U12, U4atac, and U6atac) across the completely or partially sequenced genomes of metazoan animals. Representatives of the five major spliceosomal snRNAs were found in all genomes. None of the minor splicesomal snRNAs were detected in nematodes or in the shotgun traces of Oikopleura dioica, while in all other animal genomes at most one of them is missing. Although snRNAs are present in multiple copies in most genomes, distinguishable paralogue groups are not stable over long evolutionary times, although they appear independently in several clades. In general, animal snRNA secondary structures are highly conserved, albeit, in particular, U11 and U12 in insects exhibit dramatic variations. An analysis of genomic context of snRNAs reveals that they behave like mobile elements, exhibiting very little syntenic conservation.

U6 snRNA variants isolated from the posterior silk gland of the silk moth Bombyx mori

Five U6 small nuclear RNA (snRNA) isoforms were detected and characterized from the posterior silk gland (PSG) of the silk moth Bombyx mori (Nistari strain). Using the currently accepted U6 secondary structure model as a basis for comparison, the variants were analyzed for nucleotide differences across the sequence with a focus on known functional domains. Differences were observed primarily in single-stranded areas of which sixty percent were found in the highly conserved U4-U6 binding sites. In the Nistari strain, the U6A variant was found to be approximately four times more abundant as part of high molecular weight spliceosomal complexes when compared with U6A in the total unfractionated PSG cell lysate. Additionally, the European 703 B. mori strain total cell lysate U6 snRNA was analyzed and only the dominant U6A isoform initially identified in Nistari was found. Due to U6's essential role in pre-mRNA processing, variants may modulate assemblage of the catalytic core and in doing so potentially affect the rate of splicing. Phylogenetic analysis of the U6 snRNA sequences indicate an ancient divergence of U6 from the self-splicing group II intron module and a high degree of evolutionary conservation across species possibly due to functional constraints on the gene. Using in silico analysis, 35 full-length U6 variants were observed in the recently released Whole Genome Shotgun (WGS) database of the p50T strain. The consensus sequence of these U6 genes from p50T is identical to U6A identified in the Nistari strain. Furthermore p50T variant 1, which is represented in 14 genes, is equivalent to Nistari U6A.

The silk moth Bombyx mori U1 and U2 snRNA variants are differentially expressed

Five U1 and eight U2 isoforms of the silk moth Bombyx mori exhibiting internal nucleotide differences have been previously identified and characterized in various tissues and developmental stages. In this investigation, it is demonstrated that the levels of some snRNA variants differ in egg and silk gland tissue and change during development. Qualitative and quantitative differences in the U1 and U2 variant populations were observed at three developmental points (early, middle and late) of the silk gland (SG) during the fifth instar larval stage of the silk moth. Statistical analyses of the various isoform populations across the fifth instar larval and egg stages show significant differences for some of the U1 and U2 variants. The representation of variant sequences in expressed U1 and U2 sequences (RT-PCR libraries) and in a whole-genome shotgun (WGS) assembly database was confirmed. In addition, conserved elements in the promoter 5'-flanking region of the U1 and U2 variants were identified in the WGS.

Variants of U1 small nuclear RNA assemble into spliceosomal complexes

In this study, the existence of additional U1 snRNA variants in the posterior silk gland of the Bombyx mori Nistari strain from India was investigated. Three new U1 variants were detected. One of the new isoforms (U1 SG1) was found to be preferentially assembled into high molecular weight spliceosomal complexes in comparison with the total cellular lysate RNA control. Structural and nucleotide differences were examined in these new isoforms and compared with the previously reported U1 variants. Free energy (Delta G) values for the entire U1 snRNA secondary structures as well as the individual stem/loops (I, II, III and IV) domains of the isoforms were estimated to determine their structural stability.

Hierarchical structure of donor splice site sequence multitude

An analysis of donor splice site sequence (dsss) multitude was carried out. A special method was developed to study the multitude properties. The multitude of dsss was arranged as a hierarchical tree where the "root" consists of a small group of "ancestors", VAGgtVAG (gt stands in the beginning of intron, V=C,A,G). Every ancestor has its own clan of "descendants" which represent point mutants of it. The latter are split into "progenies". The first progeny looks like point mutants of the ancestors, the second as one-point mutants of sequences of the first progeny, and so forth. This picture looks like a hierarchical phylogenetic tree. The most striking property of the multitude is the successive monotonic decrease in occurrences of the sequences from one progeny to another. It is shown that such results may be modeled by a virtual mutational process.