Synchronization of developmental processes and defense signaling by growth regulating transcription factors

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.

Transcriptomic and metabolomic profiling reveals the drought tolerance mechanism of Illicium difengpi (Schisandraceae)

Illicium difengpi (Schisandraceae), an endangered medicinal plant endemic to karst areas, is highly tolerant to drought and thus can be used as an ideal material for investigating adaptive mechanism to drought stress. The understanding of the drought tolerance of I. difengpi, especially at the molecular level, is lacking. In the present study, we aimed to clarify the molecular mechanism underlying drought tolerance in endemic I. difengpi plant in karst regions. The response characteristics of transcripts and changes in metabolite abundance of I. difengpi subjected to drought and rehydration were analyzed, the genes and key metabolites responsive to drought and rehydration were screened, and some important biosynthetic and secondary metabolic pathways were identified. A total of 231,784 genes and 632 metabolites were obtained from transcriptome and metabolome analyses, and most of the physiological metabolism in drought-treated I. difengpi plants recovered after rehydration. There were more upregulated genes than downregulated genes under drought and rehydration treatments, and rehydration treatment induced stable expression of 65.25% of genes, indicating that rehydration alleviated drought stress to some extent. Drought and rehydration treatment generated flavonoids, phenolic acids, flavonols, amino acids and their derivatives, as well as metabolites such as saccharides and alcohols in the leaves of I. difengpi plants, which alleviated the injury caused by excessive reactive oxygen species. The integration of transcriptome and metabolome analyses showed that, under drought stress, I. difengpi increased glutathione, flavonoids, polyamines, soluble sugars and amino acids, contributing to cell osmotic potential and antioxidant activity. The results show that the high drought tolerance and recovery after rehydration are the reasons for the normal growth of I. difengpi in karst mountain areas.

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.

Growth-regulating factor 15-mediated gene regulatory network enhances salt tolerance in poplar

Soil salinity is an important determinant of crop productivity and triggers salt stress response pathways in plants. The salt stress response is controlled by transcriptional regulatory networks that maintain regulatory homeostasis through combinations of transcription factor (TF)-DNA and TF-TF interactions. We investigated the transcriptome of poplar 84 K (Populus alba × Populus glandulosa) under salt stress using samples collected at 4- or 6-h intervals within 2 days of salt stress treatment. We detected 24,973 differentially expressed genes, including 2,231 TFs that might be responsive to salt stress. To explore these interactions and targets of TFs in perennial woody plants, we combined gene regulatory networks, DNA affinity purification sequencing, yeast two-hybrid-sequencing, and multi-gene association approaches. Growth-regulating factor 15 (PagGRF15) and its target, high-affinity K+ transporter 6 (PagHAK6), were identified as an important regulatory module in the salt stress response. Overexpression of PagGRF15 and PagHAK6 in transgenic lines improved salt tolerance by enhancing Na+ transport and modulating H2O2 accumulation in poplar. Yeast two-hybrid assays identified more than 420 PagGRF15-interacting proteins, including ETHYLENE RESPONSE FACTOR TFs and a zinc finger protein (C2H2) that are produced in response to a variety of phytohormones and environmental signals and are likely involved in abiotic stress. Therefore, our findings demonstrate that PagGRF15 is a multifunctional TF involved in growth, development, and salt stress tolerance, highlighting the capability of a multifaceted approach in identifying regulatory nodes in plants.

Characterization of microRNAs and Target Genes in Musa acuminata subsp. burmannicoides, var. Calcutta 4 during Interaction with Pseudocercospora musae

Endogenous microRNAs (miRNAs) are small non-coding RNAs that perform post-transcriptional regulatory roles across diverse cellular processes, including defence responses to biotic stresses. Pseudocercospora musae, the causal agent of Sigatoka leaf spot disease in banana (Musa spp.), is an important fungal pathogen of the plant. Illumina HiSeq 2500 sequencing of small RNA libraries derived from leaf material in Musa acuminata subsp. burmannicoides, var. Calcutta 4 (resistant) after inoculation with fungal conidiospores and equivalent non-inoculated controls revealed 202 conserved miRNAs from 30 miR-families together with 24 predicted novel miRNAs. Conserved members included those from families miRNA156, miRNA166, miRNA171, miRNA396, miRNA167, miRNA172, miRNA160, miRNA164, miRNA168, miRNA159, miRNA169, miRNA393, miRNA535, miRNA482, miRNA2118, and miRNA397, all known to be involved in plant immune responses. Gene ontology (GO) analysis of gene targets indicated molecular activity terms related to defence responses that included nucleotide binding, oxidoreductase activity, and protein kinase activity. Biological process terms associated with defence included response to hormone and response to oxidative stress. DNA binding and transcription factor activity also indicated the involvement of miRNA target genes in the regulation of gene expression during defence responses. sRNA-seq expression data for miRNAs and RNAseq data for target genes were validated using stem-loop quantitative real-time PCR (qRT-PCR). For the 11 conserved miRNAs selected based on family abundance and known involvement in plant defence responses, the data revealed a frequent negative correlation of expression between miRNAs and target host genes. This examination provides novel information on miRNA-mediated host defence responses, applicable in genetic engineering for the control of Sigatoka leaf spot disease.

Growth-regulating factors: conserved and divergent roles in plant growth and development and potential value for crop improvement

High yield and stress resistance are the major prerequisites for successful crop cultivation, and can be achieved by modifying plant architecture. Evolutionarily conserved growth-regulating factors (GRFs) control the growth of different tissues and organs of plants. Here, we provide a systematic overview of the expression patterns of GRF genes and the structural features of GRF proteins in different plant species. Moreover, we illustrate the conserved and divergent roles of GRFs, microRNA396 (miR396), and GRF-interacting factors (GIFs) in leaf, root, and flower development. We also describe the molecular networks involving the miR396-GRF-GIF module, and illustrate how this module coordinates with different signaling molecules and transcriptional regulators to control development of different plant species. GRFs promote leaf growth, accelerate grain filling, and increase grain size and weight. We also provide some molecular insight into how coordination between GRFs and other signaling modules enhances crop productivity; for instance, how the GRF-DELLA interaction confers yield-enhancing dwarfism while increasing grain yield. Finally, we discuss how the GRF-GIF chimera substantially improves plant transformation efficiency by accelerating shoot formation. Overall, we systematically review the conserved and divergent roles of GRFs and the miR396-GRF-GIF module in growth regulation, and also provide insights into how GRFs can be utilized to improve the productivity and nutrient content of crop plants.

Genome-wide identification, phylogeny, and expression analysis of GRF transcription factors in pineapple (Ananas comosus)

Background

Pineapple is the only commercially grown fruit crop in the Bromeliaceae family and has significant agricultural, industrial, economic, and ornamental value. GRF (growth-regulating factor) proteins are important transcription factors that have evolved in seed plants (embryophytes). They contain two conserved domains, QLQ (Gln, Leu, Gln) and WRC (Trp, Arg, Cys), and regulate multiple aspects of plant growth and stress response, including floral organ development, leaf growth, and hormone responses. The GRF family has been characterized in a number of plant species, but little is known about this family in pineapple and other bromeliads.

Main discoveries

We identified eight GRF transcription factor genes in pineapple, and phylogenetic analysis placed them into five subfamilies (I, III, IV, V, VI). Segmental duplication appeared to be the major contributor to expansion of the AcGRF family, and the family has undergone strong purifying selection during evolution. Relative to that of other gene families, the gene structure of the GRF family showed less conservation. Analysis of promoter cis-elements suggested that AcGRF genes are widely involved in plant growth and development. Transcriptome data and qRT-PCR results showed that, with the exception of AcGRF5, the AcGRFs were preferentially expressed in the early stage of floral organ development and AcGRF2 was strongly expressed in ovules. Gibberellin treatment significantly induced AcGRF7/8 expression, suggesting that these two genes may be involved in the molecular regulatory pathway by which gibberellin promotes pineapple fruit expansion.

Conclusion

AcGRF proteins appear to play a role in the regulation of floral organ development and the response to gibberellin. The information reported here provides a foundation for further study of the functions of AcGRF genes and the traits they regulate.

Genome-wide characterization and expression analysis of the growth-regulating factor family in Saccharum

BACKGROUND

Growth regulating factors (GRFs) are transcription factors that regulate diverse biological and physiological processes in plants, including growth, development, and abiotic stress. Although GRF family genes have been studied in a variety of plant species, knowledge about the identification and expression patterns of GRFs in sugarcane (Saccharum spp.) is still lacking.

RESULTS

In the present study, a comprehensive analysis was conducted in the genome of wild sugarcane (Saccharum spontaneum) and 10 SsGRF genes were identified and characterized. The phylogenetic relationship, gene structure, and expression profiling of these genes were analyzed entirely under both regular growth and low-nitrogen stress conditions. Phylogenetic analysis suggested that the 10 SsGRF members were categorized into six clusters. Gene structure analysis indicated that the SsGRF members in the same group were greatly conserved. Expression profiling demonstrated that most SsGRF genes were extremely expressed in immature tissues, implying their critical roles in sugarcane growth and development. Expression analysis based on transcriptome data and real-time quantitative PCR verification revealed that GRF1 and GRF3 were distinctly differentially expressed in response to low-nitrogen stress, which meant that they were additional participated in sugarcane stress tolerance.

CONCLUSION

Our study provides a scientific basis for the potential functional prediction of SsGRF and will be further scrutinized by examining their regulatory network in sugarcane development and abiotic stress response, and ultimately facilitating their application in cultivated sugarcane breeding.

Construction of two regulatory networks related to salt stress and lignocellulosic synthesis under salt stress based on a Populus davidiana × P. bolleana transcriptome analysis

Construction of ML-hGRN for the salt pathway in Populus davidiana × P. bolleana. Construction of ML-hGRN for the lignocellulosic pathway in Populus davidiana × P. bolleana under salt stress. Many woody plants, including Populus davidiana × P. bolleana, have made great contributions to human production and life. High salt is one of the main environmental factors that restricts the growth of poplar. This study found that high salt could induce strong biochemical changes in poplar. To detect the effect of salt treatment on gene expression, 18 libraries were sequenced on the Illumina sequencing platform. The results identified a large number of early differentially expressed genes (DEGs) and a small number of late DEGs, which indicated that most of the salt response genes of poplar were early response genes. In addition, 197 TFs, including NAC, ERF, and other TFs related to salt stress, were differentially expressed during salt treatment, which indicated that these TFs may play an important role in the salt stress response of poplar. Based on the RNA-seq analysis results, multilayered hierarchical gene regulatory networks (ML-hGRNs) of salt stress- and lignocellulosic synthesis-related DEGs were constructed using the GGM algorithm. The lignocellulosic synthesis regulatory network under salt stress revealed that lignocellulosic synthesis might play an important role in the process of salt stress resistance. Furthermore, the NAC family transcription factor PdbNAC83, which was found in the upper layer in both pathways, was selected to verify the accuracy of the ML-hGRNs. DAP-seq showed that the binding site of PdbNAC83 included a "TT(G/A)C(G/T)T" motif, and ChIP-PCR further verified that PdbNAC83 can regulate the promoters of at least six predicted downstream genes (PdbNLP2-2, PdbZFP6, PdbMYB73, PdbC2H2-like, PdbMYB93-1, PdbbHLH094) by binding to the "TT(G/A)C(G/T)T" motif, which indicates that the predicted regulatory network diagram obtained in this study is relatively accurate. In conclusion, a species-specific salt response pathway might exist in poplar, and this finding lays a foundation for further study of the regulatory mechanism of the salt stress response and provides new clues for the use of genetic engineering methods to create high-quality and highly resistant forest germplasms.

Genome-Wide Identification and Characterization of Growth Regulatory Factor Family Genes in Medicago

Growth Regulatory Factors (GRF) are plant-specific transcription factors that play critical roles in plant growth and development as well as plant tolerance against stress. In this study, a total of 16 GRF genes were identified from the genomes of Medicago truncatula and Medicago sativa. Multiple sequence alignment analysis showed that all these members contain conserved QLQ and WRC domains. Phylogenetic analysis suggested that these GRF proteins could be classified into five clusters. The GRF genes showed similar exon–intron organizations and similar architectures in their conserved motifs. Many stress-related cis-acting elements were found in their promoter region, and most of them were related to drought and defense response. In addition, analyses on microarray and transcriptome data indicated that these GRF genes exhibited distinct expression patterns in various tissues or in response to drought and salt treatments. In particular, qPCR results showed that the expression levels of gene pairs MtGRF2–MsGRF2 and MtGRF6–MsGRF6 were significantly increased under NaCl and mannitol treatments, indicating that they are most likely involved in salt and drought stress tolerance. Collectively, our study is valuable for further investigation on the function of GRF genes in Medicago and for the exploration of GRF genes in the molecular breeding of highly resistant M. sativa.

Functional conservation and divergence in plant-specific GRF gene family revealed by sequences and expression analysis

Unique to plants, growth regulatory factors (GRFs) play important roles in plant growth and reproduction. This study investigated the evolutionary and functional characteristics associated with plant growth. Using genome-wide analysis of 15 plant species, 173 members of the GRF family were identified and phylogenetically categorized into six groups. All members contained WRC and QLQ conserved domains, and the family’s expansion largely depended on segmental duplication. The promoter region of the GRF gene family mainly contained four types of cis-acting elements (light-responsive elements, development-related elements, hormone-responsive elements, and environmental stress-related elements) that are mainly related to gene expression levels. Functional divergence analysis revealed that changes in amino acid site evolution rate played a major role in the differentiation of the GRF gene family, with ten significant sites identified. Six significant sites were identified for positive selection. Moreover, the four groups of coevolutionary sites identified may play a key role in regulating the transcriptional activation of the GRF protein. Expression profiles revealed that GRF genes were generally highly expressed in young plant tissues and had tissue or organ expression specificity, demonstrating their functional conservation with distinct divergence. The results of these sequence and expression analyses are expected to provide molecular evolutionary and functional references for the plant GRF gene family.

Jasmonic Acid-Dependent MYC Transcription Factors Bind to a Tandem G-Box Motif in the YUCCA8 and YUCCA9 Promoters to Regulate Biotic Stress Responses

The indole-3-pyruvic acid pathway is the main route for auxin biosynthesis in higher plants. Tryptophan aminotransferases (TAA1/TAR) and members of the YUCCA family of flavin-containing monooxygenases catalyze the conversion of l-tryptophan via indole-3-pyruvic acid to indole-3-acetic acid (IAA). It has been described that jasmonic acid (JA) locally produced in response to mechanical wounding triggers the de novo formation of IAA through the induction of two YUCCA genes, YUC8 and YUC9. Here, we report the direct involvement of a small number of basic helix-loop-helix transcription factors of the MYC family in this process. We show that the JA-mediated regulation of the expression of the YUC8 and YUC9 genes depends on the abundance of MYC2, MYC3, and MYC4. In support of this observation, seedlings of myc knockout mutants displayed a strongly reduced response to JA-mediated IAA formation. Furthermore, transactivation assays provided experimental evidence for the binding of MYC transcription factors to a particular tandem G-box motif abundant in the promoter regions of YUC8 and YUC9, but not in the promoters of the other YUCCA isogenes. Moreover, we demonstrate that plants that constitutively overexpress YUC8 and YUC9 show less damage after spider mite infestation, thereby underlining the role of auxin in plant responses to biotic stress signals.

MicroRNAs play important roles in regulating the rapid growth of the Phyllostachys edulis culm internode

Moso bamboo (Phyllostachys edulis) is a fast-growing species with uneven growth and lignification from lower to upper segments within one internode. MicroRNAs (miRNAs) play a vital role in post-transcriptional regulation in plants. However, how miRNAs regulate fast growth in bamboo internodes is poorly understood. In this study, one moso bamboo internode was divided during early rapid growth into four segments called F4 (bottom) to F1 (upper) and these were then analysed for transcriptomes, miRNAs and degradomes. The F4 segment had a higher number of actively dividing cells as well as a higher content of auxin (IAA), cytokinin (CK) and gibberellin (GA) compared with the F1 segment. RNA-seq analysis showed DNA replication and cell division-associated genes highly expressed in F4 rather than in F1. In total, 63 miRNAs (DEMs) were identified as differentially expressed between F4 and F1. The degradome and the transcriptome indicated that many downstream transcription factors and hormonal responses genes were modulated by DEMs. Several miR-target interactions were further validated by tobacco co-infiltration. Our findings give new insights into miRNA-mediated regulatory pathways in bamboo, and will contribute to a comprehensive understanding of the molecular mechanisms governing rapid growth.

Genetic toolbox and regulatory circuits of plant-nematode associations

Plant-nematode associations are the most imperative area of study that forms the basis to understand their regulatory networks and coordinated functional aspects. Nematodes are highly parasitic organisms known so far, to cause relentless damage towards agricultural crops on a global scale. They pierce the roots of host plants and form neo-plastic feeding structures to extract out resources for their functional development. Moreover, they undergo re-differentiation within plant cells to form giant multi-nucleate feeding structures or syncytium. All these processes are facilitated by numerous transcriptomic, proteomic, metabolomic and epigenetic modifications, that regulate different biological attractions among plants and nematodes. Nevertheless, these mechanisms are quite remarkable and have been explored in the present review. Here, we have shed light on genomic as well as genetic approaches to acquire an effective understanding regarding plant-nematode associations. Transcriptomics have revealed an extensive network to unravel feeding mechanism of nematodes through gene-expression programming of target genes. Also, the regulatory circuits of epigenetic alterations through DNA-methylation, non-coding RNAs and histone modifications very well explain epigenetic profiling within plants. Since decades, research have observed many intricacies to elucidate the dynamic nature of epigenetic modulations in plant-nematode attractions. By this review, we have highlighted the functional aspects of small RNAs in inducing plant-nematode parasitism along with the putative role of miRNAs. These RNAs act as chief genetic elements to mediate the expressional changes in plants through post-transcriptional silencing of various effector proteins as well as transcriptional factors. A pragmatic role of miRNAs in modulating gene expression in nematode infection and feeding site development have also been reviewed. Hence, they have been considered master regulators for functional reprogramming the expression during establishment of feeding sites. We have also encapsulated the advancement of genome-broadened DNA-methylation and untangled the nematode mediated dynamic alterations within plant methylome along with assessing transcriptional activities of various genes and transposons. In particular, we have highlighted the role of effector proteins in stimulating epigenetic changes. Finally, we have emerged towards a molecular-based core understanding about plant-nematode associations.

Time-course transcriptomic analysis of Petunia ×hybrida leaves under water deficit stress using RNA sequencing

Water deficit limits plant growth and development, resulting in quality loss of horticultural crops. However, there is limited information on gene regulation and signaling pathways related to water deficit stress response at multiple time points. The objective of this research was to investigate global gene expression patterns under water deficit stress to provide an insight into how petunia (Petunia ×hybrida 'Mitchell Diploid') responded in the process of stress. Nine-week-old petunias were irrigated daily or placed under water stress by withholding water. Stressed plants reduced stomatal conductance after five days of water deficit, indicating they perceived stress and initiated stress response mechanisms. To analyze transcriptomic changes at the early stage of water deficit, leaf tissue samples were collected 1, 3, and 5 days after water was withheld for RNA sequencing. Under water deficit stress, 154, 3611, and 980 genes were upregulated and 41, 2806, and 253 genes were downregulated on day 1, 3, and 5, respectively. Gene Ontology analysis revealed that redox homeostasis processes through sulfur and glutathione metabolism pathways, and hormone signal transduction, especially abscisic acid and ethylene, were enriched under water deficit stress. Thirty-four transcription factor families were identified, including members of AP2/ERF, NAC, MYB-related, C2H2, and bZIP families, and TFs in AP2/ERF family was the most abundant in petunia. Interestingly, only one member of GRFs was upregulated on day 1, while most of TFs were differentially expressed on day 3 and/or 5. The transcriptome data from this research will provide valuable molecular resources for understanding the early stages of water stress-responsive networks as well as engineering petunia with enhanced water stress tolerance.

Increased water use efficiency in miR396-downregulated tomato plants

MicroRNAs regulate plant development and responses to biotic and abiotic stresses but their impact on water use efficiency (WUE) is poorly known. Increasing WUE is a major task in crop improvement programs aimed to meet the challenges posed by the reduction in water availability associated with the ongoing climatic change. We have examined the physiological and molecular response to water stress of tomato (Solanum lycopersicum L.) plants downregulated for miR396 by target mimicry. In water stress conditions, miR396-downregulated plants displayed reduced transpiration and a less then proportional decrease in the photosynthetic rate that resulted in higher WUE. The increase in WUE was associated with faster foliar accumulation of abscisic acid (ABA), with the induction of several drought-protective genes and with the activation of the jasmonic acid (JA) and γ-aminobutyric acid (GABA) pathways. We propose a model in which the downregulation of miR396 leads to the activation of a complex molecular response to water stress. This response acts synergistically with a set of leaf morphological modifications to increase stomatal closure and preserve the efficiency of the photosynthetic activity, ultimately resulting in higher WUE.

Genome-wide analysis of growth-regulating factors (GRFs) in Triticum aestivum

The Growth-Regulating Factor (GRF) family encodes a type of plant-specific transcription factor (TF). GRF members play vital roles in plant development and stress response. Although GRF family genes have been investigated in a variety of plants, they remain largely unstudied in bread wheat (Triticum aestivum L.). The present study was conducted to comprehensively identify and characterize the T. aestivum GRF (TaGRF) gene family members. We identified 30 TaGRF genes, which were divided into four groups based on phylogenetic relationship. TaGRF members within the same subgroup shared similar motif composition and gene structure. Synteny analysis suggested that duplication was the dominant reason for family member expansion. Expression pattern profiling showed that most TaGRF genes were highly expressed in growing tissues, including shoot tip meristems, stigmas and ovaries, suggesting their key roles in wheat growth and development. Further qRT-PCR analysis revealed that all 14 tested TaGRFs were significantly differentially expressed in responding to drought or salt stresses, implying their additional involvement in stress tolerance of wheat. Our research lays a foundation for functional determination of TaGRFs, and will help to promote further scrutiny of their regulatory network in wheat development and stress response.

Genome-Wide Identification of Populus Malectin/Malectin-Like Domain-Containing Proteins and Expression Analyses Reveal Novel Candidates for Signaling and Regulation of Wood Development

Malectin domain (MD) is a ligand-binding protein motif of pro- and eukaryotes. It is particularly abundant in Viridiplantae, where it occurs as either a single (MD, PF11721) or tandemly duplicated domain (PF12819) called malectin-like domain (MLD). In herbaceous plants, MD- or MLD-containing proteins (MD proteins) are known to regulate development, reproduction, and resistance to various stresses. However, their functions in woody plants have not yet been studied. To unravel their potential role in wood development, we carried out genome-wide identification of MD proteins in the model tree species black cottonwood (Populus trichocarpa), and analyzed their expression and co-expression networks. P. trichocarpa had 146 MD genes assigned to 14 different clades, two of which were specific to the genus Populus. 87% of these genes were located on chromosomes, the rest being associated with scaffolds. Based on their protein domain organization, and in agreement with the exon-intron structures, the MD genes identified here could be classified into five superclades having the following domains: leucine-rich repeat (LRR)-MD-protein kinase (PK), MLD-LRR-PK, MLD-PK (CrRLK1L), MLD-LRR, and MD-Kinesin. Whereas the majority of MD genes were highly expressed in leaves, particularly under stress conditions, eighteen showed a peak of expression during secondary wall formation in the xylem and their co-expression networks suggested signaling functions in cell wall integrity, pathogen-associated molecular patterns, calcium, ROS, and hormone pathways. Thus, P. trichocarpa MD genes having different domain organizations comprise many genes with putative foliar defense functions, some of which could be specific to Populus and related species, as well as genes with potential involvement in signaling pathways in other tissues including developing wood.

Epigenetic Mechanisms in Nematode-Plant Interactions

Epigenetic mechanisms play fundamental roles in regulating numerous biological processes in various developmental and environmental contexts. Three highly interconnected epigenetic control mechanisms, including small noncoding RNAs, DNA methylation, and histone modifications, contribute to the establishment of plant epigenetic profiles. During the past decade, a growing body of experimental work has revealed the intricate, diverse, and dynamic roles that epigenetic modifications play in plant-nematode interactions. In this review, I summarize recent progress regarding the functions of small RNAs in mediating plant responses to infection by cyst and root-knot nematodes, with a focus on the functions of microRNAs. I also recapitulate recent advances in genome-wide DNA methylation analysis and discuss how cyst nematodes induce extensive and dynamic changes in the plant methylome that impact the transcriptional activity of genes and transposable elements. Finally, the potential role of nematode effector proteins in triggering such epigenome changes is discussed.

Mutations in the miR396 binding site of the growth-regulating factor gene VvGRF4 modulate inflorescence architecture in grapevine

Bunch rot caused by Botrytis cinerea infections is a notorious problem in grapevine cultivation. To produce high quality fruits, grapevine plants are treated with fungicides, which is cost intensive and harmful to the environment. Conversely, loose cluster bunches show a considerably enhanced physical resilience to bunch diseases. With the aim to identify genetic determinants that modulate the development of bunch architecture, we have compared loose and compact 'Pinot noir' clones. Loose cluster architecture was found to be correlated with increased berry size, elongated rachis and elongated pedicels. Using transcriptome analysis in combination with whole genome sequencing, we have identified a growth-regulating factor gene, VvGRF4, upregulated and harbours heterozygous mutations in the loose cluster clones. At late stages of inflorescence development, the mRNA pools of loose cluster clones contain predominantly mRNAs derived from the mutated alleles, which are resistant to miR396 degradation. Expression of the VvGRF4 gene and its mutated variants in Arabidopsis demonstrates that it promotes pedicel elongation. Taken together, VvGRF4 modulates bunch architecture in grapevine 'Pinot noir' clones. This trait can be introduced into other cultivars using marker-assisted breeding or CRISPR-Cas9 technology. Related growth-regulating factors or other genes of the same pathway may have similar functions.

A role for Arabidopsis growth-regulating factors 1 and 3 in growth-stress antagonism

Growth-regulating factors (GRFs) belong to a small family of transcription factors that are highly conserved in plants. GRFs regulate many developmental processes and plant responses to biotic and abiotic stimuli. Despite the importance of GRFs, a detailed mechanistic understanding of their regulatory functions is still lacking. In this study, we used ChIP sequencing (ChIP-seq) to identify genome-wide binding sites of Arabidopsis GRF1 and GRF3, and correspondingly their direct downstream target genes. RNA-sequencing (RNA-seq) analysis revealed that GRF1 and GRF3 regulate the expression of a significant number of the identified direct targets. The target genes unveiled broad regulatory functions of GRF1 and GRF3 in plant growth and development, phytohormone biosynthesis and signaling, and the cell cycle. Our analyses also revealed that clock core genes and genes with stress- and defense-related functions are most predominant among the GRF1- and GRF3-bound targets, providing insights into a possible role for these transcription factors in mediating growth-defense antagonism and integrating environmental stimuli into developmental programs. Additionally, GRF1 and GRF3 target molecular nodes of growth-defense antagonism and modulate the levels of defense- and development-related hormones in opposite directions. Taken together, our results point to GRF1 and GRF3 as potential key determinants of plant fitness under stress conditions.

The miRNA-mRNA Networks Involving Abnormal Energy and Hormone Metabolisms Restrict Tillering in a Wheat Mutant dmc

Tillers not only determine plant architecture but also influence crop yield. To explore the miRNA regulatory network restraining tiller development in a dwarf-monoculm wheat mutant (dmc) derived from Guomai 301 (wild type, WT), we employed miRNome and transcriptome integrative analysis, real-time qRT-PCR, histochemistry, and determinations of the key metabolites and photosynthesis parameters. A total of 91 differentially expressed miRNAs (DEMs) were identified between dmc and WT. Among them, 40 key DEMs targeted 45 differentially expressed genes (DEGs) including the key DEGs encode growth-regulating factors (GRF), auxin response factors (ARF), and other proteins involved in the metabolisms of hormones and carbohydrates, etc. Compared with WT, both the chlorophyll contents and the photosynthesis rate were lower in dmc. The contents of glucose, sucrose, fructose, and maltose were lower in dmc. The contents of auxin (IAA) and zeatin (ZA) were significantly lower, but gibberellin (GA) was significantly higher in the tiller tissues of dmc. This research demonstrated that the DEMs regulating hormone and carbohydrate metabolisms were important causes for dmc to not tiller. A primary miRNA-mRNA regulatory model for dmc tillering was established. The lower photosynthesis rate, insufficient energy, and abnormal hormone metabolisms restrict tillering in dmc.

Biological roles and an evolutionary sketch of the GRF-GIF transcriptional complex in plants

GROWTH-REGULATING FACTORs (GRFs) are sequence-specific DNA-binding transcription factors that regulate various aspects of plant growth and development. GRF proteins interact with a transcription cofactor, GRF-INTERACTING FACTOR (GIF), to form a functional transcriptional complex. For its activities, the GRF-GIF duo requires the SWITCH2/SUCROSE NONFERMENTING2 chromatin remodeling complex. One of the most conspicuous roles of the duo is conferring the meristematic potential on the proliferative and formative cells during organogenesis. GRF expression is post-transcriptionally down-regulated by microRNA396 (miR396), thus constructing the GRF-GIF-miR396 module and fine-tuning the duo’s action. Since the last comprehensive review articles were published over three years ago, many studies have added further insight into its action and elucidated new biological roles. The current review highlights recent advances in our understanding of how the GRF-GIF-miR396 module regulates plant growth and development. In addition, I revise the previous view on the evolutionary origin of the GRF gene family.

BjuWRR1, a CC-NB-LRR gene identified in Brassica juncea, confers resistance to white rust caused by Albugo candida

BjuWRR1, a CNL-type R gene, was identified from an east European gene pool line of Brassica juncea and validated for conferring resistance to white rust by genetic transformation. White rust caused by the oomycete pathogen Albugo candida is a significant disease of crucifer crops including Brassica juncea (mustard), a major oilseed crop of the Indian subcontinent. Earlier, a resistance-conferring locus named AcB1-A5.1 was mapped in an east European gene pool line of B. juncea-Donskaja-IV. This line was tested along with some other lines of B. juncea (AABB), B. rapa (AA) and B. nigra (BB) for resistance to six isolates of A. candida collected from different mustard growing regions of India. Donskaja-IV was found to be completely resistant to all the tested isolates. Sequencing of a BAC spanning the locus AcB1-A5.1 showed the presence of a single CC-NB-LRR protein encoding R gene. The genomic sequence of the putative R gene with its native promoter and terminator was used for the genetic transformation of a susceptible Indian gene pool line Varuna and was found to confer complete resistance to all the isolates. This is the first white rust resistance-conferring gene described from Brassica species and has been named BjuWRR1. Allelic variants of the gene in B. juncea germplasm and orthologues in the Brassicaceae genomes were studied to understand the evolutionary dynamics of the BjuWRR1 gene.

Conserved miR396b-GRF Regulation Is Involved in Abiotic Stress Responses in Pitaya (Hylocereus polyrhizus)

MicroRNA396 (miR396) is a conserved microRNA family that targets growth-regulating factors (GRFs), which play significant roles in plant growth and stress responses. Available evidence justifies the idea that miR396-targeted GRFs have important functions in many plant species; however, no genome-wide analysis of the pitaya (Hylocereus polyrhizus) miR396 gene has yet been reported. Further, its biological functions remain elusive. To uncover the regulatory roles of miR396 and its targets, the hairpin sequence of pitaya miR396b and the open reading frame (ORF) of its target, HpGRF6, were isolated from pitaya. Phylogenetic analysis showed that the precursor miR396b (MIR396b) gene of plants might be clustered into three major groups, and, generally, a more recent evolutionary relationship in the intra-family has been demonstrated. The sequence analysis indicated that the binding site of hpo-miR396b in HpGRF6 is located at the conserved motif which codes the conserved "RSRKPVE" amino acid in the Trp-Arg-Cys (WRC) region. In addition, degradome sequencing analysis confirmed that four GRFs (GRF1, c56908.graph_c0; GRF4, c52862.graph_c0; GRF6, c39378.graph_c0 and GRF9, c54658.graph_c0) are hpo-miR396b targets that are regulated by specific cleavage at the binding site between the 10th and 11th nucleotides from the 5' terminus of hpo-miR396b. Furthermore, quantitative real-time polymerase chain reaction (qRT-PCR) analysis showed that hpo-miR396b is down-regulated when confronted with drought stress (15% polyethylene glycol, PEG), and its expression fluctuates under other abiotic stresses, i.e., low temperature (4 ± 1 °C), high temperature (42 ± 1 °C), NaCl (100 mM), and abscisic acid (ABA; 0.38 mM). Conversely, the expression of HpGRF6 showed the opposite trend to exposure to these abiotic stresses. Taken together, hpo-miR396b plays a regulatory role in the control of HpGRF6, which might influence the abiotic stress response of pitaya. This is the first documentation of this role in pitaya and improves the understanding of the molecular mechanisms underlying the tolerance to drought stress in this fruit.

Homeostasis in the soybean miRNA396-GRF network is essential for productive soybean cyst nematode infections

Heterodera glycines, the soybean cyst nematode, penetrates soybean roots and migrates to the vascular cylinder where it forms a feeding site called the syncytium. MiRNA396 (miR396) targets growth-regulating factor (GRF) genes, and the miR396-GRF1/3 module is a master regulator of syncytium development in model cyst nematode H. schachtii infection of Arabidopsis. Here, we investigated whether this regulatory system operates similarly in soybean roots and is likewise important for H. glycines infection. We found that a network involving nine MIR396 and 23 GRF genes is important for normal development of soybean roots and that GRF function is specified in the root apical meristem by miR396. All MIR396 genes are down-regulated in the syncytium during its formation phase while, specifically, 11 different GRF genes are up-regulated. The switch to the syncytium maintenance phase coincides with up-regulation of MIR396 and down-regulation of the 11 GRF genes specifically via post-transcriptional regulation by miR396. Furthermore, interference with the miR396-GRF6/8-13/15-17/19 regulatory network, through either overexpression or knockdown experiments, does not affect the number of H. glycines juveniles that enter the vascular cylinder to initiate syncytia, but specifically inhibits efficient H. glycines development to adult females. Therefore, homeostasis in the miR396-GRF6/8-13/15-17/19 regulatory network is essential for productive H. glycines infections.

Integration Analysis of Small RNA and Degradome Sequencing Reveals MicroRNAs Responsive to Dickeya zeae in Resistant Rice

Rice foot rot disease caused by the pathogen Dickeya zeae (formerly known as Erwinia chrysanthemi pv. zeae), is a newly emerging damaging bacterial disease in China and the southeast of Asia, resulting in the loss of yield and grain quality. However, the genetic resistance mechanisms mediated by miRNAs to D. zeae are unclear in rice. In the present study, 652 miRNAs including osa-miR396f predicted to be involved in multiple defense responses to D. zeae were identified with RNA sequencing. A total of 79 differentially expressed miRNAs were detected under the criterion of normalized reads ≥10, including 51 known and 28 novel miRNAs. Degradome sequencing identified 799 targets predicted to be cleaved by 168 identified miRNAs. Among them, 29 differentially expressed miRNA and target pairs including miRNA396f-OsGRFs were identified by co-expression analysis. Overexpression of the osa-miR396f precursor in a susceptible rice variety showed enhanced resistance to D. zeae, coupled with significant accumulation of transcripts of osa-miR396f and reduction of its target the Growth-Regulating Factors (OsGRFs). Taken together, these findings suggest that miRNA and targets including miR396f-OsGRFs have a role in resisting the infections by bacteria D. zeae.

Genome-wide identification and characterization of long non-coding RNAs responsive to Dickeya zeae in rice

Plant long non-coding RNA (lncRNA) is a type of newly emerging epigenetic regulator playing a critical role in plant growth, development, and biotic stress responses. However, it is unknown whether lncRNAs are involved in resistance responses between rice and Dickeya zeae, a bacterial agent causing rice foot rot disease. In this study, RNA-seq was performed to uncover the co-expression regulating networks mediated by D. zeae responsive lncRNAs and their candidate target genes. Of the 4709 lncRNAs identified, 2518 and 2191 were up- and down-regulated in response to D. zeae infection, respectively. Expression changes of 17 selected lncRNAs and their predicted targets with a potential role in defense response were investigated by qPCR. The expression levels of five lncRNAs were up-regulated while their cognate candidate target genes were down-regulated upon D. zeae infection. In addition, several lncRNAs were predicted to be target mimics of osa-miR396 and osa-miR156. These results suggest that lncRNAs might play a role in response to D. zeae infection by regulating the transcript levels of their targets and miRNAs in rice.

GROWTH-REGULATING FACTOR 9 negatively regulates arabidopsis leaf growth by controlling ORG3 and restricting cell proliferation in leaf primordia

Leaf growth is a complex process that involves the action of diverse transcription factors (TFs) and their downstream gene regulatory networks. In this study, we focus on the functional characterization of the Arabidopsis thaliana TF GROWTH-REGULATING FACTOR9 (GRF9) and demonstrate that it exerts its negative effect on leaf growth by activating expression of the bZIP TF OBP3-RESPONSIVE GENE 3 (ORG3). While grf9 knockout mutants produce bigger incipient leaf primordia at the shoot apex, rosette leaves and petals than the wild type, the sizes of those organs are reduced in plants overexpressing GRF9 (GRF9ox). Cell measurements demonstrate that changes in leaf size result from alterations in cell numbers rather than cell sizes. Kinematic analysis and 5-ethynyl-2'-deoxyuridine (EdU) incorporation assay revealed that GRF9 restricts cell proliferation in the early developing leaf. Performing in vitro binding site selection, we identified the 6-base motif 5'-CTGACA-3' as the core binding site of GRF9. By global transcriptome profiling, electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation (ChIP) we identified ORG3 as a direct downstream, and positively regulated target of GRF9. Genetic analysis of grf9 org3 and GRF9ox org3 double mutants reveals that both transcription factors act in a regulatory cascade to control the final leaf dimensions by restricting cell number in the developing leaf.

Bacillus amyloliquefaciens FZB42 represses plant miR846 to induce systemic resistance via a jasmonic acid-dependent signalling pathway

Bacillus amyloliquefaciens FZB42 is a type of plant growth-promoting rhizobacterium (PGPR) which activates induced systemic resistance (ISR) in Arabidopsis. Blocking of the synthesis of cyclic lipopeptides and 2,3-butanediol by FZB42, which have been demonstrated to be involved in the priming of ISR, results in the abolishment of the plant defence responses. To further clarify the ISR activated by PGPRs at the microRNA (miRNA) level, small RNA (sRNA) libraries from Arabidopsis leaves after root irrigation with FZB42, FZB42ΔsfpΔalsS and control were constructed and sequenced. After fold change selection, promoter analysis and target prediction, miR846-5p and miR846-3p from the same precursor were selected as candidate ISR-associated miRNAs. miR846 belongs to the non-conserved miRNAs, specifically exists in Arabidopsis and its function in the plant defence response remains unclear. The disease severity of transgenic Arabidopsis overexpressing miR846 (OEmiR846) or knockdown miR846 (STTM846) against Pseudomonas syringae DC3000 suggests that the miR846 expression level in Arabidopsis is negatively correlated with disease resistance. Moreover, miR846 in Arabidopsis Col-0 is repressed after methyl jasmonate treatment. In addition, jasmonic acid (JA) signalling-related genes are up-regulated in STTM846, and the stomatal apertures of STTM846 are also less than those in Arabidopsis Col-0 after methyl jasmonate treatment. Furthermore, the disease resistance of STTM846 transgenic Arabidopsis against Pseudomonas syringae pv. tomato DC3000 (Pst DC3000) is blocked by the addition of the JA biosynthetic inhibitor diethyldiethiocarbamic acid (DIECA). Taken together, our results suggest that B. amyloliquefaciens FZB42 inoculation suppresses miR846 expression to induce Arabidopsis systemic resistance via a JA-dependent signalling pathway.

miR396-OsGRFs Module Balances Growth and Rice Blast Disease-Resistance

Fitness cost is a common phenomenon in rice blast disease-resistance breeding. MiR396 is a highly conserved microRNA (miRNA) family targeting Growth Regulating Factor (OsGRF) genes. Mutation at the target site of miR396 in certain OsGRF gene or blocking miR396 expression leads to increased grain yield. Here we demonstrated that fitness cost can be trade-off in miR396-OsGRFs module via balancing growth and immunity against the blast fungus. The accumulation of miR396 isoforms was significantly increased in a susceptible accession, but fluctuated in a resistant accession upon infection of Magnaporthe oryzae. The transgenic lines over-expressing different miR396 isoforms were highly susceptible to M. oryzae. In contrast, overexpressing target mimicry of miR396 to block its function led to enhanced resistance to M. oryzae in addition to improved yield traits. Moreover, transgenic plants overexpressing OsGRF6, OsGRF7, OsGRF8, and OsGRF9 exhibited enhanced resistance to M. oryzae, but showed different alteration of growth. While overexpression of OsGRF7 led to defects in growth, overexpression of OsGRF6, OsGRF8, and OsGRF9 resulted in better or no significant change of yield traits. Collectively, our results indicate that miR396 negatively regulates rice blast disease- resistance via suppressing multiple OsGRFs, which in turn differentially control growth and yield. Therefore, miR396-OsGRFs could be a potential module to demolish fitness cost in rice blast disease-resistance breeding.

MicroRNA396a-5p and -3p induce tomato disease susceptibility by suppressing target genes and upregulating salicylic acid

Plants have evolved a variety of mechanisms to perceive and resist the assault of pathogens. The biotrophs, necrotrophs and hemibiotrophs are types of plant pathogens that activate diverse salicylic acid (SA) and jasmonic acid (JA) signaling pathways. In this study we showed that the expressions of miR396a-5p and -3p in Solanum lycopersicum (S. lycopersicum) were both down-regulated after infection by hemibiotroph Phytophthora infestans (P. infestans) and necrotroph Botrytis cinerea (B. cinerea) infection. Overexpression of miR396a-5p and -3p in transgenic tomato enhanced the susceptibility of S. lycopersicum to P. infestans and B. cinerea infection and the tendency to produce reactive oxygen species (ROS) under pathogen-related biotic stress. Additionally, miR396a regulated growth-regulating factor1 (GRF1), salicylic acid carboxyl methyltransferase (SAMT), glycosyl hydrolases (GH) and nucleotide-binding site-leucine-rich repeat (NBS-LRR) and down-regulated their levels. This ultimately led to inhibition of the expression of pathogenesis-related 1 (PR1), TGA transcription factors1 and 2 (TGA1 and TGA2) and JA-dependent proteinase inhibitors I and II (PI I and II), but enhanced the endogenous SA content and nonexpressor of pathogenesis-related genes 1 (NPR1) expression. Taken together, our results showed that negative regulation of target genes and their downstream genes expressions by miR396a-5p and -3p are critical for tomato abiotic stresses via affecting SA or JA signaling pathways.

Cooperative Regulatory Functions of miR858 and MYB83 during Cyst Nematode Parasitism

MicroRNAs (miRNAs) recently have been established as key regulators of transcriptome reprogramming that define cell function and identity. Nevertheless, the molecular functions of the greatest number of miRNA genes remain to be determined. Here, we report cooperative regulatory functions of miR858 and its MYB83 transcription factor target gene in transcriptome reprogramming during Heterodera cyst nematode parasitism of Arabidopsis (Arabidopsis thaliana). Gene expression analyses and promoter-GUS fusion assays documented a role of miR858 in posttranscriptional regulation of MYB83 in the Heterodera schachtii-induced feeding sites, the syncytia. Constitutive overexpression of miR858 interfered with H. schachtii parasitism of Arabidopsis, leading to reduced susceptibility, while reduced miR858 abundance enhanced plant susceptibility. Similarly, MYB83 expression increases were conducive to nematode infection because overexpression of a noncleavable coding sequence of MYB83 significantly increased plant susceptibility, whereas a myb83 mutation rendered the plants less susceptible. In addition, RNA-seq analysis revealed that genes involved in hormone signaling pathways, defense response, glucosinolate biosynthesis, cell wall modification, sugar transport, and transcriptional control are the key etiological factors by which MYB83 facilitates nematode parasitism of Arabidopsis. Furthermore, we discovered that miR858-mediated silencing of MYB83 is tightly regulated through a feedback loop that might contribute to fine-tuning the expression of more than a thousand of MYB83-regulated genes in the H. schachtii-induced syncytium. Together, our results suggest a role of the miR858-MYB83 regulatory system in finely balancing gene expression patterns during H. schachtii parasitism of Arabidopsis to ensure optimal cellular function.

Molecular Characterization and Expression Profiling of Tomato GRF Transcription Factor Family Genes in Response to Abiotic Stresses and Phytohormones

Growth regulating factors (GRFs) are plant-specific transcription factors that are involved in diverse biological and physiological processes, such as growth, development and stress and hormone responses. However, the roles of GRFs in vegetative and reproductive growth, development and stress responses in tomato (Solanum lycopersicum) have not been extensively explored. In this study, we characterized the 13 SlGRF genes. In silico analysis of protein motif organization, intron–exon distribution, and phylogenetic classification confirmed the presence of GRF proteins in tomato. The tissue-specific expression analysis revealed that most of the SlGRF genes were preferentially expressed in young and growing tissues such as flower buds and meristems, suggesting that SlGRFs are important during growth and development of these tissues. Some of the SlGRF genes were preferentially expressed in fruits at distinct developmental stages suggesting their involvement in fruit development and the ripening process. The strong and differential expression of different SlGRFs under NaCl, drought, heat, cold, abscisic acid (ABA), and jasmonic acid (JA) treatment, predict possible functions for these genes in stress responses in addition to their growth regulatory functions. Further, differential expression of SlGRF genes upon gibberellic acid (GA3) treatment indicates their probable function in flower development and stress responses through a gibberellic acid (GA)-mediated pathway. The results of this study provide a basis for further functional analysis and characterization of this important gene family in tomato.

The functions of plant small RNAs in development and in stress responses

Like metazoans, plants use small regulatory RNAs (sRNAs) to direct gene expression. Several classes of sRNAs, which are distinguished by their origin and biogenesis, exist in plants. Among them, microRNAs (miRNAs) and trans-acting small interfering RNAs (ta-siRNAs) mainly inhibit gene expression at post-transcriptional levels. In the past decades, plant miRNAs and ta-siRNAs have been shown to be essential for numerous developmental processes, including growth and development of shoots, leaves, flowers, roots and seeds, among others. In addition, miRNAs and ta-siRNAs are also involved in the plant responses to abiotic and biotic stresses, such as drought, temperature, salinity, nutrient deprivation, bacteria, virus and others. This review summarizes the roles of miRNAs and ta-siRNAs in plant physiology and development.

The Arabidopsis miR396 mediates pathogen-associated molecular pattern-triggered immune responses against fungal pathogens

MicroRNAs (miRNAs) play a pivotal role in regulating gene expression during plant development. Although a substantial fraction of plant miRNAs has proven responsive to pathogen infection, their role in disease resistance remains largely unknown, especially during fungal infections. In this study, we screened Arabidopsis thaliana lines in which miRNA activity has been reduced using artificial miRNA target mimics (MIM lines) for their response to fungal pathogens. Reduced activity of miR396 (MIM396 plants) was found to confer broad resistance to necrotrophic and hemibiotrophic fungal pathogens. MiR396 levels gradually decreased during fungal infection, thus, enabling its GRF (GROWTH-REGULATING FACTOR) transcription factor target genes to trigger host reprogramming. Pathogen resistance in MIM396 plants is based on a superactivation of defense responses consistent with a priming event during pathogen infection. Notably, low levels of miR396 are not translated in developmental defects in absence of pathogen challenge. Our findings support a role of miR396 in regulating plant immunity, and broaden our knowledge about the molecular players and processes that sustain defense priming. That miR396 modulates innate immunity without growth costs also suggests fine-tuning of miR396 levels as an effective biotechnological means for protection against pathogen infection.

Arabidopsis miR827 mediates post-transcriptional gene silencing of its ubiquitin E3 ligase target gene in the syncytium of the cyst nematode Heterodera schachtii to enhance susceptibility

MicroRNAs (miRNAs) are a major class of small non-coding RNAs with emerging functions in biotic and abiotic interactions. Here, we report on a new functional role of Arabidopsis miR827 and its NITROGEN LIMITATION ADAPTATION (NLA) target gene in mediating plant susceptibility to the beet cyst nematode Heterodera schachtii. Cyst nematodes are sedentary endoparasites that induce the formation of multinucleated feeding structures termed syncytia in the roots of host plants. Using promoter:GUS fusion assays we established that miR827 was activated in the initial feeding cells and this activation was maintained in the syncytium during all sedentary stages of nematode development. Meanwhile, the NLA target gene, which encodes an ubiquitin E3 ligase enzyme, was post-transcriptionally silenced in the syncytium to permanently suppress its activity during all nematode parasitic stages. Overexpression of miR827 in Arabidopsis resulted in hyper-susceptibility to H. schachtii. In contrast, inactivation of miR827 activity through target mimicry or by overexpression a miR827-resistant cDNA of NLA produced the opposite phenotype of reduced plant susceptibility to H. schachtii. Gene expression analysis of several pathogenesis-related genes together with Agrobacterium-mediated transient expression in Nicotiana benthamiana provided strong evidence that miR827-mediated downregulation of NLA to suppress basal defense pathways. In addition, using yeast two-hybrid screens we identified several candidates of NLA-interacting proteins that are involved in a wide range of biological processes and molecular functions, including three pathogenesis-related proteins. Taken together, we conclude that nematode-activated miR827 in the syncytium is necessary to suppress immune responses in order to establish infection and cause disease.

MicroRNA393 is involved in nitrogen-promoted rice tillering through regulation of auxin signal transduction in axillary buds

Rice tillering has an important influence on grain yield, and is promoted by nitrogen (N) fertilizer. Several genes controlling rice tillering, which are regulated by poor N supply, have been identified. However, the molecular mechanism associated with the regulation of tillering based on N supply is poorly understood. Here, we report that rice microRNA393 (OsmiR393) is involved in N-mediated tillering by decreasing auxin signal sensitivity in axillary buds. Expression analysis showed that N fertilizer causes up-regulation of OsmiR393, but down-regulation of two target genes (OsAFB2 and OsTB1). In situ expression analysis showed that OsmiR393 is highly expressed in the lateral axillary meristem. OsmiR393 overexpression mimicked N-mediated tillering in wild type Zhonghua 11 (ZH11). Mutation of OsMIR393 in ZH11 repressed N-promoted tillering, which simulated the effects of limited N, and this could not be restored by supplying N fertilizer. Western blot analysis showed that OsIAA6 was accumulated in both OsmiR393-overexpressing lines and N-treated wild type rice, but was reduced in the OsMIR393 mutant. Therefore, we deduced that N-induced OsmiR393 accumulation reduces the expression of OsTIR1 and OsAFB2, which alleviates sensitivity to auxin in the axillary buds and stabilizes OsIAA6, thereby promoting rice tillering.

Transcriptomic Analysis of Tea Plant Responding to Drought Stress and Recovery

Tea plant (Camellia sinensis) is an economically important beverage crop. Drought stress (DS) seriously limits the growth and development of tea plant, thus affecting crop yield and quality. To elucidate the molecular mechanisms of tea plant responding to DS, we performed transcriptomic analysis of tea plant during the three stages [control (CK) and during DS, and recovery (RC) after DS] using RNA sequencing (RNA-Seq). Totally 378.08 million high-quality trimmed reads were obtained and assembled into 59,674 unigenes, which were extensively annotated. There were 5,955 differentially expressed genes (DEGs) among the three stages. Among them, 3,948 and 1,673 DEGs were up-regulated under DS and RC, respectively. RNA-Seq data were further confirmed by qRT-PCR analysis. Genes involved in abscisic acid (ABA), ethylene, and jasmonic acid biosynthesis and signaling were generally up-regulated under DS and down-regulated during RC. Tea plant potentially used an exchange pathway for biosynthesis of indole-3-acetic acid (IAA) and salicylic acid under DS. IAA signaling was possibly decreased under DS but increased after RC. Genes encoding enzymes involved in cytokinin synthesis were up-regulated under DS, but down-regulated during RC. It seemed probable that cytokinin signaling was slightly enhanced under DS. In total, 762 and 950 protein kinases belonging to 26 families were differentially expressed during DS and RC, respectively. Overall, 547 and 604 transcription factor (TF) genes belonging to 58 families were induced in the DS vs. CK and RC vs. DS libraries, respectively. Most members of the 12 TF families were up-regulated under DS. Under DS, genes related to starch synthesis were down-regulated, while those related to starch decomposition were up-regulated. Mannitol, trehalose and sucrose synthesis-related genes were up-regulated under DS. Proline was probably mainly biosynthesized from glutamate under DS and RC. The mechanism by which ABA regulated stomatal movement under DS and RC was partly clarified. These results document the global and novel responses of tea plant during DS and RC. These data will serve as a valuable resource for drought-tolerance research and will be useful for breeding drought-resistant tea cultivars.

Sp-miR396a-5p acts as a stress-responsive genes regulator by conferring tolerance to abiotic stresses and susceptibility to Phytophthora nicotianae infection in transgenic tobacco

KEY MESSAGE

Overexpression of Sp-miR396a-5p in tobacco increased tolerance to salt, drought, cold stress and susceptibility to Phytophthora nicotianae infection. MicroRNA396 (miR396) is one of the conserved microRNA families in plants, and it targeted growth-regulating factors (GRFs) family. The GRF transcription factors are associated with growth and stress responses. However, the molecular mechanisms of miR396 responding to environmental stresses are elusive. The purpose of this study was to explore the function of tomato miR396a-5p (Sp-miR396a-5p) in Solanaceae responses to abiotic and biotic stresses. We showed that Sp-miR396a-5p transcript levels were up-regulated under salt and drought stresses and down-regulated after Phytophthora infestans (P. infestans) infection. Consistently, overexpression of Sp-miR396a-5p in tobacco enhanced its tolerance to salt, drought and cold stresses. Additionally, the expression of Sp-miR396a-5p was found to be down-regulated under pathogen-related biotic stress. Tobacco plants overexpressing Sp-miR396a-5p showed increased susceptibility to Phytophthora nicotianae (P. nicotianae) infection. Physiological analysis indicated that Sp-miR396a-5p overexpression enhanced osmoregulation and decreased production of reactive oxygen species (ROS). Furthermore, four Sp-miR396a-5p target genes, NtGRF1, NtGRF3, NtGRF7 and NtGRF8, were down-regulated in these plants. Our results suggested that Sp-miR396a-5p plays critical roles in both abiotic stresses through targeting NtGRF7-regulated expression of osmotic stress-responsive genes and pathogen infection via the regulatory networks of NtGRF1 and NtGRF3.

Regulation of plant growth and development by the GROWTH-REGULATING FACTOR and GRF-INTERACTING FACTOR duo

Transcription factors are key regulators of gene expression and play pivotal roles in all aspects of living organisms. Therefore, identification and functional characterization of transcription factors is a prerequisite step toward understanding life. This article reviews molecular and biological functions of the two transcription regulator families, GROWTH-REGULATING FACTOR (GRF) and GRF-INTERACTING FACTOR (GIF), which have only recently been recognized. A myriad of experimental evidence clearly illustrates that GRF and GIF are bona fide partner proteins and form a plant-specific transcriptional complex. One of the most conspicuous outcomes from this research field is that the GRF-GIF duo endows the primordial cells of vegetative and reproductive organs with a meristematic specification state, guaranteeing the supply of cells for organogenesis and successful reproduction. It has recently been shown that GIF1 proteins, also known as ANGUSTIFOLIA3, recruit chromatin remodelling complexes to target genes, and that AtGRF expression is directly activated by the floral identity factors, APETALA1 and SEPALLATA3, providing an important insight into understanding of the action of GRF-GIF. Moreover, GRF genes are extensively subjected to post-transcriptional control by microRNA396, revealing the presence of a complex regulatory circuit in regulation of plant growth and development by the GRF-GIF duo.

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

Growth-regulating factors (GRFs) are plant-specific transcription factors that were originally identified for their roles in stem and leaf development, but recent studies highlight them to be similarly important for other central developmental processes including flower and seed formation, root development, and the coordination of growth processes under adverse environmental conditions. The expression of several GRFs is controlled by microRNA miR396, and the GRF-miRNA396 regulatory module appears to be central to several of these processes. In addition, transcription factors upstream of GRFs and miR396 have been discovered, and gradually downstream target genes of GRFs are being unraveled. Here, we review the current knowledge of the biological functions performed by GRFs and survey available molecular data to illustrate how they exert their roles at the cellular level.

Stress responsive miRNAs and isomiRs in cereals

Abiotic and biotic stress conditions are vital determinants in the production of cereals, the major caloric source in human nutrition. Small RNAs, miRNAs and isomiRs are central to post-transcriptional regulation of gene expression in a variety of cellular processes including development and stress responses. Several miRNAs have been identified using new technologies and have roles in stress responses in plants, including cereals. The overall knowledge about the cereal miRNA repertoire, as well as an understanding of complex miRNA mediated mechanisms of target regulation in response to stress conditions, is far from complete. Ongoing efforts that add to our understanding of complex miRNA machinery have implications in plant response to stress conditions. Additionally, sequence variants of miRNAs (isomiRNAs or isomiRs), regulation of their expression through dissection of upstream regulatory elements, the role of Processing-bodies (P-bodies) in miRNA exerted gene regulation and yet unveiled organellar plant miRNAs are newly emerging topics, which will contribute to the elucidation of the miRNA machinery and its role in cereal tolerance against abiotic and biotic stresses.

Lipopolysaccharide perception leads to dynamic alterations in the microtranscriptome of Arabidopsis thaliana cells and leaf tissues

BACKGROUND

MicroRNAs (miRNAs) are non-coding RNA molecules which have recently emerged as important gene regulators in plants and their gene expression analysis is becoming increasingly important. miRNAs regulate gene expression at the post-transcriptional level by translational repression or target degradation of specific mRNAs and gene silencing. In order to profile the microtranscriptome of Arabidopsis thaliana leaf and callus tissues in response to bacterial lipopolysaccharide (LPS), small RNA libraries were constructed at 0 and 3 h post induction with LPS and sequenced by Illumina sequencing technology.

RESULTS

Differential regulation of subset of miRNAs in response to LPS treament was observed. Small RNA reads were mapped to the miRNA database and 358 miRNAs belonging to 49 miRNA families in the callus tissues and 272 miRNAs belonging to 40 miRNA families in the leaf tissues were identified. Moreover, target genes for all the identified miRNAs families in the leaf tissues and 44 of the 49 miRNAs families in the callus tissues were predicted. The sequencing analysis showed that in both callus and leaf tissues, various stress regulated-miRNAs were differentially expressed and real time PCR validated the expression profile of miR156, miR158, miR159, miR169, miR393, miR398, miR399 and miR408 along with their target genes.

CONCLUSION

A. thaliana callus and leaf callus tissues respond to LPS as a microbe-associated molecular pattern molecule through dynamic changes to the microtranscriptome associated with differential transcriptional regulation in support of immunity and basal resistance.

GROWTH REGULATING FACTOR5 stimulates Arabidopsis chloroplast division, photosynthesis, and leaf longevity

Arabidopsis (Arabidopsis thaliana) leaf development relies on subsequent phases of cell proliferation and cell expansion. During the proliferation phase, chloroplasts need to divide extensively, and during the transition from cell proliferation to expansion, they differentiate into photosynthetically active chloroplasts, providing the plant with energy. The transcription factor GROWTH REGULATING FACTOR5 (GRF5) promotes the duration of the cell proliferation period during leaf development. Here, it is shown that GRF5 also stimulates chloroplast division, resulting in a higher chloroplast number per cell with a concomitant increase in chlorophyll levels in 35S:GRF5 leaves, which can sustain higher rates of photosynthesis. Moreover, 35S:GRF5 plants show delayed leaf senescence and are more tolerant for growth on nitrogen-depleted medium. Cytokinins also stimulate leaf growth in part by extending the cell proliferation phase, simultaneously delaying the onset of the cell expansion phase. In addition, cytokinins are known to be involved in chloroplast development, nitrogen signaling, and senescence. Evidence is provided that GRF5 and cytokinins synergistically enhance cell division and chlorophyll retention after dark-induced senescence, which suggests that they also cooperate to stimulate chloroplast division and nitrogen assimilation. Taken together with the increased leaf size, ectopic expression of GRF5 has great potential to improve plant productivity.

Techniques for transferring host-pathogen protein interactions knowledge to new tasks

We consider the problem of building a model to predict protein-protein interactions (PPIs) between the bacterial species Salmonella Typhimurium and the plant host Arabidopsis thaliana which is a host-pathogen pair for which no known PPIs are available. To achieve this, we present approaches, which use homology and statistical learning methods called “transfer learning.” In the transfer learning setting, the task of predicting PPIs between Arabidopsis and its pathogen S. Typhimurium is called the “target task.” The presented approaches utilize labeled data i.e., known PPIs of other host-pathogen pairs (we call these PPIs the “source tasks”). The homology based approaches use heuristics based on biological intuition to predict PPIs. The transfer learning methods use the similarity of the PPIs from the source tasks to the target task to build a model. For a quantitative evaluation we consider Salmonella-mouse PPI prediction and some other host-pathogen tasks where known PPIs exist. We use metrics such as precision and recall and our results show that our methods perform well on the target task in various transfer settings. We present a brief qualitative analysis of the Arabidopsis-Salmonella predicted interactions. We filter the predictions from all approaches using Gene Ontology term enrichment and only those interactions involving Salmonella effectors. Thereby we observe that Arabidopsis proteins involved e.g., in transcriptional regulation, hormone mediated signaling and defense response may be affected by Salmonella.

Measuring information flow in cellular networks by the systems biology method through microarray data

In general, it is very difficult to measure the information flow in a cellular network directly. In this study, based on an information flow model and microarray data, we measured the information flow in cellular networks indirectly by using a systems biology method. First, we used a recursive least square parameter estimation algorithm to identify the system parameters of coupling signal transduction pathways and the cellular gene regulatory network (GRN). Then, based on the identified parameters and systems theory, we estimated the signal transductivities of the coupling signal transduction pathways from the extracellular signals to each downstream protein and the information transductivities of the GRN between transcription factors in response to environmental events. According to the proposed method, the information flow, which is characterized by signal transductivity in coupling signaling pathways and information transductivity in the GRN, can be estimated by microarray temporal data or microarray sample data. It can also be estimated by other high-throughput data such as next-generation sequencing or proteomic data. Finally, the information flows of the signal transduction pathways and the GRN in leukemia cancer cells and non-leukemia normal cells were also measured to analyze the systematic dysfunction in this cancer from microarray sample data. The results show that the signal transductivities of signal transduction pathways change substantially from normal cells to leukemia cancer cells.