A small XY chromosomal region explains sex determination in wild dioecious V. vinifera and the reversal to hermaphroditism in domesticated grapevines
© Picq et al.; licensee BioMed Central Ltd. 2014
Received: 10 April 2014
Accepted: 18 August 2014
Published: 3 September 2014
In Vitis vinifera L., domestication induced a dramatic change in flower morphology: the wild sylvestris subspecies is dioecious while hermaphroditism is largely predominant in the domesticated subsp. V. v. vinifera. The characterisation of polymorphisms in genes underlying the sex-determining chromosomal region may help clarify the history of domestication in grapevine and the evolution of sex chromosomes in plants. In the genus Vitis, sex determination is putatively controlled by one major locus with three alleles, male M, hermaphrodite H and female F, with an allelic dominance M > H > F. Previous genetic studies located the sex locus on chromosome 2. We used DNA polymorphisms of geographically diverse V. vinifera genotypes to confirm the position of this locus, to characterise the genetic diversity and traces of selection in candidate genes, and to explore the origin of hermaphroditism.
In V. v. sylvestris, a sex-determining region of 154.8 kb, also present in other Vitis species, spans less than 1% of chromosome 2. It displays haplotype diversity, linkage disequilibrium and differentiation that typically correspond to a small XY sex-determining region with XY males and XX females. In male alleles, traces of purifying selection were found for a trehalose phosphatase, an exostosin and a WRKY transcription factor, with strikingly low polymorphism levels between distant geographic regions. Both diversity and network analysis revealed that H alleles are more closely related to M than to F alleles.
Hermaphrodite alleles appear to derive from male alleles of wild grapevines, with successive recombination events allowing import of diversity from the X into the Y chromosomal region and slowing down the expansion of the region into a full heteromorphic chromosome. Our data are consistent with multiple domestication events and show traces of introgression from other Asian Vitis species into the cultivated grapevine gene pool.
KeywordsDioecy Domestication Hermaphroditism Sex chromosome Vitis vinifera L
The wild grapevine, Vitis vinifera L. subsp. sylvestris, is the wild ancestor of the domesticated grapevine V. v. vinifera,, cultivated for wine and table grape production . The genus Vitis, a monophyletic taxon of the family Vitaceae,, includes approximately sixty species present mainly in Asia and America, all of which -except the domesticated grapevine- are dioecious (male and female flowers borne on different plants) ,. During grapevine domestication, flower reproductive morphology has incurred radical modifications, with the change from dioecy to hermaphroditism in domesticated grapevines . The geographic origin of hermaphroditism development in the domesticated grapevine is still not elucidated, nor is it known whether it occurred during primary , and/or secondary domestication events believed to have occurred in geographically distinct areas around the Mediterranean ,.
Sex expression in Vitis flower is thought to be controlled by a major locus with three alleles, male M, hermaphrodite H and female F, with an M > H > F allelic dominance ,,[12-14]. Several genetic maps based on interspecific crosses have confirmed that sex determinism in the genus Vitis is under the control of a single major genomic region located on chromosome 2, close to the SSR marker VVIB23 [15-17]. Recently, a complex interspecific cross (V. vinifera x [V. riparia x V. cinerea]) was used by Fechter et al.  to narrow the location of the sex locus to a 143 kb genomic region located between positions 4.907.434 and 5.037.597 bp of chromosome 2  on the physical map of the V. vinifera reference genome sequence (PN40024 12x.0 version ). So far, the co-localisation on chromosome 2 of the sex locus in V. vinifera subsp. vinifera has been confirmed only in the genetic map of one intra-specific cross , with a recombination distance of 0.4 cM from the nearest genetic marker (VVIB23). Moreover, in the V. v. sylvestris subspecies, the sex locus localisation remains to be confirmed.
The evolution of proper sex chromosomes is quite rare in plants: indeed, approx. 40 species of flowering plants are currently known to have developed sex chromosomes and among them, half have heteromorphic sex chromosomes . A sex chromosome may start to develop in dioecious species through the suppression of recombination between male- and female-sterile mutations with complementary dominance in close proximity on a chromosome . Then, this sex-determination region would gradually grow in size, increasingly incorporating sex-linked genes and eventually evolving into heteromorphic sex chromosomes ,. Some of the processes involved in sex chromosome evolution, as the suppression of genetic recombination or the genetic degeneration of the Y chromosome, are not well understood and only the study of the sex-determining systems on different species and at different steps of evolution could provide some answers . While the sex determination locus in Vitis species was mainly studied to develop genetic markers for early sexing for breeding purposes ,, the work of Fechter et al. , evidencing a small sex-determination region, suggests that Vitis species could be good candidates to study the early steps of sex chromosome evolution.
In the present study, we explore the sequence polymorphisms near the sex locus in a genetically and geographically diverse panel of wild and domesticated grapevines, with the objectives to: i) confirm the position and boundaries of the sex locus in V. vinifera subsp. sylvestris; ii) characterise the sex region in terms of linkage disequilibrium, genetic diversity, selection signature and candidate genes; and iii) use this information to explore the geographic and genetic origin of hermaphroditism in domesticated grapevine.
Since wild grapevines carry the ancestral form of the sex locus from which the domesticated grapevine hermaphroditism derived, we first mapped sequence polymorphisms linked to the sex trait in Vitis vinifera subsp. sylvestris. Then, we compared the polymorphisms linked to the sex trait in diverse wild and domesticated grapevine populations to study the origin of hermaphroditism in domesticated grapevines.
Plant material and phenotypic trait data
The plant material consisted of 73 wild (39 females and 34 males) and 39 hermaphrodite domesticated grapevines (Additional file 1). These grapevines were chosen among 139 wild genotypes and 2.323 domesticated genotypes  to maximise both genetic diversity and geographic representation. Three genotypes from other species were also considered to represent genetic variation in the subgenus Vitis: V. balansaeana, V. coignetiae and V. monticola. The grapevines were sampled either in natural populations or from the French National Grapevine Germplasm Collection (INRA, Domaine de Vassal, France; http://www1.montpellier.inra.fr/vassal/). The genotypes considered varied according to the genetic analyses (Additional file 1). Sex phenotypes (male, female or hermaphrodite) were evaluated by observations of flower morphology repeated over several years, and coded according to the International Organization of Vine and Wine descriptors (code number OIV-151 ).
DNA was extracted from 150 mg of leaves according to the Dneasy Plant Mini Kit (Qiagen) instructions with 1% Polyvinylpyrrolidone (PVP 40.000) and 1% of β-mercaptoethanol added to the buffer AP1 to eliminate polyphenols, strong inhibitors of in-vitro enzymatic reactions abundantly present in the crude grape cell lysate.
Characteristics of the amplicons used in this study to cover the sex locus and its edges
4781551 - 4782603
Esterase/lipase/thioesterase family protein
4822617 - 4824068
SAUR family protein
4850582 - 4851997
Pentapeptide repeat protein
4861475 - 4862891
4883461 - 4884818
Soluble acid invertase
4900275 - 4901493
Trehalose-6-phosphate phosphatase (TPP)
4921838 - 4923352
Exostosin family protein
4953195 - 4954179**
Ethylene Overproducer-like 1 (ETO1)
4989467 - 4990268
WRKY transcription factor 21
Adenine phosphoribosyltransferase (APT3)
5036645 - 5037597
Phosphatidic acid phosphatase 2 (PAP2)
According to Fechter et al. , the 143 kb region of chromosome 2 (12×.0 version) between 4,907,434 and 5,050,616 bp corresponds to the female allele of the hermaphroditic Pinot Noir 40024, while the slightly different hermaphrodite allele is located on the unassigned scaffold_233 (chromosome UnRandom of the 12×.0). The 12×.0 scaffold_233 is collinear with the chromosome 2 of the 8× grape genome reference sequence ; both these assemblies display two regions which are absent from the chromosome 2 assembly of the 12×.0 reference sequence version: a region between the 3-Oxoacyl synthase III C terminal (KASIII) and the PLATZ transcriptor factor, and the adenine phosphoribosyl transferase (APT3) region . The APT3 distinguishes female individuals from male and hermaphroditic ones . A gene, the phosphatidic phosphatase 2 (PAP2), is not predicted by the Gaze annotation of the 12×.0 reference sequence version but it is annotated by the Gaze annotation of the 8× reference sequence version and confirmed by Fechter et al.  on the 12×.0 reference sequence version.
For our work, eight primer pairs out of the eleven could thus be designed using the Gaze annotation of the 12×.0 reference sequence version (Table 1, Figure 1). A primer pair (VSVV011) was developed in the PAP2 gene using the Gaze annotation of the 8×reference sequence version (Table 1). Another primer pair (VSVV010) was specifically developed to cover the region of the putative APT3 distinguishing female individuals from male and hermaphroditic ones . A last amplicon (VSVV008) was designed to amplify a gene present in the region between the KASSIII and the PLATZ transcriptor on the 12×.0 scaffold_233; the predicted protein of this gene blasts with an Ethylene Overproducer-like 1 gene (ETO1, blastx E-value = 4e-83). For the ETO1 and APT3 amplicons, the positions on the grape genome physical maps were estimated based upon a manual realignment of the unassigned scaffold_233 (chromosome UnRandom of the 12×.0) and the 8×reference sequence version respectively, on the chromosome 2 of the 12×.0 reference sequence version. As a consequence, in our work the 12×.0 positions of these two amplicons are approximate (Table 1).
All primer pairs were designed using the Primer3 software V.0.4.0 , so as to amplify stretches between 600 and 1.300 bp and cover a part of the promoter and the first exons and introns ,. Thermocycling consisted of an initial stringent cycle (94°C for 3 minutes followed by 12 cycles of 94°C for 30 seconds, from 65 to 56°C decreasing by 0.70°C at each cycle for 45 seconds, 72°C for 120 seconds) and additional 25 cycles of 94°C for 30 seconds, 56°C for 45 seconds, 72°C for 90-120 seconds. Sequencing was performed on PCR products purified using the AMPure® kit (Agencourt®, MA, USA); BigDye® Terminator v3.1 Cycle Sequencing Kit (Applied BioSistem, CA, USA) was used following the standard protocol and reaction products were purified with the CleanSEQ® kit (Agencourt) and read on a 3130×l Genetic Analyzer (Applied BioSystems). Raw sequence files (AB1 format) were imported, aligned and trimmed using the Staden software v.2.0.0 ; SNP calling was carried out manually using the Staden interface. Then, fasta files were exported and subsequently analysed in other softwares and pipelines.
Identification of sequence polymorphisms linked to the sex trait
Phenotypic sex inheritance in wild grapevines produces only male and female variants, with a ratio near to 1:1 in adult populations (even if some variation in sex phenotypes have been observed ,, in our sample only two morphs were found, M and F). The most parsimonious hypothesis we could make on sex inheritance in grape, based on previous observations, preliminary data analysis, and literature survey ,,,,, was that of a XY system, where, at the sex locus, the female is homozygous (XX) and the male is heterozygous (XY).
To map the sex locus on the genome, we first used a genetic association approach, looking for correlations between sex flower phenotypes and sequence polymorphisms in a panel of diverse wild genotypes from different geographic provenances (Additional file 1). However, the use of general or mixed linear models searching for association resulted in too many false positives (SNP that were correlated to sex but explained only a portion of the phenotypes). Thus, we used an approach similar to Siegismund , using Fisher tests to compare, for each polymorphism and for male and female wild grapevines separately, the expected and observed proportions of homozygous and heterozygous genotypes. The expected proportions were assumed to follow the Hardy–Weinberg law and were calculated from the allele frequencies observed in the entire population (sum of male and female individuals). The observed counts were the number of homozygous and heterozygous genotypes actually recorded in male and female grapevines. Indels were coded as present/absent (Additional file 3). Fisher tests were calculated with the fisher.test function of the R statistical software . We only considered sequence polymorphisms with less than 20% missing data and with a minimum allele frequency in the sample higher than 5%. A test was considered significant when the probability of deviation from the null hypothesis was inferior to a 0.05 P-value threshold adjusted by a Bonferroni correction for multiple hypotheses testing (0.05/n with n corresponding to the total number of studied polymorphisms).
Linkage disequilibrium in the sex region
To explore linkage disequilibrium between and within amplicons covering the sex region, we used the Measure.R2VS() function in the R package LDcorSV . r 2 VS is the square of each pairwise correlation corrected by both the relatedness and genetic structure of the sample . The sample considered here was composed of 18 male and 18 female individuals (Additional file 1). These 36 specimens were chosen among those with the least missing data, eliminating the most closely related individuals and equilibrating their geographic representation. The genetic structure matrix was calculated from a dataset of 20 SSRs  of all the wild genotypes in this study, using STRUCTURE software . We used the model with uncorrelated allele frequencies, admixture, and no prior population information, previously showed to be pertinent in grapevine . Ten STRUCTURE runs (each with 5 × 105 iterations and 5 × 105 replicates) for each K-level were obtained and compared to estimate group assignation stability. The most probable number of sub-populations was inferred based on both the similarity pattern among the 10 STRUCTURE replicates and Evanno's Δks statistics . The kinship matrix was obtained using ML-Relate software  with the same SSR markers and genotypes as above.
Diversity in M, F and Hhaplotypes and signature of selection
To compare the diversity of male, female and hermaphrodite alleles at the significant sex-linked amplicons (see Additional file 1 for the genotypes considered), the haplotypes were reconstructed using PHASE v2.1 with default parameter values ,. The attribution of individual haplotypes to the M, F and H groups (called hereafter haplogroups) were carried out with the help of haplotype trees (Additional file 4) built with a maximum likelihood method (PhyML 3.0 ) implemented in SeaView v4.3.3  and based on the Generalised Time-Reversible (GTR) model .
Allocation of 0, 1 or 2 female haplotypes ( F) to the hermaphrodite, male and female genotypes, according to the maximum likelihood trees, for the four sex linked amplicons
n = 22
n = 22
n = 22
n = 21
0 haplotype F
1 haplotype F
2 haplotypes F
n = 22
n = 22
n = 22
n = 18
0 haplotype F
1 haplotype F
2 haplotypes F
n = 24
n = 24
n = 24
n = 22
0 haplotype F
1 haplotype F
2 haplotypes F
Finally, we evaluated the intraspecific genetic differentiation between male, female and hermaphrodite haplogroups, and the interspecific differentiation between V. v sylvestris and Vitis species haplotypes, using the Fst statistics , with DnaSP v5 software as well. The Vitis species used for this statistics were V. balanseana, V. monticola and V. coignetiae.
Origin of the H haplotypes
Combining the haplotypes of the four sex-linked amplicons, the M, F and H macrohaplotypes were reconstructed. PHASE v2.1 was run again using a 100 burn-in period with 100 iterations with a thinning interval of 1 and 10 repeats. The algorithm was run several times, validating convergence. Then, to understand the origin of H haplotypes in the domesticated grapevine, a network analysis was carried out on the F, M and H macrohaplotypes using the median-joining method as described in Bandelt et al.  and implemented in Network v22.214.171.124. . A Star Contraction was run before running the network calculation.
Finally, the relationship between the network distances (in number of mutations) of the H haplotypes from the M haplogroup, and the geographic origin, grape use (table, wine or both), degree of domestication (ancient or modern cultivars ) and the genetic structure ancestry of the domesticated grapevines  were explored using an ANOVA.
Sequence polymorphisms linked to the sex trait
Among the significant polymorphisms, 28 perfectly fitted the XY sex determination model. For these polymorphisms, 100% of the male genotypes were heterozygous and 100% of the female genotypes were homozygous for the most frequent allele, i.e. for example males were A/T and females were A/A but never T/T (Figure 2c, Additional file 3). In hermaphrodite domesticated genotypes, these same polymorphisms were in the majority of cases in a heterozygous state (Additional file 5). These 28 polymorphisms, perfectly fitting the XY model, were only found in the VSVV006, VSVV007 and VSVV009 amplicons and 3 of them resulted in non-synonymous amino acid changes (38th, 61th and 66th polymorphism in VSVV006 or VSVV007, Additional file 3).
Moreover, 18 significant polymorphisms in VSVV006, VSVV007, VSVV009 and VSVV010 were only slightly deviating from the XY sex determination model, with all female genotypes homozygous for the most frequent allele and one or two non-heterozygous males (Additional file 3). For example, for the polymorphism 126 (crosses in Figure 2b, c) corresponding to the sex-linked indel in the second intron of the APT3 gene , all female were homozygous without the indel while 92% of male were heterozygous (23 heterozygous, one homozygous with the indel and one homozygous without it) (Additional file 3). In the VSVV008 amplicon, only one SNP was found slightly deviating from the XY sex determination model (Figure 2b and Additional file 3).
By contrast, and although few of them were found significantly departing from Hardy-Weinberg proportions (Fisher test), the polymorphisms found in VSVV002, VSVV003, VSVV004 and VSVV005, largely deviated from the XY model, particularly in male genotypes (Figure 2 and Additional file 3).
In summary, 46 significant polymorphisms in the VSVV006, VSVV007, VSVV009 and VSVV010 amplicons fitted either strictly (28) or closely (18) the XY sex-determination model. These results allowed us to define the boundaries of the sex locus at the positions 4.884.818 and 5.036.645 on chromosome 2 of the PN40024 physical map (12×.0 version). This 151,83 kb region, externally delimited by the gene fragments VSVV005 and VSVV011 contains 13 candidate genes (Figure 1 and Additional file 6).
Linkage disequilibrium in the sex region
The intra- and inter-amplicon linkage disequilibrium (LD) was estimated on a sub-sample of 18 male and 18 female wild grapevines (Additional file 1), by calculating the pairwise square correlation coefficient r 2 VS, correcting for the structure and kinship of the sample. Only sequence polymorphisms with less than 20% missing data and with a 0.2 minor allele frequency were analysed. At these thresholds, no polymorphisms were retained in the VSVV001 fragment.
Diversity of the M, F and Hhaplotypes and signature of selection
The M, F and H haplotypes for the four sex-associated amplicons (VSVV006, VSVV007, VSVV009 and VSVV010) were assigned using maximum likelihood haplotype trees. According to the XY model and the rules of dominance described for Vitis (M > H > F,,[12-14]), the haplogroup containing haplotypes from female, male and hermaphrodite genotypes was designated as the female F haplogroup (Additional file 4); it is supposed to contain the F haplotypes of FF females, MF males and HF hermaphrodites genotypes. By difference, the alternate haplotypes found in male and hermaphrodite genotypes but not present in the F haplogroup, were considered as the M and the H haplotypes respectively (Additional file 4).
Diversity statistics for wild male/female, cultivated hermaphrodite and female haplotypes groups
Wild male haplotypes
Fu and Li's D*
Zeng et al.'s E
DH test (p-value)
Domesticated hermaphrodite haplotypes
Fu and Li's D*
Wild female haplotypes
Fu and Li's D*
Domesticated female haplotypes
Fu and Li's D*
The H haplogroups showed an intermediate diversity pattern between the M and F haplogroups, but closer to the M haplogroups (Table 2). For VSVV010, the H haplogroup presented diversity patterns quite equivalent to that of M haplogroups, except for a lower haplotype diversity.
This dataset shows that the three grapevine flower sexes, male, female and hermaphrodite, could be correctly predicted in 97% of the genotypes of our geographically representative V. vinifera sample, using a few SNPs, i.e. n. 4 to 7 of VSVV007 and n. 8 of VSVV010 (identified by black circles respectively termed 1 and 2 in Figure 4a).
Male haplogroups revealed significantly negative Tajima's D, and Fu & Li's D* values for VSVV006 and VSVV009 (Table 2). For VSVV010, the Fu and Li's D* statistics were close to the significant threshold (0.10 > p-value > 0.05). For male haplogroups (Table 2), all amplicon revealed negative E values, but only VSVV006 showed a significant excess of low-frequency variants. The DH tests detected significantly positive selection on VSVV009 and VSVV010. No other sex haplogroup showed significant signature of selection.
Fst values between combinations of the four sex haplotype groups
Vitis vinifera intraspecific comparaison
Wild males vs. wild females
Wild males vs. domesticated hermaphrodites
Domesticated hermaphrodites vs. wild females
Wild females vs. domesticated females
Vitis sp. vs Vitis vinifera sylvestris
Origin of the H allele
To better understand the origin of the H alleles, we explored the relationship between the network distances of the H macrohaplotypes from the M ones (Figure 5b) and the geographic origin, use and degree of domestication  of the cultivated grapevines. None of these characteristics revealed a clear correlation with H macrohaplotypes positions in the phylogenic network. We then tried to match the network distances with the STRUCTURE groups defined in Bacilieri et al. . This work, based on 2.096 domesticated genotypes, has revealed four main genetic groups: a) wine cultivars from Western regions, b) table grape cultivars from Eastern Mediterranean, Caucasus, Middle and Far East countries, c) wine cultivars from the Balkans and East Europe, and d) a large group of cultivars with admixed genomes. Here, ANOVA analysis revealed a weak tendency (r 2 = 0.15, p = 0.09) for the Balkan and East Europe cultivars macrohaplotypes, as compared to wine Western cultivars, to be closer to the wild M macrohaplotypes.
Similarly, although pointing to different “degrees of domestication”, the 3 groups of F macrohaplotypes defined above did not show a clear geographic or genetic structure pattern that could explain group composition.
The network position of the macrohaplotypes of the two female V. monticola, V. coignetiae and the male V. balansaeana grapevines used as outgroups, were distributed coherently to their sex phenotype: both macrohaplotypes in the F macrohaplogroups for the females, one in the F macrohaplogroups and one close to the M and H macrohaplogroup ones for the male. The closest domesticated macrohaplotypes to the V. balansaeana ones belonged to two Russian cultivars : cv. Assyl kara (chASS) and cv. Ak ouzioum tpapskii (chAKO) (Figure 5).
Sex region location in Vitis vinifera subsp. sylvestris
From a locus defined by previous works on inter-specific crosses ,, 11 genes were partially sequenced on a diverse set of male and female wild grapevines. Forty-six polymorphisms in four amplicons were found perfectly or strongly linked to flower sex, allowing to locate in V. v. subsp. sylvestris a sex locus of 151.8 kb on chromosome 2, in full agreement with the 143 kb sex region defined by Fechter et al.  on a Vitis interspecific cross. Our results corroborates the dominance of the M allele over the F allele, characteristic of a XY sex-determination model, coherently with sex segregations in controlled crosses ,,. We also confirmed that the sex locus is situated downstream of SSR marker VVIB23, while previous studies, based on a lower marker density, had placed it upstream ,.
Within the 151.8 kb sex region, the polymorphisms of the centrally located VSVV008 amplicon associated only weakly with the sex trait, with one significant SNP only and no perfect M/F association. One hypothesis to explain this pattern may be that local recombination disrupted the association pattern in V. v. sylvestris. Unfortunately, in our work we were not able to unequivocally confirm the VSVV008 position within the sex locus. Actually, the PCR primers for this amplicon were designed based on the synteny between several genome sequence assemblies: the chromosome 2 of the 8×grape genome, the putative hermaphrodite allele on the unassigned scaffold_233 (12x.0) and the male V. cinerea BAC sequencing map , where VSVV008 is located as expected between VSVV007 and VSVV009. According to this information, we expected that VSVV008 would amplify only in males; however, in our V. sylvestris sample, it amplified indifferently of sex. Even if the sequence obtained or its PCR primers did not blast anywhere else in the genome some doubts still remain about the true coordinates of VSVV008; only new specifically designed experiments may help to definitely confirm the VSVV008 position.
Characterisation of the sex locus
Over the four genes linked to sex, we found a strong LD, unprecedented in Vitis vinifera, with a mean r 2 VS of 0.72 over 109.76 kb. In V. vinifera, LD has been shown to decay rapidly: in more than 200 gene fragments, Lijavetsky et al.  observed an LD decay lower than 0.2 at around 200 bp, a finding later confirmed through massive genotyping by Myles et al. . A larger LD in the sex locus, as compared to other genome regions, could be an indication of suppression of recombination, a feature typical of heteromorphic XY-like chromosomal regions .
The lowest values of haplotype diversity (H) were found in male haplotypes of wild grape, with the predominant occurrence of one major haplotype, distributed without variation over largely diverse geographic origins, from Eastern to Western European, Caucasian and North African provenances. Hermaphrodite domesticated grapes displayed haplotype diversity values higher than that in wild males, while female haplotypes had the highest values, without notable differences between wild and domesticated pools. The large Fst values between males and females confirm the clear genetic differentiation between the M and the F haplotypes. The negative significant Tajima's D and the Fu and Li's D* values in M haplotypes of VSVV006 and VSVV009 indicate an excess of rare polymorphisms, revealing purifying selection. Using V. balanseana as outgroup, the Zeng et al. E statistics and DH test  confirmed this pattern for VSVV009. For the VSVV007 M haplotypes, these statistics were negative but not significant, probably because over its 849 bp length, we found only one segregating site. Such monomorphism may be a signal of stabilising selection, in particular because our grape samples originated from very diverse geographic regions. Indeed, in grapevine, previous works evidenced a much higher variation rate, with an average of 1 SNP in 47-129 base pairs, according to the genome region and the population studied ,,. By contrast, the F haplotypes for the four sex-linked presented no significant traces of selection, suggesting that these alleles are evolving under a neutral model.
Overall, the sex region presents traits typical of a small XY non-recombining region . According to the commonly accepted model of sex chromosome evolution in plants, such a region can appear in dioecious species when recombination suppression occurs between two closely located male- and female-sterile mutations . The F allele is expected to contain a recessive, “loss-of-function” type, male sterility mutation whereas the M allele would harbour a fully-functioning male fertility allele with, at a nearby locus, a dominant female sterility mutation . In such a case, the M allele is expected to be constrained by selection against a recombination between the two sex-determining loci, since recombination may bring either total sterility, or reversion to the ancestral hermaphrodite state. The accumulation of insertions, inversions, repeated elements and chromosomal rearrangements in the X and the Y counterparts  may add to this mechanism, impeding local chromosome pairing at meiosis. Indeed, in this locus, Fechter et al.  reported the presence of additional repeated FMO elements and of a retrotransposon in the female allele, both absent from the male allele. These structural differences may help repress local recombination between M and F alleles. The suppression of the recombination may in turn be at the origin of the linkage disequilibrium, and it may as well explain part of the reduction of diversity in M alleles.
In this region, the weaker association with the sex trait in a distal (VSVV010) and, if accurately located, a central (VSVV008) genes could be a trace of some recombination events, sufficient to break the association with the sex causal genes, but not ample enough to completely blur LD traces (Figures 2 and 3). Rare recombination events could have prevented the evolution of this small sex region into a full sex chromosome in Vitis, although dioecy is supposed to have appeared in this taxon millions of years ago . Finally, if the VSVV008 is well located in the sex locus, sex determinism in Vitis might be the result of two distinct sets of mutation in two linked gene regions, one including VSVV006 and VSVV007, and the other including VSVV009. As in Fragaria virginiana Mill. , the female and male sterile mutations could be not completely linked allowing the appearance of neuter and hermaphrodite individual. Some hermaphrodite grapevines have been already observed in natural conditions, but their wild status is still uncertain today as they may be escapees from cultivation . Similarly, in the long-lived, late-flowering and disease-prone grapes, while non-flowering plants are observed both in the wild and in experimental breeding, it is very difficult to unequivocally establish whether these are neuter or just growing in flowering-limiting conditions.
The length of the small XY region in Vitis vinifera is less than 1% of the chromosome length, much shorter than the small sex region identified in papaya which covers 10% of the chromosome . In this small sex region, the flavin-containing monooxygenase (FMO) genes and the adenine phosphoribosyl transferase (APT3) have been already suggested as good functional candidates for flower sex determination in grapevine . Other candidate genes could be mentioned such as the trehalose-6-phosphate phosphatase (TPP) that controls inflorescence architecture in maize through sugar signal modification  and its direct product, the disaccharide trehalose, has a marked effect on flowering transition . The WRKY transcription factors are one of the largest families of transcriptional regulators  and one of these factors has been shown to regulate endosperm growth and cellularization in Arabidopsis. The VSVV008 amplicon was designed in a gene which the predicted protein reveals similarity with a Ethylene Overproducer-like 1 (ETO1). The Arabidopsis ETO1 protein specifically inhibits the enzyme activity of the 1-aminocyclopropane-1-carboxylate synthase (ACS) , known to be involved in sex determination in melons (Cucumis melo) . However, for the YABBY protein, the polymorphisms did not correlate with phenotypic sex, suggesting that the association found by Battilana et al.  is the result of an extended intergenic linkage disequilibrium (LD), and not a direct indication of a causal mutation.
Origin of the H allele and traces of domestication
The last objective of this study was to elucidate the origin of hermaphroditism in domesticated grapevines. Both Fst and network analysis revealed that H haplotypes are more closely related to M than to F haplotypes. Thus, the H allele of the domesticated hermaphrodite grapevines may have derived from the M allele of wild male grapevines as suggested by previous authors ,. Interestingly, in Carica papaya, hermaphrodites are also heterozygous for a Y chromosome variant (Yh), more similar to the male-determining Y than to the X ,. However, while all combinations of Y and/or Yh are lethal in Papaya, in Vitis HH genotypes do strive and set seeds, as in the case of certain domesticated grapevines such as Chardonnay, Muscat de Hambourg, Riesling or Cardinal (Truel pers. comm., Vassal INRA), which produce 100% hermaphrodite progenies. Thus, the H allele may be an M allele having lost the dominant female sterility mutation, explaining the dominance of the M allele over the H allele. This hypothesis could also explain the increase in diversity observed in the H haplotypes as compared to the M haplotypes.
Studying phylogenetic patterns among the haplotypes, we could only found a weak tendency of the H macrohaplotypes of cultivars from Eastern regions cultivars, as compared to Western cvs, to be closer to the wild M macrohaplotypes. Former studies situated the major grapevine domestication region in the Eastern part of the Mediterranean area ,, which is thus consistent with our data.
More interestingly, the network analysis showed that both the F and M/H haplogroups are each divided in subgroups. In particular, wild female macrohaplotypes are subdivided in two main groups, one closely connected to the V. balanseana F haplotype, and the other farther away from Vitis sp. females; domesticated female haplotypes are divided in three groups, the first one close to the V. balanseana group, and the other two branching as independent lineages from the main V. sylvestris haplogroup. Similarly, in the M/H group, while the small differentiation within M haplotypes allows for less discrimination in the wild haplogroup, the cultivated hermaphrodites are again divided in groups, one including Eastern varieties and the other with a Western component. The general picture obtained with the network analysis points to a genetic structure of the wild V. vinifera haplotypes, in relation with other species, supporting the hypothesis presented in Peros et al.  that two chloroplast lineages from different Asian species (V. piasezkii, V. amurensis and V. thunbergii) contributed to the emergence of wild V. vinifera populations in Europe. On the other hand, the group differentiation in the domesticated compartment, both for the F and the H haplotypes, suggests multiple domestication events, as advanced by Arroyo et al.  based on chloroplast genetic diversity. More surprisingly, we found that the H haplotype from cv. Assyl kara, a Russian cultivar, derives directly, via a series of mutations, from V. balansaena (Figure 5). In the F group, the cultivar closest to V. balanseana is also a Russian cultivar, cv. Ak ouzioum tapapskii. Based on this evidence, we can advance the hypothesis that, in the sex region, in addition to the already known contribution from V. vinifera ssp. sylvestris, domesticated grapes enclose a genetic contribution from different Asian species. It is historically known that during the Soviet Union period, Russian agricultural researchers were active in importing genetic variability from diverse Asian regions as a source of cold or disease resistance alleles . Indeed, Venuti et al.  recently showed that the Asian Vitis amurensis was used by breeders to introgress resistance genes into cultivated grapevines. However, since in our sample the cv. Assyl kara was recorded as one of the oldest traditional cultivar from North Caucasus , the introgression of a genetic contribution from Asian species into cultivated grapes may also significantly predate early 20th century breeding activities in Russia. It could well have occurred naturally through gene flow between different interfertile Vitis species followed by selection during domestication.
The very small differentiation found between the H and M haplotypes in the sex-linked amplicons and the small number of individuals studied here makes it difficult to clarify further the domestication pathway; this issue merits without doubt further exploration, reinforcing the argument of Venuti et al.  that new prospecting and collection of wild grapes and other Vitis species in the Eastern part of the domestication range are strongly needed presently.
The phylogeny position of V. balansaeana, V. coignetiae and V. monticola grapevines in our network, as well as segregation mapping in inter-specific crosses, both support a sex locus shared by all Vitis spp. ,, suggesting that the development of heteromorphic sex chromosomes is still in the very first stage of evolution in this taxon. In general, the age of a sex-determining region can be estimated from the age of the taxon in which it is found . As in the subgenus Vitis, dioecy is the ancestral condition, its sex-determining region should be at least as old as the separation of the Vitis and Muscadinia subgenera, thought to have diverged approx. 18 My ago . Other dioecious species with a sex region of approximately the same age, such as Silene latifolia, Bryona dioica or Rumex spp. , have reached the final stages of sex chromosome evolution, with either full heteromorphic sex chromosomes or very large regions encompassing hundreds of genes. Future works to fully sequence the sex locus in a larger sample of genotypes in Vitis species could contribute to understand why some dioecious plants rapidly developed specific sex chromosomes, while others did not.
In Vitis vinifera subsp. sylvestris, we confirmed a sex locus of 151,8 kb located downstream to the marker VVIB23 and displaying haplotype diversity, linkage disequilibrium and differentiation that typically correspond to a small XY sex-determining region with XY males and XX females. This small sex-determining region, spanning less than 1% of chromosome 2 and also present in other Vitis species, suggests that grapevines could be organisms of choice to study the early stages of evolution of sex chromosomes in perennial species.
Hermaphrodite alleles appear to derive from male alleles of wild grapevines, with successive recombination events allowing import of diversity from the X into the Y chromosomal region and slowing the expansion of the region into a full heteromorphic chromosome. Macrohaplotypes network patterns are consistent with a major grapevine domestication region in the Eastern part of the Mediterranean area and secondary domestication events in geographically distinct areas. Finally, we hypothesise that in the sex region some domesticated grapes enclose a genetic contribution from different Asian species. Our findings should encourage new prospections and collection of wild grapes, including other Vitis species, in the Eastern part of the domestication range.
Availability of supporting data
The sequences data sets supporting the results of this article are available in the Genbank repository, [GenBank: KJ575622-KJ57662; http://www.ncbi.nlm.nih.gov/genbank/]".
The authors thank the staff of the INRA Vassal grapevine collection, as well as Catherine Roux, Catherine Bréton, Amandine Launay for their laboratory work, and Jean-Pierre Peros and Sandrine Maurice for their useful suggestions on an earlier version of the manuscript. This work has been supported by the CNRS ATIP project “ARCHEO-VITIS”, the ANR program “FRUCTIMEDHIS” and by the PhD fellowships from the CNRS/Languedoc-Roussillon region (France) to S. Picq.
- Levadoux L: Les populations sauvages et cultivées de Vitis vinifera L. Ann Amélioration Plantes. 1956, 6: 59-117.Google Scholar
- Zohary D, Hopf M, Weiss E: Domestication of Plants in the Old World. Oxford University Press, Oxford; 2012.View ArticleGoogle Scholar
- Galet P: Dictionnaire Encyclopédique Des Cépages. Hachette Pratique, Paris; 2000.Google Scholar
- Soejima A, Wen J: Phylogenetic analysis of the grape family (Vitaceae) based on three chloroplast markers. Am J Bot. 2006, 93: 278-287. 10.3732/ajb.93.2.278.View ArticlePubMedGoogle Scholar
- Wen J, Nie Z-L, Soejima A, Meng Y: Phylogeny of Vitaceae based on the nuclear GAI1 gene sequences. Can J Bot. 2007, 85: 731-745. 10.1139/B07-071.View ArticleGoogle Scholar
- Antcliff AJ: Inheritance of sex in Vitis. Ann Amélioration Plantes. 1980, 30: 113-122.Google Scholar
- Carbonneau A: Stérilit's mâle et femelle dans le genre Vitis. I. Modélisation de leur hérédité. Agronomie. 1983, 3: 635-644. 10.1051/agro:19830704.View ArticleGoogle Scholar
- This P, Lacombe T, Thomas MR: Historical origins and genetic diversity of wine grapes. Trends Genet. 2006, 22: 511-519. 10.1016/j.tig.2006.07.008.View ArticlePubMedGoogle Scholar
- Zohary D, Spiegel-Roy P: Beginnings of fruit growing in the old world. Science. 1975, 187: 319-327. 10.1126/science.187.4174.319.View ArticlePubMedGoogle Scholar
- Grassi F, Labra M, Imazio S, Spada A, Sgorbati S, Scienza A, Sala F: Evidence of a secondary grapevine domestication centre detected by SSR analysis. Theor Appl Genet. 2003, 107: 1315-1320. 10.1007/s00122-003-1321-1.View ArticlePubMedGoogle Scholar
- Arroyo-Garcia R, Ruiz-Garcia L, Bolling L, Ocete R, Lopez MA, Arnold C, Ergul A, Soylemezoglu G, Uzun HI, Cabello F, Ibanez J, Aradhya MK, Atanassov A, Atanassov I, Balint S, Cenis JL, Costantini L, Goris-Lavets S, Grando MS, Klein BY, McGovern PE, Merdinoglu D, Pejic I, Pelsy F, Primikirios N, Risovannaya V, Roubelakis-Angelakis KA, Snoussi H, Sotiri P, Tamhankar S, et al: Multiple origins of cultivated grapevine (Vitis vinifera L. ssp sativa) based on chloroplast DNA polymorphisms. Mol Ecol. 2006, 15: 3707-3714. 10.1111/j.1365-294X.2006.03049.x.View ArticlePubMedGoogle Scholar
- Oberle GD: A genetic study of variations in floral morphology and function in cultivated forms of Vitis. N Y State Agric Exp Stn. 1938, 250: 1-63.Google Scholar
- Levadoux L: Etude de la fleur et de la sexualité chez la vigne. Ann LÉcole Natl Agric Montp. 1946, 27: 1-89.Google Scholar
- Doazan JP, Rives M: Mise au point sur le déterminisme génétique du sexe dans le genre Vitis. Ann Amélioration Plantes. 1967, 17: 105-121.Google Scholar
- Dalbó MA, Ye GN, Weeden NF, Steinkellner H, Sefc KM, Reisch BI: A gene controlling sex in grapevines placed on a molecular marker-based genetic map. Genome. 2000, 43: 333-340. 10.1139/g99-136.View ArticlePubMedGoogle Scholar
- Riaz S, Krivanek AF, Xu K, Walker MA: Refined mapping of the Pierce's disease resistance locus, PdR1, and Sex on an extended genetic map of Vitis rupestris x V-arizonica. Theor Appl Genet. 2006, 113: 1317-1329. 10.1007/s00122-006-0385-0.View ArticlePubMedGoogle Scholar
- Marguerit E, Boury C, Manicki A, Donnart M, Butterlin G, Nemorin A, Wiedemann-Merdinoglu S, Merdinoglu D, Ollat N, Decroocq S: Genetic dissection of sex determinism, inflorescence morphology and downy mildew resistance in grapevine. Theor Appl Genet. 2009, 118: 1261-1278. 10.1007/s00122-009-0979-4.View ArticlePubMedGoogle Scholar
- Fechter I, Hausmann L, Daum M, Soerensen TR, Viehoever P, Weisshaar B, Toepfer R: Candidate genes within a 143 kb region of the flower sex locus in Vitis. Mol Genet Genomics. 2012, 287: 247-259. 10.1007/s00438-012-0674-z.View ArticlePubMedGoogle Scholar
- Vitis 12X and 8X sequences and Gaze annotations from Genoscope and CRIBI. , [https://urgi.versailles.inra.fr/Species/Vitis/Genome-Browser]
- Battilana J, Lorenzi S, Moreira FM, Moreno-Sanz P, Failla O, Emanuelli F, Grando MS: Linkage mapping and molecular diversity at the flower sex locus in wild and cultivated grapevine reveal a prominent SSR haplotype in hermaphrodite plants. Mol Biotechnol. 2013, 54: 1031-1037. 10.1007/s12033-013-9657-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Ming R, Bendahmane A, Renner SS: Sex chromosomes in land plants. Annu Rev Plant Biol. 2011, 62: 485-514. 10.1146/annurev-arplant-042110-103914.View ArticlePubMedGoogle Scholar
- Charlesworth D, Charlesworth B, Marais G: Steps in the evolution of heteromorphic sex chromosomes. Heredity. 2005, 95: 118-128. 10.1038/sj.hdy.6800697.View ArticlePubMedGoogle Scholar
- Charlesworth D: Plant sex chromosome evolution. J Exp Bot. 2013, 64: 405-420. 10.1093/jxb/ers322.View ArticlePubMedGoogle Scholar
- Laucou V, Lacombe T, Dechesne F, Siret R, Bruno JP, Dessup M, Dessup T, Ortigosa P, Parra P, Roux C, Santoni S, Vares D, Peros JP, Boursiquot JM, This P: High throughput analysis of grape genetic diversity as a tool for germplasm collection management. Theor Appl Genet. 2011, 122: 1233-1245. 10.1007/s00122-010-1527-y.View ArticlePubMedGoogle Scholar
- Peros J-P, Berger G, Portemont A, Boursiquot J-M, Lacombe T: Genetic variation and biogeography of the disjunct Vitis subg. Vitis (Vitaceae). J Biogeogr. 2011, 38: 471-486. 10.1111/j.1365-2699.2010.02410.x.View ArticleGoogle Scholar
- OIV Descriptor List for Grape Varieties and Vitis Species. Organisation Internationale de la vigne et du vin, Paris, France; 2009.
- Lowe KM, Walker MA: Genetic linkage map of the interspecific grape rootstock cross Ramsey (Vitis champinii) x Riparia Gloire (Vitis riparia). Theor Appl Genet. 2006, 112: 1582-1592. 10.1007/s00122-006-0264-8.View ArticlePubMedGoogle Scholar
- Koressaar T, Remm M: Enhancements and modifications of primer design program Primer3. Bioinformatics. 2007, 23: 1289-1291. 10.1093/bioinformatics/btm091.View ArticlePubMedGoogle Scholar
- Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, Rozen SG: Primer3—new capabilities and interfaces. Nucleic Acids Res. 2012, 40: e115-10.1093/nar/gks596.PubMed CentralView ArticlePubMedGoogle Scholar
- Staden R: The staden sequence analysis package. Mol Biotechnol. 1996, 5: 233-241. 10.1007/BF02900361.View ArticlePubMedGoogle Scholar
- Siegismund HR: Disparity in population differentiation of sex-linked and autosamal variation in sibling species of the Jaera albifrons (Isopoda) Complex. J Hered. 2002, 93: 432-439. 10.1093/jhered/93.6.432.View ArticlePubMedGoogle Scholar
- R Core Team: R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2013. ISBN 3-900051-07-0, URL., [http://www.R-project.org/]Google Scholar
- Mangin B, Siberchicot A, Nicolas S, Doligez A, This P, Cierco-Ayrolles C: Novel measures of linkage disequilibrium that correct the bias due to population structure and relatedness. Heredity. 2012, 108: 285-291. 10.1038/hdy.2011.73.PubMed CentralView ArticlePubMedGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMed CentralPubMedGoogle Scholar
- Bacilieri R, Lacombe T, Le Cunff L, Di Vecchi-Staraz M, Laucou V, Genna B, Peros J-P, This P, Boursiquot J-M: Genetic structure in cultivated grapevines is linked to geography and human selection. BMC Plant Biol. 2013, 13: 25-10.1186/1471-2229-13-25.PubMed CentralView ArticlePubMedGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.View ArticlePubMedGoogle Scholar
- Kalinowski ST, Wagner AP, Taper ML: ML-RELATE: a computer program for maximum likelihood estimation of relatedness and relationship. Mol Ecol Notes. 2006, 6: 576-579. 10.1111/j.1471-8286.2006.01256.x.View ArticleGoogle Scholar
- Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68: 978-989. 10.1086/319501.PubMed CentralView ArticlePubMedGoogle Scholar
- Stephens M, Scheet P: Accounting for decay of linkage disequilibrium in haplotype inference and missing-data imputation. Am J Hum Genet. 2005, 76: 449-462. 10.1086/428594.PubMed CentralView ArticlePubMedGoogle Scholar
- Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59: 307-321. 10.1093/sysbio/syq010.View ArticlePubMedGoogle Scholar
- Gouy M, Guindon S, Gascuel O: SeaView Version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010, 27: 221-224. 10.1093/molbev/msp259.View ArticlePubMedGoogle Scholar
- Tavaré S: Some probabilistic and statistical problems in the analysis of DNA sequences. Lect Math Life Sci. 1986, 17: 57-86.Google Scholar
- Tajima F: Statistical methods to test for nucleotide mutation hypothesis by DNA polymorphism. Genetics. 1989, 123: 585-595.PubMed CentralPubMedGoogle Scholar
- Fu Y-X, Li W-H: Statistical tests of neutrality of mutations. Genetics. 1993, 133: 693-709.PubMed CentralPubMedGoogle Scholar
- Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.View ArticlePubMedGoogle Scholar
- Zeng K, Fu Y-X, Shi S, Wu C-I: Statistical tests for detecting positive selection by utilizing high-frequency variants. Genetics. 2006, 174: 1431-1439. 10.1534/genetics.106.061432.PubMed CentralView ArticlePubMedGoogle Scholar
- Hudson RR, Slatkin M, Maddison WP: Estimation of levels of gene flow from DNA sequence data. Genetics. 1992, 132: 583-589.PubMed CentralPubMedGoogle Scholar
- Holsinger KE, Weir BS: Genetics in geographically structured populations: defining, estimating and interpreting FST. Nat Rev Genet. 2009, 10: 639-650. 10.1038/nrg2611.PubMed CentralView ArticlePubMedGoogle Scholar
- Bandelt H-J, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48. 10.1093/oxfordjournals.molbev.a026036.View ArticlePubMedGoogle Scholar
- Network software. , [http://www.fluxus-engineering.com/sharenet.htm]
- Lacombe T: Contribution à l'étude de l'histoire évolutive de la vigne cultivée (Vitis vinifera L.) par l'analyse de la diversité génétique neutre et de gènes d'intérêt. PhD thesis. Université de Montpellier 2; 2012.,Google Scholar
- Lijavetzky D, Cabezas J, Ibanez A, Rodriguez V, Martinez-Zapater J: High throughput SNP discovery and genotyping in grapevine (Vitis vinifera L.) by combining a re-sequencing approach and SNPlex technology. BMC Genomics. 2007, 8: 424-10.1186/1471-2164-8-424.PubMed CentralView ArticlePubMedGoogle Scholar
- Myles S, Boyko AR, Owens CL, Brown PJ, Grassi F, Aradhya MK, Prins B, Reynolds A, Chia J-M, Ware D, Bustamante CD, Buckler ES: Genetic structure and domestication history of the grape. Proc Natl Acad Sci U S A. 2011, 108: 3530-3535. 10.1073/pnas.1009363108.PubMed CentralView ArticlePubMedGoogle Scholar
- Salmaso M, Faes G, Segala C, Stefanini M, Salakhutdinov I, Zyprian E, Toepfer R, Stella Grando M, Velasco R: Genome diversity and gene haplotypes in the grapevine (Vitis vinifera L.), as revealed by single nucleotide polymorphisms. Mol Breed. 2004, 14: 385-395. 10.1007/s11032-004-0261-z.View ArticleGoogle Scholar
- Vezzulli S, Micheletti D, Riaz S, Pindo M, Viola R, This P, Walker MA, Troggio M, Velasco R: A SNP transferability survey within the genus Vitis. BMC Plant Biol. 2008, 8: 128-10.1186/1471-2229-8-128.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang J, Na J-K, Yu Q, Gschwend AR, Han J, Zeng F, Aryal R, VanBuren R, Murray JE, Zhang W, Navajas-Pérez R, Feltus FA, Lemke C, Tong EJ, Chen C, Man Wai C, Singh R, Wang M-L, Min XJ, Alam M, Charlesworth D, Moore PH, Jiang J, Paterson AH, Ming R: Sequencing papaya X and Yh chromosomes reveals molecular basis of incipient sex chromosome evolution. Proc Natl Acad Sci. 2012, 109: 13710-13715. 10.1073/pnas.1207833109.PubMed CentralView ArticlePubMedGoogle Scholar
- Wan Y, Schwaninger H, Baldo A, Labate J, Zhong G-Y, Simon C: A phylogenetic analysis of the grape genus (Vitis L.) reveals broad reticulation and concurrent diversification during neogene and quaternary climate change. BMC Evol Biol. 2013, 13: 141-10.1186/1471-2148-13-141.PubMed CentralView ArticlePubMedGoogle Scholar
- Spigler RB, Lewers KS, Main DS, Ashman T-L: Genetic mapping of sex determination in a wild strawberry, Fragaria virginiana, reveals earliest form of sex chromosome. Heredity. 2008, 101: 507-517. 10.1038/hdy.2008.100.View ArticlePubMedGoogle Scholar
- Arnold C, Gillet F, Gobat JM: Occurrence of the wild vine Vitis vinifera ssp. silvestris in Europe. Vitis. 1998, 37: 159-170.Google Scholar
- Liu Z, Moore PH, Ma H, Ackerman CM, Ragiba M, Yu Q, Pearl HM, Kim MS, Charlton JW, Stiles JI, Zee FT, Paterson AH, Ming R: A primitive Y chromosome in papaya marks incipient sex chromosome evolution. Nature. 2004, 427: 348-352. 10.1038/nature02228.View ArticlePubMedGoogle Scholar
- Satoh-Nagasawa N, Nagasawa N, Malcomber S, Sakai H, Jackson D: A trehalose metabolic enzyme controls inflorescence architecture in maize. Nature. 2006, 441: 227-230. 10.1038/nature04725.View ArticlePubMedGoogle Scholar
- Paul MJ, Primavesi LF, Jhurreea D, Zhang Y: Trehalose metabolism and signaling. Annu Rev Plant Biol. 2008, 59: 417-441. 10.1146/annurev.arplant.59.032607.092945.View ArticlePubMedGoogle Scholar
- Rushton PJ, Somssich IE, Ringler P, Shen QJ: WRKY transcription factors. Trends Plany Sci. 2010, 15: 247-258. 10.1016/j.tplants.2010.02.006.View ArticleGoogle Scholar
- Luo M, Dennis ES, Berger F, Peacock WJ, Chaudhury A: MINISEED3 (MINI3), a WRKY family gene, and HAIKU2 (IKU2), a leucine-rich repeat (LRR) KINASE gene, are regulators of seed size in Arabidopsis. Proc Natl Acad Sci U S A. 2005, 102: 17531-17536. 10.1073/pnas.0508418102.PubMed CentralView ArticlePubMedGoogle Scholar
- Bostick M, Lochhead SR, Honda A, Palmer S, Callis J: Related to Ubiquitin 1 and 2 are redundant and essential and regulate vegetative growth, auxin signaling, and ethylene production in Arabidopsis. Plant Cell Online. 2004, 16: 2418-2432. 10.1105/tpc.104.024943.View ArticleGoogle Scholar
- Wang KLC, Yoshida H, Lurin C, Ecker JR: Regulation of ethylene gas biosynthesis by the Arabidopsis ETO1 protein. Nature. 2004, 428: 945-950. 10.1038/nature02516.View ArticlePubMedGoogle Scholar
- Boualem A, Fergany M, Fernandez R, Troadec C, Martin A, Morin H, Sari M-A, Collin F, Flowers JM, Pitrat M, Purugganan MD, Dogimont C, Bendahmane A: A conserved mutation in an ethylene biosynthesis enzyme leads to andromonoecy in melons. Science. 2008, 321: 836-838. 10.1126/science.1159023.View ArticlePubMedGoogle Scholar
- Negi SS, Olmo HP: Sex conversion in a male vitis vinifera L. by a Kinin. Science. 1966, 152: 1624-1624. 10.1126/science.152.3729.1624.View ArticlePubMedGoogle Scholar
- Yu Q, Navajas-Pérez R, Tong E, Robertson J, Moore P, Paterson A, Ming R: Recent origin of dioecious and gynodioecious Y chromosomes in papaya. Trop Plant Biol. 2008, 1: 49-57. 10.1007/s12042-007-9005-7.View ArticleGoogle Scholar
- Pistorius R: Scientists, plants and politics: a history of the plant genetic resources movement. Rome, International Plant Genetic Resources Institute; 1997.Google Scholar
- Venuti S, Copetti D, Foria S, Falginella L, Hoffmann S, Bellin D, Cindric P, Kozma P, Scalabrin S, Morgante M, Testolin R, Di Gaspero D: Historical Introgression of the Downy Mildew Resistance Gene Rpv12 from the Asian Species Vitis amurensis into Grapevine Varieties. PLoS ONE. 2013, 8 (4): e61228-10.1371/journal.pone.0061228.PubMed CentralView ArticlePubMedGoogle Scholar
- Maghradze D, Rustioni L, Turok J, Scienza A, Failla O: Caucasus and Northern Black Sea Region Ampelography. Quedlinburg, Vitis Journal of Grapevine Research, Special Issue, Julius Kühn Institut;m 2012.Google Scholar
- Moore RC, Kozyreva O, Lebel-Hardenack S, Siroky J, Hobza R, Vyskot B, Grant SR: Genetic and functional analysis of DD44, a sex-linked gene from the dioecious plant silene latifolia, provides clues to early events in sex chromosome evolution. Genetics. 2003, 163: 321-334.PubMed CentralPubMedGoogle Scholar
- Volz SM, Renner SS: Hybridization, polyploidy, and evolutionary transitions between monoecy and dioecy in Bryonia (Cucurbitaceae). Am J Bot. 2008, 95: 1297-1306. 10.3732/ajb.0800187.View ArticlePubMedGoogle Scholar
- Navajas-Pérez R, de la Herrán R, López González G, Jamilena M, Lozano R, Ruiz Rejón C, Ruiz Rejón M, Garrido-Ramos MA: The evolution of reproductive systems and sex-determining mechanisms within rumex (Polygonaceae) inferred from nuclear and chloroplastidial sequence data. Mol Biol Evol. 2005, 22: 1929-1939. 10.1093/molbev/msi186.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.