Full length transcriptomes analysis of cold-resistance of Apis cerana in Changbai Mountain during overwintering period

Apis cerana in Changbai Mountain is an ecological type of Apis cerana, which is an excellent breeding material with cold-resistant developed by long-term natural selection under the ecological conditions. However, the physiological and molecular mechanisms of Changbai Mountain population under cold stress are still unclear. In this study, the Nanopore sequencing was carried out for the transcriptome of Apis cerana in Changbai Mountain in the coldest period of overwintering, which will provide a reference to the cold-resistant mechanism. We determined 5,941 complete ORF sequences, 1,193 lncRNAs, 619 TFs, 10,866 SSRs and functional annotations of 11,599 new transcripts. Our results showed that the myosin family and the C2H2 zinc finger protein transcription factor family possibly have significant impacts on the response mechanism of cold stress during overwintering. In addition, the cold environment alters genes expression profiles in honeybees via different AS and APA mechanisms. These altered genes in Hippo, Foxo, and MARK pathways help them counter the stress of cold in overwinter period. Our results might provide clues about the response of eastern honeybees to extreme cold, and reflect the possible genetic basis of physiological changes.

The Long Read Transcriptome of Rice (Oryza sativa ssp. japonica var. Nipponbare) Reveals Novel Transcripts

BACKGROUND

High-throughput next-generation sequencing technologies offer a powerful approach to characterizing the transcriptomes of plants. Long read sequencing has been shown to support the discovery of novel isoforms of transcripts. This approach enables the generation of full-length sequences revealing splice variants that may be important in regulating gene action. Investigation of the diversity of transcripts in the rice transcriptome including splice variants was conducted using PacBio long-read sequence data to improve the annotation of the rice genome.

RESULTS

A cDNA library was prepared from RNA extracted from leaves, roots, seeds, inflorescences, and panicles of O. sativa ssp. japonica var Nipponbare and sequenced on a PacBio Sequel platform. This produced 346,190 non-redundant full-length non-chimeric reads (FLNC) resulting in 33,504 high-quality transcripts. Half of the transcripts were multi-exonic and entirely matched with the reference transcripts. However, 14,874 novel isoforms were also identified resulting predominantly from intron retention and at least one novel splice site. Intron retention was the prevalent alternative splicing event and exon skipping was the least observed. Of 73,659 splice junctions, 12,755 (17%) represented novel splice junctions with canonical and non-canonical intron boundaries. The complexity of the transcriptome was examined in detail for 19 starch synthesis-related genes, defining 276 spliced isoforms of which 94 splice variants were novel.

CONCLUSION

The data reveal the great complexity of the rice transcriptome. The novel transcripts provide new insights that may be a key input in future research to improve the annotation of the rice genome.

Combination analysis of single-molecule long-read and Illumina sequencing provides insights into the anthocyanin accumulation mechanism in an ornamental grass, Pennisetum setaceum cv. Rubrum

Combination analysis of single-molecule long-read and Illumina sequencing provide full-length transcriptome information and shed new light on the anthocyanin accumulation mechanism of Pennisetum setaceum cv. 'Rubrum'. Pennisetum setaceum cv. 'Rubrum' is an ornamental grass with purple leaves widely used in landscaping. However, the current next-generation sequencing (NGS) transcriptome information of this species is not satisfactory due to the difficulties in obtaining full-length transcripts. Furthermore, the molecular mechanisms of anthocyanin accumulation in P. setaceum have not been thoroughly studied. In this study, we used PacBio full-length transcriptome sequencing (SMRT) combined with NGS technology to build and improve the transcriptomic datasets and reveal the molecular mechanism of anthocyanin accumulation in P. setaceum cv. 'Rubrum'. Therefore, 280,413 full-length non-chimeric reads sequences were obtained using the SMRT technology. We obtained 97,450 high-quality non-redundant transcripts and identified 5352 alternative splicing events. In addition, 93,066 open reading frames (ORFs), including 57,457 full ORFs and 2910 long non-coding RNA (lncRNAs) were screened out. Furthermore, 10,795 differentially expressed genes were identified using NGS. We also explored key genes, synthesis pathways, and detected lncRNA involved in anthocyanin accumulation, providing new insights into anthocyanin accumulation in P. setaceum cv. 'Rubrum'. To our best knowledge, we provided the full-length transcriptome information of P. setaceum cv. 'Rubrum' for the first time. The results of this study will provide baseline information for gene function studies and pave the way for future P. setaceum cv. 'Rubrum' breeding projects.

Identification of genes involved in drought tolerance in seedlings of the desert grass, Psammochloa villosa (Poaceae), based on full-length isoform sequencing and de novo assembly from short reads

Psammochloa villosa is a perennial herbaceous plant that is dominant within arid regions of the Inner Mongolian Plateau and the Qinghai-Tibet Plateau in China, where it is an endemic species and exhibits strong drought tolerance and wind resistance. To study drought tolerance in P. villosa and determine its molecular basis, we simulated high and moderate drought stress in a controlled environment and then analyzed transcriptome sequences by combining long-read sequences from a representative, wild-grown individual with short reads from the treatment groups. We obtained 184,076 high-quality isoforms as a reference and 168,650 genes (91.6%), which we were able to annotate according to public databases. Ultimately, we obtained 119,005 unigenes representing the transcriptome of P. villosa under drought stress and, among these, we identified 3089 differentially expressed genes and 1484 transcription factors. Physiologically, P. villosa that was exposed to high and moderate drought stress had reduced germination rates and shorter buds but generated more chlorophyll, which is atypical under drought stress and possibly reflects an adaptation of these plants to their arid environment. We inferred that significantly upregulated genes were annotated as 'Chlorophyll a-b binding protein' and 'Light-harvesting chlorophyll-protein' among drought and control groups. Broadly, our analyses revealed that drought stress triggered many genome-level responses, especially related to mitigation of radical oxygen species (ROS), which increase in concentration under drought stress. In particular, in the high drought stress group compared with the control, GO enrichment analysis revealed a significant enrichment of upregulated genes (n = 10) involved in mitigation of oxidative stress. Similarly, using KEGG we found significant enrichment of genes in the phenylpropanoid biosynthesis pathway (11 genes), which yields phenols that scavenge ROS. We also inferred that many genes involved in metabolism of arginine and proline, which may serve as both scavengers of ROS and osmoprotectants that interact with stress response genes based on our protein-protein interaction network analysis. We verified the relative expression levels of eight genes associated with mitigation of ROS, DNA repair, and transmembrane transporter activity using qRT-PCR, and the results were consistent with our inferences from transcriptomes. This study provides insights into the genomic and physiological basis of drought tolerance in P. villosa and represents a resource for development of the species as a forage crop via molecular breeding within arid lands.

Applications of Oxford Nanopore Sequencing in Schizosaccharomyces pombe

Recent years have seen great progresses in third-generation sequencing. New commercial platforms from Oxford Nanopore Technologies (ONT) can generate ultra-long reads from single-molecule nucleic acid fragments of kilobases up to megabases, exceeding the limitation of short reads and dependency on template amplification suffered by the previous generation of sequencing technologies. Moreover, it can detect epigenetic modifications directly, as well as providing all-around field usage, being pocket-sized and low cost. It has already been applied to yeast research in many aspects, such as complete de novo genome assemblies, the phylogeny of large-brewing yeasts, gene isoform identification, and base modification detection. These applications have delivered novel insights into yeast genomic and transcriptomic analysis.

Full-length transcriptome analysis of Misgurnus anguillicaudatus

Faithful analysis of transcripts represents a critical component of the gene discovery, genomic characterization and annotation process in species without a reference genome. The present study generated 28,001 full-length transcripts of Misgurnus anguillicaudatus using the Pacific Bioscience (PacBio) single-molecule real-time (SMRT) sequencing system. A total of 77,346 consensus isoforms were identified, and 18,991 complete open reading frames (ORFs) detected. In total, 1345 lncRNAs were identified with high-confidence, with additional identification of a number of well-known transcription factors. The present study may facilitate additional exploration of the genomic signatures of M. anguillicaudatus.

Coordinated mechanisms of leaves and roots in response to drought stress underlying full-length transcriptome profiling in Vicia sativa L

BACKGROUND

Common vetch (Vicia sativa L.) is an important self-pollinating annual forage legume and is of interest for drought prone regions as a protein source to feed livestock and human consumption. However, the development and production of common vetch are negatively affected by drought stress. Plants have evolved common or distinct metabolic pathways between the aboveground and underground in response to drought stress. Little is known regarding the coordinated response of aboveground and underground tissues of common vetch to drought stress.

RESULTS

Our results showed that a total of 30,427 full-length transcripts were identified in 12 samples, with an average length of 2278.89 bp. Global transcriptional profiles of the above 12 samples were then analysed via Illumina-Seq. A total of 3464 and 3062 differentially expressed genes (DEGs) were identified in the leaves and roots, respectively. Gene Ontology (GO) enrichment analyses identified that the dehydrin genes and Δ¹-pyrroline-5-carboxylate synthase were induced for the biosynthesis of proline and water conservation. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis results indicated that the DEGs were significantly enriched in hormone signal transduction, starch and sucrose metabolism, and arginine and proline metabolism, and various drought response candidate genes were also identified. Abscisic acid (ABA; the AREB/ABF-SnRK2 pathway) regulates the activity of AMY3 and BAM1 to induce starch degradation in leaves and increase carbon export to roots, which may be associated with the drought stress responses in common vetch. Among the co-induced transcription factors (TFs), AREB/ABF, bHLH, MYB, WRKY, and AP2/ERF had divergent expression patterns and may be key in the crosstalk between leaves and roots during adaption to drought stress. In transgenic yeast, the overexpression of four TFs increased yeast tolerance to osmotic stresses.

CONCLUSION

The multipronged approach identified in the leaves and roots broadens our understanding of the coordinated mechanisms of drought response in common vetch, and further provides targets to improve drought resistance through genetic engineering.

Identification of putative genes for polyphenol biosynthesis in olive fruits and leaves using full-length transcriptome sequencing

Olive (Olea europaea) is a rich source of valuable bioactive polyphenols, which has attracted widespread interest. In this study, we combined targeted metabolome, Pacbio ISOseq transcriptome, and Illumina RNA-seq transcriptome to investigate the association between polyphenols and gene expression in the developing olive fruits and leaves. A total of 12 main polyphenols were measured, and 122 transcripts of 17 gene families, 101 transcripts of 9 gene families, and 106 transcripts of 6 gene families that encode for enzymes involved in flavonoid, oleuropein, and hydroxytyrosol biosynthesis were separately identified. Additionally, 232 alternative splicing events of 18 genes related to polyphenol synthesis were analyzed. This is the first time that the third generations of full-length transcriptome technology were used to study the gene expression pattern of olive fruits and leaves. The results of transcriptome combined with targeted metabolome can help us better understand the polyphenol biosynthesis pathways in the olive.

Full-length transcript sequencing and comparative transcriptomic analysis to evaluate the contribution of osmotic and ionic stress components towards salinity tolerance in the roots of cultivated alfalfa (Medicago sativa L.)

BACKGROUND

Alfalfa is the most extensively cultivated forage legume. Salinity is a major environmental factor that impacts on alfalfa's productivity. However, little is known about the molecular mechanisms underlying alfalfa responses to salinity, especially the relative contribution of the two important components of osmotic and ionic stress.

RESULTS

In this study, we constructed the first full-length transcriptome database for alfalfa root tips under continuous NaCl and mannitol treatments for 1, 3, 6, 12, and 24 h (three biological replicates for each time points, including the control group) via PacBio Iso-Seq. This resulted in the identification of 52,787 full-length transcripts, with an average length of 2551 bp. Global transcriptional changes in the same 33 stressed samples were then analyzed via BGISEQ-500 RNA-Seq. Totals of 8861 NaCl-regulated and 8016 mannitol-regulated differentially expressed genes (DEGs) were identified. Metabolic analyses revealed that these DEGs overlapped or diverged in the cascades of molecular networks involved in signal perception, signal transduction, transcriptional regulation, and antioxidative defense. Notably, several well characterized signalling pathways, such as CDPK, MAPK, CIPK, and PYL-PP2C-SnRK2, were shown to be involved in osmotic stress, while the SOS core pathway was activated by ionic stress. Moreover, the physiological shifts of catalase and peroxidase activity, glutathione and proline content were in accordance with dynamic transcript profiles of the relevant genes, indicating that antioxidative defense system plays critical roles in response to salinity stress.

CONCLUSIONS

Overall, our study provides evidence that the response to salinity stress in alfalfa includes both osmotic and ionic components. The key osmotic and ionic stress-related genes are candidates for future studies as potential targets to improve resistance to salinity stress via genetic engineering.

Uncovering full-length transcript isoforms of sugarcane cultivar Khon Kaen 3 using single-molecule long-read sequencing

Background

Sugarcane is an important global food crop and energy resource. To facilitate the sugarcane improvement program, genome and gene information are important for studying traits at the molecular level. Most currently available transcriptome data for sugarcane were generated using second-generation sequencing platforms, which provide short reads. The de novo assembled transcripts from these data are limited in length, and hence may be incomplete and inaccurate, especially for long RNAs.

Methods

We generated a transcriptome dataset of leaf tissue from a commercial Thai sugarcane cultivar Khon Kaen 3 (KK3) using PacBio RS II single-molecule long-read sequencing by the Iso-Seq method. Short-read RNA-Seq data were generated from the same RNA sample using the Ion Proton platform for reducing base calling errors.

Results

A total of 119,339 error-corrected transcripts were generated with the N50 length of 3,611 bp, which is on average longer than any previously reported sugarcane transcriptome dataset. 110,253 sequences (92.4%) contain an open reading frame (ORF) of at least 300 bp long with ORF N50 of 1,416 bp. The mean lengths of 5′ and 3′ untranslated regions in 73,795 sequences with complete ORFs are 1,249 and 1,187 bp, respectively. 4,774 transcripts are putatively novel full-length transcripts which do not match with a previous Iso-Seq study of sugarcane. We annotated the functions of 68,962 putative full-length transcripts with at least 90% coverage when compared with homologous protein coding sequences in other plants.

Discussion

The new catalog of transcripts will be useful for genome annotation, identification of splicing variants, SNP identification, and other research pertaining to the sugarcane improvement program. The putatively novel transcripts suggest unique features of KK3, although more data from different tissues and stages of development are needed to establish a reference transcriptome of this cultivar.