Functional analysis of B and C class floral organ genes in spinach demonstrates their role in sexual dimorphism

Gene Regulatory Network Controlling Flower Development in Spinach (Spinacia oleracea L.)

Spinach (Spinacia oleracea L.) is a dioecious, diploid, wind-pollinated crop cultivated worldwide. Sex determination plays an important role in spinach breeding. Hence, this study aimed to understand the differences in sexual differentiation and floral organ development of dioecious flowers, as well as the differences in the regulatory mechanisms of floral organ development of dioecious and monoecious flowers. We compared transcriptional-level differences between different genders and identified differentially expressed genes (DEGs) related to spinach floral development, as well as sex-biased genes to investigate the flower development mechanisms in spinach. In this study, 9189 DEGs were identified among the different genders. DEG analysis showed the participation of four main transcription factor families, MIKC_MADS, MYB, NAC, and bHLH, in spinach flower development. In our key findings, abscisic acid (ABA) and gibberellic acid (GA) signal transduction pathways play major roles in male flower development, while auxin regulates both male and female flower development. By constructing a gene regulatory network (GRN) for floral organ development, core transcription factors (TFs) controlling organ initiation and growth were discovered. This analysis of the development of female, male, and monoecious flowers in spinach provides new insights into the molecular mechanisms of floral organ development and sexual differentiation in dioecious and monoecious plants in spinach.

Gibberellins regulate masculinization through the SpGAI-SpSTM module in dioecious spinach

The sex of dioecious plants is mainly determined by genetic factors, but it can also be converted by environmental cues such as exogenous phytohormones. Gibberellic acids (GAs) are well-known inducers of flowering and sexual development, yet the pathway of gibberellin-induced sex conversion in dioecious spinach (Spinacia oleracea L.) remains elusive. Based on sex detection before and after GA3 application using T11A and SSR19 molecular markers, we confirmed and elevated the masculinization effect of GA on a single female plant through exogenous applications of GA3, showing complete conversion and functional stamens. Silencing of GIBBERELLIC ACID INSENSITIVE (SpGAI), a single DELLA family protein that is a central GA signaling repressor, results in similar masculinization. We also show that SpGAI can physically interact with the spinach KNOX transcription factor SHOOT MERISTEMLESS (SpSTM), which is a homolog of the flower meristem identity regulator STM in Arabidopsis. The silencing of SpSTM also masculinized female flowers in spinach. Furthermore, SpSTM could directly bind the intron of SpPI to repress SpPI expression in developing female flowers. Overall, our results suggest that GA induces a female masculinization process through the SpGAI-SpSTM-SpPI regulatory module in spinach. These insights may help to clarify the molecular mechanism underlying the sex conversion system in dioecious plants while also elucidating the physiological basis for the generation of unisexual flowers so as to establish dioecy in plants.

QTL analysis of femaleness in monoecious spinach and fine mapping of a major QTL using an updated version of chromosome-scale pseudomolecules

Although spinach is predominantly dioecious, monoecious plants with varying proportions of female and male flowers are also present. Recently, monoecious inbred lines with highly female and male conditions have been preferentially used as parents for F1-hybrids, rather than dioecious lines. Accordingly, identifying the loci for monoecism is an important issue for spinach breeding. We here used long-read sequencing and Hi-C technology to construct SOL_r2.0_pseudomolecule, a set of six pseudomolecules of spinach chromosomes (total length: 879.2 Mb; BUSCO complete 97.0%) that are longer and more genetically complete than our previous version of pseudomolecules (688.0 Mb; 81.5%). Three QTLs, qFem2.1, qFem3.1, and qFem6.1, responsible for monoecism were mapped to SOL_r2.0_pseudomolecule. qFem3.1 had the highest LOD score and corresponded to the M locus, which was previously identified as a determinant of monoecious expression, by genetic analysis of progeny from female and monoecious plants. The other QTLs were shown to modulate the ratio of female to male flowers in monoecious plants harboring a dominant allele of the M gene. Our findings will enable breeders to efficiently produce highly female- and male-monoecious parental lines for F1-hybrids by pyramiding the three QTLs. Through fine-mapping, we narrowed the candidate region for the M locus to a 19.5 kb interval containing three protein-coding genes and one long non-coding RNA gene. Among them, only RADIALIS-like-2a showed a higher expression in the reproductive organs, suggesting that it might play a role in reproductive organogenesis. However, there is no evidence that it is involved in the regulation of stamen and pistil initiation, which are directly related to the floral sex differentiation system in spinach. Given that auxin is involved in reproductive organ formation in many plant species, genes related to auxin transport/response, in addition to floral organ formation, were identified as candidates for regulators of floral sex-differentiation from qFem2.1 and qFem6.1.

Subgenome dominance shapes novel gene evolution in the decaploid pitcher plant Nepenthes gracilis

Subgenome dominance after whole-genome duplication generates distinction in gene number and expression at the level of chromosome sets, but it remains unclear how this process may be involved in evolutionary novelty. Here we generated a chromosome-scale genome assembly of the Asian pitcher plant Nepenthes gracilis to analyse how its novel traits (dioecy and carnivorous pitcher leaves) are linked to genomic evolution. We found a decaploid karyotype and a clear indication of subgenome dominance. A male-linked and pericentromerically located region on the putative sex chromosome was identified in a recessive subgenome and was found to harbour three transcription factors involved in flower and pollen development, including a likely neofunctionalized LEAFY duplicate. Transcriptomic and syntenic analyses of carnivory-related genes suggested that the paleopolyploidization events seeded genes that subsequently formed tandem clusters in recessive subgenomes with specific expression in the digestive zone of the pitcher, where specialized cells digest prey and absorb derived nutrients. A genome-scale analysis suggested that subgenome dominance likely contributed to evolutionary innovation by permitting recessive subgenomes to diversify functions of novel tissue-specific duplicates. Our results provide insight into how polyploidy can give rise to novel traits in divergent and successful high-ploidy lineages.

Ovule Transcriptome Analysis Discloses Deregulation of Genes and Pathways in Sexual and Apomictic Limonium Species (Plumbaginaceae)

The genus Limonium Mill. (sea lavenders) includes species with sexual and apomixis reproductive strategies, although the genes involved in these processes are unknown. To explore the mechanisms beyond these reproduction modes, transcriptome profiling of sexual, male sterile, and facultative apomictic species was carried out using ovules from different developmental stages. In total, 15,166 unigenes were found to be differentially expressed with apomictic vs. sexual reproduction, of which 4275 were uniquely annotated using an Arabidopsis thaliana database, with different regulations according to each stage and/or species compared. Gene ontology (GO) enrichment analysis indicated that genes related to tubulin, actin, the ubiquitin degradation process, reactive oxygen species scavenging, hormone signaling such as the ethylene signaling pathway and gibberellic acid-dependent signal, and transcription factors were found among differentially expressed genes (DEGs) between apomictic and sexual plants. We found that 24% of uniquely annotated DEGs were likely to be implicated in flower development, male sterility, pollen formation, pollen-stigma interactions, and pollen tube formation. The present study identifies candidate genes that are highly associated with distinct reproductive modes and sheds light on the molecular mechanisms of apomixis expression in Limonium sp.

Genome-wide identification and expression analysis of MADS-box transcription factors reveal their involvement in sex determination of hardy rubber tree (Eucommia ulmoides oliv.)

Eucommia ulmoides is a famous rubber-producing and medicinal tree species that produces unisexual flowers on separate individuals from the earliest stage of stamen/pistil primordium formation. To explore the genetic regulation pathway of sex in E. ulmoides, comprehensive genome-wide analyses and tissue-/sex-specific transcriptome comparisons of MADS-box transcription factors were performed for the first time in this work. Quantitative real-time PCR technique was employed to further validate the expression of genes that are assigned to floral organ ABCDE model. A total of 66 non-redundant E. ulmoides MADS-box (EuMADS) genes were identified, they were classified into Type I (M-type, 17 genes) and Type II (MIKC, 49 genes). Complex protein-motif composition, exon-intron structure and phytohormone-response cis-elements were detected in MIKC-EuMADS genes. Furthermore, 24 differentially-expressed EuMADS genes (DEGs) between male and female flowers, and two DEGs between male and female leaves were revealed. Amongst the 14 floral organ ABCDE model-related genes, there were 6 (A/B/C/E-class) and 5 (A/D/E-class) genes displayed male- and female-biased expression respectively. In particular, one B-class gene EuMADS39 and one A-class gene EuMADS65 were almost exclusively expressed in male trees, no matter in flower or leaf tissues. Collectively, these results suggested a critical role of MADS-box transcription factors in sex determination of E. ulmoides, which is conducive to decoding the molecular regulation mechanism of sex in E. ulmoides.

Genic male and female sterility in vegetable crops

Vegetable crops are greatly appreciated for their beneficial nutritional and health components. Hybrid seeds are widely used in vegetable crops for advantages such as high yield and improved resistance, which require the participation of male (stamen) and female (pistil) reproductive organs. Male- or female-sterile plants are commonly used for production of hybrid seeds or seedless fruits in vegetables. In this review we will focus on the types of genic male sterility and factors affecting female fertility, summarize typical gene function and research progress related to reproductive organ identity and sporophyte and gametophyte development in vegetable crops [mainly tomato (Solanum lycopersicum) and cucumber (Cucumis sativus)], and discuss the research trends and application perspectives of the sterile trait in vegetable breeding and hybrid production, in order to provide a reference for fertility-related germplasm innovation.

A convergent mechanism of sex determination in dioecious plants: Distinct sex-determining genes display converged regulation on floral B-class genes

Sex determination in dioecious plants has been broadly and progressively studied with the blooming of genome sequencing and editing techniques. This provides us with a great opportunity to explore the evolution and genetic mechanisms underlining the sex-determining system in dioecious plants. In this study, comprehensively reviewing advances in sex-chromosomes, sex-determining genes, and floral MADS-box genes in dioecious plants, we proposed a convergent model that governs plant dioecy across divergent species using a cascade regulation pathway connecting sex-determining genes and MADS-box genes e.g., B-class genes. We believe that this convergent mechanism of sex determination in dioecious plants will shed light on our understanding of gene regulation and evolution of plant dioecy. Perspectives concerning the evolutionary pathway of plant dioecy are also suggested.

Analysis of Floral Organ Development and Sex Determination in Schisandra chinensis by Scanning Electron Microscopy and RNA-Sequencing

S. chinensis is a typical monoecious plant, and the number and development of female flowers determines the yield of S. chinensis. Due to a lack of genetic information, the molecular mechanism of sex differentiation in S. chinensis remains unclear. In this study, the combination of scanning electron microscopy (SEM) and RNA sequencing (RNA-seq) was used to understand the way of sex differentiation of S. chinensis and to mine the related genes of sex determination. The result shows the development of male and female S. chinensis flowers was completed at the same time, the unisexual S. chinensis flowers did not undergo a transition stage between sexes, and sex may have been determined at an early stage in flower development. The results of the gene function analysis of the plant hormone signaling pathway and sucrose metabolism pathway suggest that auxin and JA could be the key hormones for sex differentiation in S. chinensis, and sucrose may promote pollen maturation at the later stage of male flower development. Two AGAMOUS (GAG) genes, 10 AGAMOUS-like MADS-box (AGLs) genes, and the MYB, NAC, WRKY, bHLH, and Trihelix transcription factor families may play important roles in sex determination in S. chinensis. Taken together, the present findings provide valuable genetic information on flower development and sex determination in S. chinensis.

The spinach YY genome reveals sex chromosome evolution, domestication, and introgression history of the species

Background

Spinach (Spinacia oleracea L.) is a dioecious species with an XY sex chromosome system, but its Y chromosome has not been fully characterized. Our knowledge about the history of its domestication and improvement remains limited.

Results

A high-quality YY genome of spinach is assembled into 952 Mb in six pseudo-chromosomes. By a combination of genetic mapping, Genome-Wide Association Studies, and genomic analysis, we characterize a 17.42-Mb sex determination region (SDR) on chromosome 1. The sex chromosomes of spinach evolved when an insertion containing sex determination genes occurred, followed by a large genomic inversion about 1.98 Mya. A subsequent burst of SDR-specific repeats (0.1–0.15 Mya) explains the large size of this SDR. We identify a Y-specific gene, NRT1/PTR 6.4 which resides in this insertion, as a strong candidate for the sex determination or differentiation factor. Resequencing of 112 spinach genomes reveals a severe domestication bottleneck approximately 10.87 Kya, which dates the domestication of spinach 7000 years earlier than the archeological record. We demonstrate that a strong selection signal associated with internode elongation and leaf area expansion is associated with domestication of edibility traits in spinach. We find that several strong genomic introgressions from the wild species Spinacia turkestanica and Spinacia tetrandra harbor desirable alleles of genes related to downy mildew resistance, frost resistance, leaf morphology, and flowering-time shift, which likely contribute to spinach improvement.

Conclusions

Analysis of the YY genome uncovers evolutionary forces shaping nascent sex chromosome evolution in spinach. Our findings provide novel insights about the domestication and improvement of spinach.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02633-x.

CmLHP1 proteins play a key role in plant development and sex determination in melon (Cucumis melo)

In monoecious melon (Cucumis melo), sex is determined by the differential expression of sex determination genes (SDGs) and adoption of sex-specific transcriptional programs. Histone modifications such as H3K27me3 have been previously shown to be a hallmark associated to unisexual flower development in melon; yet, no genetic approaches have been conducted for elucidating the roles of H3K27me3 writers, readers, and erasers in this process. Here we show that melon homologs to Arabidopsis LHP1, CmLHP1A and B, redundantly control several aspects of plant development, including sex expression. Cmlhp1ab double mutants displayed an overall loss and redistribution of H3K27me3, leading to a deregulation of genes involved in hormone responses, plant architecture, and flower development. Consequently, double mutants display pleiotropic phenotypes and, interestingly, a general increase of the male:female ratio. We associated this phenomenon with a general deregulation of some hormonal response genes and a local activation of male-promoting SDGs and MADS-box transcription factors. Altogether, these results reveal a novel function for CmLHP1 proteins in maintenance of monoecy and provide novel insights into the polycomb-mediated epigenomic regulation of sex lability in plants.

Development and Evolution of Unisexual Flowers: A Review

The development of unisexual flowers has been described in a large number of taxa, sampling the diversity of floral phenotypes and sexual systems observed in extant angiosperms, in studies focusing on floral ontogeny, on the evo-devo of unisexuality, or on the genetic and chromosomal bases of unisexuality. We review here such developmental studies, aiming at characterizing the diversity of ontogenic pathways leading to functionally unisexual flowers. In addition, we present for the first time and in a two-dimensional morphospace a quantitative description of the developmental rate of the sexual organs in functionally unisexual flowers, in a non-exhaustive sampling of angiosperms with contrasted floral morphologies. Eventually, recommendations are provided to help plant evo-devo researchers and botanists addressing macroevolutionary and ecological issues to more precisely select the taxa, the biological material, or the developmental stages to be investigated.

The underlying molecular conservation and diversification of dioecious flower and leaf buds provide insights into the development, dormancy breaking, flowering, and sex association of willows

As harbingers of bursting growth, flower buds and leaf buds generally show similar surface morphologies but different structural and functional changes. Dioecious plants further generate four types of Female/Male Flower/Leaf Buds (FFB, FLB, MFB, and MLB), showing a complex regulation. However, little is known about their underlying molecular mechanisms. Here, we exemplify the woody dioecious Salix linearistipularis to investigate their morphological characteristics and potential molecular mechanisms by combining cytological, physiological, phenological, and transcriptomic datasets. First, FFB and MFB have simultaneous development dynamics and so do FLB and MLB. Interestingly, FLB and MLB show very similar expression profiles preparing for photosynthesis and stress-tolerance, whereas FFB and MFB show great similarities but also striking sexual differences. Comparing flower buds and leaf buds after their revival from dormancy shows different cold- and vernalization-responsive genes (e.g. SliVRN1, SliAGL19, and SliAGL24), implying different programming processes for dormancy breaking between the buds. Moreover, except SliAP3, the expression of ABCDE model genes is consistent with their roles in the buds, suggesting a conserved mechanism of flower development between dioecious Salix and hermaphrodite Arabidopsis. Finally, considering sex-associated genes (e.g. SliCLE25, SliTPS21, and SliARR9) on Salix chromosomes and other reports, we hypothesize a dynamic model of sex determination on chromosomes 15 and 19 in the last ancestor of Salix and Populus but evolutionarily on 15 in Salix after their divergence. Together, our study provides new insights into the molecular mechanisms of dioecious four-type buds by showing the genes involved in their development, dormancy breaking, flowering, and sexual association.

The Dynamics of Flower Development in Castanea sativa Mill

The sweet chestnut tree (Castanea sativa Mill.) is one of the most significant Mediterranean tree species, being an important natural resource for the wood and fruit industries. It is a monoecious species, presenting unisexual male catkins and bisexual catkins, with the latter having distinct male and female flowers. Despite the importance of the sweet chestnut tree, little is known regarding the molecular mechanisms involved in the determination of sexual organ identity. Thus, the study of how the different flowers of C. sativa develop is fundamental to understand the reproductive success of this species and the impact of flower phenology on its productivity. In this study, a C. sativa de novo transcriptome was assembled and the homologous genes to those of the ABCDE model for floral organ identity were identified. Expression analysis showed that the C. sativa B- and C-class genes are differentially expressed in the male flowers and female flowers. Yeast two-hybrid analysis also suggested that changes in the canonical ABCDE protein-protein interactions may underlie the mechanisms necessary to the development of separate male and female flowers, as reported for the monoecious Fagaceae Quercus suber. The results here depicted constitute a step towards the understanding of the molecular mechanisms involved in unisexual flower development in C. sativa, also suggesting that the ABCDE model for flower organ identity may be molecularly conserved in the predominantly monoecious Fagaceae family.

Sex biased expression of hormone related genes at early stage of sex differentiation in papaya flowers

Sex types of papaya are controlled by a pair of nascent sex chromosomes, but molecular genetic mechanisms of sex determination and sex differentiation in papaya are still unclear. We performed comparative analysis of transcriptomic profiles of male and female floral buds at the early development stage before the initiation of reproductive organ primordia at which there is no morphological difference between male and female flowers. A total of 1734 differentially expressed genes (DEGs) were identified, of which 923 showed female-biased expression and 811 showed male-biased expression. Functional annotation revealed that genes related to plant hormone biosynthesis and signaling pathways, especially in abscisic acid and auxin pathways, were overrepresented in the DEGs. Transcription factor binding motifs, such as MYB2, GAMYB, and AP2/EREBP, were enriched in the promoters of the hormone-related DEGs, and transcription factors with those motifs also exhibited differential expression between sex types. Among these DEGs, we also identified 11 genes in the non-recombining region of the papaya sex chromosomes and 9 genes involved in stamen and carpel development. Our results suggested that sex differentiation in papaya may be regulated by multiple layers of regulation and coordination and involved transcriptional, epigenetic, and phytohormone regulation. Hormones, especially ABA and auxin, transcription factors, and genes in the non-recombination region of the sex chromosome could be involved in this process. Our findings may facilitate the elucidation of signal transduction and gene interaction in sex differentiation of unisexual flowers in papaya.

Transcriptome profiling of differentially expressed genes of male and female inflorescences in spinach (Spinacia oleracea L.)

Spinach (Spinacia oleracea L.) is commonly considered a dioecious plant with heterogametic (XY) and homogametic (XX) sex chromosomes. The characteristic is also utilized for the production of spinach hybrid seeds. However, the molecular mechanisms of sex determination in spinach are still unclear because of a lack of genomic and transcriptomic information. In this study, RNA sequencing (RNA-seq) was performed in male and female inflorescences to provide insight into the molecular basis of sex determination in spinach. Comparative transcriptome analyses showed that 2278 differentially expressed genes (DEGs) were identified between male and female inflorescences. A high correlation between the RNA-Seq and qRT-PCR validation for DEGs was observed. Among these, 182 DEGs were annotated to transcription factors including the MYB family protein, bHLH family, and MADS family, suggesting these factors might play a vital role in sex determination. Moreover, 26 DEGs related to flower development, including nine ABCE class genes, were detected. Expression analyses of hormone pathways showed that brassinosteroids may be key hormones related to sex determination in spinach. Overall, this study provides a large amount of DEGs related to sexual expression and lays a foundation for unraveling the regulatory mechanism of sex determination in spinach.

Comparative transcriptome analysis of male and female flowers in Spinacia oleracea L

Background

Dioecious spinach (Spinacia oleracea L.), a commercial and nutritional vegetable crop, serves as a model for studying the mechanisms of sex determination and differentiation in plants. However, this mechanism is still unclear. Herein, based on PacBio Iso-seq and Illumina RNA-seq data, comparative transcriptome analysis of male and female flowers were performed to explore the sex differentiation mechanism in spinach.

Results

Compared with published genome of spinach, 10,800 transcripts were newly annotated; alternative splicing, alternative polyadenylation and lncRNA were analyzed for the first time, increasing the diversity of spinach transcriptome. A total of 2965 differentially expressed genes were identified between female and male flowers at three early development stages. The differential expression of RNA splicing-related genes, polyadenylation-related genes and lncRNAs suggested the involvement of alternative splicing, alternative polyadenylation and lncRNA in sex differentiation. Moreover, 1946 male-biased genes and 961 female-biased genes were found and several candidate genes related to gender development were identified, providing new clues to reveal the mechanism of sex differentiation. In addition, weighted gene co-expression network analysis showed that auxin and gibberellin were the common crucial factors in regulating female or male flower development; however, the closely co-expressed genes of these two factors were different between male and female flower, which may result in spinach sex differentiation.

Conclusions

In this study, 10,800 transcripts were newly annotated, and the alternative splicing, alternative polyadenylation and long-noncoding RNA were comprehensively analyzed for the first time in spinach, providing valuable information for functional genome study. Moreover, candidate genes related to gender development were identified, shedding new insight on studying the mechanism of sex determination and differentiation in plant.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07277-4.

Application of cucumber mosaic virus to efficient induction and long-term maintenance of virus-induced gene silencing in spinach

Virus-induced gene silencing (VIGS) is a useful tool for functional genomics in plants. In this study, we tried to apply cucumber mosaic virus (CMV) to efficient induction of VIGS in spinach. Although VIGS for spinach had been previously developed based on two viruses (beet curly top virus and tobacco rattle virus), they still have some problems with systemic movement and long-term maintenance of VIGS in spinach. Although ordinary CMV strains infect spinach inducing distinct mosaic symptoms, using a CMV pseudorecombinant, we can modify the viral pathogenicity to attenuate viral symptoms that may mask the silencing phenotype. We here successfully demonstrated the viral ability to silence the phytoene desaturase (PDS) and the dihydroflavonol 4-reductase (DFR) genes in spinach. Because CMV could quickly induce VIGS even at 7-10 days postinoculation and the virus did not disappear even at the flowering stage, this CMV-based VIGS system would contribute to functional genomics in spinach and especially to the elucidation of molecular mechanisms for some properties unique to spinach such as plasticity of sex expression; the CMV-induced VIGS can last until the flowering stage after the virus was inoculated onto the seedling.

Potentials and Challenges of Genomics for Breeding Cannabis Cultivars

Cannabis (Cannabis sativa L.) is an influential yet controversial agricultural plant with a very long and prominent history of recreational, medicinal, and industrial usages. Given the importance of this species, we deepened some of the main challenges-along with potential solutions-behind the breeding of new cannabis cultivars. One of the main issues that should be fixed before starting new breeding programs is the uncertain taxonomic classification of the two main taxa (e.g., indica and sativa) of the Cannabis genus. We tried therefore to examine this topic from a molecular perspective through the use of DNA barcoding. Our findings seem to support a unique species system (C. sativa) based on two subspecies: C. sativa subsp. sativa and C. sativa subsp. indica. The second key issue in a breeding program is related to the dioecy behavior of this species and to the comprehension of those molecular mechanisms underlying flower development, the main cannabis product. Given the role of MADS box genes in flower identity, we analyzed and reorganized all the genomic and transcriptomic data available for homeotic genes, trying to decipher the applicability of the ABCDE model in Cannabis. Finally, reviewing the limits of the conventional breeding methods traditionally applied for developing new varieties, we proposed a new breeding scheme for the constitution of F₁ hybrids, without ignoring the indisputable contribution offered by genomics. In this sense, in parallel, we resumed the main advances in the genomic field of this species and, ascertained the lack of a robust set of SNP markers, provided a discriminant and polymorphic panel of SSR markers as a valuable tool for future marker assisted breeding programs.

A Model of Hormonal Regulation of Stamen Abortion during Pre-Meiosis of Litsea cubeba

Litsea cubeba (Lour.) Pers., a popular essential oil plant, is a dioecious species with degenerative sexual organs in both male and female individuals. Yet, the mechanism of degenerative organs development in male and female flowers is poorly understood. Here, we analyzed the morphological characters of degenerative organ development by morphological and histological observations, and determined the critical stage of abortion that occurs at pre-meiosis in male and female flowers. We also conducted RNA sequencing (RNA-seq) to understand the genetic basis of stamen abortion in female flowers. The differentially expressed genes (DEGs) were identified during the staminode development in female flowers; functional enrichment analysis revealed some important biological pathways involved the regulation of stamen abortion, including plant hormone signal transduction, phenylpropanoid biosynthesis, flavonoid biosynthesis and monoterpenoid biosynthesis. Furthermore, 15 DEGs involved in the hormone pathways were found to regulate stamen development. By HPLC-MS/MS analysis, there were a salicylic acid (SA) content peak and the gibberellin (GA) content lowest point in the abortion processes in female flowers, suggesting a vital function of hormonal processes. Co-expression network analysis further identified several hub genes that potentially played significant roles in the stamen abortion of L. cubeba. Taken together, we proposed a model involved in plant hormones pathways underlying stamen abortion during pre-meiosis in female flowers of L. cubeba.

Epigenetic aspects of floral homeotic genes in relation to sexual dimorphism in the dioecious plant Mercurialis annua

In plants, dioecy characterizes species that carry male and female flowers on separate plants and it occurs in about 6% of angiosperms; however, the molecular mechanisms that underlie dioecy are essentially unknown. The ability for sex-reversal by hormone application raises the hypothesis that the genes required for the expression of both sexes are potentially functional but are regulated by epigenetic means. In this study, proteomic analysis of nuclear proteins isolated from flower buds of females, males, and feminized males of the dioecious plant Mercurialis annua revealed differential expression of nuclear proteins that are implicated in chromatin structure and function, including floral homeotic proteins. Focusing on floral genes, we found that class B genes were mainly expressed in male flowers, while class D genes, as well as SUPERMAN-like genes, were mainly expressed in female flowers. Cytokinin-induced feminization of male plants was associated with down-regulation of male-specific genes concomitantly with up-regulation of female-specific genes. No correlation was found between the expression of class B and D genes and the changes in DNA methylation or chromatin conformation of these genes. Thus, we could not confirm DNA methylation or chromatin conformation of floral genes to be the major determinant regulating sexual dimorphisms. Instead, determination of sex in M. annua might be controlled upstream of floral genes by one or more sex-specific factors that affect hormonal homeostasis. A comprehensive model is proposed for sex-determination in M. annua.

Gene networks orchestrated by MeGI: a single-factor mechanism underlying sex determination in persimmon

Separating male and female sex organs is one of the main strategies used to maintain genetic diversity within a species. However, the genetic determinants and their regulatory mechanisms have been identified in only a few species. In dioecious persimmons, the homeodomain transcription factor, MeGI, which is the target of a Y chromosome-encoded small-RNA, OGI, can determine floral sexuality. The basic features of this system are conserved in the monoecious hexaploid Oriental persimmon, in which an additional epigenetic regulation of MeGI determines floral sexuality. The downstream regulatory pathways of MeGI remain uncharacterized. In this study, we examined transcriptomic data for male and female flowers from monoecious persimmon cultivars to unveil the gene networks orchestrated by MeGI. A network visualization and cistrome assessment suggested that class-1 KNOTTED-like homeobox (KNOX)/ovate family protein (OFP)/growth regulating factors (GRFs) and short vegetative phase (SVP) genes mediate the differences in gynoecium and androecium development between male and female flowers, respectively. The expression of these genes is directly controlled by MeGI. The gene networks also suggested that some cytokinin, auxin, and gibberellin signaling genes function cooperatively in the KNOX/OFP/GRF pathway during gynoecium differentiation. Meanwhile, SVP may repress PI expression in developing androecia. Overall, our results suggest that MeGI evolved the ability to promote gynoecium development and suppress androecium development by regulating KNOX/OFP/GRF and SVP expression levels, respectively. These insights may help to clarify the molecular mechanism underlying the production of unisexual flowers, while also elucidating the physiological background enabling a single-factor system to establish dioecy in plants.

RNA interference suppression of AGAMOUS and SEEDSTICK alters floral organ identity and impairs floral organ determinacy, ovule differentiation, and seed-hair development in Populus

The role of the floral homeotic gene AGAMOUS (AG) and its close homologues in development of anemophilous, unisexual catkins has not previously been studied. We transformed two RNA interference (RNAi) constructs, PTG and its matrix-attachment-region flanked version MPG, into the early-flowering female poplar clone 6K10 (Populus alba) to suppress the expression of its two duplicate AG orthologues. By early 2018, six out of 22 flowering PTG events and 11 out of 12 flowering MPG events showed modified floral phenotypes in a field trial in Oregon, USA. Flowers in catkins from modified events had 'carpel-inside-carpel' phenotypes. Complete disruption of seed production was observed in seven events, and sterile anther-like organs in 10 events. Events with strong co-suppression of both the two AG and two SEEDSTICK (STK) paralogues lacked both seeds and associated seed hairs. Alterations in all of the modified floral phenotypes were stable over 4 yr of study. Trees from floral-modified events did not differ significantly (P < 0.05) from nonmodified transgenic or nontransgenic controls in biomass growth or leaf morphology. AG and STK genes show strong conservation of gene function during poplar catkin development and are promising targets for genetic containment of exotic or genetically engineered trees.

Vitis Flower Sex Specification Acts Downstream and Independently of the ABCDE Model Genes

The most discriminating characteristic between the cultivated Vitis vinifera subsp. vinifera and the wild-form Vitis vinifera subsp. sylvestris is their sexual system. Flowers of cultivars are mainly hermaphroditic, whereas wild plants have female and male individuals whose flowers follow a hermaphroditic pattern during early stages of development and later develop non-functional reproductive organs. In angiosperms, the basic developmental system for floral organ identity is explained by the ABCDE model. This model postulates that regulatory gene functions work in a combinatorial way to confer organ identity in each whorl. In wild Vitis nothing is known about the function and expression profile of these genes. Here we show an overall view of the temporal and spatial expression pattern of the ABCDE genes as well as the pattern of VviSUPERMAN that establishes a boundary between the stamen and the carpel whorls, in the male, female and complete flower types. The results show a similar pattern in Vitis species suggesting that the pathway leading to unisexuality acts independently and/or downstream of B- and C- function genes.

Gender-specific expression of GIBBERELLIC ACID INSENSITIVE is critical for unisexual organ initiation in dioecious Spinacia oleracea

While unisexual flowers have evolved repeatedly throughout angiosperm families, the actual identification of sex-determining genes has been elusive, and their regulation within populations remains largely undefined. Here, we tested the mechanism of the feminization pathway in cultivated spinach (Spinacia oleracea), and investigated how this pathway may regulate alternative sexual development. We tested the effect of gibberellic acid (GA) on sex determination through exogenous applications of GA and inhibitors of GA synthesis and proteasome activity. GA concentrations in multiple tissues were estimated by enzyme-linked immunosorbent assay analysis. Gene function was investigated and pathway analysis was performed through virus-induced gene silencing. Relative gene expression levels were estimated by quantitative reverse transcription-polymerase chain reaction. Inhibition of GA production and proteasome activity feminized male flowers. However, there was no difference in GA content in tissues between males and females. We characterized a single DELLA family transcription factor gene (GIBBERELLIC ACID INSENSITIVE (SpGAI)) and observed inflorescence expression in females two-fold higher than in males. Reduction of SpGAI expression in females to male levels phenocopied exogenous GA application with respect to flower development. These results implicate SpGAI as the feminizing factor in spinach, and suggest that the feminizing pathway is epistatic to the masculinizing pathway. We present a unified model for alternative sexual development and discuss the implications for established theory.

Role of floral organ identity genes in the development of unisexual flowers of Quercus suber L

Monoecious species provide an excellent system to study the specific determinants that underlie male and female flower development. Quercus suber is a monoecious species with unisexual flowers at inception. Despite the overall importance of this and other tree species with a similar reproductive habit, little is known regarding the mechanisms involved in the development of their male and female flowers. Here, we have characterised members of the ABCDE MADS-box gene family of Q. suber. The temporal expression of these genes was found to be sex-biased. The B-class genes, in particular, are predominantly, or exclusively (in the case of QsPISTILLATA), expressed in the male flowers. Functional analysis in Arabidopsis suggests that the B-class genes have their function conserved. The identification of sex-biased gene expression plus the identification of unusual protein-protein interactions suggest that the floral organ identity of Q. suber may be under control of specific changes in the dynamics of the ABCDE model. This study constitutes a major step towards the characterisation of the mechanisms involved in reproductive organ identity in a monoecious tree with a potential contribution towards the knowledge of conserved developmental mechanisms in other species with a similar sex habit.

Rapid sex identification method of spinach (Spinacia oleracea L.) in the vegetative stage using loop-mediated isothermal amplification

A LAMP-mediated, simple and rapid method for sex identification in spinach was developed. Nutrient compositional analysis showed a higher iron content in male than female plants. Spinach (Spinacia oleracea L.) is a dioecious plant with its sex determined by the XY system. Male and female floral organs differ morphologically, but plants do not differ in the vegetative stage before flowering. PCR with Y chromosome markers has been used to determine the sex of dioecious plants before flowering. In this study, we developed a genotype-specific loop-mediated isothermal amplification (LAMP) for sex identification of individual vegetative-stage spinach plants, using primers designed for the genomic region flanked by male-specific markers. LAMP could specifically detect spinach males. The method was further modified to omit DNA purification and use just an aliquot of crude leaf extract homogenized in water. We compared the nutrient composition of males and females, finding higher amounts of iron in the males. Our method could therefore be used for rapidly discriminating male plants in the field, which is useful for efficient hybrid breeding.

The Quest for Molecular Regulation Underlying Unisexual Flower Development

The understanding of the molecular mechanisms responsible for the making of a unisexual flower has been a long-standing quest in plant biology. Plants with male and female flowers can be divided mainly into two categories: dioecious and monoecious, and both sexual systems co-exist in nature in ca of 10% of the angiosperms. The establishment of male and female traits has been extensively described in a hermaphroditic flower and requires the interplay of networks, directly and indirectly related to the floral organ identity genes including hormonal regulators, transcription factors, microRNAs, and chromatin-modifying proteins. Recent transcriptomic studies have been uncovering the molecular processes underlying the establishment of unisexual flowers and there are many parallelisms between monoecious, dioecious, and hermaphroditic individuals. Here, we review the paper entitled "Comparative transcriptomic analysis of male and female flowers of monoecious Quercus suber" published in 2014 in the Frontiers of Plant Science (volume 5 |Article 599) and discussed it in the context of recent studies with other dioecious and monoecious plants that utilized high-throughput platforms to obtain transcriptomic profiles of male and female unisexual flowers. In some unisexual flowers, the developmental programs that control organ initiation fail and male or female organs do not form, whereas in other species, organ initiation and development occur but they abort or arrest during different species-specific stages of differentiation. Therefore, a direct comparison of the pathways responsible for the establishment of unisexual flowers in different species are likely to reveal conserved modules of gene regulatory hubs involved in stamen or carpel development, as well as differences that reflect the different stages of development in which male and/or female organ arrest or loss-of-function occurs.

Rapid cloning and bioinformatic analysis of spinach Y chromosome-specific EST sequences

The genome of spinach single chromosome complement is about 1000 Mbp, which is the model material to study the molecular mechanisms of plant sex differentiation. The cytological study showed that the biggest spinach chromosome (chromosome 1) was taken as spinach sex chromosome. It had three alleles of sex-related X,X(m) and Y. Many researchers have been trying to clone the sex-determining genes and investigated the molecular mechanism of spinach sex differentiation. However,there are no successful cloned reports about these genes. A new technology combining chromosome microdissection with hybridization-specific amplification (HSA) was adopted. The spinach Y chromosome degenerate oligonucleotide primed-PCR (DOP-PCR) products were hybridized with cDNA of the male spinach flowers in florescence. The female spinach genome was taken as blocker and cDNA library specifically expressed in Y chromosome was constructed. Moreover, expressed sequence tag (EST) sequences in cDNA library were cloned, sequenced and bioinformatics was analysed. There were 63 valid EST sequences obtained in this study. The fragment size was between 53 and 486 bp. BLASTn homologous alignment indicated that 12 EST sequences had homologous sequences of nucleic acids, the rest were new sequences. BLASTx homologous alignment indicated that 16 EST sequences had homologous protein-encoding nucleic acid sequence. The spinach Y chromosome-specific EST sequences laid the foundation for cloning the functional genes, specifically expressed in spinach Y chromosome. Meanwhile, the establishment of the technology system in the research provided a reference for rapid cloning of other biological sex chromosome-specific EST sequences.

The sex-linked region in Populus tremuloides Turesson 141 corresponds to a pericentromeric region of about two million base pairs on P. trichocarpa chromosome 19

In the dioecious genus Populus, sex determination has been located to chromosome 19. However, despite a high degree of genome collinearity, various Populus species seem to differ with regard to the location of the sex-determining region on the respective chromosome and the apparent heterogametic sex. In this study, the boundaries of the recombination-suppressed, sex-linked region of the male P. tremuloides clone Turesson 141 were localised by genetic mapping using new SNP and InDel markers. The respective region seems to be located in a pericentromeric position. The corresponding P. trichocarpa genome region spans about two million bp and comprises 65 gene loci, which were bioinformatically evaluated for their potential as candidate genes for sex determination. Three putative transcription factor genes and four genes that are potentially involved in flower development processes, e.g. meristem transition from the vegetative to the reproductive phase, were identified. Populus tremuloides sequence data of the sex-linked region is required for a final search for candidate genes.

Microdissection and painting of the Y chromosome in spinach (Spinacia oleracea)

Spinach has long been used as a model for genetic and physiological studies of sex determination and expression. Although trisomic analysis from a cross between diploid and triploid plants identified the XY chromosome as the largest chromosome, no direct evidence has been provided to support this at the molecular level. In this study, the largest chromosomes of spinach from mitotic metaphase spreads were microdissected using glass needles. Degenerate oligonucleotide primed polymerase chain reaction was used to amplify the dissected chromosomes. The amplified products from the Y chromosome were identified using the male-specific marker T11A. For the first time, the largest spinach chromosome was confirmed to be a sex chromosome at the molecular level. PCR products from the isolated chromosomes were used in an in situ probe mixture for painting the Y chromosome. The fluorescence signals were mainly distributed on all chromosomes and four pair of weaker punctate fluorescence signal sites were observed on the terminal region of two pair of autosomes. These findings provide a foundation for the study of sex chromosome evolution in spinach.

Functional recapitulation of transitions in sexual systems by homeosis during the evolution of dioecy in Thalictrum

Sexual systems are highly variable in flowering plants and an important contributor to floral diversity. The ranunculid genus Thalictrum is especially well-suited to study evolutionary transitions in sexual systems. Homeotic transformation of sexual organs (stamens and carpels) is a plausible mechanism for the transition from hermaphroditic to unisexual flowers in this lineage because flowers of dioecious species develop unisexually from inception. The single-copy gene PISTILLATA (PI) constitutes a likely candidate for rapid switches between stamen and carpel identity. Here, we first characterized the expression pattern of all B class genes in the dioecious species T. dioicum. As expected, all B class orthologs are expressed in stamens from the earliest stages. Certain AP3 lineages were also expressed late in sepal development. We then tested whether orthologs of PI could potentially control sexual system transitions in Thalictrum, by knocking-down their expression in T. dioicum and the hermaphroditic species T. thalictroides. In T. dioicum, we found that ThdPI-1/2 silencing caused stamen primordia to develop into carpels, resulting in male to female flower conversions. In T. thalictroides, we found that ThtPI silencing caused stamen primordia to develop into supernumerary carpels, resulting in hermaphroditic to female flower conversions. These phenotypes illustrate the ability for homeotic mutations to bring about sudden and potentially adaptive changes by altering the function of a single gene. We propose a two-step evolutionary model where transitions from hermaphroditic to unisexual plants in Thalictrum result from two independent mutations at a B class gene locus. Our PI knockdown experiments in T. thalictroides recapitulate the second step in this model: the evolution of female plants as a result of a loss-of-function mutation in a B class gene.

Virus-induced gene silencing (VIGS) in plants: an overview of target species and the virus-derived vector systems

The analysis of gene functions in non-model plant species is often hampered by the fact that stable genetic transformation to downregulate gene expression is laborious and time-consuming, or, for some species, even not achievable. Virus-induced gene silencing (VIGS) can serve as an alternative to mutant collections or stable transgenic plants to allow the characterization of gene functions in a wide range of angiosperm species, albeit in a transient way. VIGS vector systems have been developed from both RNA and DNA plant viral sources to specifically silence target genes in plants. VIGS is nowadays widely used in plant genetics for gene knockdown due to its ease of use and the short time required to generating phenotypes. Here, we summarize successfully targeted eudicot and monocot plant species along with their specific VIGS vector systems which are already available for researchers.

Developmental and evolutionary diversity of plant MADS-domain factors: insights from recent studies

Members of the MADS-box transcription factor family play essential roles in almost every developmental process in plants. Many MADS-box genes have conserved functions across the flowering plants, but some have acquired novel functions in specific species during evolution. The analyses of MADS-domain protein interactions and target genes have provided new insights into their molecular functions. Here, we review recent findings on MADS-box gene functions in Arabidopsis and discuss the evolutionary history and functional diversification of this gene family in plants. We also discuss possible mechanisms of action of MADS-domain proteins based on their interactions with chromatin-associated factors and other transcriptional regulators.

Multiple developmental processes underlie sex differentiation in angiosperms

The production of unisexual flowers has evolved numerous times in dioecious and monoecious plant taxa. Based on repeated evolutionary origins, a great variety of developmental and genetic mechanisms underlying unisexual flower development is predicted. Here, we comprehensively review the modes of development of unisexual flowers, test potential correlations with sexual system, and end with a synthesis of the genetics and hormonal regulation of plant sex determination. We find that the stage of organ abortion in male and female flowers is temporally correlated within species and also confirm that the arrest of development does not tend to occur preferentially at a particular stage, or via a common process.

Mapping of the genes for dioecism and monoecism in Spinacia oleracea L.: evidence that both genes are closely linked

Spinach is basically a dioecious species, with occasional monoecious plants in some populations. Sexual dimorphism in dioecious spinach plants is controlled by an allelic pair termed X and Y located on the short arm of the longest chromosome (x = 6). Ten AFLP markers, closely linked to the X/Y locus, were identified using bulked segregant analysis, four of which were revealed to co-segregate with Y in the present mapping population. We mapped the AFLP markers and two known male-specific DNAs to a 13.4 cM region encompassing the locus. These markers will be the basis for positional cloning of the sex-determination gene. We also showed that a single, incompletely dominant gene is responsible for the highly staminate monoecious character. The gene was found to be located at a distance of 4.3 cM from microsatellite marker SO4, which mapped 1.6 cM from the X/Y locus. This indicates that the monoecious gene seems not to be allelic to but closely linked to the X/Y gene pair. SO4 will enable breeders to efficiently select highly male monoecious plants for preferential use as the pollen parent for hybrid seed production.