Growth-Regulating Factors (GRFs): A Small Transcription Factor Family with Important Functions in Plant Biology

miR396b/GRF6 module contributes to salt tolerance in rice

Salinity, as one of the most challenging environmental factors restraining crop growth and yield, poses a severe threat to global food security. To address the rising food demand, it is urgent to develop crop varieties with enhanced yield and greater salt tolerance by delving into genes associated with salt tolerance and high-yield traits. MiR396b/GRF6 module has previously been demonstrated to increase rice yield by shaping the inflorescence architecture. In this study, we revealed that miR396b/GRF6 module can significantly improve salt tolerance of rice. In comparison with the wild type, the survival rate of MIM396 and OE-GRF6 transgenic lines increased by 48.0% and 74.4%, respectively. Concurrent with the increased salt tolerance, the transgenic plants exhibited reduced H2O2 accumulation and elevated activities of ROS-scavenging enzymes (CAT, SOD and POD). Furthermore, we identified ZNF9, a negative regulator of rice salt tolerance, as directly binding to the promoter of miR396b to modulate the expression of miR396b/GRF6. Combined transcriptome and ChIP-seq analysis showed that MYB3R serves as the downstream target of miR396b/GRF6 in response to salt tolerance, and overexpression of MYB3R significantly enhanced salt tolerance. In conclusion, this study elucidated the potential mechanism underlying the response of the miR396b/GRF6 network to salt stress in rice. These findings offer a valuable genetic resource for the molecular breeding of high-yield rice varieties endowed with stronger salt tolerance.

Genome-wide characterization of the GRF transcription factors in potato (Solanum tuberosum L.) and expression analysis of StGRF genes during potato tuber dormancy and sprouting

Growth-regulating factors (GRFs) are transcription factors that play a pivotal role in plant growth and development. This study identifies 12 Solanum tuberosum GRF transcription factors (StGRFs) and analyzes their physicochemical properties, phylogenetic relationships, gene structures and gene expression patterns using bioinformatics. The StGRFs exhibit a length range of 266 to 599 amino acids, with a molecular weight of 26.02 to 64.52 kDa. The majority of StGRFs possess three introns. The promoter regions contain a plethora of cis-acting elements related to plant growth and development, as well as environmental stress and hormone response. All the members of the StGRF family contain conserved WRC and QLQ domains, with the sequences of these two conserved domain modules exhibiting high levels of conservation. Transcriptomic data indicates that StGRFs play a significant role in the growth and development of stamens, roots, young tubers, and other tissues or organs in potatoes. Furthermore, a few StGRFs exhibit differential expression patterns in response to Phytophthora infestans, chemical elicitors, heat, salt, and drought stresses, as well as multiple hormone treatments. The results of the expression analysis indicate that StGRF1, StGRF2, StGRF5, StGRF7, StGRF10 and StGRF12 are involved in the process of tuber sprouting, while StGRF4 and StGRF9 may play a role in tuber dormancy. These findings offer valuable insights that can be used to investigate the roles of StGRFs during potato tuber dormancy and sprouting.

Differential expression of microRNAs in response to Papaya ringspot virus infection in differentially responding genotypes of papaya (Carica papaya L.) and its wild relative

Papaya ringspot virus (PRSV) is one of the most devastating viruses of papaya that has significantly hampered papaya production across the globe. Although PRSV resistance is known in some of its wild relatives, such as Vasconcellea cauliflora and in some of the improved papaya genotypes, the molecular basis of this resistance mechanism has not been studied and understood. Plant microRNAs are an important class of small RNAs that regulate the gene expression in several plant species against the invading plant pathogens. These miRNAs are known to manifest the expression of genes involved in resistance against plant pathogens, through modulation of the plant's biochemistry and physiology. In this study we made an attempt to study the overall expression pattern of small RNAs and more specifically the miRNAs in different papaya genotypes from India, that exhibit varying levels of tolerance or resistance to PRSV. Our study found that the expression of some of the miRNAs was differentially regulated in these papaya genotypes and they had entirely different miRNA expression profile in healthy and PRSV infected symptomatic plants. This data may help in improvement of papaya cultivars for resistance against PRSV through new breeding initiatives or biotechnological approaches such as genome editing.

Global gene expression profile and functional analysis reveal the conservation of reproduction-associated gene networks in Gossypium hirsutum

KEY MESSAGE

Lastly, the bZIP gene family encompasses genes that have been reported to play a role in flower development, such as bZIP14 (FD). Notably, bZIP14 is essential for Flowering Locus T (FT) initiation of floral development in Arabidopsis (Abe et al. 2005). Cotton (Gossypium hirsutum L.) is the world's most extensively cultivated fiber crop. However, its reproductive development is poorly characterized at the molecular level. Thus, this study presents a detailed transcriptomic analysis of G. hirsutum at three different reproductive stages. We provide evidence that more than 64,000 genes are active in G. hirsutum during flower development, among which 94.33% have been assigned to functional terms and specific pathways. Gene set enrichment analysis (GSEA) revealed that the biological process categories of floral organ development, pollen exine formation, and stamen development were enriched among the genes expressed during the floral development of G. hirsutum. Furthermore, we identified putative Arabidopsis homologs involved in the G. hirsutum gene regulatory network (GRN) of pollen and flower development, including transcription factors such as WUSCHEL (WUS), INNER NO OUTER (INO), AGAMOUS-LIKE 66 (AGL66), SPOROCYTELESS/NOZZLE (SPL/NZZ), DYSFUNCTIONAL TAPETUM 1 (DYT1), ABORTED MICROSPORES (AMS), and ASH1-RELATED 3 (ASHR3), which are known crucial genes for plant reproductive success. The cotton MADS-box protein-protein interaction pattern resembles the previously described patterns for AGAMOUS (AG), SEEDSTICK (STK), SHATTERPROOF (SHP), and SEPALLATA3 (SEP3) homolog proteins from Arabidopsis. In addition to serving as a resource for comparative flower development studies, this work highlights the changes in gene expression profiles and molecular networks underlying stages that are valuable for cotton breeding improvement.

Comparative genomic studies on the TGF-β superfamily in blue whale

TGF-β supergene family has a wide range of physiological functions including cell adhesion, motility, proliferation, apoptosis, and differentiation. We systematically analyzed and characterized the TGF-β gene superfamily from the whole blue whale (Balaenoptera musculus) genome, using comparative genomic and evolutionary analysis. We identified 30 TGF-β genes and were split into two subgroups, BMP-like and TGF-like. All TGF-β proteins demonstrating a basic nature, with the exception of BMP1, BMP2, BMP10, GDF2, MSTN, and NODAL modulator, had acidic characteristics. All the blue whale (B. musculus) TGF-β proteins, excluding BMP1, are thermostable based on aliphatic index. The instability index showed all proteins except the NODAL modulator was unstable. TGF-β proteins showed a hydrophilic character, with the exception of GDF1 and INHBC. Moreover, all the detected TGF-β genes showed evolutionary conserved nature. A segmental duplication was indicated by TGF-β gene family, and the Ka/Ks ratio showed that the duplicated gene pairs were subjected to selection pressure, indicating both purifying and positive selection pressure. Two possible recombination breakpoints were also predicted. This study provides insights into the genetic characterization and evolutionary aspects of the TGF-β superfamily in blue whales (B. musculus).

MIR396-GRF/GIF enhances in planta shoot regeneration of Dendrobium catenatum

Recent studies on co-transformation of the growth regulator, TaGRF4-GIF1 chimera (Growth Regulating Factor 4-GRF Interacting Factor 1), in cultivated wheat varieties (Triticum aestivum), showed improved regeneration efficiency, marking a significant breakthrough. Here, a simple and reproducible protocol using the GRF4-GIF1 chimera was established and tested in the medicinal orchid Dendrobium catenatum, a monocot orchid species. TaGRF4-GIF1 from T. aestivum and DcGRF4-GIF1 from D. catenatum were reconstructed, with the chimeras significantly enhancing the regeneration efficiency of D. catenatum through in planta transformation. Further, mutating the microRNA396 (miR396) target sites in TaGRF4 and DcGRF4 improved regeneration efficiency. The target mimicry version of miR396 (MIM396) not only boosted shoot regeneration but also enhanced plant growth. Our methods revealed a powerful tool for the enhanced regeneration and genetic transformation of D. catenatum.

Blocking miR396 activity by overexpression MIM396 improved switchgrass tiller number and biomass yield

BACKGROUND

MicroRNA396 (miR396) plays an important role in the regulation of plant growth and development by repressing the expression level of its target growth-regulating factor (GRF) family genes. In our previous study, we found that overexpression of miR396 negatively regulated both tillering and biomass yield in switchgrass (Panicum virgatum L.). We, therefore, speculated that blocking the expression of miR396 could enhance switchgrass tillering and biomass yield. Here, we produced transgenic switchgrass plants overexpressing a target mimicry form of miR396 (MIM396) in wild type (WT) and Os-MIR319b overexpressing switchgrass plant (with higher enzymatic hydrolysis efficiency, but reduced tillering), in which the expression of miR396 was blocked. The phenotype and biological yields of these plants were analyzed.

RESULTS

Blocking miR396 to improve its target PvGRFs expression in switchgrass improved the tiller number and dry weight of transgenic plants. Further morphological analysis revealed that MIM396 plants increased the number of aerial branches and basal tillers compared to those of wild-type plants. The enzymatic efficiency of MIM396 plants was reduced; however, the total sugar production per plant was still significantly higher than that of wild-type plants due to the increase in biomass. In addition, blocking miR396 in a transgenic switchgrass plant overexpressing Os-MIR319b (TG21-Ms) significantly increased the PvGRF1/3/5 expression level and tiller number and biomass yield. The miR156-target gene PvSPL4, playing a negative role in aerial and basal buds outgrowth, showed significant downregulated in MIM396 and TG21-Ms. Those results indicate that miR396-PvGRFs, through disrupting the PvSPL4 expression, are involved in miR319-PvPCFs in regulating tiller number, at least partly.

CONCLUSIONS

MIM396 could be used as a molecular tool to improving tiller number and biomass yield in switchgrass wild type and miR319b transgenic plants. This finding may be applied to other graminaceous plants to regulate plant biological yield.

Genome-wide identification and expression pattern analysis of the GRF transcription factor family in Astragalus mongholicus

BACKGROUND

Astragalus membranaceus is a plant of the Astragalus genus, which is used as a traditional Chinese herbal medicine with extremely high medicinal and edible value. Astragalus mongholicus, as one of the representative medicinal materials with the same origin of medicine and food, has a rising market demand for its raw materials, but the quality is different in different production areas. Growth-regulating factors (GRF) are transcription factors unique to plants that play important roles in plant growth and development. Up to now, there is no report about GRF in A. mongholicus.

METHODS AND RESULTS

This study conducted a genome-wide analysis of the AmGRF gene family, identifying a total of nine AmGRF genes that were classified into subfamily V based on phylogenetic relationships. In the promoter region of the AmGRF gene, we successfully predicted cis-elements that respond to abiotic stress, growth, development, and hormone production in plants. Based on transcriptomic data and real-time quantitative polymerase chain reaction (qPCR) validation, the results showed that AmGRFs were expressed in the roots, stems, and leaves, with overall higher expression in leaves, higher expression of AmGRF1 and AmGRF8 in roots, and high expression levels of AmGRF1 and AmGRF9 in stems.

CONCLUSIONS

The results of this study provide a theoretical basis for the further exploration of the functions of AmGRFs in plant growth and development.

Application of Developmental Regulators for Enhancing Plant Regeneration and Genetic Transformation

Establishing plant regeneration systems and efficient genetic transformation techniques plays a crucial role in plant functional genomics research and the development of new crop varieties. The inefficient methods of transformation and regeneration of recalcitrant species and the genetic dependence of the transformation process remain major obstacles. With the advancement of plant meristematic tissues and somatic embryogenesis research, several key regulatory genes, collectively known as developmental regulators, have been identified. In the field of plant genetic transformation, the application of developmental regulators has recently garnered significant interest. These regulators play important roles in plant growth and development, and when applied in plant genetic transformation, they can effectively enhance the induction and regeneration capabilities of plant meristematic tissues, thus providing important opportunities for improving genetic transformation efficiency. This review focuses on the introduction of several commonly used developmental regulators. By gaining an in-depth understanding of and applying these developmental regulators, it is possible to further enhance the efficiency and success rate of plant genetic transformation, providing strong support for plant breeding and genetic engineering research.

Temporal dynamics of genetic architecture governing leaf development in Populus

Leaf development is a multifaceted and dynamic process orchestrated by a myriad of genes to shape the proper size and morphology. The dynamic genetic network underlying leaf development remains largely unknown. Utilizing a synergistic genetic approach encompassing dynamic genome-wide association study (GWAS), time-ordered gene co-expression network (TO-GCN) analyses and gene manipulation, we explored the temporal genetic architecture and regulatory network governing leaf development in Populus. We identified 42 time-specific and 18 consecutive genes that displayed different patterns of expression at various time points. We then constructed eight TO-GCNs that covered the cell proliferation, transition, and cell expansion stages of leaf development. Integrating GWAS and TO-GCN, we postulated the functions of 27 causative genes for GWAS and identified PtoGRF9 as a key player in leaf development. Genetic manipulation via overexpression and suppression of PtoGRF9 revealed its primary influence on leaf development by modulating cell proliferation. Furthermore, we elucidated that PtoGRF9 governs leaf development by activating PtoHB21 during the cell proliferation stage and attenuating PtoLD during the transition stage. Our study provides insights into the dynamic genetic underpinnings of leaf development and understanding the regulatory mechanism of PtoGRF9 in this dynamic process.

Genome-Wide Identification and Expression Analysis of Growth-Regulating Factors in Eucommia ulmoides Oliver (Du-Zhong)

The roots, stems, leaves, and seeds of Eucommia ulmoides contain a large amount of trans-polyisoprene (also known as Eu-rubber), which is considered to be an important laticiferous plant with valuable industrial applications. Eu-rubber used in industry is mainly extracted from leaves. Therefore, it is of great significance to identify genes related to regulating the leaf size of E. ulmoides. Plant growth-regulating factors (GRFs) play important roles in regulating leaf size, and their functions are highly conserved across different plant species. However, there have been very limited reports on EuGRFs until now. In this study, eight canonical EuGRFs with both QLQ and WRC domains and two putative eul-miR396s were identified in the chromosome-level genome of E. ulmoides. It is found that, unlike AtGRFs, all EuGRFs contain the miR396s binding site in the terminal of WRC domains. These EuGRFs were distributed on six chromosomes in the genome of E. ulmoides. Collinearity analysis of the E. ulmoides genome revealed that EuGRF1 and EuGRF3 exhibit collinear relationships with EuGRF2, suggesting that those three genes may have emerged via gene replication events. The collinear relationship between EuGRFs, AtGRFs, and OsGRFs showed that EuGRF5 and EuGRF8 had no collinear members in Arabidopsis and rice. Almost all EuGRFs show a higher expression level in growing and developing tissues, and most EuGRF promoters process phytohormone-response and stress-induced cis-elements. Moreover, we found the expression of EuGRFs was significantly induced by gibberellins (GA3) in three hours, and the height of E. ulmoides seedlings was significantly increased one week after GA3 treatment. The findings in this study provide potential candidate genes for further research and lay the foundation for further exploring the molecular mechanism underlying E. ulmoides development in response to GA3.

Overexpression of HbGRF4 or HbGRF4-HbGIF1 Chimera Improves the Efficiency of Somatic Embryogenesis in Hevea brasiliensis

Transgenic technology is a crucial tool for gene functional analysis and targeted genetic modification in the para rubber tree (Hevea brasiliensis). However, low efficiency of plant regeneration via somatic embryogenesis remains a bottleneck of successful genetic transformation in H. brasiliensis. Enhancing expression of GROWTH-REGULATING FACTOR 4 (GRF4)-GRF-INTERACTING FACTOR 1 (GIF1) has been reported to significantly improve shoot and embryo regeneration in multiple crops. Here, we identified endogenous HbGRF4 and HbGIF1 from the rubber clone Reyan7-33-97, the expressions of which dramatically increased along with somatic embryo (SE) production. Intriguingly, overexpression of HbGRF4 or HbGRF4-HbGIF1 markedly enhanced the efficiency of embryogenesis in two H. brasiliensis callus lines with contrasting rates of SE production. Transcriptional profiling revealed that the genes involved in jasmonic acid response were up-regulated, whereas those in ethylene biosynthesis and response as well as the S-adenosylmethionine-dependent methyltransferase activity were down-regulated in HbGRF4- and HbGRF4-HbGIF1-overexpressing H. brasiliensis embryos. These findings open up a new avenue for improving SE production in rubber tree, and help to unravel the underlying mechanisms of HbGRF4-enhanced somatic embryogenesis.

Common signatures of selection reveal target loci for breeding across soybean populations

Understanding the underlying genetic bases of yield-related selection and distinguishing these changes from genetic drift are critical for both improved understanding and future success of plant breeding. Soybean [Glycine max (L.) Merr.] is a key species for world food security, yet knowledge of the mechanism of selective breeding in soybean, such as the century-long program of artificial selection in U.S. soybean germplasm, is currently limited to certain genes and loci. Here, we identify genome-wide signatures of selection in separate populations of soybean subjected to artificial selection for increased yield by multiple breeding programs in the United States. We compared the alternative soybean breeding population (AGP) created by USDA-ARS to the conventional public soybean lines (CGP) developed at three different stages of breeding (ancestral, intermediate, and elite) to identify shared signatures of selection and differentiate these from drift. The results showed a strong selection for specific haplotypes identified by single site frequency and haplotype homozygosity methods. A set of common selection signatures was identified in both AGP and CGP that supports the hypothesis that separate breeding programs within similar environments coalesce on the fixation of the same key haplotypes. Signatures unique to each breeding program were observed. These results raise the possibility that selection analysis can allow the identification of favorable alleles to enhance directed breeding approaches.

Maize-Azospirillum brasilense interaction: accessing maize's miRNA expression under the effect of an inhibitor of indole-3-acetic acid production by the plant

MicroRNA (miRNA) is a class of non-coding RNAs. They play essential roles in plants' physiology, as in the regulation of plant development, response to biotic and abiotic stresses, and symbiotic processes. This work aimed to better understand the importance of maize's miRNA during Azospirillum-plant interaction when the plant indole-3-acetic acid (IAA) production was inhibited with yucasin, an inhibitor of the TAM/YUC pathway. Twelve cDNA libraries from a previous Dual RNA-Seq experiment were used to analyze gene expression using a combined analysis approach. miRNA coding genes (miR) and their predicted mRNA targets were identified among the differentially expressed genes. Statistical differences among the groups indicate that Azospirillum brasilense, yucasin, IAA concentration, or all together could influence the expression of several maize's miRNAs. The miRNA's probable targets were identified, and some of them were observed to be differentially expressed. Dcl4, myb122, myb22, and morf3 mRNAs were probably regulated by their respective miRNAs. Other probable targets were observed responding to the IAA level, the bacterium, or all of them. A. brasilense was able to influence the expression of some maize's miRNA, for example, miR159f, miR164a, miR169j, miR396c, and miR399c. The results allow us to conclude that the bacterium can influence directly or indirectly the expression of some of the identified mRNA targets, probably due to an IAA-independent pathway, and that they are somehow involved in the previously observed physiological effects.

A rapid and highly efficient sorghum transformation strategy using GRF4-GIF1/ternary vector system

Sorghum is an important crop for food, forage, wine and biofuel production. To enhance its transformation efficiency without negative developmental by-effects, we investigated the impact of GRF4-GIF1 chimaera and GRF5 on sorghum transformation. Both GRF4-GIF1 and GRF5 effectively improved the transformation efficiency of sorghum and accelerated the transformation process of sorghum to less than 2 months which was not observed when using BBM-WUS. As agrobacterium  effectors increase the ability of T-DNA transfer into plant cells, we checked whether ternary vector system can additively enhance sorghum transformation. The combination of GRF4-GIF1 with helper plasmid pVS1-VIR2 achieved the highest transformation efficiency, reaching 38.28%, which is 7.71-fold of the original method. Compared with BBM-WUS, overexpressing GRF4-GIF1 caused no noticeable growth defects in sorghum. We further developed a sorghum CRISPR/Cas9 gene-editing tool based on this GRF4-GIF1/ternary vector system, which achieved an average gene mutation efficiency of 41.36%, and null mutants were created in the T0 generation.

Growth-regulating factor 15-mediated vascular cambium differentiation positively regulates wood formation in hybrid poplar (Populus alba × P. glandulosa)

Introduction

Hybrid poplars are industrial trees in China. An understanding of the molecular mechanism underlying wood formation in hybrid poplars is necessary for molecular breeding. Although the division and differentiation of vascular cambial cells is important for secondary growth and wood formation, the regulation of this process is largely unclear.

Methods

In this study, mPagGRF15 OE and PagGRF15-SRDX transgenic poplars were generated to investigate the function of PagGRF15. RNA-seq and qRT-PCR were conducted to analyze genome-wide gene expression, while ChIP‒seq and ChIP-PCR were used to identified the downstream genes regulated by PagGRF15.

Results and discussion

We report that PagGRF15 from hybrid poplar (Populus alba × P. glandulosa), a growth-regulating factor, plays a critical role in the regulation of vascular cambium activity. PagGRF15 was expressed predominantly in the cambial zone of vascular tissue. Overexpression of mPagGRF15 (the mutated version of GRF15 in the miR396 target sequence) in Populus led to decreased plant height and internode number. Further stem cross sections showed that the mPagGRF15 OE plants exhibited significant changes in vascular pattern with an increase in xylem and a reduction in phloem. In addition, cambium cell files were decreased in the mPagGRF15 OE plants. However, dominant suppression of the downstream genes of PagGRF15 using PagGRF15-SRDX showed an opposite phenotype. Based on the RNA-seq and ChIP-seq results, combining qRT-PCR and ChIP-PCR analysis, candidate genes, such as WOX4b, PXY and GID1.3, were obtained and found to be mainly involved in cambial activity and xylem differentiation. Accordingly, we speculated that PagGRF15 functions as a positive regulator mediating xylem differentiation by repressing the expression of the WOX4a and PXY genes to set the pace of cambial activity. In contrast, PagGRF15 mediated the GA signaling pathway by upregulating GID1.3 expression to stimulate xylem differentiation. This study provides valuable information for further studies on vascular cambium differentiation mechanisms and genetic improvement of the specific gravity of wood in hybrid poplars.

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

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

miRNAs and genes as molecular regulators of rice grain morphology and yield

Rice is one of the most consumed crops worldwide and the genetic and molecular basis of its grain yield attributes are well understood. Various studies have identified different yield-related parameters in rice that are regulated by the microRNAs (miRNAs). MiRNAs are endogenous small non-coding RNAs that silence gene expression during or after transcription. They control a variety of biological or genetic activities in plants including growth, development and response to stress. In this review, we have summarized the available information on the genetic control of panicle architecture and grain yield (number and morphology) in rice. The miRNA nodes that are associated with their regulation are also described while focussing on the central role of miR156-SPL node to highlight the co-regulation of two master regulators that determine the fate of panicle development. Since abiotic stresses are known to negatively affect yield, the impact of abiotic stress induced alterations on the levels of these miRNAs are also discussed to highlight the potential of miRNAs for regulating crop yields.

Emerging roles of plant microRNAs during Colletotrichum spp. infection

MAIN CONCLUSION

This review provides valuable insights into plant molecular regulatory mechanisms during fungus attacks, highlighting potential miRNA candidates for future disease management. Plant defense responses to biotic stress involve intricate regulatory mechanisms, including post-transcriptional regulation of genes mediated by microRNAs (miRNAs). These small RNAs play a vital role in the plant's innate immune system, defending against viral, bacterial, and fungal attacks. Among the plant pathogenic fungi, Colletotrichum spp. are notorious for causing anthracnose, a devastating disease affecting economically important crops worldwide. Understanding the molecular machinery underlying the plant immune response to Colletotrichum spp. is crucial for developing tools to reduce production losses. In this comprehensive review, we examine the current understanding of miRNAs associated with plant defense against Colletotrichum spp. We summarize the modulation patterns of miRNAs and their respective target genes. Depending on the function of their targets, miRNAs can either contribute to host resistance or susceptibility. We explore the multifaceted roles of miRNAs during Colletotrichum infection, including their involvement in R-gene-dependent immune system responses, hormone-dependent defense mechanisms, secondary metabolic pathways, methylation regulation, and biosynthesis of other classes of small RNAs. Furthermore, we employ an integrative approach to correlate the identified miRNAs with various strategies and distinct phases of fungal infection. This study provides valuable insights into the current understanding of plant miRNAs and their regulatory mechanisms during fungus attacks.

Genome-wide molecular evolution analysis of the GRF and GIF gene families in Plantae (Archaeplastida)

BACKGROUND

Plant growth-regulating factors (GRFs) and GRF-interacting factors (GIFs) interact with each other and collectively have important regulatory roles in plant growth, development, and stress responses. Therefore, it is of great significance to explore the systematic evolution of GRF and GIF gene families. However, our knowledge and understanding of the role of GRF and GIF genes during plant evolution has been fragmentary.

RESULTS

In this study, a large number of genomic and transcriptomic datasets of algae, mosses, ferns, gymnosperms and angiosperms were used to systematically analyze the evolution of GRF and GIF genes during the evolution of plants. The results showed that GRF gene first appeared in the charophyte Klebsormidium nitens, whereas the GIF genes originated relatively early, and these two gene families were mainly expanded by segmental duplication events after plant terrestrialization. During the process of evolution, the protein sequences and functions of GRF and GIF family genes are relatively conservative. As cooperative partner, GRF and GIF genes contain the similar types of cis-acting elements in their promoter regions, which enables them to have similar transcriptional response patterns, and both show higher levels of expression in reproductive organs and tissues and organs with strong capacity for cell division. Based on protein-protein interaction analysis and verification, we found that the GRF-GIF protein partnership began to be established in pteridophytes and is highly conserved across different terrestrial plants.

CONCLUSIONS

These results provide a foundation for further exploration of the molecular evolution and biological functions of GRF and GIF genes.

AsHSP26.2, a creeping bentgrass chloroplast small heat shock protein positively regulates plant development

KEY MESSAGE

The creeping bentgrass small heat shock protein AsHSP26.2 positively regulates plant growth and is a novel candidate for use in crop genetic engineering for enhanced biomass production and grain yield. Small heat shock proteins (sHSPs), a family of proteins with high level of diversity, significantly influence plant stress tolerance and plant development. We have cloned a creeping bentgrass chloroplast-localized sHSP gene, AsHSP26.2 responsive to IAA, GA and 6-BA stimulation. Transgenic creeping bentgrass overexpressing AsHSP26.2 exhibited significantly enhanced plant growth with increased stolon number and length as well as enlarged leaf blade width and leaf sheath diameters, but inhibited leaf trichomes initiation and development in the abaxial epidermis. These phenotypes are completely opposite to those displayed in the transgenic plants overexpressing AsHSP26.8, another chloroplast sHSP26 isoform that contains additional seven amino acids (AEGQGDG) between the consensus regions III and IV (Sun et al., Plant Cell Environ 44:1769-1787, 2021). Furthermore, AsHSP26.2 overexpression altered phytohormone biosynthesis and signaling transduction, resulting in elevated auxin and gibberellins (GA) accumulation. The results obtained provide novel insights implicating the sHSPs in plant growth and development regulation, and strongly suggest AsHSP26.2 to be a novel candidate for use in crop genetic engineering for enhanced plant biomass production and grain yield.

Genome-wide identification of the GRF family in sweet orange (Citrus sinensis) and functional analysis of the CsGRF04 in response to multiple abiotic stresses

BACKGROUND

Citrus is one of the most valuable fruits worldwide and an economic pillar industry in southern China. Nevertheless, it frequently suffers from undesirable environmental stresses during the growth cycle, which severely restricts the growth, development and yield of citrus. In plants, the growth-regulating factor (GRF) family of transcription factors (TF) is extensively distributed and plays an vital part in plant growth and development, hormone response, as well as stress adaptation. However, the systematic identification and functional analysis of GRF TFs in citrus have not been reported.

RESULTS

Here, a genome-wide identification of GRF TFs was performed in Citrus sinensis, 9 members of CsGRFs were systematically identified and discovered to be scattered throughout 5 chromosomes. Subsequently, physical and chemical properties, phylogenetic relationships, structural characteristics, gene duplication events, collinearity and cis-elements of promoter were elaborately analyzed. In particular, the expression patterns of the CsGRF genes in response to multiple phytohormone and abiotic stress treatments were investigated. Predicated on this result, CsGRF04, which exhibited the most differential expression pattern under multiple phytohormone and abiotic stress treatments was screened out. Virus-induced gene silencing (VIGS) technology was utilized to obtain gene silenced plants for CsGRF04 successfully. After the three stress treatments of high salinity, low temperature and drought, the CsGRF04-VIGS lines showed significantly reduced resistance to high salinity and low temperature stresses, but extremely increased resistance to drought stress.

CONCLUSIONS

Taken together, our findings systematically analyzed the genomic characterization of GRF family in Citrus sinensis, and excavated a CsGRF04 with potential functions under multiple abiotic stresses. Our study lay a foundation for further study on the function of CsGRFs in abiotic stress and hormone signaling response.

Multifaceted roles of transcription factors during plant embryogenesis

Transcription factors (TFs) are diverse groups of regulatory proteins. Through their specific binding domains, TFs bind to their target genes and regulate their expression, therefore TFs play important roles in various growth and developmental processes. Plant embryogenesis is a highly regulated and intricate process during which embryos arise from various sources and undergo development; it can be further divided into zygotic embryogenesis (ZE) and somatic embryogenesis (SE). TFs play a crucial role in the process of plant embryogenesis with a number of them acting as master regulators in both ZE and SE. In this review, we focus on the master TFs involved in embryogenesis such as BABY BOOM (BBM) from the APETALA2/Ethylene-Responsive Factor (AP2/ERF) family, WUSCHEL and WUSCHEL-related homeobox (WOX) from the homeobox family, LEAFY COTYLEDON 2 (LEC2) from the B3 family, AGAMOUS-Like 15 (AGL15) from the MADS family and LEAFY COTYLEDON 1 (LEC1) from the Nuclear Factor Y (NF-Y) family. We aim to present the recent progress pertaining to the diverse roles these master TFs play in both ZE and SE in Arabidopsis, as well as other plant species including crops. We also discuss future perspectives in this context.

Plant-TFClass: a structural classification for plant transcription factors

Transcription factors (TFs) bind DNA at specific sequences to regulate gene expression. This universal process is achieved via their DNA-binding domain (DBD). In mammals, the vast diversity of DBD structural conformations and the way in which they contact DNA has been used to organize TFs in the TFClass hierarchical classification. However, the numerous DBD types present in plants but absent from mammalian genomes were missing from this classification. We reviewed DBD 3D structures and models available for plant TFs to classify most of the 56 recognized plant TF types within the TFClass framework. This extended classification adds eight new classes and 37 new families corresponding to DBD structures absent in mammals. Plant-TFClass provides a unique resource for TF comparison across families and organisms.

Growing Conditions and Varietal Ecologies Differently Regulates the Growth-regulating-factor (GRFs) Gene Family in Rice

Background

Growth-regulating factors (GRFs) are crucial in rice for controlling plant growth and development. Among the rice cultivation practices, aerobic methods are water efficient but result in significant yield reduction relative to non-aerobic cultivation. Therefore, mechanistic insights into aerobic rice cultivation are important for improving the aerobic performance of rice.

Objectives

This study aimed to examine the evolution of GRFs in different rice species, analyse the phenotypic differences between aerobic and non-aerobic conditions in three rice varieties, and assess the expression of GRFs in these varieties under both aerobic and non-aerobic conditions.

Materials and Methods

This study comprehensively examined the GRFs gene family in 11 rice species (Oryza barthii, Oryza brachyantha, Oryza glaberrima, Oryza glumipatula, Oryza sativa subsp. indica, Oryza longistaminata, Oryza meridionalis, Oryza nivara, Oryza punctata, Oryza rufipogon, Oryza sativa subsp. japonica) focusing on phylogenetic analysis. Additionally, the expression patterns of 12 GRFs were investigated in three distinct genotypes of O. sativa subsp. indica rice, under both non-aerobic and aerobic conditions.

Results

Three major phylogenetic clades were formed based on conserved motifs in the 123 GRFs proteins in eleven rice species. Further, novel motifs were identified especially in O. longistaminata indicative of the species level evolutionary differences in rice. Among the trait performance, the number of tillers was reduced by ~ 36% under aerobic conditions, but the reduction was found to be less in CR Dhan 201, an aerobic variety. Besides, three GRFs namely GRF3, GRF4, and GRF7 were found to be distinct in expression between aerobic and non-aerobic conditions.

Conclusion

Three GRF genes namely GRF3, GRF4, and GRF7 could be associated with the aerobic adaptation in rice.

Genome-wide identification and expression analysis of growth-regulating factors in Dendrobium officinale and Dendrobium chrysotoxum

Background

Dendrobium, one of the largest genera in Orchidaceae, is popular not only for its aesthetic appeal but for its significant medicinal value. Growth-regulating factors (GRFs) play an essential role in plant growth and development. However, there is still a lack of information about the evolution and biological function analysis of the GRF gene family among Dendrobiumspecies.

Methods

Growth-regulating factors from Dendrobium officinale Kimura et Migo and Dendrobium chrysotoxum Lindl. were identified by HMMER and BLAST. Detailed bioinformatics analysis was conducted to explore the evolution and function of GRF gene family in D. officinale and D. chrysotoxum using genomic data, transcriptome data and qRT-PCR technology.

Results

Here, we evaluated the evolution of the GRF gene family based on the genome sequences of D. officinale and D. chrysotoxum. Inferred from phylogenetic trees, the GRF genes were classified into two clades, and each clade contains three subclades. Sequence comparison analysis revealed relatively conserved gene structures and motifs among members of the same subfamily, indicating a conserved evolution of GRF genes within Dendrobiumspecies. However, considering the distribution of orthologous DoGRFs and DcGRFs, and the differences in the number of GRFs among species, we suggest that the GRF gene family has undergone different evolutionary processes. A total of 361 cis-elements were detected, with 33, 141, and 187 related to plant growth and development, stress, and hormones, respectively. The tissue-specific expression of GRFs showed that DoGRF8 may have a significant function in the stem elongation of D. officinale. Moreover, four genes were up-regulated under Methyl-jasmonic acid/methyl jasmonate (MeJA) treatment, showing that DoGRFs and DcGRFs play a crucial role in stress response. These findings provide valuable information for further investigations into the evolution and function of GRF genes in D. officinale and D. chrysotoxum.

Genome-Wide Analysis of Flax (Linum usitatissimum L.) Growth-Regulating Factor (GRF) Transcription Factors

Flax is an important cash crop globally with a variety of commercial uses. It has been widely used for fiber, oil, nutrition, feed and in composite materials. Growth regulatory factor (GRF) is a transcription factor family unique to plants, and is involved in regulating many processes of growth and development. Bioinformatics analysis of the GRF family in flax predicted 17 LuGRF genes, which all contained the characteristic QLQ and WRC domains. Equally, 15 of 17 LuGRFs (88%) are predicted to be regulated by lus-miR396 miRNA. Phylogenetic analysis of GRFs from flax and several other well-characterized species defined five clades; LuGRF genes were found in four clades. Most LuGRF gene promoters contained cis-regulatory elements known to be responsive to hormones and stress. The chromosomal locations and collinearity of LuGRF genes were also analyzed. The three-dimensional structure of LuGRF proteins was predicted using homology modeling. The transcript expression data indicated that most LuGRF family members were highly expressed in flax fruit and embryos, whereas LuGRF3, LuGRF12 and LuGRF16 were enriched in response to salt stress. Real-time quantitative fluorescent PCR (qRT-PCR) showed that both LuGRF1 and LuGRF11 were up-regulated under ABA and MeJA stimuli, indicating that these genes were involved in defense. LuGRF1 was demonstrated to be localized to the nucleus as expected for a transcription factor. These results provide a basis for further exploration of the molecular mechanism of LuGRF gene function and obtaining improved flax breeding lines.

Design principles for synthetic control systems to engineer plants

KEY MESSAGE

Synthetic control systems have led to significant advancement in the study and engineering of unicellular organisms, but it has been challenging to apply these tools to multicellular organisms like plants. The ability to predictably engineer plants will enable the development of novel traits capable of alleviating global problems, such as climate change and food insecurity. Engineering predictable multicellular phenotypes will require the development of synthetic control systems that can precisely regulate how the information encoded in genomes is translated into phenotypes. Many efficient control systems have been developed for unicellular organisms. However, it remains challenging to use such tools to study or engineer multicellular organisms. Plants are a good chassis within which to develop strategies to overcome these challenges, thanks to their capacity to withstand large-scale reprogramming without lethality. Additionally, engineered plants have great potential for solving major societal problems. Here we briefly review the progress of control system development in unicellular organisms, and how that information can be leveraged to characterize control systems in plants. Further, we discuss strategies for developing control systems designed to regulate the expression of transgenes or endogenous loci and generate dosage-dependent or discrete traits. Finally, we discuss the utility that mathematical models of biological processes have for control system deployment.

Genetic insights into the dissolution of dioecy in diploid persimmon Diospyros oleifera Cheng

BACKGROUND

Dioecy, a sexual system of single-sexual (gynoecious/androecious) individuals, is rare in flowering plants. This rarity may be a result of the frequent transition from dioecy into systems with co-sexual individuals.

RESULTS

In this study, co-sexual expression (monoecy and hermaphroditic development), previously thought to be polyploid-specific in Diospyros species, was identified in the diploid D. oleifeara historically. We characterized potential genetic mechanisms that underlie the dissolution of dioecy to monoecy and andro(gyno)monoecy, based on multiscale genome-wide investigations of 150 accessions of Diospyros oleifera. We found all co-sexual plants, including monoecious and andro(gyno)monoecious individuals, possessed the male determinant gene OGI, implying the presence of genetic factors controlling gynoecia development in genetically male D. oleifera. Importantly, discrepancies in the OGI/MeGI module were found in diploid monoecious D. oleifera compared with polyploid monoecious D. kaki, including no Kali insertion on the promoter of OGI, no different abundance of smRNAs targeting MeGI (a counterpart of OGI), and no different expression of MeGI between female and male floral buds. On the contrary, in both single- and co-sexual plants, female function was expressed in the presence of a genome-wide decrease in methylation levels, along with sexually distinct regulatory networks of smRNAs and their targets. Furthermore, a genome-wide association study (GWAS) identified a genomic region and a DUF247 gene cluster strongly associated with the monoecious phenotype and several regions that may contribute to andromonoecy.

CONCLUSIONS

Collectively, our findings demonstrate stable breakdown of the dioecious system in D. oleifera, presumably also a result of genomic features of the Y-linked region.

Genetic Basis of Grain Size and Weight in Rice, Wheat, and Barley

Grain size is a key component of grain yield in cereals. It is a complex quantitative trait controlled by multiple genes. Grain size is determined via several factors in different plant development stages, beginning with early tillering, spikelet formation, and assimilates accumulation during the pre-anthesis phase, up to grain filling and maturation. Understanding the genetic and molecular mechanisms that control grain size is a prerequisite for improving grain yield potential. The last decade has brought significant progress in genomic studies of grain size control. Several genes underlying grain size and weight were identified and characterized in rice, which is a model plant for cereal crops. A molecular function analysis revealed most genes are involved in different cell signaling pathways, including phytohormone signaling, transcriptional regulation, ubiquitin-proteasome pathway, and other physiological processes. Compared to rice, the genetic background of grain size in other important cereal crops, such as wheat and barley, remains largely unexplored. However, the high level of conservation of genomic structure and sequences between closely related cereal crops should facilitate the identification of functional orthologs in other species. This review provides a comprehensive overview of the genetic and molecular bases of grain size and weight in wheat, barley, and rice, focusing on the latest discoveries in the field. We also present possibly the most updated list of experimentally validated genes that have a strong effect on grain size and discuss their molecular function.

Specificity Testing for NGT PCR-Based Detection Methods in the Context of the EU GMO Regulations

The term new genomic techniques (NGTs) is an umbrella term used to describe a variety of techniques that can alter the genetic material of an organism and that have emerged or have been developed since 2001, when the existing genetically modified organism (GMO) legislation was adopted. The analytical framework used to detect GMOs in Europe is an established single harmonized procedure that is mandatory for the authorization of GM food and feed, thus generating a reliable, transparent, and effective labeling scheme for GMO products. However, NGT products can challenge the implementation and enforcement of the current regulatory system in the EU, relating in particular to the detection of NGT products that contain no foreign genetic material. Consequently, the current detection methods might fail to meet the minimum performance requirements. Although existing detection methods may be able to detect and quantify even small alterations in the genome, this does not necessarily confirm the distinction between products resulting from NGTs subject to the GMO legislation and other products. Therefore, this study provides a stepwise approach for the in silico prediction of PCR systems' specificity by testing a bioinformatics pipeline for amplicon and primer set searches in current genomic databases. In addition, it also empirically tested the PCR system evaluated during the in silico analysis. Two mutant genotypes produced by CRISPR-Cas9 in Arabidopsis thaliana were used as a case study. Overall, our results demonstrate that the single PCR system developed for identifying a nucleotide insertion in the grf1-3 genotype has multiple matches in the databases, which do not enable the discrimination of this mutated event. Empirical assays further support this demonstration. In contrast, the second mutated genotype, grf8-61, which contains a -3 bp deletion, did not yield any matches in the sequence variant database. However, the primer sequences were not efficient during the empirical assay. Our approach represents a first step in decision making for analytical methods for NGT detection, identification, and quantification in light of the European labeling regulations.

Molecular Mechanisms of the miR396b-GRF1 Module Underlying Rooting Regulation in Acer rubrum L

Rooting and root development in Acer rubrum have important effects on overall growth. A. rubrum does not take root easily in natural conditions. In this study, the mechanisms of the miR396b-GRF1 module underlying rooting regulation in A. rubrum were studied. The subcellular localization and transcriptional activation of miR396b and its target gene growth regulating factor 1 (GRF1) were investigated. These experiments showed that GRF1 was localized in the nucleus and had transcriptional activation activity. Functional validation experiments in transgenic plants demonstrated that overexpression of Ar-miR396b inhibited adventitious root growth, whereas overexpression of ArGRF1 increased adventitious root growth. These results help clarify the molecular regulatory mechanisms underlying adventitious root growth in A. rubrum and provide some new insights into the rooting rate in this species.

CRISPR/Cas9 mutagenesis of the Arabidopsis GROWTH-REGULATING FACTOR (GRF) gene family

Genome editing in plants typically relies on T-DNA plasmids that are mobilized by Agrobacterium-mediated transformation to deliver the CRISPR/Cas machinery. Here, we introduce a series of CRISPR/Cas9 T-DNA vectors for minimal settings, such as teaching labs. Gene-specific targeting sequences can be inserted as annealed short oligonucleotides in a single straightforward cloning step. Fluorescent markers expressed in mature seeds enable reliable selection of transgenic or transgene-free individuals using a combination of inexpensive LED lamps and colored-glass alternative filters. Testing these tools on the Arabidopsis GROWTH-REGULATING FACTOR (GRF) genes, we were able to create a collection of predicted null mutations in all nine family members with little effort. We then explored the effects of simultaneously targeting two, four and eight GRF genes on the rate of induced mutations at each target locus. In our hands, multiplexing was associated with pronounced disparities: while mutation rates at some loci remained consistently high, mutation rates at other loci dropped dramatically with increasing number of single guide RNA species, thereby preventing a systematic mutagenesis of the family.

Dynamic roles of small RNAs and DNA methylation associated with heterosis in allotetraploid cotton (Gossypium hirsutum L.)

BACKGROUND

Heterosis is a complex phenomenon wherein the hybrids outperform their parents. Understanding the underlying molecular mechanism by which hybridization leads to higher yields in allopolyploid cotton is critical for effective breeding programs. Here, we integrated DNA methylation, transcriptomes, and small RNA profiles to comprehend the genetic and molecular basis of heterosis in allopolyploid cotton at three developmental stages.

RESULTS

Transcriptome analysis revealed that numerous DEGs responsive to phytohormones (auxin and salicylic acid) were drastically altered in F1 hybrid compared to the parental lines. DEGs involved in energy metabolism and plant growth were upregulated, whereas DEGs related to basal defense were downregulated. Differences in homoeologous gene expression in F1 hybrid were greatly reduced after hybridization, suggesting that higher levels of parental expression have a vital role in heterosis. Small RNAome and methylome studies showed that the degree of DNA methylation in hybrid is higher when compared to the parents. A substantial number of allele-specific expression genes were found to be strongly regulated by CG allele-specific methylation levels. The hybrid exhibited higher 24-nt-small RNA (siRNA) expression levels than the parents. The regions in the genome with increased levels of 24-nt-siRNA were chiefly related to genes and their flanking regulatory regions, demonstrating a possible effect of these molecules on gene expression. The transposable elements correlated with siRNA clusters in the F1 hybrid had higher methylation levels but lower expression levels, which suggest that these non-additively expressed siRNA clusters, reduced the activity of transposable elements through DNA methylation in the hybrid.

CONCLUSIONS

These multi-omics data provide insights into how changes in epigenetic mechanisms and gene expression patterns can lead to heterosis in allopolyploid cotton. This makes heterosis a viable tool in cotton breeding.

A multi-level approach reveals key physiological and molecular traits in the response of two rice genotypes subjected to water deficit at the reproductive stage

Rice is more vulnerable to drought than maize, wheat, and sorghum because its water requirements remain high throughout the rice life cycle. The effects of drought vary depending on the timing, intensity, and duration of the events, as well as on the rice genotype and developmental stage. It can affect all levels of organization, from genes to the cells, tissues, and/or organs. In this study, a moderate water deficit was applied to two contrasting rice genotypes, IAC 25 and CIRAD 409, during their reproductive stage. Multi-level transcriptomic, metabolomic, physiological, and morphological analyses were performed to investigate the complex traits involved in their response to drought. Weighted gene network correlation analysis was used to identify the specific molecular mechanisms regulated by each genotype, and the correlations between gene networks and phenotypic traits. A holistic analysis of all the data provided a deeper understanding of the specific mechanisms regulated by each genotype, and enabled the identification of gene markers. Under non-limiting water conditions, CIRAD 409 had a denser shoot, but shoot growth was slower despite better photosynthetic performance. Under water deficit, CIRAD 409 was weakly affected regardless of the plant level analyzed. In contrast, IAC 25 had reduced growth and reproductive development. It regulated transcriptomic and metabolic activities at a high level, and activated a complex gene regulatory network involved in growth-limiting processes. By comparing two contrasting genotypes, the present study identified the regulation of some fundamental processes and gene markers, that drive rice development, and influence its response to water deficit, in particular, the importance of the biosynthetic and regulatory pathways for cell wall metabolism. These key processes determine the biological and mechanical properties of the cell wall and thus influence plant development, organ expansion, and turgor maintenance under water deficit. Our results also question the genericity of the antagonism between morphogenesis and organogenesis observed in the two genotypes.

Spatiotemporal miRNA and transcriptomic network dynamically regulate the developmental and senescence processes of poplar leaves

Poplar is an important afforestation and urban greening species. Poplar leaf development occurs in stages, from young to mature and then from mature to senescent; these are accompanied by various phenotypic and physiological changes. However, the associated transcriptional regulatory network is relatively unexplored. We first used principal component analysis to classify poplar leaves at different leaf positions into two stages: developmental maturity (the stage of maximum photosynthetic capacity); and the stage when photosynthetic capacity started to decline and gradually changed to senescence. The two stages were then further subdivided into five intervals by gene expression clustering analysis: young leaves, the period of cell genesis and functional differentiation (L1); young leaves, the period of development and initial formation of photosynthetic capacity (L3-L7); the period of maximum photosynthetic capacity of functional leaves (L9-L13); the period of decreasing photosynthetic capacity of functional leaves (L15-L27); and the period of senescent leaves (L29). Using a weighted co-expression gene network analysis of regulatory genes, high-resolution spatiotemporal transcriptional regulatory networks were constructed to reveal the core regulators that regulate leaf development. Spatiotemporal transcriptome data of poplar leaves revealed dynamic changes in genes and miRNAs during leaf development and identified several core regulators of leaf development, such as GRF5 and MYB5. This in-depth analysis of transcriptional regulation during leaf development provides a theoretical basis for exploring the biological basis of the transcriptional regulation of leaf development and the molecular design of breeding for delaying leaf senescence.

Genome-wide identification of auxin-responsive microRNAs in the poplar stem

BACKGROUND

Wood (secondary xylem) of forests is a material of great economic importance. Wood development is strictly controlled by both the phytohormone auxin and microRNAs (miRNAs). Currently, the regulatory mechanisms underlying wood formation by auxin-associated miRNAs remain unclear.

OBJECTIVE

This report was designed to identify auxin-responsive miRNAs during wood formation.

METHODS

Morphological observation of wood development in the poplar stems was performed under the treatment of different concentrations (0 mg/L, CK; 5 mg/L, Low; 10 mg/L, High) of indol-3-butyric acid (IBA). Using a small RNA sequencing strategy, the effect of IBA treatment on miRNAs expression was genome-widely analyzed.

RESULTS

In this study, we found that wood development of poplar was promoted by low concentration of IBA treatment but inhibited by high concentration of IBA treatment. Stringent bioinformatic analysis led to identification of 118 known and 134 novel miRNAs candidates. Sixty-nine unique developmental-related miRNAs, corresponding to 269 target genes, exhibited specific expression patterns in response to auxin, as was consistent with the influence of auxin application on wood formation. Three novel miRNAs had the most number (≥ 9) of target genes, belonging to SPL, GRF and ARF gene families. The evolutionary relationships and tissue expression patterns of 41 SPL, GRF and ARF genes in poplar were thus analyzed. Of them, four representative members and corresponding miRNAs were confirmed using RT-qPCR.

CONCLUSIONS

Our results may be helpful for a better understanding of auxin-induced regulation of wood formation in tree species.

Comparatively Evolution and Expression Analysis of GRF Transcription Factor Genes in Seven Plant Species

Growth regulatory factors (GRF) are plant-specific transcription factors that play pivotal roles in growth and various abiotic stresses regulation. However, adaptive evolution of GRF gene family in land plants are still being elucidated. Here, we performed the evolutionary and expression analysis of GRF gene family from seven representative species. Extensive phylogenetic analyses and gene structure analysis revealed that the number of genes, QLQ domain and WRC domain identified in higher plants was significantly greater than those identified in lower plants. Besides, dispersed duplication and WGD/segmental duplication effectively promoted expansion of the GRF gene family. The expression patterns of GRF gene family and target genes were found in multiple floral organs and abundant in actively growing tissues. They were also found to be particularly expressed in response to various abiotic stresses, with stress-related elements in promoters, implying potential roles in floral development and abiotic stress. Our analysis in GRF gene family interaction network indicated the similar results that GRFs resist to abiotic stresses with the cooperation of other transcription factors like GIFs. This study provides insights into evolution in the GRF gene family, together with expression patterns valuable for future functional researches of plant abiotic stress biology.

Optimised Agrobacterium-Mediated Transformation and Application of Developmental Regulators Improve Regeneration Efficiency in Melons

Melon (Cucumis melo L.) is a protected crop in China with high economic value. Agrobacterium-mediated genetic transformation is a powerful tool to improve agronomic traits and obtain elite germplasm. However, current transformation protocols in melons are inefficient and highly genotype-dependent. To improve transformation in melon, we tested different infiltration methods for Agrobacterium-mediated transformation. Among these methods, micro-brushing and sonication for 20 s, followed by vacuum infiltration at -1.0 kPa for 90 s, resulted in the strongest green fluorescent protein signal and increased the proportion of infected explants. We transformed melon with developmental regulatory genes AtGRF5, AtPLT5, AtBBM, AtWUS, AtWOX5, and AtWIND1 from Arabidopsis and estimated regeneration frequencies as the number of regenerating shoots/total number of inoculated explants in the selection medium. The overexpression of AtGRF5 and AtPLT5 in melon resulted in transformation efficiencies of 42.3% and 33% in ZHF and 45.6% and 32.9% in Z12, respectively, which were significantly higher than those of the control. AtGRF5 and AtPLT5 expression cassettes were added to CRISPR/Cas9 genome-editing vectors to obtain transgenic phytoene desaturase CmPDS knockout mutants. Using AtGRF5 or AtPLT5, multi-allelic mutations were observed at CmPDS target sites in recalcitrant melon genotypes. This strategy enables genotype-flexible transformation and promotes precise genome modification technologies in melons.

Technological Development and Application of Plant Genetic Transformation

Genetic transformation is an important strategy for enhancing plant biomass or resistance in response to adverse environments and population growth by imparting desirable genetic characteristics. Research on plant genetic transformation technology can promote the functional analysis of plant genes, the utilization of excellent traits, and precise breeding. Various technologies of genetic transformation have been continuously discovered and developed for convenient manipulation and high efficiency, mainly involving the delivery of exogenous genes and regeneration of transformed plants. Here, currently developed genetic transformation technologies were expounded and compared. Agrobacterium-mediated gene delivery methods are commonly used as direct genetic transformation, as well as external force-mediated ways such as particle bombardment, electroporation, silicon carbide whiskers, and pollen tubes as indirect ones. The regeneration of transformed plants usually involves the de novo organogenesis or somatic embryogenesis pathway of the explants. Ectopic expression of morphogenetic transcription factors (Bbm, Wus2, and GRF-GIF) can significantly improve plant regeneration efficiency and enable the transformation of some hard-to-transform plant genotypes. Meanwhile, some limitations in these gene transfer methods were compared including genotype dependence, low transformation efficiency, and plant tissue damage, and recently developed flexible approaches for plant genotype transformation are discussed regarding how gene delivery and regeneration strategies can be optimized to overcome species and genotype dependence. This review summarizes the principles of various techniques for plant genetic transformation and discusses their application scope and limiting factors, which can provide a reference for plant transgenic breeding.

Whole-genome identification and expression profiling of growth-regulating factor (GRF) and GRF-interacting factor (GIF) gene families in Panax ginseng

BACKGROUND

Panax ginseng is a perennial herb and one of the most widely used traditional medicines in China. During its long growth period, it is affected by various environmental factors. Past studies have shown that growth-regulating factors (GRFs) and GRF-interacting factors (GIFs) are involved in regulating plant growth and development, responding to environmental stress, and responding to the induction of exogenous hormones. However, GRF and GIF transcription factors in ginseng have not been reported.

RESULTS

In this study, 20 GRF gene members of ginseng were systematically identified and found to be distributed on 13 chromosomes. The ginseng GIF gene family has only ten members, which are distributed on ten chromosomes. Phylogenetic analysis divided these PgGRFs into six clades and PgGIFs into two clades. In total, 18 of the 20 PgGRFs and eight of the ten PgGIFs are segmental duplications. Most PgGRF and PgGIF gene promoters contain some hormone- and stress- related cis-regulatory elements. Based on the available public RNA-Seq data, the expression patterns of PgGRF and PgGIF genes were analysed from 14 different tissues. The responses of the PgGRF gene to different hormones (6-BA, ABA, GA3, IAA) and abiotic stresses (cold, heat, drought, and salt) were studied. The expression of the PgGRF gene was significantly upregulated under GA3 induction and three weeks of heat treatment. The expression level of the PgGIF gene changed only slightly after one week of heat treatment.

CONCLUSIONS

The results of this study may be helpful for further study of the function of PgGRF and PgGIF genes and lay a foundation for further study of their role in the growth and development of Panax ginseng.

Properties of INDETERMINATE DOMAIN Proteins from Physcomitrium patens: DNA-Binding, Interaction with GRAS Proteins, and Transcriptional Activity

INDETERMINATE DOMAIN (IDD) proteins are plant-specific transcription factors that interact with GRAS proteins, such as DELLA and SHORT ROOT (SHR), to regulate target genes. The combination of IDD and DELLA proteins regulates genes involved in gibberellic acid (GA) synthesis and GA signaling, whereas the combination of IDD with the complex of SHR and SCARECROW, another GRAS protein, regulates genes involved in root tissue formation. Previous bioinformatic research identified seven IDDs, two DELLA, and two SHR genes in Physcomitrium patens, a model organism for non-vascular plants (bryophytes), which lack a GA signaling pathway and roots. In this study, DNA-binding properties and protein-protein interaction of IDDs from P. patens (PpIDD) were analyzed. Our results showed that the DNA-binding properties of PpIDDs were largely conserved between moss and seed plants. Four PpIDDs showed interaction with Arabidopsis DELLA (AtDELLA) proteins but not with PpDELLAs, and one PpIDD showed interaction with PpSHR but not with AtSHR. Moreover, AtIDD10 (JACKDAW) interacted with PpSHR but not with PpDELLAs. Our results indicate that DELLA proteins have modified their structure to interact with IDD proteins during evolution from moss lineage to seed plants, whereas the interaction of IDD and SHR was already present in moss lineage.

Characterization of HSP70 family in watermelon (Citrullus lanatus): identification, structure, evolution, and potential function in response to ABA, cold and drought stress

Watermelon (Citrullus lanatus) as a crop with important economic value, is widely cultivated around the world. The heat shock protein 70 (HSP70) family in plant is indispensable under stress conditions. However, no comprehensive analysis of watermelon HSP70 family is reported to date. In this study, 12 ClHSP70 genes were identified from watermelon, which were unevenly located in 7 out of 11 chromosomes and divided into three subfamilies. ClHSP70 proteins were predicted to be localized primarily in cytoplasm, chloroplast, and endoplasmic reticulum. Two pairs of segmental repeats and 1 pair of tandem repeats existed in ClHSP70 genes, and ClHSP70s underwent strong purification selection. There were many abscisic acid (ABA) and abiotic stress response elements in ClHSP70 promoters. Additionally, the transcriptional levels of ClHSP70s in roots, stems, true leaves, and cotyledons were also analyzed. Some of ClHSP70 genes were also strongly induced by ABA. Furthermore, ClHSP70s also had different degrees of response to drought and cold stress. The above data indicate that ClHSP70s may be participated in growth and development, signal transduction and abiotic stress response, laying a foundation for further analysis of the function of ClHSP70s in biological processes.

Screening Candidate Genes at the Co Locus Conferring to the Columnar Growth Habit in Apple (Malus × Domestica Borkh.)

The columnar growth trait of apple (Malus × domestica Borkh.) is genetically controlled by the Columnar (Co) locus on 10 chromosomes, including several candidate genes. Except for MdCo31, other candidate genes at the Co locus are less elucidated. In this study, a strategy of step-by-step screening was adopted to select 11 candidate genes by experimental cloning, transient expression, and genetic transformation. There existed several SNPs in four genes by sequence alignment in columnar and non-columnar apples. Two genes were detected in the nucleus and three genes in the cell membrane, other genes were located in multiple cellular structures by subcellular location. Ectopic expression demonstrated that more branching occurred in MdCo38-OE by upregulating NtPIN1 and NtGA2ox and enlarged leaves in MdCo41-OE tobaccos by upregulating NtCCDs. Transcripts of MdCo38 and MdCo41 were associated with the Co genotypes in apples. The results indicate that MdCo38 and MdCo41 are involved in the columnar growth phenotype in apple, probably through altering polar auxin transport, active gibberellin levels, and strigolactone biosynthesis.

Transcriptomic Complexity of Culm Growth and Development in Different Types of Moso Bamboo

Moso bamboo is capable of both sexual and asexual reproduction during natural growth, resulting in four distinct types of culms: the bamboo shoot-culm, the seedling stem, the leptomorph rhizome, and a long-ignored culm-the outward-rhizome. Sometimes, when the outward rhizomes break through the soil, they continue to grow longitudinally and develop into a new individual. However, the roles of alternative transcription start sites (aTSS) or termination sites (aTTS) as well as alternative splicing (AS) have not been comprehensively studied for their development. To re-annotate the moso bamboo genome and identify genome-wide aTSS, aTTS, and AS in growing culms, we utilized single-molecule long-read sequencing technology. In total, 169,433 non-redundant isoforms and 14,840 new gene loci were identified. Among 1311 lncRNAs, most of which showed a positive correlation with their target mRNAs, one-third of these IncRNAs were preferentially expressed in winter bamboo shoots. In addition, the predominant AS type observed in moso bamboo was intron retention, while aTSS and aTTS events occurred more frequently than AS. Notably, most genes with AS events were also accompanied by aTSS and aTTS events. Outward rhizome growth in moso bamboo was associated with a significant increase in intron retention, possibly due to changes in the growth environment. As different types of moso bamboo culms grow and develop, a significant number of isoforms undergo changes in their conserved domains due to the regulation of aTSS, aTTS, and AS. As a result, these isoforms may play different roles than their original functions. These isoforms then performed different functions from their original roles, contributing to the transcriptomic complexity of moso bamboo. Overall, this study provided a comprehensive overview of the transcriptomic changes underlying different types of moso bamboo culm growth and development.

Integration of mRNA and miRNA Analysis Reveals the Post-Transcriptional Regulation of Salt Stress Response in Hemerocallis fulva

MicroRNAs (miRNAs) belong to non-coding small RNAs which have been shown to take a regulatory function at the posttranscriptional level in plant growth development and response to abiotic stress. Hemerocallis fulva is an herbaceous perennial plant with fleshy roots, wide distribution, and strong adaptability. However, salt stress is one of the most serious abiotic stresses to limit the growth and production of Hemerocallis fulva. To identify the miRNAs and their targets involved in the salt stress resistance, the salt-tolerant H. fulva with and without NaCl treatment were used as materials, and the expression differences of miRNAs-mRNAs related to salt-tolerance were explored and the cleavage sites between miRNAs and targets were also identified by using degradome sequencing technology. In this study, twenty and three significantly differential expression miRNAs (p-value < 0.05) were identified in the roots and leaves of H. fulva separately. Additionally, 12,691 and 1538 differentially expressed genes (DEGs) were also obtained, respectively, in roots and leaves. Moreover, 222 target genes of 61 family miRNAs were validated by degradome sequencing. Among the DE miRNAs, 29 pairs of miRNA targets displayed negatively correlated expression profiles. The qRT-PCR results also showed that the trends of miRNA and DEG expression were consistent with those of RNA-seq. A gene ontology (GO) enrichment analysis of these targets revealed that the calcium ion pathway, oxidative defense response, microtubule cytoskeleton organization, and DNA binding transcription factor responded to NaCl stress. Five miRNAs, miR156, miR160, miR393, miR166, and miR396, and several hub genes, squamosa promoter-binding-like protein (SPL), auxin response factor 12 (ARF), transport inhibitor response 1-like protein (TIR1), calmodulin-like proteins (CML), and growth-regulating factor 4 (GRF4), might play central roles in the regulation of NaCl-responsive genes. These results indicate that non-coding small RNAs and their target genes that are related to phytohormone signaling, Ca²⁺ signaling, and oxidative defense signaling pathways are involved in H. fulva's response to NaCl stress.

OsUGE1 is directly targeted by OsGRF6 to regulate root hair length in rice

KEY MESSAGE

Root hairs are required for water and nutrient acquisition in plants. Here, we report a novel mechanism that OsUGE1 is negatively controlled by OsGRF6 to regulate root hair elongation in rice. Root hairs are tubular outgrowths generated by the root epidermal cells. They effectively enlarge the soil-root contact area and play essential roles for nutrient and water absorption. Here, in this study, we demonstrated that the Oryza sativa UDP-glucose 4-epimerase 1-like (OsUGE1) negatively regulated root hair elongation and was directly targeted by Oryza sativa growth regulating factor 6 (OsGRF6). Knockout mutants of OsUGE1 using CRISPR-Cas9 technology showed longer root hairs than those of wild type. In contrast, overexpression lines of OsUGE1 displayed shorter root hair compared with those of wild type. GUS staining showed that it could specifically express in root hair. Subcellular localization analysis indicates that OsUGE1 is located in endoplasmic reticulum, nucleus and plasma membrane. More importantly, ChIP-qPCR, Yeast-one-hybrid and BiFC experiments revealed that OsGRF6 could bind to the promoter of OsUGE1. Furthermore, knockout mutants of OsGRF6 showed shorter root hair than those of wild type, and OsGRF6 dominantly expressed in root. In addition, the expression level of OsUGE1 is significantly downregulated in Osgrf6 mutant. Taken together, our study reveals a novel pathway that OsUGE1 is negatively controlled by OsGRF6 to regulate root hair elongation in rice.

Dissection of Developmental Programs and Regulatory Modules Directing Endosperm Transfer Cell and Aleurone Identity in the Syncytial Endosperm of Barley

Endosperm development in barley starts with the formation of a multinucleate syncytium, followed by cellularization in the ventral part of the syncytium generating endosperm transfer cells (ETCs) as first differentiating subdomain, whereas aleurone (AL) cells will originate from the periphery of the enclosing syncytium. Positional signaling in the syncytial stage determines cell identity in the cereal endosperm. Here, we performed a morphological analysis and employed laser capture microdissection (LCM)-based RNA-seq of the ETC region and the peripheral syncytium at the onset of cellularization to dissect developmental and regulatory programs directing cell specification in the early endosperm. Transcriptome data revealed domain-specific characteristics and identified two-component signaling (TCS) and hormone activities (auxin, ABA, ethylene) with associated transcription factors (TFs) as the main regulatory links for ETC specification. On the contrary, differential hormone signaling (canonical auxin, gibberellins, cytokinin) and interacting TFs control the duration of the syncytial phase and timing of cellularization of AL initials. Domain-specific expression of candidate genes was validated by in situ hybridization and putative protein-protein interactions were confirmed by split-YFP assays. This is the first transcriptome analysis dissecting syncytial subdomains of cereal seeds and provides an essential framework for initial endosperm differentiation in barley, which is likely also valuable for comparative studies with other cereal crops.

Comprehensive evolutionary analysis of growth-regulating factor gene family revealing the potential molecular basis under multiple hormonal stress in Gramineae crops

Growth-regulating factors (GRFs) are plant-specific transcription factors that contain two highly conserved QLQ and WRC domains, which control a range of biological functions, including leaf growth, floral organ development, and phytohormone signaling. However, knowledge of the evolutionary patterns and driving forces of GRFs in Gramineae crops is limited and poorly characterized. In this study, a total of 96 GRFs were identified from eight crops of Brachypodium distachyon, Hordeum vulgare, Oryza sativa L. ssp. indica, Oryza rufipogon, Oryza sativa L. ssp. japonica, Setaria italic, Sorghum bicolor and Zea mays. Based on their protein sequences, the GRFs were classified into three groups. Evolutionary analysis indicated that the whole-genome or segmental duplication plays an essential role in the GRFs expansion, and the GRFs were negatively selected during the evolution of Gramineae crops. The GRFs protein function as transcriptional activators with distinctive structural motifs in different groups. In addition, the expression of GRFs was induced under multiple hormonal stress, including IAA, BR, GA3, 6BA, ABA, and MeJ treatments. Specifically, OjGRF11 was significantly induced by IAA at 6 h after phytohormone treatment. Transgenic experiments showed that roots overexpressing OjGRF11 were more sensitive to IAA and affect root elongation. This study will broaden our insights into the origin and evolution of the GRF family in Gramineae crops and will facilitate further research on GRF function.

Interplay among ZF-HD and GRF transcription factors during Arabidopsis leaf development

The growth-regulating factor (GRF) family of transcriptional factors are involved in the control of leaf size and senescence, inflorescence and root growth, grain size, and plant regeneration. However, there is limited information about the genes regulated by these transcriptional factors, which are in turn responsible for their functions. Using a meta-analysis approach, we identified genes encoding Arabidopsis (Arabidopsis thaliana) zinc-finger homeodomain (ZF-HD) transcriptional factors, as potential targets of the GRFs. We further showed that GRF3 binds to the promoter of one of the members of the ZF-HD family, HOMEOBOX PROTEIN 33 (HB33), and activates its transcription. Increased levels of HB33 led to different modifications in leaf cell number and size that were dependent on its expression levels. Furthermore, we found that expression of HB33 for an extended period during leaf development increased leaf longevity. To cope with the functional redundancy among ZF-HD family members, we generated a dominant repressor version of HB33, HB33-SRDX. Expression of HB33-SRDX from HB33 regulatory regions was seedling-lethal, revealing the importance of the ZF-HD family in plant development. Misexpression of HB33-SRDX in early leaf development caused a reduction in both cell size and number. Interestingly, the loss-of-function of HB33 in lines carrying a GRF3 allele insensitive to miR396 reverted the delay in leaf senescence characteristic of these plants. Our results revealed functions for ZF-HDs in leaf development and linked them to the GRF pathway.