Genetic architecture of qualitative and quantitative Melampsora larici-populina leaf rust resistance in hybrid poplar: genetic mapping and QTL detection

In order to elucidate the genetic control of resistance to Melampsora larici-populina leaf rust in hybrid poplars, a Populus deltoides x P. trichocarpa F(1) progeny was analysed for qualitative and quantitative rust resistances. This progeny was evaluated for three components of quantitative resistance (latent period, uredinia number and uredinia size) to seven M. larici-populina strains in controlled conditions, and for one component of field susceptibility (rust colonization on the most infected leaf). One qualitative resistance locus inherited from P. deltoides, R(1), was localized on the genetic map. It segregates 1 : 1 in the F(1) progeny and is effective against four of the studied strains. QTL analysis was performed separately on R(1) and r(1) genotype subsets. An additional detection was conducted on the entire F(1) progeny for the three strains able to overcome R(1) and for MAX2. A total of nine QTLs were detected. Two had large, broad-spectrum effects. One (R(US)) is inherited from the P. trichocarpa parent; the other is inherited from P. deltoides and colocalized with R(1). Seven QTLs had only limited and specific effects. Significant interaction effects were detected mainly between the two major QTLs. Implications of these results for durable resistance breeding strategies, and possible benefits from the Populus genome sequence, are discussed.     

Dense genetic linkage maps of three Populus species (Populus deltoides, P. nigra and P. trichocarpa) based on AFLP and microsatellite markers

Populus deltoides, P. nigra, and P. trichocarpa are the most important species for poplar breeding programs worldwide. In addition, Populus has become a model for fundamental research on trees. Linkage maps were constructed for these three species by analyzing progeny of two controlled crosses sharing the same female parent, Populus deltoides cv. S9-2 x P. nigra cv. Ghoy and P. deltoides cv. S9-2 x P. trichocarpa cv. V24. The two-way pseudotestcross mapping strategy was used to construct the maps. Amplified fragment length polymorphism (AFLP) markers that segregated 1:1 were used to form the four parental maps. Microsatellites and sequence-tagged sites were used to align homoeologous groups between the maps and to merge linkage groups within the individual maps. Linkage analysis and alignment of the homoeologous groups resulted in 566 markers distributed over 19 groups for P. deltoides covering 86% of the genome, 339 markers distributed over 19 groups for P. trichocarpa covering 73%, and 369 markers distributed over 28 groups for P. nigra covering 61%. Several tests for randomness showed that the AFLP markers were randomly distributed over the genome.     

Comparative genome mapping among Populus adenopoda, P. alba, P. deltoides, P. euramericana and P. trichocarpa

Among the genus Populus, the sections Populus (white poplar), Aigeiros Duby (black poplar) and Tacamahaca Spach contain many tree species of economical and ecological important properties. Two parental maps for the inter-specific hybrid population of Populus adenopoda × P. alba (two species of Populus section) were constructed based on SSR and SRAP markers by means of a two-way pseudo-test cross mapping strategy. The same set of SSR markers developed from the P. trichocarpa (belonging to Tacamahaca section) genome which were used to construct the maps of P. deltoides and P. euramericana (two species of Aigeiros section) was chosen to analyze the genotype of the experimental population of P. adenopoda × P. alba. Using the mapped SSR markers as allelic bridges, the alignment of the white and black poplar maps to each other and to the P. trichocarpa physical map was conducted. The alignment showed high degree of marker synteny and colinearity and the closer relationship between Aigeiros and Tacamahaca sections than that of Populus and Tacamahaca. Moreover, there was evidence for the chromosomal duplication and inter-chromosomal reorganization involving some poplar linkage groups, suggesting a complicated course of fission or fusion in one of the lineages. A poplar consensus map based on the comparisons could be constructed will be useful in practical applications including marker assisted selection.     

Differential susceptibility of poplar hybrids to the aphid Chaitophorus leucomelas (Homoptera: Aphididae)

After its recent introduction to Chile, the aphid Chaitophorus leucomelas Koch is becoming a serious pest affecting commercial poplar, Populus spp., plantations. The pattern of natural infestation of C. leucomelas among poplar hybrids with different pedigrees and the aphid intrinsic rate of increase (r(m)), of C. leucomelas were assessed in the field. In most of the hybrids, aphid abundance peaked in March (late summer). Among 12 types of poplar crosses, [(P. trichocarpa Torr. & Gray x P. deltoides Bartram ex Marshall) x (P. trichocarpa x P. deltoides)] and [(P. trichocarpa x P. maximowiczii Henry) x P. maximowiczii] showed the highest and lowest aphid densities, respectively. A trend to find more aphids in branch bases was apparent. The intrinsic rate of C. leucomelas increase was higher in [(P. trichocarpa x P. deltoides) x P. deltoides] hybrids, and lower in [(P. trichocarpa x P. maximowiczii) x P. trichocarpa] hybrids. Aphid density and performance were higher in hybrids with P. deltoides parentage, whereas hybrids with P. maximowiczii parentage showed lower aphid densities and performance. Hybrids with P. nigra L. parentage, namely, [P. trichocarpa x P. nigra], also had high aphid density, but aphid performance was lower compared with hybrids with P. deltoides parentage. These results suggest that among poplar hybrids studied, susceptibility to C. leucomelas is inherited through P. deltoides, whereas resistance seems to be inherited through P. maximowiczii. Thus, P. maximowiczii hybrids are recommended for commercial or ornamental planting programs in zones where there is a high risk of aphid infestation.     

A Radiation Hybrid Map of 95 STSs Spanning Human Chromosome 13q

We have constructed a high-resolution physical map of the long arm of human chromosome 13 using a panel of 94 radiation hybrids. A comprehensive map of 95 chromosome 13-specific sequence tagged sites (STSs) spanning 13q from the presumed centromere at D13Z1 to the known telomere was obtained by multipoint maximum likelihood statistical methods. The 95 markers have an average retention frequency of 10%, with markers closer to the centromere having much greater retention frequencies (22-49%) than distal 13q markers (2-12%) The most likely radiation hybrid map localized the 95 STSs into 54 unique map positions, 34 with odds of 1000:1 or greater; the comprehensive map localized all but 17 STSs with odds exceeding 10:1. The total map length of 13q was 1302 cR₉₀₀₀ (range 6.4-94.4 cR₉₀₀₀) and a physical distance of 98 Mb, so that 1% breakage in the RH panel corresponds to 75 kb. A comparison of the comprehensive RH map to genetic maps of chromosome 13q shows identical locus orders for the common markers, with two exceptions over 1-cM distances. We discuss the possible relationships between the genetic and the radiation hybrid maps.     

Exploiting comparative mapping among Brassica species to accelerate the physical delimitation of a genic male-sterile locus (BnRf) in Brassica napus

The recessive genic male sterility (RGMS) line 9012AB has been used as an important pollination control system for rapeseed hybrid production in China. Here, we report our study on physical mapping of one male-sterile locus (BnRf) in 9012AB by exploiting the comparative genomics among Brassica species. The genetic maps around BnRf from previous reports were integrated and enriched with markers from the Brassica A7 chromosome. Subsequent collinearity analysis of these markers contributed to the identification of a novel ancestral karyotype block F that possibly encompasses BnRf. Fourteen insertion/deletion markers were further developed from this conserved block and genotyped in three large backcross populations, leading to the construction of high-resolution local genetic maps where the BnRf locus was restricted to a less than 0.1-cM region. Moreover, it was observed that the target region in Brassica napus shares a high collinearity relationship with a region from the Brassica rapa A7 chromosome. A BnRf-cosegregated marker (AT3G23870) was then used to screen a B. napus bacterial artificial chromosome (BAC) library. From the resulting 16 positive BAC clones, one (JBnB089D05) was identified to most possibly contain the BnRf (c) allele. With the assistance of the genome sequence from the Brassica rapa homolog, the 13.8-kb DNA fragment covering both closest flanking markers from the BAC clone was isolated. Gene annotation based on the comparison of microcollinear regions among Brassica napus, B. rapa and Arabidopsis showed that five potential open reading frames reside in this fragment. These results provide a foundation for the characterization of the BnRf locus and allow a better understanding of the chromosome evolution around BnRf.     

A physical map of the highly heterozygous Populus genome: integration with the genome sequence and genetic map and analysis of haplotype variation

As part of a larger project to sequence the Populus genome and generate genomic resources for this emerging model tree, we constructed a physical map of the Populus genome, representing one of the few such maps of an undomesticated, highly heterozygous plant species. The physical map, consisting of 2802 contigs, was constructed from fingerprinted bacterial artificial chromosome (BAC) clones. The map represents approximately 9.4-fold coverage of the Populus genome, which has been estimated from the genome sequence assembly to be 485 ± 10 Mb in size. BAC ends were sequenced to assist long-range assembly of whole-genome shotgun sequence scaffolds and to anchor the physical map to the genome sequence. Simple sequence repeat-based markers were derived from the end sequences and used to initiate integration of the BAC and genetic maps. A total of 2411 physical map contigs, representing 97% of all clones assigned to contigs, were aligned to the sequence assembly (JGI Populus trichocarpa , version 1.0). These alignments represent a total coverage of 384 Mb (79%) of the entire poplar sequence assembly and 295 Mb (96%) of linkage group sequence assemblies. A striking result of the physical map contig alignments to the sequence assembly was the co-localization of multiple contigs across numerous regions of the 19 linkage groups. Targeted sequencing of BAC clones and genetic analysis in a small number of representative regions showed that these co-aligning contigs represent distinct haplotypes in the heterozygous individual sequenced, and revealed the nature of these haplotype sequence differences.     

Physical mapping of the barley stem rust resistance gene rpg4

The barley stem rust resistance gene rpg4 was physically and genetically localized on two overlapping BAC clones covering an estimated 300-kb region of the long arm of barley chromosome 7(5H). Initially, our target was mapped within a 6.0-cM region between the previously described flanking markers MWG740 and ABG391. This region was then saturated by integrating new markers from several existing barley and rice maps and by using BAC libraries of barley cv. Morex and rice cv. Nipponbare. Physical/genetic distances in the vicinity of rpg4 were found to be 1.0 Mb/cM, which is lower than the average for barley (4 Mb/cM) and lower than that determined by translocation breakpoint mapping (1.8 Mb/cM). Synteny at high resolution levels has been established between the region of barley chromosome 7(5H) containing the rpg4 locus and the subtelomeric region of rice chromosome 3 between markers S16474 and E10757. This 1.7-cM segment of the rice genome was covered by two overlapping BAC clones, about 250 kb of total length. In barley the markers S16474 and E10757 genetically delimit rpg4, lying 0.6 cM distal and 0.4 cM proximal to the locus, respectively.     

Alignment of Escherichia coli K12 DNA sequences to a genomic restriction map.

We use the extensive published information describing the genome of Escherichia coli and new restriction map alignment software to align DNA sequence, genetic, and physical maps. Restriction map alignment software is used which considers restriction maps as strings analogous to DNA or protein sequences except that two values, enzyme name and DNA base address, are associated with each position on the string. The resulting alignments reveal a nearly linear relationship between the physical and genetic maps of the E. coli chromosome. Physical map comparisons with the 1976, 1980, and 1983 genetic maps demonstrate a better fit with the more recent maps. The results of these alignments are genomic kilobase coordinates, orientation and rank of the alignment that best fits the genetic data. A statistical measure based on extreme value distribution is applied to the alignments. Additional computer analyses allow us to estimate the accuracy of the published E. coli genomic restriction map, simulate rearrangements of the bacterial chromosome, and search for repetitive DNA. The procedures we used are general enough to be applicable to other genome mapping projects.     

Closure of a genetic linkage map of human chromosome 7q with centromere and telomere polymorphisms.

We have constructed a 2.4-cM resolution genetic linkage map for chromosome 7q that is bounded by centromere and telomere polymorphisms and contains 66 loci (88 polymorphic systems), 38 of which are uniquely placed with odds for order of at least 1000:1. Ten genes are included in the map and 11 markers have heterozygosities of at least 70%. This map is the first to incorporate several highly informative markers derived from a telomere YAC clone HTY146 (locus D7S427), including HTY146c3 (HET 92%). The telomere locus markers span at least 200 kb of the 7q terminus and no crossovers within the physical confines of the locus were observed in approximately 240 jointly informative meioses. The sex-equal map length is 158 cM and the largest genetic interval between uniquely localized markers in this map is 11 cM. The female and male map lengths are 181 and 133 cM, respectively. The map is based on the CEPH reference pedigrees and includes over 4000 new genotypes, our previously reported data plus 29 allele systems from the published CEPH version 5 database, and was constructed using the program package CRI-MAP. This genetic linkage map can be considered a baseline map for 7q, and will be useful for defining the extent of chromosome deletions previously reported for breast and prostate cancers, for developing additional genetic maps such as index marker and 1-cM maps, and ultimately for developing a fully integrated genetic and physical map for this chromosome.     

Restriction fragment length polymorphisms of a chloroplast photosystem II gene from poplar and their use for species identification.

In a total DNA library from the poplar clone Beaupré (Populus trichocarpa x P. deltoides) one DNA clone was found to identify restriction site polymorphisms in different poplar species. This clone represents a cpDNA gene that shows close homology to a photosystem II gene of pea and spinach coding for the D2 protein and the 44 kDa reaction centre. In Southern blot analysis this probe identified interspecific restriction site variation among the different poplar species; intraspecific variation was not detectable. As the chloroplast genome is maternally inherited in poplars this cpDNA probe was used for identification of P. nigra or P. deltoides as the seed parents of F1 hybrid trees in natural stands of western Germany.     

Molecular linkage maps of the Populus genome.

We report molecular genetic linkage maps for an interspecific hybrid population of Populus, a model system in forest-tree biology. The hybrids were produced by crosses between P. deltoides (mother) and P. euramericana (father), which is a natural hybrid of P. deltoides (grandmother) and P. nigra (grandfather). Linkage analysis from 93 of the 450 backcross progeny grown in the field for 15 years was performed using random amplified polymorphic DNAs (RAPDs), amplified fragment length polymorphisms (AFLPs), and inter-simple sequence repeats (ISSRs). Of a total of 839 polymorphic markers identified, 560 (67%) were testcross markers heterozygous in one parent but null in the other (segregating 1:1), 206 (25%) were intercross dominant markers heterozygous in both parents (segregating 3:1), and the remaining 73 (9%) were 19 non-parental RAPD markers (segregating 1:1) and 54 codominant AFLP markers (segregating 1:1:1:1). A mixed set of the testcross markers, non-parental RAPD markers, and codominant AFLP markers was used to construct two linkage maps, one based on the P. deltoides (D) genome and the other based on P. euramericana (E). The two maps showed nearly complete coverage of the genome, spanning 3801 and 3452 cM, respectively. The availability of non-parental RAPD and codominant AFLP markers as orthologous genes allowed for a direct comparison of the rate of meiotic recombination between the two different parental species. Generally, the rate of meiotic recombination was greater for males than females in our interspecific poplar hybrids. The confounded effect of sexes and species causes the mean recombination distance of orthologous markers to be 11% longer for the father (P. euramericana; interspecific hybrid) than for the mother (P. deltoides; pure species). The linkage maps constructed and the interspecific poplar hybrid population in which clonal replicates for individual genotypes are available present a comprehensive foundation for future genomic studies and quantitative trait locus (QTL) identification.     

Quantitative genetic variation of leaf size and shape in a mixed diploid and triploid population of Populus

In the interspecific cross of Populus trichocarpa x P. deltoides, unexpected simultaneous occurrence of diploid hybrids and triploid hybrids (with two alleles from the female parent and one from the male parent at each locus) led us to examine the evolutionary genetic significance of this phenomenon. As expected, leaf size and shape of the triploid progeny are closer to the female P. trichocarpa than male P. deltoides parent. Although the pure triploid progeny population did not have higher genetic variance in leaf traits than the pure diploid population, the former appears to hide much non-additive genetic variance and display strong genetic control over the phenotypic plasticity of leaf traits. It is suggested that the cryptic non-additive variance, especially epistasis, can be released when a population is disturbed by changes in the environment. A mixed diploid and triploid progeny population combines phenotypic and genetic characteristics of both pure hybrids and is considered to be of adaptive significance for populars to survive and evolve in a fluctuating environment. The significant effect due to general and specific combining ability differences at the population level suggests that the population divergence of these two species is under additive and non-additive genetic control.     

Microsatellite DNA fingerprinting, differentiation, and genetic relationships of clones, cultivars, and varieties of six poplar species from three sections of the genus Populus.

Accurate identification of Populus clones and cultivars is essential for effective selection, breeding, and genetic resource management programs. The unit of cultivation and breeding in poplars is a clone, and individual cultivars are normally represented by a single clone. Microsatellite DNA markers of 10 simple sequence repeat loci were used for genetic fingerprinting and differentiation of 96 clones/cultivars and varieties belonging to six Populus species (P. deltoides, P. nigra, P. balsamifera, P. trichocarpa, P. grandidentata, and P maximowiczii) from three sections of the genus. All 96 clones/cultivars could be uniquely fingerprinted based on their single- or multilocus microsatellite genotypes. The five P. grandidentata clones could be differentiated based on their single-locus genotypes, while six clones of P. trichocarpa and 11 clones of P. maximowiczii could be identified by their two-locus genotypes. Twenty clones of P. deltoides and 25 clones of P. nigra could be differentiated by their multilocus genotypes employing three loci, and 29 clones of P. balsamifera required the use of multilocus genotypes at five loci for their genetic fingerprinting and differentiation. The loci PTR3, PTR5, and PTR7 were found to be the most informative for genetic fingerprinting and differentiation of the clones. The mean number of alleles per locus ranged from 2.9 in P. trichocarpa or P. grandidentata to 6.0 in P. balsamifera and 11.2 in 96 clones of the six species. The mean number of observed genotypes per locus ranged from 2.4 in P. grandidentata to 7.4 in P. balsamifera and 19.6 in 96 clones of the six species. The mean number of unique genotypes per locus ranged from 1.3 in P. grandidentata to 3.9 in P. deltoides and 8.8 in 96 clones of the six species. The power of discrimination of the microsatellite DNA markers in the 96 clones ranged from 0.726 for PTR4 to 0.939 for PTR7, with a mean of 0.832 over the 10 simple sequence repeat loci. Clones/cultivars from the same species showed higher microsatellite DNA similarities than the clones from the different species. A UPGMA cluster plot constructed from the microsatellite genotypic similarities separated the 96 clones into six major groups corresponding to their species. Populus nigra var. italica clones were genetically differentiated from the P. nigra var. nigra clones. Microsatellite DNA markers could be useful in genetic fingerprinting, identification, classification, certification, and registration of clones, clultivars, and varieties as well as genetic resource management and protection of plant breeders' rights in Populus.     

A microarray-based genotyping and genetic mapping approach for highly heterozygous outcrossing species enables localization of a large fraction of the unassembled Populus trichocarpa genome sequence

Microarrays have demonstrated significant power for genome-wide analyses of gene expression, and recently have also revolutionized the genetic analysis of segregating populations by genotyping thousands of loci in a single assay. Although microarray-based genotyping approaches have been successfully applied in yeast and several inbred plant species, their power has not been proven in an outcrossing species with extensive genetic diversity. Here we have developed methods for high-throughput microarray-based genotyping in such species using a pseudo-backcross progeny of 154 individuals of Populus trichocarpa and P. deltoides analyzed with long-oligonucleotide in situ- synthesized microarray probes . Our analysis resulted in high-confidence genotypes for 719 single-feature polymorphism (SFP) and 1014 gene expression marker (GEM) candidates. Using these genotypes and an established microsatellite (SSR) framework map, we produced a high-density genetic map comprising over 600 SFPs, GEMs and SSRs. The abundance of gene-based markers allowed us to localize over 35 million base pairs of previously unplaced whole-genome shotgun (WGS) scaffold sequence to putative locations in the genome of P. trichocarpa . A high proportion of sampled scaffolds could be verified for their placement with independently mapped SSRs, demonstrating the previously un-utilized power that high-density genotyping can provide in the context of map-based WGS sequence reassembly. Our results provide a substantial contribution to the continued improvement of the Populus genome assembly, while demonstrating the feasibility of microarray-based genotyping in a highly heterozygous population. The strategies presented are applicable to genetic mapping efforts in all plant species with similarly high levels of genetic diversity.     

Exploitation of Arabidopsis-tomato synteny to construct a high-resolution map of the ovatecontaining region in tomato chromosome 2.

High-resolution genetic and physical maps were constructed for the region of chromosome 2 containing the major fruit-shape locus ovate. A total of 3,000 NIL F2 and F3 NILs derived from Lycopersicon esculentum cv. Yellow Pear (TA503) x L. pennellii (a wild tomato) were used to position ovate adjacent to the marker TG645 and flanked by markers TX700 and BA10R (a 0.03-cM interval). BAC libraries and a BIBAC library were screened with the closest marker, TG645. Genetic mapping with the ends of isolated BAC clones revealed that two BAC clones (100 and 140 kb) both contained the ovate locus. Screening of sequences from these BAC clones revealed synteny between this segment of tomato chromosome 2 and the chromosome-4 region of Arabidopsis containing the BAC clone ATAP22. Microsynteny between the two genomes was exploited to find additional markers near the ovate locus. The placement of ovate on a BAC clone will now allow cloning of this locus and, hence, may open the door to understanding the molecular basis of fruit development and also facilitate the genetic engineering of fruit-shape characteristics. This also represents the first time that microsynteny with Arabidopsis has been exploited for positional cloning purposes in a different plant family.     

Integration of physical and genetic maps in apple confirms whole-genome and segmental duplications in the apple genome

A total of 355 simple sequence repeat (SSR) markers were developed, based on expressed sequence tag (EST) and bacterial artificial chromosome (BAC)-end sequence databases, and successfully used to construct an SSR-based genetic linkage map of the apple. The consensus linkage map spanned 1143 cM, with an average density of 2.5 cM per marker. Newly developed SSR markers along with 279 SSR markers previously published by the HiDRAS project were further used to integrate physical and genetic maps of the apple using a PCR-based BAC library screening approach. A total of 470 contigs were unambiguously anchored onto all 17 linkage groups of the apple genome, and 158 contigs contained two or more molecular markers. The genetically mapped contigs spanned ～421 Mb in cumulative physical length, representing 60.0% of the genome. The sizes of anchored contigs ranged from 97 kb to 4.0 Mb, with an average of 995 kb. The average physical length of anchored contigs on each linkage group was ～24.8 Mb, ranging from 17.0 Mb to 37.73 Mb. Using BAC DNA as templates, PCR screening of the BAC library amplified fragments of highly homologous sequences from homoeologous chromosomes. Upon integrating physical and genetic maps of the apple, the presence of not only homoeologous chromosome pairs, but also of multiple locus markers mapped to adjacent sites on the same chromosome was detected. These findings demonstrated the presence of both genome-wide and segmental duplications in the apple genome and provided further insights into the complex polyploid ancestral origin of the apple.     

Primary, nonsyndromic vesicoureteric reflux and its nephropathy is genetically heterogeneous, with a locus on chromosome 1

Primary vesicoureteric reflux (VUR) affects 1%-2% of whites, and reflux nephropathy (RN) causes up to 15% of end-stage renal failure in children and adults. There is a 30-50-fold increased incidence of VUR in first-degree relatives of probands, compared with the general population. We report the results of the first genomewide search of VUR and RN; we studied seven European families whose members exhibit apparently dominant inheritance. We initially typed 387 polymorphic markers spaced, on average, at 10 cM throughout the genome; we used the GENEHUNTER program to provide parametric and nonparametric linkage analyses of affected individuals. The most positive locus spanned 20 cM on 1p13 between GATA176C01 and D1S1653 and had a nonparametric LOD score (NPL) of 5.76 (P=.0002) and a parametric LOD score of 3.16. Saturation with markers at 1-cM intervals increased the NPL to 5.94 (P=.00009). Hence, VUR maps to a locus on chromosome 1. There was evidence of genetic heterogeneity at the chromosome 1 locus, and 12 additional loci were identified genomewide, with P<.05. No significant linkage was found to 6p, where a renal and ureteric malformation locus has been reported, or to PAX2, mutations of which cause VUR in renal-coloboma syndrome. Our results support the hypothesis that VUR is a genetic disorder.     

Fine mapping and SNP analysis of positional candidates at the preeclampsia susceptibility locus (PREG1) on chromosome 2

Genome scans in Icelandic, Australian and New Zealand, and Finnish families have localized putative susceptibility loci for preeclampsia/ eclampsia to chromosome 2. The locus mapped in the Australian and New Zealand study (designated PREG1) was thought to be the same locus as that identified in the Icelandic study. In both these studies, two distinct quantitative trait locus (QTL) regions were evident on chromosome 2. Here, we describe our fine mapping of the PREG1 locus and a genetic analysis of two positional candidate genes. Twenty-five additional microsatellite markers were genotyped within the 74-cM linkage region defined by the combined Icelandic and Australian and New Zealand genome scans. The overall position and shape of the localization evidence obtained using nonparametric multipoint analysis did not change from that seen previously in our 10-cM resolution genome scan; two peaks were displayed, one on chromosome 2p at marker D2S388 (107.46 cM) and the other on chromosome 2q at 151.5 cM at marker D2S2313. Using the robust two-point linkage analysis implemented in the Analyze program, all 25 markers gave positive LOD scores with significant evidence of linkage being seen at marker D2S2313 (151.5 cM), achieving a LOD score of 3.37 under a strict diagnostic model. Suggestive evidence of linkage was seen at marker D2S388 (107.46 cM) with a LOD score of 2.22 under the general diagnostic model. Two candidate genes beneath the peak on chromosome 2p were selected for further analysis using public single nucleotide polymorphisms (SNPs) within these genes. Maximum LOD scores were obtained for an SNP in TACR1 (LOD = 3.5) and for an SNP in TCF7L1 (LOD = 3.33), both achieving genome-wide significance. However, no evidence of association was seen with any of the markers tested. These data strongly support the presence of a susceptibility gene on chromosome 2p11-12 and substantiate the possibility of a second locus on chromosome 2q23.