Evaluating methods to correct for population stratification when estimating paternity indexes

Thousands of years of Malay and Chinese population history in Indonesia and its implication on Paternity Index in DNA paternity testing

The existence of the Chinese population in the predominantly Malay population in Indonesia can be traced back thousands of years, and it has been suspected that it played an essential role in the history of the Malay population origin in Maritime South East Asia. With the fact that the Malay-Indonesian population is currently predominant compared to the Chinese population in Indonesia (Chinese-Indonesian), the selection of the origin of the STRs allele frequency panel population becomes an issue in DNA profiling, including in paternity testing. This study analyses the genetic relationship between the Chinese-Indonesian and Malay-Indonesian populations and how this affects the Paternity Index (PI) ​​calculation in paternity test cases. The study of the relationship between populations was carried out using neighbour-joining (NJ) tree analysis and multidimensional scaling (MDS) on the allele frequency panel of 19 autosomal STRs loci of Malay-Indonesian (n = 210) and Chinese-Indonesian (n = 78) populations. Four population groups were used as references: Malay-Malaysian, Filipino, Chinese, and Caucasian. An MDS analysis was also performed based on the pairwise F_ST calculation. The combined Paternity Index (CPI) calculation was carried out on 132 paternity cases from the Malay-Indonesian population with inclusive results using a panel of allele frequencies from the six populations. The pairwise F_ST MDS indicates a closer relationship between the Chinese-Indonesian and Malay-Indonesian compared to the Chinese population, which is in line with the CPIs comparison test. The outcome suggests that the alternative use of allele frequency database between Malay-Indonesian and Chinese-Indonesian for CPI calculations is not very influential. These results can also be considered in studying the extent of genetic assimilation between the two populations. In addition, these results support the robustness claim of multivariate analysis to represent phenomena that phylogenetic analyses may not be able to demonstrate, especially for massive panel data.

Isolation of novel microsatellite markers and their application for genetic diversity and parentage analyses in sika deer

Every part of the sika deer (Cervus nippon) body is valuable traditional Chinese medicine. And sika deer is the most important semi-domestic medicinal animal that is widely bred in Jilin province northeast of China. But few studies had been conducted to characterize the microsatellite markers derived from sika deer. We firstly used IlluminaHiSeq™2500 sequencing technology obtained 125Mbp genomic data of sika deer. Using microsatellite identification tool (MISA), 22,479 microsatellites were identified. From these data, 100 potential primers were selected for further polymorphic validation, finally, 76 primer pairs were successfully amplified and 29 primer pairs were found to be obvious polymorphic in 8 different individuals. Using those polymorphic microsatellite markers, we analyzed the genetic diversity of Jilin sika deer population. The mean number of alleles of the 29 loci is 9.31 based on genotyping blood DNA from 96 Jilin sika deer; The mean expected heterozygosity and polymorphic information content (PIC) value of the 29 loci is 0.72 and 0.68 respectively, and among which 26 loci are highly polymorphic (PIC>0.50). According to the electrophoretic results and PIC value of these 29 loci, 10 loci with combined paternity exclusion probabilities>99.99% were selected to use in parentage verification for 16 sika deer. All the offspring of a family could be successfully assigned to their biological father. These microsatellite markers generated in this study could greatly facilitate future studies of molecular breeding in sika deer.

The geographic mosaic of Ecuadorian Y-chromosome ancestry

Ecuadorians originated from a complex mixture of Native American indigenous people with Europeans and Africans. We analyzed Y-chromosome STRs (Y-STRs) in a sample of 415 Ecuadorians (145 using the AmpFlSTR^® Yfiler™ system [Life Technologies, USA] and 270 using the PowerPlex^®Y23 system [Promega Corp., USA]; hereafter Yfiler and PPY23, respectively) representing three main ecological continental regions of the country, namely Amazon rainforest, Andes, and Pacific coast. Diversity values are high in the three regions, and the PPY23 exhibits higher discrimination power than the Yfiler set. While summary statistics, AMOVA, and R_ST distances show low to moderate levels of population stratification, inferred ancestry derived from Y-STRs reveal clear patterns of geographic variation. The major ancestry in Ecuadorian males is European (61%), followed by an important Native American component (34%); whereas the African ancestry (5%) is mainly concentrated in the Northwest corner of the country. We conclude that classical procedures for measuring population stratification do not have the desirable sensitivity. Statistical inference of ancestry from Y-STRS is a satisfactory alternative for revealing patterns of spatial variation that would pass unnoticed when using popular statistical summary indices.

Phylogeographic and genome-wide investigations of Vietnam ethnic groups reveal signatures of complex historical demographic movements

The territory of present-day Vietnam was the cradle of one of the world’s earliest civilizations, and one of the first world regions to develop agriculture. We analyzed the mitochondrial DNA (mtDNA) complete control region of six ethnic groups and the mitogenomes from Vietnamese in The 1000 Genomes Project (1000G). Genome-wide data from 1000G (~55k SNPs) were also investigated to explore different demographic scenarios. All Vietnamese carry South East Asian (SEA) haplotypes, which show a moderate geographic and ethnic stratification, with the Mong constituting the most distinctive group. Two new mtDNA clades (M7b1a1f1 and F1f1) point to historical gene flow between the Vietnamese and other neighboring countries. Bayesian-based inferences indicate a time-deep and continuous population growth of Vietnamese, although with some exceptions. The dramatic population decrease experienced by the Cham 700 years ago (ya) fits well with the Nam tiến (“southern expansion”) southwards from their original heartland in the Red River Delta. Autosomal SNPs consistently point to important historical gene flow within mainland SEA, and add support to a main admixture event occurring between Chinese and a southern Asian ancestral composite (mainly represented by the Malay). This admixture event occurred ~800 ya, again coinciding with the Nam tiến.

Meta-Analysis of Mitochondrial DNA Variation in the Iberian Peninsula

The Iberian Peninsula has been the focus of attention of numerous studies dealing with mitochondrial DNA (mtDNA) variation, most of them targeting the control region segment. In the present study we sequenced the control region of 3,024 Spanish individuals from areas where available data were still limited. We also compiled mtDNA haplotypes from the literature involving 4,588 sequences and 28 population groups or small regions. We meta-analyzed all these data in order to shed further light on patterns of geographic variation, taking advantage of the large sample size and geographic coverage, in contrast with the atomized sampling strategy of previous work. The results indicate that the main mtDNA haplogroups show primarily clinal geographic patterns across the Iberian geography, roughly along a North-South axis. Haplogroup HV0 (where haplogroup U is nested) is more prevalent in the Franco Cantabrian region, in good agreement with previous findings that identified this area as a climate refuge during the Last Glacial Maximum (LGM), prior to a subsequent demographic re-expansion towards Central Europe and the Mediterranean. Typical sub-Saharan and North African lineages are slightly more prevalent in South Iberia, although at low frequencies; this pattern has been shaped mainly by the transatlantic slave trade and the Arab invasion of the Iberian Peninsula. The results also indicate that summary statistics that aim to measure molecular variation, or AMOVA, have limited sensitivity to detect population substructure, in contrast to patterns revealed by phylogeographic analysis. Overall, the results suggest that mtDNA variation in Iberia is substantially stratified. These patterns might be relevant in biomedical studies given that stratification is a common cause of false positives in case-control mtDNA association studies, and should be also considered when weighting the DNA evidence in forensic casework, which is strongly dependent on haplotype frequencies.

Revealing latitudinal patterns of mitochondrial DNA diversity in Chileans

The territory of Chile is particularly long and narrow, which combined with its mountainous terrain, makes it a unique scenario for human genetic studies. We obtained 995 control region mitochondrial DNA (mtDNA) sequences from Chileans representing populations living at different latitudes of the country from the North to the southernmost region. The majority of the mtDNA profiles are of Native American origin (∼88%). The remaining haplotypes are mostly of recent European origin (∼11%), and only a minor proportion is of recent African ancestry (∼1%). While these proportions are relatively uniform across the country, more structured patterns of diversity emerge when examining the variation from a phylogeographic perspective. For instance, haplogroup A2 reaches ∼9% in the North, and its frequency decreases gradually to ∼1% in the southernmost populations, while the frequency of haplogroup D (sub-haplogroups D1 and D4) follows the opposite pattern: 36% in the southernmost region, gradually decreasing to 21% in the North. Furthermore, there are remarkable signatures of founder effects in specific sub-clades of Native American (e.g. haplogroups D1j and D4p) and European (e.g. haplogroups T2b3 and K1a4a1a+195) ancestry. We conclude that the magnitude of the latitudinal differences observed in the patterns of mtDNA variation might be relevant in forensic casework.

Northern Slavs from Serbia do not show a founder effect at autosomal and Y-chromosomal STRs and retain their paternal genetic heritage

Studies on Y-chromosomal markers revealed significant genetic differentiation between Southern and Northern (Western and Eastern) Slavic populations. The northern Serbian region of Vojvodina is inhabited by Southern Slavic Serbian majority and, inter alia, Western Slavic (Slovak) and Eastern Slavic (Ruthenian) minorities. In the study, 15 autosomal STR markers were analysed in unrelated Slovaks, Ruthenians and Serbs from northern Serbia and western Slovakia. Additionally, Slovak males from Serbia were genotyped for 17 Y-chromosomal STR loci. The results were compared to data available for other Slavic populations. Genetic distances for autosomal markers revealed homogeneity between Serbs from northern Serbia and Slovaks from western Slovakia and distinctiveness of Serbian Slovaks and Ruthenians. Y-STR variation showed a clear genetic departure of the Slovaks and Ruthenians inhabiting Vojvodina from their Serbian neighbours and genetic similarity to the Northern Slavic populations of Slovakia and Ukraine. Admixture estimates revealed negligible Serbian paternal ancestry in both Northern Slavic minorities of Vojvodina, providing evidence for their genetic isolation from the Serbian majority population. No reduction of genetic diversity at autosomal and Y-chromosomal markers was found, excluding genetic drift as a reason for differences observed at autosomal STRs. Analysis of molecular variance detected significant population stratification of autosomal and Y-chromosomal microsatellites in the three Slavic populations of northern Serbia, indicating necessity for separate databases used for estimations of frequencies of autosomal and Y-chromosomal STR profiles in forensic casework. Our results demonstrate that regarding Y-STR haplotypes, Serbian Slovaks and Ruthenians fit in the Eastern European metapopulation defined in the Y chromosome haplotype reference database.

Ancestry analysis reveals a predominant Native American component with moderate European admixture in Bolivians

We have genotyped 46 Ancestry Informative Markers (AIMs) in two of the most populated areas in Bolivia, namely, La Paz (Andean region; n=105), and Chuquisaca (Sub-Andean region; n=73). Using different analytical tools, we inferred admixture proportions of these two American communities by comparing the genetic profiles with those publicly available from the CEPH (Centre d'Etude du Polymorphisme Humain) panel representing three main continental groups (Africa, Europe, and America). By way of simulations, we first evaluated the minimum sample size needed in order to obtain accurate estimates of ancestry proportions. The results indicated that sample sizes above 30 individuals could be large enough to estimate main continental ancestry proportions using the 46 AIMs panel. With the exception of a few individuals, the results also indicated that Bolivians showed a predominantly Native American ancestry with variable levels of European admixture. The proportions of ancestry were statistically different in La Paz and Chuquisaca: the Native American component was 86% and 77% (Mann-Whitney U-test: un-adjusted P-value=2.1×10(-5)), while the European ancestry was 13% and 21% (Mann-Whitney U-test: un-adjusted P-value=3.6×10(-5)), respectively. The African ancestry in Bolivians captured by the AIMs analyzed in the present study was below 2%. The inferred ancestry of Bolivians fits well with previous studies undertaken on haplotype data, indicating a major proportion of Native American lineages. The genetic differences observed in these two groups suggest that forensic genetic analysis should be better performed based on local databases built in the main Bolivian areas.