Identification and comparative analysis of the microRNA transcriptome in roots of two contrasting tobacco genotypes in response to cadmium stress

Mechanistic basis for mitigating drought tolerance by selenium application in tobacco (Nicotiana tabacum L.): a multi-omics approach

The lack of irrigation water in agricultural soils poses a significant constraint on global crop production. In-depth investigation into microRNAs (miRNAs) has been widely used to achieve a comprehensive understanding of plant defense mechanisms. However, there is limited knowledge on the association of miRNAs with drought tolerance in cigar tobacco. In this study, a hydroponic experiment was carried out to identify changes in plant physiological characteristics, miRNA expression and metabolite profile under drought stress, and examine the mitigating effects of selenium (Se) application. The shoot dry weight of drought-stressed plants was approximately half (50.3%) of that in non-stressed (control) conditions. However, plants supplied with Se attained 38.8% greater shoot dry weight as compared to plants with no Se supply under drought stress. Thirteen miRNAs were identified to be associated with drought tolerance. These included 7 known (such as nta-miR156b and nta-miR166a) and 6 novel miRNAs (such as novel-nta-miR156-5p and novel-nta-miR209-5p) with the target genes of squamosa promoter-binding-like protein 4 (SPL4), serine/threonine protein phosphatase 2A (PPP2A), cation/calcium exchanger 4-like (CCX4), extensin-1-like (EXT1) and reduced wall acetylation 2 (RWA2). Further investigation revealed that the expression levels of Ext1 and RWA2 were significantly decreased under drought stress but increased with Se addition. Moreover, key metabolites such as catechin and N-acetylneuraminic acid were identified, which may play a role in the regulation of drought tolerance. The integrated analysis of miRNA sequencing and metabolome highlighted the significance of the novel-nta-miR97-5p- LRR-RLK- catechin pathway in regulating drought tolerance. Our findings provide valuable insights into the molecular mechanisms underlying drought tolerance and Se-induced stress alleviation in cigar tobacco.

Noncoding RNAs and their roles in regulating the agronomic traits of crops

Molecular breeding is one of the most effective methods for improving the performance of crops. Understanding the genome features of crops, especially the physiological functions of individual genes, is of great importance to molecular breeding. Evidence has shown that genomes of both animals and plants transcribe numerous non-coding RNAs, which are involved in almost every aspect of development. In crops, an increasing number of studies have proven that non-coding RNAs are new genetic resources for regulating crop traits. In this review, we summarize the current knowledge of non-coding RNAs, which are potential crop trait regulators, and focus on the functions of long non-coding RNAs (lncRNAs) in determining crop grain yield, phased small-interfering RNAs (phasiRNAs) in regulating fertility, small interfering RNAs (siRNAs) and microRNAs (miRNAs) in facilitating plant immune response and disease resistance, and miRNAs mediating nutrient and metal stress. Finally, we also discuss the next-generation method for ncRNA application in crop domestication and breeding.

Genome-wide identification and analysis of the invertase gene family in tobacco (Nicotiana tabacum) reveals NtNINV10 participating the sugar metabolism

Sucrose (Suc) is directly associated with plant growth and development as well as tolerance to various stresses. Invertase (INV) enzymes played important role in sucrose metabolism by irreversibly catalyzing Suc degradation. However, genome-wide identification and function of individual members of the INV gene family in Nicotiana tabacum have not been conducted. In this report, 36 non-redundant NtINV family members were identified in Nicotiana tabacum including 20 alkaline/neutral INV genes (NtNINV1-20), 4 vacuolar INV genes (NtVINV1-4), and 12 cell wall INV isoforms (NtCWINV1-12). A comprehensive analysis based on the biochemical characteristics, the exon-intron structures, the chromosomal location and the evolutionary analysis revealed the conservation and the divergence of NtINVs. For the evolution of the NtINV gene, fragment duplication and purification selection were major factors. Besides, our analysis revealed that NtINV could be regulated by miRNAs and cis-regulatory elements of transcription factors associated with multiple stress responses. In addition, 3D structure analysis has provided evidence for the differentiation between the NINV and VINV. The expression patterns in diverse tissues and under various stresses were investigated, and qRT-PCR experiments were conducted to confirm the expression patterns. Results revealed that changes in NtNINV10 expression level were induced by leaf development, drought and salinity stresses. Further examination revealed that the NtNINV10-GFP fusion protein was located in the cell membrane. Furthermore, inhibition of the expression of NtNINV10 gene decreased the glucose and fructose in tobacco leaves. Overall, we have identified possible NtINV genes functioned in leaf development and tolerance to environmental stresses in tobacco. These findings provide a better understanding of the NtINV gene family and establish the basis for future research.

Comprehensive analysis of the carboxylesterase gene reveals that NtCXE22 regulates axillary bud growth through strigolactone metabolism in tobacco

Carboxylesterases (CXE) are a class of hydrolytic enzymes with α/β-folding domains that play a vital role in plant growth, development, stress response, and activation of herbicide-active substances. In this study, 49 Nicotiana tabacum L. CXE genes (NtCXEs) were identified using a sequence homology search. The basic characteristics, phylogenetic evolution, gene structure, subcellular location, promoter cis-elements, and gene expression patterns of the CXE family were systematically analyzed. RNA-seq data and quantitative real-time PCR showed that the expression level of CXEs was associated with various stressors and hormones; gene expression levels were significantly different among the eight tissues examined and at different developmental periods. As a new class of hormones, strigolactones (SLs) are released from the roots of plants and can control the germination of axillary buds.NtCXE7, NtCXE9, NtCXE22, and NtCXE24 were homologous to Arabidopsis SLs hydrolase AtCXE15, and changes in their expression levels were induced by topping and by GR24 (a synthetic analogue of strigolactone). Further examination revealed that NtCXE22-mutant (ntcxe22) plants generated by CRISPR-Cas9 technology had shorter bud outgrowth with lower SLs content. Validation of NtCXE22 was also performed in NtCCD8-OE plants (with fewer axillary buds) and in ntccd8 mutant plants (with more axillary buds). The results suggest that NtCXE22 may act as an efficient SLs hydrolase and affects axillary bud development, thereby providing a feasible method for manipulating endogenous SLs in crops and ornamental plants.

Genome-Wide Identification and Analysis of the Class III Peroxidase Gene Family in Tobacco (Nicotiana tabacum)

Class III peroxidases (PODs) are plant-specific enzymes that play significant roles in plant physiological processes and stress responses. However, a comprehensive analysis of the POD gene family in tobacco has not yet been conducted. In this study, 210 non-redundant POD gene members (NtPODs) were identified in tobacco (Nicotiana tabacum) and distributed unevenly throughout 24 tobacco chromosomes. Phylogenetic analysis clustered these genes into six subgroups (I-VI). Gene structure and motif analyses showed the structural and functional diversity among the subgroups. Segmental duplication and purifying selection were the main factors affecting NtPOD gene evolution. Our analyses also suggested that NtPODs might be regulated by miRNAs and cis-acting regulatory elements of transcription factors that are involved in various biological processes. In addition, the expression patterns in different tissues and under various stress treatments were investigated. The results showed that the majority of NtPODs had tissue-specific expression patterns and may be involved in many biotic and abiotic responses. qRT-PCR analyses of different tissues and stress treatments were performed to verify transcriptome patterns. Expression of a green fluorescent protein-NtPOD fusion confirmed the plasma membrane localization of NtPOD121 and NtPOD4. Furthermore, 3D structures provided evidences of membrane-bound peroxidase. These findings provide useful information to better understand the evolution of the NtPOD gene family and lay the foundation for further studies on POD gene function in tobacco.

Regulatory role of microRNAs (miRNAs) in the recent development of abiotic stress tolerance of plants

MicroRNAs (miRNAs) are a distinct groups of single-stranded non-coding, tiny regulatory RNAs approximately 20-24 nucleotides in length. miRNAs negatively influence gene expression at the post-transcriptional level and have evolved considerably in the development of abiotic stress tolerance in a number of model plants and economically important crop species. The present review aims to deliver the information on miRNA-mediated regulation of the expression of major genes or Transcription Factors (TFs), as well as genetic and regulatory pathways. Also, the information on adaptive mechanisms involved in plant abiotic stress responses, prediction, and validation of targets, computational tools, and databases available for plant miRNAs, specifically focus on their exploration for engineering abiotic stress tolerance in plants. The regulatory function of miRNAs in plant growth, development, and abiotic stresses consider in this review, which uses high-throughput sequencing (HTS) technologies to generate large-scale libraries of small RNAs (sRNAs) for conventional screening of known and novel abiotic stress-responsive miRNAs adds complexity to regulatory networks in plants. The discoveries of miRNA-mediated tolerance to multiple abiotic stresses, including salinity, drought, cold, heat stress, nutritional deficiency, UV-radiation, oxidative stress, hypoxia, and heavy metal toxicity, are highlighted and discussed in this review.

Genome-wide identification, phylogeny, and expression analysis of the bHLH gene family in tobacco (Nicotiana tabacum)

The basic helix-loop-helix (bHLH) is the second-largest TF family in plants that play important roles in plant growth, development, and stress responses. In this study, a total of 100 bHLHs were identified using Hidden Markov Model profiles in the Nicotiana tabacum genome, clustered into 15 major groups (I-XV) based on their conserved domains and phylogenetic relationships. Group VIII genes were found to be the most abundant, with 27 NtbHLH members. The expansion of NtbHLHs in the genome was due to segmental and tandem duplication. The purifying selection was found to have an important role in the evolution of NtHLHs. Subsequent qRT-PCR validation of five selected genes from transcriptome data revealed that NtbHLH3.1, NtbHLH3.2, NtbHLH24, NtbHLH50, and NtbHLH59.2 have higher expressions at 12 and 24 h in comparison to 0 h (control) of chilling stress. The validated results demonstrated that NtbHLH3.2 and NtbHLH24 genes have 3 and fivefold higher expression at 12 h and 2 and threefold higher expression at 24 h than control plant, shows high sensitivity towards chilling stress. Moreover, the co-expression of positively correlated genes of NtbHLH3.2 and NtbHLH24 confirmed their functional significance in chilling stress response. Therefore, suggesting the importance of NtbHLH3.2 and NtbHLH24 genes in exerting control over the chilling stress responses in tobacco.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01042-x.

The Relationship between Cadmium Toxicity and the Modulation of Epigenetic Traits in Plants

Elevated concentrations of heavy metals such as cadmium (Cd) have a negative impact on staple crop production due to their ability to elicit cytotoxic and genotoxic effects on plants. In order to understand the relationship between Cd stress and plants in an effort to improve Cd tolerance, studies have identified genetic mechanisms which could be important for conferring stress tolerance. In recent years epigenetic studies have garnered much attention and hold great potential in both improving the understanding of Cd stress in plants as well as revealing candidate mechanisms for future work. This review describes some of the main epigenetic mechanisms involved in Cd stress responses. We summarize recent literature and data pertaining to chromatin remodeling, DNA methylation, histone acetylation and miRNAs in order to understand the role these epigenetic traits play in cadmium tolerance. The review aims to provide the framework for future studies where these epigenetic traits may be used in plant breeding and molecular studies in order to improve Cd tolerance.

HvHOX9, a novel homeobox leucine zipper transcription factor, positively regulates aluminum tolerance in Tibetan wild barley

Aluminum (Al) toxicity is the primary limiting factor of crop production on acid soils. Tibetan wild barley germplasm is a valuable source of potential genes for breeding barley with acid and Al tolerance. We performed microRNA and RNA sequencing using wild (XZ16, Al-tolerant; XZ61, Al-sensitive) and cultivated (Dayton, Al-tolerant) barley. A novel homeobox-leucine zipper transcription factor, HvHOX9, was identified as a target gene of miR166b and functionally characterized. HvHOX9 was up-regulated by Al stress in XZ16 (but unchanged in XZ61 and Dayton) and was significantly induced only in root tip. Phylogenetic analysis showed that HvHOX9 is most closely related to wheat TaHOX9 and orthologues of HvHOX9 are present in the closest algal relatives of Zygnematophyceae. Barley stripe mosaic virus-induced gene silencing of HvHOX9 in XZ16 led to significantly increased Al sensitivity but did not affect its sensitivity to other metals and low pH. Disruption of HvHOX9 did not change Al concentration in the root cell sap, but led to more Al accumulation in root cell wall after Al exposure. Silencing of HvHOX9 decreased H+ influx after Al exposure. Our findings suggest that miR166b/HvHOX9 play a critical role in Al tolerance by decreasing root cell wall Al binding and increasing apoplastic pH for Al detoxification in the root.

Identification and expression analysis of conserved microRNAs during short and prolonged chromium stress in rice (Oryza sativa)

MicroRNAs (miRNAs) are one of the most critical epigenetic regulators of gene expression which modulate a spectrum of development and defence response processes in plants. Chromium (Cr) contamination in rice imposes a serious concern to human health as rice is used as staple food throughout the world. Although several studies have established the differential response of miRNAs in rice during heavy metal (arsenic, cadmium) and heat or cold stress, no report is available about the response of miRNAs during Cr stress. In the present study, we identified 512 and 568 known miRNAs from Cr treated and untreated samples, respectively. Expression analysis revealed that 13 conserved miRNAs (miR156, miR159, miR160, miR166, miR169, miR171, miR396, miR397, miR408, miR444, miR1883, miR2877, miR5072) depicted preferential up- or down-regulation (> 4-fold change; P value < 0.05). Target gene prediction of differentially expressed miRNAs and their functional annotation suggested the important role of miRNAs in defence and detoxification of Cr though ATP-binding cassette transporters (ABC transporters), transcription factors, heat shock proteins, auxin response, and metal ion transport. Real-time PCR analysis validated the differential expression of selected miRNAs and their putative target genes. In conclusion, our study identifies and predicts miRNA-mediated regulation of signalling pathway in rice during Cr stress.

Cyclic nucleotide gated channels (CNGCs) in plant signalling-Current knowledge and perspectives

Cell signaling is an evolutionarily conserved mechanism that responds and adapts to various internal and external factors. Generally, a signal is mediated by various signaling molecules and is transferred to a cascade of effector proteins. To date, there is significant evidence that cyclic nucleotides (cNMPs), e.g., adenosine 3',5'-cyclic monophosphate (cAMP) and guanosine 3',5'-cyclic monophosphate (cGMP), may represent important elements of many signaling pathways in plants. However, in contrast to the impressive progress made in understanding cyclic nucleotide signaling in mammalian hosts, only few studies have investigated this topic in plants. Existing evidence indicates that cNMPs participate in growth and developmental processes, as well as the response to various stresses. Once synthesized by adenylyl or guanylyl cyclases, these signals are transduced by acting through a number of cellular effectors. The regulatory effects of cNMPs in eukaryotes can be mediated via various downstream effector proteins, such as protein kinases, Exchange Protein directly Activated by cAMP (EPAC), and Cyclic Nucleotide-Gated ion Channels (CNGC). These proteins sense changes in intracellular cNMP levels and regulate numerous cellular responses. Moreover, the amplitude of cNMP levels and the duration of its signal in the cell is also governed by phosphodiesterases (PDEs), enzymes that are responsible for the breakdown of cNMPs. Data collected in recent years strongly suggest that cyclic nucleotide gated channels are the main cNMP effectors in plant cells. These channels are important cellular switches that transduce changes in intracellular concentrations of cyclic nucleotides into changes in membrane potential and ion concentrations. Structurally, these channels belong to the superfamily of pore-loop cation channels. In this review, we provide an overview of the molecular properties of CNGC structure, regulation and ion selectivity, and subcellular localization, as well as describing the signal transduction pathways in which these channels are involved. We will also summarize recent insights into the role of CNGC proteins in plant growth, development and response to stressors.

Genome-Wide Identification, Comprehensive Gene Feature, Evolution, and Expression Analysis of Plant Metal Tolerance Proteins in Tobacco Under Heavy Metal Toxicity

Plant metal tolerance proteins (MTPs) comprise a family of membrane divalent cation transporters that play essential roles in plant mineral nutrition maintenance and heavy metal stresses resistance. However, the evolutionary relationships and biological functions of MTP family in tobacco remain unclear. In the present study, 26, 13, and 12 MTPs in three main Nicotiana species (N. tabacum, N. sylvestris, and N. tomentosiformis) were identified and designated, respectively. The phylogenetic relationships, gene structures, chromosome distributions, conserved motifs, and domains of NtMTPs were systematic analyzed. According to the phylogenetic features, 26 NtMTPs were classified into three major substrate-specific groups that were Zn-cation diffusion facilitators (CDFs), Zn/Fe-CDFs, and Mn-CDFs, and seven primary groups (1, 5, 6, 7, 8, 9, and 12). All of the NtMTPs contained a modified signature sequence and the cation_efflux domain, whereas some of them also harbored the ZT_dimer. Evolutionary analysis showed that NtMTP family of N. tabacum originated from its parental genome of N. sylvestris and N. tomentosiformis, and further underwent gene loss and expanded via one segmental duplication event. Moreover, the prediction of cis-acting elements (CREs) and the microRNA target sites of NtMTP genes suggested the diverse and complex regulatory mechanisms that control NtMTPs gene expression. Expression profile analysis derived from transcriptome data and quantitative real-time reverse transcription-PCR (qRT-PCR) analysis showed that the tissue expression patterns of NtMTPs in the same group were similar but varied among groups. Besides, under heavy metal toxicity, NtMTP genes exhibited various responses in either tobacco leaves or roots. 19 and 15 NtMTPs were found to response to at least one metal ion treatment in leaves and roots, respectively. In addition, NtMTP8.1, NtMTP8.4, and NtMTP11.1 exhibited Mn transport abilities in yeast cells. These results provided a perspective on the evolution of MTP genes in tobacco and were helpful for further functional characterization of NtMTP genes.

Intrusive Growth of Phloem Fibers in Flax Stem: Integrated Analysis of miRNA and mRNA Expression Profiles

Phloem fibers are important elements of plant architecture and the target product of many fiber crops. A key stage in fiber development is intrusive elongation, the mechanisms of which are largely unknown. Integrated analysis of miRNA and mRNA expression profiles in intrusivelygrowing fibers obtained by laser microdissection from flax (Linum usitatissimum L.) stem revealed all 124 known flax miRNA from 23 gene families and the potential targets of differentially expressed miRNAs. A comparison of the expression between phloem fibers at different developmental stages, and parenchyma and xylem tissues demonstrated that members of miR159, miR166, miR167, miR319, miR396 families were down-regulated in intrusively growing fibers. Some putative target genes of these miRNA families, such as those putatively encoding growth-regulating factors, an argonaute family protein, and a homeobox-leucine zipper family protein were up-regulated in elongating fibers. miR160, miR169, miR390, and miR394 showed increased expression. Changes in the expression levels of miRNAs and their target genes did not match expectations for the majority of predicted target genes. Taken together, poorly understood intrusive fiber elongation, the key process of phloem fiber development, was characterized from a miRNA-target point of view, giving new insights into its regulation.

Evolutionary and expression analysis of CAMTA gene family in Nicotiana tabacum yielded insights into their origin, expansion and stress responses

Calmodulin-binding transcription activators (CAMTAs) represent the novel gene family of transcriptional regulators, which play important biological functions. Though, the first ever plant CAMTA gene was evidenced in Nicotiana tabacum in 2002. But, the systematic identification, origin and function of this gene family has not been performed due to the lack of reference genome information until now. Here, we identified 29 CAMTA genes in four Nicotiana species, including thirteen NtabCAMTAs, six NsylCAMTAs, and five NtomCAMTAs and NbenCAMTAs. These CAMTA families were classified into five phylogenetic groups (I-V), among which, the group-IV CAMTAs probably emerged the earliest. The NtabCAMTA family genes have diverse structures, and are randomly localized on five chromosomes and scaffolds. N. tabacum acquired 11 copies of homolog CAMATA genes from the parental genomes of N. tomentosiformis and N. sylvestris, followed by expansion through polyploidization and duplication. The NtabCAMTA genes were differentially expressed in different plant parts, and showed sensitivity towards different abiotic and biotic stresses. Co-expression network analysis revealed that some NtabCAMTA subunits interact with each other, and co-expressed. The current study is the first report presenting a comprehensive overview of Nicotiana CAMTA families, and opens a new avenue for the improvement of the cultivated tobacco.

Traversing the Links between Heavy Metal Stress and Plant Signaling

Plants confront multifarious environmental stresses widely divided into abiotic and biotic stresses, of which heavy metal stress represents one of the most damaging abiotic stresses. Heavy metals cause toxicity by targeting crucial molecules and vital processes in the plant cell. One of the approaches by which heavy metals act in plants is by over production of reactive oxygen species (ROS) either directly or indirectly. Plants act against such overdose of metal in the environment by boosting the defense responses like metal chelation, sequestration into vacuole, regulation of metal intake by transporters, and intensification of antioxidative mechanisms. This response shown by plants is the result of intricate signaling networks functioning in the cell in order to transmit the extracellular stimuli into an intracellular response. The crucial signaling components involved are calcium signaling, hormone signaling, and mitogen activated protein kinase (MAPK) signaling that are discussed in this review. Apart from signaling components other regulators like microRNAs and transcription factors also have a major contribution in regulating heavy metal stress. This review demonstrates the key role of MAPKs in synchronously controlling the other signaling components and regulators in metal stress. Further, attempts have been made to focus on metal transporters and chelators that are regulated by MAPK signaling.

Cold-responsive miRNAs and their target genes in the wild eggplant species Solanum aculeatissimum

Background

Low temperature is an important abiotic stress in plant growth and development, especially for thermophilic plants. Eggplants are thermophilic vegetables, although the molecular mechanism of their response to cold stress remains to be elucidated. MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs that play an essential role during plant development and stress responses. Although the role of many plant miRNAs in facilitating chilling tolerance has been verified, little is known about the mechanisms of eggplant chilling tolerance.

Results

Here, we used high-throughput sequencing to extract the miRNA and target genes expression profiles of Solanum aculeatissimum (S. aculeatissimum) under low temperature stress at different time periods(0 h, 2 h, 6 h, 12 h, 24 h). Differentially regulated miRNAs and their target genes were analyzed by comparing the small RNA (sRNA) and miRBase 20.0 databases using BLAST or BOWTIE, respectively. Fifty-six down-regulated miRNAs and 28 up-regulated miRNAs corresponding to 220 up-regulated mRNAs and 94 down-regulated mRNAs, respectively, were identified in S. aculeatissimum. Nine significant differentially expressed miRNAs and twelve mRNAs were identified by quantitative Real-time PCR and association analysis, and analyzed for their GO function enrichment and KEGG pathway association.

Conclusions

In summary, numerous conserved and novel miRNAs involved in the chilling response were identified using high-throughput sequencing, which provides a theoretical basis for the further study of low temperature stress-related miRNAs and the regulation of cold-tolerance mechanisms of eggplant at the miRNA level.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4341-y) contains supplementary material, which is available to authorized users.