The Making of Leaves: How Small RNA Networks Modulate Leaf Development

The Emerging Role of Non-Coding RNAs (ncRNAs) in Plant Growth, Development, and Stress Response Signaling

Plant species utilize a variety of regulatory mechanisms to ensure sustainable productivity. Within this intricate framework, numerous non-coding RNAs (ncRNAs) play a crucial regulatory role in plant biology, surpassing the essential functions of RNA molecules as messengers, ribosomal, and transfer RNAs. ncRNAs represent an emerging class of regulators, operating directly in the form of small interfering RNAs (siRNAs), microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs). These ncRNAs exert control at various levels, including transcription, post-transcription, translation, and epigenetic. Furthermore, they interact with each other, contributing to a variety of biological processes and mechanisms associated with stress resilience. This review primarily concentrates on the recent advancements in plant ncRNAs, delineating their functions in growth and development across various organs such as root, leaf, seed/endosperm, and seed nutrient development. Additionally, this review broadens its scope by examining the role of ncRNAs in response to environmental stresses such as drought, salt, flood, heat, and cold in plants. This compilation offers updated information and insights to guide the characterization of the potential functions of ncRNAs in plant growth, development, and stress resilience in future research.

Non-coding RNAs and leaf senescence: Small molecules with important roles

Non-coding RNAs (ncRNAs) are a special class of functional RNA molecules that are not translated into proteins. ncRNAs have emerged as pivotal regulators of diverse developmental processes in plants. Recent investigations have revealed the association of ncRNAs with the regulation of leaf senescence, a complex and tightly regulated developmental process. However, a comprehensive review of the involvement of ncRNAs in the regulation of leaf senescence is still lacking. This manuscript aims to summarize the molecular mechanisms underlying ncRNAs-mediated leaf senescence and the potential applications of ncRNAs to manipulate the onset and progression of leaf senescence. Various classes of ncRNAs, including microRNAs (miRNAs), small interfering RNAs (siRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs), are discussed in terms of their regulatory mechanisms in leaf senescence. Furthermore, we explore the interactions between ncRNA and the key regulators of senescence, including transcription factors as well as core components in phytohormone signaling pathways. We also discuss the possible challenges and approaches related to ncRNA-mediated leaf senescence. This review contributes to a further understanding of the intricate regulatory network involving ncRNAs in leaf senescence.

Anatomical Mechanisms of Leaf Blade Morphogenesis in Sasaella kogasensis 'Aureostriatus'

There are limited studies on the cytology of bamboo leaf development from primordium to maturity. This study delves into the leaf morphological characteristics and growth patterns of Sasaella kogasensis 'Aureostriatus' and provides a three-dimensional anatomical analysis of cell division, expansion, and degradation. Leaves on the same branch develop bottom-up, while individual leaves develop the other way around. Like bamboo shoots and culms, the leaves follow a "slow-fast-slow" growth pattern, with longitudinal growth being predominant during their development. The growth zones of individual leaves included division, elongation, and maturation zones based on the distribution of growth space. By measuring 13,303 epidermal long cells and 3293 mesophyll cells in longitudinal sections of rapidly elongating leaves, we observed that in the rapid elongation phase (S4-S5), the division zone was located in the 1-2 cm segment at the bottom of the leaf blade and maintained a constant size, continuously providing new cells for leaf elongation, whereas in the late rapid elongation phase (S6), when the length of the leaf blade was approaching that of a mature leaf, its cells at the bottom of the blade no longer divided and were replaced by the ability to elongate. Furthermore, to gain an insight into the dynamic changes in the growth of the S. kogasensis 'Aureostriatus' leaves in the lateral and periclinal directions, the width and thickness of 1459 epidermal and 2719 mesophyll cells were counted in the mid-cross section of leaves at different developmental stages. The results showed that during the early stages of development (S1-S3), young leaves maintained vigorous division in the lateral direction, while periplasmic division gradually expanded from the bottom to the top of the leaf blade and the number of cell layers stabilized at S4. The meristematic tissues on both sides of the leaf were still able to divide at S4 but the frequency of the division gradually decreased, while cell division and expansion occurred simultaneously between the veins. At S6, the cells at the leaf margins and between the veins were completely differentiated and the width of the leaf blade no longer expanded. These findings revealed changes in cell growth anisotropically during the leaf development of S. kogasensis 'Aureostriatus' and demonstrated that leaf elongation was closely related to the longitudinal expansion of epidermal cells and proliferative growth of mesophyll cells, whereas the cell division of meristematic tissues and expansion of post-divisional cells contributed to the increases in blade width and thickness. The presented framework will facilitate a further exploration of the molecular regulatory mechanisms of leaf development in S. kogasensis 'Aureostriatus' and provide relevant information for developmental and taxonomic studies of bamboo plants.

Genome-wide identification of GROWTH-REGULATING FACTORs in Liriodendron chinense and functional characterization of LcGRF2 in leaf size regulation

GROWTH-REGULATING FACTORs (GRFs) play a pivotal role in the regulation of leaf size in plants and have been widely reported in plants. However, their specific functions in leaf size regulation in Liriodendron chinense remains unclear. Therefore, in this study, we identified GRF genes on a genome-wide scale in L. chinense to characterize the roles of LcGRFs in regulating leaf size. A total of nine LcGRF genes were identified, and these genes exhibited weak expression in mature leaves but strong expression in shoot apex. Notably, LcGRF2 exhibited the highest expression level in the shoot apex of L. chinense. Further RT-qPCR assay revealed that the expression level of LcGRF2 gradually decreased along with the leaf development process, and also displayed a gradient along the leaf proximo-distal and medio-lateral axes. Furthermore, overexpression of LcGRF2 in Arabidopsis thaliana resulted in increased leaf size, and significantly up-regulated the expression of genes involved in cell division like AtCYCD3;1, AtKNOLLE, and AtCYCB1;1, indicating that LcGRF2 may influence leaf size by promoting cell proliferation. This work contributes to a better understanding of the roles and molecular mechanisms of LcGRFs in the regulation of leaf size in L. chinense.

Coordination of floral and fiber development in cotton (Gossypium) by hormone- and flavonoid-signalling associated regulatory miRNAs

Various plant development activities and stress responses are tightly regulated by various microRNAs (miRNA) and their target genes, or transcription factors in a spatiotemporal manner. Here, to exemplify how flowering-associated regulatory miRNAs synchronize their expression dynamics during floral and fiber development in cotton, constitutive expression diminution transgenic lines of auxin-signaling regulatory Gh-miR167 (35S-MIM167) were developed through target mimicry approach. 'Moderate' (58% to 80%)- and 'high' (> 80%)-Gh-miR167 diminution mimic lines showed dosage-dependent developmental deformities in anther development, pollen maturation, and fruit (= boll) formation. Cross pollination of 'moderate' 35S-MIM167 mimic lines with wild type (WT) plant partially restored boll formation and emergence of fiber initials on the ovule surface. Gh-miR167 diminution favored organ-specific transcription biases in miR159, miR166 as well as miR160, miR164, and miR172 along with their target genes during anther and petal development, respectively. Similarly, accumulative effect of percent Gh-miR167 diminution, cross regulation of its target ARF6/8 genes, and temporal mis-expression of hormone signaling- and flavonoid biosynthesis-associated regulatory miRNAs at early fiber initiation stage caused irregular fiber formation. Spatial and temporal transcription proportions of regulatory miRNAs were also found crucial for the execution of hormone- and flavonoid-dependent progression of floral and fiber development. These observations discover how assorted regulatory genetic circuits get organized in response to Gh-miR167 diminution and converge upon ensuing episodes of floral and fiber development in cotton.

Independent Responses of Photosynthesis and Plant Morphology to Alterations of PIF Proteins and Light-Dependent MicroRNA Contents in Arabidopsis thaliana pif Mutants Grown under Lights of Different Spectral Compositions

The effects of the quality of light on the content of phytochrome interacting factors (PIFs) such as PIF3, PIF4 and PIF5, as well as the expression of various light-dependent microRNAs, in adult Arabidopsis thaliana pif mutant plants (pif4, pif5, pif3pif5, pif4pif5, pif3pif4pif5) were studied. We demonstrate that under blue light, the pif4 mutant had maximal expression of most of the studied microRNAs (miR163, miR319, miR398, miR408, miR833) when the PIF4 protein in plants was reduced. This finding indicates that the PIF4 protein is involved in the downregulation of this group of microRNAs. This assumption is additionally confirmed by the fact that under the RL spectrum in pif5 mutants, practically the same miRNAs decrease expression against the background of an increase in the amount of PIF4 protein. Unlike the WT and other mutants, the pif4 mutant responded to the BL spectrum not only by activating the expression of light-dependent miRNAs, but also by a significant increase in the expression of transcription factors and key light signalling genes. These molecular reactions do not affect the activity of photosynthesis but may be involved in the formation of a light quality-dependent phenotype.

Ecotype-specific blockage of tasiARF production by two different RNA viruses in Arabidopsis

Arabidopsis thaliana is one of the most studied model organisms of plant biology with hundreds of geographical variants called ecotypes. One might expect that this enormous genetic variety could result in differential response to pathogens. Indeed, we observed previously that the Bur ecotype develops much more severe symptoms (upward curling leaves and wavy leaf margins) upon infection with two positive-strand RNA viruses of different families (turnip vein-clearing virus, TVCV, and turnip mosaic virus, TuMV). To find the genes potentially responsible for the ecotype-specific response, we performed a differential expression analysis of the mRNA and sRNA pools of TVCV and TuMV-infected Bur and Col plants along with the corresponding mock controls. We focused on the genes and sRNAs that showed an induced or reduced expression selectively in the Bur virus samples in both virus series. We found that the two ecotypes respond to the viral infection differently, yet both viruses selectively block the production of the TAS3-derived small RNA specimen called tasiARF only in the virus-infected Bur plants. The tasiARF normally forms a gradient through the adaxial and abaxial parts of the leaf (being more abundant in the adaxial part) and post-transcriptionally regulates ARF4, a major leaf polarity determinant in plants. The lack of tasiARF-mediated silencing could lead to an ectopically expressed ARF4 in the adaxial part of the leaf where the misregulation of auxin-dependent signaling would result in an irregular growth of the leaf blade manifesting as upward curling leaf and wavy leaf margin. QTL mapping using Recombinant Inbred Lines (RILs) suggests that the observed symptoms are the result of a multigenic interaction that allows the symptoms to develop only in the Bur ecotype. The particular nature of genetic differences leading to the ecotype-specific symptoms remains obscure and needs further study.

Genetic dissection of morphological variation in rosette leaves and leafy heads in cabbage (Brassica oleracea var. capitata)

Correlations between morphological traits of cabbage rosette leaves and heads were found. Genome-wide association studies of these traits identified 50 robust quantitative trait loci in multiple years. Half of these loci affect both organs. Cabbage (Brassica oleracea var. capitata) is an economically important vegetable crop cultivated worldwide. Cabbage plants go through four vegetative stages: seedling, rosette, folding and heading. Rosette leaves are the largest leaves of cabbage plants and provide most of the energy needed to produce the leafy head. To understand the relationship and the genetic basis of leaf development and leafy head formation, 308 cabbage accessions were scored for rosette leaf and head traits in three-year field trials. Significant correlations were found between morphological traits of rosette leaves and heads, namely leaf area with the head area, height and width, and leaf width with the head area and head height, when heads were harvested at a fixed number of days after sowing. Fifty robust quantitative trait loci (QTLs) for rosette leaf and head traits distributed over all nine chromosomes were identified with genome-wide association studies. All these 50 loci were identified in multiple years and generally affect multiple traits. Twenty-five of the QTL were associated with both rosette leaf and leafy head traits. We discuss thirteen candidate genes identified in these QTL that are expressed in heading leaves, with an annotation related to auxin and other phytohormones, leaf development, and leaf polarity that likely play a role in leafy head development or rosette leaf expansion.

Parasitic plant small RNA analyses unveil parasite-specific signatures of microRNA retention, loss, and gain

Parasitism is a successful life strategy that has evolved independently in several families of vascular plants. The genera Cuscuta and Orobanche represent examples of the two profoundly different groups of parasites: one parasitizing host shoots and the other infecting host roots. In this study, we sequenced and described the overall repertoire of small RNAs from Cuscuta campestris and Orobanche aegyptiaca. We showed that C. campestris contains a number of novel microRNAs (miRNAs) in addition to a conspicuous retention of miRNAs that are typically lacking in other Solanales, while several typically conserved miRNAs seem to have become obsolete in the parasite. One new miRNA appears to be derived from a horizontal gene transfer event. The exploratory analysis of the miRNA population (exploratory due to the absence of a full genomic sequence for reference) from the root parasitic O. aegyptiaca also revealed a loss of a number of miRNAs compared to photosynthetic species from the same order. In summary, our study shows partly similar evolutionary signatures in the RNA silencing machinery in both parasites. Our data bear proof for the dynamism of this regulatory mechanism in parasitic plants.

Identification and analysis of HD-Zip genes involved in the leaf development of Liriodendron chinense using multidimensional analysis

Homeodomain-leucine zipper (HD-Zip) proteins are plant-specific transcription factors that play important roles in different biological processes, especially leaf development. However, no studies to date have identified the HD-Zip genes in Liriodendron chinense nor characterized their functions. We identified the HD-Zip genes in L. chinense by analysing the phylogeny, chromosome location, structure, conserved motif, cis-regulatory elements, synteny, post-transcriptional regulation and expression patterns of these genes during leaf development. A total of 36 LcHD-Zip genes were identified and divided into four subfamilies (HD-Zip I to IV). Synteny analysis revealed that segmental duplication was the main force driving the expansion of LcHD-Zip genes. These 36 LcHD-Zip genes exhibited 11 different expression patterns. Pattern 1, 2, 3, 4, 6, 7, 8 and 9 genes may play important roles in leaf development, such as leaf initiation, leaf polarity establishment, leaf shape development, phytohormone-mediated leaf growth and leaf epidermal structure formation. Four HD-Zip III genes were targeted by microRNAs (miRNAs), and the miR165/166a-HD-Zip regulatory module formed regulated leaf initiation and leaf polarity establishment. Overall, LcHD-Zip genes play key roles in leaf development of L. chinense. This work provides a foundation for the functional verification of HD-Zip genes identified in this study.

Identification of Crucial Amino Acids in Begomovirus C4 Proteins Involved in the Modulation of the Severity of Leaf Curling Symptoms

Begomoviruses frequently inflict upward or downward leaf curling symptoms on infected plants, leading to severe economic damages. Knowledge of the underlying mechanism controlling the leaf curling severity may facilitate the development of alternative disease management strategies. In this study, through genomic recombination between Ageratum yellow vein virus Nan-Tou strain (AYVV-NT) and Tomato leaf curl virus Tai-Chung Strain (TLCV-TC), which caused upward and downward leaf curling on Nicotiana benthamiana, respectively, it was found that the coding region of C4 protein might be involved in the determination of leaf curling directions. Sequence comparison and mutational analysis revealed that the cysteine and glycine at position 8 and 14 of AYVV-TC C4 protein, respectively, are involved in the modulation of leaf curling symptoms. Cross-protection assays further demonstrated that N. benthamiana inoculated with AYVV-carrying mutations of the aforementioned amino acids exhibited attenuated leaf curling symptoms under the challenge of wild-type AYVV-NT. Together, these findings revealed a new function of begomovirus C4 proteins involved in the modulation of leaf curling severity during symptom formation and suggested potential applications for managing viral diseases through manipulating the symptoms.

The Roles of microRNA-Long Non-coding RNA-mRNA Networks in the Regulation of Leaf and Flower Development in Liriodendron chinense

The leaf and the flower are vital plant organs owing to their roles in photosynthesis and reproduction. Long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and transcription factors (TFs) are very important to the development of these organs. Liriodendron chinense is a common ornamental tree species in southern China with an unusual leaf shape and tulip-like flowers. The genetic mechanisms underlying leaf and flower development in L. chinense and the miRNA-lncRNA-TF regulatory networks are poorly studied. Through the integration and analysis of different types of sequencing data, we identified the miRNA-lncRNA-TF regulatory networks that were related to leaf and flower development. These networks contained 105 miRNAs, 258 lncRNAs, 393 TFs, and 22 endogenous target mimics. Notably, lch-lnc7374-miR156h-SPL3 and lch-lnc7374-miR156j-SPL9 were potential regulators of stamen and pistil development in L. chinense, respectively. miRNA-lncRNA-mRNA regulatory networks were shown to impact anther development, male and female fertility, and petal color by regulating the biosynthesis of phenylpropanoid metabolites. Phenylpropanoid metabolite biosynthesis genes and TFs that were targeted by miRNAs and lncRNAs were differentially expressed in the leaf and flower. Moreover, RT-qPCR analysis confirmed 22 differentially expressed miRNAs, among which most of them showed obvious leaf or flower specificity; miR157a-SPL and miR160a-ARF module were verified by using RLM-RACE, and these two modules were related to leaf and flower development. These findings provide insight into the roles of miRNA-lncRNA-mRNA regulatory networks in organ development and function in L. chinense, and will facilitate further investigation into the regulatory mechanisms of leaf and flower development in L. chinense.

Transgenic Chrysanthemum indicum overexpressing cin-miR396a exhibits altered plant development and reduced salt and drought tolerance

The conserved microRNA396 (miR396) is involved in growth, development, and abiotic stress responses in a variety of plants by regulating target genes. Here, we obtained transgenic Chrysanthemum indicum (C. indicum) overexpressing the cin-miR396a gene. The transgenic plants (TGs) had longer internodes and fewer epidermal hairs in contrast with the wild-type (WT) control. cin-miR396a overexpression in C. indicum reduced salt tolerance and drought tolerance. After salt and drought stress compared with WT plants, the transgenic C. indicum exhibited a relative decrease in leaf water content, and the leaf free proline content, also exhibited a relative increase, in the leaf conductivity and leaf Malondialdehyde content, while the total chlorophyll content did not differ significantly from WT, and the Na⁺/K⁺ ratio in the roots of transgenic C. indicum increased after salt stress. We also identified two target genes of cin-miR396a, CiGRF1 and CiGRF5, whose expression was induced by salt and drought treatments and suppressed in transgenic C. indicum. Taken together, our results reveal a unique role for the regulatory module of miR396a-GRFs in C. indicum development and response to abiotic stresses. cin-miR396a plays a negative regulatory role in C. indicum in response to salt and drought stresses.

MicroRNA775 regulates intrinsic leaf size and reduces cell wall pectin levels by targeting a galactosyltransferase gene in Arabidopsis

Plants possess unique primary cell walls made of complex polysaccharides that play critical roles in determining intrinsic cell and organ size. How genes responsible for synthesizing and modifying the polysaccharides in the cell wall are regulated by microRNAs (miRNAs) to control plant size remains largely unexplored. Here we identified 23 putative cell wall-related miRNAs, termed as CW-miRNAs, in Arabidopsis thaliana and characterized miR775 as an example. We showed that miR775 post-transcriptionally silences GALT9, which encodes an endomembrane-located galactosyltransferase belonging to the glycosyltransferase 31 family. Over-expression of miR775 and deletion of GALT9 led to significantly enlarged leaf-related organs, primarily due to increased cell size. Monosaccharide quantification, confocal Raman imaging, and immunolabeling combined with atomic force microscopy revealed that the MIR775A-GALT9 circuit modulates pectin levels and the elastic modulus of the cell wall. We also showed that MIR775A is directly repressed by the transcription factor ELONGATED HYPOCOTYL5 (HY5). Genetic analysis confirmed that HY5 is a negative regulator of leaf size that acts through the HY5-MIR775A-GALT9 repression cascade to control pectin levels. These findings demonstrate that miR775-regulated cell wall remodeling is an integral determinant of intrinsic leaf size in A. thaliana. Studying other CW-miRNAs would provide more insights into cell wall biology.

Identification and molecular characterization of miRNAs and their target genes associated with seed development through small RNA sequencing in chickpea

Multiple studies have attempted to dissect the molecular mechanism underlying seed development in chickpea (Cicer arietinum L.). These studies highlight the need to focus on the role of miRNAs in regulating storage protein accumulation in seeds. Therefore, a total of 8,856,691 short-read sequences were generated from a small RNA library of developing chickpea seeds and were analyzed using miRDeep-P to identify 74 known and 26 novel miRNA sequences. Known miRNAs were classified into 22 miRNA families with miRNA156 family being most abundant. Of the 26 putative novel miRNAs identified, only 22 could be experimentally validated using stem loop end point PCR. Differential expression analyses led to the identification of known as well as novel miRNAs that could regulate various stages of chickpea seed development. In silico target prediction revealed several important target genes and transcription factors like SPL, mediator of RNA Polymerase II transcription subunit 12, aspartic proteinase and NACs, which were further validated by real-time PCR analysis. A comparative expression analysis in chickpea genotypes with contrasting seed protein content revealed one known (Car-miR156h) and two novel miRNA (Car-novmiR7 and Car-novmiR23) candidates to be highly expressed in the LPC (low protein content) chickpea genotypes, targets of which are known to regulate seed storage protein accumulation. Therefore, this study provides a useful resource in the form of miRNA and their targets which can be further utilized to understand and manipulate various regulatory mechanisms involved in seed development with the overall aim of improving yield and nutrition attributes in chickpea.

The Polycomb group methyltransferase StE(z)2 and deposition of H3K27me3 and H3K4me3 regulate the expression of tuberization genes in potato

Polycomb repressive complex (PRC) group proteins regulate various developmental processes in plants by repressing target genes via H3K27 trimethylation, and they function antagonistically with H3K4 trimethylation mediated by Trithorax group proteins. Tuberization in potato has been widely studied, but the role of histone modifications in this process is unknown. Recently, we showed that overexpression of StMSI1, a PRC2 member, alters the expression of tuberization genes in potato. As MSI1 lacks histone-modification activity, we hypothesized that this altered expression could be caused by another PRC2 member, StE(z)2, a potential H3K27 methyltransferase in potato. Here, we demonstrate that a short-day photoperiod influences StE(z)2 expression in the leaves and stolons. StE(z)2 overexpression alters plant architecture and reduces tuber yield, whereas its knockdown enhances yield. ChIP-sequencing using stolons induced by short-days indicated that several genes related to tuberization and phytohormones, such as StBEL5/11/29, StSWEET11B, StGA2OX1, and StPIN1 carry H3K4me3 or H3K27me3 marks and/or are StE(z)2 targets. Interestingly, we observed that another important tuberization gene, StSP6A, is targeted by StE(z)2 in leaves and that it has increased deposition of H3K27me3 under long-day (non-induced) conditions compared to short days. Overall, our results show that StE(z)2 and deposition of H3K27me3 and/or H3K4me3 marks might regulate the expression of key tuberization genes in potato.

Photocontrol of Axillary Bud Outgrowth by MicroRNAs: Current State-of-the-Art and Novel Perspectives Gained From the Rosebush Model

Shoot branching is highly dependent on environmental factors. While many species show some light dependence for branching, the rosebush shows a strict requirement for light to allow branching, making this species an excellent model to further understand how light impinges on branching. Here, in the first part, we provide a review of the current understanding of how light may modulate the complex regulatory network of endogenous factors like hormones (SL, IAA, CK, GA, and ABA), nutrients (sugar and nitrogen), and ROS to control branching. We review the regulatory contribution of microRNAs (miRNAs) to branching in different species, highlighting the action of such evolutionarily conserved factors. We underline some possible pathways by which light may modulate miRNA-dependent regulation of branching. In the second part, we exploit the strict light dependence of rosebush for branching to identify putative miRNAs that could contribute to the photocontrol of branching. For this, we first performed a profiling of the miRNAs expressed in early light-induced rosebush buds and next tested whether they were predicted to target recognized regulators of branching. Thus, we identified seven miRNAs (miR156, miR159, miR164, miR166, miR399, miR477, and miR8175) that could target nine genes (CKX1/6, EXPA3, MAX4, CYCD3;1, SUSY, 6PFK, APX1, and RBOHB1). Because these genes are affecting branching through different hormonal or metabolic pathways and because expression of some of these genes is photoregulated, our bioinformatic analysis suggests that miRNAs may trigger a rearrangement of the regulatory network to modulate branching in response to light environment.

Genome-Wide Screening and Characterization of Non-Coding RNAs in Coffea canephora

Coffea canephora grains are highly traded commodities worldwide. Non-coding RNAs (ncRNAs) are transcriptional products involved in genome regulation, environmental responses, and plant development. There is not an extensive genome-wide analysis that uncovers the ncRNA portion of the C. canephora genome. This study aimed to provide a curated characterization of six ncRNA classes in the Coffea canephora genome. For this purpose, we employed a combination of similarity-based and structural-based computational approaches with stringent curation. Candidate ncRNA loci had expression evidence analyzed using sRNA-seq libraries. We identified 7455 ncRNA loci (6976 with transcriptional evidence) in the C. canephora genome. This comprised of total 115 snRNAs, 1031 snoRNAs, 92 miRNA precursors, 602 tRNAs, 72 rRNAs, and 5064 lncRNAs. For miRNAs, we identified 159 putative high-confidence targets. This study was the most extensive genomic catalog of curated ncRNAs in the Coffea genus. This data might help elaborating more robust hypotheses in future comparative genomic studies as well as gene regulation and genome dynamics, helping to understand the molecular basis of domestication, environmental adaptation, resistance to pests and diseases, and coffee productivity.

Identification and Characterization of Stress-Responsive TAS3-Derived TasiRNAs in Melon

Small interfering RNAs (siRNA) are key regulators of gene expression that play essential roles in diverse biological processes. Trans-acting siRNAs (tasiRNAs) are a class of plant-endogenous siRNAs that lead the cleavage of nonidentical transcripts. TasiRNAs are usually involved in fine-tuning development. However, increasing evidence supports that tasiRNAs may be involved in stress response. Melon is a crop of great economic importance extensively cultivated in semiarid regions frequently exposed to changing environmental conditions that limit its productivity. However, knowledge of the precise role of siRNAs in general, and of tasiRNAs in particular, in regulating the response to adverse environmental conditions is limited. Here, we provide the first comprehensive analysis of computationally inferred melon-tasiRNAs responsive to two biotic (viroid-infection) and abiotic (cold treatment) stress conditions. We identify two TAS3-loci encoding to length (TAS3-L) and short (TAS3-S) transcripts. The TAS candidates predicted from small RNA-sequencing data were characterized according to their chromosome localization and expression pattern in response to stress. The functional activity of cmTAS genes was validated by transcript quantification and degradome assays of the tasiRNA precursors and their predicted targets. Finally, the functionality of a representative cmTAS3-derived tasiRNA (TAS3-S) was confirmed by transient assays showing the cleavage of ARF target transcripts.

Integrated Analysis of Small RNA, Transcriptome and Degradome Sequencing Provides New Insights into Floral Development and Abscission in Yellow Lupine (Lupinus luteus L.)

The floral development in an important legume crop yellow lupine (Lupinus luteus L., Taper cv.) is often affected by the abscission of flowers leading to significant economic losses. Small non-coding RNAs (sncRNAs), which have a proven effect on almost all developmental processes in other plants, might be of key players in a complex net of molecular interactions regulating flower development and abscission. This study represents the first comprehensive sncRNA identification and analysis of small RNA, transcriptome and degradome sequencing data in lupine flowers to elucidate their role in the regulation of lupine generative development. As shedding in lupine primarily concerns flowers formed at the upper part of the inflorescence, we analyzed samples from extreme parts of raceme separately and conducted an additional analysis of pedicels from abscising and non-abscising flowers where abscission zone forms. A total of 394 known and 28 novel miRNAs and 316 phased siRNAs were identified. In flowers at different stages of development 59 miRNAs displayed differential expression (DE) and 46 DE miRNAs were found while comparing the upper and lower flowers. Identified tasiR-ARFs were DE in developing flowers and were strongly expressed in flower pedicels. The DEmiR-targeted genes were preferentially enriched in the functional categories related to carbohydrate metabolism and plant hormone transduction pathways. This study not only contributes to the current understanding of how lupine flowers develop or undergo abscission but also holds potential for research aimed at crop improvement.

Roles of Small RNAs in Virus-Plant Interactions

Small RNAs (sRNAs), including microRNAs (miRNAs) and short interfering RNAs (siRNAs), are non-coding but powerful RNA molecules of 20–30 nucleotides in length. sRNAs play crucial regulatory roles in diverse plant biological processes. Recently, many studies on sRNAs have been reported. We summarize new findings of sRNAs in virus-plant interactions to accelerate the function analysis of sRNAs. The main content of this review article includes three parts: virus-responsive sRNAs, function analysis of sRNAs in virus pathogenicity or host resistance, and some sRNAs-mediated underlying mechanisms in virus-plant interactions. New findings of sRNAs deepen our understanding about sRNAs’ roles, which might contribute to the design of novel control measures against plant viruses.

Overexpression of OsAGO1b Induces Adaxially Rolled Leaves by Affecting Leaf Abaxial Sclerenchymatous Cell Development in Rice

BACKGROUND

ARGONAUTE 1 (AGO1) proteins can recruit small RNAs to regulate gene expression, involving several growth and development processes in Arabidopsis. Rice genome contains four AGO1 genes, OsAGO1a to OsAGO1d. However, the regulatory functions to rice growth and development of each AGO1 gene are still unknown.

RESULTS

We obtained overexpression and RNAi transgenic lines of each OsAGO1 gene. However, only up- and down-regulation of OsAGO1b caused multiple abnormal phenotypic changes in rice, indicating that OsAGO1b is the key player in rice growth and organ development compared with other three OsAGO1s. qRT-PCR assays showed that OsAGO1b was almost unanimously expressed in leaves at different developmental stages, and strongly expressed in spikelets at S1 to S3 stages. OsAGO1b is a typical AGO protein, and co-localized in both the nucleus and cytoplasm simultaneously. Overexpression of OsAGO1b caused adaxially rolled leaves and a series of abnormal phenotypes, such as the reduced tiller number and plant height. Knockdown lines of OsAGO1b showed almost normal leaves, but the seed setting percentage was significantly reduced accompanied by the disturbed anther patterning and reduced pollen fertility. Further anatomical observation revealed that OsAGO1b overexpression plants showed the partially defective development of sclerenchymatous cells on the abaxial side of leaves. In situ hybridization showed OsAGO1b mRNA was uniformly accumulated in P1 to P3 primordia without polarity property, suggesting OsAGO1b did not regulate the adaxial-abaxial polarity development directly. The expression levels of several genes related to leaf polarity development and vascular bundle differentiation were observably changed. Notably, the accumulation of miR166 and TAS3-siRNA was decreased, and their targeted OSHBs and OsARFs were significantly up-regulated. The mRNA distribution patterns of OSHB3 and OsARF4 in leaves remained almost unchanged between ZH11 and OsAGO1b overexpression lines, but their expression levels were enhanced at the regions of vascular bundles and sclerenchymatous cell differentiation.

CONCLUSIONS

In summary, we demonstrated OsAGO1b is the leading player among four OsAGO1s in rice growth and development. We propose that OsAGO1b may regulate the abaxial sclerenchymatous cell differentiation by affecting the expression of OSHBs in rice.

The interaction between miR160 and miR165/166 in the control of leaf development and drought tolerance in Arabidopsis

MicroRNAs (miRNAs) are a class of non-coding RNAs that play important roles in plant development and abiotic stresses. To date, studies have mainly focused on the roles of individual miRNAs, however, a few have addressed the interactions among multiple miRNAs. In this study, we investigated the interplay and regulatory circuit between miR160 and miR165/166 and its effect on leaf development and drought tolerance in Arabidopsis using Short Tandem Target Mimic (STTM). By crossing STTM160 Arabidopsis with STTM165/166, we successfully generated a double mutant of miR160 and miR165/166. The double mutant plants exhibited a series of compromised phenotypes in leaf development and drought tolerance in comparison to phenotypic alterations in the single STTM lines. RNA-seq and qRT-PCR analyses suggested that the expression levels of auxin and ABA signaling genes in the STTM-directed double mutant were compromised compared to the two single mutants. Our results also suggested that miR160-directed regulation of auxin response factors (ARFs) contribute to leaf development via auxin signaling genes, whereas miR165/166- mediated HD-ZIP IIIs regulation confers drought tolerance through ABA signaling. Our studies further indicated that ARFs and HD-ZIP IIIs may play opposite roles in the regulation of leaf development and drought tolerance that can be further applied to other crops for agronomic traits improvement.

A Resource for Inactivation of MicroRNAs Using Short Tandem Target Mimic Technology in Model and Crop Plants

microRNAs (miRNAs) are endogenous small non-coding RNAs that bind to mRNAs and target them for cleavage and/or translational repression, leading to gene silencing. We previously developed short tandem target mimic (STTM) technology to deactivate endogenous miRNAs in Arabidopsis. Here, we created hundreds of STTMs that target both conserved and species-specific miRNAs in Arabidopsis, tomato, rice, and maize, providing a resource for the functional interrogation of miRNAs. We not only revealed the functions of several miRNAs in plant development, but also demonstrated that tissue-specific inactivation of a few miRNAs in rice leads to an increase in grain size without adversely affecting overall plant growth and development. RNA-seq and small RNA-seq analyses of STTM156/157 and STTM165/166 transgenic plants revealed the roles of these miRNAs in plant hormone biosynthesis and activation, secondary metabolism, and ion-channel activity-associated electrophysiology, demonstrating that STTM technology is an effective approach for studying miRNA functions. To facilitate the study and application of STTM transgenic plants and to provide a useful platform for storing and sharing of information about miRNA-regulated gene networks, we have established an online Genome Browser (https://blossom.ffr.mtu.edu/designindex2.php) to display the transcriptomic and miRNAomic changes in STTM-induced miRNA knockdown plants.

Mechanisms Underlying the Environmentally Induced Plasticity of Leaf Morphology

The primary function of leaves is to provide an interface between plants and their environment for gas exchange, light exposure and thermoregulation. Leaves have, therefore a central contribution to plant fitness by allowing an efficient absorption of sunlight energy through photosynthesis to ensure an optimal growth. Their final geometry will result from a balance between the need to maximize energy uptake while minimizing the damage caused by environmental stresses. This intimate relationship between leaf and its surroundings has led to an enormous diversification in leaf forms. Leaf shape varies between species, populations, individuals or even within identical genotypes when those are subjected to different environmental conditions. For instance, the extent of leaf margin dissection has, for long, been found to inversely correlate with the mean annual temperature, such that Paleobotanists have used models based on leaf shape to predict the paleoclimate from fossil flora. Leaf growth is not only dependent on temperature but is also regulated by many other environmental factors such as light quality and intensity or ambient humidity. This raises the question of how the different signals can be integrated at the molecular level and converted into clear developmental decisions. Several recent studies have started to shed the light on the molecular mechanisms that connect the environmental sensing with organ-growth and patterning. In this review, we discuss the current knowledge on the influence of different environmental signals on leaf size and shape, their integration as well as their importance for plant adaptation.