Male lineages in Brazilian populations and performance of haplogroup prediction tools

Detection of a novel 16.3 variant allele at locus DYS533 in R1b males inhabiting southern South America: A 19-nucleotide insertion explains its origin based on Sanger sequencing results

We typed 1541 Y-STR haplotypes from reference samples along forensic casework investigations. In three haplotypes, we detected a variant allele designed as 16.3 at locus DYS533. This was confirmed by amplification using two commercial kits. Sanger sequencing revealing a novel motif corresponding to [TATC]₁₂ repeats with a 19-bp insertion in the flanking upstream region. We propose its origin as an insertion at - 9.1 upstream of the repeat motifs. We searched other local databases and found this allele in various geographical areas of Argentina and neighbouring countries. The haplotypes share a common core of 10 Y-STRs (DYS389-I/13; DYS389-II/30; DYS19/14; DYS481/22; DYS438/12; DYS437/16; DYS635/23; DYS392/13; DYS393/13; GATA H4/11) and belong to the R1b haplogroup. This 16.3 allele is restricted to southern South America, which allows us to propose a local and relatively recent origin. The sequence described herein constitutes a novelty that could be considered in future criteria for the nomenclature of STRs based on massively parallel sequencing.

A forensic population database in El Salvador: 58 STRs and 94 SNPs

We have genotyped the 58 STRs (27 autosomal, 24 Y-STRs and 7 X-STRs) and 94 autosomal SNPs in Illumina ForenSeq™ Primer Mix A in a sample of 248 men and 143 women from El Salvador, Central America. Regional division (Centro, Oriente, Occidente) showed in almost all cases F_ST values not significantly different from 0, and further analyses were applied only to the undivided, country-wide population. The overall random match probability (RMP) decreased from 6.79 × 10^-31 in length-based genotypes in the 27 autosomal STRs to 1.47 × 10^-34 in repeat-sequence based genotypes. Combining the autosomal loci in this set, RMP reaches 2.97 × 10^-70. In a population genetic analysis, El Salvador showed the lowest F_ST values with US Hispanics both for autosomal and X-STRs; however, it was much closer to Native Americans for the latter than for the former, in accordance with the well-known gender-biased admixture that created most Latin American populations.

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.

The Y-STR landscape of coastal southeastern Han: Forensic characteristics, haplotype analyses, mutation rates, and population genetics

The Y-STR landscape of Coastal Southeastern Han (CSEH) living in Chinese southeast areas (including Guangdong, Fujian, and Zhejiang provinces) is still unclear. We investigated 62 Y-STR markers in a reasonably large number of 1021 unrelated males and 1027 DNA-confirmed father-son pairs to broaden the genetic backgrounds of CSEH. In total, 85 null alleles, 121 off-ladder alleles, and 95 copy number variants were observed, and 1012 distinct haplotypes were determined with the overall HD and DC values of 0.999974 and 0.9912. We observed 369 mutations in 76 099 meiotic transfers, and the average estimated Y-STR mutation rate was 4.85 × 10^-3 (95% CI, 4.4 × 10^-3 -5.4 × 10^-3 ). The Spearman correlation analyses indicated that GD values (R² = 0.6548) and average allele sizes (R² = 0.5989) have positive correlations with Y-STR mutation rates. Our RM Y-STR set including 8 candidate RM Y-STRs, of which DYS534, DYS630, and DYS713 are new candidates in CSEH, distinguished 18.52% of father-son pairs. This study also clarified the population structures of CSEH which isolated in population-mixed South China relatively. The strategy, SM Y-STRs for familial searching and RM Y-STRs for individual identification regionally, could be applicable based on enough knowledge of the Y-STR mutability of different populations.

Insights From Y-STRs: Forensic Characteristics, Genetic Affinities, and Linguistic Classifications of Guangdong Hakka and She Groups

Guangdong province is situated in the south of China with a population size of 113.46 million. Hakka is officially recognized as a branch of Han Chinese, and She is the official minority group in mainland China. There are approximately 25 million Hakka people who mainly live in the East and North regions of China, while there are only 0.7 million She people. The genetic characterization and forensic parameters of these two groups are poorly defined (She) or still need to be explored (Hakka). In this study, we have genotyped 475 unrelated Guangdong males (260 Hakka and 215 She) with Promega PowerPlex^® Y23 System. A total of 176 and 155 different alleles were observed across all 23 Y-STRs for Guangdong Hakka (with a range of allele frequencies from 0.0038 to 0.7423) and Guangdong She (0.0047–0.8605), respectively. The gene diversity ranged from 0.4877 to 0.9671 (Guangdong Hakka) and 0.3277–0.9526 (Guangdong She), while the haplotype diversities were 0.9994 and 0.9939 for Guangdong Hakka and Guangdong She, with discrimination capacity values of 0.8885 and 0.5674, respectively. With reference to geographical and linguistic scales, the phylogenetic analyses showed us that Guangdong Hakka has a close relationship with Southern Han, and the genetic pool of Guangdong Hakka was influenced by surrounding Han populations. The predominant haplogroups of the Guangdong She group were O2-M122 and O2a2a1a2-M7, while Guangdong She clustered with other Tibeto-Burman language-speaking populations (Guizhou Tujia and Hunan Tujia), which shows us that the Guangdong She group is one of the branches of Tibeto-Burman populations and the Huonie dialect of She languages may be a branch of Tibeto-Burman language families.

The Y chromosome ancestry marker R1b1b2: a surrogate of the SARS-CoV-2 population affinity

Individual and population susceptibilities to disease remain a murky area of investigation, clouded by past bias based on ideological differences and wars. The current SARS-CoV-2 pandemic, the largest in living memory, brought this matter to forefront as the disparity in disease burden became apparent. A timeline analysis of the pandemic revealed the presence of country clusters that display a marked preponderance of disease among populations carrying the ancestry marker R1b1b2, notably associated with both infection and mortality. This marker is a relic of past human expansions from western Asia and subsequently Europe and the rest of the world, which may have been accompanied by peculiar biological events rendering these populations vulnerable to SARS-CoV-2.

Predicting haplogroups using a versatile machine learning program (PredYMaLe) on a new mutationally balanced 32 Y-STR multiplex (CombYplex): Unlocking the full potential of the human STR mutation rate spectrum to estimate forensic parameters

We developed a new mutationally well-balanced 32 Y-STR multiplex (CombYplex) together with a machine learning (ML) program PredYMaLe to assess the impact of STR mutability on haplogourp prediction, while respecting forensic community criteria (high DC/HD). We designed CombYplex around two sub-panels M1 and M2 characterized by average and high-mutation STR panels. Using these two sub-panels, we tested how our program PredYmale reacts to mutability when considering basal branches and, moving down, terminal branches. We tested first the discrimination capacity of CombYplex on 996 human samples using various forensic and statistical parameters and showed that its resolution is sufficient to separate haplogroup classes. In parallel, PredYMaLe was designed and used to test whether a ML approach can predict haplogroup classes from Y-STR profiles. Applied to our kit, SVM and Random Forest classifiers perform very well (average 97 %), better than Neural Network (average 91 %) and Bayesian methods (< 90 %). We observe heterogeneity in haplogroup assignation accuracy among classes, with most haplogroups having high prediction scores (99-100 %) and two (E1b1b and G) having lower scores (67 %). The small sample sizes of these classes explain the high tendency to misclassify the Y-profiles of these haplogroups; results were measurably improved as soon as more training data were added. We provide evidence that our ML approach is a robust method to accurately predict haplogroups when it is combined with a sufficient number of markers, well-balanced mutation rate Y-STR panels, and large ML training sets. Further research on confounding factors (such as CNV-STR or gene conversion) and ideal STR panels in regard to the branches analysed can be developed to help classifiers further optimize prediction scores.

Mutational data and population profiling of 23 Y-STRs in three Brazilian populations

Y-chromosomal STRs are important markers in forensic genetics, due to some peculiar characteristics. The absence of recombination makes them a useful tool to infer kinship in complex cases involving distant paternal relatives, or to infer paternal bio-geographic ancestry. The presence of a single copy, being transmitted from father to son, allow tracing mutational events in Y-STRs without ambiguity. For the statistical interpretation of forensic evidences based on Y-STR profiles, it is necessary to have estimates on both mutation rates and haplotype frequencies. In this work, 407 father-son duos from São Paulo and Rio de Janeiro states and 204 unrelated individuals from Manaus were analyzed. Haplotype frequencies and mutation rates for the Y-STRs from the PowerPlex Y23 commercial kit were estimated. Thirty-six mutations were observed in 15 of the 22 Y-STRs analyzed, for an average mutation rate of 3.84 × 10^-3 (95 % CI 2.69 × 10^-3 to 5.32 × 10^-3). All mutations in GAAA repeats occurred in alleles with 13 or more uninterrupted units. Mutations in GATA repeats were observed in alleles with 9-17 uninterrupted units. An analysis carried out in different father's age groups showed an increase of 2.48 times the mutation rate in the age group of 40-50 years, when compared to the 20-30 age group, in agreement with the described for autosomal STRs. A high haplotype diversity was found in the three Brazilian populations. Pairwise genetic distance analysis (F_ST) showed no significant differences between the three populations in this study, which were also close to populations with strong European influence. The highest distances among the Brazilian populations were with São Gabriel da Cachoeira, which has a high Native American ancestry.

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.