Y-chromosome lineages in Cabo Verde Islands witness the diverse geographic origin of its first male settlers

The Y-chromosome haplogroup composition of the population of the Cabo Verde Archipelago was profiled by using 32 single-nucleotide polymorphism markers and compared with potential source populations from Iberia, west Africa, and the Middle East. According to the traditional view, the major proportion of the founding population of Cabo Verde was of west African ancestry with the addition of a minor fraction of male colonizers from Europe. Unexpectedly, more than half of the paternal lineages (53.5%) of Cabo Verdeans clustered in haplogroups I, J, K, and R1, which are characteristic of populations of Europe and the Middle East, while being absent in the probable west African source population of Guiné-Bissau. Moreover, a high frequency of J* lineages in Cabo Verdeans relates them more closely to populations of the Middle East and probably provides the first genetic evidence of the legacy of the Jews. In addition, the considerable proportion (20.5%) of E3b(xM81) lineages indicates a possible gene flow from the Middle East or northeast Africa, which, at least partly, could be ascribed to the Sephardic Jews. In contrast to the predominance of west African mitochondrial DNA haplotypes in their maternal gene pool, the major west African Y-chromosome lineage E3a was observed only at a frequency of 15.9%. Overall, these results indicate that gene flow from multiple sources and various sex-specific patterns have been important in the formation of the genomic diversity in the Cabo Verde islands.An erratum to this article can be found at     

Mitochondrial portraits of the Madeira and Açores archipelagos witness different genetic pools of its settlers

We have studied the matrilineal genetic composition of the Madeira and Açores north Atlantic archipelagos, which were settled by the Portuguese in the 15th century. Both archipelagos, and particularly Madeira, were involved in a complex commercial network established by the Portuguese, which included the trading of slaves across the Atlantic. One hundred and fifty-five mtDNAs sampled from the Madeira and 179 from the Açores archipelagos were analysed for the hypervariable segment I (HVS-I), and for haplogroup-diagnostic coding-region RFLPs. The different settlement histories of both groups of islands are well reflected in their present day mtDNA pool. Although both archipelagos show identical diversity values, they are clearly different in their haplogroup content. Madeira displays a stronger sub-Saharan imprint, with haplogroups L1–L3 constituting about 13% of the lineages. Also, the relative frequencies of L sub-clusters in Madeira and mainland Portugal suggests that, at least in part, African presence in Madeira can be attributed to a direct gene flow from West Africa and not via Portugal. A comparison of the genetic composition of these two archipelagos with the Canary Islands, specially taking into account that their European source population was essentially from the Iberian Peninsula, testifies the stronger impact of the North African U6 cluster in the Canaries. This group is present in Madeira at a moderate frequency, but very reduced in the Açores. Nevertheless the recorded introduction of Canary native Guanches, who are characterized by the presence of particular sub-clade U6b1, has left no detectable imprints in the present day population of Madeira.     

Biogeographic,climatic and spatial drivers differentially affect α-, β- and γ-diversities on oceanic archipelagos

Island biogeographic studies traditionally treat single islands as units of analysis. This ignores the fact that most islands are spatially nested within archipelagos. Here, we took a fundamentally different approach and focused on entire archipelagos using species richness of vascular plants on 23 archipelagos worldwide and their 174 constituent islands. We assessed differential effects of biogeographic factors (area, isolation, age, elevation), current and past climate (temperature, precipitation, seasonality, climate change velocity) and intra-archipelagic spatial structure (archipelago area, number of islands, area range, connectivity, environmental volume, inter-island distance) on plant diversity. Species diversity of each archipelago (γ) was additively partitioned into α, β, nestedness and replacement β-components to investigate the relative importance of environmental and spatial drivers. Multiple regressions revealed strong effects of biogeography and climate on α and γ, whereas spatial factors, particularly number of islands, inter-island distance and area range, were key to explain β. Structural equation models additionally suggested that γ is predominantly determined by indirect abiotic effects via its components, particularly β. This highlights that β and the spatial arrangement of islands are essential to understand insular ecology and evolution. Our methodological framework can be applied more widely to other taxa and archipelago-like systems, allowing new insights into biodiversity origin and maintenance.     

Mitochondrial DNA patterns in the Macaronesia islands: Variation within and among archipelagos

Macaronesia covers four Atlantic archipelagos: the Azores, Madeira, the Canary Islands, and the Cape Verde islands. When discovered by Europeans in the 15th century, only the Canaries were inhabited. Historical reports highlight the impact of Iberians on settlement in Macaronesia. Although important differences in their settlement are documented, its influence on their genetic structures and relationships has yet to be ascertained. In this study, the hypervariable region I (HVRI) sequence and coding region polymorphisms of mitochondrial DNA (mtDNA) in 623 individuals from the Azores (120) and Canary Islands (503) were analyzed. Combined with published data, these give a total of 1,542 haplotypes from Macaronesia and 1,067 from the Iberian Peninsula. The results obtained indicate that Cape Verde is the most distinctive archipelago, with an mtDNA pool composed almost exclusively of African lineages. However, the other archipelagos present an mtDNA profile dominated by the presence of West‐Eurasian mtDNA haplogroups with African lineages present in varying proportions. Moreover, no signs of integration of typical Canarian U6 lineages in the other archipelagos were detected. The four Macaronesia archipelagos currently have differentiated genetic profiles, and the Azores present the highest intra‐archipelago differentiation and the lowest values of diversity. The analyses performed show that the present‐day genetic profile of the Macaronesian archipelagos was mainly determined by the initial process of settlement and further microdifferentiation probably as a consequence of the small population size of some islands. Moreover, contacts between archipelagos seem to have had a low impact on the mtDNA genetic pool of each archipelago. Am J Phys Anthropol, 2010. © 2009 Wiley‐Liss, Inc.     

Deep Rooting In-Situ Expansion of mtDNA Haplogroup R8 in South Asia

Background

The phylogeny of the indigenous Indian-specific mitochondrial DNA (mtDNA) haplogroups have been determined and refined in previous reports. Similar to mtDNA superhaplogroups M and N, a profusion of reports are also available for superhaplogroup R. However, there is a dearth of information on South Asian subhaplogroups in particular, including R8. Therefore, we ought to access the genealogy and pre-historic expansion of haplogroup R8 which is considered one of the autochthonous lineages of South Asia.

Methodology/Principal Findings

Upon screening the mtDNA of 5,836 individuals belonging to 104 distinct ethnic populations of the Indian subcontinent, we found 54 individuals with the HVS-I motif that defines the R8 haplogroup. Complete mtDNA sequencing of these 54 individuals revealed two deep-rooted subclades: R8a and R8b. Furthermore, these subclades split into several fine subclades. An isofrequency contour map detected the highest frequency of R8 in the state of Orissa. Spearman''s rank correlation analysis suggests significant correlation of R8 occurrence with geography.

Conclusions/Significance

The coalescent age of newly-characterized subclades of R8, R8a (15.4±7.2 Kya) and R8b (25.7±10.2 Kya) indicates that the initial maternal colonization of this haplogroup occurred during the middle and upper Paleolithic period, roughly around 40 to 45 Kya. These results signify that the southern part of Orissa currently inhabited by Munda speakers is likely the origin of these autochthonous maternal deep-rooted haplogroups. Our high-resolution study on the genesis of R8 haplogroup provides ample evidence of its deep-rooted ancestry among the Orissa (Austro-Asiatic) tribes.     

Population Differentiation of Southern Indian Male Lineages Correlates with Agricultural Expansions Predating the Caste System

Previous studies that pooled Indian populations from a wide variety of geographical locations, have obtained contradictory conclusions about the processes of the establishment of the Varna caste system and its genetic impact on the origins and demographic histories of Indian populations. To further investigate these questions we took advantage that both Y chromosome and caste designation are paternally inherited, and genotyped 1,680 Y chromosomes representing 12 tribal and 19 non-tribal (caste) endogamous populations from the predominantly Dravidian-speaking Tamil Nadu state in the southernmost part of India. Tribes and castes were both characterized by an overwhelming proportion of putatively Indian autochthonous Y-chromosomal haplogroups (H-M69, F-M89, R1a1-M17, L1-M27, R2-M124, and C5-M356; 81% combined) with a shared genetic heritage dating back to the late Pleistocene (10–30 Kya), suggesting that more recent Holocene migrations from western Eurasia contributed <20% of the male lineages. We found strong evidence for genetic structure, associated primarily with the current mode of subsistence. Coalescence analysis suggested that the social stratification was established 4–6 Kya and there was little admixture during the last 3 Kya, implying a minimal genetic impact of the Varna (caste) system from the historically-documented Brahmin migrations into the area. In contrast, the overall Y-chromosomal patterns, the time depth of population diversifications and the period of differentiation were best explained by the emergence of agricultural technology in South Asia. These results highlight the utility of detailed local genetic studies within India, without prior assumptions about the importance of Varna rank status for population grouping, to obtain new insights into the relative influences of past demographic events for the population structure of the whole of modern India.     

Population genomics and conservation management of a declining tropical rodent

Conservation management is improved by incorporating information about the spatial distribution of population genetic diversity into planning strategies. Northern Australia is the location of some of the world’s most severe ongoing declines of endemic mammal species, yet we have little genetic information from this regional mammal assemblage to inform a genetic perspective on conservation assessment and planning. We used next-generation sequencing data from remnant populations of the threatened brush-tailed rabbit-rat (Conilurus penicillatus) to compare patterns of genomic diversity and differentiation across the landscape and investigate standardised hierarchical genomic diversity metrics to better understand brush-tailed rabbit-rat population genomic structure. We found strong population structuring, with high levels of differentiation between populations (F_ST = 0.21–0.78). Two distinct genomic lineages between the Tiwi Islands and mainland are also present. Prioritisation analysis showed that one population in both lineages would need to be conserved to retain at least ~80% of alleles for the species. Analysis of standardised genomic diversity metrics showed that approximately half of the total diversity occurs among lineages (δ = 0.091 from grand total γ = 0.184). We suggest that a focus on conserving remnant island populations may not be appropriate for the preservation of species-level genomic diversity and adaptive potential, as these populations represent a small component of the total diversity and a narrow subset of the environmental conditions in which the species occurs. We also highlight the importance of considering both genomic and ecological differentiation between source and receiving populations when considering translocations for conservation purposes.Subject terms: Ecological genetics, Population genetics, Conservation biology, Biogeography     

Evidence that disease-induced population decline changes genetic structure and alters dispersal patterns in the Tasmanian devil

Infectious disease has been shown to be a major cause of population declines in wild animals. However, there remains little empirical evidence on the genetic consequences of disease-mediated population declines, or how such perturbations might affect demographic processes such as dispersal. Devil facial tumour disease (DFTD) has resulted in the rapid decline of the Tasmanian devil, Sarcophilus harrisii, and threatens to cause extinction. Using 10 microsatellite DNA markers, we compared genetic diversity and structure before and after DFTD outbreaks in three Tasmanian devil populations to assess the genetic consequences of disease-induced population decline. We also used both genetic and demographic data to investigate dispersal patterns in Tasmanian devils along the east coast of Tasmania. We observed a significant increase in inbreeding (F_IS pre/post-disease −0.030/0.012, P<0.05; relatedness pre/post-disease 0.011/0.038, P=0.06) in devil populations after just 2–3 generations of disease arrival, but no detectable change in genetic diversity. Furthermore, although there was no subdivision apparent among pre-disease populations (θ=0.005, 95% confidence interval (CI) −0.003 to 0.017), we found significant genetic differentiation among populations post-disease (θ=0.020, 0.010–0.027), apparently driven by a combination of selection and altered dispersal patterns of females in disease-affected populations. We also show that dispersal is male-biased in devils and that dispersal distances follow a typical leptokurtic distribution. Our results show that disease can result in genetic and demographic changes in host populations over few generations and short time scales. Ongoing management of Tasmanian devils must now attempt to maintain genetic variability in this species through actions designed to reverse the detrimental effects of inbreeding and subdivision in disease-affected populations.     

Low Genetic Diversity and Strong Geographical Structure of the Critically Endangered White-Headed Langur (Trachypithecus leucocephalus) Inferred from Mitochondrial DNA Control Region Sequences

Many Asian colobine monkey species are suffering from habitat destruction and population size decline. There is a great need to understand their genetic diversity, population structure and demographic history for effective species conservation. The white-headed langur (Trachypithecus leucocephalus) is a Critically Endangered colobine species endemic to the limestone karst forests in southwestern China. We analyzed the mitochondrial DNA (mtDNA) control region sequences of 390 fecal samples from 40 social groups across the main distribution areas, which represented one-third of the total extant population. Only nine haplotypes and 10 polymorphic sites were identified, indicating remarkably low genetic diversity in the species. Using a subset of 77 samples from different individuals, we evaluated genetic variation, population structure, and population demographic history. We found very low values of haplotype diversity (h = 0.570 ± 0.056) and nucleotide diversity (π = 0.00323 ± 0.00044) in the hypervariable region I (HVRI) of the mtDNA control region. Distribution of haplotypes displayed marked geographical pattern, with one population (Chongzuo, CZ) showing a complete lack of genetic diversity (having only one haplotype), whereas the other population (Fusui, FS) having all nine haplotypes. We detected strong population genetic structure among habit patches (Φ _ST = 0.375, P < 0.001). In addition, the Mantel test showed a significant correlation between the pairwise genetic distances and geographical distances among social groups in FS (correlation coefficient = 0.267, P = 0.003), indicting isolation-by-distance pattern of genetic divergence in the mtDNA sequences. Analyses of demographic history suggested an overall stable historical population size and modest population expansion in the last 2,000 years. Our results indicate different genetic diversity and possibly distinct population history for different local populations, and suggest that CZ and FS should be considered as one evolutionarily significant unit (ESU) and two management units (MUs) pending further investigation using nuclear markers.     

Genetic structure of the population of Cabo Verde (West Africa): Evidence of substantial european admixture

The population of Cabo Verde was founded in the fifteenth century (1462), on the basis of slaves brought from the West African coast and a few Europeans, mainly from Portugal. The polymorphism of six red cell enzymes (ADA, AK1, ALAD, ESD, GLO1, and PGD) and ten plasma proteins (AHSG, BF, F13A, F13B, GC, HP, ORM, PLG, TBG, and TF) was studied in a sample of 268 individuals from Cabo Verde (West Africa). There is no statistical evidence of genetic heterogeneity between the two groups of islands which constitute the archipelago, Barlavento and Sotavento. The gene frequency distribution observed in Cabo Verde differs, in many markers, from that of West African populations, suggesting an important European influence. The proportion of Caucasian genes in the population of Cabo Verde has been calculated to be M = 0.3634 ± 0.0510, and the considerable dispersion of the locus-specific admixture estimates seems to indicate random drift has also played a role in the evolution of the allele frequencies in the archipelago. Partition of the variance of the mean estimate in evolutionary and sampling variance shows the evolutionary variance is more than ten times higher than the sampling variance. When dendrograms are constructed on the basis of different genetic distances, the population of Cabo Verde clusters with Afro-Americans, forming a different group from the populations of the African continent. This is interpreted as a consequence of the importance of Caucasian admixture both in Afro-Americans and in the population of Cabo Verde. © 1995 Wiley-Liss, Inc.     

Population Genetics of the S?o Tomé Caecilian (Gymnophiona: Dermophiidae: Schistometopum thomense) Reveals Strong Geographic Structuring

Islands provide exciting opportunities for exploring ecological and evolutionary mechanisms. The oceanic island of São Tomé in the Gulf of Guinea exhibits high diversity of fauna including the endemic caecilian amphibian, Schistometopum thomense. Variation in pigmentation, morphology and size of this taxon over its c. 45 km island range is extreme, motivating a number of taxonomic, ecological, and evolutionary hypotheses to explain the observed diversity. We conducted a population genetic study of S. thomense using partial sequences of two mitochondrial DNA genes (ND4 and 16S), together with morphological examination, to address competing hypotheses of taxonomic or clinal variation. Using Bayesian phylogenetic analysis and Spatial Analysis of Molecular Variance, we found evidence of four geographic clades, whose range and approximated age (c. 253 Kya – 27 Kya) are consistent with the spread and age of recent volcanic flows. These clades explained 90% of variation in ND4 (φ_CT = 0.892), and diverged by 4.3% minimum pairwise distance at the deepest node. Most notably, using Mismatch Distributions and Mantel Tests, we identified a zone of population admixture that dissected the island. In the northern clade, we found evidence of recent population expansion (Fu''s F_s = −13.08 and Tajima''s D = −1.80) and limited dispersal (Mantel correlation coefficient = 0.36, p = 0.01). Color assignment to clades was not absolute. Paired with multinomial regression of chromatic data, our analyses suggested that the genetic groups and a latitudinal gradient together describe variation in color of S. thomense. We propose that volcanism and limited dispersal ability are the likely proximal causes of the observed genetic structure. This is the first population genetic study of any caecilian and demonstrates that these animals have deep genetic divisions over very small areas in accordance with previous speculations of low dispersal abilities.     

The Plasmodium vivax Merozoite Surface Protein 3β Sequence Reveals Contrasting Parasite Populations in Southern and Northwestern Thailand

Background

Malaria control efforts have a significant impact on the epidemiology and parasite population dynamics. In countries aiming for malaria elimination, malaria transmission may be restricted to limited transmission hot spots, where parasite populations may be isolated from each other and experience different selection forces. Here we aim to examine the Plasmodium vivax population divergence in geographically isolated transmission zones in Thailand.

Methodology

We employed the P. vivax merozoite surface protein 3β (PvMSP3β) as a molecular marker for characterizing P. vivax populations based on the extensive diversity of this gene in Southeast Asian parasite populations. To examine two parasite populations with different transmission levels in Thailand, we obtained 45 P. vivax isolates from Tak Province, northwestern Thailand, where the annual parasite incidence (API) was more than 2%, and 28 isolates from Yala and Narathiwat Provinces, southern Thailand, where the API was less than 0.02%. We sequenced the PvMSP3β gene and examined its genetic diversity and molecular evolution between the parasite populations.

Principal Findings

Of 58 isolates containing single PvMSP3β alleles, 31 sequence types were identified. The overall haplotype diversity was 0.77±0.06 and nucleotide diversity 0.0877±0.0054. The northwestern vivax malaria population exhibited extensive haplotype diversity (HD) of PvMSP3β (HD = 1.0). In contrast, the southern parasite population displayed a single PvMSP3β allele (HD = 0), suggesting a clonal population expansion. This result revealed that the extent of allelic diversity in P. vivax populations in Thailand varies among endemic areas.

Conclusion

Malaria parasite populations in a given region may vary significantly in genetic diversity, which may be the result of control and influenced by the magnitude of malaria transmission intensity. This is an issue that should be taken into account for the implementation of P. vivax control measures such as drug policy and vaccine development.     

Region‐wide and ecotype‐specific differences in demographic histories of threespine stickleback populations,estimated from whole genome sequences

We analysed 81 whole genome sequences of threespine sticklebacks from Pacific North America, Greenland and Northern Europe, representing 16 populations. Principal component analysis of nuclear SNPs grouped populations according to geographical location, with Pacific populations being more divergent from each other relative to European and Greenlandic populations. Analysis of mitogenome sequences showed Northern European populations to represent a single phylogeographical lineage, whereas Greenlandic and particularly Pacific populations showed admixture between lineages. We estimated demographic history using a genomewide coalescence with recombination approach. The Pacific populations showed gradual population expansion starting >100 Kya, possibly reflecting persistence in cryptic refuges near the present distributional range, although we do not rule out possible influence of ancient admixture. Sharp population declines ca. 14–15 Kya were suggested to reflect founding of freshwater populations by marine ancestors. In Greenland and Northern Europe, demographic expansion started ca. 20–25 Kya coinciding with the end of the Last Glacial Maximum. In both regions, marine and freshwater populations started to show different demographic trajectories ca. 8–9 Kya, suggesting that this was the time of recolonization. In Northern Europe, this estimate was surprisingly late, but found support in subfossil evidence for presence of several freshwater fish species but not sticklebacks 12 Kya. The results demonstrate distinctly different demographic histories across geographical regions with potential consequences for adaptive processes. They also provide empirical support for previous assumptions about freshwater populations being founded independently from large, coherent marine populations, a key element in the Transporter Hypothesis invoked to explain the widespread occurrence of parallel evolution across freshwater stickleback populations.     

The Q2 Mitochondrial Haplogroup in Oceania

Many details surrounding the origins of the peoples of Oceania remain to be resolved, and as a step towards this we report seven new complete mitochondrial genomes from the Q2a haplogroup, from Papua New Guinea, Fiji and Kiribati. This brings the total to eleven Q2 genomes now available. The Q haplogroup (that includes Q2) is an old and diverse lineage in Near Oceania, and is reasonably common; within our sample set of 430, 97 are of the Q haplogroup. However, only 8 are Q2, and we report 7 here. The tree with all complete Q genomes is proven to be minimal. The dating estimate for the origin of Q2 (around 35 Kya) reinforces the understanding that humans have been in Near Oceania for tens of thousands of years; nevertheless the Polynesian maternal haplogroups remain distinctive. A major focus now, with regard to Polynesian ancestry, is to address the differences and timing of the ‘Melanesian’ contribution to the maternal and paternal lineages as people moved further and further into Remote Oceania. Input from other fields such as anthropology, history and linguistics is required for a better understanding and interpretation of the genetic data.     

Gene Pool Structure of Eastern Ukrainians as Inferred from the Y-Chromosome Haplogroups

Y chromosomes from representative sample of Eastern Ukrainians (94 individuals) were analyzed for composition and frequencies of haplogroups, defined by 11 biallelic loci located in non-recombining part of the chromosome (SRY1532, YAP, 92R7, DYF155S2, 12f2, Tat, M9, M17, M25,M89, andM56). In the Ukrainian gene, pool six haplogroups were revealed: E, F (including G and I), J, N3, P, and R1a1. These haplogroups were earlier detected in a study of Y-chromosome diversity on the territory of Europe as a whole. The major haplogroup in the Ukrainian gene pool, haplogroup R1a1 (earlier designated HG3), accounted for about 44% of all Y chromosomes in the sample examined. This haplogroup is thought to mark the migration patterns of the early Indo-Europeans and is associated with the distribution of the Kurgan archaeological culture. The second major haplogroup is haplogroup F (21.3%), which is a combination of the lineages differing by the time of appearance. Haplogroup P found with the frequency of 9.6%, represents the genetic contribution of the population originating from the ancient autochthonous population of Europe. Haplogroups J and E (11.7 and 4.2%, respectively) mark the migration patterns of the Middle-Eastern agriculturists during the Neolithic. The presence of the N3 lineage (9.6%) is likely explained by a contribution of the assimilated Finno–Ugric tribes. The data on the composition and frequencies of Y-chromosome haplogroups in the sample studied substantially supplement the existing picture of the male lineage distribution in the Eastern Slav population.     

Genetic diversity and structure of Rhododendron meddianum,a plant species with extremely small populations

Rhododendron meddianum is a critically endangered species with important ornamental value and is also a plant species with extremely small populations. In this study, we used double digest restriction-site-associated DNA sequencing (ddRAD) technology to assess the genetic diversity, genetic structure and demographic history of the three extant populations of R. meddianum. Analysis of SNPs indicated that R. meddianum populations have a high genetic diversity (π = 0.0772 ± 0.0024, H_E = 0.0742 ± 0.002). Both F_ST values (0.1582–0.2388) and AMOVA showed a moderate genetic differentiation among the R. meddianum populations. Meanwhile, STRUCTURE, PCoA and NJ trees indicated that the R. meddianum samples were clustered into three distinct genetic groups. Using the stairway plot, we found that R. meddianum underwent a population bottleneck about 70,000 years ago. Furthermore, demographic models of R. meddianum and its relative, Rhododendron cyanocarpum, revealed that these species diverged about 3.05 (2.21–5.03) million years ago. This divergence may have been caused by environmental changes that occurred after the late Pliocene, e.g., the Asian winter monsoon intensified, leading to a drier climate. Based on these findings, we recommend that R. meddianum be conserved through in situ, ex situ approaches and that its seeds be collected for germplasm.     

Statistical tests for detecting positive selection by utilizing high-frequency variants

By comparing the low-, intermediate-, and high-frequency parts of the frequency spectrum, we gain information on the evolutionary forces that influence the pattern of polymorphism in population samples. We emphasize the high-frequency variants on which positive selection and negative (background) selection exhibit different effects. We propose a new estimator of θ (the product of effective population size and neutral mutation rate), θ_L, which is sensitive to the changes in high-frequency variants. The new θ_L allows us to revise Fay and Wu's H-test by normalization. To complement the existing statistics (the H-test and Tajima's D-test), we propose a new test, E, which relies on the difference between θ_L and Watterson's θ_W. We show that this test is most powerful in detecting the recovery phase after the loss of genetic diversity, which includes the postselective sweep phase. The sensitivities of these tests to (or robustness against) background selection and demographic changes are also considered. Overall, D and H in combination can be most effective in detecting positive selection while being insensitive to other perturbations. We thus propose a joint test, referred to as the DH test. Simulations indicate that DH is indeed sensitive primarily to directional selection and no other driving forces.     

The Expanded mtDNA Phylogeny of the Franco-Cantabrian Region Upholds the Pre-Neolithic Genetic Substrate of Basques

The European genetic landscape has been shaped by several human migrations occurred since Paleolithic times. The accumulation of archaeological records and the concordance of different lines of genetic evidence during the last two decades have triggered an interesting debate concerning the role of ancient settlers from the Franco-Cantabrian region in the postglacial resettlement of Europe. Among the Franco-Cantabrian populations, Basques are regarded as one of the oldest and more intriguing human groups of Europe. Recent data on complete mitochondrial DNA genomes focused on macrohaplogroup R0 revealed that Basques harbor some autochthonous lineages, suggesting a genetic continuity since pre-Neolithic times. However, excluding haplogroup H, the most representative lineage of macrohaplogroup R0, the majority of maternal lineages of this area remains virtually unexplored, so that further refinement of the mtDNA phylogeny based on analyses at the highest level of resolution is crucial for a better understanding of the European prehistory. We thus explored the maternal ancestry of 548 autochthonous individuals from various Franco-Cantabrian populations and sequenced 76 mitogenomes of the most representative lineages. Interestingly, we identified three mtDNA haplogroups, U5b1f, J1c5c1 and V22, that proved to be representative of Franco-Cantabria, notably of the Basque population. The seclusion and diversity of these female genetic lineages support a local origin in the Franco-Cantabrian area during the Mesolithic of southwestern Europe, ∼10,000 years before present (YBP), with signals of expansions at ∼3,500 YBP. These findings provide robust evidence of a partial genetic continuity between contemporary autochthonous populations from the Franco-Cantabrian region, specifically the Basques, and Paleolithic/Mesolithic hunter-gatherer groups. Furthermore, our results raise the current proportion (≈15%) of the Franco-Cantabrian maternal gene pool with a putative pre-Neolithic origin to ≈35%, further supporting the notion of a predominant Paleolithic genetic substrate in extant European populations.     

Geographic Patterns of Genetic Variation in a Broadly Distributed Marine Vertebrate: New Insights into Loggerhead Turtle Stock Structure from Expanded Mitochondrial DNA Sequences

Previous genetic studies have demonstrated that natal homing shapes the stock structure of marine turtle nesting populations. However, widespread sharing of common haplotypes based on short segments of the mitochondrial control region often limits resolution of the demographic connectivity of populations. Recent studies employing longer control region sequences to resolve haplotype sharing have focused on regional assessments of genetic structure and phylogeography. Here we synthesize available control region sequences for loggerhead turtles from the Mediterranean Sea, Atlantic, and western Indian Ocean basins. These data represent six of the nine globally significant regional management units (RMUs) for the species and include novel sequence data from Brazil, Cape Verde, South Africa and Oman. Genetic tests of differentiation among 42 rookeries represented by short sequences (380 bp haplotypes from 3,486 samples) and 40 rookeries represented by long sequences (∼800 bp haplotypes from 3,434 samples) supported the distinction of the six RMUs analyzed as well as recognition of at least 18 demographically independent management units (MUs) with respect to female natal homing. A total of 59 haplotypes were resolved. These haplotypes belonged to two highly divergent global lineages, with haplogroup I represented primarily by CC-A1, CC-A4, and CC-A11 variants and haplogroup II represented by CC-A2 and derived variants. Geographic distribution patterns of haplogroup II haplotypes and the nested position of CC-A11.6 from Oman among the Atlantic haplotypes invoke recent colonization of the Indian Ocean from the Atlantic for both global lineages. The haplotypes we confirmed for western Indian Ocean RMUs allow reinterpretation of previous mixed stock analysis and further suggest that contemporary migratory connectivity between the Indian and Atlantic Oceans occurs on a broader scale than previously hypothesized. This study represents a valuable model for conducting comprehensive international cooperative data management and research in marine ecology.