Haplotype block and superblock structures of the alpha1-adrenergic receptor genes reveal echoes from the chromosomal past

A significant proportion of the human genome is contained within haplotype blocks across which pairwise linkage disequilibrium (LD) is very high. However, LD is also often high between markers at more remote distances, and within different haplotype blocks. Here, we evaluate the origins of haplotype block structure in the three genes for alpha1 adrenergic receptors (alpha1-AR) in the human genome ( ADRA1A, ADRA1B and ADRA1D) by genotyping dense single-nucleotide polymorphism (SNP) marker maps, and show that LD signals between distant markers are due to the presence of extended haplotype superblocks in individuals with ancient chromosomes which have escaped historic recombination. ARs mediate the physiological effects of epinephrine and norepinephrine, and are targets of many therapeutic drugs. This work has identified haplotype backgrounds of alpha1-AR missense variants, haplotype block structures in US Caucasians and African Americans, and haplotype tag SNPs for each block, and we present strong evidence for ancient haplotype block superstructure at these genes which has been partially disrupted by recombination, and evidence for reinstatement of linkage disequilibrium by subsequent recombination events. ADRA1A is comprised of four haplotype blocks in US Caucasians, while in African Americans Block 1 is split. ADRA1B has four blocks in US Caucasians, but in African Americans only the first two blocks are present. ADRA1D has two blocks in US Caucasians, and the first block is replaced by two smaller blocks in African Americans. For both ADRA1A and ADRA1B, haplotype superstructures may represent a novel, higher-level hierarchy in the human genome, which may reduce redundancy of testing by further aggregation of genotype data.Electronic Supplementary Material Supplementary material is available in the online version of this article at Communicated by W. R. McCombie     

Optimal haplotype structure for linkage disequilibrium-based fine mapping of quantitative trait loci using identity by descent

A linkage disequilibrium-based method for fine mapping quantitative trait loci (QTL) has been described that uses similarity between individuals' marker haplotypes to determine if QTL alleles are identical by descent (IBD) to model covariances among individuals' QTL alleles for a mixed linear model. Mapping accuracy with this method was found to be sensitive to the number of linked markers that was included in the haplotype when fitting the model at a putative position of the QTL. The objective of this study was to determine the optimal haplotype structure for this IBD-based method for fine mapping a QTL in a previously identified QTL region. Haplotypes consisting of 1, 2, 4, 6, or all 10 available markers were fit as a "sliding window" across the QTL region under ideal and nonideal simulated population conditions. It was found that using haplotypes of 4 or 6 markers as a sliding "window" resulted in the greatest mapping accuracy under nearly all conditions, although the true IBD state at a putative QTL position was most accurately predicted by IBD probabilities obtained using all markers. Using 4 or 6 markers resulted in greater discrimination of IBD probabilities between positions while maintaining sufficient accuracy of IBD probabilities to detect the QTL. Fitting IBD probabilities on the basis of a single marker resulted in the worst mapping accuracy under all conditions because it resulted in poor accuracy of IBD probabilities. In conclusion, for fine mapping using IBD methods, marker information must be used in a manner that results in sensitivity of IBD probabilities to the putative position of the QTL while maintaining sufficient accuracy of IBD probabilities to detect the QTL. Contrary to expectation, use of haplotypes of 4-6 markers to derive IBD probabilities, rather than all available markers, best fits these criteria. Thus for populations similar to those simulated here, optimal mapping accuracy for this IBD-based fine-mapping method is obtained with a haplotype structure including a subset of all available markers.     

Recombination hotspots and population structure in Plasmodium falciparum

Understanding the influences of population structure, selection, and recombination on polymorphism and linkage disequilibrium (LD) is integral to mapping genes contributing to drug resistance or virulence in Plasmodium falciparum. The parasite's short generation time, coupled with a high cross-over rate, can cause rapid LD break-down. However, observations of low genetic variation have led to suggestions of effective clonality: selfing, population admixture, and selection may preserve LD in populations. Indeed, extensive LD surrounding drug-resistant genes has been observed, indicating that recombination and selection play important roles in shaping recent parasite genome evolution. These studies, however, provide only limited information about haplotype variation at local scales. Here we describe the first (to our knowledge) chromosome-wide SNP haplotype and population recombination maps for a global collection of malaria parasites, including the 3D7 isolate, whose genome has been sequenced previously. The parasites are clustered according to continental origin, but alternative groupings were obtained using SNPs at 37 putative transporter genes that are potentially under selection. Geographic isolation and highly variable multiple infection rates are the major factors affecting haplotype structure. Variation in effective recombination rates is high, both among populations and along the chromosome, with recombination hotspots conserved among populations at chromosome ends. This study supports the feasibility of genome-wide association studies in some parasite populations.     

Microsatellite based linkage disequilibrium analyses reveal <Emphasis Type="Italic">Saltol</Emphasis> haplotype fragmentation and identify novel QTLs for seedling stage salinity tolerance in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

A set of 84 diverse rice genotypes were assessed for seedling stage salt tolerance and their genetic diversity using 41 polymorphic SSR markers comprising of 19 Saltol QTL linked and 22 random markers. Phenotypic screening under hydroponics identified three indica landraces (Badami, Shah Pasand and Pechi Badam), two Oryza rufipogon accessions (NKSWR2 and NKSWR17) and one each of Basmati rice (Seond Basmati) and japonica cultivars (Tompha Khau) as salt tolerant, having similar tolerance as of Pokkali and FL478. Among the salt tolerant genotypes, biomass showed positive correlation with shoot fresh weight and negative association with root and shoot Na⁺ content. The results indicated repression of Na⁺ loading within the tolerant plants. Linkage disequilibrium (LD) of the Saltol linked markers was weak, suggestive of high fragmentation of Pokkali haplotype, a result of evolutionary active recombination events. Poor haplotype structure of the Saltol region, may reduce its usefulness in marker assisted breeding programmes, if the target foreground markers chosen are wide apart. LD mapping identified eight robust marker-trait associations (QTLs), of which RM10927 was found linked to root and shoot Na⁺ content and RM10871 with shoot Na⁺/K⁺ ratio. RM271 on chromosome 10, an extra Saltol marker, was found associated to root Na⁺/K⁺ ratio. This marker showed a distinct allele among O. rufipogon accessions. There were also other novel loci detected on chromosomes 2, 5 and 10 influencing salt tolerance in the tested germplasm. Although Saltol remained as the key locus, the role of other genomic regions cannot be neglected in tailoring seedling stage salt tolerance in rice.     

Influence of genetic background and heterozygosity on meiotic recombination in Arabidopsis thaliana.

Plant breeding relies on genetic variability generated by meiotic recombination. Control of recombination frequencies is not yet possible, but would significantly extend the options for plant-breeding strategies. A prerequisite would be variability of recombination frequencies. In this study, 15 transgenic kanamycin (KR) and hygromycin (HR) resistance gene insertions mapping to the five Arabidopsis thaliana chromosomes were used as genetic markers. Recombination frequencies were determined from the frequencies of resistance phenotypes within populations segregating for linked KR and HR markers. Recombination frequencies of marker pairs were compared among these four ecotypes, among F1s in both reciprocal forms derived from these ecotypes, and between F1s and their parent lines. On average, the recombination frequencies in F1 crosses were substantially higher (up to 2-fold) than in the homozygous parental ecotypes. A strong negative correlation between genetic similarities of ecotypes and recombination frequencies was detected for two adjacent marker pairs located on the long arm of chromosome 3, but not for marker pairs in other genomic regions. Our results suggest that heterozygosity influences recombination in plant breeding, and cannot be ignored in genetic mapping of genomes.     

Deriving haplotypes through recombination and gene conversion pathways

Retracing the trajectories of past genetic events is crucial to understand the structure of the genome, both in individuals and across populations. A haplotype describes a string of polymorphic sites along a DNA segment. Haplotype diversity is due to mutations creating new variants, and to recombinations and gene conversions that mix and redistribute these variants among individual chromosomes in populations. A number of studies have revealed a relatively simple pattern of haplotype diversity in the human genome, dominated by a few common haplotypes representing founder ancestral ones. New haplotypes are usually rare and have a limited geographic distribution. We propose a method to derive a new haplotype from a set of putative ancestral haplotypes, once mutations in place, through minimal recombination and gene conversion pathways. We describe classes of pathways that represent the whole set of minimal pathways leading to a new haplotype. We show that obtaining this set of pathways can be represented as a problem of finding "secondary structures" of minimum energy. We present a polynomial algorithm solving this folding problem.     

The structure of linkage disequilibrium at the DBH locus strongly influences the magnitude of association between diallelic markers and plasma dopamine beta-hydroxylase activity

There is currently a great deal of interest in using linkage disequilibrium (LD) mapping to locate both disease and quantitative-trait loci on a genomewide scale. Recent findings suggest that much of the human genome is organized in discrete "blocks" of low haplotype diversity, but the utility of such blocks in identifying genes influencing complex traits is not yet known and must ultimately be tested empirically through use of real data. We recently identified a putative functional polymorphism (-1021C-->T) in the 5' upstream region of the DBH gene that accounted for 35%-52% of the total phenotypic variance in plasma dopamine beta-hydroxylase (DBH) activity in samples from three distinct populations. In the present study, we genotyped 11 diallelic markers at the DBH locus surrounding -1021C-->T in 386 unrelated individuals of European origin. We identified a single 10-kb block containing -1021C-->T, in which four haplotypes comprised 93% of the observed chromosomes. Only markers within the block were highly associated with phenotype (P< or =2.2 x 10(-10)), with one exception. In general, association with phenotype was strongly correlated with the degree of LD between each marker and -1021C-->T. Of four LD measures assessed, d(2) was the best predictor of this relationship. Had one attempted to map quantitative-trait loci for plasma DBH activity on a genomewide basis without prior knowledge of candidate regions and not included (by chance) markers within this haplotype block, the DBH locus might have been missed entirely. These results provide a direct example of the potential value of constructing a haplotype map of the human genome prior to embarking on large-scale association studies.     

Recombination Hotspots and Population Structure in Plasmodium falciparum

Understanding the influences of population structure, selection, and recombination on polymorphism and linkage disequilibrium (LD) is integral to mapping genes contributing to drug resistance or virulence in Plasmodium falciparum. The parasite's short generation time, coupled with a high cross-over rate, can cause rapid LD break-down. However, observations of low genetic variation have led to suggestions of effective clonality: selfing, population admixture, and selection may preserve LD in populations. Indeed, extensive LD surrounding drug-resistant genes has been observed, indicating that recombination and selection play important roles in shaping recent parasite genome evolution. These studies, however, provide only limited information about haplotype variation at local scales. Here we describe the first (to our knowledge) chromosome-wide SNP haplotype and population recombination maps for a global collection of malaria parasites, including the 3D7 isolate, whose genome has been sequenced previously. The parasites are clustered according to continental origin, but alternative groupings were obtained using SNPs at 37 putative transporter genes that are potentially under selection. Geographic isolation and highly variable multiple infection rates are the major factors affecting haplotype structure. Variation in effective recombination rates is high, both among populations and along the chromosome, with recombination hotspots conserved among populations at chromosome ends. This study supports the feasibility of genome-wide association studies in some parasite populations.     

Durability of marker-quantitative trait loci haplotypes in structured populations

Given the relative ease of identifying genetic markers linked to QTL (compared to finding the loci themselves), it is natural to ask whether linked markers can be used to address questions concerning the contemporary dynamics and recent history of the QTL. In particular, can a marker allele found associated with a QTL allele in a QTL mapping study be used to track population dynamics or the history of the QTL allele? For this strategy to succeed, the marker-QTL haplotype must persist in the face of recombination over the relevant time frame. Here we investigate the dynamics of marker-QTL haplotype frequencies under recombination, population structure, and divergent selection to assess the potential utility of linked markers for a population genetic study of QTL. For two scenarios, described as "secondary contact" and "novel allele," we use both deterministic and stochastic methods to describe the influence of gene flow between habitats, the strength of divergent selection, and the genetic distance between a marker and the QTL on the persistence of marker-QTL haplotypes. We find that for most reasonable values of selection on a locus (s < or = 0.5) and migration (m > 1%) between differentially selected populations, haplotypes of typically spaced markers (5 cM) and QTL do not persist long enough (>100 generations) to provide accurate inference of the allelic state at the QTL.     

Comparing linkage disequilibrium-based methods for fine mapping quantitative trait loci

Recently, a method for fine mapping quantitative trait loci (QTL) using linkage disequilibrium was proposed to map QTL by modeling covariance between individuals, due to identical-by-descent (IBD) QTL alleles, on the basis of the similarity of their marker haplotypes under an assumed population history. In the work presented here, the advantage of using marker haplotype information for fine mapping QTL was studied by comparing the IBD-based method with 10 markers to regression on a single marker, a pair of markers, or a two-locus haplotype under alternative population histories. When 10 markers were genotyped, the IBD-based method estimated the position of the QTL more accurately than did single-marker regression in all populations. When 20 markers were genotyped for regression, as single-marker methods do not require knowledge of haplotypes, the mapping accuracy of regression in all populations was similar to or greater than that of the IBD-based method using 10 markers. Thus for populations similar to those simulated here, the IBD-based method is comparable to single-marker regression analysis for fine mapping QTL.     

Assessing the performance of the haplotype block model of linkage disequilibrium

Several recent studies have suggested that linkage disequilibrium (LD) in the human genome has a fundamentally "blocklike" structure. However, thus far there has been little formal assessment of how well the haplotype block model captures the underlying structure of LD. Here we propose quantitative criteria for assessing how blocklike LD is and apply these criteria to both real and simulated data. Analyses of several large data sets indicate that real data show a partial fit to the haplotype block model; some regions conform quite well, whereas others do not. Some improvement could be obtained by genotyping higher marker densities but not by increasing the number of samples. Nonetheless, although the real data are only moderately blocklike, our simulations indicate that, under a model of uniform recombination, the structure of LD would actually fit the block model much less well. Simulations of a model in which much of the recombination occurs in narrow hotspots provide a much better fit to the observed patterns of LD, suggesting that there is extensive fine-scale variation in recombination rates across the human genome.     

Demography,recombination hotspot intensity,and the block structure of linkage disequilibrium

BACKGROUND: Effective gene mapping based on genetic association data will require detailed knowledge of patterns of linkage disequilibrium (LD) in human populations. It has been recently suggested that linkage disequilibrium in humans may be organized in a block-like structure, with islands of high LD separated by regions of rapid breakdown of LD due to recombination hotspots. The experimental data to date, however, are limited, and fundamental questions remain about the implications of recombination rate heterogeneity. Here, we use computer simulations to evaluate how such heterogeneity influences patterns of LD, and we develop formal criteria to assess whether the patterns are functionally block like in the context of association mapping.RESULTS: Our analyses suggest that, even in models of extreme recombination rate heterogeneity, some human populations will have a functionally block-like structure to the pattern of LD, but others will not, depending on their precise demographic histories. In fact, for many models, we find that, following an LD-generating event, populations may move through discrete phases that can be functionally described as pre-block, block, and post-block. An analysis of observed and expected patterns of LD surrounding hotspots within the MHC Class II region confirms these theoretical expectations.CONCLUSIONS: Even if highly punctuated patterns of recombination are the rule, patterns of LD are still likely to show differences among populations and among genomic regions that are of practical importance in the design of genetic association studies. The notion that the average extent of LD is a useful concept for the design of association studies must be abandoned in light of the experimental and theoretical evidence.     

Distribution of genetic variation underlying adult migration timing in steelhead of the Columbia River basin

Fish migrations are energetically costly, especially when moving between freshwater and saltwater, but are a viable strategy for Pacific salmon and trout (Oncorhynchus spp.) due to the advantageous resources available at various life stages. Anadromous steelhead (O. mykiss) migrate vast distances and exhibit variation for adult migration phenotypes that have a genetic basis at candidate genes known as greb1L and rock1. We examined the distribution of genetic variation at 13 candidate markers spanning greb1L, intergenic, and rock1 regions versus 226 neutral markers for 113 populations (n = 9,471) of steelhead from inland and coastal lineages in the Columbia River. Patterns of population structure with neutral markers reflected genetic similarity by geographic region as demonstrated in previous studies, but candidate markers clustered populations by genetic variation associated with adult migration timing. Mature alleles for late migration had the highest frequency overall in steelhead populations throughout the Columbia River, with only 9 of 113 populations that had a higher frequency of premature alleles for early migration. While a single haplotype block was evident for the coastal lineage, we identified multiple haplotype blocks for the inland lineage. The inland lineage had one haplotype block that corresponded to candidate markers within the greb1L gene and immediately upstream in the intergenic region, and the second block only contained candidate markers from the intergenic region. Haplotype frequencies had similar patterns of geographic distribution as single markers, but there were distinct differences in frequency between the two haplotype blocks for the inland lineage. This may represent multiple recombination events that differed between lineages where phenotypic differences exist between freshwater entry versus arrival timing as indicated by Micheletti et al. (2018a). Redundancy analyses were used to model environmental effects on allelic frequencies of candidate markers, and significant variables were migration distance, temperature, isothermality, and annual precipitation. This study improves our understanding of the spatial distribution of genetic variation underlying adult migration timing in steelhead as well as associated environmental factors and has direct conservation and management implications.     

Haplotype block structure is conserved across mammals

Genetic variation in genomes is organized in haplotype blocks, and species-specific block structure is defined by differential contribution of population history effects in combination with mutation and recombination events. Haplotype maps characterize the common patterns of linkage disequilibrium in populations and have important applications in the design and interpretation of genetic experiments. Although evolutionary processes are known to drive the selection of individual polymorphisms, their effect on haplotype block structure dynamics has not been shown. Here, we present a high-resolution haplotype map for a 5-megabase genomic region in the rat and compare it with the orthologous human and mouse segments. Although the size and fine structure of haplotype blocks are species dependent, there is a significant interspecies overlap in structure and a tendency for blocks to encompass complete genes. Extending these findings to the complete human genome using haplotype map phase I data reveals that linkage disequilibrium values are significantly higher for equally spaced positions in genic regions, including promoters, as compared to intergenic regions, indicating that a selective mechanism exists to maintain combinations of alleles within potentially interacting coding and regulatory regions. Although this characteristic may complicate the identification of causal polymorphisms underlying phenotypic traits, conservation of haplotype structure may be employed for the identification and characterization of functionally important genomic regions.     

Comparison of Linkage Disequilibrium in Elite European Maize Inbred Lines using AFLP and SSR Markers

Application of association mapping to plant breeding populations has the potential to revolutionize plant genetics. The main objectives of this study were to (i) investigate the extent and genomic distribution of linkage disequilibrium (LD) between pairs of amplified fragment length polymorphism (AFLP) markers, (ii) compare these results with those obtained with simple sequence repeat (SSR) markers, and (iii) compare the usefulness of AFLP and SSR markers for genomewide association mapping in plant breeding populations. We examined LD in a cross-section of 72 European elite inbred lines genotyped with 452 AFLP and 93 SSR markers. LD was significant (p < 0.05) for about 15% of the AFLP marker pairs and for about 49% of the SSR marker pairs in each of the two germplasm groups, flint and dent. In both germplasm groups the ratio of linked to unlinked loci pairs in LD was higher for AFLPs than for SSRs. The observation of LD due to linkage for both marker types suggested that genome-wide association mapping should be possible using either AFLPs or SSRs. The results of our study indicated that SSRs should be favored over AFLPs but the opposite applies to populations with a long history of recombination.     

The impact of Yangtze River discharge, ocean currents and historical events on the biogeographic pattern of Cellana toreuma along the China coast

Aim

Genetic data were used to measure the phylogeographic distribution of the limpet, Cellana toreuma along the China coast in order to acsertain impacts of historic events, ocean currents and especially freshwater discharge from the Yangtze River on the connectivity of intertidal species with limited larval dispersal capability.

Methodology/Principal Findings

Genetic variation in 15 populations of C. toreuma (n = 418), ranging from the Yellow Sea (YS), East China Sea (ECS) and South China Sea (SCS), were determined from partial mitochondrial cytochrome c oxidase subunit I gene. Genetic diversity and divergence based on haplotype frequencies were analyzed using CONTRIB, and AMOVA was used to examine genetic population structure. Historic demographic expansions were evaluated from both neutrality tests and mismatch distribution tests. Among the 30 haplotypes identified, a dominant haplotype No. 1 (H1) existed in all the populations, and a relatively abundant private haplotype (H2) in YS. Pairwise F_ST values between YS and the other two groups were relatively high and the percentage of variation among groups was 10.9%.

Conclusions

The high nucleotide and gene diversity in the YS, with large pairwise genetic distances and relatively high percentages of variation among groups, suggests that this group was relatively isolated from ECS and SCS. This is likely driven by historic events, ocean currents, and demographic expansion. We propose that freshwater discharge from the Yangtze River, which may act as physical barrier limiting the southward dispersal of larvae from northern populations, is especially important in determining the separation of the YS group from the rest of the Chinese populations of C. toreuma.     

A nearest-neighboring-end algorithm for genetic mapping

MOTIVATION: High-throughput methods are beginning to make possible the genotyping of thousands of loci in thousands of individuals, which could be useful for tightly associating phenotypes to candidate loci. Current mapping algorithms cannot handle so many data without building hierarchies of framework maps. RESULTS: A version of Kruskal's minimum spanning tree algorithm can solve any genetic mapping problem that can be stated as marker deletion from a set of linkage groups. These include backcross, recombinant inbred, haploid and double-cross recombinational populations, in addition to conventional deletion and radiation hybrid populations. The algorithm progressively joins linkage groups at increasing recombination fractions between terminal markers, and attempts to recognize and correct erroneous joins at peaks in recombination fraction. The algorithm is O (mn3) for m individuals and n markers, but the mean run time scales close to mn2. It is amenable to parallel processing and has recovered true map order in simulations of large backcross, recombinant inbred and deletion populations with up to 37,005 markers. Simulations were used to investigate map accuracy in response to population size, allelic dominance, segregation distortion, missing data and random typing errors. It produced accurate maps when marker distribution was sufficiently uniform, although segregation distortion could induce translocated marker orders. The algorithm was also used to map 1003 loci in the F7 ITMI population of bread wheat, Triticum aestivum L. emend Thell., where it shortened an existing standard map by 16%, but it failed to associate blocks of markers properly across gaps within linkage groups. This was because it depends upon the rankings of recombination fractions at individual markers, and is susceptible to sampling error, typing error and joint selection involving the terminal markers of nearly finished linkage groups. Therefore, the current form of the algorithm is useful mainly to improve local marker ordering in linkage groups obtained in other ways. AVAILABILITY: The source code and supplemental data are http://www.iubio.bio.indiana.edu/soft/molbio/qtl/flipper/ CONTACT: ccrane@purdue.edu.     

Localization of the familial Mediterranean fever gene (FMF) to a 250-kb interval in non-Ashkenazi Jewish founder haplotypes. The French FMF Consortium.

Chromosome 16p13.3 harbors a gene (MEF) associated with familial Mediterranean fever (FMF), a recessive disease very common in populations of Mediterranean ancestry. In the course of positional cloning of MEF, we genotyped 26 non-Ashkenazi Jewish FMF pedigrees (310 meioses) with 15 microsatellite markers, most of which were recently developed by Généthon. Identification of recombination events in the haplotypes allowed narrowing of the MEF interval to a region between D16S3124 (telomeric) and D16S475 (centromeric). Two markers, D16S3070 and D16S3275, a microsatellite marker isolated from a YAC that also contains D16S3070, showed no recombination with the disease. Linkage disequilibrium and haplotype analysis highlighted the existence of a founder haplotype in our population. The core ancestral alleles were present in 71% of MEF-bearing chromosomes at loci D16S3070 and D16S3275. Furthermore, identification of historical crossing-over events in these pedigrees indicated that MEF is located between these two loci, which are both contained in a 250-kb genomic fragment.     

Analysis of genome-wide structure, diversity and fine mapping of Mendelian traits in traditional and village chickens

Extensive phenotypic variation is a common feature among village chickens found throughout much of the developing world, and in traditional chicken breeds that have been artificially selected for traits such as plumage variety. We present here an assessment of traditional and village chicken populations, for fine mapping of Mendelian traits using genome-wide single-nucleotide polymorphism (SNP) genotyping while providing information on their genetic structure and diversity. Bayesian clustering analysis reveals two main genetic backgrounds in traditional breeds, Kenyan, Ethiopian and Chilean village chickens. Analysis of linkage disequilibrium (LD) reveals useful LD (r(2) ≥ 0.3) in both traditional and village chickens at pairwise marker distances of ~10 Kb; while haplotype block analysis indicates a median block size of 11-12 Kb. Association mapping yielded refined mapping intervals for duplex comb (Gga 2:38.55-38.89 Mb) and rose comb (Gga 7:18.41-22.09 Mb) phenotypes in traditional breeds. Combined mapping information from traditional breeds and Chilean village chicken allows the oocyan phenotype to be fine mapped to two small regions (Gga 1:67.25-67.28 Mb, Gga 1:67.28-67.32 Mb) totalling ~75 Kb. Mapping the unmapped earlobe pigmentation phenotype supports previous findings that the trait is sex-linked and polygenic. A critical assessment of the number of SNPs required to map simple traits indicate that between 90 and 110K SNPs are required for full genome-wide analysis of haplotype block structure/ancestry, and for association mapping in both traditional and village chickens. Our results demonstrate the importance and uniqueness of phenotypic diversity and genetic structure of traditional chicken breeds for fine-scale mapping of Mendelian traits in the species, with village chicken populations providing further opportunities to enhance mapping resolutions.