Joshua M. Akey

Estimating African American Admixture Proportions by Use of Population-Specific Alleles

We analyzed the European genetic contribution to 10 populations of African descent in the United States (Maywood, Illinois; Detroit; New York; Philadelphia; Pittsburgh; Baltimore; Charleston, South Carolina; New Orleans; and Houston) and in Jamaica, using nine autosomal DNA markers. These markers either are population-specific or show frequency differences >45% between the parental populations and are thus especially informative for admixture. European genetic ancestry ranged from 6.8% (Jamaica) to 22.5% (New Orleans). The unique utility of these markers is reflected in the low variance associated with these admixture estimates (SEM 1.3%-2.7%). We also estimated the male and female European contribution to African Americans, on the basis of informative mtDNA (haplogroups H and L) and Y Alu polymorphic markers. Results indicate a sex-biased gene flow from Europeans, the male contribution being substantially greater than the female contribution. mtDNA haplogroups analysis shows no evidence of a significant maternal Amerindian contribution to any of the 10 populations. We detected significant nonrandom association between two markers located 22 cM apart (FY-null and AT3), most likely due to admixture linkage disequilibrium created in the interbreeding of the two parental populations. The strength of this association and the substantial genetic distance between FY and AT3 emphasize the importance of admixed populations as a useful resource for mapping traits with different prevalence in two parental populations.

Trans-acting regulatory variation in Saccharomyces cerevisiae and the role of transcription factors

Natural genetic variation can cause significant differences in gene expression, but little is known about the polymorphisms that affect gene regulation. We analyzed regulatory variation in a cross between laboratory and wild strains of Saccharomyces cerevisiae. Clustering and linkage analysis defined groups of coregulated genes and the loci involved in their regulation. Most expression differences mapped to trans-acting loci. Positional cloning and functional assays showed that polymorphisms in GPA1 and AMN1 affect expression of genes involved in pheromone response and daughter cell separation, respectively. We also asked whether particular classes of genes were more likely to contain trans-regulatory polymorphisms. Notably, transcription factors showed no enrichment, and trans-regulatory variation seems to be broadly dispersed across classes of genes with different molecular functions.

Population history and natural selection shape patterns of genetic variation in 132 genes.

Identifying regions of the human genome that have been targets of natural selection will provide important insights into human evolutionary history and may facilitate the identification of complex disease genes. Although the signature that natural selection imparts on DNA sequence variation is difficult to disentangle from the effects of neutral processes such as population demographic history, selective and demographic forces can be distinguished by analyzing multiple loci dispersed throughout the genome. We studied the molecular evolution of 132 genes by comprehensively resequencing them in 24 African-Americans and 23 European-Americans. We developed a rigorous computational approach for taking into account multiple hypothesis tests and demographic history and found that while many apparent selective events can instead be explained by demography, there is also strong evidence for positive or balancing selection at eight genes in the European-American population, but none in the African-American population. Our results suggest that the migration of modern humans out of Africa into new environments was accompanied by genetic adaptations to emergent selective forces. In addition, a region containing four contiguous genes on Chromosome 7 showed striking evidence of a recent selective sweep in European-Americans. More generally, our results have important implications for mapping genes underlying complex human diseases.

Haplotypes vs single marker linkage disequilibrium tests: what do we gain?

The genetic dissection of complex diseases represents a formidable challenge for modern human genetics. Recently, it has been suggested that linkage disequilibrium (LD) based methods will be a powerful approach for delineating complex disease genes. Most proposed LD test statistics search for association between a single marker and a putative trait locus. However, the power of a single marker association test may suffer because LD information contained in flanking markers is ignored. Intuitively, haplotypes (which can be regarded as a collection of ordered markers) may be more powerful than individual, unorganised markers. In this study, we derive the analytical tools based on standard chi-square statistics to directly investigate and compare the power between multilocus haplotypes and single marker LD tests. More specifically, novel formulas are obtained in order to calculate expected haplotype frequencies of unlimited size. This study demonstrates that the use of haplotypes can significantly improve the power and robustness of mapping disease genes. Additionally, we detail how the power of haplotype based association tests are affected by important population genetic parameters such as the genetic distance between markers and disease locus, mode of disease inheritance, age of trait causing mutation, frequency of associated marker allele, and level of initial LD. Finally, published data from the Hereditary Hemochromatosis disease region is used to illustrate the utility of haplotypes.

Distribution of recombination crossovers and the origin of haplotype blocks: The interplay of population history, recombination, and mutation

Recent studies suggest that haplotypes are arranged into discrete blocklike structures throughout the human genome. Here, we present an alternative haplotype block definition that assumes no recombination within each block but allows for recombination between blocks, and we use it to study the combined effects of demographic history and various population genetic parameters on haplotype block characteristics. Through extensive coalescent simulations and analysis of published haplotype data on chromosome 21, we find that (1) the combined effects of population demographic history, recombination, and mutation dictate haplotype block characteristics and (2) haplotype blocks can arise in the absence of recombination hot spots. Finally, we provide practical guidelines for designing and interpreting studies investigating haplotype block structure.

Y-Chromosome Evidence for a Northward Migration of Modern Humans into Eastern Asia during the Last Ice Age

The timing and nature of the arrival and the subsequent expansion of modern humans into eastern Asia remains controversial. Using Y-chromosome biallelic markers, we investigated the ancient human-migration patterns in eastern Asia. Our data indicate that southern populations in eastern Asia are much more polymorphic than northern populations, which have only a subset of the southern haplotypes. This pattern indicates that the first settlement of modern humans in eastern Asia occurred in mainland Southeast Asia during the last Ice Age, coinciding with the absence of human fossils in eastern Asia, 50,000-100,000 years ago. After the initial peopling, a great northward migration extended into northern China and Siberia.

The genomic distribution of population substructure in four populations using 8,525 autosomal SNPs

Understanding the nature of evolutionary relationships among persons and populations is important for the efficient application of genome science to biomedical research. We have analysed 8,525 autosomal single nucleotide polymorphisms (SNPs) in 84 individuals from four populations: African-American, European-American, Chinese and Japanese. Individual relationships were reconstructed using the allele sharing distance and the neighbour-joining tree making method. Trees show clear clustering according to population, with the root branching from the African-American clade. The African-American cluster is much less star-like than European-American and East Asian clusters, primarily because of admixture. Furthermore, on the East Asian branch, all ten Chinese individuals cluster together and all ten Japanese individuals cluster together. Using positional information, we demonstrate strong correlations between inter-marker distance and both locus-specific FST (the proportion of total variation due to differentiation) levels and branch lengths. Chromosomal maps of the distribution of locus-specific branch lengths were constructed by combining these data with other published SNP markers (total of 33,704 SNPs). These maps clearly illustrate a non-uniform distribution of human genetic substructure, an instructional and useful paradigm for education and research.

Telomere Length as a Quantitative Trait: Genome-Wide Survey and Genetic Mapping of Telomere Length-Control Genes in Yeast

Telomere length-variation in deletion strains of Saccharomyces cerevisiae was used to identify genes and pathways that regulate telomere length. We found 72 genes that when deleted confer short telomeres, and 80 genes that confer long telomeres relative to those of wild-type yeast. Among identified genes, 88 have not been previously implicated in telomere length control. Genes that regulate telomere length span a variety of functions that can be broadly separated into telomerase-dependent and telomerase-independent pathways. We also found 39 genes that have an important role in telomere maintenance or cell proliferation in the absence of telomerase, including genes that participate in deoxyribonucleotide biosynthesis, sister chromatid cohesion, and vacuolar protein sorting. Given the large number of loci identified, we investigated telomere lengths in 13 wild yeast strains and found substantial natural variation in telomere length among the isolates. Furthermore, we crossed a wild isolate to a laboratory strain and analyzed telomere length in 122 progeny. Genome-wide linkage analysis among these segregants revealed two loci that account for 30%–35% of telomere length-variation between the strains. These findings support a general model of telomere length-variation in outbred populations that results from polymorphisms at a large number of loci. Furthermore, our results laid the foundation for studying genetic determinants of telomere length-variation and their roles in human disease.

Large-scale SNP analysis reveals clustered and continuous patterns of human genetic variation

Understanding the distribution of human genetic variation is an important foundation for research into the genetics of common diseases. Some of the alleles that modify common disease risk are themselves likely to be common and, thus, amenable to identification using gene-association methods. A problem with this approach is that the large sample sizes required for sufficient statistical power to detect alleles with moderate effect make gene-association studies susceptible to false-positive findings as the result of population stratification [1, 2]. Such type I errors can be eliminated by using either family-based association tests or methods that sufficiently adjust for population stratification [3–5]. These methods require the availability of genetic markers that can detect and, thus, control for sources of genetic stratification among populations. In an effort to investigate population stratification and identify appropriate marker panels, we have analysed 11,555 single nucleotide polymorphisms in 203 individuals from 12 diverse human populations. Individuals in each population cluster to the exclusion of individuals from other populations using two clustering methods. Higher-order branching and clustering of the populations are consistent with the geographic origins of populations and with previously published genetic analyses. These data provide a valuable resource for the definition of marker panels to detect and control for population stratification in population-based gene identification studies. Using three US resident populations (European-American, African-American and Puerto Rican), we demonstrate how such studies can proceed, quantifying proportional ancestry levels and detecting significant admixture structure in each of these populations.

Interaction between the melanocortin-1 receptor and P genes contributes to inter-individual variation in skin pigmentation phenotypes in a Tibetan population

Abstract. The melanocortin-1 receptor (MC1R) andP gene product are two important components of the human pigmentary system that have been shown to be associated with red hair/fair skin and cause type II oculocutaneous albinism, respectively. However, their contribution to inter-individual variation at the population level is not well defined. To this end, we genotyped 3 single nucleotide polymorphisms (SNPs) in theMC1R gene (Arg67Gln, Gln163Arg, Val92Met) and 2 SNPs in theP gene (IVS13–15 and Gly780Gly) in 184 randomly ascertained Tibetan subjects, whose skin color was measured as a quantitative trait by reflective spectroscopy. Single locus analyses failed to demonstrate an association between any of the 5 SNPs and skin pigmentation. However, when an epistatic model was applied to the data, a significant gene-gene interaction was identified between Val92Met in MC1R and IVS13–15 in theP gene (F=2.43;P=0.0105). We also discuss the possible mechanisms of how gene interactions arise in signal transduction pathways.