Endogenous small-noncoding RNAs and their roles in chilling response and stress acclimation in Cassava

LncRNAs exert indispensable roles in orchestrating the interaction among diverse noncoding RNAs and enrich the regulatory network of plant growth and its adaptive environmental stress response

With the advent of advanced sequencing technologies, non-coding RNAs (ncRNAs) are increasingly pivotal and play highly regulated roles in the modulation of diverse aspects of plant growth and stress response. This includes a spectrum of ncRNA classes, ranging from small RNAs to long non-coding RNAs (lncRNAs). Notably, among these, lncRNAs emerge as significant and intricate components within the broader ncRNA regulatory networks. Here, we categorize ncRNAs based on their length and structure into small RNAs, medium-sized ncRNAs, lncRNAs, and circle RNAs. Furthermore, the review delves into the detailed biosynthesis and origin of these ncRNAs. Subsequently, we emphasize the diverse regulatory mechanisms employed by lncRNAs that are located at various gene regions of coding genes, embodying promoters, 5'UTRs, introns, exons, and 3'UTR regions. Furthermore, we elucidate these regulatory modes through one or two concrete examples. Besides, lncRNAs have emerged as novel central components that participate in phase separation processes. Moreover, we illustrate the coordinated regulatory mechanisms among lncRNAs, miRNAs, and siRNAs with a particular emphasis on the central role of lncRNAs in serving as sponges, precursors, spliceosome, stabilization, scaffolds, or interaction factors to bridge interactions with other ncRNAs. The review also sheds light on the intriguing possibility that some ncRNAs may encode functional micropeptides. Therefore, the review underscores the emergent roles of ncRNAs as potent regulatory factors that significantly enrich the regulatory network governing plant growth, development, and responses to environmental stimuli. There are yet-to-be-discovered roles of ncRNAs waiting for us to explore.

Non-coding RNAs (ncRNAs) in plant: Master regulators for adapting to extreme temperature conditions

Unusual daily temperature fluctuations caused by climate change and climate variability adversely impact agricultural crop production. Since plants are immobile and constantly receive external environmental signals, such as extreme high (heat) and low (cold) temperatures, they have developed complex molecular regulatory mechanisms to cope with stressful situations to sustain their natural growth and development. Among these mechanisms, non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), small-interfering RNAs (siRNAs), and long-non-coding RNAs (lncRNAs), play a significant role in enhancing heat and cold stress tolerance. This review explores the pivotal findings related to miRNAs, siRNAs, and lncRNAs, elucidating how they functionally regulate plant adaptation to extreme temperatures. In addition, this review addresses the challenges associated with uncovering these non-coding RNAs and understanding their roles in orchestrating heat and cold tolerance in plants.

High throughput sRNA sequencing revealed gene regulatory role mediated by pathogen-derived small RNAs during Sri Lankan Cassava Mosaic Virus infection in Cassava

Small RNA (sRNA) mediated gene regulation during Sri Lankan Cassava Mosaic Virus (SLCMV) infection was studied from the Indian Cassava Cultivar H226. Our study generated high throughput sRNA dataset of 23.64 million reads from the control and SLCMV infected H226 leaf libraries. mes-miR9386 was detected as the most prominent miRNA expressed in control and infected leaf. Among the differentially expressed miRNAs, mes-miR156, mes- miR395 and mes-miR535a/b showed significant down regulation in the infected leaf. Genome-wide analysis of the three small RNA profiles revealed critical role of virus-derived small RNAs (vsRNAs) from the infected leaf tissues of H226. The vsRNAs were mapped to the bipartite SLCMV genome and high expression of siRNAs generated from the virus genomic region encoding AV1/AV2 genes in the infected leaf pointed towards the susceptibility of H226 cultivars to SLCMV. Furthermore, the sRNA reads mapped to the antisense strand of the SLCMV ORFs was higher than the sense strand. These vsRNAs were potential to target key host genes involved in virus interaction such as aldehyde dehydrogenase, ADP-ribosylation factor1 and ARF1-like GTP-binding proteins. The sRNAome-assisted analysis also revealed the origin of virus-encoded miRNAs from the SLCMV genome in the infected leaf. These virus-derived miRNAs were predicted to have hair-pin like secondary structures, and have different isoforms. Moreover, our study revealed that the pathogen sRNAs play a critical role in the infection process in H226 plants.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03494-2.

Antisense Transcription in Plants: A Systematic Review and an Update on cis-NATs of Sugarcane

Initially, natural antisense transcripts (NATs, natRNAs, or asRNAs) were considered repressors; however, their functions in gene regulation are diverse. Positive, negative, or neutral correlations to the cognate gene expression have been noted. Although the first studies were published about 50 years ago, there is still much to be investigated regarding antisense transcripts in plants. A systematic review of scientific publications available in the Web of Science databases was conducted to contextualize how the studying of antisense transcripts has been addressed. Studies were classified considering three categories: “Natural antisense” (208), artificial antisense used in “Genetic Engineering” (797), or “Natural antisense and Genetic Engineering”-related publications (96). A similar string was used for a systematic search in the NCBI Gene database. Of the 1132 antisense sequences found for plants, only 0.8% were cited in PubMed and had antisense information confirmed. This value was the lowest when compared to fungi (2.9%), bacteria (2.3%), and mice (54.1%). Finally, we present an update for the cis-NATs identified in Saccharum spp. Of the 1413 antisense transcripts found in different experiments, 25 showed concordant expressions, 22 were discordant, 1264 did not correlate with the cognate genes, and 102 presented variable results depending on the experiment.

Small RNA and Degradome Deep Sequencing Reveals the Roles of microRNAs in Peanut (Arachis hypogaea L.) Cold Response

Cold stress is a major environmental factor that affects plant growth and development, as well as fruit postharvest life and quality. MicroRNAs (miRNAs) are a class of non-coding small RNAs that play crucial roles in various abiotic stresses. Peanuts (Arachis hypogaea L.), one of the most important grain legumes and source of edible oils and proteins, are cultivated in the semi-arid tropical and subtropical regions of the world. To date, there has been no report on the role of miRNAs in the response to cold stress in cultivated peanuts. In this study, we profiled cold-responsive miRNAs in peanuts using deep sequencing in cold-sensitive (WQL20) alongside a tolerant variety (WQL30). A total of 407 known miRNAs and 143 novel peanut-specific miRNAs were identified. The expression of selected known and novel miRNAs was validated by northern blotting and six known cold-responsive miRNAs were revealed. Degradome sequencing identified six cold-responsive miRNAs that regulate 12 target genes. The correlative expression patterns of several miRNAs and their target genes were further validated using qRT-PCR. Our data showed that miR160-ARF, miR482-WDRL, miR2118-DR, miR396-GRF, miR162-DCL, miR1511-SRF, and miR1511-SPIRAL1 modules may mediate cold stress responses. Transient expression analysis in Nicotiana benthamiana found that miR160, miR482, and miR2118 may play positive roles, and miR396, miR162, and miR1511 play negative roles in the regulation of peanut cold tolerance. Our results provide a foundation for understanding miRNA-dependent cold stress response in peanuts. The characterized correlations between miRNAs and their response to cold stress could serve as markers in breeding programs or tools for improving cold tolerance of peanuts.

Auxin response factors in plant adaptation to drought and salinity stress

Salinity and drought stresses affect plant growth worldwide and limit crop production. Auxin is crucial in regulating plants' salinity and drought stress adaptative response. As a chemical messenger, auxin influences gene expression through a family of functionally distinct transcription factors, the DNA-binding AUXIN RESPONSE FACTORS (ARFs). Various studies have revealed the important roles of ARFs in regulating drought and salinity stress responses in plants. Different ARFs regulate soluble sugar content, promote root development, and maintain chlorophyll content under drought and saline stress conditions to help plants adapt to these stresses. The functional characterization of ARFs pertaining to the regulation of drought and salinity stress responses is still in its infancy. Interestingly, the small RNA-mediated post-transcriptional regulation of ARF expression has been shown to influence plant responses to both stresses. The current knowledge on the diverse roles of ARFs in conferring specificity to auxin-mediated drought and salinity stress responses has not been reviewed to date. In this review, we summarize the recent research concerning the role of ARFs in response to drought and salinity stresses: gene expression patterns, functional characterization, and post-transcriptional regulation under drought and salinity stresses. We have also reviewed the modulation of ARF expression by other molecular regulators in the context of drought and salt stress signaling.

TRANS-ACTING SIRNA3-derived short interfering RNAs confer cleavage of mRNAs in rice

Plant TRANS-ACTING SIRNA3 (TAS3)-derived short interfering RNAs (siRNAs) include tasiR-AUXIN RESPONSE FACTORs (ARFs), which are functionally conserved in targeting ARF genes, and a set of non-tasiR-ARF siRNAs, which have rarely been studied. In this study, TAS3 siRNAs were systematically characterized in rice (Oryza sativa). Small RNA sequencing results showed that an overwhelming majority of TAS3 siRNAs belong to the non-tasiR-ARF group, while tasiR-ARFs occupy a diminutive fraction. Phylogenetic analysis of TAS3 genes across dicot and monocot plants revealed that the siRNA-generating regions were highly conserved in grass species, especially in the Oryzoideae. Target genes were identified for not only tasiR-ARFs but also non-tasiR-ARF siRNAs by analyzing rice Parallel Analysis of RNA Ends datasets, and some of these siRNA-target interactions were experimentally confirmed using tas3 mutants generated by genome editing. Consistent with the de-repression of target genes, phenotypic alterations were observed for mutants in three TAS3 loci in comparison to wild-type rice. The regulatory role of ribosomes in the TAS3 siRNA-target interactions was further revealed by the fact that TAS3 siRNA-mediated target cleavage, in particular tasiR-ARFs targeting ARF2/3/14/15, occurred extensively in rice polysome samples. Altogether, our study sheds light into TAS3 genes in plants and expands our knowledge about rice TAS3 siRNA-target interactions.

Identification and Analysis of Small Interfering RNAs Associated With Heat Stress in Flowering Chinese Cabbage Using High-Throughput Sequencing

Endogenous small interfering RNAs (siRNAs) are substantial gene regulators in eukaryotes and play key functions in plant development and stress tolerance. Among environmental factors, heat is serious abiotic stress that severely influences the productivity and quality of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). However, how siRNAs are involved in regulating gene expression during heat stress is not fully understood in flowering Chinese cabbage. Combining bioinformatical and next-generation sequencing approaches, we identified heat-responsive siRNAs in four small RNA libraries of flowering Chinese cabbage using leaves collected at 0, 1, 6, and 12 h after a 38°C heat-stress treatment; 536, 816, and 829 siRNAs exhibited substantial differential expression at 1, 6, and 12 h, respectively. Seventy-five upregulated and 69 downregulated differentially expressed siRNAs (DE-siRNAs) were common for the three time points of heat stress. We identified 795 target genes of DE-siRNAs, including serine/threonine-protein kinase SRK2I, CTR1-like, disease resistance protein RML1A-like, and RPP1, which may play a role in regulating heat tolerance. Gene ontology showed that predictive targets of DE-siRNAs may have key roles in the positive regulation of biological processes, organismal processes, responses to temperature stimulus, signaling, and growth and development. These novel results contribute to further understanding how siRNAs modulate the expression of their target genes to control heat tolerance in flowering Chinese cabbage.

The Dynamism of Transposon Methylation for Plant Development and Stress Adaptation

Plant development processes are regulated by epigenetic alterations that shape nuclear structure, gene expression, and phenotypic plasticity; these alterations can provide the plant with protection from environmental stresses. During plant growth and development, these processes play a significant role in regulating gene expression to remodel chromatin structure. These epigenetic alterations are mainly regulated by transposable elements (TEs) whose abundance in plant genomes results in their interaction with genomes. Thus, TEs are the main source of epigenetic changes and form a substantial part of the plant genome. Furthermore, TEs can be activated under stress conditions, and activated elements cause mutagenic effects and substantial genetic variability. This introduces novel gene functions and structural variation in the insertion sites and primarily contributes to epigenetic modifications. Altogether, these modifications indirectly or directly provide the ability to withstand environmental stresses. In recent years, many studies have shown that TE methylation plays a major role in the evolution of the plant genome through epigenetic process that regulate gene imprinting, thereby upholding genome stability. The induced genetic rearrangements and insertions of mobile genetic elements in regions of active euchromatin contribute to genome alteration, leading to genomic stress. These TE-mediated epigenetic modifications lead to phenotypic diversity, genetic variation, and environmental stress tolerance. Thus, TE methylation is essential for plant evolution and stress adaptation, and TEs hold a relevant military position in the plant genome. High-throughput techniques have greatly advanced the understanding of TE-mediated gene expression and its associations with genome methylation and suggest that controlled mobilization of TEs could be used for crop breeding. However, development application in this area has been limited, and an integrated view of TE function and subsequent processes is lacking. In this review, we explore the enormous diversity and likely functions of the TE repertoire in adaptive evolution and discuss some recent examples of how TEs impact gene expression in plant development and stress adaptation.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

Plant Non-Coding RNAs: Origin, Biogenesis, Mode of Action and Their Roles in Abiotic Stress

As sessile species, plants have to deal with the rapidly changing environment. In response to these environmental conditions, plants employ a plethora of response mechanisms that provide broad phenotypic plasticity to allow the fine-tuning of the external cues related reactions. Molecular biology has been transformed by the major breakthroughs in high-throughput transcriptome sequencing and expression analysis using next-generation sequencing (NGS) technologies. These innovations have provided substantial progress in the identification of genomic regions as well as underlying basis influencing transcriptional and post-transcriptional regulation of abiotic stress response. Non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), short interfering RNAs (siRNAs), and long non-coding RNAs (lncRNAs), have emerged as essential regulators of plants abiotic stress response. However, shared traits in the biogenesis of ncRNAs and the coordinated cross-talk among ncRNAs mechanisms contribute to the complexity of these molecules and might play an essential part in regulating stress responses. Herein, we highlight the current knowledge of plant microRNAs, siRNAs, and lncRNAs, focusing on their origin, biogenesis, modes of action, and fundamental roles in plant response to abiotic stresses.

Genomic and Transcriptomic Analysis Identified Novel Putative Cassava lncRNAs Involved in Cold and Drought Stress

Long non-coding RNAs (lncRNAs) play important roles in the regulation of complex cellular processes, including transcriptional and post-transcriptional regulation of gene expression relevant for development and stress response, among others. Compared to other important crops, there is limited knowledge of cassava lncRNAs and their roles in abiotic stress adaptation. In this study, we performed a genome-wide study of ncRNAs in cassava, integrating genomics- and transcriptomics-based approaches. In total, 56,840 putative ncRNAs were identified, and approximately half the number were verified using expression data or previously known ncRNAs. Among these were 2229 potential novel lncRNA transcripts with unmatched sequences, 250 of which were differentially expressed in cold or drought conditions, relative to controls. We showed that lncRNAs might be involved in post-transcriptional regulation of stress-induced transcription factors (TFs) such as zinc-finger, WRKY, and nuclear factor Y gene families. These findings deepened our knowledge of cassava lncRNAs and shed light on their stress-responsive roles.

Integrated annotations and analyses of small RNA-producing loci from 47 diverse plants

Plant endogenous small RNAs (sRNAs) are important regulators of gene expression. There are two broad categories of plant sRNAs: microRNAs (miRNAs) and endogenous short interfering RNAs (siRNAs). MicroRNA loci are relatively well-annotated but compose only a small minority of the total sRNA pool; siRNA locus annotations have lagged far behind. Here, we used a large data set of published and newly generated sRNA sequencing data (1333 sRNA-seq libraries containing more than 20 billion reads) and a uniform bioinformatic pipeline to produce comprehensive sRNA locus annotations of 47 diverse plants, yielding more than 2.7 million sRNA loci. The two most numerous classes of siRNA loci produced mainly 24- and 21-nucleotide (nt) siRNAs, respectively. Most often, 24-nt-dominated siRNA loci occurred in intergenic regions, especially at the 5′-flanking regions of protein-coding genes. In contrast, 21-nt-dominated siRNA loci were most often derived from double-stranded RNA precursors copied from spliced mRNAs. Genic 21-nt-dominated loci were especially common from disease resistance genes, including from a large number of monocots. Individual siRNA sequences of all types showed very little conservation across species, whereas mature miRNAs were more likely to be conserved. We developed a web server where our data and several search and analysis tools are freely accessible.

Deep-sequencing of Solanum commersonii small RNA libraries reveals riboregulators involved in cold stress response

Among wild species used in potato breeding, Solanum commersonii displays the highest tolerance to low temperatures under both acclimated (ACC) and non-acclimated (NACC) conditions. It is also the first wild potato relative with a known whole genome sequence. Recent studies have shown that abiotic stresses induce changes in the expression of many small non-coding RNA (sncRNA). We determined the small non-coding RNA (sncRNAome) of two clones of S. commersonii contrasting in their cold response phenotypes via smRNAseq. Differential analysis provided evidence that expression of several miRNAs changed in response to cold stress conditions. Conserved miR408a and miR408b changed their expression under NACC conditions, whereas miR156 and miR169 were differentially expressed only under ACC conditions. We also report changes in tasiRNA and secondary siRNA expression under both stress conditions. Our results reveal possible roles of sncRNA in the regulatory networks associated with tolerance to low temperatures and provide useful information for a more strategic use of genomic resources in potato breeding.

A perspective on crassulacean acid metabolism photosynthesis evolution of orchids on different continents: Dendrobium as a case study

Members of the Orchidaceae, one of the largest families of flowering plants, evolved the crassulacean acid metabolism (CAM) photosynthesis strategy. It is thought that CAM triggers adaptive radiation into new niche spaces, yet very little is known about its origin and diversification on different continents. Here, we assess the prevalence of CAM in Dendrobium, which is one of the largest genera of flowering plants and found in a wide range of environments, from the high altitudes of the Himalayas to relatively arid habitats in Australia. Based on phylogenetic time trees, we estimated that CAM, as determined by δ 13C values less negative than -20.0‰, evolved independently at least eight times in Dendrobium. The oldest lineage appeared in the Asian clade during the middle Miocene, indicating the origin of CAM was associated with a pronounced climatic cooling that followed a period of aridity. Divergence of the four CAM lineages in the Asian clade appeared to be earlier than divergence of those in the Australasian clade. However, CAM species in the Asian clade are much less diverse (25.6%) than those in the Australasian clade (57.9%). These findings shed new light on CAM evolutionary history and the aridity levels of the paleoclimate on different continents.

Cassava AGPase genes and their encoded proteins are different from those of other plants

Cassava AGPase and AGPase genes have some unique characteristics. ADP-glucose pyrophosphorylase (AGPase) is a rate-limiting enzyme for starch synthesis. In this study, cassava AGPase genes (MeAGP) were analyzed based on six cultivars and one wild species. A total of seven MeAGPs was identified, including four encoding AGPase large subunits (MeAGPLs 1, 2, 3 and 4) and three encoding AGPase small subunits (MeAGPSs 1, 2 and 3). The copy number of MeAGPs varied in cassava germplasm materials. There were 14 introns for MeAGPLs 1, 2 and 3, 13 introns for MeAGPL4, and 8 introns for other three MeAGPSs. Multiple conservative amino acid sequence motifs were found in the MeAGPs. There were differences in amino acids at binding sites of substrates and regulators among different MeAGP subunits and between MeAGPs and a potato AGPase small subunit (1YP2:B). MeAGPs were all located in chloroplasts. MeAGP expression was not only associated with gene copy number and types/combinations, regions and levels of the DNA methylation but was also affected by environmental factors with the involvement of various transcription factors in multiple regulation networks and in various cis-elements in the gene promoter regions. The MeAGP activity also changed with environmental conditions and had potential differences among the subunits. Taken together, MeAGPs differ in number from those of Arabidopsis, potato, maize, banana, sweet potato, and tomato.

Genome-Wide Identification of Putative MicroRNAs in Cassava (Manihot esculenta Crantz) and Their Functional Landscape in Cellular Regulation

MicroRNAs are small noncoding RNAs, involved in the regulation of many cellular processes in plants. Hundreds of miRNAs have been identified in cassava by various techniques, yet these identifications were constrained by a lack of miRNA templates and the narrow range of conditions in transcriptome study. In this research, we conducted genome-wide analysis identification, whereby miRNAs from cassava genome were thoroughly screened using bioinformatics approach independent of predefined templates and studied conditions. Our work provided a catalog of putative mature miRNAs and explored the landscape of miRNAome in cassava. These putative miRNAs were validated using statistical analysis as well as available cassava expression data. We showed that the crowded locations of cassava miRNAs are consistent with other plants and animals and hypothesized to have the same evolutionary origin. At least 10 conserved miRNAs were identified in cassava based on the comparative study of miRNA conservation. Finally, investigation of miRNAs and target gene relationships enabled us to envisage the complexities of cellular regulatory systems modulated at posttranscriptional level.

Identification and characterization of drought-responsive CC-type glutaredoxins from cassava cultivars reveals their involvement in ABA signalling

BACKGROUND

CC-type glutaredoxins (GRXs) are plant-specific glutaredoxin, play regulatory roles in response of biotic and abiotic stress. However, it is not clear whether the CC-type GRXs are involve in drought response in cassava (Manihot esculenta), an important tropical tuber root crop.

RESULTS

Herein, genome-wide analysis identified 18 CC-type GRXs in the cassava genome, of which six (namely MeGRXC3, C4, C7, C14, C15, and C18) were induced by drought stress in leaves of two cassava cultivars Argentina 7 (Arg7) and South China 124 (SC124). Exogenous abscisic acid (ABA) application induced the expression of all the six CC-type GRXs in leaves of both Arg7 and SC124 plants. Overexpression of MeGRXC15 in Arabidopsis (Col-0) increases tolerance of ABA on the sealed agar plates, but results in drought hypersensitivity in soil-grown plants. The results of microarray assays show that MeGRXC15 overexpression affected the expression of a set of transcription factors which involve in stress response, ABA, and JA/ET signalling pathway. The results of protein interaction analysis show that MeGRXC15 can interact with TGA5 from Arabidopsis and MeTGA074 from cassava.

CONCLUSIONS

CC-type glutaredoxins play regulatory roles in cassava response to drought possibly through ABA signalling pathway.

High-throughput sequencing of small RNAs revealed the diversified cold-responsive pathways during cold stress in the wild banana (Musa itinerans)

BACKGROUND

Cold stress is one of the most severe abiotic stresses affecting the banana production. Although some miRNAs have been identified, little is known about the role of miRNAs in response to cold stress in banana, and up to date, there is no report about the role of miRNAs in the response to cold stress in the plants of the cultivated or wild bananas.

RESULT

Here, a cold-resistant line wild banana (Musa itinerans) from China was used to profile the cold-responsive miRNAs by RNA-seq during cold stress. Totally, 265 known mature miRNAs and 41 novel miRNAs were obtained. Cluster analysis of differentially expressed (DE) miRNAs indicated that some miRNAs were specific for chilling or 0 °C treated responses, and most of them were reported to be cold-responsive; however, some were seldom reported to be cold-responsive in response to cold stress, e.g., miR395, miR408, miR172, suggesting that they maybe play key roles in response to cold stress. The GO and KEGG pathway enrichment analysis of DE miRNAs targets indicated that there existed diversified cold-responsive pathways, and miR172 was found likely to play a central coordinating role in response to cold stress, especially in the regulation of CK2 and the circadian rhythm. Finally, qPCR assays indicated the related targets were negatively regulated by the tested DE miRNAs during cold stress in the wild banana.

CONCLUSIONS

In this study, the profiling of miRNAs by RNA-seq in response to cold stress in the plants of the wild banana (Musa itinerans) was reported for the first time. The results showed that there existed diversified cold-responsive pathways, which provided insight into the roles of miRNAs during cold stress, and would be helpful for alleviating cold stress and cold-resistant breeding in bananas.

Auxin Response Factors (ARFs) are potential mediators of auxin action in tomato response to biotic and abiotic stress (Solanum lycopersicum)

Survival biomass production and crop yield are heavily constrained by a wide range of environmental stresses. Several phytohormones among which abscisic acid (ABA), ethylene and salicylic acid (SA) are known to mediate plant responses to these stresses. By contrast, the role of the plant hormone auxin in stress responses remains so far poorly studied. Auxin controls many aspects of plant growth and development, and Auxin Response Factors play a key role in the transcriptional activation or repression of auxin-responsive genes through direct binding to their promoters. As a mean to gain more insight on auxin involvement in a set of biotic and abiotic stress responses in tomato, the present study uncovers the expression pattern of SlARF genes in tomato plants subjected to biotic and abiotic stresses. In silico mining of the RNAseq data available through the public TomExpress web platform, identified several SlARFs as responsive to various pathogen infections induced by bacteria and viruses. Accordingly, sequence analysis revealed that 5' regulatory regions of these SlARFs are enriched in biotic and abiotic stress-responsive cis-elements. Moreover, quantitative qPCR expression analysis revealed that many SlARFs were differentially expressed in tomato leaves and roots under salt, drought and flooding stress conditions. Further pointing to the putative role of SlARFs in stress responses, quantitative qPCR expression studies identified some miRNA precursors as potentially involved in the regulation of their SlARF target genes in roots exposed to salt and drought stresses. These data suggest an active regulation of SlARFs at the post-transcriptional level under stress conditions. Based on the substantial change in the transcript accumulation of several SlARF genes, the data presented in this work strongly support the involvement of auxin in stress responses thus enabling to identify a set of candidate SlARFs as potential mediators of biotic and abiotic stress responses.

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.

MicroRNA-like RNAs from the same miRNA precursors play a role in cassava chilling responses

MicroRNAs (miRNAs) are known to play important roles in various cellular processes and stress responses. MiRNAs can be identified by analyzing reads from high-throughput deep sequencing. The reads realigned to miRNA precursors besides canonical miRNAs were initially considered as sequencing noise and ignored from further analysis. Here we reported a small-RNA species of phased and half-phased miRNA-like RNAs different from canonical miRNAs from cassava miRNA precursors detected under four distinct chilling conditions. They can form abundant multiple small RNAs arranged along precursors in a tandem and phased or half-phased fashion. Some of these miRNA-like RNAs were experimentally confirmed by re-amplification and re-sequencing, and have a similar qRT-PCR detection ratio as their cognate canonical miRNAs. The target genes of those phased and half-phased miRNA-like RNAs function in process of cell growth metabolism and play roles in protein kinase. Half-phased miR171d.3 was confirmed to have cleavage activities on its target gene P-glycoprotein 11, a broad substrate efflux pump across cellular membranes, which is thought to provide protection for tropical cassava during sharp temperature decease. Our results showed that the RNAs from miRNA precursors are miRNA-like small RNAs that are viable negative gene regulators and may have potential functions in cassava chilling responses.

High throughput deep sequencing reveals the important roles of microRNAs during sweetpotato storage at chilling temperature

Sweetpotato (Impomoea batatas L.) is a globally important economic food crop with a potential of becoming a bioenergy and pharmaceutical crop. Thus, studying the molecular mechanism of tuberous root development and storage is very important. However, not too much progress has been made in this field. In this study, we employed the next generation high-throughput deep sequencing technology to sequence all small RNAs and degradome of sweetpotato for systematically investigating sweetpotato response to chilling stress during storage. A total of 190 known microRNAs (miRNAs) and 191 novel miRNAs were identified, and 428 transcripts were targeted by 184 identified miRNAs. More importantly, we identified 26 miRNAs differentially expressed between chilling stress and control conditions. The expression of these miRNAs and their targets was also confirmed by qRT-PCR. Integrated analysis of small RNAs and degradome sequencing reveals that miRNA-mediated SA signaling, ABA-dependent, and ROS response pathways are involved in sweetpotato root response to chilling stress during storage.

Identification of drought-responsive miRNAs and physiological characterization of tea plant (Camellia sinensis L.) under drought stress

BACKGROUND

Drought stress is one of the major natural challenges in the main tea-producing regions of China. The tea plant (Camellia sinensis) is a traditional beverage plant whose growth status directly affects tea quality. Recent studies have revealed that microRNAs (miRNAs) play key functions in plant growth and development. Although some miRNAs have been identified in C. sinensis, little is known about their roles in the drought stress response of tea plants.

RESULTS

Physiological characterization of Camellia sinensis 'Tieguanyin' under drought stress showed that the malondialdehyde concentration and electrical conductivity of leaves of drought-stressed plants increased when the chlorophyll concentration decreased under severe drought stress. We sequenced four small-RNA (sRNA) libraries constructed from leaves of plants subjected to four different treatments, normal water supply (CK); mild drought stress (T1); moderate drought stress (T2) and severe drought stress (T3). A total of 299 known mature miRNA sequences and 46 novel miRNAs were identified. Gene Ontology enrichment analysis revealed that most of the differentially expressed-miRNA target genes were related to regulation of transcription. Kyoto Encyclopedia of Genes and Genomes analysis revealed that the most highly enriched pathways under drought stress were D-alanine metabolism, sulfur metabolism, and mineral absorption pathways. Real-time quantitative PCR (qPCR) was used to validate the expression patterns of 21 miRNAs (2 up-regulated and 19 down-regulated under drought stress). The observed co-regulation of the miR166 family and their targets ATHB-14-like and ATHB-15-like indicate the presence of negative feedback regulation in miRNA pathways.

CONCLUSIONS

Analyses of drought-responsive miRNAs in tea plants showed that most of differentially expressed-miRNA target genes were related to regulation of transcription. The results of study revealed that the expressions of phase-specific miRNAs vary with morphological, physiological, and biochemical changes. These findings will be useful for research on drought resistance and provide insights into the mechanisms of drought adaptation and resistance in C. sinensis.

Small RNAs in regulating temperature stress response in plants

Due to global climate change, temperature stress has become one of the primary causes of crop losses worldwide. Much progress has been made in unraveling the complex stress response mechanisms in plants, particularly in the identification of temperature stress responsive protein-coding genes. Recently discovered microRNAs (miRNAs) and endogenous small-interfering RNAs (siRN As) have also been demonstrated as important players in plant temperature stress response. Using high-throughput sequencing, many small RNAs, especially miRNAs, have been identified to be triggered by cold or heat. Subsequently, several studies have shown an important functional role for these small RNAs in cold or heat tolerance. These findings greatly broaden our understanding of endogenous small RNAs in plant stress response control. Here, we highlight new findings regarding the roles of miRNAs and siRNAs in plant temperature stress response and acclimation. We also review the current understanding of the regulatory mechanisms of small RNAs in temperature stress response, and explore the outlook for the use of these small RNAs in molecular breeding for improvement of temperature stress tolerance in plants.

Non-coding RNAs and Their Roles in Stress Response in Plants

Eukaryotic genomes encode thousands of non-coding RNAs (ncRNAs), which play crucial roles in transcriptional and post-transcriptional regulation of gene expression. Accumulating evidence indicates that ncRNAs, especially microRNAs (miRNAs) and long ncRNAs (lncRNAs), have emerged as key regulatory molecules in plant stress responses. In this review, we have summarized the current progress on the understanding of plant miRNA and lncRNA identification, characteristics, bioinformatics tools, and resources, and provided examples of mechanisms of miRNA- and lncRNA-mediated plant stress tolerance.

Genome-wide characterization and expression analysis enables identification of abiotic stress-responsive MYB transcription factors in cassava (Manihot esculenta)

The myeloblastosis (MYB) transcription factor superfamily is the largest transcription factor family in plants, playing different roles during stress response. However, abiotic stress-responsive MYB transcription factors have not been systematically studied in cassava (Manihot esculenta), an important tropical tuber root crop. In this study, we used a genome-wide transcriptome analysis to predict 299 putative MeMYB genes in the cassava genome. Under drought and cold stresses, many MeMYB genes exhibited different expression patterns in cassava leaves, indicating that these genes might play a role in abiotic stress responses. We found that several stress-responsive MeMYB genes responded to abscisic acid (ABA) in cassava leaves. We characterize four MeMYBs, namely MeMYB1, MeMYB2, MeMYB4, and MeMYB9, as R2R3-MYB transcription factors. Furthermore, RNAi-driven repression of MeMYB2 resulted in drought and cold tolerance in transgenic cassava. Gene expression assays in wild-type and MeMYB2-RNAi cassava plants revealed that MeMYB2 may affect other MeMYBs as well as MeWRKYs under drought and cold stress, suggesting crosstalk between MYB and WRKY family genes under stress conditions in cassava.

SRNAome and degradome sequencing analysis reveals specific regulation of sRNA in response to chilling injury in tomato fruit

Plant genomes encode diverse small RNA classes that function in distinct gene-silencing pathways. To elucidate the intricate regulation of microRNAs (miRNAs) and endogenous small-interfering RNAs (siRNAs) in response to chilling injury in tomato fruit, the deep sequencing and bioinformatic methods were combined to decipher the small RNAs landscape in the control and chilling-injured groups. Except for the known miRNAs and ta-siRNAs, 85 novel miRNAs and 5 ta-siRNAs members belonging to 3 TAS families (TAS5, TAS9 and TAS10) were identified, 34 putative phased small RNAs and 740 cis/trans-natural antisense small-interfering RNAs (nat-siRNAs) were also found in our results which enriched the tomato small RNAs repository. A large number of genes targeted by those miRNAs and siRNAs were predicted to be involved in the chilling injury responsive process and five of them were verified via degradome sequencing. Based on the above results, a regulatory model that comprehensively reveals the relationships between the small RNAs and their targets was set up. This work provides a foundation for further study of the regulation of miRNAs and siRNAs in the plant in response to chilling injury.

Genome-wide identification and functional prediction of cold and/or drought-responsive lncRNAs in cassava

Cold and drought stresses seriously affect cassava (Manihot esculenta) plant growth and yield. Recently, long noncoding RNAs (lncRNAs) have emerged as key regulators of diverse cellular processes in mammals and plants. To date, no systematic screening of lncRNAs under abiotic stress and their regulatory roles in cassava has been reported. In this study, we present the first reference catalog of 682 high-confidence lncRNAs based on analysis of strand-specific RNA-seq data from cassava shoot apices and young leaves under cold, drought stress and control conditions. Among them, 16 lncRNAs were identified as putative target mimics of cassava known miRNAs. Additionally, by comparing with small RNA-seq data, we found 42 lncNATs and sense gene pairs can generate nat-siRNAs. We identified 318 lncRNAs responsive to cold and/or drought stress, which were typically co-expressed concordantly or discordantly with their neighboring genes. Trans-regulatory network analysis suggested that many lncRNAs were associated with hormone signal transduction, secondary metabolites biosynthesis, and sucrose metabolism pathway. The study provides an opportunity for future computational and experimental studies to uncover the functions of lncRNAs in cassava.

Arabidopsis thaliana siRNA biogenesis mutants have the lower frequency of homologous recombination

Small interfering RNAs (siRNAs) are involved in the regulation of plant development and response to stress. We have previously shown that mutants impaired in Dicer-like 2 (DCL2), DCL3 and DCL4, RDR2, RDR6 and NPRD1 are partially impaired in their response to stress and dcl2 and dcl3 plants are also impaired in transgenerational response to stress, including changes in homologous recombination frequency (HRF). Here, we have analyzed genome stability of dcl2, dcl3, dcl4, dcl2 dcl3, dcl2 dcl3 dcl4 and rdr6 mutants by measuring the non-induced and the stress-induced recombination frequency. We found that all mutants had the lower spontaneous HRF. The analysis of strand breaks showed that all tested Arabidopsis mutants had a higher level of spontaneous strand breaks, suggesting that the lower HRF is not due to the unusually low level of breaks. Exposure to methyl methane sulfonate (MMS) resulted in an increase in the level of strand breaks in wild-type plants and a decrease in mutants. All mutants had the higher methylation of cytosines at CpG sites under non-induced conditions. Exposure to MMS resulted in a decrease in methylation level in wild-type plants and an increase in methylation in all dcl mutants. The expression of several DNA repair genes was altered in dcl4 plants under non-induced and induced conditions. Our data suggest that siRNA biogenesis may be essential for the maintenance of the genome stability and stress response in Arabidopsis.

Unveiling the Micronome of Cassava (Manihot esculenta Crantz)

MicroRNAs (miRNAs) are an important class of endogenous non-coding single-stranded small RNAs (21-24 nt in length), which serve as post-transcriptional negative regulators of gene expression in plants. Despite the economic importance of Manihot esculenta Crantz (cassava) only 153 putative cassava miRNAs (from multiple germplasm) are available to date in miRBase (Version 21), and identification of a number of miRNAs from the cassava EST database have been limited to comparisons with Arabidopsis. In this study, mature sequences of all known plant miRNAs were used as a query for homologous searches against cassava EST and GSS databases, and additional identification of novel and conserved miRNAs were gleaned from next generation sequencing (NGS) of two cassava landraces (T200 from southern Africa and TME3 from West Africa) at three different stages post explant transplantation and acclimatization. EST and GSS derived data revealed 259 and 32 miRNAs in cassava, and one of the miRNA families (miR2118) from previous studies has not been reported in cassava. NGS data collectively displayed expression of 289 conserved miRNAs in leaf tissue, of which 230 had not been reported previously. Of the 289 conserved miRNAs identified in T200 and TME3, 208 were isomiRs. Thirty-nine novel cassava-specific miRNAs of low abundance, belonging to 29 families, were identified. Thirty-eight (98.6%) of the putative new miRNAs identified by NGS have not been previously reported in cassava. Several miRNA targets were identified in T200 and TME3, highlighting differential temporal miRNA expression between the two cassava landraces. This study contributes to the expanding knowledge base of the micronome of this important crop.

Integrated RNA-Seq and sRNA-Seq Analysis Identifies Chilling and Freezing Responsive Key Molecular Players and Pathways in Tea Plant (Camellia sinensis)

Tea [Camellia sinensis (L) O. Kuntze, Theaceae] is one of the most popular non-alcoholic beverages worldwide. Cold stress is one of the most severe abiotic stresses that limit tea plants’ growth, survival and geographical distribution. However, the genetic regulatory network and signaling pathways involved in cold stress responses in tea plants remain unearthed. Using RNA-Seq, DGE and sRNA-Seq technologies, we performed an integrative analysis of miRNA and mRNA expression profiling and their regulatory network of tea plants under chilling (4℃) and freezing (-5℃) stress. Differentially expressed (DE) miRNA and mRNA profiles were obtained based on fold change analysis, miRNAs and target mRNAs were found to show both coherent and incoherent relationships in the regulatory network. Furthermore, we compared several key pathways (e.g., ‘Photosynthesis’), GO terms (e.g., ‘response to karrikin’) and transcriptional factors (TFs, e.g., DREB1b/CBF1) which were identified as involved in the early chilling and/or freezing response of tea plants. Intriguingly, we found that karrikins, a new group of plant growth regulators, and β-primeverosidase (BPR), a key enzyme functionally relevant with the formation of tea aroma might play an important role in both early chilling and freezing response of tea plants. Quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) analysis further confirmed the results from RNA-Seq and sRNA-Seq analysis. This is the first study to simultaneously profile the expression patterns of both miRNAs and mRNAs on a genome-wide scale to elucidate the molecular mechanisms of early responses of tea plants to cold stress. In addition to gaining a deeper insight into the cold resistant characteristics of tea plants, we provide a good case study to analyse mRNA/miRNA expression and profiling of non-model plant species using next-generation sequencing technology.

Potential functions of microRNAs in starch metabolism and development revealed by miRNA transcriptome profiling of cassava cultivars and their wild progenitor

BACKGROUND

MicroRNAs (miRNAs) are small (approximately 21 nucleotide) non-coding RNAs that are key post-transcriptional gene regulators in eukaryotic organisms. More than 100 cassava miRNAs have been identified in a conservation analysis and a repertoire of cassava miRNAs have also been characterised by next-generation sequencing (NGS) in recent studies. Here, using NGS, we profiled small non-coding RNAs and mRNA genes in two cassava cultivars and their wild progenitor to identify and characterise miRNAs that are potentially involved in plant growth and starch biosynthesis.

RESULTS

Six small RNA and six mRNA libraries from leaves and roots of the two cultivars, KU50 and Arg7, and their wild progenitor, W14, were subjected to NGS. Analysis of the sequencing data revealed 29 conserved miRNA families and 33 new miRNA families. Together, these miRNAs potentially targeted a total of 360 putative target genes. Whereas 16 miRNA families were highly expressed in cultivar leaves, another 13 miRNA families were highly expressed in storage roots of cultivars. Co-expression analysis revealed that the expression level of some targets had negative relationship with their corresponding miRNAs in storage roots and leaves; these targets included MYB33, ARF10, GRF1, RD19, APL2, NF-YA3 and SPL2, which are known to be involved in plant development, starch biosynthesis and response to environmental stimuli.

CONCLUSION

The identified miRNAs, target mRNAs and target gene ontology annotation all shed light on the possible functions of miRNAs in Manihot species. The differential expression of miRNAs between cultivars and their wild progenitor, together with our analysis of GO annotation and confirmation of miRNA: target pairs, might provide insight into know the differences between wild progenitor and cultivated cassava.