SMRT sequencing of full-length transcriptome of seagrasses Zostera japonica

Diverse RNA viruses discovered in multiple seagrass species

Seagrasses are marine angiosperms that form highly productive and diverse ecosystems. These ecosystems, however, are declining worldwide. Plant-associated microbes affect critical functions like nutrient uptake and pathogen resistance, which has led to an interest in the seagrass microbiome. However, despite their significant role in plant ecology, viruses have only recently garnered attention in seagrass species. In this study, we produced original data and mined publicly available transcriptomes to advance our understanding of RNA viral diversity in Zostera marina, Zostera muelleri, Zostera japonica, and Cymodocea nodosa. In Z. marina, we present evidence for additional Zostera marina amalgavirus 1 and 2 genotypes, and a complete genome for an alphaendornavirus previously evidenced by an RNA-dependent RNA polymerase gene fragment. In Z. muelleri, we present evidence for a second complete alphaendornavirus and near complete furovirus. Both are novel, and, to the best of our knowledge, this marks the first report of a furovirus infection naturally occurring outside of cereal grasses. In Z. japonica, we discovered genome fragments that belong to a novel strain of cucumber mosaic virus, a prolific pathogen that depends largely on aphid vectoring for host-to-host transmission. Lastly, in C. nodosa, we discovered two contigs that belong to a novel virus in the family Betaflexiviridae. These findings expand our knowledge of viral diversity in seagrasses and provide insight into seagrass viral ecology.

RNA Transcript Diversity in Neuromuscular Research

Three decades since the Human Genome Project began, scientists have now identified more then 25,000 protein coding genes in the human genome. The vast majority of the protein coding genes (>  90%) are multi-exonic, with the coding DNA being interrupted by intronic sequences, which are removed from the pre-mRNA transcripts before being translated into proteins, a process called splicing maturation. Variations in this process, i.e. by exon skipping, intron retention, alternative 5' splice site (5'ss), 3' splice site (3'ss), or polyadenylation usage, lead to remarkable transcriptome and proteome diversity in human tissues. Given its critical biological importance, alternative splicing is tightly regulated in a tissue- and developmental stage-specific manner. The central nervous system and skeletal muscle are amongst the tissues with the highest number of differentially expressed alternative exons, revealing a remarkable degree of transcriptome complexity. It is therefore not surprising that splicing mis-regulation is causally associated with a myriad of neuromuscular diseases, including but not limited to amyotrophic lateral sclerosis (ALS), spinal muscular atrophy (SMA), Duchenne muscular dystrophy (DMD), and myotonic dystrophy type 1 and 2 (DM1, DM2). A gene's transcript diversity has since become an integral and an important consideration for drug design, development and therapy. In this review, we will discuss transcript diversity in the context of neuromuscular diseases and current approaches to address splicing mis-regulation.

Single-molecule real-time sequencing of the full-length transcriptome of Halophila beccarii

Ecologically, Halophila beccarii Asch. is considered as a colonizing or a pioneer seagrass species and a “tiny but mighty” seagrass species, since it may recover quickly from disturbance generally. The use of transcriptome technology can provide a better understanding of the physiological processes of seagrasses. To date, little is known about the genome and transcriptome information of H. beccarii. In this study, we used single molecule real-time (SMRT) sequencing to obtain full-length transcriptome data and characterize the transcriptome structure. A total of 11,773 of the 15,348 transcripts were successfully annotated in seven databases. In addition, 1573 long non-coding RNAs, 8402 simple sequence repeats and 2567 transcription factors were predicted in all the transcripts. A GO analysis showed that 5843 transcripts were divided into three categories, including biological process (BP), cellular component (CC) and molecular function (MF). In these three categories, metabolic process (1603 transcripts), protein-containing complex (515 transcripts) and binding (3233 transcripts) were the primary terms in BP, CC, and MF, respectively. The major types of transcription factors were involved in MYB-related and NF-YB families. To the best of our knowledge, this is the first report of the transcriptome of H. beccarii using SMRT sequencing technology.

Overexpression of Zostera japonica 14-3-3 gene ZjGRF1 enhances the resistance of transgenic Arabidopsis to copper stress

BACKGROUND

Copper is both a nutrient essential for plant growth and a pollutant. In recent decades, with the rapid development of industrial and agricultural production, copper has been used more and more widely, and its consumption has also increased rapidly. Excessive soil copper contents induce phytotoxicity, affecting plant growth, development and yields. Moreover, copper can accumulate in crops and enter the food chain through enrichment, harming human health.

METHODS AND RESULTS

In this study, Arabidopsis wild-type (WT) and Zostera japonica 14-3-3 gene ZjGRF1 overexpression lines were used to explore the physiological function and molecular mechanism of ZjGRF1 in Arabidopsis in the copper stress response. Under copper stress, compared with WT plants, transgenic ZjGRF1 Arabidopsis plants exhibited less inhibition of root growth and development and had higher fresh weights. Under copper stress, the soluble sugar and soluble protein contents in transgenic ZjGRF1 Arabidopsis plants were significantly higher than those in WT plants, while the superoxide dismutase (SOD), peroxidase and catalase (CAT) activities were significantly higher than those in WT plants. Additionally, the malonaldehyde content of transgenic plants was significantly lower than that of WT plants. Furthermore, qRT-PCR results showed that under copper stress, the SOD, CAT1 and HMA5 expression levels in transgenic ZjGRF1 Arabidopsis plants were significantly higher than those in WT plants, while COPT1 expression was significantly lower than that in WT plants.

CONCLUSIONS

ZjGRF1 enhanced the copper stress resistance of Arabidopsis by maintaining high antioxidant enzyme activity, increasing copper efflux and reducing copper uptake under copper stress.

Full-Length Transcriptome Maps of Reef-Building Coral Illuminate the Molecular Basis of Calcification, Symbiosis, and Circa-Dian Genes

Coral transcriptomic data largely rely on short-read sequencing, which severely limits the understanding of coral molecular mechanisms and leaves many important biological questions unresolved. Here, we sequence the full-length transcriptomes of four common and frequently dominant reef-building corals using the PacBio Sequel II platform. We obtain information on reported gene functions, structures, and expression profiles. Among them, a comparative analysis of biomineralization-related genes provides insights into the molecular basis of coral skeletal density. The gene expression profiles of the symbiont Symbiodiniaceae are also isolated and annotated from the holobiont sequence data. Finally, a phylogenetic analysis of key circadian clock genes among 40 evolutionarily representative species indicates that there are four key members in early metazoans, including cry genes; Clock or Npas2; cyc or Arntl; and tim, while per, as the fifth member, occurs in Bilateria. In summary, this work provides a foundation for further work on the manipulation of skeleton production or symbiosis to promote the survival of these important organisms.

De novo assembly and annotation of the transcriptome of the endangered seagrass Zostera capensis: Insights from differential gene expression under thermal stress

Seagrasses are important marine ecosystem engineers but anthropogenic impacts and climate change have led to numerous population declines globally. In South Africa, Zostera capensis is endangered due to fragmented populations and heavy anthropogenic pressures on estuarine ecosystems that house the core of the populations. Addressing questions of how pressures such as climate change affect foundational species, including Z. capensis are crucial to supporting their conservation and underpin restoration efforts. Here we use ecological transcriptomics to study key functional responses of Z. capensis through quantification of gene expression after thermal stress and present the first reference transcriptome of Z. capensis. Four de novo reference assemblies (Trinity, IDBA-tran, RNAspades, SOAPdenovo) filtered through the EvidentialGene pipeline resulted in 153,755 transcripts with a BUSCO score of 66.1% for completeness. Differential expression analysis between heat stressed (32 °C for three days) and pre-warming plants identified genes involved in photosynthesis, oxidative stress, translation, metabolic and biosynthetic processes in the Z. capensis thermal stress response. This reference transcriptome is a significant contribution to the limited available genomic resources for Z. capensis and represents a vital tool for addressing questions around the species restoration and potential functional responses to warming marine environments.

Exploring epitranscriptomics for crop improvement and environmental stress tolerance

Climate change and stressful environmental conditions severely hamper crop growth, development and yield. Plants respond to environmental perturbations, through their plasticity provided by key-genes, governed at post-/transcriptional levels. Gene-regulation in plants is a multilevel process controlled by diverse cellular entities that includes transcription factors (TF), epigenetic regulators and non-coding RNAs beside others. There are successful studies confirming the role of epigenetic modifications (DNA-methylation/histone-modifications) in gene expression. Recent years have witnessed emergence of a highly specialized field the "Epitranscriptomics". Epitranscriptomics deals with investigating post-transcriptional RNA chemical-modifications present across the life forms that change structural, functional and biological characters of RNA. However, deeper insights on of epitranscriptomic modifications, with >140 types known so far, are to be understood fully. Researchers have identified epitranscriptome marks (writers, erasers and readers) and mapped the site-specific RNA modifications (m6A, m⁵C, 3' uridylation, etc.) responsible for fine-tuning gene expression in plants. Simultaneous advancement in sequencing platforms, upgraded bioinformatic tools and pipelines along with conventional labelled techniques have further given a statistical picture of these epitranscriptomic modifications leading to their potential applicability in crop improvement and developing climate-smart crops. We present herein the insights on epitranscriptomic machinery in plants and how epitranscriptome and epitranscriptomic modifications underlying plant growth, development and environmental stress responses/adaptations. Third-generation sequencing technology, advanced bioinformatics tools and databases being used in plant epitranscriptomics are also discussed. Emphasis is given on potential exploration of epitranscriptome engineering for crop-improvement and developing environmental stress tolerant plants covering current status, challenges and future directions.

Overexpression of Zostera japonica heat shock protein gene ZjHsp70 enhances the thermotolerance of transgenic Arabidopsis

BACKGROUND

Heat shock protein 70s (Hsp70s) are major members of the heat shock protein family and play a variety of roles to protect plants against stress. Plant Hsp70s are a conserved and widely expressed family of heat shock proteins. They have two main functional regions: N-terminal nucleic acid binding region and C-terminal substrate binding region.

METHODS AND RESULTS

In this study, we cloned the Hsp70 gene of Zostera japonica (ZjHsp70) based on the sequence obtained by transcriptome sequencing. The transcriptional levels of ZjHsp70 increased significantly at 1 h after heat treatment. ZjHsp70 was located in the cytoplasm and nucleus. The overexpression of ZjHsp70 in Arabidopsis resulted in increased heat tolerance, lower contents of malondialdehyde and higher antioxidant enzyme activity than in the wild type. ZjHsp70 may achieve this goal by maintaining highly active antioxidant enzymes.

CONCLUSIONS

We show that ZjHsp70 can improve plant heat tolerance by maintaining high antioxidant enzyme activity under high temperature stress. This study provided a basis to study the role of ZjHsp70 in thermotolerance in more detail.

SMRT Sequencing Reveals Candidate Genes and Pathways With Medicinal Value in Cipangopaludina chinensis

Cipangopaludina chinensis is an economically important aquatic snail with high medicinal value. However, molecular biology research on C. chinensis is limited by the lack of a reference genome, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. chinensis. Herein, we conducted the first full-length transcriptome analysis of C. chinensis using PacBio single-molecule real-time (SMRT) sequencing technology. We identified a total of 26,312 unigenes with an average length of 2,572 bp, of which the largest number of zf-c2h2 transcription factor families (120,18.24%) were found, and also observed that the majority of the 8,058 SSRs contained 4-7 repeat units, which provided data for subsequent work on snail genetics Subsequently, 91.86% (24,169) of the genes were successfully annotated to the four major databases, while the highest homology was observed with Pomacea canaliculata. Functional annotation revealed that the majority of transcripts were enriched in metabolism, signal transduction and Immune-related pathways, and several candidate genes involved in drug metabolism and immune response were identified (e.g., CYP1A1, CYP2J, CYP2U1, GST, ,PIK3, PDE3A, PRKAG). This study lays a foundation for future molecular biology research and provides a reference for studying genes associated with the medicinal value of C. chinensis.

SMRT sequencing of full-length transcriptome and gene expression analysis in two chemical types of Pogostemon cablin (Blanco) Benth

BACKGROUND

Pogostemon cablin (Blanco) Benth. also called patchouli, is a traditional medicinal and aromatic plant that grows mainly in Southeast Asia and China. In China, P. cablin is divided into two chemical types: the patchouliol-type and the pogostone-type. Patchouliol-type patchouli usually grow taller, with thicker stems and bigger leaves, and produce more aromatic oil.

METHODS

To better understand the genetic differences between the two chemical types that contribute to their differences in morphology and biosynthetic capabilities, we constructed de novo transcriptomes from both chemical types using the Pacific Biosciences (PacBio) Sequel platform and performed differential expression analysis of multiple tissues using Illumina short reads.

RESULTS

In this study, using single-molecule real-time (SMRT) long-read sequencing, we obtained 22.07 GB of clean data and 134,647 nonredundant transcripts from two chemical types. Additionally, we identified 126,576 open reading frames (ORFs), 100,638 coding sequences (CDSs), 4,106 long noncoding RNAs (lncRNAs) and 6,829 transcription factors (TFs) from two chemical types of P. cablin. We adopted PacBio and Illumina sequencing to identify differentially expressed transcripts (DEGs) in three tissues of the two chemical types. More DEGs were observed in comparisons of different tissues collected from the same chemical type relative to comparisons of the same tissue collected from different chemical types. Furthormore, using KEGG enrichment analysis of DEGs, we found that the most enriched biosynthetic pathways of secondary metabolites of the two chemical types were "terpenoid backbone biosynthesis", "phenylpropanoid biosynthesis", "plant hormone signal transduction", "sesquiterpenoid and triterpenoid biosynthesis", "ubiquinone and other terpenoid-quinone biosynthesis", "flavonoid biosynthesis", and "flavone and flavonol biosynthesis". However, the main pathways of the patchouliol-type also included "diterpene biosynthesis" and "monoterpene biosynthesis". Additionally, by comparing the expression levels of the three tissues verified by qRT-PCR, more DEGs in the roots were upregulated in the mevalonate (MVA) pathway in the cytoplasm, but more DEGs in the leaves were upregulated in the methylerythritol phosphate (MEP) pathway in the plastid, both of which are important pathways for terpenoids biosynthesis. These findings promote the study of further genome annotation and transcriptome research in P. cablin.

Sensitive and accurate analysis of gene expression signatures enabled by oligonucleotide-labelled cDNA

High-throughput RNA sequencing offers a comprehensive analysis of transcriptome complexity originated from regulatory events, such as differential gene expression, alternative polyadenylation and others, and allows the increase in diagnostic capacity and precision. For gene expression profiling applications that do not specifically require information on alternative splicing events, the mRNA 3' termini counting approach is a cost-effective alternative to whole transcriptome sequencing. Here, we report MTAS-seq (mRNA sequencing via terminator-assisted synthesis) - a novel RNA-seq library preparation method directed towards mRNA 3' termini. We demonstrate the specific enrichment for 3'-terminal regions by simple and quick single-tube protocol with built-in molecular barcoding to enable accurate estimation of transcript abundance. To achieve that, we synthesized oligonucleotide-modified dideoxynucleotides which enable the generation of cDNA libraries at the reverse transcription step. We validated the performance of MTAS-seq on well-characterized reference bulk RNA and further tested it with eukaryotic cell lysates.

Overexpression of seagrass DnaJ gene ZjDjB1 enhances the thermotolerance of transgenic arabidopsis thaliana

Seagrass meadows are one of the most important marine resources that grow along the coast. They provide habitat and a food source for animals. They also protect the coast, fix sediment and purify seawater. In the current period of global climate change, anomalies in coastal water temperatures are increasing. A sudden increase in water temperature owing to a heat wave can have a profound effect on seagrass. Zostera japonica is a type of intertidal seagrasses, which is exposed to the air at low tide. High temperatures in the summer often lead to a decline in seagrass meadows. DnaJ proteins, also known as J proteins, are a family of conserved chaperone proteins. They are designated as J proteins because they contain a highly conserved J domain. They function as chaperones of heat shock proteins in organisms. In this study, the role of DnaJ protein (ZjDjB1) of Z. japonica under heat stress was studied. ZjDjB1 was localized to the cytoplasm and nucleus. The overexpression of ZjDjB1 in Arabidopsis thaliana results in an increase in thermotolerance and a decrease in the accumulation of reactive oxygen species and also a reduction in membrane damage. ZjDjB1 may achieve this goal by maintaining a low activity of proteolytic enzymes.

Full Transcriptome Analysis of Callus Suspension Culture System of Bletilla striata

Background

Bletilla striata has been widely used in the pharmacology industry. To effectively produce the secondary metabolites through suspension cultured cells of B. striata, it is important to exploring the full-length transcriptome data and the genes related to cell growth and chemical producing of all culture stages. We applied a combination of Real-Time Sequencing of Single Molecule (SMRT) and second-generation sequencing (SGS) to generate the complete and full-length transcriptome of B. striata suspension cultured cells.

Methods

The B. striata transcriptome was formed in de novo way by using PacBio isoform sequencing (Iso-Seq) on a pooled RNA sample derived from 23 samples of 10 culture stages, to explore the potential for capturing full-length transcript isoforms. All unigenes were obtained after splicing, assembling, and clustering, and corrected by the SGS results. The obtained unigenes were compared with the databases, and the functions were annotated and classified.

Results and conclusions

A total of 100,276 high-quality full-length transcripts were obtained, with an average length of 2530 bp and an N50 of 3302 bp. About 52% of total sequences were annotated against the Gene Ontology, 53,316 unigenes were hit by KOG annotations and divided into 26 functional categories, 80,020 unigenes were mapped by KEGG annotations and clustered into 363 pathways. Furthermore, 15,133 long-chain non-coding RNAs (lncRNAs) were detected. And 68,996 coding sequences were identified based on SSR analysis, among which 31 pairs of primers selected at random were amplified and obtained stable bands. In conclusion, our results provide new full-length transcriptome data and genetic resources for identifying growth and metabolism-related genes, which provide a solid foundation for further research on its growth regulation mechanisms and genetic engineering breeding mechanisms of B. striata.