The role of selection in the evolution of human mitochondrial genomes

More reliable estimates of divergence times in Pan using complete mtDNA sequences and accounting for population structure

Here, we report the sequencing and analysis of eight complete mitochondrial genomes of chimpanzees (Pan troglodytes) from each of the three established subspecies (P. t. troglodytes, P. t. schweinfurthii and P. t. verus) and the proposed fourth subspecies (P. t. ellioti). Our population genetic analyses are consistent with neutral patterns of evolution that have been shaped by demography. The high levels of mtDNA diversity in western chimpanzees are unlike those seen at nuclear loci, which may reflect a demographic history of greater female to male effective population sizes possibly owing to the characteristics of the founding population. By using relaxed-clock methods, we have inferred a timetree of chimpanzee species and subspecies. The absolute divergence times vary based on the methods and calibration used, but relative divergence times show extensive uniformity. Overall, mtDNA produces consistently older times than those known from nuclear markers, a discrepancy that is reduced significantly by explicitly accounting for chimpanzee population structures in time estimation. Assuming the human-chimpanzee split to be between 7 and 5 Ma, chimpanzee time estimates are 2.1-1.5, 1.1-0.76 and 0.25-0.18 Ma for the chimpanzee/bonobo, western/(eastern + central) and eastern/central chimpanzee divergences, respectively.

Population expansion in the North African late Pleistocene signalled by mitochondrial DNA haplogroup U6

BACKGROUND

The archaeology of North Africa remains enigmatic, with questions of population continuity versus discontinuity taking centre-stage. Debates have focused on population transitions between the bearers of the Middle Palaeolithic Aterian industry and the later Upper Palaeolithic populations of the Maghreb, as well as between the late Pleistocene and Holocene.

RESULTS

Improved resolution of the mitochondrial DNA (mtDNA) haplogroup U6 phylogeny, by the screening of 39 new complete sequences, has enabled us to infer a signal of moderate population expansion using Bayesian coalescent methods. To ascertain the time for this expansion, we applied both a mutation rate accounting for purifying selection and one with an internal calibration based on four approximate archaeological dates: the settlement of the Canary Islands, the settlement of Sardinia and its internal population re-expansion, and the split between haplogroups U5 and U6 around the time of the first modern human settlement of the Near East.

CONCLUSIONS

A Bayesian skyline plot placed the main expansion in the time frame of the Late Pleistocene, around 20 ka, and spatial smoothing techniques suggested that the most probable geographic region for this demographic event was to the west of North Africa. A comparison with U6's European sister clade, U5, revealed a stronger population expansion at around this time in Europe. Also in contrast with U5, a weak signal of a recent population expansion in the last 5,000 years was observed in North Africa, pointing to a moderate impact of the late Neolithic on the local population size of the southern Mediterranean coast.

High-throughput sequencing of complete human mtDNA genomes from the Philippines

Because of the time and cost associated with Sanger sequencing of complete human mtDNA genomes, practically all evolutionary studies have screened samples first to define haplogroups and then either selected a few samples from each haplogroup, or many samples from a particular haplogroup of interest, for complete mtDNA genome sequencing. Such biased sampling precludes many analyses of interest. Here, we used high-throughput sequencing platforms to generate, rapidly and inexpensively, 109 complete mtDNA genome sequences from random samples of individuals from three Filipino groups, including one Negrito group, the Mamanwa. We obtained on average ∼55-fold coverage per sequence, with <1% missing data per sequence. Various analyses attest to the accuracy of the sequences, including comparison to sequences of the first hypervariable segment of the control region generated by Sanger sequencing; patterns of nucleotide substitution and the distribution of polymorphic sites across the genome; and the observed haplogroups. Bayesian skyline plots of population size change through time indicate similar patterns for all three Filipino groups, but sharply contrast with such plots previously constructed from biased sampling of complete mtDNA genomes, as well as with an artificially constructed sample of sequences that mimics the biased sampling. Our results clearly demonstrate that the high-throughput sequencing platforms are the methodology of choice for generating complete mtDNA genome sequences.

Genomic-scale capture and sequencing of endogenous DNA from feces

Genomic-level analyses of DNA from non-invasive sources would facilitate powerful conservation and evolutionary studies in natural populations of endangered and otherwise elusive species. However, the typical low quantity and poor quality of DNA that is extracted from non-invasive samples have generally precluded such work. Here we apply a modified DNA capture protocol that, when used in combination with massively-parallel sequencing technology, facilitates efficient and highly-accurate resequencing of megabases of specified nuclear genomic regions from fecal DNA samples. We validated our approach by comparing genetic variants identified from corresponding fecal and blood DNA samples of six western chimpanzees (Pan troglodytes verus) across more than 1.5 megabases of chromosome 21, chromosome X, and the complete mitochondrial genome. Our results suggest that it is now feasible to conduct genomic studies in natural populations for which constraints on invasive sampling have otherwise long been a barrier. The data we collected also provided an opportunity to examine western chimpanzee genetic diversity at unprecedented scale. Despite high mitochondrial genome diversity (π = 0.585%), western chimpanzees have a low ratio (0.42) of X chromosomal (π = 0.034%) to autosomal (chromosome 21 π = 0.081%) sequence diversity, a pattern that may reflect an unusual demographic history of this subspecies.

Genetic architecture of metabolic rate: environment specific epistasis between mitochondrial and nuclear genes in an insect

The extent to which mitochondrial DNA (mtDNA) variation is involved in adaptive evolutionary change is currently being reevaluated. In particular, emerging evidence suggests that mtDNA genes coevolve with the nuclear genes with which they interact to form the energy producing enzyme complexes in the mitochondria. This suggests that intergenomic epistasis between mitochondrial and nuclear genes may affect whole-organism metabolic phenotypes. Here, we use crossed combinations of mitochondrial and nuclear lineages of the seed beetle Callosobruchus maculatus and assay metabolic rate under two different temperature regimes. Metabolic rate was affected by an interaction between the mitochondrial and nuclear lineages and the temperature regime. Sequence data suggests that mitochondrial genetic variation has a role in determining the outcome of this interaction. Our genetic dissection of metabolic rate reveals a high level of complexity, encompassing genetic interactions over two genomes, and genotype × genotype × environment interactions. The evolutionary implications of these results are twofold. First, because metabolic rate is at the root of life histories, our results provide insights into the complexity of life-history evolution in general, and thermal adaptation in particular. Second, our results suggest a mechanism that could contribute to the maintenance of nonneutral mtDNA polymorphism.

Potential pitfalls in MitoChip detected tumor-specific somatic mutations: a call for caution when interpreting patient data

BACKGROUND

Several investigators have employed high throughput mitochondrial sequencing array (MitoChip) in clinical studies to search mtDNA for markers linked to cancers. In consequence, a host of somatic mtDNA mutations have been identified as linked to different types of cancers. However, closer examination of these data show that there are a number of potential pitfalls in the detection tumor-specific somatic mutations in clinical case studies, thus urging caution in the interpretation of mtDNA data to the patients. This study examined mitochondrial sequence variants demonstrated in cancer patients, and assessed the reliability of using detected patterns of polymorphisms in the early diagnosis of cancer.

METHODS

Published entire mitochondrial genomes from head and neck, adenoid cystic carcinoma, sessile serrated adenoma, and lung primary tumor from clinical patients were examined in a phylogenetic context and compared with known, naturally occurring mutations which characterize different populations.

RESULTS

The phylogenetic linkage analysis of whole arrays of mtDNA mutations from patient cancerous and non-cancerous tissue confirmed that artificial recombination events occurred in studies of head and neck, adenoid cystic carcinoma, sessile serrated adenoma, and lung primary tumor. Our phylogenetic analysis of these tumor and control leukocyte mtDNA haplotype sequences shows clear cut evidence of mixed ancestries found in single individuals.

CONCLUSIONS

Our study makes two prescriptions: both in the clinical situation and in research 1. more care should be taken in maintaining sample identity and 2. analysis should always be undertaken with respect to all the data available and within an evolutionary framework to eliminate artifacts and mix-ups.

The origins and evolution of genetic disease risk in modern humans

Patterns and risks of human disease have evolved. In this article, I review evidence regarding the importance of recent adaptive evolution, positive selection, and genomic conflicts in shaping the genetic and phenotypic architectures of polygenic human diseases. Strong recent selection in human populations can create and maintain genetically based disease risk primarily through three processes: increased scope for dysregulation from recent human adaptations, divergent optima generated by intraspecific genomic conflicts, and transient or stable deleterious by-products of positive selection caused by antagonistic pleiotropy, ultimately due to trade-offs at the levels of molecular genetics, development, and physiology. Human disease due to these processes appears to be concentrated in three sets of phenotypes: cognition and emotion, reproductive traits, and life-history traits related to long life-span. Diverse, convergent lines of evidence suggest that a small set of tissues whose pleiotropic patterns of gene function and expression are under especially strong selection-brain, placenta, testis, prostate, breast, and ovary-has mediated a considerable proportion of disease risk in modern humans.

Mitochondrial haplogroup N1a phylogeography, with implication to the origin of European farmers

BACKGROUND

Tracing the genetic origin of central European farmer N1a lineages can provide a unique opportunity to assess the patterns of the farming technology spread into central Europe in the human prehistory. Here, we have chosen twelve N1a samples from modern populations which are most similar with the farmer N1a types and performed the complete mitochondrial DNA genome sequencing analysis. To assess the genetic and phylogeographic relationship, we performed a detailed survey of modern published N1a types from Eurasian and African populations.

RESULTS

The geographic origin and expansion of farmer lineages related N1a subclades have been deduced from combined analysis of 19 complete sequences with 166 N1a haplotypes. The phylogeographic analysis revealed that the central European farmer lineages have originated from different sources: from eastern Europe, local central Europe, and from the Near East via southern Europe.

CONCLUSIONS

The results obtained emphasize that the arrival of central European farmer lineages did not occur via a single demic diffusion event from the Near East at the onset of the Neolithic spread of agriculture into Europe. Indeed these results indicate that the Neolithic transition process was more complex in central Europe and possibly the farmer N1a lineages were a result of a 'leapfrog' colonization process.

Maternally inherited susceptibility to cancer

Tumor microenvironment promotes mtDNA mutations. A number of these mutations will affect cell metabolism and increase cell survival. These mutations are positively selected and contribute to other tumor features, such as extracellular matrix remodeling and angiogenic processes, thus favoring metastases. Like somatic mutations, although with less marked effects, some mtDNA population polymorphisms will affect OXPHOS function, cell metabolism, and homeostasis. Thus, they could behave as inherited susceptibility factors for cancer. However, in addition to epidemiological evidence, other more direct clues are required. The cybrid approach can help to clarify this issue.

Large-scale mtDNA screening reveals a surprising matrilineal complexity in east Asia and its implications to the peopling of the region

In order to achieve a thorough coverage of the basal lineages in the Chinese matrilineal pool, we have sequenced the mitochondrial DNA (mtDNA) control region and partial coding region segments of 6,093 mtDNAs sampled from 84 populations across China. By comparing with the available complete mtDNA sequences, 194 of those mtDNAs could not be firmly assigned into the available haplogroups. Completely sequencing 51 representatives selected from these unclassified mtDNAs identified a number of novel lineages, including five novel basal haplogroups that directly emanate from the Eurasian founder nodes (M and N). No matrilineal contribution from the archaic hominid was observed. Subsequent analyses suggested that these newly identified basal lineages likely represent the genetic relics of modern humans initially peopling East Asia instead of being the results of gene flow from the neighboring regions. The observation that most of the newly recognized mtDNA lineages have already differentiated and show the highest genetic diversity in southern China provided additional evidence in support of the Southern Route peopling hypothesis of East Asians. Specifically, the enrichment of most of the basal lineages in southern China and their rather ancient ages in Late Pleistocene further suggested that this region was likely the genetic reservoir of modern humans after they entered East Asia.

Detecting heteroplasmy from high-throughput sequencing of complete human mitochondrial DNA genomes

Heteroplasmy, the existence of multiple mtDNA types within an individual, has been previously detected by using mostly indirect methods and focusing largely on just the hypervariable segments of the control region. Next-generation sequencing technologies should enable studies of heteroplasmy across the entire mtDNA genome at much higher resolution, because many independent reads are generated for each position. However, the higher error rate associated with these technologies must be taken into consideration to avoid false detection of heteroplasmy. We used simulations and phiX174 sequence data to design criteria for accurate detection of heteroplasmy with the Illumina Genome Analyzer platform, and we used artificial mixtures and replicate data to test and refine the criteria. We then applied these criteria to mtDNA sequence reads for 131 individuals from five Eurasian populations that had been generated via a parallel tagged approach. We identified 37 heteroplasmies at 10% frequency or higher at 34 sites in 32 individuals. The mutational spectrum does not differ between heteroplasmic mutations and polymorphisms in the same individuals, but the relative mutation rate at heteroplasmic mutations is significantly higher than that estimated for all mutable sites in the human mtDNA genome. Moreover, there is also a significant excess of nonsynonymous mutations observed among heteroplasmies, compared to polymorphism data from the same individuals. Both mutation-drift and negative selection influence the fate of heteroplasmies to determine the polymorphism spectrum in humans. With appropriate criteria for avoiding false positives due to sequencing errors, next-generation technologies can provide novel insights into genome-wide aspects of mtDNA heteroplasmy.

Whole-genome genetic diversity in a sample of Australians with deep Aboriginal ancestry

Australia was probably settled soon after modern humans left Africa, but details of this ancient migration are not well understood. Debate centers on whether the Pleistocene Sahul continent (composed of New Guinea, Australia, and Tasmania) was first settled by a single wave followed by regional divergence into Aboriginal Australian and New Guinean populations (common origin) or whether different parts of the continent were initially populated independently. Australia has been the subject of relatively few DNA studies even though understanding regional variation in genomic structure and diversity will be important if disease-association mapping methods are to be successfully evaluated and applied across populations. We report on a genome-wide investigation of Australian Aboriginal SNP diversity in a sample of participants from the Riverine region. The phylogenetic relationship of these Aboriginal Australians to a range of other global populations demonstrates a deep common origin with Papuan New Guineans and Melanesians, with little evidence of substantial later migration until the very recent arrival of European colonists. The study provides valuable and robust insights into an early and important phase of human colonization of the globe. A broader survey of Australia, including diverse geographic sample populations, will be required to fully appreciate the continent's unique population history and consequent genetic heritage, as well as the importance of both to the understanding of health issues.

Evidence for sub-haplogroup h5 of mitochondrial DNA as a risk factor for late onset Alzheimer's disease

BACKGROUND

Alzheimer's Disease (AD) is the most common neurodegenerative disease and the leading cause of dementia among senile subjects. It has been proposed that AD can be caused by defects in mitochondrial oxidative phosphorylation. Given the fundamental contribution of the mitochondrial genome (mtDNA) for the respiratory chain, there have been a number of studies investigating the association between mtDNA inherited variants and multifactorial diseases, however no general consensus has been reached yet on the correlation between mtDNA haplogroups and AD.

METHODOLOGY/PRINCIPAL FINDINGS

We applied for the first time a high resolution analysis (sequencing of displacement loop and restriction analysis of specific markers in the coding region of mtDNA) to investigate the possible association between mtDNA-inherited sequence variation and AD in 936 AD patients and 776 cognitively assessed normal controls from central and northern Italy. Among over 40 mtDNA sub-haplogroups analysed, we found that sub-haplogroup H5 is a risk factor for AD (OR=1.85, 95% CI:1.04-3.23) in particular for females (OR=2.19, 95% CI:1.06-4.51) and independently from the APOE genotype. Multivariate logistic regression revealed an interaction between H5 and age. When the whole sample is considered, the H5a subgroup of molecules, harboring the 4336 transition in the tRNAGln gene, already associated to AD in early studies, was about threefold more represented in AD patients than in controls (2.0% vs 0.8%; p=0.031), and it might account for the increased frequency of H5 in AD patients (4.2% vs 2.3%). The complete re-sequencing of the 56 mtDNAs belonging to H5 revealed that AD patients showed a trend towards a higher number (p=0.052) of sporadic mutations in tRNA and rRNA genes when compared with controls.

CONCLUSIONS

Our results indicate that high resolution analysis of inherited mtDNA sequence variation can help in identifying both ancient polymorphisms defining sub-haplogroups and the accumulation of sporadic mutations associated with complex traits such as AD.

Human genetic differentiation across the Strait of Gibraltar

BACKGROUND

The Strait of Gibraltar is a crucial area in the settlement history of modern humans because it represents a possible connection between Africa and Europe. So far, genetic data were inconclusive about the fact that this strait constitutes a barrier to gene flow, as previous results were highly variable depending on the genetic locus studied. The present study evaluates the impact of the Gibraltar region in reducing gene flow between populations from North-Western Africa and South-Western Europe, by comparing formally various genetic loci. First, we compute several statistics of population differentiation. Then, we use an original simulation approach in order to infer the most probable evolutionary scenario for the settlement of the area, taking into account the effects of both demography and natural selection at some loci.

RESULTS

We show that the genetic patterns observed today in the region of the Strait of Gibraltar may reflect an ancient population genetic structure which has not been completely erased by more recent events such as Neolithic migrations. Moreover, the differences observed among the loci (i.e. a strong genetic boundary revealed by the Y-chromosome polymorphism and, at the other extreme, no genetic differentiation revealed by HLA-DRB1 variation) across the strait suggest specific evolutionary histories like sex-mediated migration and natural selection. By considering a model of balancing selection for HLA-DRB1, we here estimate a coefficient of selection of 2.2% for this locus (although weaker in Europe than in Africa), which is in line with what was estimated from synonymous versus non-synonymous substitution rates. Selection at this marker thus appears strong enough to leave a signature not only at the DNA level, but also at the population level where drift and migration processes were certainly relevant.

CONCLUSIONS

Our multi-loci approach using both descriptive analyses and Bayesian inferences lead to better characterize the role of the Strait of Gibraltar in the evolution of modern humans. We show that gene flow across the Strait of Gibraltar occurred at relatively high rates since pre-Neolithic times and that natural selection and sex-bias migrations distorted the demographic signal at some specific loci of our genome.

Genetic heritage and native identity of the Seaconke Wampanoag tribe of Massachusetts

The name "Wampanoag" means "Eastern People" or "People of the First Light" in the local dialect of the Algonquian language. Once extensively populating the coastal lands and neighboring islands of the eastern United States, the Wampanoag people now consist of two federally recognized tribes, the Aquinnah and Mashpee, the state-recognized Seaconke Wampanoag tribe, and a number of bands and clans in present-day southern Massachusetts. Because of repeated epidemics and conflicts with English colonists, including King Philip's War of 1675-76, and subsequent colonial laws forbidding tribal identification, the Wampanoag population was largely decimated, decreasing in size from as many as 12,000 individuals in the 16th century to less than 400, as recorded in 1677. To investigate the influence of the historical past on its biological ancestry and native cultural identity, we analyzed genetic variation in the Seaconke Wampanoag tribe. Our results indicate that the majority of their mtDNA haplotypes belongs to West Eurasian and African lineages, thus reflecting the extent of their contacts and interactions with people of European and African descent. On the paternal side, Y-chromosome analysis identified a range of Native American, West Eurasian, and African haplogroups in the population, and also surprisingly revealed the presence of a paternal lineage that appears at its highest frequencies in New Guinea and Melanesia. Comparison of the genetic data with genealogical and historical information allows us to reconstruct the tribal history of the Seaconke Wampanoag back to at least the early 18th century.

Genetic diversity among ancient Nordic populations

Using established criteria for work with fossil DNA we have analysed mitochondrial DNA from 92 individuals from 18 locations in Denmark ranging in time from the Mesolithic to the Medieval Age. Unequivocal assignment of mtDNA haplotypes was possible for 56 of the ancient individuals; however, the success rate varied substantially between sites; the highest rates were obtained with untouched, freshly excavated material, whereas heavy handling, archeological preservation and storage for many years influenced the ability to obtain authentic endogenic DNA. While the nucleotide diversity at two locations was similar to that among extant Danes, the diversity at four sites was considerably higher. This supports previous observations for ancient Britons. The overall occurrence of haplogroups did not deviate from extant Scandinavians, however, haplogroup I was significantly more frequent among the ancient Danes (average 13%) than among extant Danes and Scandinavians (∼2.5%) as well as among other ancient population samples reported. Haplogroup I could therefore have been an ancient Southern Scandinavian type “diluted” by later immigration events. Interestingly, the two Neolithic samples (4,200 YBP, Bell Beaker culture) that were typed were haplogroup U4 and U5a, respectively, and the single Bronze Age sample (3,300–3,500 YBP) was haplogroup U4. These two haplogroups have been associated with the Mesolithic populations of Central and Northern Europe. Therefore, at least for Southern Scandinavia, our findings do not support a possible replacement of a haplogroup U dominated hunter-gatherer population by a more haplogroup diverse Neolithic Culture.

Analysis of Canis mitochondrial DNA demonstrates high concordance between the control region and ATPase genes

Background

Phylogenetic studies of wild Canis species have relied heavily on the mitochondrial DNA control region (mtDNA CR) to infer species relationships and evolutionary lineages. Previous analyses of the CR provided evidence for a North American evolved eastern wolf (C. lycaon), that is more closely related to red wolves (C. rufus) and coyotes (C. latrans) than grey wolves (C. lupus). Eastern wolf origins, however, continue to be questioned. Therefore, we analyzed mtDNA from 89 wolves and coyotes across North America and Eurasia at 347 base pairs (bp) of the CR and 1067 bp that included the ATPase6 and ATPase8 genes. Phylogenies and divergence estimates were used to clarify the evolutionary history of eastern wolves, and regional comparisons of nonsynonomous to synonomous substitutions (dN/dS) at the ATPase6 and ATPase8 genes were used to elucidate the potential role of selection in shaping mtDNA geographic distribution.

Results

We found high concordance across analyses between the mtDNA regions studied. Both had a high percentage of variable sites (CR = 14.6%; ATP = 9.7%) and both phylogenies clustered eastern wolf haplotypes monophyletically within a North American evolved lineage apart from coyotes. Divergence estimates suggest the putative red wolf sequence is more closely related to coyotes (D_xyCR = 0.01982 ± 0.00494 SD; D_xyATP = 0.00332 ± 0.00097 SD) than the eastern wolf sequences (D_xyCR = 0.03047 ± 0.00664 SD; D_xyATP = 0.00931 ± 0.00205 SD). Neutrality tests on both genes were indicative of the population expansion of coyotes across eastern North America, and dN/dS ratios suggest a possible role for purifying selection in the evolution of North American lineages. dN/dS ratios were higher in European evolved lineages from northern climates compared to North American evolved lineages from temperate regions, but these differences were not statistically significant.

Conclusions

These results demonstrate high concordance between coding and non-coding regions of mtDNA, and provide further evidence that the eastern wolf possessed distinct mtDNA lineages prior to recent coyote introgression. Purifying selection may have influenced North American evolved Canis lineages, but detection of adaptive selection in response to climate is limited by the power of current statistical tests. Increased sampling and development of alternative analytical tools will be necessary to disentangle demographic history from processes of natural selection.

Mitochondrial function, mitochondrial DNA and ageing: a reappraisal

The impressive performance of the research in biology of mitochondrion has greatly improved our knowledge on the functions of this organelle and highlighted the influence its functioning has on numerous human phenotypes. In particular, many studies have focused on the involvement of mitochondrion function (and dysfunction) in human ageing. To date, the literature in this specific field of mitochondrial biology is so vast that it is often difficult to properly put new data and new findings in the right context. The present paper aims to review the findings of the last few years in order to outline a general framework to understand how mitochondria can affect ageing and how ageing affects mitochondria.

Malaria parasite sequences from chimpanzee support the co-speciation hypothesis for the origin of virulent human malaria (Plasmodium falciparum)

Phylogenetic analyses of the mitochondrial cytochrome b (cytb), apicoplast caseinolytic protease C (clpC), and 18S rRNA sequences of Plasmodium isolates from chimpanzees along with those of the virulent human malaria parasite P. falciparum showed that the common chimpanzee (Pan troglodytes) malaria parasites, assigned by Rich et al. (2009) to P. reichenowi, constitute a paraphyletic assemblage. The assumption that P. falciparum diverged from P. reichenowi as recently as 5000-50,000 years ago would require a rate of synonymous substitution/site/year in cytb and clpC on the order of 10(-5)-10(-6), several orders of magnitude higher than any known from eukaryotic organelle genomes, and would imply an unrealistically recent timing of the most recent common ancestor of P. falciparum mitochondrial genomes. The available data are thus most consistent with the hypothesis that P. reichenowi (in the strict sense) and P. falciparum co-speciated with their hosts about 5-7 million years ago.

The archaeogenetics of Europe

A new timescale has recently been established for human mitochondrial DNA (mtDNA) lineages, making mtDNA at present the most informative genetic marker system for studying European prehistory. Here, we review the new chronology and compare mtDNA with Y-chromosome patterns, in order to summarize what we have learnt from archaeogenetics concerning five episodes over the past 50,000 years which significantly contributed to the settlement history of Europe: the pioneer colonisation of the Upper Palaeolithic, the Late Glacial re-colonisation of the continent from southern refugia after the Last Glacial Maximum, the postglacial re-colonization of deserted areas after the Younger Dryas cold snap, the arrival of Near Easterners with an incipient Neolithic package, and the small-scale migrations along continent-wide economic exchange networks beginning with the Copper Age. The available data from uniparental genetic systems have already transformed our view of the prehistory of Europe, but our knowledge of these processes remains limited. Nevertheless, their legacy remains as sedimentary layers in the gene pool of modern Europeans, and our understanding of them will improve substantially when more mtDNAs are completely sequenced, the Y chromosome more thoroughly analysed, and haplotype blocks of the autosomal genome become amenable to phylogeographic studies.

A method to correct for the effects of purifying selection on genealogical inference

Accurate reconstruction of the divergence times among individuals is an essential step toward inferring population parameters from genetic data. However, our ability to reconstruct accurate genealogies is often thwarted by the evolutionary forces we hope to detect, most prominently natural selection. Here, I demonstrate that purifying selection acting at many linked sites can systematically bias current methods of genealogical reconstruction, and I present a new method that corrects for this bias by allowing a class of sites to have a time-dependent rate. The parameters influencing the time dependency can be estimated from the data, allowing for a general method to detect the presence of selected sites and correcting for their distortion of the apparent mutation rate. The method works well under a variety of scenarios, including gamma-distributed selection coefficients as well as entirely neutral evolution. I also compare the performance of the new method to relaxed clock models, and I demonstrate the method on a data set from the mitochondrion of the North Atlantic whale-"louse" Cyamus ovalis.

The trans-Saharan slave trade - clues from interpolation analyses and high-resolution characterization of mitochondrial DNA lineages

Background

A proportion of 1/4 to 1/2 of North African female pool is made of typical sub-Saharan lineages, in higher frequencies as geographic proximity to sub-Saharan Africa increases. The Sahara was a strong geographical barrier against gene flow, at least since 5,000 years ago, when desertification affected a larger region, but the Arab trans-Saharan slave trade could have facilitate enormously this migration of lineages. Till now, the genetic consequences of these forced trans-Saharan movements of people have not been ascertained.

Results

The distribution of the main L haplogroups in North Africa clearly reflects the known trans-Saharan slave routes: West is dominated by L1b, L2b, L2c, L2d, L3b and L3d; the Center by L3e and some L3f and L3w; the East by L0a, L3h, L3i, L3x and, in common with the Center, L3f and L3w; while, L2a is almost everywhere. Ages for the haplogroups observed in both sides of the Saharan desert testify the recent origin (holocenic) of these haplogroups in sub-Saharan Africa, claiming a recent introduction in North Africa, further strengthened by the no detection of local expansions.

Conclusions

The interpolation analyses and complete sequencing of present mtDNA sub-Saharan lineages observed in North Africa support the genetic impact of recent trans-Saharan migrations, namely the slave trade initiated by the Arab conquest of North Africa in the seventh century. Sub-Saharan people did not leave traces in the North African maternal gene pool for the time of its settlement, some 40,000 years ago.

Using mitochondrial DNA to test the hypothesis of a European post-glacial human recolonization from the Franco-Cantabrian refuge

It has been proposed that the distribution patterns and coalescence ages found in Europeans for mitochondrial DNA (mtDNA) haplogroups V, H1 and H3 are the result of a post-glacial expansion from a Franco-Cantabrian refuge that recolonized central and northern areas. In contrast, in this refined mtDNA study of the Cantabrian Cornice that contributes 413 partial and 9 complete new mtDNA sequences, including a large Basque sample and a sample of Asturians, no experimental evidence was found to support the human refuge-expansion theory. In fact, all measures of gene diversity point to the Cantabrian Cornice in general and the Basques in particular, as less polymorphic for V, H1 and H3 than other southern regions in Iberia or in Central Europe. Genetic distances show the Cantabrian Cornice is a very heterogeneous region with significant local differences. The analysis of several minor subhaplogroups, based on complete sequences, also suggests different focal expansions over a local and peninsular range that did not affect continental Europe. Furthermore, all detected clinal trends show stronger longitudinal than latitudinal profiles. In Northern Iberia, it seems that the highest diversity values for some haplogroups with Mesolithic coalescence ages are centred on the Mediterranean side, including Catalonia and South-eastern France.

Evolution of the nuclear-encoded cytochrome oxidase subunits in vertebrates

Vertebrate mitochondrial cytochrome c oxidase (COX) possesses 10 nuclear-encoded subunits. Six subunits have paralogs in mammals, but the origins and distribution of isoforms among vertebrates have not been analyzed. We used Bayesian phylogenetic analysis to interpret the origins of each subunit, inferring the roles of gene and genome duplications. The paralogous ancestries of five genes were identical throughout the major vertebrate taxa: no paralogs of COX6c and COX7c, two paralogs of COX4 and COX6a, and three paralogs of COX7a. Two genes had an extra copy in teleosts (COX5a, COX5b), and three genes had additional copies in mammals (COX6b, COX7b, COX8). Focusing on early vertebrates, we examined structural divergence and explored transcriptional profiles across zebrafish tissues. Quantitative transcript profiles revealed dramatic differences in transcript abundance for different subunits. COX7b and COX4 transcripts were typically present at very low levels, whereas COX5a and COX8 were in vast excess in all tissues. For genes with paralogs, two general patterns emerged. For COX5a and COX8, there was ubiquitous expression of one paralog, with the other paralog in lower abundance in all tissues. COX4 and COX6a shared a distinct expression pattern, with one paralog dominant in brain and gills and the other in muscles. The isoform profiles in combination with phylogenetic analyses show that vertebrate COX isoform patterns are consistent with the hypothesis that early whole genome duplications in basal vertebrates governed the isoform repertoire in modern fish and tetrapods, though more recent lineage-specific gene/genome duplications also play a role in select subunits.

eCOMPAGT integrates mtDNA: import, validation and export of mitochondrial DNA profiles for population genetics, tumour dynamics and genotype-phenotype association studies

Background

Mitochondrial DNA (mtDNA) is widely being used for population genetics, forensic DNA fingerprinting and clinical disease association studies. The recent past has uncovered severe problems with mtDNA genotyping, not only due to the genotyping method itself, but mainly to the post-lab transcription, storage and report of mtDNA genotypes.

Description

eCOMPAGT, a system to store, administer and connect phenotype data to all kinds of genotype data is now enhanced by the possibility of storing mtDNA profiles and allowing their validation, linking to phenotypes and export as numerous formats. mtDNA profiles can be imported from different sequence evaluation programs, compared between evaluations and their haplogroup affiliations stored. Furthermore, eCOMPAGT has been improved in its sophisticated transparency (support of MySQL and Oracle), security aspects (by using database technology) and the option to import, manage and store genotypes derived from various genotyping methods (SNPlex, TaqMan, and STRs). It is a software solution designed for project management, laboratory work and the evaluation process all-in-one.

Conclusions

The extended mtDNA version of eCOMPAGT was designed to enable error-free post-laboratory data handling of human mtDNA profiles. This software is suited for small to medium-sized human genetic, forensic and clinical genetic laboratories. The direct support of MySQL and the improved database security options render eCOMPAGT a powerful tool to build an automated workflow architecture for several genotyping methods. eCOMPAGT is freely available at http://dbis-informatik.uibk.ac.at/ecompagt.

Compensatory evolution in mitochondrial tRNAs navigates valleys of low fitness

A long-standing controversy in evolutionary biology is whether or not evolving lineages can cross valleys on the fitness landscape that correspond to low-fitness genotypes, which can eventually enable them to reach isolated fitness peaks. Here we study the fitness landscapes traversed by switches between different AU and GC Watson-Crick nucleotide pairs at complementary sites of mitochondrial transfer RNA stem regions in 83 mammalian species. We find that such Watson-Crick switches occur 30-40 times more slowly than pairs of neutral substitutions, and that alleles corresponding to GU and AC non-Watson-Crick intermediate states segregate within human populations at low frequencies, similar to those of non-synonymous alleles. Substitutions leading to a Watson-Crick switch are strongly correlated, especially in mitochondrial tRNAs encoded on the GT-nucleotide-rich strand of the mitochondrial genome. Using these data we estimate that a typical Watson-Crick switch involves crossing a fitness valley of a depth of about 10(-3) or even about 10(-2), with AC intermediates being slightly more deleterious than GU intermediates. This compensatory evolution must proceed through rare intermediate variants that never reach fixation. The ubiquitous nature of compensatory evolution in mammalian mitochondrial tRNAs and other molecules implies that simultaneous fixation of two alleles that are individually deleterious may be a common phenomenon at the molecular level.

SNaPshot minisequencing to resolve mitochondrial macro-haplogroups found in Africa

African mitochondrial DNA (mtDNA) haplogroups are divided into seven macro-haplogroups (L0'1'2'3'4'5'6), while the rest of the world's lineages are classified as subgroups of macro-haplogroups M, N and R. The most common approach to characterizing mtDNA variation is the sequencing of hypervariable segments I and II of the non-coding control region of the molecule. Given the higher mutation rate within the control region compared with the coding regions of the molecule, recurrent mutations in the former can sometimes hide possible phylogenetic structure. The incorporation of haplogroup-defining coding region mutations has helped in overcoming this limitation. By judiciously selecting 14 coding region SNPs and incorporating them into a multiplex minisequencing assay we were able to resolve mtDNA sequences from some sub-Saharan African populations into ten macro-haplogroups (L0-L6, M, N and R). We tested the efficacy of the panel by screening 699 individuals, consisting mostly of Khoe-San, Bantu speakers and individuals with mixed ancestries (Coloreds) and found no inconsistencies compared with hypervariable segment sequencing results. The panel provided a fast and efficient means of classifying mtDNA into the ten mitochondrial macro-haplogroups and provided a reliable screening to distinguish African from non-African-derived mtDNA lineages.

Explaining the imperfection of the molecular clock of hominid mitochondria

The molecular clock of mitochondrial DNA has been extensively used to date various genetic events. However, its substitution rate among humans appears to be higher than rates inferred from human-chimpanzee comparisons, limiting the potential of interspecies clock calibrations for intraspecific dating. It is not well understood how and why the substitution rate accelerates. We have analyzed a phylogenetic tree of 3057 publicly available human mitochondrial DNA coding region sequences for changes in the ratios of mutations belonging to different functional classes. The proportion of non-synonymous and RNA genes substitutions has reduced over hundreds of thousands of years. The highest mutation ratios corresponding to fast acceleration in the apparent substitution rate of the coding sequence have occurred after the end of the Last Ice Age. We recalibrate the molecular clock of human mtDNA as 7990 years per synonymous mutation over the mitochondrial genome. However, the distribution of substitutions at synonymous sites in human data significantly departs from a model assuming a single rate parameter and implies at least 3 different subclasses of sites. Neutral model with 3 synonymous substitution rates can explain most, if not all, of the apparent molecular clock difference between the intra- and interspecies levels. Our findings imply the sluggishness of purifying selection in removing the slightly deleterious mutations from the human as well as the Neandertal and chimpanzee populations. However, for humans, the weakness of purifying selection has been further exacerbated by the population expansions associated with the out-of Africa migration and the end of the Last Ice Age.

Complete mitochondrial DNA sequences provide new insights into the Polynesian motif and the peopling of Madagascar

More than a decade of mitochondrial DNA (mtDNA) studies have given the 'Polynesian motif' renowned status as a marker for tracing the late-Holocene expansion of Austronesian speaking populations. Despite considerable research on the Polynesian motif in Oceania, there has been little equivalent work on the western edge of its expansion - leaving major issues unresolved regarding the motif's evolutionary history. This has also led to considerable uncertainty regarding the settlement of Madagascar. In this study, we assess mtDNA variation in 266 individuals from three Malagasy ethnic groups: the Mikea, Vezo, and Merina. Complete mtDNA genome sequencing reveals a new variant of the Polynesian motif in Madagascar; two coding region mutations define a Malagasy-specific sub-branch. This newly defined 'Malagasy motif' occurs at high frequency in all three ethnic groups (13-50%), and its phylogenetic position, geographic distribution, and estimated age all support a recent origin, but without conclusively identifying a specific source region. Nevertheless, the haplotype's limited diversity, similar to those of other mtDNA haplogroups found in our Malagasy groups, best supports a small number of initial settlers arriving to Madagascar through the same migratory process. Finally, the discovery of this lineage provides a set of new polymorphic positions to help localize the Austronesian ancestors of the Malagasy, as well as uncover the origin and evolution of the Polynesian motif itself.

ReseqChip: automated integration of multiple local context probe data from the MitoChip array in mitochondrial DNA sequence assembly

Background

The Affymetrix MitoChip v2.0 is an oligonucleotide tiling array for the resequencing of the human mitochondrial (mt) genome. For each of 16,569 nucleotide positions of the mt genome it holds two sets of four 25-mer probes each that match the heavy and the light strand of a reference mt genome and vary only at their central position to interrogate all four possible alleles. In addition, the MitoChip v2.0 carries alternative local context probes to account for known mtDNA variants. These probes have been neglected in most studies due to the lack of software for their automated analysis.

Results

We provide ReseqChip, a free software that automates the process of resequencing mtDNA using multiple local context probes on the MitoChip v2.0. ReseqChip significantly improves base call rate and sequence accuracy. ReseqChip is available at http://code.open-bio.org/svnweb/index.cgi/bioperl/browse/bioperl-live/trunk/Bio/Microarray/Tools/.

Conclusions

ReseqChip allows for the automated consolidation of base calls from alternative local mt genome context probes. It thereby improves the accuracy of resequencing, while reducing the number of non-called bases.

Fluorescent duplex allele-specific PCR and amplicon melting for rapid homogeneous mtDNA haplogroup H screening and sensitive mixture detection

BACKGROUND

For large scale studies aiming at a better understanding of mitochondrial DNA (mtDNA), sequence variation in particular mt haplogroups (hgs) and population structure, reliable low-cost high-throughput genotyping assays are needed. Furthermore, methods facilitating sensitive mixture detection and relative quantification of allele proportions are indispensable for the study of heteroplasmy, mitochondrial sequence evolution, and mitochondrial disorders. Here the properties of a homogeneous competitive duplex allele specific PCR (ARMS) assay were scrutinized in the light of these requirements.

METHODOLOGY/PRINCIPAL FINDINGS

A duplex ARMS assay amplifying either the ancestral mtDNA 2706G allele (non-hg H samples) or the derived 7028C allele (hg H samples) in the presence of SYBR Green fluorescent reporter dye was developed and characterized. Product detection, allele calling, and hg inference were based on the amplicon-characteristic melting-point temperatures obtained with on-line post-PCR fluorescent dissociation curve analysis (DCA). The analytical window of the assay covered at least 5 orders of magnitude of template DNA input with a detection limit in the low picogram range of genomic DNA. A set of forensically relevant test specimens was analyzed successfully. The presence of mtDNA mixtures was detected over a broad range of input DNA amounts and mixture ratios, and the estimation of allele proportions in samples with known total mtDNA content was feasible with limitations. A qualified DNA analyst successfully analyzed approximately 2,200 DNA extracts within three regular working days, without using robotic lab-equipment. By performing the amplification on-line, the assay also facilitated absolute mtDNA quantification.

CONCLUSIONS

Although this assay was developed just for a particular purpose, the approach is general in that it is potentially suitable in a broad variety of assay-layouts for many other applications, including the analysis of mixtures. Homogeneous ARMS-DCA is a valuable tool for large-volume studies targeting small numbers of single nucleotide polymorphisms (SNPs).

A new deep branch of eurasian mtDNA macrohaplogroup M reveals additional complexity regarding the settlement of Madagascar

Background

Current models propose that mitochondrial DNA macrohaplogroups M and N evolved from haplogroup L3 soon after modern humans left Africa. Increasingly, however, analysis of isolated populations is filling in the details of, and in some cases challenging, aspects of this general model.

Results

Here, we present the first comprehensive study of three such isolated populations from Madagascar: the Mikea hunter-gatherers, the neighbouring Vezo fishermen, and the Merina central highlanders (n = 266). Complete mitochondrial DNA genome sequences reveal several unresolved lineages, and a new, deep branch of the out-of-Africa founder clade M has been identified. This new haplogroup, M23, has a limited global distribution, and is restricted to Madagascar and a limited range of African and Southwest Asian groups.

Conclusions

The geographic distribution, phylogenetic placement and molecular age of M23 suggest that the colonization of Madagascar was more complex than previously thought.

Gene genealogies strongly distorted by weakly interfering mutations in constant environments

Neutral nucleotide diversity does not scale with population size as expected, and this "paradox of variation" is especially severe for animal mitochondria. Adaptive selective sweeps are often proposed as a major cause, but a plausible alternative is selection against large numbers of weakly deleterious mutations subject to Hill-Robertson interference. The mitochondrial genealogies of several species of whale lice (Amphipoda: Cyamus) are consistently too short relative to neutral-theory expectations, and they are also distorted in shape (branch-length proportions) and topology (relative sister-clade sizes). This pattern is not easily explained by adaptive sweeps or demographic history, but it can be reproduced in models of interference among forward and back mutations at large numbers of sites on a nonrecombining chromosome. A coalescent simulation algorithm was used to study this model over a wide range of parameter values. The genealogical distortions are all maximized when the selection coefficients are of critical intermediate sizes, such that Muller's ratchet begins to turn. In this regime, linked neutral nucleotide diversity becomes nearly insensitive to N. Mutations of this size dominate the dynamics even if there are also large numbers of more strongly and more weakly selected sites in the genome. A genealogical perspective on Hill-Robertson interference leads directly to a generalized background-selection model in which the effective population size is progressively reduced going back in time from the present.

Mitochondrial genome evidence reveals successful Late Paleolithic settlement on the Tibetan Plateau

Due to its numerous environmental extremes, the Tibetan Plateau--the world's highest plateau--is one of the most challenging areas of modern human settlement. Archaeological evidence dates the earliest settlement on the plateau to the Late Paleolithic, while previous genetic studies have traced the colonization event(s) to no earlier than the Neolithic. To explore whether the genetic continuity on the plateau has an exclusively Neolithic time depth, we studied mitochondrial DNA (mtDNA) genome variation within 6 regional Tibetan populations sampled from Tibet and neighboring areas. Our results confirm that the vast majority of Tibetan matrilineal components can trace their ancestry to Epipaleolithic and Neolithic immigrants from northern China during the mid-Holocene. Significantly, we also identified an infrequent novel haplogroup, M16, that branched off directly from the Eurasian M founder type. Its nearly exclusive distribution in Tibetan populations and ancient age (>21 kya) suggest that M16 may represent the genetic relics of the Late Paleolithic inhabitants on the plateau. This partial genetic continuity between the Paleolithic inhabitants and the contemporary Tibetan populations bridges the results and inferences from archaeology, history, and genetics.

Mitochondrial control region sequences from an African American population sample

Entire mitochondrial control region data were generated for 248 African American individuals, which had been previously typed for 15 autosomal STRs [J.M. Butler, R. Schoske, P.M. Vallone, J.W. Redman, M.C. Kline, Allele frequencies for 15 autosomal STR loci on U.S. Caucasian, African American, and Hispanic populations, J. Forensic Sci. 48 (2003) 908-911].

High mitochondrial DNA stability in B-cell chronic lymphocytic leukemia

BACKGROUND

Chronic Lymphocytic Leukemia (CLL) leads to progressive accumulation of lymphocytes in the blood, bone marrow, and lymphatic tissues. Previous findings have suggested that the mtDNA could play an important role in CLL.

METHODOLOGY/PRINCIPAL FINDINGS

The mitochondrial DNA (mtDNA) control-region was analyzed in lymphocyte cell DNA extracts and compared with their granulocyte counterpart extract of 146 patients suffering from B-Cell CLL; B-CLL (all recruited from the Basque country). Major efforts were undertaken to rule out methodological artefacts that would render a high false positive rate for mtDNA instabilities and thus lead to erroneous interpretation of sequence instabilities. Only twenty instabilities were finally confirmed, most of them affecting the homopolymeric stretch located in the second hypervariable segment (HVS-II) around position 310, which is well known to constitute an extreme mutational hotspot of length polymorphism, as these mutations are frequently observed in the general human population. A critical revision of the findings in previous studies indicates a lack of proper methodological standards, which eventually led to an overinterpretation of the role of the mtDNA in CLL tumorigenesis.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that mtDNA instability is not the primary causal factor in B-CLL. A secondary role of mtDNA mutations cannot be fully ruled out under the hypothesis that the progressive accumulation of mtDNA instabilities could finally contribute to the tumoral process. Recommendations are given that would help to minimize erroneous interpretation of sequencing results in mtDNA studies in tumorigenesis.

Human mitochondrial haplogroup H: the highest VO2max consumer--is it a paradox?

Mitochondrial background has been demonstrated to influence maximal oxygen uptake (VO(2max), in mLkg(-1)min(-1)), but this genetic influence can be compensated for by regular exercise. A positive correlation among electron transport chain (ETC) coupling, ATP and reactive oxygen species (ROS) production has been established, and mitochondrial variants have been reported to show differences in their ETC performance. In this study, we examined in detail the VO(2max) differences found among mitochondrial haplogroups. We recruited 81 healthy male Spanish Caucasian individuals and determined their mitochondrial haplogroup. Their VO(2max) was determined using incremental cycling exercise (ICE). VO(2max) was lower in J than in non-J haplogroup individuals (P=0.04). The H haplogroup was responsible for this difference (VO(2max); J vs. H; P=0.008) and this group also had significantly higher mitochondrial oxidative damage (mtOD) than the J haplogroup (P=0.04). In agreement with these results, VO(2max) and mtOD were positively correlated (P=0.01). Given that ROS production is the major contributor to mtOD and consumes four times more oxygen per electron than the ETC, our results strongly suggest that ROS production is responsible for the higher VO(2max) found in the H variant. These findings not only contribute to a better understanding of the mechanisms underneath VO(2max), but also help to explain some reported associations between mitochondrial haplogroups and mtOD with longevity, sperm motility, premature aging and susceptibility to different pathologies.

High mitogenomic evolutionary rates and time dependency

Using entire modern and ancient mitochondrial genomes of Adélie penguins (Pygoscelis adeliae) that are up to 44000 years old, we show that the rates of evolution of the mitochondrial genome are two to six times greater than those estimated from phylogenetic comparisons. Although the rate of evolution at constrained sites, including nonsynonymous positions and RNAs, varies more than twofold with time (between shallow and deep nodes), the rate of evolution at synonymous sites remains the same. The time-independent neutral evolutionary rates reported here would be useful for the study of recent evolutionary events.

Updating phylogeny of mitochondrial DNA macrohaplogroup m in India: dispersal of modern human in South Asian corridor

To construct maternal phylogeny and prehistoric dispersals of modern human being in the Indian sub continent, a diverse subset of 641 complete mitochondrial DNA (mtDNA) genomes belonging to macrohaplogroup M was chosen from a total collection of 2,783 control-region sequences, sampled from 26 selected tribal populations of India. On the basis of complete mtDNA sequencing, we identified 12 new haplogroups - M53 to M64; redefined/ascertained and characterized haplogroups M2, M3, M4, M5, M6, M8′C′Z, M9, M10, M11, M12-G, D, M18, M30, M33, M35, M37, M38, M39, M40, M41, M43, M45 and M49, which were previously described by control and/or coding-region polymorphisms. Our results indicate that the mtDNA lineages reported in the present study (except East Asian lineages M8′C′Z, M9, M10, M11, M12-G, D ) are restricted to Indian region.The deep rooted lineages of macrohaplogroup ‘M’ suggest in-situ origin of these haplogroups in India. Most of these deep rooting lineages are represented by multiple ethnic/linguist groups of India. Hierarchical analysis of molecular variation (AMOVA) shows substantial subdivisions among the tribes of India (Fst = 0.16164). The current Indian mtDNA gene pool was shaped by the initial settlers and was galvanized by minor events of gene flow from the east and west to the restricted zones. Northeast Indian mtDNA pool harbors region specific lineages, other Indian lineages and East Asian lineages. We also suggest the establishment of an East Asian gene in North East India through admixture rather than replacement.

Traces of sub-Saharan and Middle Eastern lineages in Indian Muslim populations

Islam is the second most practiced religion in India, next to Hinduism. It is still unclear whether the spread of Islam in India has been only a cultural transformation or is associated with detectable levels of gene flow. To estimate the contribution of West Asian and Arabian admixture to Indian Muslims, we assessed genetic variation in mtDNA, Y-chromosomal and LCT/MCM6 markers in 472, 431 and 476 samples, respectively, representing six Muslim communities from different geographical regions of India. We found that most of the Indian Muslim populations received their major genetic input from geographically close non-Muslim populations. However, low levels of likely sub-Saharan African, Arabian and West Asian admixture were also observed among Indian Muslims in the form of L0a2a2 mtDNA and E1b1b1a and J(*)(xJ2) Y-chromosomal lineages. The distinction between Iranian and Arabian sources was difficult to make with mtDNA and the Y chromosome, as the estimates were highly correlated because of similar gene pool compositions in the sources. In contrast, the LCT/MCM6 locus, which shows a clear distinction between the two sources, enabled us to rule out significant gene flow from Arabia. Overall, our results support a model according to which the spread of Islam in India was predominantly cultural conversion associated with minor but still detectable levels of gene flow from outside, primarily from Iran and Central Asia, rather than directly from the Arabian Peninsula.