Integration of small RNAs and transcriptome sequencing uncovers a complex regulatory network during vernalization and heading stages of orchardgrass (Dactylis glomerata L.)

Transcriptome Profiling Provides Insights into the Early Development of Tiller Buds in High- and Low-Tillering Orchardgrass Genotypes

Orchardgrass (Dactylis glomerata L.) is among the most economically important perennial cool-season grasses, and is considered an excellent hay, pasture, and silage crop in temperate regions worldwide. Tillering is a vital feature that dominates orchardgrass regeneration and biomass yield. However, transcriptional dynamics underlying early-stage bud development in high- and low-tillering orchardgrass genotypes are unclear. Thus, this study assessed the photosynthetic parameters, the partially essential intermediate biomolecular substances, and the transcriptome to elaborate the early-stage profiles of tiller development. Photosynthetic efficiency and morphological development significantly differed between high- (AKZ-NRGR667) and low-tillering genotypes (D20170203) at the early stage after tiller formation. The 206.41 Gb of high-quality reads revealed stage-specific differentially expressed genes (DEGs), demonstrating that signal transduction and energy-related metabolism pathways, especially photosynthetic-related processes, influence tiller induction and development. Moreover, weighted correlation network analysis (WGCNA) and functional enrichment identified distinctively co-expressed gene clusters and four main regulatory pathways, including chlorophyll, lutein, nitrogen, and gibberellic acid (GA) metabolism pathways. Therefore, photosynthesis, carbohydrate synthesis, nitrogen efficient utilization, and phytohormone signaling pathways are closely and intrinsically linked at the transcriptional level. These findings enhance our understanding of tillering in orchardgrass and perennial grasses, providing a new breeding strategy for improving forage biomass yield.

Transcriptome Analysis Reveals the Potential Molecular Mechanisms of Tiller Bud Development in Orchardgrass

Tillering is a special type of branching and one of the important contributors to the yield of cereal crops. Strigolactone and sucrose play a vital role in controlling tiller formation, but their mechanism has not been elucidated completely in most crops. Orchardgrass (Dactylis glomerata L.) is an important perennial forage with prominent tillering ability among crops. To date, the mechanism of tillering in orchardgrass is still largely unknown. Therefore, we performed a transcriptome and miRNA analysis to reveal the potential RNA mechanism of tiller formation under strigolactone and sucrose treatment in orchardgrass. Our results found that D3, COL5, NCED1, HXK7, miRNA4393-z, and miRNA531-z could be key factors to control tiller bud development in orchardgrass. In addition, strigolactones might affect the ABA biosynthesis pathway to regulate the tiller bud development of orchardgrass, which may be related to the expression changes in miRNA4393-z, NCED1, and D10. miRNA531-z could be involved in the interaction of strigolactones and sucrose in regulating tillering. These results will be further used to clarify the potential mechanism of tillering for breeding new high-tillering and high-production orchardgrass varieties and beneficial to improving the production and reproduction of crops.

The cotton miR530-SAP6 module activated by systemic acquired resistance mediates plant defense against Verticillium dahliae

Many cotton miRNAs in root responding to Verticillium dahliae infection have been identified. Conversely, the miRNAs in leaf distantly responding to this fungal infection from roots via systemic acquired resistance (SAR) remain to be explored. Here, we constructed two groups of leaf sRNA libraries in cotton treated with V. dahliae via root-dipped method at 7- and 10-day post inoculation. Analysis of high-throughput sRNA sequencing identified 75 known and 379 novel miRNAs, of which 41 miRNAs significantly differentially expressed in fungal treatment plant leaves compared to the mock treatment at two time points. Then we characterized the cotton miR530-SAP6 module as a representative in the distant response to V. dahliae infection in roots. Based on degradome data and a luciferase (LUC) fusion reporter analysis, ghr-miR530 directedly cleaved GhSAP6 mRNA during the post-transcriptional process. Silencing of ghr-miR530 increased plant defense to this fungus, while its overexpression attenuated plant resistance. In link with ghr-miR530 function, the knockdown of GhSAP6 also decreased the plant resistance, resulting from down-regulation of SA-relative gene expression including GhNPR1 and GhPR1. In all, these results demonstrated that there are numerous miRNAs in leaf distantly responding to V. dahliae infection in roots mediate plant immunity.

Coordinating Diverse Functions of miRNA and lncRNA in Fleshy Fruit

Non-coding RNAs play vital roles in the diverse biological processes of plants, and they are becoming key topics in horticulture research. In particular, miRNAs and long non-coding RNAs (lncRNAs) are receiving increased attention in fruit crops. Recent studies in horticulture research provide both genetic and molecular evidence that miRNAs and lncRNAs regulate biological function and stress responses during fruit development. Here, we summarize multiple regulatory modules of miRNAs and lncRNAs and their biological roles in fruit sets and stress responses, which would guide the development of molecular breeding techniques on horticultural crops.

Comparative transcriptomic analysis of the gene expression and underlying molecular mechanism of submergence stress response in orchardgrass roots

INTRODUCTION

Submergence stress creates a hypoxic environment. Roots are the first plant organ to face these low-oxygen conditions, which causes damage and affects the plant growth and yield. Orchardgrass (Dactylis glomerata L.) is one of the most important cold-season forage grasses globally. However, their submergence stress-induced gene expression and the underlying molecular mechanisms of orchardgrass roots are still unknown.

METHODS

Using the submergence-tolerant 'Dianbei' and submergence-sensitive 'Anba', the transcriptomic analysis of orchardgrass roots at different time points of submergence stress (0 h, 8 h, and 24 h) was performed.

RESULTS

We obtained 118.82Gb clean data by RNA-Seq. As compared with the control, a total of 6663 and 9857 differentially expressed genes (DEGs) were detected in Dianbei, while 7894 and 11215 DEGs were detected in Anba at 8 h and 24 h post-submergence-stress, respectively. Gene Ontology (GO) enrichment analysis obtained 986 terms, while Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis obtained 123 pathways. Among them, the DEGs in plant hormones, mitogen-activated protein kinase (MAPK) and Ca²⁺ signal transduction were significantly differentially expressed in Dianbei, but not in Anba.

DISCUSSION

This study was the first to molecularly elucidate the submergence stress tolerance in the roots of two orchardgrass cultivars. These findings not only enhanced our understanding of the orchardgrass submergence tolerance, but also provided a theoretical basis 36 for the cultivation of submergence-tolerant forage varieties.

Large-scale analyses of angiosperm Flowering Locus T genes reveal duplication and functional divergence in monocots

FLOWERING LOCUS T (FT) are well-known key genes for initiating flowering in plants. Delineating the evolutionary history and functional diversity of FT genes is important for understanding the diversification of flowering time and how plants adapt to the changing surroundings. We performed a comprehensive phylogenetic analysis of FT genes in 47 sequenced flowering plants and the 1,000 Plant Transcriptomes (1KP) database with a focus on monocots, especially cereals. We revealed the evolutionary history of FT genes. The FT genes in monocots can be divided into three clades (I, II, and III), whereas only one monophyletic group was detected in early angiosperms, magnoliids, and eudicots. Multiple rounds of whole-genome duplications (WGD) events followed by gene retention contributed to the expansion and variation of FT genes in monocots. Amino acid sites in the clade II and III genes were preferentially under high positive selection, and some sites located in vital domain regions are known to change functions when mutated. Clade II and clade III genes exhibited high variability in important regions and functional divergence compared with clade I genes; thus, clade I is more conserved than clade II and III. Genes in clade I displayed higher expression levels in studied organs and tissues than the clade II and III genes. The co-expression modules showed that some of the FT genes might have experienced neofunctionalization and subfunctionalization, such as the acquisition of environmental resistance. Overall, FT genes in monocots might form three clades by the ancient gene duplication, and each clade was subsequently subjected to different selection pressures and amino acid substitutions, which eventually led to different expression patterns and functional diversification. Our study provides a global picture of FT genes' evolution in monocots, paving a road for investigating FT genes' function in future.

Integrating physiology, genetics, and transcriptome to decipher a new thermo-sensitive and light-sensitive virescent leaf gene mutant in cucumber

Chloroplasts are the material basis of photosynthesis, and temperature and light severely affect chloroplast development and thus influence photosynthetic efficiency. This study identified a spontaneous virescent leaf mutant, SC311Y, whose cotyledons and true leaves were yellow and gradually turned green. However, temperature and light affected the process of turning green. In addition, this mutant (except at the seedling stage) had ruffled leaves with white stripes, sterile males, and poorly fertile female flowers. Genetic characteristics analysis revealed that the recessive gene controlled the virescent leaf. Two F₂ populations mapped v-3 to the interval of 33.54-35.66 Mb on chromosome 3. In this interval, BSA-Seq, RNA-Seq, and cDNA sequence analyses revealed only one nonsynonymous mutation in the Csa3G042730 gene, which encoded the RNA exosome supercomplex subunit resurrection1 (RST1). Csa3G042730 was predicted to be the candidate gene controlling the virescent leaf, and the candidate gene may regulate chloroplast development by regulating plastid division2 (PDV2). A transcriptome analysis showed that different factors caused the reduced chlorophyll and carotenoid content in the mutants. To our knowledge, this study is the first report of map-based cloning related to virescent leaf, male-sterile, and chloroplast RNA regulation in cucumber. The results could accelerate the study of the RNA exosome supercomplex for the dynamic regulation of chloroplast RNA.

Genome-wide identification, phylogenetic analysis, and expression analysis of the SPL gene family in orchardgrass (Dactylis glomerata L.)

SPL (SQUAMOSA promoter binding protein-like) is a plant-specific transcription factor family that contains the conserved SBP domain, which plays a vital role in the vegetative-to-reproductive phase transition, flowering development and regulation, tillering/branching, and stress responses. Although the SPL family has been identified and characterized in various plant species, limited information about it has been obtained in orchardgrass, which is a critical forage crop worldwide. In this study, 17 putative DgSPL genes were identified among seven chromosomes, and seven groups that share similar gene structures and conserved motifs were determined by phylogenetic analysis. Of these, eight genes have potential target sites for miR156. cis-Element and gene ontology annotation analysis indicated DgSPLs may be involved in regulating development and abiotic stress responses. The expression patterns of eight DgSPL genes at five developmental stages, in five tissues, and under three stress conditions were determined by RNA-seq and qRT-PCR. These assays indicated DgSPLs are involved in vegetative-to-reproductive phase transition, floral development, and stress responses. The transient expression analysis in tobacco and heterologous expression assays in yeast indicated that miR156-targeted DG1G01828.1 and DG0G01071.1 are nucleus-localized proteins, that may respond to drought, salt, and heat stress. Our study represents the first systematic analysis of the SPL family in orchardgrass. This research provides a comprehensive assessment of the DgSPL family, which lays the foundation for further examination of the role of miR156/DgSPL in regulating development and stress responses in forages grasses.

Blocking miR530 Improves Rice Resistance, Yield, and Maturity

MicroRNAs fine-tune plant growth and resistance against multiple biotic and abiotic stresses. The trade-off between biomass and resistance can penalize crop yield. In this study, we have shown that rice miR530 regulates blast disease resistance, yield, and growth period. While the overexpression of miR530 results in compromised blast disease resistance, reduced grain yield, and late maturity, blocking miR530 using a target mimic (MIM530) leads to enhanced resistance, increased grain yield, and early maturity. Further study revealed that the accumulation of miR530 was decreased in both leaves and panicles along with the increase of age. Such expression patterns were accordant with the enhanced resistance from seedlings to adult plants, and the grain development from panicle formation to fully-filled seeds. Divergence analysis of miR530 precursor with upstream 1,000-bp promoter sequence in 11 rice species revealed that miR530 was diverse in Oryza sativa japonica and O. sativa indica group, which was consistent with the different accumulation of miR530 in japonica accessions and indica accessions. Altogether, our results indicate that miR530 coordinates rice resistance, yield, and maturity, thus providing a potential regulatory module for breeding programs aiming to improve yield and disease resistance.

Small RNA and Transcriptome Sequencing Reveals miRNA Regulation of Floral Thermogenesis in Nelumbo nucifera

The sacred lotus (Nelumbo nucifera Gaertn.) can produce heat autonomously and maintain a relatively stable floral chamber temperature for several days when blooming. Floral thermogenesis is critical for flower organ development and reproductive success. However, the regulatory role of microRNA (miRNA) underlying floral thermogenesis in N. nucifera remains unclear. To comprehensively understand the miRNA regulatory mechanism of thermogenesis, we performed small RNA sequencing and transcriptome sequencing on receptacles from five different developmental stages. In the present study, a total of 172 known miRNAs belonging to 39 miRNA families and 126 novel miRNAs were identified. Twenty-nine thermogenesis-related miRNAs and 3024 thermogenesis-related mRNAs were screened based on their expression patterns. Of those, seventeen differentially expressed miRNAs (DEMs) and 1765 differentially expressed genes (DEGs) had higher expression during thermogenic stages. The upregulated genes in the thermogenic stages were mainly associated with mitochondrial function, oxidoreductase activity, and the energy metabolism process. Further analysis showed that miR156_2, miR395a_5, miR481d, and miR319p may play an important role in heat-producing activity by regulating cellular respiration-related genes. This study provides comprehensive miRNA and mRNA expression profile of receptacle during thermogenesis in N. nucifera, which advances our understanding on the regulation of floral thermogenesis mediated by miRNA.

Genome-wide identification, structural analysis and expression profiles of GRAS gene family in orchardgrass

The GRAS gene family is a family of transcription factors that regulates plant growth and development. Despite being well-studied in many plant species, little is known about this gene family in orchardgrass (Dactylis glomerata L.), one of the top four economically important perennial forage grasses cultivated worldwide. We identified 46 GRAS genes in orchardgrass and analyzed their characteristics by phylogenetic, gene structural, motifs and expression patterns analysis. The phylogenetic analysis of eight species revealed that DgGRAS family had the evolutional conservation and closer homology relationship with the GRAS family of rice, barley and Brachypodium distachyon. Moreover, 46 DgGRAS proteins were divided into eight subfamilies based on the tree topology and rice or Arabidopsis classification, and LISCL subfamily was the largest one. Besides, we found that the motif 15 may be unique to the orchardgrass LISCL subfamily, and the motif 6 and motif 17 had indispensable functions in the orchardgrass LISCL subfamily. We further analyzed the expression profiles of DgGRAS genes at mature and seeding stage. And we found that DgGRAS17 played an important role in the growth and development no matter what stage it was at. DgGRAS5, DgGRAS28, DgGRAS31, DgGRAS42 and DgGRAS44 got involved in processes of the growth and development at seeding stage instead of mature stage. These results indicated that the major expression patterns and detailed functions of the DgGRAS genes varied with developmental stages. Taken together, this is the first systematic analysis of the GRAS gene family in the orchardgrass genome and the results provide insights into the potential functions of GRAS genes.