Genome-wide analysis of plant miRNA action clarifies levels of regulatory dynamics across developmental contexts

PmiR-Select® - a computational approach to plant pre-miRNA identification in genomes

Precursors of microRNAs (pre-miRNAs) are less used in silico to mine miRNAs. This study developed PmiR-Select® based on covariance models (CMs) to identify new pre-miRNAs, detecting conserved secondary structural features across RNA sequences and eliminating the redundancy. The pipeline preceded PmiR-Select® filtered 20% plant pre-miRNAs (from 38589 to 8677) from miRBase. The second filter reduced pre-miRNAs by 7% (from 8677 to 8045) through length limit to pre-miRNAs (70-300 nt) and miRNAs (20-24 nt). The 80% redundancy threshold was statistically the best, eliminating 55% pre-miRNAs (from 8045 to 3608). Angiosperms retained the highest number of pre-miRNAs and their families (2981 and 2202), followed by gymnosperms (362 and 271), bryophytes (183 and 119), and algae (82 and 78). Thirty-seven conserved pre-miRNA families happened among plant land clades, but none with algae. The PmiR-Select® was applied to the rice genome, producing 8536 pre-miRNAs from 36 families. The 80% redundancy threshold retained 3% pre-miRNAs (n = 264) from 36 families, valuable experimental and computational research resources. 14% (n = 1216) of 8536 were new pre-miRNAs from 19 new families in rice. Only 16 new sequences from six families overlapped (39 to 54% identities) with rice pre-miRNAs and five species on miRBase. The validation against mature miRNAs identified 8086 pre-miRNAs from 13 families. Eleven ones have already been recorded, but two new and abundant pre-miRNAs [miR437 (n = 296) and miR1435 (n = 725)] scattered in all 12-rice chromosomes. PmiR-Select® identified pre-miRNAs, decreased the redundancy, and discovered new miRNAs. These findings pave the way to delineating benchtop and computational experiments.

Quantitative and spatially resolved detection of multiplexed microRNA from plant tissue via hybridization to hydrogel-bound DNA probes in nanoliter well arrays

Understanding complex regulatory networks in plant systems requires elucidating the roles of various gene regulators under a spatial landscape. MicroRNA are key regulators that impart high information value through their tissue specificity and stability when using expression patterns for evaluating network outcomes. However, current techniques that utilize spatial multiplexing and quantitation of microRNA are limited to primarily mammalian systems. Here, we present a method to spatially resolve and quantify multiple endogenous microRNA in situ using ethanol fixed, paraffin embedded model plant species. This method utilizes target-specific microRNA capture along with universal ligating and labelling, all within functionalized hydrogel posts containing DNA probes in nanoliter well arrays. We demonstrate the platform's multiplexing capabilities through analyzing three endogenous microRNA in Arabidopsis thaliana rosettes which provide useful answers to fundamental plant growth and development from the unique expression patterns. The spatial tissue technique is also validated using non-spatial small RNA assays to demonstrate the versatility of the well array platform. Our new platform expands the toolkit of spatial omics technologies for plants.

A diffusible small-RNA-based Turing system dynamically coordinates organ polarity

The formation of a flat and thin leaf presents a developmentally challenging problem, requiring intricate regulation of adaxial-abaxial (top-bottom) polarity. The patterning principles controlling the spatial arrangement of these domains during organ growth have remained unclear. Here we show that this regulation in Arabidopsis thaliana is achieved by an organ-autonomous Turing reaction-diffusion system centred on mobile small RNAs. The data illustrate how Turing dynamics transiently instructed by prepatterned information is sufficient to self-sustain properly oriented polarity in a dynamic, growing organ, presenting intriguing parallels to left-right patterning in the vertebrate embryo. Computational modelling demonstrates that this self-organizing system continuously adapts to coordinate the robust planar polarity of a flat leaf while affording flexibility to generate the tissue patterns of evolutionarily diverse organ shapes. Our findings identify a small-RNA-based Turing network as a dynamic regulator of organ polarity that accounts for leaf shape diversity at the level of the individual organ, plant or species.

Plant virology: an RNA treasure trove

Key principles pertaining to RNA biology not infrequently have their origins in plant virology. Examples have arisen from studies on viral RNA-intrinsic properties and the infection process from gene expression, replication, movement, and defense evasion to biotechnological applications. Since RNA is at the core of the central dogma in molecular biology, how plant virology assisted in the reinforcement or adaptations of this concept, while at other instances shook up elements of the doctrine, is discussed. Moreover, despite the negative effects of viral diseases in agriculture worldwide, plant viruses can be considered a scientific treasure trove. Today they remain tools of discovery for biotechnology, studying evolution, cell biology, and host-microbe interactions.

Delineating the tissue-mediated drought stress governed tuning of conserved miR408 and its targets in rice

Engineering drought tolerance in rice needs to focus on regulators that enhance tolerance while boosting plant growth and vigor. The present study delineated the concealed function and tissue-mediated interplay of the miR408/target module in imparting drought stress tolerance in rice. The plant miR408 family comprises three dominant mature forms (21 nt), including a distinct monocot variant (F-7 with 5' C) and is divided into six groups. miR408 majorly cleaves genes belonging to the blue copper protein in addition to several other species-specific targets in plants. Comparative sequence analysis in 4726 rice accessions identified 22 sequence variants (SNP and InDELs) in its promoter (15) and pre-miR408 region. Haplotype analysis of the sequence variants indicated eight haplotypes (three: Japonica-specific and five: Indica-specific) of the miR408 promoter. In drought-tolerant Nagina 22, miR408 follows flag leaf preferential expression. Under drought conditions, its levels are upregulated in flag leaf and roots which seems to be regulated by a differential fraction of methylated cytosines (mCs) in the precursor region. The active pool of miR408 regulated targets under control and drought conditions is impacted by the tissue type. Comparative expression analysis of the miR408/target module under different sets of conditions features 83 targets exhibiting antagonistic expression in rice, out of which 12 genes, including four PLANTACYANINS (OsUCL6, 7, 9 and 30), PIRIN, OsLPR1, OsCHUP1, OsDOF12, OsBGLU1, glycine-rich cell wall gene, OsDUT, and OsERF7, are among the high confidence targets. Further, overexpression of MIR408 in drought-sensitive rice cultivar (PB1) leads to the massive enhancement of vegetative growth in rice with improved ETR and Y(II) and enhanced dehydration stress tolerance. The above results suggest that miR408 is likely to act as a positive regulator of growth and vigor, as well as dehydration stress, making it a potential candidate for engineering drought tolerance in rice.

A Wox3-patterning module organizes planar growth in grass leaves and ligules

Grass leaves develop from a ring of primordial initial cells within the periphery of the shoot apical meristem, a pool of organogenic stem cells that generates all of the organs of the plant shoot. At maturity, the grass leaf is a flattened, strap-like organ comprising a proximal supportive sheath surrounding the stem and a distal photosynthetic blade. The sheath and blade are partitioned by a hinge-like auricle and the ligule, a fringe of epidermally derived tissue that grows from the adaxial (top) leaf surface. Together, the ligule and auricle comprise morphological novelties that are specific to grass leaves. Understanding how the planar outgrowth of grass leaves and their adjoining ligules is genetically controlled can yield insight into their evolutionary origins. Here we use single-cell RNA-sequencing analyses to identify a 'rim' cell type present at the margins of maize leaf primordia. Cells in the leaf rim have a distinctive identity and share transcriptional signatures with proliferating ligule cells, suggesting that a shared developmental genetic programme patterns both leaves and ligules. Moreover, we show that rim function is regulated by genetically redundant Wuschel-like homeobox3 (WOX3) transcription factors. Higher-order mutations in maize Wox3 genes greatly reduce leaf width and disrupt ligule outgrowth and patterning. Together, these findings illustrate the generalizable use of a rim domain during planar growth of maize leaves and ligules, and suggest a parsimonious model for the homology of the grass ligule as a distal extension of the leaf sheath margin.

Discretizing the cellular bases of plant morphogenesis: Emerging properties from subcellular and noisy patterning

A fundamental question in biology is how multicellular organisms robustly shape their organs. In the past decade, much progress has been made not just in identifying biochemical and biophysical factors underpinning morphogenesis, but also in analyzing their spatio-temporal dynamics. A remarkable outcome of such analyses is that morphogenesis involves high levels of heterogeneity and fluctuations at local scales. Although this could be considered as white noise to be averaged over time, there is increasing evidence that these heterogeneities and fluctuations are instructive cues for development. In this review, we highlight some of the new questions that such heterogeneities raise for plant morphogenesis. We also investigate their effects across scales, focusing on how subcellular heterogeneities contribute to organ shape robustness and evolvability.

MicroRNAs balance growth and salt stress responses in sweet sorghum

Salt stress is one of the major causes of reduced crop production, limiting agricultural development globally. Plants have evolved with complex systems to maintain the balance between growth and stress responses, where signaling pathways such as hormone signaling play key roles. Recent studies revealed that hormones are modulated by microRNAs (miRNAs). Previously, two sweet sorghum (Sorghum bicolor) inbred lines with different salt tolerance were identified: the salt-tolerant M-81E and the salt-sensitive Roma. The levels of endogenous hormones in M-81E and Roma varied differently under salt stress, showing a different balance between growth and stress responses. miRNA and degradome sequencing showed that the expression of many upstream transcription factors regulating signal transduction and hormone-responsive genes was directly induced by differentially expressed miRNAs, whose levels were very different between the two sweet sorghum lines. Furthermore, the effects of representative miRNAs on salt tolerance in sorghum were verified through a transformation system mediated by Agrobacterium rhizogenes. Also, miR-6225-5p reduced the level of Ca²⁺ in the miR-6225-5p-overexpressing line by inhibiting the expression of the Ca²⁺ uptake gene SbGLR3.1 in the root epidermis and affected salt tolerance in sorghum. This study provides evidence for miRNA-mediated growth and stress responses in sweet sorghum.

Overexpression of circ PTK2 suppresses the progression of nonalcoholic fatty liver disease via the miR-200c/SIK2/PI3K/Akt axis

Excessive hepatic lipid accumulation is involved in the pathogenesis of nonalcoholic fatty liver disease (NAFLD). A previous study showed that the circular RNA (circRNA) PTK2 was significantly downregulated in NAFLD mice. However, the detailed function of circ PTK2 in NAFLD remains unclear. A high-fat diet (HFD) was used to establish a mouse model of NAFLD, and free fatty acid (FFA) treatment was used to establish an in vitro model of NAFLD. Oil red O staining was used to evaluate lipid accumulation. The pathological changes in mice were observed by HE staining. Western blotting and RT-qPCR were applied to assess protein and mRNA levels, respectively. A dual luciferase reporter assay and RIP were used to explore the relationship among circ PTK2, miR-200c and SIK2. Circ PTK2 and SIK2 were downregulated and miR-200c was upregulated in NAFLD. Upregulation of circ PTK2 reversed lipid accumulation in FFA-treated HepG2 cells. Moreover, circ PTK2 bound to miR-200c, and SIK2 was identified as the direct target of miR-200c. Moreover, the miR-200c inhibitor-induced decrease in lipid accumulation was reversed by SIK2 knockdown. Furthermore, the impact of circ PTK2 overexpression on PI3K/Akt signaling was partially reversed by SIK2 silencing. Circ PTK2 overexpression alleviates NAFLD development via the miR-200c/SIK2/PI3K/Akt axis. Thus, our work might provide new methods for NAFLD treatment.

pH-activated DNA nanomachine for miRNA-21 imaging to accurately identify cancer cell

MicroRNA (miRNA) imaging has been employed to distinguish cancer cells from normal cells by exploiting the overexpression of miRNA in cancer. Inspired by the acidic extracellular tumor microenvironment, we designed a pH-activated DNA nanomachine to enable the specific detection of cancer cells using miRNA imaging. The DNA nanomachine was engineered by assembling two hairpins (Y1 and Y2) onto the surface of a ZIF-8 metal-organic framework (MOF), which decomposed under acidic conditions to release the adsorbed DNA hairpin molecules in situ. The released hairpins were captured by the target miRNA-21 and underwent catalytic hairpin assembly amplification between Y1 and Y2. The detection limit for miRNA assays using the DNA nanomachine was determined to be 27 pM, which is low enough for sensitive detection in living cells. Living cell imaging of miRNA-21 further corroborated the application of the DNA nanomachine in the identification of cancer cell.

The control of compound inflorescences: insights from grasses and legumes

A major challenge in biology is to understand how organisms have increased developmental complexity during evolution. Inflorescences, with remarkable variation in branching systems, are a fitting model to understand architectural complexity. Inflorescences bear flowers that may become fruits and/or seeds, impacting crop productivity and species fitness. Great advances have been achieved in understanding the regulation of complex inflorescences, particularly in economically and ecologically important grasses and legumes. Surprisingly, a synthesis is still lacking regarding the common or distinct principles underlying the regulation of inflorescence complexity. Here, we synthesize the similarities and differences in the regulation of compound inflorescences in grasses and legumes, and propose that the emergence of novel higher-order repetitive modules is key to the evolution of inflorescence complexity.

AGL15 Promotion of Somatic Embryogenesis: Role and Molecular Mechanism

Plants have amazing regenerative properties with single somatic cells, or groups of cells able to give rise to fully formed plants. One means of regeneration is somatic embryogenesis, by which an embryonic structure is formed that "converts" into a plantlet. Somatic embryogenesis has been used as a model for zygotic processes that are buried within layers of maternal tissues. Understanding mechanisms of somatic embryo induction and development are important as a more accessible model for seed development. We rely on seed development not only for most of our caloric intake, but also as a delivery system for engineered crops to meet agricultural challenges. Regeneration of transformed cells is needed for this applied work as well as basic research to understand gene function. Here we focus on a MADS-domain transcription factor, AGAMOUS-Like15 (AGL15) that shows a positive correlation between accumulation levels and capacity for somatic embryogenesis. We relate AGL15 function to other transcription factors, hormones, and epigenetic modifiers involved in somatic embryo development.

Integration of small RNA, degradome, and transcriptome sequencing data illustrates the mechanism of low phosphorus adaptation in Camellia oleifera

Phosphorus (P) is an indispensable macronutrient for plant growth and development, and it is involved in various cellular biological activities in plants. Camellia oleifera is a unique high-quality woody oil plant that grows in the hills and mountains of southern China. However, the available P content is deficient in southern woodland soil. Until now, few studies focused on the regulatory functions of microRNAs (miRNAs) and their target genes under low inorganic phosphate (Pi) stress. In this study, we integrated small RNA, degradome, and transcriptome sequencing data to investigate the mechanism of low Pi adaptation in C. oleifera. We identified 40,689 unigenes and 386 miRNAs by the deep sequencing technology and divided the miRNAs into four different groups. We found 32 miRNAs which were differentially expressed under low Pi treatment. A total of 414 target genes of 108 miRNAs were verified by degradome sequencing. Gene ontology (GO) functional analysis of target genes found that they were related to the signal response to the stimulus and transporter activity, indicating that they may respond to low Pi stress. The integrated analysis revealed that 31 miRNA-target pairs had negatively correlated expression patterns. A co-expression regulatory network was established based on the profiles of differentially expressed genes. In total, three hub genes (ARF22, WRKY53, and SCL6), which were the targets of differentially expressed miRNAs, were discovered. Our results showed that integrated analyses of the small RNA, degradome, and transcriptome sequencing data provided a valuable basis for investigating low Pi in C. oleifera and offer new perspectives on the mechanism of low Pi tolerance in woody oil plants.