Defeating numts: semi-pure mitochondrial DNA from eggs and simple purification methods for field-collected wildlife tissues.

Mitochondrial DNA (mtDNA) continues to play a pivotal role in phylogeographic, phylogenetic, and population genetic studies. PCR amplification with mitochondrial primers often yields ambiguous sequences, in part because of the co-amplification of nuclear copies of mitochondrial genes (numts) and true mitochondrial heteroplasmy arising from mutations, hybridization with paternal leakage, gene duplications, and recombination. Failing to detect numts or to distinguish the origin of such homologous sequences results in the incorrect interpretation of data. However, few studies obtain purified mtDNA to confirm the mitochondrial origin of the first reference sequences for a species. Here, we demonstrate the importance and ease of obtaining semi-pure mtDNA from wildlife tissues, preserved under various typical field conditions, and investigate the success of 3 commercial extraction kits, cesium-chloride gradient mtDNA purification, long-template PCR amplification, cloning, and more species-specific degenerate primers. Using more detailed avian examples, we illustrate that unfertilized or undeveloped eggs provide the purest sources of mtDNA; that kits provide an alternative to cesium-chloride gradient methods; and that long-template PCR, cloning, and degenerate primers cannot be used to produce reliable mitochondrial reference sequences, but can be powerful tools when used in conjunction with purified mtDNA stocks to distinguish numts from true heteroplasmy.     

Near panmixia at the distribution‐wide scale but evidence of genetic differentiation in a geographically isolated population of the Terek Sandpiper Xenus cinereus

Populations from different parts of a species range may vary in their genetic structure, variation and dynamics. Geographically isolated populations or those located at the periphery of the range may differ from those located in the core of the range. Such peripheral populations may harbour genetic variation important for the adaptive potential of the species. We studied the distribution‐wide population genetic structure of the Terek Sandpiper Xenus cinereus using 13 microsatellite loci and the mitochondrial DNA (mtDNA) control region. In addition, we estimated whether genetic variation changes from the core towards the edge of the breeding range. We used the results to evaluate the management needs of the sampled populations. Distribution‐wide genetic structure was negligible; the only population that showed significant genetic differentiation was the geographically isolated Dnieper River basin population in Eastern Europe. The genetic variation of microsatellites decreased towards the edge of the distribution, supporting the abundant‐centre hypotheses in which the core area of the distribution preserves the most genetic variation; however, no such trend could be seen with mtDNA. Overall genetic variation was low and there were signs of past population contractions followed by expansion; this pattern is found in most northern waders. The current effective population size (N_e) is large, and therefore global conservation measures are not necessary. However, the marginal Dnieper River population needs to be considered its own management unit. In addition, the Finnish population warrants conservation actions due to its extremely small size and degree of isolation from the main range, which makes it vulnerable to genetic depletion.     

Recombination sequences in plant mitochondrial genomes: diversity and homologies to known mitochondrial genes.

Several plant mitochondrial genomes contain repeated sequences that are postulated to be sites of homologous intragenomic recombination (1-3). In this report, we have used filter hybridizations to investigate sequence relationships between the cloned mitochondrial DNA (mtDNA) recombination repeats from turnip, spinach and maize and total mtDNA isolated from thirteen species of angiosperms. We find that strong sequence homologies exist between the spinach and turnip recombination repeats and essentially all other mitochondrial genomes tested, whereas a major maize recombination repeat does not hybridize to any other mtDNA. The sequences homologous to the turnip repeat do not appear to function in recombination in any other genome, whereas the spinach repeat hybridizes to reiterated sequences within the mitochondrial genomes of wheat and two species of pokeweed that do appear to be sites of recombination. Thus, although intragenomic recombination is a widespread phenomenon in plant mitochondria, it appears that different sequences either serve as substrates for this function in different species, or else surround a relatively short common recombination site which does not cross-hybridize under our experimental conditions. Identified gene sequences from maize mtDNA were used in heterologous hybridizations to show that the repeated sequences implicated in recombination in turnip and spinach/pokeweed/wheat mitochondria include, or are closely linked to genes for subunit II of cytochrome c oxidase and 26S rRNA, respectively. Together with previous studies indicating that the 18S rRNA gene in wheat mtDNA is contained within a recombination repeat (3), these results imply an unexpectedly frequent association between recombination repeats and plant mitochondrial genes.     

Genetic differentiation across a latitudinal gradient in two co-occurring butterfly species: revealing population differences in a context of climate change

Genetic differentiation within a species' range is determined by natural selection, genetic drift, and gene flow. Selection and drift enhance genetic differences if populations are sufficiently isolated, while gene flow precludes differentiation and local adaptation. Over large geographical areas, these processes can create a variety of scenarios, ranging from admixture to a high degree of population differentiation. Genetic differences among populations may signal functional differences within a species' range, potentially leading to population or ecotype-specific responses to global change. We investigated differentiation within the geographical range of two butterfly species along a broad latitudinal gradient. This gradient is the primary axis of climatic variation, and many ecologists expect populations at the poleward edge of this gradient to expand under climate change. Our study species inhabit a shared ecosystem and differ in body size and resource specialization; both also find their poleward range limit on an island. We find evidence for divergence of peripheral populations from the core in both taxa, suggesting the potential for genetic distinctiveness at the leading edge of climate change. We also find differences between the species in the extent of peripheral differentiation with the smaller and more specialized species showing greater population divergence (microsatellites and mtDNA) and reduced gene flow (mtDNA). Finally, gene flow estimates in both species differed strongly between two marker types. These findings suggest caution in assuming that populations are invariant across latitude and thus will respond as a single ecotype to climatic change.     

The amount and integrity of mtDNA in maize decline with development

In maize and other grasses there is a developmental gradient from the meristematic cells at the base of the stalk to the differentiated cells at the leaf tip. This gradient presents an opportunity to investigate changes in mitochondrial DNA (mtDNA) that accompany growth under light and dark conditions, as done previously for plastid DNA. Maize mtDNA was analyzed by DAPI-DNA staining of individual mitochondria, gel electrophoresis/blot hybridization, and real-time qPCR. Both the amount and integrity of the mtDNA were found to decline with development. There was a 20-fold decline in mtDNA copy number per cell from the embryo to the light-grown leaf blade. The amount of DNA per mitochondrial particle was greater in dark-grown leaf blade (24 copies, on average) than in the light (2 copies), with some mitochondria lacking any detectable DNA. Three factors that influence the demise of mtDNA during development are considered: (1) the decision to either repair or degrade mtDNA molecules that are damaged by the reactive oxygen species produced as byproducts of respiration; (2) the generation of ATP by photophosphorylation in chloroplasts, reducing the need for respiratory-competent mitochondria; and (3) the shift in mitochondrial function from energy-generating respiration to photorespiration during the transition from non-green to green tissue.     

Distribution of mitochondrial DNA variation in lake sturgeon (Adpenser fulvescens) from the Moose River basin, Ontario, Canada

We analysed mitochondrial DNA (mtDNA) variation of lake sturgeon ( Adpenser fulvescens ) from the Moose River basin. Our objective was to address various proximate and ultimate factors which may influence the distribution of lake sturgeon mtDNA haplotype lineages in this watershed. The lake sturgeon sampled were characterized by only two mtDNA hapiotypes based on a restriction fragment length polymorphism analysis with 40 restriction endonucleases and direct sequencing of 275 nucleotides in the mtDNA control region. We detected no heterogeneity in the mtDNA haplotype frequencies of lake sturgeon captured from different sites within rivers including those separated by major hydroelectric installations. However, lake sturgeon from one tributary had significantly different haplotype frequencies than those from other tributaries suggesting that they composed a discrete genetic stock. These results suggest that gene flow among most sites is significant and is an important factor affecting the distribution of mtDNA variation in this species. The genetic structuring and diversity are discussed in relation to lake sturgeon management and conservation.     

Three divergent mitochondrial genomes from California populations of the copepod Tigriopus californicus

Previous work on the harpacticoid copepod Tigriopus californicus has focused on the extensive population differentiation in three mtDNA protein coding genes (COXI, COXII, Cytb). In order to get a more complete understanding of mtDNA evolution in this species, we sequenced three complete mitochondrial genomes (one from each of three California populations) and compared them to two published mtDNA genomes from an Asian congener, Tigriopus japonicus. Several features of the mtDNA genome appear to be conserved within the genus: 1) the unique order of the protein coding genes, rRNA genes and most of the tRNA genes, 2) the genome is compact, varying between 14.3 and 14.6 kb, and 3) all genes are encoded on the same strand of the mtDNA. Within T. californicus, extremely high levels of nucleotide divergence (>20%) are observed across much of the mitochondrial genome. Inferred amino acid sequences of the proteins encoded in the mtDNAs also show high levels of divergence; at the extreme, the three ND3 variants in T. californicus showed >25% amino acid substitutions, compared with <3% amino acid divergence at the previously studied COXI locus. Unusual secondary structures make functional assignments of some tRNAs difficult. The only apparent tRNA(trp) in these genomes completely overlaps the 5' end of the 16S rRNA in all three T. californicus mtDNAs. Although not previously noted, this feature is also conserved in T. japonicus mtDNAs; whether this sequence is processed into a functional tRNA has not been determined. The putative control region contains a duplicated segment of different length (from 88 to 155 bp) in each of the T. californicus sequences. In each case, the duplicated segments are not tandem repeats; despite their different lengths, the distance between the start of the first and the start of the second repeat is conserved (520 bp). The functional significance, if any, of this repeat structure remains unknown.     

TFAM knockdown-triggered mtDNA-nucleoid aggregation and a decrease in mtDNA copy number induce the reorganization of nucleoid populations and mitochondria-associated ER-membrane contacts

The correct organization of mitochondrial DNA (mtDNA) in nucleoids and the contacts of mitochondria with the ER play an important role in maintaining the mitochondrial genome distribution within the cell. Mitochondria-associated ER membranes (MAMs) consist of interacting proteins and lipids located in the outer mitochondrial membrane and ER membrane, forming a platform for the mitochondrial inner membrane-associated genome replication factory as well as connecting the nucleoids with the mitochondrial division machinery. We show here that knockdown of a core component of mitochondrial nucleoids, TFAM, causes changes in the mitochondrial nucleoid populations, which subsequently impact ER-mitochondria membrane contacts. Knockdown of TFAM causes a significant decrease in the copy number of mtDNA as well as aggregation of mtDNA nucleoids. At the same time, it causes significant upregulation of the replicative TWNK helicase in the membrane-associated nucleoid fraction. This is accompanied by a transient elevation of MAM proteins, indicating a rearrangement of the linkage between ER and mitochondria triggered by changes in mitochondrial nucleoids. Reciprocal knockdown of the mitochondrial replicative helicase TWNK causes a decrease in mtDNA copy number and modifies mtDNA membrane association, however, it does not cause nucleoid aggregation and considerable alterations of MAM proteins in the membrane-associated fraction. Our explanation is that the aggregation of mitochondrial nucleoids resulting from TFAM knockdown triggers a compensatory mechanism involving the reorganization of both mitochondrial nucleoids and MAM. These results could provide an important insight into pathological conditions associated with impaired nucleoid organization or defects of mtDNA distribution.     

Losing your edge: climate change and the conservation value of range‐edge populations

Populations occurring at species' range edges can be locally adapted to unique environmental conditions. From a species' perspective, range‐edge environments generally have higher severity and frequency of extreme climatic events relative to the range core. Under future climates, extreme climatic events are predicted to become increasingly important in defining species' distributions. Therefore, range‐edge genotypes that are better adapted to extreme climates relative to core populations may be essential to species' persistence during periods of rapid climate change. We use relatively simple conceptual models to highlight the importance of locally adapted range‐edge populations (leading and trailing edges) for determining the ability of species to persist under future climates. Using trees as an example, we show how locally adapted populations at species' range edges may expand under future climate change and become more common relative to range‐core populations. We also highlight how large‐scale habitat destruction occurring in some geographic areas where many species range edge converge, such as biome boundaries and ecotones (e.g., the arc of deforestation along the rainforest‐cerrado ecotone in the southern Amazonia), can have major implications for global biodiversity. As climate changes, range‐edge populations will play key roles in helping species to maintain or expand their geographic distributions. The loss of these locally adapted range‐edge populations through anthropogenic disturbance is therefore hypothesized to reduce the ability of species to persist in the face of rapid future climate change.     

Mitochondrial Mutations in Subjects with Psychiatric Disorders

A considerable body of evidence supports the role of mitochondrial dysfunction in psychiatric disorders and mitochondrial DNA (mtDNA) mutations are known to alter brain energy metabolism, neurotransmission, and cause neurodegenerative disorders. Genetic studies focusing on common nuclear genome variants associated with these disorders have produced genome wide significant results but those studies have not directly studied mtDNA variants. The purpose of this study is to investigate, using next generation sequencing, the involvement of mtDNA variation in bipolar disorder, schizophrenia, major depressive disorder, and methamphetamine use. MtDNA extracted from multiple brain regions and blood were sequenced (121 mtDNA samples with an average of 8,800x coverage) and compared to an electronic database containing 26,850 mtDNA genomes. We confirmed novel and rare variants, and confirmed next generation sequencing error hotspots by traditional sequencing and genotyping methods. We observed a significant increase of non-synonymous mutations found in individuals with schizophrenia. Novel and rare non-synonymous mutations were found in psychiatric cases in mtDNA genes: ND6, ATP6, CYTB, and ND2. We also observed mtDNA heteroplasmy in brain at a locus previously associated with schizophrenia (T16519C). Large differences in heteroplasmy levels across brain regions within subjects suggest that somatic mutations accumulate differentially in brain regions. Finally, multiplasmy, a heteroplasmic measure of repeat length, was observed in brain from selective cases at a higher frequency than controls. These results offer support for increased rates of mtDNA substitutions in schizophrenia shown in our prior results. The variable levels of heteroplasmic/multiplasmic somatic mutations that occur in brain may be indicators of genetic instability in mtDNA.     

Role of phenotypic plasticity and population differentiation in adaptation to novel environmental conditions

Species can adapt to new environmental conditions either through individual phenotypic plasticity, intraspecific genetic differentiation in adaptive traits, or both. Wild emmer wheat, Triticum dicoccoides, an annual grass with major distribution in Eastern Mediterranean region, is predicted to experience in the near future, as a result of global climate change, conditions more arid than in any part of the current species distribution. To understand the role of the above two means of adaptation, and the effect of population range position, we analyzed reaction norms, extent of plasticity, and phenotypic selection across two experimental environments of high and low water availability in two core and two peripheral populations of this species. We studied 12 quantitative traits, but focused primarily on the onset of reproduction and maternal investment, which are traits that are closely related to fitness and presumably involved in local adaptation in the studied species. We hypothesized that the population showing superior performance under novel environmental conditions will either be genetically differentiated in quantitative traits or exhibit higher phenotypic plasticity than the less successful populations. We found the core population K to be the most plastic in all three trait categories (phenology, reproductive traits, and fitness) and most successful among populations studied, in both experimental environments; at the same time, the core K population was clearly genetically differentiated from the two edge populations. Our results suggest that (1) two means of successful adaptation to new environmental conditions, phenotypic plasticity and adaptive genetic differentiation, are not mutually exclusive ways of achieving high adaptive ability; and (2) colonists from some core populations can be more successful in establishing beyond the current species range than colonists from the range extreme periphery with conditions seemingly closest to those in the new environment.     

猪线粒体控制区微卫星(mtMs)变异及其分布特征

猪线粒体DNA长度因D-loop中mtMs序列重复数变异而不同。为阐明猪mtMs变异及其在不同群体中的分布特征,本研究对4个藏猪群体、八眉猪、烟台黑猪及3个引入猪群体共164个个体mtDNA D-loop全序进行测定,并与GenBank中发表的相关序列进行比对,分析了其中的重复单元核苷酸变异、重复次数及其在各群体间的分布规律。结果表明,因重复单元“5’-TGCGTACACG-3’”第2~4和10位上碱基变异,猪mtMs形成了以其为核心序列的多种重复单元(R^A^R^G)组成的复合结构,重复数介于3~47之间。单一重复(R^A)组成的mtMs结构在各群体中表现出分布优势,而大多数复合结构(R^AR^B)分布在国外选育群体中。藏猪的mtMs复合结构多达8个,R^AR^CR^E为其特有,另有5个与八眉猪共享。中国野猪和韩国地方猪种也有其特有mtMs结构。丰富的mtMs变异和特有结构与长期适应进化有关,本研究为地方猪种遗传资源及适应性研究提供了标记工具和理论基础。     

Comprehensive, rapid and sensitive detection of sequence variants of human mitochondrial tRNA genes.

In the present study, a comprehensive, rapid and sensitive method for screening sequence variation of the human mitochondrial tRNA genes has been developed. For this purpose, the denaturing gradient gel electrophoresis (DGGE) technique has been appropriately modified for simultaneous mutation analysis of a large number of samples and adapted so as to circumvent the problems caused by the anomalous electrophoretic behavior of DNA fragments encoding tRNA genes. Eighteen segments of mitochondrial DNA (mtDNA), each containing a single uniform melting domain, were selected to cover all tRNA-encoding regions using the computer program MELT94. All 18 segments were simultaneously analyzed by electrophoresis through a single broad range denaturing gradient gel under rigorously defined conditions, which prevent band broadening and other migration abnormalities from interfering with detection of sequence variants. All base substitutions tested, which include six natural mutations and 14 artificially introduced ones, have been detected successfully in the present study. Several types of evidence strongly suggest that the anomalous behavior in DGGE of tRNA gene-containing mtDNA fragments reflects their tendency to form temporary or stable alternative secondary structures under semi-denaturing conditions. The high sensitivity of the method, which can detect as low as 10% of mutant mtDNA visually, makes it valuable for the analysis of heteroplasmic mutations.     

Denaturing gradient gel method for mapping single base changes in human mitochondrial DNA.

A denaturing gradient gel electrophoresis (DGGE) method is described that detects even single base pair changes in mitochondrial DNA (mtDNA). In this method, restriction fragments of mtDNA are electrophoresed in a urea/formamide gradient gel at 60 degrees C. Migration distance of each mtDNA fragment in the gel depends on melting behavior which reflects base composition. Fragments are located by Southern blotting with specific mtDNA probes. With just four carefully chosen restriction enzymes and as little as 50-100 ng of mtDNA, the method covers almost the entire human mitochondrial genome. To demonstrate the method, human mtDNA was analyzed. In six normal individuals, DGGE revealed melting behavior polymorphisms (MBPs) in mtDNA fragments that were not detected by restriction fragment length polymorphism (RFLP) analysis in agarose gels. Another individual, shown to have a melting behavior polymorphism in the cytochrome b coding region, was studied in detail. By mapping, the mutation was deduced to lie between nt 14905 and 15370. The affected fragment was amplified by PCR and sequenced. Specific base changes were identified in the region predicted by the gel result. This method will be especially useful as a diagnostic tool in mitochondrial disease for rapid localization of mtDNA mutations to specific regions of the genome, but DGGE also could complement RFLP analysis as a more sensitive method to follow maternal lineage in human and animal populations in a variety of research fields.     

The mitochondrial genome of Malus domestica and the import-driven hypothesis of mitochondrial genome expansion in seed plants

Mitochondrial genomes of spermatophytes are the largest of all organellar genomes. Their large size has been attributed to various factors; however, the relative contribution of these factors to mitochondrial DNA (mtDNA) expansion remains undetermined. We estimated their relative contribution in Malus domestica (apple). The mitochondrial genome of apple has a size of 396 947 bp and a one to nine ratio of coding to non-coding DNA, close to the corresponding average values for angiosperms. We determined that 71.5% of the apple mtDNA sequence was highly similar to sequences of its nuclear DNA. Using nuclear gene exons, nuclear transposable elements and chloroplast DNA as markers of promiscuous DNA content in mtDNA, we estimated that approximately 20% of the apple mtDNA consisted of DNA sequences imported from other cell compartments, mostly from the nucleus. Similar marker-based estimates of promiscuous DNA content in the mitochondrial genomes of other species ranged between 21.2 and 25.3% of the total mtDNA length for grape, between 23.1 and 38.6% for rice, and between 47.1 and 78.4% for maize. All these estimates are conservative, because they underestimate the import of non-functional DNA. We propose that the import of promiscuous DNA is a core mechanism for mtDNA size expansion in seed plants. In apple, maize and grape this mechanism contributed far more to genome expansion than did homologous recombination. In rice the estimated contribution of both mechanisms was found to be similar.     

A replication-linked mutational gradient drives somatic mutation accumulation and influences germline polymorphisms and genome composition in mitochondrial DNA

Mutations in mitochondrial DNA (mtDNA) cause maternally inherited diseases, while somatic mutations are linked to common diseases of aging. Although mtDNA mutations impact health, the processes that give rise to them are under considerable debate. To investigate the mechanism by which de novo mutations arise, we analyzed the distribution of naturally occurring somatic mutations across the mouse and human mtDNA obtained by Duplex Sequencing. We observe distinct mutational gradients in G→A and T→C transitions delimited by the light-strand origin and the mitochondrial Control Region (mCR). The gradient increases unequally across the mtDNA with age and is lost in the absence of DNA polymerase γ proofreading activity. In addition, high-resolution analysis of the mCR shows that important regulatory elements exhibit considerable variability in mutation frequency, consistent with them being mutational ‘hot-spots’ or ‘cold-spots’. Collectively, these patterns support genome replication via a deamination prone asymmetric strand-displacement mechanism as the fundamental driver of mutagenesis in mammalian DNA. Moreover, the distribution of mtDNA single nucleotide polymorphisms in humans and the distribution of bases in the mtDNA across vertebrate species mirror this gradient, indicating that replication-linked mutations are likely the primary source of inherited polymorphisms that, over evolutionary timescales, influences genome composition during speciation.     

The telomeres of the linear mitochondrial DNA of tetrahymena thermophila consist of 53 bp tandem repeats

We have cloned and sequenced the telomeric DNA of the linear mitochondrial DNA (mtDNA) of T. thermophila BVII. The mtDNA telomeres consist of a 53 bp sequence tandemly repeated from 4 to 30 times, with most molecules having 15 ± 4 repetitions. The previously recognized terminal heterogeneity of the mtDNA is completely accounted for by the variability in the number of repeats. The 53 bp repeat does not resemble known telomeric DNA in sequence, repeat size, or number of repetitions. The termini occur at heterogeneous positions within the 53 bp repeat. The junction of the telomeric repeat with the internal DNA is at a different position within the telomeric repeat on each end of the mtDNA. We propose a model for the maintenance of the mtDNA ends involving unequal homologous recombination.     

Mitochondrial DNA sequence variation is associated with free-living activity energy expenditure in the elderly

The decline in activity energy expenditure underlies a range of age-associated pathological conditions, neuromuscular and neurological impairments, disability, and mortality. The majority (90%) of the energy needs of the human body are met by mitochondrial oxidative phosphorylation (OXPHOS). OXPHOS is dependent on the coordinated expression and interaction of genes encoded in the nuclear and mitochondrial genomes. We examined the role of mitochondrial genomic variation in free-living activity energy expenditure (AEE) and physical activity levels (PAL) by sequencing the entire (~16.5 kilobases) mtDNA from 138 Health, Aging, and Body Composition Study participants. Among the common mtDNA variants, the hypervariable region 2 m.185G>A variant was significantly associated with AEE (p=0.001) and PAL (p=0.0005) after adjustment for multiple comparisons. Several unique nonsynonymous variants were identified in the extremes of AEE with some occurring at highly conserved sites predicted to affect protein structure and function. Of interest is the p.T194M, CytB substitution in the lower extreme of AEE occurring at a residue in the Qi site of complex III. Among participants with low activity levels, the burden of singleton variants was 30% higher across the entire mtDNA and OXPHOS complex I when compared to those having moderate to high activity levels. A significant pooled variant association across the hypervariable 2 region was observed for AEE and PAL. These results suggest that mtDNA variation is associated with free-living AEE in older persons and may generate new hypotheses by which specific mtDNA complexes, genes, and variants may contribute to the maintenance of activity levels in late life.     

Genetic structure of pastoral and farmer populations in the African Sahel

Traditional pastoralists survive in few places in the world. They can still be encountered in the African Sahel, where annual alternations of dry and wet seasons force them to continual mobility. Little is known about the genetic structure of these populations. We present here the population distribution of 312 hypervariable segment I mitochondrial DNA (mtDNA) and 364 Y-short tandem repeat haplotypes in both farmer and pastoralist groups from the Lake Chad Basin and the West African Sahel. We show that the majority of pastoral populations (represented in the African Sahel by the Fulani nomads) fail to show significant departure from neutrality for mtDNA as evidenced by Fu's Fs statistics and exhibit lower levels of intrapopulation diversity measures for mtDNA when contrasted with farmers. These differences were not observed for the Y chromosome. Furthermore, analyses of molecular variance and population distributions of the mtDNA haplotypes show more heterogeneity in the sedentary groups than in the pastoralists. On the other hand, pastoralists retain a signature of a wide phylogenetic distance contributing to their male gene pool, whereas in at least some of the farmer populations, a founder effect and/or drift might have led to the presence of a single major lineage. Interestingly, these observations are in contrast with those recorded in Central Asia, where similar comparisons of farmer and pastoral groups have recently been carried out. We can conclude that in Africa, there have been no substantial mating exchanges between the Fulani pastoralists coming to the Lake Chad Basin from the West African Sahel and their farmer neighbors. At the same time, we suggest that the emergence of pastoralism might be an earlier and/or a demographically more important event than the introduction of sedentary agriculture, at least in this part of Africa.