Association mapping with single-feature polymorphisms

We develop methods for exploiting "single-feature polymorphism" data, generated by hybridizing genomic DNA to oligonucleotide expression arrays. Our methods enable the use of such data, which can be regarded as very high density, but imperfect, polymorphism data, for genomewide association or linkage disequilibrium mapping. We use a simulation-based power study to conclude that our methods should have good power for organisms like Arabidopsis thaliana, in which linkage disequilibrium is extensive, the reason being that the noisiness of single-feature polymorphism data is more than compensated for by their great number. Finally, we show how power depends on the accuracy with which single-feature polymorphisms are called.     

Detecting strong positive selection in the genome

New statistical tests have been developed in the past decade that enable us to infer evidence of recent strong positive selection from genome-wide data on single-nucleotide polymorphism and to localize the targets of selection in the genome. Based on these tests, past demographic events that led to distortions of the site-frequency spectrum of variation can be distinguished from selection, in particular if linkage disequilibrium is taken into account. These methods have been successfully applied to species from which complete sequence information and polymorphism data are available, including Drosophila melanogaster, humans, and several plant species. To make full use of the available data, however, the tests that were primarily designed for panmictic populations need to be extended to spatially structured populations.     

Relative roles of mutation and recombination in generating allelic polymorphism at an MHC class II locus in Peromyscus maniculatus

The MHC class II loci encoding cell surface antigens exhibit extremely high allelic polymorphism. There is considerable uncertainty in the literature over the relative roles of recombination and de novo mutation in generating this diversity. We studied class II sequence diversity and allelic polymorphism in two populations of Peromyscus maniculatus, which are among the most widespread and abundant mammals of North America. We find that intragenic recombination (or gene conversion) has been the predominant mode for the generation of allelic polymorphism in this species, with the amount of population recombination per base pair exceeding mutation by at least an order of magnitude during the history of the sample. Despite this, patchwork motifs of sites with high linkage disequilibrium are observed. This does not appear to be consistent with the much larger amount of recombination versus mutation in the history of the sample, unless the recombination rate is highly non-uniform over the sequence or selection maintains certain sites in linkage disequilibrium. We conclude that selection is most likely to be responsible for preserving sequence motifs in the presence of abundant recombination.     

Linkage disequilibrium in the neurofibromatosis 1 (NF1) region: implications for gene mapping.

To test the usefulness of linkage disequilibrium for gene mapping, we compared physical distances and linkage disequilibrium among eight RFLPs in the neurofibromatosis 1 (NF1) region. Seven of the polymorphisms span most of the NF1 gene, while the remaining polymorphism lies approximately 70 kb 3' to a stop codon in exon 49. By using Centre d'Etude du Polymorphisme Humain (CEPH) kindreds, 91-110 unrelated parents were genotyped. A high degree of disequilibrium is maintained among the seven intragenic polymorphisms (r > .82, P < 10(-7)), even though they are separated by as much as 340 kb. The 3' polymorphism is only 68 kb distal to the next polymorphism, but disequilibrium between the 3' polymorphism and all others is comparatively low (magnitude of 4 < .33, P values .27-.001). This result was replicated in three sets of unrelated kindreds: the Utah CEPH families, the non-Utah CEPH families, and an independent set of NF1 families. Trigenic, quadrigenic, three-locus, and four-locus disequilibrium measures were also estimated. There was little evidence of higher-order linkage disequilibrium. As expected for a disease with multiple mutations, no disequilibrium was observed between the disease gene and any of the RFLPs. The observed pattern of high disequilibrium within the gene and a loss of disequilibrium 3' to the stop codon could have implications for gene mapping studies. These are discussed, and guidelines for linkage disequilibrium studies are suggested.     

Molecular evolution of two linked genes, Est-6 and Sod, in Drosophila melanogaster.

We have obtained 15 sequences of Est-6 from a natural population of Drosophila melanogaster to test whether linkage disequilibrium exists between Est-6 and the closely linked Sod, and whether natural selection may be involved. An early experiment with allozymes had shown linkage disequilibrium between these two loci, while none was detected between other gene pairs. The Sod sequences for the same 15 haplotypes were obtained previously. The two genes exhibit similar levels of nucleotide polymorphism, but the patterns are different. In Est-6, there are nine amino acid replacement polymorphisms, one of which accounts for the S-F allozyme polymorphism. In Sod, there is only one replacement polymorphism, which corresponds to the S-F allozyme polymorphism. The transversion/transition ratio is more than five times larger in Sod than in Est-6. At the nucleotide level, the S and F alleles of Est-6 make up two allele families that are quite different from each other, while there is relatively little variation within each of them. There are also two families of alleles in Sod, one consisting of a subset of F alleles, and the other consisting of another subset of F alleles, designed F(A), plus all the S alleles. The Sod F(A) and S alleles are completely or nearly identical in nucleotide sequence, except for the replacement mutation that accounts for the allozyme difference. The two allele families have independent evolutionary histories in the two genes. There are traces of statistically significant linkage disequilibrium between the two genes that, we suggest, may have arisen as a consequence of selection favoring one particular sequence at each locus.     

Characterization of single nucleotide polymorphism markers for the green sea turtle (Chelonia mydas)

We present data on 29 new single nucleotide polymorphism assays for the green sea turtle, Chelonia mydas. DNA extracts from 39 green turtles were used for two methods of single nucleotide polymorphism discovery. The first approach employed an amplified fragment length polymorphism technique. The second technique screened a microsatellite library. Allele-specific amplification assays were developed for high-throughput single nucleotide polymorphism genotyping and tested on two Pacific C. mydas nesting populations. Observed heterozygosities ranged from 0 to 0.95 for a Hawaiian population and from 0 to 0.85 for a Galapagos population. Each of the populations had one locus out of Hardy-Weinberg equilibrium, SSCM2b and SSCM5 for Hawaii and Galapagos, respectively. No loci showed significant genotypic linkage disequilibrium across an expanded set of four Pacific nesting populations. However, two loci, SSCM4 and SSCM10b showed linkage disequilibrium across three populations indicating possible association.     

Hitchhiking and recombination in birds: evidence from Mhc-linked and unlinked loci in Red-winged Blackbirds (Agelaius phoeniceus)

Hitchhiking phenomena and genetic recombination have important consequences for a variety of fields for which birds are model species, yet we know virtually nothing about naturally occurring rates of recombination or the extent of linkage disequilibrium in birds. We took advantage of a previously sequenced cosmid clone from Red-winged Blackbirds (Agelaius phoeniceus) bearing a highly polymorphic Mhc class II gene, Agph-DABI, to measure the extent of linkage disequilibrium across approximately 40 kb of genomic DNA and to determine whether non-coding nucleotide diversity was elevated as a result of physical proximity to a target of balancing selection. Application of coalescent theory predicts that the hitchhiking effect is enhanced by the larger effective population size of blackbirds compared with humans, despite the presumably higher rates of recombination in birds. We surveyed sequence polymorphism at three Mhc-linked loci occurring 1.5-40 kb away from Agph-DAB1 and found that nucleotide diversity was indistinguishable from that found at three presumably unlinked, non-coding introns (beta-actin intron 2, beta-fibrinogen intron 7 and rhodopsin intron 2). Linkage disequilibrium as measured by Lewontin's D' was found only across a few hundred base pairs within any given locus, and was not detectable among any Mhc-linked loci. Estimated rates of the per site recombination rate p derived from three different analytical methods suggest that the amounts of recombination in blackbirds are up to two orders of magnitude higher than in humans, a discrepancy that cannot be explained entirely by the higher effective population size of blackbirds relative to humans. In addition, the ratio of the number of estimated recombination events per mutation frequently exceeds 1, as in Drosophila, again much higher than estimates in humans. Although the confidence limits of the blackbird estimates themselves span an order of magnitude, these data suggest that in blackbirds the hitchhiking effect for this region is negligible and may imply that the per site per individual recombination rate is high, resembling those of Drosophila more than those of humans.     

Genetic variation at the matrix metalloproteinase-9 locus on chromosome 20q12.2–13.1

Allelic association methods are better suited than linkage analysis for mapping of susceptibility genes that confer modest increases in risk in complex diseases. In both family- and population-based association studies, it is very useful to have prior knowledge of all sequence variants and the degree of linkage disequilibrium in a candidate gene region. In this study, we scanned sequence variants in a 2.2-kb promoter sequence and all 13 exons (totalling 3.3 kb) of the matrix metalloproteinase-9 gene, which is associated with coronary heart disease and a candidate for other diseases involving connective tissue remodelling, such as cancer metastasis. The sequences had a total of ten variable sites, four in the promoter, five in the coding region (three of which alter the amino acid encoded) and one in the 3' untranslated sequence. Sequence inspection suggests that some of the variants will have a functional impact on either level of expression or enzymatic activity. Tight linkage disequilibrium was detected between variants across the entire length of the gene (approximately 9 kb), and frequencies of different haplotypes were determined. The data provide an essential tool for studies of the possible contribution of genetic variation at the matrix metalloproteinase-9 locus to genetically determined susceptibility to a number of important diseases. The results also provide experimental data on the extent of linkage disequilibrium in the general population, which is yet to be resolved.     

Population structure and its effect on haplotype diversity and linkage disequilibrium surrounding the xa5 locus of rice (Oryza sativa L.)

To assess the usefulness of linkage disequilibrium mapping in an autogamous, domesticated species, we have characterized linkage disequilibrium in the candidate region for xa5, a recessive gene conferring race-specific resistance to bacterial blight in rice. This trait and locus have good mapping information, a tractable phenotype, and available sequence data, but no cloned gene. We sampled 13 short segments from the 70-kb candidate region in 114 accessions of Oryza sativa. Five additional segments were sequenced from the adjacent 45-kb region in resistant accessions to estimate the distance at which linkage disequilibrium decays. The data show significant linkage disequilibrium between sites 100 kb apart. The presence of the xa5 resistant reaction in two ecotypes and in accessions with different haplotypes in the candidate region may indicate multiple origins or genetic heterogeneity for resistance. In addition, genetic differentiation between ecotypes emphasizes the need for controlling for population structure in the design of linkage disequilibrium studies in rice.     

Effects of 5' regulatory-region polymorphisms on paraoxonase-gene (PON1) expression

Human HDL-associated paraoxonase (PON1) hydrolyzes a number of toxic organophosphorous compounds and reduces oxidation of LDLs and HDLs. These properties of PON1 account for its ability to protect against pesticide poisonings and atherosclerosis. PON1 also hydrolyzes a number of lactone and cyclic-carbonate drugs. Among individuals in a population, PON1 levels vary widely. We previously identified three polymorphisms in the PON1 regulatory region that affect expression levels in cultured human hepatocytes. In this study, we determined the genotypes of three regulatory-region polymorphisms for 376 white individuals and examined their effect on plasma-PON1 levels, determined by rates of phenylacetate hydrolysis. The -108 polymorphism had a significant effect on PON1-activity level, whereas the -162 polymorphism had a lesser effect. The -909 polymorphism, which is in linkage disequilibrium with the other sites, appears to have little or no independent effect on PON1-activity level in vivo. Other studies have found that the L55M polymorphism in the PON1-coding region is associated with differences in both PON1-mRNA and PON1-activity levels. The results presented here indicate that the L55M effect of lowered activity is not due to the amino acid change but is, rather, largely due to linkage disequilibrium with the -108 regulatory-region polymorphism. The codon 55 polymorphism marginally appeared to account for 15.3% of the variance in PON1 activity, but this dropped to 5% after adjustments for the effects of the -108 and Q192R polymorphisms were made. The -108C/T polymorphism accounted for 22.8% of the observed variability in PON1-expression levels, which was much greater than that attributable to the other PON1 polymorphisms. We also identified four sequence differences in the 3' UTR of the PON1 mRNA.     

Isolation and characterization of microsatellite loci from Macadamia

Thirty‐three microsatellite loci were isolated for the Australian rainforest tree Macadamia integrifolia. Genotyping across a test panel of 43 commercial cultivars generated an average polymorphic information content of 0.480. Five loci showed no polymorphism across cultivars. Significant linkage disequilibrium was detected in 10 pairwise comparisons, including two pairs of loci identified from the same clone sequence. The 33 microsatellite loci represent a significant tool for genome mapping and population genetic studies.     

Comparison of genetic models for analysing the effects of a PvuII polymorphism in the oestrogen receptor 1 (ESR1) gene on prolificacy in an Iberian x Meishan pig population

The effect of a previously reported PvuII polymorphism in oestrogen receptor 1 (ESR1) was analysed in an F(2) population of Iberian x Meishan pigs. We tested three hypotheses: (1) that a causal mutation was fixed in the parental populations, (2) that a causal mutation existed that was in complete linkage disequilibrium with the alleles of the PvuII polymorphism and (3) that a causal mutation existed in linkage disequilibrium within each parental population. The third model was the most plausible based on the available data. ESR1 alleles displayed different patterns of linkage disequilibrium with the causal mutation in each of the parental populations and the PvuII polymorphism was clearly not the causal mutation. As a consequence, the use of the ESR1 mutation for selection must be evaluated for a particular pig population before it is applied.     

Genetic variation in Native Americans, inferred from Latino SNP and resequencing data

Analyses of genetic polymorphism data have the potential to be highly informative about the demographic history of Native American populations, but due to a combination of historical and political factors, there are essentially no autosomal sequence polymorphism data from any Native American group. However, there are many resequencing studies involving Latinos, whose genomes contain segments inherited from their Native American ancestors. In this study, we introduce a new method for estimating local ancestry across the genomes of admixed individuals and show how this method, along with dense genotyping and targeted resequencing, can be used to assay genetic variation in ancestral Native American groups. We analyze roughly 6 Mb of resequencing data from 22 Mexican Americans to provide the first large-scale view of sequence level variation in Native Americans. We observe low levels of diversity and high levels of linkage disequilibrium in the Native American-derived sequences, consistent with a recent severe population bottleneck associated with the initial peopling of the Americas. Using two different computational approaches, one novel, we estimate that this bottleneck occurred roughly 12.5 Kya; when uncertainty in the estimation process is taken into account, our results are consistent with archeological estimates for the colonization of the Americas.     

Linkage disequilibrium and linkage analysis of the glucose-6-phosphatase gene

Recent studies have indicated that the four most common mutations account for 78% of mutant alleles in the glucose-6-phosphatase (G6Pase) gene. A significant fraction of mutant alleles remain unidentified. Thus, informative polymorphic markers are necessary for linkage analysis in carrier testing and prenatal diagnosis in families where mutations can not be identified. The common mutations appear to be ethnic-specific, suggesting that the individual mutations may have a common founder. With the recent discovery of the nucleotide 1176 polymorphism, we have studied whether these mutations are in linkage disequilibrium with the polymorphism. The results of polymerase chain reaction/allele-specific oligonucleotide analysis show that nucleotide 1176 C is in linkage disequilibrium with mutations R83 C and R83H, and with the splicing mutation 727G→T. The 1176 T polymorphism is in linkage disequilibrium with 459insTA. A GT repeat polymorphism has also been found. However, its heterozygosity is low. The 1176 nucleotide polymorphic marker can be used in carrier and prenatal diagnosis of GSD1a families that have unidentified mutations and are informative for this marker. Received: 27 January 1998 / Accepted: 17 April 1998     

PAH 399 GTA(Val)→GTT(Val), a new silent mutation found in the Chinese

Summary A silent mutation or sequence polymorphism, an A to T substitution at codon 399 in exon 11 of the phenylalanine hydroxylase (PAH) gene has been identified by DNA sequence analysis in the Chinese. The frequencies of this new mutation in normal and abnormal (phenylketonuria; PKU) genes are 0.005 and 0.09, respectively, based on the analyses of 100 apparently normal individuals and 39 PKU patients, as demonstrated by DNA amplification with polymerase chain reaction (PCR) and oligonucleotide hybridization methods. The results suggest that there is linkage disequilibrium between this polymorphism and PKU mutations in the PAH gene; approximately 10% of defect PAH alleles in the Chinese population may be identified with this sequence polymorphic marker.     

A haplotype-based 'haplotype relative risk' approach to detecting allelic associations.

A novel variation of the Haplotype Relative Risk (HRR) of Rubinstein et al. [Hum Immunol 1981;3:384] is proposed, in order to glean increased information about linkage disequilibrium or allelic associations by analyzing haplotype-based data rather than genotypic data. It is shown that statistical tests based on our design give much higher power than those based on the original HRR approach. Several additional nonparametric tests based on the same data are analyzed, and power is computed for each of them. Further, parametric likelihood methods are applied to testing linkage equilibrium, and estimating delta, the coefficient of linkage disequilibrium, from the same data.     

Linkage disequilibrium among commonly genotyped SNP variants detected from bull sequence,

Genomic prediction utilizing causal variants could increase selection accuracy above that achieved with SNPs genotyped by currently available arrays used for genomic selection. A number of variants detected from sequencing influential sires are likely to be causal, but noticeable improvements in prediction accuracy using imputed sequence variant genotypes have not been reported. Improvement in accuracy of predicted breeding values may be limited by the accuracy of imputed sequence variants. Using genotypes of SNPs on a high‐density array and non‐synonymous SNPs detected in sequence from influential sires of a multibreed population, results of this examination suggest that linkage disequilibrium between non‐synonymous and array SNPs may be insufficient for accurate imputation from the array to sequence. In contrast to 75% of array SNPs being strongly correlated to another SNP on the array, less than 25% of the non‐synonymous SNPs were strongly correlated to an array SNP. When correlations between non‐synonymous and array SNPs were strong, distances between the SNPs were greater than separation that might be expected based on linkage disequilibrium decay. Consistently near‐perfect whole‐genome linkage disequilibrium between the full array and each non‐synonymous SNP within the sequenced bulls suggests that whole‐genome approaches to infer sequence variants might be more accurate than imputation based on local haplotypes. Opportunity for strong linkage disequilibrium between sequence and array SNPs may be limited by discrepancies in allele frequency distributions, so investigating alternate genotyping approaches and panels providing greater chances of frequency‐matched SNPs strongly correlated to sequence variants is also warranted. Genotypes used for this study are available from https://www.animalgenome.org/repository/pub/ ;USDA2017.0519/.     

Linkage disequilibrium in the North American Holstein population

Linkage disequilibrium was estimated using 7119 single nucleotide polymorphism markers across the genome and 200 animals from the North American Holstein cattle population. The analysis of maternally inherited haplotypes revealed strong linkage disequilibrium ( r ²   >   0.8) in genomic regions of ∼50 kb or less. While linkage disequilibrium decays as a function of genomic distance, genomic regions within genes showed greater linkage disequilibrium and greater variation in linkage disequilibrium compared with intergenic regions. Identification of haplotype blocks could characterize the most common haplotypes. Although maximum haplotype block size was over 1 Mb, mean block size was 26–113 kb by various definitions, which was larger than that observed in humans (∼10 kb). Effective population size of the dairy cattle population was estimated from linkage disequilibrium between single nucleotide polymorphism marker pairs in various haplotype ranges. Rapid reduction of effective population size of dairy cattle was inferred from linkage disequilibrium in recent generations. This result implies a loss of genetic diversity because of the high rate of inbreeding and high selection intensity in dairy cattle. The pattern observed in this study indicated linkage disequilibrium in the current dairy cattle population could be exploited to refine mapping resolution. Changes in effective population size during past generations imply a necessity of plans to maintain polymorphism in the Holstein population.     

Linkage disequilibrium: what history has to tell us

Linkage disequilibrium has become important in the context of gene mapping. We argue that to understand the pattern of association between alleles at different loci, and of DNA sequence polymorphism in general, it is useful first to consider the underlying genealogy of the chromosomes. The stochastic process known as the coalescent is a convenient way to model such genealogies, and in this paper we set out the theory behind the coalescent and its implications for understanding linkage disequilibrium.