De novo transcriptome sequencing and analysis of genes related to salt stress response in Glehnia littoralis

Screening and identification of multiple abiotic stress responsive candidate genes based on hybrid-sequencing in Vicia sativa

Common vetch is an important leguminous forage for both livestock fodder and green manure and has a tremendous latent capacity in a sustainable agroecosystem. In the present study, a comprehensive transcriptome analysis of the aboveground leaves and underground roots of common vetch under multiple abiotic stress treatments, including NaCl, drought, cold, and cold drought, was performed using hybrid-sequencing technology, i. e. single-molecule real-time sequencing technology (SMRT) and supplemented by next-generation sequencing (NGS) technology. A total of 485,038 reads of insert (ROIs) with a mean length of 2606 bp and 228,261 full-length nonchimeric (FLNC) reads were generated. After deduplication, 39,709 transcripts were generated. Of these transcripts, we identified 1059 alternative splicing (AS) events, 17,227 simple sequence repeats (SSRs), and 1647 putative transcription factors (TFs). Furthermore, 640 candidates long noncoding RNAs (lncRNAs) and 28,256 complete coding sequences (CDSs) were identified. In gene annotation analyses, a total of 38,826 transcripts (97.78%) were annotated in eight public databases. Finally, seven multiple abiotic stress-responsive candidate genes were obtained through gene expression, annotation information, and protein-protein interaction (PPI) networks. Our research not only enriched the structural information of FL transcripts in common vetch, but also provided useful information for exploring the molecular mechanism of multiple abiotic stress tolerance between aboveground and underground tissues in common vetch and related legumes.

Physiological Measurements and Transcriptome Survey Reveal How Semi-mangrove Clerodendrum inerme Tolerates Saline Adversity

Salinity adversity has been a major environmental stressor for plant growth and reproduction worldwide. Semi-mangrove Clerodendrum inerme, a naturally salt-tolerant plant, can be studied as a successful example to understand the biological mechanism of saline resistance. Since it is a sophisticated and all-round scale process for plants to react to stress, our greenhouse study interpreted the response of C. inerme to salt challenge in the following aspects: morphology, osmotic protectants, ROS production and scavenging, ion homeostasis, photosynthetic efficiency, and transcriptome reprogramming. The results drew an overview picture to illustrate the tolerant performance of C. inerme from salt acclimatization (till medium NaCl level, 0.3 mol/L) to salinity stress (high NaCl level, 0.5 mol/L). The overall evaluation leads to a conclusion that the main survival strategy of C. inerme is globally reshaping metabolic and ion profiles to adapt to saline adversity. These findings uncover the defense mechanism by which C. inerme moderates its development rate to resist the short- and long-term salt adversity, along with rebalancing the energy allocation between growth and stress tolerance.

Integrative RNA-Seq analysis of Capsicum annuum L.-Phytophthora capsici L. pathosystem reveals molecular cross-talk and activation of host defence response

Chili pepper (Capsicum annuum L.) is economically one of the most important spice. But, it's productivity is highly affected by the pathogen, Phytophthora capsici L. Our current understanding of the molecular mechanisms associated with the defence response in C. annuum-P. capsici pathosystem is limited. The current study used RNA-seq technology to dissect the genes associated with defence response against P. capsici infection in two contrasting landraces, i.e. GojamMecha_9086 (Resistant) and Dabat_80045 (Susceptible) exposed to P. capsici infection. The transcriptomes from four leaf samples (RC, RI, SC and SI) of chili pepper resulted in a total of 118,879 assembled transcripts along with 52,384 pooled unigenes. The enrichment analysis of the transcripts indicated 23 different KEGG pathways under five main categories. Out of 774 and 484 differentially expressed genes (DEGs) of two landraces (under study), respectively, 57 and 29 DEGs were observed as associated with defence responses against P. capsici infection in RC vs. RI and SC vs. SI leaf samples, respectively. qRT-PCR analysis of six randomly selected genes validated the results of Illumina NextSeq500 sequencing. A total of 58 transcription factor families (bHLH most abundant) and 2095 protein families (Protein kinase most abundant) were observed across all the samples with maximum hits in RI and SI samples. Expression analysis revealed differential regulation of genes associated with defence and signalling response with shared coordination of molecular function, cellular component and biological processing. The results presented here would enhance our present understanding of the defence response in chili pepper against P. capsici infection, which the molecular breeders could utilize to develop resistant chili genotypes.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01122-y.

Characterization of the Glehnia littoralis Non-specific Phospholipase C Gene GlNPC3 and Its Involvement in the Salt Stress Response

Glehnia littoralis is a medicinal halophyte that inhabits sandy beaches and has high ecological and commercial value. However, the molecular mechanism of salt adaptation in G. littoralis remains largely unknown. Here, we cloned and identified a non-specific phospholipase C gene (GlNPC3) from G. littoralis, which conferred lipid-mediated signaling during the salt stress response. The expression of GlNPC3 was induced continuously by salt treatment. Overexpression of GlNPC3 in Arabidopsis thaliana increased salt tolerance compared to wild-type (WT) plants. GlNPC3-overexpressing plants had longer roots and higher fresh and dry masses under the salt treatment. The GlNPC3 expression pattern revealed that the gene was expressed in most G. littoralis tissues, particularly in roots. The subcellular localization of GlNPC3 was mainly at the plasma membrane, and partially at the tonoplast. GlNPC3 hydrolyzed common membrane phospholipids, such as phosphotidylserine (PS), phosphoethanolamine (PE), and phosphocholine (PC). In vitro enzymatic assay showed salt-induced total non-specific phospholipase C (NPC) activation in A. thaliana GlNPC3-overexpressing plants. Plant lipid profiling showed a significant change in the membrane-lipid composition of A. thaliana GlNPC3-overexpressing plants compared to WT after the salt treatment. Furthermore, downregulation of GlNPC3 expression by virus-induced gene silencing in G. littoralis reduced the expression levels of some stress-related genes, such as SnRK2, P5SC5, TPC1, and SOS1. Together, these results indicated that GlNPC3 and GlNPC3-mediated membrane lipid change played a positive role in the response of G. littoralis to a saline environment.

Identification of differentially expressed genes and pathways in isonuclear kenaf genotypes under salt stress

Salinity is a potential abiotic stress and globally threatens crop productivity. However, the molecular mechanisms underlying salt stress tolerance with respect to cytoplasmic effect, gene expression, and metabolism pathway under salt stress have not yet been reported in isonuclear kenaf genotypes. To fill this knowledge gap, growth, physiological, biochemical, transcriptome, and cytoplasm changes in kenaf cytoplasmic male sterile (CMS) line (P3A) and its iso-nuclear maintainer line (P3B) under 200 mM sodium chloride (NaCl) stress and control conditions were analyzed. Salt stress significantly reduced leaf soluble protein, soluble sugars, proline, chlorophyll content, antioxidant enzymatic activity, and induced oxidative damage in terms of higher MDA contents in both genotypes. The reduction of these parameters resulted in a reduced plant growth compared with control. However, P3A was relatively more tolerant to salt stress than P3B. This tolerance of P3A was further confirmed by improved physio-biochemical traits under salt stress conditions. Transcriptome analysis showed that 4256 differentially expressed genes (DEGs) between the two genotypes under salt stress were identified. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that photosynthesis, photosynthesis antenna-protein, and plant hormone signal transduction pathways might be associated with the improved NaCl stress tolerance in P3A. Conclusively, P3A cytoplasmic male sterile could be a potential salt-tolerant material for future breeding program of kenaf and can be used for phytoremediation of salt-affected soils. These data provide a valuable resource on the cytoplasmic effect of salt-responsive genes in kenaf and salt stress tolerance in kenaf.

Illuminating the Plant Rhabdovirus Landscape through Metatranscriptomics Data

Rhabdoviruses infect a large number of plant species and cause significant crop diseases. They have a negative-sense, single-stranded unsegmented or bisegmented RNA genome. The number of plant-associated rhabdovirid sequences has grown in the last few years in concert with the extensive use of high-throughput sequencing platforms. Here, we report the discovery of 27 novel rhabdovirus genomes associated with 25 different host plant species and one insect, which were hidden in public databases. These viral sequences were identified through homology searches in more than 3000 plant and insect transcriptomes from the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) using known plant rhabdovirus sequences as the query. The identification, assembly and curation of raw SRA reads resulted in sixteen viral genome sequences with full-length coding regions and ten partial genomes. Highlights of the obtained sequences include viruses with unique and novel genome organizations among known plant rhabdoviruses. Phylogenetic analysis showed that thirteen of the novel viruses were related to cytorhabdoviruses, one to alphanucleorhabdoviruses, five to betanucleorhabdoviruses, one to dichorhaviruses and seven to varicosaviruses. These findings resulted in the most complete phylogeny of plant rhabdoviruses to date and shed new light on the phylogenetic relationships and evolutionary landscape of this group of plant viruses. Furthermore, this study provided additional evidence for the complexity and diversity of plant rhabdovirus genomes and demonstrated that analyzing SRA public data provides an invaluable tool to accelerate virus discovery, gain evolutionary insights and refine virus taxonomy.

Discovery of putative novel viruses in the transcriptomes of endangered plant species native to India and China

Viruses are abundant entities that infect almost every living organism. In recent years, Next Generation Sequencing coupled with bioinformatic analyses is widely adopted for identification of known and unknown viruses in a plant sample. In the present study, nine putative novel viruses were discovered from public domain transcriptome datasets of five endangered plant species by de novo assembly of reads using CLC and SPAdes followed by BLAST analysis. Of the identified viruses, ten coding-complete and five partial genomic segments were recovered. Based on phylogeny and BLAST analysis, the identified viruses were putatively assigned to various plant viral genera except dactylorhiza hatagirea benylike virus that probably represents a new group of plant virus. The methodology followed can be adopted for the discovery of novel viruses in plant species with little genomic information. Viral genome sequences recovered in the study will serve as a valuable resource for further characterization of identified viruses.

Investigation of an Antioxidative System for Salinity Tolerance in Oenanthe javanica

Abiotic stress, such as drought and salinity, severely affect the growth and yield of many plants. Oenanthe javanica (commonly known as water dropwort) is an important vegetable that is grown in the saline-alkali soils of East Asia, where salinity is the limiting environmental factor. To study the defense mechanism of salt stress responses in water dropwort, we studied two water dropwort cultivars, V11E0022 and V11E0135, based on phenotypic and physiological indexes. We found that V11E0022 were tolerant to salt stress, as a result of good antioxidant defense system in the form of osmolyte (proline), antioxidants (polyphenols and flavonoids), and antioxidant enzymes (APX and CAT), which provided novel insights for salt-tolerant mechanisms. Then, a comparative transcriptomic analysis was conducted, and Gene Ontology (GO) analysis revealed that differentially expressed genes (DEGs) involved in the carbohydrate metabolic process could reduce oxidative stress and enhance energy production that can help in adaptation against salt stress. Similarly, lipid metabolic processes can also enhance tolerance against salt stress by reducing the transpiration rate, H₂O₂, and oxidative stress. Furthermore, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that DEGs involved in hormone signals transduction pathway promoted the activities of antioxidant enzymes and reduced oxidative stress; likewise, arginine and proline metabolism, and flavonoid pathways also stimulated the biosynthesis of proline and flavonoids, respectively, in response to salt stress. Moreover, transcription factors (TFs) were also identified, which play an important role in salt stress tolerance of water dropwort. The finding of this study will be helpful for crop improvement under salt stress.

De novo transcriptome sequencing and analysis of salt-, alkali-, and drought-responsive genes in Sophora alopecuroides

Background

Salinity, alkalinity, and drought stress are the main abiotic stress factors affecting plant growth and development. Sophora alopecuroides L., a perennial leguminous herb in the genus Sophora, is a highly salt-tolerant sand-fixing pioneer species distributed mostly in Western Asia and northwestern China. Few studies have assessed responses to abiotic stress in S. alopecuroides. The transcriptome of the genes that confer stress-tolerance in this species has not previously been sequenced. Our objective was to sequence and analyze this transcriptome.

Results

Twelve cDNA libraries were constructed in triplicate from mRNA obtained from Sophora alopecuroides for the control and salt, alkali, and drought treatments. Using de novo assembly, 902,812 assembled unigenes were generated, with an average length of 294 bp. Based on similarity searches, 545,615 (60.43%) had at least one significant match in the Nr, Nt, Pfam, KOG/COG, Swiss-Prot, and GO databases. In addition, 1673 differentially expressed genes (DEGs) were obtained from the salt treatment, 8142 from the alkali treatment, and 17,479 from the drought treatment. A total of 11,936 transcription factor genes from 82 transcription factor families were functionally annotated under salt, alkali, and drought stress, these include MYB, bZIP, NAC and WRKY family members. DEGs were involved in the hormone signal transduction pathway, biosynthesis of secondary metabolites and antioxidant enzymes; this suggests that these pathways or processes may be involved in tolerance towards salt, alkali, and drought stress in S. alopecuroides.

Conclusion

Our study first reported transcriptome reference sequence data in Sophora alopecuroides, a non-model plant without a reference genome. We determined digital expression profile and discovered a broad survey of unigenes associated with salt, alkali, and drought stress which provide genomic resources available for Sophora alopecuroides.

Selection and Validation of Reference Genes for Normalisation of Gene Expression in Glehnia littoralis

Glehnia littoralis is an important medicinal halophyte—the dried root of which is used as Chinese herbal medicine. However, the use, selection and stability of reference genes are rarely verified in studies of G. littoralis, which hampers investigation of its salt tolerance and metabolism. In this study, we selected 13 candidate reference genes from the transcriptome data of G. littoralis—serine/threonine-protein phosphatase PP2A (PP2A), polyubiquitin 10 (UBQ10), actin (ACT), elongation factor 1-α (EF1-α), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), α-tubulin (α-TUB), β-tubulin (β-TUB), polypyrimidine tract-binding protein 1 (PTBP1), expressed protein 1 (EXP1), expressed protein 2 (EXP2), TIP41-like (TIP41), SAND family (SAND), and cyclophilin 2 (CYP2), and used qRT-PCR to analyse their expression levels in roots of G. littoralis treated with NaCl, polyethylene glycol (PEG), abscisic acid (ABA), and methyl jasmonate (MeJA), as well as in various organs of G. littoralis. The ΔCt, geNorm, NormFinder, and BestKeeper algorithms were used to assess the expression stability of the candidate reference genes and the results were then used to generate a comprehensive rank list with the RankAggreg R package. The most stable reference genes for normalisation were EXP1 and PP2A in response to NaCl, EXP2 and PP2A in response to ABA, CYP2 and α-TUB in response to MeJA, and ACT and EXP1 in the PEG and the organ subsets. GAPDH, β-TUB, and UBQ10 exhibited low stability and so were unsuitable for normalisation. This study is the first systematic analysis of candidate reference genes in G. littoralis and will facilitate further investigation of normalisation of gene expression in G. littoralis.

RNA-Seq analysis of Clerodendrum inerme (L.) roots in response to salt stress

Background

Clerodendrum inerme (L.) Gaertn, a halophyte, usually grows on coastal beaches as an important mangrove plant. The salt-tolerant mechanisms and related genes of this species that respond to short-term salinity stress are unknown for us. The de novo transcriptome of C. inerme roots was analyzed using next-generation sequencing technology to identify genes involved in salt tolerance and to better understand the response mechanisms of C. inerme to salt stress.

Results

Illumina RNA-sequencing was performed on root samples treated with 400 mM NaCl for 0 h, 6 h, 24 h, and 72 h to investigate changes in C. inerme in response to salt stress. The de novo assembly identified 98,968 unigenes. Among these unigenes, 46,085 unigenes were annotated in the NCBI non-redundant protein sequences (NR) database, 34,756 sequences in the Swiss-Prot database and 43,113 unigenes in the evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG) database. 52 Gene Ontology (GO) terms and 31 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were matched to those unigenes. Most differentially expressed genes (DEGs) related to the GO terms “single-organism process”, “membrane” and “catalytic activity” were significantly enriched while numerous DEGs related to the plant hormone signal transduction pathway were also significantly enriched. The detection of relative expression levels of 9 candidate DEGs by qRT-PCR were basically consistent with fold changes in RNA sequencing analysis, demonstrating that transcriptome data can accurately reflect the response of C. inerme roots to salt stress.

Conclusions

This work revealed that the response of C. inerme roots to saline condition included significant alteration in response of the genes related to plant hormone signaling. Besides, our findings provide numerous salt-tolerant genes for further research to improve the salt tolerance of functional plants and will enhance research on salt-tolerant mechanisms of halophytes.