Spliced-leader RNA trans splicing in a chordate, Oikopleura dioica, with a compact genome

Extreme genome scrambling in marine planktonic Oikopleura dioica cryptic species

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicate Oikopleura dioica around the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes, Drosophila, or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation in O. dioica suggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures.

Ribosomal Protein Cluster Organization in Asgard Archaea

It has been proposed that the superphylum of Asgard Archaea may represent a historical link between the Archaea and Eukarya. Following the discovery of the Archaea, it was soon appreciated that archaeal ribosomes were more similar to those of Eukarya rather than Bacteria. Coupled with other eukaryotic-like features, it has been suggested that the Asgard Archaea may be directly linked to eukaryotes. However, the genomes of Bacteria and non-Asgard Archaea generally organize ribosome-related genes into clusters that likely function as operons. In contrast, eukaryotes typically do not employ an operon strategy. To gain further insight into conservation of the r-protein genes, the genome order of conserved ribosomal protein (r-protein) coding genes was identified in 17 Asgard genomes (thirteen complete genomes and four genomes with less than 20 contigs) and compared with those found previously in non-Asgard archaeal and bacterial genomes. A universal core of two clusters of 14 and 4 cooccurring r-proteins, respectively, was identified in both the Asgard and non-Asgard Archaea. The equivalent genes in the E. coli version of the cluster are found in the S10 and spc operons. The large cluster of 14 r-protein genes (uS19-uL22-uS3-uL29-uS17 from the S10 operon and uL14-uL24-uL5-uS14-uS8-uL6-uL18-uS5-uL30-uL15 from the spc operon) occurs as a complete set in the genomes of thirteen Asgard genomes (five Lokiarchaeotes, three Heimdallarchaeotes, one Odinarchaeote, and four Thorarchaeotes). Four less conserved clusters with partial bacterial equivalents were found in the Asgard. These were the L30e (str operon in Bacteria) cluster, the L18e (alpha operon in Bacteria) cluster, the S24e-S27ae-rpoE1 cluster, and the L31e, L12..L1 cluster. Finally, a new cluster referred to as L7ae was identified. In many cases, r-protein gene clusters/operons are less conserved in their organization in the Asgard group than in other Archaea. If this is generally true for nonribosomal gene clusters, the results may have implications for the history of genome organization. In particular, there may have been an early transition to or from the operon approach to genome organization. Other nonribosomal cellular features may support different relationships. For this reason, it may be important to consider ribosome features separately.

Neighbours matter: Effects of genomic organization on gene expression plasticity in response to environmental stresses during biological invasions

Gene expression regulation has been widely recognized as an important molecular mechanism underlying phenotypic plasticity in environmental adaptation. However, it remains largely unexplored on the effects of genomic organization on gene expression plasticity under environmental stresses during biological invasions. Here, we use an invasive model ascidian, Ciona robusta, to investigate how genomic organization affects gene expression in response to salinity stresses during range expansions. Our study showed that neighboring genes were co-expressed and approximately 30% of stress responsive genes were physically clustered on chromosomes. Such coordinated expression was substantially affected by the physical distance and orientation of genes. Interestingly, the overall expression correlation of neighboring genes was significantly decreased under high salinity stresses, illustrating that the co-expression regulation could be disrupted by salinity challenges. Furthermore, the clustering of genes was associated with their function constraints and expression patterns - operon genes enriched in gene expression machinery had the highest transcriptional activity and expression stability. Notably, our analyses showed that the tail-to-tail genes, mainly involved in biological functions related to phosphorylation, homeostatic process, and ion transport, exhibited higher intrinsic expression variability and greater response to salinity challenges. Altogether, the results obtained here provide new insights into the effects of gene organization on gene expression plasticity under environmental challenges, hence improving our knowledge on mechanisms of rapid environmental adaptation during biological invasions.

Developmental biology of the larvacean Oikopleura dioica: Genome resources, functional screening, and imaging

The larvacean Oikopleura dioica is a cosmopolitan planktonic chordate and is closely related to vertebrates. It is characterized by a tadpole-shaped morphology with notochord flanked by muscle in the tail and brain on the dorsal side, a short life cycle of five days, a compact genome of approximately 56 Mb, a simple and transparent body with a small number of cells (~4000 in functional juveniles), invariant embryonic cell lineages, and fast development that ensures complete morphogenesis and organ formation 10 h after fertilization. With these features, this marine chordate is a promising and advantageous animal model in which genetic manipulation is feasible. In this review, we introduce relevant resources and modern techniques that have been developed: (1) Genome and transcriptomes. Oikopleura dioica has the smallest genome among non-parasitic metazoans. Its genome databases have been generated using three geographically distant O. dioica populations, and several intra-species sequence differences are becoming evident; (2) Functional genetic knockdown techniques. Comprehensive screening of genes is feasible using ovarian microinjection and double-strand DNA-induced gene knockdown; and (3) Live imaging of embryos and larvae. Application of these techniques has uncovered novel aspects of development, including meiotic cell arrest, left-right patterning, epidermal cell patterning, and mouth formation involving the connection of ectoderm and endoderm sheets. Oikopleura dioca has become very useful for developmental and evolutionary studies in chordates.

Computational analysis of spliced leader trans-splicing in the regenerative flatworm Macrostomum lignano reveals its prevalence in conserved and stem cell related genes

In eukaryotes, trans-splicing is a process of nuclear pre-mRNA maturation where two different RNA molecules are joined together by the spliceosomal machinery utilizing mechanisms similar to cis-splicing. In diverse taxa of lower eukaryotes, spliced leader (SL) trans-splicing is the most frequent type of trans-splicing, when the same sequence derived from short small nuclear RNA molecules, called SL RNAs, is attached to the 5’ ends of different non-processed pre-mRNAs. One of the functions of SL trans-splicing is processing polycistronic pre-mRNA molecules transcribed from operons, when several genes are transcribed as one pre-mRNA molecule. However, only a fraction of trans-spliced genes reside in operons, suggesting that SL trans-splicing must also have some other, less understood functions. Regenerative flatworms are informative model organisms which hold the keys to understand the mechanism of stem cell regulation and specialization during regeneration and homeostasis. Their ability to regenerate is fueled by the division and differentiation of the adult somatic stem cell population called neoblasts. Macrostomum lignano is a flatworm model organism where substantial technological advances have been achieved in recent years, including the development of transgenesis. Although a large fraction of genes in M. lignano were estimated to be SL trans-spliced, SL trans-splicing was not studied in detail in M. lignano before. Here, we performed the first comprehensive study of SL trans-splicing in M. lignano. By reanalyzing the existing genome and transcriptome data of M. lignano, we estimate that 30 % of its genes are SL trans-spliced, 15 % are organized in operons, and almost 40 % are both SL trans-spliced and in operons. We annotated and characterized the sequence of SL RNA and characterized conserved cis- and SL transsplicing motifs. Finally, we found that a majority of SL trans-spliced genes are evolutionarily conserved and significantly over-represented in neoblast-specific genes. Our findings suggest an important role of SL trans-splicing in the regulation and maintenance of neoblasts in M. lignano.

Monodopsis and Vischeria Genomes Shed New Light on the Biology of Eustigmatophyte Algae

Members of eustigmatophyte algae, especially Nannochloropsis and Microchloropsis, have been tapped for biofuel production owing to their exceptionally high lipid content. Although extensive genomic, transcriptomic, and synthetic biology toolkits have been made available for Nannochloropsis and Microchloropsis, very little is known about other eustigmatophytes. Here we present three near-chromosomal and gapless genome assemblies of Monodopsis strains C73 and C141 (60 Mb) and Vischeria strain C74 (106 Mb), which are the sister groups to Nannochloropsis and Microchloropsis in the order Eustigmatales. These genomes contain unusually high percentages of simple repeats, ranging from 12% to 21% of the total assembly size. Unlike Nannochloropsis and Microchloropsis, long interspersed nuclear element repeats are abundant in Monodopsis and Vischeria and might constitute the centromeric regions. We found that both mevalonate and nonmevalonate pathways for terpenoid biosynthesis are present in Monodopsis and Vischeria, which is different from Nannochloropsis and Microchloropsis that have only the latter. Our analysis further revealed extensive spliced leader trans-splicing in Monodopsis and Vischeria at 36-61% of genes. Altogether, the high-quality genomes of Monodopsis and Vischeria not only serve as the much-needed outgroups to advance Nannochloropsis and Microchloropsis research, but also shed new light on the biology and evolution of eustigmatophyte algae.

Screening a Spliced Leader-Based Symbiodinium microadriaticum cDNA Library Using the Yeast-Two Hybrid System Reveals a Hemerythrin-Like Protein as a Putative SmicRACK1 Ligand

The dinoflagellate Symbiodiniaceae family plays a central role in the health of the coral reef ecosystem via the symbiosis that establishes with its inhabiting cnidarians and supports the host metabolism. In the last few decades, coral reefs have been threatened by pollution and rising temperatures which have led to coral loss. These events have raised interest in studying Symbiodiniaceae and their hosts; however, progress in understanding their metabolism, signal transduction pathways, and physiology in general, has been slow because dinoflagellates present peculiar characteristics. We took advantage of one of these peculiarities; namely, the post-transcriptional addition of a Dino Spliced Leader (Dino-SL) to the 5' end of the nuclear mRNAs, and used it to generate cDNA libraries from Symbiodinium microadriaticum. We compared sequences from two Yeast-Two Hybrid System cDNA Libraries, one based on the Dino-SL sequence, and the other based on the SMART technology (Switching Mechanism at 5' end of RNA Transcript) which exploits the template switching function of the reverse transcriptase. Upon comparison of the performance of both libraries, we obtained a significantly higher yield, number and length of sequences, number of transcripts, and better 5' representation from the Dino-SL based library than from the SMART library. In addition, we confirmed that the cDNAs from the Dino-SL library were adequately expressed in the yeast cells used for the Yeast-Two Hybrid System which resulted in successful screening for putative SmicRACK1 ligands, which yielded a putative hemerythrin-like protein.

Telomere-to-telomere assembly of the genome of an individual Oikopleura dioica from Okinawa using Nanopore-based sequencing

Background

The larvacean Oikopleura dioica is an abundant tunicate plankton with the smallest (65–70 Mbp) non-parasitic, non-extremophile animal genome identified to date. Currently, there are two genomes available for the Bergen (OdB3) and Osaka (OSKA2016) O. dioica laboratory strains. Both assemblies have full genome coverage and high sequence accuracy. However, a chromosome-scale assembly has not yet been achieved.

Results

Here, we present a chromosome-scale genome assembly (OKI2018_I69) of the Okinawan O. dioica produced using long-read Nanopore and short-read Illumina sequencing data from a single male, combined with Hi-C chromosomal conformation capture data for scaffolding. The OKI2018_I69 assembly has a total length of 64.3 Mbp distributed among 19 scaffolds. 99% of the assembly is contained within five megabase-scale scaffolds. We found telomeres on both ends of the two largest scaffolds, which represent assemblies of two fully contiguous autosomal chromosomes. Each of the other three large scaffolds have telomeres at one end only and we propose that they correspond to sex chromosomes split into a pseudo-autosomal region and X-specific or Y-specific regions. Indeed, these five scaffolds mostly correspond to equivalent linkage groups in OdB3, suggesting overall agreement in chromosomal organization between the two populations. At a more detailed level, the OKI2018_I69 assembly possesses similar genomic features in gene content and repetitive elements reported for OdB3. The Hi-C map suggests few reciprocal interactions between chromosome arms. At the sequence level, multiple genomic features such as GC content and repetitive elements are distributed differently along the short and long arms of the same chromosome.

Conclusions

We show that a hybrid approach of integrating multiple sequencing technologies with chromosome conformation information results in an accurate de novo chromosome-scale assembly of O. dioica’s highly polymorphic genome. This genome assembly opens up the possibility of cross-genome comparison between O. dioica populations, as well as of studies of chromosomal evolution in this lineage.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07512-6.

SLIDR and SLOPPR: flexible identification of spliced leader trans-splicing and prediction of eukaryotic operons from RNA-Seq data

BACKGROUND

Spliced leader (SL) trans-splicing replaces the 5' end of pre-mRNAs with the spliced leader, an exon derived from a specialised non-coding RNA originating from elsewhere in the genome. This process is essential for resolving polycistronic pre-mRNAs produced by eukaryotic operons into monocistronic transcripts. SL trans-splicing and operons may have independently evolved multiple times throughout Eukarya, yet our understanding of these phenomena is limited to only a few well-characterised organisms, most notably C. elegans and trypanosomes. The primary barrier to systematic discovery and characterisation of SL trans-splicing and operons is the lack of computational tools for exploiting the surge of transcriptomic and genomic resources for a wide range of eukaryotes.

RESULTS

Here we present two novel pipelines that automate the discovery of SLs and the prediction of operons in eukaryotic genomes from RNA-Seq data. SLIDR assembles putative SLs from 5' read tails present after read alignment to a reference genome or transcriptome, which are then verified by interrogating corresponding SL RNA genes for sequence motifs expected in bona fide SL RNA molecules. SLOPPR identifies RNA-Seq reads that contain a given 5' SL sequence, quantifies genome-wide SL trans-splicing events and predicts operons via distinct patterns of SL trans-splicing events across adjacent genes. We tested both pipelines with organisms known to carry out SL trans-splicing and organise their genes into operons, and demonstrate that (1) SLIDR correctly detects expected SLs and often discovers novel SL variants; (2) SLOPPR correctly identifies functionally specialised SLs, correctly predicts known operons and detects plausible novel operons.

CONCLUSIONS

SLIDR and SLOPPR are flexible tools that will accelerate research into the evolutionary dynamics of SL trans-splicing and operons throughout Eukarya and improve gene discovery and annotation for a wide range of eukaryotic genomes. Both pipelines are implemented in Bash and R and are built upon readily available software commonly installed on most bioinformatics servers. Biological insight can be gleaned even from sparse, low-coverage datasets, implying that an untapped wealth of information can be retrieved from existing RNA-Seq datasets as well as from novel full-isoform sequencing protocols as they become more widely available.

Conservative route to genome compaction in a miniature annelid

The causes and consequences of genome reduction in animals are unclear because our understanding of this process mostly relies on lineages with often exceptionally high rates of evolution. Here, we decode the compact 73.8-megabase genome of Dimorphilus gyrociliatus, a meiobenthic segmented worm. The D. gyrociliatus genome retains traits classically associated with larger and slower-evolving genomes, such as an ordered, intact Hox cluster, a generally conserved developmental toolkit and traces of ancestral bilaterian linkage. Unlike some other animals with small genomes, the analysis of the D. gyrociliatus epigenome revealed canonical features of genome regulation, excluding the presence of operons and trans-splicing. Instead, the gene-dense D. gyrociliatus genome presents a divergent Myc pathway, a key physiological regulator of growth, proliferation and genome stability in animals. Altogether, our results uncover a conservative route to genome compaction in annelids, reminiscent of that observed in the vertebrate Takifugu rubripes.

This study reports the genome of the miniature segmented annelid Dimorphilus gyrociliatus and reveals no drastic changes in genome architecture and regulation, unlike other cases of genome miniaturization.

Phylogenetic Analyses of Glycosyl Hydrolase Family 6 Genes in Tunicates: Possible Horizontal Transfer

Horizontal gene transfer (HGT) is the movement of genetic material between different species. Although HGT is less frequent in eukaryotes than in bacteria, several instances of HGT have apparently shaped animal evolution. One well-known example is the tunicate cellulose synthase gene, CesA, in which a gene, probably transferred from bacteria, greatly impacted tunicate evolution. A Glycosyl Hydrolase Family 6 (GH6) hydrolase-like domain exists at the C-terminus of tunicate CesA, but not in cellulose synthases of other organisms. The recent discovery of another GH6 hydrolase-like gene (GH6-1) in tunicate genomes further raises the question of how tunicates acquired GH6. To examine the probable origin of these genes, we analyzed the phylogenetic relationship of GH6 proteins in tunicates and other organisms. Our analyses show that tunicate GH6s, the GH6-1 gene, and the GH6 part of the CesA gene, form two independent, monophyletic gene groups. We also compared their sequence signatures and exon splice sites. All tunicate species examined have shared splice sites in GH6-containing genes, implying ancient intron acquisitions. It is likely that the tunicate CesA and GH6-1 genes existed in the common ancestor of all extant tunicates.

SLFinder, a pipeline for the novel identification of splice-leader sequences: a good enough solution for a complex problem

BACKGROUND

Spliced Leader trans-splicing is an important mechanism for the maturation of mRNAs in several lineages of eukaryotes, including several groups of parasites of great medical and economic importance. Nevertheless, its study across the tree of life is severely hindered by the problem of identifying the SL sequences that are being trans-spliced.

RESULTS

In this paper we present SLFinder, a four-step pipeline meant to identify de novo candidate SL sequences making very few assumptions regarding the SL sequence properties. The pipeline takes transcriptomic de novo assemblies and a reference genome as input and allows the user intervention on several points to account for unexpected features of the dataset. The strategy and its implementation were tested on real RNAseq data from species with and without SL Trans-Splicing.

CONCLUSIONS

SLFinder is capable to identify SL candidates with good precision in a reasonable amount of time. It is especially suitable for species with unknown SL sequences, generating candidate sequences for further refining and experimental validation.

Trans-splicing of mRNAs links gene transcription to translational control regulated by mTOR

BACKGROUND

In phylogenetically diverse organisms, the 5' ends of a subset of mRNAs are trans-spliced with a spliced leader (SL) RNA. The functions of SL trans-splicing, however, remain largely enigmatic.

RESULTS

We quantified translation genome-wide in the marine chordate, Oikopleura dioica, under inhibition of mTOR, a central growth regulator. Translation of trans-spliced TOP mRNAs was suppressed, consistent with a role of the SL sequence in nutrient-dependent translational control of growth-related mRNAs. Under crowded, nutrient-limiting conditions, O. dioica continued to filter-feed, but arrested growth until favorable conditions returned. Upon release from unfavorable conditions, initial recovery was independent of nutrient-responsive, trans-spliced genes, suggesting animal density sensing as a first trigger for resumption of development.

CONCLUSION

Our results are consistent with a proposed role of trans-splicing in the coordinated translational down-regulation of nutrient-responsive genes under growth-limiting conditions.

Combinatorial chromatin dynamics foster accurate cardiopharyngeal fate choices

During embryogenesis, chromatin accessibility profiles control lineage-specific gene expression by modulating transcription, thus impacting multipotent progenitor states and subsequent fate choices. Subsets of cardiac and pharyngeal/head muscles share a common origin in the cardiopharyngeal mesoderm, but the chromatin landscapes that govern multipotent progenitors competence and early fate choices remain largely elusive. Here, we leveraged the simplicity of the chordate model Ciona to profile chromatin accessibility through stereotyped transitions from naive Mesp+ mesoderm to distinct fate-restricted heart and pharyngeal muscle precursors. An FGF-Foxf pathway acts in multipotent progenitors to establish cardiopharyngeal-specific patterns of accessibility, which govern later heart vs. pharyngeal muscle-specific expression profiles, demonstrating extensive spatiotemporal decoupling between early cardiopharyngeal enhancer accessibility and late cell-type-specific activity. We found that multiple cis-regulatory elements, with distinct chromatin accessibility profiles and motif compositions, are required to activate Ebf and Tbx1/10, two key determinants of cardiopharyngeal fate choices. We propose that these 'combined enhancers' foster spatially and temporally accurate fate choices, by increasing the repertoire of regulatory inputs that control gene expression, through either accessibility and/or activity.

Eukaryotic Acquisition of a Bacterial Operon

Operons are a hallmark of bacterial genomes, where they allow concerted expression of functionally related genes as single polycistronic transcripts. They are rare in eukaryotes, where each gene usually drives expression of its own independent messenger RNAs. Here, we report the horizontal operon transfer of a siderophore biosynthesis pathway from relatives of Escherichia coli into a group of budding yeast taxa. We further show that the co-linearly arranged secondary metabolism genes are expressed, exhibit eukaryotic transcriptional features, and enable the sequestration and uptake of iron. After transfer, several genetic changes occurred during subsequent evolution, including the gain of new transcription start sites that were sometimes within protein-coding sequences, acquisition of polyadenylation sites, structural rearrangements, and integration of eukaryotic genes into the cluster. We conclude that the genes were likely acquired as a unit, modified for eukaryotic gene expression, and maintained by selection to adapt to the highly competitive, iron-limited environment.

Heterodera glycines utilizes promiscuous spliced leaders and demonstrates a unique preference for a species-specific spliced leader over C. elegans SL1

Spliced leader trans-splicing (SLTS) plays a part in the maturation of pre-mRNAs in select species across multiple phyla but is particularly prevalent in Nematoda. The role of spliced leaders (SL) within the cell is unclear and an accurate assessment of SL occurrence within an organism is possible only after extensive sequencing data are available, which is not currently the case for many nematode species. SL discovery is further complicated by an absence of SL sequences from high-throughput sequencing results due to incomplete sequencing of the 5'-ends of transcripts during RNA-seq library preparation, known as 5'-bias. Existing datasets and novel methodology were used to identify both conserved SLs and unique hypervariable SLs within Heterodera glycines, the soybean cyst nematode. In H. glycines, twenty-one distinct SL sequences were found on 2,532 unique H. glycines transcripts. The SL sequences identified on the H. glycines transcripts demonstrated a high level of promiscuity, meaning that some transcripts produced as many as nine different individual SL-transcript combinations. Most uniquely, transcriptome analysis revealed that H. glycines is the first nematode to demonstrate a higher SL trans-splicing rate using a species-specific SL over well-conserved Caenorhabditis elegans SL-like sequences.

Oikopleura dioica: An Emergent Chordate Model to Study the Impact of Gene Loss on the Evolution of the Mechanisms of Development

The urochordate Oikopleura dioica is emerging as a nonclassical animal model in the field of evolutionary developmental biology (a.k.a. evo-devo) especially attractive for investigating the impact of gene loss on the evolution of mechanisms of development. This is because this organism fulfills the requirements of an animal model (i.e., has a simple and accessible morphology, a short generation time and life span, and affordable culture in the laboratory and amenable experimental manipulation), but also because O. dioica occupies a key phylogenetic position to understand the diversification and origin of our own phylum, the chordates. During its evolution, O. dioica genome has suffered a drastic process of compaction, becoming the smallest known chordate genome, a process that has been accompanied by exacerbating amount of gene losses. Interestingly, however, despite the extensive gene losses, including entire regulatory pathways essential for the embryonic development of other chordates, O. dioica retains the typical chordate body plan. This unexpected situation led to the formulation of the so-called inverse paradox of evo-devo, that is, when a genetic diversity is able to maintain a phenotypic unity. This chapter reviews the biological features of O. dioica as a model animal, along with the current data on the evolution of its genes and genome. We pay special attention to the numerous examples of gene losses that have taken place during the evolution of this unique animal model, which is helping us to understand to which the limits of evo-devo can be pushed off.

Characterization of spliced leader trans-splicing in a photosynthetic rhizarian amoeba, Paulinella micropora, and its possible role in functional gene transfer

Paulinella micropora is a rhizarian thecate amoeba, belonging to a photosynthetic Paulinella species group that has a unique organelle termed chromatophore, whose cyanobacterial origin is distinct from that of plant and algal chloroplasts. Because acquisition of the chromatophore was quite a recent event compared with that of the chloroplast ancestor, the Paulinella species are thought to be model organisms for studying the early process of primary endosymbiosis. To obtain insight into how endosymbiotically transferred genes acquire expression competence in the host nucleus, here we analyzed the 5′ end sequences of the mRNAs of P. micropora MYN1 strain with the aid of a cap-trapper cDNA library. As a result, we found that mRNAs of 27 genes, including endosymbiotically transferred genes, possessed the common 5′ end sequence of 28–33 bases that were posttranscriptionally added by spliced leader (SL) trans-splicing. We also found two subtypes of SL RNA genes encoded by the P. micropora MYN1 genome. Differing from the other SL trans-splicing organisms that usually possess poly(A)-less SL RNAs, this amoeba has polyadenylated SL RNAs. In this study, we characterize the SL trans-splicing of this unique organism and discuss the putative merits of SL trans-splicing in functional gene transfer and genome evolution.

Distinct core promoter codes drive transcription initiation at key developmental transitions in a marine chordate

Background

Development is largely driven by transitions between transcriptional programs. The initiation of transcription at appropriate sites in the genome is a key component of this and yet few rules governing selection are known. Here, we used cap analysis of gene expression (CAGE) to generate bp-resolution maps of transcription start sites (TSSs) across the genome of Oikopleura dioica, a member of the closest living relatives to vertebrates.

Results

Our TSS maps revealed promoter features in common with vertebrates, as well as striking differences, and uncovered key roles for core promoter elements in the regulation of development. During spermatogenesis there is a genome-wide shift in mode of transcription initiation characterized by a novel core promoter element. This element was associated with > 70% of male-specific transcription, including the use of cryptic internal promoters within operons. In many cases this led to the exclusion of trans-splice sites, revealing a novel mechanism for regulating which mRNAs receive the spliced leader. Binding of the cell cycle regulator, E2F1, is enriched at the TSS of maternal genes in endocycling nurse nuclei. In addition, maternal promoters lack the TATA-like element found in zebrafish and have broad, rather than sharp, architectures with ordered nucleosomes. Promoters of ribosomal protein genes lack the highly conserved TCT initiator. We also report an association between DNA methylation on transcribed gene bodies and the TATA-box.

Conclusions

Our results reveal that distinct functional promoter classes and overlapping promoter codes are present in protochordates like in vertebrates, but show extraordinary lineage-specific innovations. Furthermore, we uncover a genome-wide, developmental stage-specific shift in the mode of TSS selection. Our results provide a rich resource for the study of promoter structure and evolution in Metazoa.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4504-5) contains supplementary material, which is available to authorized users.

Nonprotein-Coding RNAs as Regulators of Development in Tunicates

Tunicates, or urochordates, are a group of small marine organisms that are found widely throughout the seas of the world. As most plausible sister group of the vertebrates, they are of utmost importance for a comprehensive understanding of chordate evolution; hence, they have served as model organisms for many aspects of the developmental biology. Current genomic analysis of tunicates indicates that their genomes evolved with a fast rate not only at the level of nucleotide substitutions but also in terms of genomic organization. The latter involves genome reduction, rearrangements, as well as the loss of some important coding and noncoding RNA (ncRNAs) elements and even entire genomic regions that are otherwise well conserved. These observations are largely based on evidence from comparative genomics resulting from the analysis of well-studied gene families such as the Hox genes and their noncoding elements. In this chapter, the focus lies on the ncRNA complement of tunicates, with a particular emphasis on microRNAs, which have already been studied extensively for other animal clades. MicroRNAs are known as important regulators of key genes in animal development, and they are intimately related to the increase morphological complexity in higher metazoans. Here we review the discovery, evolution, and genome organization of the miRNA repertoire, which has been drastically reduced and restructured in tunicates compared to the chordate ancestor. Known functions of microRNAs as regulators of development in tunicates are a central topic. For instance, we consider the role of miRNAs as regulators of the muscle development and their importance in the regulation of the differential expression during the oral siphon regeneration. Beyond microRNAs, we touch upon the functions of some other ncRNAs such as yellow crescent RNA, moRNAs, RMST lncRNAs, or spliced-leader (SL) RNAs, which have diverse functions associated with the embryonic development, neurogenesis, and mediation of mRNA stability in general.

Genomic and Transcriptomic Analysis Reveals Spliced Leader Trans-Splicing in Cryptomonads

Spliced leader trans-splicing (SLTS) is a poorly understood mechanism that is found in a diversity of eukaryotic lineages. In SLTS, a short RNA sequence is added near the 5′ ends of the transcripts of protein-coding genes by a modified spliceosomal reaction. Available data suggest that SLTS has evolved many times, and might be more likely to evolve in animals. That SLTS might be more likely to evolve in the context of the generally complex transcriptomes characteristic of animals suggests the possibility that SLTS functions in gene regulation or transcriptome diversification, however no general novel function for SLTS is known. Here, I report SLTS in a lineage of cellularly complex unicellular eukaryotes. Cryptomonads are a group of eukaryotic algae that acquired photosynthetic capacity by secondary endosymbiosis of a red alga, and that retain a reduced copy of the nucleus of the engulfed alga. I estimate that at least one-fifth of genes in the model cryptomonad Guillardia theta and its relative Hanusia phi undergo SLTS. I show that hundreds of genes in G. theta generate alternative transcripts by SLTS at alternative sites, however I find little evidence for alternative protein production by alternative SLTS splicing. Interestingly, I find no evidence for substantial operon structure in the G. theta genome, in contrast to previous findings in other lineages with SLTS. These results extend SLTS to another major group of eukaryotes, and heighten the mystery of the evolution of SLTS and its association with cellular and transcriptomic complexity.

Sex-specific chromatin landscapes in an ultra-compact chordate genome

Background

In multicellular organisms, epigenome dynamics are associated with transitions in the cell cycle, development, germline specification, gametogenesis and inheritance. Evolutionarily, regulatory space has increased in complex metazoans to accommodate these functions. In tunicates, the sister lineage to vertebrates, we examine epigenome adaptations to strong secondary genome compaction, sex chromosome evolution and cell cycle modes.

Results

Across the 70 MB Oikopleura dioica genome, we profiled 19 histone modifications, and RNA polymerase II, CTCF and p300 occupancies, to define chromatin states within two homogeneous tissues with distinct cell cycle modes: ovarian endocycling nurse nuclei and mitotically proliferating germ nuclei in testes. Nurse nuclei had active chromatin states similar to other metazoan epigenomes, with large domains of operon-associated transcription, a general lack of heterochromatin, and a possible role of Polycomb PRC2 in dosage compensation. Testis chromatin states reflected transcriptional activity linked to spermatogenesis and epigenetic marks that have been associated with establishment of transgenerational inheritance in other organisms. We also uncovered an unusual chromatin state specific to the Y-chromosome, which combined active and heterochromatic histone modifications on specific transposable elements classes, perhaps involved in regulating their activity.

Conclusions

Compacted regulatory space in this tunicate genome is accompanied by reduced heterochromatin and chromatin state domain widths. Enhancers, promoters and protein-coding genes have conserved epigenomic features, with adaptations to the organization of a proportion of genes in operon units. We further identified features specific to sex chromosomes, cell cycle modes, germline identity and dosage compensation, and unusual combinations of histone PTMs with opposing consensus functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0110-4) contains supplementary material, which is available to authorized users.

Co-evolution of SNF spliceosomal proteins with their RNA targets in trans-splicing nematodes

Although the mechanism of pre-mRNA splicing has been well characterized, the evolution of spliceosomal proteins is poorly understood. The U1A/U2B″/SNF family (hereafter referred to as the SNF family) of RNA binding spliceosomal proteins participates in both the U1 and U2 small interacting nuclear ribonucleoproteins (snRNPs). The highly constrained nature of this system has inhibited an analysis of co-evolutionary trends between the proteins and their RNA binding targets. Here we report accelerated sequence evolution in the SNF protein family in Phylum Nematoda, which has allowed an analysis of protein:RNA co-evolution. In a comparison of SNF genes from ecdysozoan species, we found a correlation between trans-splicing species (nematodes) and increased phylogenetic branch lengths of the SNF protein family, with respect to their sister clade Arthropoda. In particular, we found that nematodes (~70-80 % of pre-mRNAs are trans-spliced) have experienced higher rates of SNF sequence evolution than arthropods (predominantly cis-spliced) at both the nucleotide and amino acid levels. Interestingly, this increased evolutionary rate correlates with the reliance on trans-splicing by nematodes, which would alter the role of the SNF family of spliceosomal proteins. We mapped amino acid substitutions to functionally important regions of the SNF protein, specifically to sites that are predicted to disrupt protein:RNA and protein:protein interactions. Finally, we investigated SNF's RNA targets: the U1 and U2 snRNAs. Both are more divergent in nematodes than arthropods, suggesting the RNAs have co-evolved with SNF in order to maintain the necessarily high affinity interaction that has been characterized in other species.

A functional difference between native and horizontally acquired genes in bdelloid rotifers

The form of RNA processing known as SL trans-splicing involves the transfer of a short conserved sequence, the spliced leader (SL), from a noncoding SL RNA to the 5' ends of mRNA molecules. SL trans-splicing occurs in several animal taxa, including bdelloid rotifers (Rotifera, Bdelloidea). One striking feature of these aquatic microinvertebrates is the large proportion of foreign genes, i.e. those acquired by horizontal gene transfer from other organisms, in their genomes. However, whether such foreign genes behave similarly to native genes has not been tested in bdelloids or any other animal. We therefore used a combination of experimental and computational methods to examine whether transcripts of foreign genes in bdelloids were SL trans-spliced, like their native counterparts. We found that many foreign transcripts contain SLs, use similar splice acceptor sequences to native genes, and are able to undergo alternative trans-splicing. However, a significantly lower proportion of foreign mRNAs contains SL sequences than native transcripts. This demonstrates a novel functional difference between foreign and native genes in bdelloids and suggests that SL trans-splicing is not essential for the expression of foreign genes, but is acquired during their domestication.

DNA interference: DNA-induced gene silencing in the appendicularian Oikopleura dioica

RNA interference is widely employed as a gene-silencing system in eukaryotes for host defence against invading nucleic acids. In response to invading double-stranded RNA (dsRNA), mRNA is degraded in sequence-specific manner. So far, however, DNA interference (DNAi) has been reported only in plants, ciliates and archaea, and has not been explored in Metazoa. Here, we demonstrate that linear double-stranded DNA promotes both sequence-specific transcription blocking and mRNA degradation in developing embryos of the appendicularian Oikopleura dioica. Introduced polymerase chain reaction (PCR) products or linearized plasmids encoding Brachyury induced tail malformation and mRNA degradation. This malformation was also promoted by DNA fragments of the putative 5'-flanking region and intron without the coding region. PCR products encoding Zic-like1 and acetylcholine esterase also induced loss of sensory organ and muscle acetylcholinesterase activity, respectively. Co-injection of mRNA encoding EGFP and mCherry, and PCR products encoding these fluorescent proteins, induced sequence-specific decrease in the green or red fluorescence, respectively. These results suggest that O. dioica possesses a defence system against exogenous DNA and RNA, and that DNA fragment-induced gene silencing would be mediated through transcription blocking as well as mRNA degradation. This is the first report of DNAi in Metazoa.

Genome-wide identification and characterization of transcription start sites and promoters in the tunicate Ciona intestinalis

The tunicate Ciona intestinalis, an invertebrate chordate, has recently emerged as a powerful model organism for gene regulation analysis. However, few studies have been conducted to identify and characterize its transcription start sites (TSSs) and promoters at the genome-wide level. Here, using TSS-seq, we identified TSSs at the genome-wide scale and characterized promoters in C. intestinalis. Specifically, we identified TSS clusters (TSCs), high-density regions of TSS-seq tags, each of which appears to originate from an identical promoter. TSCs were found not only at known TSSs but also in other regions, suggesting the existence of many unknown transcription units in the genome. We also identified candidate promoters of 79 ribosomal protein (RP) genes, each of which had the major TSS in a polypyrimidine tract and showed a sharp TSS distribution like human RP gene promoters. Ciona RP gene promoters, however, did not appear to have typical TATA boxes, unlike human RP gene promoters. In Ciona non-RP promoters, two pyrimidine-purine dinucleotides, CA and TA, were frequently used as TSSs. Despite the absence of CpG islands, Ciona TATA-less promoters showed low expression specificity like CpG-associated human TATA-less promoters. By using TSS-seq, we also predicted trans-spliced gene TSSs and found that their downstream regions had higher G+T content than those of non-trans-spliced gene TSSs. Furthermore, we identified many putative alternative promoters, some of which were regulated in a tissue-specific manner. Our results provide valuable information about TSSs and promoter characteristics in C. intestinalis and will be helpful in future analysis of transcriptional regulation in chordates.

Transcriptomic profiling of gene expression and RNA processing during Leishmania major differentiation

Protozoan parasites of the genus Leishmania are the etiological agents of leishmaniasis, a group of diseases with a worldwide incidence of 0.9–1.6 million cases per year. We used RNA-seq to conduct a high-resolution transcriptomic analysis of the global changes in gene expression and RNA processing events that occur as L. major transforms from non-infective procyclic promastigotes to infective metacyclic promastigotes. Careful statistical analysis across multiple biological replicates and the removal of batch effects provided a high quality framework for comprehensively analyzing differential gene expression and transcriptome remodeling in this pathogen as it acquires its infectivity. We also identified precise 5′ and 3′ UTR boundaries for a majority of Leishmania genes and detected widespread alternative trans-splicing and polyadenylation. An investigation of possible correlations between stage-specific preferential trans-splicing or polyadenylation sites and differentially expressed genes revealed a lack of systematic association, establishing that differences in expression levels cannot be attributed to stage-regulated alternative RNA processing. Our findings build on and improve existing expression datasets and provide a substantially more detailed view of L. major biology that will inform the field and potentially provide a stronger basis for drug discovery and vaccine development efforts.

Trans-splicing in metazoans: A link to translational control?

The trans-splicing of a spliced-leader RNA to a subset of mRNAs is a phenomenon that occurs in many species, including Caenorhabditis elegans, and yet the driving force for its evolution in disparate groups of animals remains unclear. Polycistronic mRNA resulting from the transcription of operons is resolved via trans-splicing, but operons comprise only a sub-set of trans-spliced genes. Using the marine chordate, Oikopleura dioica, we recently tested the hypothesis that metazoan operons accelerate recovery from growth arrest. We found no supporting evidence for this in O. dioica. Instead we found a striking relationship between trans-splicing and maternal mRNA in O. dioica, C. elegans and the ascidian, Ciona intestinalis. Furthermore, in O. dioica and C. elegans, we found evidence to suggest a role for mTOR signaling in the translational control of growth-related, trans-spliced maternal mRNAs. We propose that this may be a mechanism for adjusting egg number in response to nutrient levels in these species.

Maternal and zygotic transcriptomes in the appendicularian, Oikopleura dioica: novel protein-encoding genes, intra-species sequence variations, and trans-spliced RNA leader

RNA sequencing analysis was carried out to characterize egg and larval transcriptomes in the appendicularian, Oikopleura dioica, a planktonic chordate, which is characterized by rapid development and short life cycle of 5 days, using a Japanese population of the organism. De novo transcriptome assembly matched with 16,423 proteins corresponding to 95.4% of the protein-encoding genes deposited in the OikoBase, the genome database of the Norwegian population. Nucleotide and amino acid sequence identities between the Japanese and Norwegian O. dioica were estimated to be around 91.0 and 94.8%, respectively. We discovered 175 novel protein-encoding genes: 144 unigenes were common to both the Japanese and Norwegian populations, whereas 31 unigenes were not found in the OikoBase genome reference. Among the total 12,311 unigenes, approximately 63% were detected in egg-stage RNAs, whereas 99% were detected in larval stage RNAs; 3772 genes were up-regulated, and 1336 genes were down-regulated more than four-fold in the larvae. Gene ontology analyses characterized gene activities in these two developmental stages. We found a messenger RNA (mRNA) 5' trans-spliced leader, which was observed in 40.8% of the total unique transcripts. It showed preferential linkage to adenine at the 5' ends of the downstream exons. Trans-splicing was observed more frequently in egg mRNAs compared with larva-specific mRNAs.

Trans-splicing and operons in metazoans: translational control in maternally regulated development and recovery from growth arrest

Polycistronic mRNAs transcribed from operons are resolved via the trans-splicing of a spliced-leader (SL) RNA. Trans-splicing also occurs at monocistronic transcripts. The phlyogenetically sporadic appearance of trans-splicing and operons has made the driving force(s) for their evolution in metazoans unclear. Previous work has proposed that germline expression drives operon organization in Caenorhabditis elegans, and a recent hypothesis proposes that operons provide an evolutionary advantage via the conservation of transcriptional machinery during recovery from growth arrested states. Using a modified cap analysis of gene expression protocol we mapped sites of SL trans-splicing genome-wide in the marine chordate Oikopleura dioica. Tiled microarrays revealed the expression dynamics of trans-spliced genes across development and during recovery from growth arrest. Operons did not facilitate recovery from growth arrest in O. dioica. Instead, we found that trans-spliced transcripts were predominantly maternal. We then analyzed data from C. elegans and Ciona intestinalis and found that an enrichment of trans-splicing and operon gene expression in maternal mRNA is shared between all three species, suggesting that this may be a driving force for operon evolution in metazoans. Furthermore, we found that the majority of known terminal oligopyrimidine (TOP) mRNAs are trans-spliced in O. dioica and that the SL contains a TOP-like motif. This suggests that the SL in O. dioica confers nutrient-dependent translational control to trans-spliced mRNAs via the TOR-signaling pathway. We hypothesize that SL-trans-splicing provides an evolutionary advantage in species that depend on translational control for regulating early embryogenesis, growth and oocyte production in response to nutrient levels.

Operons are a conserved feature of nematode genomes

The organization of genes into operons, clusters of genes that are co-transcribed to produce polycistronic pre-mRNAs, is a trait found in a wide range of eukaryotic groups, including multiple animal phyla. Operons are present in the class Chromadorea, one of the two main nematode classes, but their distribution in the other class, the Enoplea, is not known. We have surveyed the genomes of Trichinella spiralis, Trichuris muris, and Romanomermis culicivorax and identified the first putative operons in members of the Enoplea. Consistent with the mechanism of polycistronic RNA resolution in other nematodes, the mRNAs produced by genes downstream of the first gene in the T. spiralis and T. muris operons are trans-spliced to spliced leader RNAs, and we are able to detect polycistronic RNAs derived from these operons. Importantly, a putative intercistronic region from one of these potential enoplean operons confers polycistronic processing activity when expressed as part of a chimeric operon in Caenorhabditis elegans. We find that T. spiralis genes located in operons have an increased likelihood of having operonic C. elegans homologs. However, operon structure in terms of synteny and gene content is not tightly conserved between the two taxa, consistent with models of operon evolution. We have nevertheless identified putative operons conserved between Enoplea and Chromadorea. Our data suggest that operons and "spliced leader" (SL) trans-splicing predate the radiation of the nematode phylum, an inference which is supported by the phylogenetic profile of proteins known to be involved in nematode SL trans-splicing.

Retroduplication and loss of parental genes is a mechanism for the generation of intronless genes in Ciona intestinalis and Ciona savignyi

Tunicates, the sister clade of vertebrates, have miniature genomes and numerous intronless genes compared to other animals. It is still unclear how the tunicates acquired such a large number of intronless genes. Here, we analyzed sequences and intron-exon organizations of homologous genes from two closely related tunicates, Ciona intestinalis and Ciona savignyi. We found seven cases in which ancestral introns of a gene were completely lost in a species after their divergence. In four cases, both the intronless copy and the intron-containing copy were present in the genome, indicating that the intronless copy was generated by retroduplication. In the other three cases, the intron-containing copy was absent, implying it was lost after retroduplication. This result suggests that retroduplication and loss of parental genes is a major mechanism for the accumulation of intronless genes in tunicates.

Evolutionary genomics of fast evolving tunicates

Tunicates have been extensively studied because of their crucial phylogenetic location (the closest living relatives of vertebrates) and particular developmental plan. Recent genome efforts have disclosed that tunicates are also remarkable in their genome organization and molecular evolutionary patterns. Here, we review these latter aspects, comparing the similarities and specificities of two model species of the group: Oikopleura dioica and Ciona intestinalis. These species exhibit great genome plasticity and Oikopleura in particular has undergone a process of extreme genome reduction and compaction that can be explained in part by gene loss, but is mostly due to other mechanisms such as shortening of intergenic distances and introns, and scarcity of mobile elements. In Ciona, genome reorganization was less severe being more similar to the other chordates in several aspects. Rates and patterns of molecular evolution are also peculiar in tunicates, being Ciona about 50% faster than vertebrates and Oikopleura three times faster. In fact, the latter species is considered as the fastest evolving metazoan recorded so far. Two processes of increase in evolutionary rates have taken place in tunicates. One of them is more extreme, and basically restricted to genes encoding regulatory proteins (transcription regulators, chromatin remodeling proteins, and metabolic regulators), and the other one is less pronounced but affects the whole genome. Very likely adaptive evolution has played a very significant role in the first, whereas the functional and/or evolutionary causes of the second are less clear and the evidence is not conclusive. The evidences supporting the incidence of increased mutation and less efficient negative selection are presented and discussed.

Molecular cloning and characterization of SL3: a stem cell-specific SL RNA from the planarian Schmidtea mediterranea

Spliced leader (SL) trans-splicing is a biological phenomenon, common among many metazoan taxa, consisting in the transfer of a short leader sequence from a small SL RNA to the 5' end of a subset of pre-mRNAs. While knowledge of the biochemical mechanisms driving this process has accumulated over the years, the functional consequences of such post-transcriptional event at the organismal level remain unclear. In addition, the fact that functional analyses have been undertaken mainly in trypanosomes and nematodes leaves a somehow fragmented picture of the possible biological significance and evolution of SL trans-splicing in eukaryotes. Here, we analyzed the spatial expression of SL RNAs in the planarian flatworm Schmidtea mediterranea, with the goal of identifying novel developmental paradigms for the study of trans-splicing in metazoans. Besides the previously identified SL1 and SL2, S. mediterranea expresses a third SL RNA described here as SL3. While, SL1 and SL2 are collectively expressed in a broad range of planarian cell types, SL3 is highly enriched in a subset of the planarian stem cells engaged in regenerative responses. Our findings provide new opportunities to study how trans-splicing may regulate the phenotype of a cell.

To grow or not to grow: nutritional control of development during Caenorhabditis elegans L1 arrest

It is widely appreciated that larvae of the nematode Caenorhabditis elegans arrest development by forming dauer larvae in response to multiple unfavorable environmental conditions. C. elegans larvae can also reversibly arrest development earlier, during the first larval stage (L1), in response to starvation. "L1 arrest" (also known as "L1 diapause") occurs without morphological modification but is accompanied by increased stress resistance. Caloric restriction and periodic fasting can extend adult lifespan, and developmental models are critical to understanding how the animal is buffered from fluctuations in nutrient availability, impacting lifespan. L1 arrest provides an opportunity to study nutritional control of development. Given its relevance to aging, diabetes, obesity and cancer, interest in L1 arrest is increasing, and signaling pathways and gene regulatory mechanisms controlling arrest and recovery have been characterized. Insulin-like signaling is a critical regulator, and it is modified by and acts through microRNAs. DAF-18/PTEN, AMP-activated kinase and fatty acid biosynthesis are also involved. The nervous system, epidermis, and intestine contribute systemically to regulation of arrest, but cell-autonomous signaling likely contributes to regulation in the germline. A relatively small number of genes affecting starvation survival during L1 arrest are known, and many of them also affect adult lifespan, reflecting a common genetic basis ripe for exploration. mRNA expression is well characterized during arrest, recovery, and normal L1 development, providing a metazoan model for nutritional control of gene expression. In particular, post-recruitment regulation of RNA polymerase II is under nutritional control, potentially contributing to a rapid and coordinated response to feeding. The phenomenology of L1 arrest will be reviewed, as well as regulation of developmental arrest and starvation survival by various signaling pathways and gene regulatory mechanisms.

The candidate histocompatibility locus of a Basal chordate encodes two highly polymorphic proteins

The basal chordate Botryllus schlosseri undergoes a natural transplantation reaction governed by a single, highly polymorphic locus called the fuhc. Our initial characterization of this locus suggested it encoded a single gene alternatively spliced into two transcripts: a 555 amino acid–secreted form containing the first half of the gene, and a full-length, 1008 amino acid transmembrane form, with polymorphisms throughout the ectodomain determining outcome. We have now found that the locus encodes two highly polymorphic genes which are separated by a 227 bp intergenic region: first, the secreted form as previously described, and a second gene encoding a 531 amino acid membrane-bound gene containing three extracellular immunoglobulin domains. While northern blotting revealed only these two mRNAs, both PCR and mRNA-seq detect a single capped and polyadenylated transcript that encodes processed forms of both genes linked by the intergenic region, as well as other transcripts in which exons of the two genes are spliced together. These results might suggest that the two genes are expressed as an operon, during which both genes are co-transcribed and then trans-spliced into two separate messages. This type of transcriptional regulation has been described in tunicates previously; however, the membrane-bound gene does not encode a typical Splice Leader (SL) sequence at the 5′ terminus that usually accompanies trans-splicing. Thus, the presence of stable transcripts encoding both genes may suggest a novel mechanism of regulation, or conversely may be rare but stable transcripts in which the two mRNAs are linked due to a small amount of read-through by RNA polymerase. Both genes are highly polymorphic and co-expressed on tissues involved in histocompatibility. In addition, polymorphisms on both genes correlate with outcome, although we have found a case in which it appears that the secreted form may be major allorecognition determinant.

Transcriptomic study reveals widespread spliced leader trans-splicing, short 5'-UTRs and potential complex carbon fixation mechanisms in the euglenoid Alga Eutreptiella sp

Eutreptiella are an evolutionarily unique and ecologically important genus of microalgae, but they are poorly understood with regard to their genomic make-up and expression profiles. Through the analysis of the full-length cDNAs from a Eutreptiella species, we found a conserved 28-nt spliced leader sequence (Eut-SL, ACACUUUCUGAGUGUCUAUUUUUUUUCG) was trans-spliced to the mRNAs of Eutreptiella sp. Using a primer derived from Eut-SL, we constructed four cDNA libraries under contrasting physiological conditions for 454 pyrosequencing. Clustering analysis of the ∼1.9×10(6) original reads (average length 382 bp) yielded 36,643 unique transcripts. Although only 28% of the transcripts matched documented genes, this fraction represents a functionally very diverse gene set, suggesting that SL trans-splicing is likely ubiquitous in this alga's transcriptome. The mRNAs of Eutreptiella sp. seemed to have short 5'- untranslated regions, estimated to be 21 nucleotides on average. Among the diverse biochemical pathways represented in the transcriptome we obtained, carbonic anhydrase and genes known to function in the C4 pathway and heterotrophic carbon fixation were found, posing a question whether Eutreptiella sp. employs multifaceted strategies to acquire and fix carbon efficiently. This first large-scale transcriptomic dataset for a euglenoid uncovers many potential novel genes and overall offers a valuable genetic resource for research on euglenoid algae.

Genome-wide computational identification of bicistronic mRNA in humans

Mammalian bicistronic mRNA is a recently discovered mammalian gene structure. Several reported cases of mammalian bicistronic mRNA indicated that genes of this structure play roles in some important biological processes. However, a genome-wide computational identification of bicistronic mRNA in mammalian genome, such as human genome, is still lacking. Here we used a comparative genomics approach to identify the frequency of human bicistronic mRNA. We then validated the result by using a new support vector machine (SVM) model. We identified 43 human bicistronic mRNAs in 30 distinct genes. Our literature analysis shows that our method recovered 100 % (6/6) of the previously known bicistronic mRNAs which had been experimentally confirmed by other groups. Our graph theory-based analysis and GO analysis indicated that human bicistronic mRNAs are prone to produce different yet closely functionally related proteins. In addition, we also described and analyzed three different mechanisms of ORF fusion. Our method of identifying bicistronic mRNAs in human genome provides a model for the computational identification of characteristic gene structures in mammalian genomes. We anticipate that our data will facilitate further molecular characterization and functional study of human bicistronic mRNA.

Transcriptional enhancers in ascidian development

The study of cis-regulatory DNAs that control developmental gene expression is integral to the modeling of comprehensive genomic regulatory networks for embryogenesis. Ascidian embryos provide a unique opportunity for the analysis of cis-regulatory DNAs with cellular resolution in the context of a simple but typical chordate body plan. Here, we review landmark studies that have laid the foundations for the study of transcriptional enhancers, among other cis-regulatory DNAs, and their roles in ascidian development. The studies using ascidians of the Ciona genus have capitalized on a unique electroporation technique that permits the simultaneous transfection of hundreds of fertilized eggs, which develop rapidly and express transgenes with little mosaicism. Current studies using the ascidian embryo benefit from extensively annotated genomic resources to characterize transcript models in silico. The search for functional noncoding sequences can be guided by bioinformatic analyses combining evolutionary conservation, gene coexpression, and combinations of overrepresented short-sequence motifs. The power of the transient transfection assays has allowed thorough dissection of numerous cis-regulatory modules, which provided insights into the functional constraints that shape enhancer architecture and diversification. Future studies will benefit from pioneering stable transgenic lines and the analysis of chromatin states. Whole genome expression, functional and DNA binding data are being integrated into comprehensive genomic regulatory network models of early ascidian cell specification with a single-cell resolution that is unique among chordate model systems.

TEMPERATURE-DEPENDENT STRUCTURAL VARIABILITY OF RNAs: SPLICED LEADER RNAs AND THEIR EVOLUTIONARY HISTORY

The structures attained by RNA molecules depend not only on their sequence but also on environmental parameters such as their temperature. So far, this effect has been largely neglected in bioinformatics studies. Here, we show that structural comparisons can be facilitated and more coherent structural models can be obtained when differences in environmental parameters are taken into account. We re-evaluate the secondary structures of the spliced leader (SL) RNAs from the seven eukaryotic phyla in which SL RNA trans-splicing has been described. Adjusting structure prediction to the natural growth temperatures and considering energetically similar secondary structures, we observe striking similarities among Euglenida, Kinetoplastida, Dinophyceae, Cnidaria, Rotifera, Nematoda, Platyhelminthes, and Tunicata that cannot be explained easily by the independent innovation of SL RNAs in each of these phyla. Supplementary Table is available at .

Metazoan operons accelerate recovery from growth-arrested states

Existing theories explain why operons are advantageous in prokaryotes, but their occurrence in metazoans is an enigma. Nematode operon genes, typically consisting of growth genes, are significantly upregulated during recovery from growth-arrested states. This expression pattern is anticorrelated to nonoperon genes, consistent with a competition for transcriptional resources. We find that transcriptional resources are initially limiting during recovery and that recovering animals are highly sensitive to any additional decrease in transcriptional resources. We provide evidence that operons become advantageous because, by clustering growth genes into operons, fewer promoters compete for the limited transcriptional machinery, effectively increasing the concentration of transcriptional resources and accelerating recovery. Mathematical modeling reveals how a moderate increase in transcriptional resources can substantially enhance transcription rate and recovery. This design principle occurs in different nematodes and the chordate C. intestinalis. As transition from arrest to rapid growth is shared by many metazoans, operons could have evolved to facilitate these processes.

Evolutionary crossroads in developmental biology: the tunicates

The tunicates, or urochordates, constitute a large group of marine animals whose recent common ancestry with vertebrates is reflected in the tadpole-like larvae of most tunicates. Their diversity and key phylogenetic position are enhanced, from a research viewpoint, by anatomically simple and transparent embryos, compact rapidly evolving genomes, and the availability of powerful experimental and computational tools with which to study these organisms. Tunicates are thus a powerful system for exploring chordate evolution and how extreme variation in genome sequence and gene regulatory network architecture is compatible with the preservation of an ancestral chordate body plan.

Analysis of a botryllid enriched-full-length cDNA library: insight into the evolution of spliced leader trans-splicing in tunicates

In some animals, mRNA may be modified after transcription by the addition of a 5' spliced leader sequence. This is known as spliced leader (SL) trans-splicing, and is of uncertain function and evolutionary origin. Here, we report the identification of SL trans-splicing in the colonial ascidian Botryllus schlosseri. Combining our own expressed sequence tag (EST) data with additional data from GenBank, we identify the dominant spliced leader sequence and show it to be similar to that of other ascidians and to that of Oikopleura dioica, a basally diverging tunicate. Gene Ontology analysis of B. schlosseri ESTs with and without a 5' spliced leader shows that genes encoding ribosomal proteins tend not to be trans-spliced, a character shared with the ascidian Ciona intestinalis. We also examine individual cases of genes that produce mRNAs that are SL trans-spliced in B. schlosseri but not in C. intestinalis. We conclude that SL trans-splicing evolved early in the tunicate lineage and shows stability over considerable evolutionary time. However, SL trans-splicing may be gained or lost in individual genes.