A cucurbit androecy gene reveals how unisexual flowers develop and dioecy emerges

Genome assembly of the basket willow, Salix viminalis, reveals earliest stages of sex chromosome expansion

BACKGROUND

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

RESULTS

Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

CONCLUSIONS

Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

Genome assembly of the basket willow, Salix viminalis, reveals earliest stages of sex chromosome expansion

BACKGROUND

Sex chromosomes have evolved independently multiple times in eukaryotes and are therefore considered a prime example of convergent genome evolution. Sex chromosomes are known to emerge after recombination is halted between a homologous pair of chromosomes, and this leads to a range of non-adaptive modifications causing gradual degeneration and gene loss on the sex-limited chromosome. However, the proximal causes of recombination suppression and the pace at which degeneration subsequently occurs remain unclear.

RESULTS

Here, we use long- and short-read single-molecule sequencing approaches to assemble and annotate a draft genome of the basket willow, Salix viminalis, a species with a female heterogametic system at the earliest stages of sex chromosome emergence. Our single-molecule approach allowed us to phase the emerging Z and W haplotypes in a female, and we detected very low levels of Z/W single-nucleotide divergence in the non-recombining region. Linked-read sequencing of the same female and an additional male (ZZ) revealed the presence of two evolutionary strata supported by both divergence between the Z and W haplotypes and by haplotype phylogenetic trees. Gene order is still largely conserved between the Z and W homologs, although the W-linked region contains genes involved in cytokinin signaling regulation that are not syntenic with the Z homolog. Furthermore, we find no support across multiple lines of evidence for inversions, which have long been assumed to halt recombination between the sex chromosomes.

CONCLUSIONS

Our data suggest that selection against recombination is a more gradual process at the earliest stages of sex chromosome formation than would be expected from an inversion and may result instead from the accumulation of transposable elements. Our results present a cohesive understanding of the earliest genomic consequences of recombination suppression as well as valuable insights into the initial stages of sex chromosome formation and regulation of sex differentiation.

A single gene underlies the dynamic evolution of poplar sex determination

Although hundreds of plant lineages have independently evolved dioecy (that is, separation of the sexes), the underlying genetic basis remains largely elusive¹. Here we show that diverse poplar species carry partial duplicates of the ARABIDOPSIS RESPONSE REGULATOR 17 (ARR17) orthologue in the male-specific region of the Y chromosome. These duplicates give rise to small RNAs apparently causing male-specific DNA methylation and silencing of the ARR17 gene. CRISPR-Cas9-induced mutations demonstrate that ARR17 functions as a sex switch, triggering female development when on and male development when off. Despite repeated turnover events, including a transition from the XY system to a ZW system, the sex-specific regulation of ARR17 is conserved across the poplar genus and probably beyond. Our data reveal how a single-gene-based mechanism of dioecy can enable highly dynamic sex-linked regions and contribute to maintaining recombination and integrity of sex chromosomes.

Integrative genome-wide analysis reveals the role of WIP proteins in inhibition of growth and development

In cucurbits, CmWIP1 is a master gene controlling sex determination. To bring new insight in the function of CmWIP1, we investigated two Arabidopsis WIP transcription factors, AtWIP1/TT1 and AtWIP2/NTT. Using an inducible system we showed that WIPs are powerful inhibitor of growth and inducer of cell death. Using ChIP-seq and RNA-seq we revealed that most of the up-regulated genes bound by WIPs display a W-box motif, associated with stress signaling. In contrast, the down-regulated genes contain a GAGA motif, a known target of polycomb repressive complex. To validate the role of WIP proteins in inhibition of growth, we expressed AtWIP1/TT1 in carpel primordia and obtained male flowers, mimicking CmWIP1 function in melon. Using other promoters, we further demonstrated that WIPs can trigger growth arrest of both vegetative and reproductive organs. Our data supports an evolutionary conserved role of WIPs in recruiting gene networks controlling growth and adaptation to stress.

Maria V Gomez Roldan, Farhaj Izhaq et al. study the role of WIP transcription factors in Arabidopsis thaliana. Using an inducible system, they generate ChIP-seq and RNA-seq data that reveal how WIP proteins repress development and growth. Through ectopic expression of WIP in Arabidopsis carpels, they obtain male flower, mimicking sex determination in melon.

A Bridge to the World

Zhi-Hong Xu is a plant physiologist who studied botany at Peking University (1959–1965). He joined the Shanghai Institute of Plant Physiology (SIPP), Chinese Academy of Sciences (CAS), as a graduate student in 1965. He recalls what has happened for the institute, during the Cultural Revolution, and he witnessed the spring of science eventually coming to China. Xu was a visiting scholar at the John Innes Institute and in the Department of Botany at Nottingham University in the United Kingdom (1979–1981). He became deputy director of SIPP in 1983 and director in 1991; he also chaired the State Key Laboratory of Plant Molecular Genetics SIPP (1988–1996). He worked as a visiting scientist in the Institute of Molecular and Cell Biology, National University of Singapore, for three months each year (1989–1992). He served as vice president of CAS (1992–2002) and as president of Peking University (1999–2008). Over these periods he was heavily involved in the design and implementation of major scientific projects in life sciences and agriculture in China. He is an academician of CAS and member of the Academy of Sciences for the Developing World. His scientific contributions mainly cover plant tissue culture, hormone mechanism in development, as well as plant developmental response to environment. Xu, as a scientist and leader who has made an impact in the community, called up a lot of excellent young scientists returning to China. His efforts have promoted the fast development of China's plant and agricultural sciences.

Pathways to sex determination in plants: how many roads lead to Rome?

The presence of thousands of independent origins of dioecy in angiosperms provides a unique opportunity to address the parallel evolution of the molecular pathways underlying unisexual flowers. Recent progress towards identifying sex determination genes has identified hormone response pathways, mainly associated with cytokinin and ethylene response pathways, as having been recruited multiple times independently to control unisexuality. Moreover, transcriptomics has begun to identify commonalities among intermediate sections of signal transduction pathways. These recent advances set the stage for development of a comparative evolutionary development research program to identify the shared and unique aspects of the genetic pathways of unisexual flower development in angiosperms.

Gynoecy instability in cucumber (Cucumis sativus L.) is due to unequal crossover at the copy number variation-dependent Femaleness (F) locus

Cucumber, Cucumis sativus is an important vegetable crop, and gynoecy has played a critical role in yield increase of hybrid cucumber production. Cucumber has a unique genetic system for gynoecious sex expression, which is determined by the copy number variation (CNV)-based, dominant, and dosage-dependent femaleness (F) locus. However, this gynoecy expression system seems unstable since monecious plants could often be found in F-dependent gynoecious cucumber inbreds. We hypothesized that gynoecy instability (gynoecy loss) may be due to unequal crossing over (UCO) during meiosis among repeat units of the CNV. In this study, using high throughput genome resequencing, fiber-FISH and genomic qPCR analyses, we first confirmed and refined the structure of the F locus, which was a CNV of a 30.2-kb tandem repeat. Gynoecious plants contained three genes: CsACS1, CsACS1G, and CsMYB, of which CsACS1G is a duplication of CsACS1 but with a recombinant distal promoter that may contribute to gynoecy sex expression. In two large populations from self-pollinated gynoecious inbred lines, 'gynoecy loss' mutants were identified with similar mutation rates (~0.12%). We show that these monecious mutants have lost CsACS1G. In addition, we identified gynoecious lines in natural populations that carry two copies of CSACS1G. We proposed a model to explain gynoecy instability in F-dependent cucumbers, which is caused by UCO among CSACS1/G units during meiosis. The findings present a convincing case that the phenotypic variation of an economically important trait is associated with the dynamic changes of copy numbers at the F locus. This work also has important implications in cucumber breeding.

De novo genome assembly and Hi-C analysis reveal an association between chromatin architecture alterations and sex differentiation in the woody plant Jatropha curcas

BACKGROUND

Chromatin architecture is an essential factor regulating gene transcription in different cell types and developmental phases. However, studies on chromatin architecture in perennial woody plants and on the function of chromatin organization in sex determination have not been reported.

RESULTS

Here, we produced a chromosome-scale de novo genome assembly of the woody plant Jatropha curcas with a total length of 379.5 Mb and a scaffold N50 of 30.7 Mb using Pacific Biosciences long reads combined with genome-wide chromosome conformation capture (Hi-C) technology. Based on this high-quality reference genome, we detected chromatin architecture differences between monoecious and gynoecious inflorescence buds of Jatropha. Differentially expressed genes were significantly enriched in the changed A/B compartments and topologically associated domain regions and occurred preferentially in differential contact regions between monoecious and gynoecious inflorescence buds. Twelve differentially expressed genes related to flower development or hormone synthesis displayed significantly different genomic interaction patterns in monoecious and gynoecious inflorescence buds. These results demonstrate that chromatin organization participates in the regulation of gene transcription during the process of sex differentiation in Jatropha.

CONCLUSIONS

We have revealed the features of chromatin architecture in perennial woody plants and investigated the possible function of chromatin organization in Jatropha sex differentiation. These findings will facilitate understanding of the regulatory mechanisms of sex determination in higher plants.

Love is in the air: ethylene and sex determination in Cucurbita pepo

This article comments on:

García A, Aguado E, Martínez C, Loska D, Beltrán S, Valenzuela JL, Garrido D, Jamilena M. 2019. The ethylene receptors CpETR1A and CpETR2B cooperate in the control of sex determination in Cucurbita pepo. Journal of Experimental Botany 70, 154–167.

The ethylene receptors CpETR1A and CpETR2B cooperate in the control of sex determination in Cucurbita pepo

High-throughput screening of an ethyl methanesulfonate-generated mutant collection of Cucurbita pepo using the ethylene triple-response test resulted in the identification of two semi-dominant ethylene-insensitive mutants: etr1a and etr2b. Both mutations altered sex determination mechanisms, promoting conversion of female into bisexual or hermaphrodite flowers, and monoecy into andromonoecy, thereby delaying the transition to female flowering and reducing the number of pistillate flowers per plant. The mutations also altered the growth rate and maturity of petals and carpels in pistillate flowers, lengthening the time required for flowers to reach anthesis, as well as stimulating the growth rate of ovaries and the parthenocarpic development of fruits. Whole-genome sequencing allowed identification of the causal mutation of the phenotypes as two missense mutations in the coding region of CpETR1A and CpETR2B, each one corresponding to one of the duplicates of ethylene receptor genes highly homologous to Arabidopsis ETR1 and ETR2. The phenotypes of homozygous and heterozygous single- and double-mutant plants indicated that the two ethylene receptors cooperate in the control of the ethylene response. The level of ethylene insensitivity, which was determined by the strength of each mutant allele and the dose of wild-type and mutant etr1a and etr2b alleles, correlated with the degree of phenotypic changes in the mutants.

Molecularly tagged genes and quantitative trait loci in cucumber with recommendations for QTL nomenclature

Cucumber, Cucumis sativus L. (2n = 2x = 14), is an important vegetable crop worldwide. It was the first specialty crop with a publicly available draft genome. Its relatively small, diploid genome, short life cycle, and self-compatible mating system offers advantages for genetic studies. In recent years, significant progress has been made in molecular mapping, and identification of genes and QTL responsible for key phenotypic traits, but a systematic review of the work is lacking. Here, we conducted an extensive literature review on mutants, genes and QTL that have been molecularly mapped or characterized in cucumber. We documented 81 simply inherited trait genes or major-effect QTL that have been cloned or fine mapped. For each gene, detailed information was compiled including chromosome locations, allelic variants and associated polymorphisms, predicted functions, and diagnostic markers that could be used for marker-assisted selection in cucumber breeding. We also documented 322 QTL for 42 quantitative traits, including 109 for disease resistances against seven pathogens. By alignment of these QTL on the latest version of cucumber draft genomes, consensus QTL across multiple studies were inferred, which provided insights into heritable correlations among different traits. Through collaborative efforts among public and private cucumber researchers, we identified 130 quantitative traits and developed a set of recommendations for QTL nomenclature in cucumber. This is the first attempt to systematically summarize, analyze and inventory cucumber mutants, cloned or mapped genes and QTL, which should be a useful resource for the cucurbit research community.

A unique chromosome translocation disrupting ClWIP1 leads to gynoecy in watermelon

To understand sex determination in watermelon (Citrullus lanatus), a spontaneous gynoecious watermelon mutant, XHBGM, was selected from the monoecious wild type XHB. Using map-based cloning, resequencing and fluorescence in situ hybridization analysis, a unique chromosome translocation between chromosome 2 and chromosome 3 was found in XHBGM. Based on the breakpoint location in chromosome 2, a putative C2H2 zinc finger transcription factor gene, ClWIP1 (gene ID Cla008537), an orthologue of the melon gynoecy gene CmWIP1, was disrupted. Using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated system 9 to edit ClWIP1, we obtained gynoecious watermelon lines. Functional studies showed that ClWIP1 is expressed specifically in carpel primordia and is related to the abortion of carpel primordia in early floral development. To identify the cellular and metabolic processes associated with ClWIP1, we compared the shoot apex transcriptomes of two gynoecious mutants and their corresponding wild types. Transcriptome analysis showed that differentially expressed genes related to the ethylene and cytokinin pathways were upregulated in the gynoecious mutants. This study explores the molecular mechanism of sex determination in watermelon and provides a theoretical and technical basis for breeding elite gynoecious watermelon lines.

Default Sex and Single Gene Sex Determination in Dioecious Plants

A well-established hypothesis for the evolution of dioecy involves two genes linked at a sex-determining region (SDR). Recently there has been increased interest in possible single gene sex determination. Work in Populus has finally provided direct experimental evidence for single gene sex determination in plants using CRISPR-Cas9 to knock out a single gene and convert individuals from female to male. In poplar, the feminizing factor popARR17 acts as a "master regulator", analogous to the mammalian masculinizing factor SRY. The production of fully functional males from females by a simple single gene knockout is experimental evidence that an antagonistic male-determining factor does not exist in Populus. Mammals have a "default sex" (female), as do poplar trees (Populus), although the default sex in poplars is male. The occurrence of single gene sex determination with a default sex may be much commoner in plants than hitherto expected, especially when dioecy evolved via monoecy. The master regulator does not even need to be at the SDR (although it may be). In most poplars the feminizing factor popARR17 is not at the SDR, but instead a negative regulator of it. So far there is little information on how high-level regulators are connected to floral phenotype. A model is presented of how sex-determining genes could lead to different floral morphologies via MADS-box floral developmental genes.

Genetic architecture of fruit size and shape variation in cucurbits: a comparative perspective

The Cucurbitaceae family hosts many economically important fruit vegetables (cucurbits) such as cucumber, melon, watermelon, pumpkin/squash, and various gourds. The cucurbits are probably best known for the diverse fruit sizes and shapes, but little is known about their genetic basis and molecular regulation. Here, we reviewed the literature on fruit size (FS), shape (FSI), and fruit weight (FW) QTL identified in cucumber, melon, and watermelon, from which 150 consensus QTL for these traits were inferred. Genome-wide survey of the three cucurbit genomes identified 253 homologs of eight classes of fruit or grain size/weight-related genes cloned in Arabidopsis, tomato, and rice that encode proteins containing the characteristic CNR (cell number regulator), CSR (cell size regulator), CYP78A (cytochrome P450), SUN, OVATE, TRM (TONNEAU1 Recruiting Motif), YABBY, and WOX domains. Alignment of the consensus QTL with candidate gene homologs revealed widespread structure and function conservation of fruit size/shape gene homologs in cucurbits, which was exemplified with the fruit size/shape candidate genes CsSUN25-26-27a and CsTRM5 in cucumber, CmOFP1a in melon, and ClSUN25-26-27a in watermelon. In cucurbits, the andromonoecy (for 1-aminocyclopropane-1-carboxylate synthase) and the carpel number (for CLAVATA3) loci are known to have pleiotropic effects on fruit shape, which may complicate identification of fruit size/shape candidate genes in these regions. The present work illustrates the power of comparative analysis in understanding the genetic architecture of fruit size/shape variation, which may facilitate QTL mapping and cloning for fruit size-related traits in cucurbits. The limitations and perspectives of this approach are also discussed.

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.

Anatomy and RNA-Seq reveal important gene pathways regulating sex differentiation in a functionally Androdioecious tree, Tapiscia sinensis

BACKGROUND

Gametogenesis is a key step in the production of ovules or pollen in higher plants. The sex-determination aspects of gametogenesis have been well characterized in the model plant Arabidopsis. However, little is known about this process in androdioecious plants. Tapiscia sinensis Oliv. is a functionally androdioecious tree, with both male and hermaphroditic individuals. Hermaphroditic flowers (HFs) are female-fertile flowers that can produce functional pollen and set fruits. However, compared with male flowers (MFs), the pollen viability and number of pollen grains per flower are markedly reduced in HFs. MFs are female-sterile flowers that fail to set fruit and that eventually drop.

RESULTS

Compared with HF, a notable cause of MF female sterility in T. sinensis is when the early gynoecium meristem is disrupted. During the early stage of HF development (stage 6), the ring meristem begins to form as a ridge around the center of the flower. At this stage, the internal fourth-whorl organ is stem-like rather than carpelloid in MF. A total of 52,945 unigenes were identified as transcribed in MF and HF. A number of differentially expressed genes (DEGs) and metabolic pathways were detected as involved in the development of the gynoecium, especially the ovule, carpel and style. At the early gynoecium development stage, DEGs were shown to function in the metabolic pathways regulating ethylene biosynthesis and signal transduction (upstream regulator), auxin, cytokinin transport and signalling, and sex determination (or flower meristem identity).

CONCLUSIONS

Pathways for the female sterility model were initially proposed to shed light on the molecular mechanisms of gynoecium development at early stages in T. sinensis.

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.

A High-Quality Melon Genome Assembly Provides Insights into Genetic Basis of Fruit Trait Improvement

Accurate reference genomes have become indispensable tools for characterization of genetic and functional variations. Here we generated a high-quality assembly of the melon Payzawat using a combination of short-read sequencing, single-molecule real-time sequencing, Hi-C, and a high-density genetic map. The final 12 chromosome-level scaffolds cover ∼94.13% of the estimated genome (398.57 Mb). Compared with the published DHL92 genome, our assembly exhibits a 157-fold increase in contig length and remarkable improvements in the assembly of centromeres and telomeres. Six genes within STHQF12.4 on pseudochromosome 12, identified from whole-genome comparison between Payzawat and DHL92, may explain a considerable proportion of the skin thickness. In addition, our population study showed that melon domesticated at multiple times from whole-genome perspective and melons in China are introduced from different routes. Selective sweeps underlying the genes related to desirable traits, haplotypes of alleles associated with agronomic traits, and the variants from resequencing data enable efficient breeding.

•

Provides a high-quality assembly for melon genome

•

Explains a considerable proportion of epidermis thickness

•

Melons in China are introduced from different routes

•

Haplotypes of alleles associated with agronomic traits enable efficient breeding

The gynoecious CmWIP1 transcription factor interacts with CmbZIP48 to inhibit carpel development

In angiosperms, sex determination leads to development of unisexual flowers. In Cucumis melo, development of unisexual male flowers results from the expression of the sex determination gene, CmWIP1, in carpel primordia. To bring new insight on the molecular mechanisms through which CmWIP1 leads to carpel abortion in male flowers, we used the yeast two-hybrid approach to look for CmWIP1-interacting proteins. We found that CmWIP1 physically interacts with an S2 bZIP transcription factor, CmbZIP48. We further determined the region mediating the interaction and showed that it involves the N-terminal part of CmWIP1. Using laser capture microdissection coupled with quantitative real-time gene expression analysis, we demonstrated that CmWIP1 and CmbZIP48 share a similar spatiotemporal expression pattern, providing the plant organ context for the CmWIP1-CmbZIP48 protein interaction. Using sex transition mutants, we demonstrated that the expression of the male promoting gene CmWIP1 correlates with the expression of CmbZIP48. Altogether, our data support a model in which the coexpression and the physical interaction of CmWIP1 and CmbZIP48 trigger carpel primordia abortion, leading to the development of unisexual male flowers.

Developmental mechanisms involved in the diversification of flowers

We all appreciate the fantastic diversity of flowers. How flowers diversified, however, remains largely enigmatic. The mechanisms underlying the diversification of flowers are complex because the overall appearance of a flower is determined by many factors, such as the shape and size of its receptacle, and the arrangement, number, type, shape and colour of floral organs. Modifications of the developmental trajectories of a flower and its components, therefore, can lead to the generation of new floral types. In this Review, by summarizing the recent progress in studying the initiation, identity determination, morphogenesis and maturation of floral organs, we present our current understanding of the mechanisms underlying the diversification of flowers.

Flower bud proteome reveals modulation of sex-biased proteins potentially associated with sex expression and modification in dioecious Coccinia grandis

BACKGROUND

Dioecy is an important sexual system wherein, male and female flowers are borne on separate unisexual plants. Knowledge of sex-related differences can enhance our understanding in molecular and developmental processes leading to unisexual flower development. Coccinia grandis is a dioecious species belonging to Cucurbitaceae, a family well-known for diverse sexual forms. Male and female plants have 22A + XY and 22A + XX chromosomes, respectively. Previously, we have reported a gynomonoecious form (22A + XX) of C. grandis bearing morphologically hermaphrodite flowers (GyM-H) and female flowers (GyM-F). Also, we have showed that foliar spray of AgNO₃ on female plant induces morphologically hermaphrodite bud development (Ag-H) despite the absence of Y-chromosome.

RESULTS

To identify sex-related differences, total proteomes from male, female, GyM-H and Ag-H flower buds at early and middle stages of development were analysed by label-free proteomics. Protein search against the cucumber protein sequences (Phytozome) as well as in silico translated C. grandis flower bud transcriptome database, resulted in the identification of 2426 and 3385 proteins (FDR ≤ 1%), respectively. The latter database was chosen for further analysis as it led to the detection of higher number of proteins. Identified proteins were annotated using BLAST2GO pipeline. SWATH-MS-based comparative abundance analysis between Female_Early_vs_Male_Early, Ag_Early_vs_Female_Early, GyM-H_Middle_vs_Male_Middle and Ag_Middle_vs_ Male_Middle led to the identification of 650, 1108, 905 and 805 differentially expressed proteins, respectively, at fold change ≥1.5 and P ≤ 0.05. Ethylene biosynthesis-related candidates as highlighted in protein interaction network were upregulated in female buds compared to male buds. AgNO₃ treatment on female plant induced proteins related to pollen development in Ag-H buds. Additionally, a few proteins governing pollen germination and tube growth were highly enriched in male buds compared to Ag-H and GyM-H buds.

CONCLUSION

Overall, current proteomic analysis provides insights in the identification of key proteins governing dioecy and unisexual flower development in cucurbitaceae, the second largest horticultural family in terms of economic importance. Also, our results suggest that the ethylene-mediated stamen inhibition might be conserved in dioecious C. grandis similar to its monoecious cucurbit relatives. Further, male-biased proteins associated with pollen germination and tube growth identified here can help in understanding pollen fertility.

Morphological, Transcriptomic and Hormonal Characterization of Trimonoecious and Subandroecious Pumpkin (Cucurbita maxima) Suggests Important Roles of Ethylene in Sex Expression

Sex expression is a complex process, and in-depth knowledge of its mechanism in pumpkin is important. In this study, young shoot apices at the one-true-leaf stage and 10-leaf stage in Cucurbita maxima trimonoecious line ‘2013–12’ and subandroecious line ‘9–6’ were collected as materials, and transcriptome sequencing was performed using an Illumina HiSeq^TM 2000 System. 496 up-regulated genes and 375 down-regulated genes were identified between shoot apices containing mostly male flower buds and only female flower buds. Based on gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the differentially expressed genes were mainly enriched in the ethylene and auxin synthesis and signal transduction pathways. In addition, shoot apices at the 4-leaf stage were treated with the ethylene-releasing agent 2-chloroethylphosphonic acid (Ethrel), aminoethoxyvinyl glycine (AVG), AgNO₃ and indoleacetic acid (IAA). The number of female flowers up to node 20 on the main stem of ‘2013–12’ increased significantly after Ethrel and IAA treatment and decreased significantly after AVG and AgNO₃ treatment. The female flowers in ‘9–6’ showed slight changes after treatment with the exogenous chemicals. The expression of key genes in ethylene synthesis and signal transduction (CmaACS7, CmaACO1, CmaETR1 and CmaEIN3) was determined using quantitative RT-PCR, and the expression of these four genes was positively correlated with the number of female flowers in ‘2013–12’. The variations in gene expression, especially that of CmaACS7, after chemical treatment were small in ‘9–6’. From stage 1 (S1) to stage 7 (S7) of flower development, the expression of CmaACS7 in the stamen was much lower than that in the ovary, stigma and style. These transcriptome data and chemical treatment results indicated that IAA might affect pumpkin sex expression by inducing CmaACS7 expression and indirectly affecting ethylene production, and the ethylene synthesis and signal transduction pathways play crucial roles in pumpkin flower sex expression. A possible reason for the differences in sex expression between pumpkin lines ‘2013–12’ and ‘9–6’ was proposed based on the key gene expression. Overall, these transcriptome data and chemical treatment results suggest important roles for ethylene in pumpkin sex expression.

Comparative transcriptome analysis reveals new molecular pathways for cucumber genes related to sex determination

Transcriptome data and qPCR analysis revealed new insight into genes regulatory mechanism related to cucumber sex determination. Cucumber (Cucumis sativus L.) is an economically important crop cultivated worldwide. Enhancing the genomic resources for cucumber may enable the regulation of traits relevant to crop productivity and quality. Sequencing technologies and bioinformatics tools provide opportunities for the development of such resources. The aims of this study were to identify and characterize the genes involved in sex determination and flower morphogenesis in cucumber isogenic lines that differed regarding flower sex type. We obtained transcripts for 933 genes related to shoot apex development, among which 310 were differentially expressed genes (DEGs) among the male, female, and hermaphroditic lines. We performed gene ontology and molecular network analyses and explored the DEGs related to already known processes like: hormone synthesis and signaling, lipid and sugar metabolism; and also newly discovered processes related to cell wall, membrane, and cytoskeleton modifications; ion homeostasis which appears to be important for ethylene perception and signaling, and genes expression mediated by transcription factors related to floral organ identities. We proposed a new model of regulatory mechanism network of sex development in cucumber. Our results may be useful for clarifying the molecular genetics and the functional mechanisms underlying the sex determination processes.

Molecular basis of cucumber fruit domestication

Cucumber (Cucumis sativus L.) is an economically important vegetable crop that is cultivated worldwide. Compared to the wild ancestor bearing small, bitter and seedy fruit, domesticated cucumbers exhibit significant variation in fruit appearance, size and flavor. Understanding the molecular basis of domestication related traits can provide insights into fruit evolution and make crop breeding more efficient. Here we review recent advances in relating to the genetic basis of fruit morphological traits (femaleness, fruit spine, wart, size, color and carpel development) and organoleptic features (bitterness) during cucumber domestication.

Gene Interactions Regulating Sex Determination in Cucurbits

The family Cucurbitaceae includes many economically important crops, such as cucumber (Cucumis sativus), melon (Cucumis melo), watermelon (Citrullus lanatus), and zucchini (Cucurbita pepo), which share homologous gene pathways that control similar phenotypes. Sex determination is a research hotspot associated with yield and quality, and the genes involved are highly orthologous and conserved in cucurbits. In the field, six normal sex types have been categorized according to the distribution of female, male, or bisexual flowers in a given plant. To date, five orthologous genes involved in sex determination have been cloned, and their various combinations and expression patterns can explain all the identified sex types. In addition to genetic mechanisms, ethylene controls sex expression in this family. Two ethylene signaling components have been identified recently, which will help us to explore the ethylene signaling-mediated interactions among sex-related genes. This review discusses recent advances relating to the mechanism of sex determination in cucurbits and the prospects for research in this area.

Sex and the flower - developmental aspects of sex chromosome evolution

Background

The evolution of dioecious plants is occasionally accompanied by the establishment of sex chromosomes: both XY and ZW systems have been found in plants. Structural studies of sex chromosomes are now being followed up by functional studies that are gradually shedding light on the specific genetic and epigenetic processes that shape the development of separate sexes in plants.

Scope

This review describes sex determination diversity in plants and the genetic background of dioecy, summarizes recent progress in the investigation of both classical and emerging model dioecious plants and discusses novel findings. The advantages of interspecies hybrids in studies focused on sex determination and the role of epigenetic processes in sexual development are also overviewed.

Conclusions

We integrate the genic, genomic and epigenetic levels of sex determination and stress the impact of sex chromosome evolution on structural and functional aspects of plant sexual development. We also discuss the impact of dioecy and sex chromosomes on genome structure and expression.

Integrative Analysis of Transcriptomic and Methylomic Data in Photoperiod-Dependent Regulation of Cucumber Sex Expression

The cucumber (Cucumis sativus) is characterized by its diversity and seasonal plasticity in sexual type. A long day length condition significantly decreased the cucumber female flower ratio by 17.7-52.9%, and the effect of photoperiod treatment is more significant under low temperature than under high temperature. Transcriptome analysis indicates that the photoperiod treatment preferentially significantly influenced flower development processes, particularly MADS-box genes in shoot apices. The long-day treatment resulted in predominantly transposable element (TE)- and gene-associated CHH-types of DNA methylation changes. Nevertheless, there was significant enrichment of CG- and CHG-types of DNA methylation changes nearing transcription start sites (TSSs)/transcription end sites (TESs) and gene bodies, respectively. Predominantly negative association between differentially methylated regions (DMRs) and differentially expressed genes (DEGs) were observed which implied epiregulation of DEGs. Two MADS-box genes that were significantly downregulated by long photoperiod showed significant hypermethylation in promoter regions that is essentially TE-rich. This study indicates MADS-box genes which are partially regulated by promoter methylation state may mediate photoperiod-dependent regulation of cucumber sex expression.

Temperature and photoperiod changes affect cucumber sex expression by different epigenetic regulations

BACKGROUND

Cucumbers (Cucumis sativus) are known for their plasticity in sex expression. DNA methylation status determines gene activity but is susceptible to environmental condition changes. Thus, DNA methylation-based epigenetic regulation may at least partially account for the instability of cucumber sex expression. Do temperature and photoperiod that are the two most important environmental factors have equal effect on cucumber sex expression by similar epigenetic regulation mechanism? To answer this question, we did a two-factor experiment of temperature and photoperiod and generated methylome and transcriptome data from cucumber shoot apices.

RESULTS

The seasonal change in the femaleness of a cucumber core germplasm collection was investigated over five consecutive years. As a result, 71.3% of the 359 cucumber accessions significantly decreased their femaleness in early autumn when compared with spring. High temperature and long-day photoperiod treatments, which mimic early autumn conditions, are both unfavorable for female flower formation, and temperature is the predominant factor. High temperatures and long-day treatments both predominantly resulted in hypermethylation compared to demethylation, and temperature effect was decisive. The targeted cytosines shared in high-temperature and long-day photoperiod treatment showed the same change in DNA methylation level. Moreover, differentially expressed TEs (DETs) and the predicted epiregulation sites were clustered across chromosomes, and importantly, these sites were reproducible among different treatments. Essentially, the photoperiod treatment preferentially and significantly influenced flower development processes, while temperature treatment produced stronger responses from phytohormone-pathway-related genes. Cucumber AGAMOUS was likely epicontrolled exclusively by photoperiod while CAULIFLOWER A and CsACO₃ were likely epicontrolled by both photoperiod and temperature.

CONCLUSIONS

Seasonal change of sex expression is a germplasm-wide phenomenon in cucumbers. High temperature and long-day photoperiod might have the same effect on the methylome via the same mechanism of gene-TE interaction but resulted in different epicontrol sites that account for different mechanisms between temperature- and photoperiod-dependent sex expression changes.

Comprehensive Analysis of Cucumber Gibberellin Oxidase Family Genes and Functional Characterization of CsGA20ox1 in Root Development in Arabidopsis

Cucumber (Cucumis sativus L.) is an important vegetable crop worldwide and gibberellins (GAs) play important roles in the regulation of cucumber developmental and growth processes. GA oxidases (GAoxs), which are encoded by different gene subfamilies, are particularly important in regulating bioactive GA levels by catalyzing the later steps in the biosynthetic pathway. Although GAoxs are critical enzymes in GA synthesis pathway, little is known about GAox genes in cucumber, in particular about their evolutionary relationships, expression profiles and biological function. In this study, we identified 17 GAox genes in cucumber genome and classified them into five subfamilies based on a phylogenetic tree, gene structures, and conserved motifs. Synteny analysis indicated that the tandem duplication or segmental duplication events played a minor role in the expansion of cucumber GA2ox, GA3ox and GA7ox gene families. Comparative syntenic analysis combined with phylogenetic analysis provided deep insight into the phylogenetic relationships of CsGAox genes and suggested that protein homology CsGAox are closer to AtGAox than OsGAox. In addition, candidate transcription factors BBR/BPC (BARLEY B RECOMBINANT/BASIC PENTACYSTEINE) and GRAS (GIBBERELLIC ACID-INSENSITIVE, REPRESSOR of GAI, and SCARECROW) which may directly bind promoters of CsGAox genes were predicted. Expression profiles derived from transcriptome data indicated that some CsGAox genes, especially CsGA20ox1, are highly expressed in seedling roots and were down-regulated under GA₃ treatment. Ectopic over-expression of CsGA20ox1 in Arabidopsis significantly increased primary root length and lateral root number. Taken together, comprehensive analysis of CsGAoxs would provide a basis for understanding the evolution and function of the CsGAox family.

Transcriptomic insights into the blue light-induced female floral sex expression in cucumber (Cucumis sativus L.)

In cucurbitaceous crops, sex differentiation of flower buds is a crucial developmental process that directly affects fruit yield. Here we showed that the induction of female flower was the highest in the blue light-treated monoecious cucumber plants compared with that in other light qualities (white, green and red). High-throughput RNA-Seq analysis of the shoot apexes identified a total of 74 differently-expressed genes (DEGs), in which 52 up-regulated and 22 down-regulated under the blue light compared with that in white light. The DEGs were mainly involved in metabolic pathways, biosynthesis of secondary metabolites, plant hormone signal transduction, starch and sucrose metabolism and phenylpropanoid biosynthesis. While the  ethylene and gibberellins synthesis and signaling related genes were down-regulated, the abscisic acid and auxin signal transduction pathways were up-regulated by the blue light treatment. Furthermore, the blue light treatment up-regulated the transcription of genes relating to photosynthesis, starch and sucrose metabolism. Meanwhile, the blue light suppressed the GA3 concentration but promoted the concentrations of auxin and photosynthetic pigments. Taken together, the results suggest that the blue light-induced female floral sex expression is closely associated with the blue light-induced changes in abscisic acid, auxin, gibberellins, photosynthesis, starch and sucrose metabolism pathways, which is potentially different from the traditional ethylene-dependent pathway.

Ethylene responsive factor ERF110 mediates ethylene-regulated transcription of a sex determination-related orthologous gene in two Cucumis species

In plants, unisexual flowers derived from developmental sex determination form separate stamens and pistils that facilitate cross pollination. In cucumber and melon, ethylene plays a key role in sex determination. Six sex determination-related genes have been identified in ethylene biosynthesis in these Cucumis species. The interactions among these genes are thought to involve ethylene signaling; however, the underlying mechanism of regulation remains unknown. In this study, hormone treatment and qPCR assays were used to confirm expression of these sex determination-related genes in cucumber and melon is ethylene sensitive. RNA-Seq analysis subsequently helped identify the ethylene responsive factor (ERF) gene, CsERF110, related to ethylene signaling and sex determination. CsERF110 and its melon ortholog, CmERF110, shared a conserved AP2/ERF domain and showed ethylene-sensitive expression. Yeast one-hybrid and ChIP-PCR assays further indicated that CsERF110 bound to at least two sites in the promoter fragment of CsACS11, while transient transformation analysis showed that CsERF110 and CmERF110 enhance CsACS11 and CmACS11 promoter activity, respectively. Taken together, these findings suggest that CsERF110 and CmERF110 respond to ethylene signaling, mediating ethylene-regulated transcription of CsACS11 and CmACS11 in cucumber and melon, respectively. Furthermore, the mechanism involved in its regulation is thought to be conserved in these two Cucumis species.

One Hundred Ways to Invent the Sexes: Theoretical and Observed Paths to Dioecy in Plants

Dioecy, the presence of male and female flowers on separate individuals, is both widespread and uncommon within flowering plants, with only a few percent of dioecious species spread across most major phylogenetic taxa. It is therefore safe to assume that dioecy evolved independently in these different groups, which allows us to ask questions regarding the molecular and developmental mechanisms underlying these independent transitions to dioecy. We start this review by examining the problem from the standpoint of a genetic engineer trying to develop dioecy, discuss various potential solutions, and compare them to models proposed in the past and based on genetic and evolutionary considerations. Next, we present recent information regarding candidate sex determinants in three species, acquired using newly established genomic approaches. Although such specific information is still scarce, it is slowly becoming apparent that various genes or pathways can be altered to evolve dioecy.

A RAD-Based Genetic Map for Anchoring Scaffold Sequences and Identifying QTLs in Bitter Gourd (Momordica charantia)

Genetic mapping is a basic tool necessary for anchoring assembled scaffold sequences and for identifying QTLs controlling important traits. Though bitter gourd (Momordica charantia) is both consumed and used as a medicinal, research on its genomics and genetic mapping is severely limited. Here, we report the construction of a restriction site associated DNA (RAD)-based genetic map for bitter gourd using an F2 mapping population comprising 423 individuals derived from two cultivated inbred lines, the gynoecious line 'K44' and the monoecious line 'Dali-11.' This map comprised 1,009 SNP markers and spanned a total genetic distance of 2,203.95 cM across the 11 linkage groups. It anchored a total of 113 assembled scaffolds that covered about 251.32 Mb (85.48%) of the 294.01 Mb assembled genome. In addition, three horticulturally important traits including sex expression, fruit epidermal structure, and immature fruit color were evaluated using a combination of qualitative and quantitative data. As a result, we identified three QTL/gene loci responsible for these traits in three environments. The QTL/gene gy/fffn/ffn, controlling sex expression involved in gynoecy, first female flower node, and female flower number was detected in the reported region. Particularly, two QTLs/genes, Fwa/Wr and w, were found to be responsible for fruit epidermal structure and white immature fruit color, respectively. This RAD-based genetic map promotes the assembly of the bitter gourd genome and the identified genetic loci will accelerate the cloning of relevant genes in the future.

Transcriptome Analysis of Litsea cubeba Floral Buds Reveals the Role of Hormones and Transcription Factors in the Differentiation Process

BACKGROUND

Litsea cubeba (Lour.) Pers. is an important economic plant that is rich in valuable essential oil. The essential oil is often used as a raw material for perfumes, food additives, insecticides and bacteriostats. Most of the essential oil is contained in the fruit, and the quantity and quality of fruit are dependent on the flowers. To explore the molecular mechanism of floral bud differentiation, high-throughput RNA sequencing was used to detect differences in the gene expression of L. cubeba female and male floral buds at three differentiation stages.

RESULTS

This study obtained 160.88 Gbp of clean data that were assembled into 100,072 unigenes, and a total of 38,658 unigenes were annotated. A total of 27,521 simple sequence repeats (SSRs) were identified after scanning the assembled transcriptome, and the mono-nucleotide repeats were predominant, followed by di-nucleotide and tri-nucleotide repeats. A total of 12,559 differentially expressed genes (DEGs) were detected from the female (F) and male (M) floral bud comparisons. The gene ontology (GO) databases revealed that these DEGs were primarily contained in "metabolic processes", "cellular processes", and "single-organism processes". The Kyoto Encyclopedia of Genes and Genomes (KEGG) databases suggested that the DEGs belonged to "plant hormone signal transduction" and accounted for a relatively large portion in all of these comparisons. We analyzed the expression level of plant hormone-related genes and detected the contents of several relevant plant hormones in different stages. The results revealed that the dynamic changes in each hormone content were almost consistent with the expression levels of relevant genes. The transcription factors selected from the DEGs were analyzed. Most DEGs of MADS-box were upregulated and most DEGs of bZIP were downregulated. The expression trends of the DEGs were nearly identical in female and male floral buds, and qRT-PCR analysis revealed consistency with the transcriptome data.

CONCLUSIONS

We sequenced and assembled a high-quality L. cubeba floral bud transcriptome, and the data appeared to be well replicated (n = 3) over three developmental time points during flower development. Our study explored the changes in the contents of several plant hormones during floral bud differentiation using biochemical and molecular biology techniques, and the changes in expression levels of several flower development related transcription factors. These results revealed the role of these factors (i.e., hormones and transcription factors) and may advance our understanding of their functions in flower development in L. cubeba.

CRISPR/Cas9-Mediated Mutagenesis of Four Putative Symbiosis Genes of the Tropical Tree Parasponia andersonii Reveals Novel Phenotypes

Parasponia represents five fast-growing tropical tree species in the Cannabaceae and is the only plant lineage besides legumes that can establish nitrogen-fixing nodules with rhizobium. Comparative analyses between legumes and Parasponia allows identification of conserved genetic networks controlling this symbiosis. However, such studies are hampered due to the absence of powerful reverse genetic tools for Parasponia. Here, we present a fast and efficient protocol for Agrobacterium tumefaciens-mediated transformation and CRISPR/Cas9 mutagenesis of Parasponia andersonii. Using this protocol, knockout mutants are obtained within 3 months. Due to efficient micro-propagation, bi-allelic mutants can be studied in the T0 generation, allowing phenotypic evaluation within 6 months after transformation. We mutated four genes - PanHK4, PanEIN2, PanNSP1, and PanNSP2 - that control cytokinin, ethylene, or strigolactone hormonal networks and that in legumes commit essential symbiotic functions. Knockout mutants in Panhk4 and Panein2 displayed developmental phenotypes, namely reduced procambium activity in Panhk4 and disturbed sex differentiation in Panein2 mutants. The symbiotic phenotypes of Panhk4 and Panein2 mutant lines differ from those in legumes. In contrast, PanNSP1 and PanNSP2 are essential for nodule formation, a phenotype similar as reported for legumes. This indicates a conserved role for these GRAS-type transcriptional regulators in rhizobium symbiosis, illustrating the value of Parasponia trees as a research model for reverse genetic studies.

GLABROUS (CmGL) encodes a HD-ZIP IV transcription factor playing roles in multicellular trichome initiation in melon

Map-based cloning identified CmGL that encodes a HD-ZIP type IV transcription factor that controls multicellular trichome initiation in melon. Trichomes are small hairs covering the aerial parts of plants that originate from the epidermal cells, which can protect plants against the damage by insects and pathogens. The regulatory pathway of unicellular trichomes has been well studied in the model plant Arabidopsis. Little is known about the genetic control and regulation of trichome development in melon (Cucumis melo L.) which has multicellular trichomes. In this study, we identified a melon mutant, cmgl, which showed completely glabrous on all aerial organs. A bulked segregant analysis was conducted to identify polymorphic markers for linkage analysis in a population with 256 F₂ plants, which allowed to locate the cmgl locus in melon chromosome VIII. Next-generation sequencing-aided marker discovery and fine mapping in a large population with 1536 F₂ plants narrowed the candidate gene region to 12 kb that harbored only one candidate gene for cmgl, which encoded a class IV homeodomain-associated leucine zipper transcription factor. Four SNPs in the coding region of the CmGL gene were identified between the two parental lines; a single base substitution from C to A resulted in a premature termination codon and a truncated protein in the cmgl. The SNP was converted into a dCAPS marker, which showed co-segregation in the F₂ population and 564 melon accessions. Result of this study will be helpful for better understanding of genetic control of trichome development in melon and marker-assisted selection in developing new cultivars.

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.

Differential Gene Expression Caused by the F and M Loci Provides Insight Into Ethylene-Mediated Female Flower Differentiation in Cucumber

In cucumber (Cucumis sativus L.), the differentiation and development of female flowers are important processes that directly affect the fruit yield and quality. Sex differentiation is mainly controlled by three ethylene synthase genes, F (CsACS1G), M (CsACS2), and A (CsACS11). Thus, ethylene plays a key role in the sex differentiation in cucumber. The "one-hormone hypothesis" posits that F and M regulate the ethylene levels and initiate female flower development in cucumber. Nonetheless, the precise molecular mechanism of this process remains elusive. To investigate the mechanism by which F and M regulate the sex phenotype, three cucumber near-isogenic lines, namely H34 (FFmmAA, hermaphroditic), G12 (FFMMAA, gynoecious), and M12 (ffMMAA, monoecious), with different F and M loci were generated. The transcriptomic analysis of the apical shoots revealed that the expression of the B-class floral homeotic genes, CsPI (Csa4G358770) and CsAP3 (Csa3G865440), was immensely suppressed in G12 (100% female flowers) but highly expressed in M12 (∼90% male flowers). In contrast, CAG2 (Csa1G467100), which is an AG-like C-class floral homeotic gene, was specifically highly expressed in G12. Thus, the initiation of female flowers is likely to be caused by the downregulation of B-class and upregulation of C-class genes by ethylene production in the floral primordium. Additionally, CsERF31, which was highly expressed in G12, showed temporal and spatial expression patterns similar to those of M and responded to the ethylene-related chemical treatments. The biochemical experiments further demonstrated that CsERF31 could directly bind the promoter of M and promote its expression. Thus, CsERF31 responded to the ethylene signal derived from F and mediated the positive feedback regulation of ethylene by activating M expression, which offers an extended "one-hormone hypothesis" of sex differentiation in cucumber.

De novo transcriptome assembly from flower buds of dioecious, gynomonoecious and chemically masculinized female Coccinia grandis reveals genes associated with sex expression and modification

BACKGROUND

Coccinia grandis (ivy gourd), is a dioecious member of Cucurbitaceae having heteromorphic sex chromosomes. Chromosome constitution of male and female plants of C. grandis is 22A + XY and 22A + XX respectively. Earlier we showed that a unique gynomonoecious form of C. grandis (22A + XX) also exists in nature bearing morphologically hermaphrodite flowers (GyM-H). Additionally, application of silver nitrate (AgNO₃) on female plants induces stamen development leading to the formation of morphologically hermaphrodite flowers (Ag-H) despite the absence of Y-chromosome. Due to the unavailability of genome sequence and the slow pace at which sex-linked genes are identified, sex expression and modification in C. grandis are not well understood.

RESULTS

We have carried out a comprehensive RNA-Seq study from early-staged male, female, GyM-H, and Ag-H as well as middle-staged male and GyM-H flower buds. A de novo transcriptome was assembled using Trinity and annotated by BLAST2GO and Trinotate pipelines. The assembled transcriptome consisted of 467,233 'Trinity Transcripts' clustering into 378,860 'Trinity Genes'. Female_Early_vs_Male_Early, Ag_Early_vs_Female_Early, and GyM-H_Middle_vs_Male_Middle comparisons exhibited 35,694, 3574, and 14,954 differentially expressed transcripts respectively. Further, qRT-PCR analysis of selected candidate genes validated digital gene expression profiling results. Interestingly, ethylene response-related genes were found to be upregulated in female buds compared to male buds. Also, we observed that AgNO₃ treatment suppressed ethylene responses in Ag-H flowers by downregulation of ethylene-responsive transcription factors leading to stamen development. Further, GO terms related to stamen development were enriched in early-staged male, GyM-H, and Ag-H buds compared to female buds supporting the fact that stamen growth gets arrested in female flowers.

CONCLUSIONS

Suppression of ethylene responses in both male and Ag-H compared to female buds suggests a probable role of ethylene in stamen suppression similar to monoecious cucurbits such as melon and cucumber. Also, pollen fertility associated GO terms were depleted in middle-staged GyM-H buds compared to male buds indicating the necessity of Y-chromosome for pollen fertility. Overall, this study would enable identification of new sex-biased genes for further investigation of stamen arrest, pollen fertility, and AgNO₃-mediated sex modification.

The asparagus genome sheds light on the origin and evolution of a young Y chromosome

Sex chromosomes evolved from autosomes many times across the eukaryote phylogeny. Several models have been proposed to explain this transition, some involving male and female sterility mutations linked in a region of suppressed recombination between X and Y (or Z/W, U/V) chromosomes. Comparative and experimental analysis of a reference genome assembly for a double haploid YY male garden asparagus (Asparagus officinalis L.) individual implicates separate but linked genes as responsible for sex determination. Dioecy has evolved recently within Asparagus and sex chromosomes are cytogenetically identical with the Y, harboring a megabase segment that is missing from the X. We show that deletion of this entire region results in a male-to-female conversion, whereas loss of a single suppressor of female development drives male-to-hermaphrodite conversion. A single copy anther-specific gene with a male sterile Arabidopsis knockout phenotype is also in the Y-specific region, supporting a two-gene model for sex chromosome evolution.

Several models have been proposed to explain the emergence of sex chromosomes. Here, through comparative genomics and mutant analysis, Harkess et al. show that linked but separate genes on the Y chromosome are responsible for sex determination in Asparagus, supporting a two-gene model for sex chromosome evolution.

The Evolution of Sex Chromosomes and Dosage Compensation in Plants

Plant sex chromosomes can be vastly different from those of the few historical animal model organisms from which most of our understanding of sex chromosome evolution is derived. Recently, we have seen several advancements from studies on green algae, brown algae, and land plants that are providing a broader understanding of the variable ways in which sex chromosomes can evolve in distant eukaryotic groups. Plant sex-determining genes are being identified and, as expected, are completely different from those in animals. Species with varying levels of differentiation between the X and Y have been found in plants, and these are hypothesized to be representing different stages of sex chromosome evolution. However, we are also finding that sex chromosomes can remain morphologically unchanged over extended periods of time. Where degeneration of the Y occurs, it appears to proceed similarly in plants and animals. Dosage compensation (a phenomenon that compensates for the consequent loss of expression from the Y) has now been documented in a plant system, its mechanism, however, remains unknown. Research has also begun on the role of sex chromosomes in sexual conflict resolution, and it appears that sex-biased genes evolve similarly in plants and animals, although the functions of these genes remain poorly studied. Because the difficulty in obtaining sex chromosome sequences is increasingly being overcome by methodological developments, there is great potential for further discovery within the field of plant sex chromosome evolution.

Flower Development and Sex Determination between Male and Female Flowers in Vernicia fordii

Vernicia fordii is a monoecious and diclinous species with male and female flowers on the same inflorescence. Low female to male flower ratio is one of the main reasons for low yield in this species. However, little is known of its floral development and sex determination. Here, according to the results of scanning electron microscopy and histological analysis, the floral development of V. fordii was divided into 12 stages and the first morphological divergence between the male and female flowers was found to occur at stage 7. The male flowers are always unisexual, but the female flowers present bisexual characteristics, with sterile stamen (staminode) restricted to pre-meiosis of mother sporogenous cells and cell death occurring at later development stages. To further elucidate the molecular mechanism underling sex determination at the divergence stage for male and female flowers, comparative transcriptome analysis was performed. In total, 56,065 unigenes were generated and 608 genes were differentially expressed between male and female flowers, among which 310 and 298 DEGs (differentially expressed genes) showed high expression levels in males and females, respectively. The transcriptome data showed that the sexual dimorphism of female flowers was affected by jasmonic acid, transcription factors, and some genes related to the floral meristem activity. Ten candidate genes showed consistent expression in the qRT-PCR validation and DEGs data. In this study, we provide developmental characterization and transcriptomic information for better understanding of the development of unisexual flowers and the regulatory networks underlying the mechanism of sex determination in V. fordii, which would be helpful in the molecular breeding of V. fordii to improve the yield output.

The quest for epigenetic regulation underlying unisexual flower development in Cucumis melo

BACKGROUND

Melon (Cucumis melo) is an important vegetable crop from the Cucurbitaceae family and a reference model specie for sex determination, fruit ripening and vascular fluxes studies. Nevertheless, the nature and role of its epigenome in gene expression regulation and more specifically in sex determination remains largely unknown.

RESULTS

We have investigated genome wide H3K27me3 and H3K9ac histone modifications and gene expression dynamics, in five melon organs. H3K9ac and H3K27me3 were mainly distributed along gene-rich regions and constrained to gene bodies. H3K9ac was preferentially located at the TSS, whereas H3K27me3 distributed uniformly from TSS to TES. As observed in other species, H3K9ac and H3K27me3 correlated with high and low gene expression levels, respectively. Comparative analyses of unisexual flowers pointed out sex-specific epigenetic states of TFs involved in ethylene response and flower development. Chip-qPCR analysis of laser dissected carpel and stamina primordia, revealed sex-specific histone modification of MADS-box genes. Using sex transition mutants, we demonstrated that the female promoting gene, CmACS11, represses the expression of the male promoting gene CmWIP1 via deposition of H3K27me3.

CONCLUSIONS

Our findings reveal the organ-specific landscapes of H3K9ac and H3K27me3 in melon. Our results also provide evidence that the sex determination genes recruit histone modifiers to orchestrate unisexual flower development in monoecious species.

The association of changes in DNA methylation with temperature-dependent sex determination in cucumber

Cucumber (Cucumis sativus L.) is characterized by its diverse and flexible sexual types. Here, we evaluated the effect of low temperature (LT) exposure on cucumber femaleness under short-day conditions. Shoot apices were subjected to whole-genome bisulfate sequencing (WGBS), mRNA-seq, and sRNA-seq. The results showed that temperature had a substantial and global impact on transposable element (TE)-related small RNA-directed DNA methylation (RdDM) mechanisms, resulting in large amounts of CHH-type cytosine demethylation. In the cucumber genome, TEs are common in regions near genes that are also subject to DNA demethylation. TE-gene interactions showed very strong reactions to LT treatment, as nearly 80% of the differentially methylated regions (DMRs) were distributed in genic regions. Demethylation near genes led to the co-ordinated expression of genes and TEs. More importantly, genome-wide de novo methylation changes also resulted in small amounts of CG- and CHG-type DMRs. Methylation changes in CG-DMRs located <600 bp from the transcription start and end sites (TSSs/TESs) negatively correlated with transcription changes in differentially expressed genes (DEGs), probably indicating epiregulation. Ethylene is called the 'sex hormone' of cucumbers. We observed the up-regulation of ethylene biosynthesis-related CsACO3 and the down-regulation of an Arabidopsis RAP2.4-like ethylene-responsive (AP2/ERF) transcription factor, demonstrating the inferred epiregulation. Our study characterized the response of the apex methylome to LT and predicted the possible epiregulation of temperature-dependent sex determination (TSD) in cucumber.

Draft genome sequence of bitter gourd (Momordica charantia), a vegetable and medicinal plant in tropical and subtropical regions

Bitter gourd (Momordica charantia) is an important vegetable and medicinal plant in tropical and subtropical regions globally. In this study, the draft genome sequence of a monoecious bitter gourd inbred line, OHB3-1, was analyzed. Through Illumina sequencing and de novo assembly, scaffolds of 285.5 Mb in length were generated, corresponding to ∼84% of the estimated genome size of bitter gourd (339 Mb). In this draft genome sequence, 45,859 protein-coding gene loci were identified, and transposable elements accounted for 15.3% of the whole genome. According to synteny mapping and phylogenetic analysis of conserved genes, bitter gourd was more related to watermelon (Citrullus lanatus) than to cucumber (Cucumis sativus) or melon (C. melo). Using RAD-seq analysis, 1507 marker loci were genotyped in an F₂ progeny of two bitter gourd lines, resulting in an improved linkage map, comprising 11 linkage groups. By anchoring RAD tag markers, 255 scaffolds were assigned to the linkage map. Comparative analysis of genome sequences and predicted genes determined that putative trypsin-inhibitor and ribosome-inactivating genes were distinctive in the bitter gourd genome. These genes could characterize the bitter gourd as a medicinal plant.

On the rarity of dioecy in flowering plants

Dioecy, the coexistence of separate male and female individuals in a population, is a rare but phylogenetically widespread sexual system in flowering plants. While research has concentrated on why and how dioecy evolves from hermaphroditism, the question of why dioecy is rare, despite repeated transitions to it, has received much less attention. Previous phylogenetic and theoretical studies have suggested that dioecy might be an evolutionary dead end. However, recent research indicates that the phylogenetic support for this hypothesis is attributable to a methodological bias and that there is no evidence for reduced diversification in dioecious angiosperms. The relative rarity of dioecy thus remains a puzzle. Here, we review evidence for the hypothesis that dioecy might be rare not because it is an evolutionary dead end, but rather because it easily reverts to hermaphroditism. We review what is known about transitions between hermaphroditism and dioecy, and conclude that there is an important need to consider more widely the possibility of transitions away from dioecy, both from an empirical and a theoretical point of view, and by combining tools from molecular evolution and insights from ecology.

Identification of RAN1 orthologue associated with sex determination through whole genome sequencing analysis in fig (Ficus carica L.)

With the aim of identifying sex determinants of fig, we generated the first draft genome sequence of fig and conducted the subsequent analyses. Linkage analysis with a high-density genetic map established by a restriction-site associated sequencing technique, and genome-wide association study followed by whole-genome resequencing analysis identified two missense mutations in RESPONSIVE-TO-ANTAGONIST1 (RAN1) orthologue encoding copper-transporting ATPase completely associated with sex phenotypes of investigated figs. This result suggests that RAN1 is a possible sex determinant candidate in the fig genome. The genomic resources and genetic findings obtained in this study can contribute to general understanding of Ficus species and provide an insight into fig’s and plant’s sex determination system.

An SNP-based saturated genetic map and QTL analysis of fruit-related traits in Zucchini using Genotyping-by-sequencing

Background

Cucurbita pepo is a cucurbit with growing economic importance worldwide. Zucchini morphotype is the most important within this highly variable species. Recently, transcriptome and Simple Sequence Repeat (SSR)- and Single Nucleotide Polymorphism (SNP)-based medium density maps have been reported, however further genomic tools are needed for efficient molecular breeding in the species. Our objective is to combine currently available complete transcriptomes and the Zucchini genome sequence with high throughput genotyping methods, mapping population development and extensive phenotyping to facilitate the advance of genomic research in this species.

Results

We report the Genotyping-by-sequencing analysis of a RIL population developed from the inter subspecific cross Zucchini x Scallop (ssp. pepo x ssp. ovifera). Several thousands of SNP markers were identified and genotyped, followed by the construction of a high-density linkage map based on 7,718 SNPs (average of 386 markers/linkage group) covering 2,817.6 cM of the whole genome, which is a great improvement with respect to previous maps. A QTL analysis was performed using phenotypic data obtained from the RIL population from three environments. In total, 48 consistent QTLs for vine, flowering and fruit quality traits were detected on the basis of a multiple-environment analysis, distributed in 33 independent positions in 15 LGs, and each QTL explained 1.5–62.9% of the phenotypic variance. Eight major QTLs, which could explain greater than 20% of the phenotypic variation were detected and the underlying candidate genes identified.

Conclusions

Here we report the first SNP saturated map in the species, anchored to the physical map. Additionally, several consistent QTLs related to early flowering, fruit shape and length, and rind and flesh color are reported as well as candidate genes for them. This information will enhance molecular breeding in C. pepo and will assist the gene cloning underlying the studied QTLs, helping to reveal the genetic basis of the studied processes in squash.

Electronic supplementary material

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