Profiling of drought-responsive microRNA and mRNA in tomato using high-throughput sequencing

Effect of long-term deficit irrigation on tomato and goji berry quality: from fruit composition to in vitro bioaccessibility of carotenoids

Drought is a persistent challenge for horticulture, affecting various aspects of fruit development and ultimately fruit quality, but the effect on nutritional value has been under-investigated. Here, fruit quality was studied on six tomato genotypes and one goji cultivar under deficit irrigation (DI), from fruit composition to in vitro bioaccessibility of carotenoids. For both species, DI concentrated most health-related metabolites in fresh fruit. On a dry mass basis, DI increased total phenolic and sugar concentration, but had a negative or insignificant impact on fruit ascorbic acid, organic acid, and alcohol-insoluble matter contents. DI also reduced total carotenoids content in tomato (-18.7% on average), especially β-carotene (-32%), but not in goji berry DW (+15.5% and +19.6%, respectively). DI reduced the overall in vitro bioaccessibility of carotenoids to varying degrees depending on the compound and plant species. Consequently, mixed micelles produced by digestion of fruits subjected to DI contained either the same or lesser quantities of carotenoids, even though fresh fruits could contain similar or higher quantities. Thus, DI effects on fruit composition were species and genotype dependent, but an increase in the metabolite concentration did not necessarily translate into greater bioaccessibility potentially due to interactions with the fruit matrix.

Controlling water deficiency as an abiotic stress factor to improve tomato nutritional and flavour quality

Water deficit (WD) irrigation techniques to improve water use efficiency have been rapidly developed. However, the effect of WD irrigation on tomato quality has not been sufficiently studied. Here, we investigated the effects of varying water irrigation levels [T1-T4: 80%, 65%, 55%, and 45% of maximum field moisture capacity (FMC)] and full irrigation (CK: 90% of maximum FMC) on tomato fruits from the mature-green to red-ripening stages, to compare the nutritional and flavour qualities of the resulting tomatoes. The proline, aspartic, malic, citric, and ascorbic acid contents increased, phenylalanine and glutamic acid contents decreased, and the total amino and organic acid contents increased by 18.91% and 26.12%, respectively, in T2-treated fruits. Furthermore, the T2-treated fruits exhibited higher K and P contents alongside improved characteristic aromas. These findings provide novel insights for further improvements in tomato quality while also developing water-saving irrigation techniques.

Candidate regulators of drought stress in tomato revealed by comparative transcriptomic and proteomic analyses

Drought is among the most common abiotic constraints of crop growth, development, and productivity. Integrating different omics approaches offers a possibility for deciphering the metabolic pathways and fundamental mechanisms involved in abiotic stress tolerance. Here, we explored the transcriptional and post-transcriptional changes in drought-stressed tomato plants using transcriptomic and proteomic profiles to determine the molecular dynamics of tomato drought stress responses. We identified 22467 genes and 5507 proteins, among which the expression of 3765 genes and 294 proteins was significantly changed under drought stress. Furthermore, the differentially expressed genes (DEGs) and differentially abundant proteins (DAPs) showed a good correlation (0.743). The results indicated that integrating different omics approaches is promising in exploring the multilayered regulatory mechanisms of plant drought resistance. Gene ontology (GO) and pathway analysis identified several GO terms and pathways related to stress resistance, including response to stress, abiotic stimulus, and oxidative stress. The plant hormone abscisic acid (ABA) plays pivotal roles in response to drought stress, ABA-response element binding factor (AREB) is a key positive regulator of ABA signaling. Moreover, our analysis indicated that drought stress increased the abscisic acid (ABA) content, which activated AREB1 expression to regulate the expression of TAS14, GSH-Px-1, and Hsp, ultimately improving tomato drought resistance. In addition, the yeast one-hybrid assay demonstrated that the AREB1 could bind the Hsp promoter to activate Hsp expression. Thus, this study involved a full-scale analysis of gene and protein expression in drought-stressed tomato, deepening the understanding of the regulatory mechanisms of the essential drought-tolerance genes in tomato.

The miRNA-mRNA Regulatory Modules of Pinus massoniana Lamb. in Response to Drought Stress

Masson pine (Pinus massoniana Lamb.) is a major fast-growing woody tree species and pioneer species for afforestation in barren sites in southern China. However, the regulatory mechanism of gene expression in P. massoniana under drought remains unclear. To uncover candidate microRNAs, their expression profiles, and microRNA-mRNA interactions, small RNA-seq was used to investigate the transcriptome from seedling roots under drought and rewatering in P. massoniana. A total of 421 plant microRNAs were identified. Pairwise differential expression analysis between treatment and control groups unveiled 134, 156, and 96 differential expressed microRNAs at three stages. These constitute 248 unique microRNAs, which were subsequently categorized into six clusters based on their expression profiles. Degradome sequencing revealed that these 248 differentially expressed microRNAs targeted 2069 genes. Gene Ontology enrichment analysis suggested that these target genes were related to translational and posttranslational regulation, cell wall modification, and reactive oxygen species scavenging. miRNAs such as miR482, miR398, miR11571, miR396, miR166, miRN88, and miRN74, along with their target genes annotated as F-box/kelch-repeat protein, 60S ribosomal protein, copper-zinc superoxide dismutase, luminal-binding protein, S-adenosylmethionine synthase, and Early Responsive to Dehydration Stress may play critical roles in drought response. This study provides insights into microRNA responsive to drought and rewatering in Masson pine and advances the understanding of drought tolerance mechanisms in Pinus.

Comparative Genome-Wide Analysis of MicroRNAs and Their Target Genes in Roots of Contrasting Indica Rice Cultivars under Reproductive-Stage Drought

Recurrent occurrence of drought stress in varying intensity has become a common phenomenon in the present era of global climate change, which not only causes severe yield losses but also challenges the cultivation of rice. This raises serious concerns for sustainable food production and global food security. The root of a plant is primarily responsible to perceive drought stress and acquire sufficient water for the survival/optimal growth of the plant under extreme climatic conditions. Earlier studies reported the involvement/important roles of microRNAs (miRNAs) in plants' responses to environmental/abiotic stresses. A number (738) of miRNAs is known to be expressed in different tissues under varying environmental conditions in rice, but our understanding of the role, mode of action, and target genes of the miRNAs are still elusive. Using contrasting rice [IR-64 (reproductive-stage drought sensitive) and N-22 (drought-tolerant)] cultivars, imposed with terminal (reproductive-stage) drought stress, we demonstrate differential expression of 270 known and 91 novel miRNAs in roots of the contrasting rice cultivars in response to the stress. Among the known miRNAs, osamiR812, osamiR166, osamiR156, osamiR167, and osamiR396 were the most differentially expressed miRNAs between the rice cultivars. In the root of N-22, 18 known and 12 novel miRNAs were observed to be exclusively expressed, while only two known (zero novels) miRNAs were exclusively expressed in the roots of IR-64. The majority of the target gene(s) of the miRNAs were drought-responsive transcription factors playing important roles in flower, grain development, auxin signaling, root development, and phytohormone-crosstalk. The novel miRNAs identified in this study may serve as good candidates for the genetic improvement of rice for terminal drought stress towards developing climate-smart rice for sustainable food production.

Transcriptomic analysis reveals the gene regulatory networks involved in leaf and root response to osmotic stress in tomato

Introduction

Tomato (Solanum lycopersicum L.) is a major horticultural crop that is cultivated worldwide and is characteristic of the Mediterranean agricultural system. It represents a key component of the diet of billion people and an important source of vitamins and carotenoids. Tomato cultivation in open field often experiences drought episodes, leading to severe yield losses, since most modern cultivars are sensitive to water deficit. Water stress leads to changes in the expression of stress-responsive genes in different plant tissues, and transcriptomics can support the identification of genes and pathways regulating this response.

Methods

Here, we performed a transcriptomic analysis of two tomato genotypes, M82 and Tondo, in response to a PEG-mediated osmotic treatment. The analysis was conducted separately on leaves and roots to characterize the specific response of these two organs.

Results

A total of 6,267 differentially expressed transcripts related to stress response was detected. The construction of gene co-expression networks defined the molecular pathways of the common and specific responses of leaf and root. The common response was characterized by ABA-dependent and ABA-independent signaling pathways, and by the interconnection between ABA and JA signaling. The root-specific response concerned genes involved in cell wall metabolism and remodeling, whereas the leaf-specific response was principally related to leaf senescence and ethylene signaling. The transcription factors representing the hubs of these regulatory networks were identified. Some of them have not yet been characterized and can represent novel candidates for tolerance.

Discussion

This work shed new light on the regulatory networks occurring in tomato leaf and root under osmotic stress and set the base for an in-depth characterization of novel stress-related genes that may represent potential candidates for improving tolerance to abiotic stress in tomato.

Genes related to cell wall metabolisms are targeted by miRNAs in immature tomato fruits under drought stress

The metabolism of tomato fruits changes when plants experience drought stress. In this study, we investigated changes in microRNA (miRNA) abundance and detected 32 miRNAs whose expression changes in fruit. The candidate target genes for each miRNA were predicted from the differentially expressed genes identified by transcriptome analysis at the same fruit maturation stage. The predicted targeted genes were related to cell wall metabolisms, response to pathogens, and plant hormones. Among these, we focused on cell wall metabolism-related genes and performed a dual luciferase assay to assess the targeting of their mRNAs by their predicted miRNA. As a result, sly-miR10532 and sly-miR7981e suppress the expression of mRNAs of galacturonosyltransferase-10 like encoding the main enzyme of pectin biosynthesis, while sly-miR171b-5p targets β-1,3-glucosidase mRNAs involved in glucan degradation. These results will allow the systematic characterization of miRNA and their target genes in the tomato fruit under drought stress conditions.

Characterization of Tomato Brown Rugose Fruit Virus (ToBRFV) Detected in Czech Republic

Tomato is the most consumed vegetable in the world. The tomato brown rugose fruit virus (ToBRFV) is an important destructive virus that damages tomatoes and peppers with significant economic impact. The detection and characterization of this important viral pathogen were evaluated at the molecular and morphological level. The viral isolate was purified and inoculated on tomato and pepper plants. Small RNAs were sequenced in both plants and the profiles were compared. The complete genome of the isolate was obtained, and microRNA (miRNA) profiles were unveiled by small RNA sequencing. Symptoms caused by the isolate were also described and the morphology of the isolate was observed by transmission electron microscopy. Our results contribute to further understanding of the role of miRNAs in ToBRFV pathogenesis, which may be crucial for understanding disease symptom development in tomatoes and peppers.

Comparison of Tomato Transcriptomic Profiles Reveals Overlapping Patterns in Abiotic and Biotic Stress Responses

Until a few years ago, many studies focused on the transcriptomic response to single stresses. However, tomato cultivations are often constrained by a wide range of biotic and abiotic stress that can occur singularly or in combination, and several genes can be involved in the defensive mechanism response. Therefore, we analyzed and compared the transcriptomic responses of resistant and susceptible genotypes to seven biotic stresses (Cladosporium fulvum, Phytophthora infestans, Pseudomonas syringae, Ralstonia solanacearum, Sclerotinia sclerotiorum, Tomato spotted wilt virus (TSWV) and Tuta absoluta) and five abiotic stresses (drought, salinity, low temperatures, and oxidative stress) to identify genes involved in response to multiple stressors. With this approach, we found genes encoding for TFs, phytohormones, or participating in signaling and cell wall metabolic processes, participating in defense against various biotic and abiotic stress. Moreover, a total of 1474 DEGs were commonly found between biotic and abiotic stress. Among these, 67 DEGs were involved in response to at least four different stresses. In particular, we found RLKs, MAPKs, Fasciclin-like arabinogalactans (FLAs), glycosyltransferases, genes involved in the auxin, ET, and JA pathways, MYBs, bZIPs, WRKYs and ERFs genes. Detected genes responsive to multiple stress might be further investigated with biotechnological approaches to effectively improve plant tolerance in the field.

Recent allopolyploidy alters Spartina microRNA expression in response to xenobiotic-induced stress

Environmental contamination by xenobiotics represents a major threat for natural ecosystems and public health. In response, xenobiotic detoxification is a fundamental trait of organisms for developmental plasticity and stress tolerance, but the underlying molecular mechanisms remain poorly understood in plants. To decipher this process, we explored the consequences of allopolyploidy on xenobiotic tolerance in the genus Spartina Schreb. Specifically, we focused on microRNAs (miRNAs) owing to their central function in the regulation of gene expression patterns, including responses to stress. Small RNA-Seq was conducted on the parents S. alterniflora and S. maritima, their F1 hybrid S. x townsendii and the allopolyploid S. anglica under phenanthrene-induced stress (phe), a model Polycyclic Aromatic Hydrocarbon (PAH) compound. Differentially expressed miRNAs in response to phe were specifically identified within species. In complement, the respective impacts of hybridization and genome doubling were detected, through changes in miRNA expression patterns between S. x townsendii, S. anglica and the parents. The results support the impact of allopolyploidy in miRNA-guided regulation of plant response to phe. In total, we identified 17 phe-responsive miRNAs in Spartina among up-regulated MIR156 and down-regulated MIR159. We also describe novel phe-responsive miRNAs as putative Spartina-specific gene expression regulators in response to stress. Functional validation using Arabidopsis (L.) Heynh. T-DNA lines inserted in homologous MIR genes was performed, and the divergence of phe-responsive miRNA regulatory networks between Arabidopsis and Spartina was discussed.

miRNAs and lncRNAs in tomato: Roles in biotic and abiotic stress responses

Plants are continuously exposed to various biotic and abiotic stresses in the natural environment. To cope with these stresses, they have evolved a multitude of defenses mechanisms. With the rapid development of genome sequencing technologies, a large number of non-coding RNA (ncRNAs) have been identified in tomato, like microRNAs (miRNAs) and long non-coding RNAs (lncRNAs). Recently, more and more evidence indicates that many ncRNAs are involved in plant response to biotic and abiotic stresses in tomato. In this review, we summarize recent updates on the regulatory roles of ncRNAs in tomato abiotic/biotic responses, including abiotic (high temperature, drought, cold, salinization, etc.) and biotic (bacteria, fungi, viruses, insects, etc.) stresses. Understanding the molecular mechanisms mediated by ncRNAs in response to these stresses will help us to clarify the future directions for ncRNA research and resistance breeding in tomato.

Integrative analysis of microRNAs and mRNAs reveals the regulatory networks of triterpenoid saponin metabolism in Soapberry (Sapindus mukorossi Gaertn.)

Triterpenoid saponin are important secondary metabolites and bioactive constituents of soapberry (Sapindus mukorossi Gaertn.) and are widely used in medicine and toiletry products. However, little is known about the roles of miRNAs in the regulation of triterpenoid saponin biosynthesis in soapberry. In this study, a total of 3036 miRNAs were identified, of which 1372 miRNAs were differentially expressed at different stages of pericarp development. Important KEGG pathways, such as terpenoid backbone biosynthesis, sesquiterpenoid and triterpenoid biosynthesis, and basal transcription factors were highlighted, as well the roles of some key miRNAs, such as ath-miR5021, han-miR3630-3p, and ppe-miR858, which may play important roles in regulating triterpenoid saponin biosynthesis. In addition, 58 miRNAs might participate in saponin biosynthesis pathways by predicting the targets of those miRNAs to 53 saponin biosynthesis structural genes. And 75 miRNAs were identified to potentially play vital role in saponin accumulation by targeting transcript factor genes, bHLH, bZIP, ERF, MYB, and WRKY, respectively, which are candidate regulatory genes in the pathway of saponin biosynthesis. The results of weighted gene coexpression network analysis (WGCNA) suggested that two saponin-specific miRNA modules and 10 hub miRNAs may participate in saponin biosynthesis. Furthermore, multiple miRNA-mRNA regulatory networks potentially involved in saponin biosynthesis were generated, e.g., ath-miR5021-SmIDI2/SmGPS5/SmbAS1/SmCYP71D-3/SmUGT74G-2, han-miR3630-3p-SmCYP71A-14/SmbHLH54/SmMYB135/SmWRKY32, and ppe-miR858-SmMYB5/SmMYB32. qRT-PCR analysis validated the expression patterns of nine miRNAs and 12 corresponding target genes. This study represents the first comprehensive analysis of miRNAs in soapberry and lays the foundation for further understanding of miRNA-based regulation in triterpenoid saponin biosynthesis.

Spatiotemporal dynamics of the tomato fruit transcriptome under prolonged water stress

Water availability influences all aspects of plant growth and development; however, most studies of plant responses to drought have focused on vegetative organs, notably roots and leaves. Far less is known about the molecular bases of drought acclimation responses in fruits, which are complex organs with distinct tissue types. To obtain a more comprehensive picture of the molecular mechanisms governing fruit development under drought, we profiled the transcriptomes of a spectrum of fruit tissues from tomato (Solanum lycopersicum), spanning early growth through ripening and collected from plants grown under varying intensities of water stress. In addition, we compared transcriptional changes in fruit with those in leaves to highlight different and conserved transcriptome signatures in vegetative and reproductive organs. We observed extensive and diverse genetic reprogramming in different fruit tissues and leaves, each associated with a unique response to drought acclimation. These included major transcriptional shifts in the placenta of growing fruit and in the seeds of ripe fruit related to cell growth and epigenetic regulation, respectively. Changes in metabolic and hormonal pathways, such as those related to starch, carotenoids, jasmonic acid, and ethylene metabolism, were associated with distinct fruit tissues and developmental stages. Gene coexpression network analysis provided further insights into the tissue-specific regulation of distinct responses to water stress. Our data highlight the spatiotemporal specificity of drought responses in tomato fruit and indicate known and unrevealed molecular regulatory mechanisms involved in drought acclimation, during both vegetative and reproductive stages of development.

Transcriptomic and splicing changes underlying tomato responses to combined water and nutrient stress

Tomato is a horticultural crop of high economic and nutritional value. Suboptimal environmental conditions, such as limited water and nutrient availability, cause severe yield reductions. Thus, selection of genotypes requiring lower inputs is a goal for the tomato breeding sector. We screened 10 tomato varieties exposed to water deficit, low nitrate or a combination of both. Biometric, physiological and molecular analyses revealed different stress responses among genotypes, identifying T270 as severely affected, and T250 as tolerant to the stresses applied. Investigation of transcriptome changes caused by combined stress in roots and leaves of these two genotypes yielded a low number of differentially expressed genes (DEGs) in T250 compared to T270, suggesting that T250 tailors changes in gene expression to efficiently respond to combined stress. By contrast, the susceptible tomato activated approximately one thousand and two thousand genes in leaves and roots respectively, indicating a more generalized stress response in this genotype. In particular, developmental and stress-related genes were differentially expressed, such as hormone responsive factors and transcription factors. Analysis of differential alternative splicing (DAS) events showed that combined stress greatly affects the splicing landscape in both genotypes, highlighting the important role of AS in stress response mechanisms. In particular, several stress and growth-related genes as well as transcription and splicing factors were differentially spliced in both tissues. Taken together, these results reveal important insights into the transcriptional and post-transcriptional mechanisms regulating tomato adaptation to growth under reduced water and nitrogen inputs.

Deeper genomic insights into tomato CLE genes repertoire identify new active peptides

Background

In eukaryotes, cell-to-cell communication relies on the activity of small signaling peptides. In plant genomes, many hundreds of genes encode for such short peptide signals. However, only few of them are functionally characterized and due to the small gene size and high sequence variability, the comprehensive identification of such peptide-encoded genes is challenging. The CLAVATA3 (CLV3)/EMBRYO SURROUNDING REGION-RELATED (CLE) gene family encodes for short peptides that have a role in plant meristem maintenance, vascular patterning and responses to environment. The full repertoire of CLE genes and the role of CLE signaling in tomato (Solanum lycopersicum)- one of the most important crop plants- has not yet been fully studied.

Results

By using a combined approach, we performed a genome-wide identification of CLE genes using the current tomato genome version SL 4.0. We identified 52 SlCLE genes, including 37 new non annotated before. By analyzing publicly available RNAseq datasets we could confirm the expression of 28 new SlCLE genes. We found that SlCLEs are often expressed in a tissue-, organ- or condition-specific manner. Our analysis shows an interesting gene diversification within the SlCLE family that seems to be a result of gene duplication events. Finally, we could show a biological activity of selected SlCLE peptides in the root growth arrest that was SlCLV2-dependent.

Conclusions

Our improved combined approach revealed 37 new SlCLE genes. These findings are crucial for better understanding of the CLE signaling in tomato. Our phylogenetic analysis pinpoints the closest homologs of Arabidopsis CLE genes in tomato genome and can give a hint about the function of newly identified SlCLEs. The strategy described here can be used to identify more precisely additional short genes in plant genomes. Finally, our work suggests that the mechanism of root-active CLE peptide perception is conserved between Arabidopsis and tomato. In conclusion, our work paves the way to further research on the CLE-dependent circuits modulating tomato development and physiological responses.

Comprehensive Genome-Wide Analysis and Expression Pattern Profiling of PLATZ Gene Family Members in Solanum Lycopersicum L. under Multiple Abiotic Stresses

PLATZ (plant AT-rich sequence and zinc-binding) family proteins with two conserved zinc-dependent DNA-binding motifs are transcription factors specific to the plant kingdom. The functions of PLATZ proteins in growth, development, and adaptation to multiple abiotic stresses have been investigated in various plant species, but their role in tomato has not been explored yet. In the present work, 20 non-redundant Solanum lycopersicum PLATZ (SlPLATZ) genes with three segmentally duplicated gene pairs and four tandemly duplicated gene pairs were identified on eight tomato chromosomes. The comparative modeling and gene ontology (GO) annotations of tomato PLATZ proteins indicated their probable roles in defense response, transcriptional regulation, and protein metabolic processes as well as their binding affinity for various ligands, including nucleic acids, peptides, and zinc. SlPLATZ10 and SlPLATZ17 were only expressed in 1 cm fruits and flowers, respectively, indicating their preferential involvement in the development of these organs. The expression of SlPLATZ1, SlPLATZ12, and SlPLATZ19 was up- or down-regulated following exposure to various abiotic stresses, whereas that of SlPLATZ11 was induced under temperature stresses (i.e., cold and heat stress), revealing their probable function in the abiotic stress tolerance of tomato. Weighted gene co-expression network analysis corroborated the aforementioned findings by spotlighting the co-expression of several stress-associated genes with SlPLATZ genes. Confocal fluorescence microscopy revealed the localization of SlPLATZ−GFP fusion proteins in the nucleus, hinting at their functions as transcription factors. These findings provide a foundation for a better understanding of the structure and function of PLATZ genes and should assist in the selection of potential candidate genes involved in the development and abiotic stress adaptation in tomato.

Tomato MicroRNAs and Their Functions

MicroRNAs (miRNAs) define an essential class of non-coding small RNAs that function as posttranscriptional modulators of gene expression. They are coded by MIR genes, several hundreds of which exist in the genomes of Arabidopsis and rice model plants. The functional analysis of Arabidopsis and rice miRNAs indicate that their miRNAs regulate a wide range of processes including development, reproduction, metabolism, and stress. Tomato serves as a major model crop for the study of fleshy fruit development and ripening but until recently, information on the identity of its MIR genes and their coded miRNAs was limited and occasionally contradictory. As a result, the majority of tomato miRNAs remained uncharacterized. Recently, a comprehensive annotation of tomato MIR genes has been carried out by several labs and us. In this review, we curate and organize the resulting partially overlapping MIR annotations into an exhaustive and non-redundant atlas of tomato MIR genes. There are 538 candidate and validated MIR genes in the atlas, of which, 169, 18, and 351 code for highly conserved, Solanaceae-specific, and tomato-specific miRNAs, respectively. Furthermore, a critical review of functional studies on tomato miRNAs is presented, highlighting validated and possible functions, creating a useful resource for future tomato miRNA research.

Identification of microRNAs responsive to arbuscular mycorrhizal fungi in Panicum virgatum (switchgrass)

BACKGROUND

MicroRNAs (miRNAs) are important post-transcriptional regulators involved in the control of a range of processes, including symbiotic interactions in plants. MiRNA involvement in arbuscular mycorrhizae (AM) symbiosis has been mainly studied in model species, and our study is the first to analyze global miRNA expression in the roots of AM colonized switchgrass (Panicum virgatum), an emerging biofuel feedstock. AM symbiosis helps plants gain mineral nutrition from the soil and may enhance switchgrass biomass production on marginal lands. Our goals were to identify miRNAs and their corresponding target genes that are controlling AM symbiosis in switchgrass.

RESULTS

Through genome-wide analysis of next-generation miRNA sequencing reads generated from switchgrass roots, we identified 122 mature miRNAs, including 28 novel miRNAs. By comparing miRNA expression profiles of AM-inoculated and control switchgrass roots, we identified 15 AM-responsive miRNAs across lowland accession "Alamo", upland accession "Dacotah", and two upland/lowland F₁ hybrids. We used degradome sequencing to identify target genes of the AM-responsive miRNAs revealing targets of miRNAs residing on both K and N subgenomes. Notably, genes involved in copper ion binding were targeted by downregulated miRNAs, while upregulated miRNAs mainly targeted GRAS family transcription factors.

CONCLUSION

Through miRNA analysis and degradome sequencing, we revealed that both upland and lowland switchgrass genotypes as well as upland-lowland hybrids respond to AM by altering miRNA expression. We demonstrated complex GRAS transcription factor regulation by the miR171 family, with some miR171 family members being AM responsive while others remained static. Copper miRNA downregulation was common amongst the genotypes tested and we identified superoxide dismutases and laccases as targets, suggesting that these Cu-miRNAs are likely involved in ROS detoxification and lignin deposition, respectively. Other prominent targets of the Cu miRNAs were blue copper proteins. Overall, the potential effect of AM colonization on lignin deposition pathways in this biofuel crop highlights the importance of considering AM and miRNA in future biofuel crop development strategies.

A Comprehensive Evaluation of Effects on Water-Level Deficits on Tomato Polyphenol Composition, Nutritional Quality and Antioxidant Capacity

Tomatoes have high nutritional value and abundant bioactive compounds. Moderate water deficit irrigation alters metabolic levels of fruits, improving composition and quality. We investigated the effects of water deficit (T1, T2, T3, and T4) treatments and adequate irrigation (CK) on tomato polyphenol composition, antioxidant capacity, and nutritional quality. Compared with CK, the total flavonoid content increased by 33.66% and 44.73% in T1 and T2, and total phenols increased by 57.64%, 72.22%, and 55.78% in T1, T2, and T3, respectively. The T2 treatment significantly enhanced antioxidant’ capacities (ABTS, HSRA, FRAP, and DPPH). There were multiple groups of significant or extremely significant positive correlations between polyphenol components and antioxidant activity. For polyphenols and antioxidant capacity, the classification models divided the treatments: CK and T4 and T1−T3. The contents of soluble solids, soluble protein, vitamin C, and soluble sugar of the treatment groups were higher than those of CK. The soluble sugar positively correlated with sugar−acid ratios. In the PCA-based model, T3 in the first quadrant indicated the best treatment in terms of nutritional quality. Overall, comprehensive rankings using principal component analysis (PCA) revealed T2 > T1 > T3 > T4 > CK. Therefore, the T2 treatment is a suitable for improving quality and antioxidant capacity. This study provides novel insights into improving water-use efficiency and quality in the context of water scarcity worldwide.

Drought tolerance improvement in Solanum lycopersicum: an insight into "OMICS" approaches and genome editing

Solanum lycopersicum (tomato) is an internationally acclaimed vegetable crop that is grown worldwide. However, drought stress is one of the most critical challenges for tomato production, and it is a crucial task for agricultural biotechnology to produce drought-resistant cultivars. Although breeders have done a lot of work on the tomato to boost quality and quantity of production and enhance resistance to biotic and abiotic stresses, conventional tomato breeding approaches have been limited to improving drought tolerance because of the intricacy of drought traits. Many efforts have been made to better understand the mechanisms involved in adaptation and tolerance to drought stress in tomatoes throughout the years. "Omics" techniques, such as genomics, transcriptomics, proteomics, and metabolomics in combination with modern sequencing technologies, have tremendously aided the discovery of drought-responsive genes. In addition, the availability of biotechnological tools, such as plant transformation and the recently developed genome editing system for tomatoes, has opened up wider opportunities for validating the function of drought-responsive genes and the generation of drought-tolerant varieties. This review highlighted the recent progresses for tomatoes improvement against drought stress through "omics" and "multi-omics" technologies including genetic engineering. We have also discussed the roles of non-coding RNAs and genome editing techniques for drought stress tolerance improvement in tomatoes.

An Epigenetic Alphabet of Crop Adaptation to Climate Change

Crop adaptation to climate change is in a part attributed to epigenetic mechanisms which are related to response to abiotic and biotic stresses. Although recent studies increased our knowledge on the nature of these mechanisms, epigenetics remains under-investigated and still poorly understood in many, especially non-model, plants, Epigenetic modifications are traditionally divided into two main groups, DNA methylation and histone modifications that lead to chromatin remodeling and the regulation of genome functioning. In this review, we outline the most recent and interesting findings on crop epigenetic responses to the environmental cues that are most relevant to climate change. In addition, we discuss a speculative point of view, in which we try to decipher the “epigenetic alphabet” that underlies crop adaptation mechanisms to climate change. The understanding of these mechanisms will pave the way to new strategies to design and implement the next generation of cultivars with a broad range of tolerance/resistance to stresses as well as balanced agronomic traits, with a limited loss of (epi)genetic variability.

Metabolome and Transcriptome Analyses of Cucurbitacin Biosynthesis in Luffa (Luffa acutangula)

Cucurbitacins are extremely bitter compounds mainly present in Cucurbitaceae, where Luffa belongs. However, there is no comprehensive analysis of cucurbitacin biosynthesis in Luffa fruit. Therefore, this study analyzed bitter (WM709) and non-bitter (S1174) genotypes of Luffa to reveal the underlying mechanism of cucurbitacin biosynthesis by integrating metabolome and transcriptome analyses. A total of 422 metabolites were detected, including vitamins, essential amino acids, antioxidants, and antitumor substances. Of these, 131 metabolites showed significant differences between bitter (WM709) and non-bitter (S1174) Luffa fruits. The levels of isocucurbitacin B, cucurbitacin D, 23,24-dihydro cucurbitacin E, cucurbitacin F were significantly higher in bitter than in non-bitter Luffa. Transcriptome analysis showed that Bi, cytochromes P450s (CYP450s), and acyltransferase (ACT) of the cucurbitacin biosynthesis pathway, were significantly up-regulated. Moreover, drought stress and abscisic acid (ABA) activated genes of the cucurbitacin biosynthesis pathway. Furthermore, dual-luciferase reporter and yeast one-hybrid assays demonstrated that ABA-response element binding factor 1 (AREB1) binds to the Bi promoter to activate Bi expression. Comparative analysis of the Luffa and cucumber genomes showed that Bi, CYP450s, and ACT are located in the conserved syntenic loci, and formed a cucurbitacin biosynthesis cluster. This study provides important insights into major genes and metabolites of the cucurbitacin biosynthetic pathway, deepening the understanding of regulatory mechanisms of cucurbitacin biosynthesis in Luffa.

Research on Cancer Molecular Typing Based on High-Throughput Sequencing Technology

This paper studies the role of high-throughput measurement technology in cancer molecular typing. Based on the Dendrix algorithm, the model proposed in this paper selects the gene replication time as an inherent attribute that affects the frequency of gene mutations and adds it to the model. After setting the size of the gene set, compared with the Dendrix algorithm, the model does not need to delete the gene set that has been found in the process of searching the pathway, and it can find more driving pathway gene sets. Based on the high coverage and high exclusivity of the driving gene set in the pathway and the influence of gene covariates, this paper constructs an adaptive multiobjective optimization model. In order to overcome the problem of gene mutation heterogeneity, this model introduces gene covariates as the weight of gene mutation frequency so that the model is adaptive to each gene. The analysis of the research results shows the reliability of high-throughput sequencing technology.

Integrated transcriptome and small RNA sequencing analyses reveal a drought stress response network in Sophora tonkinensis

BACKGROUND

Sophora tonkinensis Gagnep is a traditional Chinese medical plant that is mainly cultivated in southern China. Drought stress is one of the major abiotic stresses that negatively impacts S. tonkinensis growth. However, the molecular mechanisms governing the responses to drought stress in S. tonkinensis at the transcriptional and posttranscriptional levels are not well understood.

RESULTS

To identify genes and miRNAs involved in drought stress responses in S. tonkinensis, both mRNA and small RNA sequencing was performed in root samples under control, mild drought, and severe drought conditions. mRNA sequencing revealed 66,476 unigenes, and the differentially expressed unigenes (DEGs) were associated with several key pathways, including phenylpropanoid biosynthesis, sugar metabolism, and quinolizidine alkaloid biosynthesis pathways. A total of 10 and 30 transcription factors (TFs) were identified among the DEGs under mild and severe drought stress, respectively. Moreover, small RNA sequencing revealed a total of 368 miRNAs, including 255 known miRNAs and 113 novel miRNAs. The differentially expressed miRNAs and their target genes were involved in the regulation of plant hormone signal transduction, the spliceosome, and ribosomes. Analysis of the regulatory network involved in the response to drought stress revealed 37 differentially expressed miRNA-mRNA pairs.

CONCLUSION

This is the first study to simultaneously profile the expression patterns of mRNAs and miRNAs on a genome-wide scale to elucidate the molecular mechanisms of the drought stress responses of S. tonkinensis. Our results suggest that S. tonkinensis implements diverse mechanisms to modulate its responses to drought stress.

Drought and crop yield

Episodes of water shortage occur in most agricultural regions of the world. Their durations and intensities increase, and their seasonal timing alters with changing climate. During the ontogenic cycle of crop plants, each development stage, such as seed germination, seedling establishment, vegetative root and shoot growth, flowering, pollination and seed and fruit development, is specifically sensitive to dehydration. Desiccation threatens yield and leads to specific patterns, depending on the type of crop plant and the harvested plant parts, e.g. leafy vegetables, tubers, tap roots or fruits. This review summarizes the effects of drought stress on crop plants and relates the dehydration-dependent yield penalty to the harvested organ and tissue. The control of shoot transpiration and the reorganization of root architecture are of core importance for maintaining proper plant water relationships. Upon dehydration, the provision and partitioning of assimilates and the uptake and distribution of nutrients define remaining growth activity. Domestication of crops by selection for high yield under high input has restricted the genetic repertoire for achieving drought stress tolerance. Introgression of suitable alleles from wild relatives into commercial cultivars might improve the ability to grow with less water. Future research activities should focus more on field studies in order to generate more realistic improvements to crops. Robotic field phenotyping should be integrated into genetic mapping for the identification of relevant traits.

The Effect of Low Temperature Stress on the Leaves and MicroRNA Expression of Potato Seedlings

In recent years, with the wanton destruction of the ecological environment by humans and the frequent occurrence of extreme bad weather, many places that should have been warm and blooming in spring have instead experienced the phenomenon of the “April blizzard,” which has seriously affected China's crops, especially spring potato production in most areas. Potato cultivars, especially potato seedlings, are sensitive to frost, and low temperature frost has become one of the most important abiotic stresses affecting potato production. Potato cold tolerance is regulated by a complex gene network. Although some low temperature resistant microRNAs have been identified, little is known about the role of miRNAs in response to low temperature stress in potato. Therefore, the objective of this study is to clarify the influence of low temperature stress on the miRNA expression of potato by comparing the expression differences of miRNA in potato which was treated with different low temperatures. For the study, 307 known miRNAs belonging to 73 small RNA families and 211 novel miRNAs were obtained. When the temperature decreased, the number of both known and novel miRNA decreased, and the minimum temperature was −2°C. Most of the miRNAs respond to low temperature, drought, and disease stress; some conserved miRNAs were first found to respond to low temperature stress in potato, such as stu-miR530, stu-miR156d, and stu-miR167b. The Gene Ontology, Kyoto Encyclopedia of Genes, and Genomes pathway enrichment analysis of 442 different expression miRNAs target genes indicated that there existed diversified low temperature responsive pathways, but Abscisic Acid was found likely to play a central coordinating role in response to low temperature stress in many metabolism pathways. Quantitative real-time PCR assays indicated that the related targets were negatively regulated by the tested different expression miRNAs during low temperature stress. The results indicated that miRNAs may play an important coordination role in response to low temperature stress in many metabolic pathways by regulating abscisic acid and gibberellin, which provided insight into the roles of miRNAs during low temperature stress and would be helpful for alleviating low temperature stress and promoting low temperature resistant breeding in potatoes.

miRNAs play critical roles in response to abiotic stress by modulating cross-talk of phytohormone signaling

One of the most interesting signaling molecules that regulates a wide array of adaptive stress responses in plants are the micro RNAs (miRNAs) that are a unique class of non-coding RNAs constituting novel mechanisms of post-transcriptional gene regulation. Recent studies revealed the role of miRNAs in several biotic and abiotic stresses by regulating various phytohormone signaling pathways as well as by targeting a number of transcription factors (TFs) and defense related genes. Phytohormones are signal molecules modulating the plant growth and developmental processes by regulating gene expression. Studies concerning miRNAs in abiotic stress response also show their vital roles in abiotic stress signaling. Current research indicates that miRNAs may act as possible candidates to create abiotic stress tolerant crop plants by genetic engineering. Yet, the detailed mechanism governing the dynamic expression networks of miRNAs in response to stress tolerance remains unclear. In this review, we provide recent updates on miRNA-mediated regulation of phytohormones combating various stress and its role in adaptive stress response in crop plants.

Identification and analysis of novel microRNAs provide insights to reproductive capacity of the cultured Asian yellow pond turtle Mauremys mutica

The crucial roles of miRNAs in regulating animal growth, development, and disease resistance have been extensively reported, but their roles in relation to the reproductive capacity of aquatic animals (numbers of eggs laid and hatchlings), especially reptiles, remain unclear. In this study, high-throughput sequencing technology was used to screen miRNAs related to reproductive capacity based on the construction of a cDNA library of ovaries from higher-fecundity (HF) and lower-fecundity (LF) M. mutica. The results showed that 15,767,494 (93.98%) and 14,137,621 (94.17%) high-quality reads were obtained from the HF and LF groups, respectively. We screened 131 miRNAs that were differentially expressed between the HF and LF groups, of which 78 were upregulated and 53 were downregulated compared with the M. mutica reference genome. GO and KEGG pathway enrichment analyses of the target genes of differentially expressed miRNAs revealed significant differences in the enrichment frequencies of genes associated with ATP binding and proteolysis between the HF and LF groups, while the tricarboxylic acid cycle, glucagon signaling pathway and vitamin B6 metabolic pathway were shown to potentially help determine reproductive capacity. Ten miRNAs were verified by qRT-PCR to confirm the reliability and accuracy of the sequencing results, and a miRNA-mRNA target gene interaction network was constructed. These results will further our understanding of the regulatory mechanism of miRNAs in regards to turtle reproductive capacity.

Chemical and structural quality traits during postharvest ripening regulated by chromosome segments from a wild relative of tomato Solanum pennellii IL4-2 and IL5-1

Tomato is usually harvested at an early ripening stage with high firmness suitable for storage and transportation but lacks many quality parameters such as sugars, organic acids, and phenolics. In a recent study, we have selected introgression lines (ILs) IL4-2 and IL5-1, developed from a cross between the Solanum pennellii and the Solanum lycopersicum M82, that exhibit differentiated postharvest shelf-life characteristics in the fruit compared to M82 and the rest of the ILs. Here, we first structurally and biochemically characterized IL4-2, IL5-1, and their parent M82 to decipher the cell wall mechanistic difference between soft (IL4-2) and firm (IL5-1) lines at two postharvest ripening periods. Generally, IL4-2 had more active cell wall modifications in terms of ripening-related gene expression, water-soluble pectin, and cell wall structure under the microscope, which probably makes this line softer than IL5-1. We also evaluated these lines based on commercial quality parameters, sugars, phenolics, organic, and amino acids to gain insight into their commercial and functional quality and reveal noticeable differences. In summary, the contribution of the S. pennellii IL5-1 and IL4-2 to the shelf life of the tomato was structurally characterized, and the component differences meeting the quality criteria were revealed.

High-throughput miRNA deep sequencing in response to drought stress in sugarcane

Drought is a major factor which reduces cane growth and productivity. In the present study, we sequenced drought susceptible (V1) and drought tolerant (V2) sugarcane varieties using high-throughput miRNA deep sequencing method to study the regulation of gene expression by miRNAs during drought stress in sugarcane. A total of 1224 conserved miRNAs which belong to 89 miRNA families were identified and 38% of the differentially regulated miRNAs were common for both varieties. Additionally 435 novel miRNAs were also identified from four small RNA libraries. We identified 145 miRNAs that were differentially expressed in susceptible variety (V1-31) and 143 miRNAs differentially expressed in the tolerant variety (V2-31). Target prediction revealed that the genes mainly encoded transcription factors, proteins, phosphatase and kinases involved in signal transduction pathways, integral component of membrane and inorganic ion transport metabolism, enzymes involved in carbohydrate transport and metabolism and drought-stress-related proteins involved in defense mechanisms. Pathway analysis of targets revealed that "General function prediction only" was the most significant pathway observed in both tolerant and susceptible genotypes followed by "signal transduction mechanisms". Functional annotation of the transcripts revealed genes like calcium-dependent protein kinase, respiratory burst oxidase, caffeic acid 3-O-methyltransferase, peroxidase, calmodulin, glutathione S-transferase and transcription factors like MYB, WRKY that are involved in drought tolerant pathways. qRT-PCR was used to verify the expression levels of miRNAs and their potential targets obtained from RNA sequencing results.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-021-02857-x.

Evolution of plasticity in production and transgenerational inheritance of small RNAs under dynamic environmental conditions

In a changing environment, small RNAs (sRNAs) play an important role in the post-transcriptional regulation of gene expression and can vary in abundance depending on the conditions experienced by an individual (phenotypic plasticity) and its parents (non-genetic inheritance). Many sRNAs are unusual in that they can be produced in two ways, either using genomic DNA as the template (primary sRNAs) or existing sRNAs as the template (secondary sRNAs). Thus, organisms can evolve rapid plastic responses to their current environment by adjusting the amplification rate of sRNA templates. sRNA levels can also be transmitted transgenerationally by the direct transfer of either sRNAs or the proteins involved in amplification. Theory is needed to describe the selective forces acting on sRNA levels, accounting for the dual nature of sRNAs as regulatory elements and templates for amplification and for the potential to transmit sRNAs and their amplification agents to offspring. Here, we develop a model to study the dynamics of sRNA production and inheritance in a fluctuating environment. We tested the selective advantage of mutants capable of sRNA-mediated phenotypic plasticity within resident populations with fixed levels of sRNA transcription. Even when the resident was allowed to evolve an optimal constant rate of sRNA production, plastic amplification rates capable of responding to environmental conditions were favored. Mechanisms allowing sRNA transcripts or amplification agents to be inherited were favored primarily when parents and offspring face similar environments and when selection acts before the optimal level of sRNA can be reached within the organism. Our study provides a clear set of testable predictions for the evolution of sRNA-related mechanisms of phenotypic plasticity and transgenerational inheritance.

MicroRNAs: Potential Targets for Developing Stress-Tolerant Crops

Crop yield is challenged every year worldwide by changing climatic conditions. The forecasted climatic scenario urgently demands stress-tolerant crop varieties to feed the ever-increasing global population. Molecular breeding and genetic engineering approaches have been frequently exploited for developing crops with desired agronomic traits. Recently, microRNAs (miRNAs) have emerged as powerful molecules, which potentially serve as expression markers during stress conditions. The miRNAs are small non-coding endogenous RNAs, usually 20-24 nucleotides long, which mediate post-transcriptional gene silencing and fine-tune the regulation of many abiotic- and biotic-stress responsive genes in plants. The miRNAs usually function by specifically pairing with the target mRNAs, inducing their cleavage or repressing their translation. This review focuses on the exploration of the functional role of miRNAs in regulating plant responses to abiotic and biotic stresses. Moreover, a methodology is also discussed to mine stress-responsive miRNAs from the enormous amount of transcriptome data available in the public domain generated using next-generation sequencing (NGS). Considering the functional role of miRNAs in mediating stress responses, these molecules may be explored as novel targets for engineering stress-tolerant crop varieties.

Increased water use efficiency in miR396-downregulated tomato plants

MicroRNAs regulate plant development and responses to biotic and abiotic stresses but their impact on water use efficiency (WUE) is poorly known. Increasing WUE is a major task in crop improvement programs aimed to meet the challenges posed by the reduction in water availability associated with the ongoing climatic change. We have examined the physiological and molecular response to water stress of tomato (Solanum lycopersicum L.) plants downregulated for miR396 by target mimicry. In water stress conditions, miR396-downregulated plants displayed reduced transpiration and a less then proportional decrease in the photosynthetic rate that resulted in higher WUE. The increase in WUE was associated with faster foliar accumulation of abscisic acid (ABA), with the induction of several drought-protective genes and with the activation of the jasmonic acid (JA) and γ-aminobutyric acid (GABA) pathways. We propose a model in which the downregulation of miR396 leads to the activation of a complex molecular response to water stress. This response acts synergistically with a set of leaf morphological modifications to increase stomatal closure and preserve the efficiency of the photosynthetic activity, ultimately resulting in higher WUE.

Genome-Wide Identification of Key Candidate microRNAs and Target Genes Associated with Peanut Drought Tolerance

Peanut is an important crash crop worldwide, and it is often threatened by drought stress due to unexpected extreme weather events. In this work, NH5 and FH18 were selected as drought-tolerant and drought-sensitive varieties, respectively. Comparison of their physiological responses revealed that NH5 showed less wilting, higher relative water content and lower water loss rate of detached leaves, lower electrolyte leakage, and stronger antioxidant ability under drought stress than did FH18. Based on comparative transcriptomic analysis, 5376 differentially expressed mRNAs were commonly identified in the two varieties, and 2993 genes specifically changed in the drought-tolerant variety and were mainly enriched in photosynthesis-antenna proteins and photosynthetic pathways. Furthermore, 73 microRNAs (miRNAs) were differentially expressed in the drought tolerance variety specifically under drought stress; of these, two key candidate miRNAs, novel miR_416 and novel miR_73, were identified, and the majority of their target genes were enriched in phenylpropanoid biosynthesis, linoleic acid metabolism, and cutin, suberine, and wax biosynthesis. This study lays the foundation for the analysis of the molecular mechanism of drought tolerance and promotes the genetic improvement of peanut drought tolerance.

Genome-wide investigation of proline transporter (ProT) gene family in tomato: Bioinformatics and expression analyses in response to drought stress

Proline has various functions in plants, such as growth, development and stress response to biotic and abiotic factors. Therefore, proline accumulation and transport are vital for crop production in higher quality and quantity. The present study addresses genome-wide identification and bioinformatics analyses of tomato (Solanum lycopersicum) proline transporter (ProT) genes and their expression profiles under drought stress. The analyses indicated four novel ProT genes (SlProTs) in the tomato genome and their protein lengths ranged from 439 to 452 amino acid residues. All SlProTs contained a PF01490 (transmembrane amino acid transporter protein) domain and seven exons, and they had a basic pI. The phylogeny analysis proved that monocot-dicot divergence was not present and the SlProT proteins were distinct from the ProT proteins in monocots and Arabidopsis. Based on the digital expression analysis, SlProT1 and SlProT2 genes seemed to be more active than the others in response to abiotic stress conditions. However, detected by RT-qPCR, the expression levels of all SlProT genes under drought stress were similar. The promotor analyses of SlProT genes revealed that they contained many transcription factors binding sites in cis-elements, such as MYB, Dof, Hox, bZIP, bHLH, AP2/ERF and WRKY. Finally, our findings could contribute to the understanding of SlProT genes and proline metabolism in plants.

Genome-wide identification and characterization of long non-coding RNAs related to grain yield in foxtail millet [Setaria italica (L.) P. Beauv.]

Background

Long noncoding RNAs (lncRNAs) have been reported to play critical roles in diverse growth and development processes in plants. However, the systematic identification and characterization of lncRNAs in foxtail millet is nearly blank.

Results

In this study, we performed high-throughput sequencing of young spikelets from four foxtail millet varieties in different yield levels at booting stage. As a result, a total of 12,378 novel lncRNAs were identified, and 70 were commonly significantly differentially expressed in comparisons between high-yield varieties and conventional varieties, suggesting that they involved in yield formation and regulation in foxtail millet. Functional analysis revealed that among the 70 significantly differentially expressed lncRNAs, 67 could transcriptionally modulate target genes in cis and in trans. Moreover, 18 lncRNAs related to grain yield in foxtail millet were predicted to function as miRNA target mimics and regulate gene expression by competing for the interaction between miRNAs and their target mRNAs.

Conclusion

Our results will provide materials for elucidation of the molecular mechanisms of lncRNAs participate in yield regulation, and will contribute to high yield foxtail millet breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07272-9.

Identification of conserved genes linked to responses to abiotic stresses in leaves among different plant species

As a consequence of global climate change, certain stress factors that have a negative impact on crop productivity such as heat, cold, drought and salinity are becoming increasingly prevalent. We conducted a meta-analysis to identify genes conserved across plant species involved in (1) general abiotic stress conditions, and (2) specific and unique abiotic stress factors (drought, salinity, extreme temperature) in leaf tissues. We collected raw data and re-analysed eight RNA-Seq studies using our previously published bioinformatic pipeline. A total of 68 samples were analysed. Gene set enrichment analysis was performed using MapMan and PageMan whereas DAVID (Database for Annotation, Visualisation and Integrated Discovery) was used for metabolic process enrichment analysis. We identified of a total of 5122 differentially expressed genes when considering all abiotic stresses (3895 were upregulated and 1227 were downregulated). Jasmonate-related genes were more commonly upregulated by drought, whereas gibberellin downregulation was a key signal for drought and heat. In contrast, cold stress clearly upregulated genes involved in ABA (abscisic acid), cytokinin and gibberellins. A gene (non-phototrophic hypocotyl) involved in IAA (indoleacetic acid) response was induced by heat. Regarding secondary metabolism, as expected, MVA pathway (mevalonate pathway), terpenoids and alkaloids were generally upregulated by all different stresses. However, flavonoids, lignin and lignans were more repressed by heat (cinnamoyl coA reductase 1 and isopentenyl pyrophosphatase). Cold stress drastically modulated genes involved in terpenoid and alkaloids. Relating to transcription factors, AP2-EREBP, MADS-box, WRKY22, MYB, homoebox genes members were significantly modulated by drought stress whereas cold stress enhanced AP2-EREBPs, bZIP members, MYB7, BELL 1 and one bHLH member. C2C2-CO-LIKE, MADS-box and a homeobox (HOMEOBOX3) were mostly repressed in response to heat. Gene set enrichment analysis showed that ubiquitin-mediated protein degradation was enhanced by heat, which unexpectedly repressed glutaredoxin genes. Cold stress mostly upregulated MAP kinases (mitogen-activated protein kinase). Findings of this work will allow the identification of new molecular markers conserved across crops linked to major genes involved in quantitative agronomic traits affected by different abiotic stress.

Novel targets for engineering Physostegia chlorotic mottle and tomato brown rugose fruit virus-resistant tomatoes: in silico prediction of tomato microRNA targets

Background

Physostegia chlorotic mottle virus (PhCMoV; genus: Alphanucleorhabdovirus, family: Rhabdoviridae) and tomato brown rugose fruit virus (ToBRFV; genus: Tobamovirus, family: Virgaviridae) are newly emerging plant viruses that have a dramatic effect on tomato production. Among various known virus-control strategies, RNAi-mediated defence has shown the potential to protect plants against various pathogens including viral infections. Micro(mi)RNAs play a major role in RNAi-mediated defence.

Methods

Using in silico analyses, we investigated the possibility of tomato-encoded miRNAs (TomiRNA) to target PhCMoV and ToBRFV genomes using five different algorithms, i.e., miRanda, RNAhybrid, RNA22, Tapirhybrid and psRNATarget.

Results

The results revealed that 14 loci on PhCMoV and 10 loci on ToBRFV can be targeted by the TomiRNAs based on the prediction of at least three algorithms. Interestingly, one TomiRNA, miR6026, can target open reading frames from both viruses, i.e., the phosphoprotein encoding gene of PhCMoV, and the two replicase components of ToBRFV. There are currently no commercially available PhCMoV- or ToBRFV-resistant tomato varieties, therefore the predicted data provide useful information for the development of PhCMoV- and ToBFRV-resistant tomato plants.

Global mRNA and microRNA expression dynamics in response to anthracnose infection in sorghum

BACKGROUND

Anthracnose is a damaging disease of sorghum caused by the fungal pathogen Colletotrichum sublineolum. Genome-wide mRNA and microRNA (miRNA) profiles of resistant and susceptible sorghum genotypes were studied to understand components of immune responses, and fungal induced miRNA and target gene networks.

RESULTS

A total of 18 mRNA and 12 miRNA libraries from resistant and susceptible sorghum lines were sequenced prior to and after inoculation with C. sublineolum. Significant differences in transcriptomes of the susceptible and resistant genotypes were observed with dispersion distance and hierarchical cluster tree analyses. Of the total 33,032 genes predicted in the sorghum genome, 19,593 were induced by C. sublineolum, and 15,512 were differentially expressed (DEGs) between the two genotypes. The resistant line was marked by significant reprogramming of the transcriptome at 24 h post inoculation (hpi), and a decrease at 48 hpi, whereas the susceptible line displayed continued changes in gene expression concordant with elevated fungal growth in the susceptible genotype. DEGs encode proteins implicated in diverse functions including photosynthesis, synthesis of tetrapyrrole, carbohydrate and secondary metabolism, immune signaling, and chitin binding. Genes encoding immune receptors, MAPKs, pentatricopeptide repeat proteins, and WRKY transcription factors were induced in the resistant genotype. In a parallel miRNA profiling, the susceptible line displayed greater number of differentially expressed miRNAs than the resistant line indicative of a widespread suppression of gene expression. Interestingly, we found 75 miRNAs, including 36 novel miRNAs, which were differentially expressed in response to fungal inoculation. The expression of 50 miRNAs was significantly different between resistant and susceptible lines. Subsequently, for 35 differentially expressed miRNAs, the corresponding 149 target genes were identified. Expression of 56 target genes were significantly altered after inoculation, showing inverse expression with the corresponding miRNAs.

CONCLUSIONS

We provide insights into genome wide dynamics of mRNA and miRNA profiles, biological and cellular processes underlying host responses to fungal infection in sorghum. Resistance is correlated with early transcriptional reprogramming of genes in various pathways. Fungal induced genes, miRNAs and their targets with a potential function in host responses to anthracnose were identified, opening avenues for genetic dissection of resistance mechanisms.

Comparative Study of Pine Reference Genomes Reveals Transposable Element Interconnected Gene Networks

Sequencing the giga-genomes of several pine species has enabled comparative genomic analyses of these outcrossing tree species. Previous studies have revealed the wide distribution and extraordinary diversity of transposable elements (TEs) that occupy the large intergenic spaces in conifer genomes. In this study, we analyzed the distribution of TEs in gene regions of the assembled genomes of Pinus taeda and Pinus lambertiana using high-performance computing resources. The quality of draft genomes and the genome annotation have significant consequences for the investigation of TEs and these aspects are discussed. Several TE families frequently inserted into genes or their flanks were identified in both species’ genomes. Potentially important sequence motifs were identified in TEs that could bind additional regulatory factors, promoting gene network formation with faster or enhanced transcription initiation. Node genes that contain many TEs were observed in multiple potential transposable element-associated networks. This study demonstrated the increased accumulation of TEs in the introns of stress-responsive genes of pines and suggests the possibility of rewiring them into responsive networks and sub-networks interconnected with node genes containing multiple TEs. Many such regulatory influences could lead to the adaptive environmental response clines that are characteristic of naturally spread pine populations.

Stress induced dynamic adjustment of conserved miR164:NAC module

Aims including the rationale

Salinity and drought are the two major stresses limiting the productivity of economically important crops such as Glycine max (soybean). The incidence of these stresses during the pod development stages affects the quality and quantity of seeds, which compromise the yield of soybean. The miR164:NAC module has been shown to play a critical role in regulating the response to salt and drought stress in several plant species. However, biological role of miR164:NAC module in salt stress in soybean is not fully understood.

Methods

In this study, we identified 215 salt responsive miRNAs, using miScript miRNA array with a sensitive and a tolerant soybean genotype, William82 and INCASoy36, respectively. The targets of these salt regulated miRNAs were searched in the degradome datasets.

Key results

It was found that four salt stress deregulated miRNAs targeted the NAC transcription factor and among these miR164k and miR408d showed antagonistic expression in the two soybean genotypes. The expression of miR164k was higher in salt tolerant INCASoy36 as compared to salt sensitive William82, under unstressed conditions. However under salt stress, miR164k was downregulated in INCASoy36 (-2.65 fold), whereas it was upregulated in William82 (4.68 fold). A transient co-expression assay validated that gma-miR164k directs the cleavage of GmNAC1 transcript. Bioinformatics analysis revealed that the regulation of NAC transcription factor family by members of miR164 family is conserved across many species. The dynamic expression profiles of miR164 and NAC-TFs were captured in different tissues of rice, tobacco, and two soybean genotypes under drought and salt stress conditions.

Main conclusion

Collectively, our results suggest that genetically determined dynamic modulation of the conserved miR164:NAC-TF module may play an important role in determining the adaptive response of plants to stress.

Integrated Analysis of Small RNA, Transcriptome, and Degradome Sequencing Reveals the Water-Deficit and Heat Stress Response Network in Durum Wheat

Water-deficit and heat stress negatively impact crop production. Mechanisms underlying the response of durum wheat to such stresses are not well understood. With the new durum wheat genome assembly, we conducted the first multi-omics analysis with next-generation sequencing, providing a comprehensive description of the durum wheat small RNAome (sRNAome), mRNA transcriptome, and degradome. Single and combined water-deficit and heat stress were applied to stress-tolerant and -sensitive Australian genotypes to study their response at multiple time-points during reproduction. Analysis of 120 sRNA libraries identified 523 microRNAs (miRNAs), of which 55 were novel. Differentially expressed miRNAs (DEMs) were identified that had significantly altered expression subject to stress type, genotype, and time-point. Transcriptome sequencing identified 49,436 genes, with differentially expressed genes (DEGs) linked to processes associated with hormone homeostasis, photosynthesis, and signaling. With the first durum wheat degradome report, over 100,000 transcript target sites were characterized, and new miRNA-mRNA regulatory pairs were discovered. Integrated omics analysis identified key miRNA-mRNA modules (particularly, novel pairs of miRNAs and transcription factors) with antagonistic regulatory patterns subject to different stresses. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis revealed significant roles in plant growth and stress adaptation. Our research provides novel and fundamental knowledge, at the whole-genome level, for transcriptional and post-transcriptional stress regulation in durum wheat.

High-Throughput Sequencing and Expression Analysis Suggest the Involvement of Pseudomonas putida RA-Responsive microRNAs in Growth and Development of Arabidopsis

Beneficial soil microorganisms largely comprise of plant growth-promoting rhizobacteria (PGPR), which adhere to plant roots and facilitate their growth and development. Pseudomonas putida (RA) strain MTCC5279 is one such PGPR that exhibits several characteristics of plant growth promotion, such as P-solubilization, and siderophores and IAA production. Plant–PGPR interactions are very complex phenomena, and essentially modulate the expression of numerous genes, consequently leading to changes in the physiological, biochemical, cellular and molecular responses of plants. Therefore, in order to understand the molecular bases of plant–PGPR interactions, we carried out the identification of microRNAs from the roots of Arabidopsis upon P. putida RA-inoculation, and analyses of their expression. MicroRNAs (miRNAs) are 20- to 24-nt non-coding small RNAs known to regulate the expression of their target genes. Small RNA sequencing led to the identification of 293 known and 67 putative novel miRNAs, from the control and RA-inoculated libraries. Among these, 15 known miRNAs showed differential expression upon RA-inoculation in comparison to the control, and their expressions were corroborated by stem-loop quantitative real-time PCR. Overall, 28,746 and 6931 mRNAs were expected to be the targets of the known and putative novel miRNAs, respectively, which take part in numerous biological, cellular and molecular processes. An inverse correlation between the expression of RA-responsive miRNAs and their target genes also strengthened the crucial role of RA in developmental regulation. Our results offer insights into the understanding of the RA-mediated modulation of miRNAs and their targets in Arabidopsis, and pave the way for the further exploitation and characterization of candidate RA-responsive miRNA(s) for various crop improvement strategies directed towards plant sustainable growth and development.

Pathogenesis related protein-1 (PR-1) genes in tomato (Solanum lycopersicum L.): Bioinformatics analyses and expression profiles in response to drought stress

The pathogenesis-related protein 1 (PR-1) gene family play important roles in the plant metabolism in response to biotic and abiotic stresses. The present study aimed genome-wide identification and bioinformatics analyses of PR-1 genes in tomato (Solanum lycopersicum L.). The analyses resulted in the identification of 13 novel SlPR-1 genes, each of which produce a protein belonging to the CAP superfamily (PF00188). The KEGG annotation analyses revealed that the SlPR-1 proteins functioned in the environmental information processing (09130). The expression patterns of the PR-1 genes and some stress-related physiological parameters were investigated in Fusarium oxysporum sensitive and tolerant tomato varieties under drought stress. The drought stress leaded upregulation of all SlPR-1 genes, reaching up to 50 folds. The results indicate that the SlPR-1 genes play active roles in response to drought. This is the first study exhibiting the expression profiles of SlPR-1 genes under an abiotic stress, drought, in tomato.

A novel strigolactone-miR156 module controls stomatal behaviour during drought recovery

miR156 is a conserved microRNA whose role and induction mechanisms under stress are poorly known. Strigolactones are phytohormones needed in shoots for drought acclimation. They promote stomatal closure ABA-dependently and independently; however, downstream effectors for the former have not been identified. Linkage between miR156 and strigolactones under stress has not been reported. We compared ABA accumulation and sensitivity as well as performances of wt and miR156-overexpressing (miR156-oe) tomato plants during drought. We also quantified miR156 levels in wt, strigolactone-depleted and strigolactone-treated plants, exposed to drought stress. Under irrigated conditions, miR156 overexpression and strigolactone treatment led to lower stomatal conductance and higher ABA sensitivity. Exogenous strigolactones were sufficient for miR156 accumulation in leaves, while endogenous strigolactones were required for miR156 induction by drought. The "after-effect" of drought, by which stomata do not completely re-open after rewatering, was enhanced by both strigolactones and miR156. The transcript profiles of several miR156 targets were altered in strigolactone-depleted plants. Our results show that strigolactones act as a molecular link between drought and miR156 in tomato, and identify miR156 as a mediator of ABA-dependent effect of strigolactones on the after-effect of drought on stomata. Thus, we provide insights into both strigolactone and miR156 action on stomata.

Post-transcriptional regulation of several biological processes involved in latex production in Hevea brasiliensis

Background

Small RNAs modulate plant gene expression at both the transcriptional and post-transcriptional level, mostly through the induction of either targeted DNA methylation or transcript cleavage, respectively. Small RNA networks are involved in specific plant developmental processes, in signaling pathways triggered by various abiotic stresses and in interactions between the plant and viral and non-viral pathogens. They are also involved in silencing maintenance of transposable elements and endogenous viral elements. Alteration in small RNA production in response to various environmental stresses can affect all the above-mentioned processes. In rubber trees, changes observed in small RNA populations in response to trees affected by tapping panel dryness, in comparison to healthy ones, suggest a shift from a transcriptional to a post-transcriptional regulatory pathway. This is the first attempt to characterise small RNAs involved in post-transcriptional silencing and their target transcripts in Hevea.

Methods

Genes producing microRNAs (MIR genes) and loci producing trans-activated small interfering RNA (ta-siRNA) were identified in the clone PB 260 re-sequenced genome. Degradome libraries were constructed with a pool of total RNA from six different Hevea tissues in stressed and non-stressed plants. The analysis of cleaved RNA data, associated with genomics and transcriptomics data, led to the identification of transcripts that are affected by 20–22 nt small RNA-mediated post-transcriptional regulation. A detailed analysis was carried out on gene families related to latex production and in response to growth regulators.

Results

Compared to other tissues, latex cells had a higher proportion of transcript cleavage activity mediated by miRNAs and ta-siRNAs. Post-transcriptional regulation was also observed at each step of the natural rubber biosynthesis pathway. Among the genes involved in the miRNA biogenesis pathway, our analyses showed that all of them are expressed in latex. Using phylogenetic analyses, we show that both the Argonaute and Dicer-like gene families recently underwent expansion. Overall, our study underlines the fact that important biological pathways, including hormonal signalling and rubber biosynthesis, are subject to post-transcriptional silencing in laticifers.

Bioinformatics Resources for Plant Abiotic Stress Responses: State of the Art and Opportunities in the Fast Evolving -Omics Era

Abiotic stresses are among the principal limiting factors for productivity in agriculture. In the current era of continuous climate changes, the understanding of the molecular aspects involved in abiotic stress response in plants is a priority. The rise of -omics approaches provides key strategies to promote effective research in the field, facilitating the investigations from reference models to an increasing number of species, tolerant and sensitive genotypes. Integrated multilevel approaches, based on molecular investigations at genomics, transcriptomics, proteomics and metabolomics levels, are now feasible, expanding the opportunities to clarify key molecular aspects involved in responses to abiotic stresses. To this aim, bioinformatics has become fundamental for data production, mining and integration, and necessary for extracting valuable information and for comparative efforts, paving the way to the modeling of the involved processes. We provide here an overview of bioinformatics resources for research on plant abiotic stresses, describing collections from -omics efforts in the field, ranging from raw data to complete databases or platforms, highlighting opportunities and still open challenges in abiotic stress research based on -omics technologies.

The Tomato Translational Landscape Revealed by Transcriptome Assembly and Ribosome Profiling

Recent applications of translational control in Arabidopsis (Arabidopsis thaliana) highlight the potential power of manipulating mRNA translation for crop improvement. However, to what extent translational regulation is conserved between Arabidopsis and other species is largely unknown, and the translatome of most crops remains poorly studied. Here, we combined de novo transcriptome assembly and ribosome profiling to study global mRNA translation in tomato (Solanum lycopersicum) roots. Exploiting features corresponding to active translation, we discovered widespread unannotated translation events, including 1,329 upstream open reading frames (uORFs) within the 5' untranslated regions of annotated coding genes and 354 small ORFs (sORFs) among unannotated transcripts. uORFs may repress translation of their downstream main ORFs, whereas sORFs may encode signaling peptides. Besides evolutionarily conserved sORFs, we uncovered 96 Solanaceae-specific sORFs, revealing the importance of studying translatomes directly in crops. Proteomic analysis confirmed that some of the unannotated ORFs generate stable proteins in planta. In addition to defining the translatome, our results reveal the global regulation by uORFs and microRNAs. Despite diverging over 100 million years ago, many translational features are well conserved between Arabidopsis and tomato. Thus, our approach provides a high-throughput method to discover unannotated ORFs, elucidates evolutionarily conserved and unique translational features, and identifies regulatory mechanisms hidden in a crop genome.

Expression of miR159 Is Altered in Tomato Plants Undergoing Drought Stress

In a scenario of global climate change, water scarcity is a major threat for agriculture, severely limiting crop yields. Therefore, alternatives are urgently needed for improving plant adaptation to drought stress. Among them, gene expression reprogramming by microRNAs (miRNAs) might offer a biotechnologically sound strategy. Drought-responsive miRNAs have been reported in many plant species, and some of them are known to participate in complex regulatory networks via their regulation of transcription factors involved in water stress signaling. We explored the role of miR159 in the response of Solanum lycopersicum Mill. plants to drought stress by analyzing the expression of sly-miR159 and its target SlMYB transcription factor genes in tomato plants of cv. Ailsa Craig grown in deprived water conditions or in response to mechanical damage caused by the Colorado potato beetle, a devastating insect pest of Solanaceae plants. Results showed that sly-miR159 regulatory function in the tomato plants response to distinct stresses might be mediated by differential stress-specific MYB transcription factor targeting. sly-miR159 targeting of SlMYB33 transcription factor transcript correlated with accumulation of the osmoprotective compounds proline and putrescine, which promote drought tolerance. This highlights the potential role of sly-miR159 in tomato plants’ adaptation to water deficit conditions.

Identification of Stress Associated microRNAs in Solanum lycopersicum by High-Throughput Sequencing

Tomato (Solanum lycopersicum) is one of the most important crops around the world and also a model plant to study response to stress. High-throughput sequencing was used to analyse the microRNA (miRNA) profile of tomato plants undergoing five biotic and abiotic stress conditions (drought, heat, P. syringae infection, B. cinerea infection, and herbivore insect attack with Leptinotarsa decemlineata larvae) and one chemical treatment with a plant defence inducer, hexanoic acid. We identified 104 conserved miRNAs belonging to 37 families and we predicted 61 novel tomato miRNAs. Among those 165 miRNAs, 41 were stress-responsive. Reverse transcription quantitative PCR (RT-qPCR) was used to validate high-throughput expression analysis data, confirming the expression profiles of 10 out of 11 randomly selected miRNAs. Most of the differentially expressed miRNAs were stress-specific, except for sly-miR167c-3p upregulated in B. cinerea and P. syringae infection, sly-newmiR26-3p upregulated in drought and Hx treatment samples, and sly-newmiR33-3p, sly-newmiR6-3p and sly-newmiR8-3p differentially expressed both in biotic and abiotic stresses. From mature miRNAs sequences of the 41 stress-responsive miRNAs 279 targets were predicted. An inverse correlation between the expression profiles of 4 selected miRNAs (sly-miR171a, sly-miR172c, sly-newmiR22-3p and sly-miR167c-3p) and their target genes (Kinesin, PPR, GRAS40, ABC transporter, GDP and RLP1) was confirmed by RT-qPCR. Altogether, our analysis of miRNAs in different biotic and abiotic stress conditions highlight the interest to understand the functional role of miRNAs in tomato stress response as well as their putative targets which could help to elucidate plants molecular and physiological adaptation to stress.