BackgroundWith 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.ResultsTo 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.ConclusionsOur 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.

BackgroundTea 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.ResultsUsing 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.ConclusionsOur 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.

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.

The discovery of N 6 -methyldeoxyadenine (6mA) across eukaryotes led to a search for additional epigenetic mechanisms. However, some studies have highlighted confounding factors that challenge the prevalence of 6mA in eukaryotes. We developed a metagenomic method to quantitatively deconvolve 6mA events from a genomic DNA sample into species of interest, genomic regions, and sources of contamination. Applying this method, we observed high-resolution 6mA deposition in two protozoa. We found that commensal or soil bacteria explained the vast majority of 6mA in insect and plant samples. We found no evidence of high abundance of 6mA in Drosophila , Arabidopsis , or humans. Plasmids used for genetic manipulation, even those from Dam methyltransferase mutant Escherichia coli , could carry abundant 6mA, confounding the evaluation of candidate 6mA methyltransferases and demethylases. On the basis of this work, we advocate for a reassessment of 6mA in eukaryotes.

Molluscs form their shells out of CaCO 3 and a matrix of biomacromolecules. Understanding the role of matrices may shed some light on the mechanism of biomineralization. Here, a 1401-bp full-length cDNA sequence encoding a novel matrix protein was cloned from the mantle of the bivalve oyster, Pinctada fucata. The deduced protein (Prisilkin-39), which has a molecular mass of 39.3 kDa and an isoelectric point of 8.83, was fully characterized, and its role in biomineralization was demonstrated using both in vivo and in vitro crystal growth assays. Prisilkin-39 is a highly repetitive protein with an unusual composition of Gly, Tyr, and Ser residues. Expression of Prisilkin-39 was localized to columnar epithelial cells of the mantle edge, corresponding to the calcitic prismatic layer formation. Immunostaining in situ and immunodetection in vitro revealed the presence of a characteristic pattern of Prisilkin-39 in the organic sheet and in sheaths around the prisms. Prisilkin-39 binds tightly with chitin, an insoluble polysaccharide that forms the highly structured framework of the shell. Antibody injection in vivo resulted in dramatic morphological deformities in the inner shell surface structure, where large amounts of CaCO 3 were deposited in an uncontrolled manner. Moreover, Prisilkin-39 strictly prohibited the precipitation of aragonite in vitro. Taken together, Prisilkin-39 is the first protein shown to have dual function, involved both in the chitinous framework building and in crystal growth regulation during the prismatic layer mineralization. These observations may extend our view on the rare group of basic matrices and their functions during elaboration of the molluscan shell.

Mechanical stress plays a key role in regulating cartilage degradation in osteoarthritis (OA). The aim of this study was to evaluate the effects and mechanisms of mechanical stress on articular cartilage. A total of 80 male Sprague‐Dawley rats were randomly divided into eight groups ( n = 10 for each group): control group (CG), OA group (OAG), and CG or OAG subjected to low‐, moderate‐, or high‐intensity treadmill exercise (CL, CM, CH, OAL, OAM, and OAH, respectively). Chondrocytes were obtained from the knee joints of rats; they were cultured on Bioflex 6‐well culture plates and subjected to different durations of cyclic tensile strain (CTS) with or without exposure to interleukin‐1β (IL‐1β). The results of the histological score, immunohistochemistry, enzyme‐linked immunosorbent assay, and western‐blot analyses indicated that there were no differences between CM and CG, but OAM showed therapeutic effects compared with OAG. However, CH and OAH experienced more cartilage damage than CG and OAG, respectively. CTS had no therapeutic effects on collagen II of normal chondrocytes, which is consistent with findings after treadmill exercise. However, CTS for 4 hr could alleviate the chondrocyte damage induced by IL‐1β by activating AMP‐activated protein kinase (AMPK) phosphorylation and suppressing nuclear translocation of nuclear factor (NF)‐κB p65. Our findings indicate that mechanical stress had no therapeutic effects on normal articular cartilage and chondrocytes; mechanical stress only caused damage with excessive stimulation. Still, moderate biomechanical stress could reduce sensitization to the inflammatory response of articular cartilage and chondrocytes through the AMPK/NF‐κB signaling pathway.

The Wnt signaling pathway is crucial for tissue morphogenesis, participating in cellular behavior changes, notably during the process of convergent-extension. Interactions between Wnt-secreting and receiving cells during convergent-extension remain elusive. We investigated the role and genetic interactions of Wnt ligands and their trafficking factors Wls, Gpc4 and Frzb in the context of palate morphogenesis in zebrafish. We describe that the chaperon Wls and its ligands Wnt9a and Wnt5b are expressed in the ectoderm, whereas juxtaposed chondrocytes express Frzb and Gpc4. Using wls, gpc4, frzb, wnt9a and wnt5b mutants, we genetically dissected the Wnt signals operating between secreting ectoderm and receiving chondrocytes. Our analysis delineates that non-canonical Wnt signaling is required for cell intercalation, and that wnt5b and wnt9a are required for palate extension in the anteroposterior and transverse axes, respectively.