Evolutionary conservation of microRNA regulatory programs in plant flower development

Comprehensive, Genome-Wide Identification and Expression Analyses of Phenylalanine Ammonia-Lyase Family under Abiotic Stresses in Brassica oleracea

Phenylalanine ammonia-lyase (PAL) acts as the rate-limiting enzyme for anthocyanin biosynthesis through the phenylpropanoid pathway, a crucial component of plant secondary metabolism. The PAL gene family plays a crucial role in plants' defense and stress responses, but its in silico identification and expression analyses in Brassica oleracea under different abiotic stresses remain unexplored. In this study, nine BolPAL, seven BrPAL, four AtPAL, and seventeen BnPAL genes were obtained from the genomes of B. oleracea, Brassica rapa, Arabidopsis thaliana, and Brassica napus, respectively. Segmental duplication and purifying selection are the causes of the BolPAL gene's amplification and evolution. The BolPAL genes with comparable intron-exon architectures and motifs were grouped together in the same clade. Three categories comprised the cis-regulatory elements: abiotic stressors, phytohormones, and light. According to the results of the qRT-PCR experiments, the majority of the BolPAL genes were expressed highly under MeJA, a low temperature, and a high temperature, and they were downregulated under ABA. Under white light (100 µmol m-2 s-1) with 50, 100, or 150 µmol m-2 s-1 far-red (FR), only a small number of the PAL genes were expressed at 50 and 100 µmol m-2 s-1 FR, while the majority of the PAL genes were slightly elevated at 150 µmol m-2 s-1 FR. This work offers a theoretical foundation for molecular breeding research to investigate the role of BolPAL genes and their role in anthocyanin biosynthesis.

Widespread changes to the translational landscape in a maize microRNA biogenesis mutant

MicroRNAs are short, non-coding RNAs that repress gene expression in both plants and animals and have diverse functions related to growth, development, and stress responses. The ribonuclease, DICER-LIKE1 (DCL1) is required for two steps in plant miRNA biogenesis: cleavage of the primary miRNAs (pri-miRNAs) to release a hairpin structure, called the precursor miRNA (pre-miRNA) and cleavage of the pre-miRNA to generate the miRNA/miRNA* duplex. The mature miRNA guides the RNA-induced silencing complex to target RNAs with complementary sequences, resulting in translational repression and/or RNA cleavage of target mRNAs. However, the relative contribution of translational repression versus mRNA degradation by miRNAs remains unknown at the genome-level in crops, especially in maize. The maize fuzzy tassel (fzt) mutant contains a hypomorphic mutation in DCL1 resulting in broad developmental defects. While most miRNAs are reduced in fzt, the levels of miRNA-targeted mRNAs are not dramatically increased, suggesting that translational regulation by miRNAs may be common. To gain insight into the repression mechanism of plant miRNAs, we combined ribosome profiling and RNA-sequencing to globally survey miRNA activities in maize. Our data indicate that translational repression contributes significantly to regulation of most miRNA targets and that approximately one-third of miRNA targets are regulated primarily at the translational level. Surprisingly, ribosomes appear altered in fzt mutants suggesting that DCL1 may also have a role in ribosome biogenesis. Thus, DICER-LIKE1 shapes the translational landscape in plants through both miRNA-dependent and miRNA-independent mechanisms.

Recent advancements in the physiological, genetic, and genomic research on Rhododendrons for trait improvement

High species diversity, hybridization potential, broad geographical dispersal range and ornamental characteristics (i.e., attractive size, shape, structure, flowers, and evergreen) have fetched a good international market for Rhododendron. However, most species are restricted to specific geographic areas due to their habitat specificity in acidic soil and cold climates, resulting many species being classified under threat categories of the IUCN. In this review, advances in research on Rhododendron for improvement to floral display quality and stress resistance have been described. The low genetic barrier among species has created opportunities for extensive hybridization and ploidy alteration for introducing quality and adaptive traits during the development of new varieties. Recent technological advances have supported investigations into the mechanism of flower development, as well as cold tolerance and pathogen resistance mechanisms in the Rhododendron. However, most of the species have limited adaptability to drought, line-tolerance, pathogen resistance, and high-temperature conditions and this resistance ability present in few species largely remains unexplored. Additionally, the available genetic diversity and genomic information on species, and possibilities for their application in molecular breeding have been summarized. Overall, genomic resource data are scarce in the majority of the members of this genus. Finally, various research gaps such as genetic mapping of quality traits, understanding the molecular mechanism of quality-related traits and genomic assortment in Rhododendron members have been discussed in the future perspective section.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04006-6.

Integrated analysis of miRNAs, transcriptome and phytohormones in the flowering time regulatory network of tea oil camellia

Camellia oleifera is a crucial cash crop in the southern region of China. Timely flowering is a crucial characteristic for maximizing crop productivity. Nevertheless, the cold temperature and wet weather throughout the fall and winter seasons in South China impact the timing of flowering and the yield produced by C. oleifera. This study examined the miRNAs, transcriptomes, and phytohormones that are part of the flowering time regulatory networks in distinct varieties of C. oleifera (Sep, Oct, and Nov). This study provides evidence that phytohormones significantly impact the timing of flowering in C. oleifera leaves. There is a positive correlation between the accumulation variations of zeatin (cZ), brassinolide (BL), salicylic acid (SA), 1-amino cyclopropane carboxylic acid (ACC), and jasmonic acid (JA) and flowering time. This means that blooming occurs earlier when the quantity of these substances in leaves increases. Abscisic acid (ABA), trans-zeatin-riboside (tZR), dihydrozeatin (dh-Z), and IP (N6-Isopentenyladenine) exhibit contrasting effects. Furthermore, both miR156 and miR172 play a crucial function in regulating flowering time in C. oleifera leaves by modulating the expression of SOC1, primarily through the miR156-SPL and miR172-AP2 pathways. These findings establish a strong basis for future research endeavors focused on examining the molecular network associated with the flowering period of C. oleifera and controlling flowering time management through external treatments.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01473-2.

Genome-Wide Identification and Characterization of miRNAs and Natural Antisense Transcripts Show the Complexity of Gene Regulatory Networks for Secondary Metabolism in Aristolochia contorta

Aristolochia contorta Bunge is an academically and medicinally important plant species. It belongs to the magnoliids, with an uncertain phylogenetic position, and is one of the few plant species lacking a whole-genome duplication (WGD) event after the angiosperm-wide WGD. A. contorta has been an important traditional Chinese medicine material. Since it contains aristolochic acids (AAs), chemical compounds with nephrotoxity and carcinogenicity, the utilization of this plant has attracted widespread attention. Great efforts are being made to increase its bioactive compounds and reduce or completely remove toxic compounds. MicroRNAs (miRNAs) and natural antisense transcripts (NATs) are two classes of regulators potentially involved in metabolism regulation. Here, we report the identification and characterization of 223 miRNAs and 363 miRNA targets. The identified miRNAs include 51 known miRNAs belonging to 20 families and 172 novel miRNAs belonging to 107 families. A negative correlation between the expression of miRNAs and their targets was observed. In addition, we identified 441 A. contorta NATs and 560 NAT-sense transcript (ST) pairs, of which 12 NATs were targets of 13 miRNAs, forming 18 miRNA-NAT-ST modules. Various miRNAs and NATs potentially regulated secondary metabolism through the modes of miRNA-target gene-enzyme genes, NAT-STs, and NAT-miRNA-target gene-enzyme genes, suggesting the complexity of gene regulatory networks in A. contorta. The results lay a solid foundation for further manipulating the production of its bioactive and toxic compounds.

Mulberry-derived miR168a downregulates BmMthl1 to promote physical development and fecundity in silkworms

Plant-derived miRNAs and their interactions with host organisms are considered important factors in regulating host physiological processes. In this study, we investigated the interaction between the silkworm, an oligophagous insect, and its primary food source, mulberry, to determine whether mulberry-derived miRNAs can penetrate silkworm cells and regulate their functions. Our results demonstrated that miR168a from mulberry leaves enters the silkworm hemolymph and binds to the silkworm Argonaute1 BmAGO1, which is transported via vesicles secreted by silkworm cells to exert its regulatory functions. In vivo and in vitro functional studies revealed that miR168a targets the mRNA of silkworm G protein-coupled receptor, BmMthl1, thereby inhibiting its expression and activating the JNK-FoxO pathway. This activation reduces oxidative stress responses, prolongs the lifespan of silkworms, and improves their reproductive capacity. These findings highlight the challenges of replacing mulberry leaves with alternative protein sources and provide a foundation for developing silkworm germplasms suitable for factory rearing.

Genome-Wide Identification and Expression Analyses of the FAR1/FHY3 Gene Family Provide Insight into Inflorescence Development in Maize

As transcription factors derived from transposase, FAR-RED IMPAIRED RESPONSE1 (FAR1) and its homolog FHY3 play crucial roles in the regulation of light signaling and various stress responses by coordinating the expression of downstream target genes. Despite the extensive investigation of the FAR1/FHY3 family in Arabidopsis thaliana and other species, a comprehensive examination of these genes in maize has not been conducted thus far. In this study, we employed a genomic mining approach to identify 16 ZmFAR1 genes in the maize inbred line B73, which were further classified into five subgroups based on their phylogenetic relationships. The present study characterized the predicted polypeptide sequences, molecular weights, isoelectric points, chromosomal distribution, gene structure, conserved motifs, subcellular localizations, phylogenetic relationships, and cis-regulatory elements of all members belonging to the ZmFAR1 family. Furthermore, the tissue-specific expression of the 16 ZmFAR1 genes was analyzed using RNA-seq, and their expression patterns under far-red light conditions were validated in the ear and tassel through qRT-qPCR. The observed highly temporal and spatial expression patterns of these ZmFAR1 genes were likely associated with their specific functional capabilities under different light conditions. Further analysis revealed that six ZmFAR1 genes (ZmFAR1-1, ZmFAR1-10, ZmFAR1-11, ZmFAR1-12, ZmFAR1-14, and ZmFAR1-15) exhibited a response to simulated shading treatment and actively contributed to the development of maize ears. Through the integration of expression quantitative trait loci (eQTL) analyses and population genetics, we identified the presence of potential causal variations in ZmFAR1-14 and ZmFAR1-9, which play a crucial role in regulating the kernel row number and kernel volume weight, respectively. In summary, this study represents the initial identification and characterization of ZmFAR1 family members in maize, uncovering the functional variation in candidate regulatory genes associated with the improvement of significant agronomic traits during modern maize breeding.

MiRNAs differentially expressed in vegetative and reproductive organs of Marchantia polymorpha - insights into their expression pattern, gene structures and function

MicroRNAs regulate gene expression affecting a variety of plant developmental processes. The evolutionary position of Marchantia polymorpha makes it a significant model to understand miRNA-mediated gene regulatory pathways in plants. Previous studies focused on conserved miRNA-target mRNA modules showed their critical role in Marchantia development. Here, we demonstrate that the differential expression of conserved miRNAs among land plants and their targets in selected organs of Marchantia additionally underlines their role in regulating fundamental developmental processes. The main aim of this study was to characterize selected liverwort-specific miRNAs, as there is a limited knowledge on their biogenesis, accumulation, targets, and function in Marchantia. We demonstrate their differential accumulation in vegetative and generative organs. We reveal that all liverwort-specific miRNAs examined are encoded by independent transcriptional units. MpmiR11737a, MpmiR11887 and MpmiR11796, annotated as being encoded within protein-encoding genes, have their own independent transcription start sites. The analysis of selected liverwort-specific miRNAs and their pri-miRNAs often reveal correlation in their levels, suggesting transcriptional regulation. However, MpmiR11796 shows a reverse correlation to its pri-miRNA level, suggesting post-transcriptional regulation. Moreover, we identify novel targets for selected liverwort-specific miRNAs and demonstrate an inverse correlation between their expression and miRNA accumulation. In the case of one miRNA precursor, we provide evidence that it encodes two functional miRNAs with two independent targets. Overall, our research sheds light on liverwort-specific miRNA gene structure, provides new data on their biogenesis and expression regulation. Furthermore, identifying their targets, we hypothesize the potential role of these miRNAs in early land plant development and functioning.

Small RNA profiling reveals that an ovule-specific microRNA, cja-miR5179, targets a B-class MADS-box gene in Camellia japonica

BACKGROUND AND AIMS

The functional specialization of microRNA and its target genes is often an important factor in the establishment of spatiotemporal patterns of gene expression that are essential to plant development and growth. In different plant lineages, understanding the functional conservation and divergence of microRNAs remains to be explored.

METHODS

To identify small regulatory RNAs underlying floral patterning, we performed a tissue-specific profiling of small RNAs in various floral organs from single and double flower varieties (flowers characterized by multiple layers of petals) in Camellia japonica. We identified cja-miR5179, which belongs to a deeply conserved microRNA family that is conserved between angiosperms and basal plants but frequently lost in eudicots. We characterized the molecular function of cja-miR5179 and its target - a B-function MADS-box gene - through gene expression analysis and transient expression assays.

KEY RESULTS

We showed that cja-miR5179 is exclusively expressed in ovule tissues at the early stage of floral development. We found that cja-miR5179 targets the coding sequences of a DEFICIENS-like B-class gene (CjDEF) mRNA, which is located in the K motif of the MADS-box domain; and the target sites of miR5179/MADS-box were consistent in Camellia and orchids. Furthermore, through a petal transient-expression assay, we showed that the BASIC PENTACYSTEINE proteins bind to the GA-rich motifs in the cja-miR5179 promoter region and suppresses its expression.

CONCLUSIONS

We propose that the regulation between miR5179 and a B-class MADS-box gene in C. japonica has a deep evolutionary origin before the separation of monocots and dicots. During floral development of C. japonica, cja-miR5179 is specifically expressed in the ovule, which may be required for the inhibition of CjDEF function. This work highlights the evolutionary conservation as well as functional divergence of small RNAs in floral development.

MIR159 regulates multiple aspects of stamen and carpel development and requires dissection and delimitation of differential downstream regulatory network for manipulating fertility traits

Unravelling genetic networks regulating developmental programs are key to devising and implementing genomics assisted trait modification strategies. It is crucial to understand the role of small RNAs, and the basis of their ability to modify traits. MIR159 has been previously reported to cause defects in anther development in Arabidopsis; however, the complete spectrum and basis of the defects remained unclear. The present study was therefore undertaken to comprehensively investigate the role of miR159 from Brassica juncea in modulating vegetative and reproductive traits. Owing to the polyploid nature of Brassica, paralogous and homeologous copies of MIR159A, MIR159B, and, MIR159C were identified and analysis of the precursor uncovered extensive structural and sequence variation. The MIR159 locus with mature miR159 with perfect target complimentarily with MYB65, was cloned from Brassica juncea var. Varuna for functional characterization by generating constitutively over-expressing lines in Arabidopsis thaliana Col-0. Apart from statistically significant difference in multiple vegetative traits, drastic differences were observed in stamen and pistil. Over-expression of miR159a led to shortening of filament length and loss of tetradynamous condition. Anthers were apiculate, with improper lobe formation, and unsynchronized cellular growth between connective tissue and another lobe development. Analysis revealed arrested meiosis/cytokinesis in microspores, and altered lignin deposition pattern in endothecial walls thus affecting anther dehiscence. In the gynoecium, flaccid, dry stigmatic papillae, and large embryo sac in the female gametophyte was observed. Over-expression of miR159a thus severely affected pollination and seed-set. Analysis of the transcriptome data revealed components of regulatory networks of anther and carpel developmental pathway, and lignin metabolism that are affected. Expression analysis allowed us to position the miR159a-MYB65 module in the genetic network of stamen development, involved in pollen-grain maturation; in GA-mediated regulation of stamen development, and in lignin metabolism. The study, on one hand indicates role of miR159a-MYB65 in regulating multiple aspects of reproductive organ development that can be manipulated for trait modification, but also raises several unaddressed questions such as relationship between miR159a and male-meiosis, miR159a and filament elongation for future investigations. Accession numbers: KC204951-KC204960. Project number PRJNA1035268.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01377-7.

The miR167-OsARF12 module regulates rice grain filling and grain size downstream of miR159

Grain weight and quality are always determined by grain filling. Plant microRNAs have drawn attention as key targets for regulation of grain size and yield. However, the mechanisms that underlie grain size regulation remain largely unclear because of the complex networks that control this trait. Our earlier studies demonstrated that suppressed expression of miR167 (STTM/MIM167) substantially increased grain weight. In a field test, the yield increased up to 12.90%-21.94% because of a significantly enhanced grain filling rate. Here, biochemical and genetic analyses revealed the regulatory effects of miR159 on miR167 expression. Further analysis indicated that OsARF12 is the major mediator by which miR167 regulates rice grain filling. Overexpression of OsARF12 produced grain weight and grain filling phenotypes resembling those of STTM/MIM167 plants. Upon in-depth analysis, we found that OsARF12 activates OsCDKF;2 expression by directly binding to the TGTCGG motif in its promoter region. Flow cytometry analysis of young panicles from OsARF12-overexpressing plants and examination of cell number in cdkf;2 mutants verified that OsARF12 positively regulates grain filling and grain size by targeting OsCDKF;2. Moreover, RNA sequencing results suggested that the miR167-OsARF12 module is involved in the cell development process and hormone pathways. OsARF12-overexpressing plants and cdkf;2 mutants exhibited enhanced and reduced sensitivity to exogenous auxin and brassinosteroid (BR) treatment, confirming that targeting of OsCDKF;2 by OsARF12 mediates auxin and BR signaling. Our results reveal that the miR167-OsARF12 module works downstream of miR159 to regulate rice grain filling and grain size via OsCDKF;2 by controlling cell division and mediating auxin and BR signals.

Zinc oxide nanoparticles (ZnONPs) influenced seed development, grain quality, and remobilization by affecting the transcription of microRNA 171 (miR171), miR156, NAM, and SUT genes in wheat (Triticum aestivum): a biological advantage and risk assessment study

Limited studies have been conducted on the role of microRNAs (miRs) and transcription factors in regulating plant cell responses to nanoparticles. This study attempted to address whether the foliar application of zinc oxide nanoparticles (ZnONPs; 0, 10, 25, and 50 mgL^-1) can affect miRs, gene expression, and wheat grain quality. The seedlings were sprayed with ZnONPs (0, 10, 25, and 50 mgL^-1) or bulk counterpart (BZnO) five times at 72 h intervals. The application of ZnONPs at 10 mgL^-1 increased the number of spikelets and seed weight, while the nano-supplement at 50 mgL^-1 was accompanied by severe restriction on developing spikes and grains. ZnONPs, in a dose-dependent manner, transcriptionally influenced miR156 and miR171. The expression of miR171 showed a similar trend to that of miR156. The ZnONPs at optimum concentration upregulated the NAM transcription factor and sucrose transporter (SUT) at transcriptional levels. However, the transcription of both NAM and SUT genes displayed a downward trend in response to the toxic dose of ZnONPs (50 mgL^-1). Utilization of ZnONPs increased proline and total soluble phenolic content. Monitoring the accumulation of carbohydrates, including fructan, glucose, fructose, and sucrose, revealed that ZnONPs at 10 mgL^-1 modified the source/sink communication and nutrient remobilization. The molecular and physiological data revealed that the expression of miR156 and miR171 is tightly linked to seed grain development, remobilization of carbohydrates, and genes involved in nutrient transportation. This study establishes a novel strategy for obtaining higher yields in crops. This biological risk assessment investigation also displays the potential hazard of applying ZnONPs at the flowering developmental phase.

Transgene Bioconfinement: Don't Flow There

The adoption of genetically engineered (GE) crops has led to economic and environmental benefits. However, there are regulatory and environmental concerns regarding the potential movement of transgenes beyond cultivation. These concerns are greater for GE crops with high outcrossing frequencies to sexually compatible wild relatives and those grown in their native region. Newer GE crops may also confer traits that enhance fitness, and introgression of these traits could negatively impact natural populations. Transgene flow could be lessened or prevented altogether through the addition of a bioconfinement system during transgenic plant production. Several bioconfinement approaches have been designed and tested and a few show promise for transgene flow prevention. However, no system has been widely adopted despite nearly three decades of GE crop cultivation. Nonetheless, it may be necessary to implement a bioconfinement system in new GE crops or in those where the potential of transgene flow is high. Here, we survey such systems that focus on male and seed sterility, transgene excision, delayed flowering, as well as the potential of CRISPR/Cas9 to reduce or eliminate transgene flow. We discuss system utility and efficacy, as well as necessary features for commercial adoption.

Identification and analysis of miRNAs differentially expressed in male and female Trichosanthes kirilowii maxim

Trichosanthes kirilowii Maxim. (TK) is a dioecious plant in the Cucurbitaceae family of which different sexes have separate medicinal uses. We used Illumina high-throughput sequencing technology to sequence miRNAs from male and female flower buds of TK. We performed bioinformatics analysis, miRNA identification, and target gene prediction on the data obtained from sequencing, and association analysis was performed in combination with the results of a previous transcriptome sequencing study. As a result, there were 80 differentially expressed miRNAs (DESs) between the female and male plants (48 upregulated and 32 downregulated in female plants). Moreover, 27 novel miRNAs in DESs were predicted to have 282 target genes, and 51 known miRNAs were predicted to have 3418 target genes. By establishing a regulatory network between miRNAs and target genes, 12 core genes were screened, including 7 miRNAs and 5 target genes. Among them, tkmiR157a-5p, tkmiR156c, tkmiR156_2, and tkmiR156k_2 jointly target the regulation of tkSPL18 and tkSPL13B. These two target genes are specifically expressed in male and female plants, respectively, and are involved in the biosynthesis process of BR, which is closely related to the sex differentiation process of TK. The identification of these miRNAs will provide a reference for the analysis of the sex differentiation mechanism of TK.

Plant non-coding RNAs function in pollen development and male sterility

Male sterility is classified as either cytoplasmic male sterility (CMS) or genic male sterility (GMS). Generally, CMS involves mitochondrial genomes interacting with the nuclear genome, while GMS is caused by nuclear genes alone. Male sterility is regulated by multilevel mechanisms in which non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and phased small interfering RNAs (phasiRNAs), which have been proven to be critical elements. The development of high-throughput sequencing technology offers new opportunities to evaluate the genetic mechanism of ncRNAs in plant male sterility. In this review, we summarize the critical ncRNAs that regulate gene expression in ways dependent on or independent of hormones, which involve the differentiation of the stamen primordia, degradation of the tapetum, formation of microspores, and the release of pollen. In addition, the key mechanisms of the miRNA-lncRNA-mRNA interaction networks mediating male sterility in plants are elaborated. We present a different perspective on exploring the ncRNA-mediated regulatory pathways that control CMS in plants and create male-sterile lines through hormones or genome editing. A refined understanding of the ncRNA regulatory mechanisms in plant male sterility for the development of new sterile lines would be conducive to improve hybridization breeding.

AGO1 and HSP90 buffer different genetic variants in Arabidopsis thaliana

Argonaute 1 (AGO1), the principal protein component of microRNA-mediated regulation, plays a key role in plant growth and development. AGO1 physically interacts with the chaperone HSP90, which buffers cryptic genetic variation in plants and animals. We sought to determine whether genetic perturbation of AGO1 in Arabidopsis thaliana would also reveal cryptic genetic variation, and if so, whether AGO1-dependent loci overlap with those dependent on HSP90. To address these questions, we introgressed a hypomorphic mutant allele of AGO1 into a set of mapping lines derived from the commonly used Arabidopsis strains Col-0 and Ler. Although we identified several cases in which AGO1 buffered genetic variation, none of the AGO1-dependent loci overlapped with those buffered by HSP90 for the traits assayed. We focused on 1 buffered locus where AGO1 perturbation uncoupled the traits days to flowering and rosette leaf number, which are otherwise closely correlated. Using a bulk segregant approach, we identified a nonfunctional Ler hua2 mutant allele as the causal AGO1-buffered polymorphism. Introduction of a nonfunctional hua2 allele into a Col-0 ago1 mutant background recapitulated the Ler-dependent ago1 phenotype, implying that coupling of these traits involves different molecular players in these closely related strains. Taken together, our findings demonstrate that even though AGO1 and HSP90 buffer genetic variation in the same traits, these robustness regulators interact epistatically with different genetic loci, suggesting that higher-order epistasis is uncommon. Plain Language Summary Argonaute 1 (AGO1), a key player in plant development, interacts with the chaperone HSP90, which buffers environmental and genetic variation. We found that AGO1 buffers environmental and genetic variation in the same traits; however, AGO1-dependent and HSP90-dependent loci do not overlap. Detailed analysis of a buffered locus found that a nonfunctional HUA2 allele decouples days to flowering and rosette leaf number in an AGO1-dependent manner, suggesting that the AGO1-dependent buffering acts at the network level.

The Integrated mRNA and miRNA Approach Reveals Potential Regulators of Flowering Time in Arundina graminifolia

Orchids are among the most precious flowers in the world. Regulation of flowering time is one of the most important targets to enhance their ornamental value. The beauty of Arundina graminifolia is its year-round flowering, although the molecular mechanism of this flowering ability remains masked. Therefore, we performed a comprehensive assessment to integrate transcriptome and miRNA sequencing to disentangle the genetic regulation of flowering in this valuable species. Clustering analyses provided a set of molecular regulators of floral transition and floral morphogenesis. We mined candidate floral homeotic genes, including FCA, FPA, GI, FT, FLC, AP2, SOC1, SVP, GI, TCP, and CO, which were targeted by a variety of miRNAs. MiR11091 targeted the highest number of genes, including candidate regulators of phase transition and hormonal control. The conserved miR156-miR172 pathway of floral time regulation was evident in our data, and we found important targets of these miRNAs in the transcriptome. Moreover, endogenous hormone levels were determined to decipher the hormonal control of floral buds in A. graminifolia. The qRT-PCR analysis of floral and hormonal integrators validated the transcriptome expression. Therefore, miRNA-mediated mining of candidate genes with hormonal regulation forms the basis for comprehending the complex regulatory network of perpetual flowering in precious orchids. The findings of this study can do a great deal to broaden the breeding programs for flowering time manipulation of orchids.

Identification of differentially expressed miRNAs and their target genes in response to brassinolide treatment on flowering of tree peony (Paeonia ostii)

Tree peony is a famous flower plant in China, but the short and concentrated flowering period limits its ornamental value and economic value. Brassinolide (BR) plays an important role in plant growth and development including flowering. There have been a large number of reports on the molecular aspects of the flowering process, but the genetic mechanism that was responsible for miRNA-guided regulation of tree peony is almost unclear. In this study, the leaves of tree peony cultivar, 'Feng Dan', were sprayed with different concentrations of BR, and the obvious bloom delay was found at the treatment with BR 50 μg/L. The small RNA sequencing and transcriptome sequencing were performed on the petals of tree peony under an untreated control (CK) and the treatment with BR 50 μg/L during four consecutive flowering development stages. A total of 22 known miRNAs belonging to 12 families were identified and 84 novel miRNAs were predicted. Combined with transcriptome data, a total of 376 target genes were predicted for the 18 differentially expressed known miRNAs and 177 target genes were predicted for the 23 differentially expressed novel miRNAs. Additionally, the potential miRNAs and their target genes were identified, including miR156b targeting SPL, miR172a_4 targeting AP2 and four novel miRNAs targeting SPA1, and revealed that they might affect the flowering time in tree peony. Collectively, these results would provide a theoretical basis for further analysis of miRNA-guided regulation on flowering period in tree peony.

Characteristics and potential functional effects of long insertions in Asian butternuts

Background

Structural variants (SVs) play important roles in adaptation evolution and species diversification. Especially, in plants, many phenotypes of response to the environment were found to be associated with SVs. Despite the prevalence and significance of SVs, long insertions remain poorly detected and studied in all but model species.

Results

We used whole-genome resequencing of paired reads from 80 Asian butternuts to detect long insertions and further analyse their characteristics and potential functional effects. By combining of mapping-based and de novo assembly-based methods, we obtained a multiple related species pangenome representing higher taxonomic groups. We obtained 89,312 distinct contigs totaling 147,773,999 base pair (bp) of new sequences, of which 347 were putative long insertions placed in the reference genome. Most of the putative long insertions appeared in multiple species; in contrast, only 62 putative long insertions appeared in one species, which may be involved in the response to the environment. 65 putative long insertions fell into 61 distinct protein-coding genes involved in plant development, and 105 putative long insertions fell into upstream of 106 distinct protein-coding genes involved in cellular respiration. 3,367 genes were annotated in 2,606 contigs. We propose PLAINS (https://github.com/CMB-BNU/PLAINS.git), a streamlined, comprehensive pipeline for the prediction and analysis of long insertions using whole-genome resequencing.

Conclusions

Our study lays down an important foundation for further whole-genome long insertion studies, allowing the investigation of their effects by experiments.

The miR156-Targeted SQUAMOSA PROMOTER BINDING PROTEIN (PmSBP) Transcription Factor Regulates the Flowering Time by Binding to the Promoter of SUPPRESSOR OF OVEREXPRESSION OF CO1 (PmSOC1) in Prunus mume

Prunus mume, a famous perennial ornamental plant and fruit tree in Asia, blooms in winter or early spring in the Yangtze River area. The flowering time directly determines its ornamental and economic value, so it is of great significance to study the molecular mechanism of flowering time. SQUAMOSA PROMOTER BINDING PROTEIN (SBP), often regulated by miR156, is an important flowering regulator, although its function is unknown in P. mume. Here, 11 miR156 precursors were analyzed and located in five chromosomes of the P. mume genome. The expression pattern showed that PmSBP1/6 was negatively correlated with miR156. The promoters of PmSBP1/6 were specifically expressed in the apical meristem. Overexpression of PmSBP1/6 in tobacco promoted flowering and changed the length ratio of pistil and stamen. Moreover, PmSBP1 also affected the number and vitality of pollen and reduced the fertility of transgenic tobacco. Furthermore, ectopic expression of PmSBP1/6 caused up-regulated expression of endogenous SUPPRESSOR OF OVEREXPRESSION OF CO1 (NtSOC1). The yeast-one hybrid assay showed that PmSBP1 was bonded to the promoters of PmSOC1s. In conclusion, a miR156-PmSBP1-PmSOC1s pathway was formed to participate in the regulation of flowering time in P. mume, which provided references for the molecular mechanism of flowering time regulation and molecular breeding of P. mume.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

The molecular mechanism of microRNA duplex selectivity of Arabidopsis ARGONAUTE10

Small RNAs (sRNAs), including microRNAs (miRNAs) and small interfering RNAs (siRNAs), are essential gene regulators for plant and animal development. The loading of sRNA duplexes into the proper ARGONAUTE (AGO) protein is a key step to forming a functional silencing complex. In Arabidopsis thaliana, the specific loading of miR166/165 into AGO10 (AtAGO10) is critical for the maintenance of the shoot apical meristem, the source of all shoot organs, but the mechanism by which AtAGO10 distinguishes miR166/165 from other cellular miRNAs is not known. Here, we show purified AtAGO10 alone lacks loading selectivity towards miR166/165 duplexes. However, phosphate and HSP chaperone systems reshape the selectivity of AtAGO10 to its physiological substrates. A loop in the AtAGO10 central cleft is essential for recognizing specific mismatches opposite the guide strand 3' region in miR166/165 duplexes. Replacing this loop with the equivalent loop from Homo sapiens AGO2 (HsAGO2) changes AtAGO10 miRNA loading behavior such that 3' region mismatches are ignored and mismatches opposite the guide 5' end instead drive loading, as in HsAGO2. Thus, this study uncovers the molecular mechanism underlying the miR166/165 selectivity of AtAGO10, essential for plant development, and provides new insights into how miRNA duplex structures are recognized for sRNA sorting.

Identification of miRNAs and Target Genes at Key Stages of Sexual Differentiation in Androdioecious Osmanthus fragrans

Androdioecy is the crucial transition state in the evolutionary direction of hermaphroditism to dioecy, however, the molecular mechanisms underlying the formation of this sex system remain unclear. While popular in China for its ornamental and cultural value, Osmanthus fragrans has an extremely rare androdioecy breeding system, meaning that there are both male and hermaphroditic plants in a population. To unravel the mechanisms underlying the formation of androdioecy, we performed small RNA sequencing studies on male and hermaphroditic O. fragrans. A total of 334 miRNAs were identified, of which 59 were differentially expressed. Functional categorization revealed that the target genes of differentially expressed miRNAs were mainly involved in the biological processes of reproductive development and the hormone signal transduction pathway. We speculated that the miRNA160, miRNA167, miRNA393 and miRNA396 families may influence the sex differentiation in O. fragrans. Overall, our study is the first exploration of miRNAs in the growth and development process of O. fragrans, and is also the first study of androdioecious plants from the miRNA sequencing perspective. The analysis of miRNAs and target genes that may be involved in the sex differentiation process lay a foundation for the ultimate discovery of the androdioecious molecular mechanism in O. fragrans.

miR3633a-GA3ox2 Module Conducts Grape Seed-Embryo Abortion in Response to Gibberellin

Seedlessness is one of the important quality and economic traits favored by grapevine consumers, which are mainly affected by phytohormones, especially gibberellin (GA). GA is widely utilized in seedless berry production and could effectively induce grape seed embryo abortion. However, the molecular mechanism underlying this process, like the role of RNA silencing in the biosynthesis pathway of GA remains elusive. Here, Gibberellin 3-β dioxygenase2 (GA3ox2) as the last key enzyme in GA biosynthesis was predicated as a potential target gene for miR3633a, and two of them were identified as a GA response in grape berries. We also analyzed the promoter regions of genes encoding GA biosynthesis and found the hormone-responsive elements to regulate grape growth and development. The cleavage interaction between VvmiR3633a and VvGA3ox2 was validated by RLM-RACE and the transient co-transformation technique in tobacco in vivo. Interestingly, during GA-induced grape seed embryo abortion, exogenous GA promoted the expression of VvmiR3633a, thereby mainly repressing the level of VvGA3ox2 in seed embryos. We also observed a negative correlation between down-regulated VvGA20ox2/VvGA3ox2 and up-regulated VvGA2ox3/VvGA2ox1, of which GA inactivation was greater than GA synthesis, inhibited active GA content, accompanied by the reduction of VvSOD and VvCAT expression levels and enzymatic activities. These series of changes might be the main causes of grape seed embryo abortion. In conclusion, we have preliminarily drawn a schematic mode of GA-mediated VvmiR3633a and related genes regulatory network during grape seed abortion induced by exogenous GA. Our findings provide novel insights into the GA-responsive roles of the VvmiR3633a-VvGA3ox2 module in the modulation of grape seed-embryo abortion, which has implications for the molecular breeding of high-quality seedless grape berries.

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.

Evolutionary Analysis and Functional Identification of Clock-Associated PSEUDO-RESPONSE REGULATOR (PRRs) Genes in the Flowering Regulation of Roses

Pseudo-response regulators (PRRs) are the important genes for flowering in roses. In this work, clock PRRs were genome-wide identified using Arabidopsis protein sequences as queries, and their evolutionary analyses were deliberated intensively in Rosaceae in correspondence with angiosperms species. To draw a comparative network and flow of clock PRRs in roses, a co-expression network of flowering pathway genes was drawn using a string database, and their functional analysis was studied by silencing using VIGS and protein-to-protein interaction. We revealed that the clock PRRs were significantly expanded in Rosaceae and were divided into three major clades, i.e., PRR5/9 (clade 1), PRR3/7 (clade 2), and TOC1/PRR1 (clade 3), based on their phylogeny. Within the clades, five clock PRRs were identified in Rosa chinensis. Clock PRRs had conserved RR domain and shared similar features, suggesting the duplication occurred during evolution. Divergence analysis indicated the role of duplication events in the expansion of clock PRRs. The diverse cis elements and interaction of clock PRRs with miRNAs suggested their role in plant development. Co-expression network analysis showed that the clock PRRs from Rosa chinensis had a strong association with flowering controlling genes. Further silencing of RcPRR1b and RcPRR5 in Rosa chinensis using VIGS led to earlier flowering, confirming them as negative flowering regulators. The protein-to-protein interactions between RcPRR1a/RcPRR5 and RcCO suggested that RcPRR1a/RcPRR5 may suppress flowering by interfering with the binding of RcCO to the promoter of RcFT. Collectively, these results provided an understanding of the evolutionary profiles as well as the functional role of clock PRRs in controlling flowering in roses.

Unravelling the Role of Epigenetic Modifications in Development and Reproduction of Angiosperms: A Critical Appraisal

Epigenetics are the heritable changes in gene expression patterns which occur without altering DNA sequence. These changes are reversible and do not change the sequence of the DNA but can alter the way in which the DNA sequences are read. Epigenetic modifications are induced by DNA methylation, histone modification, and RNA-mediated mechanisms which alter the gene expression, primarily at the transcriptional level. Such alterations do control genome activity through transcriptional silencing of transposable elements thereby contributing toward genome stability. Plants being sessile in nature are highly susceptible to the extremes of changing environmental conditions. This increases the likelihood of epigenetic modifications within the composite network of genes that affect the developmental changes of a plant species. Genetic and epigenetic reprogramming enhances the growth and development, imparts phenotypic plasticity, and also ensures flowering under stress conditions without changing the genotype for several generations. Epigenetic modifications hold an immense significance during the development of male and female gametophytes, fertilization, embryogenesis, fruit formation, and seed germination. In this review, we focus on the mechanism of epigenetic modifications and their dynamic role in maintaining the genomic integrity during plant development and reproduction.

The RpTOE1-RpFT Module Is Involved in Rejuvenation during Root-Based Vegetative Propagation in Robinia pseudoacacia

Vegetative propagation is an important method of reproduction and rejuvenation in horticulture and forestry plants with a long lifespan. Although substantial juvenile clones have been obtained through the vegetative propagation of ornamental plants, the molecular factors that regulate rejuvenation during vegetative propagation are largely unknown. Here, root sprouting and root cutting of Robinia pseudoacacia were used as two vegetative propagation methods. From two consecutive years of transcriptome data from rejuvenated seedlings and mature trees, one gene module and one miRNA module were found to be specifically associated with rejuvenation during vegetative propagation through weighted gene co-expression network analysis (WGCNA). In the gene module, a transcription factor-encoding gene showed high expression during vegetative propagation, and it was subsequently named RpTOE1 through homology analysis. Heterologous overexpression of RpTOE1 in wild-type Arabidopsis and toe1 toe2 double mutants prolonged the juvenile phase. The qRT-PCR results predicted RpFT to be a downstream gene that was regulated by RpTOE1. Further investigation of the protein-DNA interactions using yeast one-hybrid, electrophoretic mobility shift, and dual luciferase reporter assays confirmed that RpTOE1 negatively regulated RpFT by binding directly to the TOE binding site (TBS)-like motif on its promoter. On the basis of these results, we showed that the high expression of RpTOE1 during vegetative propagation and its inhibition of RpFT played a key role in the phase reversal of R. pseudoacacia.

Temporal Dynamic Transcriptome Landscape Reveals Regulatory Network During the Early Differentiation of Female Strobilus Buds in Ginkgo biloba

Reproductive bud differentiation is one of the most critical events for the reproductive success of seed plants. Yet, our understanding of genetic basis remains limited for the development of the reproductive organ of gymnosperms, namely, unisexual strobilus or cone, leaving its regulatory network largely unknown for strobilus bud differentiation. In this study, we analyzed the temporal dynamic landscapes of genes, long non-coding RNAs (lncRNAs), and microRNAs (miRNAs) during the early differentiation of female strobilus buds in Ginkgo biloba based on the whole transcriptome sequencing. Results suggested that the functions of three genes, i.e., Gb_19790 (GbFT), Gb_13989 (GinNdly), and Gb_16301 (AG), were conserved in both angiosperms and gymnosperms at the initial differentiation stage. The expression of genes, lncRNAs, and miRNAs underwent substantial changes from the initial differentiation to the enlargement of ovule stalk primordia. Besides protein-coding genes, 364 lncRNAs and 15 miRNAs were determined to be functional. Moreover, a competing endogenous RNA (ceRNA) network comprising 10,248 lncRNA-miRNA-mRNA pairs was identified, which was highly correlated with the development of ovulate stalk primordia. Using the living fossil ginkgo as the study system, this study not only reveals the expression patterns of genes related to flowering but also provides novel insights into the regulatory networks of lncRNAs and miRNAs, especially the ceRNA network, paving the way for future studies concerning the underlying regulation mechanisms of strobilus bud differentiation.

An Integrated Analysis of Transcriptome and miRNA Sequencing Provides Insights into the Dynamic Regulations during Flower Morphogenesis in Petunia

Published genome sequences can facilitate multiple genome sequencing studies of flower development, which can serve as the basis for later analysis of variation in flower phenotypes. To identify potential regulators related to flower morphology, we captured dynamic expression patterns under five different developmental stages of petunia flowers, a popular bedding plant, using transcriptome and miRNA sequencing. The significant transcription factor (TF) families, including MYB, MADS, and bHLH, were elucidated. MADS-box genes exhibited co-expression patterns with BBR-BPC, GATA, and Dof genes in different modules according to a weighted gene co-expression network analysis. Through miRNA sequencing, a total of 45 conserved and 26 novel miRNAs were identified. According to GO and KEGG enrichment analysis, the carbohydrate metabolic process, photosynthesis, and phenylalanine metabolism were significant at the transcriptomic level, while the response to hormone pathways was significantly enriched by DEmiR-targeted genes. Finally, an miRNA–RNA network was constructed, which suggested the possibility of novel miRNA-mediated regulation pathways being activated during flower development. Overall, the expression data in the present study provide novel insights into the developmental gene regulatory network facilitated by TFs, miRNA, and their target genes.

Investigation of Roles of TaTALE Genes during Development and Stress Response in Bread Wheat

The three amino acid loop extension (TALE) genes of the homeobox superfamily are responsible for numerous biological functions in plants. Herein, we identified a total of 72 TaTALE genes in the allohexaploid genome of bread wheat (Triticum aestivum L.) and performed a comprehensive investigation for gene and protein structural properties, phylogeny, expression patterns, and multilevel gene regulations. The identified TaTALE proteins were further classified into two groups, TaBLHs and TaKNOXs, which were tightly clustered into the phylogeny. The negative Ka/Ks ratio of duplicated genes suggested purifying selection pressure with confined functional divergence. Various signature domains and motifs were found conserved in both groups of proteins. The occurrence of diverse cis-regulatory elements and modulated expression during various developmental stages and in the presence of abiotic (heat, drought, salt) and two different fungal stresses suggested their roles in development and stress response, as well. The interaction of TaTALEs with the miRNAs and other development-related homeobox proteins also suggested their roles in growth and development and stress response. The present study revealed several important aspects of TaTALEs that will be useful in further functional validation of these genes in future studies.

The Evolution of euAPETALA2 Genes in Vascular Plants: From Plesiomorphic Roles in Sporangia to Acquired Functions in Ovules and Fruits

The field of evolutionary developmental biology can help address how morphological novelties evolve, a key question in evolutionary biology. In Arabidopsis thaliana, APETALA2 (AP2) plays a role in the development of key plant innovations including seeds, flowers, and fruits. AP2 belongs to the AP2/ETHYLENE RESPONSIVE ELEMENT BINDING FACTOR family which has members in all viridiplantae, making it one of the oldest and most diverse gene lineages. One key subclade, present across vascular plants is the euAPETALA2 (euAP2) clade, whose founding member is AP2. We reconstructed the evolution of the euAP2 gene lineage in vascular plants to better understand its impact on the morphological evolution of plants, identifying seven major duplication events. We also performed spatiotemporal expression analyses of euAP2/TOE3 genes focusing on less explored vascular plant lineages, including ferns, gymnosperms, early diverging angiosperms and early diverging eudicots. Altogether, our data suggest that euAP2 genes originally contributed to spore and sporangium development, and were subsequently recruited to ovule, fruit and floral organ development. Finally, euAP2 protein sequences are highly conserved; therefore, changes in the role of euAP2 homologs during development are most likely due to changes in regulatory regions.

Dynamic expression of miRNAs and functional analysis of target genes involved in the response to male sterility of the wheat line YS3038

Thermosensitive cytoplasmic male sterile (TCMS) lines play an important role in wheat breeding, heterosis utilization, and germplasm innovation. MicroRNAs (miRNAs) can regulate the expression level of target genes by inhibiting the translation of these genes. YS3038 is a wheat TCMS line. In this study, the fertility conversion mechanism of YS3038 was studied by examining the abortion characteristics of YS3038, the regulation pattern of miRNAs and the target genes of miRNAs in YS3038. MiRNA-seq was performed on three important stages of YS3038 under sterile and fertile conditions. Then, the clean reads were aligned with some databases to filter other ncRNAs and repeats. The known miRNAs and novel miRNAs were predicted by sequence comparison with known miRNAs from miRbase. Differential expression of miRNAs between different stages and between different fertile conditions was analyzed, and functional analysis of target genes with opposite expression patterns as those of the miRNAs was conducted. The Ubisch bodies and microspores of sterile anthers were covered with filamentous materials. The degradation of the tapetum cells, the chloroplast structure of endothecium cells, and the microspore structure were abnormal. Microspore development was hindered from the late uninucleate stage to the binucleate stage. Twenty, 52, and 68 differentially expressed miRNAs (DEmiRs) were identified at the early uninucleate, late uninucleate, and binucleate stages, respectively, and there were 0, 7, and 72 differentially expressed target genes (DETGs), respectively, at these three stages. At the binucleate stage, 29 DEmiRs had 41 target mRNAs in total, and the expression patterns of the 41 target mRNAs were opposite to those of the 29 miRNAs. Fifteen significantly enriched KEGG pathways were associated with the 41 target mRNAs. Leucine-rich repeat receptor-like kinases (LRR-RLKs) play important roles in plant developmental and physiological processes. Some studies have shown that the expression of LRR-RLKs is related to the differentiation of microsporocytes and tapetum cells and to male sterility. An LRR-RLK (TaeRPK) gene was silenced by the barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) method, and the seed setting rates of the TaeRPK-silenced plants (3.51%) were significantly lower than those of the negative control plants (88.78%) (P < 0.01). Thus, the TaeRPK gene is likely to be involved in the fertility conversion of YS3038.

Analysis of miRNA-mediated regulation of flowering induction in Lilium × formolongi

BACKGROUND

MicroRNAs play pivotal roles in plant vegetative phase change and flowering induction via integrating into multiple flowering pathways. Lilium × formolongi is an important ornamental lily cultivar that can flower within one year after sowing. However, it remains unresolved how miRNA-mediated regulation networks contribute to the L. × formolongi characteristics of a short vegetative growth period and rapid flowering.

RESULTS

In this study, the small RNA libraries and one degradome library were constructed for L. × formolongi during vegetative growth and flowering initiation, and 366 conserved miRNAs and 32 novel miRNAs were identified. Additionally, 84 miRNAs were significantly differentially expressed during development. A total of 396 targets of 185 miRNAs were identified and validated through degradome sequencing. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that functions of the targets were top enriched in the cold and cadmium ion responses, pentose phosphate pathway and carbon fixation in photosynthetic organisms. Furthermore, among 23 differentially expressed miRNA-target pairs, the miR156s-LfSPL2, miR172a-LfAP2 and miR164a-LfNAC pairs as well as miR159a-LfSPL2 were found to be relevant to flowering based on the correlation analysis of expression profiles in the miRNA libraries, degradome and transcriptome. A coexpression regulatory network focused on differentially expressed pairs was also constructed by WGCNA, and 14 miRNAs were considered putative key miRNAs during vegetative development and flowering induction. miR156a/ d/ e showed particularly strong relationships with other miRNAs in the coexpression network.

CONCLUSIONS

This study provides cues for the further exploration of the regulatory mechanisms of short vegetative development and flowering in L. × formolongi.

Secrets of the MIR172 family in plant development and flowering unveiled

In plants, conserved microRNAs (miRNAs) tend to be encoded by gene families with multiple members. Two recent studies interrogated the functions of the 5-member MIR172 family in Arabidopsis and revealed complexities and intricacies of gene regulatory networks underlying floral transition.

Systematic analyses of the MIR172 family members of Arabidopsis define their distinct roles in regulation of APETALA2 during floral transition

MicroRNAs (miRNAs) play important roles in regulating flowering and reproduction of angiosperms. Mature miRNAs are encoded by multiple MIRNA genes that can differ in their spatiotemporal activities and their contributions to gene regulatory networks, but the functions of individual MIRNA genes are poorly defined. We functionally analyzed the activity of all 5 Arabidopsis thaliana MIR172 genes, which encode miR172 and promote the floral transition by inhibiting the accumulation of APETALA2 (AP2) and APETALA2-LIKE (AP2-LIKE) transcription factors (TFs). Through genome editing and detailed confocal microscopy, we show that the activity of miR172 at the shoot apex is encoded by 3 MIR172 genes, is critical for floral transition of the shoot meristem under noninductive photoperiods, and reduces accumulation of AP2 and TARGET OF EAT2 (TOE2), an AP2-LIKE TF, at the shoot meristem. Utilizing the genetic resources generated here, we show that the promotion of flowering by miR172 is enhanced by the MADS-domain TF FRUITFULL, which may facilitate long-term silencing of AP2-LIKE transcription, and that their activities are partially coordinated by the TF SQUAMOSA PROMOTER-BINDING-LIKE PROTEIN 15. Thus, we present a genetic framework for the depletion of AP2 and AP2-LIKE TFs at the shoot apex during floral transition and demonstrate that this plays a central role in floral induction.

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

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

A Composite Analysis of Flowering Time Regulation in Lettuce

Plants undergo profound physiological changes when transitioning from vegetative to reproductive growth. These changes affect crop production, as in the case of leafy vegetables. Lettuce is one of the most valuable leafy vegetable crops in the world. Past genetic studies have identified multiple quantitative trait loci (QTLs) that affect the timing of the floral transition in lettuce. Extensive functional molecular studies in the model organism Arabidopsis provide the opportunity to transfer knowledge to lettuce to explore the mechanisms through which genetic variations translate into changes in flowering time. In this review, we integrated results from past genetic and molecular studies for flowering time in lettuce with orthology and functional inference from Arabidopsis. This summarizes the basis for all known genetic variation underlying the phenotypic diversity of flowering time in lettuce and how the genetics of flowering time in lettuce projects onto the established pathways controlling flowering time in plants. This comprehensive overview reveals patterns across experiments as well as areas in need of further study. Our review also represents a resource for developing cultivars with delayed flowering time.

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

Backgrounds

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

Results

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

Conclusions

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

Supplementary Information

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

Genome-wide analysis of PmTCP4 transcription factor binding sites by ChIP-Seq during pistil abortion in Japanese apricot

The TCP4 transcription factor plays an important role in plant growth and development, especially in flower development. PmTCP4 is involved in the process of pistil abortion in Japanese apricot, but its molecular mechanism, particularly the DNA binding sites and co-regulatory genes, are quite unknown. Therefore, to identify the genome-wide binding sites of PmTCP4 transcription factors and their co-regulatory genes, chromatin immunoprecipitation sequencing (ChIP-Seq) was carried out. ChIP-Seq data produced the maximum enriched peaks in two Japanese apricot cultivars 'Daqiandi' (DQD) and 'Longyan' (LY), which showed that the majority of DNA-protein interactions are relevant and have a significant function in binding sites. Moreover, 720 and 251 peak-associated genes regulated by PmTCP4 were identified in DQD and LY, respectively, and most of them were involved in the flower and pistil development process. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that photosynthesis and oxidative phosphorylation were the most enriched pathways in both cultivars and all identified genes related to these pathways were down-regulated. This study will provide a reference for a better understanding of the PmTCP4 regulatory mechanism during pistil abortion in Japanese apricot.

MiR396 is involved in plant response to vernalization and flower development in Agrostis stolonifera

MicroRNA396 (miR396) has been demonstrated to regulate flower development by targeting growth-regulating factors (GRFs) in annual species. However, its role in perennial grasses and its potential involvement in flowering time control remain unexplored. Here we report that overexpression of miR396 in a perennial species, creeping bentgrass (Agrostis stolonifera L.), alters flower development. Most significantly, transgenic (TG) plants bypass the vernalization requirement for flowering. Gene expression analysis reveals that miR396 is induced by long-day (LD) photoperiod and vernalization. Further study identifies VRN1, VRN2, and VRN3 homologs whose expression patterns in wild-type (WT) plants are similar to those observed in wheat and barley during transition from short-day (SD) to LD, and SD to cold conditions. However, compared to WT controls, TG plants overexpressing miR396 exhibit significantly enhanced VRN1 and VRN3 expression, but repressed VRN2 expression under SD to LD conditions without vernalization, which might be associated with modified expression of methyltransferase genes. Collectively, our results unveil a potentially novel mechanism by which miR396 suppresses the vernalization requirement for flowering which might be related to the epigenetic regulation of VRN genes and provide important new insight into critical roles of a miRNA in regulating vernalization-mediated transition from vegetative to reproductive growth in monocots.

Computational Identification and Comparative Analysis of Conserved miRNAs and Their Putative Target Genes in the Juglans regia and J. microcarpa Genomes

MicroRNAs (miRNAs) are important factors for the post-transcriptional regulation of protein-coding genes in plants and animals. They are discovered either by sequencing small RNAs or computationally. We employed a sequence-homology-based computational approach to identify conserved miRNAs and their target genes in Persian (English) walnut, Juglans regia, and its North American wild relative, J. microcarpa. A total of 119 miRNA precursors (pre-miRNAs) were detected in the J. regia genome and 121 in the J. microcarpa genome and miRNA target genes were predicted and their functional annotations were performed in both genomes. In the J. regia genome, 325 different genes were targets; 87.08% were regulated by transcript cleavage and 12.92% by translation repression. In the J. microcarpa genome, 316 different genes were targets; 88.92% were regulated by transcript cleavage and 11.08% were regulated by translation repression. Totals of 1.3% and 2.0% of all resistance gene analogues (RGA) and 2.7% and 2.6% of all transcription factors (TFs) were regulated by miRNAs in the J. regia and J. microcarpa genomes, respectively. Juglans genomes evolved by a whole genome duplication (WGD) and consist of eight pairs of fractionated homoeologous chromosomes. Within each pair, the chromosome that has more genes with greater average transcription also harbors more pre-miRNAs and more target genes than its homoeologue. While only minor differences were detected in pre-miRNAs between the J. regia and J. microcarpa genomes, about one-third of the pre-miRNA loci were not conserved between homoeologous chromosome within each genome. Pre-miRNA and their corresponding target genes showed a tendency to be collocated within a subgenome.

OsmiR167a-targeted auxin response factors modulate tiller angle via fine-tuning auxin distribution in rice

Rice tiller angle determines plant growth density and further contributes grain production. Although a few genes have been characterized to regulate tiller angle in rice, the molecular mechanism underlying the control of tiller angle via microRNA is poorly understood. Here, we report that rice tiller angle is controlled by OsmiR167a-targeted auxin response factors OsARF12, OsARF17 and OsARF25. In the overexpression of OsMIR167a plants, the expression of OsARF12, OsARF17 and OsARF25 was severely repressed and displayed larger tiller angle as well as the osarf12/osarf17 and osarf12/ osarf25 plants. In addition, those plants showed compromised abnormal auxin distribution and less sensitive to gravity. We also demonstrate that OsARF12, OsARF17 and OsARF25 function redundantly and might be involved in HSFA2D and LAZY1-dependent asymmetric auxin distribution pathway to control rice tiller angle. Our results reveal that OsmiR167a represses its targets, OsARF12, OsARF17 and OsARF25, to control rice tiller angle by fine-tuning auxin asymmetric distribution in shoots.

TIAs pathway genes and associated miRNA identification in Vinca minor: supporting aspidosperma and eburnamine alkaloids linkage via transcriptomic analysis

V. minor contains monomeric eburnamine-type of indole alkaloids having utilization as a neuro-medicinal plant. The biosynthetic pathway studies using miRNAs has been the focal point for plant genomic research in recent years and this technique is utilized to get an insight into a possible pathway level study in V. minor as understanding of genes in this prized medicinal plant is meagrely understood. The de novo transcriptomic analysis using Illumina Next gen sequencing has been performed in glasshouse shifted plant and transformed roots to elucidate the possible non confirmed steps of terpenoid indole alkaloids (TIAs) pathway in V. minor. A putative TIA pathway is elucidated in the study including twelve possible TIAs biosynthetic genes. The specific miRNA associated with TIAs pathway were identified and their roles were discussed for the first time in V. minor. The comparative analysis of transcriptomic data of glasshouse shifted plant and transformed roots showed that the raw reads of transformed roots were higher (83,740,316) compared to glasshouse shifted plant (67,733,538). The EST-SSR prediction showed the maximum common repeats among glasshouse shifted plant and transformed roots, although small variation was found in trinucleotide repeats restricted to glasshouse shifted plant. The study reveals overall 37 miRNAs which were observed to be true and can have a role in pathway as they can regulate the growth and alkaloid production. The identification of putative pathway genes plays an important role in establishing linkage between Aspidosperma and Eburnamine alkaloids.

Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

Background

Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant.

Results

Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes.

Conclusions

Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

Ectopic expression of miRNA172 in tomato (Solanum lycopersicum) reveals novel function in fruit development through regulation of an AP2 transcription factor

BACKGROUND

MicroRNAs (miRNAs) are short non-coding RNAs that can influence gene expression via diverse mechanisms. Tomato is a fruit widely consumed for its flavor, culinary attributes, and high nutritional quality. Tomato fruit are climacteric and fleshy, and their ripening is regulated by endogenous and exogenous signals operating through a coordinated genetic network. Much research has been conducted on mechanisms of tomato fruit ripening, but the roles of miRNA-regulated repression/expression of specific regulatory genes are not well documented.

RESULTS

In this study, we demonstrate that miR172 specifically targets four SlAP2 transcription factor genes in tomato. Among them, SlAP2a was repressed by the overexpression of SlmiR172, manifesting in altered flower morphology, development and accelerated ripening. miR172 over-expression lines specifically repressed SlAP2a, enhancing ethylene biosynthesis, fruit color and additional ripening characteristics. Most previously described ripening-regulatory genes, including RIN-MADS, NR, TAGL1 and LeHB-1 were not influenced by miR172 while CNR showed altered expression.

CONCLUSIONS

Tomato fruit ripening is directly influenced by miR172 targeting of the APETALA2 transcription factor, SlAP2a, with minimal influence over additional known ripening-regulatory genes. miR172a-guided SlAP2a expression provides insight into another layer of genetic control of ripening and a target for modifying the quality and nutritional value of tomato and possibly other fleshy fruit crops.

The Critical Role of miRNAs in Regulation of Flowering Time and Flower Development

Flowering is an important biological process for plants that ensures reproductive success. The onset of flowering needs to be coordinated with an appropriate time of year, which requires tight control of gene expression acting in concert to form a regulatory network. MicroRNAs (miRNAs) are non-coding RNAs known as master modulators of gene expression at the post-transcriptional level. Many different miRNA families are involved in flowering-related processes such as the induction of floral competence, floral patterning, and the development of floral organs. This review highlights the diverse roles of miRNAs in controlling the flowering process and flower development, in combination with potential biotechnological applications for miRNAs implicated in flower regulation.