Optimization of de novo transcriptome assembly from next-generation sequencing data

Hepatic transcriptome profiling of largemouth bass (Micropterus salmoides Lacépède) injected with Flavobacterium covae or lipopolysaccharide

Flavobacterium covae (columnaris) is the most detrimental bacterial disease affecting the largemouth bass (Micropterus salmoides Lacépède) aquaculture industry. In the current study, fish received an intraperitoneal injection of either 1× PBS (100 μL), LPS in PBS (100 μL, 10 μg/mL), or F. covae (100 μL, 2.85 × 1011 CFU/mL) to simulate immunological challenges. After 24 h post-injection, liver tissue from the control and treated groups were then collected for transcriptome analysis. Results of the Gene Ontology (GO) and KEGG pathway analyses for the F. covae and LPS-injected groups found differentially expressed genes (DEGs) enriched primarily in toll-like receptors (TLRs), cytokine-cytokine receptors, complement and coagulation cascades, and the PPAR signalling pathways. This suggests that the liver immune system is enhanced by these five combined pathways. Additionally, the DEGs TLR5, MYD88, and IL-1 were significantly upregulated in F. covae and LPS-injected fish compared to the controls, whereas IL-8 was downregulated. The upregulation of TLR5 was unexpected as F. covae lacks flagellin, the protein that binds to TLR5. Additionally, it is unknown whether the TLR5 is upregulated by LPS. Further research into the upregulation of TLR5 is warranted. These results provide insight into immune responses and associated pathways contributing to the immune system in the liver during columnaris infection and induced response to LPS in largemouth bass.

Roast: a tool for reference-free optimization of supertranscriptome assemblies

BACKGROUND

Transcriptomic studies involving organisms for which reference genomes are not available typically start by generating de novo transcriptome or supertranscriptome assembly from the raw RNA-seq reads. Assembling a supertranscriptome is, however, a challenging task due to significantly varying abundance of mRNA transcripts, alternative splicing, and sequencing errors. As a result, popular de novo supertranscriptome assembly tools generate assemblies containing contigs that are partially-assembled, fragmented, false chimeras or have local mis-assemblies leading to decreased assembly accuracy. Commonly available tools for assembly improvement rely primarily on running BLAST using closely related species making their accuracy and reliability conditioned on the availability of the data for closely related organisms.

RESULTS

We present ROAST, a tool for optimization of supertranscriptome assemblies that uses paired-end RNA-seq data from Illumina sequencing platform to iteratively identify and fix assembly errors solely using the error signatures generated by RNA-seq alignment tools including soft-clips, unexpected expression coverage, and reads with mates unmapped or mapped on a different contig to identify and fix various supertranscriptome assembly errors without performing BLAST searches against other organisms. Evaluation results using simulated as well as real datasets show that ROAST significantly improves assembly quality by identifying and fixing various assembly errors.

CONCLUSION

ROAST provides a reference-free approach to optimizing supertranscriptome assemblies highlighting its utility in refining de novo supertranscriptome assemblies of non-model organisms.

TransIntegrator: capture nearly full protein-coding transcript variants via integrating Illumina and PacBio transcriptomes

Genes have the ability to produce transcript variants that perform specific cellular functions. However, accurately detecting all transcript variants remains a long-standing challenge, especially when working with poorly annotated genomes or without a known genome. To address this issue, we have developed a new computational method, TransIntegrator, which enables transcriptome-wide detection of novel transcript variants. For this, we determined 10 Illumina sequencing transcriptomes and a PacBio full-length transcriptome for consecutive embryo development stages of amphioxus, a species of great evolutionary importance. Based on the transcriptomes, we employed TransIntegrator to create a comprehensive transcript variant library, namely iTranscriptome. The resulting iTrancriptome contained 91 915 distinct transcript variants, with an average of 2.4 variants per gene. This substantially improved current amphioxus genome annotation by expanding the number of genes from 21 954 to 38 777. Further analysis manifested that the gene expansion was largely ascribed to integration of multiple Illumina datasets instead of involving the PacBio data. Moreover, we demonstrated an example application of TransIntegrator, via generating iTrancriptome, in aiding accurate transcriptome assembly, which significantly outperformed other hybrid methods such as IDP-denovo and Trinity. For user convenience, we have deposited the source codes of TransIntegrator on GitHub as well as a conda package in Anaconda. In summary, this study proposes an affordable but efficient method for reliable transcriptomic research in most species.

Metabolomics and transcriptomics analyses for characterizing the alkaloid metabolism of Chinese jujube and sour jujube fruits

Introduction

Jujube is an important economic forest tree whose fruit is rich in alkaloids. Chinese jujube (Ziziphus jujuba Mill.) and sour jujube (Ziziphus spinosa Hu.) are the two most important species of the jujube genus. However, the mechanisms underlying the synthesis and metabolism of alkaloids in jujube fruits remain poorly understood.

Methods

In this study, the fruits of Ziziphus jujuba 'Hupingzao' and Ziziphus spinosa 'Taigusuanzao' in different harvest stages were used as test materials, we first integrated widely targeted metabolomics and transcriptomics analyses to elucidate the metabolism of alkaloids of jujube fruits.

Results

In the metabolomics analysis, 44 alkaloid metabolites were identified in 4 samples, 3 of which were unique to sour jujube fruit. The differential alkaloid metabolites (DAMs) were more accumulated in sour jujube than in Chinese jujube; further, they were more accumulated in the white ripening stage than in the red stage. DAMs were annotated to 12 metabolic pathways. Additionally, transcriptomics data revealed 259 differentially expressed genes (DEGs) involved in alkaloid synthesis and metabolism. By mapping the regulatory networks of DAMs and DEGs, we screened out important metabolites and 11 candidate genes.

Discussion

This study preliminarily elucidated the molecular mechanism of jujube alkaloid synthesis. The candidate genes regulated the synthesis of key alkaloid metabolites, but the specific regulation mechanism is unclear. Taken together, our results provide insights into the metabolic networks of alkaloid synthesis in Chinese jujube and sour jujube fruits at different harvest stages, thereby providing a theoretical reference for further research on the regulatory mechanism of jujube alkaloids and their development and utilization.

Optimizing an efficient ensemble approach for high-quality de novo transcriptome assembly of Thymus daenensis

Non-erroneous and well-optimized transcriptome assembly is a crucial prerequisite for authentic downstream analyses. Each de novo assembler has its own algorithm-dependent pros and cons to handle the assembly issues and should be specifically tested for each dataset. Here, we examined efficiency of seven state-of-art assemblers on ~ 30 Gb data obtained from mRNA-sequencing of Thymus daenensis. In an ensemble workflow, combining the outputs of different assemblers associated with an additional redundancy-reducing step could generate an optimized outcome in terms of completeness, annotatability, and ORF richness. Based on the normalized scores of 16 benchmarking metrics, EvidentialGene, BinPacker, Trinity, rnaSPAdes, CAP3, IDBA-trans, and Velvet-Oases performed better, respectively. EvidentialGene, as the best assembler, totally produced 316,786 transcripts, of which 235,730 (74%) were predicted to have a unique protein hit (on uniref100), and also half of its transcripts contained an ORF. The total number of unique BLAST hits for EvidentialGene was approximately three times greater than that of the worst assembler (Velvet-Oases). EvidentialGene could even capture 17% and 7% more average BLAST hits than BinPacker and Trinity. Although BinPacker and CAP3 produced longer transcripts, the EvidentialGene showed a higher collinearity between transcript size and ORF length. Compared with the other programs, EvidentialGene yielded a higher number of optimal transcript sets, further full-length transcripts, and lower possible misassemblies. Our finding corroborates that in non-model species, relying on a single assembler may not give an entirely satisfactory result. Therefore, this study proposes an ensemble approach of accompanying EvidentialGene pipelines to acquire a superior assembly for T. daenensis.

Elucidating the Mesocarp Drupe Transcriptome of Açai (Euterpe oleracea Mart.): An Amazonian Tree Palm Producer of Bioactive Compounds

Euterpe oleracea palm, endemic to the Amazon region, is well known for açai, a fruit violet beverage with nutritional and medicinal properties. During E. oleracea fruit ripening, anthocyanin accumulation is not related to sugar production, contrarily to grape and blueberry. Ripened fruits have a high content of anthocyanins, isoprenoids, fibers, and proteins, and are poor in sugars. E. oleracea is proposed as a new genetic model for metabolism partitioning in the fruit. Approximately 255 million single-end-oriented reads were generated on an Ion Proton NGS platform combining fruit cDNA libraries at four ripening stages. The de novo transcriptome assembly was tested using six assemblers and 46 different combinations of parameters, a pre-processing and a post-processing step. The multiple k-mer approach with TransABySS as an assembler and Evidential Gene as a post-processer have shown the best results, with an N50 of 959 bp, a read coverage mean of 70x, a BUSCO complete sequence recovery of 36% and an RBMT of 61%. The fruit transcriptome dataset included 22,486 transcripts representing 18 Mbp, of which a proportion of 87% had significant homology with other plant sequences. Approximately 904 new EST-SSRs were described, and were common and transferable to Phoenix dactylifera and Elaeis guineensis, two other palm trees. The global GO classification of transcripts showed similar categories to that in P. dactylifera and E. guineensis fruit transcriptomes. For an accurate annotation and functional description of metabolism genes, a bioinformatic pipeline was developed to precisely identify orthologs, such as one-to-one orthologs between species, and to infer multigenic family evolution. The phylogenetic inference confirmed an occurrence of duplication events in the Arecaceae lineage and the presence of orphan genes in E. oleracea. Anthocyanin and tocopherol pathways were annotated entirely. Interestingly, the anthocyanin pathway showed a high number of paralogs, similar to in grape, whereas the tocopherol pathway exhibited a low and conserved gene number and the prediction of several splicing forms. The release of this exhaustively annotated molecular dataset of E. oleracea constitutes a valuable tool for further studies in metabolism partitioning and opens new great perspectives to study fruit physiology with açai as a model.

Normalized Workflow to Optimize Hybrid De Novo Transcriptome Assembly for Non-Model Species: A Case Study in Lilium ledebourii (Baker) Boiss

A high-quality transcriptome is required to advance numerous bioinformatics workflows. Nevertheless, the effectuality of tools for de novo assembly and real precision assembled transcriptomes looks somewhat unexplored, particularly for non-model organisms with complicated (very long, heterozygous, polyploid) genomes. To disclose the performance of various transcriptome assembly programs, this study built 11 single assemblies and analyzed their performance on some significant reference-free and reference-based criteria. As well as to reconfirm the outputs of benchmarks, 55 BLAST were performed and compared using 11 constructed transcriptomes. Concisely, normalized benchmarking demonstrated that Velvet–Oases suffer from the worst results, while the EvidentialGene strategy can provide the most comprehensive and accurate transcriptome of Lilium ledebourii (Baker) Boiss. The BLAST results also confirmed the superiority of EvidentialGene, so it could capture even up to 59% more (than Velvet–Oases) unique gene hits. To promote assembly optimization, with the help of normalized benchmarking, PCA and AHC, it is emphasized that each metric can only provide part of the transcriptome status, and one should never settle for just a few evaluation criteria. This study supplies a framework for benchmarking and optimizing the efficiency of assembly approaches to analyze RNA-Seq data and reveals that selecting an inefficient assembly strategy might result in less identification of unique gene hits.

Salivary and Intestinal Transcriptomes Reveal Differential Gene Expression in Starving, Fed and Trypanosoma cruzi-Infected Rhodnius neglectus

Rhodnius neglectus is a potential vector of Trypanosoma cruzi (Tc), the causative agent of Chagas disease. The salivary glands (SGs) and intestine (INT) are actively required during blood feeding. The saliva from SGs is injected into the vertebrate host, modulating immune responses and favoring feeding for INT digestion. Tc infection significantly alters the physiology of these tissues; however, studies that assess this are still scarce. This study aimed to gain a better understanding of the global transcriptional expression of genes in SGs and INT during fasting (FA), fed (FE), and fed in the presence of Tc (FE + Tc) conditions. In FA, the expression of transcripts related to homeostasis maintenance proteins during periods of stress was predominant. Therefore, the transcript levels of Tret1-like and Hsp70Ba proteins were increased. Blood appeared to be responsible for alterations found in the FE group, as most of the expressed transcripts, such as proteases and cathepsin D, were related to digestion. In FE + Tc group, there was a decreased expression of blood processing genes for insect metabolism (e.g., Antigen-5 precursor, Pr13a, and Obp), detoxification (Sult1) in INT and acid phosphatases in SG. We also found decreased transcriptional expression of lipocalins and nitrophorins in SG and two new proteins, pacifastin and diptericin, in INT. Several transcripts of unknown proteins with investigative potential were found in both tissues. Our results also show that the presence of Tc can change the expression in both tissues for a long or short period of time. While SG homeostasis seems to be re-established on day 9, changes in INT are still evident. The findings of this study may be used for future research on parasite-vector interactions and contribute to the understanding of food physiology and post-meal/infection in triatomines.

Transcriptional expression changes during compensatory plasticity in the terminal ganglion of the adult cricket Gryllus bimaculatus

BACKGROUND

Damage to the adult central nervous system often leads to long-term disruptions in function due to the limited capacity for neurological recovery. The central nervous system of the Mediterranean field cricket, Gryllus bimaculatus, shows an unusual capacity for compensatory plasticity, most obviously in the auditory system and the cercal escape system. In both systems, unilateral sensory disruption leads the central circuitry to compensate by forming and/or strengthening connections with the contralateral sensory organ. While this compensatory plasticity in the auditory system relies on robust dendritic sprouting and novel synapse formation, the compensatory plasticity in the cercal escape circuitry shows little obvious dendritic sprouting and instead may rely on shifts in excitatory and inhibitory synaptic strength.

RESULTS

In order to better understand what types of molecular pathways might underlie this compensatory shift in the cercal system, we used a multiple k-mer approach to assemble a terminal ganglion transcriptome that included ganglia collected one, three, and 7 days after unilateral cercal ablation in adult, male animals. We performed differential expression analysis using EdgeR and DESeq2 and examined Gene Ontologies to identify candidates potentially involved in this plasticity. Enriched GO terms included those related to the ubiquitin-proteosome protein degradation system, chromatin-mediated transcriptional pathways, and the GTPase-related signaling system.

CONCLUSION

Further exploration of these GO terms will provide a clearer picture of the processes involved in compensatory recovery of the cercal escape system in the cricket and can be compared and contrasted with the distinct pathways that have been identified upon deafferentation of the auditory system in this same animal.

Transcriptome Signature of Immune Cells Post Reovirus Treatment in KRAS Mutated Colorectal Cancer

Purpose

Reovirus propagates with high efficiency in KRAS mutated colorectal cancer (CRC). About 45–50% of CRC patients possess a KRAS mutation. Oncolytic reovirus treatment in combination with chemotherapy was tested in patients possessing KRAS mutated metastatic CRC. This study evaluates the biological responses to reovirus treatment by determining the gene expression patterns in RAS-related signaling pathways.

Methods

Reovirus was administered as a 60-min intravenous infusion for 5 consecutive days every 28 days, at a tissue culture infective dose (TCID₅₀) of 3×10¹⁰. Peripheral blood mononuclear cells (PBMCs) were isolated from whole-blood pre- and post-reovirus administration at 48 hr, day-8, and day-15. Clariom_D_Human_Assay was used to determine the expression of vital genes compared to pre-reovirus treatment by RNA sequencing. Using exported sample signals, ΔΔCt method was used to analyze the fold changes of genes within seven gene pathways. Significance was calculated by students-two-tail-t-test. Hierarchical clustering dendrogram was constructed by calculating Pearson’s correlation coefficients.

Results

As compared to the control, SOS1[48 hr; 2.49X], RRAS [48 hr; 2.24X], PIK3CB [D8, D15; 2.27X, 3.16X], MIR 16–2 [D15; 1.70X], CHORDC1 [48 hr, D15; 1.89X, 4.54X], RTN4 [48 hr; 4.66X], FAM96A [48 hr; 4.54X], NFKB [D8, D15; 19.0X, 1.42X], CASP8 [D8, D15; 2.11X, 1.77X], and CASP9 [D8; 1.45X] are upregulated post-reovirus. NOS3 [D15; 0.61X], SYNE1 [D8, D15; 0.78X, 0.71X], ANGPT1 [D8; 0.62X], VEGFB [48 hr, D8, D15; 0.44X, 0.28X, 0.28X], JUN [D15; 0.69X], and IGF2 [D8; 0.73X] are downregulated post-reovirus. Fold change values were significant [p<0.05].

Conclusion

This study highlights reovirus as a novel treatment option for KRAS mutated CRC and showcases its effect on the expression of crucial genes.

Transcriptome analysis of hepatopancreas from different living states oriental river prawn (Macrobrachium nipponense) in response to hypoxia

As an important economical freshwater prawn, Macrobrachium nipponense has difficulty with adapting to hypoxia. In this study, comparative transcriptome analysis was used for the first time to explore the differences between different living states of Macrobrachium nipponense under hypoxia. A total of 94.22 Gb clean reads were obtained and assembled into 54,688 unigenes. A total of 224, 266, and 750 differently expressed genes were found in the comparison of the control and death groups, the control and moribund groups, and the control and survived groups, respectively. Three signal pathways closely related to hypoxia were found by enriching of the signal pathways in three comparison groups. In addition, much attention was focused on the differential genes in these pathways. Oxidative stress related genes, such as 70 kDa heat shock protein, phosphoenolpyruvate carboxykinase and cyclooxygenase were differentially expressed in different comparisons. After comparing with previous studies, cyclooxygenase was found to be an important hypoxia-related gene that is fully involved in the hypoxic response. Interestingly, two new genes with no Nr annotation were found in this manuscript. This manuscript will enrich our understanding of oxidative stress response to hypoxia and provide a theoretical basis for the subsequent solution of apoptosis caused by hypoxia.

Constructing a de novo transcriptome and a reference proteome for the bivalve Scrobicularia plana: Comparative analysis of different assembly strategies and proteomic analysis

Scrobicularia plana is a coastal and estuarine bivalve widely used in ecotoxicological studies. However, the underlying molecular mechanisms for S. plana pollutant responses are hardly known due to the lack of molecular databases. Thus, in this study we present a holistic approach to assess a robust reference transcriptome and proteome of this clam. A mixture of control and metal-exposed individuals was used for mRNA isolation. Four sets of high quality filtered preprocessed reads were generated (two quality scores and two sequenced lengths) and assembled with Mira, Ray and Trinity algorithms. The sixty-four generated assemblies were refined, filtered and evaluated for their proteomic quality. Eight assemblies presented top Detonate scores but one was selected due to its compactness and biological representation, which was generated: (i) from the highest quality dataset (Q20L100), (ii) using Trinity algorithm with all k-mers (AtKa), (iii) removing redundancy by CD-HIT (RR80), and (iv) filtering out poor contigs (F), that was subsequently named Q20L100AtKaRR80F. S. plana proteomic analysis revealed 10,017 peptide groups that corresponded to 2066 proteins with a wide coverage of molecular functions and biological processes, confirming the strength of the database generated.

How to Become an Apomixis Model: The Multifaceted Case of Paspalum

In the past decades, the grasses of the Paspalum genus have emerged as a versatile model allowing evolutionary, genetic, molecular, and developmental studies on apomixis as well as successful breeding applications. The rise of such an archetypal system progressed through integrative phases, which were essential to draw conclusions based on solid standards. Here, we review the steps adopted in Paspalum to establish the current body of knowledge on apomixis and provide model breeding programs for other agronomically important apomictic crops. In particular, we discuss the need for previous detailed cytoembryological and cytogenetic germplasm characterization; the establishment of sexual and apomictic materials of identical ploidy level; the development of segregating populations useful for inheritance analysis, positional mapping, and epigenetic control studies; the development of omics data resources; the identification of key molecular pathways via comparative gene expression studies; the accurate molecular characterization of genomic loci governing apomixis; the in-depth functional analysis of selected candidate genes in apomictic and model species; the successful building of a sexual/apomictic combined breeding scheme.

BHap: a novel approach for bacterial haplotype reconstruction

MOTIVATION

The bacterial haplotype reconstruction is critical for selecting proper treatments for diseases caused by unknown haplotypes. Existing methods and tools do not work well on this task, because they are usually developed for viral instead of bacterial populations.

RESULTS

In this study, we developed BHap, a novel algorithm based on fuzzy flow networks, for reconstructing bacterial haplotypes from next generation sequencing data. Tested on simulated and experimental datasets, we showed that BHap was capable of reconstructing haplotypes of bacterial populations with an average F1 score of 0.87, an average precision of 0.87 and an average recall of 0.88. We also demonstrated that BHap had a low susceptibility to sequencing errors, was capable of reconstructing haplotypes with low coverage and could handle a wide range of mutation rates. Compared with existing approaches, BHap outperformed them in terms of higher F1 scores, better precision, better recall and more accurate estimation of the number of haplotypes.

AVAILABILITY AND IMPLEMENTATION

The BHap tool is available at http://www.cs.ucf.edu/∼xiaoman/BHap/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Combining independent de novo assemblies to optimize leaf transcriptome of Persian walnut

Transcriptome resources can facilitate to increase yield and quality of walnuts. Finding the best transcriptome assembly has not been the subject of walnuts research as yet. This research generated 240,179,782 reads from 11 walnut leaves according to cDNA libraries. The reads provided a complete de novo transcriptome assembly. Fifteen different transcriptome assemblies were constructed from five different well-known assemblers used in scientific literature with different k-mer lengths (Bridger, BinPacker, SOAPdenovo-Trans, Trinity and SPAdes) as well as two merging approaches (EvidentialGene and Transfuse). Based on the four quality metrics of assembly, the results indicated an efficiency in the process of merging the assemblies after being generated by de novo assemblers. Finally, EvidentialGene was recognized as the best assembler for the de novo assembly of the leaf transcriptome in walnut. Among a total number of 183,191 transcripts which were generated by EvidentialGene, there were 109,413 transcripts capable of protein potential (59.72%) and 104,926 were recognized as ORFs (57.27%). In addition, 79,185 transcripts were predicted to exist with at least one hit to the Pfam database. A number of 3,931 transcription factors were identified by BLAST searching against PlnTFDB. Furthermore, 6,591 of the predicted peptide sequences contained signaling peptides, while 92,704 contained transmembrane domains. Comparison of the assembled transcripts with transcripts of the walnut and published genome assembly for the 'Chandler' cultivar using the BLAST algorithm led to identify a total number of 27,304 and 19,178 homologue transcripts, respectively. De novo transcriptomes in walnut leaves can be developed for the future studies in functional genomics and genetic studies of walnuts.

Effect of vitamin D3 on the immunomodulation of head kidney after Edwardsiella ictaluri challenge in yellow catfish (Pelteobagrus fulvidraco)

Edwardsiella ictaluri (E. ictaluri) causes severe infections in yellow catfish (Pelteobagrus fulvidraco), which leads to a massive loss in the aquaculture industry especially in catfish commercial production. Previous studies have confirmed that vitamin D₃ is essential in immune regulation in mammals. Based on next-generation sequencing, this study explored the immunomodulatory effects of dietary vitamin D₃ on the head kidney of yellow catfish after E. ictaluri challenge. Current results showed that increasing the content of dietary vitamin D₃ within the experimental concentration range (1120IU/kg-16600IU/kg) could reduce the mortality of the yellow catfish after E. ictaluri challenge. Results of the next-generation sequencing showed that dietary vitamin D₃ regulates the immune mechanism of the head kidney mainly through three pathways i.e. negative regulation of interferon-β production, negative regulation of interleukin-6 production and neutrophil chemotaxis. Proteins HSPA8, MAP4K4 and MRC1 may be involved in vitamin D₃-mediated immunoregulation in the head kidney. qPCR results showed that increasing the content of dietary vitamin D₃ can improve the immune function of the yellow catfish by down-regulating ifn-β and pro-inflammatory factors tnf-α, il1-β, il-6, il-8 and up-regulating the anti-inflammatory factor il-10. The above results indicated that dietary addition of vitamin D₃ regulated the immune response in head kidney of yellow catfish and helped the fish to resist the negative effects of infection by E. ictaluri in a dose-dependent manner.

De novo transcriptome analysis of halotolerant bacterium Staphylococcus sp. strain P-TSB-70 isolated from East coast of India: In search of salt stress tolerant genes

In the present study, we identified salt stress tolerant genes from the marine bacterium Staphylococcus sp. strain P-TSB-70 through transcriptome sequencing. In favour of whole-genome transcriptome profiling of Staphylococcus sp. strain P-TSB-70 (GenBank Accn. No. KP117091) which tolerated upto 20% NaCl stress, the strain was cultured in the laboratory condition with 20% NaCl stress. Transcriptome analyses were performed by SOLiD4.0 sequencing technology from which 10280 and 9612 transcripts for control and treated, respectively, were obtained. The coverage per base (CPB) statistics were analyzed for both the samples. Gene ontology (GO) analysis has been categorized at varied graph levels based on three primary ontology studies viz. cellular components, biological processes, and molecular functions. The KEGG analysis of the assembled transcripts using KAAS showed presumed components of metabolic pathways which perhaps implicated in diverse metabolic pathways responsible for salt tolerance viz. glycolysis/gluconeogenesis, oxidative phosphorylation, glutathione metabolism, etc. further involving in salt tolerance. Overall, 90 salt stress tolerant genes were identified as of 186 salt-related transcripts. Several genes have been found executing normally in the TCA cycle pathway, integral membrane proteins, generation of the osmoprotectants, enzymatic pathway associated with salt tolerance. Recognized genes fit diverse groups of salt stress genes viz. abc transporter, betaine, sodium antiporter, sodium symporter, trehalose, ectoine, and choline, that belong to different families of genes involved in the pathway of salt stress. The control sample of the bacterium showed elevated high proportion of transcript contigs (29%) while upto 20% salt stress treated sample of the bacterium showed a higher percentage of transcript contigs (31.28%). A total of 1,288 and 1,133 transcript contigs were measured entirely as novel transcript contigs in both control and treated samples, respectively. The structure and function of 10 significant salt stress tolerant genes of Staphylococcus sp. have been analyzed in this study. The information acquired in the present study possibly used to recognize and clone the salt stress tolerant genes and support in developing the salt stress-tolerant plant varieties to expand the agricultural productivity in the saline system.

Analysis of Metabolites and Gene Expression Changes Relative to Apricot (Prunus armeniaca L.) Fruit Quality During Development and Ripening

Apricot (Prunus armeniaca L.) is a valuable worldwide agronomical crop, with a delicious fruit highlighted as a functional food with both nutritional and bioactive properties, remarkably beneficial to human health. Apricot fruit ripening is a coordinated developmental process which requires change in the expression of hundreds to thousands of genes to modify many biochemical and physiological processes arising from quality characteristics in ripe fruit. In addition, enhancing fruit and nutraceutical quality is one of the central objectives to be improved in the new varieties developed by breeding programs. In this study we analyzed the contents of main metabolites linked to the nutraceutical value of apricot fruits, together with the most important pomological characteristics and biochemical contents of fruit during the ripening process in two contrasted apricot genotypes. Additionally, the gene expression changes were analyzed using RNA-Seq and real time qPCR. Results showed that genes with differential expression in the biosynthetic pathways, such as phenylpropanoids, flavonoids, starch and sucrose and carotenoid metabolism, could be possible candidates as molecular markers of fruit quality characteristics for fruit color and soluble solid content. The gene involves in carotenoid metabolism carotenoid cleavage dioxygenase 4, and the gene sucrose synthase in starch and sucrose metabolism were identified as candidate genes in the ripening process for white skin ground color and flesh color and high soluble sugar content. The application of these candidate genes on marker-assisted selection in apricot breeding programs may contribute to the early selection of high-quality fruit genotypes with suitable nutraceutical values.

Microarray-Based Quality Assessment as a Supporting Criterion for de novo Transcriptome Assembly Selection

RNA-Sequencing and de novo assembly have enabled the analysis of species with non-available reference transcriptomes, although intrinsic features (biological and technical) induce errors in the reconstruction. A strategy to resolve these errors consists of varying assembling process parameters to generate multiple reconstructions. However, the best assembly selection remains a challenge. Quantitative metrics for quality assessment have been inconsistent when compared with pertinent references. In this paper, a criterion for supporting assembly selection based on mapping DNA microarray hybridized probes to assembly sets is proposed. Mouse and fruit fly RNA-Seq datasets were assembled with standard de novo procedures. Quality assessment was estimated using quantitative metrics and the proposed criterion. The assembly that best mapped to the available reference transcriptomes of these model species provided the highest quality assembly. The hybridized probes identified the best assemblies, whereas quantitative metrics remained inconsistent. For example, subtle probe mapping difference of 0.25 percent, but statistically significant (ANOVA, p < 0.05), enabled the assembly selection that led to identify 3,719 more contigs and led to 1,049 further mapped contigs to the mouse reference transcriptome. The microarray data availability for non-model species makes the proposed criterion suitable for quality assessment of multiple de novo assembly strategies.

Comparative Analysis of Strategies for De Novo Transcriptome Assembly in Prokaryotes: Streptomyces clavuligerus as a Case Study

The performance of software tools for de novo transcriptome assembly greatly depends on the selection of software parameters. Up to now, the development of de novo transcriptome assembly for prokaryotes has not been as remarkable as that for eukaryotes. In this contribution, Rockhopper2 was used to perform a comparative transcriptome analysis of Streptomyces clavuligerus exposed to diverse environmental conditions. The study focused on assessing the incidence of software parameters on software performance for the identification of differentially expressed genes as a final goal. For this, a statistical optimization was performed using the Transrate Assembly Score (TAS). TAS was also used for evaluating the software performance and for comparing it with related tools, e.g., Trinity. Transcriptome redundancy and completeness were also considered for this analysis. Rockhopper2 and Trinity reached a TAS value of 0.55092 and 0.58337, respectively. Trinity assembles transcriptomes with high redundancy, with 55.6% of transcripts having some duplicates. Additionally, we observed that the total number of differentially expressed genes (DEG) and their annotation greatly depends on the method used for removing redundancy and the tools used for transcript quantification. To our knowledge, this is the first work aimed at assessing de novo assembly software for prokaryotic organisms.

Utilization of Tissue Ploidy Level Variation in de Novo Transcriptome Assembly of Pinus sylvestris

Compared to angiosperms, gymnosperms lag behind in the availability of assembled and annotated genomes. Most genomic analyses in gymnosperms, especially conifer tree species, rely on the use of de novo assembled transcriptomes. However, the level of allelic redundancy and transcript fragmentation in these assembled transcriptomes, and their effect on downstream applications have not been fully investigated. Here, we assessed three assembly strategies for short-reads data, including the utility of haploid megagametophyte tissue during de novo assembly as single-allele guides, for six individuals and five different tissues in Pinus sylvestris. We then contrasted haploid and diploid tissue genotype calls obtained from the assembled transcriptomes to evaluate the extent of paralog mapping. The use of the haploid tissue during assembly increased its completeness without reducing the number of assembled transcripts. Our results suggest that current strategies that rely on available genomic resources as guidance to minimize allelic redundancy are less effective than the application of strategies that cluster redundant assembled transcripts. The strategy yielding the lowest levels of allelic redundancy among the assembled transcriptomes assessed here was the generation of SuperTranscripts with Lace followed by CD-HIT clustering. However, we still observed some levels of heterozygosity (multiple gene fragments per transcript reflecting allelic redundancy) in this assembled transcriptome on the haploid tissue, indicating that further filtering is required before using these assemblies for downstream applications. We discuss the influence of allelic redundancy when these reference transcriptomes are used to select regions for probe design of exome capture baits and for estimation of population genetic diversity.

Transcriptome Sequencing Reveals Regulatory Mechanisms of Taxol Synthesis in Taxus wallichiana var. Mairei

Taxol is one of the most potent and effective anticancer drugs and is originally isolated from Taxus species. To investigate the specific regulatory mechanisms of taxol synthesis in Taxus wallichiana var. mairei, RNA-seq was conducted to reveal the differences in transcriptional levels between wild type (WT) and “Jinxishan” (JXS), a cultivar selected from a population of Taxus mairei that shows about 3-fold higher taxol content in the needles than WT. Our results indicated that high expressions of the genes taxadienol acetyltransferase (TAT), taxadiene 5-alpha hydroxylase (T5H), 5-alpha-taxadienol-10-beta-hydroxylase (T10OH), and 2-debenzoyl-7,13-diacetylbaccatin III-2-O-benzoyl-transferase (DBBT), which catalyze a series of key acetylation and hydroxylation steps, are the main cause of high taxol content in JXS. Moreover, in the present study, the activation of jasmonic acid (JA) signal transduction and its crosstalk with gibberellin (GA), auxin, and ethylene (ET) explained the elevation of differentially expressed genes (DEGs) from the taxol biosynthesis pathway. This also indicates that taxol biosynthesis in T. mairei is associated with the balance of cell development and defense. TF-encoding (transcriptional factor) genes, represented by the ethylene-responsive transcription factor (ERF), basic/helix-loop-helix (bHLH), MYB, and WRKY families, were detected as differentially expressed between JXS and WT, further indicating that the regulation of hormone signaling on taxol biosynthesis genes was mediated by transcription factors (TFs). To our knowledge, this is the first study to illustrate the regulatory mechanisms of taxol synthesis in a new cultivar of T. mairei with a high taxol content in its needles. These transcriptome data provide reasonable explanations for the variation of taxol content between WT and JXS.

Biosynthesis and Signal Transduction of ABA, JA, and BRs in Response to Drought Stress of Kentucky Bluegrass

Kentucky bluegrass (KB, Poa pratensis) is one of the most widely used cool-season turfgrass species, but it is sensitive to drought stress. Molecular studies in KB are hindered by its large and complex genome structure. In this study, a comparative transcriptomic study was conducted between a short and long period of water deficiency. Three transcriptome libraries were constructed and then sequenced by using leaf RNA samples of plants at 0, 2, and 16 h after PEG6000 treatment. A total of 199,083 differentially expressed genes (DEGs) were found. The Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation revealed that DEGs were enriched in "Plant hormone signal transduction" and "MAPK signaling pathway-Plant". Some key up-regulated genes, including PYL, JAZ, and BSK, were involved in hormone signaling transduction of abscisic acid, jasmonic acid, and brassinosteroid and possibly these genes play important roles in coping with drought stress in KB. Furthermore, our results showed that the concentrations of ABA, JA and BR increased significantly with the extension of the drought period. The specific DEGs encoding functional proteins, kinase and transcription factors, could be valuable information for genetic manipulation to promote drought tolerance of KB in the future.

Biosynthesis and Signal Transduction of ABA, JA, and BRs in Response to Drought Stress of Kentucky Bluegrass

Kentucky bluegrass (KB, Poa pratensis) is one of the most widely used cool-season turfgrass species, but it is sensitive to drought stress. Molecular studies in KB are hindered by its large and complex genome structure. In this study, a comparative transcriptomic study was conducted between a short and long period of water deficiency. Three transcriptome libraries were constructed and then sequenced by using leaf RNA samples of plants at 0, 2, and 16 h after PEG6000 treatment. A total of 199,083 differentially expressed genes (DEGs) were found. The Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation revealed that DEGs were enriched in "Plant hormone signal transduction" and "MAPK signaling pathway-Plant". Some key up-regulated genes, including PYL, JAZ, and BSK, were involved in hormone signaling transduction of abscisic acid, jasmonic acid, and brassinosteroid and possibly these genes play important roles in coping with drought stress in KB. Furthermore, our results showed that the concentrations of ABA, JA and BR increased significantly with the extension of the drought period. The specific DEGs encoding functional proteins, kinase and transcription factors, could be valuable information for genetic manipulation to promote drought tolerance of KB in the future.

A high-throughput SNP discovery strategy for RNA-seq data

Background

Single nucleotide polymorphisms (SNP) have been applied as important molecular markers in genetics and breeding studies. The rapid advance of next generation sequencing (NGS) provides a high-throughput means of SNP discovery. However, SNP development is limited by the availability of reliable SNP discovery methods. Especially, the optimum assembler and SNP caller for accurate SNP prediction from next generation sequencing data are not known.

Results

Herein we performed SNP prediction based on RNA-seq data of peach and mandarin peel tissue under a comprehensive comparison of two paired-end read lengths (125 bp and 150 bp), five assemblers (Trinity, IDBA, oases, SOAPdenovo, Trans-abyss) and two SNP callers (GATK and GBS). The predicted SNPs were compared with the authentic SNPs identified via PCR amplification followed by gene cloning and sequencing procedures. A total of 40 and 240 authentic SNPs were presented in five anthocyanin biosynthesis related genes in peach and in nine carotenogenic genes in mandarin. Putative SNPs predicted from the same RNA-seq data with different strategies led to quite divergent results. The rate of false positive SNPs was significantly lower when the paired-end read length was 150 bp compared with 125 bp. Trinity was superior to the other four assemblers and GATK was substantially superior to GBS due to a low rate of missing authentic SNPs. The combination of assembler Trinity, SNP caller GATK, and the paired-end read length 150 bp had the best performance in SNP discovery with 100% accuracy both in peach and in mandarin cases. This strategy was applied to the characterization of SNPs in peach and mandarin transcriptomes.

Conclusions

Through comparison of authentic SNPs obtained by PCR cloning strategy and putative SNPs predicted from different combinations of five assemblers, two SNP callers, and two paired-end read lengths, we provided a reliable and efficient strategy, Trinity-GATK with 150 bp paired-end read length, for SNP discovery from RNA-seq data. This strategy discovered SNP at 100% accuracy in peach and mandarin cases and might be applicable to a wide range of plants and other organisms.

Electronic supplementary material

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

Usability of reference-free transcriptome assemblies for detection of differential expression: a case study on Aethionema arabicum dimorphic seeds

Background

RNA-sequencing analysis is increasingly utilized to study gene expression in non-model organisms without sequenced genomes. Aethionema arabicum (Brassicaceae) exhibits seed dimorphism as a bet-hedging strategy – producing both a less dormant mucilaginous (M⁺) seed morph and a more dormant non-mucilaginous (NM) seed morph. Here, we compared de novo and reference-genome based transcriptome assemblies to investigate Ae. arabicum seed dimorphism and to evaluate the reference-free versus -dependent approach for identifying differentially expressed genes (DEGs).

Results

A de novo transcriptome assembly was generated using sequences from M⁺ and NM Ae. arabicum dry seed morphs. The transcripts of the de novo assembly contained 63.1% complete Benchmarking Universal Single-Copy Orthologs (BUSCO) compared to 90.9% for the transcripts of the reference genome. DEG detection used the strict consensus of three methods (DESeq2, edgeR and NOISeq). Only 37% of 1533 differentially expressed de novo assembled transcripts paired with 1876 genome-derived DEGs. Gene Ontology (GO) terms distinguished the seed morphs: the terms translation and nucleosome assembly were overrepresented in DEGs higher in abundance in M⁺ dry seeds, whereas terms related to mRNA processing and transcription were overrepresented in DEGs higher in abundance in NM dry seeds. DEGs amongst these GO terms included ribosomal proteins and histones (higher in M⁺), RNA polymerase II subunits and related transcription and elongation factors (higher in NM). Expression of the inferred DEGs and other genes associated with seed maturation (e.g. those encoding late embryogenesis abundant proteins and transcription factors regulating seed development and maturation such as ABI3, FUS3, LEC1 and WRI1 homologs) were put in context with Arabidopsis thaliana seed maturation and indicated that M⁺ seeds may desiccate and mature faster than NM. The 1901 transcriptomic DEG set GO-terms had almost 90% overlap with the 2191 genome-derived DEG GO-terms.

Conclusions

Whilst there was only modest overlap of DEGs identified in reference-free versus -dependent approaches, the resulting GO analysis was concordant in both approaches. The identified differences in dry seed transcriptomes suggest mechanisms underpinning previously identified contrasts between morphology and germination behaviour of M⁺ and NM seeds.

Electronic supplementary material

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

Usability of reference-free transcriptome assemblies for detection of differential expression: a case study on Aethionema arabicum dimorphic seeds

BACKGROUND

RNA-sequencing analysis is increasingly utilized to study gene expression in non-model organisms without sequenced genomes. Aethionema arabicum (Brassicaceae) exhibits seed dimorphism as a bet-hedging strategy - producing both a less dormant mucilaginous (M⁺) seed morph and a more dormant non-mucilaginous (NM) seed morph. Here, we compared de novo and reference-genome based transcriptome assemblies to investigate Ae. arabicum seed dimorphism and to evaluate the reference-free versus -dependent approach for identifying differentially expressed genes (DEGs).

RESULTS

A de novo transcriptome assembly was generated using sequences from M⁺ and NM Ae. arabicum dry seed morphs. The transcripts of the de novo assembly contained 63.1% complete Benchmarking Universal Single-Copy Orthologs (BUSCO) compared to 90.9% for the transcripts of the reference genome. DEG detection used the strict consensus of three methods (DESeq2, edgeR and NOISeq). Only 37% of 1533 differentially expressed de novo assembled transcripts paired with 1876 genome-derived DEGs. Gene Ontology (GO) terms distinguished the seed morphs: the terms translation and nucleosome assembly were overrepresented in DEGs higher in abundance in M⁺ dry seeds, whereas terms related to mRNA processing and transcription were overrepresented in DEGs higher in abundance in NM dry seeds. DEGs amongst these GO terms included ribosomal proteins and histones (higher in M⁺), RNA polymerase II subunits and related transcription and elongation factors (higher in NM). Expression of the inferred DEGs and other genes associated with seed maturation (e.g. those encoding late embryogenesis abundant proteins and transcription factors regulating seed development and maturation such as ABI3, FUS3, LEC1 and WRI1 homologs) were put in context with Arabidopsis thaliana seed maturation and indicated that M⁺ seeds may desiccate and mature faster than NM. The 1901 transcriptomic DEG set GO-terms had almost 90% overlap with the 2191 genome-derived DEG GO-terms.

CONCLUSIONS

Whilst there was only modest overlap of DEGs identified in reference-free versus -dependent approaches, the resulting GO analysis was concordant in both approaches. The identified differences in dry seed transcriptomes suggest mechanisms underpinning previously identified contrasts between morphology and germination behaviour of M⁺ and NM seeds.

Sexual dimorphism and heightened conditional expression in a sexually selected weapon in the Asian rhinoceros beetle

Among the most dramatic examples of sexual selection are the weapons used in battles between rival males over access to females. As with ornaments of female choice, the most "exaggerated" sexually selected weapons vary from male to male more widely than other body parts (hypervariability), and their growth tends to be more sensitive to nutritional state or physiological condition compared with growth of other body parts ("heightened" conditional expression). Here, we use RNAseq analysis to build on recent work exploring these mechanisms in the exaggerated weapons of beetles, by examining patterns of differential gene expression in exaggerated (head and thorax horns) and non-exaggerated (wings, genitalia) traits in the Asian rhinoceros beetle, Trypoxylus dichotomus. Our results suggest that sexually dimorphic expression of weaponry involves large-scale changes in gene expression, relative to other traits, while nutrition-driven changes in gene expression in these same weapons are less pronounced. However, although fewer genes overall were differentially expressed in high- vs. low-nutrition individuals, the number of differentially expressed genes varied predictably according to a trait's degree of condition dependence (head horn > thorax horn > wings > genitalia). Finally, we observed a high degree of similarity in direction of effects (vectors) for subsets of differentially expressed genes across both sexually dimorphic and nutritionally responsive growth. Our results are consistent with a common set of mechanisms governing sexual size dimorphism and condition dependence.

Characterization of Laguncularia racemosa transcriptome and molecular response to oil pollution

Mangroves are ecosystems of economic and ecological importance. Laguncularia racemosa (Combretaceae), popularly known as white mangrove, is a species that greatly contributes to the community structure of neotropical and West African mangrove forests. Despite the significance of these ecosystems, they have been destroyed by oil spills that can cause yellowing of leaves, increased sensitivity to other stresses and death of trees. However, the molecular response of plants to oil stress is poorly known. In this work, Illumina reads were de novo assembled into 46,944 transcripts of L. racemosa roots and leaves, including putative isoform variants. In addition to improving the genomic information available for mangroves, the L. racemosa assembled transcriptome allowed us to identify reference genes to normalize quantitative real-time PCR (qPCR) expression data from oil-stressed mangrove plants, which were used in RNASeq validation. The analysis of expression changes induced by the oil exposure revealed 310 and 286 responsive transcripts of leaves and roots, respectively, mainly up-regulated. Enriched GO categories related to chloroplasts and photosynthesis were found among both leaf and root oil-responsive transcripts, while "response to heat" and "response to hypoxia" were exclusively enriched in leaves and roots, respectively. The comparison of L. racemosa 12-h-oil-stressed leaf expression profile to previous Arabidopsis heat-stress studies and co-expression evidence also pointed to similarities between the heat and oil responses, in which the HSP-coding genes seem to play a key role. A subset of the L. racemosa oil-responsive root genes exhibited similar up-regulation profiles to their Arabidopsis homologs involved in hypoxia responses, including the HRA1 and LBD41 TF-coding genes. Genes linked to the ethylene pathway such as those coding for ERF TFs were also modulated during the L. racemosa root response to oil stress. Taken together, these results show that oil contamination affects photosynthesis, protein metabolism, hypoxia response and the ethylene pathway in L. racemosa 12-h-oil-exposed leaves and roots.

Transcriptional profiling reveals that a MYB transcription factor MsMYB4 contributes to the salinity stress response of alfalfa

MYB transcription factors are important regulators of the plant response to abiotic stress. Their participation in the salinity stress of the key forage legume species alfalfa (Medicago sativa) was investigated here by comparing the transcriptomes of the two cultivars Dryland (DL) and Sundory (SD), which differed with respect to their ability to tolerate salinity stress. When challenged by the stress, DL plants were better able than SD ones to scavenge reactive oxygen species. A large number of genes encoding transcription regulators, signal transducers and proteins involved in both primary and secondary metabolism were differentially transcribed in the two cultivars, especially when plants were subjected to salinity stress. The set of induced genes included 17 MYB family of transcription factors, all of which were subsequently isolated. The effect of constitutively expressing these genes on the salinity tolerance expressed by Arabidopsis thaliana was investigated. The introduction of MsMYB4 significantly increased the plants’ salinity tolerance in an abscisic acid-dependent manner. A sub-cellular localization experiment and a transactivation assay indicated that MsMYB4 was deposited in the nucleus and was able to activate transcription in yeast. Based on this information, we propose that the MsMYB4 products is likely directly involved in alfalfa’s response to salinity stress.

Overexpression of Lilium formosanumMADS-box (LFMADS) Causing Floral Defects While Promoting Flowering in Arabidopsis thaliana, Whereas Only Affecting Floral Transition Time in Nicotiana tabacum

The Formosa lily (Lilium formosanum) is one of the most common horticultural species in Taiwan. To explore gene regulation involved in this species, we used transcriptome analysis to generate PH-FB (mixed floral buds) and PH-LF (mature leaves) datasets. Combination of the PH-FB and PH-LF constructed a de novo assembly of the ALL dataset, including 18,041 contigs and 23,807 unigenes by Nr, GO, COG, and KEGG databases. The differential gene expression (DGE) analysis revealed 9937 genes were upregulated while 10,383 genes were downregulated in the developing floral buds compared to mature leaves. Seven putative genes (LFMADS1 to 7) encoding floral organ identity proteins were selected for further analysis. LFMADS1-6 genes were specifically expressed in the floral organ, while LFMADS7 in the floral buds and mature leaves. Phylogenetic analysis revealed that LFMADS1-3 is classified into B-class, LFMADS4 into C-class, LFMADS5 into D-class, and LFMADS6-7 into E-class, respectively. LFMADS-GFP fusion proteins appeared to localize in the nucleus, supporting their roles as transcription factors (TFs). Overexpression of the LFMADS2, LFMADS4, and LFMADS6 genes in Arabidopsis resulted in early flowering and floral defect, however, only early flowering in transgenic tobacco was observed. Highly expressed floral integrator genes, including AtFT, AtLFY, and AtFUL in transgenic Arabidopsis and NtFUL and NtSOC1 in transgenic tobacco, resulted in early flowering phenotype through qRT-PCR analysis. Yeast two-hybrid analysis suggested that LFMADSs may form higher order complexes with the B-, C-, D, and/or E-class proteins to determine the floral organ identity. Furthermore, E-class LFMADS proteins may function as a glue to mediate and strengthen the protein-protein interactions. Therefore, our de novo datasets would provide information for investigating other differentially expressed candidate transcripts. In addition, functional conservation of LFMADSs appears to be vital in floral transition and floral organ identity.

A comprehensive analysis of the Korean fir (Abies koreana) genes expressed under heat stress using transcriptome analysis

Korean fir (Abies koreana), a rare species endemic to South Korea, is sensitive to climate change. Here, we used next-generation massively parallel sequencing technology and de novo transcriptome assembly to gain a comprehensive overview of the Korean fir transcriptome under heat stress. Sequencing control and heat-treated samples of Korean fir, we obtained more than 194,872,650 clean reads from each sample. After de novo assembly and quantitative assessment, 42,056 unigenes were generated with an average length of 908 bp. In total, 6,401 differentially expressed genes were detected, of which 2,958 were up-regulated and 3,443 down-regulated, between the heat-treated and control samples. A gene ontology analysis of these unigenes revealed heat-stress-related terms, such as "response to stimulus". Further, in depth analysis revealed 204 transcription factors and 189 Hsps as differentially expressed. Finally, 12 regulated candidate genes associated with heat stress were examined using quantitative real-time PCR (qRT-PCR). In this study, we present the first comprehensive characterisation of Korean fir subjected to heat stress using transcriptome analysis. It provides an important resource for future studies of Korean fir with the objective of identifying heat stress tolerant lines.

IDP-denovo: de novo transcriptome assembly and isoform annotation by hybrid sequencing

Motivation

In the past years, the long read (LR) sequencing technologies, such as Pacific Biosciences and Oxford Nanopore Technologies, have been demonstrated to substantially improve the quality of genome assembly and transcriptome characterization. Compared to the high cost of genome assembly by LR sequencing, it is more affordable to generate LRs for transcriptome characterization. That is, when informative transcriptome LR data are available without a high-quality genome, a method for de novo transcriptome assembly and annotation is of high demand.

Results

Without a reference genome, IDP-denovo performs de novo transcriptome assembly, isoform annotation and quantification by integrating the strengths of LRs and short reads. Using the GM12878 human data as a gold standard, we demonstrated that IDP-denovo had superior sensitivity of transcript assembly and high accuracy of isoform annotation. In addition, IDP-denovo outputs two abundance indices to provide a comprehensive expression profile of genes/isoforms. IDP-denovo represents a robust approach for transcriptome assembly, isoform annotation and quantification for non-model organism studies. Applying IDP-denovo to a non-model organism, Dendrobium officinale, we discovered a number of novel genes and novel isoforms that were not reported by the existing annotation library. These results reveal the high diversity of gene isoforms in D.officinale, which was not reported in the existing annotation library.

Availability and implementation

The dataset of Dendrobium officinale used/analyzed during the current study has been deposited in SRA, with accession code SRP094520. IDP-denovo is available for download at www.healthcare.uiowa.edu/labs/au/IDP-denovo/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Identification of genes associated with stress tolerance in moth bean [Vigna aconitifolia (Jacq.) Marechal], a stress hardy crop

Moth bean is the most drought and heat tolerant cultigens among Asian Vigna. We performed comparative transcriptome analysis of moth bean cultivar "Marumoth" under control and stress condition. De novo transcriptome assembly was carried out by using Velvet followed by Oases softwares. Differential expression analyses, SSR identification and validation and mapping of pathways and transcription factors were conducted. A total of 179,979 and 201,888 reads were generated on Roche 454 platform and 48,617,205 and 45,449,053 reads were generated on ABI Solid platform for the control and stressed samples. Combined assembly from Roche and ABI Solid platforms generated 16,090 and 15,096 transcripts for control and stressed samples. We found 1287 SSRs and 5606 transcripts involved in 179 pathways. The 55 transcription factor families represented 19.42% of total mothbean transcripts. In expression profiling, ten transcripts were found to be up-regulated and 41 down-regulated while 490 showed no major change under moisture stress condition. Stress inducible genes like Catalase, Cyt P450 monooxygenase, heat shock proteins (HSP 90 and HSP 70), oxidoreductase, protein kinases, dehydration responsive protein (DRP), universal stress protein and ferridoxin NADH oxidoreductase genes were up-regulated in stressed sample. Genes which might be involved in moisture stress tolerance in moth bean were identified and these might be useful for stress tolerance breeding in moth bean and other related crops.

Evaluating the Performance of De Novo Assembly Methods for Venom-Gland Transcriptomics

Venom-gland transcriptomics is a key tool in the study of the evolution, ecology, function, and pharmacology of animal venoms. In particular, gene-expression variation and coding sequences gained through transcriptomics provide key information for explaining functional venom variation over both ecological and evolutionary timescales. The accuracy and usefulness of inferences made through transcriptomics, however, is limited by the accuracy of the transcriptome assembly, which is a bioinformatic problem with several possible solutions. Several methods have been employed to assemble venom-gland transcriptomes, with the Trinity assembler being the most commonly applied among them. Although previous evidence of variation in performance among assembly software exists, particularly regarding recovery of difficult-to-assemble multigene families such as snake venom metalloproteinases, much work to date still employs a single assembly method. We evaluated the performance of several commonly used de novo assembly methods for the recovery of both nontoxin transcripts and complete, high-quality venom-gene transcripts across eleven snake and four scorpion transcriptomes. We varied k-mer sizes used by some assemblers to evaluate the impact of k-mer length on transcript recovery. We showed that the recovery of nontoxin transcripts and toxin transcripts is best accomplished through different assembly software, with SDT at smaller k-mer lengths and Trinity being best for nontoxin recovery and a combination of SeqMan NGen and a seed-and-extend approach implemented in Extender as the best means of recovering a complete set of toxin transcripts. In particular, Extender was the only means tested capable of assembling multiple isoforms of the diverse snake venom metalloproteinase family, while traditional approaches such as Trinity recovered at most one metalloproteinase transcript. Our work demonstrated that traditional metrics of assembly performance are not predictive of performance in the recovery of complete and high quality toxin genes. Instead, effective venom-gland transcriptomic studies should combine and quality-filter the results of several assemblers with varying algorithmic strategies.

Optimization of Assembly Pipeline may Improve the Sequence of the Chloroplast Genome in Quercus spinosa

Obtaining chloroplast (cp) genome sequence is necessary for studying physiological roles in plants. However, it is difficult to use traditional sequencing methods to get cp genome sequences because of the complex procedures of preparing templates. With the advent of next-generation sequencing technology, massive genome sequences can be produced. Thus, a good pipeline to assemble next-generation sequence reads with optimized k-mer length is essential to get whole cp genome sequences. Moreover, adjustment of other parameters is also very important, especially for the assembly of the cp genome. In this study, we developed a pipeline to generate the cp genome for Quercus spinosa. When Quercus rubra was used as a reference, we achieved coverage of 97.75% after optimizing k-mer length as well as other parameters. The efficiency of the pipeline makes it a useful method for cp genome construction in plants. It also provides great perspective on the analysis of cp genome characteristics and evolution.

Transcriptome sequencing and annotation of the halophytic microalga Dunaliella salina

The unicellular green alga Dunaliella salina is well adapted to salt stress and contains compounds (including β-carotene and vitamins) with potential commercial value. A large transcriptome database of D. salina during the adjustment, exponential and stationary growth phases was generated using a high throughput sequencing platform. We characterized the metabolic processes in D. salina with a focus on valuable metabolites, with the aim of manipulating D. salina to achieve greater economic value in large-scale production through a bioengineering strategy. Gene expression profiles under salt stress verified using quantitative polymerase chain reaction (qPCR) implied that salt can regulate the expression of key genes. This study generated a substantial fraction of D. salina transcriptional sequences for the entire growth cycle, providing a basis for the discovery of novel genes. This first full-scale transcriptome study of D. salina establishes a foundation for further comparative genomic studies.

The aquatic animals' transcriptome resource for comparative functional analysis

Background

Aquatic animals have great economic and ecological importance. Among them, non-model organisms have been studied regarding eco-toxicity, stress biology, and environmental adaptation. Due to recent advances in next-generation sequencing techniques, large amounts of RNA-seq data for aquatic animals are publicly available. However, currently there is no comprehensive resource exist for the analysis, unification, and integration of these datasets. This study utilizes computational approaches to build a new resource of transcriptomic maps for aquatic animals. This aquatic animal transcriptome map database dbATM provides de novo assembly of transcriptome, gene annotation and comparative analysis of more than twenty aquatic organisms without draft genome.

Results

To improve the assembly quality, three computational tools (Trinity, Oases and SOAPdenovo-Trans) were employed to enhance individual transcriptome assembly, and CAP3 and CD-HIT-EST software were then used to merge these three assembled transcriptomes. In addition, functional annotation analysis provides valuable clues to gene characteristics, including full-length transcript coding regions, conserved domains, gene ontology and KEGG pathways. Furthermore, all aquatic animal genes are essential for comparative genomics tasks such as constructing homologous gene groups and blast databases and phylogenetic analysis.

Conclusion

In conclusion, we establish a resource for non model organism aquatic animals, which is great economic and ecological importance and provide transcriptomic information including functional annotation and comparative transcriptome analysis. The database is now publically accessible through the URL http://dbATM.mbc.nctu.edu.tw/.

Electronic supplementary material

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

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

De Novo Transcriptome Assembly in Polyploid Species

In the absence of a reference genome, the ultimate goal of a de novo transcriptome assembly is to accurately and comprehensively reconstruct the set of messenger RNA transcripts represented in the sample. Non-reference assembly of the transcriptome of polyploid species poses a particular challenge because of the presence of homeologs that are difficult to disentangle at the sequence level. This is especially true for hexaploid oats, which have three highly similar subgenomes, two of which are thought to be nearly identical. Under these circumstances, most software packages and established pipelines encounter difficulties in rendering an accurate transcriptome because they are typically developed, refined, and tested for diploid organisms. We present a protocol for transcriptome assembly in oats that can be extended both to other polyploids and species with highly duplicated genomes.

Transcriptome sequencing reveals high isoform diversity in the ant Formica exsecta

Transcriptome resources for social insects have the potential to provide new insight into polyphenism, i.e., how divergent phenotypes arise from the same genome. Here we present a transcriptome based on paired-end RNA sequencing data for the ant Formica exsecta (Formicidae, Hymenoptera). The RNA sequencing libraries were constructed from samples of several life stages of both sexes and female castes of queens and workers, in order to maximize representation of expressed genes. We first compare the performance of common assembly and scaffolding software (Trinity, Velvet-Oases, and SOAPdenovo-trans), in producing de novo assemblies. Second, we annotate the resulting expressed contigs to the currently published genomes of ants, and other insects, including the honeybee, to filter genes that have annotation evidence of being true genes. Our pipeline resulted in a final assembly of altogether 39,262 mRNA transcripts, with an average coverage of >300X, belonging to 17,496 unique genes with annotation in the related ant species. From these genes, 536 genes were unique to one caste or sex only, highlighting the importance of comprehensive sampling. Our final assembly also showed expression of several splice variants in 6,975 genes, and we show that accounting for splice variants affects the outcome of downstream analyses such as gene ontologies. Our transcriptome provides an outstanding resource for future genetic studies on F. exsecta and other ant species, and the presented transcriptome assembly can be adapted to any non-model species that has genomic resources available from a related taxon.

De novo transcriptome sequencing of Isaria cateniannulata and comparative analysis of gene expression in response to heat and cold stresses

Isaria cateniannulata is a very important and virulent entomopathogenic fungus that infects many insect pest species. Although I. cateniannulata is commonly exposed to extreme environmental temperature conditions, little is known about its molecular response mechanism to temperature stress. Here, we sequenced and de novo assembled the transcriptome of I. cateniannulata in response to high and low temperature stresses using Illumina RNA-Seq technology. Our assembly encompassed 17,514 unigenes (mean length = 1,197 bp), in which 11,445 unigenes (65.34%) showed significant similarities to known sequences in NCBI non-redundant protein sequences (Nr) database. Using digital gene expression analysis, 4,483 differentially expressed genes (DEGs) were identified after heat treatment, including 2,905 up-regulated genes and 1,578 down-regulated genes. Under cold stress, 1,927 DEGs were identified, including 1,245 up-regulated genes and 682 down-regulated genes. The expression patterns of 18 randomly selected candidate DEGs resulting from quantitative real-time PCR (qRT-PCR) were consistent with their transcriptome analysis results. Although DEGs were involved in many pathways, we focused on the genes that were involved in endocytosis: In heat stress, the pathway of clathrin-dependent endocytosis (CDE) was active; however at low temperature stresses, the pathway of clathrin-independent endocytosis (CIE) was active. Besides, four categories of DEGs acting as temperature sensors were observed, including cell-wall-major-components-metabolism-related (CWMCMR) genes, heat shock protein (Hsp) genes, intracellular-compatible-solutes-metabolism-related (ICSMR) genes and glutathione S-transferase (GST). These results enhance our understanding of the molecular mechanisms of I. cateniannulata in response to temperature stresses and provide a valuable resource for the future investigations.

Convergent evolution of Y chromosome gene content in flies

Sex-chromosomes have formed repeatedly across Diptera from ordinary autosomes, and X-chromosomes mostly conserve their ancestral genes. Y-chromosomes are characterized by abundant gene-loss and an accumulation of repetitive DNA, yet the nature of the gene repertoire of fly Y-chromosomes is largely unknown. Here we trace gene-content evolution of Y-chromosomes across 22 Diptera species, using a subtraction pipeline that infers Y genes from male and female genome, and transcriptome data. Few genes remain on old Y-chromosomes, but the number of inferred Y-genes varies substantially between species. Young Y-chromosomes still show clear evidence of their autosomal origins, but most genes on old Y-chromosomes are not simply remnants of genes originally present on the proto-sex-chromosome that escaped degeneration, but instead were recruited secondarily from autosomes. Despite almost no overlap in Y-linked gene content in different species with independently formed sex-chromosomes, we find that Y-linked genes have evolved convergent gene functions associated with testis expression. Thus, male-specific selection appears as a dominant force shaping gene-content evolution of Y-chromosomes across fly species.

While X-chromosome gene content tends to be conserved, Y-chromosome evolution is dynamic and difficult to reconstruct. Here, Mahajan and Bachtrog use a subtraction pipeline to identify Y-linked genes in 22 Diptera species, revealing patterns of Y-chromosome gene-content evolution.

The first transcriptome of Italian wall lizard, a new tool to infer about the Island Syndrome

Some insular lizards show a high degree of differentiation from their conspecific mainland populations, like Licosa island lizards, which are described as affected by Reversed Island Syndrome (RIS). In previous works, we demonstrated that some traits of RIS, as melanization, depend on a differential expression of gene encoding melanocortin receptors. To better understand the basis of syndrome, and providing raw data for future investigations, we generate the first de novo transcriptome of the Italian wall lizard. Comparing mainland and island transcriptomes, we link differences in life-traits to differential gene expression. Our results, taking together testis and brain sequences, generated 275,310 and 269,885 transcripts, 18,434 and 21,606 proteins in Gene Ontology annotation, for mainland and island respectively. Variant calling analysis identified about the same number of SNPs in island and mainland population. Instead, through a differential gene expression analysis we found some putative genes involved in syndrome more expressed in insular samples like Major Histocompatibility Complex class I, Immunoglobulins, Melanocortin 4 receptor, Neuropeptide Y and Proliferating Cell Nuclear Antigen.

De novo assembly of Phlomis purpurea after challenging with Phytophthora cinnamomi

BACKGROUND

Phlomis plants are a source of biological active substances with potential applications in the control of phytopathogens. Phlomis purpurea (Lamiaceae) is autochthonous of southern Iberian Peninsula and Morocco and was found to be resistant to Phytophthora cinnamomi. Phlomis purpurea has revealed antagonistic effect in the rhizosphere of Quercus suber and Q. ilex against P. cinnamomi. Phlomis purpurea roots produce bioactive compounds exhibiting antitumor and anti-Phytophthora activities with potential to protect susceptible plants. Although these important capacities of P. purpurea have been demonstrated, there is no transcriptomic or genomic information available in public databases that could bring insights on the genes underlying this anti-oomycete activity.

RESULTS

Using Illumina technology we obtained a de novo assembly of P. purpurea transcriptome and differential transcript abundance to identify putative defence related genes in challenged versus non-challenged plants. A total of 1,272,600,000 reads from 18 cDNA libraries were merged and assembled into 215,739 transcript contigs. BLASTX alignment to Nr NCBI database identified 124,386 unique annotated transcripts (57.7%) with significant hits. Functional annotation identified 83,550 out of 124,386 unique transcripts, which were mapped to 141 pathways. 39% of unigenes were assigned GO terms. Their functions cover biological processes, cellular component and molecular functions. Genes associated with response to stimuli, cellular and primary metabolic processes, catalytic and transporter functions were among those identified. Differential transcript abundance analysis using DESeq revealed significant differences among libraries depending on post-challenge times. Comparative cyto-histological studies of P. purpurea roots challenged with P. cinnamomi zoospores and controls revealed specific morphological features (exodermal strips and epi-cuticular layer), that may provide a constitutive efficient barrier against pathogen penetration. Genes involved in cutin biosynthesis and in exodermal Casparian strips formation were up-regulated.

CONCLUSIONS

The de novo assembly of transcriptome using short reads for a non-model plant, P. purpurea, revealed many unique transcripts useful for further gene expression, biological function, genomics and functional genomics studies. The data presented suggest a combination of a constitutive resistance and an increased transcriptional response from P. purpurea when challenged with the pathogen. This knowledge opens new perspectives for the understanding of defence responses underlying pathogenic oomycete/plant interaction upon challenge with P. cinnamomi.

De novo transcriptome assembly for the spiny mouse (Acomys cahirinus)

Spiny mice of the genus Acomys display several unique physiological traits, including menstruation and scar-free wound healing; characteristics that are exceedingly rare in mammals, and of considerable interest to the scientific community. These unique attributes, and the potential for spiny mice to accurately model human diseases, are driving increased use of this genus in biomedical research, however little genetic information is accessible for this species. This project aimed to generate a draft transcriptome for the Common spiny mouse (Acomys cahirinus). Illumina sequencing of RNA from 15 organ types (male and female) produced 451 million, 150 bp paired-end reads (92.4Gbp). An extensive survey of de novo transcriptome assembly approaches using Trinity, SOAPdenovo-Trans, and Oases at multiple kmer lengths was conducted, producing 50 single-kmer assemblies from this dataset. Non-redundant transcripts from all assemblies were merged into a meta-assembly using the EvidentialGene tr2aacds pipeline, producing the largest gene catalogue to date for Acomys cahirinus. This study provides the first detailed characterization of the spiny mouse transcriptome. It validates use of the EvidentialGene tr2aacds pipeline in mammals to augment conventional de novo assembly approaches, and provides a valuable scientific resource for further investigation into the unique physiological characteristics inherent in the genus Acomys.

Phenotypic and Comparative Transcriptome Analysis of Different Ploidy Plants in Dendrocalamus latiflorus Munro

Elucidating the differences in gene expression profiles of plants with different ploidy levels and how they affect phenotypic traits is vital to allow genetic improvement of plants such as Ma bamboo (Dendrocalamus latiflorus Munro). We previously obtained triploid (2n = 3X = 36), hexaploid (2n = 6X = 72), and dodecaploid (2n = 12X = 144) Ma bamboo plants from embryogenic callus by anther culturing. Phenotypic differences between these plants appeared to be correlated with differences in ploidy. Here, we performed transcriptome profiling and sequencing of anther-regenerated plants and F1 seedlings of different ploidy levels using RNA-Seq technology. Pair-wise comparisons of the four resulting libraries revealed 8,396 differentially expressed genes. These differentially expressed genes were annotated, functionally classified, and partially validated. We found that the chromosome doubling led to substantially up- or down-regulation of genes that were involved in cell growth and differentiation; the polyploidy levels altered the anatomical, physiological and growth characteristics, such as leaf thickness, fusoid cell and stomatal size, shoot number, photosynthesis and respiration rate and so on. Additionally, two candidate genes, EXPB3 and TCP with potenitial regulatory roles in cell division and differentiation, were identified through gene coexpresseion network analysis. These results highlight the significance of potential applications of polyploidy, and provide valuable information for the genetic breeding of bamboo species.

The Euphausia superba transcriptome database, SuperbaSE: An online, open resource for researchers

Antarctic krill (Euphausia superba) is a crucial component of the Southern Ocean ecosystem, acting as the major link between primary production and higher trophic levels with an annual predator demand of up to 470 million tonnes. It also acts as an ecosystem engineer, affecting carbon sequestration and recycling iron and nitrogen, and has increasing importance as a commercial product in the aquaculture and health industries. Here we describe the creation of a de novo assembled head transcriptome for E. superba. As an example of its potential as a molecular resource, we relate its exploitation in identifying and characterizing numerous genes related to the circadian clock in E. superba, including the major components of the central feedback loop. We have made the transcriptome openly accessible for a wider audience of ecologists, molecular biologists, evolutionary geneticists, and others in a user-friendly format at SuperbaSE, hosted at http://www.krill.le.ac.uk.