TaEPFL1, an EPIDERMAL PATTERNING FACTOR-LIKE (EPFL) secreted peptide gene, is required for stamen development in wheat

StEPF2 and StEPFL9 Play Opposing Roles in Regulating Stomatal Development and Drought Tolerance in Potato (Solanum tuberosum L.)

Stomata are essential for photosynthesis and water-use efficiency in plants. When expressed in transgenic Arabidopsis thaliana plants, the potato (Solanum tuberosum) proteins EPIDERMAL PATTERNING FACTOR 2 (StEPF2) and StEPF-LIKE9 (StEPFL9) play antagonistic roles in regulating stomatal density. Little is known, however, about how these proteins regulate stomatal development, growth, and response to water deficit in potato. Transgenic potato plants overexpressing StEPF2 (E2 plants) or StEPFL9 (ST plants) were generated, and RT-PCR and Western blot analyses were used to select two lines overexpressing each gene. E2 plants showed reduced stomatal density, whereas ST plants produced excessive stomata. Under well-watered conditions, ST plants displayed vigorous growth with improved leaf gas exchange and also showed increased biomass/yields compared with non-transgenic and E2 plants. E2 plants maintained lower H2O2 content and higher levels of stomatal conductance and photosynthetic capacity than non-transgenic and ST plants, which resulted in higher water-use efficiency and biomass/yields during water restriction. These results suggest that StEPF2 and StEPFL9 functioned in pathways regulating stomatal development. These genes are thus promising candidates for use in future breeding programs aimed at increasing potato water-use efficiency and yield under climate change scenarios.

Comparative Genomics Analysis of the Populus Epidermal Pattern Factor (EPF) Family Revealed Their Regulatory Effects in Populus euphratica Stomatal Development

Drought stress seriously threatens plant growth. The improvement of plant water use efficiency (WUE) and drought tolerance through stomatal regulation is an effective strategy for coping with water shortages. Epidermal patterning factor (EPF)/EPF-like (EPFL) family proteins regulate stomatal formation and development in plants and thus contribute to plant stress adaptation. Here, our analysis revealed the presence of 14 PeEPF members in the Populus euphratica genome, which exhibited a relatively conserved gene structure with 1-3 introns. Subcellular localisation prediction revealed that 9 PeEPF members were distributed in the chloroplasts of P. euphratica, and 5 were located extracellularly. Phylogenetic analysis indicated that PeEPFs can be divided into three clades, with genes within the same clade revealing a relatively conserved structure. Furthermore, we observed the evolutionary conservation of PeEPFs and AtEPF/EPFLs in certain domains, which suggests their conserved function. The analysis of cis-acting elements suggested the possible involvement of PeEPFs in plant response to multiple hormones. Transcriptomic analysis revealed considerable changes in the expression level of PeEPFs during treatment with polyethylene glycol and abscisic acid. The overexpression of PeEPF2 resulted in low stomatal density in transgenetic lines. These findings provide a basis for gaining insights into the function of PeEPFs in response to abiotic stress.

Genome-wide identification and expression analysis of EPF/EPFL gene family in Populus trichocarpa

The Epidermal Patterning Factor/EPF-like (EPF/EPFL) family encodes a specific type of secreted protein in plants and plays an important role in plant growth and development, especially in the process of morphogenesis. To investigate the characteristics of EPF/EPFL gene family members and their regulatory functions in stomatal development of Populus trichocarpa, a total of 15 EPF/EPFL family genes were identified. Then the gene structure, chromosome location, phylogenetic relationship, protein conserved domain and gene expression profile were analyzed. According to phylogenetic analysis, PtEPF/EPFL can be classified into four groups. The gene structure and protein conservation motifs within the EPF family indicate the high conservation of the PtEPF/EPFL sequence. The promoter region of PtEPF/EPFL was found to contain cis-elements in response to stress and plant hormones. In addition, RT-qPCR results indicated that the PtEPF/EPFL have a differentially expressed in different tissues. Under drought stress treatment, a substantial upregulation was observed in the majority of PtEPF/EPFL members, suggesting their potential involvement in drought response. These results provide a theoretical basis for future exploration of the characteristics and functions of more PtEPF/EPFL genes.

Genome-wide analyses of member identification, expression pattern, and protein-protein interaction of EPF/EPFL gene family in Gossypium

BACKGROUND

Epidermal patterning factor / -like (EPF/EPFL) gene family encodes a class of cysteine-rich secretory peptides, which are widelyfound in terrestrial plants.Multiple studies has indicated that EPF/EPFLs might play significant roles in coordinating plant development and growth, especially as the morphogenesis processes of stoma, awn, stamen, and fruit skin. However, few research on EPF/EPFL gene family was reported in Gossypium.

RESULTS

We separately identified 20 G. raimondii, 24 G. arboreum, 44 G. hirsutum, and 44 G. barbadense EPF/EPFL genes in the 4 representative cotton species, which were divided into four clades together with 11 Arabidopsis thaliana, 13 Oryza sativa, and 17 Selaginella moellendorffii ones based on their evolutionary relationships. The similar gene structure and common motifs indicated the high conservation among the EPF/EPFL members, while the uneven distribution in chromosomes implied the variability during the long-term evolutionary process. Hundreds of collinearity relationships were identified from the pairwise comparisons of intraspecifc and interspecific genomes, which illustrated gene duplication might contribute to the expansion of cotton EPF/EPFL gene family. A total of 15 kinds of cis-regulatory elements were predicted in the promoter regions, and divided into three major categories relevant to the biological processes of development and growth, plant hormone response, and abiotic stress response. Having performing the expression pattern analyses with the basic of the published RNA-seq data, we found most of GhEPF/EPFL and GbEPF/EPFL genes presented the relatively low expression levels among the 9 tissues or organs, while showed more dramatically different responses to high/low temperature and salt or drought stresses. Combined with transcriptome data of developing ovules and fibers and quantitative Real-time PCR results (qRT-PCR) of 15 highly expressed GhEPF/EPFL genes, it could be deduced that the cotton EPF/EPFL genes were closely related with fiber development. Additionally, the networks of protein-protein interacting among EPF/EPFLs concentrated on the cores of GhEPF1 and GhEPF7, and thosefunctional enrichment analyses indicated that most of EPF/EPFLs participate in the GO (Gene Ontology) terms of stomatal development and plant epidermis development, and the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways of DNA or base excision repair.

CONCLUSION

Totally, 132 EPF/EPFL genes were identified for the first time in cotton, whose bioinformatic analyses of cis-regulatory elements and expression patterns combined with qRT-PCR experiments to prove the potential functions in the biological processes of plant growth and responding to abiotic stresses, specifically in the fiber development. These results not only provide comprehensive and valuable information for cotton EPF/EPFL gene family, but also lay solid foundation for screening candidate EPF/EPFL genes in further cotton breeding.

Genome-wide identification and phylogenetic and expression pattern analyses of EPF/EPFL family genes in the Rye (Secale cereale L.)

Rye (Secale cereale L.) is one of the major cereal crop species in the Triticeae family and is known to be most tolerant to diverse abiotic stresses, such as cold, heat, osmotic, and salt stress. The EPIDERMAL PATTERNING FACTOR (EPF) and EPF-LIKE (EPFL) families of small secreted peptides act to regulate many aspects of plant growth and development; however, their functions are not widely characterized in rye. In this study, we identified 12 ScEPF/EPFL genes, which can be divided into six groups and are evenly distributed on six rye chromosomes. Further examination of the gene structure and protein conservation motifs of EPF/EPFL family members demonstrated the high conservation of the ScEPF/EPFL sequence. Interactions between ScEPF/EPFL proteins and promoters containing hormone- and stress-responsive cis-acting elements suggest that the regulation of ScEPF/EPFL expression is complex. Expression profiling analyses revealed that ScEPF/EPFL genes exhibited tissue-specific expression patterns. Notably, ScEPFL1,ScEPFL7, ScEPFL9, and ScEPFL10 displayed significantly higher expression levels in spikelets compared to other tissues. Moreover, fluorescence quantification experiments demonstrated that these genes exhibited distinct expression patterns in response to various stress conditions, suggesting that each gene plays a unique role in stress signaling pathways. Our research findings provide a solid basis for further investigation into the functions of ScEPF/EPFLs. Furthermore, these genes can serve as potential candidates for breeding stress-resistant rye varieties and improving production yields.

A transcriptional atlas identifies key regulators and networks for the development of spike tissues in barley

Grain number and size determine grain yield in crops and are closely associated with spikelet fertility and grain filling in barley (Hordeum vulgare). Abortion of spikelet primordia within individual barley spikes causes a 30%-50% loss in the potential number of grains during development from the awn primordium stage to the tipping stage, after that grain filling is the primary factor regulating grain size. To identify transcriptional signatures associated with spike development, we use a six-rowed barley cultivar (Morex) to develop a spatiotemporal transcriptome atlas containing 255 samples covering 17 stages and 5 positions along the spike. We identify several fundamental regulatory networks, in addition to key regulators of spike development and morphology. Specifically, we show HvGELP96, encoding a GDSL domain-containing protein, as a regulator of spikelet fertility and grain number. Our transcriptional atlas offers a powerful resource to answer fundamental questions in spikelet development and degeneration in barley.

Sex-Related Gene Network Revealed by Transcriptome Differentiation of Bisexual and Unisexual Flowers of Orchid Cymbidium tortisepalum

Despite extensive research on orchid reproductive strategies, the genetic studies of sex differentiation in the orchid family are still lacking. In this study, we compared three sexual phenotypes of Cymbidium tortisepalum bisexual flowers as well as female and male unisexual mutants. Through comparative transcriptomes, we analyzed the sex-biased differentially expressed genes (DEGs) and gene co-expression networks of sex organs (gynostemium and ovary) among them, identified the candidate genes of sex differentiation, and validated their expression by qRT-PCR. The C. tortisepalum unisexual mutants with degenerated phenotypes were compared to the bisexual plants with respect to both the flower organs and plant morphologies. Totally, 12,145, 10,789, and 14,447 genes were uniquely expressed in the female, male, and hermaphrodite sex organs, respectively. A total of 4291 sex-biased DEGs were detected among them, with 871, 2867, and 1937 DEGs in the comparisons of bisexual vs. female, bisexual vs. male, and male vs. female flowers, respectively. Two co-expressed network modules, with 81 and 419 genes were tightly correlated with female sexual traits, while two others with 265 and 135 genes were highly correlated with male sexual traits. Two female-biased hub genes (CtSDR3b and CtSDR3b-like) nested in the female modules, the homologs of maize sex determinant tasselseed2, may control the feminization of C. tortisepalum. At the same time, two male-biased hub genes (CtYAB2 and CtYAB5) nested in the male modules, the homologs of grape sex determinant VviYABBY3, may control the androphany of C. tortisepalum. This study discovered the molecular regulation networks and proposed a model for orchid sex differentiation, therefore providing for the first time the genetic basis for the sex separation in the orchid family.

Genetic Dissection of Spike Productivity Traits in the Siberian Collection of Spring Barley

Barley (Hordeum vulgare L.) is one of the most commonly cultivated cereals worldwide. Its local varieties can represent a valuable source of unique genetic variants useful for crop improvement. The aim of this study was to reveal loci contributing to spike productivity traits in Siberian spring barley and to develop diagnostic DNA markers for marker-assisted breeding programs. For this purpose we conducted a genome-wide association study using a panel of 94 barley varieties. In total, 64 SNPs significantly associated with productivity traits were revealed. Twenty-three SNP markers were validated by genotyping in an independent sample set using competitive allele-specific PCR (KASP). Finally, fourteen markers associated with spike productivity traits on chromosomes 2H, 4H and 5H can be suggested for use in breeding programs.

The TaTCP4/10-B1 cascade regulates awn elongation in wheat (Triticum aestivum L.)

Awns are important morphological markers for wheat and exert a strong physiological effect on wheat yield. The awn elongation suppressor B1 has recently been cloned through association and linkage analysis in wheat. However, the mechanism of awn inhibition centered around B1 remains to be clarified. Here, we identified an allelic variant in the coding region of B1 through analysis of re-sequencing data; this variant causes an amino acid substitution and premature termination, resulting in a long-awn phenotype. Transcriptome analysis indicated that B1 inhibited awn elongation by impeding cytokinin- and auxin-promoted cell division. Moreover, B1 directly repressed the expression of TaRAE2 and TaLks2, whose orthologs have been reported to promote awn development in rice or barley. More importantly, we found that TaTCP4 and TaTCP10 synergistically inhibited the expression of B1, and a G-to-A mutation in the B1 promoter attenuated its inhibition by TaTCP4/10. Taken together, our results reveal novel mechanisms of awn development and provide genetic resources for trait improvement in wheat.

Small EPIDERMAL PATTERNING FACTOR-LIKE2 peptides regulate awn development in rice

The EPIDERMAL PATTERNING FACTOR (EPF) and EPF-LIKE (EPFL) family of small secreted peptides act to regulate many aspects of plant growth and development; however, their functions are not widely characterized in rice (Oryza sativa). Here, we used clustered regularly interspaced short palindromic repeats (CRISPR/Cas9) technology to individually knockout each of 11 EPF/EPFL genes in the rice cultivar Kasalath. Loss of function of most OsEPF/EPFL genes generated no obvious phenotype alteration, while disruption of OsEPFL2 in Kasalath caused a short or no awn phenotype and reduced grain size. OsEPFL2 is strongly expressed in the young panicle, consistent with a role in regulating awn and grain development. Haplotype analysis indicated that OsEPFL2 can be classified into six major haplotypes. Nucleotide diversity and genetic differentiation analyses suggested that OsEPFL2 was positively selected during the domestication of rice. Our work to systematically investigate the function of EPF/EPFL peptides demonstrates that different members of the same gene family have been independently selected for their ability to regulate a similar biological function and provides perspective on rice domestication.

Co-expression network analysis of genes and networks associated with wheat pistillody

Crop male sterility has great value in theoretical research and breeding application. HTS-1, whose stamens transformed into pistils or pistil-like structures, is an important male sterility material selecting from Chinese Spring three-pistil (CSTP) wheat. However the molecular mechanism of pistillody development in HTS-1 remains a mystery. RNA-seq data of 11 wheat tissues were obtained from the National Center for Biotechnology Information (NCBI), including the stamens of CSTP and the pistils and pistillodic stamen of HTS-1. The Salmon program was utilized to quantify the gene expression levels of the 11 wheat tissues; and gene quantification results were normalized by transcripts per million (TPM). In total, 58,576 genes were used to construct block-wise network by co-expression networks analysis (WGCNA) R package. We obtained all of modules significantly associated with the 11 wheat tissues. AgriGO V2.0 was used to do Gene Ontology (GO) enrichment analysis; and genes and transcription factors (TFs) in these significant modules about wheat pistillody development were identified from GO enrichment results. Basic local alignment search tool (BLAST) was used to align HTS-1 proteins with the published pistillody-related proteins and TFs. Genes about wheat pistillody development were analyzed and validated by qRT-PCR. The MEturquoise, MEsaddlebrown, MEplum, MEcoral1, MElightsteelblue1, and MEdarkslateblue modules were significantly corelated to pistillodic stamen (correlation p < 0.05). Moreover, 206 genes related to carpel development (GO:0048440) or gynoecium development (GO:0048467) were identified only in the MEturquoise module by Gene Ontology (GO) analysis, and 42 of 206 genes were hub genes in MEturquoise module. qRT-PCR results showed that 38 of the 42 hub genes had highly expressed in pistils and pistillodic stamens than in stamens. A total of 15 pistillody development-related proteins were validated by BLAST. Transcription factors (TFs) were also analyzed in the MEturquoise module, and 618 TFs were identified. In total, 56 TFs from 11 families were considered to regulate the development of pistillodic stamen. The co-expression network showed that six of HB and three of BES1 genes were identified in 42 hub genes. This indicated that TFs played important roles in wheat pistillody development. In addition, there were 11 of ethylene-related genes connected with TFs or hub genes, suggesting the important roles of ethylene-related genes in pistillody development. These results provide important insights into the molecular interactions underlying pistillody development.

Comparative transcriptome analysis indicates conversion of stamens into pistil-like structures in male sterile wheat (Triticum aestivum L.) with Aegilops crassa cytoplasm

BACKGROUND

Aegilops crassa cytoplasm is an important source for investigating cytoplasmic male sterility (CMS). Moreover, the stamens of line C303A exhibit a high degree of pistillody, turning almost white. However, the molecular mechanism that underlies pistillody in C303A remains unclear. Therefore, to obtain a better understanding of pistillody in C303A, the phenotypic and cytological features of C303A were observed to identify the key stage for the homeotic transformation of stamens into pistil-like structures. Transcriptome profiles were determined for stamens using Illumina RNA sequencing.

RESULTS

Morphological observations of the CMS wheat line with Aegilops crassa cytoplasm C303A showed that the pistils developed normally, but the stamens were ultimately aborted and they released no pollen when mature. According to paraffin section observations, the stamens began to transform into pistils or pistil-like structures in the binucleate stage (BNS). Therefore, the stamens were collected from line C303A and its maintainer 303B in the BNS for transcriptome sequencing. In total, 20,444 wheat genes were determined as differentially expressed in C303A and 303B stamens, with 10,283 upregulated and 10,161 downregulated genes. Gene Ontology enrichment analyses showed that most of the differentially expressed genes (DEGs) were annotated with GO terms comprising metabolic process, cell, cellular process, catalytic activity, and cell part. Analysis based on the Kyoto Encyclopedia of Genes and Genomes database showed that the enriched DEGs were mainly associated with energy metabolism. We also found several essential genes that may contribute to pistillody in C303A. These findings suggest that disrupted energy metabolism and reactive oxygen metabolism induce pistillody and eventually lead to abortion in C303A.

CONCLUSION

We determined the complex transcriptome profiles for C303A stamens and demonstrated that disrupted energy metabolism and class B MADS-box genes are related to pistillody. These findings may facilitate future studies of the mechanistic response of the wheat stamen and pollen development in CMS.