A human genome diversity cell line panel

A mutation in KIR3DS1 that results in truncation and lack of cell surface expression

The KIR gene cluster exhibits a high degree of polymorphism in terms of gene content as well as allelic polymorphism, and data suggest that it is evolving rapidly. The KIR3DL1 locus is one of the most polymorphic loci within this cluster and is unique in that it encodes an activating receptor KIR3DS1, as well as multiple inhibitory KIR3DL1 allotypes. Because KIR3DS1 has been implicated in a number of diseases, we tested for the presence of KIR3DS1 variants that might affect its expression and activating capacity. Preliminary FACS analysis indicated that indeed some individuals with the KIR3DS1 allele showed no cell surface expression of the molecule. Sequencing analysis identified a variant with a complex deletion/substitution mutation in exon 4 (which encodes the D1 extracellular domain), resulting in a premature stop codon. We subsequently genotyped 3,960 unrelated individuals and determined the frequencies of this allele across geographically distinct world populations. The data indicate that the null KIR3DS1 allele is uncommon, arose on a single haplotype, and spread across geographically distinct populations.

Variation in estimated recombination rates across human populations

Recently it has been reported that recombination hotspots appear to be highly variable between humans and chimpanzees, and there is evidence for between-person variability in hotspots, and evolutionary transience. To understand the nature of variation in human recombination rates, it is important to describe patterns of variability across populations. Direct measurement of recombination rates remains infeasible on a large scale, and population-genetic approaches can be imprecise, and are affected by demographic history. Reports to date have suggested broad similarity in recombination rates at large genomic scales and across human populations. Here, we examine recombination rate estimates at a finer population and genomic scale: 28 worldwide populations and 107 SNPs in a 1 Mb stretch of chromosome 22q. We employ analysis of variance of recombination rate estimates, corrected for differences in effective population size using genome-wide microsatellite mutation rate estimates. We find substantial variation in fine-scale rates between populations, but reduced variation within continental groups. All effects examined (SNP-pair, region, population and interactions) were highly significant. Adjustment for effective population size made little difference to the conclusions. Observed hotspots tended to be conserved across populations, albeit at varying intensities. This holds particularly for populations from the same region, and also to a considerable degree across geographical regions. However, some hotspots appear to be population-specific. Several results from studies on the population history of humans are in accordance with our analysis. Our results suggest that between-population variation in DNA sequences may underly recombination rate variation.

A compass that points to lupus: genetic studies on type I interferon pathway

It was more than 20 years ago that patients with systemic lupus erythematosus (SLE) were first reported to display elevated serum levels of type I interferon (IFN). Since then, extensive studies revealed a crucial role for type I IFN in SLE pathogenesis. The current model proposes that small increase of type I IFN production by plasmacytoid dendritic cells (pDCs) is sufficient to induce unabated activation of immature peripheral DCs. IFN-matured DCs select and activate autoreactive T cells and B cells, rather than deleting them, resulting in peripheral tolerance breakdown, a characteristic feature of SLE. Furthermore, immune complexes provide an amplification loop to pDCs for further IFN production. In the past 5 years, high-throughput technologies such as expression profiling and single-nucleotide polymorphism (SNP) typing established the role of altered type I IFN system in SLE, and a detailed picture of its molecular mechanisms is beginning to emerge. In this review, we discuss two major lines of genetics studies on type I IFN pathway related to human SLE: (1) expression profiling of IFN-responsive genes and (2) disease-associated SNPs of IFN-related genes, especially IRF5 (IFN-regulatory factor 5). Lastly, we discuss how such genetic alterations in type I IFN pathway fit in the current model of SLE pathogenesis.

Population stratification of a common APOBEC gene deletion polymorphism

The APOBEC3 gene family plays a role in innate cellular immunity inhibiting retroviral infection, hepatitis B virus propagation, and the retrotransposition of endogenous elements. We present a detailed sequence and population genetic analysis of a 29.5-kb common human deletion polymorphism that removes the APOBEC3B gene. We developed a PCR-based genotyping assay, characterized 1,277 human diversity samples, and found that the frequency of the deletion allele varies significantly among major continental groups (global FST = 0.2843). The deletion is rare in Africans and Europeans (frequency of 0.9% and 6%), more common in East Asians and Amerindians (36.9% and 57.7%), and almost fixed in Oceanic populations (92.9%). Despite a worldwide frequency of 22.5%, analysis of data from the International HapMap Project reveals that no single existing tag single nucleotide polymorphism may serve as a surrogate for the deletion variant, emphasizing that without careful analysis its phenotypic impact may be overlooked in association studies. Application of haplotype-based tests for selection revealed potential pitfalls in the direct application of existing methods to the analysis of genomic structural variation. These data emphasize the importance of directly genotyping structural variation in association studies and of accurately resolving variant breakpoints before proceeding with more detailed population-genetic analysis.

Perforin gene mutations in patients with acquired aplastic anemia

Perforin is a cytolytic protein expressed mainly in activated cytotoxic lymphocytes and natural killer cells. Inherited perforin mutations account for 20% to 40% of familial hemophagocytic lymphohistiocytosis, a fatal disease of early childhood characterized by the absence of functional perforin. Aplastic anemia, the paradigm of immune-mediated bone marrow failure syndromes, is characterized by hematopoietic stem cell destruction by activated T cells and Th1 cytokines. We examined whether mutations in the perforin gene occurred in acquired aplastic anemia. Three nonsynonymous PRF1 mutations among 5 unrelated patients were observed. Four of 5 patients with the mutations showed some hemophagocytosis in the bone marrow at diagnosis. Perforin protein levels in these patients were very low or absent, and perforin granules were completely absent. Natural killer (NK) cell cytotoxicity from these patients was significantly decreased. Our data suggest that PRF1 genetic alterations help explain the aberrant proliferation and activation of cytotoxic T cells and may represent genetic risk factors for bone marrow failure.

PediDraw: a web-based tool for drawing a pedigree in genetic counseling

Background

Drawing a pedigree is a prerequisite in genetic counseling. In medical records, a pedigree is useful to document the family history of the patient. Drawing a pedigree is also necessary in collecting genetic resources for medical research such as positional cloning. Currently, most pedigrees are drawn by hand or by drawing software. Due to the special requirements in a standardized pedigree, generating a pedigree by these methods is usually time-consuming and requires professionals. This limits the usage of a pedigree as demanded in remote diagnosis or online counseling from the counselees to send an electronic pedigree.

Results

We developed an online pedigree drawing tool, PediDraw, which enables users to generate pedigrees after inputting the family information step-by-step on web. It outputs a pedigree or table to present a family history to the counselors.

Conclusion

PediDraw is a user-friendly web-based drawing tool. It is accessible via Internet.

Using population mixtures to optimize the utility of genomic databases: linkage disequilibrium and association study design in India

When performing association studies in populations that have not been the focus of large-scale investigations of haplotype variation, it is often helpful to rely on genomic databases in other populations for study design and analysis - such as in the selection of tag SNPs and in the imputation of missing genotypes. One way of improving the use of these databases is to rely on a mixture of database samples that is similar to the population of interest, rather than using the single most similar database sample. We demonstrate the effectiveness of the mixture approach in the application of African, European, and East Asian HapMap samples for tag SNP selection in populations from India, a genetically intermediate region underrepresented in genomic studies of haplotype variation.

CYP2D6 worldwide genetic variation shows high frequency of altered activity variants and no continental structure

BACKGROUND AND OBJECTIVE

CYP2D6, a member of the cytochrome P450 superfamily, is responsible for the metabolism of about 25% of the commonly prescribed drugs. Its activity ranges from complete deficiency to excessive activity, potentially causing toxicity of medication or therapeutic failure with recommended drug dosages. This study aimed to describe the CYP2D6 diversity at the global level.

METHODS

A total of 1060 individuals belonging to 52 worldwide-distributed populations were genotyped at 12 highly informative variable sites, as well as for gene deletion and duplications. Phenotypes were predicted on the basis of haplotype combinations.

RESULTS AND CONCLUSIONS

Our study shows that (i) CYP2D6 diversity is far greater within than between populations and groups thereof, (ii) null or low-activity variants occur at high frequencies in various areas of the world, (iii) linkage disequilibrium is lowest in Africa and highest in the Americas. Patterns of variation, within and among populations, are similar to those observed for other autosomal markers (e.g. microsatellites and protein polymorphisms), suggesting that the diversity observed at the CYP2D6 locus reflects the same factors affecting variation at random genome markers.

Extreme individual marker F(ST )values do not imply population-specific selection in humans: the NRG1 example

Extreme population differentiation, as measured by the F(ST) value, has been suggested as an indicator of recent population-specific positive selection. Elevated F(ST) values indicating high differentiation between continental groups were previously reported on a linkage disequilibrium region in the Neuregulin 1 gene, a gene which has been associated to schizophrenia. In the present study we show evidence that high F(ST) values may not necessarily imply the action of selection, in particular positive selection, neither globally nor regionally, using the example of the NRG1 gene.

Association of CYP1B1 germ line mutations with hepatocyte nuclear factor 1alpha-mutated hepatocellular adenoma

Biallelic somatic mutations of TCF1 coding for hepatocyte nuclear factor 1alpha (HNF1alpha) are found in 50% of the hepatocellular adenoma (HCA) cases usually associated with oral contraception. In rare cases, HNF1alpha germ line mutations could also predispose to familial adenomatosis. In order to identify new genetic factors predisposing to HNF1alpha-mutated HCA, we searched for mutations in genes involved in the metabolism of estrogen. For 10 genes (CYP1A1, CYP1A2, CYP3A4, CYP3A5, COMT, UGT2B7, NQO1, GSTM1, GSTP1, and GSTT1), we did not find mutations nor differences in the allele distribution among 32 women presenting HNF1alpha-mutated adenomas compared with 58 controls. In contrast, we identified a CYP1B1 germ line heterozygous mutation in 4 of 32 women presenting HNF1alpha-mutated adenomas compared with none in 58 controls. We confirmed these results with the identification of four additional CYP1B1 mutations in a second series of 26 cases. No mutations were found in the control group, which was extended to 98 individuals, and only a known rare genetic variant was observed in two controls (P = 0.0003). We did an ethoxyresorufin O-deethylase assay to evaluate the functional consequence of the CYP1B1 mutations. We found reduced enzymatic activity in each CYP1B1 variant. In addition, an E229K CYP1B1 mutation was found in a woman with a germ line HNF1alpha mutation in a familial adenomatosis context. In this large family, all three patients with adenomatosis bore both HNF1 and CYP1B1 germ line mutations. In conclusion, our data suggested that CYP1B1 germ line-inactivating mutations might increase the incidence of HCA in women with HNF1alpha mutations.

A human derived SSADH coding variant is replacing the ancestral allele shared with primates

BACKGROUND

A growing number of reports describe markers with high frequencies of the ancestral alleles in Africa, contrasting with high frequencies and possibly fixation of derived variants out of Africa. Such a pattern can be explained by either neutral or non-neutral processes.

AIM

The study examined worldwide frequencies of two non-synonymous variants in NAD(+)-dependent succinic semialdehyde dehydrogenase (SSADH), in a search for possible signatures of natural selection favouring the derived alleles.

SUBJECTS AND METHODS

The typing of 1574 subjects were compiled, representing 60 populations from all continents. SSADH haplotype frequencies were correlated across 52 populations to those of 260 single nucleotide polymorphism (SNP) markers deposited in the CEPH database and of markers reported to be under positive Darwinian selection.

RESULTS

In the world population, the c.538C variant is proceeding to replace the ancestral c.538T, shared with primates. The overall population differentiation is within the normal range. A significant correlation was also found between the frequencies of the derived alleles in SSADH and Microcephalin (MCPH1), which showed concerted changes worldwide and, at least in Asian populations, also on a restricted geographical scale.

CONCLUSION

The analysis of robust correlations based on a large panel of populations is potentially able to identify clusters of genomic regions or genes showing co-evolution of the frequencies of derived alleles.

Two novel point mutations in the long-range SHH enhancer in three families with triphalangeal thumb and preaxial polydactyly

Spatio-temporal expression of sonic hedgehog (SHH) is driven by a regulatory element (ZRS) that lies 1 Mb upstream from SHH. Point mutations within the highly conserved ZRS have been described in the hemimelic extra toes mouse and in four families with preaxial polydactyly [Lettice et al., 2003]. Four North American Caucasian families were identified with autosomal dominant triphalangeal thumb. DNA from 20 affected and 36 unaffected family members was evaluated by sequence analysis of a 774-bp highly conserved ZRS contained within LMBR1 intron 5. Mutations within ZRS were identified in three of four families. In pedigree A and C, a novel A/G transition was identified near the 5' end of ZRS at bp 739 that segregated with disease or carrier status. Pedigree A, described previously [Dobbs et al., 2000], is a large family with 19 affected members who exhibit a milder phenotype, including predominantly triphalangeal thumbs and low penetrance (82%) relative to other families. Pedigree C is a small family with two affected family members with triphalangeal thumb, and one affected with both triphalangeal thumb and preaxial polydactyly. A novel C/G mutation at bp 621 was identified in pedigree B that segregated with the disease in all four affected individuals who manifested both preaxial polydactyly and triphalangeal thumb. Both mutations alter putative Cdx transcription factor binding sites. Mutations within ZRS appear to be a common cause of familial triphalangeal thumb and preaxial polydactyly. A genotype/phenotype correlate is suggested by pedigree A, whose mutation lies near the 5' end of ZRS; this family demonstrates a higher rate of nonpenetrance and milder phenotype. However, modifier genes may be contributing to the milder phenotype in this family.

Insights into modern disease from our distant evolutionary past

An EMBO workshop entitled 'Human Evolution and Disease' was held recently (6-9 December 2006, Hyderabad, India) where 141 scientists from many disciplines came together to discuss recent studies of human variation, origins and dispersal, natural selection and disease susceptibility. The meeting tackled the subject of human evolution and disease from the different perspectives of archaeology, linguistics, genetics and genomics based on both new and publicly available data sets. In this report, we highlight the latest fashion crazes in the discipline, in particular, the use of large public data sets and new methods to analyse modern human variation and the links between human evolution and disease susceptibility.

Impaired FGF signaling contributes to cleft lip and palate

Nonsyndromic cleft lip and palate (NS CLP) is a complex birth defect resulting from a combination of genetic and environmental factors. Several members of the FGF and FGFR families are expressed during craniofacial development and can rarely harbor mutations that result in human clefting syndromes. We hypothesized that disruptions in this pathway might also contribute to NS CLP. We sequenced the coding regions and performed association testing on 12 genes (FGFR1, FGFR2, FGFR3, FGF2, FGF3, FGF4, FGF7, FGF8, FGF9, FGF10, FGF18, and NUDT6) and used protein structure analyses to predict the function of amino acid variants. Seven likely disease-causing mutations were identified, including: one nonsense mutation (R609X) in FGFR1, a de novo missense mutation (D73H) in FGF8, and other missense variants in FGFR1, FGFR2, and FGFR3. Structural analysis of FGFR1, FGFR2, and FGF8 variants suggests that these mutations would impair the function of the proteins, albeit through different mechanisms. Genotyping of SNPs in the genes found associations between NS CLP and SNPs in FGF3, FGF7, FGF10, FGF18, and FGFR1. The data suggest that the FGF signaling pathway may contribute to as much as 3-5% of NS CLP and will be a consideration in the clinical management of CLP.

The population history of the Xibe in northern China: a comparison of autosomal, mtDNA and Y-chromosomal analyses of migration and gene flow

The Xibe population originated in northeastern China, but migrated to northwestern China in 1764-6. We have used a combination of published autosomal and Y-chromosomal data, together with newly derived mtDNA HVSI sequences, to investigate the extent to which genetic data can reveal the geographical origin of this population and the level of gene flow after the migration. We find that mtDNA data are uninformative, but that both autosomal and Y-chromosomal data indicate a northeastern origin, and that the Y data suggest 28% subsequent male-mediated gene flow into the population. We thus conclude that an appropriate combination of genetic data and analytical methods can reveal even recent and local events in the history of a population. In the Chinese samples examined, the combination of a northeastern autosomal background with a northwestern Y chromosome is indicative of a Xibe individual.

Generalized analysis of molecular variance

Many studies in the fields of genetic epidemiology and applied population genetics are predicated on, or require, an assessment of the genetic background diversity of the individuals chosen for study. A number of strategies have been developed for assessing genetic background diversity. These strategies typically focus on genotype data collected on the individuals in the study, based on a panel of DNA markers. However, many of these strategies are either rooted in cluster analysis techniques, and hence suffer from problems inherent to the assignment of the biological and statistical meaning to resulting clusters, or have formulations that do not permit easy and intuitive extensions. We describe a very general approach to the problem of assessing genetic background diversity that extends the analysis of molecular variance (AMOVA) strategy introduced by Excoffier and colleagues some time ago. As in the original AMOVA strategy, the proposed approach, termed generalized AMOVA (GAMOVA), requires a genetic similarity matrix constructed from the allelic profiles of individuals under study and/or allele frequency summaries of the populations from which the individuals have been sampled. The proposed strategy can be used to either estimate the fraction of genetic variation explained by grouping factors such as country of origin, race, or ethnicity, or to quantify the strength of the relationship of the observed genetic background variation to quantitative measures collected on the subjects, such as blood pressure levels or anthropometric measures. Since the formulation of our test statistic is rooted in multivariate linear models, sets of variables can be related to genetic background in multiple regression-like contexts. GAMOVA can also be used to complement graphical representations of genetic diversity such as tree diagrams (dendrograms) or heatmaps. We examine features, advantages, and power of the proposed procedure and showcase its flexibility by using it to analyze a wide variety of published data sets, including data from the Human Genome Diversity Project, classical anthropometry data collected by Howells, and the International HapMap Project.

Humans exhibit great genetic diversity. Understanding the factors that contribute to and sustain this diversity is an important research area. Not only can such understanding shed light on human origins, but it can also assist in the discovery of genes and genetic factors that contribute to debilitating diseases. Statistical analysis methods that can facilitate the identification of factors contributing to or associated with human genetic diversity are growing in number as new high-throughput molecular genetic assays and technologies are developed. We consider the use of an analysis method termed generalized analysis of molecular variance (GAMOVA), which builds off of previously proposed analysis methods for testing hypotheses about the factors associated with genetic background diversity. We apply the method in a wide variety of settings and show that it is both flexible and powerful. GAMOVA has great potential to assist in population-based human genetic studies, as it can be used to address questions such as: Is a sample of affected cases and unaffected controls from a homogeneous population, or is there evidence of heterogeneity that could affect the results of an association study? Is there reason to believe that the ancestry of a set of individuals influences the traits that they have?

PVRL1 variants contribute to non-syndromic cleft lip and palate in multiple populations

Poliovirus Receptor Like-1 (PVRL1) is a member of the immunoglobulin super family that acts in the initiation and maintenance of epithelial adherens junctions and is mutated in the cleft lip and palate/ectodermal dysplasia 1 syndrome (CLPED1, OMIM #225000). In addition, a common non-sense mutation in PVRL1 was discovered more often among non-syndromic sporadic clefting cases in Northern Venezuela in a previous case-control study. The present work sought to ascertain the role of PVRL1 in the sporadic forms of orofacial clefting in multiple populations. Multiple rare and common variants from all three splice isoforms were initially ascertained by sequencing 92 Iowan and 86 Filipino cases and CEPH controls. Using a family-based analysis to examine these variants, the common glycine allele of the G361V coding variant was significantly overtransmitted among all orofacial clefting phenotypes (P = 0.005). This represented G361V genotyping from over 800 Iowan, Danish, and Filipino families. Among four rare amino acid changes found within the V1 and C1 domains, S112T and T131A were found adjacent to critical amino acid positions within the V1 variable domain, regions previously shown to mediate cell-to-cell and cell-to-virus adhesion. The T131A variant was not found in over 1,300 non-affected control samples although the alanine is found in other species. The serine of the S112T variant position is conserved across all known PVRL1 sequences. Together these data suggest that both rare and common mutations within PVRL1 make a minor contribution to disrupting the initiation and regulation of cell-to-cell adhesion and downstream morphogenesis of the embryonic face.

The genetical bandwidth mapping: a spatial and graphical representation of population genetic structure based on the Wombling method

Characterizing the spatial variation of allele frequencies in a population has a wide range of applications in population genetics. This article introduces a new nonparametric method, which provides a two-dimensional representation of a structural parameter called the genetical bandwidth, which describes genetic structure around arbitrary spatial locations in a study area. This parameter corresponds to the shortest distance to areas of significant allele variation, and its computation is based on the Womble's systemic function. A simulation study and application to data sets taken from the literature give evidence that the method is particularly demonstrative when the fine-scale structure is stronger than the large-scale structure, and that it is generally able to locate genetic boundaries or clines precisely.

Worldwide analysis of multiple microsatellites: Language diversity has a detectable influence on DNA diversity

Previous studies of the correlations between the languages spoken by human populations and the genes carried by the members of those populations have been limited by the small amount of genetic markers available and by approximations in the treatment of linguistic data. In this study we analyzed a large collection of polymorphic microsatellite loci (377), distributed on all autosomes, and used Ruhlen's linguistic classification, to investigate the relative roles of geography and language in shaping the distribution of human DNA diversity at a worldwide scale. For this purpose, we performed three different kinds of analysis: (i) we partitioned genetic variances at three hierarchical levels of population subdivision according to language group by means of a molecular analysis of variance (AMOVA); (ii) we quantified by a series of Mantel's tests the correlation between measures of genetic and linguistic differentiation; and (iii) we tested whether linguistic differences are increased across known zones of increased genetic change between populations. Genetic differences appear to more closely reflect geographic than linguistic differentiation. However, our analyses show that language differences also have a detectable effect on DNA diversity at the genomic level, above and beyond the effects of geographic distance.

Human Genetic Diversity and the Epidemiology of Parasitic and Other Transmissible Diseases

This paper aims to review human genetic studies that are generally poorly known by parasitologists and scientists working on other pathogenic agents. The key proposals of this paper are as follows: (i) human susceptibility to transmissible diseases may often have a complex, multigenic background; (ii) recent discoveries indicate that major genomic rearrangements may be involved, possibly more so than DNA sequence; (iii) it is crucial to have a general population genetics framework of the human species based on neutral/historical markers to analyse reliably genetic susceptibility to infectious diseases; and (iv) the population level is a key factor. Ethnic diversity, a highly adaptive genetically driven phenotypic diversity, is possibly a valuable source for exploring human genetic susceptibility to transmissible diseases, since different populations have been exposed to drastically different geographic/climatic environments and different pathogens and vectors for tens of thousands of years. Studies dealing with human genetic susceptibility to transmissible diseases have mostly been based on the hypothesis that this factor is driven by only one or a few genes, and considered the individual more than the population level. Two different approaches have been developed for identifying the genes involved: (i) candidate genes and (ii) blind association studies (linkage analysis), screening the genome with a large number of high-resolution markers. Some loci involved in susceptibility to leishmaniosis, malaria and schistosomosis, for example, have already been identified. South American trypanosomosis (Chagas disease) is reviewed in detail to show the methodological problems of this classical approach. Current knowledge on the general impact of transmissible diseases on human genetic diversity, mainly HLA polymorphism, and the hopes raised by recent major international programmes such as the Human Genome Project (HGP), Human Genome Diversity Project (HGDP), International Human Haplotype Map Project (Hap Map) and extended databases, networks and networks of networks will also be reviewed.

Signatures of strong population differentiation shape extended haplotypes across the human CD28, CTLA4, and ICOS costimulatory genes

The three members of the costimulatory receptor family, CD28, CTLA-4, and ICOS, have complementary effects on T cell activation, and their balance controls the overall outcome of immune and autoimmune responses. They are encoded in a short genomic interval, and overall activity may result from interplay between allelic variants at each locus. With multiethnic DNA panels that represent a wide spectrum of human populations, we demonstrate long-range linkage disequilibrium among the three genes. A large fraction of the variation found in the locus can be explained by the presence of extended haplotypes encompassing variants at CD28, CTLA4, and the ICOS promoter. There are unusual differences in the distribution of some variants and haplotypes between geographic regions. The differences may reflect demographic events and/or the adaptation to diverse environmental and microbial challenges encountered in the course of human migrations and will be important to consider when interpreting association to immune/autoimmune responsiveness.

Structural variation on the short arm of the human Y chromosome: recurrent multigene deletions encompassing Amelogenin Y

Structural polymorphism is increasingly recognized as a major form of human genome variation, and is particularly prevalent on the Y chromosome. Assay of the Amelogenin Y gene (AMELY) on Yp is widely used in DNA-based sex testing, and sometimes reveals males who have interstitial deletions. In a collection of 45 deletion males from 12 populations, we used a combination of sequence-tagged site mapping, and binary-marker and Y-short tandem repeat haplotyping to understand the structural basis of this variation. Of the 45 deletion males, 41 carry indistinguishable deletions, 3.0-3.8 Mb in size. Breakpoint mapping strongly implicates a mechanism of non-allelic homologous recombination between the proximal major array of TSPY gene-containing repeats, and a single distal copy of TSPY; this is supported by the estimation of TSPY copy number in deleted and non-deleted males. The remaining four males carry three distinct non-recurrent deletions (2.5-4.0 Mb), which may be due to non-homologous mechanisms. Haplotyping shows that TSPY-mediated deletions have arisen seven times independently in the sample. One instance, represented by 30 chromosomes mostly of Indian origin within haplogroup J2e1*/M241, has a time-to-most-recent-common-ancestor of approximately 7700+/-1300 years. In addition to AMELY, deletion males all lack the genes PRKY and TBL1Y, and the rarer deletion classes also lack PCDH11Y. The persistence and expansion of deletion lineages, together with direct phenotypic evidence, suggests that absence of these genes has no major deleterious effects.

Analysis of IFT74 as a candidate gene for chromosome 9p-linked ALS-FTD

BACKGROUND

A new locus for amyotrophic lateral sclerosis--frontotemporal dementia (ALS-FTD) has recently been ascribed to chromosome 9p.

METHODS

We identified chromosome 9p segregating haplotypes within two families with ALS-FTD (F476 and F2) and undertook mutational screening of candidate genes within this locus.

RESULTS

Candidate gene sequencing at this locus revealed the presence of a disease segregating stop mutation (Q342X) in the intraflagellar transport 74 (IFT74) gene in family 476 (F476), but no mutation was detected within IFT74 in family 2 (F2). While neither family was sufficiently informative to definitively implicate or exclude IFT74 mutations as a cause of chromosome 9-linked ALS-FTD, the nature of the mutation observed within F476 (predicted to truncate the protein by 258 amino acids) led us to sequence the open reading frame of this gene in a large number of ALS and FTD cases (n = 420). An additional sequence variant (G58D) was found in a case of sporadic semantic dementia. I55L sequence variants were found in three other unrelated affected individuals, but this was also found in a single individual among 800 Human Diversity Gene Panel samples.

CONCLUSION

Confirmation of the pathogenicity of IFT74 sequence variants will require screening of other chromosome 9p-linked families.

Standardized subsets of the HGDP-CEPH Human Genome Diversity Cell Line Panel, accounting for atypical and duplicated samples and pairs of close relatives

The HGDP-CEPH Human Genome Diversity Cell Line Panel is a widely-used resource for studies of human genetic variation. Here, pairs of close relatives that have been included in the panel are identified. Together with information on atypical and duplicated samples, the inferred relative pairs suggest standardized subsets of the panel for use in future population-genetic studies.

Is urbanization scrambling the genetic structure of human populations? A case study

Recent population expansion and increased migration linked to urbanization are assumed to be eroding the genetic structure of human populations. We investigated change in population structure over three generations by analysing both demographic and mitochondrial DNA (mtDNA) data from a random sample of 2351 men from 22 Iranian populations. Potential changes in genetic diversity (theta) and genetic distance (F(ST)) over the last three generations were analysed by assigning mtDNA sequences to populations based on the individual's place of birth or that of their mother or grandmother. Despite the fact that several areas included cities of over one million inhabitants, we detected no change in genetic diversity, and only a small decrease in population structure, except in the capital city (Tehran), which was characterized by massive immigration, increased theta and a large decrease in F(ST) over time. Our results suggest that recent erosion of human population structure might not be as important as previously thought, except in some large conurbations, and this clearly has important implications for future sampling strategies.

A chromosomal rearrangement hotspot can be identified from population genetic variation and is coincident with a hotspot for allelic recombination

Insights into the origins of structural variation and the mutational mechanisms underlying genomic disorders would be greatly improved by a genomewide map of hotspots of nonallelic homologous recombination (NAHR). Moreover, our understanding of sequence variation within the duplicated sequences that are substrates for NAHR lags far behind that of sequence variation within the single-copy portion of the genome. Perhaps the best-characterized NAHR hotspot lies within the 24-kb-long Charcot-Marie-Tooth disease type 1A (CMT1A)-repeats (REPs) that sponsor deletions and duplications that cause peripheral neuropathies. We investigated structural and sequence diversity within the CMT1A-REPs, both within and between species. We discovered a high frequency of retroelement insertions, accelerated sequence evolution after duplication, extensive paralogous gene conversion, and a greater than twofold enrichment of SNPs in humans relative to the genome average. We identified an allelic recombination hotspot underlying the known NAHR hotspot, which suggests that the two processes are intimately related. Finally, we used our data to develop a novel method for inferring the location of an NAHR hotspot from sequence variation within segmental duplications and applied it to identify a putative NAHR hotspot within the LCR22 repeats that sponsor velocardiofacial syndrome deletions. We propose that a large-scale project to map sequence variation within segmental duplications would reveal a wealth of novel chromosomal-rearrangement hotspots.

Transferability of tag SNPs in genetic association studies in multiple populations

A general question for linkage disequilibrium-based association studies is how power to detect an association is compromised when tag SNPs are chosen from data in one population sample and then deployed in another sample. Specifically, it is important to know how well tags picked from the HapMap DNA samples capture the variation in other samples. To address this, we collected dense data uniformly across the four HapMap population samples and eleven other population samples. We picked tag SNPs using genotype data we collected in the HapMap samples and then evaluated the effective coverage of these tags in comparison to the entire set of common variants observed in the other samples. We simulated case-control association studies in the non-HapMap samples under a disease model of modest risk, and we observed little loss in power. These results demonstrate that the HapMap DNA samples can be used to select tags for genome-wide association studies in many samples around the world.

A worldwide survey of haplotype variation and linkage disequilibrium in the human genome

Recent genomic surveys have produced high-resolution haplotype information, but only in a small number of human populations. We report haplotype structure across 12 Mb of DNA sequence in 927 individuals representing 52 populations. The geographic distribution of haplotypes reflects human history, with a loss of haplotype diversity as distance increases from Africa. Although the extent of linkage disequilibrium (LD) varies markedly across populations, considerable sharing of haplotype structure exists, and inferred recombination hotspot locations generally match across groups. The four samples in the International HapMap Project contain the majority of common haplotypes found in most populations: averaging across populations, 83% of common 20-kb haplotypes in a population are also common in the most similar HapMap sample. Consequently, although the portability of tag SNPs based on the HapMap is reduced in low-LD Africans, the HapMap will be helpful for the design of genome-wide association mapping studies in nearly all human populations.

The genetic structure of human populations studied through short insertion-deletion polymorphisms

In a landmark study Rosenberg et al. (2002) analyzed human genome diversity with 377 microsatellites in the HGDP-CEPH Genome Diversity Panel and reported that the populations were structured into five geographical regions: America, Sub-Saharan Africa, East Asia, Oceania and a cluster composed of Europe, the Middle East and Central Asia. They also observed that the within-population component accounted for 93-95%, and that the among-regions portion was only 3.6%, of the total genetic variance. We have also studied the HGDP-CEPH Diversity Panel (1,064 individuals from 52 populations) with a set of 40 biallelic slow-evolving short insertion-deletion polymorphisms (indels). We confirmed the partition of worldwide diversity into five genetic clusters that correspond to major geographic regions. Using the indels we have also disclosed an among-regions component of genetic variance considerably larger (12.1%) than had been estimated using microsatellites. Our study demonstrates that a set of 40 well-chosen biallelic markers is sufficient for the characterization of human population structure at the global level.

Bayesian clustering using hidden Markov random fields in spatial population genetics

We introduce a new Bayesian clustering algorithm for studying population structure using individually geo-referenced multilocus data sets. The algorithm is based on the concept of hidden Markov random field, which models the spatial dependencies at the cluster membership level. We argue that (i) a Markov chain Monte Carlo procedure can implement the algorithm efficiently, (ii) it can detect significant geographical discontinuities in allele frequencies and regulate the number of clusters, (iii) it can check whether the clusters obtained without the use of spatial priors are robust to the hypothesis of discontinuous geographical variation in allele frequencies, and (iv) it can reduce the number of loci required to obtain accurate assignments. We illustrate and discuss the implementation issues with the Scandinavian brown bear and the human CEPH diversity panel data set.

Identifying the environmental factors that determine the genetic structure of populations

The study of population genetic structure is a fundamental problem in population biology because it helps us obtain a deeper understanding of the evolutionary process. One of the issues most assiduously studied in this context is the assessment of the relative importance of environmental factors (geographic distance, language, temperature, altitude, etc.) on the genetic structure of populations. The most widely used method to address this question is the multivariate Mantel test, a nonparametric method that calculates a correlation coefficient between a dependent matrix of pairwise population genetic distances and one or more independent matrices of environmental differences. Here we present a hierarchical Bayesian method that estimates F(ST) values for each local population and relates them to environmental factors using a generalized linear model. The method is demonstrated by applying it to two data sets, a data set for a population of the argan tree and a human data set comprising 51 populations distributed worldwide. We also carry out a simulation study to investigate the performance of the method and find that it can correctly identify the factors that play a role in the structuring of genetic diversity under a wide range of scenarios.

A geographically explicit genetic model of worldwide human-settlement history

Currently available genetic and archaeological evidence is generally interpreted as supportive of a recent single origin of modern humans in East Africa. However, this is where the near consensus on human settlement history ends, and considerable uncertainty clouds any more detailed aspect of human colonization history. Here, we present a dynamic genetic model of human settlement history coupled with explicit geographical distances from East Africa, the likely origin of modern humans. We search for the best-supported parameter space by fitting our analytical prediction to genetic data that are based on 52 human populations analyzed at 783 autosomal microsatellite markers. This framework allows us to jointly estimate the key parameters of the expansion of modern humans. Our best estimates suggest an initial expansion of modern humans approximately 56,000 years ago from a small founding population of approximately 1,000 effective individuals. Our model further points to high growth rates in newly colonized habitats. The general fit of the model with the data is excellent. This suggests that coupling analytical genetic models with explicit demography and geography provides a powerful tool for making inferences on human-settlement history.

Impact of human population history on distributions of individual-level genetic distance

Summaries of human genomic variation shed light on human evolution and provide a framework for biomedical research. Variation is often summarised in terms of one or a few statistics (eg F(ST) and gene diversity). Now that multilocus genotypes for hundreds of autosomal loci are available for thousands of individuals, new approaches are applicable. Recently, trees of individuals and other clustering approaches have demonstrated the power of an individual-focused analysis. We propose analysing the distributions of genetic distances between individuals. Each distribution, or common ancestry profile (CAP), is unique to an individual, and does not require a priori assignment of individuals to populations. Here, we consider a range of models of population history and, using coalescent simulation, reveal the potential insights gained from a set of CAPs. Information lies in the shapes of individual profiles--sometimes captured by variance of individual CAPs--and the variation across profiles. Analysis of short tandem repeat genotype data for over 1,000 individuals from 52 populations is consistent with dramatic differences in population histories across human groups.

Inference of ancestry: constructing hierarchical reference populations and assigning unknown individuals

The ability to infer personal genetic ancestry is being increasingly utilised in certain medical and forensic situations. Herein, the unsupervised Bayesian clustering algorithms structure, is employed to analyse 377 autosomal short tandem repeats typed on 1,056 individuals from the Centre d'Etude du Polymorphisme Humain Human Diversity Panel. Individuals of known geographical origin were hierarchically classified into a framework of increasingly homogeneous clusters to serve as reference populations into which individuals of unknown ancestry can be assigned. The groupings were characterised by the geographical affinities of cluster members and the accuracy of these procedures was verified using several genetic indices. Fine-scale substructure was detectable beyond the broad population level classifications that previously have been explored in this dataset. Metrics indicated that within certain lines, the strongest structuring signals were detected at the leaves of the hierarchy where lineage-specific groupings were identified. The accuracy of unknown assignment was assessed at each level of the hierarchy using a 'leave one out' strategy in which each individual was stripped of cluster membership and then re-assigned using the supervised Bayesian clustering algorithm implemented in GeneClass2. Although most clusters at all levels of resolution experienced highly accurate assignment, a decline was observed in the finer levels due to the mixed membership characteristics of some individuals. The parameters defined by this study allowed for assignment of unknown individuals to genetically defined clusters with measured likelihood. Shared ancestry data can then be inferred for the unknown individual.

Analysis of nucleotide diversity of NAT2 coding region reveals homogeneity across Native American populations and high intra-population diversity

N-acetyltransferase 2 (NAT2), an important enzyme in clinical pharmacology, metabolizes antibiotics such as isoniazid and sulfamethoxazole, and catalyzes the transformation of aromatic and heterocyclic amines from the environment and diet into carcinogenic intermediates. Polymorphisms in NAT2 account for variability in the acetylator phenotype and the pharmacokinetics of metabolized drugs. Native Americans, settled in rural areas and large cities of Latin America, are under-represented in pharmacogenetics studies; therefore, we sequenced the coding region of NAT2 in 456 chromosomes from 13 populations from the Americas, and two from Siberia, detecting nine substitutions and 11 haplotypes. Variants *4 (37%), *5B (23%) and *7B (24%) showed high frequencies. Average frequencies of fast, intermediate and slow acetylators across Native Americans were 18, 56 and 25%, respectively. NAT2 intra-population genetic diversity for Native Americans is higher than East Asians and similar to the rest of the world, and NAT2 variants are homogeneously distributed across native populations of the continent.

Genetic evidence for the role of loci at 19q13 in cleft lip and palate

BACKGROUND

Clefts of the lip and palate are common birth defects, affecting approximately 1 in 700 births worldwide. The aetiology of clefting is complex, with multiple genetic and environmental influences.

METHODS

Genotype based linkage disequilibrium analysis was conducted using the family based association test (FBAT) and the likelihood ratio test (LRT). We also carried out direct sequencing of the PVR and PVRL2 candidate genes based on their homology to PVRL1, a gene shown previously to cause Margarita Island clefting. Participants included 434 patients with cleft lip with or without cleft palate or cleft palate only and their mothers from eight countries in South America, 205 nuclear triads (father-mother-affected child) from Iowa, 541 nuclear triads from Denmark, and 100 patients with cleft lip and palate from the Philippines.

RESULTS

An allelic variant in the PVR gene showed statistically significant association with both South American and Iowa populations (p = 0.0007 and p = 0.0009, respectively). Direct sequencing of PVR and PVRL2 yielded 26 variants, including two rare amino acid changes, one in each gene, which were not seen in controls.

CONCLUSIONS

We found an association between a common variant in a gene at 19q and isolated clefting in two heterogeneous populations. However, it is unclear from our data if rare variants in PVR and PVRL2 are sufficient to cause clefting in isolation.

Homozygous mutations in LPIN2 are responsible for the syndrome of chronic recurrent multifocal osteomyelitis and congenital dyserythropoietic anaemia (Majeed syndrome)

BACKGROUND

Majeed syndrome is an autosomal recessive, autoinflammatory disorder characterised by chronic recurrent multifocal osteomyelitis and congenital dyserythropoietic anaemia. The objectives of this study were to map, identify, and characterise the Majeed syndrome causal gene and to speculate on its function and role in skin and bone inflammation.

METHODS

Six individuals with Majeed syndrome from two unrelated families were identified for this study. Homozygosity mapping and parametric linkage analysis were employed for the localisation of the gene responsible for Majeed syndrome. Direct sequencing was utilised for the identification of mutations within the genes contained in the region of linkage. Expression studies and in silico characterisation of the identified causal gene and its protein were carried out.

RESULTS

The phenotype of Majeed syndrome includes inflammation of the bone and skin, recurrent fevers, and dyserythropoietic anaemia. The clinical picture of the six affected individuals is briefly reviewed. The gene was mapped to a 5.5 cM interval (1.8 Mb) on chromosome 18p. Examination of genes in this interval led to the identification of homozygous mutations in LPIN2 in affected individuals from the two families. LPIN2 was found to be expressed in almost all tissues. The function of LPIN2 and its role in inflammation remains unknown.

CONCLUSIONS

We conclude that homozygous mutations in LPIN2 result in Majeed syndrome. Understanding the aberrant immune response in this condition will shed light on the aetiology of other inflammatory disorders of multifactorial aetiology including isolated chronic recurrent multifocal osteomyelitis, Sweet syndrome, and psoriasis.

Signatures of demographic history and natural selection in the human major histocompatibility complex Loci

Many lines of evidence show that several HLA loci have experienced balancing selection. However, distinguishing among demographic and selective explanations for patterns of variation observed with HLA genes remains a challenge. In this study we address this issue using data from a diverse set of human populations at six classical HLA loci and, employing a comparative genomics approach, contrast results for HLA loci to those for non-HLA markers. Using a variety of analytic methods, we confirm and extend evidence for selection acting on several HLA loci. We find that allele frequency distributions for four of the six HLA loci deviate from neutral expectations and show that this is unlikely to be explained solely by demographic factors. Other features of HLA variation are explained in part by demographic history, including decreased heterozygosity and increased LD for populations at greater distances from Africa and a similar apportionment of genetic variation for HLA loci compared to putatively neutral non-HLA loci. On the basis of contrasts among different HLA loci and between HLA and non-HLA loci, we conclude that HLA loci bear detectable signatures of both natural selection and demographic history.

Association between the transforming growth factor alpha gene and nonsyndromic oral clefts: a HuGE review

Transforming growth factor alpha (TGFA) is a well-characterized mammalian growth factor. Since the first report of an association between DNA sequence variants at the TGFA genetic locus and nonsyndromic oral clefts, 47 studies have been carried out, producing conflicting results. In this review, the author synthesizes findings from published reports on the association between the TGFA gene and clefting in humans. Bias, lack of statistical power, and genuine population diversity can explain the diverse results. In the aggregate, TGFA is probably a genetic modifier of clefting in humans, which is consistent with the oligogenic model suggested for nonsyndromic oral clefts.

Spread of an inactive form of caspase-12 in humans is due to recent positive selection

The human caspase-12 gene is polymorphic for the presence or absence of a stop codon, which results in the occurrence of both active (ancestral) and inactive (derived) forms of the gene in the population. It has been shown elsewhere that carriers of the inactive gene are more resistant to severe sepsis. We have now investigated whether the inactive form has spread because of neutral drift or positive selection. We determined its distribution in a worldwide sample of 52 populations and resequenced the gene in 77 individuals from the HapMap Yoruba, Han Chinese, and European populations. There is strong evidence of positive selection from low diversity, skewed allele-frequency spectra, and the predominance of a single haplotype. We suggest that the inactive form of the gene arose in Africa approximately 100-500 thousand years ago (KYA) and was initially neutral or almost neutral but that positive selection beginning approximately 60-100 KYA drove it to near fixation. We further propose that its selective advantage was sepsis resistance in populations that experienced more infectious diseases as population sizes and densities increased.

Complete sequence data support lack of balancing selection on PRNP in a natural Chinese population

The M129V mutation in the human prion protein gene (PRNP) is the primary site linked to susceptibility to prion diseases in humans. The heterozygous state of this allele has been proven to be more resistant to prion-related diseases such as Creutzfeldt-Jakob disease. Based on a study of the common genetic variations, it has been proposed that balancing selection has played a critical role in shaping the distribution of 129V in worldwide populations. Kreitman and Di-Rienzo [Kreitman and Di-Rienzo (2004)Trend Genet 20:300-304] challenged this hypothesis by pointing out that the exclusion of polymorphisms with low frequency may introduce an ascertainment bias and, in turn, lead to a wrong conclusion. By studying sequence variations in samples from the Human Genome Diversity Project, Soldevila et al. [Soldevila et al. (2005) Trends Genet 21:389-391] showed that this ascertainment bias does exist. We argue that the evidence presented by Soldevila et al. may be compromised by the samples tested since they were selected from a population that may be substructured. In this study, we re-evaluated the hypothesis of balancing selection in a natural Chinese population using a much longer segment encompassing the entire genomic region of the PRNP gene (15 kb). We showed that the pattern of genetic variation in PRNP is not consistent with the presence of balancing selection in this gene.

Y-chromosome analysis confirms highly sex-biased dispersal and suggests a low male effective population size in bonobos (Pan paniscus)

Dispersal is a rare event that is difficult to observe in slowly maturing, long-lived wild animal species such as the bonobo. In this study we used sex-linked (mitochondrial DNA sequence and Y-chromosome microsatellite) markers from the same set of individuals to estimate the magnitude of difference in effective dispersal between the sexes and to investigate the long-term demographic history of bonobos. We sampled 34 males from four distinct geographical areas across the bonobo distribution range. As predicted for a female-dispersing species, we found much higher levels of differentiation among local bonobo populations based upon Y-chromosomal than mtDNA genetic variation. Specifically, almost all of the Y-chromosomal variation distinguished populations, while nearly all of the mtDNA variation was shared between populations. Furthermore, genetic distance correlated with geographical distance for mtDNA but not for the Y chromosome. Female bonobos have a much higher migration rate and/or effective population size as compared to males, and the estimate for the mitochondrial TMRCA (time to most recent common ancestor) was approximately 10 times greater than the estimate for the Y chromosome (410,000 vs. 40,000-45,000). For humans the difference is merely a factor of two, suggesting a more stable demographic history in bonobos in comparison to humans.

Deciphering host migrations and origins by means of their microbes

Mitochondrial DNA and microsatellite sequences are powerful genetic markers for inferring the genealogy and the population genetic structure of animals but they have only limited resolution for organisms that display low genetic variability due to recent strong bottlenecks. An alternative source of data for deciphering migrations and origins in genetically uniform hosts can be provided by some of their microbes, if their evolutionary history correlates closely with that of the host. In this review, we first discuss how a variety of viruses, and the bacterium Helicobacter pylori, can be used as genetic tracers for one of the most intensively studied species, Homo sapiens. Then, we review statistical problems and limitations that affect the calculation of particular population genetic parameters for these microbes, such as mutation rates, with particular emphasis on the effects of recombination, selection and mode of transmission. Finally, we extend the discussion to other host-parasite systems and advocate the adoption of an integrative approach to both sampling and analysis.

The portability of tagSNPs across populations: a worldwide survey

In the search for common genetic variants that contribute to prevalent human diseases, patterns of linkage disequilibrium (LD) among linked markers should be considered when selecting SNPs. Genotyping efficiency can be increased by choosing tagging SNPs (tagSNPs) in LD with other SNPs. However, it remains to be seen whether tagSNPs defined in one population efficiently capture LD in other populations; that is, how portable tagSNPs are. Indeed, tagSNP portability is a challenge for the applicability of HapMap results. We analyzed 144 SNPs in a 1-Mb region of chromosome 22 in 1055 individuals from 38 worldwide populations, classified into seven continental groups. We measured tagSNP portability by choosing three reference populations (to approximate the three HapMap populations), defining tagSNPs, and applying them to other populations independently on the availability of information on the tagSNPs in the compared population. We found that tagSNPs are highly informative in other populations within each continental group. Moreover, tagSNPs defined in Europeans are often efficient for Middle Eastern and Central/South Asian populations. TagSNPs defined in the three reference populations are also efficient for more distant and differentiated populations (Oceania, Americas), in which the impact of their special demographic history on the genetic structure does not interfere with successfully detecting the most common haplotype variation. This high degree of portability lends promise to the search for disease association in different populations, once tagSNPs are defined in a few reference populations like those analyzed in the HapMap initiative.

Testing association between LRRK2 and Parkinson's disease and investigating linkage disequilibrium

BACKGROUND

We and others recently identified the gene underlying PARK8 linked Parkinson's disease (PD). This gene, LRRK2, contains mutations that cause an autosomal dominant PD, including a mutation, G2019S, which is the most common PD causing mutation identified to date. Common genetic variability in genes that contain PD causing mutations has previously been implicated as a risk factor for typical sporadic disease.

METHODS

We undertook a case-control association analysis of LRRK2 in two independent European PD cohorts using 31 tagging single nucleotide polymorphisms (tSNPs) and five potentially functional SNPs. To assess the structure of this locus in different populations, we have performed linkage disequilibrium (LD) analysis using these variants in a human diversity panel.

RESULTS

We show that common genetic variability in LRRK2 is not associated with risk for PD in the European populations studied here. We also show inter-population variability in the strength of LD across this locus.

CONCLUSIONS

To our knowledge this is the first comprehensive analysis of common variability within LRRK2 as a risk factor for PD.

Polarity and temporality of high-resolution y-chromosome distributions in India identify both indigenous and exogenous expansions and reveal minor genetic influence of Central Asian pastoralists

Although considerable cultural impact on social hierarchy and language in South Asia is attributable to the arrival of nomadic Central Asian pastoralists, genetic data (mitochondrial and Y chromosomal) have yielded dramatically conflicting inferences on the genetic origins of tribes and castes of South Asia. We sought to resolve this conflict, using high-resolution data on 69 informative Y-chromosome binary markers and 10 microsatellite markers from a large set of geographically, socially, and linguistically representative ethnic groups of South Asia. We found that the influence of Central Asia on the pre-existing gene pool was minor. The ages of accumulated microsatellite variation in the majority of Indian haplogroups exceed 10,000-15,000 years, which attests to the antiquity of regional differentiation. Therefore, our data do not support models that invoke a pronounced recent genetic input from Central Asia to explain the observed genetic variation in South Asia. R1a1 and R2 haplogroups indicate demographic complexity that is inconsistent with a recent single history. Associated microsatellite analyses of the high-frequency R1a1 haplogroup chromosomes indicate independent recent histories of the Indus Valley and the peninsular Indian region. Our data are also more consistent with a peninsular origin of Dravidian speakers than a source with proximity to the Indus and with significant genetic input resulting from demic diffusion associated with agriculture. Our results underscore the importance of marker ascertainment for distinguishing phylogenetic terminal branches from basal nodes when attributing ancestral composition and temporality to either indigenous or exogenous sources. Our reappraisal indicates that pre-Holocene and Holocene-era--not Indo-European--expansions have shaped the distinctive South Asian Y-chromosome landscape.

A microarray system for genotyping 150 single nucleotide polymorphisms in the coding region of human mitochondrial DNA

We established a genotyping system for a panel of 150 SNPs in the coding regions of mitochondrial DNA based on multiplex tag-array minisequencing. We show the feasibility of this system for simultaneous identification of individuals and prediction of the geographical origin of the mitochondrial DNA population lineage of the sample donors by genotyping the panel of SNPs in 265 samples representing nine different populations from Africa, Europe, and Asia. Nearly 40,000 genotypes were produced in the study, with an overall genotyping success rate of 95% and accuracy close to 100%. The gene diversity value of the panel of 150 SNPs was 0.991, compared to 0.995 for sequencing 500 nucleotides of the hypervariable regions I and II of mtDNA. For 17 individuals with identical sequences in the hypervariable regions of mtDNA, our panel of SNPs increased the power of discrimination. We observed 144 haplotypes that correspond to previously determined mitochondrial "haplogroups," and they allowed prediction of the origin of the maternal population lineage of 97% of the analyzed samples.