Identification of novel transcripts with differential dorso-ventral expression in Xenopus gastrula using serial analysis of gene expression

Cell lineage-guided mass spectrometry reveals increased energy metabolism and reactive oxygen species in the vertebrate organizer

Molecular understanding of the vertebrate Organizer, a tissue center critical for inductive signaling during gastrulation, has so far been mostly limited to transcripts and a few proteins, the latter due to limitations in detection and sensitivity. The Spemann-Mangold Organizer (SMO) in the South African Clawed Frog (X. laevis), a popular model of development, has long been known to be the origin of signals that pattern the mesoderm and central nervous system. Molecular screens of the SMO have identified several genes responsible for the ability of the SMO to establish the body axis. Nonetheless, a comprehensive study of proteins and metabolites produced specifically in the SMO and their functional roles has been lacking. Here, we pioneer a deep discovery proteomic and targeted metabolomic screen of the SMO in comparison to the remainder of the embryo using high-resolution mass spectrometry (HRMS). Quantification of ~4,600 proteins and a panel of targeted metabolites documented differential expression for 460 proteins and multiple intermediates of energy metabolism in the SMO. Upregulation of oxidative phosphorylation and redox regulatory proteins gave rise to elevated oxidative stress and an accumulation of reactive oxygen species in the SMO. Imaging experiments corroborated these findings, discovering enrichment of hydrogen peroxide in the SMO. Chemical perturbation of the redox gradient perturbed mesoderm involution during early gastrulation. HRMS expands the bioanalytical toolbox of cell and developmental biology, providing previously unavailable information on molecular classes to challenge and refine our classical understanding of the Organizer and its function during early patterning of the embryo.

Transcriptome analysis of the almond moth, Cadra cautella, female abdominal tissues and identification of reproduction control genes

BACKGROUND

The almond moth, Cadra cautella is a destructive pest of stored food commodities including dates that causes severe economic losses for the farming community worldwide. To date, no genetic information related to the molecular mechanism/strategies of its reproduction is available. Thus, transcriptome analysis of C. cautella female abdominal tissues was performed via next-generation sequencing (NGS) to recognize the genes responsible for reproduction.

RESULTS

The NGS was performed with an Illumina Hiseq 2000 sequencer (Beijing Genomics Institute: BGI). From the transcriptome data, 9,804,804,120 nucleotides were generated and their assemblage resulted in 62,687 unigenes. The functional annotation analyses done by different databases, annotated, 27,836 unigenes in total. The transcriptome data of C. cautella female abdominal tissue was submitted to the National Center for Biotechnology Information (accession no: PRJNA484692). The transcriptome analysis yielded several genes responsible for C. cautella reproduction including six Vg gene transcripts. Among the six Vg gene transcripts, only one was highly expressed with 3234.95 FPKM value (fragments per kilobase per million mapped reads) that was much higher than that of the other five transcripts. Higher differences in the expression level of the six Vg transcripts were confirmed by running the RT-PCR using gene specific primers, where the expression was observed only in one transcript it was named as the CcVg.

CONCLUSIONS

This is the first study to explore C. cautella reproduction control genes and it might be supportive to explore the reproduction mechanism in this pest at the molecular level. The NGS based transcriptome pool is valuable to study the functional genomics and will support to design biotech-based management strategies for C. cautella.

Spemann organizer transcriptome induction by early beta-catenin, Wnt, Nodal, and Siamois signals in Xenopus laevis

The earliest event in Xenopus development is the dorsal accumulation of nuclear β-catenin under the influence of cytoplasmic determinants displaced by fertilization. In this study, a genome-wide approach was used to examine transcription of the 43,673 genes annotated in the Xenopus laevis genome under a variety of conditions that inhibit or promote formation of the Spemann organizer signaling center. Loss of function of β-catenin with antisense morpholinos reproducibly reduced the expression of 247 mRNAs at gastrula stage. Interestingly, only 123 β-catenin targets were enriched on the dorsal side and defined an early dorsal β-catenin gene signature. These genes included several previously unrecognized Spemann organizer components. Surprisingly, only 3 of these 123 genes overlapped with the late Wnt signature recently defined by two other groups using inhibition by Dkk1 mRNA or Wnt8 morpholinos, which indicates that the effects of β-catenin/Wnt signaling in early development are exquisitely regulated by stage-dependent mechanisms. We analyzed transcriptome responses to a number of treatments in a total of 46 RNA-seq libraries. These treatments included, in addition to β-catenin depletion, regenerating dorsal and ventral half-embryos, lithium chloride treatment, and the overexpression of Wnt8, Siamois, and Cerberus mRNAs. Only some of the early dorsal β-catenin signature genes were activated at blastula whereas others required the induction of endomesoderm, as indicated by their inhibition by Cerberus overexpression. These comprehensive data provide a rich resource for analyzing how the dorsal and ventral regions of the embryo communicate with each other in a self-organizing vertebrate model embryo.

Towards the bridging of molecular genetics data across Xenopus species

BACKGROUND

The clawed African frog Xenopus laevis has been one of the main vertebrate models for studies in developmental biology. However, for genetic studies, Xenopus tropicalis has been the experimental model of choice because it shorter life cycle and due to a more tractable genome that does not result from genome duplication as in the case of X. laevis. Today, although still organized in a large number of scaffolds, nearly 85% of X. tropicalis and 89% of X. laevis genomes have been sequenced. There is expectation for a comparative physical map that can be used as a Rosetta Stone between X. laevis genetic studies and X. tropicalis genomic research.

RESULTS

In this work, we have mapped using coarse-grained alignment the 18 chromosomes of X. laevis, release 9.1, on the 10 reference scaffolds representing the haploid genome of X. tropicalis, release 9.0. After validating the mapping with theoretical data, and estimating reference averages of genome sequence identity, 37 to 44% between the two species, we have carried out a synteny analysis for 2,112 orthologous genes. We found that 99.6% of genes are in the same organization.

CONCLUSIONS

Taken together, our results make possible to establish the correspondence between 62 and 65.5% of both genomes, percentage of identity, synteny and automatic annotation of transcripts of both species, providing a new and more comprehensive tool for comparative analysis of these two species, by allowing to bridge molecular genetics data among them.

A robust method for transcript quantification with RNA-seq data

The advent of high throughput RNA-seq technology allows deep sampling of the transcriptome, making it possible to characterize both the diversity and the abundance of transcript isoforms. Accurate abundance estimation or transcript quantification of isoforms is critical for downstream differential analysis (e.g., healthy vs. diseased cells) but remains a challenging problem for several reasons. First, while various types of algorithms have been developed for abundance estimation, short reads often do not uniquely identify the transcript isoforms from which they were sampled. As a result, the quantification problem may not be identifiable, i.e., lacks a unique transcript solution even if the read maps uniquely to the reference genome. In this article, we develop a general linear model for transcript quantification that leverages reads spanning multiple splice junctions to ameliorate identifiability. Second, RNA-seq reads sampled from the transcriptome exhibit unknown position-specific and sequence-specific biases. We extend our method to simultaneously learn bias parameters during transcript quantification to improve accuracy. Third, transcript quantification is often provided with a candidate set of isoforms, not all of which are likely to be significantly expressed in a given tissue type or condition. By resolving the linear system with LASSO, our approach can infer an accurate set of dominantly expressed transcripts while existing methods tend to assign positive expression to every candidate isoform. Using simulated RNA-seq datasets, our method demonstrated better quantification accuracy and the inference of dominant set of transcripts than existing methods. The application of our method on real data experimentally demonstrated that transcript quantification is effective for differential analysis of transcriptomes.

Characterization of small RNAs in X. tropicalis gastrulae

Here, we report and characterize deep sequencing data and bioinformatics analysis of small RNAs from X. tropicalis gastrula. A total of 17,553,124 reads with perfect match to the genome derived from 2,616,053 unique sequences were identified. Seventy-seven percent of theses sequences were not found in previous reports from X. tropicalis oocytes and somatic tissues. Bioinformatics analyses indicate that a large fraction of the small RNAs are PIWI-interacting RNAs. Up to 23.9% of small RNAs mapped to transposable elements and 27% to genic regions. Most of the abundant transposon-derived small RNAs are found in oocyte and gastrula libraries, suggesting that transposons also need to be silenced during early embryonic development. Importantly, novel clusters of piRNAs whose expression is activated after zygotic transcription begins were identified in the genome of X. tropicalis. Additionally, miRNAs were also identified and many of them are not present in oocytes, suggesting that miRNA expression is stage-specific. To the best of our knowledge, this is the first high throughput data release and bioinformatics characterization of small RNAs during Xenopus early embryonic development.

Characterization of small RNAs in Xenopus tropicalis gastrulae

Here, we report and characterize deep sequencing data and bioinformatics analysis of small RNAs from Xenopus tropicalis gastrula. A total of 17,553,124 reads with perfect match to the genome derived from 2,616,053 unique sequences were identified. Seventy-seven percent of theses sequences were not found in previous reports from X. tropicalis oocytes and somatic tissues. Bioinformatics analyses indicate that a large fraction of the small RNAs are PIWI-interacting RNAs. Up to 23.9% of small RNAs mapped to transposable elements and 27% to genic regions. Most of abundant transposable derived small RNAs are found in oocyte and gastrula libraries, suggesting that transposon needs to be silenced also during early development. Additionally, miRNAs were identified and many of them are not present in oocytes, suggesting that miRNA expression is stage specific. To the best of our knowledge, this is the first high throughput data release and bioinformatics characterization of small RNAs during Xenopus development.

Analysis of transcriptome differences between resistant and susceptible strains of the citrus red mite Panonychus citri (Acari: Tetranychidae)

BACKGROUND

The citrus red mite is a worldwide citrus pest and a common sensitizing allergen of asthma and rhinitis. It has developed strong resistance to many registered acaricides, However, the molecular mechanisms of resistance remain unknown. we therefore used next generation sequencing technology to investigate the global transcriptomes between resistant strains and susceptible strains.

RESULTS

We obtained 34,159, 30,466 and 32,217 unigenes by assembling the SS reads, RS reads and SS&RS reads respectively. There are total 17,581 annotated unigenes from SS&RS reads by BLAST searching databases of nr, the Clusters of Orthologous Groups (COGs) and Kyoto Encyclopedia of Genes and Genomes (KEGG) with an E-value ≤ 1e-5, in which 7,075 unigenes were annotated in the COG database, 12, 712 unigenes were found in the KEGG database and 3,812 unigenes were assigned to Gene ontology (GO). Moreover, 2,701 unigenes were judged to be the differentially expressed genes (DEGs) based on the uniquely mapped reads. There are 219 pathways in all annotated unigenes and 198 pathways in DEGs that mapped to the KEGG database. We identified 211 metabolism genes and target genes related to general insecticide resistance such as P450 and Cytochrome b, and further compared their differences between RS and SS. Meanwhile, we identified 105 and 194 genes related to growth and reproduction, respectively, based on the mode of action of Hexythiazox. After further analyses, we found variation in sequences but not in gene expression related to mite growth and reproduction between different strains.

CONCLUSION

To our knowledge, this is the first comparative transcriptome study to discover candidate genes involved in phytophagous mite resistance. This study identified differential unigenes related to general pesticide resistance and organism growth and reproduction in P. citri. The assembled, annotated transcriptomes provide a valuable genomic resource for further understanding the molecular basis of resistance mechanisms.

Expression of DNA transposable elements during nervous system development

Transposable elements (retrotransposons and DNA transposons) comprise a large proportion of animal genomes, for example 20% in D. melanogaster, 36% in X. tropicalis and 45% in humans. After invading a new genome, the transposable element increases its copy number and subsequently accumulates mutations. These may eventually result in inactive copies. Until recent days transposons have been considered “junk” DNA and no clear function have been assigned for this important amount of information on genomes.

Expression of transposable elements in neural tissues during Xenopus development

Transposable elements comprise a large proportion of animal genomes. Transposons can have detrimental effects on genome stability but also offer positive roles for genome evolution and gene expression regulation. Proper balance of the positive and deleterious effects of transposons is crucial for cell homeostasis and requires a mechanism that tightly regulates their expression. Herein we describe the expression of DNA transposons of the Tc1/mariner superfamily during Xenopus development. Sense and antisense transcripts containing complete Tc1-2_Xt were detected in Xenopus embryos. Both transcripts were found in zygotic stages and were mainly localized in Spemann's organizer and neural tissues. In addition, the Tc1-like elements Eagle, Froggy, Jumpy, Maya, Xeminos and TXr were also expressed in zygotic stages but not oocytes in X. tropicalis. Interestingly, although Tc1-2_Xt transcripts were not detected in Xenopus laevis embryos, transcripts from other two Tc1-like elements (TXr and TXz) presented a similar temporal and spatial pattern during X. laevis development. Deep sequencing analysis of Xenopus tropicalis gastrulae showed that PIWI-interacting RNAs (piRNAs) are specifically derived from several Tc1-like elements. The localized expression of Tc1-like elements in neural tissues suggests that they could play a role during the development of the Xenopus nervous system.

Digital gene expression analysis of two life cycle stages of the human-infective parasite, Trypanosoma brucei gambiense reveals differentially expressed clusters of co-regulated genes

Background

The evolutionarily ancient parasite, Trypanosoma brucei, is unusual in that the majority of its genes are regulated post-transcriptionally, leading to the suggestion that transcript abundance of most genes does not vary significantly between different life cycle stages despite the fact that the parasite undergoes substantial cellular remodelling and metabolic changes throughout its complex life cycle. To investigate this in the clinically relevant sub-species, Trypanosoma brucei gambiense, which is the causative agent of the fatal human disease African sleeping sickness, we have compared the transcriptome of two different life cycle stages, the potentially human-infective bloodstream forms with the non-human-infective procyclic stage using digital gene expression (DGE) analysis.

Results

Over eleven million unique tags were generated, producing expression data for 7360 genes, covering 81% of the genes in the genome. Compared to microarray analysis of the related T. b. brucei parasite, approximately 10 times more genes with a 2.5-fold change in expression levels were detected. The transcriptome analysis revealed the existence of several differentially expressed gene clusters within the genome, indicating that contiguous genes, presumably from the same polycistronic unit, are co-regulated either at the level of transcription or transcript stability.

Conclusions

DGE analysis is extremely sensitive for detecting gene expression differences, revealing firstly that a far greater number of genes are stage-regulated than had previously been identified and secondly and more importantly, this analysis has revealed the existence of several differentially expressed clusters of genes present on what appears to be the same polycistronic units, a phenomenon which had not previously been observed in microarray studies. These differentially regulated clusters of genes are in addition to the previously identified RNA polymerase I polycistronic units of variant surface glycoproteins and procyclin expression sites, which encode the major surface proteins of the parasite. This raises a number of questions regarding the function and regulation of the gene clusters that clearly warrant further study.