MiR396 regulatory network and its expression during grain development in wheat

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.

BpGRP1 acts downstream of BpmiR396c/BpGRF3 to confer salt tolerance in Betula platyphylla

Glycine-rich RNA-binding proteins (GRPs) have been implicated in the responses of plants to environmental stresses, but the function of GRP genes involved in salt stress and the underlying mechanism remain unclear. In this study, we identified BpGRP1 (glycine-rich RNA-binding protein), a Betula platyphylla gene that is induced under salt stress. The physiological and molecular responses to salt tolerance were investigated in both BpGRP1-overexpressing and suppressed conditions. BpGRF3 (growth-regulating factor 3) was identified as a regulatory factor upstream of BpGRP1. We demonstrated that overexpression of BpGRF3 significantly increased the salt tolerance of birch, whereas the grf3-1 mutant exhibited the opposite effect. Further analysis revealed that BpGRF3 and its interaction partner, BpSHMT, function upstream of BpGRP1. We demonstrated that BpmiR396c, as an upstream regulator of BpGRF3, could negatively regulate salt tolerance in birch. Furthermore, we uncovered evidence showing that the BpmiR396c/BpGRF3 regulatory module functions in mediating the salt response by regulating the associated physiological pathways. Our results indicate that BpmiR396c regulates the expression of BpGRF3, which plays a role in salt tolerance by targeting BpGRP1.

Endoplasmic reticulum stress-responsive microRNAs are involved in the regulation of abiotic stresses in wheat

KEY MESSAGE

ER stress-responsive miRNAs, tae-miR164, tae-miR2916, and tae-miR396e-5p, are essential in response to abiotic stress. Investigating ER stress-responsive miRNAs is necessary to improve plant tolerance to environmental stress. MicroRNAs (miRNAs) play vital regulatory roles in plant responses to environmental stress. Recently, the endoplasmic reticulum (ER) stress pathway, an essential signalling pathway in plants in response to adverse conditions, has been widely studied in model plants. However, miRNAs associated with ER stress response remain largely unknown. Using high-throughput sequencing, three ER stress-responsive miRNAs, tae-miR164, tae-miR2916, and tae-miR396e-5p were identified, and their target genes were confirmed. These three miRNAs and their target genes actively responded to dithiothreitol, polyethylene glycol, salt, heat, and cold stresses. Furthermore, in some instances, the expression patterns of the miRNAs and their corresponding target genes were contrasting. Knockdown of tae-miR164, tae-miR2916, or tae-miR396e-5p using a barley stripe mosaic virus-based miRNA silencing system substantially enhanced the tolerance of wheat plants to drought, salt, and heat stress. Under conditions involving these stresses, inhibiting the miR164 function by using the short tandem target mimic approach in Arabidopsis thaliana resulted in phenotypes consistent with those of miR164-silenced wheat plants. Correspondingly, overexpression of tae-miR164 in Arabidopsis resulted in a decreased tolerance to drought stress and, to some extent, a decrease in tolerance to salt and high temperature. These results revealed that tae-miR164 plays a negative regulatory role in wheat/Arabidopsis in response to drought, salt, and heat stress. Taken together, our study provides new insights into the regulatory role of ER stress-responsive miRNAs in abiotic stress responses.

Screening of Reference miRNA of Different Early- and Late-Flowering Tree Peony Varieties

miRNA plays an important role in plant growth and development and in response to various stresses. Quantitative real-time PCR (qRT-PCR) technology is often used to detect the expression level of miRNAs and genes by comparing with reference genes. In order to screen out the optimal reference miRNAs in different tree peony varieties, the petals of 42 different early- and late-flowering tree peony varieties were used as experimental materials, and geNorm, NormFinder, Bestkeeper, and RefFinder software were used to evaluate the stability of 16 candidate reference miRNAs. The results showed that the average Ct values of all candidate reference miRNAs were between 15.34 ± 0.29 and 32.64 ± 0.38. The optimal number of reference miRNAs was four, which were PsPC-5p-19095, PsPC-3p-51259, PsmiR159a, and PsPC-3p-6660 in geNorm. The stability of PsPC-3p-6660 was the highest in the analysis results of NormFinder software. Among the analysis results of Bestkeeper software, PsMIR319-p5 has the highest stability. Among the results of comprehensive evaluation and analysis of several software using RefFinder, the candidate reference miRNA with the highest stability was PsPC-3p-6660. When PsPC-3p-6660 was used as the reference miRNA, the expression of PomiR171 and PomiR414 in response to different flowering times of tree peony was relatively stable in 42 tree peony varieties, indicating that PsPC-3p-6660 was stable and reliable. The results of this study provide a reference miRNA for studying the expression changes of miRNA in different tree peony varieties and further exploring the regulatory mechanism of miRNA in different peony varieties.

PeGRF6-PeGIF1 complex regulates cell proliferation in the leaf of Phalaenopsis equestris

Phalaenopsis equestris is an ornamental plant with very large leaves. In this study, we identified genes related to the regulation of leaf development in Phalaenopsis and explored their mechanism of action. Sequence alignment and phylogenetic analyses revealed that PeGRF6 in the PeGRF family of P. equestris has similarities with the Arabidopsis genes AtGRF1 and AtGRF2, which are known to be involved in the regulation of leaf development. Among the PeGRFs, PeGRF6 was continuously and stably expressed at various stages of leaf development. The functions of PeGRF6 and of its complex formed with PeGIF1 in leaf development were verified by virus-induced gene silencing (VIGS) technology. The results show that the PeGRF6-PeGIF1 complex forms in the nucleus and positively regulates leaf cell proliferation via influencing cell size. Interestingly, VIGS suppression of PeGRF6 resulted in anthocyanin accumulation in Phalaenopsis leaves. Analyses of the regulatory mechanism of the miR396-PeGRF6 model based on the P. equestris small RNA library constructed here suggested that PeGRF6 transcripts are cleaved by Peq-miR396. These results show that, compared with PeGRF6 or PeGIF1 alone, the PeGRF6-PeGIF1 complex plays a more important role in the leaf development of Phalaenopsis, possibly by regulating the expression of cell cycle-related genes.

Molecular prospective on the wheat grain development

Wheat grain development is an important biological process to determine grain yield and quality, which is controlled by the interplay of genetic, epigenetic, and environmental factors. Wheat grain development has been extensively characterized at the phenotypic and genetic levels. The advent of innovative molecular technologies allows us to characterize genes, proteins, and regulatory factors involved in wheat grain development, which have enhanced our understanding of the wheat seed development process. However, wheat is an allohexaploid with a large genome size, the molecular mechanisms underlying the wheat grain development have not been well understood as those in diploids. Understanding grain development, and how it is regulated, is of fundamental importance for improving grain yield and quality through conventional breeding or genetic engineering. Herein, we review the current discoveries on the molecular mechanisms underlying wheat grain development. Notably, only a handful of genes that control wheat grain development have, thus far, been well characterized, their interplay underlying the grain development remains elusive. The synergistic network-integrated genomics and epigenetics underlying wheat grain development and how the subgenome divergence dynamically and precisely regulates wheat grain development are unknown.

Genome‑wide identification and characterization of miR396 family members and their target genes GRF in sorghum (Sorghum bicolor (L.) moench)

MicroRNAs (miRNAs) widely participate in plant growth and development. The miR396 family, one of the most conserved miRNA families, remains poorly understood in sorghum. To reveal the evolution and expression pattern of Sbi-miR396 gene family in sorghum, bioinformatics analysis and target gene prediction were performed on the sequences of the Sbi-miR396 gene family members. The results showed that five Sbi-miR396 members, located on chromosomes 4, 6, and 10, were identified at the whole-genome level. The secondary structure analysis showed that the precursor sequences of all five Sbi-miR396 potentially form a stable secondary stem-loop structure, and the mature miRNA sequences were generated on the 5' arm of the precursors. Sequence analysis identified the mature sequences of the five sbi-miR396 genes were high identity, with differences only at the 1st, 9th and 21st bases at the 5' end. Phylogenetic analysis revealed that Sbi-miR396a, Sbi-miR396b, and Sbi-miR396c were clustered into Group I, and Sbi-miR396d and Sbi-miR396e were clustered into Group II, and all five sbi-miR396 genes were closely related to those of maize and foxtail millet. Expression analysis of different tissue found that Sbi-miR396d/e and Sbi-miR396a/b/c were preferentially and barely expressed, respectively, in leaves, flowers, and panicles. Target gene prediction indicates that the growth-regulating factor family members (SbiGRF1/2/3/4/5/6/7/8/10) were target genes of Sbi-miR396d/e. Thus, Sbi-miR396d/e may affect the growth and development of sorghum by targeting SbiGRFs. In addition, expression analysis of different tissues and developmental stages found that all Sbi-miR396 target genes, SbiGRFs, were barely expressed in leaves, root and shoot, but were predominantly expressed in inflorescence and seed development stage, especially SbiGRF1/5/8. Therefore, inhibition the expression of sbi-miR396d/e may increase the expression of SbiGRF1/5/8, thereby affecting floral organ and seed development in sorghum. These findings provide the basis for studying the expression of the Sbi-mir396 family members and the function of their target genes.

Multiple Functions of MiRNAs in Brassica napus L

The worldwide climate changes every year due to global warming, waterlogging, drought, salinity, pests, and pathogens, impeding crop productivity. Brassica napus is one of the most important oil crops in the world, and rapeseed oil is considered one of the most health-beneficial edible vegetable oils. Recently, miRNAs have been found and confirmed to control the expression of targets under disruptive environmental conditions. The mechanism is through the formation of the silencing complex that mediates post-transcriptional gene silencing, which pairs the target mRNA and target cleavage and/or translation inhibition. However, the functional role of miRNAs and targets in B. napus is still not clarified. This review focuses on the current knowledge of miRNAs concerning development regulation and biotic and abiotic stress responses in B. napus. Moreover, more strategies for miRNA manipulation in plants are discussed, along with future perspectives, and the enormous amount of transcriptome data available provides cues for miRNA functions in B. napus. Finally, the construction of the miRNA regulatory network can lead to the significant development of climate change-tolerant B. napus through miRNA manipulation.

Analysis of the MIR396 gene family and the role of MIR396b in regulating fiber length in cotton

Cotton fiber is one of the most important natural raw materials in the world textile industry. Improving fiber yield and quality has always been the main goal. MicroRNAs, as typical small noncoding RNAs, could affect fiber length during different stages of fiber development. Based on differentially expressed microRNA in the two interspecific backcross inbred lines (BILs) with a significant difference in fiber length, we identified the miR396 gene family in the two tetraploid cotton genomes and found MIR396b_D13 as the functional precursor to produce mature miR396 during the fiber elongation stage. Among 46 target genes regulated by miR396b, the GROWTH-REGULATING FACTOR 5 gene (GRF5, Gh_A10G0492) had a differential expression level in the two BILs during fiber elongation stage. The expression patterns indicated that the miR396b-GRF5 regulatory module has a critical role in fiber development. Furthermore, virus-induced gene silencing (VIGS) of miR396b significantly produced longer fiber than the wild type, and the expression level of GRF5 showed the reverse trends of the miR396b expression level. The analysis of co-expression network for the GRF5 gene suggested that a cytochrome P450 gene functions as an allene oxide synthase (Gh_D06G0089, AOS), which plays a critical role in jasmonate biosynthetic pathway. In conclusion, our results revealed that the miR396b-GRF5 module has a critical role in fiber development. These findings provide a molecular foundation for fiber quality improvement in the future.

Identification of long non-coding RNAs and microRNAs involved in anther development in the tropical Camellia oleifera

BACKGROUND

Explored the molecular science of anther development is important for improving productivity and overall yield of crops. Although the role of regulatory RNAs, including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), in regulating anther development has been established, their identities and functions in Camellia oleifera, an important industrial crop, have yet not been clearly explored. Here, we report the identification and characterization of genes, lncRNAs and miRNAs during three stages of the tropical C. oleifera anther development by single-molecule real-time sequencing, RNA sequencing and small RNA sequencing, respectively.

RESULTS

These stages, viz. the pollen mother cells stage, tetrad stage and uninucleate pollen stage, were identified by analyzing paraffin sections of floral buds during rapid expansion periods. A total of 18,393 transcripts, 414 putative lncRNAs and 372 miRNAs were identified, of which 5,324 genes, 115 lncRNAs, and 44 miRNAs were differentially accumulated across three developmental stages. Of these, 44 and 92 genes were predicted be regulated by 37 and 30 differentially accumulated lncRNAs and miRNAs, respectively. Additionally, 42 differentially accumulated lncRNAs were predicted as targets of 27 miRNAs. Gene ontology enrichment indicated that potential target genes of lncRNAs were enriched in photosystem II, regulation of autophagy and carbohydrate phosphatase activity, which are essential for anther development. Functional annotation of genes targeted by miRNAs indicated that they are relevant to transcription and metabolic processes that play important roles in microspore development. An interaction network was built with 2 lncRNAs, 6 miRNAs and 10 mRNAs. Among these, miR396 and miR156 family were up-regulated, while their targets, genes (GROWTH REGULATING FACTORS and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes) and lncRNAs, were down-regulated. Further, the trans-regulated targets of these lncRNAs, like wall-associated kinase2 and phosphomannose isomerase1, are involved in pollen wall formation during anther development.

CONCLUSIONS

This study unravels lncRNAs, miRNAs and miRNA-lncRNA-mRNA networks involved in development of anthers of the tropical C. oleifera lays a theoretical foundation for further elucidation of regulatory roles of lncRNAs and miRNAs in anther development.

Integrated mRNA and Small RNA Sequencing Reveals a microRNA Regulatory Network Associated with Starch Biosynthesis in Lotus (Nelumbo nucifera Gaertn.) Rhizomes

Internode starch biosynthesis is one of the most important traits in lotus rhizome because of its relation to crop productivity. Understanding the microRNA (miRNA) and mRNA expression profiles related to lotus internode starch biosynthesis would help develop molecular improvement strategies, but they are not yet well-investigated. To identify genes and miRNAs involved in internode starch biosynthesis, the cDNA and small RNA libraries of Z6-1, Z6-2, and Z6-3 were sequenced, and their expression were further studied. Through combined analyses of transcriptome data and small RNA sequencing data, a complex co-expression regulatory network was constructed, in which 20 miRNAs could modulate starch biosynthesis in different internodes by tuning the expression of 10 target genes. QRT-PCR analysis, transient co-expression experiment and dual luciferase assay comprehensively confirmed that NnumiR396a down-regulated the expression of NnSS2 and ultimately prevents the synthesis of amylopectin, and NnumiR396b down-regulated the expression of NnPGM2 and ultimately prevents the synthesis of total starch. Our results suggest that miRNAs play a critical role in starch biosynthesis in lotus rhizome, and that miRNA-mediated networks could modulate starch biosynthesis in this tissue. These results have provided important insights into the molecular mechanism of starch biosynthesis in developing lotus rhizome.

N6-Methyladenosine dynamic changes and differential methylation in wheat grain development

More methylation changes occur in late interval than in early interval of wheat seed development with protein and the starch synthesis-related pathway enriched in the later stages. Wheat seed development is a critical process to determining wheat yield and quality, which is controlled by genetics, epigenetics and environments. The N6-methyladenosine (m⁶A) modification is a reversible and dynamic process and plays regulatory role in plant development and stress responses. To better understand the role of m⁶A in wheat grain development, we characterized the m⁶A modification at 10 day post-anthesis (DPA), 20 DPA and 30 DPA in wheat grain development. m⁶A-seq identified 30,615, 30,326, 27,676 high confidence m⁶A peaks from the 10DPA, 20DPA, and 30DPA, respectively, and enriched at 3'UTR. There were 29,964, 29,542 and 26,834 unique peaks identified in AN0942_10d, AN0942_20d and AN0942_30d. One hundred and forty-two genes were methylated by m⁶A throughout seed development, 940 genes methylated in early grain development (AN0942_20d vs AN0942_10d), 1542 genes in late grain development (AN0942_30d vs AN0942_20d), and 1190 genes between early and late development stage (AN0942_30d vs AN0942_10d). KEGG enrichment analysis found that protein-related pathways and the starch synthesis-related pathway were significantly enriched in the later stages of seed development. Our results provide novel knowledge on m⁶A dynamic changes and its roles in wheat grain development.

Plant RNA-mediated gene regulatory network

Not all transcribed RNAs are protein-coding RNAs. Many of them are non-protein-coding RNAs in diverse eukaryotes. However, some of them seem to be non-functional and are resulted from spurious transcription. A lot of non-protein-coding transcripts have a significant function in the translation process. Gene expressions depend on complex networks of diverse gene regulatory pathways. Several non-protein-coding RNAs regulate gene expression in a sequence-specific system either at the transcriptional level or post-transcriptional level. They include a significant part of the gene expression regulatory network. RNA-mediated gene regulation machinery is evolutionarily ancient. They well-evolved during the evolutionary time and are becoming much more complex than had been expected. In this review, we are trying to summarizing the current knowledge in the field of RNA-mediated gene silencing.

Identification and profiling of conserved microRNAs in different developmental stages of crown imperial (Fritillaria imperialis L.) using high-throughput sequencing

BACKGROUND

Novel strategies for improvement of ornamental plants and their properties relay on miRNA control of differential plant gene expression modulation. Still, in response to the same abiotic stresses, some conserved miRNA families show different expression patterns in different plant species. In parallel, the use of deep sequencing technologies reveals new levels of complexity of regulatory networks in plants through identification of new miRNAs.

METHODS AND RESULTS

Fritillaria imperialis plants were collected from their natural habitats in Koohrang, Chaharmahal va Bakhtiari, Iran. Several tissues including stamen, pistil, petal, sepal, leaf, stem, bulb and fruit were collected during three developmental stages (stem elongation, flower development and seed head stages). Using RNAseq and qRT-PCR approach, this research revealed 21 conserved miRNAs, matching 15 miRNA families, in Fritilaria imperialis.

CONCLUSIONS

The expression of seven conserved miRNAs (Fim-miR156b, Fim-miR159, Fim-miR166a-5p, Fim-miR169d-5p, Fim-miR171c, Fim-miR393 and Fim-miR396e-3p) was further investigated in different tissues and three developmental stages, suggesting different roles for these miRNAs during growth and development of crown imperial. Gained knowledge from this research can open the door to find efficient ways to secure crown imperial survival, preservation and utilization and if proven useful may be applied in other plant species as well.

Small RNAs, Degradome, and Transcriptome Sequencing Provide Insights into Papaya Fruit Ripening Regulated by 1-MCP

As an inhibitor of ethylene receptors, 1-methylcyclopropene (1-MCP) can delay the ripening of papaya. However, improper 1-MCP treatment will cause a rubbery texture in papaya. Understanding of the underlying mechanism is still lacking. In the present work, a comparative sRNA analysis was conducted after different 1-MCP treatments and identified a total of 213 miRNAs, of which 44 were known miRNAs and 169 were novel miRNAs in papaya. Comprehensive functional enrichment analysis indicated that plant hormone signal pathways play an important role in fruit ripening. Through the comparative analysis of sRNAs and transcriptome sequencing, a total of 11 miRNAs and 12 target genes were associated with the ethylene and auxin signaling pathways. A total of 1741 target genes of miRNAs were identified by degradome sequencing, and nine miRNAs and eight miRNAs were differentially expressed under the ethylene and auxin signaling pathways, respectively. The network regulation diagram of miRNAs and target genes during fruit ripening was drawn. The expression of 11 miRNAs and 12 target genes was verified by RT-qPCR. The target gene verification showed that cpa-miR390a and cpa-miR396 target CpARF19-like and CpERF RAP2-12-like, respectively, affecting the ethylene and auxin signaling pathways and, therefore, papaya ripening.