Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

The Andromonoecious Sex Determination Gene Predates the Separation of Cucumis and Citrullus Genera

  • Adnane Boualem,

    Affiliation INRA, Institute of Plant Sciences Paris-Saclay (IPS2), CNRS, Université Paris-Sud, Bâtiment 630, 91405, Orsay, France

  • Afef Lemhemdi,

    Affiliation INRA, Institute of Plant Sciences Paris-Saclay (IPS2), CNRS, Université Paris-Sud, Bâtiment 630, 91405, Orsay, France

  • Marie-Agnes Sari,

    Affiliation CNRS, UMR 8601, Laboratoire de Chimie et Biochimie Pharmacologiques et Toxicologiques, Université René Descartes, 75006, Paris, France

  • Sarah Pignoly,

    Affiliation INRA, UR 1052, Unité de Génétique et d’Amélioration des Fruits et Légumes, BP 94, 84143, Montfavet, France

  • Christelle Troadec,

    Affiliation INRA, Institute of Plant Sciences Paris-Saclay (IPS2), CNRS, Université Paris-Sud, Bâtiment 630, 91405, Orsay, France

  • Fadi Abou Choucha,

    Affiliation INRA, Institute of Plant Sciences Paris-Saclay (IPS2), CNRS, Université Paris-Sud, Bâtiment 630, 91405, Orsay, France

  • Ilknur Solmaz,

    Affiliation Department of Horticulture, Faculty of Agriculture, University of Cukurova, Adana, Turkey

  • Nebahat Sari,

    Affiliation Department of Horticulture, Faculty of Agriculture, University of Cukurova, Adana, Turkey

  • Catherine Dogimont,

    Affiliation INRA, UR 1052, Unité de Génétique et d’Amélioration des Fruits et Légumes, BP 94, 84143, Montfavet, France

  • Abdelhafid Bendahmane

    bendahm@evry.inra.fr

    Affiliation INRA, Institute of Plant Sciences Paris-Saclay (IPS2), CNRS, Université Paris-Sud, Bâtiment 630, 91405, Orsay, France

Abstract

Understanding the evolution of sex determination in plants requires the cloning and the characterization of sex determination genes. Monoecy is characterized by the presence of both male and female flowers on the same plant. Andromonoecy is characterized by plants carrying both male and bisexual flowers. In watermelon, the transition between these two sexual forms is controlled by the identity of the alleles at the A locus. We previously showed, in two Cucumis species, melon and cucumber, that the transition from monoecy to andromonoecy results from mutations in 1-aminocyclopropane-1-carboxylic acid synthase (ACS) gene, ACS-7/ACS2. To test whether the ACS-7/ACS2 function is conserved in cucurbits, we cloned and characterized ClACS7 in watermelon. We demonstrated co-segregation of ClACS7, the homolog of CmACS-7/CsACS2, with the A locus. Sequence analysis of ClACS7 in watermelon accessions identified three ClACS7 isoforms, two in andromonoecious and one in monoecious lines. To determine whether the andromonoecious phenotype is due to a loss of ACS enzymatic activity, we expressed and assayed the activity of the three protein isoforms. Like in melon and cucumber, the isoforms from the andromonoecious lines showed reduced to no enzymatic activity and the isoform from the monoecious line was active. Consistent with this, the mutations leading andromonoecy were clustered in the active site of the enzyme. Based on this, we concluded that active ClACS7 enzyme leads to the development of female flowers in monoecious lines, whereas a reduction of enzymatic activity yields hermaphrodite flowers. ClACS7, like CmACS-7/CsACS2 in melon and cucumber, is highly expressed in carpel primordia of buds determined to develop carpels and not in male flowers. Based on this finding and previous investigations, we concluded that the monoecy gene, ACS7, likely predated the separation of the Cucumis and Citrullus genera.

Introduction

The majority of angiosperms are hermaphroditic producing exclusively bisexual flowers. Sex determination is a developmental evolutionary process that leads to unisexual flowers in 10% of the species [1]. Monoecious species (5–6%) exhibit male and female flowers on the same plant. Dioecious species (5–6%) have separate male and female individuals. There has been concerted effort in the last few decades to pinpoint the master genes controlling sex determination and to understand the evolutionary forces that drive transitions between sexual morphs [2,3]. In these investigations, several plant families have emerged as model systems [47]. The choice of the Cucurbitaceae is justified by the wide spread of sex among the species of this plant family. Out of 800 investigated species, 460 are monoecious and 340 are dioecious. Throughout evolution of Cucurbitaceae, there have been also numerous shifts between monoecy and other sexual morphs, which makes this plant family a practical model to investigate the molecular mechanisms controlling sex determination.

In monoecious Cucurbitaceae species, such as Cucumis melo, floral primordia are initially bisexual with sex determination occurring by the selective developmental arrest of either the stamen or the carpel, resulting in unisexual flowers [8]. Cloning and characterization of natural mutations leading to transition from monoecy to different sexual morphs in melon and cucumber has brought new insight to the monoecy sex determination pathway. Three master sex genes, Monoecious, Androecious and Gynoecious genes, interact to explain the gender and the position of the flower on the plant [912]. In this model the expression of the carpel inhibitor, CmWIP1, encoded by the gynoecious gene, is dependent on non-expression of CmACS11, encoded by the androecious gene. The expression of the stamina inhibitor, CmACS-7, encoded by the monoecious gene is dependent on non-expression of CmWIP1. In monoecious and andromonoecious plants, male flowers develop because of non-expression of CmACS11 that permits the expression of CmWIP1. Female flowers develop on the branches because of expression of CmACS11 that represses the expression of CmWIP1, thus releasing the expression of CmACS-7 that inhibits the stamina development. If nonfunctional CmACS-7 is expressed, hermaphrodite instead of female flowers develop. Androecious plants are obtained by loss of function of CmACS11 that leads to expression of CmWIP1 in all the flowers of the plants. Gynoecious plants are obtained by inactivation of CmWIP1 function. Hermaphrodite plants are obtained by inactivation of CmWIP1 and CmACS-7 [912].

Although sex determination occurred many times independently during evolution [1,13,14], one would expect that related species use common sexual mechanisms. We previously showed that in melon and cucumber, two species that diverged more than 10 million years ago (Mya) share the same pathway leading to monoecy [10]. For instance the melon and the cucumber monoecious gene are orthologs and encode for the rate-limiting enzyme in ethylene biosynthesis, the 1-aminocyclopropane-1-carboxylic acid synthase (ACS). In both species, the ACS genes, referred to as CmACS-7 in melon and CsACS2 in cucumber, are expressed in carpel primordia, and loss of enzymatic activity leads to stamen development [10].

Similarly, the melon and the cucumber androecious genes that encode for an ACS are also orthologs [11]. In both species, the ACS genes, referred to as CmACS11 in melon and CsACS11 in cucumber, are expressed in carpel primordia, and loss of enzymatic activity leads to stamen development [11].

In melon and cucumber, two species of the Cucumis genus, our previous studies [9,10] have shown that the sexual transition from monoecy to andromonoecy results from loss-of-function mutations in CmACS-7 and CsACS2, respectively. In watermelon, a monoecious cucurbit, the recessive a locus controls the andromonoecious sexual phenotype and acts as CmACS-7 and CsACS2 [15]. Recently, a QTL analysis has shown that the major QTL controlling the watermelon andromonoecy is located on chromosome 3 close to the Cla011230 gene encoding an ACC synthase enzyme [16].

In the present study, we investigated whether components of the monoecy pathway are conserved between the Cucumis and Citrullus genus that diverged 20 Mya [17]. To test this hypothesis, we cloned the andromonoecious (a) gene from Citrullus lanatus (watermelon) and characterized natural mutations leading to sexual transition between monoecy and andromonoecy. We showed that, like in melon and cucumber, an ACS gene, ClACS7, inhibits stamina development in female flowers and loss of enzymatic activity leads to stamen development. ClACS7 is syntenic to CmACS-7 in melon and CsACS2 in cucumber. Based on this finding, we concluded that the ACS7 function is ancestral to the divergence of the Cucumis and the Citrullus genus.

Materials and Methods

Plant material and segregating population

The genetic linkage analysis of the A locus and ClACS7 was carried out using a segregating population obtained from a cross between a monoecious accession, Citrullus lanatus var. citroides, PI 296341, and an andromonoecious Citrullus lanatus var. lanatus variety Halep Karasi. The resultant F1 was backcrossed with Halep Karasi.

Plants were grown in the greenhouse in summer at Montfavet, France, under standard agronomic conditions and evaluated for the sex of their flowers at each node along the main stem and side branches.

Plant genotyping

To identify plants carrying recombination events, plant DNA was extracted from each individual of the backcross segregating population. Determination of the sexual phenotype was performed for all plants. For genotyping, a 556 bp DNA fragment of the ClACS7 gene was amplified with primers listed in S1 Table and sequenced. Polymorphisms within this fragment were used to determine recombinant alleles.

Sequence analysis

Multiple sequence alignment of full-length protein sequences of ACS was performed using the ClustalW (http://www.ebi.ac.uk/Tools/clustalw2). Phylogenetic trees were constructed by MEGA6 (http://www.megasoftware.net/index.html) based on the Neighbor-Joining method. The three-dimensional structure predictions were generated using the Geno3D server (http://geno3d-pbil.ibcp.fr) and the visualization was carried out using the Chimera server (http://www.cgl.ucsf.edu/chimera).

Quantitative RT-PCR experiments

Total RNA was extracted from frozen leaves, stems and flowers with the Trizol reagent (Invitrogen, France). First strand cDNA synthesis and primer design were performed as decribed in [11]. Primers sequences are listed in S1 Table. To check specificity of the designed primers, all amplicons were sequenced and blasted against NCBI database. Quantitative RT-PCR reactions were performed as described in [11]. A negative control without cDNA, technical replicates on three independent synthesis of cDNA (derived from the same RNA sample), and three independent biological experiments were performed in all cases. To compare data from different PCR runs and cDNA samples, CT values for ClACS7 were normalized to the CT value of ClActin2 (primers shown in S1 Table). ClACS7 relative expressions were determined as described in [9].

In situ hybridization

ClACS7 in situ hybridisation was performed as described in [18]. Primers used for this experiment are listed in S1 Table.

Expression and purification of the recombinant ClACS7 isoforms

The BL21(DE3)pLYSS E. coli strain {F- ompT hsdSB(rB- m rB) gal dcm (DE3) pLysS (CmR)} was used for the expression of the enzyme. The plasmid pET15b (Novagen, Germany) encodes T7 promoter, and confers ampicillin resistance. S-Adenosyl Methionine (SAM), Pyridoxal 5’phosphate (PLP) and 5’Adenylic Acid Deaminase from Aspergillus (Deaminase), were purchased from Sigma (France).

The ClACS7 enzymes expression and purification were performed as described in [9], using BL21(DE3)pLYS cells expressing the pClACS7 constructs.

Enzyme activity assays

Enzyme activity was determined by monitoring the MTA formation by differential spectroscopy recorded on an Uvikon 940 spectrophotometer (Biotek-Kontron Instruments, France) as according to [19]. Namely, specific activity measurement were performed on 3 different enzyme preparations by incubating 60 μM SAM in 100 mM KPhos buffer pH 8.5 (0.2 ml) and deaminase (8 μg for 0.2 ml) in the presence of 5 μM PLP in each quartz cuvette for 3 min at 25°C after addition of the purified enzyme (1 to 2 μg). The conversion of the MTA produced by ACC into an inosine derivative was monitored at 265 nm (Δε = –7740 M-1.cm-1) and the specific activity was expressed as nmol of MTA formed per min per mg of protein.

Results

Cla011230/ClACS7 is the watermelon orthologue of CmACS-7/CsACS2

To test whether the genetic determinant controlling monoecy to andromonoecy sexual transition is conserved between the Citrullus and Cucumis genus and whether the watermelon ortholog of CmACS-7/CsACS2 could be encoded by the watermelon A gene, we isolated all the ACS genes from melon, cucumber and watermelon. Both cucurbit genomes contain 8 ACS genes and phylogenetic analysis determined that Cla011230, hereafter ClACS7, CmACS-7 and CsACS2 define a monophyletic clade (Fig 1A). ClACS7 encodes an ACC synthase of 444 amino acids that share 92% amino acid sequence identity with CmACS-7 and CsACS2 and all the polymorphic residues are located in non-conserved positions (S1 Fig). Sequence analysis shows that ClACS7, CmACS-7 and CsACS2 are highly synthenic and share conserved gene structure with 3 exons and 2 introns (Fig 1B and 1C). Based on this, we concluded that ClACS7 is most likely the watermelon ortholog of CmACS-7 and CsACS2.

thumbnail
Fig 1. Phylogeny relationship of cucurbit ACS protein family and syntenic and gene structure analysis of watermelon, ClACS7, melon, CmACS-7, and cucumber, CsACS2, genes.

(A) Using the ClACS7 protein sequence as the query in tBLASTn searches, ACS protein family members were identified from watermelon, cucumber and melon. Protein sequences were aligned and used to generate the neighbor-joining phylogenetic tree with 1,000 bootstrap replicates. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test are shown next to the branches. ClACS7 clade is indicated in red. (B) Syntenic segment between the cucumber chromosome 1 (top), the melon chromosome 2 (middle) and the watermelon chromosome 3 (bottom) harboring CsACS2 (Csa1M580750), CmACS-7 (MELO3C015444) and ClACS7 (Cla011230). Broad horizontal arrows indicated genes and their transcriptional orientations. Syntenic genes are joined by light purple lines. (C) Schematic diagram of CsACS2, CmACS-7 and ClACS7 gene structure. The numbers, in bps, indicate the size of the exons (filled boxes) and the introns (black lines).

https://doi.org/10.1371/journal.pone.0155444.g001

ClACS7 co-segregates with the watermelon andromonoecious (a) locus

To map ClACS7 relative to the andromonoecious (a) locus, we first sequenced the ClACS7 gene in four watermelon lines, Sugar Baby, Charleston Grey, PI 296341 which are monoecious (A) harboring male and female flowers(Fig 2A and 2B), and Halep Karasi, an andromonoecious (a) line harboring male and hermaphrodite flowers (Fig 2A and 2C). PI 296341 and Halep Karasi were used to generate the mapping population. Among the four accessions, we observed 19 DNA polymorphisms in the introns and 7 in the exons (Fig 2D). Out of the 7 polymorphisms in the coding region, six were silent and only the SNP at nucleotide position 1477 produced a cysteine to tryptophan amino acid substitution at position 364 in the protein (C364W; Fig 2D). The identified C364W substitution was mapped relative to (a) sex locus (Fig 2E). In this experiment, a perfect co-segregation of ClACS7 with the (a) locus was obtained in a segregating population of 78 backcross plants (Fig 2E). Based on the functional similarity between monoecy versus andromonoecy in watermelon, melon and cucumber, the high sequence identity between ClACS7, CmACS-7 and CsACS2 and the co-segregation of the A gene with ClACS7, we concluded that it is likely that the A gene encodes for ClACS7.

thumbnail
Fig 2. ClACS7 sequence variation and co-segregation analysis.

(A-C) watermelon flower types from monoecious and andromonoecious plants. Monoecious watermelon bears male (A) and female (B) flowers whereas andromonoecious bears male (A) and hermaphrodite (C) flowers. sg, stigma; st, stamen. (D) ClACS7 sequence variation in monoecious and andromonoecious watermelon lines used in the co-segregation analysis. Polymorphic sites indicated below the gene model were positioned relative to the first nucleotide of the start codon. The red colored nucleotide indicates the C to G transition leading to the amino acid substitution C364W. Genotypes of the parental lines, PI 296341 and Halep Karasi are AAGyGy and aaGyGy, respectively. A and a genotypes on the left indicate whether the accessions harbour the dominant A or the recessive a alleles at the A locus. (E) Linkage analysis of andromonoecious (A) locus with ClACS7. The observed sexual phenotype, the number of plants, the genotype at the SNP 1477 and the genotype-phenotype linkage frequency are indicated.

https://doi.org/10.1371/journal.pone.0155444.g002

Loss of ClACS7 enzymatic activity leads to andromonoecy

To understand the basis of allelic differences in the A gene, we compared the sequence of ClACS7 in 43 watermelon accessions (Fig 3) selected based on their sexual types. All tested monoecious accessions harbouring the dominant A allele presented the same ClACS7 protein sequence. In contrast, analysis of the sequence obtained from the andromonoecious accessions, harbouring the recessive a allele, revealed three missense mutations, C364W, N264K and L395F leading to two ClACS7 isoforms, C364W and N264K_L395F (Fig 3). When the closest ClACS7 homologous sequences from different plants were aligned, we observed that the residues C364 and L395 are conserved across seed plants whereas N264 is located in a non-conserved region of the protein (Fig 4A) [20]. Crystallographic studies have determined that C364 and L395 residues are located in close proximity to the active site (Fig 4B–4D) [21]. The L395 residue is located in the conserved box 7 that contains residues involved in the binding of the enzyme substrate SAM (Fig 4B–4D) [21]. Surprisingly, the L395 is located at 3 amino acids of the S399 and the S399L mutation was previously described as a natural andromonoecious m allele in cucumber [10].

thumbnail
Fig 3. Analysis of ClACS7 DNA sequence variation in 43 watermelon accessions representing all the sexual morphs of the species.

Polymorphic sites indicated below the gene model were positioned relative to the first nucleotide of the start codon. The red boxes point out the nucleotide polymorphisms leading to the missense mutations, N264K, C364W and L395F leading to andromonoecious sexual morph. Presence of male, female or hermaphrodite flowers in the accessions is indicated by (x).

https://doi.org/10.1371/journal.pone.0155444.g003

thumbnail
Fig 4. Molecular characterization of ClACS7 alleles.

(A) Amino acid alignments of ClACS7 isoform with homologous proteins from Cm (Cucumis melo), Cs (Cucumis sativus), Vv (Vitis vinifera), At (Arabidopsis thaliana), Le (Lycopersicum esculentum), Ph (Petunia hybrida), Zm (Zea mays) and Pg (Picea glauca). Numbers above the alignment indicate the amino acid positions along the ClACS7 protein. Box 7 indicates a conserved domain in ACS proteins. N264K, C364W and L395F amino acid polymorphisms are shown above the alignment. (B-D) Model of the three-dimensional structure of ClACS7. (B) Superposition of the tomato ACS 3D crystal structure (1IAY.pdb) indicated in blue and the 3D model of ClACS7 indicated in pink. The three-dimensional structure were generated using the Geno3D server (http://geno3d-pbil.ibcp.fr) and the visualization was carried out using the Chimera server (http://www.cgl.ucsf.edu/chimera). (C and D) Zoom in on the active site of ClACS7. The residues that are polymorphic in the different ClACS7 isoforms are depicted in red. The cofactor PLP is indicated by purple ball and stick. The competitive inhibitor AVG (aminoethoxyvinylglycine) is indicated by yellow ball and stick. (E and F) Biochemical characterization of ClACS7, W364, F395 and K264F395 enzyme isoforms. (E) Enzymatic activity of ClACS7 (black bars), W364 (pink bars), F395 (orange bars) and K264F395 (blue bars) protein forms. Specific activities were measured on dialyzed enzymes in the presence of 60 μM SAM and various PLP concentrations. (F) Initial velocity of the ClACS7 isoforms measured at 50 mM of PLP. Experimental data were fit into a line following Michaelis-Menten equation.

https://doi.org/10.1371/journal.pone.0155444.g004

To determine whether hermaphrodite flowers in andromonoecious plants are due to loss of a ClACS7 enzymatic activity, we expressed the different forms of the protein as six-histidine tag (His6-ClACS7) fusion proteins in Escherichia coli. The purified proteins were assessed for enzymatic activity in vitro by monitoring 5’-methylthioadenosine (MTA) formation at different PLP concentrations. The K264F395 isoform was found to be totally inactive (Fig 4E and 4F, Table 1). In the presence of high concentrations of PLP (100μM), the W364 and the ClACS7 isoforms share similar biochemical constants (Table 1). However, at a physiological PLP concentration [9], the W364 isoform has a reduced enzymatic activity; its Vmax is approximately 25% of the Vmax of the ClACS7 isoform (Fig 4E, Table 1). For the K264F395 protein form, to determine whether the amino acid substitution of the conserved residue L395 is responsible for the inactivation of the enzyme, we expressed and assessed the enzymatic activity of the F395 isoform. Like the W364 isoform, at a physiological PLP concentration, the F395 isoform has a reduced enzymatic activity whereas at 100μM PLP, the F395 and ClACS7 isoforms have similar biochemical constants (Fig 4E, Table 1). The fact that the L395F amino acid substitution does not lead to a complete ClACS7 enzymatic inactivation suggests that the N264K substitution could also reduce enzymatic activity and thus the complete loss of enzymatic activity observed in the K264F395 isoform is due to a synergic effect of the N264K and L395F substitutions.

thumbnail
Table 1. The kinetic parameters, Km and Vm, measured at 5μM and 100μM PLP for ClACS7, W364, K264F395 and F395 protein isoforms.

https://doi.org/10.1371/journal.pone.0155444.t001

Based on this and because the ClACS7 isoform from the monoecious lines was active and isoforms from the andromonoecious lines showed reduced to no enzymatic activity, we concluded that the loss of ClACS7 enzyme activity is likely the cause of the development of hermaphrodite flowers in watermelon and that ethylene production by ClACS7 is necessary for the developmental inhibition of the male organs in the female flowers of monoecious watermelon. This conclusion is reinforced by recent experiments showing that monoecious watermelons treated with an ethylene biosynthesis inhibitor, the 1-aminoethoxyvinylglycine (AVG), develop hermaphrodite flowers [22]. Moreover, these data show that, in watermelon, the andromonoecy evolved from monoecy, independently in wild C. lanatus var. citroides accessions and in cultivated C. lanatus var. lanatus.

ClACS7 is expressed in carpel primordia of female and hermaphrodite flowers

Quantitative real-time PCR analysis demonstrates that ClACS7 is highly expressed in female and hermaphrodite flower buds and it is weak or undetectable in buds determined to develop into male flowers. No ClACS7 expression was detected in vegetative tissues (Fig 5A). To investigate in more detail the spatial expression pattern of ClACS7, we carried out in situ hybridizations. ClACS7 mRNA was detected in female and hermaphrodite flower buds at stages where buds are not morphologically distinguishable. ClACS7 mRNA was strongly localised in the carpel primordia (Fig 5C and 5E). In the male flowers of the monoecious or andromonoecious watermelon lines, no ClACS7 transcript was detected (Fig 5B and 5D).

thumbnail
Fig 5. Expression patterns of ClACS7 in watermelon flowers.

(A) Quantitative RT-PCR of ClACS7 in reproductive organs, male or female flower buds from monoecious watermelon and male or hermaphrodite flower buds from andromonoecious watermelon, collected at stage 4, and in the vegetative tissues, leaves (L) and stems (S); shown are the mean ± standard deviation (s.d.) of three biological replicates. (B-E) ClACS7 in situ expression in flower buds at stages 4 in male (B) and female (C) flowers of monoecious watermelon and male (D) and hermaphrodite (E) flowers of andromonoecious watermelon. C, Carpel; St, Stamen; P, Petal; S, Sepal. Bars, 100 μm.

https://doi.org/10.1371/journal.pone.0155444.g005

Since ClACS7 is expressed in both female and hermaphrodite flowers (Fig 5) and because the loss of ClACS7 activity accounts for the functional variation between the monoecious and andromonoecious watermelon lines (Fig 4), we concluded that ClACS7-mediated ethylene production in the carpel primordia prevents the development of the stamina in female flowers and is also not required for carpel development.

Moreover, since the watermelon ClACS7, the melon CmACS-7 and the cucumber CsACS2 have the same expression pattern and since the loss of ClACS7, CmACS-7 and CsACS2 enzymatic activities lead to the sexual transition from monoecy to andromonoecy, we concluded that ClACS7, CmACS-7 and CsACS2 are regulated by mechanisms common to the Citrullus and Cucumis genus and that ethylene signalling and response mechanisms leading to inhibition of stamina development in female flowers of monoecious plants are conserved in the two genera, Citrullus and Cucumis.

Discussion

Most of the species of the Cucurbitales order have unisexual flowers and are mostly monoecious (>75% of the species) bearing both male and female flowers on the same plant. Andromonoecy, a major sexual system in which plants carry both male and perfect bisexual flowers, has evolved independently numerous times [23] and is found in ~ 4,000 species in 33 angiosperm families [24]. Investigations on the genetic inheritance of the different sexual morphs in the Cucurbitaceae plant family have described that the andromonoecious trait is controlled by a recessive monogenic locus [5,18,25,26]. In our previous studies, we have demonstrated that the sexual transition from monoecy to andromonoecy, in melon and cucumber, is controlled by the common gene CmACS-7/CsACS2 [9,10]. In this context, the genetic determinant controlling the monoecy is likely to have been established before the speciation of melon and cucumber, two species of the genus Cucumis that have diverged over 10 My ago [17].

In this work, we described the isolation and characterization of the watermelon, Citrullus lanatus, andromonoecious (a) locus and showed that the (a) locus encodes for a rate-limiting enzyme in ethylene biosynthesis, ClACS7, the orthologous gene of the melon CmACS-7 and cucumber CsACS2. Sequence analysis of monoecious and andromonoecious watermelon accessions revealed two ClACS7 isoforms in the andromonoecious lines, the C364W and the N264K_L395F. Based on the kinetic parameters of the ClACS7 protein isoforms, we concluded that active ClACS7 enzyme leads to the development of female flowers in monoecious lines, whereas a reduction of enzymatic activity yields hermaphrodite flowers in andromonoecious lines. Expression pattern analysis showed that ClACS7 is specifically expressed in carpel primordia of female and hermaphrodite flowers and no expression was detected in male flowers and the vegetative tissues, leaves and stems. On the basis of enzymatic activity and gene expression studies and as previously described in melon and cucumber [9,10], we concluded that ClACS7-mediated ethylene production in the carpel primordia prevents the development of the stamina in female flowers and is also not required for carpel development.

Altogether, the conserved gene structure, the genomic synteny, the conserved expression patterns and the fact that ClACS7, CmACS-7 and CsACS2 loss-of-function activities lead to hermaphrodite flower development suggest that ClACS7, CmACS-7 and CsACS2 derive from a common ancestral gene that controls monoecy before the divergence of the Citrullus and Cucumis genera, 20 My ago [17].

Recently, a study on sex determination in squash, Cucurbita pepo, has linked the partial andromonoecy phenotype, i.e. an incomplete conversion of female into hermaphrodite flowers, to the CpACS27A gene, an ACS gene homologous to CmACS-7 and CsACS2 [27]. The similar expression pattern in the carpels of carpel-bearing flowers (female and hermaphrodite) suggest that ClACS7, CmACS-7, CsACS2 and CpACS27A share common regulatory mechanisms along evolution and that monoecy is controlled by an ancestral gene that precedes the Cucumis, Citrullus and Cucurbita genus divergence 30 My ago [17].

Supporting Information

S1 Fig. Amino acid alignments of ClACS7 with CsACS2 and CmACS-7 homologous proteins from cucumber and melon, respectively.

Box 1 to Box 7 indicate conserved domains in ACS proteins. Non-conserved amino acids between ClACS7, CsACS2 and CmACS-7 proteins are highlighted in grey.

https://doi.org/10.1371/journal.pone.0155444.s001

(PDF)

Acknowledgments

We thank P. Audigier and N. Giovinazzo for plant care and V. Rittener for technical assistance. We thank M. Chatterjee for helpful discussions on the manuscript. We thank J. Fauve for providing watermelon andromonoecious accessions. This work was supported by BAP department in INRA, the grants Program Saclay Plant Sciences (SPS, ANR-10-LABX-40), ANR MELODY (ANR-11-BSV7-0024) and the European Research Council (ERC-SEXYPARTH).

Author Contributions

Conceived and designed the experiments: A Boualem NS CD A Bendahmane. Performed the experiments: A Boualem AL MAS SP CT FAC IS. Analyzed the data: A Boualem AL MAS SP. Contributed reagents/materials/analysis tools: A Boualem FAC IS NS CD A Bendahmane. Wrote the paper: A Boualem A Bendahmane.

References

  1. 1. Renner SS (2014) The relative and absolute frequencies of angiosperm sexual systems: dioecy, monoecy, gynodioecy, and an updated online database. Am J Bot 101: 1588–1596. pmid:25326608
  2. 2. Diggle PK, Di Stilio VS, Gschwend AR, Golenberg EM, Moore RC, Russell JR, et al. (2011) Multiple developmental processes underlie sex differentiation in angiosperms. Trends Genet 27: 368–376. pmid:21962972
  3. 3. Ainsworth C, Parker J, Buchanan-Wollaston V (1998) Sex determination in plants. Curr Top Dev Biol 38: 167–223. pmid:9399079
  4. 4. Kubicki B (1969a) Investigations on sex determination in cucumber (Cucumis sativus L.). V. Genes controlling intensity of femaleness. Genetica Polonica 10: 69–85.
  5. 5. Kubicki B (1969b) Investigations on sex determination in cucumber (Cucumis sativus L.). VII. Andromonoecism and hermaphroditism. Genetica Polonica 10: 101–120.
  6. 6. Poole CF, Grimball PC (1939) Inheritance of new sex forms in Cucumis melo L. Journal of Heredity 30: 21–25.
  7. 7. Malepszy S, Niemirowicz-Szczytt K (1991) Sex determination in cucumber (Cucumis sativus) as a model system for molecular biology. Plant Science 80: 39–47.
  8. 8. Bai SL, Peng YB, Cui JX, Gu HT, Xu LY, Li YQ, et al. (2004) Developmental analyses reveal early arrests of the spore-bearing parts of reproductive organs in unisexual flowers of cucumber (Cucumis sativus L.). Planta 220: 230–240. pmid:15290297
  9. 9. Boualem A, Fergany M, Fernandez R, Troadec C, Martin A, Morin H, et al. (2008) A conserved mutation in an ethylene biosynthesis enzyme leads to andromonoecy in melons. Science 321: 836–838. pmid:18687965
  10. 10. Boualem A, Troadec C, Kovalski I, Sari MA, Perl-Treves R, Bendahmane A (2009) A conserved ethylene biosynthesis enzyme leads to andromonoecy in two cucumis species. PLoS One 4: e6144. pmid:19578542
  11. 11. Boualem A, Troadec C, Camps C, Lemhemdi A, Morin H, Sari MA, et al. (2015) A cucurbit androecy gene reveals how unisexual flowers develop and dioecy emerges. Science 350: 688–691. pmid:26542573
  12. 12. Martin A, Troadec C, Boualem A, Rajab M, Fernandez R, Morin H, et al. (2009) A transposon-induced epigenetic change leads to sex determination in melon. Nature 461: 1135–1138. pmid:19847267
  13. 13. Barrett SC (2002) The evolution of plant sexual diversity. Nat Rev Genet 3: 274–284. pmid:11967552
  14. 14. Lewis D (1942) The evolution of sex in flowering plants. Biological Review of the Cambridge Philosophical Society 17: 46–67.
  15. 15. Rosa JT (1928) The inheritance of flowers types in Cucumis and Citrullus. Hilgardia 3: 233–250.
  16. 16. Prothro J, Abdel-Haleem H, Bachlava E, White V, Knapp S, McGregor C (2013) Quantitative trait loci associated with sex expression in an inter-subspecific watermelon population. J Amer Soc Hort Sci 138: 125–130.
  17. 17. Schaefer H, Heibl C, Renner SS (2009) Gourds afloat: a dated phylogeny reveals an Asian origin of the gourd family (Cucurbitaceae) and numerous oversea dispersal events. Proc Biol Sci 276: 843–851. pmid:19033142
  18. 18. Nikovics K, Blein T, Peaucelle A, Ishida T, Morin H, Aida M, et al. (2006) The balance between the MIR164A and CUC2 genes controls leaf margin serration in Arabidopsis. Plant Cell 18: 2929–2945. pmid:17098808
  19. 19. White MF, Vasquez J, Yang SF, Kirsch JF (1994) Expression of apple 1-aminocyclopropane-1-carboxylate synthase in Escherichia coli: kinetic characterization of wild-type and active-site mutant forms. Proc Natl Acad Sci U S A 91: 12428–12432. pmid:7809054
  20. 20. Ralph SG, Hudgins JW, Jancsik S, Franceschi VR, Bohlmann J (2007) Aminocyclopropane carboxylic acid synthase is a regulated step in ethylene-dependent induced conifer defense. Full-length cDNA cloning of a multigene family, differential constitutive, and wound- and insect-induced expression, and cellular and subcellular localization in spruce and Douglas fir. Plant Physiol 143: 410–424. pmid:17122070
  21. 21. Huai Q, Xia Y, Chen Y, Callahan B, Li N, Ke H (2001) Crystal structures of 1-aminocyclopropane-1-carboxylate (ACC) synthase in complex with aminoethoxyvinylglycine and pyridoxal-5'-phosphate provide new insight into catalytic mechanisms. J Biol Chem 276: 38210–38216. pmid:11431475
  22. 22. Manzano S, Martinez C, Garcia JM, Megias Z, Jamilena M (2014) Involvement of ethylene in sex expression and female flower development in watermelon (Citrullus lanatus). Plant Physiol Biochem 85: 96–104. pmid:25463265
  23. 23. Miller JS, Diggle PK (2003) Diversification of andromonoecy in Solanum section Lasiocarpa (Solanaceae): the roles of phenotypic plasticity and architecture. American Journal of Botany 90: 707–715. pmid:21659166
  24. 24. Yampolsky E, Yampolsky H (1922) Distribution of sex forms in phanerogamic flora. Bibliographia Genetica 3: 1–62.
  25. 25. Poole CF, Grimball PC (1944) Interaction of sex, shape and weight genes in watermelon. J Agric Res 71: 533–552.
  26. 26. Ji GJ, Zhang J, Gong GY, Shi JT, Zhang HY, Ren Y, et al. (2015) Inheritance of sex forms in watermelon (Citrullus lanatus). Scientia Horticulturae 193: 367–373.
  27. 27. Martinez C, Manzano S, Megias Z, Barrera A, Boualem A, Garrido D, et al. (2014) Molecular and functional characterization of CpACS27A gene reveals its involvement in monoecy instability and other associated traits in squash (Cucurbita pepo L.). Planta 239: 1201–1215. pmid:24595516