- Research article
- Open Access
The powdery mildew resistance gene REN1 co-segregates with an NBS-LRR gene cluster in two Central Asian grapevines
© Coleman et al; licensee BioMed Central Ltd. 2009
- Received: 27 June 2009
- Accepted: 30 December 2009
- Published: 30 December 2009
Grape powdery mildew is caused by the North American native pathogen Erysiphe necator. Eurasian Vitis vinifera varieties were all believed to be susceptible. REN1 is the first resistance gene naturally found in cultivated plants of Vitis vinifera.
REN1 is present in 'Kishmish vatkana' and 'Dzhandzhal kara', two grapevines documented in Central Asia since the 1920's. These cultivars have a second-degree relationship (half sibs, grandparent-grandchild, or avuncular), and share by descent the chromosome on which the resistance allele REN1 is located. The REN1 interval was restricted to 1.4 cM using 38 SSR markers distributed across the locus and the segregation of the resistance phenotype in two progenies of collectively 461 offspring, derived from either resistant parent. The boundary markers delimit a 1.4-Mbp sequence in the PN40024 reference genome, which contains 27 genes with known functions, 2 full-length coiled-coil NBS-LRR genes, and 9 NBS-LRR pseudogenes. In the REN1 locus of PN40024, NBS genes have proliferated through a mixture of segmental duplications, tandem gene duplications, and intragenic recombination between paralogues, indicating that the REN1 locus has been inherently prone to producing genetic variation. Three SSR markers co-segregate with REN1, the outer ones confining the 908-kb array of NBS-LRR genes. Kinship and clustering analyses based on genetic distances with susceptible cultivars representative of Central Asian Vitis vinifera indicated that 'Kishmish vatkana' and 'Dzhandzhal kara' fit well into local germplasm. 'Kishmish vatkana' also has a parent-offspring relationship with the seedless table grape 'Sultanina'. In addition, the distant genetic relatedness to rootstocks, some of which are derived from North American species resistant to powdery mildew and have been used worldwide to guard against phylloxera since the late 1800's, argues against REN1 being infused into Vitis vinifera from a recent interspecific hybridisation.
The REN1 gene resides in an NBS-LRR gene cluster tightly delimited by two flanking SSR markers, which can assist in the selection of this DNA block in breeding between Vitis vinifera cultivars. The REN1 locus has multiple layers of structural complexity compared with its two closely related paralogous NBS clusters, which are located some 5 Mbp upstream and 4 Mbp downstream of the REN1 interval on the same chromosome.
- Powdery Mildew
- Powdery Mildew Resistance
- Cinnamyl Alcohol Dehydrogenase
- Table Grape
- REN1 Region
The repertoire of grapevine resistance genes against fungal diseases has recently embraced a novel member, the first naturally present in a genome of the cultivated species Vitis vinifera. Grapevines were domesticated in the Mediterranean basin, in Western and Central Asia, and then cultivated in all temperate regions around the world [2, 3]. Western Eurasian cultivars all succumb to powdery mildew outbreaks in unprotected vineyards, following the worldwide spread of the pathogen from the American continent since the second half of the 19th century . Russian breeders sought after sources of resistance in Central Asian varieties, which kindled hopes of finding resistance genes in cultivated vines (reviewed in ). 'Kishmish vatkana' is a seedless grape originally found at the beginning of the 20th century in the Uzbek oasis of Shakhrisabz, about 70 km south of Samarqand, and then maintained in local germplasm collections. 'Kishmish vatkana' is able to mount a post-penetration reaction against Erysiphe necator Schwein, a trait controlled by the dominant gene REN1. REN1 is a novel resistance gene located on chromosome 13, different from those that had been previously found in wild North American grapes and their interspecific hybrids [6–12].
REN1 or other resistance genes might be present in other accessions of V. vinifera grown in Armenia, Moldova, Russia, Georgia, and Uzbekistan (reviewed in ). Another cultivar called 'Dzhandzhal kara', also written as 'Djandjal Kara', "Karadzhandal', and 'Janjal kara' in Western literature [13, 14] and identified with the synonym 'Semiz kara', has been utilised by breeders as a source of powdery mildew resistance [14, 15]. 'Dzhandzhal kara' was first reported in the Uzbek district of Urgut, about 40 km southeast of Samarqand, and a few more stocks also existed in Kitab, within a 10-km range from Shakhrisabz where 'Kishmish vatkana' was also found. 'Dzhandzhal kara' is cultivated for the production of table grapes. 'Dzhandzhal kara' has normal hermaphrodite flowers that develop viable embryos and yield seeded fruit. 'Kishmish vatkana' is regarded as a stenospermocarpic cultivar, because its ovules degenerate soon after fertilisation at the initial stages of embryogenesis. 'Kishmish vatkana' and 'Dzhandzhal kara' survive natural powdery mildew infections better than any V. vinifera cultivar. Under high disease pressure, E. necator occasionally escapes the plant's defence reaction, which leads to microcolony formation and limited conidiation in both cultivars, with 'Dzhandzhal kara' allowing slightly more hyphal growth than 'Kishmish vatkana' (P. Kozma, unpublished results). Adaptive evolution of V. vinifera under the selective pressure of the pathogen was generally disregarded, because this crop is vegetatively propagated, and plants in the wild have probably undergone few mating generations since the arrival of the pathogen to the Eurasian continent. Thus, one might suspect a recent hybridisation with a resistant non-V. vinifera ancestor accidentally present in the area of cultivation, such as a North American rootstock or a wild grape from East Asia brought to the fertile oases of Central Asia along the ancient trade routes to the west . 'Kishmish vatkana' and 'Dzhandzhal kara' are unique within the species V. vinifera in their resistance to powdery mildew, but they are typical among domesticated varieties in all other respects, particularly when compared with varieties grown in the Uzbek oases where 'Kishmish vatkana' and 'Dzhandzhal kara' were found, such as 'Nimrang' and 'Katta kurgan'. 'Nimrang' and 'Katta kurgan' have been imported, utilised, and documented in Western and in Far Eastern literature since the middle of the last century [17–21]. Their pure V. vinifera origin has never been called into question: instead they are credited with some of the most distinguishing features of the V. vinifera subsp. sativa group orientalis subgroup antasiatica [22, 23]. Similarly, 'Kishmish vatkana' and 'Dzhandzhal kara' conform to many morphological descriptors that differentiate V. vinifera varieties from wild relatives, and thus do not appear to be the result of a recent genome introgression from feral grapevines.
If the REN1 gene had not moved into the V. vinifera germplasm from outside, the questions remain as to which functional class this gene belongs to and how it has evolved into functionality. Histological analyses have shown that REN1 does not prevent cell entry in 'Kishmish vatkana' but it inhibits the progression of pathogen invasion once haustoria have penetrated the host cells . This mode of action bears a strong resemblance to that of host resistance genes. Many plant resistances to powdery mildews are controlled by NBS-LRR genes, which all confer post penetration resistance , including Run1 in the muscadine grape . NBS-LRR genes have proliferated through reiterate tandem duplication in plant genomes, forming gene arrays prone to random shuffling of DNA segments. Genes in this class evolve rapidly, with their clustered organisation facilitating some mutated members to express novel receptors that may match the effectors of a virulent pathogen. In vegetatively propagated crops, R-genes are more likely to develop in wild populations or at the outer margins of the field crops, where plants sexually propagate and generate genetic novelty, some of which might provide a competitive advantage.
In this paper, we show (1) the presence of REN1 in a second Central Asian variety named 'Dzhandzhal kara', which has some distinctive viticultural traits with respect to 'Kishmish vatkana' (e.g. seeded vs. seedless table grape); (2) the precise location of the REN1 locus using recombinants from 461 offspring in two separate groups of progeny, one from each resistant variety crossed with susceptible grapevines; (3) gene content, structural organisation, and dynamics of evolution of the homologous 1.4-Mb REN1 region in the sequenced grapevine PN40024; (4) kinship and population analyses for uncovering the genetic background of 'Kishmish vatkana' and 'Dzhandzhal kara'.
Plant material and evaluation of powdery mildew resistance
'Kishmish vatkana' and the progeny derived from its cross with 'Nimrang' (04-10/nnn, where nnn is the identifier of each seedling) are described in . 'Dzhandzhal kara' was kindly provided by P. Cindric', University of Novi Sad, Serbia. 'Dzhandzhal kara' is one of the grandparents of the segregating population derived from the cross ('Dzhandzhal kara' × 'Lasta') × ('Katta kurgan' × 'Perlette'), coded 07-12/nnn. Lasta is a downy mildew resistant descendent of 'Villard blanc'. The 'Dzhandzhal kara' × 'Lasta' offspring named SK91-4/27, later used for further crossing, was selected at the University of Novi Sad, Serbia, with the aim of combining downy and powdery mildew resistances with table grape traits. The other grandparents are the table varieties 'Katta kurgan' and 'Perlette', a traditional Central Asian grapevine and a recently bred V. vinifera intraspecific hybrid, respectively. Except for 'Dzhandzhal kara', all grandparents are susceptible to powdery mildew, and 'Villard blanc' is the only non-V. vinifera ancestor in the pedigree of the 07-12 population. All grandparents were analysed in order to trace the origin of each allele segregating in the progeny. The 07-12 progeny were evaluated for powdery mildew resistance using the same procedure previously adopted for the 04-10 progeny .
Development and mapping of markers in the REN1 region
A search for microsatellites around the REN1 region was performed within the homologous interval of the PN40024 genome sequence using a modified version of the Sputnik repeat-finder [25, 26]. SSR primers were designed using Primer 3. The uniqueness of their annealing sites was confirmed by a BlastN search of the primer sequences against the PN40024 genome sequence. Genomic DNA was extracted from young leaves. Segregation analysis was performed as reported in . Local linkage maps were constructed using CarthaGene .
Gene annotation and structural analyses of the REN1 locus
Gene models and repetitive elements were predicted using GAZE, Repeat Masker, Tandem Repeats Finder, RepeatScout, and ReAS . GAZE gene models were annotated by a BlastX search of the NCBI protein database, and assigned to the predicted function of their best match. The alignment of RNA-Seq short reads to the genomic DNA sequence was used to find further putative coding regions outside of the GAZE models .
A representative set of NBS-LRR proteins from Arabidopsis and rice [30, 31] was used for a tBLASTn search of the REN1 region in order to uncover gene fragments and pseudogenes. The matching regions and the corresponding NBS models were numbered progressively based on their genome position. Each sequence displaying similarity supported by an E value of < e-10 was extended on each side until the next predicted gene, and manually annotated using a gene prediction consensus obtained with GenScan, FGENESH, GENEMARK, and Geneid. Amino acid alignments to other plant CC-NBS-LRR proteins guided the choice of the most likely gene model when gene finders provided conflicting predictions. The deduced amino acid sequences were aligned to reference plant CC-NBS-LRR proteins using the module Blast 2 Sequences. Protein motifs were predicted using Pfam, SMART, and COILS. Paralogous NBS genes of the NBS members found in the REN1 interval were identified by the best reciprocal match, and retained at the thresholds of E values < e-10 and significant identity over > 100 bp. Amino acid alignments were carried out using ClustalX. Similarity trees were constructed using the Neighbor Joining (NJ) algorithm implemented in MEGA, with 1000 bootstrap replicates. Comparison of duplicated regions and drawing of DNA dot plots were performed using Dotter. Global DNA alignments were done using LAGAN in a window of 100 bp with a minimum identity of 70%, and visualised with VISTA. Nucleotide positions along chr13 refer to the version of the grapevine genome assembly deposited under the accession number FN597036.
Kinship and population analyses
Genotypic data were produced using a set of 24 SSR markers, covering all chromosomes and including a standard set of 6 markers for homogenisation with the European Vitis Database . Genomic DNA was PCR-amplified using 200 μM of each dNTP, 1 U of Taq polymerase, and 2.5 pmol of each primer at an annealing temperature of 55°C with 30 PCR cycles. PCR products were run on a Beckman Coulter CEQ 8000 capillary sequencer and alleles were sized with the CEQ 8000 Genetic Analysis software, Fragment Analysis Module.
Kinship analysis was performed using KINGROUP v2 . Allele frequencies at 19 unlinked SSR loci (one per chromosome) were calculated from a sample of 16 cultivars representative of the gene pool of V. vinifera ssp. sativa convarietas orientalis and pontica, which include grapes from Central Asia and from regions around the Black Sea. Pairwise relatedness (r), which is a measure of the proportion of alleles in common that are identical by descent, and associated p-values were estimated for each possible pair of cultivars using the algorithm from . Then, a likelihood ratio (LR), with associated p-values, was calculated for each possible relationship category between the pairs, using a hypothesized relationship as the primary hypothesis, and the next closest relationship category as the null hypothesis.
Clustering was performed using NTSYS. For the 16 accessions used in the kinship analysis, genetic similarity was calculated using a simple match coefficient and clustering was performed with the unweighted pair-group method with arithmetic mean (UPGMA) and NJ algorithms. For the extended germplasm comparison, genotypic data at six markers recommended by the European Vitis Database were obtained from public databases [35, 36] and literature. Similarity was calculated based on a distance coefficient (DIST) in order to handle the lack of harmonisation of allelic sizes from different databases. Principal component analysis (PCA) was based on standardised covariance of genetic distances between 'Kishmish vatkana', 'Dzhandzhal kara', and alleles of 216 other accessions. PCA was performed using GenAlEx .
Genetic determinant of PM resistance in 'Dzhandzhal kara'
Powdery mildew resistance is transmitted by 'Dzhandzhal kara' as a Mendelian trait. The inheritance of PM resistance was investigated using an available cross population, hereafter referred to as the 07-12 population, where 'Dzhandzhal kara' is one of the grandparents, the parents being ('Dzhandzhal kara' × 'Lasta') × ('Katta kurgan' × 'Perlette'). Each of the grandparents is heterogeneously heterozygous and homozygous at a genome-wide level, while they are Rr, rr, rr, and rr for the PM resistance phenotypic trait, respectively. The parents of the 07-12 population are Rr and rr, respectively. Thus, the 07-12 segregating population is formally an intercross or cross population, according to the nomenclature given by popular mapping software such as Carthagene and JoinMap. However, this kind of population is the best approximation for a classical testcross in the genetic analysis of a heterozygous Rr trait in outcrossing species, in the absence of full homozygous lines.
The inheritance of PM resistance was studied in 151 offspring of the 07-12 population by scoring the presence or absence of visible symptoms. Of these, 71 siblings showed masses of conidiophores detectable to the naked eye, while 80 siblings showed no symptoms during the entire vegetative season. The observed segregation ratio fits a monogenic control of the trait (1:1, χ2 = 0.12), as previously observed for 'Kishmish vatkana' .
Test of independence for six markers, found to be linked to REN1 in 'Kishmish vatkana', and phenotypes of powdery mildew resistance/susceptibility in a ('Dzhandzhal kara' × 'Lasta') × ('Katta kurgan' × 'Perlette') F1 family.
χ2 test of independence
Marker saturation of the REN1 interval in 'Kishmish vatkana'
Restriction of the REN1 interval based on crossovers in 'Kishmish vatkana' and 'Dzhandzhal kara' × 'Lasta'
The haplotypes of the recombinant individuals that progressively restricted the interval down to scaffold_47 and to the 1.4 Mbp region within scaffold_47 were visually inspected and compared with their corresponding phenotype [see Additional file 2]. Genotyping of 310 offspring of 'Kishmish vatkana' revealed four recombination events that occurred most closely to REN1 in the germ cells that gave rise to the individuals 04-10/48 and 04-10/248, which both have a one-point crossover between markers s47_11 and s47_14, and in the germ cells that gave rise to the individuals 04-10/3 and 04-10/337, whose chromosome (chr) 13 tetrads recombined between sc47_20 and sc47_22, on the opposite side of REN1. The interval between sc47_11 and sc47_22 corresponds to 1.7 cM in the 'Kishmish vatkana' genetic map, and to a physical distance of 2.2 Mbp in PN40024. Despite being sampled from a smaller family of 151 offspring, the most informative recombinants in the 'Dzhandzhal kara' × 'Lasta' descendents originated from crossovers closer to REN1 on both sides of the gene. The two closest recombination events occurred in the germ cells that gave rise to the offspring 07-12/82, which bore a crossover between the markers sc47_14 and SC8_0071_014 on the upper side of the gene, and the offspring 07-12/139, which bore a crossover between sc47_18 and sc47_20 downstream of the gene. The interval between markers sc47_14 and sc47_20 is 1.4-cM long in the 'Dzhandzhal kara' × 'Lasta' genetic map, and corresponds to 1.4 Mbp in the PN40024 genome assembly.
The REN1 region in PN40024: genes and structural features
The REN1 region from marker sc47_14 to sc47_20 spans 1,440,612 bp in the homologous interval of the PN40024 genome. Within this region there are 27 genes predicted by GAZE and BlastX that putatively encompass 15 different functions, and 11 other gene models that encode full-length or partial NBS-LRR disease resistance genes with a coiled-coil domain in each N-terminus [see Additional file 3 and Additional file 4]. The hypothetical protein products of 11 other gene models had no significant match in NCBI. The average gene density is one gene every 28.8 kb. Gene distribution is quite even and complementary to the repetitive fraction, which overall covers 29% of the interval.
As for other genes with putative functions that might be physiologically relevant to disease resistance, the most abundant gene category was represented by 9 cinnamyl alcohol dehydrogenases (CADs) involved in lignin biosynthesis, two WD40-repeat proteins, a kinase-like protein, and a Rab1/RabD small GTPase [see Additional file 3 for the complete list of genes in the region]. Among these, 6 CADs, 2 DNA gyrases, a dihydropyrimidinase, a kinase-like protein, an NADH-dehydrogenase, a vacuolar protein, and a tRNAval synthetase showed evidence of constitutive expression in plantlets of PN40024, as assessed by deep RNA-Seq .
NBS paralogues in the REN1 region: local and ectopic duplications
The role of segmental duplications and/or small-scale gene duplication between and within gene clusters was assessed through three complementary approaches: dot plot analyses, global alignments of the entire paralogous segments (from 400 kb to 1.4 Mb), and percentage of conserved sequences in pairwise alignments of the 10-kb gene environment around each NBS-LRR gene.
There are 261.7-kb of conserved non-coding sequences in the chr13:11.8..12.2 Mb cluster (65% of the total length) that are shared with the REN1 interval, at an average identity of 93.9%. Physical expansion of the REN1 interval can be largely attributed to segmental local duplications of DNA blocks containing NBS and CAD genes and also by colonisation of genic and intergenic space by TEs (Figure 3). Overall, these events caused an approximate three-fold expansion of the REN1 region with respect to the paralogous locus chr13:11.8..12.2 Mb. Regions of extended colinearity are evident by dot plot comparison in the distal part of the two clusters, from immediately upstream of NBSk to the end of the REN1 interval, as well as at and downstream of the NBSc location [see Additional file 5]. Large duplications did not occur in the chr13:11.8..12.2 Mb cluster [see Additional file 5]. The relationship between the REN1 interval and the lower NBS cluster (chr13:21.5..21.9 Mb) is more limited: they share only 3.9% of non-coding sequences, but gene order and orientation of two progenitor NBS paralogues and a CAD gene in chr13:21.5..21.9 Mb are conserved, including the paralogue of NBSg, the only NBS with forward orientation in the REN1 interval [see Additional file 5]. Similar levels of conserved sequences (2.6%) and dot plot colinearity are recognisable between the chr13:11.8..12.2 Mb cluster and the chr13:21.5..21.9 Mb cluster [see Additional file 5].
Based on the level of sequence conservation in non-coding regions over the entirety of these intervals, and the nucleotide similarity in the 10-kb NBS gene environment between and within clusters and their best reciprocal matches, the most parsimonious model of evolution, which explains the present-day organisation of the REN1 locus in PN40024, is as follows. Three progenitor NBS-LRR genes, structured in a proto-cluster in an ancestral genome, generated the current genomic structure of the REN1 locus and its paralogous loci through several episodes of ectopic and local duplications (Figure 3, ancientNBS-1 to -3). A segment containing two progenitor NBS sequences (ancientNBS-3 and the oppositely oriented ancientNBS-2) and the intervening CAD, duplicated ectopically, giving rise to the common ancestor of the REN1 locus and the chr13:11.8..12.2 Mb cluster. The third progenitor NBS sequence, that is ancientNBS-1 or one of its six locally duplicated copies, duplicated and moved ectopically to the present-day location of NBS353 (chr13:7.4 Mb), where it has remained as a dispersed gene. The founder array composed of ancientNBS-1, -2, and -3 locally proliferated through tandem duplications, involving small blocks of DNA (<10 kb) around each of the three founder genes, and generated a somewhat ordered pattern of gene distribution at the present-day chr13:21.5..21.9 Mb cluster (Figure 3 on top) [see Additional file 5].
The DNA block containing an intact progenitor of NBS358, an intact paralogue descended from the oppositely oriented ancientNBS-3, a CAD, intact progenitors of NBS359 and NBS360, and distally the genes tRNAval synthetase and vacuolar-sorting-protein-family-36, were the progenitor sequences that ectopically duplicated giving rise to the chr13:11.8..12.2 Mb cluster and to the segment that subsequently evolved into the REN1 locus.
The founder NBS genes of the REN1 locus were NBSc (which is paralogous to NBS358), NBSg, NBSj (which is paralogous to NBS360), and NBSk (which is paralogous to NBS359). Comparison between NBSc, NBSj, and NBSk corroborated the hypothesis that NBSj (or its paralogue NBS 360) originated by a sequence exchange, which brought together the CC-NBS domains of NBSc (or NBS358) and the LRR of NBSk (or NBS359), far beyond the boundary of the coding region on either side. The hybrid position of NBSj and NBS360 in the NJ tree of Figure 5, between the pairs NBSc/NBS358 and NBSk/NBS359, provides additional support for NBSj and NBS360 being chimeric copies generated by intergenic recombination.
Five gene islands populated by duplicated NBS and CAD sequences were distinguishable by dot plot self-comparison of the REN1 interval (Figure 3), [see Additional file 5]. In order to understand the relative role of single-gene tandem duplications and larger duplications in generating this expansion, the percentage of conserved sequence between each pair of 10-kb sequence around each NBS-LRR (Figure 4) was compared in the context of the percentage of conserved sequence between gene islands, which ranged in size between 118 and 186 kb. Gene islands 2-5 share 63.8 kb with 93.8% nucleotide identity, which represent 42.4% of the total length. Gene islands 3-4 share 48.5 kb (40.3% of the total length) with 94.3% nucleotide identity. Gene island 1 is more similar to 2 than to any other (90.6% identity across 64.4 kb, 39% of the total length of the pairwise comparison). The founder NBS genes are localised in three of the five gene islands (NBSc-1, NBSj-4, and NBSk-5). Gene island 2 appears to have originated from a duplication of gene island 5, which might have involved as much as 63.8 kb. Within gene island 2, NBSf has the highest similarity with NBSk, located in gene island 5, and their sequence conservation extends several kilobases outside of the coding regions (Figure 4). The other two NBS genes in gene island 2 (NBSd and NBSe) appear to have arisen in this region later by small scale duplications (see below). Gene island 3 contains NBSi, which is the most similar paralogue of NBSj, a founder NBS in gene island 4 (Figure 4). Gene island 1 is the most heterogeneous, and it has been shaped by independent small-scale duplications. NBSa is the paralogue of NBSc within the same gene island. NBSb shares the highest identity with NBSj in gene island 4, but the overlapping region does not extend significantly beyond their coding regions (Figure 4). NBSb appears to have later generated NBSe by a gene duplication, which invaded gene island 2 together with a further duplication of NBSf into NBSd. NBSg is the only founder NBS that did not undergo single gene duplication and did not remain involved in larger segmental duplications.
CAD genes in the REN1 interval have proliferated in a similar way [see Additional file 6]. CAD8 and CAD9 in gene island 5 are the most similar paralogues of the two CAD genes present in the chr13:11.8..12.2 Mb cluster. Within the REN1 interval, CAD1 and CAD2, upstream of all NBS genes, originated from a duplication of CAD8 and CAD9, after the insertion of a LINE in CAD8, in addition to the LINE shared with the progenitor CAD in the chr13:11.8..12.2 Mb cluster. CAD4 and CAD7 share high identity with CAD8, while CAD 5 and 6 show partial overlap with CAD8 in their coding sequences.
Kinship and population analyses: genealogy of 'Kishmish vatkana' and 'Dzhandzhal kara'
Clustering methods applied to outcrossing diploids result in an average estimation of the two component haploid genomes present in each individual, which might obscure the contribution of outlying alleles in hybrids with a single feral parent. Utilising the parent-offspring information for 'Kishmish vatkana', the haploid genome that contributed the resistant REN1 allele could be separated from the haploid genome contributed by 'Sultanina'. This allowed the in silico construction of a di-haploid genotype which contained those alleles of 'Kishmish vatkana' that were not inherited from 'Sultanina'. The resulting virtual genotype was used for clustering with the same grape varieties as described above, in order to shed light on the phylogenetic relatedness of the resistant genotype (Figure 8). Though the virtual di-haploid did not become separated from the Central Asian group, it clustered more distantly from the Uzbek susceptible cultivars, and more tightly to grapevine varieties from Georgia and Armenia. Interestingly, this region is considered one of the centres of origin for domesticated grapes and is home to the highest genetic diversity of V. vinifera around the world. The gene pool analysed in this study encompassed 134 alleles at 24 loci. Only 3 alleles in 'Dzhandzhal kara' (6.2%) and 4 in 'Kishmish vatkana' (8.3%), including one allele at the REN1-linked locus SC8_0071_014, were not found in the reference gene pool.
In order to more widely sample the genetic variability represented in the grape germplasm, a principal component analysis was performed using the 6 markers recommended by the European Vitis Database and genotypes of 218 accessions including cultivars of V. vinifera, rootstocks, hybrids, and wild species [see Additional file 9]. Twenty three accessions are included in Cluster 1 [see Additional file 9], 21 of which are Western European cultivars, and two (Orion and Phönix) are interspecific hybrids with large parts of their genomes inherited from V. vinifera backcrosses. Cluster 2 encompasses 36 accessions: two (41B and Léon Millot) contain 50% of a non-V. vinifera genome, and seven are interspecific hybrids with most of the genome inherited from V. vinifera, while the remainder have a predominantly Western European origin. Cluster 3 is composed of Western European V. vinifera with a few hybrids intermixed. Cluster 4 encompasses the broadest heterogeneity, including North American species and hybrids used as rootstocks. 'Kishmish vatkana' and 'Dzhandzhal kara' grouped together in Cluster 5, along with other V. vinifera varieties of the proles orientalis and pontica. Cluster 6 includes only 11 varieties, including 'Katta kurgan', 'Sultanina', and 'Monukka'.
In this paper, we present evidence that the REN1 resistance allele, originally identified in 'Kishmish vatkana' , is also present in the genome of a second Central Asian grapevine cultivar named 'Dzhandzhal kara'. 'Kishmish vatkana' and 'Dzhandzhal kara' are estimated to share 30% of their genome by descent, in addition to sharing the resistant haplotype at the REN1 locus. The presence of 'Kishmish vatkana' and 'Dzhandzhal kara' was documented by Russian breeders after World War I, within a few kilometres in the Kitab-Shakhrisabz oasis, the fertile northeastern corner of the Kashkadarya River basin, south of Samarqand. 'Dzhandzhal kara' plants were also found in the district of Urgut, about 50 km to the northeast of the Kitab-Shakhrisabz oasis, on the western edges of the Pamir Mountains. While 'Kishmish vatkana' has been neglected by viticulturists and has been maintained only in germplasm collections, 'Dzhandzhal kara' has gained a wide reputation for fresh fruit production in several countries of the former Soviet Union, and it has been firmly classified as a pure V. vinifera in the international scientific literature.
By identifying two grapevines that share REN1 by descent, more information has been gained about the genealogy of these resistant grapevines. A likelihood-based kinship analysis, which predicts parent-offspring and second-degree relationships without prior knowledge of the parents , provided evidence that one of the parents of 'Kishmish vatkana' is 'Sultanina', also known as 'Thompson Seedless'. This parent-offspring relationship is not surprising, since both cultivars produce stenospermocarpic fruit, and 'Sultanina' has been widely cultivated for centuries in Uzbekistan and Northern Afghanistan, with the local synonym of 'Ak-Kishmish'. Since 'Kishmish vatkana' but not 'Sultanina' has a second degree relationship with 'Dzhandzhal kara', and the resistant and seedless 'Kishmish vatkana' displays a combination of biological features that are uncoupled in the parents (the pollen parent 'Sultanina' does not produce viable embryos and it is susceptible to powdery mildew, the unknown resistant parent must have been able to produce seeds), the direction of the relationship could also be established. It is not known whether the cross was made deliberately or occurred by chance in a promiscuous vineyard. However, the fact that 'Kishmish vatkana' is poorly documented in chronicles and historical records, in spite of its attractiveness as a seedless and powdery mildew resistant table grape, suggests that this hybridisation occurred rather recently.
'Dzhandzhal kara' has a longer recorded history and was already disseminated worldwide by the early 1970's, which might suggest that the seedling 'Dzhandzhal kara' is more ancient than 'Kishmish vatkana'. A grape accession catalogued as Vitis vinifera subsp. vinifera 'Karadzhandal' in a US collection in 1965 is currently maintained at the USDA-ARS National Germplasm Repository in Davis, California. A variety named 'Kara dzandzal' was mentioned in a trial conducted in Japan in 1969 on 74 varieties imported from the former USSR, and recommended for breeding in a strict list of five top cultivars, including the typical Central Asian table grapes 'Katta kurgan' and 'Nimrang' . A variety spelled 'Djandjal Kara' was included in a privileged set of accessions collected in Australian germplasm repositories and used in a seminal work on cultivar DNA fingerprinting in the 1990's .
Second-degree relationships could not be identified between 'Dzhandzhal kara' and any of the few other Central Asian grapes analysed in this paper, except for 'Kishmish vatkana'. However, the alleles carried by 'Dzhandzhal kara' at 19 SSR loci sampled over all chromosomes and the 'Kishmish vatkana' alleles at the same loci not shared with 'Sultanina' seem not to be foreign to the gene pool of V. vinifera ssp. sativa proles pontica and orientalis, which include domesticated cultivars from Central Asia and around the Caspian Sea. Thus, the kinship data, in addition to the phenotypic features of these plants, argue against close relatedness to exotic germplasm and the introduction of REN1 from outside of V. vinifera, such as East Asian genotypes or rootstocks of North American origin. Populations of Vitis vinifera ssp. sylvestris, the grapevines that climb trees in natural environments, extend far into Tajikistan and the Western Himalayas. The possibility of an infusion of genes, including resistance alleles if they exist, from indigenous and undomesticated grapevines into Central Asian varieties can not be excluded. The existence of continuous gene flow between natural and cultivated compartments of V. vinifera and the frequency with which these grapevines spontaneously hybridise has been recently estimated in a few Western European sites .
An inquiry further back into the genealogy of the powdery mildew-resistant genotypes could provide additional information regarding the origin of REN1. This might be accomplished by DNA analysis of more exhaustive collections of the Asian V. vinifera germplasm. Grape germplasm is very rich with indigenous and introduced accessions in the geographical region where the REN1 resistant genotypes have been found. Viticulture in fertile oases of northern Afghanistan and Uzbekistan was founded by Greek colonists in the fourth century B.C. and had already flourished by the second century B.C., at the time when General Zhang Qian brought viticulture back to China from his expeditions in the far provinces of the Western Han realm and in the modern-day Uzbekistan . The geographical location of Uzbek oases on the harsh northern routes of the Silk Road accelerated the introduction of varieties traded from commercial hubs in the eastern Mediterranean basin (such as 'Sultanina', which is thought to have originated in Iran or Turkey), varieties from the Black and Caspian Seas, and grapes collected in the Xinjiang province of modern China and the Turfan depression where V. vinifera cultivars have been cultivated for over 2,000 years [see Additional file 10]. All of this transiting material might have mixed and cross-pollinated with locally domesticated grapes or indigenous forms of V. vinifera subsp. sylvestris. In particular, Samarqand and its surroundings, where 'Kishmish vatkana' and 'Dzhandzhal kara' were discovered, was an important hub for caravans that headed south crossing the Kashkadarya valley to Termez, the gate to northern Afghanistan, and to the junctions with the southern trade routes [see Additional file 10].
The restriction of the genetic locus in the resistant genotypes and the investigation of the genes contained in the homologous interval of the PN40024 reference sequence provided clues to the possible molecular nature of the REN1 gene and to the structural organisation of the locus. Functional candidates in the interval belong to two gene categories: NBS-LRR genes and cinnamyl alcohol dehydrogenases. NBS-LRRs are by far the most abundant class of plant resistance genes, in particular in the context of host-resistance and in response to biotrophs. The presence of this gene category in the REN1 locus is in agreement with the type of histological responses observed in inoculated leaves of 'Kishmish vatkana' . The necrosis of penetrated host cells resembles the hypersensitive response (HR) triggered by a gene-for-gene type of plant-pathogen recognition. By contrast, it is not known if the presence of CAD genes in the REN1 locus is accidental or if it might have some biological significance. Cinnamyl alcohol dehydrogenase catalyzes the last step in the synthesis of the monomeric precursors of lignin. Different CAD classes are involved specifically in wood development or in responses to biotic stresses.
CAD inhibition partially suppresses papilla-associated cell wall appositions and jeopardises penetration resistance in powdery mildew-cereal interactions [44, 45]. CADs are also involved in post-penetration resistance. Intact CAD activity is necessary for expressing an epidermal cell HR in several plant-pathogen interactions [46, 47]. This observation has been explained by the role of CADs in phenylpropanoid biosynthesis products derived from cinnamyl alcohol, such as phenolic compounds and free radicals, which are integral to oxidative processes and lignin cell wall encrustations associated with HR. In grapevines, CAD transcript levels were constitutively higher in the powdery mildew resistant V. aestivalis 'Norton' than in the susceptible 'Cabernet Sauvignon' . It is noteworthy that evolution has frequently favoured physical clustering of genes belonging to widely different pathways but with related or complementary functions in defence responses (reviewed in ). The accumulation of favourable alleles of physically linked NBS-LRR and CAD genes that trigger and contribute to exerting HR, respectively, would provide an evolutionary advantage to the individual that inherited that block of DNA.
NBS-LRR and CAD genes and pseudogenes in the REN1 locus were generated through several rounds of duplication events, which appear more chaotic than duplications in the most closely related paralogous loci of the grape genome. This DNA block appears inherently prone to structural rearrangements, which makes it somewhat similar to the powdery mildew Mla resistance locus in barley , the Rps1-k locus in soybean , and Rp1 complex in maize . In addition to segmental duplications and extensive invasion of transposable elements, NBS genes in the REN1 interval also evolved through intragenic recombination. The gene NBSj displayed a sharp variation of nucleotide identity immediately downstream of the MHDV site, the 3' border of the NBS domain, when compared to its paralogues NBSc and NBSk. NBSj is highly similar to the 5' end of NBSc and to the 3' end of NBSk, but substantially dissimilar to these in the complementary parts. The level of identity with either paralogue extends homogeneously beyond each end of the coding region, indicating that NBSj is more likely the result of intragenic recombination between tandemly arrayed paralogues rather than the result of domain-specific accumulation of advantageous mutations. Unequal crossing-over and gene conversion are events of nonreciprocal exchange between similar sequences, which can generate hybrid genes with novel functions by domain swapping (reviewed in ). Intragenic recombination frequently occurs between NBS-LRR genes and promotes diversification by reassorting domains between functional paralogues and by recruiting fragments from pseudogenes that may have accumulated mutations in the absence of selective pressure [53, 54]. Either intrachromosomal gene conversions or unequal crossing-over rearranged NBS-LRR paralogues in Arabidopsis thaliana and at the maize Rp1 locus , respectively. Both events are favoured by close physical vicinity between the partners, by high sequence similarity, and by direct orientation of the gene pair , all features that occur in the REN1 locus. In some instances, it was proven that reassortment of domains by unequal crossing over generate chimeric NBS-LRR variants associated with novel race specificities (reviewed in ).
The structural organisation of the 1.4-Mbp sequence across the REN1 locus in PN40024 bears the marks of a dynamic evolution, which suggests that this genomic region has an inherent capacity to generate high DNA variation at the population level. To pinpoint the gene that confers the REN1 phenotype and determine how the structure of the resistant haplotype differs across the REN1 locus from that in the PN40024 reference sequence, cloning and sequencing of large insert DNA BAC clones from 'Kishmish vatkana' will be initiated. The functional identification of the effective REN1 gene and its sequence information will shed light on the precise DNA variation that brought about the powdery mildew resistance phenotype in V. vinifera.
We have provided evidence that the resistance gene REN1, which prevents grape powdery mildew, is shared by descent between two varieties of V. vinifera, which have been documented in Central Asia since the 1920's. In the Eurasian continent, outbreaks of powdery mildew have occurred since 1845, after the accidental introduction of the pathogen into Western Europe. The discovery of two resistant accessions that are unique in their powdery mildew resistance among the several thousands of susceptible V. vinifera varieties provides a model for studying the molecular evolution of a resistance gene, active against a recently introduced pathogen, in a host plant species previously thought to be entirely susceptible.
The analysis of the structural organisation of the REN1 locus in the reference grapevine genome sequence has provided clues to the functional categories of candidate genes present in the genetic interval, and revealed a chaotic expansion of NBS-LRR genes within the REN1 locus, compared to its most closely related paralogous loci in the grapevine genome. Tandemly duplicated paralogues in the REN1 locus and other repetitive features in the region may have created a highly unstable region prone to meiotic mispairing, which has accelerated duplications and reshuffling of domains between gene copies, as witnessed by at least one case of a chimeric NBS-LRR generated by intragenic recombination between adjacent paralogues. These events are likely to generate new functional variants, in particular when they occur within or immediately upstream of coding regions and recombine the promoter regions or functional domains originally present in separated gene copies. This scenario may have occurred in NBSj, which bears the marks of a sequence exchange that brought together the 5' promoter region, the coiled-coil domain, and the complete nucleotide-binding-site of the paralogue NBSc with the LRR region of the paralogue NBSk.
The authors thank Petar Cindric', University of Novi Sad, Serbia, for the cross 'Dzhandzhal kara' × 'Lasta', for selecting their offspring SK91-4/27 which was used as a parent of the 07-12 mapping population used in this study, and for kindly providing the original 'Dzhandzhal kara' plant used as a grandparent of the 07-12 family. We also thank Alberto Stefan, Michelle Bowe, and John Heywood, for help in bioinformatics analyses, technical assistance, and kinship analysis, respectively. This research was supported by funds from the Italian Ministry of Agriculture (Vitis Genome Analysis, VIGNA project), from the Regional Government of Friuli Venezia Giulia (Grape Breeding Project), USDA CSREES (Grant #32008-38901-193670), and the Missouri State University Graduate College.
- Hoffmann S, Di Gaspero G, Kovács L, Howard S, Kiss E, Galbács Zs, Testolin R, Kozma P: Resistance to Erisiphe necator in the grapevine 'Kishmish vatkana' is controlled by a single locus through restriction of hyphal growth. Theor Appl Genet. 2008, 116: 427-438. 10.1007/s00122-007-0680-4.View ArticlePubMedGoogle Scholar
- Arroyo-García R, Ruiz-García L, Bolling L, Ocete R, López MA, Arnold C, Ergul A, Söylemezoğlu G, Uzun HI, Cabello F, Ibáñez 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, This P, Troshin L, Malpica JM, Lefort F, Martinez-Zapater JM: 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
- This P, Lacombe T, Thomas M: Historical origin and genetic diversity of wine grapes. Trends Genet. 2006, 22: 511-519. 10.1016/j.tig.2006.07.008.View ArticlePubMedGoogle Scholar
- Pearson RC: Powdery mildew. Compendium of grape diseases. Edited by: Pearson RC, Goheen AC. 1988, St Paul: APS Press, 9-11.Google Scholar
- Vojtovic KA: Vospriimchivost sortov vinograda k oidiumu (Powdery mildew susceptibility of grapevine cultivars). Novüje kompleksno - ustojchevüje stolovüje sorta vinograda i metodü ih polichenija. (New complex resistant table grape cultivars and methods for breeding). 1987, Kishinev: Kisinev Kartja Moldovenjaske, 42-46.Google Scholar
- Dalbó MA, Ye GN, Weeden NF, Wilcox WF, Reisch BI: Marker-assisted selection for powdery mildew resistance in grapes. J Amer Soc Hort Sci. 2001, 26: 83-89.Google Scholar
- Pauquet J, Bouquet A, This P, Adam-Blondon A-F: Establishment of a local map of AFLP markers around the powdery mildew resistance gene Run1 in grapevine and assessment of their usefulness for marker assisted selection. Theor Appl Genet. 2001, 103: 1201-1210. 10.1007/s001220100664.View ArticleGoogle Scholar
- Donald TM, Pellerone F, Adam-Blondon AF, Bouquet A, Thomas MR, Dry IB: Identification of resistance gene analogs linked to a powdery mildew resistance locus in grapevine. Theor Appl Genet. 2002, 104: 610-618. 10.1007/s00122-001-0768-1.View ArticlePubMedGoogle Scholar
- Barker CL, Donald T, Pauquet J, Ratnaparkhe MB, Bouquet A, Adam-Blondon AF, Thomas MR, Dry I: Genetic and physical mapping of the grapevine powdery mildew resistance gene, Run1, using a bacterial artificial chromosome library. Theor Appl Genet. 2005, 111: 370-377. 10.1007/s00122-005-2030-8.View ArticlePubMedGoogle Scholar
- Fischer BM, Salakhutdinov I, Akkurt M, Eibach R, Edwards KJ, Töpfer R, Zyprian EM: Quantitative trait locus analysis of fungal disease resistance factors on a molecular map of grapevine. Theor Appl Genet. 2004, 108: 501-515. 10.1007/s00122-003-1445-3.View ArticlePubMedGoogle Scholar
- Akkurt M, Welter L, Maul E, Töpfer R, Zyprian E: Development of SCAR markers linked to powdery mildew (Uncinula necator) resistance in grapevine (Vitis vinifera L. and Vitis sp.). Mol Breed. 2006, 19: 103-111. 10.1007/s11032-006-9047-9.View ArticleGoogle Scholar
- Welter LJ, Göktürk-Baydar N, Akkurt M, Maul E, Eibach R, Töpfer R, Zyprian EM: Genetic mapping and localization of quantitative trait loci affecting fungal disease resistance and leaf morphology in grapevine (Vitis vinifera L). Mol Breed. 2007, 20: 359-374. 10.1007/s11032-007-9097-7.View ArticleGoogle Scholar
- Thomas MR, Cain P, Scott NS: DNA typing of grapevines: A universal methodology and database for describing cultivars and evaluating genetic relatedness. Plant Mol Biol. 1994, 25: 939-949. 10.1007/BF00014668.View ArticlePubMedGoogle Scholar
- Korbuly J: Results of breeding for resistance to winter frosts and different pathogens using Vitis amurensis. Acta Hort. 2000, 528: 551-557.View ArticleGoogle Scholar
- Filippenko IM, Stin LT: Szort evropejskogo vida Vitis vinifera L. Dzhandzhal kara usztojcsiv k oidiumu. (Vitis vinifera L. 'Dzhandzhal kara' is resistant to powdery mildew). Bull Nauchn Inform CGL IV Michurina. 1977, 25: 57-58.Google Scholar
- Wan Y, Schwaniniger H, He P, Wang Y: Comparison of resistance to powdery mildew and downy mildew in Chinese wild grapes. Vitis. 2007, 46: 132-136.Google Scholar
- Winkler AJ, Kliewer WM, Lider LA, Cook JA: General Viticulture. 1974, Berkeley: University of California PressGoogle Scholar
- Kurihara A, Kishi T, Yamane H, Nagata K: Studies on the lately introduced grape varieties from U.S.S.R. Bull Fruit Tree Res Stn E. 1976, 1: 1-9.Google Scholar
- Heuertz M, Goryslavets S, Hausman J-F, Risovanna V: Characterization of grapevine accessions from Ukraine using microsatellite markers. Am J Enol Vitic. 2008, 59: 169-178.Google Scholar
- Le Cunff L, Fournier-Level A, Laucou V, Vezzulli S, Lacombe T, Adam-Blondon AF, Boursiquot JM, This P: Construction of nested genetic core collections to optimize the exploitation of natural diversity in Vitis vinifera L. subsp. sativa. BMC Plant Biol. 2008, 8: 31-10.1186/1471-2229-8-31.PubMed CentralView ArticlePubMedGoogle Scholar
- Obón de Castro C, Rivera D, Carreño E, Alcaraz F, Palazón JA: Seed morphology of Vitis vinifera and its relationship to ecogeographical groups and chlorotypes. Acta Hort. 2008, 799: 51-59.View ArticleGoogle Scholar
- Negrul AM: The genetic basis of grape breeding. Bull Appl Bot Leningrad Ser VIII. 1936, 6: 1-149.Google Scholar
- Troshin L, Nedov P, Litvak A, Guzun N: Improvement of Vitis vinifera sativa DC. taxonomy. Vitis, Proceedings 5th Intl Symp Grape Breeding. 1990, 37-43. SpecialGoogle Scholar
- Halterman D, Zhou F, Wei F, Wise RP, Schulze-Lefert P: The MLA6 coiled-coil, NBS-LRR protein confers AvrMla6-dependent resistance specificity to Blumeria graminis f. sp. hordei in barley and wheat. Plant J. 2001, 25: 335-348. 10.1046/j.1365-313x.2001.00982.x.View ArticlePubMedGoogle Scholar
- Sputnik. . [http://espressosoftware.com/sputnik/index.html]
- Morgante M, Hanafey M, Powell W: Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nat Genet. 2002, 30: 194-200. 10.1038/ng822.View ArticlePubMedGoogle Scholar
- de Givry S, Bouchez M, Chabrier P, Milan D, Schiex T: CarthaGene: multipopulation integrated genetic and radiation hybrid mapping. Bioinformatics. 2005, 21: 1703-1704. 10.1093/bioinformatics/bti222.View ArticlePubMedGoogle Scholar
- Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, Vezzi A, Legeai F, Hugueney P, Dasilva C, Horner D, Mica E, Jublot D, Poulain J, Bruyère C, Billault A, Segurens B, Gouyvenoux M, Ugarte E, Cattonaro F, Anthouard V, Vico V, Del Fabbro C, Alaux M, Di Gaspero G, Dumas V, Felice N, Paillard S, Juman I, Moroldo M, Scalabrin S, Canaguier A, Le Clainche I, Malacrida G, Durand E, Pesole G, Laucou V, Chatelet P, Merdinoglu D, Delledonne M, Pezzotti M, Lecharny A, Scarpelli C, Artiguenave F, Pè ME, Valle G, Morgante M, Caboche M, Adam-Blondon AF, Weissenbach J, Quétier F, Wincker P: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449: 463-468. 10.1038/nature06148.View ArticlePubMedGoogle Scholar
- Denoeud F, Aury JM, Da Silva C, Noel B, Rogier O, Delledonne M, Morgante M, Valle G, Wincker P, Scarpelli C, Jaillon O, Artiguenave F: Annotating genomes with massive-scale RNA sequencing. Genome Biol. 2008, 9: R175-10.1186/gb-2008-9-12-r175.PubMed CentralView ArticlePubMedGoogle Scholar
- Meyers B, Kozik A, Griego A, Hanhui K, Michelmore RW: Genome-wide analysis of NBS-LRR-encoding genes in Arabidopsis. Plant Cell. 2003, 15: 809-834. 10.1105/tpc.009308.PubMed CentralView ArticlePubMedGoogle Scholar
- Monosi B, Wisser RJ, Pennill L, Hulbert SH: Full-genome analysis of resistance gene homologues in rice. Theor Appl Genet. 2004, 109: 1423-1447. 10.1007/s00122-004-1758-x.View ArticleGoogle Scholar
- European Vitis Database. . [http://www.genres.de/eccdb/vitis/]
- Konovalov DA, Manning C, Henshaw MT: KINGROUP: a program for pedigree relationship reconstruction and kin group assignments using genetic markers. Mol Ecol Notes. 2004, 4: 779-782. 10.1111/j.1471-8286.2004.00796.x.View ArticleGoogle Scholar
- Wang J: An estimator of pairwise relatedness using molecular markers. Genetics. 2002, 160: 1203-1215.PubMed CentralPubMedGoogle Scholar
- Grape Microsatellite Collection. . [http://meteo.iasma.it/genetica/gmc.html]
- Swiss Vitis Microsatellite Database. . [http://www1.unine.ch/svmd/]
- Peakall R, Smouse PE: GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006, 6: 288-295. 10.1111/j.1471-8286.2005.01155.x.View ArticleGoogle Scholar
- Cipriani G, Di Gaspero G, Canaguier A, Jusseaumes J, Tassin J, Lemainque A, Roux C, Adam-Blondon A-F, Testolin R: Molecular linkage maps: Strategies, resources and achievements. Grapes. Edited by: Zapater MM, Adam-Blondon AF. 2009, Enfield: Sciences Publishers; [Chittaranjan Kole (Series Editor): Genomics of fruit and vegetables crops]Google Scholar
- Blouin MS: DNA-based methods for pedigree reconstruction and kinship analysis in natural populations. Trends Ecol Evol. 2003, 18: 503-511. 10.1016/S0169-5347(03)00225-8.View ArticleGoogle Scholar
- Vitis International Variety Catalogue. . [http://www.vivc.bafz.de]
- Vouillamoz JF, Grando MS: Genealogy of wine grape cultivars: 'Pinot' is related to 'Syrah'. Heredity. 2006, 97: 102-110. 10.1038/sj.hdy.6800842.View ArticlePubMedGoogle Scholar
- Di Vecchi-Staraz M, Laucou V, Bruno G, Lacombe T, Gerber S, Bourse T, Boselli M, This P: Low level of pollen-mediated gene flow from cultivated to wild grapevine: consequences for the evolution of the endangered subspecies Vitis vinifera L. subsp. silvestris. J Hered. 2009, 100: 66-75. 10.1093/jhered/esn084.View ArticlePubMedGoogle Scholar
- McGovern PE: Ancient wine: The search for the origins of viniculture. 2003, Princeton: Princeton University PressGoogle Scholar
- Barakat A, Bagniewska-Zadworna A, Choi A, Plakkat U, Diloreto DS, Yellanki P, Carlson JE: The cinnamyl alcohol dehydrogenase gene family in Populus: phylogeny, organization, and expression. BMC Plant Biol. 2009, 9: 26-10.1186/1471-2229-9-26.PubMed CentralView ArticlePubMedGoogle Scholar
- Bhuiyan NM, Selvaraj G, Wei Y, King J: Gene expression profiling and silencing reveal that monolignol biosynthesis plays a critical role in penetration defence in wheat against powdery mildew invasion. J Exp Bot. 2009, 60: 509-521. 10.1093/jxb/ern290.PubMed CentralView ArticlePubMedGoogle Scholar
- Zeyen RJ, Bushnell WR, Carver TLW, Robbins MP, Clark TA, Boyles DA, Vance CP: Inhibiting phenylalanine ammonia-lyase and cinnamyl-alcohol dehydrogenase suppresses Mla1 (HR) but not Mlo5 (non-HR) barley powdery mildew resistances. Physiol Mol Plant Pathol. 1995, 47: 119-140. 10.1006/pmpp.1995.1047.View ArticleGoogle Scholar
- Prats E, Martínez F, Rojas-Molina MM, Rubiales D: Differential effects of phenylalanine ammonia lyase, cinnamyl alcohol dehydrogenase, and energetic metabolism inhibition on resistance of appropriate host and nonhost cereal-rust interactions. Phytopathology. 2007, 97: 1578-1583. 10.1094/PHYTO-97-12-1578.View ArticlePubMedGoogle Scholar
- Fung R, Qiu W, Su Y, Schachtman D, Huppert K, Fekete C, Kovács L: Gene expression variation in grapevine species Vitis vinifera L. and Vitis aestivalis Michx. Genet Res Crop Evol. 2007, 54: 1541-1553. 10.1007/s10722-006-9146-9.View ArticleGoogle Scholar
- Wei F, Wing RA, Wise RP: Genome dynamics and evolution of the Mla (powdery mildew) resistance locus in barley. Plant Cell. 2002, 14: 1903-1917. 10.1105/tpc.002238.PubMed CentralView ArticlePubMedGoogle Scholar
- Gao H, Bhattacharyya MK: The soybean-Phytophthora resistance locus Rps1-k encompasses coiled coil-nucleotide binding-leucine rich repeat-like genes and repetitive sequences. BMC Plant Biol. 2008, 8: 29-10.1186/1471-2229-8-29.PubMed CentralView ArticlePubMedGoogle Scholar
- Sun Q, Collins NC, Ayliffe M, Smith SM, Drake J, Pryor T, Hulbert SH: Recombination between paralogues at the Rp1 rust resistance locus in maize. Genetics. 2001, 158: 423-438.PubMed CentralPubMedGoogle Scholar
- Chen JM, Cooper DN, Chuzhanova N, Férec C, Patrinos GP: Gene conversion: mechanisms, evolution and human disease. Nat Rev Genet. 2007, 8: 762-775. 10.1038/nrg2193.View ArticlePubMedGoogle Scholar
- Mondragon-Palomino M, Gaut BS: Gene conversion and the evolution of three leucine-rich repeat gene families in Arabidopsis thaliana. Mol Biol Evol. 2005, 22: 2444-2456. 10.1093/molbev/msi241.View ArticlePubMedGoogle Scholar
- Kuang H, Caldwell KS, Meyers BC, Michelmore RW: Frequent sequence exchanges between homologs of RPP8 in Arabidopsis are not necessarily associated with genomic proximity. Plant J. 2008, 54: 69-80. 10.1111/j.1365-313X.2008.03408.x.View ArticlePubMedGoogle Scholar
- Xu S, Clark T, Zheng H, Vang S, Li R, Wong GK-S, Wang J, Zheng X: Gene conversion in the rice genome. BMC Genomics. 2008, 9: 93-10.1186/1471-2164-9-93.PubMed CentralView ArticlePubMedGoogle Scholar
- Cipriani G, Marrazzo MT, Di Gaspero G, Pfeiffer A, Morgante M, Testolin R: A set of microsatellite markers with long core repeat optimized for grape (Vitis spp.) genotyping. BMC Plant Biol. 2008, 8: 127-10.1186/1471-2229-8-127.PubMed CentralView ArticlePubMedGoogle Scholar
- Sefc KM, Regner F, Turetschek E, Glössl J, Steinkellner H: Identification of microsatellite sequences in Vitis riparia and their applicability for genotyping of different Vitis species. Genome. 1999, 42: 367-373. 10.1139/gen-42-3-367.View ArticlePubMedGoogle Scholar
- Bowers JE, Dangl GS, Vignani R, Meredith CP: Isolation and characterization of new polymorphic simple sequence repeat loci in grape (Vitis vinifera L.). Genome. 1996, 39: 628-633. 10.1139/g96-080.View ArticlePubMedGoogle Scholar
- Bowers JE, Dangl GS, Meredith CP: Development and characterization of additional microsatellite DNA markers for grape. Am J Enol Vitic. 1999, 50: 243-246.Google Scholar
- Thomas MR, Scott NS: Microsatellite repeats in grapevine reveal DNA polymorphisms when analysed as sequence-tagged sites (STSs). Theor Appl Genet. 1993, 86: 985-990.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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.