Tamarix microRNA Profiling Reveals New Insight into Salt Tolerance

Molecular mechanism of flower colour formation in Rhododendron simsii Planchon revealed by integration of microRNAome and RNAomics

Systems-wide understanding of gene expression profile regulating flower colour formation in Rhododendron simsii Planchon is insufficient. In this research, integration analysis of ribonucleic acid (RNA)omics and microRNAome were performed to reveal the molecular mechanism of flower colour formation in three R. simsii varieties with red, pink and crimson flowers, respectively. Totally, 3129, 5755 and 5295 differentially expressed gene (DEG)s were identified through comparative transcriptome analysis between 'Red variety' and 'Pink variety' (1507 up-regulated and 1622 down-regulated), 'Red variety' and 'Crimson variety' (2148 up-regulated 3607 down-regulated), as well as 'Pink variety' and 'Crimson variety' (2089 up-regulated and 3206 down-regulated), which were involved in processes of 'catalytic activity', 'binding', 'metabolic process' and 'cellular process', as well as pathways of 'metabolic pathways', 'biosynthesis of secondary metabolites', 'plant-pathogen interaction' and 'phenylpropanoid biosynthesis'. A total of 215 miRNAs, containing 153 known miRNAs belonging to 57 families and 62 novel miRNA, were involved in flower colour formation. In particular, 55 miRNAs were significantly differently expressed. Based on miRNA-mRNA regulatory network, ath-miR5658 could affect the synthesis of pelargonidin, cyanidin and delphinidin through downregulating accumulation of anthocyanidin 3-O-glucosyltransferase; ath-miR868-3p could regulate isoflavonoid biosynthesis through downregulating expression of CYP81E1/E7; ath-miR156g regulated the expression of flavonoid 3',5'-hydroxylase; and ath-miR829-5p regulated flavonol synthasein flavonoid biosynthesis process. This research will provide important roles in breeding new varieties with rich flower colour.

Small RNA and Degradome Sequencing Reveal Roles of miRNAs in the Petal Color Fading of Malus Crabapple

The Malus crabapple is an important woody ornamental plant. The fading of petals during its development significantly affects their ornamental value. Petal color is related to anthocyanin content and miRNAs play an important role in the post-transcriptional regulation of anthocyanin synthesis. However, the mechanisms underlying miRNA regulation of petal fading have rarely been studied. Transcriptome and small RNA sequencing of petals from the blooming phases of Malus. 'Indian Summer' varieties S1 (small bud), S2 (initial-flowering), and S3 (late-flowering) allowed us to identify 230 known miRNAs and 17 novel miRNAs, including 52 differentially expressed miRNAs which targeted 494 genes and formed 823 miRNA-target pairs. Based on the target gene annotation results, miRNA-target pairs were screened that may be involved in the fading process of Malus crabapple petals through three different pathways: anthocyanin synthesis, transport, and degradation, involving mcr-miR858-MYB1\MYB5 and mcr-miR396-McCHI inhibiting anthocyanin synthesis; mcr-miR167, mcr-miR390, mcr-miR535, and mcr-miR858 inhibiting anthocyanin transport from the cytoplasm to the vacuole by targeting ABC transporter genes (ABCB, ABCC, ABCD, and ABCG); and mcr-miR398 targeting the superoxide dismutase genes (CZSOD2 and CCS) to accelerate anthocyanin degradation. These findings offer a novel approach to understanding the mechanism of petal fading and serve as a reference for other plants with floral fading.

A comprehensive analysis of transcriptomic data for comparison of plants with different photosynthetic pathways in response to drought stress

The main factor leading to a decrease in crop productivity is abiotic stresses, particularly drought. Plants with C4 and CAM photosynthesis are better adapted to drought-prone areas than C3 plants. Therefore, it is beneficial to compare the stress response of plants with different photosynthetic pathways. Since most crops are C3 and C4 plants, this study focused on conducting an RNA-seq meta-analysis to investigate and compare how C3 and C4 plants respond to drought stress at the gene expression level in their leaves. Additionally, the accuracy of the meta-analysis results was confirmed with RT-qPCR. Based on the functional enrichment and network analysis, hub genes related to ribosomal proteins and photosynthesis were found to play a potential role in stress response. Moreover, our findings suggest that the low abundant amino acid degradation pathway, possibly through providing ATP source for the TCA cycle, in both groups of plants and the activation of the OPPP pathway in C4 plants, through providing the electron source required by this plant, can help to improve drought stress tolerance.

Halophytes as new model plant species for salt tolerance strategies

Soil salinity is becoming a growing issue nowadays, severely affecting the world's most productive agricultural landscapes. With intersecting and competitive challenges of shrinking agricultural lands and increasing demand for food, there is an emerging need to build resilience for adaptation to anticipated climate change and land degradation. This necessitates the deep decoding of a gene pool of crop plant wild relatives which can be accomplished through salt-tolerant species, such as halophytes, in order to reveal the underlying regulatory mechanisms. Halophytes are generally defined as plants able to survive and complete their life cycle in highly saline environments of at least 200-500 mM of salt solution. The primary criterion for identifying salt-tolerant grasses (STGs) includes the presence of salt glands on the leaf surface and the Na⁺ exclusion mechanism since the interaction and replacement of Na⁺ and K⁺ greatly determines the survivability of STGs in saline environments. During the last decades or so, various salt-tolerant grasses/halophytes have been explored for the mining of salt-tolerant genes and testing their efficacy to improve the limit of salt tolerance in crop plants. Still, the utility of halophytes is limited due to the non-availability of any model halophytic plant system as well as the lack of complete genomic information. To date, although Arabidopsis (Arabidopsis thaliana) and salt cress (Thellungiella halophila) are being used as model plants in most salt tolerance studies, these plants are short-lived and can tolerate salinity for a shorter duration only. Thus, identifying the unique genes for salt tolerance pathways in halophytes and their introgression in a related cereal genome for better tolerance to salinity is the need of the hour. Modern technologies including RNA sequencing and genome-wide mapping along with advanced bioinformatics programs have advanced the decoding of the whole genetic information of plants and the development of probable algorithms to correlate stress tolerance limit and yield potential. Hence, this article has been compiled to explore the naturally occurring halophytes as potential model plant species for abiotic stress tolerance and to further breed crop plants to enhance salt tolerance through genomic and molecular tools.

Genome-wide identification of lncRNAs, miRNAs, mRNAs and their regulatory networks involved in tapping panel dryness in rubber tree (Hevea brasiliensis)

Noncoding RNAs (ncRNAs) play pivotal roles in various biological processes in plants. However, the role of ncRNAs in tapping panel dryness (TPD) of rubber tree (Hevea brasiliensis Muell. Arg.) is largely unknown. Here, the whole transcriptome analyses of bark tissues from healthy and TPD trees were performed to identify differentially expressed long ncRNAs (DELs), microRNAs/miRNAs (DEMs), genes (DEGs) and their regulatory networks involved in TPD. A total of 263 DELs, 174 DEMs and 1574 DEGs were identified in the bark of TPD tree compared with that of healthy tree. Kyoto Encyclopedia of Genes and Genomes analysis revealed that most of the DEGs and targets of DELs and DEMs were mainly enriched in metabolic pathways, biosynthesis of secondary metabolites and plant hormone signal transduction. Additionally, the majority of DEGs and DELs related to rubber biosynthesis were downregulated in TPD trees. Furthermore, 98 DEGs and 44 DELs were targeted by 54 DEMs, 190 DEGs were identified as putative targets of 56 DELs, and 2 and 44 DELs were predicted as precursors and endogenous target mimics of 2 and 6 DEMs, respectively. Based on these, the DEL-DEM-DEG regulatory network involved in TPD was constructed, and 13 hub DELs, 3 hub DEMs and 2 hub DEGs were identified. The results provide novel insights into the regulatory roles of ncRNAs underlying TPD and lay a foundation for future functional characterization of long ncRNAs, miRNAs and genes involved in TPD in rubber tree.

Meta-Analysis of Common and Differential Transcriptomic Responses to Biotic and Abiotic Stresses in Arabidopsis thaliana

Environmental stresses adversely affect crop growth and yield, resulting in major losses to plants. These stresses occur simultaneously in nature, and we therefore conducted a meta-analysis in this study to identify differential and shared genes, pathways, and transcriptomic mechanisms involved in Arabidopsis response to biotic and abiotic stresses. The results showed a total of 436/21 significant up-/downregulated differentially expressed genes (DEGs) in response to biotic stresses, while 476 and 71 significant DEGs were respectively up- and downregulated in response to abiotic stresses in Arabidopsis thaliana. In addition, 21 DEGs (2.09%) were commonly regulated in response to biotic and abiotic stresses. Except for WRKY45 and ATXTH22, which were respectively up-/down- and down-/upregulated in response to biotic and abiotic stresses, other common DEGs were upregulated in response to all biotic and abiotic treatments. Moreover, the transcription factors (TFs) bHLH, MYB, and WRKY were the common TFs in response to biotic and abiotic stresses. In addition, ath-miR414 and ath-miR5658 were identified to be commonly expressed in response to both biotic and abiotic stresses. The identified common genes and pathways during biotic and abiotic stresses may provide potential candidate targets for the development of stress resistance breeding programs and for the genetic manipulation of crop plants.

Identification, evolution, expression analysis of phospholipase D (PLD) gene family in tea (Camellia sinensis)

Phospholipase D (PLD) (EC 3.1.4.4) plays important roles in plants growth, development, and response to environmental stresses. Tea plant (Camellia sinensis) is the most important non-alcoholic beverage in the world with health benefits, but tea production decreases in response to environmental stresses such as cold and drought. Therefore, a genome-wide analysis of the C. sinensis PLD gene family (CsPLDs) was carried out. In the current study, identification, evolutionary relationship, duplication, selection pressure, gene structure, promoter analysis, transcript-targeted miRNA, and simple sequence repeat markers prediction, RNA-seq data analysis, and three-dimensional structure of the CsPLDs have been investigated using bioinformatics tools. 15 PLDs were identified from the tea genome which belongs to five groups, including CsPLDα, CsPLDβ, CsPLDδ, CsPLDε, and CsPLDζ. Both segmental and tandem duplications have occurred in the CsPLD gene family. Ka/Ks ratio for the most duplicated pair genes was less than 1 which implies negative selection to conserve their function during the tea evolution. 68 cis-elements have been found in CsPLDs indicating the contribution of these genes in response to environmental stresses. Likewise, 72 SSR loci and 96 miRNA molecules in 14 and 15 CsPLDs have been detected. According to RNA-seq data, the highest expression in all tissues under various abiotic stresses was related to CsPLDα1. Besides, a three-dimensional structure of the CsPLDα1 was evaluated to better understand its biological activity. This research provides comprehensive information that could be useful in future studies to develop stress-tolerant tea.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01007-0.

ARF family identification in Tamarix chinensis reveals the salt responsive expression of TcARF6 targeted by miR167

Auxin response factors (ARFs) are important transcription factors (TFs) that are differentially expressed in response to various abiotic stresses. The important roles of ARFs and small RNA-ARF pathways in mediating plant growth and stress responses have emerged in several recent studies. However, no studies on the involvement of ARFs in tamarisk trees, which are resistant to salinity, have been conducted. In this study, systematic analysis revealed 12 TcARF genes belonging to five different groups in Tamarix chinensis. The microRNA response elements of miR160, which belongs to group I and miR167, which belongs to group III, were conserved in terms of their location and sequence. Moreover, digital gene expression profiles suggested that a potential miR167 target gene, TcARF6, was rapidly expressed in response to salt stress. Cloning of TcARF6 revealed that TcARF6 could be an activation TF with a glutamine-rich region and expression pattern analysis revealed that the expression of TcARF6 was significantly downregulated specifically in the roots. A significant negative correlation in the expression pattern of tch-miR167/TcARF6 indicated that this module may play a key role in the response to salt stress. Overall, these results provide basic information on the posttranscriptional regulation of TcARF6 for future investigations of the T. chinensis salt-stress response.

Long noncoding RNAs and miRNAs regulating terpene and tartaric acid biosynthesis in rose-scented geranium

This study aimed to explore the noncoding RNAs, which have emerged as key regulatory molecules in biological processes, in rose-scented geranium. We analyzed RNA-seq data revealing 26 784 long noncoding RNAs (lncRNAs) and 871 miRNAs in rose-scented geranium. A total of 466 lncRNAs were annotated using different plant lncRNA public databases. Furthermore, 372 lncRNAs and 99 miRNAs were detected that target terpene and tartarate biosynthetic pathways. An interactome, comprising of lncRNAs, miRNAs, and mRNAs, was constructed that represents a noncoding RNA regulatory network of the target mRNAs. Real-time quantitative PCR expression validation was done for selected lncRNAs involved in the regulation of terpene and tartaric acid pathways. This study provides the first insights into the regulatory functioning of noncoding RNAs in rose-scented geranium.