Uncovering small RNA-mediated responses to cold stress in a wheat thermosensitive genic male-sterile line by deep sequencing

Identification and Analysis of lncRNA and circRNA Related to Wheat Grain Development

The role of lncRNA and circRNA in wheat grain development is still unclear. The objectives of this study were to characterize the lncRNA and circRNA in the wheat grain development and to construct the interaction network among lncRNA, circRNA, and their target miRNA to propose a lncRNA-circRNA-miRNA module related to wheat grain development. Full transcriptome sequencing on two wheat varieties (Annong 0942 and Anke 2005) with significant differences in 1000-grain weight at 10 d (days after pollination), 20 d, and 30 d of grain development were conducted. We detected 650, 736, and 609 differentially expressed lncRNA genes, and 769, 1054, and 1062 differentially expressed circRNA genes in the grains of 10 days, 20 days and 30 days after pollination between Annong 0942 and Anke 2005, respectively. An analysis of the lncRNA-miRNA and circRNA-miRNA targeting networks reveals that circRNAs exhibit a more complex and extensive interaction network in the development of cereal grains and the formation of grain shape. Central to these interactions are tae-miR1177, tae-miR1128, and tae-miR1130b-3p. In contrast, lncRNA genes only form a singular network centered around tae-miR1133 and tae-miR5175-5p when comparing between varieties. Further analysis is conducted on the underlying genes of all target miRNAs, we identified TaNF-YB1 targeted by tae-miR1122a and TaTGW-7B targeted by miR1130a as two pivotal regulatory genes in the development of wheat grains. The quantitative real-time PCR (qRT-PCR) and dual-luciferase reporter assays confirmed the target regulatory relationships between miR1130a-TaTGW-7B and miR1122a-TaNF-YB1. We propose a network of circRNA and miRNA-mediated gene regulation in the development of wheat grains.

Genome-wide expression analysis of novel heat-responsive microRNAs and their targets in contrasting wheat genotypes at reproductive stage under terminal heat stress

Introduction

Heat stress at terminal stage of wheat is critical and leads to huge yield losses worldwide. microRNAs (miRNAs) play significant regulatory roles in gene expression associated with abiotic and biotic stress at the post-transcriptional level.

Methods

In the present study, we carried out a comparative analysis of miRNAs and their targets in flag leaves as well as developing seeds of heat tolerant (RAJ3765) and heat susceptible (HUW510) wheat genotypes under heat stress and normal conditions using small RNA and degradome sequencing.

Results and discussion

A total of 84 conserved miRNAs belonging to 35 miRNA families and 93 novel miRNAs were identified in the 8 libraries. Tae-miR9672a-3p, tae-miR9774, tae-miR9669-5p, and tae-miR5048-5p showed the highest expression under heat stress. Tae-miR9775, tae-miR9662b-3p, tae-miR1120a, tae-miR5084, tae-miR1122a, tae-miR5085, tae-miR1118, tae-miR1130a, tae-miR9678-3p, tae-miR7757-5p, tae-miR9668-5p, tae-miR5050, tae-miR9652-5p, and tae-miR9679-5p were expressed only in the tolerant genotype, indicating their role in heat tolerance. Comparison between heat-treated and control groups revealed that 146 known and 57 novel miRNAs were differentially expressed in the various tissues. Eight degradome libraries sequence identified 457 targets of the differentially expressed miRNAs. Functional analysis of the targets indicated their involvement in photosynthesis, spliceosome, biosynthesis of nucleotide sugars and protein processing in the endoplasmic reticulum, arginine and proline metabolism and endocytosis.

Conclusion

This study increases the number of identified and novel miRNAs along with their roles involved in heat stress response in contrasting genotypes at two developing stages of wheat.

The roles of non-coding RNAs in male reproductive development and abiotic stress responses during this unique process in flowering plants

Successful male reproductive development is the guarantee for sexual reproduction of flowering plants. Male reproductive development is a complicated and multi-stage process that integrates physiological processes and adaptation and tolerance to a myriad of environmental stresses. This well-coordinated process is governed by genetic and epigenetic machineries. Non-coding RNAs (ncRNAs) play pleiotropic roles in the plant growth and development. The identification, characterization and functional analysis of ncRNAs and their target genes have opened a new avenue for comprehensively revealing the regulatory network of male reproductive development and its response to environmental stresses in plants. This review briefly addresses the types, origin, biogenesis and mechanisms of ncRNAs in plants, highlights important updates on the roles of ncRNAs in regulating male reproductive development and emphasizes the contribution of ncRNAs, especially miRNAs and lncRNAs, in responses to abiotic stresses during this unique process in flowering plants.

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

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

The plant noncoding transcriptome: a versatile environmental sensor

Plant noncoding RNA transcripts have gained increasing attention in recent years due to growing evidence that they can regulate developmental plasticity. In this review article, we comprehensively analyze the relationship between noncoding RNA transcripts in plants and their response to environmental cues. We first provide an overview of the various noncoding transcript types, including long and small RNAs, and how the environment modulates their performance. We then highlight the importance of noncoding RNA secondary structure for their molecular and biological functions. Finally, we discuss recent studies that have unveiled the functional significance of specific long noncoding transcripts and their molecular partners within ribonucleoprotein complexes during development and in response to biotic and abiotic stress. Overall, this review sheds light on the fascinating and complex relationship between dynamic noncoding transcription and plant environmental responses, and highlights the need for further research to uncover the underlying molecular mechanisms and exploit the potential of noncoding transcripts for crop resilience in the context of global warming.

Functional role of microRNA in the regulation of biotic and abiotic stress in agronomic plants

The increasing demand for food is the result of an increasing population. It is crucial to enhance crop yield for sustainable production. Recently, microRNAs (miRNAs) have gained importance because of their involvement in crop productivity by regulating gene transcription in numerous biological processes, such as growth, development and abiotic and biotic stresses. miRNAs are small, non-coding RNA involved in numerous other biological functions in a plant that range from genomic integrity, metabolism, growth, and development to environmental stress response, which collectively influence the agronomic traits of the crop species. Additionally, miRNA families associated with various agronomic properties are conserved across diverse plant species. The miRNA adaptive responses enhance the plants to survive environmental stresses, such as drought, salinity, cold, and heat conditions, as well as biotic stresses, such as pathogens and insect pests. Thus, understanding the detailed mechanism of the potential response of miRNAs during stress response is necessary to promote the agronomic traits of crops. In this review, we updated the details of the functional aspects of miRNAs as potential regulators of various stress-related responses in agronomic plants.

Non-coding RNAs-mediated environmental surveillance determines male fertility in plants

Plants are continuously exposed to environmental stresses leading to significant yield losses. With the changing climatic conditions, the intensity and duration of these stresses are expected to increase, posing a severe threat to crop productivity worldwide. Male gametogenesis is one of the most sensitive developmental stages. Exposure to environmental stresses during this stage leads to male sterility and yield loss. Elucidating the underlying molecular mechanism of environment-affected male sterility is essential to address this challenge. High-throughput RNA sequencing studies, loss-of-function phenotypes of sRNA biogenesis genes and functional genomics studies with non-coding RNAs have started to unveil the roles of small RNAs, long non-coding RNAs and the complex regulatory interactions between them in regulating male fertility under different growth regimes. Here, we discuss the current understanding of the non-coding RNA-mediated environmental stress surveillance and regulation of male fertility in plants. The candidate ncRNAs emerging from these studies can be leveraged to generate environment-sensitive male sterile lines for hybrid breeding or mitigate the impact of climate change on male fertility, as the situation demands.

Microspore-expressed SCULP1 is required for p-coumaroylation of sporopollenin, exine integrity, and pollen development in wheat

Sporopollenin is one of the most structurally sophisticated and chemically recalcitrant biopolymers. In higher plants, sporopollenin is the dominant component of exine, the outer wall of pollen grains, and contains covalently linked phenolics that protect the male gametes from harsh environments. Although much has been learned about the biosynthesis of sporopollenin precursors in the tapetum, the nutritive cell layer surrounding developing microspores, little is known about how the biopolymer is assembled on the microspore surface. We identified SCULP1 (SKS clade universal in pollen) as a seed plant conserved clade of the multicopper oxidase family. We showed that SCULP1 in common wheat (Triticum aestivum) is specifically expressed in the microspore when sporopollenin assembly takes place, localized to the developing exine, and binds p-coumaric acid in vitro. Through genetic, biochemical, and 3D reconstruction analyses, we demonstrated that SCULP1 is required for p-coumaroylation of sporopollenin, exine integrity, and pollen viability. Moreover, we found that SCULP1 accumulation is compromised in thermosensitive genic male sterile wheat lines and its expression partially restored exine integrity and male fertility. These findings identified a key microspore protein in autonomous sporopollenin polymer assembly, thereby laying the foundation for elucidating and engineering sporopollenin biosynthesis.

Research progress on plant noncoding RNAs in response to low-temperature stress

Low temperature (LT) is an important factor limiting plant growth and distribution. Plants have evolved sophisticated adaptive mechanisms to cope with hypothermia. RNA silencing is the orchestrator of these cellular responses. RNA silencing, which modifies gene expression through noncoding RNAs (ncRNAs), is a strategy used by plants to combat environmental stress. ncRNAs, which have very little protein-coding capacity, work by binding reverse complementary endogenous transcripts. In plants, ncRNAs include small non-coding RNAs (sncRNAs), medium-sized non-coding RNAs (mncRNAs), and long non-coding RNAs (lncRNAs). Apart from describing the biogenesis of different ncRNAs (miRNAs, siRNAs, and lncRNAs), we thoroughly discuss the functions of these ncRNAs during cold acclimation. Two major classes of sncRNAs, microRNAs and siRNAs, play essential regulatory roles in cold response processes through the posttranscriptional gene silencing (PTGS) pathway or transcriptional gene silencing (TGS) pathway. Microarray or transcriptome sequencing analysis can reveal a large number of cold-responsive miRNAs in plants. In this review, the cold-response patterns of miRNAs verified by Northern blotting or quantitative PCR in Arabidopsis thaliana, rice, and many other important crops are discussed. The detailed molecular mechanisms of several miRNAs in Arabidopsis (miR397, miR408, miR402, and miR394) and rice (Osa-miR156, Osa-miR319, and Osa-miR528) that regulate plant cold resistance are elucidated. In addition, the regulatory mechanism of the lncRNA SVALKA in the cold signaling pathway is explained in detail. Finally, we present the challenges for understanding the roles of small ncRNAs in cold signal transduction.

Comparative transcriptome analysis provides insight into the important pathways and key genes related to the pollen abortion in the thermo-sensitive genic male sterile line 373S in Brassica napus L

The thermo-sensitive genic male sterility (TGMS) system plays a key role in the production of two-line hybrids in rapeseed (Brassica napus). To uncover key cellular events and genetic regulation associated with TGMS, a combined study using cytological methods and RNA-sequencing analysis was conducted for the rapeseed TGMS line 373S. Cytological studies showed that microspore cytoplasm of 373S plants was condensed, the microspore nucleus was degraded at an early stage, the exine was irregular, and the tapetum developed abnormally, eventually leading to male sterility. RNA-sequencing analysis identified 430 differentially expressed genes (298 upregulated and 132 downregulated) between the fertile and sterile samples. Gene ontology analysis demonstrated that the most highly represented biological processes included sporopollenin biosynthetic process, pollen exine formation, and extracellular matrix assembly. Kyoto encyclopedia of genes and genomes analysis indicated that the enriched pathways included amino acid metabolism, carbohydrate metabolism, and lipid metabolism. Moreover, 26 transcript factors were identified, which may be associated with abnormal tapetum degeneration and exine formation. Subsequently, 19 key genes were selected, which are considered to regulate pollen development and even participate in pollen exine formation. Our results will provide important insight into the molecular mechanisms underlying TGMS in rapeseed.

Comparative Transcriptome Analysis Reveals Hormone Signal Transduction and Sucrose Metabolism Related Genes Involved in the Regulation of Anther Dehiscence in Photo-Thermo-Sensitive Genic Male Sterile Wheat

Anther dehiscence is an important process to release pollen and then is a critical event in pollination. In the wheat photo-thermo-sensitive genic male sterility (PTGMS) line, pollen cannot release from anther since the anther cannot dehisce during anther dehiscence stage in a sterile condition. In this study, we carried out RNA-sequencing to analyze the transcriptome of one wheat PTGMS line BS366 during anther dehiscence under fertile and sterile conditions to explore the mechanism. We identified 6306 differentially expressed genes (DEGs). Weighted gene co-expression network analysis (WGCNA) and KEGG analysis showed that DEGs were mainly related to “hormone signal transduction pathway” and “starch and sucrose metabolism”. We identified 35 and 23 DEGs related hormone signal transduction and sucrose metabolism, respectively. Compared with conventional wheat Jing411, there were some changes in the contents of hormones, including JA, IAA, BR, ABA and GA3, and sucrose, during three anther dehiscence stages in the sterile condition in BS366. We performed qRT-PCR to verify the expression levels of some critical DEGs of the hormone signaling pathway and the starch and sucrose metabolism pathway. The results showed disparate expression patterns of the critical DEGs of the hormone signaling pathway and the starch and sucrose metabolism pathway in different conditions, suggesting these genes may be involved in the regulation of the anther dehiscence in BS366. Finally, we conducted a hypothesis model to reveal the regulation pathway of hormones and sucrose on anther dehiscence. The information provided new clues to the molecular mechanisms of anther dehiscence in wheat and improved wheat hybrid breeding.

Integrated Analysis of Microarray, Small RNA, and Degradome Datasets Uncovers the Role of MicroRNAs in Temperature-Sensitive Genic Male Sterility in Wheat

Temperature-sensitive genic male sterile (TGMS) line Beijing Sterility 366 (BS366) has been utilized in hybrid breeding for a long time, but the molecular mechanism underlying male sterility remains unclear. Expression arrays, small RNA, and degradome sequencing were used in this study to explore the potential role of miRNA in the cold-induced male sterility of BS366. Microspore observation showed defective cell plates in dyads and tetrads and shrunken microspores at the vacuolated stage. Differential regulation of Golgi vesicle transport, phragmoplast formation, sporopollenin biosynthesis, pollen exine formation, and lipid metabolism were observed between cold and control conditions. Pollen development was significantly represented in the 352 antagonistic miRNA-target pairs in the integrated analysis of miRNA and mRNA profiles. The specific cleavage of ARF17 and TIR1 by miR160 and miR393 were found in the cold-treated BS366 degradome, respectively. Thus, the cold-mediated miRNAs impaired cell plate formation through repression of Golgi vesicle transport and phragmoplast formation. The repressed expression of ARF17 and TIR1 impaired pollen exine formation. The results of this study will contribute to our understanding of the roles of miRNAs in male sterility in wheat.

Comprehensive molecular evaluation of the histone methyltransferase gene family and their important roles in two-line hybrid wheat

BACKGROUND

Histone methylation usually plays important roles in plant development through post-translational regulation and may provide a new visual field for heterosis. The histone methyltransferase gene family has been identified in various plants, but its members and functions in hybrid wheat related in heterosis is poorly studied.

RESULTS

In this study, 175 histone methyltransferase (HMT) genes were identified in wheat, including 152 histone lysine methyltransferase (HKMT) genes and 23 protein arginine methyltransferase (PRMT) genes. Gene structure analysis, physicochemical properties and subcellular localization predictions of the proteins, exhibited the adequate complexity of this gene family. As an allohexaploid species, the number of the genes (seven HKMTs orthologous groups and four PRMTs orthologous groups) in wheat were about three times than those in diploids and showed certain degrees of conservation, while only a small number of subfamilies such as ASH-like and Su-(var) subfamilies have expanded their members. Transcriptome analysis showed that HMT genes were mainly expressed in the reproductive organs. Expression analysis showed that some TaHMT genes with different trends in various hybrid combinations may be regulated by lncRNAs with similar expression trends. Pearson correlation analysis of the expression of TaHMT genes and two yield traits indicated that four DEGs may participate in the yield heterosis of two-line hybrid wheat. ChIP-qPCR results showed that the histone modifications (H3K4me3, H3K36me3 and H3K9ac) enriched in promoter regions of three TaCCA1 genes which are homologous to Arabidopsis heterosis-related CCA1/LHY genes. The higher expression levels of TaCCA1 in F₁ than its parents are positive with these histone modifications. These results showed that histone modifications may play important roles in wheat heterosis.

CONCLUSIONS

Our study identified characteristics of the histone methyltransferase gene family and enhances the understanding of the evolution and function of these members in allohexaploid wheat. The causes of heterosis of two-line hybrid wheat were partially explained from the perspective of histone modifications.

Genome wide Identification and Characterization of Wheat GH9 Genes Reveals Their Roles in Pollen Development and Anther Dehiscence

Glycoside hydrolase family 9 (GH9) is a key member of the hydrolase family in the process of cellulose synthesis and hydrolysis, playing important roles in plant growth and development. In this study, we investigated the phenotypic characteristics and gene expression involved in pollen fertility conversion and anther dehiscence from a genomewide level. In total, 74 wheat GH9 genes (TaGH9s) were identified, which were classified into Class A, Class B and Class C and unevenly distributed on chromosomes. We also investigated the gene duplication and reveled that fragments and tandem repeats contributed to the amplification of TaGH9s. TaGH9s had abundant hormone-responsive elements and light-responsive elements, involving JA–ABA crosstalk to regulate anther development. Ten TaGH9s, which highly expressed stamen tissue, were selected to further validate their function in pollen fertility conversion and anther dehiscence. Based on the cell phenotype and the results of the scanning electron microscope at the anther dehiscence period, we found that seven TaGH9s may target miRNAs, including some known miRNAs (miR164 and miR398), regulate the level of cellulose by light and phytohormone and play important roles in pollen fertility and anther dehiscence. Finally, we proposed a hypothesis model to reveal the regulation pathway of TaGH9 on fertility conversion and anther dehiscence. Our study provides valuable insights into the GH9 family in explaining the male sterility mechanism of the wheat photo-thermo-sensitive genetic male sterile (PTGMS) line and generates useful male sterile resources for improving wheat hybrid breeding.

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

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

The Characters of Non-Coding RNAs and Their Biological Roles in Plant Development and Abiotic Stress Response

Plant growth and development are greatly affected by the environment. Many genes have been identified to be involved in regulating plant development and adaption of abiotic stress. Apart from protein-coding genes, more and more evidence indicates that non-coding RNAs (ncRNAs), including small RNAs and long ncRNAs (lncRNAs), can target plant developmental and stress-responsive mRNAs, regulatory genes, DNA regulatory regions, and proteins to regulate the transcription of various genes at the transcriptional, posttranscriptional, and epigenetic level. Currently, the molecular regulatory mechanisms of sRNAs and lncRNAs controlling plant development and abiotic response are being deeply explored. In this review, we summarize the recent research progress of small RNAs and lncRNAs in plants, focusing on the signal factors, expression characters, targets functions, and interplay network of ncRNAs and their targets in plant development and abiotic stress responses. The complex molecular regulatory pathways among small RNAs, lncRNAs, and targets in plants are also discussed. Understanding molecular mechanisms and functional implications of ncRNAs in various abiotic stress responses and development will benefit us in regard to the use of ncRNAs as potential character-determining factors in molecular plant breeding.

Global survey of alternative splicing and gene modules associated with fertility regulation in a thermosensitive genic male sterile wheat

Thermosensitive genic male sterile (TGMS) wheat lines are the core of two-line hybrid systems. Understanding the mechanism that regulates male sterility in TGMS wheat lines is helpful for promoting wheat breeding. Several studies have obtained information regarding the mechanisms associated with male sterility at the transcriptional level, but it is not clear how the post-transcriptional process of alternative splicing might contribute to controlling male sterility. In this study, we performed genome-wide analyses of alternative splicing during the meiosis stage in TGMS line BS366 using PacBio and RNA-Seq hybrid sequencing. Cytological observations indicated that cytoskeleton assembly in pollen cells, calcium deposition in pollen and tapetal cells, and vesicle transport in tapetal cells were deficient in BS366. According to our cytological findings, 49 differentially spliced genes were isolated. Moreover, 25 long non-coding RNA targets and three bHLH transcription factors were identified. Weighted gene co-expression network analysis detected four candidate differentially spliced genes that had strong co-relation with the seed setting percentage, which is the direct representation of male sterility in BS366. In this study, we obtained comprehensive data regarding the alternative splicing-mediated regulation of male sterility in TGMS wheat. The candidates identified may provide the molecular basis for an improved understanding of male sterility.

Temporal expression study of miRNAs in the crown tissues of winter wheat grown under natural growth conditions

BACKGROUND

Winter wheat requires prolonged exposure to low temperature to initiate flowering (vernalization). Shoot apical meristem of the crown is the site of cold perception, which produces leaf primordia during vegetative growth before developing into floral primordia at the initiation of the reproductive phase. Although many essential genes for winter wheat cold acclimation and floral initiation have been revealed, the importance of microRNA (miRNA) meditated post-transcriptional regulation in crowns is not well understood. To understand the potential roles of miRNAs in crown tissues, we performed a temporal expression study of miRNAs in crown tissues at the three-leaf stage, winter dormancy stage, spring green-up stage, and jointing stage of winter wheat grown under natural growth conditions.

RESULTS

In total, 348 miRNAs belonging to 298 miRNA families, were identified in wheat crown tissues. Among them, 92 differentially expressed miRNAs (DEMs) were found to be significantly regulated from the three-leaf stage to the jointing stage. Most of these DEMs were highly expressed at the three-leaf stage and winter dormancy stage, and then declined in later stages. Six DEMs, including miR156a-5p were markedly induced during the winter dormancy stage. Eleven DEMs, including miR159a.1, miR390a-5p, miR393-5p, miR160a-5p, and miR1436, were highly expressed at the green-up stage. Twelve DEMs, such as miR172a-5p, miR394a, miR319b-3p, and miR9676-5p were highly induced at the jointing stage. Moreover, 14 novel target genes of nine wheat or Pooideae-specific miRNAs were verified using RLM-5' RACE assay. Notably, six mTERFs and two Rf1 genes, which are associated with mitochondrial gene expression, were confirmed as targets of three wheat-specific miRNAs.

CONCLUSIONS

The present study not only confirmed the known miRNAs associated with phase transition and floral development, but also identified a number of wheat or Pooideae-specific miRNAs critical for winter wheat cold acclimation and floral development. Most importantly, this study provided experimental evidence that miRNA could regulate mitochondrial gene expression by targeting mTERF and Rf1 genes. Our study provides valuable information for further exploration of the mechanism of miRNA mediated post-transcriptional regulation during winter wheat vernalization and inflorescent initiation.

The Dynamism of Transposon Methylation for Plant Development and Stress Adaptation

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

Comprehensive analysis of formin gene family highlights candidate genes related to pollen cytoskeleton and male fertility in wheat (Triticum aestivum L.)

Background

Formin, a highly conserved multi-domain protein, interacts with microfilaments and microtubules. Although specifically expressed formin genes in anthers are potentially significant in research on male sterility and hybrid wheat breeding, similar reports in wheat, especially in thermo-sensitive genic male sterile (TGMS) wheat, remain elusive.

Results

Herein, we systematically characterized the formin genes in TGMS wheat line BS366 named TaFormins (TaFHs) and predicted their functions in inducing stress response. In total, 25 TaFH genes were uncovered, majorly localized in 2A, 2B, and 2D chromosomes. According to the neighbor-joining (NJ) method, all TaFH proteins from wheat and other plants clustered in 6 sub-groups (A-F). The modeled 3D structures of TaFH1-A/B, TaFH2-A/B, TaFH3-A/B and TaFH3-B/D were validated. And different numbers of stress and hormone-responsive regulatory elements in their 1500 base pair promoter regions were contained in the TaFH genes copies. TaFHs had specific temporal and spatial expression characteristics, whereby TaFH1, TaFH4, and TaFH5 were expressed highly in the stamen of BS366. Besides, the accumulation of TaFHs was remarkably lower in a low-temperature sterile condition (Nanyang) than fertile condition (Beijing), particularly at the early stamen development stage. The pollen cytoskeleton of BS366 was abnormal in the three stages under sterile and fertile environments. Furthermore, under different stress levels, TaFHs expression could be induced by drought, salt, abscisic acid (ABA), salicylic acid (SA), methyl jasmonate (MeJA), indole-3-acetic acid (IAA), polyethylene glycol (PEG), and low temperature. Some miRNAs, including miR167, miR1120, and miR172, interacts with TaFH genes; thus, we constructed an interaction network between microRNAs, TaFHs, phytohormone responses, and distribution of cytoskeleton to reveal the regulatory association between upstream genes of TaFH family members and sterile.

Conclusions

Collectively, this comprehensive analysis provides novel insights into TaFHs and miRNA resources for wheat breeding. These findings are, therefore, valuable in understanding the mechanism of TGMS fertility conversion in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07878-7.

Genome-wide profiling of long noncoding RNAs involved in wheat spike development

BACKGROUND

Long noncoding RNAs (lncRNAs) have been shown to play important roles in the regulation of plant growth and development. Recent transcriptomic analyses have revealed the gene expression profiling in wheat spike development, however, the possible regulatory roles of lncRNAs in wheat spike morphogenesis remain largely unclear.

RESULTS

Here, we analyzed the genome-wide profiling of lncRNAs during wheat spike development at six stages, and identified a total of 8,889 expressed lncRNAs, among which 2,753 were differentially expressed lncRNAs (DE lncRNAs) at various developmental stages. Three hundred fifteen differentially expressed cis- and trans-regulatory lncRNA-mRNA pairs comprised of 205 lncRNAs and 279 genes were predicted, which were found to be mainly involved in the stress responses, transcriptional and enzymatic regulations. Moreover, the 145 DE lncRNAs were predicted as putative precursors or target mimics of miRNAs. Finally, we identified the important lncRNAs that participate in spike development by potentially targeting stress response genes, TF genes or miRNAs.

CONCLUSIONS

This study outlines an overall view of lncRNAs and their possible regulatory networks during wheat spike development, which also provides an alternative resource for genetic manipulation of wheat spike architecture and thus yield.

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.

TarDB: an online database for plant miRNA targets and miRNA-triggered phased siRNAs

Background

In plants, microRNAs (miRNAs) are pivotal regulators of plant development and stress responses. Different computational tools and web servers have been developed for plant miRNA target prediction; however, in silico prediction normally contains false positive results. In addition, many plant miRNA target prediction servers lack information for miRNA-triggered phased small interfering RNAs (phasiRNAs). Creating a comprehensive and relatively high-confidence plant miRNA target database is much needed.

Results

Here, we report TarDB, an online database that collects three categories of relatively high-confidence plant miRNA targets: (i) cross-species conserved miRNA targets; (ii) degradome/PARE (Parallel Analysis of RNA Ends) sequencing supported miRNA targets; (iii) miRNA-triggered phasiRNA loci. TarDB provides a user-friendly interface that enables users to easily search, browse and retrieve miRNA targets and miRNA initiated phasiRNAs in a broad variety of plants. TarDB has a comprehensive collection of reliable plant miRNA targets containing previously unreported miRNA targets and miRNA-triggered phasiRNAs even in the well-studied model species. Most of these novel miRNA targets are relevant to lineage-specific or species-specific miRNAs. TarDB data is freely available at http://www.biosequencing.cn/TarDB.

Conclusions

In summary, TarDB serves as a useful web resource for exploring relatively high-confidence miRNA targets and miRNA-triggered phasiRNAs in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07680-5.

Aux/IAA14 Regulates microRNA-Mediated Cold Stress Response in Arabidopsis Roots

The phytohormone auxin and microRNA-mediated regulation of gene expressions are key regulators of plant growth and development at both optimal and under low-temperature stress conditions. However, the mechanistic link between microRNA and auxin in regulating plant cold stress response remains elusive. To better understand the role of microRNA (miR) in the crosstalk between auxin and cold stress responses, we took advantage of the mutants of Arabidopsis thaliana with altered response to auxin transport and signal. Screening of the mutants for root growth recovery after cold stress at 4 °C revealed that the auxin signaling mutant, solitary root 1 (slr1; mutation in Aux/IAA14), shows a hypersensitive response to cold stress. Genome-wide expression analysis of miRs in the wild-type and slr1 mutant roots using next-generation sequencing revealed 180 known and 71 novel cold-responsive microRNAs. Cold stress also increased the abundance of 26-31 nt small RNA population in slr1 compared with wild type. Comparative analysis of microRNA expression shows significant differential expression of 13 known and 7 novel miRs in slr1 at 4 °C compared with wild type. Target gene expression analysis of the members from one potential candidate miR, miR169, revealed the possible involvement of miR169/NF-YA module in the Aux/IAA14-mediated cold stress response. Taken together, these results indicate that SLR/IAA14, a transcriptional repressor of auxin signaling, plays a crucial role in integrating miRs in auxin and cold responses.

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

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

Uncloaking lncRNA-meditated gene expression as a potential regulator of CMS in cotton (Gossypium hirsutum L.)

Cytoplasmic male sterility is a well-proven mechanism for cotton hybrid production. Long non-coding RNAs belong to a class of transcriptional regulators that function in multiple biological processes. The cDNA libraries from the flower buds of the cotton CGMS, it's restorer (Rf) and maintainer lines were sequenced using high throughput NGS technique. A total of 1531 lncRNAs showed significant differential expression patterns between these three lines. Functional analysis of the co-expression network of lncRNA-mRNA using gene ontology vouchsafes that, lncRNAs play a crucial role in cytoplasmic male sterility and fertility restoration through pollen development, INO80 complex, development of anther wall tapetum, chromatin remodeling, and histone modification. Additionally, 94 lncRNAs were identified as putative precursors of 49 miRNAs. qRT-PCR affirms the concordance of expression pattern to RNA-seq data. These findings divulge the lncRNA driven miRNA-mediated regulation of gene expression profiling superintended for a better understanding of the CMS mechanisms of cotton.

Molecular genetic analyses of abiotic stress responses during plant reproductive development

Plant responses to abiotic stresses during vegetative growth have been extensively studied for many years. Daily environmental fluctuations can have dramatic effects on plant vegetative growth at multiple levels, resulting in molecular, cellular, physiological, and morphological changes. Plants are even more sensitive to environmental changes during reproductive stages. However, much less is known about how plants respond to abiotic stresses during reproduction. Fortunately, recent advances in this field have begun to provide clues about these important processes, which promise further understanding and a potential contribution to maximize crop yield under adverse environments. Here we summarize information from several plants, focusing on the possible mechanisms that plants use to cope with different types of abiotic stresses during reproductive development, and present a tentative molecular portrait of plant acclimation during reproductive stages. Additionally, we discuss strategies that plants use to balance between survival and productivity, with some comparison among different plants that have adapted to distinct environments.

PAREameters: a tool for computational inference of plant miRNA-mRNA targeting rules using small RNA and degradome sequencing data

MicroRNAs (miRNAs) are short, non-coding RNAs that modulate the translation-rate of messenger RNAs (mRNAs) by directing the RNA-induced silencing complex to sequence-specific targets. In plants, this typically results in cleavage and subsequent degradation of the mRNA. Degradome sequencing is a high-throughput technique developed to capture cleaved mRNA fragments and thus can be used to support miRNA target prediction. The current criteria used for miRNA target prediction were inferred on a limited number of experimentally validated A. thaliana interactions and were adapted to fit these specific interactions; thus, these fixed criteria may not be optimal across all datasets (organisms, tissues or treatments). We present a new tool, PAREameters, for inferring targeting criteria from small RNA and degradome sequencing datasets. We evaluate its performance using a more extensive set of experimentally validated interactions in multiple A. thaliana datasets. We also perform comprehensive analyses to highlight and quantify the differences between subsets of miRNA-mRNA interactions in model and non-model organisms. Our results show increased sensitivity in A. thaliana when using the PAREameters inferred criteria and that using data-driven criteria enables the identification of additional interactions that further our understanding of the RNA silencing pathway in both model and non-model organisms.

MicroRNA expression patterns unveil differential expression of conserved miRNAs and target genes against abiotic stress in safflower

Environmental stresses influence the growth and development of plants by influencing patterns of gene expression. Different regulators control gene expression, including transcription factors (TFs) and microRNAs. MicroRNAs (miRNAs: ~21 nucleotides long) are encoded by miRNA genes transcribed by RNA polymerase II (RNP-II) and play key roles in plant development and physiology. There is little knowledge currently available on miRNAs and their function in response to environmental stresses in safflower. To obtain more information on safflower miRNAs, we initially used a comparative genomics approach and succeeded in identifying 126 miRNAs belonging to 29 conserved families, along with their target genes. In this study, we investigated the expression profiles of seven conserved miRNAs related to drought, salinity, heat, and Cd stress in the leaf and root organs using qRT-PCR, for the first time. Gene Ontology (GO) analysis found that target genes of miRNAs are often TFs such as AP2/ERF and HD-ZIP as well as NAC domain-containing proteins. Expression analyses confirmed that miRNAs can play a vital role in keeping safflower stress-tolerant. Differential expression of miR156, miR162, miR164, miR166, miR172, miR398, and miR408 regulate the expression of their respective target genes. These genes activate several pathways leading to physiological and biochemical responses to abiotic stresses. Some conserved miRNAs were regulated by abiotic stresses. Our finding provides valuable information to understand miRNAs in relation to different abiotic stresses in safflower.

Identification and characterization of miRNAs involved in cold acclimation of zebrafish ZF4 cells

MicroRNAs (miRNAs) play vital roles in various biological processes under multiple stress conditions by leading to mRNA cleavage or translational repression. However, the detailed roles of miRNAs in cold acclimation in fish are still unclear. In the present study, high-throughput sequencing was performed to identify miRNAs from 6 small RNA libraries from the zebrafish embryonic fibroblast ZF4 cells under control (28°C, 30 days) and cold-acclimation (18°C, 30 days) conditions. A total of 414 miRNAs, 349 known and 65 novel, were identified. Among those miRNAs, 24 (19 known and 5 novel) were up-regulated, and 23 (9 known and 14 novel) were down-regulated in cold acclimated cells. The Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses indicated that the target genes of known differentially expressed miRNAs (DE-miRNA) are involved in cold acclimation by regulation of phosphorylation, cell junction, intracellular signal transduction, ECM-receptor interaction and so on. Moreover, both miR-100-3p inhibitor and miR-16b mimics could protect ZF4 cells under cold stress, indicating the involvement of miRNA in cold acclimation. Further study showed that miR-100-3p and miR-16b could regulate inversely the expression of their target gene (atad5a, cyp2ae1, lamp1, rilp, atxn7, tnika, btbd9), and that overexpression of miR-100-3p disturbed the early embryonic development of zebrafish. In summary, the present data show that miRNAs are closely involved in cold acclimation in zebrafish ZF4 cells and provide information for further understanding of the roles of miRNAs in cold acclimation in fish.

Comprehensive Atlas of Wheat (Triticum aestivum L.) AUXIN RESPONSE FACTOR Expression During Male Reproductive Development and Abiotic Stress

AUXIN RESPONSE FACTOR (ARF) proteins regulate a wide range of signaling pathways, from general plant growth to abiotic stress responses. Here, we performed a genome-wide survey in wheat (Triticum aestivum) and identified 69 TaARF members that formed 24 homoeologous groups. Phylogenetic analysis clustered TaARF genes into three clades, similar to ARF genes in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa). Structural characterization suggested that ARF gene structure and domain composition are well conserved between plant species. Expression profiling revealed diverse patterns of TaARF transcript levels across a range of developmental stages, tissues, and abiotic stresses. A number of TaARF genes shared similar expression patterns and were preferentially expressed in anthers. Moreover, our systematic analysis identified three anther-specific TaARF genes (TaARF8, TaARF9, and TaARF21) whose expression was significantly altered by low temperature in thermosensitive genic male-sterile (TGMS) wheat; these TaARF genes are candidates to participate in the cold-induced male sterility pathway, and offer potential applications in TGMS wheat breeding and hybrid seed production. Moreover, we identified putative functions for a set of TaARFs involved in responses to abscisic acid and abiotic stress. Overall, this study characterized the wheat ARF gene family and generated several hypotheses for future investigation of ARF function during anther development and abiotic stress.

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

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

The Major Factors Causing the Microspore Abortion of Genic Male Sterile Mutant NWMS1 in Wheat (Triticum aestivum L.)

Male sterility is a valuable trait for genetic research and production application of wheat (Triticum aestivum L.). NWMS1, a novel typical genic male sterility mutant, was obtained from Shengnong 1, mutagenized with ethyl methane sulfonate (EMS). Microstructure and ultrastructure observations of the anthers and microspores indicated that the pollen abortion of NWMS1 started at the early uninucleate microspore stage. Pollen grain collapse, plasmolysis, and absent starch grains were the three typical characteristics of the abnormal microspores. The anther transcriptomes of NWMS1 and its wild type Shengnong 1 were compared at the early anther development stage, pollen mother cell meiotic stage, and binucleate microspore stage. Several biological pathways clearly involved in abnormal anther development were identified, including protein processing in endoplasmic reticulum, starch and sucrose metabolism, lipid metabolism, and plant hormone signal transduction. There were 20 key genes involved in the abnormal anther development, screened out by weighted gene co-expression network analysis (WGCNA), including SKP1B, BIP5, KCS11, ADH3, BGLU6, and TIFY10B. The results indicated that the defect in starch and sucrose metabolism was the most important factor causing male sterility in NWMS1. Based on the experimental data, a primary molecular regulation model of abnormal anther and pollen developments in mutant NWMS1 was established. These results laid a solid foundation for further research on the molecular mechanism of wheat male sterility.

Blocked synthesis of sporopollenin and jasmonic acid leads to pollen wall defects and anther indehiscence in genic male sterile wheat line 4110S at high temperatures

Environment-sensitive genic male sterility is a valid tool for hybrid production and hybrid breeding, but there are no previous reports of the molecular mechanism of fertility conversion. In this study, RNA-seq, phenotypic and cytological observations, and physiological indexes were applied to analyze thermo-sensitive genic male sterility line 4110S under different temperature conditions to explore the fertility transformation mechanism. In total, 3420 differentially expressed genes (DEGs) were identified comprising 2331 upregulated genes and 1089 downregulated genes. The DEGs were apparently distributed among 54 Gene Ontology functional groups. The phenylpropanoid, long-chain fatty acid, and jasmonic acid (JA) biosynthesis pathways were related to male sterility, where their downregulation blocked the synthesis of sporopollenin and JA. Phenotypic and cytological analyses showed that pollen wall defects and anther indehiscence at high temperatures induced sterility. Moreover, enzyme-linked immunosorbent assay results indicated that the abundance of JA was lower in 4110S under restrictive conditions (high temperature) than permissive conditions (low temperature). A possible regulated network of pathways associated with male sterility was suggested. These results provided insights into the molecular mechanism of fertility conversion in the thermosensitive male sterility system.

Comparative profiling of roots small RNA expression and corresponding gene ontology and pathway analyses for low- and high-cadmium-accumulating genotypes of wheat in response to cadmium stress

MicroRNAs (miRNAs) participate in multiple biological processes in plant. Cd accumulation ability differs among varieties in wheat, but little is known about miRNAs and their function in Cd accumulation of wheat under Cd stress. Therefore, the present study detected small RNAs responsible for differential Cd accumulation between two contrasting wheat genotypes (low-Cd accumulation one L17 and high-Cd accumulation one H17) to identify novel targets to further study Cd stress in wheat. Under normal conditions, 139 miRNAs were differentially expressed between L17 and H17, while this value reached 142 after Cd exposure. For Cd-induced DEMs, total 25 miRNAs were differentially expressed in L17 after Cd treatment, while, 70 Cd-induced DEMs were found in H17. Moreover, GO analysis revealed that target genes of DEMs related to lipid biosynthetic process and chlorophyll binding are uniquely enriched in L17, target genes of DEMs related to ribosome biogenesis and sucrose alpha-glucosidase activity are uniquely enriched in H17. By pathway analysis, target genes of DEMs related to PI3K-Akt signaling pathway was specifically enriched in L17, target genes of DEMs related to carbohydrate digestion and absorption pathway was uniquely enriched in H17. In addition, miRNA-gene co-expression showed that tae-miR9774 was uniquely expressed between L17Cd and L17CK, while tae-miR398 was specially expressed between H17Cd and H17CK. Our results suggested that Cd-accumulating ability of L17 and H17 varied from the expression of induced unique miRNA, such as expression of tae-miR-9774 and tae-miR-398. Our study not only provide the foundation for further exploring the miRNAs-induced molecular mechanisms of Cd accumulation in wheat but also supply novel strategies for improving phytoremediation ability of food plants through genetic engineering.

PAREsnip2: a tool for high-throughput prediction of small RNA targets from degradome sequencing data using configurable targeting rules

Small RNAs (sRNAs) are short, non-coding RNAs that play critical roles in many important biological pathways. They suppress the translation of messenger RNAs (mRNAs) by directing the RNA-induced silencing complex to their sequence-specific mRNA target(s). In plants, this typically results in mRNA cleavage and subsequent degradation of the mRNA. The resulting mRNA fragments, or degradome, provide evidence for these interactions, and thus degradome analysis has become an important tool for sRNA target prediction. Even so, with the continuing advances in sequencing technologies, not only are larger and more complex genomes being sequenced, but also degradome and associated datasets are growing both in number and read count. As a result, existing degradome analysis tools are unable to process the volume of data being produced without imposing huge resource and time requirements. Moreover, these tools use stringent, non-configurable targeting rules, which reduces their flexibility. Here, we present a new and user configurable software tool for degradome analysis, which employs a novel search algorithm and sequence encoding technique to reduce the search space during analysis. The tool significantly reduces the time and resources required to perform degradome analysis, in some cases providing more than two orders of magnitude speed-up over current methods.

Identification of Differentially Expressed miRNAs between a Wheat K-type Cytoplasmic Male Sterility Line and Its Near-Isogenic Restorer Line

K-type cytoplasmic male sterility (KCMS) lines were ideal material for three-line hybrid wheat system due to the major role in hybrid wheat production. In this study, the morphology of developing microspore and mature pollen was compared between a KCMS line and its near-isogenic restorer line (KCMS-NIL). The most striking difference is that the microspore was unable to develop into tricellular pollen in the KCMS line. MicroRNA plays vital roles in flowering and gametophyte development. Small RNA sequencing identified a total of 274 known and 401 novel miRNAs differentially expressed between two lines or two developmental stages. Most of miRNAs with high abundance were differentially expressed at the uninucleate stage, and their expression level recovered or remained at the binucleate stage. Further degradome sequencing identified target genes which were mainly enriched in transcription regulation, phytohormone signaling and RNA degradation pathways. Combining with the transcriptome data, a correlation was found between the abnormal anther development, such as postmeiotic mitosis cessation, deformative pollen wall and the chromosome condensation of the vegetative cell, and the alterations in the related miRNA and their targets expression profiles. According to the correlation and pathway analysis, we propose a hypothetic miRNA-mediated network for the control of KCMS restoration.

Comparative profile analysis reveals differentially expressed microRNAs regulate anther and pollen development in kenaf cytoplasmic male sterility line

Cytoplasmic male sterility (CMS) is advantageous in extensive crop breeding and represents a perfect model for understanding anther and pollen development research. MicroRNAs (miRNAs) play key roles in regulating various biological processes. However, the miRNA-mediated regulatory network in kenaf CMS occurrence remains largely unknown. In the present study, a comparative deep sequencing approach was used to investigate the miRNAs and their roles in regulating anther and pollen development during CMS occurrence. We identified 283 known and 46 new candidate miRNAs in kenaf anther. A total of 67 differentially expressed miRNAs (DEMs) were discovered between CMS and its maintainer line. Among them, 40 and 27 miRNAs were up- and downregulated, respectively. These 67 DEMs were predicted to target 189 genes. Validation of DEMs and putative target genes were confirmed by using real-time quantitative PCR. In addition, a potential miRNA-mediated regulatory network, which mainly involves the auxin signaling pathway, signal transduction, glycolysis and energy metabolism, gene expression, transmembrane transport, protein modification and metabolism, and floral development, that mediates anther development during CMS occurrence was proposed. Taken together, our findings provide a better understanding of the molecular mechanism of miRNA regulation in pollen development and CMS occurrence in kenaf.

Lineage-Specific Evolved MicroRNAs Regulating NB-LRR Defense Genes in Triticeae

Disease resistance genes encoding proteins with nucleotide binding sites and Leucine-Rich Repeat (NB-LRR) domains include many members involved in the effector-triggered immunity pathway in plants. The transcript levels of these defense genes are negatively regulated by diverse microRNAs (miRNAs) in angiosperms and gymnosperms. In wheat, using small RNA expression datasets and degradome datasets, we identified five miRNA families targeting NB-LRR defense genes in monocots, some of which arose in the Triticeae species era. These miRNAs regulate different types of NB-LRR genes, most of them with coil-coiled domains, and trigger the generation of secondary small interfering RNAs (siRNA) as a phased pattern in the target site regions. In addition to acting in response to biotic stresses, they are also responsive to abiotic stresses such as heat, drought, salt, and light stress. Their copy number and expression variation in Triticeae suggest a rapid birth and death frequency. Altogether, non-conserved miRNAs as conserved transcriptional regulators in gymnosperms and angiosperms regulating the disease resistance genes displayed quick plasticity including the variations of sequences, gene copy number, functions, and expression level, which accompanied with NB-LRR genes may be tune-regulated to plants in natural environments with various biotic and abiotic stresses.

Genome-Wide Analysis of the miRNA-mRNAs Network Involved in Cold Tolerance in Populus simonii × P. nigra

Background: Cold tolerance is important for plants’ geographical distribution and survival in extreme seasonal variations of climate. However, Populus simonii × P. nigra shows wide adaptability and strong cold resistance. Transcriptional and post-transcriptional regulation of cold-responsive genes is crucial for cold tolerance in plants. To understand the roles of regulatory RNAs under cold induction in Populus simonii × P. nigra, we constructed cDNA and small RNA libraries from leaf buds treated or not with −4 °C for 8 h for analysis. Results: Through high-throughput sequencing and differential expression analysis, 61 miRNAs and 1229 DEGs were identified under cold induction condition in Populus simonii × P. nigra. The result showed that miR167a, miR1450, miR319a, miR395b, miR393a-5p, miR408-5p, and miR168a-5p were downregulated, whereas transcription level of miR172a increased under the cold treatment. Thirty-one phased-siRNA were also obtained (reads ≥ 4) and some of them proceeded from TAS3 loci. Analysis of the differentially expressed genes (DEGs) showed that transcription factor genes such as Cluster-15451.2 (putative MYB), Cluster-16493.29872 (putative bZIP), Cluster-16493.29175 (putative SBP), and Cluster-1378.1 (putative ARF) were differentially expressed in cold treated and untreated plantlets of Populus simonii × P. nigra. Integrated analysis of miRNAs and transcriptome showed miR319, miR159, miR167, miR395, miR390, and miR172 and their target genes, including MYB, SBP, bZIP, ARF, LHW, and ATL, were predicted to be involved in ARF pathway, SPL pathway, DnaJ related photosystem II, and LRR receptor kinase, and many of them are known to resist chilling injury. Particularly, a sophisticated regulatory model including miRNAs, phasiRNAs, and targets of them was set up. Conclusions: Integrated analysis of miRNAs and transcriptome uncovered the complicated regulation of the tolerance of cold in Populus simonii × P. nigra. MiRNAs, phasiRNAs, and gene-encoded transcription factors were characterized at a whole genome level and their expression patterns were proved to be complementary. This work lays a foundation for further research of the pathway of sRNAs and regulatory factors involved in cold tolerance.

Identification of miRNAs and Their Response to Cold Stress in Astragalus Membranaceus

Astragalus membranaceus is an important medicinal plant widely cultivated in East Asia. MicroRNAs (miRNAs) are endogenous regulatory molecules that play essential roles in plant growth, development, and the response to environmental stresses. Cold is one of the key environmental factors affecting the yield and quality of A. membranaceus, and miRNAs may mediate the gene regulation network under cold stress in A. membranaceus. To identify miRNAs and reveal their functions in cold stress response in A. membranaceus, small RNA sequencing was conducted followed by bioinformatics analysis, and quantitative real time PCR (qRT-PCR) analysis was performed to profile the expression of miRNAs under cold stress. A total of 168 conserved miRNAs belonging to 34 families and 14 putative non-conserved miRNAs were identified. Many miRNA targets were predicted and these targets were involved in diversified regulatory and metabolic pathways. By using qRT-PCR, 27 miRNAs were found to be responsive to cold stress, including 4 cold stress-induced and 17 cold-repressed conserved miRNAs, and 6 cold-induced non-conserved miRNAs. These cold-responsive miRNAs probably mediate the response to cold stress by regulating development, hormone signaling, defense, redox homeostasis, and secondary metabolism in A. membranaceus. These cold-corresponsive miRNAs may be used as the candidate genes in further molecular breeding for improving cold tolerance of A. membranaceus.

Development of species specific putative miRNA and its target prediction tool in wheat (Triticum aestivum L.)

MicroRNA are 20-24 nt, non-coding, single stranded molecule regulating traits and stress response. Tissue and time specific expression limits its detection, thus is major challenge in their discovery. Wheat has limited 119 miRNAs in MiRBase due to limitation of conservation based methodology where old and new miRNA genes gets excluded. This is due to origin of hexaploid wheat by three successive hybridization, older AA, BB and younger DD subgenome. Species specific miRNA prediction (SMIRP concept) based on 152 thermodynamic features of training dataset using support vector machine learning approach has improved prediction accuracy to 97.7%. This has been implemented in TamiRPred ( http://webtom.cabgrid.res.in/tamirpred ). We also report highest number of putative miRNA genes (4464) of wheat from whole genome sequence populated in database developed in PHP and MySQL. TamiRPred has predicted 2092 (>45.10%) additional miRNA which was not predicted by miRLocator. Predicted miRNAs have been validated by miRBase, small RNA libraries, secondary structure, degradome dataset, star miRNA and binding sites in wheat coding region. This tool can accelerate miRNA polymorphism discovery to be used in wheat trait improvement. Since it predicts chromosome-wise miRNA genes with their respective physical location thus can be transferred using linked SSR markers. This prediction approach can be used as model even in other polyploid crops.

Genome-wide identification and analysis of the COI gene family in wheat (Triticum aestivum L.)

Background

COI (CORONATINE INSENSITIVE), an F-box component of the Skp1-Cullin-F-box protein (SCF^COI1) ubiquitin E3 ligase, plays important roles in the regulation of plant growth and development. Recent studies have shown that COIs are involved in pollen fertility. In this study, we identified and characterized COI genes in the wheat genome and analyzed expression patterns under abiotic stress.

Results

A total of 18 COI candidate sequences for 8 members of COI gene family were isolated in wheat (Triticum aestivum L.). Phylogenetic and structural analyses showed that these COI genes could be divided into seven distinct subfamilies. The COI genes showed high expression in stamens and glumes. The qRT-PCR results revealed that wheat COIs were involved in several abiotic stress responses and anther/glume dehiscence in the photoperiod-temperature sensitive genic male sterile (PTGMS) wheat line BS366.

Conclusions

The structural characteristics and expression patterns of the COI gene family in wheat as well as the stress-responsive and differential tissue-specific expression profiles of each TaCOI gene were examined in PTGMS wheat line BS366. In addition, we examined SA- and MeJA-induced gene expression in the wheat anther and glume to investigate the role of COI in the JA signaling pathway, involved in the regulation of abnormal anther dehiscence in the PTGMS wheat line. The results of this study contribute novel and detailed information about the TaCOI gene family in wheat and could be used as a benchmark for future studies of the molecular mechanisms of PTGMS in other crops.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5116-9) contains supplementary material, which is available to authorized users.

Cellular and molecular characteristics of pollen abortion in chrysanthemum cv. Kingfisher

Microspore degeneration at the tetrad stage is associated with tapetum degeneration retardation. Some genes and proteins related to cell senescence and death are the key factors for pollen abortion. Chrysanthemum (Chrysanthemum morifolium) is a major floriculture crop in the world, but pollen contamination is an urgent problem to be solved in chrysanthemum production. C. morifolium 'Kingfisher' is a chrysanthemum cultivar that does not contain any pollen in mature anthers, thus it is a very important material for developing chrysanthemum without pollen contamination. However, the mechanism of its pollen abortion remains unclear. In this study, the cellular and molecular mechanisms of 'Kingfisher' pollen abortion were investigated using transmission electron microscopy, RNA sequencing, isobaric tags for relative and absolute quantitation, and bioinformatics. It was found that the meiosis of microspore mother cells was normal before the tetrad stage, the microspores began to degenerate at the tetrad stage, and no microspores were observed in the anthers after the tetrad stage. In addition, transcriptomic and proteomic analyses showed that some genes and proteins related to cell senescence and death were identified to be implicated in chrysanthemum pollen abortion. These results indicated that the tetrad stage was the main period of pollen abortion, and the genes and proteins related to cell senescence and death contributed to pollen abortion. These add to our understanding of chrysanthemum pollen abortion and will be helpful for development of flowers without pollen contamination in the future.

Small RNA and degradome sequencing used to elucidate the basis of tolerance to salinity and alkalinity in wheat

BACKGROUND

Soil salinity and/or alkalinity impose a major constraint over crop yield and quality. An understanding of the molecular basis of the plant response to these stresses could inform the breeding of more tolerant varieties. The bread wheat cultivar SR3 exhibits an enhanced level of salinity tolerance, while SR4 is distinguished by its superior tolerance of alkalinity.

RESULTS

The small RNA and degradome sequencing was used to explore the miRNAs and corresponding targets associated with the superior stress tolerance of the SR lines. An examination of the small RNA content of these two closely related lines revealed the presence of 98 known and 219 novel miRNA sequences. Degradome libraries were constructed in order to identify the targets of the miRNAs, leading to the identification of 58 genes targeted by 26 of the known miRNAs and 549 targeted by 65 of the novel ones. The function of two of the stress-responsive miRNAs was explored using virus-induced gene silencing.

CONCLUSIONS

This analysis indicated that regulation mediated by both auxin and epigenetic modification can be important in determining both salinity and alkalinity tolerance, while jasmonate signaling and carbohydrate metabolism are important for salinity tolerance, as is proton transport for alkalinity tolerance.

Uncovering key small RNAs associated with gametocidal action in wheat

Gametocidal (Gc) chromosomes can kill gametes that lack them by causing chromosomal breakage to ensure their preferential transmission, and they have been exploited in genetic breeding. The present study investigated the possible roles of small RNAs (sRNAs) in Gc action. By sequencing two small RNA libraries from the anthers of Triticum aestivum cv. Chinese Spring (CS) and the Chinese Spring-Gc 3C chromosome monosomic addition line (CS-3C), we identified 239 conserved and 72 putative novel miRNAs, including 135 differentially expressed miRNAs. These miRNAs were predicted to target multiple genes with various molecular functions relevant to the features of Gc action, including sterility and genome instability. The transgenic overexpression of miRNA, which was up-regulated in CS-3C, reduced rice fertility. The CS-3C line exhibited a genome-wide reduction in 24 nt siRNAs compared with that of the CS line, particularly in transposable element (TE) and repetitive DNA sequences. Corresponding to this reduction, the bisulfite sequencing analysis of four retro-TE sequences showed a decrease in CHH methylation, typical of RNA-directed DNA methylation (RdDM). These results demonstrate that both miRNA-directed regulation of gene expression and siRNA-directed DNA methylation of target TE loci could play a role in Gc action.

Transcriptome analysis reveals new microRNAs-mediated pathway involved in anther development in male sterile wheat

Background

337S is a novel bi-pole-photo-thermo-sensitive genic male sterile line in wheat, and sensitive to both long day length/high temperature and short day length/low temperature condition. Although the regulatory function of MicroRNAs (miRNAs) in reproductive development has been increasingly studied, their roles in pre-meiotic and meiotic cells formation of plants have not been clearly explored. Here, we explored the roles of miRNAs in regulating male sterility of 337S at short day length/low temperature condition.

Results

Small RNA sequencing and degradome analyses were employed to identify miRNAs and their targets in the 337S whose meiotic cells collapsed rapidly during male meiotic prophase, resulting in failure of meiosis at SL condition. A total of 102 unique miRNAs were detected. Noticeably, the largest miRNA family was MiR1122. The target CCR4-associated factor 1 (CAF1) of miR2275, a subunit of the Carbon Catabolite Repressed 4-Negative on TATA-less (CCR4-NOT) complex, contributes to the process of early meiosis, and was first identified here. Further studies showed that the expression of several pivotal anther-related miRNAs was altered in 337S at SL condition, especially tae-miR1127a, which may be related to male sterility of 337S. Here, we also identified a new member of SWI/SNF factors SWI/SNF-related matrix-associated actin-dependent regulator of chromatin subfamily A, member 3-like 3 (SMARCA3L3) targeted by tae-miR1127a, whose function might be involved in faithful progression of meiosis in male reproductive cells.

Conclusion

The miRNA-target interactions of tae-miR2275-CAF1 and tae-miR1127a-SMARCA3L3 might be involved in regulating male fertility in 337S. Our results also implied that multiple roles for SMARCA3L3 and CAF1 in DNA repair and transcriptional regulation jointly orchestrated a tight and orderly system for maintaining chromatin and genome integrity during meiosis.

Electronic supplementary material

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

Genome-Wide Identification of the Alba Gene Family in Plants and Stress-Responsive Expression of the Rice Alba Genes

Architectural proteins play key roles in genome construction and regulate the expression of many genes, albeit the modulation of genome plasticity by these proteins is largely unknown. A critical screening of the architectural proteins in five crop species, viz., Oryza sativa, Zea mays, Sorghum bicolor, Cicer arietinum, and Vitis vinifera, and in the model plant Arabidopsis thaliana along with evolutionary relevant species such as Chlamydomonas reinhardtii, Physcomitrella patens, and Amborella trichopoda, revealed 9, 20, 10, 7, 7, 6, 1, 4, and 4 Alba (acetylation lowers binding affinity) genes, respectively. A phylogenetic analysis of the genes and of their counterparts in other plant species indicated evolutionary conservation and diversification. In each group, the structural components of the genes and motifs showed significant conservation. The chromosomal location of the Alba genes of rice (OsAlba), showed an unequal distribution on 8 of its 12 chromosomes. The expression profiles of the OsAlba genes indicated a distinct tissue-specific expression in the seedling, vegetative, and reproductive stages. The quantitative real-time PCR (qRT-PCR) analysis of the OsAlba genes confirmed their stress-inducible expression under multivariate environmental conditions and phytohormone treatments. The evaluation of the regulatory elements in 68 Alba genes from the 9 species studied led to the identification of conserved motifs and overlapping microRNA (miRNA) target sites, suggesting the conservation of their function in related proteins and a divergence in their biological roles across species. The 3D structure and the prediction of putative ligands and their binding sites for OsAlba proteins offered a key insight into the structure-function relationship. These results provide a comprehensive overview of the subtle genetic diversification of the OsAlba genes, which will help in elucidating their functional role in plants.

Regulation of low temperature stress in plants by microRNAs

Low temperature is one of the most common environmental stresses that seriously affect the growth and development of plants. However, plants have the plasticity in their defence mechanisms enabling them to tolerate and, sometimes, even survive adverse environmental conditions. MicroRNAs (miRNAs) are small non-coding RNAs, approximately 18-24 nucleotides in length, and are being increasingly recognized as regulators of gene expression at the post-transcriptional level and have the ability to influence a broad range of biological processes. There is growing evidence in the literature that reprogramming of gene expression mediated through miRNAs is a major defence mechanism in plants enabling them to respond to stresses. To date, numerous studies have established the importance of miRNA-based regulation of gene expression under low temperature stress. Individual miRNAs can modulate the expression of multiple mRNA targets, and, therefore, the manipulation of a single miRNA has the potential to affect multiple biological processes. Numerous functional studies have attempted to identify the miRNA-target interactions and have elaborated the role of several miRNAs in cold-stress regulation. This review summarizes the current understanding of miRNA-mediated modulation of the expression of key genes as well as genetic and regulatory pathways, involved in low temperature stress responses in plants.

Identification and molecular characterization of a trans-acting small interfering RNA producing locus regulating leaf rust responsive gene expression in wheat (Triticum aestivum L.)

A novel leaf rust responsive ta-siRNA-producing locus was identified in wheat showing similarity to 28S rRNA and generated four differentially expressing ta-siRNAs by phasing which targeted stress responsive genes. Trans-acting-small interfering RNAs (Ta-siRNAs) are plant specific molecules generally involved in development and are also stress responsive. Ta-siRNAs identified in wheat till date are all responsive to abiotic stress only. Wheat cultivation is severely affected by rusts and leaf rust particularly affects grain filling. This study reports a novel ta-siRNA producing locus (TAS) in wheat which is a segment of 28S ribosomal RNA but shows differential expression during leaf rust infestation. Four small RNA libraries prepared from wheat Near Isogenic Lines were treated with leaf rust pathogen and compared with untreated controls. A TAS with the ability to generate four ta-siRNAs by phasing events was identified along with the microRNA TamiR16 as the phase initiator. The targets of the ta-siRNAs included α-gliadin, leucine rich repeat, trans-membrane proteins, glutathione-S-transferase, and fatty acid desaturase among others, which are either stress responsive genes or are essential for normal growth and development of plants. Expression of the TAS, its generated ta-siRNAs, and their target genes were profiled at five different time points after pathogen inoculation of susceptible and resistant wheat isolines and compared with mock-inoculated controls. Comparative analysis of expression unveiled differential and reciprocal relationship as well as discrete patterns between susceptible and resistant isolines. The expression profiles of the target genes of the identified ta-siRNAs advocate more towards effector triggered susceptibility favouring pathogenesis. The study helps in discerning the functions of wheat genes regulated by ta-siRNAs in response to leaf rust.