Digging deeper into East African human Y chromosome lineages

UYSD: a novel data repository accessible via public website for worldwide population frequencies of Y-SNP haplogroups

For decades, there has been scientific interest in the variation and geographic distribution of paternal lineages associated with the human Y chromosome. However, the relevant data have been dispersed across numerous publications, making it difficult to consolidate. Additionally, understanding the relationships between different variants, and the tools used to analyze them, have evolved over time, further complicating efforts to harmonize this information. The Universal Y-SNP Database (UYSD) marks a substantial advancement by providing a comprehensive and accessible platform for Y-SNP and haplogroup data from populations around the world. UYSD harmonizes diverse datasets into a unified repository, facilitating the exploration of global Y-chromosomal variation. The platform handles data generated with both high- and low-throughput technology and is compatible with the automated analysis software tool, Yleaf v3. Key functionalities include the ability to: i) visualize haplogroup distributions on an interactive world map, ii) estimate haplogroup frequencies in geographic regions with sparse data through interpolation, and iii) display detailed phylogenetic trees of Y-chromosomal haplogroups. Currently, UYSD encompasses data from over 6,600 males across 27 populations. This dataset largely aligns with known global Y-haplogroup patterns, but also reveals unexplored finer-scale geographic variations. While the present dataset is largely European-centered, UYSD is designed for ongoing expansion by the scientific community, aiming to include more global data and higher-resolution population sequencing data. The platform thus offers valuable insights into human genetic diversity and migration patterns, serving several fields of research such as: human population genetics, genetic anthropology, ancient DNA analysis and forensic genetics.

Application of Targeted Y-Chromosomal Capture Enrichment to Increase the Resolution of Native American Haplogroup Q

Y-chromosomal haplogroups and the Y-SNPs defining them are relevant for the exploration of male lineages, inference of paternal ancestry, and reconstruction of migration pathways, to name a few. Currently, over 300,000 Y-SNPs have been reported, defining 20 main haplogroups. However, ascertainment bias in the investigations has led to some haplogroups being overlooked, which hinders a representative depiction of certain populations and their migration events. For migration pattern analyses of the first settlers of the Americas, the Native American main founding lineage Q-M3 needs to be further investigated to allow clear genetic differentiation of individuals of different ethnogeographic origins. To increase the resolution within this haplogroup, a total of 7.45 Mb of the Y chromosome of 59 admixed South Americans of haplogroup Q was targeted for sequencing using hybridization capture enrichment. Data were combined with 218 publicly available sequences of Central and South Americans of haplogroup Q. After rigorous data processing, variants not meeting the quality criteria were excluded and 4128 reliable Y-SNPs were reported. A total of 2224 Y-SNPs had previously unknown positions in the phylogenetic tree, and 1291 of these are novel. The phylogenetic relationships between the Y-SNPs were established using the software SNPtotree in order to report a redesigned phylogenetic tree containing 300 branches, defined by 3400 Y-SNPs. The new tree introduces 117 previously undescribed branches and is the most comprehensive phylogenetic tree of the Native American haplogroup Q lineages to date. The 214 sequences were assigned to 135 different low- to high-resolution branches, while in the previous phylogenetic tree, only 195 sequences could be sorted into 14 low-resolution branches with the same quality criteria. The improved genetic differentiation of subhaplogroup Q-M3 has a great potential to resolve migration patterns of Native Americans.

Impact of patrilocality on contrasting patterns of paternal and maternal heritage in Central-West Africa

Despite their ancient past and high diversity, African populations are the least represented in human population genetic studies. In this study, uniparental markers (mtDNA and Y chromosome) were used to investigate the impact of sociocultural factors on the genetic diversity and inter-ethnolinguistic gene flow in the three major Nigerian groups: Hausa (n = 89), Yoruba (n = 135) and Igbo (n = 134). The results show a distinct history from the maternal and paternal perspectives. The three Nigerian groups present a similar substrate for mtDNA, but not for the Y chromosome. The two Niger-Congo groups, Yoruba and Igbo, are paternally genetically correlated with populations from the same ethnolinguistic affiliation. Meanwhile, the Hausa is paternally closer to other Afro-Asiatic populations and presented a high diversity of lineages from across Africa. When expanding the analyses to other African populations, it is observed that language did not act as a major barrier to female-mediated gene flow and that the differentiation of paternal lineages is better correlated with linguistic than geographic distances. The results obtained demonstrate the impact of patrilocality, a common and well-established practice in populations from Central-West Africa, in the preservation of the patrilineage gene pool and in the affirmation of identity between groups.

The Ancestry of Eastern Paraguay: A Typical South American Profile with a Unique Pattern of Admixture

Immigrants from diverse origins have arrived in Paraguay and produced important demographic changes in a territory initially inhabited by indigenous Guarani. Few studies have been performed to estimate the proportion of Native ancestry that is still preserved in Paraguay and the role of females and males in admixture processes. Therefore, 548 individuals from eastern Paraguay were genotyped for three marker sets: mtDNA, Y-SNPs and autosomal AIM-InDels. A genetic homogeneity was found between departments for each set of markers, supported by the demographic data collected, which showed that only 43% of the individuals have the same birthplace as their parents. The results show a sex-biased intermarriage, with higher maternal than paternal Native American ancestry. Within the native mtDNA lineages in Paraguay (87.2% of the total), most haplogroups have a broad distribution across the subcontinent, and only few are concentrated around the Paraná River basin. The frequency distribution of the European paternal lineages in Paraguay (92.2% of the total) showed a major contribution from the Iberian region. In addition to the remaining legacy of the colonial period, the joint analysis of the different types of markers included in this study revealed the impact of post-war migrations on the current genetic background of Paraguay.

Insights into the Middle Eastern paternal genetic pool in Tunisia: high prevalence of T-M70 haplogroup in an Arab population

To obtain refreshed insights into the paternal lineages of Tunisian populations, Y-chromosome diversity was assessed in two populations belonging to an Arab genealogical lineage, Kairouan and Wesletia, as well as in four Tunisian Andalusian populations, Testour, Slouguia, Qalaat-El-Andalous and El Alia. The Arabs from Kairouan revealed 73.47% of E-M81 and close affinities with Berber groups, indicating they are likely arabized Berbers, clearly differentiated from the Arabs from Wesletia, who harbored the highest frequency (71.8%) of the Middle Eastern component ever observed in North Africa. In the Tunisian Andalusians, the North African component largely prevailed, followed by the Middle Eastern contribution. Global comparative analysis highlighted the heterogeneity of Tunisian populations, among which, as a whole, dominated a set of lineages ascribed to be of autochthonous Berber origin (71.67%), beside a component of essentially Middle Eastern extraction (18.35%), and signatures of Sub-Saharan (5.2%), European (3.45%) and Asiatic (1.33%) contributions. The remarkable frequency of T-M70 in Wesletia (17.4%) prompted to refine its phylogeographic analysis, allowing to confirm its Middle Eastern origin, though signs of local evolution in Northern Africa were also detected. Evidence was clear on the ancient introduction of T lineages into the region, probably since Neolithic times associated to spread of agriculture.

Traditional food plants of the upper Aswa River catchment of northern Uganda-a cultural crossroads

BACKGROUND

In the parkland agroforestry system of northern Uganda, smallholder farming households rely on a diversity of plant species to fulfil their nutritional requirements, many of which also serve a range of medicinal, cultural, and livelihood functions. The purpose of the study was to assemble an inventory of indigenous plant species used as food in four districts within the Aswa River catchment of northern Uganda, and to document their utilization and management by rural communities.

METHODS

From July 1999 to August 2000, a series of 61 community-based focus group discussions on the utilization of plant biodiversity were conducted in the vernacular language at 34 sites in four districts of northern Uganda, with participation by key informants self-selected on basis of their technical knowledge and personal interest. Of these, 232 respondents subsequently contributed to a collection of herbarium specimens, which were submitted to the Makerere University Herbarium for identification. On receipt of each specimen collected, a structured interview was conducted to document the botanical, ecological, seasonal, and alimentary attributes of each identified taxon, and details of its processing and utilization by the community from which it was obtained. The data analysis was undertaken during 2019 and 2020, including statistical tests to assess the relative importance of the cited taxa using the Relative Importance Index (RI), and to determine the similarity of edible plant use between the four cultures using the Jaccard Index of similarity (JI).

RESULTS

Key informant interviews yielded 1347 use reports (URs) for 360 identified specimens of 88 indigenous edible plant species. The data describes patterns of use of indigenous edible plants of four cultures of the Aswa River catchment of northern Uganda. RI scores ranged from 0.93 to 0.11, with fruit trees occupying the top 25 taxa (RI 0.45 and above). Jaccard similarity scores ranged from 25.8% between Lango and Acholi, to 15.8% between Acholi and Ethur, indicating that cultural factors appear to be more significant than shared ancestry as determinants of cultural similarity of plant use.

CONCLUSIONS

The data constitute an inventory of on-farm plant species, including cultivated, semi-cultivated, and wild plants, integrated into a parkland agroforestry system in which useful trees and other plant species are sustained and managed under cultivation. Agricultural and on-farm plant biodiversity may be seen as a food security resource, and a nutritional buffer against increasing risks and stressors on low-input smallholder agriculture. Further studies should assess the intra-species biodiversity of these resources, with respect to farmer-valued traits and vernacular (folk) classification systems.

Middle eastern genetic legacy in the paternal and maternal gene pools of Chuetas

Chuetas are a group of descendants of Majorcan Crypto-Jews (Balearic Islands, Spain) who were socially stigmatized and segregated by their Majorcan neighbours until recently; generating a community that, although after the seventeenth century no longer contained Judaic religious elements, maintained strong group cohesion, Jewishness consciousness, and endogamy. Collective memory fixed 15 surnames as a most important defining element of Chueta families. Previous studies demonstrated Chuetas were a differentiated population, with a considerable proportion of their original genetic make-up. Genetic data of Y-chromosome polymorphism and mtDNA control region showed, in Chuetas’ paternal lineages, high prevalence of haplogroups J2-M172 (33%) and J1-M267 (18%). In maternal lineages, the Chuetas hallmark is the presence of a new sub-branching of the rare haplogroup R0a2m as their modal haplogroup (21%). Genetic diversity in both Y-chromosome and mtDNA indicates the Chueta community has managed to avoid the expected heterogeneity decrease in their gene pool after centuries of isolation and inbreeding. Moreover, the composition of their uniparentally transmitted lineages demonstrates a remarkable signature of Middle Eastern ancestry—despite some degree of host admixture—confirming Chuetas have retained over the centuries a considerable degree of ancestral genetic signature along with the cultural memory of their Jewish origin.

A Y-chromosomal survey of Ecuador's multi-ethnic population reveals new insights into the tri-partite population structure and supports an early Holocene age of the rare Native American founder lineage C3-MPB373

Ecuador is a multiethnic and pluricultural country with a complex history defined by migration and admixture processes. The present study aims to increase our knowledge on the Ecuadorian Native Amerindian groups and the unique South American Y-chromosome haplogroup C3-MPB373 through the analysis of up to 23 Y-chromosome STRs (Y-STRs) and several Y-SNPs in a sample of 527 Ecuadorians from 7 distinct populations and geographic areas, including Kichwa and non-Kichwa Native Amerindians, Mestizos and Afro-Ecuadorians. Our results reveal the presence of C3-MPB373 both in the Amazonian lowland Kichwa with frequencies up to 28 % and, for the first time, in notable proportions in Kichwa populations from the Ecuadorian highlands. The substantially higher frequencies of C3-MPB373 in the Amazonian lowlands found in Kichwa and Waorani individuals suggest a founder effect in that area. Notably, estimates for the time to the most recent common ancestor (TMRCA) in the range of 7.2-9.0 kya point to an ancient origin of the haplogroup and suggest an early Holocene expansion of C3-MPB373 into South America. Finally, the pairwise genetic distances (R_ST) separate the Kichwa Salasaka from all the other Native Amerindian and Ecuadorian groups, indicating a so far hidden diversity among the Kichwa-speaking populations and suggesting a more southern origin of this population. In sum, our study provides a more in-depth knowledge of the male genetic structure of the multiethnic Ecuadorian population, as well as a valuable reference dataset for forensic use.

New insights on intercontinental origins of paternal lineages in Northeast Brazil

Background

The current Brazilian population is the product of centuries of admixture between intercontinental founding groups. Although previous results have revealed a heterogeneous distribution of mitochondrial lineages in the Northeast region, the most targeted by foreign settlers during the sixteenth century, little is known about the paternal ancestry of this particular population. Considering historical records have documented a series of territorial invasions in the Northeast by various European populations, we aimed to characterize the male lineages found in Brazilian individuals in order to discover to what extent these migrations have influenced the present-day gene pool. Our approach consisted of employing four hierarchical multiplex assays for the investigation of 45 unique event polymorphisms in the non-recombining portion of the Y-chromosome of 280 unrelated men from several Northeast Brazilian states.

Results

Primary multiplex results allowed the identification of six major haplogroups, four of which were screened for downstream SNPs and enabled the observation of 19 additional lineages. Results reveal a majority of Western European haplogroups, among which R1b-S116* was the most common (63.9%), corroborating historical records of colonizations by Iberian populations. Nonetheless, F_ST genetic distances show similarities between Northeast Brazil and several other European populations, indicating multiple origins of settlers. Regarding Native American ancestry, our findings confirm a strong sexual bias against such haplogroups, which represented only 2.5% of individuals, highly contrasting previous results for maternal lineages. Furthermore, we document the presence of several Middle Eastern and African haplogroups, supporting a complex historical formation of this population and highlighting its uniqueness among other Brazilian regions.

Conclusions

We performed a comprehensive analysis of the major Y-chromosome lineages that form the most dynamic migratory region from the Brazilian colonial period. This evidence suggests that the ongoing entry of European, Middle Eastern, and African males in the Brazilian Northeast, since at least 500 years, was significantly responsible for the present-day genetic architecture of this population.

Forensic Y-SNP analysis beyond SNaPshot: High-resolution Y-chromosomal haplogrouping from low quality and quantity DNA using Ion AmpliSeq and targeted massively parallel sequencing

Y-chromosomal haplogroups assigned from male-specific Y-chromosomal single nucleotide polymorphisms (Y-SNPs) allow paternal lineage identification and paternal bio-geographic ancestry inference, both being relevant in forensic genetics. However, most previously developed forensic Y-SNP tools did not provide Y haplogroup resolution on the high level needed in forensic applications, because the limited multiplex capacity of the DNA technologies used only allowed the inclusion of a relatively small number of Y-SNPs. In a proof-of-principle study, we recently demonstrated that high-resolution Y haplogrouping is feasible via two AmpliSeq PCR analyses and simultaneous massively parallel sequencing (MPS) of 530 Y-SNPs allowing the inference of 432 Y-haplogroups. With the current study, we present a largely improved Y-SNP MPS lab tool that we specifically designed for the analysis of low quality and quantity DNA often confronted with in forensic DNA analysis. Improvements include i) Y-SNP marker selection based on the "minimal reference phylogeny for the human Y chromosome" (PhyloTree Y), ii) strong increase of the number of targeted Y-SNPs allowing many more Y haplogroups to be inferred, iii) focus on short amplicon length enabling successful analysis of degraded DNA, and iv) combination of all amplicons in a single AmpliSeq PCR and simultaneous sequencing allowing single DNA aliquot use. This new MPS tool simultaneously analyses 859 Y-SNPs and allows inferring 640 Y haplogroups. Preliminary forensic developmental validation testing revealed that this tool performs highly accurate, is sensitive and robust. We also provide a revised software tool for analysing the sequencing data produced by the new MPS lab tool including final Y haplogroup assignment. We envision the tools introduced here for high-resolution Y-chromosomal haplogrouping to determine a man's paternal lineage and/or paternal bio-geographic ancestry to become widely used in forensic Y-chromosome DNA analysis and other applications were Y haplogroup information from low quality / quantity DNA samples is required.

Paternal portrait of populations of the middle Magdalena River region (Tolima and Huila, Colombia): New insights on the peopling of Central America and northernmost South America

The valley of the Magdalena River is one of the main population pathways in Colombia. The gene pool and spatial configuration of human groups in this territory have been outlined throughout three historical stages: the Native pre-Hispanic world, Spanish colonization, and XIX century migrations. This research was designed with the goal of characterizing the diversity and distribution pattern of Y-chromosome lineages that are currently present in the Tolima and Huila departments (middle Magdalena River region). Historic cartography was used to identify the main geographic sites where the paternal lineages belonging to this area have gathered. Twelve municipalities were chosen, and a survey that included genealogical information was administered. Samples collected from 83 male volunteers were analyzed for 48 Y-SNPs and 17 Y-STRs. The results showed a highly diverse region characterized by the presence of 16 sublineages within the major clades R, Q, J, G, T and E and revealed that 93% (n = 77) of haplotypes were different. Among these haplogroups, European-specific R1b-M269 lineages were the most representative (57.83%), with six different subhaplogroups and 43 unique haplotypes. Native American paternal ancestry was also detected based on the presence of the Q1a2-M3*(xM19, M194, M199) and Q1a2-M346*(xM3) lineages. Interestingly, all Q1a2-M346*(xM3) samples (n = 7, with five different haplotypes) carried allele six at the DYS391 locus. This allele has a worldwide frequency of 0.169% and was recently associated with a new Native subhaplogroup. An in-depth phylogenetic analysis of these samples suggests the Tolima and Huila region to be the principal area in all Central and South America where this particular Native lineage is found. This lineage has been present in the region for at least 1,809 (+/- 0,5345) years.

Genetic structure and sex-biased gene flow in the history of southern African populations

Objectives

We investigated the genetic history of southern African populations with a special focus on their paternal history. We reexamined previous claims that the Y‐chromosome haplogroup E1b1b (E‐M293) was brought to southern Africa by pastoralists from eastern Africa, and investigated patterns of sex‐biased gene flow in southern Africa.

Materials and methods

We analyzed previously published complete mtDNA genome sequences and ∼900 kb of NRY sequences from 23 populations from Namibia, Botswana, and Zambia, as well as haplogroup frequencies from a large sample of southern African populations and 23 newly genotyped Y‐linked STR loci for samples assigned to haplogroup E1b1b.

Results

Our results support an eastern African origin for Y‐chromosome haplogroup E1b1b (E‐M293); however, its current distribution in southern Africa is not strongly associated with pastoralism, suggesting more complex demographic events and/or changes in subsistence practices in this region. The Bantu expansion in southern Africa had a notable genetic impact and was probably a rapid, male‐dominated expansion. Our finding of a significant increase in the intensity of the sex‐biased gene flow from north to south may reflect changes in the social dynamics between Khoisan and Bantu groups over time.

Conclusions

Our study shows that the population history of southern Africa has been complex, with different immigrating groups mixing to different degrees with the autochthonous populations. The Bantu expansion led to heavily sex‐biased admixture as a result of interactions between Khoisan females and Bantu males, with a geographic gradient which may reflect changes in the social dynamics between Khoisan and Bantu groups over time.

The peopling of the last Green Sahara revealed by high-coverage resequencing of trans-Saharan patrilineages

BACKGROUND

Little is known about the peopling of the Sahara during the Holocene climatic optimum, when the desert was replaced by a fertile environment.

RESULTS

In order to investigate the role of the last Green Sahara in the peopling of Africa, we deep-sequence the whole non-repetitive portion of the Y chromosome in 104 males selected as representative of haplogroups which are currently found to the north and to the south of the Sahara. We identify 5,966 mutations, from which we extract 142 informative markers then genotyped in about 8,000 subjects from 145 African, Eurasian and African American populations. We find that the coalescence age of the trans-Saharan haplogroups dates back to the last Green Sahara, while most northern African or sub-Saharan clades expanded locally in the subsequent arid phase.

CONCLUSIONS

Our findings suggest that the Green Sahara promoted human movements and demographic expansions, possibly linked to the adoption of pastoralism. Comparing our results with previously reported genome-wide data, we also find evidence for a sex-biased sub-Saharan contribution to northern Africans, suggesting that historical events such as the trans-Saharan slave trade mainly contributed to the mtDNA and autosomal gene pool, whereas the northern African paternal gene pool was mainly shaped by more ancient events.

Forensic use of Y-chromosome DNA: a general overview

The male-specific part of the human Y chromosome is widely used in forensic DNA analysis, particularly in cases where standard autosomal DNA profiling is not informative. A Y-chromosomal gene fragment is applied for inferring the biological sex of a crime scene trace donor. Haplotypes composed of Y-chromosomal short tandem repeat polymorphisms (Y-STRs) are used to characterise paternal lineages of unknown male trace donors, especially suitable when males and females have contributed to the same trace, such as in sexual assault cases. Y-STR haplotyping applied in crime scene investigation can (i) exclude male suspects from involvement in crime, (ii) identify the paternal lineage of male perpetrators, (iii) highlight multiple male contributors to a trace, and (iv) provide investigative leads for finding unknown male perpetrators. Y-STR haplotype analysis is employed in paternity disputes of male offspring and other types of paternal kinship testing, including historical cases, as well as in special cases of missing person and disaster victim identification involving men. Y-chromosome polymorphisms are applied for inferring the paternal bio-geographic ancestry of unknown trace donors or missing persons, in cases where autosomal DNA profiling is uninformative. In this overview, all different forensic applications of Y-chromosome DNA are described. To illustrate the necessity of forensic Y-chromosome analysis, the investigation of a prominent murder case is described, which initiated two changes in national forensic DNA legislation both covering Y-chromosome use, and was finally solved via an innovative Y-STR dragnet involving thousands of volunteers after 14 years. Finally, expectations for the future of forensic Y-chromosome DNA analysis are discussed.

Male Lineages in Brazil: Intercontinental Admixture and Stratification of the European Background

The non-recombining nature of the Y chromosome and the well-established phylogeny of Y-specific Single Nucleotide Polymorphisms (Y-SNPs) make them useful for defining haplogroups with high geographical specificity; therefore, they are more apt than the Y-STRs to detect population stratification in admixed populations from diverse continental origins. Different Y-SNP typing strategies have been described to address issues of population history and movements within geographic territories of interest. In this study, we investigated a set of 41 Y-SNPs in 1217 unrelated males from the five Brazilian geopolitical regions, aiming to disclose the genetic structure of male lineages in the country. A population comparison based on pairwise FST genetic distances did not reveal statistically significant differences in haplogroup frequency distributions among populations from the different regions. The genetic differences observed among regions were, however, consistent with the colonization history of the country. The sample from the Northern region presented the highest Native American ancestry (8.4%), whereas the more pronounced African contribution could be observed in the Northeastern population (15.1%). The Central-Western and Southern samples showed the higher European contributions (95.7% and 93.6%, respectively). The Southeastern region presented significant European (86.1%) and African (12.0%) contributions. The subtyping of the most frequent European lineage in Brazil (R1b1a-M269) allowed differences in the genetic European background of the five Brazilian regions to be investigated for the first time.

Y chromosome diversity in a linguistic isolate (Mirandese, NE Portugal)

OBJECTIVES

The purpose of this study was to genetically characterize the male lineages of people who speak Mirandese, an interesting case of a linguistic relict that can still be found in the municipality of Miranda do Douro, NE Portugal. This region lies within the area of the Leonese dialects, which are remnants of the Romance dialects spoken in the Kingdom of Leon currently grouped in the Astur-Leonese linguistic continuum. We intended to disclose affinities with surrounding populations, namely from Spain where the Astur-Leonese is also spoken.

METHODS

Eighty-eight unrelated males (58 from Miranda and 30 from Bragança, the broad Portuguese region where Miranda is located) were genotyped with the combined use of 17 Y chromosome short tandem repeats (Y-STRs) and a high resolution Y chromosome single nucleotide polymorphism (Y-SNPs) strategy. Moreover, 236 males from Miranda and neighboring regions, previously classified as R-M269, were also genotyped.

RESULTS

R-P312 was the most frequent haplogroup in the Mirandese, followed by J-12f2.1 and T-M70. The male lineages J-12f2.1 and T-M70 were also well represented, and both were shared with descendants of Sephardic Jews. No signs of diversity reduction were detected.

CONCLUSIONS

Mirandese speakers display a Y chromosome gene pool that shows a subtle differentiation from neighboring populations, mainly attributable to the assimilation of lineages ascribed to be of Jewish ancestry. Although not revealing signs of geographic/linguistic isolation, no clear affinities with other Astur-Leonese populations were detected. The results suggest that in Miranda language sharing is not accompanied by significant gene flow between populations from both sides of the political border. Am. J. Hum. Biol. 28:671-680, 2016. © 2016 Wiley Periodicals, Inc.

Mosaic maternal ancestry in the Great Lakes region of East Africa

The Great Lakes lie within a region of East Africa with very high human genetic diversity, home of many ethno-linguistic groups usually assumed to be the product of a small number of major dispersals. However, our knowledge of these dispersals relies primarily on the inferences of historical, linguistics and oral traditions, with attempts to match up the archaeological evidence where possible. This is an obvious area to which archaeogenetics can contribute, yet Uganda, at the heart of these developments, has not been studied for mitochondrial DNA (mtDNA) variation. Here, we compare mtDNA lineages at this putative genetic crossroads across 409 representatives of the major language groups: Bantu speakers and Eastern and Western Nilotic speakers. We show that Uganda harbours one of the highest mtDNA diversities within and between linguistic groups, with the various groups significantly differentiated from each other. Despite an inferred linguistic origin in South Sudan, the data from the two Nilotic-speaking groups point to a much more complex history, involving not only possible dispersals from Sudan and the Horn but also large-scale assimilation of autochthonous lineages within East Africa and even Uganda itself. The Eastern Nilotic group also carries signals characteristic of West-Central Africa, primarily due to Bantu influence, whereas a much stronger signal in the Western Nilotic group suggests direct West-Central African ancestry. Bantu speakers share lineages with both Nilotic groups, and also harbour East African lineages not found in Western Nilotic speakers, likely due to assimilating indigenous populations since arriving in the region ~3000 years ago.

Phylogeographic Refinement and Large Scale Genotyping of Human Y Chromosome Haplogroup E Provide New Insights into the Dispersal of Early Pastoralists in the African Continent

Haplogroup E, defined by mutation M40, is the most common human Y chromosome clade within Africa. To increase the level of resolution of haplogroup E, we disclosed the phylogenetic relationships among 729 mutations found in 33 haplogroup DE Y-chromosomes sequenced at high coverage in previous studies. Additionally, we dissected the E-M35 subclade by genotyping 62 informative markers in 5,222 samples from 118 worldwide populations. The phylogeny of haplogroup E showed novel features compared with the previous topology, including a new basal dichotomy. Within haplogroup E-M35, we resolved all the previously known polytomies and assigned all the E-M35* chromosomes to five new different clades, all belonging to a newly identified subhaplogroup (E-V1515), which accounts for almost half of the E-M35 chromosomes from the Horn of Africa. Moreover, using a Bayesian phylogeographic analysis and a single nucleotide polymorphism-based approach we localized and dated the origin of this new lineage in the northern part of the Horn, about 12 ka. Time frames, phylogenetic structuring, and sociogeographic distribution of E-V1515 and its subclades are consistent with a multistep demic spread of pastoralism within north-eastern Africa and its subsequent diffusion to subequatorial areas. In addition, our results increase the discriminative power of the E-M35 haplogroup for use in forensic genetics through the identification of new ancestry-informative markers.

Admixture and genetic diversity distribution patterns of non-recombining lineages of Native American ancestry in Colombian populations

Genetic diversity of present American populations results from very complex demographic events involving different types and degrees of admixture. Through the analysis of lineage markers such as mtDNA and Y chromosome it is possible to recover the original Native American haplotypes, which remained identical since the admixture events due to the absence of recombination. However, the decrease in the effective population sizes and the consequent genetic drift effects suffered by these populations during the European colonization resulted in the loss or under-representation of a substantial fraction of the Native American lineages. In this study, we aim to clarify how the diversity and distribution of uniparental lineages vary with the different demographic characteristics (size, degree of isolation) and the different levels of admixture of extant Native groups in Colombia. We present new data resulting from the analyses of mtDNA whole control region, Y chromosome SNP haplogroups and STR haplotypes, and autosomal ancestry informative insertion-deletion polymorphisms in Colombian individuals from different ethnic and linguistic groups. The results demonstrate that populations presenting a high proportion of non-Native American ancestry have preserved nevertheless a substantial diversity of Native American lineages, for both mtDNA and Y chromosome. We suggest that, by maintaining the effective population sizes high, admixture allowed for a decrease in the effects of genetic drift due to Native population size reduction and thus resulting in an effective preservation of the Native American non-recombining lineages.

Simultaneous analysis of hundreds of Y-chromosomal SNPs for high-resolution paternal lineage classification using targeted semiconductor sequencing

SNPs from the non-recombining part of the human Y chromosome (Y-SNPs) are informative to classify paternal lineages in forensic, genealogical, anthropological, and evolutionary studies. Although thousands of Y-SNPs were identified thus far, previous Y-SNP multiplex tools target only dozens of markers simultaneously, thereby restricting the provided Y-haplogroup resolution and limiting their applications. Here, we overcome this shortcoming by introducing a high-resolution multiplex tool for parallel genotyping-by-sequencing of 530 Y-SNPs using the Ion Torrent PGM platform, which allows classification of 432 worldwide Y haplogroups. Contrary to previous Y-SNP multiplex tools, our approach covers branches of the entire Y tree, thereby maximizing the paternal lineage classification obtainable. We used a default DNA input amount of 10 ng per reaction but preliminary sensitivity testing revealed positive results from as little as 100 pg input DNA. Furthermore, we demonstrate that sample pooling using barcodes is feasible, allowing increased throughput for lower per-sample costs. In addition to the wetlab protocol, we provide a software tool for automated data quality control and haplogroup classification. The unique combination of ultra-high marker density and high sensitivity achievable from low amounts of potentially degraded DNA makes this new multiplex tool suitable for a wide range of Y-chromosome applications.

Static and moving frontiers: the genetic landscape of Southern African Bantu-speaking populations

A consensus on Bantu-speaking populations being genetically similar has emerged in the last few years, but the demographic scenarios associated with their dispersal are still a matter of debate. The frontier model proposed by archeologists postulates different degrees of interaction among incoming agropastoralist and resident foraging groups in the presence of "static" and "moving" frontiers. By combining mitochondrial DNA and Y chromosome data collected from several southern African populations, we show that Bantu-speaking populations from regions characterized by a moving frontier developing after a long-term static frontier have larger hunter-gatherer contributions than groups from areas where a static frontier was not followed by further spatial expansion. Differences in the female and male components suggest that the process of assimilation of the long-term resident groups into agropastoralist societies was gender biased. Our results show that the diffusion of Bantu languages and culture in Southern Africa was a process more complex than previously described and suggest that the admixture dynamics between farmers and foragers played an important role in shaping the current patterns of genetic diversity.

Male lineage strata of Brazilian population disclosed by the simultaneous analysis of STRs and SNPs

Brazil has a large territory divided in five geographical regions harboring highly diverse populations that resulted from different degrees and modes of admixture between Native Americans, Europeans and Africans. In this study, a sample of 605 unrelated males was genotyped for 17 Y-STRs and 46 Y-SNPs aiming a deep characterization of the male gene pool of Rio de Janeiro and its comparison with other Brazilian populations. High values of Y-STR haplotype diversity (0.9999±0.0001) and Y-SNP haplogroup diversity (0.7589±0.0171) were observed. Population comparisons at both haplotype and haplogroup levels showed significant differences between Brazilian South Eastern and Northern populations that can be explained by differences in the proportion of African and Native American Y chromosomes. Statistical significant differences between admixed urban samples from the five regions of Brazil were not previously detected at haplotype level based on smaller size samples from South East, which emphasizes the importance of sample size to detected population stratification for an accurate interpretation of profile matches in kinship and forensic casework. Although not having an intra-population discrimination power as high as the Y-STRs, the Y-SNPs are more powerful to disclose differences in admixed populations. In this study, the combined analysis of these two types of markers proved to be a good strategy to predict population sub-structure, which should be taken into account when delineating forensic database strategies for Y chromosome haplotypes.

Association between Y haplogroups and autosomal AIMs reveals intra-population substructure in Bolivian populations

For the correct evaluation of the weight of genetic evidence in a forensic context, databases must reflect the structure of the population, with all possible groups being represented. Countries with a recent history of admixture between strongly differentiated populations are usually highly heterogeneous and sub-structured. Bolivia is one of these countries, with a high diversity of ethnic groups and different levels of admixture (among Native Americans, Europeans and Africans) across the territory. For a better characterization of the male lineages in Bolivia, 17 Y-STR and 42 Y-SNP loci were genotyped in samples from La Paz and Chuquisaca. Only European and Native American Y-haplogroups were detected, and no sub-Saharan African chromosomes were found. Significant differences were observed between the two samples, with a higher frequency of European lineages in Chuquisaca than in La Paz. A sample belonging to haplogroup Q1a3a1a1-M19 was detected in La Paz, in a haplotype background different from those previously found in Argentina. This result supports an old M19 North-south dispersion in South America, possibly via two routes. When comparing the ancestry of each individual assessed through his Y chromosome with the one estimated using autosomal AIMs, (a) increased European ancestry in individuals with European Y chromosomes and (b) higher Native American ancestry in the carriers of Native American Y-haplogroups were observed, revealing an association between autosomal and Y-chromosomal markers. The results of this study demonstrate that a sub-structure does exist in Bolivia at both inter- and intrapopulation levels, a fact which must be taken into account in the evaluation of forensic genetic evidence.

Colombia's racial crucible: Y chromosome evidence from six admixed communities in the Department of Bolivar

OBJECTIVES

To determine the African, European and Native-American paternal contributions in genetic samples from the Department of Bolivar (Colombia) with the aims of establishing (1) possible population substructures, and (2) the proportion of biological African heritage in admixed populations of European, Amerindian, and African descent.

METHODS

Y-SNPs were typed in samples from six communities, including Palenque (renowned for its African linguistic and cultural heritage).

RESULTS

Findings reveal a high diversity of Y-haplogroups. With the exception of Palenque, the sum of European male lineages uniformly exceeded 57%. In Palenque, African lineages accounted for 57.7% of its chromosomes, with European male lineages constituting a mere 38.5%. In Pinillos, a significant proportion (23.8%) of the chromosomes belongs to the Native American haplogroup Q1a3a*-M3. Genetic differentiation analyses reveal significant divergences in most pairwise comparisons among the Bolivar municipalities, and the same holds between Bolivar and other South American populations.

CONCLUSIONS

Heterogeneous patterns of admixture reveal a genetic substructure within the Department of Bolivar. On the paternal side, five out of the six communities studied exhibit a predominantly European gene pool. The exception is Palenque, where European input (38%) is more significant than we had expected.

AMY-tree: an algorithm to use whole genome SNP calling for Y chromosomal phylogenetic applications

Background

Due to the rapid progress of next-generation sequencing (NGS) facilities, an explosion of human whole genome data will become available in the coming years. These data can be used to optimize and to increase the resolution of the phylogenetic Y chromosomal tree. Moreover, the exponential growth of known Y chromosomal lineages will require an automatic determination of the phylogenetic position of an individual based on whole genome SNP calling data and an up to date Y chromosomal tree.

Results

We present an automated approach, ‘AMY-tree’, which is able to determine the phylogenetic position of a Y chromosome using a whole genome SNP profile, independently from the NGS platform and SNP calling program, whereby mistakes in the SNP calling or phylogenetic Y chromosomal tree are taken into account. Moreover, AMY-tree indicates ambiguities within the present phylogenetic tree and points out new Y-SNPs which may be phylogenetically relevant. The AMY-tree software package was validated successfully on 118 whole genome SNP profiles of 109 males with different origins. Moreover, support was found for an unknown recurrent mutation, wrong reported mutation conversions and a large amount of new interesting Y-SNPs.

Conclusions

Therefore, AMY-tree is a useful tool to determine the Y lineage of a sample based on SNP calling, to identify Y-SNPs with yet unknown phylogenetic position and to optimize the Y chromosomal phylogenetic tree in the future. AMY-tree will not add lineages to the existing phylogenetic tree of the Y-chromosome but it is the first step to analyse whole genome SNP profiles in a phylogenetic framework.

mtDNA variation in East Africa unravels the history of Afro-Asiatic groups

East Africa (EA) has witnessed pivotal steps in the history of human evolution. Due to its high environmental and cultural variability, and to the long-term human presence there, the genetic structure of modern EA populations is one of the most complicated puzzles in human diversity worldwide. Similarly, the widespread Afro-Asiatic (AA) linguistic phylum reaches its highest levels of internal differentiation in EA. To disentangle this complex ethno-linguistic pattern, we studied mtDNA variability in 1,671 individuals (452 of which were newly typed) from 30 EA populations and compared our data with those from 40 populations (2970 individuals) from Central and Northern Africa and the Levant, affiliated to the AA phylum. The genetic structure of the studied populations--explored using spatial Principal Component Analysis and Model-based clustering--turned out to be composed of four clusters, each with different geographic distribution and/or linguistic affiliation, and signaling different population events in the history of the region. One cluster is widespread in Ethiopia, where it is associated with different AA-speaking populations, and shows shared ancestry with Semitic-speaking groups from Yemen and Egypt and AA-Chadic-speaking groups from Central Africa. Two clusters included populations from Southern Ethiopia, Kenya and Tanzania. Despite high and recent gene-flow (Bantu, Nilo-Saharan pastoralists), one of them is associated with a more ancient AA-Cushitic stratum. Most North-African and Levantine populations (AA-Berber, AA-Semitic) were grouped in a fourth and more differentiated cluster. We therefore conclude that EA genetic variability, although heavily influenced by migration processes, conserves traces of more ancient strata.

Y-STR diversity and sex-biased gene flow among Caribbean populations

In the present study, we report, for the first time, the allele and haplotype frequencies of 17 Y-STR (Y-filer) loci in the populations of Haiti, Jamaica and the Bahamas (Abaco, Eleuthera, Exuma, Grand Bahama, Long Island and New Providence). This investigation was undertaken to assess the paternal genetic structure of the abovementioned Caribbean islands. A total of 607 different haplotypes were identified among the 691 males examined, of which 537 (88.5%) were unique. Haplotype diversities (HD) ranged from 0.989 in Long Island to 1.000 in Grand Bahama, with limited haplotype sharing observed among these Caribbean collections. Discriminatory capacity (DC) values were also high, ranging from 79.1% to 100% in Long Island and Grand Bahama, respectively, illustrating the capacity of this set of markers to differentiate between patrilineal related individuals within each population. Phylogenetic comparison of the Bahamian, Haitian and Jamaican groups with available African, European, East Asian and Native American populations reveals strong genetic ties with the continental African collections, a finding that corroborates our earlier work using autosomal STR and Y-chromosome binary markers. In addition, various degrees of sex-biased gene flow exhibiting disproportionately higher European paternal (as compared to autosomal) influences were detected in all Caribbean islands genotyped except for Abaco and Eleuthera. We attribute the presence or absence of asymmetric gene flow to unique, island specific demographic events and family structures.

Molecular dissection of the basal clades in the human Y chromosome phylogenetic tree

One hundred and forty-six previously detected mutations were more precisely positioned in the human Y chromosome phylogeny by the analysis of 51 representative Y chromosome haplogroups and the use of 59 mutations from literature. Twenty-two new mutations were also described and incorporated in the revised phylogeny. This analysis made it possible to identify new haplogroups and to resolve a deep trifurcation within haplogroup B2. Our data provide a highly resolved branching in the African-specific portion of the Y tree and support the hypothesis of an origin in the north-western quadrant of the African continent for the human MSY diversity.

The genetic landscape of Equatorial Guinea and the origin and migration routes of the Y chromosome haplogroup R-V88

Human Y chromosomes belonging to the haplogroup R1b1-P25, although very common in Europe, are usually rare in Africa. However, recently published studies have reported high frequencies of this haplogroup in the central-western region of the African continent and proposed that this represents a 'back-to-Africa' migration during prehistoric times. To obtain a deeper insight into the history of these lineages, we characterised the paternal genetic background of a population in Equatorial Guinea, a Central-West African country located near the region in which the highest frequencies of the R1b1 haplogroup in Africa have been found to date. In our sample, the large majority (78.6%) of the sequences belong to subclades in haplogroup E, which are the most frequent in Bantu groups. However, the frequency of the R1b1 haplogroup in our sample (17.0%) was higher than that previously observed for the majority of the African continent. Of these R1b1 samples, nine are defined by the V88 marker, which was recently discovered in Africa. As high microsatellite variance was found inside this haplogroup in Central-West Africa and a decrease in this variance was observed towards Northeast Africa, our findings do not support the previously hypothesised movement of Chadic-speaking people from the North across the Sahara as the explanation for these R1b1 lineages in Central-West Africa. The present findings are also compatible with an origin of the V88-derived allele in the Central-West Africa, and its presence in North Africa may be better explained as the result of a migration from the south during the mid-Holocene.

Afghanistan's ethnic groups share a Y-chromosomal heritage structured by historical events

Afghanistan has held a strategic position throughout history. It has been inhabited since the Paleolithic and later became a crossroad for expanding civilizations and empires. Afghanistan's location, history, and diverse ethnic groups present a unique opportunity to explore how nations and ethnic groups emerged, and how major cultural evolutions and technological developments in human history have influenced modern population structures. In this study we have analyzed, for the first time, the four major ethnic groups in present-day Afghanistan: Hazara, Pashtun, Tajik, and Uzbek, using 52 binary markers and 19 short tandem repeats on the non-recombinant segment of the Y-chromosome. A total of 204 Afghan samples were investigated along with more than 8,500 samples from surrounding populations important to Afghanistan's history through migrations and conquests, including Iranians, Greeks, Indians, Middle Easterners, East Europeans, and East Asians. Our results suggest that all current Afghans largely share a heritage derived from a common unstructured ancestral population that could have emerged during the Neolithic revolution and the formation of the first farming communities. Our results also indicate that inter-Afghan differentiation started during the Bronze Age, probably driven by the formation of the first civilizations in the region. Later migrations and invasions into the region have been assimilated differentially among the ethnic groups, increasing inter-population genetic differences, and giving the Afghans a unique genetic diversity in Central Asia.

Encephalocele in Uganda: ethnic distinctions in lesion location, endoscopic management of hydrocephalus, and survival in 110 consecutive children

OBJECT

This study characterizes the first clinical series of encephalocele (EC) from East or Central Africa, and is the largest reported from the African continent. The authors explored survival, the efficacy of primary endoscopic management of associated hydrocephalus, and ethnic differences in EC location.

METHODS

One hundred ten consecutive children presented to CURE Children's Hospital of Uganda for treatment of EC over a 9-year period. Clinical data, including patient demographic information, birth date, lesion type (sincipital, parietal, or occipital), operative data, and subsequent course had been entered prospectively into a clinical database. Home visits to update the status of those lost to follow-up were done when possible. With appropriate institutional approvals, the database was reviewed for this retrospective study. Two-tailed probability values calculated using the Fisher exact test were used to assess the significance of differences among groups, with p < 0.05 being considered significant. The Kaplan-Meier method was used for analysis of survival and treatment success probabilities.

RESULTS

There were 53 (48%) occipital, 33 (30%) sincipital, and 24 (22%) parietal lesions. Occipital lesions were significantly more common among children of Bantu origin (p = 0.02). Nilotes demonstrated a roughly equal distribution among sincipital, parietal, and occipital locations. The female/male ratio was 1.2, with no difference between EC types (range 1.0-1.4, p = 0.6-0.8). Of 110 patients, 108 (98%) underwent surgical repair at a median age of 1 month (mean 15.7 months), whereas 2 had treatment for hydrocephalus only. Wound revision was required in 13% of cases. Surgery-related mortality was 3%. One-year and 5-year survival rates were 87% (95% CI 0.79-0.93) and 61% (95% CI 0.51-0.70), respectively. Hydrocephalus required treatment in 32%, and was equally common among the 3 EC types. Thirteen patients were treated with combined endoscopic third ventriculostomy/choroid plexus cauterization (ETV/CPC) and 2 with ETV alone, whereas 18 patients received primary shunt placement. Predicted treatment success at 1 year was 79% for ETV or ETV/CPC (95% CI 0.50-0.94) and 47% for shunt placement (95% CI 0.24-0.71).

CONCLUSIONS

Analysis of this first EC series from this region suggests that sincipital lesions are 3 times more common in East than in West Africa. Occipital lesions predominate in patients of Bantu origin, but not among those of Nilotic descent. Hydrocephalus incidence was equally common among different EC types, and endoscopic treatment was more successful (79%) than shunting (47%) at 1 year. The 5-year mortality rate was similar to that for infants with myelomeningocele in Uganda, and more than twice that for their unaffected peers.

Y-chromosomal variation in sub-Saharan Africa: insights into the history of Niger-Congo groups

Technological and cultural innovations as well as climate changes are thought to have influenced the diffusion of major language phyla in sub-Saharan Africa. The most widespread and the richest in diversity is the Niger-Congo phylum, thought to have originated in West Africa ∼ 10,000 years ago (ya). The expansion of Bantu languages (a family within the Niger-Congo phylum) ∼ 5,000 ya represents a major event in the past demography of the continent. Many previous studies on Y chromosomal variation in Africa associated the Bantu expansion with haplogroup E1b1a (and sometimes its sublineage E1b1a7). However, the distribution of these two lineages extends far beyond the area occupied nowadays by Bantu-speaking people, raising questions on the actual genetic structure behind this expansion. To address these issues, we directly genotyped 31 biallelic markers and 12 microsatellites on the Y chromosome in 1,195 individuals of African ancestry focusing on areas that were previously poorly characterized (Botswana, Burkina Faso, Democratic Republic of Congo, and Zambia). With the inclusion of published data, we analyzed 2,736 individuals from 26 groups representing all linguistic phyla and covering a large portion of sub-Saharan Africa. Within the Niger-Congo phylum, we ascertain for the first time differences in haplogroup composition between Bantu and non-Bantu groups via two markers (U174 and U175) on the background of haplogroup E1b1a (and E1b1a7), which were directly genotyped in our samples and for which genotypes were inferred from published data using linear discriminant analysis on short tandem repeat (STR) haplotypes. No reduction in STR diversity levels was found across the Bantu groups, suggesting the absence of serial founder effects. In addition, the homogeneity of haplogroup composition and pattern of haplotype sharing between Western and Eastern Bantu groups suggests that their expansion throughout sub-Saharan Africa reflects a rapid spread followed by backward and forward migrations. Overall, we found that linguistic affiliations played a notable role in shaping sub-Saharan African Y chromosomal diversity, although the impact of geography is clearly discernible.

Y-chromosome R-M343 African lineages and sickle cell disease reveal structured assimilation in Lebanon

We have sought to identify signals of assimilation of African male lines in Lebanon by exploring the association of sickle cell disease (SCD) in Lebanon with Y-chromosome haplogroups that are informative of the disease origin and its exclusivity to the Muslim community. A total of 732 samples were analyzed, including 33 SCD patients from Lebanon genotyped for 28 binary markers and 19 short tandem repeats on the non-recombinant segment of the Y chromosome. Genetic organization was identified using populations known to have influenced the genetic structure of the Lebanese population, in addition to African populations with high incidence of SCD. Y-chromosome haplogroup R-M343 sub-lineages distinguish between sub-Saharan African and Lebanese Y chromosomes. We detected a limited penetration of SCD into Lebanese R-M343 carriers, restricted to Lebanese Muslims. We suggest that this penetration brought the sickle cell gene along with the African R-M343, probably with the Saharan caravan slave trade.