A Genome-Wide Hybrid Incompatibility Landscape between Caenorhabditis briggsae and C. nigoni

Systematic characterization of ẖybrid incompatibility (HI) between related species remains the key to understanding speciation. The genetic basis of HI has been intensively studied in Drosophila species, but remains largely unknown in other species, including nematodes, which is mainly due to the lack of a sister species with which C. elegans can mate and produce viable progeny. The recent discovery of a C. briggsae sister species, C. nigoni, has opened up the possibility of dissecting the genetic basis of HI in nematode species. However, the paucity of dominant and visible marker prevents the efficient mapping of HI loci between the two species. To elucidate the genetic basis of speciation in nematode species, we first generated 96 chromosomally integrated GFP markers in the C. briggsae genome and mapped them into the defined locations by PCR and Next-Generation Sequencing (NGS). Aided by the marker, we backcrossed the GFP-associated C. briggsae genomic fragments into C. nigoni for at least 15 generations and produced 111 independent introgressions. The introgression fragments cover most of the C. briggsae genome. We finally dissected the patterns of HI by scoring the embryonic lethality, larval arrest, sex ratio and male sterility for each introgression line, through which we identified pervasive HI loci and produced a genome-wide landscape of HI between the two nematode species, the first of its type for any non-Drosophila species. The HI data not only provided insights into the genetic basis of speciation, but also established a framework for the possible cloning of HI loci between the two nematode species. Furthermore, the data on hybrids confirmed Haldane’s rule and suggested the presence of a large X effect in terms of fertility between the two species. Importantly, this work opens a new avenue for studying speciation genetics between nematode species and allows parallel comparison of the HI with that in Drosophila and other species.     

Genetic variation for postzygotic reproductive isolation between Caenorhabditis briggsae and Caenorhabditis sp. 9

The process of speciation is key to the origins of biodiversity, and yet the Caenorhabditis nematode model system has contributed little to this topic. Genetic studies of speciation in the genus are now feasible, owing to crosses between the recently discovered Caenorhabditis sp. 9 and the well-known C. briggsae producing fertile F(1) hybrid females. We dissected patterns of postzygotic reproductive isolation between these species by crossing eight isogenic strains of C. briggsae reciprocally with six strains of C. sp. 9. We determined that overall patterns of reproductive isolation are robust across these genetic backgrounds. However, we also quantified significant heritable variation within each species for interspecific hybrid incompatibilities for total adult progeny, egg-to-adult viability, and the percentage of male progeny. This demonstrates that intraspecific variation for interspecific hybrid incompatibility occurs despite extensive, albeit incomplete, reproductive isolation. Therefore, this emerging general phenomenon of variable reproductive isolation is not restricted to highly interfertile, early-stage incipient species, but also applies to species in the latest stages of the speciation process. Furthermore, we confirm Haldane's rule and demonstrate strongly asymmetric parent-of-origin effects (Darwin's corollary) that consistently manifest more extremely when hermaphroditic C. briggsae serves as maternal parent. These findings highlight Caenorhabditis as an emerging system for understanding the genetics of general patterns of reproductive isolation.     

Heterogeneity and concordance in locus‐specific differentiation and introgression between species of towhees

The maintenance or breakdown of reproductive isolation is an observable outcome of secondary contact between species. In cases where hybrids beyond the F1 are formed, the representation of each species' ancestry can vary dramatically among genomic regions. This genomic heterogeneity in ancestry and introgression can offer insight into evolutionary processes, particularly if introgression is compared in multiple hybrid zones. Similarly, considerable heterogeneity exists across the genome in the extent to which populations and species have diverged, reflecting the combined effects of different evolutionary processes on genetic variation. We studied hybridization across two hybrid zones of two phenotypically well‐differentiated bird species in Mexico (Pipilo maculatus and P. ocai), to investigate genomic heterogeneity in differentiation and introgression. Using genotyping‐by‐sequencing (GBS) and hierarchical Bayesian models, we genotyped 460 birds at over 41 000 single nucleotide polymorphism (SNP) loci. We identified loci exhibiting extreme introgression relative to the genome‐wide expectation using a Bayesian genomic cline model. We also estimated locus‐specific F_ST and identified loci with exceptionally high genetic divergence between the parental species. We found some concordance of locus‐specific introgression in the two independent hybrid zones (6–20% of extreme loci shared across zones), reflecting areas of the genome that experience similar gene flow when the species interact. Additionally, heterogeneity in introgression and divergence across the genome revealed another subset of loci under the influence of locally specific factors. These results are consistent with a history in which reproductive isolation has been influenced by a common set of loci in both hybrid zones, but where local environmental and stochastic factors also lead to genomic differentiation.     

Caenorhabditis briggsae recombinant inbred line genotypes reveal inter-strain incompatibility and the evolution of recombination

The nematode Caenorhabditis briggsae is an emerging model organism that allows evolutionary comparisons with C. elegans and exploration of its own unique biological attributes. To produce a high-resolution C. briggsae recombination map, recombinant inbred lines were generated from reciprocal crosses between two strains and genotyped at over 1,000 loci. A second set of recombinant inbred lines involving a third strain was also genotyped at lower resolution. The resulting recombination maps exhibit discrete domains of high and low recombination, as in C. elegans, indicating these are a general feature of Caenorhabditis species. The proportion of a chromosome's physical size occupied by the central, low-recombination domain is highly correlated between species. However, the C. briggsae intra-species comparison reveals striking variation in the distribution of recombination between domains. Hybrid lines made with the more divergent pair of strains also exhibit pervasive marker transmission ratio distortion, evidence of selection acting on hybrid genotypes. The strongest effect, on chromosome III, is explained by a developmental delay phenotype exhibited by some hybrid F2 animals. In addition, on chromosomes IV and V, cross direction-specific biases towards one parental genotype suggest the existence of cytonuclear epistatic interactions. These interactions are discussed in relation to surprising mitochondrial genome polymorphism in C. briggsae, evidence that the two strains diverged in allopatry, the potential for local adaptation, and the evolution of Dobzhansky-Muller incompatibilities. The genetic and genomic resources resulting from this work will support future efforts to understand inter-strain divergence as well as facilitate studies of gene function, natural variation, and the evolution of recombination in Caenorhabditis nematodes.     

Isolation and mapping of DNA probes within the linkage group I gene cluster of Caenorhabditis elegans

We have isolated probes for DNA polymorphisms across the linkage group I gene cluster in Caenorhabditis elegans, using Tc1-linkage selection. The probes detect strain polymorphism between the wild-type strains of var. Bristol and var. Bergerac. As a result of mapping the sites hP4, hP5, hP6, hP7, hP9, and sPl, more than 1000 kilobases (kb) of cloned cosmid DNA has been positioned on the genetic map. We found there is more DNA per map unit in the center of the gene cluster than expected on the basis of the genomic average. Furthermore, the amount is not constant across the entire region but reaches a peak in the hP9 unc-13 interval. To find the coding regions, we examined DNA cross-homology between two species, Caenorhabditis elegans and Caenorhabditis briggsae. Approximately one-third of the DNA in the hP5 hP9 interval was examined for coding regions and 21 sequences were identified within 318 kb of DNA.     

Differential patterns of introgression across the X chromosome in a hybrid zone between two species of house mice

A complete understanding of the speciation process requires the identification of genomic regions and genes that confer reproductive barriers between species. Empirical and theoretical research has revealed two important patterns in the evolution of reproductive isolation in animals: isolation typically arises as a result of disrupted epistatic interactions between multiple loci and these disruptions map disproportionately to the X chromosome. These patterns suggest that a targeted examination of natural gene flow between closely related species at X-linked markers with known positions would provide insight into the genetic basis of speciation. We take advantage of the existence of genomic data and a well-documented European zone of hybridization between two species of house mice, Mus domesticus and M. musculus, to conduct such a survey. We evaluate patterns of introgression across the hybrid zone for 13 diagnostic X-linked loci with known chromosomal positions using a maximum likelihood model. Interlocus comparisons clearly identify one locus with reduced introgression across the center of the hybrid zone, pinpointing a candidate region for reproductive isolation. Results also reveal one locus with high frequencies of M. domesticus alleles in populations on the M. musculus side of the zone, suggesting the possibility that positive selection may act to drive the spread of alleles from one species on to the genomic background of the other species. Finally, cline width and cline center are strongly positively correlated across the X chromosome, indicating that gene flow of the X chromosome may be asymmetrical. This study highlights the utility of natural populations of hybrids for mapping speciation genes and suggests that the middle of the X chromosome may be important for reproductive isolation between species of house mice.     

Bayesian estimation of genomic clines

We developed a Bayesian genomic cline model to study the genetic architecture of adaptive divergence and reproductive isolation between hybridizing lineages. This model quantifies locus‐specific patterns of introgression with two cline parameters that describe the probability of locus‐specific ancestry as a function of genome‐wide admixture. ‘Outlier’ loci with extreme patterns of introgression relative to most of the genome can be identified. These loci are potentially associated with adaptive divergence or reproductive isolation. We simulated genetic data for admixed populations that included neutral introgression, as well as loci that were subject to directional, epistatic or underdominant selection, and analysed these data using the Bayesian genomic cline model. Under many demographic conditions, underdominance or directional selection had detectable and predictable effects on cline parameters, and ‘outlier’ loci were greatly enriched for genetic regions affected by selection. We also analysed previously published genetic data from two transects through a hybrid zone between Mus domesticus and M. musculus. We found considerable variation in rates of introgression across the genome and particularly low rates of introgression for two X‐linked markers. There were similarities and differences in patterns of introgression between the two transects, which likely reflects a combination of stochastic variability because of genetic drift and geographic variation in the genetic architecture of reproductive isolation. By providing a robust framework to quantify and compare patterns of introgression among genetic regions and populations, the Bayesian genomic cline model will advance our understanding of the genetics of reproductive isolation and the speciation process.     

Clustered organization of reproductive genes in the C. elegans genome

Defining the forces that sculpt genome organization is fundamental for understanding the origin, persistence, and diversification of species. The genomic sequences of the nematodes Caenorhabditis elegans and Caenorhabditis briggsae provide an excellent opportunity to explore the dynamics of chromosome evolution. Extensive chromosomal rearrangement has accompanied divergence from their common ancestor, an event occurring roughly 100 million years ago (Mya); yet, morphologically, these species are nearly indistinguishable and both reproduce primarily by self-fertilization. Here, we show that genes expressed during spermatogenesis (sperm genes) are nonrandomly distributed across the C. elegans genome into three large clusters located on two autosomes. In addition to sperm genes, these chromosomal regions are enriched for genes involved in the hermaphrodite sperm/oocyte switch and in the reception of sperm signals that control fertilization. Most loci are present in single copy, suggesting that cluster formation is largely due to gene aggregation and not to tandem duplication. Comparative mapping indicates that the C. briggsae genome differs dramatically from the C. elegans genome in clustering. Because clustered genes have a direct role in reproduction and thus fitness, their aggregated pattern might have been shaped by natural selection, perhaps as hermaphroditism evolved.     

Patterns of nucleotide polymorphism distinguish temperate and tropical wild isolates of Caenorhabditis briggsae

Caenorhabditis briggsae provides a natural comparison species for the model nematode C. elegans, given their similar morphology, life history, and hermaphroditic mode of reproduction. Despite C. briggsae boasting a published genome sequence and establishing Caenorhabditis as a model genus for genetics and development, little is known about genetic variation across the geographic range of this species. In this study, we greatly expand the collection of natural isolates and characterize patterns of nucleotide variation for six loci in 63 strains from three continents. The pattern of polymorphisms reveals differentiation between C. briggsae strains found in temperate localities in the northern hemisphere from those sampled near the Tropic of Cancer, with diversity within the tropical region comparable to what is found for C. elegans in Europe. As in C. elegans, linkage disequilibrium is pervasive, although recombination is evident among some variant sites, indicating that outcrossing has occurred at a low rate in the history of the sample. In contrast to C. elegans, temperate regions harbor extremely little variation, perhaps reflecting colonization and recent expansion of C. briggsae into northern latitudes. We discuss these findings in relation to their implications for selection, demographic history, and the persistence of self-fertilization.     

Smoke without fire: most reported cases of intron gain in nematodes instead reflect intron losses

Identification of recently gained spliceosomal introns would provide crucial evidence in the continuing debate concerning the age and evolutionary significance of introns. A previously published genomic analysis reported to have identified 122 introns that had been gained since the divergence of the nematodes Caenorhabidits elegans and Caenorhabditis briggsae approximately 100 MYA. However, using newly available genomic sequence from additional Caenorhabditis species, we show that 74% (60/81) of the reported gains in C. elegans are present in a C. briggsae relative. This pattern indicates that these introns represent losses in C. briggsae, not gains in C. elegans. In addition, 61% (25/41) of the reported gains in C. briggsae are present in the more distant C. briggsae relative, in a pattern suggesting that additional reported gains in C. elegans and/or C. briggsae may in fact represent unrecognized losses. These results underscore the dominance of intron loss over intron gain in recent eukaryotic evolution, the pitfalls associated with parsimony in inferring intron gains, and the importance of genomic sequencing of clusters of closely related species for drawing accurate inferences about genome evolution.     

Molecular evolution and quantitative variation for chemosensory behaviour in the nematode genus Caenorhabditis

Caenorhabditis elegans is a model organism in biology, yet despite the tremendous information generated from genetic, genomic and functional analyses, C. elegans has rarely been used to address questions in ecological genetics. Here, we analyse genetic variation for chemosensory behaviour, an ecologically important trait that is also genetically well characterized, at both the phenotypic and molecular levels within three species of the genus Caenorhabditis. We show that the G-protein ODR-3 plays an important role in chemosensory avoidance behaviour and identify orthologues of odr-3 in C. briggsae and C. remanei. Both quantitative genetic analysis of chemosensory behaviour and molecular population genetic analysis of odr-3 show that there is little genetic variation among a worldwide collection of isolates of the primarily selfing C. elegans, whereas there is substantially more variation within a single population of the outcrossing C. remanei. Although there are a large number of substitutions at silent sites within odr-3 among the three species, molecular evolution at the protein level is extremely conserved, suggesting that odr-3 plays an important role in cell signalling during chemosensation and/or neuronal cilia development in C. remanei and in C. briggsae as it does in C. elegans. Our results suggest that C. remanei may be a more suitable subject for ecological and evolutionary genetic studies than C. elegans.     

Polymorphism and evolution of vulval precursor cell lineages within two nematode genera, Caenorhabditis and Oscheius

BACKGROUND: The cell lineage of nematodes is mostly invariant for a given species, but varies between species. One can thus wonder how a cell lineage varies during evolution. We have started a microevolutionary approach within two genera by observing lineage variations of vulval precursor cells in different natural nematode populations of the same and closely related species. RESULTS: In Caenorhabditis elegans, the P3.p cell lineage is variable within a genetically homogeneous population and polymorphic between wild strains. Irrespective of its division pattern, P3.p is competent to form vulval tissue in different C. elegans strains, whereas it is not competent in C. briggsae. In Oscheius sp. 1, P4.p and P8.p lineages are strongly polymorphic. Within each genus, these intraspecies polymorphisms in cell lineages are amplified between closely related species. In Oscheius sp. 1, the large polymorphisms in P4.p and P8.p lineages allowed us to undertake a genetic analysis of the variation between two pairs of strains. Multiple loci are involved in cell lineage differences, and variation at one locus appears to have a relatively strong effect. In addition to these large lineage variations in cells that do not normally contribute to the vulva, we find minor variations (errors) in vulval lineages, which represent the precision level of the vulval-patterning process and point to a selection pressure for maintenance of a large vulval equivalence group. CONCLUSIONS: Polymorphisms in vulval cell lineage are found within a given nematode species, and could be instrumental in explaining evolutionary variations between closely related species.     

GENOME‐WIDE INVESTIGATION OF REPRODUCTIVE ISOLATION IN EXPERIMENTAL LINEAGES AND NATURAL SPECIES OF NEUROSPORA: IDENTIFYING CANDIDATE REGIONS BY MICROARRAY‐BASED GENOTYPING AND MAPPING

Inherent incompatibilities between genetic components from genomes of different species may cause intrinsic reproductive isolation. In evolution experiments designed to instigate speciation in laboratory populations of the filamentous fungus Neurospora, we previously discovered a pair of incompatibility loci (dfe and dma) that interact negatively to cause severe defects in sexual reproduction. Here we show that the dfe‐dma incompatibility also is a significant cause of genetic isolation between two naturally occurring species of Neurospora (N. crassa and N. intermedia). The strong incompatibility interaction has a simple genetic basis (two biallelic loci) and antagonistic epistasis occurs between heterospecific alleles only, consistent with the Dobzhansky–Muller model of genic incompatibility. We developed microarray‐based, restriction‐site associated DNA (RAD) markers that identified ～1500 polymorphisms between the genomes of the two species, and constructed the first interspecific physical map of Neurospora. With this new mapping resource, the approximate genomic locations of the incompatibility loci were determined using three different approaches: genome scanning, bulk‐segregant analyses, and introgression. These population, quantitative, and classical genetics methods concordantly identified two candidate regions, narrowing the search for each incompatibility locus to only ～2% of the nuclear genome. This study demonstrates how advances in high‐throughput, genome‐wide genotyping can be applied to mapping reproductive isolation genes and speciation research.     

Genetic mapping of species boundaries in Louisiana irises using IRRE retrotransposon display markers

Genetic mapping studies provide insight into the pattern and extent of genetic incompatibilities affecting hybridization between closely related species. Genetic maps of two species of Louisiana Irises, Iris fulva and I. brevicaulis, were constructed from transposon-based molecular markers segregating in reciprocal backcross (BC1) interspecific hybrids and used to investigate genomic patterns of species barriers inhibiting introgression. Linkage mapping analyses indicated very little genetic incompatibility between I. fulva and I. brevicaulis in the form of map regions exhibiting transmission ratio distortion, and this was confirmed using a Bayesian multipoint mapping analysis. These results demonstrate the utility of transposon-based marker systems for genetic mapping studies of wild plant species and indicate that the genomes of I. fulva and I. brevicaulis are highly permeable to gene flow and introgression from one another via backcrossing.     

Hakuna Nematoda: genetic and phenotypic diversity in African isolates of Caenorhabditis elegans and C. briggsae

Caenorhabditis elegans and C. briggsae have many parallels in terms of morphology, life history and breeding system. Both species also share similar low levels of molecular diversity, although the global sampling of natural populations has been limited and geographically biased. In this study, we describe the first cultured isolates of C. elegans and C. briggsae from sub-Saharan Africa. We characterize these samples for patterns of nucleotide polymorphism and vulva precursor cell lineage, and conduct a series of hybrid crosses in C. briggsae to test for genetic incompatibilities. The distribution of genetic diversity confirms a lack of geographic structure to C. elegans sequences but shows genetic differentiation of C. briggsae into three distinct clades that may correspond to three latitudinal ranges. Despite low levels of molecular diversity, we find considerable variation in cell division frequency in African C. elegans for the P3.p vulva precursor cell, and in African C. briggsae for P4.p, a variation that was not previously observed in this species. Hybrid crosses did not reveal major incompatibilities between C. briggsae strains from Africa and elsewhere, and there was some evidence of inbreeding depression. These new African isolates suggest that important ecological factors may be shaping the patterns of diversity in C. briggsae, and that despite many similarities between C. elegans and C. briggsae, there may be more subtle differences in their natural histories than previously appreciated.     

Maximum likelihood techniques for the mapping and analysis of quantitative trait loci with the aid of genetic markers

A method is presented to estimate the biometric parameters of a quantitative trait locus linked to a genetic marker when both loci are segregating in the F-2 generation of a cross between two inbred lines. The method, which assumes underlying normal distributions, is a combination of maximum likelihood and moments methods and uses the statistics of the genetic marker genotype samples for the quantitative trait to estimate the recombination frequency between the two loci and the means and variances of the genotypes of the quantitative trait locus. With this method, the genetic parameters of a locus affecting plant height linked to an electrophoretic marker for esterase were accurately estimated from a sample of 1596 F-2 progeny of a cross between two species of Lycopersicon (tomato). Linkage distance between the two loci was 38 map units and the effect of the quantitative trait locus was 1.6 phenotypic standard deviation units. Accurate estimates of the genetic parameters and linkage distance for populations of 2000 individuals simulated with a segregating codominant locus with an effect of 1.63 standard deviations linked to a genetic marker with .2 recombination were also derived by this method. The method is not effective in distinguishing between complete and partial linkage in samples of only 500 individuals or for quantitative loci with effects less than a phenotypic standard deviation. The method is more effective for codominant than for dominant loci.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Genome-wide associations between hybrid sterility QTL and marker transmission ratio distortion

Marker transmission ratio distortion (TRD) in genetic mapping populations is frequently ascribed to selection against allelic combinations that cause hybrid incompatibility. Accordingly, genomic regions of TRD should be nonrandomly associated (colocated) with loci that underlie hybrid incompatibility. To directly test this hypothesis, we evaluated the genome-wide qualitative and quantitative agreement between chromosomal regions exhibiting marker TRD and those known to contain hybrid incompatibility quantitative trait locus (QTL). Incompatibility data came from a near-isogenic line (NIL) analysis of pollen and seed sterility in a cross between two Solanum (formerly Lycopersicon) species. We assessed (1) whether these incompatibility loci are colocated with markers that show significant TRD in two earlier generations preceding these introgression lines and (2) whether the magnitude of marker distortion quantitatively matches the estimated strength of selection against each incompatibility locus. We found evidence that TRD regions are chromosomally colocated with hybrid incompatibility loci more frequently than is expected by chance: pollen sterility QTLs were most closely associated with distorted heterozygote frequencies in later-generation backcrosses. Nonetheless, there was no evidence for an association between TRD and seed sterility and little evidence of a quantitative association between the magnitude of marker TRD and the fitness effects of heterospecific alleles at each chromosomal location. We propose and test a model (the "dance partner" model) to explain several cases where regions of TRD are not associated with hybrid incompatibility loci. Under this model, some NILs containing greater than one heterospecific introgression may not express hybrid incompatibility phenotypes because they carry both appropriate genetic dance partners required for a fully functional interaction. Accordingly, negative interactions expressed in earlier backcross generations are masked in these double-introgression NILs. Based on this model, we identify the location of several new putative pairwise interactors underlying hybrid incompatibility in this species cross.     

Transmission ratio distortion in intraspecific hybrids of Mimulus guttatus: implications for genomic divergence

We constructed a genetic linkage map between two divergent populations of Mimulus guttatus. We genotyped an F(2) mapping population (N = 539) at 154 AFLP, microsatellite, and gene-based markers. A framework map was constructed consisting of 112 marker loci on 14 linkage groups with a total map length of 1518 cM Kosambi. Nearly half of all markers (48%) exhibited significant transmission ratio distortion (alpha = 0.05). By using a Bayesian multipoint mapping method and visual inspection of significantly distorted markers, we detected 12 transmission ratio distorting loci (TRDL) throughout the genome. The high degree of segregation distortion detected in this intraspecific map indicates substantial genomic divergence that perhaps suggests genomic incompatibilities between these two populations. We compare the pattern of transmission ratio distortion in this map to an interspecific map constructed between M. guttatus and M. nasutus. A similar level of segregation distortion is detected in both maps. Collinear regions between maps are compared to determine if there are shared genetic patterns of non-Mendelian segregation distortion within and among Mimulus species.     

Comparative peptidomics of Caenorhabditis elegans versus C. briggsae by LC-MALDI-TOF MS

Neuropeptides are important signaling molecules that function in cell-cell communication as neurotransmitters or hormones to orchestrate a wide variety of physiological conditions and behaviors. These endogenous peptides can be monitored by high throughput peptidomics technologies from virtually any tissue or organism. The neuropeptide complement of the soil nematode Caenorhabditis elegans has been characterized by on-line two-dimensional liquid chromatography and quadrupole time-of-flight tandem mass spectrometry (2D-nanoLC Q-TOF MS/MS). Here, we use an alternative peptidomics approach combining liquid chromatography (LC) with matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry to map the peptide content of C. elegans and another Caenorhabditis species, Caenorhabditis briggsae. This study allows a better annotation of neuropeptide-encoding genes from the C. briggsae genome and provides a promising basis for further evolutionary comparisons.     

Genetic and physical mapping of the grapevine powdery mildew resistance gene, Run1, using a bacterial artificial chromosome library

Resistance to grapevine powdery mildew is controlled by Run1, a single dominant gene present in the wild grapevine species, Muscadinia rotundifolia, but absent from the cultivated species, Vitis vinifera. Run1 has been introgressed into V. vinifera using a pseudo-backcross strategy, and genetic markers have previously been identified that are linked to the resistance locus. Here we describe the construction of comprehensive genetic and physical maps spanning the resistance locus that will enable future positional cloning of the resistance gene. Physical mapping was performed using a bacterial artificial chromosome (BAC) library constructed using genomic DNA extracted from a resistant V. vinifera individual carrying Run1 within an introgression. BAC contig assembly has enabled 20 new genetic markers to be identified that are closely linked to Run1, and the position of the resistance locus has been refined, locating the gene between the simple sequence repeat (SSR) marker, VMC4f3.1, and the BAC end sequence-derived marker, CB292.294. This region contains two multigene families of resistance gene analogues (RGA). A comparison of physical and genetic mapping data indicates that recombination is severely repressed in the vicinity of Run1, possibly due to divergent sequence contained within the introgressed fragment from M. rotundifolia that carries the Run1 gene.