Rice osa-miR171c Mediates Phase Change from Vegetative to Reproductive Development and Shoot Apical Meristem Maintenance by Repressing Four OsHAM Transcription Factors

Functions and Regulation of HAM Family Genes in Meristems During Gametophyte and Sporophyte Generations

A fascinating feature of land plants is their ability to continually initiate new tissues and organs throughout their lifespan, driven by a pool of pluripotent stem cells located in meristems. In seed plants, various types of meristems are initiated and maintained during the sporophyte generation, while their gametophytes lack meristems and rely on sporophyte tissues for growth. In contrast, seed-free vascular plants, such as ferns, develop meristems during both the sporophyte and gametophyte generations, allowing for the independent growth of both generations. Recent findings have highlighted both conserved and lineage-specific roles of the HAIRY MERISTEM (HAM) family of GRAS-domain transcriptional regulators in various meristems throughout the land plant lifecycle. Here, we review and discuss how HAM genes maintain meristem indeterminacy in both sporophytes and gametophytes, with a focus on studies performed in two model species: the flowering plant Arabidopsis thaliana and the fern Ceratopteris richardii. Additionally, we summarize the crucial and tightly regulated functions of the microRNA171 (miR171)-HAM regulatory modules, which define HAM spatial patterns and activities during meristem development across various meristem identities in land plants.

The role of Ancestral MicroRNAs in grass inflorescence development

Plant inflorescences are complex, highly diverse structures whose morphology is determined in meristems that form during reproductive development. Inflorescence structure influences flower formation, and consequently grain number, and yield in crops. Correct inflorescence and flower development require tight control of gene expression via complex interplay between regulatory networks. MicroRNAs (miRNAs) have emerged as fundamental modulators of gene expression at the transcriptional and/or post-transcriptional level in plant inflorescence development. First discovered more than three decades ago, miRNAs have proved to be revolutionary in advancing our mechanistic understanding of gene expression. This review highlights current knowledge of downstream target genes and pathways of some highly conserved miRNAs that regulate the maintenance, identity, and activity of inflorescence and floral meristems in economically and agriculturally important grass species, including rice (Oryza sativa), maize (Zea mays), barley (Hordeum vulgare), and wheat (Triticum aestivum). Furthermore, we summarize emerging regulatory networks of miRNAs and their targets to suggest new avenues and strategies for application of miRNAs as a tool to enhance crop yield and performance.

Genome-Wide Identification of the Rehmannia glutinosa miRNA Family and Exploration of Their Expression Characteristics Caused by the Replant Disease Formation-Related Principal Factor

Background/Objectives: Rehmannia glutinosa, a highly valuable medicinal plant in China, is encountering severe replant disease. Replant disease represents a complex stress driven by multiple principal factors (RDFs), including allelochemicals, microbes, and their interactions. miRNAs are recognized as key regulators of plant response to stresses; however, their specific roles within RDFs are not entirely clear. Methods: This study builds six RDF treatments, comprising R. glutinosa continuously planted (SP), normally planted (NP), and NP treated with ferulic acid (FA), Fusarium oxysporum (FO), and a combination of FA with FO (FAFO). sRNA-seq technology was used to identify crucial miRNAs in response to diverse RDFs. Results: In total, 30 sRNA datasets were generated from the SP, NP, FA, FO, and FAFO samples. A total of 160 known and 41 novel miRNAs (RgmiRNAs) were identified in the R. glutinosa genome based on the sRNA database. Abundance analysis revealed that RgmiRNAs in SP exhibited a distinct expression profile in comparison with others. Of these, 124, 86, 86, and 90 RgmiRNAs were differentially expressed in SP, FA, FO, and FAFO compared with NP. Target analysis indicated that RgmiRNAs downregulated in both SP and RDFs impede the organism growth of R. glutinosa. RgmiRNAs upregulated in SP can disrupt root formation and nutrient metabolism, in which, two RgmiR398 were uniquely expressed in SP. It was confirmed to target RgCSD genes. The expression patterns of RgmiR398 and RgCSD indicated that replant disease induces the oxidative damage of R. glutinosa through RgmiR398. Conclusions:RgmiRNA profiling under RDFs provides a theoretical basis for the further clarification of RgmiRNA function in replant disease.

Unraveling the Molecular Mechanisms by Which the miR171b-SCL6 Module Regulates Maturation in Lilium

Lilium is one of the most widely cultivated ornamental bulbous plants in the world. Although research has shown that variable temperature treatments can accelerate the development process from vegetative to reproductive growth in Lilium, the molecular regulation mechanisms of this development are not clear. In this study, Lbr-miR171b and its target gene, LbrSCL6, were selected and validated using transgenic functional verification, subcellular localization, and transcriptional activation. This study also investigated the differential expression of Lbr-miR171b and LbrSCL6 in two temperature treatment groups (25 °C and 15 °C). Lbr-miR171b expression significantly increased after the temperature change, whereas that of LbrSCL6 exhibited the opposite trend. Through in situ hybridization experiments facilitated by the design of hybridization probes targeting LbrSCL6, a reduction in LbrSCL6 expression was detected following variable temperature treatment at 15 °C. The transgenic overexpression of Lbr-miR171b in plants promoted the phase transition, while LbrSCL6 overexpression induced a delay in the phase transition. In addition, LbrWOX4 interacted with LbrSCL6 in yeast two-hybrid and bimolecular fluorescence complementation assays. In conclusion, these results explain the molecular regulatory mechanisms governing the phase transition in Lilium.

GhHAM regulates GoPGF-dependent gland development and contributes to broad-spectrum pest resistance in cotton

Cotton is a globally cultivated crop, producing 87% of the natural fiber used in the global textile industry. The pigment glands, unique to cotton and its relatives, serve as a defense structure against pests and pathogens. However, the molecular mechanism underlying gland formation and the specific role of pigment glands in cotton's pest defense are still not well understood. In this study, we cloned a gland-related transcription factor GhHAM and generated the GhHAM knockout mutant using CRISPR/Cas9. Phenotypic observations, transcriptome analysis, and promoter-binding experiments revealed that GhHAM binds to the promoter of GoPGF, regulating pigment gland formation in cotton's multiple organs via the GoPGF-GhJUB1 module. The knockout of GhHAM significantly reduced gossypol production and increased cotton's susceptibility to pests in the field. Feeding assays demonstrated that more than 80% of the cotton bollworm larvae preferred ghham over the wild type. Furthermore, the ghham mutants displayed shorter cell length and decreased gibberellins (GA) production in the stem. Exogenous application of GA3 restored stem cell elongation but not gland formation, thereby indicating that GhHAM controls gland morphogenesis independently of GA. Our study sheds light on the functional differentiation of HAM proteins among plant species, highlights the significant role of pigment glands in influencing pest feeding preference, and provides a theoretical basis for breeding pest-resistant cotton varieties to address the challenges posed by frequent outbreaks of pests.

Overexpression of Larch SCL6 Inhibits Transitions from Vegetative Meristem to Inflorescence and Flower Meristem in Arabidopsis thaliana (L.) Heynh

SCARECROW-LIKE6 (SCL6) plays a role in the formation and maintenance of the meristem. In Larix kaempferi (Lamb.) Carr., an important afforestation tree species in China, SCL6 (LaSCL6) has two alternative splicing variants-LaSCL6-var1 and LaSCL6-var2-which are regulated by microRNA171. However, their roles are still unclear. In this study, LaSCL6-var1 and LaSCL6-var2 were transformed into the Arabidopsis thaliana (L.) Heynh. genome, and the phenotypic characteristics of transgenic A. thaliana, including the germination percentage, root length, bolting time, flower and silique formation times, inflorescence axis length, and branch and silique numbers, were analyzed to reveal their functions. It was found that LaSCL6-var1 and LaSCL6-var2 overexpression shortened the root length by 41% and 31%, respectively, and increased the inflorescence axis length. Compared with the wild type, the bolting time in transgenic plants was delayed by approximately 2-3 days, the first flower and silique formation times were delayed by approximately 3-4 days, and the last flower and silique formation times were delayed by about 5 days. Overall, the life cycle in transgenic plants was prolonged by approximately 5 days. These results show that LaSCL6 overexpression inhibited the transitions from the vegetative meristem to inflorescence meristem and from the flower meristem to meristem arrest in A. thaliana, revealing the roles of LaSCL6-var1 and LaSCL6-var2 in the fate transition and maintenance of the meristem.

Network of GRAS transcription factors in plant development, fruit ripening and stress responses

The plant-specific family of GRAS transcription factors has been wide implicated in the regulation of transcriptional reprogramming associated with a diversity of biological functions ranging from plant development processes to stress responses. Functional analyses of GRAS transcription factors supported by in silico structural and comparative analyses are emerging and clarifying the regulatory networks associated with their biological roles. In this review, a detailed analysis of GRAS proteins' structure and biochemical features as revealed by recent discoveries indicated how these characteristics may impact subcellular location, molecular mechanisms, and function. Nomenclature issues associated with GRAS classification into different subfamilies in diverse plant species even in the presence of robust genomic resources are discussed, in particular how it affects assumptions of biological function. Insights into the mechanisms driving evolution of this gene family and how genetic and epigenetic regulation of GRAS contributes to subfunctionalization are provided. Finally, this review debates challenges and future perspectives on the application of this complex but promising gene family for crop improvement to cope with challenges of environmental transition.

Negative regulators of grain yield and mineral contents in rice: potential targets for CRISPR-Cas9-mediated genome editing

Rice is a major global staple food crop, and improving its grain yield and nutritional quality has been a major thrust research area since last decades. Yield and nutritional quality are complex traits which are controlled by multiple signaling pathways. Sincere efforts during past decades of research have identified several key genetic and molecular regulators that governed these complex traits. The advent of clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein 9 (Cas9)-mediated gene knockout approaches has accelerated the development of improved varieties; however, finding out target gene with negative regulatory function in particular trait without giving any pleiotropic effect remains a challenge. Here, we have reviewed past and recent literature and identified important negative regulators of grain yield and mineral contents which could be potential targets for CRISPR-Cas9-mediated gene knockout. Additionally, we have also compiled a list of microRNAs (miRNAs), which target positive regulators of grain yield, plant stress tolerance, and grain mineral contents. Knocking out these miRNAs could help to increase expression of such positive regulators and thus improve the plant trait. The knowledge presented in this review would help to further accelerate the CRISPR-Cas9-mediated trait improvement in rice.

Genome-wide identification and expression analysis of the GRAS family under low-temperature stress in bananas

Introduction

GRAS, named after GAI, RGA, and SCR, is a class of plant-specific transcription factors family that plays a crucial role in growth and development, signal transduction, and various stress responses.

Methods

To understand the biological functions of the banana GRAS gene family, a genome-wide identification and bioinformatics analysis of the banana GRAS gene family was performed based on information from the M. acuminata, M. balbisiana, and M. itinerans genomic databases.

Result

In the present study, we identified 73 MaGRAS, 59 MbGRAS, and 58 MiGRAS genes in bananas at the whole-genome scale, and 56 homologous genes were identified in the three banana genomes. Banana GRASs can be classified into 10 subfamilies, and their gene structures revealed that most banana GRAS gDNAs lack introns. The promoter sequences of GRASs had a large number of cis-acting elements related to plant growth and development, phytohormone, and adversity stress responsiveness. The expression pattern of seven key members of MaGRAS response to low-temperature stress and different tissues was also examined by quantitative reverse transcription polymerase chain reaction (qRT-PCR). The microRNAs-MaGRASs target prediction showed perfect complementarity of seven GRAS genes with the five mac-miRNAs. The expression of all seven genes was lowest in roots, and the expression of five genes was highest in leaves during low-temperature stress. The expression of MaSCL27-2, MaSCL27-3, and MaSCL6-1 was significantly lower under low-temperature stress compared to the control, except for MaSCL27-2, which was slightly higher than the 28°C control at 4 h. The expression of MaSCL27-2, MaSCL27-3, and MaSCL6-1 dropped to the lowest levels at 24 h, 12 h, and 4 h, respectively. The MaSCL27-4 and MaSCL6-2 expression was intermittently upregulated, rising to the highest expression at 24h, while the expression of MaSCL22 was less variable, remaining at the control level with small changes.

Discussion

In summary, it is tentatively hypothesized that the GRAS family has an important function in low-temperature stress in bananas. This study provides a theoretical basis for further analyzing the function of the banana GRAS gene and the resistance of bananas to cold temperatures.

Role of microRNA miR171 in plant development

MicroRNAs (miRNAs) are endogenous non-coding small RNA with 19-24 nucleotides (nts) in length, which play an essential role in regulating gene expression at the post-transcriptional level. As one of the first miRNAs found in plants, miR171 is a typical class of conserved miRNAs. The miR171 sequences among different species are highly similar, and the vast majority of them have both "GAGCCG" and "CAAUAU" fragments. In addition to being involved in plant growth and development, hormone signaling and stress response, miR171 also plays multiple and important roles in plants through interactions with microbe and other small-RNAs. The miRNA functions by regulating the expression of target genes. Most of miR171's target genes are in the GRAS gene family, but also include some NSP, miRNAs, lncRNAs, and other genes. This review is intended to summarize recent updates on miR171 regarding its function in plant life and hopefully provide new ideas for understanding miR171 function and regulatory mechanisms.

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.

Comparative protein analysis of two maize genotypes with contrasting tolerance to low temperature

BACKGROUND

Low temperature (LT) stress is one of the major environmental stress factors affecting the growth and yield of maize (Zea mays L.). Hence, it is important to unravel the molecular mechanisms behind LT stress tolerance to improve molecular breeding in LT tolerant genotypes. In the present study, two maize genotypes viz. Gurez local from Kashmir Himalaya and tropical grown GM6, were dissected for their LT stress response in terms of accumulation of differentially regulated proteins (DRPs). Leaf proteome analysis at three-leaf stage of maize seedlings subjected to LT stress of 6 °C for a total of 12 h duration was performed using two dimensional gel electrophoresis (2D-PAGE) followed by subsequent identification of the proteins involved.

RESULTS

After MALDI-TOF (Matrix-assisted laser desorption/ionization-time of flight) and bioinformatics analysis, 19 proteins were successfully identified in Gurez local, while as 10 proteins were found to get successful identification in GM6. The interesting observations from the present investigation is the identification of three novel proteins viz. threonine dehydratase biosynthetic chloroplastic, thylakoidal processing peptidase 1 chloroplastic, and nodulin-like protein, whose role in abiotic stress tolerance, in general, and LT stress, in particular, has not been reported so far. It is important to highlight here that most of LT responsive proteins including the three novel proteins were identified from Gurez local only, owing to its exceptional LT tolerance. From the protein profiles, obtained in both genotypes immediately after LT stress perception, it was inferred that stress responsive protein accumulation and their expression fashion help the Gurez local in seedling establishment and withstand unfavorable conditions as compared to GM6. This was inferred from the findings of pathway enrichment analysis like regulation of seed growth, timing of floral transition, lipid glycosylation, and aspartate family amino acid catabolic processes, besides other key stress defense mechanisms. However, in GM6, metabolic pathways enriched were found to be involved in more general processes including cell cycle DNA replication and regulation of phenylpropanoid metabolism. Furthermore, majority of the qRT-PCR results of the selected proteins demonstrated positive correlation between protein levels and transcript abundance, thereby strengthening our findings.

CONCLUSIONS

In conclusion, our findings reported majority of the identified proteins in Gurez local exhibiting up-regulated pattern under LT stress as compared to GM6. Furthermore, three novel proteins induced by LT stress were found in Gurez local, requiring further functional validation. Therefore, our results offer more insights for elucidating the molecular networks mediating LT stress tolerance in maize.

miR171 modulates induction of somatic embryogenesis in citrus callus

Overexpression of miR171 restored SE competence in the recalcitrant citrus callus, and inhibition of miR171 function weakened SE competence in the strongly embryogenic citrus callus. Somatic embryogenesis (SE) is an important way of in vitro regeneration for plants. For perennial woody crops such as citrus, embryogenic callus is usually induced from unfertilized aborted ovules and widely used in biotechnology aided breeding. However, SE capacity always declines in callus during subculture, which makes regeneration difficult and hinders the application of biotechnology. We previously found that miR171 may be a regulator of SE in citrus, based on the abundant expression of csi-miR171c in the embryogenic callus and during SE of citrus. Here, we report that miR171 promotes SE and is required for SE in citrus. Overexpression of miR171 restored SE competence in the recalcitrant callus of 'Guoqing No.1' Satsuma mandarin (G1), whereas inhibition of miR171 function by Short Tandem Target Mimic (STTM) weakened SE competence in the strongly embryogenic callus of 'Valencia' sweet orange (V). The comparative transcriptomic analysis in miR171 overexpressed callus line (OE) and the wild type callus (WT) indicated that overexpression of miR171 decreased the expression level of its SCARECROW-LIKE (CsSCL) targets, and activated stress response related biological processes and metabolic processes that are required for cell differentiation. However, CsSCLs were up-regulated in the OE callus during SE induction process, which activated the cell division and developmental processes that are required for embryogenesis progress. Our results validate the function of miR171 in regulation of SE and reveal the biological responses provoked by miR171 in citrus that may promote SE.

Function and Regulation of microRNA171 in Plant Stem Cell Homeostasis and Developmental Programing

MicroRNA171 (miR171), a group of 21-nucleotide single-strand small RNAs, is one ancient and conserved microRNA family in land plants. This review focuses on the recent progress in understanding the role of miR171 in plant stem cell homeostasis and developmental patterning, and the regulation of miR171 by developmental cues and environmental signals. Specifically, miR171 regulates shoot meristem activity and phase transition through repressing the HAIRYMERISTEM (HAM) family genes. In the model species Arabidopsis, miR171 serves as a short-range mobile signal, which initiates in the epidermal layer of shoot meristems and moves downwards within a limited distance, to pattern the apical-basal polarity of gene expression and drive stem cell dynamics. miR171 levels are regulated by light and various abiotic stresses, suggesting miR171 may serve as a linkage between environmental factors and cell fate decisions. Furthermore, miR171 family members also demonstrate both conserved and lineage-specific functions in land plants, which are summarized and discussed here.

Functional analysis of the eTM-miR171-SCL6 module regulating somatic embryogenesis in Lilium pumilum DC. Fisch

Somatic embryogenesis (SE) is of great significance in Lilium bulb production, germplasm preservation and genetic improvement. miRNAs are important regulators of plant growth and development at the transcriptional level. Previous research by our group has shown that lpu-miR171 and its target gene SCARECROW-LIKE 6 (SCL6) play an important regulatory role in lily SE, and we predicted and identified that endogenous target mimics (eTMs) can regulate lpu-miR171. However, the associated mechanism and internal regulatory network are not yet clear. In the present study, lpu-miR171 was used as an entry point to explore the regulatory network between its upstream eTMs and its downstream target gene LpSCL6, as well as to identify the mechanism of this regulatory network in Lilium SE. Tobacco transient transformation confirmed that miRNA171 significantly inhibited the expression of LpSCL6. On this basis, the Lilium stable genetic transformation system was used to demonstrate that silencing lpu-miR171a and lpu-miR171b and overexpressing LpSCL6-II and LpSCL6-I promoted starch accumulation in calli and the expression of key cell cycle genes, thus providing energy to meet preconditions for SE and accelerate the formation and development of Lilium somatic embryos. LpSCL6-II and LpSCL6-I are nuclear proteins with self-activation activity in yeast cells. In addition, we confirmed in Lilium that lpu-eTM171 is the eTM of lpu-miR171 that binds lpu-miR171 to prevent cleavage of the target gene LpSCL6, thereby promoting SE. Therefore, the present study established a new mechanism whereby the eTM-miR171-SCL6 module regulates SE in Lilium pumilum DC. Fisch. and provided new insights clarifying the mechanism of SE.

Rice microRNA171f/SCL6 module enhances drought tolerance by regulation of flavonoid biosynthesis genes

Plants have evolved sophisticated defense systems to enhance drought tolerance. These include the microRNA (miRNA) group of small noncoding RNAs that act as post-transcriptional regulators; however, details of the mechanisms by which they confer drought tolerance are not well understood. Here, we show that osa-MIR171f, a member of osa-MIR171 gene family, is mainly expressed in response to drought stress and regulates the transcript levels of SCARECROW-LIKE6-I (SCL6-I) and SCL6-II in rice (Oryza sativa). The SCL6 genes are known to be involved in shoot branching and flag leaf morphology. Osa-MIR171f-overexpressing (osa-MIR171f-OE) transgenic plants showed reduced drought symptoms compared with non-transgenic (NT) control plants under both field drought and polyethylene glycol (PEG)-mediated dehydration stress conditions. Transcriptome analysis of osa-MIR171f-OE plants and osa-mir171f-knockout (K/O) lines generated by clustered regularly interspaced short palindromic repeats (CRISPR/Cas9) revealed that osa-mature-miR171a-f (osa-miR171) regulates the expression of flavonoid biosynthesis genes, consequently leading to drought tolerance. This upregulation in the osa-MIR171f-OE plants, which did not occur in NT control plants, was observed under both normal and drought conditions. Our findings indicate that osa-miR171 plays a role in drought tolerance by regulating SCL6-I and SCL6-II transcript levels.

Genome-wide in silico analysis of long intergenic non-coding RNAs from rice peduncles at the heading stage

Long intergenic non-coding RNAs (lincRNAs) belong to the category of long non-coding RNAs (lncRNAs), originated from intergenic regions, which do not code for proteins. LincRNAs perform prominent role in regulation of gene expression during plant development and stress response by directly interacting with DNA, RNA, or proteins, or triggering production of small RNA regulatory molecules. Here, we identified 2973 lincRNAs and investigated their expression dynamics during peduncle elongation in two Indian rice cultivars, Pokkali and Swarna, at the time of heading. Differential expression analysis revealed common and cultivar-specific expression patterns, which we utilized to infer the lincRNA candidates with potential involvement in peduncle elongation and panicle exsertion. Their putative targets were identified using in silico prediction methods followed by pathway mapping and literature-survey based functional analysis. Further, to infer the mechanism of action, we identified the lincRNAs which potentially act as miRNA precursors or target mimics.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01059-2.

miR156/157 Targets SPLs to Regulate Flowering Transition, Plant Architecture and Flower Organ Size in Petunia

miR156/157 plays multiple pivotal roles during plant growth and development. In this study, we identified 11 miR156- and 5 miR157-encoding loci from the genome of Petunia axillaris and Petunia inflata, designated as PaMIR0156/157s and PiMIR0156/157s, respectively. Real-time quantitative reverse transcription PCR (qRT-PCR) analysis indicated that PhmiR156/157 was expressed predominantly in cotyledons, germinating seeds, flower buds, young fruits and seedlings. PhmiR156/157 levels declined in shoot apical buds and leaves of petunia before flowering as the plant ages; moreover, the temporal expression patterns of most miR156/157-targeted PhSPLs were complementary to that of PhmiR156/157. Ectopic expression of PhMIR0157a in Arabidopsis and petunia resulted in delayed flowering, dwarf plant stature, increased branches and reduced organ size. However, PhMIR0156f-overexpressing Arabidopsis and petunia plants showed only delayed flowering. In addition, downregulation of PhmiR156/157 level by overexpressing STTM156/157 led to taller plants with less branches, longer internodes and precocious flowering. qRT-PCR analysis indicated that PhmiR156/157 modulates these traits mainly by downregulating their PhSPL targets and subsequently decreasing the expression of flowering regulatory genes. Our results demonstrate that the PhmiR156/157-PhSPL module has conserved but also divergent functions in growth and development, which will help us decipher the genetic basis for the improvement of flower transition, plant architecture and organ development in petunia.

Genome-Wide Analysis of the IQM Gene Family in Rice (Oryza sativa L.)

Members of the IQM (IQ-Motif Containing) gene family are involved in plant growth and developmental processes, biotic and abiotic stress response. To systematically analyze the IQM gene family and their expression profiles under diverse biotic and abiotic stresses, we identified 8 IQM genes in the rice genome. In the current study, the whole genome identification and characterization of OsIQMs, including the gene and protein structure, genome localization, phylogenetic relationship, gene expression and yeast two-hybrid were performed. Eight IQM genes were classified into three subfamilies (I-III) according to the phylogenetic analysis. Gene structure and protein motif analyses showed that these IQM genes are relatively conserved within each subfamily of rice. The 8 OsIQM genes are distributed on seven out of the twelve chromosomes, with three IQM gene pairs involved in segmental duplication events. The evolutionary patterns analysis revealed that the IQM genes underwent a large-scale event within the last 20 to 9 million years. In addition, quantitative real-time PCR analysis of eight OsIQMs genes displayed different expression patterns at different developmental stages and in different tissues as well as showed that most IQM genes were responsive to PEG, NaCl, jasmonic acid (JA), abscisic acid (ABA) treatment, suggesting their crucial roles in biotic, and abiotic stress response. Additionally, a yeast two-hybrid assay showed that OsIQMs can interact with OsCaMs, and the IQ motif of OsIQMs is required for OsIQMs to combine with OsCaMs. Our results will be valuable to further characterize the important biological functions of rice IQM genes.

The GRAS gene family and its roles in seed development in litchi (Litchi chinensis Sonn)

BACKGROUND

The GRAS gene family plays crucial roles in multiple biological processes of plant growth, including seed development, which is related to seedless traits of litchi (Litchi chinensis Sonn.). However, it hasn't been fully identified and analyzed in litchi, an economic fruit tree cultivated in subtropical regions.

RESULTS

In this study, 48 LcGRAS proteins were identified and termed according to their chromosomal location. LcGRAS proteins can be categorized into 14 subfamilies through phylogenetic analysis. Gene structure and conserved domain analysis revealed that different subfamilies harbored various motif patterns, suggesting their functional diversity. Synteny analysis revealed that the expansion of the GRAS family in litchi may be driven by their tandem and segmental duplication. After comprehensively analysing degradome data, we found that four LcGRAS genes belong to HAM subfamily were regulated via miR171-mediated degradation. The various expression patterns of LcGRAS genes in different tissues uncovered they were involved in different biological processes. Moreover, the different temporal expression profiles of LcGRAS genes between abortive and bold seed indicated some of them were involved in maintaining the normal development of the seed.

CONCLUSION

Our study provides comprehensive analyses on GRAS family members in litchi, insight into a better understanding of the roles of GRAS in litchi development, and lays the foundation for further investigations on litchi seed development.

MADS-box transcription factors MADS11 and DAL1 interact to mediate the vegetative-to-reproductive transition in pine

The reproductive transition is an important event that is crucial for plant survival and reproduction. Relative to the thorough understanding of the vegetative phase transition in angiosperms, a little is known about this process in perennial conifers. To gain insight into the molecular basis of the regulatory mechanism in conifers, we used temporal dynamic transcriptome analysis with samples from seven different ages of Pinus tabuliformis to identify a gene module substantially associated with aging. The results first demonstrated that the phase change in P. tabuliformis occurred as an unexpectedly rapid transition rather than a slow, gradual progression. The age-related gene module contains 33 transcription factors and was enriched in genes that belong to the MADS (MCMl, AGAMOUS, DEFICIENS, SRF)-box family, including six SOC1-like genes and DAL1 and DAL10. Expression analysis in P. tabuliformis and a late-cone-setting P. bungeana mutant showed a tight association between PtMADS11 and reproductive competence. We then confirmed that MADS11 and DAL1 coordinate the aging pathway through physical interaction. Overexpression of PtMADS11 and PtDAL1 partially rescued the flowering of 35S::miR156A and spl1,2,3,4,5,6 mutants in Arabidopsis (Arabidopsis thaliana), but only PtMADS11 could rescue the flowering of the ft-10 mutant, suggesting PtMADS11 and PtDAL1 play different roles in flowering regulatory networks in Arabidopsis. The PtMADS11 could not alter the flowering phenotype of soc1-1-2, indicating it may function differently from AtSOC1 in Arabidopsis. In this study, we identified the MADS11 gene in pine as a regulatory mediator of the juvenile-to-adult transition with functions differentiated from the angiosperm SOC1.

Genome Identification and Expression Analysis of GRAS Family Related to Development, Hormone and Pathogen Stress in Brachypodium distachyon

GRAS transcription factors are widely present in the plant kingdom and play important roles in regulating multiple plant physiological processes. Brachypodium distachyon is a model for grasses for researching plant-pathogen interactions. However, little is known about the BdGRAS family genes involved in plant response to biotic stress. In this study, we identified 63 genes of the GRAS family in B. distachyon. The phylogenetic analysis showed that BdGRAS genes were divided into ten subfamilies and unevenly distributed on five chromosomes. qRT-PCR results showed that the BdGRAS family genes were involved in the growth and development of B. distachyon. Moreover, the expression of the HAM subfamily genes of BdGRAS changed during the interaction between B. distachyon and Magnaporthe oryzae. Interestingly, BdGRAS31 in the HAM subfamily was regulated by miR171 after inoculation with M. oryzae. These results provide insight into the potential functions of the BdGRAS family in disease resistance.

Conservation and diversification of HAIRY MERISTEM gene family in land plants

The shoot apical meristems (SAMs) of land plants are crucial for plant growth and organ formation. In several angiosperms, the HAIRY MERISTEM (HAM) genes function as key regulators that control meristem development and stem cell homeostasis. To date, the origin and evolutionary history of the HAM family in land plants remains unclear. Potentially shared and divergent functions of HAM family members from angiosperms and non-angiosperms are also not known. In constructing a comprehensive phylogeny of the HAM family, we show that HAM proteins are widely present in land plants and that HAM proteins originated prior to the divergence of bryophytes. The HAM family was duplicated in a common ancestor of angiosperms, leading to two distinct groups: type I and type II. Type-II HAM members are widely present in angiosperms, whereas type-I HAM members were independently lost in different orders of monocots. Furthermore, HAM members from angiosperms and non-angiosperms (including bryophytes, lycophytes, ferns and gymnosperms) are able to replace the role of the type-II HAM genes in Arabidopsis, maintaining established SAMs and promoting the initiation of new stem cell niches. Our results uncover the conserved functions of HAM family members and reveal the conserved regulatory mechanisms underlying HAM expression patterning in meristems, providing insight into the evolution of key stem cell regulators in land plants.

A De Novo Transcriptome Assembly of Ceratopteris richardii Provides Insights into the Evolutionary Dynamics of Complex Gene Families in Land Plants

As the closest extant sister group to seed plants, ferns are an important reference point to study the origin and evolution of plant genes and traits. One bottleneck to the use of ferns in phylogenetic and genetic studies is the fact that genome-level sequence information of this group is limited, due to the extreme genome sizes of most ferns. Ceratopteris richardii (hereafter Ceratopteris) has been widely used as a model system for ferns. In this study, we generated a transcriptome of Ceratopteris, through the de novo assembly of the RNA-seq data from 17 sequencing libraries that are derived from two sexual types of gametophytes and five different sporophyte tissues. The Ceratopteris transcriptome, together with 38 genomes and transcriptomes from other species across the Viridiplantae, were used to uncover the evolutionary dynamics of orthogroups (predicted gene families using OrthoFinder) within the euphyllophytes and identify proteins associated with the major shifts in plant morphology and physiology that occurred in the last common ancestors of euphyllophytes, ferns, and seed plants. Furthermore, this resource was used to identify and classify the GRAS domain transcriptional regulators of many developmental processes in plants. Through the phylogenetic analysis within each of the 15 GRAS orthogroups, we uncovered which GRAS family members are conserved or have diversified in ferns and seed plants. Taken together, the transcriptome database and analyses reported here provide an important platform for exploring the evolution of gene families in land plants and for studying gene function in seed-free vascular plants.

HAM Gene Family and Shoot Meristem Development

Land plants develop highly diversified shoot architectures, all of which are derived from the pluripotent stem cells in shoot apical meristems (SAMs). As sustainable resources for continuous organ formation in the aboveground tissues, SAMs play an important role in determining plant yield and biomass production. In this review, we summarize recent advances in understanding one group of key regulators - the HAIRY MERISTEM (HAM) family GRAS domain proteins - in shoot meristems. We highlight the functions of HAM family members in dictating shoot stem cell initiation and proliferation, the signaling cascade that shapes HAM expression domains in shoot meristems, and the conservation and diversification of HAM family members in land plants. We also discuss future directions that potentially lead to a more comprehensive view of the HAM gene family and stem cell homeostasis in land plants.

Profiling of MicroRNAs Involved in Mepiquat Chloride-Mediated Inhibition of Internode Elongation in Cotton (Gossypium hirsutum L.) Seedlings

Mepiquat chloride (MC) is the most important plant growth retardant that is widely used in cotton (Gossypium hirsutum L.) production to suppress excessive vegetative growth and improve plant architecture. MicroRNAs (miRNAs) are important gene expression regulators that control plant growth and development. However, miRNA-mediated post-transcriptional regulation in MC-induced growth inhibition remains unclear. In this study, the dynamic expression profiles of miRNAs responsive to MC in cotton internodes were investigated. A total of 508 known miRNAs belonging to 197 families and five novel miRNAs were identified. Among them, 104 miRNAs were differentially expressed at 48, 72, or 96 h post MC treatment compared with the control (0 h); majority of them were highly conserved miRNAs. The number of differentially expressed miRNAs increased with time after treatment. The expression of 14 known miRNAs was continuously suppressed, whereas 12 known miRNAs and one novel miRNA were continuously induced by MC. The expression patterns of the nine differentially expressed miRNAs were verified using qRT-PCR. The targets of the known and novel miRNAs were predicted. Four conserved and six novel targets were validated using the RLM-5' RACE assay. This study revealed that miRNAs play crucial regulatory roles in the MC-induced inhibition of internode elongation. It can improve our understanding of post-transcriptional gene regulation in MC-mediated growth inhibition and could potentially facilitate the breeding of dwarf cotton.

miRNA expression profiling and zeatin dynamic changes in a new model system of in vivo indirect regeneration of tomato

Callus formation and adventitious shoot differentiation could be observed on the cut surface of completely decapitated tomato plants. We propose that this process can be used as a model system to investigate the mechanisms that regulate indirect regeneration of higher plants without the addition of exogenous hormones. This study analyzed the patterns of trans-zeatin and miRNA expression during in vivo regeneration of tomato. Analysis of trans-zeatin revealed that the hormone cytokinin played an important role in in vivo regeneration of tomato. Among 183 miRNAs and 1168 predicted target genes sequences identified, 93 miRNAs and 505 potential targets were selected based on differential expression levels for further characterization. Expression patterns of six miRNAs, including sly-miR166, sly-miR167, sly-miR396, sly-miR397, novel 156, and novel 128, were further validated by qRT-PCR. We speculate that sly-miR156, sly-miR160, sly-miR166, and sly-miR397 play major roles in callus formation of tomato during in vivo regeneration by regulating cytokinin, IAA, and laccase levels. Overall, our microRNA sequence and target analyses of callus formation during in vivo regeneration of tomato provide novel insights into the regulation of regeneration in higher plants.

Genetic Basis of Tiller Dynamics of Rice Revealed by Genome-Wide Association Studies

A tiller number is the key determinant of rice plant architecture and panicle number and consequently controls grain yield. Thus, it is necessary to optimize the tiller number to achieve the maximum yield in rice. However, comprehensive analyses of the genetic basis of the tiller number, considering the development stage, tiller type, and related traits, are lacking. In this study, we sequence 219 Korean rice accessions and construct a high-quality single nucleotide polymorphism (SNP) dataset. We also evaluate the tiller number at different development stages and heading traits involved in phase transitions. By genome-wide association studies (GWASs), we detected 20 significant association signals for all traits. Five signals were detected in genomic regions near known candidate genes. Most of the candidate genes were involved in the phase transition from vegetative to reproductive growth. In particular, HD1 was simultaneously associated with the productive tiller ratio and heading date, indicating that the photoperiodic heading gene directly controls the productive tiller ratio. Multiple linear regression models of lead SNPs showed coefficients of determination (R²) of 0.49, 0.22, and 0.41 for the tiller number at the maximum tillering stage, productive tiller number, and productive tiller ratio, respectively. Furthermore, the model was validated using independent japonica rice collections, implying that the lead SNPs included in the linear regression model were generally applicable to the tiller number prediction. We revealed the genetic basis of the tiller number in rice plants during growth, By GWASs, and formulated a prediction model by linear regression. Our results improve our understanding of tillering in rice plants and provide a basis for breeding high-yield rice varieties with the optimum the tiller number.

Genetic Basis of Tiller Dynamics of Rice Revealed by Genome-Wide Association Studies

A tiller number is the key determinant of rice plant architecture and panicle number and consequently controls grain yield. Thus, it is necessary to optimize the tiller number to achieve the maximum yield in rice. However, comprehensive analyses of the genetic basis of the tiller number, considering the development stage, tiller type, and related traits, are lacking. In this study, we sequence 219 Korean rice accessions and construct a high-quality single nucleotide polymorphism (SNP) dataset. We also evaluate the tiller number at different development stages and heading traits involved in phase transitions. By genome-wide association studies (GWASs), we detected 20 significant association signals for all traits. Five signals were detected in genomic regions near known candidate genes. Most of the candidate genes were involved in the phase transition from vegetative to reproductive growth. In particular, HD1 was simultaneously associated with the productive tiller ratio and heading date, indicating that the photoperiodic heading gene directly controls the productive tiller ratio. Multiple linear regression models of lead SNPs showed coefficients of determination (R²) of 0.49, 0.22, and 0.41 for the tiller number at the maximum tillering stage, productive tiller number, and productive tiller ratio, respectively. Furthermore, the model was validated using independent japonica rice collections, implying that the lead SNPs included in the linear regression model were generally applicable to the tiller number prediction. We revealed the genetic basis of the tiller number in rice plants during growth, By GWASs, and formulated a prediction model by linear regression. Our results improve our understanding of tillering in rice plants and provide a basis for breeding high-yield rice varieties with the optimum the tiller number.

Genome-Wide Analysis of the GRAS Gene Family in Barley (Hordeum vulgare L.)

The GRAS (named after first three identified proteins within this family, GAI, RGA, and SCR) family contains plant-specific genes encoding transcriptional regulators that play a key role in gibberellin (GA) signaling, which regulates plant growth and development. Even though GRAS genes have been characterized in some plant species, little research is known about the GRAS genes in barley (Hordeum vulgare L.). In this study, we observed 62 GRAS members from barley genome, which were grouped into 12 subgroups by using phylogenomic analysis together with the GRAS genes from Arabidopsis (Arabidopsis thaliana), maize (Zea mays), and rice (Oryza sativa). Chromosome localization and gene structure analysis suggested that duplication events and abundant presence of intronless genes might account for the massive expansion of GRAS gene family in barley. The analysis of RNA-seq data indicates the expression pattern of GRAS genes in various tissues at different stages in barley. Noteworthy, our qRT-PCR analysis showed the expression of 18 candidate GRAS genes abundantly in the developing inflorescence, indicating their potential roles in the barley inflorescence development and reproduction. Collectively, our evolutionary and expression analysis of GRAS family are useful for future functional characterization of GA signaling in barley and agricultural improvement.

Genome-wide survey and expression analyses of the GRAS gene family in Brassica napus reveals their roles in root development and stress response

Genome-wide identification, classification, expression analyses, and functional characterization of GRAS genes in oil crop, Brassica napus, indicate their importance in root development and stress response. GRAS proteins are a plant-specific transcription factor gene family involved in tissues development and stress response. We classified 87 putative GRAS genes in the Brassica napus genome (BnGRASs) into 13 subfamilies by phylogenetic analysis. The C-terminal GRAS domains of Brassica napus (B. napus) proteins were less conserved among subfamilies, but were conserved within each subfamily. A series of analyses revealed that 89.7% of the BnGRASs did not have intron insertions, and 24 specific-motifs were found at the N-terminal. A highly conserved microRNA 171 (miRNA171) target was observed specifically in the HAM subfamily across land plants. A total of 868 pairs of interaction proteins were predicted, the primary of which were transcription factors involved in transcriptional regulation and signal transduction. Integrated comparative analysis of GRAS genes across 26 species of algae, mosses, ferns, gymnosperms, and angiosperms revealed that this gene family originated in early mosses and was classified into 19 subfamilies, 14 of which may have originated prior to bryophyte evolution. RNA-Seq analysis demonstrated that most BnGRASs were widely expressed in different tissues/organs at different stages in B. napus, and 24 BnGRASs were highly/specifically expressed in roots. Results from a qRT-PCR analysis suggested that two BnGRASs belonging to SCR and LISCL subfamilies potentially have important roles in the stress response of roots.

Overexpression of SlGRAS7 Affects Multiple Behaviors Leading to Confer Abiotic Stresses Tolerance and Impacts Gibberellin and Auxin Signaling in Tomato

Abiotic stresses remain the key environmental issues that reduce plant development and therefore affect crop production. Transcription factors, such as the GRAS family, are involved in various functions of abiotic stresses and plant growth. The GRAS family of tomato (Solanum lycopersicum), SlGRAS7, is described in this study. We produced overexpressing SlGARS7 plants to learn more about the GRAS transcription factors. Plants overexpressing SlGARS7 (SlGRAS7-OE) showed multiple phenotypes related to many behaviors, including plant height, root and shoot length, and flowering time. We observed that many genes in the SlGRAS7-OE seedlings that are associated with auxin and gibberellin (GA) are downregulated and have altered sensitivity to GA₃/IAA. SlGRAS7 was upregulated during abiotic stresses following treatment with sodium chloride (NaCl) and D-mannitol in the wild-type (WT) tomato. Tomato plants overexpressing SlGRAS7 showed more resistance to drought and salt stress comparison with WT. Our study of SlGRAS7 in tomato demonstrates how GRAS showed an integrative role, improving resistance to abiotic stresses and enhancing gibberellin/auxin signaling through reproductive as well as vegetative processes.

miRNAs associated with auxin signaling, stress response, and cellular activities mediate adventitious root formation in apple rootstocks

Adventitious root (AR) formation is essential for the vegetative propagation of apple rootstocks. miRNAs play a significant role in regulating AR development, however, large-scale transcriptomic data on miRNA mediated AR formation in apple rootstocks is lacking. Therefore, in order to identify the molecular mechanisms underlying AR formation in 'M9-T337' apple rootstocks, transcriptomic changes occurring during key time points of AR formation (0, 3, and 16 days) were analyzed using high-throughput sequencing with a focus on miRNAs. A total of 84 known miRNAs and 56 novel miRNAs have differentially expressed were identified. Additionally, a total of 88 target genes of known miRNAs and 76 target genes of novel miRNAs were identified by degradome sequencing. The expression levels of the miRNAs and target genes were quantified by RT-qPCR. Results indicate that miRNAs and their target genes are associated with auxin signal-related (miR160 and miR390), stress response-related (miR398, miR395 and miR408), cell fate transformation-, proliferation- and enlargement-related (miR171, miR156, miR166, miR319 and miR396). These all involve pathways that participate in AR formation in 'M9-T337' apple rootstock. In addition, hormones (AUX, CTK, GA3, BR, JA, and ABA) are also involved in regulating AR formation. The candidate genes belonging to pathways associated with AR formation exhibited significantly higher expression levels, providing evidence that they may be involved in the regulation of AR development. The collective results of the present study indicate that the developmental process associated with AR formation in apple rootstock is extremely complex. The known and novel miRNAs and target genes that were identified by high-throughput and degradome sequencing, respectively, provide a framework for the future analysis of miRNAs associated with AR development in apple rootstocks, and provide new information that can be used to better understand AR development in woody plants.

Current Advances in Molecular Mechanisms and Physiological Basis of Panicle Degeneration in Rice

Panicle degeneration, also known as panicle abortion, is a serious defect and causes heavy losses to reproductive yield in cereals. Several mutants have been reported to display the phenotype of spikelet abortion in rice. Recent findings have resulted in significant breakthroughs, but comprehensive understanding about the molecular pathways and physiological basis of panicle degeneration still remain a dilemma. In this review, we have summarized all the responsible genes and mechanisms underlying the panicle development with a special focus on degeneration. Here, we hypothesized a model by using knowledge and coherent logic in order to understand the molecular regulation of panicle degeneration. In addition to this, we included all the previous discoveries, schools of thoughts, ancient working theories, and crosstalk of phytohormones and provided new insights for future studies.

A genome-wide association study using a Vietnamese landrace panel of rice (Oryza sativa) reveals new QTLs controlling panicle morphological traits

CONTEXT

Yield improvement is an important issue for rice breeding. Panicle architecture is one of the key components of rice yield and exhibits a large diversity. To identify the morphological and genetic determinants of panicle architecture, we performed a detailed phenotypic analysis and a genome-wide association study (GWAS) using an original panel of Vietnamese landraces.

RESULTS

Using a newly developed image analysis tool, morphological traits of the panicles were scored over two years: rachis length; primary, secondary and tertiary branch number; average length of primary and secondary branches; average length of internode on rachis and primary branch. We observed a high contribution of spikelet number and secondary branch number per panicle to the overall phenotypic diversity in the dataset. Twenty-nine stable QTLs associated with seven traits were detected through GWAS over the two years. Some of these QTLs were associated with genes already implicated in panicle development. Importantly, the present study revealed the existence of new QTLs associated with the spikelet number, secondary branch number and primary branch number traits.

CONCLUSIONS

Our phenotypic analysis of panicle architecture variation suggests that with the panel of samples used, morphological diversity depends largely on the balance between indeterminate vs. determinate axillary meristem fate on primary branches, supporting the notion of differences in axillary meristem fate between rachis and primary branches. Our genome-wide association study led to the identification of numerous genomic sites covering all the traits studied and will be of interest for breeding programs aimed at improving yield. The new QTLs detected in this study provide a basis for the identification of new genes controlling panicle development and yield in rice.

Ectopic Overexpression of bol-miR171b Increases Chlorophyll Content and Results in Sterility in Broccoli ( Brassica oleracea L var. italica)

MiR171 plays pleiotropic roles in the growth and development of several plant species. However, the mechanism underlying the miR171-mediated regulation of organ development in broccoli remains unknown. In this study, bol-miR171b was characterized and found to be differentially expressed in various broccoli organs. The ectopic overexpression of bol-miR171b in Arabidopsis affected the leaf and silique development of transgenic lines. In particular, the chlorophyll content of leaves from overexpressed bol-miR171b transgenic Arabidopsis was higher than that of the vector controls. The fertility and seed yield of Arabidopsis with overexpressed bol-miR171b were markedly lower than those of the vector controls. Similarly, overexpressed bol-miR171b transgenic broccoli exhibited dark green leaves with high chlorophyll content, and nearly all of the flowers were sterile. These results demonstrated that overexpression of bol-miR171b could increase the chlorophyll content of transgenic plants. Degradome sequencing was conducted to identify the targets of bol-miR171b. Two members of the GRAS gene family, BolSCL6 and BolSCL27, were cleaved by bol-miR171b-3p in broccoli. In addition to the genes targeted by bol-miR171b-3p, adenylylsulfate reductase 3 ( APSR3), which played important roles in plant sulfate assimilation and reduction, was speculated to be cleaved by bol-miR171b-5p, suggesting that the star sequence of bol-miR171b may also have functions in broccoli. Comparative transcriptome analysis further revealed that the genes involved in chloroplast development and sulfate homeostasis should participate in the bol-miR171b -mediated regulatory network. Taken together, these findings provided new insights into the function and regulation of bol-miR171b in broccoli and indicated the potential of bol-miR171b as a small RNA molecule that increased leaf chlorophyll in plants by genetic engineering.

Combined analysis and miRNA expression profiles of the flowering related genes in common wild rice (oryza rufipogon Griff.)

Common wild rice (Oryza rufipogon Griff.) is the most closely related ancestral species to Asian cultivated rice (Oryza sativa L.). It contains various valuable traits with regard to tolerance to cold, drought and salinity, flowering diversity and many quantitative trait loci with agronomic important traits. Flowering is one of the most important agronomic traits. However, flowering-related transcriptome and how to be regulated by miRNAs have not been estimated in O.rufipogon. To identify how the genes and miRNAs regulating flowering in O.rufipogon. Three O.rufipogon RNA libraries, two vegetative stages (CWRT-V1 and CWRT-V2) and one flowering stage (CWRT-F2) were constructed using leaves tissue and sequenced using Illumina deep sequencing. 27,405, 27,333, 28,979 unique genes were obtained after mapping to the reference genome from CWRT-V1, CWRT-V2 and CWRT-F2, respectively. Then differentially expressed genes (DEGs) were screened and got 1419 unique genes are likely to involve in flower development. Detailed information showed that MADS box and floral meristem identity genes, such as MADS 1, MADS14, Hd1 are involved in common wild rice. Then, combined analysis of miRNA and mRNA expression profiles was performed. Twenty three known miRNA-mRNA pairs and five new candidates were presented an anti-correlationship. Interestingly, 12 miRNAs were negatively correlated with 20 mRNAs encoding flowering-related proteins, indicating that miRNAs regulated target genes to promote flowering in CWRT-F2 group. The results provided here genomic resources for flowering related genes and how these flowering genes were regulated by miRNAs in common wild rice.

Under a New Light: Regulation of Light-Dependent Pathways by Non-coding RNAs

The biological relevance of non-protein coding RNAs in the regulation of critical plant processes has been firmly established in recent years. This has been mostly achieved with the discovery and functional characterization of small non-coding RNAs, such as small interfering RNAs and microRNAs (miRNAs). However, recent next-generation sequencing techniques have widened our view of the non-coding RNA world, which now includes long non-coding RNAs (lncRNAs). Small and lncRNAs seem to diverge in their biogenesis and mode of action, but growing evidence highlights their relevance in developmental processes and in responses to particular environmental conditions. Light can affect MIRNA gene transcription, miRNA biogenesis, and RNA-induced silencing complex (RISC) activity, thus controlling not only miRNA accumulation but also their biological function. In addition, miRNAs can mediate several light-regulated processes. In the lncRNA world, few reports are available, but they already indicate a role in the regulation of photomorphogenesis, cotyledon greening, and photoperiod-regulated flowering. In this review, we will discuss how light controls MIRNA gene expression and the accumulation of their mature forms, with a particular emphasis on those miRNAs that respond to different light qualities and are conserved among species. We will also address the role of small non-coding RNAs, particularly miRNAs, and lncRNAs in the regulation of light-dependent pathways. We will mainly focus on the recent progress done in understanding the interconnection between these non-coding RNAs and photomorphogenesis, circadian clock function, and photoperiod-dependent flowering.

Identification and characterization of Prunus persica miRNAs in response to UVB radiation in greenhouse through high-throughput sequencing

Background

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression of target mRNAs involved in plant growth, development, and abiotic stress. As one of the most important model plants, peach (Prunus persica) has high agricultural significance and nutritional values. It is well adapted to be cultivated in greenhouse in which some auxiliary conditions like temperature, humidity, and UVB etc. are needed to ensure the fruit quality. However, little is known about the genomic information of P. persica under UVB supplement. Transcriptome and expression profiling data for this species are therefore important resources to better understand the biological mechanism of seed development, formation and plant adaptation to environmental change. Using a high-throughput miRNA sequencing, followed by qRT-PCR tests and physiological properties determination, we identified the responsive-miRNAs under low-dose UVB treatment and described the expression pattern and putative function of related miRNAs and target genes in chlorophyll and carbohydrate metabolism.

Results

A total of 164 known peach miRNAs belonging to 59 miRNA families and 109 putative novel miRNAs were identified. Some of these miRNAs were highly conserved in at least four other plant species. In total, 1794 and 1983 target genes for known and novel miRNAs were predicted, respectively. The differential expression profiles of miRNAs between the control and UVB-supplement group showed that UVB-responsive miRNAs were mainly involved in carbohydrate metabolism and signal transduction. UVB supplement stimulated peach to synthesize more chlorophyll and sugars, which was verified by qRT-PCR tests of related target genes and metabolites’ content measurement.

Conclusion

The high-throughput sequencing data provided the most comprehensive miRNAs resource available for peach study. Our results identified a series of differentially expressed miRNAs/target genes that were predicted to be low-dose UVB-responsive. The correlation between transcriptional profiles and metabolites contents in UVB supplement groups gave novel clues for the regulatory mechanism of miRNAs in Prunus. Low-dose UVB supplement could increase the chlorophyll and sugar (sorbitol) contents via miRNA-target genes and therefore improve the fruit quality in protected cultivation of peaches.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12864-017-4347-5) contains supplementary material, which is available to authorized users.

Genome-wide characterization of GRAS family genes in Medicago truncatula reveals their evolutionary dynamics and functional diversification

The GRAS gene family is a large plant-specific family of transcription factors that are involved in diverse processes during plant development. Medicago truncatula is an ideal model plant for genetic research in legumes, and specifically for studying nodulation, which is crucial for nitrogen fixation. In this study, 59 MtGRAS genes were identified and classified into eight distinct subgroups based on phylogenetic relationships. Motifs located in the C-termini were conserved across the subgroups, while motifs in the N-termini were subfamily specific. Gene duplication was the main evolutionary force for MtGRAS expansion, especially proliferation of the LISCL subgroup. Seventeen duplicated genes showed strong effects of purifying selection and diverse expression patterns, highlighting their functional importance and diversification after duplication. Thirty MtGRAS genes, including NSP1 and NSP2, were preferentially expressed in nodules, indicating possible roles in the process of nodulation. A transcriptome study, combined with gene expression analysis under different stress conditions, suggested potential functions of MtGRAS genes in various biological pathways and stress responses. Taken together, these comprehensive analyses provide basic information for understanding the potential functions of GRAS genes, and will facilitate further discovery of MtGRAS gene functions.

Identification of protoplast-isolation responsive microRNAs in Citrus reticulata Blanco by high-throughput sequencing

Protoplast isolation is a stress-inducing process, during which a variety of physiological and molecular alterations take place. Such stress response affects the expression of totipotency of cultured protoplasts. MicroRNAs (miRNAs) play important roles in plant growth, development and stress responses. However, the underlying mechanism of miRNAs involved in the protoplast totipotency remains unclear. In this study, high-throughput sequencing technology was used to sequence two populations of small RNA from calli and callus-derived protoplasts in Citrus reticulata Blanco. A total of 67 known miRNAs from 35 families and 277 novel miRNAs were identified. Among these miRNAs, 18 known miRNAs and 64 novel miRNAs were identified by differentially expressed miRNAs (DEMs) analysis. The expression patterns of the eight DEMs were verified by qRT-PCR. Target prediction showed most targets of the miRNAs were transcription factors. The expression levels of half targets showed a negative correlation to those of the miRNAs. Furthermore, the physiological analysis showed high levels of antioxidant activities in isolated protoplasts. In short, our results indicated that miRNAs may play important roles in protoplast-isolation response.

Altered accumulation of osa-miR171b contributes to rice stripe virus infection by regulating disease symptoms

Viral infection affects the pattern of plant miRNA expression. It has been presumed that reduction of miR171 and several other miRNAs influences viral symptoms in plants. We here experimentally demonstrate the association of osa-miR171b with rice stripe virus (RSV) symptoms in rice. Inhibition of osa-miR171b caused stunting with reduced chlorophyll content in leaves similar to viral symptoms. Overexpression of osa-miR171b by an artificial miRNA extended vegetative growth and enhanced chlorophyll accumulation in leaves. Tillers were thicker, and panicles were longer with more spikelets in plants overexpressing osa-miR171b than in controls, but there were no differences in tiller numbers. Targets of osa-miR171b, OsSCL6-IIa, OsSCL6-IIb, and OsSCL6-IIc, were respectively up- and down-regulated in plants where osa-miR171b was inhibited or overexpressed. In plants overexpressing osa-miR171b, five positive regulators for heading development, Ehd1, Ehd2, Ehd3, Ehd4, and Hd3a were up-regulated, while the negative regulator Ghd7 was down-regulated. Plants overexpressing osa-miR171b were less susceptible to RSV and virus symptoms were attenuated. Taken together, the results reveal that a reduction of osa-miR171b in RSV-infected rice contributes to RSV symptoms, and provide more insight into the roles of osa-miR171b in rice.

Short tandem target mimic rice lines uncover functions of miRNAs in regulating important agronomic traits

Improvements in plant agricultural productivity are urgently needed to reduce the dependency on limited natural resources and produce enough food for a growing world population. Human intervention over thousands of years has improved the yield of important crops; however, it is increasingly difficult to find new targets for genetic improvement. MicroRNAs (miRNAs) are promising targets for crop improvement, but their inactivation is technically challenging and has hampered functional analyses. We have produced a large collection of transgenic short tandem target mimic (STTM) lines silencing 35 miRNA families in rice as a resource for functional studies and crop improvement. Visual assessment of field-grown miRNA-silenced lines uncovered alterations in many valuable agronomic traits, including plant height, tiller number, and grain number, that remained stable for up to five generations. We show that manipulation of miR398 can increase panicle length, grain number, and grain size in rice. In addition, we discovered additional agronomic functions for several known miRNAs, including miR172 and miR156. Our collection of STTM lines thus represents a valuable resource for functional analysis of rice miRNAs, as well as for agronomic improvement that can be readily transferred to other important food crops.

Overexpression of a tomato miR171 target gene SlGRAS24 impacts multiple agronomical traits via regulating gibberellin and auxin homeostasis

In Arabidopsis, the miR171-GRAS module has been clarified as key player in meristem maintenance. However, the knowledge about its role in fruit crops like tomato (Solanum lycopersicum) remains scarce. We previously identified tomato SlGRAS24 as a target gene of Sly-miR171. To study the role of this probable transcription factor, we generated transgenic tomato plants underexpressing SlGRAS24, overexpressing SlGRAS24, overexpressing Sly-miR171 and expressing β-glucuronidase (GUS) under the SlGRAS24 promoter (proSlGRAS24-GUS). Plants overexpressing SlGRAS24 (SlGRAS24-OE) had pleiotropic phenotypes associated with multiple agronomical traits including plant height, flowering time, leaf architecture, lateral branch number, root length, fruit set and development. Many GA/auxin-related genes were down-regulated and altered responsiveness to exogenous IAA/NAA or GA₃ application was observed in SlGRAS24-OE seedlings. Moreover, compromised fruit set and development in SlGRAS24-OE was also observed. These newly identified phenotypes for SlGRAS24 homologs in tomato were later proved to be caused by impaired pollen sacs and fewer viable pollen grains. At anthesis, the comparative transcriptome results showed altered expression of genes involved in pollen development and hormone signalling. Taken together, our data demonstrate that SlGRAS24 participates in a series of developmental processes through modulating gibberellin and auxin signalling, which sheds new light on the involvement of hormone crosstalk in tomato development.

miRNA and Degradome Sequencing Reveal miRNA and Their Target Genes That May Mediate Shoot Growth in Spur Type Mutant "Yanfu 6"

The spur-type growth habit in apple trees is characterized by short internodes, increased number of fruiting spurs, and compact growth that promotes flowering and facilitates management practices, such as pruning. The molecular mechanisms responsible for regulating spur-type growth have not been elucidated. In the present study, miRNAs and the expression of their potential target genes were evaluated in shoot tips of "Nagafu 2" (CF) and spur-type bud mutation "Yanfu 6" (YF). A total of 700 mature miRNAs were identified, including 202 known apple miRNAs and 498 potential novel miRNA candidates. A comparison of miRNA expression in CF and YF revealed 135 differentially expressed genes, most of which were downregulated in YF. YF also had lower levels of GA, ZR, IAA, and ABA hormones, relative to CF. Exogenous applications of GA promoted YF shoot growth. Based on the obtained results, a regulatory network involving plant hormones, miRNA, and their potential target genes is proposed for the molecular mechanism regulating the growth of YF. miRNA164, miRNA166, miRNA171, and their potential targets, and associated plant hormones, appear to regulate shoot apical meristem (SAM) growth. miRNA159, miRNA167, miRNA396, and their potential targets, and associated plant hormones appear to regulate cell division and internode length. This study provides a foundation for further studies designed to elucidate the mechanism underlying spur-type apple architecture.

Expression Variations of miRNAs and mRNAs in Rice (Oryza sativa)

Differences in expression levels are an important source of phenotypic variation within and between populations. MicroRNAs (miRNAs) are key players in post-transcriptional gene regulation that are important for plant development and stress responses. We surveyed expression variation of miRNAs and mRNAs of six accessions from two rice subspecies Oryza sativa L. ssp. indica and Oryza sativa L. ssp. japonica using deep sequencing. While more than half (53.7%) of the mature miRNAs exhibit differential expression between grains and seedlings of rice, only 11.0% show expression differences between subspecies, with an additional 2.2% differentiated for the development-by-subspecies interaction. Expression variation is greater for lowly conserved miRNAs than highly conserved miRNAs, whereas the latter show stronger negative correlation with their targets in expression changes between subspecies. Using a permutation test, we identified 51 miRNA–mRNA pairs that correlate negatively or positively in expression level among cultivated rice. Genes involved in various metabolic processes and stress responses are enriched in the differentially expressed genes between rice indica and japonica subspecies. Our results indicate that stabilizing selection is the major force governing miRNA expression in cultivated rice, albeit positive selection may be responsible for much of the between-subspecies expression divergence.

Over-Expression of GmGIa-Regulated Soybean miR172a Confers Early Flowering in Transgenic Arabidopsis thaliana

Flowering is a pivotal event in the life cycle of plants. miR172 has been widely confirmed to play critical roles in flowering time control by regulating its target gene expression in Arabidopsis. However, the role of its counterpart in soybean remains largely unclear. In the present study, we found that the gma-miR172a was regulated by a GIGANTEA ortholog, GmGIa, in soybean through miRNA metabolism. The expression analysis revealed that gma-miR172a has a pattern of diurnal rhythm expression and its abundance increased rapidly as plants grew until the initiation of flowering phase in soybean. One target gene of gma-miR172a, Glyma03g33470, was predicted and verified using a modified RLM 5'-RACE (RNA ligase-mediated rapid amplification of 5' cDNA ends) assay. Overexpression of gma-miR172a exhibited an early flowering phenotype and the expression of FT, AP1 and LFY were simultaneously increased in gma-miR172a-transgenic Arabidopsis plants, suggesting that the early flowering phenotype was associated with up-regulation of these genes. The overexpression of the gma-miR172a-resistant version of Glyma03g33470 weakened early flowering phenotype in the toe1 mutant of Arabidopsis. Taken together, our results suggested that gma-miR172a played an important role in GmGIa-mediated flowering by repressing Glyma03g33470, which in turn increased the expression of FT, AP1 and LFY to promote flowering in soybean.

Ready, aim, shoot: stem cell regulation of the shoot apical meristem

Plant shoot meristems contain stem cells that are continuously renewed to replenish cells that exit and differentiate during lateral organ formation. Complex cell-to-cell signaling systems balance division and differentiation. These center on ligand-receptor networks, hormone pathways, and transcriptional regulators that function in an integrated manner. In this review, we aim to highlight new findings in shoot stem cell regulation across species.