Development and validation of YARN: A novel SE-400 MPS kit for East Asian paternal lineage analysis

Y-chromosomal short tandem repeat polymorphisms (Y-STRs) and Y-chromosomal single nucleotide polymorphisms (Y-SNPs) are valuable genetic markers used in paternal lineage identification and population genetics. Currently, there is a lack of an effective panel that integrates Y-STRs and Y-SNPs for studying paternal lineages, particularly in East Asian populations. Hence, we developed a novel Y-chromosomal targeted panel called YARN (Y-chromosome Ancestry and Region Network) based on multiplex PCR and a single-end 400 massive parallel sequencing (MPS) strategy, consisting of 44 patrilineage Y-STRs and 260 evolutionary Y-SNPs. A total of 386 reactions were validated for the effectiveness and applicability of YARN according to SWGDAM validation guidelines, including sensitivity (with a minimum input gDNA of 0.125 ng), mixture identification (ranging from 1:1-1:10), PCR inhibitor testing (using substances such as 50 μM hematin, 100 μM hemoglobin, 100 μM humic acid, and 2.5 mM indigo dye), species specificity (successfully distinguishing humans from other animals), repeatability study (achieved 100% accuracy), and concordance study (with 99.91% accuracy for 1121 Y-STR alleles). Furthermore, we conducted a pilot study using YARN in a cohort of 484 Han Chinese males from Huaiji County, Zhaoqing City, Guangdong, China (GDZQHJ cohort). In this cohort, we identified 52 different Y-haplogroups and 73 different surnames. We found weak to moderate correlations between the Y-haplogroups, Chinese surnames, and geographical locations of the GDZQHJ cohort (with λ values ranging from 0.050 to 0.340). However, when we combined two different categories into a new independent variable, we observed stronger correlations (with λ values ranging from 0.617 to 0.754). Overall, the YARN panel, which combines Y-STR and Y-SNP genetic markers, meets forensic DNA quality assurance guidelines and holds potential for East Asian geographical origin inference and paternal lineage analysis.

Multiple founding paternal lineages inferred from the newly-developed 639-plex Y-SNP panel suggested the complex admixture and migration history of Chinese people

BACKGROUND

Non-recombining regions of the Y-chromosome recorded the evolutionary traces of male human populations and are inherited haplotype-dependently and male-specifically. Recent whole Y-chromosome sequencing studies have identified previously unrecognized population divergence, expansion and admixture processes, which promotes a better understanding and application of the observed patterns of Y-chromosome genetic diversity.

RESULTS

Here, we developed one highest-resolution Y-chromosome single nucleotide polymorphism (Y-SNP) panel targeted for uniparental genealogy reconstruction and paternal biogeographical ancestry inference, which included 639 phylogenetically informative SNPs. We genotyped these loci in 1033 Chinese male individuals from 33 ethnolinguistically diverse populations and identified 256 terminal Y-chromosomal lineages with frequency ranging from 0.0010 (singleton) to 0.0687. We identified six dominant common founding lineages associated with different ethnolinguistic backgrounds, which included O2a2b1a1a1a1a1a1a1-M6539, O2a1b1a1a1a1a1a1-F17, O2a2b1a1a1a1a1b1a1b-MF15397, O2a2b2a1b1-A16609, O1b1a1a1a1b2a1a1-F2517, and O2a2b1a1a1a1a1a1-F155. The AMOVA and nucleotide diversity estimates revealed considerable differences and high genetic diversity among ethnolinguistically different populations. We constructed one representative phylogenetic tree among 33 studied populations based on the haplogroup frequency spectrum and sequence variations. Clustering patterns in principal component analysis and multidimensional scaling results showed a genetic differentiation between Tai-Kadai-speaking Li, Mongolic-speaking Mongolian, and other Sinitic-speaking Han Chinese populations. Phylogenetic topology inferred from the BEAST and Network relationships reconstructed from the popART further showed the founding lineages from culturally/linguistically diverse populations, such as C2a/C2b was dominant in Mongolian people and O1a/O1b was dominant in island Li people. We also identified many lineages shared by more than two ethnolinguistically different populations with a high proportion, suggesting their extensive admixture and migration history.

CONCLUSIONS

Our findings indicated that our developed high-resolution Y-SNP panel included major dominant Y-lineages of Chinese populations from different ethnic groups and geographical regions, which can be used as the primary and powerful tool for forensic practice. We should emphasize the necessity and importance of whole sequencing of more ethnolinguistically different populations, which can help identify more unrecognized population-specific variations for the promotion of Y-chromosome-based forensic applications.

Development and evaluations of the ancestry informative markers of the VISAGE Enhanced Tool for Appearance and Ancestry

The VISAGE Enhanced Tool for Appearance and Ancestry (ET) has been designed to combine markers for the prediction of bio-geographical ancestry plus a range of externally visible characteristics into a single massively parallel sequencing (MPS) assay. We describe the development of the ancestry panel markers used in ET, and the enhanced analyses they provide compared to previous MPS-based forensic ancestry assays. As well as established autosomal single nucleotide polymorphisms (SNPs) that differentiate sub-Saharan African, European, East Asian, South Asian, Native American, and Oceanian populations, ET includes autosomal SNPs able to efficiently differentiate populations from Middle East regions. The ability of the ET autosomal ancestry SNPs to distinguish Middle East populations from other continentally defined population groups is such that characteristic patterns for this region can be discerned in genetic cluster analysis using STRUCTURE. Joint cluster membership estimates showing individual co-ancestry that signals North African or East African origins were detected, or cluster patterns were seen that indicate origins from central and Eastern regions of the Middle East. In addition to an augmented panel of autosomal SNPs, ET includes panels of 85 Y-SNPs, 16 X-SNPs and 21 autosomal Microhaplotypes. The Y- and X-SNPs provide a distinct method for obtaining extra detail about co-ancestry patterns identified in males with admixed backgrounds. This study used the 1000 Genomes admixed African and admixed American sample sets to fully explore these enhancements to the analysis of individual co-ancestry. Samples from urban and rural Brazil with contrasting distributions of African, European, and Native American co-ancestry were also studied to gauge the efficiency of combining Y- and X-SNP data for this purpose. The small panel of Microhaplotypes incorporated in ET were selected because they showed the highest levels of haplotype diversity amongst the seven population groups we sought to differentiate. Microhaplotype data was not formally combined with single-site SNP genotypes to analyse ancestry. However, the haplotype sequence reads obtained with ET from these loci creates an effective system for de-convoluting two-contributor mixed DNA. We made simple mixture experiments to demonstrate that when the contributors have different ancestries and the mixture ratios are imbalanced (i.e., not 1:1 mixtures) the ET Microhaplotype panel is an informative system to infer ancestry when this differs between the contributors.

Distribution of Y-chromosome haplogroups in Serbian population groups originating from historically and geographically significant distinct parts of the Balkan Peninsula

Our study enrolled 1200 Serbian males originating from three geographical regions in the Balkan Peninsula inhabited by Serbs: present-day Serbia, regions of Old Herzegovina and Kosovo and Metohija. These samples were genotyped using the combination of 23 Y-chromosomal short tandem repeats (Y-STRs) loci and 17 Ychromosomal single nucleotide polymorphisms (Y-SNPs) loci for the haplotype and haplogroup analysis in order to characterize in detail Y chromosome flow in the recent history. Serbia's borders have changed through history, forcing Serbs constantly to migrate to different regions of Balkan Peninsula. The most significant migration waves in the recent history towards present-day Serbia occurred from the regions of Old- Herzegovina and Kosovo and Metohija that lie in the south-west/south. High haplotype diversity and discrimination capacity were observed in all three datasets, with the highest number of unique haplotypes (381) and discrimination capacity (0.97) detected in the samples originating from the present-day Serbia. Haplogroup composition didn't differ significantly among datasets, with three dominant haplogroups (I-M170, E-P170 and R-M198), and haplogroup I-M170 being the most frequent in all three datasets. Haplogroup E-P170 was the second most dominant in the dataset originating from geographical region of Kosovo and Metohija, whereas haplogroup R-M198 was the second most prevalent in the dataset from historical region of Old Herzegovina. Based on the phylogenetic three for haplogroup I constructed within this study, haplogroup I2a1-P37.2 was the most dominant within all three datasets, especially in the dataset from historical region of Old Herzegovina, where 182 out of 400 samples were derived for SNP P37.2. Genetic distances between three groups of samples, evaluated by the Fst and Rst statistical values, and further visualized through multidimensional scaling plot, showed great genetic similarity between datasets from Old Herzegovina and present-day Serbia. Genetic difference in the haplogroup distribution and frequency between datasets from historical region of Old Herzegovina and from geographical region of Kosovo and Metohija was confirmed with highest Fst and Rst vaules. In this study we have distinguished genetic structure, diversity and haplogroup frequencies within 1200 Serbian males from three datasets, relationships among them as well as with other Balkan and European populations, which is useful for studying recent demographic history.

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.

Genetic reconstruction and phylogenetic analysis by 193 Y-SNPs and 27 Y-STRs in a Chinese Yi ethnic group

Yi is the seventh-largest ethnic group in China and features mountainous regional characteristics. The Liangshan Yi Autonomous Prefecture is the largest Yi agglomeration with isolated geographical conditions, profoundly impeding genetic communication. Here, we investigated 427 unrelated males of Liangshan from 193 Y-chromosome single nucleotide polymorphisms (Y-SNPs) and 27 Y-chromosome short tandem repeats (Y-STRs) to reveal the genetic structure and paternal phylogeny of the group. The haplogroup diversity reached 0.9169 with 46 different subhaplogroups by 193 Y-SNPs, and the haplotype diversity reached 0.9999 by 27 Y-STR loci. Multidimensional scaling (MDS), N-J tree, and Network were constructed to decipher and visualize the genetic relations between Yi and worldwide groups. Our results revealed: (1) the Network by Y-STRs and Y-SNPs showed the haplogroup D1a1a-M15 in the Liangshan Yi population was a ramification of Tibetan groups' expansion from west to east on the plateau; (2) the haplogroup distribution and the mismatch mutation analysis indicated the haplogroup O2a2b1a1a1a4a2-Z25929 of Liangshan Yi derived from manifold Southeast Asian immigrants; (3) a high-resolution Y-SNPs panel is vital to depict accurate paternal derivations and build an integrated and refining genetic structure of ethnic groups.

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.

Massively parallel sequencing of sex-chromosomal STRs in Saudi Arabia reveals patrilineage-associated sequence variants

Massively parallel sequencing (MPS) of forensic STRs has the potential to reveal additional allele diversity compared to conventional capillary electrophoresis (CE) typing strategies, but population studies are currently relatively few in number. The Verogen ForenSeq™ DNA Signature Prep Kit includes both Y-STRs and X-STRs among its targeted loci, and here we report the sequences of these loci, analysed using Verogen's ForenSeq™ Universal Analysis Software (UAS) v1.3 and STRait Razor v3.0, in a representative sample of 89 Saudi Arabian males. We identified 56 length variants (equivalent to CE alleles) and 75 repeat sequence sub-variants across the six X-STRs analysed; equivalent figures for the set of 24 Y-STRs were 147 and 192 respectively. We also observed two flanking sequence variants for the X-, and six for the Y-STRs. Recovery of sequence data and concordance with CE data (where available) across the tested loci was good, though rare flanking variation affected interpretation and allele calling at DYF387S1 and DXS7132. Examination of flanking sequences of the Y-STRs revealed five SNPs (L255, M4790, BY7692, Z16708 and S17543) previously shown to define specific haplogroups by Y-chromosome sequencing. These define Y-haplogroups in 62 % of our sample, a proportion that increases to 91 % when haplogroup-associated repeat-sequence motifs are also considered. A population-level comparison of the Saudi Arabian X-STRs with a global sample showed our dataset to be part of a large cluster of populations of West Eurasian and Middle Eastern origin.

Development and validation of a custom panel including 183 Y-SNPs for Chinese Y-chromosomal haplogroups dissection using a MALDI-TOF MS system

Y-chromosome SNP haplogroups exhibit geographic structuring in many populations around the world. Therefore, Y-chromosome haplogroups have been widely used to infer paternal biogeographical ancestry and high-resolution paternal lineage classification. In the present study, we designed a customized panel containing 183 Y-SNPs based on previous studies and evaluated the genotyping performance and repeatability, concordance, sensitivity, and ability of analysing case-type samples using a MALDI-TOF MS platform. The average call rate for duplicate typing of any one SNP in the panel was 97.0%. In the concordance and accuracy study, the results of haplogroup designation obtained from MALDI-TOF MS platform were fully consistent with those obtained from the next-generation sequencing (NGS) platform. The optimal amount of template DNA in the PCR seemed to be 10 ng. However, if less DNA (≥156.25 pg) was available, it was still possible to obtain meaningful haplogroup information. For the application part, this panel could be applied for the detection of blood, semen, and buccal swabs samples. Particularly, blood stain on FTA card samples could be dissected by direct PCR amplification on the MALDI-TOF MS platform. Besides, 371 unrelated male individuals from four Chinese ethnic groups (Han, Hui, Mongolian, and Kazak) were detected using this panel. Total 78 terminal haplogroups were found and the haplogroup diversity was 0.933576. The results demonstrate that this panel could be an accurate, fast, and cost-effective method for database construction where the amount of sample material is less of a concern and when the cost of the assay is taken into consideration.

Y Haplogroup Diversity of the Dominican Republic: Reconstructing the Effect of the European Colonization and the Trans-Atlantic Slave Trades

The Dominican Republic is one of the two countries on the Hispaniola island, which is part of the Antilles. Hispaniola was affected by the European colonization and massive deportation of African slaves since the XVI century and these events heavily shaped the genetic composition of the present-day population. To shed light about the effect of the European rules, we analyzed 92 single nucleotide polymorphisms on the Y chromosome in 182 Dominican individuals from three different locations. The Dominican Y haplogroup composition was characterized by an excess of northern African/European lineages (59%), followed by the African clades (38%), whereas the Native-American lineages were rare (3%). The comparison with the mitochondrial DNA variability, dominated by African clades, revealed a sex-biased admixture pattern, in line with the colonial society dominated by European men. When other Caribbean and non-Caribbean former colonies were also considered, we noted a difference between territories under a Spanish rule (like the Dominican Republic) and British/French rule, with the former characterized by an excess of European Y lineages reflecting the more permissive Iberian legislation about mixed people and slavery. Finally, we analyzed the distribution in Africa of the Dominican lineages with a putative African origin, mainly focusing on central and western Africa, which were the main sources of African slaves. We found that most (83%) of the African lineages observed in Santo Domingo have a central African ancestry, suggesting that most of the slaves were deported from regions.

Ethnic fragmentation and degree of urbanization strongly affect the discrimination power of Y-STR haplotypes in central Sahel

Y chromosome short tandem repeats (Y-STRs) are commonly used to identify male lineages for investigative and judicial purposes and could represent the only source of male-specific genetic information from unbalanced female-male mixtures. The Yfiler Plus multiplex, which includes twenty conventional and seven rapidly-mutating Y-STRs, represents the most discriminating patrilineal system commercially available to date. Over the past five years, this multiplex has been used to analyze several Eurasian populations, with a reported discrimination capacity (DC) approaching or corresponding to the highest possible value. However, despite the inclusion of rapidly mutating Y-STRs, extensive haplotype sharing was still reported for some African populations due to a number of different factors affecting the effective population size. In the present study, we analyzed 27 Y-STRs included in the Yfiler Plus multiplex and 82 Y-SNPs in central Sahel (northern Cameroon and western Chad), an African region characterized by a strong ethnic fragmentation and linguistic diversity. We evaluated the effects of population sub-structuring on genetic diversity by stratifying a sample composed of 431 males according to their ethnicity (44 different ethnic groups) and urbanization degree (four villages and four towns). Overall, we observed a low discrimination capacity (DC = 0.90), with 71 subjects (16.5 %) sharing 27 Y-STR haplotypes. Haplotype sharing was essentially limited to subjects with the same binary haplogroup, coming from the same location and belonging to the same ethnic group. Haplotype sharing was much higher in rural areas (average DC = 0.83) than urban settlements (average DC = 0.96) with a significant correlation between DC and census size (r = 0.89; p = 0.003). Notably, we found that genetic differentiation between villages from the same country (Φ_ST = 0.14) largely exceeded that found among countries (Φ_ST = 0.02). These findings have important implications for the choice of the appropriate reference population database to evaluate the statistical relevance of forensic Y-haplotype matches.

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.

The paternal genetic structure of Jingpo and Dai in southwest China

Yunnan province harbours substantial genetic, cultural and linguistic diversity, with the largest number of Aborigines in China, but the relationship among these Aborigines remains enigmatic. This study genotyped 45 Y chromosomal single nucleotide polymorphisms (SNPs) of 500 males from two aboriginal cross-border populations, Jingpo and Dai, from Dehong, Yunnan. It is reported that Haplogroup O2a2b1a1-M117 is the dominant lineage in both Jingpo and Dai. The Jingpo people show affinity with Tibeto-Burman speaking populations with a relatively high frequency of Haplogroup D-M174, and the Dai people are generally genetically similar with Tai-Kadai speaking populations with high frequencies of Haplogroup O1a-M119 and O1b1a1a-M95, which is consistent with their language classification.

Identification and diversity of Y-chromosome haplotypes in Qinghai yak populations

The aim of the present study was to perform a preliminary analysis of the characterization and diversity of Y-chromosome haplotypes/haplogroups in yak of Qinghai Province, China. A total of 322 male yaks from nine populations belonging to three officially recognized breeds (Gaoyuan, Huanhu and Datong) were sampled. Animals were genotyped using six previously reported Y-SNPs present in the SRY, USP9Y, UTY, AMELY and OFD1Y genes and four new Y-SNPs in the OFD1Y gene (g.569A>C, g.578A>C, g.608G>T and g.653G>C) identified in this study. Seven Y-chromosome haplotypes (H1-H7) were identified according to the combination of the 10 Y-SNPs. H1, H2 and H6 were the most common and shared haplotypes across all yak populations/breeds. Private haplotypes H3 and H7 were detected in the Datong breed; H4 in Guoleimude, Qumalai, Qilian, Tianjun and Ganglong populations; and H5 in Qumalai of Gaoyuan breed. Haplotype clustering and network analyses inferred two haplogroups, Y1 and Y2, indicating two divergent lineages of paternal origins of Qinghai yak. The analysis of molecular variance showed a significant difference among individuals (P < 0.0001) with more than 93% of the total genetic variation present within populations, suggesting a weak genetic structure among Qinghai yak populations. The overall Y-haplotype diversity was 0.538 ± 0.028, showing a relatively high diversity in Qinghai yak. The Gaoyuan and Datong breeds had similar haplotype diversities (0.547 ± 0.030 and 0.553 ± 0.083, respectively), which were higher than that of the Huanhu breed (0.441 ± 0.098). Our results support the conservation and sustainable use of unique yak genetic resources in Qinghai.

Effective resolution of the Y chromosome sublineages of the Iberian haplogroup R1b-DF27 with forensic purposes

Single-nucleotide polymorphisms (SNPs) found within the non-recombining region of the Y chromosome (NRY) represent a powerful tool in forensic genetics for inferring the paternal ancestry of a vestige and complement the determination of biogeographical origin in combination with other markers like AIMs. In the present study, we introduce a panel of 15 Y-SNPs for a fine-resolution subtyping of the haplogroup R1b-DF27, in a single minisequencing reaction. This is the first minisequencing panel that allows a fine subtyping of R1b-DF27, which displays high frequencies in Iberian and Iberian-influenced populations. This panel includes subhaplogroups of DF27 that display moderate geographical differentiation, of interest to link a sample with a specific location of the Iberian Peninsula or with Iberian ancestry. Conversely, part of the intricacy of a new minisequencing panel is to have all the included variants available to test the effectiveness of the analysis method. We have overcome the absence of the least common variants through site-directed mutagenesis. Overall, the results show that our panel is a robust and effective method for subtyping R1b-DF27 lineages from a minimal amount of DNA, and its high resolution enables to improve male lineage discrimination in Iberian and Southwest European descent individuals. The small length of the amplicons and its reproducibility makes this assay suitable for forensic and population genetics purposes.

Genetic ancestry of families of putative Inka descent

This study focuses on the descendants of the royal Inka family. The Inkas ruled Tawantinsuyu, the largest pre-Columbian empire in South America, which extended from southern Colombia to central Chile. The origin of the royal Inkas is currently unknown. While the mummies of the Inka rulers could have been informative, most were destroyed by Spaniards and the few remaining disappeared without a trace. Moreover, no genetic studies have been conducted on present-day descendants of the Inka rulers. In the present study, we analysed uniparental DNA markers in 18 individuals predominantly from the districts of San Sebastian and San Jerónimo in Cusco (Peru), who belong to 12 families of putative patrilineal descent of Inka rulers, according to documented registries. We used single-nucleotide polymorphisms and short tandem repeat (STR) markers of the Y chromosome (Y-STRs), as well as mitochondrial DNA D-loop sequences, to investigate the paternal and maternal descent of the 18 alleged Inka descendants. Two Q-M3* Y-STR clusters descending from different male founders were identified. The first cluster, named AWKI-1, was associated with five families (eight individuals). By contrast, the second cluster, named AWKI-2, was represented by a single individual; AWKI-2 was part of the Q-Z19483 sub-lineage that was likely associated with a recent male expansion in the Andes, which probably occurred during the Late Intermediate Period (1000-1450 AD), overlapping the Inka period. Concerning the maternal descent, different mtDNA lineages associated with each family were identified, suggesting a high maternal gene flow among Andean populations, probably due to changes in the last 1000 years.

Defining Y-SNP variation among the Flemish population (Western Europe) by full genome sequencing

Y-chromosomal single nucleotide polymorphisms (Y-SNPs) represent a powerful tool in forensic research and casework, especially for inferring paternal ancestry of unknown perpetrators and unidentified bodies. However, the wealth of recently discovered Y-SNPs, the 'jungle' of different evolutionary lineage trees and nomenclatures, and the lack of population-wide data of many phylogenetically mapped Y-SNPs, limits the use of Y-SNPs in routine forensic approaches. Recently, a concise reference phylogeny of the human Y chromosome, the 'Minimal Reference Y-tree', was introduced aiming to provide a stable phylogeny with optimal global discrimination capacity by including the most resolving Y-SNPs. Here, we obtained a representative sample of 270 whole-genome sequences (WGS) to grasp the Y-SNP variation within the autochthonous Flemish population (Belgium, Western Europe) according to this reference Y-tree. The high quality of the Y-SNP calling was guaranteed for the WGS sample as well as its representativeness for the Flemish population based on the comparison of the main haplogroup frequencies with those from earlier studies on Flanders and the Netherlands. The 270 Flemish Y chromosomes were assigned to 98 different sub-haplogroups of the Minimal Reference Y-tree, showing its high potential of discrimination and confirming the spectrum of evolutionary lineages within Western Europe in general and within Flanders in particular. The full database with all Y-SNP calls of the Flemish sample is public available for future updates including forensic and population genetic studies. New initiatives to categorise Y-SNP variation in other populations according to the reference phylogeny of the Y chromosome are highly encouraged for forensic applications. Recommendations to realise such future population sample sets are discussed based on this study.

The Genetic History of Peruvian Quechua-Lamistas and Chankas: Uniparental DNA Patterns among Autochthonous Amazonian and Andean Populations

This study focuses on the genetic history of the Quechua‐Lamistas, inhabitants of the Lamas Province in the San Martin Department, Peru, who speak their own distinct variety of the Quechua family of languages. It has been suggested that different pre‐Columbian ethnic groups from the Peruvian Amazonia, like the Motilones or “shaven heads”, assimilated the Quechua language and then formed the current native population of Lamas. However, many Quechua‐Lamistas claim to be direct descendants of the Chankas, a famous pre‐Columbian indigenous group that escaped from Inca rule in the Andes. To investigate the Quechua‐Lamistas and Chankas’ ancestries, we compared uniparental genetic profiles (17 STRs of Q‐M3 Y‐chromosome and mtDNA complete control region haplotypes) among autochthonous Amazonian and Andean populations from Peru, Bolivia and Ecuador. The phylogeographic and population genetic analyses indicate a fairly heterogeneous ancestry for the Quechua‐Lamistas, while they are closely related to their neighbours who speak Amazonian languages, presenting no direct relationships with populations from the region where the ancient Chankas lived. On the other hand, the genetic profiles of self‐identified Chanka descendants living in Andahuaylas (located in the Apurimac Department, Peru, in the Central Andes) were closely related to those living in Huancavelica and the assumed Chanka Confederation area before the Inca expansion.

Genetic population study of Y-chromosome markers in Benin and Ivory Coast ethnic groups

Ninety-six single nucleotide polymorphisms (SNPs) and seventeen short tandem repeat (STRs) were investigated on the Y-chromosome of 288 unrelated healthy individuals from populations in Benin (Bariba, Yoruba, and Fon) and the Ivory Coast (Ahizi and Yacouba). We performed a multidimensional scaling analysis based on FST and RST genetic distances using a large extensive database of sub-Saharan African populations. There is more genetic homogeneity in Ivory Coast populations compared with populations from Benin. Notably, the Beninese Yoruba are significantly differentiated from neighbouring groups, but also from the Yoruba from Nigeria (FST>0.05; P<0.01). The Y-chromosome dataset presented here provides new valuable data to understand the complex genetic diversity and human male demographic events in West Africa.

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.

Testing Central and Inner Asian admixture among contemporary Hungarians

Historically, the Carpathian Basin was the final destination for many nomadic peoples who migrated westward from Inner and Central Asia towards Europe. Proto-Hungarians (Steppe Magyars) were among those who came from the East, the Eurasian Steppe in the early middle ages. In order to detect the paternal genetic contribution from nomadic Steppe tribes, we tested 966 samples from Central Asian (Uzbekistan, Kazakhstan), Inner Asian (Mongolians and Buryats in Mongolia) and Hungarian-speaking European (Hungarian, Sekler and Csango) populations. We constructed median-joining networks of certain haplogroups in Hungarian-speaking European, and Altaic-speaking Central and Inner Asian populations. We estimated that the possible paternal genetic contribution from the above described populations among contemporary Hungarian speaking populations ranged between 5% and 7.4%. It is lowest among Hungarians from Hungary (5.1%), while higher among Hungarian-speaking groups in Romania, notably Sekler (7.4%) and Csango (6.3%). However, these results represent only an upper limit. Actual Central/Inner Asian admixture might be somewhat lower as some of the related lineages may have come from a common third source. The main haplogroups responsible for the Central/Inner Asian admixture among Hungarians are J2*-M172 (xM47, M67, M12), J2-L24, R1a-Z93; Q-M242 and E-M78. Earlier studies showed very limited Uralic genetic influence among Hungarians, and based on the present study, Altaic/Turkic genetic contribution is also not significant, although significantly higher than the Uralic one. The conclusion of this study is that present-day Hungarian speakers are genetically very similar to neighbouring populations, isolated Hungarian speaking groups having relatively higher presence of Central and Inner Asian genetic elements. At the same time, the reliable historical and genetic conclusions require an extension of the study to a significantly larger database with deep haplogroup resolution, including ancient DNA data.

Towards a consensus Y-chromosomal phylogeny and Y-SNP set in forensics in the next-generation sequencing era

Currently, several different Y-chromosomal phylogenies and haplogroup nomenclatures are presented in scientific literature and at conferences demonstrating the present diversity in Y-chromosomal phylogenetic trees and Y-SNP sets used within forensic and anthropological research. This situation can be ascribed to the exponential growth of the number of Y-SNPs discovered due to mostly next-generation sequencing (NGS) studies. As Y-SNPs and their respective phylogenetic positions are important in forensics, such as for male lineage characterization and paternal bio-geographic ancestry inference, there is a need for forensic geneticists to know how to deal with these newly identified Y-SNPs and phylogenies, especially since these phylogenies are often created with other aims than to carry out forensic genetic research. Therefore, we give here an overview of four categories of currently used Y-chromosomal phylogenies and the associated Y-SNP sets in scientific research in the current NGS era. We compare these categories based on the construction method, their advantages and disadvantages, the disciplines wherein the phylogenetic tree can be used, and their specific relevance for forensic geneticists. Based on this overview, it is clear that an up-to-date reduced tree with a consensus Y-SNP set and a stable nomenclature will be the most appropriate reference resource for forensic research. Initiatives to reach such an international consensus are therefore highly recommended.

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.

Maternal admixture and population structure in Mexican-Mestizos based on mtDNA haplogroups

The maternal ancestry (mtDNA) has important applications in different research fields, such as evolution, epidemiology, identification, and human population history. This is particularly interesting in Mestizos, which constitute the main population in Mexico (∼93%) resulting from post-Columbian admixture between Spaniards, Amerindians, and African slaves, principally. Consequently, we conducted minisequencing analysis (SNaPshot) of 11 mitochondrial single-nucleotide polymorphisms in 742 Mestizos of 10 populations from different regions in Mexico. The predominant maternal ancestry was Native American (92.9%), including Haplogroups A, B, C, and D (47, 23.7, 15.9, and 6.2%, respectively). Conversely, European and African ancestries were less frequent (5.3 and 1.9%, respectively). The main characteristics of the maternal lineages observed in Mexican-Mestizos comprised the following: 1) contrasting geographic gradient of Haplogroups A and C; 2) increase of European lineages toward the Northwest; 3) low or absent, but homogeneous, African ancestry throughout the Mexican territory; 4) maternal lineages in Mestizos roughly represent the genetic makeup of the surrounding Amerindian groups, particularly toward the Southeast, but not in the North and West; 5) continuity over time of the geographic distribution of Amerindian lineages in Mayas; and 6) low but significant maternal population structure (FST  = 2.8%; P = 0.0000). The average ancestry obtained from uniparental systems (mtDNA and Y-chromosome) in Mexican-Mestizos was correlated with previous ancestry estimates based on autosomal systems (genome-wide single-nucleotide polymorphisms and short tandem repeats). Finally, the comparison of paternal and maternal lineages provided additional information concerning the gender bias admixture, mating patterns, and population structure in Mestizos throughout the Mexican territory.

Increasing phylogenetic resolution still informative for Y chromosomal studies on West-European populations

Many Y-chromosomal lineages which are defined in the latest phylogenetic tree of the human Y chromosome by the Y Chromosome Consortium (YCC) in 2008 are distributed in (Western) Europe due to the fact that a large number of phylogeographic studies focus on this area. Therefore, the question arises whether newly discovered polymorphisms on the Y chromosome will still be interesting to study Western Europeans on a population genetic level. To address this question, the West-European region of Flanders (Belgium) was selected as study area since more than 1000 Y chromosomes from this area have previously been genotyped at the highest resolution of the 2008 YCC-tree and coupled to in-depth genealogical data. Based on these data the temporal changes of the population genetic pattern over the last centuries within Flanders were studied and the effects of several past gene flow events were identified. In the present study a set of recently reported novel Y-SNPs were genotyped to further characterize all those Flemish Y chromosomes that belong to haplogroups G, R-M269 and T. Based on this extended Y-SNP set the discrimination power increased drastically as previous large (sub-)haplogroups are now subdivided in several non-marginal groups. Next, the previously observed population structure within Flanders appeared to be the result of different gradients of independent sub-haplogroups. Moreover, for the first time within Flanders a significant East-West gradient was observed in the frequency of two R-M269 lineages, and this gradient is still present when considering the current residence of the DNA donors. Our results thus suggest that an update of the Y-chromosomal tree based on new polymorphisms is still useful to increase the discrimination power based on Y-SNPs and to study population genetic patterns in more detail, even in an already well-studied region such as Western Europe.

Updating the Y-chromosomal phylogenetic tree for forensic applications based on whole genome SNPs

The Y-chromosomal phylogenetic tree has a wide variety of important forensic applications and therefore it needs to be state-of-the-art. Nevertheless, since the last 'official' published tree many publications reported additional Y-chromosomal lineages and other phylogenetic topologies. Therefore, it is difficult for forensic scientists to interpret those reports and use an up-to-date tree and corresponding nomenclature in their daily work. Whole genome sequencing (WGS) data is useful to verify and optimise the current phylogenetic tree for haploid markers. The AMY-tree software is the first open access program which analyses WGS data for Y-chromosomal phylogenetic applications. Here, all published information is collected in a phylogenetic tree and the correctness of this tree is checked based on the first large analysis of 747 WGS samples with AMY-tree. The obtained result is one phylogenetic tree with all peer-reviewed reported Y-SNPs without the observed recurrent and ambiguous mutations. Nevertheless, the results showed that currently only the genomes of a limited set of Y-chromosomal (sub-)haplogroups is available and that many newly reported Y-SNPs based on WGS projects are false positives, even with high sequencing coverage methods. This study demonstrates the usefulness of AMY-tree in the process of checking the quality of the present Y-chromosomal tree and it accentuates the difficulties to enlarge this tree based on only WGS methods.