The role of small RNAs in abiotic stress

MicroRNAs in Plants Development and Stress Resistance

Plant growth and development are governed by a rigorously timed sequence of ontogenetic programmes. MicroRNAs (miRNAs), a class of short noncoding RNAs, function as master regulators of gene expression by targeting mRNAs for cleavage or direct translational inhibition at the posttranscriptional level in eukaryotes. Numerous miRNA molecules that control significant agronomic properties in plants have been found. On the one hand, miRNAs target transcription factors (TFs) to determine plant structure, such as root development, internode elongation, leaf morphogenesis, sex determination and nutrient transition. On the other hand, miRNAs alter expression levels to adapt to biological and abiotic stresses, including fungi, bacteria, viruses, drought, waterlogging, high temperature, low temperature, salinity, nutrient deficiencies, heavy metals and other abiotic stresses. To fully understand the role of miRNAs in plants, we review the regulatory role of miRNAs in plant development and stress resistance. Beyond that, we propose that the novel miRNA in review can be effectively further studied with artificial miRNA (amiRNA) or short tandem target mimics (STTM) and miRNA delivery in vitro can be used to improve crop yield and agricultural sustainability.

Evaluation of zma-miR408 and its target genes function on maize (Zea mays) leaf growth response to cold stress by VIGS-based STTM approach

miR408 is a conserved plant miRNA family that is known to regulate genes involved in copper metabolism. However, the function of miR408 in maize leaf growth regulation under cold stress isn't defined. In this study, endogenous maize miR408 was transiently silenced by using virus-induced gene silencing (VIGS) combined with short tandem target mimic (STTM) approaches. To this end, STTM-miR408a/b was designed, synthesized, and applied to maize seedlings. Subsequently, STTM-miR408a/b (STTM) and mock-treated (M) seedlings were subjected to cold stress (C) and the growth response of the seedlings was monitored. Finally, STTM-miR408a/b-treatment successfully downregulated the expression of endogenous mir408a/b and upregulated their putative targets Basic Blue Protein (BBP) and Blue Copper Protein (BCP) antagonistically in the STTM and STTM_C groups compared to M and M_C groups. On the other hand, their putative target Laccase (LAC22) gene was upregulated in the STTM group compared to the M group, but there were no significant expression differences between the M_C and STTM_C groups. The elongation rate of the STTM-miR408a/b-treated second and third leaves was reduced by 10% and 19% resulting in 19% and 11% shortening, respectively. Furthermore, the activity of catalase (CAT) and glutathione reductase (GR) was decreased by 57% in STTM, M_C, and STTM_C, and 29% and 28% in the M_C and STTM_C groups and ascorbate peroxidase (APX) was increased by 15% in M_C and STTM_C groups, respectively. These findings illuminated the maize leaf growth response to cold via regulation of expression of miR408 and its target genes and antioxidant system.

Small Interfering RNAs as Critical Regulators of Plant Life Process: New Perspectives on Regulating the Transcriptomic Machinery

Small interfering RNAs (siRNAs) are a distinct class of regulatory RNAs in plants and animals. Gene silencing by small interfering RNAs is one of the fundamental mechanisms for regulating gene expression. siRNAs are critical regulators during developmental processes. siRNAs have similar structures and functions to small RNAs but are derived from double-stranded RNA and may be involved in directing DNA methylation of target sequences. siRNAs are a less well-studied class than the miRNA group, and researchers continue to identify new classes of siRNAs that appear at specific developmental stages and in particular tissues, revealing a more complex mode of siRNA action than previously thought. This review characterizes the siRNA classes and their biogenesis process and focuses on presenting their known functions in the regulation of plant development and responses to biotic and abiotic stresses. The review also highlights the exciting potential for future research in this field, proposing methods for detecting plant siRNAs and a bioinformatic pathway for identifying siRNAs and their functions.

MYB transcription factors, their regulation and interactions with non-coding RNAs during drought stress in Brassica juncea

BACKGROUND

Brassica juncea (L.) Czern is an important oilseed crop affected by various abiotic stresses like drought, heat, and salt. These stresses have detrimental effects on the crop's overall growth, development and yield. Various Transcription factors (TFs) are involved in regulation of plant stress response by modulating expression of stress-responsive genes. The myeloblastosis (MYB) TFs is one of the largest families of TFs associated with various developmental and biological processes such as plant growth, secondary metabolism, stress response etc. However, MYB TFs and their regulation by non-coding RNAs (ncRNAs) in response to stress have not been studied in B. juncea. Thus, we performed a detailed study on the MYB TF family and their interactions with miRNAs and Long non coding RNAs.

RESULTS

Computational investigation of genome and proteome data presented a comprehensive picture of the MYB genes and their protein architecture, including intron-exon organisation, conserved motif analysis, R2R3 MYB DNA-binding domains analysis, sub-cellular localization, protein-protein interaction and chromosomal locations. Phylogenetically, BjuMYBs were further classified into different subclades on the basis of topology and classification in Arabidopsis. A total of 751 MYBs were identified in B. juncea corresponding to 297 1R-BjuMYBs, 440 R2R3-BjuMYBs, 12 3R-BjuMYBs, and 2 4R-BjuMYBs types. We validated the transcriptional profiles of nine selected BjuMYBs under drought stress through RT-qPCR. Promoter analysis indicated the presence of drought-responsive cis-regulatory components. Additionally, the miRNA-MYB TF interactions was also studied, and most of the microRNAs (miRNAs) that target BjuMYBs were involved in abiotic stress response and developmental processes. Regulatory network analysis and expression patterns of lncRNA-miRNA-MYB indicated that selected long non-coding RNAs (lncRNAs) acted as strong endogenous target mimics (eTMs) of the miRNAs regulated expression of BjuMYBs under drought stress.

CONCLUSIONS

The present study has established preliminary groundwork of MYB TFs and their interaction with ncRNAs in B. juncea and it will help in developing drought- tolerant Brassica crops.

Investigating the genetic basis of maize ear characteristics: a comprehensive genome-wide study utilizing high-throughput phenotypic measurement method and system

The morphology of maize ears plays a critical role in the breeding of new varieties and increasing yield. However, the study of traditional ear-related traits alone can no longer meet the requirements of breeding. In this study, 20 ear-related traits, including size, shape, number, and color, were obtained in 407 maize inbred lines at two sites using a high-throughput phenotypic measurement method and system. Significant correlations were found among these traits, particularly the novel trait ear shape (ES), which was correlated with traditional traits: kernel number per row and kernel number per ear. Pairwise comparison tests revealed that the inbred lines of tropical-subtropical were significantly different from other subpopulations in row numbers per ear, kernel numbers per ear, and ear color. A genome-wide association study identified 275, 434, and 362 Single nucleotide polymorphisms (SNPs) for Beijing, Sanya, and best linear unbiased prediction scenarios, respectively, explaining 3.78% to 24.17% of the phenotypic variance. Furthermore, 58 candidate genes with detailed functional descriptions common to more than two scenarios were discovered, with 40 genes being associated with color traits on chromosome 1. After analysis of haplotypes, gene expression, and annotated information, several candidate genes with high reliability were identified, including Zm00001d051328 for ear perimeter and width, zma-MIR159f for ear shape, Zm00001d053080 for kernel width and row number per ear, and Zm00001d048373 for the blue color channel of maize kernels in the red-green-blue color model. This study emphasizes the importance of researching novel phenotypic traits in maize by utilizing high-throughput phenotypic measurements. The identified genetic loci enrich the existing genetic studies related to maize ears.

Role of non-coding RNAs against salinity stress in Oryza species: Strategies and challenges in analyzing miRNAs, tRFs and circRNAs

Salinity is an imbalanced concentration of mineral salts in the soil or water that causes yield loss in salt-sensitive crops. Rice plant is vulnerable to soil salinity stress at seedling and reproductive stages. Different non-coding RNAs (ncRNAs) post-transcriptionally regulate different sets of genes during different developmental stages under varying salinity tolerance levels. While microRNAs (miRNAs) are well known small endogenous ncRNAs, tRNA-derived RNA fragments (tRFs) are an emerging class of small ncRNAs derived from tRNA genes with a demonstrated regulatory role, like miRNAs, in humans but unexplored in plants. Circular RNA (circRNA), another ncRNA produced by back-splicing events, acts as target mimics by preventing miRNAs from binding with their target mRNAs, thereby reducing the miRNA's action upon its target. Same may hold true between circRNAs and tRFs. Hence, the work done on these ncRNAs was reviewed and no reports were found for circRNAs and tRFs under salinity stress in rice, either at seedling or reproductive stages. Even the reports on miRNAs are restricted to seedling stage only, in spite of severe effects on rice crop production due to salt stress during reproductive stage. Moreover, this review sheds light on strategies to predict and analyze these ncRNAs in an effective manner.

Genome-Wide Analysis of microRNAs and Their Target Genes in Dongxiang Wild Rice (Oryza rufipogon Griff.) Responding to Salt Stress

Rice (Oryza sativa) is a staple food for more than half of the world's population, and its production is critical for global food security. Moreover, rice yield decreases when exposed to abiotic stresses, such as salinity, which is one of the most detrimental factors for rice production. According to recent trends, as global temperatures continue to rise due to climate change, more rice fields may become saltier. Dongxiang wild rice (Oryza rufipogon Griff., DXWR) is a progenitor of cultivated rice and has a high tolerance to salt stress, making it useful for studying the regulatory mechanisms of salt stress tolerance. However, the regulatory mechanism of miRNA-mediated salt stress response in DXWR remains unclear. In this study, miRNA sequencing was performed to identify miRNAs and their putative target genes in response to salt stress in order to better understand the roles of miRNAs in DXWR salt stress tolerance. A total of 874 known and 476 novel miRNAs were identified, and the expression levels of 164 miRNAs were found to be significantly altered under salt stress. The stem-loop quantitative real-time PCR (qRT-PCR) expression levels of randomly selected miRNAs were largely consistent with the miRNA sequencing results, suggesting that the sequencing results were reliable. The gene ontology (GO) analysis indicated that the predicted target genes of salt-responsive miRNAs were involved in diverse biological pathways of stress tolerance. This study contributes to our understanding of DXWR salt tolerance mechanisms regulated by miRNAs and may ultimately improve salt tolerance in cultivated rice breeding using genetic methods in the future.

Comprehensive, integrative genomic analysis of microRNA expression profiles in different tissues of two wheat cultivars with different traits

Wheat is one of the most important food sources on Earth. MicroRNAs (miRNA) play important roles in wheat productivity. To identify wheat miRNAs, we constructed and sequenced sRNA libraries from leaves and roots of two wheat cultivars (RAC875 and Kukri) with many different traits. Given that available miRNA wheat complement in the plant-specific database PmiREN ( https://pmiren.com ) does not include root tissues and root-associated miRNAs might thus be missing, we performed first the prediction of novel miRNAs using the sRNAbench tool. We found a total of 150 putatively novel miRNA genes with expression of both arms from 289 unique mature sequences and nearly 30% of all miRNA reads in roots corresponded to novel miRNAs. In contrast, this figure in leaves dropped to under 3%, confirming the undersampling of roots in the complement of known miRNAs. By using 120 publicly available wheat datasets, 598 Zea mays small RNA libraries, 64 plant species genomes, wheat degradome library, and functional enrichment analysis, a subset of novel miRNAs were confirmed as bona-fide miRNAs. Of the total 605 miRNAs identified in this study inclusive of 316 known miRNAs, 528 miRNAs were shared by both cultivars, 429 miRNAs were shared by both root tissues and 329 miRNAs were shared by both leaf tissues. In addition, 32 miRNAs were specific to Kukri while 45 miRNAs were specific to RAC875. These miRNAs had diverse functions, such as regulation of gene transcription, protein translation, energy metabolism, and cell cycle progression. Our data provide a genome-wide miRNA expression profile in these two wheat cultivars and help functional studies of wheat genomics.

Reprogramming of Fundamental miRNA and Gene Expression during the Barley-Piriformospora indica Interaction

The interactions between plants and microorganisms, which are widely present in the microbial-dominated rhizosphere, have been studied. This association is highly beneficial to the organisms involved, as plants benefit soil microorganisms by providing them with metabolites, while microorganisms promote plant growth and development by promoting nutrient uptake and/or protecting the plant from biotic and abiotic stresses. Piriformospora indica, an endophytic fungus of Sebacinales, colonizes the roots of a wide range of host plants and establishes various benefits for the plants. In this work, an interaction between barley and the P. indica was established to elucidate microRNA (miRNA)-based regulatory changes in miRNA profiles and gene expression that occurred during the symbiosis. Growth promotion and vigorous root development were confirmed in barley colonized by P. indica. The genome-wide expression profile analysis of miRNAs in barley root showed that 7,798,928, 6,418,039 and 7,136,192 clean reads were obtained from the libraries of mock, 3 dai and 7 dai roots, respectively. Sequencing of the barley genome yielded in 81 novel miRNA and 450 differently expressed genes (DEGs). Additionally, 11, 24, 6 differentially expressed microRNAs (DEMs) in barley were found in the three comparison groups, including 3 dai vs. mock, 7 dai vs. mock and 7 dai vs. 3 dai, respectively. The predicted target genes of these miRNAs are mainly involved in transcription, cell division, auxin signal perception and transduction, photosynthesis and hormone stimulus. Transcriptome analysis of P. indica identified 667 and 594 differentially expressed genes (DEG) at 3 dai and 7 dai. Annotation and GO (Gene Ontology) analysis indicated that the DEGs with the greatest changes were concentrated in oxidoreductase activity, ion transmembrane transporter activity. It implies that reprogramming of fundamental miRNA and gene expression occurs both in barley and P. indica. Analysis of global changes in miRNA profiles of barley colonized with P. indica revealed that several putative endogenous barley miRNAs expressed upon colonization belonging to known micro RNA families involved in growth and developmental regulation.

ASRmiRNA: Abiotic Stress-Responsive miRNA Prediction in Plants by Using Machine Learning Algorithms with Pseudo K-Tuple Nucleotide Compositional Features

MicroRNAs (miRNAs) play a significant role in plant response to different abiotic stresses. Thus, identification of abiotic stress-responsive miRNAs holds immense importance in crop breeding programmes to develop cultivars resistant to abiotic stresses. In this study, we developed a machine learning-based computational method for prediction of miRNAs associated with abiotic stresses. Three types of datasets were used for prediction, i.e., miRNA, Pre-miRNA, and Pre-miRNA + miRNA. The pseudo K-tuple nucleotide compositional features were generated for each sequence to transform the sequence data into numeric feature vectors. Support vector machine (SVM) was employed for prediction. The area under receiver operating characteristics curve (auROC) of 70.21, 69.71, 77.94 and area under precision-recall curve (auPRC) of 69.96, 65.64, 77.32 percentages were obtained for miRNA, Pre-miRNA, and Pre-miRNA + miRNA datasets, respectively. Overall prediction accuracies for the independent test set were 62.33, 64.85, 69.21 percentages, respectively, for the three datasets. The SVM also achieved higher accuracy than other learning methods such as random forest, extreme gradient boosting, and adaptive boosting. To implement our method with ease, an online prediction server “ASRmiRNA” has been developed. The proposed approach is believed to supplement the existing effort for identification of abiotic stress-responsive miRNAs and Pre-miRNAs.

Identification and verification of seed development related miRNAs in kernel almond by small RNA sequencing and qPCR

Many studies have investigated the role of miRNAs on the yield of various plants, but so far, no report is available on the identification and role of miRNAs in fruit and seed development of almonds. In this study, preliminary analysis by high-throughput sequencing of short RNAs of kernels from the crosses between almond cultivars 'Sefid' × 'Mamaee' (with small and large kernels, respectively) and 'Sefid' × 'P. orientalis' (with small kernels) showed that the expressions of several miRNAs such as Pdu-miR395a-3p, Pdu-miR8123-5p, Pdu-miR482f, Pdu-miR6285, and Pdu-miR396a were significantly different. These miRNAs targeted genes encoding different proteins such as NYFB-3, SPX1, PGSIP3 (GUX2), GH3.9, and BEN1. The result of RT-qPCR revealed that the expression of these genes showed significant differences between the crosses and developmental stages of the seeds, suggesting that these genes might be involved in controlling kernel size because the presence of these miRNAs had a negative effect on their target genes. Pollen source can influence kernel size by affecting hormonal signaling and metabolic pathways through related miRNAs, a phenomenon known as xenia.

Comprehensive Mechanism of Gene Silencing and Its Role in Plant Growth and Development

Gene silencing is a negative feedback mechanism that regulates gene expression to define cell fate and also regulates metabolism and gene expression throughout the life of an organism. In plants, gene silencing occurs via transcriptional gene silencing (TGS) and post-transcriptional gene silencing (PTGS). TGS obscures transcription via the methylation of 5' untranslated region (5'UTR), whereas PTGS causes the methylation of a coding region to result in transcript degradation. In this review, we summarized the history and molecular mechanisms of gene silencing and underlined its specific role in plant growth and crop production.

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.

The role of differentially expressed salivary microRNA in oral squamous cell carcinoma. A systematic review

OBJECTIVE

This study aims to systematically review the role of differentially expressed microRNA (miRNA) in saliva as potential biomarkers in oral cancer patients.

DESIGN

PubMed, Scopus and EBSCO online data bases were used as well as manual searching to extract studies from January 2008 up to October 2020.

RESULTS

A total of 14 studies that met the eligibility criteria were included. All selected studies were of case-control type. A total of 25 differentially expressed miRNAs were identified. Thirteen of these miRNAs (Let-7a, miR 27, miR 34, miR 92, miR 124, miR 125a, miR 136, miR139 miR 145, miR 146a, miR 200a, miR 205 and miR 375) were downregulated and other twelve (miR 9, miR 21, miR 31, miR 122, miR 134, miR 184, miR 191, miR 196a, miR 196b, miR 412, miR 512 and miR 8392) were upregulated. Four miRNAs were evaluated in more than one study (miR21, miR31, miR125 and miR 200).

CONCLUSION

According to these results, salivary miRNA can aid in diagnosis and prognosis of oral squamous cell carcinoma (OSCC). However, controlled clinical trials with a large sample size are required to validate the differentially expressed miRNAs of the present review.

Pi-starvation induced transcriptional changes in barley revealed by a comprehensive RNA-Seq and degradome analyses

BACKGROUND

Small RNAs (sRNAs) are 20-30 nt regulatory elements which are responsible for plant development regulation and participate in many plant stress responses. Insufficient inorganic phosphate (Pi) concentration triggers plant responses to balance the internal Pi level.

RESULTS

In this study, we describe Pi-starvation-responsive small RNAs and transcriptome changes in barley (Hordeum vulgare L.) using Next-Generation Sequencing (NGS) RNA-Seq data derived from three different types of NGS libraries: (i) small RNAs, (ii) degraded RNAs, and (iii) functional mRNAs. We find that differentially and significantly expressed miRNAs (DEMs, Bonferroni adjusted p-value < 0.05) are represented by 15 molecules in shoot and 13 in root; mainly various miR399 and miR827 isomiRs. The remaining small RNAs (i.e., those without perfect match to reference sequences deposited in miRBase) are considered as differentially expressed other sRNAs (DESs, p-value Bonferroni correction < 0.05). In roots, a more abundant and diverse set of other sRNAs (DESs, 1796 unique sequences, 0.13% from the average of the unique small RNA expressed under low-Pi) contributes more to the compensation of low-Pi stress than that in shoots (DESs, 199 unique sequences, 0.01%). More than 80% of differentially expressed other sRNAs are up-regulated in both organs. Additionally, in barley shoots, up-regulation of small RNAs is accompanied by strong induction of two nucleases (S1/P1 endonuclease and 3'-5' exonuclease). This suggests that most small RNAs may be generated upon nucleolytic cleavage to increase the internal Pi pool. Transcriptomic profiling of Pi-starved barley shoots identifies 98 differentially expressed genes (DEGs). A majority of the DEGs possess characteristic Pi-responsive cis-regulatory elements (P1BS and/or PHO element), located mostly in the proximal promoter regions. GO analysis shows that the discovered DEGs primarily alter plant defense, plant stress response, nutrient mobilization, or pathways involved in the gathering and recycling of phosphorus from organic pools.

CONCLUSIONS

Our results provide comprehensive data to demonstrate complex responses at the RNA level in barley to maintain Pi homeostasis and indicate that barley adapts to Pi-starvation through elicitation of RNA degradation. Novel P-responsive genes were selected as putative candidates to overcome low-Pi stress in barley plants.

Genome-Wide Identification of Copper Stress-Regulated and Novel MicroRNAs in Mulberry Leaf

Copper (Cu) is an essential trace element for plant growth and development. It is widely involved in respiration, photosynthesis, pollen formation, and other biological processes. Therefore, low or excessive copper causes damage to plants. Mulberry is an essential perennial economic tree. At present, research on the abiotic stress responses in mulberry is mainly focused on the identification of resistant germplasm resources and cloning of resistant genes. In contrast, studies on the resistance function of microRNAs and the regulatory gene responses to stress are rare. In this study, small RNA libraries (control and copper stressed) were constructed from mulberry leaf RNA. High-throughput sequencing and screening were employed, a total of 65 known miRNAs and 78 predicted novel mature miRNAs were identified, among which 40 miRNAs were differentially expressed under copper stress. Subsequently, expression patterns were verified for 14 miRNAs by real-time fluorescence quantitative PCR (qPCR). The target genes of miRNAs were validated by 5' RLM-RACE. Our results provide the bases for further study on the molecular mechanism of copper stress regulation in mulberry.

Identification, evolution, expression, and docking studies of fatty acid desaturase genes in wheat (Triticum aestivum L.)

Backgrounds

Fatty acid desaturases (FADs) introduce a double bond into the fatty acids acyl chain resulting in unsaturated fatty acids that have essential roles in plant development and response to biotic and abiotic stresses. Wheat germ oil, one of the important by-products of wheat, can be a good alternative for edible oils with clinical advantages due to the high amount of unsaturated fatty acids. Therefore, we performed a genome-wide analysis of the wheat FAD gene family (TaFADs).

Results

68 FAD genes were identified from the wheat genome. Based on the phylogenetic analysis, wheat FADs clustered into five subfamilies, including FAB2, FAD2/FAD6, FAD4, DES/SLD, and FAD3/FAD7/FAD8. The TaFADs were distributed on chromosomes 2A-7B with 0 to 10 introns. The Ka/Ks ratio was less than one for most of the duplicated pair genes revealed that the function of the genes had been maintained during the evolution. Several cis-acting elements related to hormones and stresses in the TaFADs promoters indicated the role of these genes in plant development and responses to environmental stresses. Likewise, 72 SSRs and 91 miRNAs in 36 and 47 TaFADs have been identified. According to RNA-seq data analysis, the highest expression in all developmental stages and tissues was related to TaFAB2.5, TaFAB2.12, TaFAB2.15, TaFAB2.17, TaFAB2.20, TaFAD2.1, TaFAD2.6, and TaFAD2.8 genes while the highest expression in response to temperature stress was related to TaFAD2.6, TaFAD2.8, TaFAB2.15, TaFAB2.17, and TaFAB2.20. Furthermore, docking simulations revealed several residues in the active site of TaFAD2.6 and TaFAD2.8 in close contact with the docked oleic acid that could be useful in future site-directed mutagenesis studies to increase the catalytic efficiency of them and subsequently improve agronomic quality and tolerance of wheat against environmental stresses.

Conclusions

This study provides comprehensive information that can lead to the detection of candidate genes for wheat genetic modification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07199-1.

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.

The lre-miR159a-LrGAMYB pathway mediates resistance to grey mould infection in Lilium regale

Grey mould is one of the most determinative factors of lily growth and plays a major role in limiting lily productivity. MicroRNA159 (miR159) is a highly conserved microRNA in plants, and participates in the regulation of plant development and stress responses. Our previous studies revealed that lre-miR159a participates in the response of Lilium regale to Botrytis elliptica according to deep sequencing analyses; however, the response mechanism remains unknown. Here, lre-miR159a and its target LrGAMYB gene were isolated from L. regale. Transgenic Arabidopsis overexpressing lre-MIR159a exhibited larger leaves and smaller necrotic spots on inoculation with Botrytis than those of wild-type and overexpressing LrGAMYB plants. The lre-MIR159a overexpression also led to repressed expression of two targets of miR159, AtMYB33 and AtMYB65, and enhanced accumulation of hormone-related genes, including AtPR1, AtPR2, AtNPR1, AtPDF1.2, and AtLOX for both the jasmonic acid and salicylic acid pathways. Moreover, lower levels of H₂ O₂ and O2- were observed in lre-MIR159a transgenic Arabidopsis, which reduced the damage from reactive oxygen species accumulation. Taken together, these results indicate that lre-miR159a positively regulates resistance to grey mould by repressing the expression of its target LrGAMYB gene and activating a defence response.

Identification and expression analysis of miRNAs and elucidation of their role in salt tolerance in rice varieties susceptible and tolerant to salinity

Soil salinization is a serious problem for cultivation of rice, as among cereals rice is the most salt sensitive crop, and more than 40% of the total agricultural land amounting to approximately 80 million ha the world over is salt affected. Salinity affects a plant in a varieties of ways, including ion toxicity, osmotic stress and oxidative damage. Since miRNAs occupy the top place in biochemical events determining a trait, understanding their role in salt tolerance is highly desirable, which may allow introduction of the trait in the rice cultivars of choice through biotechnological interventions. High throughput sequencing of sRNAs in the root and shoot tissues of the seedlings of the control and NaCl treated Pokkali, a salt-tolerant rice variety, identified 75 conserved miRNAs and mapped 200 sRNAs to the rice genome as novel miRNAs. Expression of nine novel miRNAs and two conserved miRNAs were confirmed by Northern blotting. Several of both conserved and novel miRNAs that expressed differentially in root and/or shoot tissues targeted transcription factors like AP2/EREBP domain protein, ARF, NAC, MYB, NF-YA, HD-Zip III, TCP and SBP reported to be involved in salt tolerance or in abiotic stress tolerance in general. Most of the novel miRNAs expressed in the salt tolerant wild rice Oryza coarctata, suggesting conservation of miRNAs in taxonomically related species. One of the novel miRNAs, osa-miR12477, also targeted L-ascorbate oxidase (LAO), indicating build-up of oxidative stress in the plant upon salt treatment, which was confirmed by DAB staining. Thus, salt tolerance might involve miRNA-mediated regulation of 1) cellular abundance of the hormone signaling components like EREBP and ARF, 2) synthesis of abiotic stress related transcription factors, and 3) antioxidative component like LAO for mitigation of oxidative damage. The study clearly indicated importance of osa-miR12477 regulated expression of LAO in salt tolerance in the plant.

Stress-responsive miRNAome of Glycine max (L.) Merrill: molecular insights and way forward

MAIN CONCLUSION

Analysis of stress-associated miRNAs of Glycine max (L.) Merrill reveals wider ramifications of small RNA-mediated (conserved and legume-specific miRNAs) gene regulatory foot prints in molecular adaptive responses. MicroRNAs (miRNAs) are indispensable components of gene regulatory mechanism of plants. Soybean is a crop of immense commercial potential grown worldwide for its edible oil and soy meal. Intensive research efforts, using the next generation sequencing and bioinformatics techniques, have led to the identification and characterization of numerous small RNAs, especially microRNAs (miRNAs), in soybean. Furthermore, studies have unequivocally demonstrated the significance of miRNAs during the developmental processes and various stresses in soybean. In this review, we summarize the current state of understanding of miRNA-based abiotic and biotic stress responses in soybean. In addition, the molecular insights gained from the stress-related soybean miRNAs have been compared to the miRNAs of other crops, especially legumes, and the core commonalities have been highlighted, though differences among them were not ignored. Nature of response of soybean-derived conserved miRNAs during various stresses was also analyzed to gain deeper insights regarding sRNAome-based defense responses. This review further provides way forward in legume small RNA transcriptomics based on the adaptive responses of soybean and other legume-derived miRNAs.

Identification of copper (Cu) stress-responsive grapevine microRNAs and their target genes by high-throughput sequencing

MicroRNAs (miRNAs) are a class of single-stranded non-coding small RNAs (sRNAs) that are 20-24 nucleotides (nt) in length. Extensive studies have indicated that miRNAs play important roles in plant growth, development and stress responses. With more copper (Cu) and copper containing compounds used as bactericides and fungicides in plants, Cu stress has become one of the serious environmental problems that affect plant growth and development. In order to uncover the hidden response mechanisms of Cu stress, two small RNA libraries were constructed from Cu-treated and water-treated (Control) leaves of 'Summer Black' grapevine. Following high-throughput sequencing and filtering, a total of 158 known and 98 putative novel miRNAs were identified in the two libraries. Among these, 100 known and 47 novel miRNAs were identified as differentially expressed under Cu stress. Subsequently, the expression patterns of nine Cu-responsive miRNAs were validated by quantitative real-time PCR (qRT-PCR). There existed some consistency in expression levels of Cu-responsive miRNAs between high throughput sequencing and qRT-PCR assays. The targets prediction of miRNAs indicates that miRNA may regulate some transcription factors, including AP2, SBP, NAC, MYB and ARF during Cu stress. The target genes for two known and two novel miRNAs showed specific cleavage sites at the 10th and/or 11th nucleotide from the 5'-end of the miRNA corresponding to their miRNA complementary sequences. The findings will lay the foundation for exploring the role of the regulation of miRNAs in response to Cu stress and provide valuable gene information for breeding some Cu-tolerant grapevine cultivars.

Identification of a redox-dependent regulatory network of miRNAs and their targets in wheat

Reactive oxygen species and antioxidants have an important role in the regulation of plant growth and development under both optimal and stress conditions. In this study, we investigate a possible redox control of miRNAs in wheat (Triticum aestivum ssp. aestivum). Treatment of seedlings with 10 mM H2O2 via the roots for 24 h resulted in decreased glutathione content, increased half-cell reduction potential of the glutathione disulphide/glutathione redox pair, and greater ascorbate peroxidase activity compared to the control plants. These changes were accompanied by alterations in the miRNA transcript profile, with 70 miRNAs being identified with at least 1.5-fold difference in their expression between control and treated (0, 3, 6 h) seedlings. Degradome sequencing identified 86 target genes of these miRNAs, and 6722 possible additional target genes were identified using bioinformatics tools. The H2O2-responsiveness of 1647 target genes over 24 h of treatment was also confirmed by transcriptome analysis, and they were mainly found to be related to the control of redox processes, transcription, and protein phosphorylation and degradation. In a time-course experiment (0-24 h of treatment) a correlation was found between the levels of glutathione, other antioxidants, and the transcript levels of the H2O2-responsive miRNAs and their target mRNAs. This relationship together with bioinformatics modelling of the regulatory network indicated glutathione-related redox control of miRNAs and their targets, which allows the adjustment of the metabolism to changing environmental conditions.

The miRNAome of ramie (Boehmeria nivea L.): identification, expression, and potential roles of novel microRNAs in regulation of cadmium stress response

BACKGROUND

MicroRNAs (miRNAs) regulate numerous crucial abiotic stress processes in plants. However, information is limited on their involvement in cadmium (Cd) stress response and tolerance mechanisms in plants, including ramie (Boehmeria nivea L.) that produces a number of economic valuable as an important natural fibre crop and an ideal crop for Cd pollution remediation.

RESULTS

Four small RNA libraries of Cd-stressed and non-stressed leaves and roots of ramie were constructed. Using small RNA-sequencing, 73 novel miRNAs were identified. Genome-wide expression analysis revealed that a set of miRNAs was differentially regulated in response to Cd stress. In silico target prediction identified 426 potential miRNA targets that include several uptake or transport factors for heavy metal ions. The reliability of small RNA sequencing and the relationship between the expression levels of miRNAs and their target genes were confirmed by quantitative PCR (q-PCR). We showed that the expression patterns of miRNAs obtained by q-PCR were consistent with those obtained from small RNA sequencing. Moreover, we demonstrated that the expression of six randomly selected target genes was inversely related to that of their corresponding miRNAs, indicating that the miRNAs regulate Cd stress response in ramie.

CONCLUSIONS

This study enriches the number of Cd-responsive miRNAs and lays a foundation for the elucidation of the miRNA-mediated regulatory mechanism in ramie during Cd stress.

Regulation through MicroRNAs in Response to Low-Energy N+ Ion Irradiation in Oryza sativa

MicroRNAs (miRNAs) are a non-coding regulatory RNAs that play significant roles in plant growth and development, especially in the stress response. Low-energy ion radiation, a type of environmental stress, can cause multiple biological effects. To understand the roles of miRNAs in response to low-energy N⁺ ion radiation in Oryza sativa, high-throughput sequencing of small RNAs was carried out to detect the expression of miRNAs in the shoots of the rice after 2 × 10¹⁷ N+/cm² irradiation. The differentially expressed 28 known miRNAs were identified, 17 of these identified miRNAs were validated by real-time quantitative fluorescent PCR (q-PCR), including 9 up-regulated miRNAs (miR1320-3p, miR1320-5p, miR156b-3p, miR156c-5p, miR156c-3p/g-3p, miR1561-5p, miR398b and miR6250) and 8 down-regulated miRNAs (miR156a/e/i, miR156k, miR160f-5p, miR166j-5p, miR1846e and miR399d). In addition, 45 novel radiationresponsive miRNAs were predicted, and 8 of them were verified by q-PCR. The target genes of radiation-responsive miRNAs were predicted and gene function enrichment analysis was performed with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). The expression of 9 targets of 4 known miRNA families (miR156, miR399, miR1320 and miR398) and 2 targets of 2 novel miRNAs were quantified by q-PCR, and a strong negative regulation relation between miRNAs and their targets were observed. Those targets including SQUAMOSA promoterbinding-like protein (SPL) genes, copper/zinc superoxide dismutase (Cu/Zn-SOD), copper chaperone for SOD (CCS1) and electron transporter/ heat-shock protein binding protein (HSP), which are involved in growth and defense against various stresses, especially associated with reactive oxygen species (ROS) scavenging. This work provides important information for understanding the ROS generation and elimination mechanisms closely related to miRNAs in rice seedlings after low-energy N⁺ radiation exposure.

Ascophyllum nodosum extract mitigates salinity stress in Arabidopsis thaliana by modulating the expression of miRNA involved in stress tolerance and nutrient acquisition

Ascophyllum nodosum extract (ANE) contains bioactive compounds that improve the growth of Arabidopsis in experimentally-induced saline conditions; however, the molecular mechanisms through which ANE elicits tolerance to salinity remain largely unexplored. Micro RNAs (miRNAs) are key regulators of gene expression, playing crucial roles in plant growth, development, and stress tolerance. Next generation sequencing of miRNAs from leaves of control Arabidopsis and from plants subjected to three treatments (ANE, NaCl and ANE+NaCl) was used to identify ANE-responsive miRNA in the absence and presence of saline conditions. Differential gene expression analysis revealed that ANE had a strong effect on miRNAs expression in both conditions. In the presence of salinity, ANE tended to reduce the up-regulation or the down-regulation trend induced caused by NaCl in miRNAs such as ath-miR396a-5p, ath-miR399, ath-miR2111b and ath-miR827. To further uncover the effects of ANE, the expression of several target genes of a number of ANE-responsive miRNAs was analyzed by qPCR. NaCl, but not ANE, down-regulated miR396a-5p, which negatively regulated the expression of AtGRF7 leading to a higher expression of AtDREB2a and AtRD29 in the presence of ANE+NaCl, as compared to ANE alone. ANE+NaCl initially reduced and then enhanced the expression of ath-miR169g-5p, while the expression of the target genes AtNFYA1 and ATNFYA2, known to be involved in the salinity tolerance mechanism, was increased as compared to ANE or to NaCl treatments. ANE and ANE+NaCl modified the expression of ath-miR399, ath-miR827, ath-miR2111b, and their target genes AtUBC24, AtWAK2, AtSYG1 and At3g27150, suggesting a role of ANE in phosphate homeostasis. In vivo and in vitro experiments confirmed the improved growth of Arabidopsis in presence of ANE, in saline conditions and in phosphate-deprived medium, further substantiating the influence of ANE on a variety of essential physiological processes in Arabidopsis including salinity tolerance and phosphate uptake.

Perspectives on Structural, Physiological, Cellular, and Molecular Responses to Desiccation in Resurrection Plants

Resurrection plants possess a unique ability to counteract desiccation stress. Desiccation tolerance (DT) is a very complex multigenic and multifactorial process comprising a combination of physiological, morphological, cellular, genomic, transcriptomic, proteomic, and metabolic processes. Modification in the sugar composition of the hemicellulosic fraction of the cell wall is detected during dehydration. An important change is a decrease of glucose in the hemicellulosic fraction during dehydration that can reflect a modification of the xyloglucan structure. The expansins might also be involved in cell wall flexibility during drying and disrupt hydrogen bonds between polymers during rehydration of the cell wall. Cleavages by xyloglucan-modifying enzymes release the tightly bound xyloglucan-cellulose network, thus increasing cell wall flexibility required for cell wall folding upon desiccation. Changes in hydroxyproline-rich glycoproteins (HRGPs) such as arabinogalactan proteins (AGPs) are also observed during desiccation and rehydration processes. It has also been observed that significant alterations in the process of photosynthesis and photosystem (PS) II activity along with changes in the antioxidant enzyme system also increased the cell wall and membrane fluidity resulting in DT. Similarly, recent data show a major role of ABA, LEA proteins, and small regulatory RNA in regulating DT responses. Current progress in "-omic" technologies has enabled quantitative monitoring of the plethora of biological molecules in a high throughput routine, making it possible to compare their levels between desiccation-sensitive and DT species. In this review, we present a comprehensive overview of structural, physiological, cellular, molecular, and global responses involved in desiccation tolerance.

The miRNAome dynamics during developmental and metabolic reprogramming of tomato root infected with potato cyst nematode

Cyst-forming plant-parasitic nematodes are pests threatening many crops. By means of their secretions cyst nematodes induce the developmental and metabolic reprogramming of host cells that lead to the formation of a syncytium, which is the sole food source for growing nematodes. The in depth micro RNA (miRNA) dynamics in the syncytia induced by Globodera rostochiensis in tomato roots was studied. The miRNAomes were obtained from syncytia covering the early and intermediate developmental stages, and were the subject of differential expression analysis. The expression of 1235 miRNAs was monitored. The fold change (log₂FC) ranged from -7.36 to 8.38, indicating that this transcriptome fraction was very variable. Moreover, we showed that the DE (differentially expressed) miRNAs do not fully overlap between the selected time points, suggesting infection stage specific regulation by miRNA. The correctness of RNA-seq expression profiling was confirmed by qRT-PCR (quantitative Real Time Polymerase Chain Reaction) for seven miRNA species. Down- and up-regulated miRNA species, including their isomiRs, were further used to identify their potential targets. Among them there are a large number of transcription factors linked to different aspects of plant development belonging to gene families, such as APETALA2 (AP2), SQUAMOSA (MADS-box), MYB, GRAS, and AUXIN RESPONSE FACTOR (ARF). The substantial portion of potential target genes belong to the NB-LRR and RLK (RECEPTOR-LIKE KINASE) families, indicating the involvement of miRNA mediated regulation in defense responses. We also collected the evidence for target cleavage in the case of 29 miRNAs using one of three alternative methods: 5' RACE (5' Rapid Amplification of cDNA Ends), a search of tasiRNA within our datasets, and the meta-analysis of tomato degradomes in the GEO (Gene Expression Omnibus) database. Eight target transcripts showed a negative correlation with their respective miRNAs at two or three time points. These results indicate a large regulatory potential for miRNAs in tuning the development and defense responses.

Chinese Wheat Mosaic Virus-Induced Gene Silencing in Monocots and Dicots at Low Temperature

Virus-induced gene silencing (VIGS) is an important tool for functional genomics studies in plants. With this method, it is possible to target most endogenous genes and downregulate the messenger RNA (mRNA) in a sequence-specific manner. Chinese wheat mosaic virus (CWMV) has a bipartite, single-strand positive RNA genome, and can infect both wheat and Nicotiana benthamiana, and the optimal temperature for systemic infection in plants is 17°C. To assess the potential of the virus as a vector for gene silencing at low temperature, a fragment of the N. benthamiana or wheat phytoene desaturase (PDS) gene was expressed from a modified CWMV RNA2 clone and the resulting photo bleaching in infected plants was used as a reporter for silencing. Downregulation of PDS mRNA was also measured by quantitative reverse-transcriptase polymerase chain reaction (RT-qPCR). In experiments using fragments of PDS ranging from 500 to 1500 nucleotides, insert length influenced the stability and the efficiency of VIGS. The CWMV induced silencing system was also used to suppress miR165/166 and miR3134a through expression of miRNA target mimics. The relative expression levels of mature miR165/166 and miR3134a decreased whereas the transcript levels of their target genes increased. Interestingly, we also found the CWMV-induced silencing system was more efficient compare with the vector based on Barley stripe mosaic virus (BSMV) or Foxtail mosaic virus (FoMV) in wheat or the vector based on TRV in N. benthamiana at 17°C. In summary, the CWMV vector is effective in silencing endogenous genes and miRNAs at 17°C, thereby providing a powerful tool for gene function analysis in both N. benthamiana and wheat at low temperature.

Identification and characterization of durum wheat microRNAs in leaf and root tissues

MicroRNAs are a class of post-transcriptional regulators of plant developmental and physiological processes and responses to environmental stresses. Here, we present the study regarding the annotation and characterization of MIR genes conducted in durum wheat. We characterized the miRNAome of leaf and root tissues at tillering stage under two environmental conditions: irrigated with 100% (control) and 55% of evapotranspiration (early water stress). In total, 90 microRNAs were identified, of which 32 were classified as putative novel and species-specific miRNAs. In addition, seven microRNA homeologous groups were identified in each of the two genomes of the tetraploid durum wheat. Differential expression analysis highlighted a total of 45 microRNAs significantly differentially regulated in the pairwise comparisons leaf versus root. The miRNA families, miR530, miR395, miR393, miR5168, miR396 and miR166, miR171, miR319, and miR167, were the most expressed in leaves in comparison to roots. Putative microRNA targets were predicted for both five and three prime sequences derived from the stem-loop of the MIR gene. Gene ontology analysis showed significant overrepresented gene categories in microRNA targets belonging to transcription factors, phenylpropanoids, oxydases, and lipid binding-protein. This work represents one of the first genome wide characterization of MIR genes in durum wheat, identifying leaf and root tissue-specific microRNAs. This genomic identification of microRNAs together with the analysis of their expression profiles is a well-accepted starting point leading to a better comprehension of the role of MIR genes in the genus Triticum.

Stress-responsive pathways and small RNA changes distinguish variable developmental phenotypes caused by MSH1 loss

BACKGROUND

Proper regulation of nuclear-encoded, organelle-targeted genes is crucial for plastid and mitochondrial function. Among these genes, MutS Homolog 1 (MSH1) is notable for generating an assortment of mutant phenotypes with varying degrees of penetrance and pleiotropy. Stronger phenotypes have been connected to stress tolerance and epigenetic changes, and in Arabidopsis T-DNA mutants, two generations of homozygosity with the msh1 insertion are required before severe phenotypes begin to emerge. These observations prompted us to examine how msh1 mutants contrast according to generation and phenotype by profiling their respective transcriptomes and small RNA populations.

RESULTS

Using RNA-seq, we analyze pathways that are associated with MSH1 loss, including abiotic stresses such as cold response, pathogen defense and immune response, salicylic acid, MAPK signaling, and circadian rhythm. Subtle redox and environment-responsive changes also begin in the first generation, in the absence of strong phenotypes. Using small RNA-seq we further identify miRNA changes, and uncover siRNA trends that indicate modifications at the chromatin organization level. In all cases, the magnitude of changes among protein-coding genes, transposable elements, and small RNAs increases according to generation and phenotypic severity.

CONCLUSION

Loss of MSH1 is sufficient to cause large-scale regulatory changes in pathways that have been individually linked to one another, but rarely described all together within a single mutant background. This study enforces the recognition of organelles as critical integrators of both internal and external cues, and highlights the relationship between organelle and nuclear regulation in fundamental aspects of plant development and stress signaling. Our findings also encourage further investigation into potential connections between organelle state and genome regulation vis-á-vis small RNA feedback.

Involvement of Small RNAs in Phosphorus and Sulfur Sensing, Signaling and Stress: Current Update

Plants require several essential mineral nutrients for their growth and development. These nutrients are required to maintain physiological processes and structural integrity in plants. The root architecture has evolved to absorb nutrients from soil and transport them to other parts of the plant. Nutrient deficiency affects several physiological and biological processes in plants and leads to reduction in crop productivity and yield. To compensate this adversity, plants have developed adaptive mechanisms to enhance the acquisition, conservation, and mobilization of these nutrients under deficient or adverse conditions. In addition, plants have evolved an intricate nexus of complex signaling cascades, which help in nutrient sensing and uptake as well as to maintain nutrient homeostasis. In recent years, small non-coding RNAs such as micro RNAs (miRNAs) and endogenous small interfering RNAs have emerged as important component in regulating plant stress responses. A set of these small RNAs (sRNAs) have been implicated in regulating various processes involved in nutrient uptake, assimilation, and deficiency. In response to phosphorus (P) and sulphur (S) deficiencies, role of sRNAs, miR395 and miR399, have been identified to be instrumental; however, many more miRNAs might be involved in regulating the plant response to these nutrient stresses. These sRNAs modulate expression of target genes in response to P and S deficiencies and regulate their uptake and utilization for proper growth and development of the plant. This review summarizes the current understanding of uptake, sensing, and signaling of P and S and highlights the regulatory role of sRNAs in adaptive responses to these nutrient stresses in plants.

Identification of Salt Tolerance-related microRNAs and Their Targets in Maize (Zea mays L.) Using High-throughput Sequencing and Degradome Analysis

To identify the known and novel microRNAs (miRNAs) and their targets that are involved in the response and adaptation of maize (Zea mays) to salt stress, miRNAs and their targets were identified by a combined analysis of the deep sequencing of small RNAs (sRNA) and degradome libraries. The identities were confirmed by a quantitative expression analysis with over 100 million raw reads of sRNA and degradome sequences. A total of 1040 previously known miRNAs were identified from four maize libraries, with 762 and 726 miRNAs derived from leaves and roots, respectively, and 448 miRNAs that were common between the leaves and roots. A total of 37 potential new miRNAs were selected based on the same criteria in response to salt stress. In addition to known miR167 and miR164 species, novel putative miR167 and miR164 species were also identified. Deep sequencing of miRNAs and the degradome [with quantitative reverse transcription polymerase chain reaction (qRT-PCR) analyses of their targets] showed that more than one species of novel miRNA may play key roles in the response to salinity in maize. Furthermore, the interaction between miRNAs and their targets may play various roles in different parts of maize in response to salinity.

Characterization of miRNAs responsive to exogenous ethylene in grapevine berries at whole genome level

MicroRNAs (miRNAs) are critical regulators of various biological and metabolic processes of plants. Numerous miRNAs and their functions have been identified and analyzed in many plants. However, till now, the involvement of miRNAs in the response of grapevine berries to ethylene has not been reported yet. Here, Solexa technology was employed to deeply sequence small RNA libraries constructed from grapevine berries treated with and without ethylene. A total of 124 known and 78 novel miRNAs were identified. Among these miRNAs, 162 miRNAs were clearly responsive to ethylene, with 55 downregulated, 59 upregulated, and 14 unchanged miRNAs detected only in the control. The other 35 miRNAs responsive to ethylene were induced by ethylene and detected only in the ethylene-treated grapevine materials. Expression analysis of 27 conserved and 26 novel miRNAs revealed that 13 conserved and 18 novel ones were regulated by ethylene during the whole development of grapevine berries. High-throughput sequencing and qRT-PCR assays revealed consistent results on the expression results of ethylene-responsive miRNAs. Moreover, 90 target genes for 34 novel miRNAs were predicted, most of which were involved in responses to various stresses, especially like exogenous ethylene treatment. The identified miRNAs may be mainly involved in grapevine berry development and response to various environmental conditions.

Generation of different sizes and classes of small RNAs in barley is locus, chromosome and/or cultivar-dependent

Background

Various small RNA (sRNA) sizes and varieties have been identified, but their relationship as well as relationship with their origins and allocations have not been well understood or investigated.

Results

By comparing sRNAs generated from two barley cultivars, Golden Promise (GP) and Pallas, we identified that the generation of different sizes and types of sRNAs in barley was locus-, chromosome- and/or cultivar-dependent. 20-nt sRNAs mainly comprising miRNAs and chloroplast-derived sRNAs were significantly over-expressed in Pallas vs. GP on chromosomes 3H and 6H. MiRNAs-enriched 21-nt sRNAs were significantly over-expressed in Pallas vs. GP only on chromosome 4H. On chromosome 5H this size of sRNAs was significantly under-expressed in Pallas, so were 22-nt sRNAs mainly comprising miRNAs and repeat-derived sRNAs. 24-nt sRNAs mostly derived from repeats were evenly distributed in all chromosomes and expressed similarly between GP and Pallas. Unlike other sizes of sRNAs, 24-nt sRNAs were little conserved in other plant species. Abundant sRNAs were mostly generated from 3’ terminal regions of chromosome 1H and 5’ terminal regions of chromosome 5H. Over-expressed miRNAs in GP vs. Pallas primarily function in stress responses and iron-binding.

Conclusions

Our study indicates that 23−24-nt sRNAs may be linked to repressive chromatin modifications and function in genome stability while 20−21-nt sRNAs may be important for the cultivar specificity. This study provides a novel insight into the mechanism of sRNA expression and function in barley.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3023-5) contains supplementary material, which is available to authorized users.

Time-Course Small RNA Profiling Reveals Rice miRNAs and Their Target Genes in Response to Rice Stripe Virus Infection

It has been known that many microRNAs (miRNAs) are involved in the regulation for the plant development and defense mechanism by regulating the expression of the target gene. Several previous studies has demonstrated functional roles of miRNAs in antiviral defense mechanisms. In this study, we employed high-throughput sequencing technology to identify rice miRNAs upon rice stripe virus (RSV) infection at three different time points. Six libraries from mock and RSV-infected samples were subjected for small RNA sequencing. Bioinformatic analyses revealed 374 known miRNAs and 19 novel miRNAs. Expression of most identified miRNAs was not dramatically changed at 3 days post infection (dpi) and 7 dpi by RSV infection. However, many numbers of miRNAs were up-regulated in mock and RSV-infected samples at 15 dpi by RSV infection. Moreover, expression profiles of identified miRNAs revealed that only few numbers of miRNAs were strongly regulated by RSV infection. In addition, 15 resistance genes were targets of six miRNAs suggesting that those identified miRNAs and 15 NBS-LRR resistance genes might be involved in RSV infection. Taken together, our results provide novel insight into the dynamic expression profiles of rice miRNAs upon RSV infection and clues for the understanding of the regulatory roles of miRNAs via time-course.

Comparison of miRNAs and Their Targets in Seed Development between Two Maize Inbred Lines by High-Throughput Sequencing and Degradome Analysis

MicroRNAs (miRNAs) play an important role in plant growth, development, and response to environment. For identifying and comparing miRNAs and their targets in seed development between two maize inbred lines (i.e. PH6WC and PH4CV), two sRNAs and two degradome libraries were constructed. Through high-throughput sequencing and miRNA identification, 55 conserved and 24 novel unique miRNA sequences were identified in two sRNA libraries; moreover, through degradome sequencing and analysis, 137 target transcripts corresponding to 38 unique miRNA sequences were identified in two degradome libraries. Subsequently, 16 significantly differentially expressed miRNA sequences were verified by qRT-PCR, in which 9 verified sequences obviously target 30 transcripts mainly involved with regulation in flowering and development in embryo. Therefore, the results suggested that some miRNAs (e.g. miR156, miR171, miR396 and miR444) related reproductive development might differentially express in seed development between the PH6WC and PH4CV maize inbred lines in this present study.

Nicotine biosynthesis is regulated by two more layers: Small and long non-protein-coding RNAs

In recent years, many small RNAs and long non-protein-coding RNAs (lncRNAs) have been identified and characterized. They have been proved to play essential regulatory roles in gene expression in both primary and secondary metabolisms. In nature, many plants produce alkaloids. However, there are only few reports on the involvement of non-coding RNAs in alkaloid biosynthesis. Nicotine is major alkaloid in tobacco plants. Its biosynthesis and regulation in tobacco (Nicotiana tabacum) have been well studied; and major structural genes involved in the nicotine biosynthesis and transcriptional regulators related to its biosynthesis have been identified and characterized. In our recent studies, we identified a microRNA (nta-miRX27) and also a lncRNA (nta-eTMX27) as an endogenous target mimicry (eTM) in tobacco targeting the nicotine biosynthesis key gene QPT2 encoding quinolinate phosphoribosyltransferase (QPT) and thereby regulating the nicotine content. Their regulatory pattern leads us to conclude that nicotine biosynthesis is regulated by 2 more layers besides previously known mechanisms. Future study on the relationship between the non-coding RNAs and transcription factors in nicotine biosynthesis was discussed in this article.

Increasing abscisic acid levels by immunomodulation in barley grains induces precocious maturation without changing grain composition

Abscisic acid (ABA) accumulates in seeds during the transition to the seed filling phase. ABA triggers seed maturation, storage activity, and stress signalling and tolerance. Immunomodulation was used to alter the ABA status in barley grains, with the resulting transgenic caryopses responding to the anti-ABA antibody gene expression with increased accumulation of ABA. Calculation of free versus antibody-bound ABA reveals large excess of free ABA, increasing signficantly in caryopses from 10 days after fertilization. Metabolite and transcript profiling in anti-ABA grains expose triggered and enhanced ABA-functions such as transcriptional up-regulation of sucrose-to-starch metabolism, storage protein synthesis and ABA-related signal transduction. Thus, enhanced ABA during transition phases induces precocious maturation but negatively interferes with growth and development. Anti-ABA grains display broad constitutive gene induction related to biotic and abiotic stresses. Most of these genes are ABA- and/or stress-inducible, including alcohol and aldehyde dehydrogenases, peroxidases, chaperones, glutathione-S-transferase, drought- and salt-inducible proteins. Conclusively, ABA immunomodulation results in precocious ABA accumulation that generates an integrated response of stress and maturation. Repression of ABA signalling, occurring in anti-ABA grains, potentially antagonizes effects caused by overshooting production. Finally, mature grain weight and composition are unchanged in anti-ABA plants, although germination is somewhat delayed. This indicates that anti-ABA caryopses induce specific mechanisms to desensitize ABA signalling efficiently, which finally yields mature grains with nearly unchanged dry weight and composition. Such compensation implicates the enormous physiological and metabolic flexibilities of barley grains to adjust effects of unnaturally high ABA amounts in order to ensure and maintain proper grain development.

Nice to meet you: genetic, epigenetic and metabolic controls of plant perception of beneficial associative and endophytic diazotrophic bacteria in non-leguminous plants

A wide range of rhizosphere diazotrophic bacteria are able to establish beneficial associations with plants, being able to associate to root surfaces or even endophytically colonize plant tissues. In common, both associative and endophytic types of colonization can result in beneficial outcomes to the plant leading to plant growth promotion, as well as increase in tolerance against biotic and abiotic stresses. An intriguing question in such associations is how plant cell surface perceives signals from other living organisms, thus sorting pathogens from beneficial ones, to transduce this information and activate proper responses that will finally culminate in plant adaptations to optimize their growth rates. This review focuses on the recent advances in the understanding of genetic and epigenetic controls of plant-bacteria signaling and recognition during beneficial associations with associative and endophytic diazotrophic bacteria. Finally, we propose that "soil-rhizosphere-rhizoplane-endophytes-plant" could be considered as a single coordinated unit with dynamic components that integrate the plant with the environment to generate adaptive responses in plants to improve growth. The homeostasis of the whole system should recruit different levels of regulation, and recognition between the parties in a given environment might be one of the crucial factors coordinating these adaptive plant responses.

A cotton miRNA is involved in regulation of plant response to salt stress

The present study functionally identified a new microRNA (microRNA ovual line 5, miRNVL5) with its target gene GhCHR from cotton (Gossypium hirsutum). The sequence of miRNVL5 precursor is 104 nt long, with a well developed secondary structure. GhCHR contains two DC1 and three PHD Cys/His-rich domains, suggesting that GhCHR encodes a zinc-finger domain-containing transcription factor. miRNVL5 and GhCHR express at various developmental stages of cotton. Under salt stress (50-400 mM NaCl), miRNVL5 expression was repressed, with concomitant high expression of GhCHR in cotton seedlings. Ectopic expression of GhCHR in Arabidopsis conferred salt stress tolerance by reducing Na(+) accumulation in plants and improving primary root growth and biomass. Interestingly, Arabidopsis constitutively expressing miRNVL5 showed hypersensitivity to salt stress. A GhCHR orthorlous gene At2g44380 from Arabidopsis that can be cleaved by miRNVL5 was identified by degradome sequencing, but no confidential miRNVL5 homologs in Arabidopsis have been identified. Microarray analysis of miRNVL5 transgenic Arabidopsis showed six downstream genes (CBF1, CBF2, CBF3, ERF4, AT3G22920, and AT3G49200), which were induced by salt stress in wild-type but repressed in miRNVL5-expressing Arabidopsis. These results indicate that miRNVL5 is involved in regulation of plant response to salt stress.

Comparison of small RNA profiles in Nicotiana benthamiana and Solanum lycopersicum infected by polygonum ringspot tospovirus reveals host-specific responses to viral infection

Viral small RNAs (vsRNAs) are one of the key elements involved in RNA silencing-based defense against viruses in plants. We analyzed the vsRNA profiles in Nicotiana benthamiana and Solanum lycopersicum infected by polygonum ringspot virus (PolRSV) (Tospovirus, Bunyaviridae). VsRNAs were abundant in both hosts, but a different size profile was observed, with an abundance peak at 21 in N. benthamiana and at 22 nt in tomato. VsRNAs mapping to the PolRSV L genomic segment were under-represented in both hosts, while S and M segments were differentially and highly targeted in N. benthamiana and tomato, respectively. Differences in preferential targeting of single ORFs were observed, with over-representation of NSs ORF-derived reads in N. benthamiana. Intergenic regions (IGRs)-mapping vsRNAs were under-represented, while enrichment of vsRNAs reads mapping to the NSs positive sense strand was observed in both hosts. Comparison with a previous study on tomato spotted wilt virus (TSWV) under the same experimental conditions, showed that the relative accumulation of PolRSV-specific and endogenous sRNAs was similar to the one observed for silencing suppressor-deficient TSWV strains, suggesting possible different properties of PolRSV NSs silencing suppressor compared to that of TSWV.

Genome-Wide Analysis of MicroRNA Responses to the Phytohormone Abscisic Acid in Populus euphratica

MicroRNA (miRNA) is a type of non-coding small RNA with a regulatory function at the posttranscriptional level in plant growth development and in response to abiotic stress. Previous studies have not reported on miRNAs responses to the phytohormone abscisic acid (ABA) at a genome-wide level in Populus euphratica, a model tree for studying abiotic stress responses in woody plants. Here we analyzed the miRNA response to ABA at a genome-wide level in P. euphratica utilizing high-throughput sequencing. To systematically perform a genome-wide analysis of ABA-responsive miRNAs in P. euphratica, nine sRNA libraries derived from three groups (control, treated with ABA for 1 day and treated with ABA for 4 days) were constructed. Each group included three libraries from three individual plantlets as biological replicate. In total, 151 unique mature sequences belonging to 75 conserved miRNA families were identified, and 94 unique sequences were determined to be novel miRNAs, including 56 miRNAs with miRNA(*) sequences. In all, 31 conserved miRNAs and 31 novel miRNAs response to ABA significantly differed among the groups. In addition, 4132 target genes were predicted for the conserved and novel miRNAs. Confirmed by real-time qPCR, expression changes of miRNAs were inversely correlated with the expression profiles of their putative targets. The Populus special or novel miRNA-target interactions were predicted might be involved in some biological process related stress tolerance. Our analysis provides a comprehensive view of how P. euphratica miRNA respond to ABA, and moreover, different temporal dynamics were observed in different ABA-treated libraries.