Background: With the fast advances in nextgen sequencing technology, high-throughput RNA sequencing has emerged as a powerful and cost-effective way for transcriptome study. De novo assembly of transcripts provides an important solution to transcriptome analysis for organisms with no reference genome. However, there lacked understanding on how the different variables affected assembly outcomes, and there was no consensus on how to approach an optimal solution by selecting software tool and suitable strategy based on the properties of RNA-Seq data. Results: To reveal the performance of different programs for transcriptome assembly, this work analyzed some important factors, including k-mer values, genome complexity, coverage depth, directional reads, etc. Seven program conditions, four single k-mer assemblers (SK: SOAPdenovo, ABySS, Oases and Trinity) and three multiple k-mer methods (MK: SOAPdenovo-MK, trans-ABySS and Oases-MK) were tested. While small and large k-mer values performed better for reconstructing lowly and highly expressed transcripts, respectively, MK strategy worked well for almost all ranges of expression quintiles. Among SK tools, Trinity performed well across various conditions but took the longest running time. Oases consumed the most memory whereas SOAPdenovo required the shortest runtime but worked poorly to reconstruct full-length CDS. ABySS showed some good balance between resource usage and quality of assemblies. Conclusions: Our work compared the performance of publicly available transcriptome assemblers, and analyzed important factors affecting de novo assembly. Some practical guidelines for transcript reconstruction from short-read RNA-Seq data were proposed. De novo assembly of C. sinensis transcriptome was greatly improved using some optimized methods.

Background: Tea is the most popular non-alcoholic health beverage in the world. The tea plant (Camellia sinensis (L.) O. Kuntze) needs to undergo a cold acclimation process to enhance its freezing tolerance in winter. Changes that occur at the molecular level in response to low temperatures are poorly understood in tea plants. To elucidate the molecular mechanisms of cold acclimation, we employed RNA-Seq and digital gene expression (DGE) technologies to the study of genome-wide expression profiles during cold acclimation in tea plants. Results: Using the Illumina sequencing platform, we obtained approximately 57.35 million RNA-Seq reads. These reads were assembled into 216,831 transcripts, with an average length of 356 bp and an N50 of 529 bp. In total, 1,770 differentially expressed transcripts were identified, of which 1,168 were up-regulated and 602 down-regulated. These include a group of cold sensor or signal transduction genes, cold-responsive transcription factor genes, plasma membrane stabilization related genes, osmosensing-responsive genes, and detoxification enzyme genes. DGE and quantitative RT-PCR analysis further confirmed the results from RNA-Seq analysis. Pathway analysis indicated that the “carbohydrate metabolism pathway” and the “calcium signaling pathway” might play a vital role in tea plants’ responses to cold stress. Conclusions: Our study presents a global survey of transcriptome profiles of tea plants in response to low, non-freezing temperatures and yields insights into the molecular mechanisms of tea plants during the cold acclimation process. It could also serve as a valuable resource for relevant research on cold-tolerance and help to explore the cold-related genes in improving the understanding of low-temperature tolerance and plant-environment interactions.

Representing a basal branch of arachnids, scorpions are known as ‘living fossils’ that maintain an ancient anatomy and are adapted to have survived extreme climate changes. Here we report the genome sequence of Mesobuthus martensii, containing 32,016 protein-coding genes, the most among sequenced arthropods. Although M. martensii appears to evolve conservatively, it has a greater gene family turnover than the insects that have undergone diverse morphological and physiological changes, suggesting the decoupling of the molecular and morphological evolution in scorpions. Underlying the long-term adaptation of scorpions is the expansion of the gene families enriched in basic metabolic pathways, signalling pathways, neurotoxins and cytochrome P450, and the different dynamics of expansion between the shared and the scorpion lineage-specific gene families. Genomic and transcriptomic analyses further illustrate the important genetic features associated with prey, nocturnal behaviour, feeding and detoxification. The M. martensii genome reveals a unique adaptation model of arthropods, offering new insights into the genetic bases of the living fossils.

Major secondary metabolites, including flavonoids, caffeine, and theanine, are important components of tea products and are closely related to the taste, flavor, and health benefits of tea. Secondary metabolite biosynthesis in Camellia sinensis is differentially regulated in different tissues during growth and development. Until now, little was known about the expression patterns of genes involved in secondary metabolic pathways or their regulatory mechanisms. This study aimed to generate expression profiles for C. sinensis tissues and to build a gene regulation model of the secondary metabolic pathways. RNA sequencing was performed on 13 different tissue samples from various organs and developmental stages of tea plants, including buds and leaves of different ages, stems, flowers, seeds, and roots. A total of 43.7 Gbp of raw sequencing data were generated, from which 347,827 unigenes were assembled and annotated. There were 46,693, 8446, 3814, 10,206, and 4948 unigenes specifically expressed in the buds and leaves, stems, flowers, seeds, and roots, respectively. In total, 1719 unigenes were identified as being involved in the secondary metabolic pathways in C. sinensis, and the expression patterns of the genes involved in flavonoid, caffeine, and theanine biosynthesis were characterized, revealing the dynamic nature of their regulation during plant growth and development. The possible transcription factor regulation network for the biosynthesis of flavonoid, caffeine, and theanine was built, encompassing 339 transcription factors from 35 families, namely bHLH, MYB, and NAC, among others. Remarkably, not only did the data reveal the possible critical check points in the flavonoid, caffeine, and theanine biosynthesis pathways, but also implicated the key transcription factors and related mechanisms in the regulation of secondary metabolite biosynthesis. Our study generated gene expression profiles for different tissues at different developmental stages in tea plants. The gene network responsible for the regulation of the secondary metabolic pathways was analyzed. Our work elucidated the possible cross talk in gene regulation between the secondary metabolite biosynthetic pathways in C. sinensis. The results increase our understanding of how secondary metabolic pathways are regulated during plant development and growth cycles, and help pave the way for genetic selection and engineering for germplasm improvement.

The hydrolytic deamination of adenosine to inosine (A-to-I editing) in precursor mRNA induces variable gene products at the post-transcription level. How and to what extent A-to-I RNA editing diversifies transcriptome is not fully characterized in the evolution, and very little is known about the selective constraints that drive the evolution of RNA editing events. Here we present a study on A-to-I RNA editing, by generating a global profile of A-to-I editing for a phylogeny of seven Drosophila species, a model system spanning an evolutionary timeframe of approximately 45 million years. Of totally 9281 editing events identified, 5150 (55.5%) are located in the coding sequences (CDS) of 2734 genes. Phylogenetic analysis places these genes into 1,526 homologous families, about 5% of total gene families in the fly lineages. Based on conservation of the editing sites, the editing events in CDS are categorized into three distinct types, representing events on singleton genes (type I), and events not conserved (type II) or conserved (type III) within multi-gene families. While both type I and II events are subject to purifying selection, notably type III events are positively selected, and highly enriched in the components and functions of the nervous system. The tissue profiles are documented for three editing types, and their critical roles are further implicated by their shifting patterns during holometabolous development and in post-mating response. In conclusion, three A-to-I RNA editing types are found to have distinct evolutionary dynamics. It appears that nervous system functions are mainly tested to determine if an A-to-I editing is beneficial for an organism. The coding plasticity enabled by A-to-I editing creates a new class of binary variations, which is a superior alternative to maintain heterozygosity of expressed genes in a diploid mating system. One prevalent form of RNA editing is the deamination of adenosines (A-to-I editing) in the precursor mRNA molecules, pertaining to most organisms in the metazoan lineage. While examples of A-to-I editing on critical genes have been known for years, it has not been fully characterized how A-to-I editing shapes the transcriptome and proteome in the evolution. To understand how A-to-I editing affects genes’ evolution and how itself is constrained by selection, we generated a global profile of A-to-I editing for a phylogeny of seven fly species, a model system representing an evolutionary timeframe of about 45 million years. We are focused on 5150 editing sites (of totally 9281 identified) located in the coding region of 2734 genes. Our analysis revealed the evolution dynamics of A-to-I editing sites and functional specificity of targeted genes. The shifting patterns of A-to-I editing are documented during holometabolous development and in post-mating response in flies. This work points to the important roles of regulated RNA editing in animal development and offers new insight into the evolution of A-to-I editing events and their harboring genes.