Combining different mRNA capture methods to analyze the transcriptome: analysis of the Xenopus laevis transcriptome

A program of successive gene expression in mouse one-cell embryos

At the moment of union in fertilization, sperm and oocyte are transcriptionally silent. The ensuing onset of embryonic transcription (embryonic genome activation [EGA]) is critical for development, yet its timing and profile remain elusive in any vertebrate species. We here dissect transcription during EGA by high-resolution single-cell RNA sequencing of precisely synchronized mouse one-cell embryos. This reveals a program of embryonic gene expression (immediate EGA [iEGA]) initiating within 4 h of fertilization. Expression during iEGA produces canonically spliced transcripts, occurs substantially from the maternal genome, and is mostly downregulated at the two-cell stage. Transcribed genes predict regulation by transcription factors (TFs) associated with cancer, including c-Myc. Blocking c-Myc or other predicted regulatory TF activities disrupts iEGA and induces acute developmental arrest. These findings illuminate intracellular mechanisms that regulate the onset of mammalian development and hold promise for the study of cancer.

A cap 0-dependent mRNA capture method to analyze the yeast transcriptome

Analysis of the protein coding transcriptome by the RNA sequencing requires either enrichment of the desired fraction of coding transcripts or depletion of the abundant non-coding fraction consisting mainly of rRNA. We propose an alternative mRNA enrichment strategy based on the RNA-binding properties of the human IFIT1, an antiviral protein recognizing cap 0 RNA. Here, we compare for Saccharomyces cerevisiae an IFIT1-based mRNA pull-down with yeast targeted rRNA depletion by the RiboMinus method. IFIT1-based RNA capture depletes rRNA more effectively, producing high quality RNA-seq data with an excellent coverage of the protein coding transcriptome, while depleting cap-less transcripts such as mitochondrial or some non-coding RNAs. We propose IFIT1 as a cost effective and versatile tool to prepare mRNA libraries for a variety of organisms with cap 0 mRNA ends, including diverse plants, fungi and eukaryotic microbes.

The initiation of mammalian embryonic transcription: to begin at the beginning

Gamete (sperm and oocyte) genomes are transcriptionally silent until embryonic genome activation (EGA) following fertilization. EGA in humans had been thought to occur around the eight-cell stage, but recent findings suggest that it is triggered in one-cell embryos, by fertilization. Phosphorylation and other post-translational modifications during fertilization may instate transcriptionally favorable chromatin and activate oocyte-derived transcription factors (TFs) to initiate EGA. Expressed genes lay on cancer-associated pathways and their identities predict upregulation by MYC and other cancer-associated TFs. One interpretation of this is that the onset of EGA, and the somatic cell trajectory to cancer, are mechanistically related: cancer initiates epigenetically. We describe how fertilization might be linked to the initiation of EGA and involve distinctive processes recapitulated in cancer.

Surviving in the Brine: A Multi-Omics Approach for Understanding the Physiology of the Halophile Fungus Aspergillus sydowii at Saturated NaCl Concentration

Although various studies have investigated osmoadaptations of halophilic fungi to saline conditions, only few analyzed the fungal mechanisms occurring at saturated NaCl concentrations. Halophilic Aspergillus sydowii is a model organism for the study of molecular adaptations of filamentous fungi to hyperosmolarity. For the first time a multi-omics approach (i.e., transcriptomics and metabolomics) was used to compare A. sydowii at saturated concentration (5.13 M NaCl) to optimal salinity (1 M NaCl). Analysis revealed 1,842 genes differentially expressed of which 704 were overexpressed. Most differentially expressed genes were involved in metabolism and signal transduction. A gene ontology multi-scale network showed that ATP binding constituted the main network node with direct interactions to phosphorelay signal transduction, polysaccharide metabolism, and transferase activity. Free amino acids significantly decreased and amino acid metabolism was reprogrammed at 5.13 M NaCl. mRNA transcriptional analysis revealed upregulation of genes involved in methionine and cysteine biosynthesis at extreme water deprivation by NaCl. No modifications of membrane fatty acid composition occurred. Upregulated genes were involved in high-osmolarity glycerol signal transduction pathways, biosynthesis of β-1,3-glucans, and cross-membrane ion transporters. Downregulated genes were related to the synthesis of chitin, mannose, cell wall proteins, starvation, pheromone synthesis, and cell cycle. Non-coding RNAs represented the 20% of the total transcripts with 7% classified as long non-coding RNAs (lncRNAs). The 42% and 69% of the total lncRNAs and RNAs encoding transcription factors, respectively, were differentially expressed. A network analysis showed that differentially expressed lncRNAs and RNAs coding transcriptional factors were mainly related to the regulation of metabolic processes, protein phosphorylation, protein kinase activity, and plasma membrane composition. Metabolomic analyses revealed more complex and unknown metabolites at saturated NaCl concentration than at optimal salinity. This study is the first attempt to unravel the molecular ecology of an ascomycetous fungus at extreme water deprivation by NaCl (5.13 M). This work also represents a pioneer study to investigate the importance of lncRNAs and transcriptional factors in the transcriptomic response to high NaCl stress in halophilic fungi.

Human embryonic genome activation initiates at the one-cell stage

In human embryos, the initiation of transcription (embryonic genome activation [EGA]) occurs by the eight-cell stage, but its exact timing and profile are unclear. To address this, we profiled gene expression at depth in human metaphase II oocytes and bipronuclear (2PN) one-cell embryos. High-resolution single-cell RNA sequencing revealed previously inaccessible oocyte-to-embryo gene expression changes. This confirmed transcript depletion following fertilization (maternal RNA degradation) but also uncovered low-magnitude upregulation of hundreds of spliced transcripts. Gene expression analysis predicted embryonic processes including cell-cycle progression and chromosome maintenance as well as transcriptional activators that included cancer-associated gene regulators. Transcription was disrupted in abnormal monopronuclear (1PN) and tripronuclear (3PN) one-cell embryos. These findings indicate that human embryonic transcription initiates at the one-cell stage, sooner than previously thought. The pattern of gene upregulation promises to illuminate processes involved at the onset of human development, with implications for epigenetic inheritance, stem-cell-derived embryos, and cancer.

•

Gene expression initiates at the one-cell stage in human embryos

•

Expression is of low magnitude but remains elevated until the eight-cell stage

•

Upregulated transcripts are spliced and correspond to embryonic processes

•

Upregulation is disrupted in morphologically abnormal one-cell embryos

Global approaches for profiling transcription initiation

Transcription start site (TSS) selection influences transcript stability and translation as well as protein sequence. Alternative TSS usage is pervasive in organismal development, is a major contributor to transcript isoform diversity in humans, and is frequently observed in human diseases including cancer. In this review, we discuss the breadth of techniques that have been used to globally profile TSSs and the resulting insights into gene regulation, as well as future prospects in this area of inquiry.

Comparisons among rainbow trout, Oncorhynchus mykiss, populations of maternal transcript profile associated with egg viability

Background

Transcription is arrested in the late stage oocyte and therefore the maternal transcriptome stored in the oocyte provides nearly all the mRNA required for oocyte maturation, fertilization, and early cleavage of the embryo. The transcriptome of the unfertilized egg, therefore, has potential to provide markers for predictors of egg quality and diagnosing problems with embryo production encountered by fish hatcheries. Although levels of specific transcripts have been shown to associate with measures of egg quality, these differentially expressed genes (DEGs) have not been consistent among studies. The present study compares differences in select transcripts among unfertilized rainbow trout eggs of different quality based on eyeing rate, among 2 year classes of the same line (A1, A2) and a population from a different hatchery (B). The study compared 65 transcripts previously reported to be differentially expressed with egg quality in rainbow trout.

Results

There were 32 transcripts identified as DEGs among the three groups by regression analysis. Group A1 had the most DEGs, 26; A2 had 15, 14 of which were shared with A1; and B had 12, 7 of which overlapped with A1 or A2. Six transcripts were found in all three groups, dcaf11, impa2, mrpl39_like, senp7, tfip11 and uchl1.

Conclusions

Our results confirmed maternal transcripts found to be differentially expressed between low- and high-quality eggs in one population of rainbow trout can often be found to overlap with DEGs in other populations. The transcripts differentially expressed with egg quality remain consistent among year classes of the same line. Greater similarity in dysregulated transcripts within year classes of the same line than among lines suggests patterns of transcriptome dysregulation may provide insight into causes of decreased viability within a hatchery population. Although many DEGs were identified, for each of the genes there is considerable variability in transcript abundance among eggs of similar quality and low correlations between transcript abundance and eyeing rate, making it highly improbable to predict the quality of a single batch of eggs based on transcript abundance of just a few genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07773-1.

High resolution biosensor to test the capping level and integrity of mRNAs

5′ Cap structures are ubiquitous on eukaryotic mRNAs, essential for post-transcriptional processing, translation initiation and stability. Here we describe a biosensor designed to detect the presence of cap structures on mRNAs that is also sensitive to mRNA degradation, so uncapped or degraded mRNAs can be detected in a single step. The biosensor is based on a chimeric protein that combines the recognition and transduction roles in a single molecule. The main feature of this sensor is its simplicity, enabling semi-quantitative analyses of capping levels with minimal instrumentation. The biosensor was demonstrated to detect the capping level on several in vitro transcribed mRNAs. Its sensitivity and dynamic range remained constant with RNAs ranging in size from 250 nt to approximately 2700 nt and the biosensor was able to detect variations in the capping level in increments of at least 20%, with a limit of detection of 2.4 pmol. Remarkably, it also can be applied to more complex analytes, such mRNA vaccines and mRNAs transcribed in vivo. This biosensor is an innovative example of a technology able to detect analytically challenging structures such as mRNA caps. It could find application in a variety of scenarios, from quality analysis of mRNA-based products such as vaccines to optimization of in vitro capping reactions.

Molecular asymmetry in the cephalochordate embryo revealed by single-blastomere transcriptome profiling

Studies in various animals have shown that asymmetrically localized maternal transcripts play important roles in axial patterning and cell fate specification in early embryos. However, comprehensive analyses of the maternal transcriptomes with spatial information are scarce and limited to a handful of model organisms. In cephalochordates (amphioxus), an early branching chordate group, maternal transcripts of germline determinants form a compact granule that is inherited by a single blastomere during cleavage stages. Further blastomere separation experiments suggest that other transcripts associated with the granule are likely responsible for organizing the posterior structure in amphioxus; however, the identities of these determinants remain unknown. In this study, we used high-throughput RNA sequencing of separated blastomeres to examine asymmetrically localized transcripts in two-cell and eight-cell stage embryos of the amphioxus Branchiostoma floridae. We identified 111 and 391 differentially enriched transcripts at the 2-cell stage and the 8-cell stage, respectively, and used in situ hybridization to validate the spatial distribution patterns for a subset of these transcripts. The identified transcripts could be categorized into two major groups: (1) vegetal tier/germ granule-enriched and (2) animal tier/anterior-enriched transcripts. Using zebrafish as a surrogate model system, we showed that overexpression of one animal tier/anterior-localized amphioxus transcript, zfp665, causes a dorsalization/anteriorization phenotype in zebrafish embryos by downregulating the expression of the ventral gene, eve1, suggesting a potential function of zfp665 in early axial patterning. Our results provide a global transcriptomic blueprint for early-stage amphioxus embryos. This dataset represents a rich platform to guide future characterization of molecular players in early amphioxus development and to elucidate conservation and divergence of developmental programs during chordate evolution.

Parallel poly(A) homo- and hetero-duplex formation detection with an adapted DNA nanoswitch technique

Polyriboadenylic [poly(rA)] strands of sufficient length form parallel double helices in acidic and/or ammonium-containing conditions. Poly(rA) duplexes in acidic conditions are held together by A⁺-A⁺ base-pairing also involving base interactions with the phosphate backbone. Traditional UV-melting studies of parallel poly(A) duplexes have typically examined homo-duplex formation of a single nucleic acid species in solution. We have adapted a technique utilizing a DNA nanoswitch that detects interaction of two different strands either with similar or differing lengths or modifications. Our method detected parallel duplex formation as a function of length, chemical modifications, and pH, and at a sensitivity that required over 100-fold less concentration of sample than prior UV-melting methods. While parallel polyriboadenylic acid and poly-2'-O-methyl-adenylic acid homo-duplexes formed, we did not detect homo-duplexes of polydeoxyriboadenylic acid strands or poly-locked nucleic acid (LNA)-adenylic strands. Importantly however, a poly-locked nucleic acid (LNA)-adenylic strand, as well as a poly-2'-O-methyl-adenylic strand, formed a hetero-duplex with a polyriboadenylic strand. Overall, our work validates a new tool for studying parallel duplexes and reveals fundamental properties of poly(A) parallel duplex formation. Parallel duplexes may find use in DNA nanotechnology and in molecular biology applications such as a potential poly(rA) tail capture tool as an alternative to traditional oligo(dT) based purification.

Genome-wide analysis identifies cis-acting elements regulating mRNA polyadenylation and translation during vertebrate oocyte maturation

Most cells change patterns of gene expression through transcriptional regulation. In contrast, oocytes are transcriptionally silent and regulate mRNA poly(A) tail length to control protein production. However, the genome-wide relationship of poly(A) tail changes to mRNA translation during vertebrate oocyte maturation is not known. We used Tail-seq and polyribosome analysis to measure poly(A) tail and translational changes during oocyte maturation in Xenopus laevis We identified large-scale poly(A) and translational changes during oocyte maturation, with poly(A) tail length changes preceding translational changes. Proteins important for completion of the meiotic divisions and early development exhibited increased polyadenylation and translation during oocyte maturation. A family of U-rich sequence elements was enriched near the polyadenylation signal of polyadenylated and translationally activated mRNAs. We propose that changes in mRNA polyadenylation are a conserved mechanism regulating protein expression during vertebrate oocyte maturation and that these changes are controlled by a spatial code of cis-acting sequence elements.

Functional characterization of the translation initiation factor eIF4E of Echinococcus granulosus

The eukaryotic initiation factor 4E (eIF4E) specifically recognizes the 5' mRNA cap, a rate-limiting step in the translation initiation process. Although the 7-methylguanosine cap (MMGcap) is the most common 5' cap structure in eukaryotes, the trans-splicing process that occurs in several organism groups, including nematodes and flatworms, leads to the addition of a trimethylguanosine cap (TMGcap) to some RNA transcripts. In some helminths, eIF4E can have a dual capacity to bind both MMGcap and TMGcap. In the present work, we evaluated the distribution of eIF4E protein sequences in platyhelminths and we showed that only one gene coding for eIF4E is present in most parasitic flatworms. Based on this result, we cloned the Echinococcus granulosus cDNA sequence encoding eIF4E in Escherichia coli, expressed the recombinant eIF4E as a fusion protein to GST, and tested its ability to capture mRNAs through the 5' cap using pull-down assay and qPCR. Our results indicate that the recombinant eIF4E was able to bind preferentially 5'-capped mRNAs compared with rRNAs from total RNA preparations of E. granulosus. By qPCR, we observed an enrichment in MMG-capped mRNA compared with TMG-capped mRNAs among Eg-eIF4E-GST pull-down RNAs. Eg-eIF4E structural model using the Schistosoma mansoni eIF4E as template showed to be well preserved with only a few differences between chemically similar amino acid residues at the binding sites. These data showed that E. granulosus eIF4E can be used as a potential tool to study full-length 5'-capped mRNA, besides being a potential drug target against parasitic flatworms.

Transcriptome analysis of egg viability in rainbow trout, Oncorhynchus mykiss

Background

Maternal transcripts are accumulated in the oocyte during oogenesis to provide for protein synthesis from oocyte maturation through early embryonic development, when nuclear transcription is silenced. The maternal mRNAs have short poly(A) tails after undergoing post-transcriptional processing necessary for stabilizing them for storage. The transcripts undergo cytoplasmic polyadenylation when they are to be translated. Transcriptome analyses comparing total mRNA and elongated poly(A) mRNA content among eggs of different quality can provide insight into molecular mechanisms affecting egg developmental competence in rainbow trout. The present study used RNA-seq to compare transcriptomes of unfertilized eggs of rainbow trout females yielding different eyeing rates, following rRNA removal and poly(A) retention for construction of the libraries.

Results

The percentage of embryos to reach the 32-cell stage at 24 h post fertilization was significantly correlated to family eyeing rate, indicating that inviable embryos were developmentally compromised before zygotic genome activation. RNA sequencing identified 2 differentially expressed transcripts (DETs) from total mRNA sequencing comparing females with low-quality (< 5% eyeing), medium-quality (30–50% eyeing), and high-quality (> 80% eyeing) eggs. In contrast, RNA sequencing from poly(A) captured transcripts identified 945 DETs between low- and high-quality eggs, 1012 between low- and medium-quality eggs, and only 2 between medium- and high-quality eggs. The transcripts of mitochondrial genes were enriched with polyadenylated transcript sequencing and they were significantly reduced in low-quality eggs. Similarly, mitochondrial DNA was reduced in low-quality eggs compared with medium- and high-quality eggs. The functional gene analysis classified the 945 DETs between low- and high-quality eggs into 31 functional modules, many of which were related to ribosomal and mitochondrial functions. Other modules involved transcription, translation, cell division, apoptosis, and immune responses.

Conclusions

Our results indicate that differences in egg quality may be derived from differences in maternal nuclear transcript activation and cytoplasmic polyadenylation before ovulation, as opposed to accumulation and storage of maternal nuclear transcripts during oogenesis. Transcriptome comparisons suggest low-quality eggs suffered from impaired oxidative phosphorylation and translation. The DETs identified in this study provide insight into developmental competence in rainbow trout eggs.

Electronic supplementary material

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

Analysis of Transcriptome and Epitranscriptome in Plants Using PacBio Iso-Seq and Nanopore-Based Direct RNA Sequencing

Nanopore sequencing from Oxford Nanopore Technologies (ONT) and Pacific BioSciences (PacBio) single-molecule real-time (SMRT) long-read isoform sequencing (Iso-Seq) are revolutionizing the way transcriptomes are analyzed. These methods offer many advantages over most widely used high-throughput short-read RNA sequencing (RNA-Seq) approaches and allow a comprehensive analysis of transcriptomes in identifying full-length splice isoforms and several other post-transcriptional events. In addition, direct RNA-Seq provides valuable information about RNA modifications, which are lost during the PCR amplification step in other methods. Here, we present a comprehensive summary of important applications of these technologies in plants, including identification of complex alternative splicing (AS), full-length splice variants, fusion transcripts, and alternative polyadenylation (APA) events. Furthermore, we discuss the impact of the newly developed nanopore direct RNA-Seq in advancing epitranscriptome research in plants. Additionally, we summarize computational tools for identifying and quantifying full-length isoforms and other co/post-transcriptional events and discussed some of the limitations with these methods. Sequencing of transcriptomes using these new single-molecule long-read methods will unravel many aspects of transcriptome complexity in unprecedented ways as compared to previous short-read sequencing approaches. Analysis of plant transcriptomes with these new powerful methods that require minimum sample processing is likely to become the norm and is expected to uncover novel co/post-transcriptional gene regulatory mechanisms that control biological outcomes during plant development and in response to various stresses.

Identification of critical sex-biased genes in Andrias davidianus by de novo transcriptome

The Chinese giant salamander Andrias davidianus is a protected amphibian with high nutritional and economic value. Understanding its sex determination mechanism is important for improving culture techniques and sex control in breeding. However, little information on the characterization of critical genes involved in sex is available. Herein, sequencing of ovary and test produced 40,783,222 and 46,128,902 raw reads, respectively, which were jointly assembled into 80,497 unigenes. Of these, 36,609 unigenes were annotated, of which 8907 were female-biased and 10,385 were male-biased. Several sex-related pathways were observed, including the Wnt signaling pathway. After elevated temperature and estrogen exposure, neomale and neofemale specimens were identified by a female-specific marker for the first time. RT-qPCR analysis showed the expression profile of ten selected sex-biased genes to be exhibited consistently in male and neomale and in female and neofemale, with the exception of the Amh and TfIIIa genes. Results suggested that these genes may play important roles in A. davidianus sex determination and gonad development. This provides a basis for further investigation of the molecular mechanisms of sex determination in amphibians.

Poly A tail length analysis of in vitro transcribed mRNA by LC-MS

The 3'-polyadenosine (poly A) tail of in vitro transcribed (IVT) mRNA was studied using liquid chromatography coupled to mass spectrometry (LC-MS). Poly A tails were cleaved from the mRNA using ribonuclease T1 followed by isolation with dT magnetic beads. Extracted tails were then analyzed by LC-MS which provided tail length information at single-nucleotide resolution. A 2100-nt mRNA with plasmid-encoded poly A tail lengths of either 27, 64, 100, or 117 nucleotides was used for these studies as enzymatically added poly A tails showed significant length heterogeneity. The number of As observed in the tails closely matched Sanger sequencing results of the DNA template, and even minor plasmid populations with sequence variations were detected. When the plasmid sequence contained a discreet number of poly As in the tail, analysis revealed a distribution that included tails longer than the encoded tail lengths. These observations were consistent with transcriptional slippage of T7 RNAP taking place within a poly A sequence. The type of RNAP did not alter the observed tail distribution, and comparison of T3, T7, and SP6 showed all three RNAPs produced equivalent tail length distributions. The addition of a sequence at the 3' end of the poly A tail did, however, produce narrower tail length distributions which supports a previously described model of slippage where the 3' end can be locked in place by having a G or C after the poly nucleotide region. Graphical abstract Determination of mRNA poly A tail length using magnetic beads and LC-MS.

Identification and Characterization of Mitotic Spindle-Localized Transcripts

RNAs associate with the mitotic spindle in a variety of organisms, where they can spatially regulate protein production, ensure their proper segregation during cell division, or perform translation-independent roles in spindle formation. The identification of spindle-associated RNAs is an important first step in understanding the biological consequences of this phenomenon. In this chapter, we describe a method to use Xenopus laevis egg extracts to assemble and isolate mitotic spindles and to identify the spindle-associated RNAs. The method described here can be used in combination with immunodepletions, the addition of inhibitors, or other perturbations to investigate factors that affect RNA localization to the spindle. Finally, we describe a method to assess the consequences of ablating RNA in the extract on spindle formation.

Combined ultra-low input mRNA and whole-genome sequencing of human embryonic stem cells

Background

Next Generation Sequencing has proven to be an exceptionally powerful tool in the field of genomics and transcriptomics. With recent development it is nowadays possible to analyze ultra-low input sample material down to single cells. Nevertheless, investigating such sample material often limits the analysis to either the genome or transcriptome. We describe here a combined analysis of both types of nucleic acids from the same sample material.

Methods

The method described enables the combined preparation of amplified cDNA as well as amplified whole-genome DNA from an ultra-low input sample material derived from a sub-colony of in-vitro cultivated human embryonic stem cells. cDNA is prepared by the application of oligo-dT coupled magnetic beads for mRNA capture, first strand synthesis and 3’-tailing followed by PCR. Whole-genome amplified DNA is prepared by Phi29 mediated amplification. Illumina sequencing is applied to short fragment libraries prepared from the amplified samples.

Results

We developed a protocol which enables the combined analysis of the genome as well as the transcriptome by Next Generation Sequencing from ultra-low input samples. The protocol was evaluated by sequencing sub-colony structures from human embryonic stem cells containing 150 to 200 cells. The method can be adapted to any available sequencing system.

Conclusions

To our knowledge, this is the first report where sub-colonies of human embryonic stem cells have been analyzed both at the genomic as well as transcriptome level. The method of this proof of concept study may find useful practical applications for cases where only a limited number of cells are available, e.g. for tissues samples from biopsies, tumor spheres, circulating tumor cells and cells from early embryonic development. The results we present demonstrate that a combined analysis of genomic DNA and messenger RNA from ultra-low input samples is feasible and can readily be applied to other cellular systems with limited material available.

De novo Transcriptome Assemblies of Rana (Lithobates) catesbeiana and Xenopus laevis Tadpole Livers for Comparative Genomics without Reference Genomes

In this work we studied the liver transcriptomes of two frog species, the American bullfrog (Rana (Lithobates) catesbeiana) and the African clawed frog (Xenopus laevis). We used high throughput RNA sequencing (RNA-seq) data to assemble and annotate these transcriptomes, and compared how their baseline expression profiles change when tadpoles of the two species are exposed to thyroid hormone. We generated more than 1.5 billion RNA-seq reads in total for the two species under two conditions as treatment/control pairs. We de novo assembled these reads using Trans-ABySS to reconstruct reference transcriptomes, obtaining over 350,000 and 130,000 putative transcripts for R. catesbeiana and X. laevis, respectively. Using available genomics resources for X. laevis, we annotated over 97% of our X. laevis transcriptome contigs, demonstrating the utility and efficacy of our methodology. Leveraging this validated analysis pipeline, we also annotated the assembled R. catesbeiana transcriptome. We used the expression profiles of the annotated genes of the two species to examine the similarities and differences between the tadpole liver transcriptomes. We also compared the gene ontology terms of expressed genes to measure how the animals react to a challenge by thyroid hormone. Our study reports three main conclusions. First, de novo assembly of RNA-seq data is a powerful method for annotating and establishing transcriptomes of non-model organisms. Second, the liver transcriptomes of the two frog species, R. catesbeiana and X. laevis, show many common features, and the distribution of their gene ontology profiles are statistically indistinguishable. Third, although they broadly respond the same way to the presence of thyroid hormone in their environment, their receptor/signal transduction pathways display marked differences.

The calcium-dependent ribonuclease XendoU promotes ER network formation through local RNA degradation

In both Xenopus laevis egg extract and human cells, an increase in cytosolic calcium activates the endogenous ribonuclease XendoU/hEndoU, which localizes to the ER, promotes RNA cleavage and RNP removal, and induces ER network assembly.

How cells shape and remodel organelles in response to cellular signals is a poorly understood process. Using Xenopus laevis egg extract, we found that increases in cytosolic calcium lead to the activation of an endogenous ribonuclease, XendoU. A fraction of XendoU localizes to the endoplasmic reticulum (ER) and is required for nuclear envelope assembly and ER network formation in a catalysis-dependent manner. Using a purified vesicle fusion assay, we show that XendoU functions on the surface of ER membranes to promote RNA cleavage and ribonucleoprotein (RNP) removal. Additionally, RNA removal from the surface of vesicles by RNase treatment leads to increased ER network formation. Using human tissue culture cells, we found that hEndoU localizes to the ER, where it promotes the formation of ER tubules in a catalysis-dependent manner. Together, these results demonstrate that calcium-activated removal of RNA from membranes by XendoU promotes and refines ER remodeling and the formation of tubular ER.