Extensive geographical and social structure in the paternal lineages of Saudi Arabia revealed by analysis of 27 Y-STRs

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.

Genetic Analysis of 27 Y-STR Haplotypes in 11 Iranian Ethnic Groups

BACKGROUND

The study of Y-chromosomal variations provides valuable insights into male susceptibility in certain diseases like cardiovascular disease (CVD). In this study, we analyzed paternal lineage in different Iranian ethnic groups, not only to identify developing medical etiology, but also to pave the way for gender-specific targeted strategies and personalized medicine in medical genetic research studies.

METHODS

The diversity of eleven Iranian ethnic groups was studied using 27 Y-chromosomal short tandem repeat (Y-STR) haplotypes from Y-filer® Plus kit. Analysis of molecular variance (AMOVA) based on pair-wise RST along with multidimensional scaling (MDS) calculation and Network phylogenic analysis was employed to quantify the differences between 503 unrelated individuals from each ethnicity.

RESULTS

Results from AMOVA calculation confirmed that Gilaks and Azeris showed the largest genetic distance (RST=0.35434); however, Sistanis and Lurs had the smallest considerable genetic distance (RST=0.00483) compared to other ethnicities. Although Azeris had a considerable distance from other ethnicities, they were still close to Turkmens. MDS analysis of ethnic groups gave the indication of lack of similarity between different ethnicities. Besides, network phylogenic analysis demonstrated insignificant clustering between samples.

CONCLUSION

The AMOVA analysis results explain that the close distance of Azeris and Turkmens may be the effect of male-dominant expansions across Central Asia that contributed to historical and demographics of populations in the region. Insignificant differences in network analysis could be the consequence of high mutation events that happened in the Y-STR regions over the years. Considering the ethnic group affiliations in medical research, our results provided an understanding and characterization of Iranian male population for future medical and population genetics studies.

The Qatari population's genetic structure and gene flow as revealed by the Y chromosome

The Y-chromosome has been widely used in forensic genetic applications and human population genetic studies due to its uniparental origins. A large database on the Qatari population was created for comparison with other databases from the Arabian Peninsula, the Middle East, and Africa. We provide a study of 23 Y-STR loci included in PowerPlex Y23 (Promega, USA) that were genotyped to produce haplotypes in 379 unrelated males from Qatar, a country at the crossroads of migration patterns. Overall, the most polymorphic locus provided by the Promega kit was DYS458, with a genetic diversity value of 0.85 and a haplotype diversity of 0.998924. Athey's Haplogroup Predictor tool was used to predict haplogroups from Y-STR haplotypes in the Qatari population. In a median-joining network, the haplogroup J1 predominance (49%) in Qatar generated a star-like expansion cluster. The graph of population Q-matrix was developed using Y-STR data from 38 Middle Eastern and 97 African populations (11,305 individuals), and it demonstrated a stronger sub-grouping of countries within each ethnic group and showed the effect of Arabs on the indigenous Berbers of North Africa. The estimated migration rate between the Qatari and other Arabian populations was inferred using Bayesian coalescence theory in the Migrate-n program. According to the Gene Flow study, the main migration route was from Yemen to Kuwait through Qatar. Our research, using the PowerPlex Y23 database, shows the importance of gene diversity, as well as regional and social structuring, in determining the utility of demographic and forensic databases.

Allelic and haplotypic polymorphisms and paternal genetic analysis of Chinese Shaanxi Han population utilizing a multiplex Y-STR set

BACKGROUND

The analysis of Y chromosomal genetic markers is of great significance in human genetic fields related to male individuals. The Han nationality is the most populous ethnic group. It is critical to investigate the Y-chromosome short tandem repeat (Y-STR) genetic informativeness of Han nationalities in different Chinese regions in order to gain a comprehensive understanding of their paternal genetic relationships and origin.

AIM

To assess the allelic and haplotypic polymorphisms of the novel AGCU Y SUPP STR amplification system containing seven Y-STRs in the maximal dataset of the Y-STR Haplotype Reference Database (YHRD) and 17 newly included Y-STRs, and explore the genetic relationships among the Shaanxi Han population and 12 reference populations from China.

SUBJECTS AND METHODS

A total sample of 220 Han male subjects were obtained from the Shaanxi Province, China, and genotyped by the novel AGCU Y SUPP STR amplification system. Multiplex population genetic analyses derived from the same 16 Y-STR loci were carried out among the Shaanxi Han population and 12 reference populations from China.

RESULTS

The gene diversities (GD) ranged from the maximum value of 0.9609 (DYS385a,b) to the minimum value of 0.5441 (DYS531). Besides, 217 distinct haplotypes were detected wholly in 220 individuals, of which 214 (98.62%) were exclusive. The entire haplotype diversity (HD) and discrimination capacity (DC) were 0.9999 and 0.9864, respectively, while the haplotype match probability (HMP) was 0.0045. Among the reference populations, the obtained results of population genetic analyses revealed that the Shaanxi Han population had the largest genetic distance with the Guangxi Yao group, but the smallest genetic distance with the Hunan Tujia group.

CONCLUSIONS

These Y-STR loci in the AGCU Y SUPP STR amplification system were of high genetic polymorphisms and the amplification system could be used as a prospective complementary tool for forensic application and paternal genetics in the Shaanxi Han population.

Population genetic study of 17 Y-STR Loci of the Sorani Kurds in the Province of Sulaymaniyah, Iraq

Background

The Kurds as an ethnic group are believed to be a combination of earlier Indo-European tribes who migrated and inhabited a mountainous area thousands of years ago. However, as it is difficult to describe the precise history of their origin, it is necessary to investigate their population relationship with other geographical and ethnic groups.

Results

Seventeen Short Tandem Repeat markers on the Y chromosome (Y-STR) included in the AmpFLSTR™ Yfiler™ PCR Amplification Kit (Thermo Fisher Scientific, USA) were used to type DNA samples from the Sorani (Central) Kurdish population in Sulaymaniyah province. One hundred fifty-seven haplotypes were obtained from 162 unrelated male individuals. The highest and lowest gene diversities were DYS385a/b (GD = 0.848) and DYS392 (GD = 0.392), respectively. The haplotypes were used to predict the most likely haplogroups in the Sulaymaniyah population.

Conclusion

Haplogroup prediction indicated predominance (28%) of subclade J2 (44/157) in the Sorani Kurds, northeast of Iraq. The pairwise genetic distance results showed that the Kurdish group clustered along with Asian populations, whereas the furthest countries were Europeans and Africans.

The Impact of Leukemia on the Detection of Short Tandem Repeat (STR) Markers

INTRODUCTION

Short tandem repeats (STRs) have been used for various identity typing methods worldwide. They have high discrimination power in human identification in forensics, paternity testing, missed personal identification, genetic diseases, and gene regulatory functions. They have also been used to detect and monitor the stability of diseases, including various types of cancer. This study aimed to investigate the impact of leukemia on the detection and stability of STR markers.

METHODS

DNA was isolated from 30 participants (15 with chronic myeloid leukemia( CML) and 15 healthy controls) and used to amplify STR markers using specific primers.

RESULTS

We found that the blood of those with leukemia had more 9.3 and 9 alleles at the tyrosine hydroxylase 1 (TH01) marker than the blood of the healthy control samples. The results of this study will help researchers understand leukemia's effect on the detection and stability of STR markers in leukemic patients compared to healthy individuals.

CONCLUSION

Our results demonstrate that STR markers could become useful in genetic studies of leukemia cases.

Influence of Adiposity on the Gut Microbiota Composition of Arab Women: A Case-Control Study

Recent evidence has suggested that the gut microbiota is a possible risk factor for obesity. However, limited evidence is available on the association between the gut microbiota composition and obesity markers in the Middle-Eastern region. We aimed to investigate the association between gut microbiota and obesity markers in a case-control study including 92 Saudi women aged 18-25 years, including participants with obesity (case, n = 44) and with normal weight (control, n = 48). Anthropometric, body composition, and biochemical data were collected. The whole-genome shotgun technique was used to analyze the gut microbiota. The Shannon alpha and Bray-Curtis beta diversity were determined. The microbial alpha diversity was significantly associated with only the waist-to-hip ratio (WHR) (p-value = 0.04), while the microbial beta diversity was significantly associated with body mass index (p-value = 0.048), %body fat (p-value = 0.018), and WHR (p-value = 0.050). Specific bacteria at different taxonomic levels, such as Bacteroidetes and Synergistetes, were positively associated with different obesity markers. Alistipes was higher in the control group compared with the case group. The results highlight the association of the gut microbiota with obesity and suggest that the gut microbiota of Saudi women is associated with specific obesity markers. Future studies are needed to determine the role of the identified strains in the metabolism of individuals with obesity.

Forensic genetic value of 27 Y-STR loci (Y-Filer® Plus) in the South African population

South Africa has one of the highest rape statistics in the world, with an average of 117 rapes reported daily. Y-STR genotyping is becoming a popular tool in the analysis of DNA evidence collected after a crime of a sexual nature has been committed, but has yet to be implemented in South Africa's forensic laboratories. This study aimed to investigate the forensic value of the 27 Yfiler™ Plus loci in the South African population. A total of 271 samples from the African, Asian/Indian, Mixed Ancestry¹, and Caucasian populations at the University of the Free State in Bloemfontein, South Africa were amplified and analysed using ThermoFisher Scientific's Yfiler™ Plus PCR Amplification kit. Of the 271 samples, 261 were identified to be unique, with an overall discrimination capacity of 98.15%. Discrimination capacities ranged from 91.67% for the Asian/Indian population to 100% for the Mixed Ancestry population. The haplotype diversity across the four populations is 0.9999, with an average gene diversity across all loci of 0.717. The forensic parameters estimated in this study provide evidence for the potential use of the commercial Yfiler™ Plus PCR amplification kit in a forensic application in South Africa.

Genetic Relationship Among the Kazakh People Based on Y-STR Markers Reveals Evidence of Genetic Variation Among Tribes and Zhuz

Ethnogenesis of Kazakhs took place in Central Asia, a region of high genetic and cultural diversity. Even though archaeological and historical studies have shed some light on the formation of modern Kazakhs, the process of establishment of hierarchical socioeconomic structure in the Steppe remains contentious. In this study, we analyzed haplotype variation at 15 Y-chromosomal short-tandem-repeats obtained from 1171 individuals from 24 tribes representing the three socio-territorial subdivisions (Senior, Middle and Junior zhuz) in Kazakhstan to comprehensively characterize the patrilineal genetic architecture of the Kazakh Steppe. In total, 577 distinct haplotypes were identified belonging to one of 20 haplogroups; 16 predominant haplogroups were confirmed by SNP-genotyping. The haplogroup distribution was skewed towards C2-M217, present in all tribes at a global frequency of 51.9%. Despite signatures of spatial differences in haplotype frequencies, a Mantel test failed to detect a statistically significant correlation between genetic and geographic distance between individuals. An analysis of molecular variance found that ∼8.9% of the genetic variance among individuals was attributable to differences among zhuzes and ∼20% to differences among tribes within zhuzes. The STRUCTURE analysis of the 1164 individuals indicated the presence of 20 ancestral groups and a complex three-subclade organization of the C2-M217 haplogroup in Kazakhs, a result supported by the multidimensional scaling analysis. Additionally, while the majority of the haplotypes and tribes overlapped, a distinct cluster of the O2 haplogroup, mostly of the Naiman tribe, was observed. Thus, firstly, our analysis indicated that the majority of Kazakh tribes share deep heterogeneous patrilineal ancestries, while a smaller fraction of them are descendants of a founder paternal ancestor. Secondly, we observed a high frequency of the C2-M217 haplogroups along the southern border of Kazakhstan, broadly corresponding to both the path of the Mongolian invasion and the ancient Silk Road. Interestingly, we detected three subclades of the C2-M217 haplogroup that broadly exhibits zhuz-specific clustering. Further study of Kazakh haplotypes variation within a Central Asian context is required to untwist this complex process of ethnogenesis.

A reference database of forensic autosomal and gonosomal STR markers in the Tigray population of Ethiopia

Allele frequencies of 21 autosomal STR markers (AmpF/STR GlobalFiler) and haplotype frequencies of 27 Y- and 12 X-STR markers (AmpF/STR YFiler Plus and Investigator Argus X-12, respectively) were investigated in the Tigray population of Ethiopia, representing the main population group in the Tigray regional state of Ethiopia and neighboring Eritrea. For autosomal STR allele frequencies, the average random match probability in the Tigray sample was 2.1 × 10^-27. The average locus by locus F_ST distance calculated comparing autosomal STR allele frequencies from Tigray and from a broad regional reference dataset currently available for the Horn of Africa was 0.003. The Tigray male sample displayed high Y-STR diversity, with complete individualization of haplotypes using the AmpF/STR YFiler Plus panel. Analysis of molecular variance did not detect significant heterogeneity between Y-STR haplotypes observed in the present study and those previously reported in the literature for other Tigray population samples from Ethiopia and Eritrea. Study of the X-STR landscape in Tigray evidenced several distinctive features including: the molecular characterization of a novel null allele at locus DXS10146 with frequency > 1%; allele dependency between loci within linkage groups I and III; significant differences in haplotype distribution compared to other Horn of Africa populations, that should be taken into account in kinship analysis. The collected data can be used as a reference STR database by local forensic genetics services and in genetic identification procedures of victims of human trafficking in the Mediterranean Sea, which frequently involve individuals originating from the Horn of Africa.

Genetic characterisation for Yan'an Han population in Northern Shaanxi Province, China, via 38Y-STRs using Yfiler™ Platinum

BACKGROUND

Haplotype/allele frequency data of Y-chromosomal STR loci in ethnically diverse populations are essential for forensics, anthropology and genealogy. However, genetic structure and forensic characterisation of the Chinese Han population residing in Yan'an, in the Northern Shaanxi Province, remain unclear.

AIM

To assess forensic efficiency for 38 Y-Chromosomal STR loci in Yan'an Han population and reveal the population genetic relationships between Yan'an Han and other populations at a nationwide and worldwide level.

SUBJECTS AND METHODS

719 healthy unrelated males were genotyped using the Yfiler™ Platinum system. Haplotype/allele frequencies and forensic parameters were calculated. Multi-dimensional scaling plots (MDS) and neighbor-joining (N-J) tree were used to explore the population structure based on the pairwise gene distances (Rst).

RESULTS

A total of 707 haplotypes were identified, among which 697 unique haplotypes were observed (98.59%). The overall haplotype diversity (HD) and discrimination capacity (DC) were 0.9999 and 0.9833, respectively. Comprehensive population comparisons showed Yan'an Han is genetically closer to linguistically similar populations in China, and more related to East Asian populations around the world.

CONCLUSION

The present results give a unique insight into the Yan'an Han population via the set of 38 Y-STRs, which can be used for forensic practice and human genetics research.

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.

Population data and genetic diversity analysis of 17 Y-STR loci in Saudi population

Background

The Y chromosome polymorphism has been widely studied for human migrations, population genetics, forensic applications, and paternity analysis. However, studies regarding genetic lineage and population genetic structure of the Y chromosome in different regions of Saudi Arabia are limited.

Aim

This study aimed to analyze the distribution of Y chromosome haplotypes in a sample of 125 native Saudi males from different geographic regions of Saudi Arabia and compare to previously published Y chromosome haplotype data from Saudi Arabia and some neighboring Arab populations.

Materials and methods

Buccal swabs were collected from 125 healthy unrelated native Saudi males from different geographic regions of Saudi Arabia. Genomic DNA was extracted by Chelex®100; 17 Y-STR loci were amplified using the AmpFℓlSTR Yfiler PCR amplification kit and detected on the 3130 Genetic AnalyzerTM. Allele frequency and gene diversity were calculated with online tool STRAF. The Saudi population data were compared with the neighboring populations using pairwise genetic distances and associated probability values were calculated using the Y Chromosome Haplotype Reference Database Website (YHRD) software.

Results and conclusion

One hundred six YSTR haplotypes and 102 YSTR alleles (excluding 4 null alleles) were identified having a discrimination capacity (DC) of 85.8%. The highest haplotype diversity (HD) and gene diversity (GD) were observed at the loci DYS 458 (0.817) and DYS385b (0.807), respectively. According to our results, the Iraqi and Qena (Egypt) populations appeared to have closer relatedness to the Saudi population as compared with Yemen. The UAE and Kuwait populations showed the same degree of relatedness to the Saudi population followed by Bahrain. On the contrary, the Adnanit and Qahtanit populations of Jordan demonstrated low genetic distance from the Saudi population. In short, studying a population sample of pure Saudi ethnicity enabled us to identify a unique set of haplotypes which may help in establishing genetic relatedness between Saudi and the neighboring Arab populations. The present paper, therefore, highlights the importance of ensuring ethnic originality of the study sample while conducting population genetics studies.

The Tooth Size Discrepancy among Orthodontic Patients and Normal Occlusion Individuals from Saudi Arabia: A Three-Dimensional Scan Analysis of Diagnostic Casts

Background/Purpose:

Tooth size discrepancy is one of the causative factors of malocclusion. This study aimed to establish the TDS among various malocclusion classes and normal occlusion subjects in a Southern Saudi population.

Materials and Methods:

The study casts of orthodontic patients from Southern Saudi Arabia (mean age: 19.6 years ± 2.45; n = 120) were randomly selected and grouped into eight equal classes based on Angle's classification and gender. An additional 30 study casts, devoid of malocclusion, having excellent intercuspation and class I occlusion, were grouped into two controls (n = 15) based on gender. The study casts were three-dimensional scanned to measure mesiodistal widths of all the teeth. The calculated anterior ratios (AR) and overall ratio (OR) were statistically analyzed with analysis of variance and t-tests.

Results:

There were no significant differences in “OR” and “AR” between the genders (P > 0.05) and among the malocclusion and control subgroups (P > 0.05). The mean “OR” (92.01 ± 0.18) and “AR” (78.60 ± 0.27) of the malocclusion group were significantly higher than that of Bolton's ratios (P > 0.05). The “AR” of the control group was significantly higher than Bolton's standards (P = 0.048). However, “OR” was no different (P = 0.105). Malocclusion patients displayed a discrepancy (±2 standard deviation) in “AR” of 22.5% and 'OR' of 6.7% from Bolton's mean (BM). Similarly, the control group displayed a discrepancy in “AR” of 20% and “OR” of 10% from BM.

Conclusion:

The mean “OR” and “AR” of the Southern Saudi population showed no sexual dimorphism and no significant difference among various malocclusion and control subgroups. The “AR” of the malocclusion and control subgroups did not comply with Bolton's standards.

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.

Population genetic diversity in an Iraqi population and gene flow across the Arabian Peninsula

Y-STRs have emerged as important forensic and population genetic markers for human identification and population differentiation studies. Therefore, population databases for these markers have been developed for almost all major populations around the world. The Iraqi population encompasses several ethnic groups that need to be genetically characterised and evaluated for possible substructures. Previous studies on the Iraqi population based on Y-STR markers were limited by a restricted number of markers. A larger database for Iraqi Arab population needed to be developed to help study and compare the population with other Middle Eastern populations. Twenty-three Y-STR loci included in the PowerPlex Y23 (Promega, Madison, WI, USA) were typed in 254 males from the Iraqi Arab population. Global and regional Y-STR analysis demonstrated regional genetic continuity among the populations of Iraq, the Arabian Peninsula and the Middle East. The Iraqi Arab haplotypes were used to allocate samples to their most likely haplogroups using Athey’s Haplogroup Predictor tool. Prediction indicated predominance (36.6%) of haplogroup J1 in Iraqi Arabs. The migration rate between other populations and the Iraqis was inferred using coalescence theory in the Migrate-n program. Y-STR data were used to test different out-of-Africa migration models as well as more recent migrations within the Arabian Peninsula. The migration models demonstrated that gene flow to Iraq began from East Africa, with the Levantine corridor the most probable passageway out of Africa. The data presented here will enrich our understanding of genetic diversity in the region and introduce a PowerPlex Y23 database to the forensic community.

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.

Genetic structure and forensic characterisation of 36 Y-chromosomal STR loci in Hmong-Mien-speaking Miao population

BACKGROUND

Y-chromosomal short tandem repeats (Y-STRs) are widely used in paternity identification, pedigree investigation, and human population genetic history.

AIM

To investigate the Y-STR polymorphisms in a typical Miao population, and explore the genetic differentiation between the Miao population and reference groups.

SUBJECTS AND METHODS

We detected 36 Y-STRs genotyping in 455 unrelated Miao individuals from Guizhou province, and analysed genetic differentiation between the Miao population and 76 reference groups.

RESULTS

A total of 369 alleles were obtained, and the allele frequencies ranged from 0.0022 to 0.9802. In addition, the haplotype diversity, random match probability, and discrimination capacity values were 0.99997, 0.0022, and 0.9934, respectively. Moreover, the genetic relationships between Guizhou Miao and 76 ethnic populations showed that the population stratification was almost consistent with geographic distribution and language-family.

CONCLUSIONS

The 36 Y-STR loci in this study have good polymorphism distributions in the Guizhou Miao population, and therefore would be a useful tool in forensic identification and male parentage testing and even pedigree investigation.

Genetic Reconstruction and Forensic Analysis of Chinese Shandong and Yunnan Han Populations by Co-Analyzing Y Chromosomal STRs and SNPs

Y chromosomal short tandem repeats (Y-STRs) have been widely harnessed for forensic applications, such as pedigree source searching from public security databases and male identification from male–female mixed samples. For various populations, databases composed of Y-STR haplotypes have been built to provide investigating leads for solving difficult or cold cases. Recently, the supplementary application of Y chromosomal haplogroup-determining single-nucleotide polymorphisms (SNPs) for forensic purposes was under heated debate. This study provides Y-STR haplotypes for 27 markers typed by the Yfiler^™ Plus kit and Y-SNP haplogroups defined by 24 loci within the Y-SNP Pedigree Tagging System for Shandong Han (n = 305) and Yunnan Han (n = 565) populations. The genetic backgrounds of these two populations were explicitly characterized by the analysis of molecular variance (AMOVA) and multi-dimensional scaling (MDS) plots based on 27 Y-STRs. Then, population comparisons were conducted by observing Y-SNP allelic frequencies and Y-SNP haplogroups distribution, estimating forensic parameters, and depicting distribution spectrums of Y-STR alleles in sub-haplogroups. The Y-STR variants, including null alleles, intermedia alleles, and copy number variations (CNVs), were co-listed, and a strong correlation between Y-STR allele variants (“DYS518~.2” alleles) and the Y-SNP haplogroup QR-M45 was observed. A network was reconstructed to illustrate the evolutionary pathway and to figure out the ancestral mutation event. Also, a phylogenetic tree on the individual level was constructed to observe the relevance of the Y-STR haplotypes to the Y-SNP haplogroups. This study provides the evidence that basic genetic backgrounds, which were revealed by both Y-STR and Y-SNP loci, would be useful for uncovering detailed population differences and, more importantly, demonstrates the contributing role of Y-SNPs in population differentiation and male pedigree discrimination.

Geographical structuring and low diversity of paternal lineages in Bahrain shown by analysis of 27 Y-STRs

We have determined the distribution of Y-chromosomal haplotypes and predicted haplogroups in the ethnically diverse Kingdom of Bahrain, a small archipelago in the Arabian Gulf. Paternal population structure within Bahrain was investigated using the 27 Y-STRs (short tandem repeats) in the Yfiler Plus kit to generate haplotypes from 562 unrelated Bahraini males, sub-divided into four geographical regions—Northern, Capital, Southern and Muharraq. Yfiler Plus provided a significant improvement over the 17-locus Yfiler kit in discrimination capacity (from 77% to 87.5% overall), but discrimination capacity differed widely between regions from 98.4% in Muharraq to 75.2% in the Northern region, an unusually low value possibly resulting from recent rapid population expansion. Clusters of closely related male lineages were seen, with only 79.4% of donors displaying unique haplotypes and 59% of instances of shared haplotypes occurring within, rather than between, regions. Haplogroup prediction indicated diverse origins of the population with a predominance of haplogroups J2 and J1, both typical of the Arabian Peninsula, but also haplogroups such as B2 and E1b1a likely originating in Africa, and H, L and R2 likely indicative of migration from South Asia. Haplogroup frequencies differed significantly between regions, with J2 significantly more common in the Northern region compared with the Southern, possibly due to differential settlement by Baharna and Arabs. Our study shows that paternal lineage population structure can exist even over small geographical scales, and that highly discriminating genetic tools are required where rapid expansions have occurred within tightly bounded populations.

Genome-Wide Selection Scan in an Arabian Peninsula Population Identifies a TNKS Haplotype Linked to Metabolic Traits and Hypertension

Despite the extreme and varying environmental conditions prevalent in the Arabian Peninsula, it has experienced several waves of human migrations following the out-of-Africa diaspora. Eventually, the inhabitants of the peninsula region adapted to the hot and dry environment. The adaptation and natural selection that shaped the extant human populations of the Arabian Peninsula region have been scarcely studied. In an attempt to explore natural selection in the region, we analyzed 662,750 variants in 583 Kuwaiti individuals. We searched for regions in the genome that display signatures of positive selection in the Kuwaiti population using an integrative approach in a conservative manner. We highlight a haplotype overlapping TNKS that showed strong signals of positive selection based on the results of the multiple selection tests conducted (integrated Haplotype Score, Cross Population Extended Haplotype Homozygosity, Population Branch Statistics, and log-likelihood ratio scores). Notably, the TNKS haplotype under selection potentially conferred a fitness advantage to the Kuwaiti ancestors for surviving in the harsh environment while posing a major health risk to present-day Kuwaitis.

Massively parallel sequencing of autosomal STRs and identity-informative SNPs highlights consanguinity in Saudi Arabia

While many studies have been undertaken of Middle Eastern populations using autosomal STR profiling by capillary electrophoresis, little has so far been published from this region on the forensic use of massively parallel sequencing (MPS). Here, we carried out MPS of 27 autosomal STRs and 91 identity-informative SNPs (iiSNPs) with the Verogen ForenSeq™ DNA Signature Prep Kit on a representative sample of 89 Saudi Arabian males, and analysed the resulting sequence data using Verogen's ForenSeq Universal Analysis Software (UAS) v1.3 and STRait Razor v3.0. This revealed sequence variation in the composition of complex STR arrays, and SNPs in their flanking regions, which raised the number of STR alleles from 238 distinct length variants to 357 sequence sub-variants. Similarly, between one and three additional polymorphic sites were observed within the amplicons of 37 of the 91 iiSNPs, forming up to six microhaplotypes per locus. These further enhance discrimination compared to the biallelic target SNP data presented by the primary UAS interface. In total, we observed twenty-two STR alleles previously unrecognised in the STRait Razor v3.0 default allele list, along with nine SNPs flanking target iiSNPs that were not highlighted by the UAS. Sequencing reduced the STR-based random match probability (RMP) from 2.62E-30 to 3.49E-34, and analysis of the iiSNP microhaplotypes reduced RMP from 9.97E-37 to 6.83E-40. The lack of significant linkage disequilibrium between STRs and target iiSNPs allowed the two marker types to be combined using the product rule, yielding a RMP of 2.39E-73. Evidence of consanguinity was apparent from both marker types. While TPOX was the only locus displaying a significant deviation from Hardy-Weinberg equilibrium, 23 out of 27 STRs and 63 out of 91 iiSNPs showed fewer than expected heterozygotes, demonstrating an overall homozygote excess probably reflecting the high frequency of first-cousin marriages in Saudi Arabia. We placed our data in a global context by considering the same markers in the Human Genome Diversity Panel (HGDP), revealing that the Saudi sample was typical of Middle Eastern populations, with a higher level of inbreeding than is seen in most European, African and Central/South Asian populations, correlating with known patterns of endogamy. Given reduced levels of diversity within endogamous groups, the ability to combine the discrimination power of both STRs and SNPs offers significant benefits in the analysis of forensic evidence in Saudi Arabia and the Middle East region more generally.

Demand for larger Y-STR reference databases in ethnic melting-pot countries: Argentina as a test case

The Y chromosome behaves as a single locus. Its genetic information is useful in forensic casework, deficiency kinship testing, and population genetics studies. Continuous increases of loci number within commercial kits forced modification of worldwide reference databases. In Pan American countries, like Argentina, diverse parental ethnic groups contributed to the extant admixed urban populations. We report 509 additional haplotypes of 23 Y-STRs from donors inhabiting urban areas of six Argentinean provinces: Buenos Aires, Santiago del Estero, Santa Cruz, Rio Negro, Santa Fe, and Formosa. To better understand the demographic landscape of the admixed urban paternal lineages, structural analysis was performed using published data from other Argentinean provinces. AMOVA by Rst distance and inferred haplogroups by two predictive online software methods based on haplotypes yielded complementary results with respect to detected population structure, probably due to the different proportions of the Native American Q3-M3 haplogroup in the studied samples. This situation, which is common to most North, Meso, and South American countries, underscores the need for the additional step of typing specific SNPs for haplogroup diagnosis. We propose organizing Y-STR haplotype reference databases according to the most frequent haplogroups detected in a given admixed population.

Recombination hotspots in an extended human pseudoautosomal domain predicted from double-strand break maps and characterized by sperm-based crossover analysis

The human X and Y chromosomes are heteromorphic but share a region of homology at the tips of their short arms, pseudoautosomal region 1 (PAR1), that supports obligate crossover in male meiosis. Although the boundary between pseudoautosomal and sex-specific DNA has traditionally been regarded as conserved among primates, it was recently discovered that the boundary position varies among human males, due to a translocation of ~110 kb from the X to the Y chromosome that creates an extended PAR1 (ePAR). This event has occurred at least twice in human evolution. So far, only limited evidence has been presented to suggest this extension is recombinationally active. Here, we sought direct proof by examining thousands of gametes from each of two ePAR-carrying men, for two subregions chosen on the basis of previously published male X-chromosomal meiotic double-strand break (DSB) maps. Crossover activity comparable to that seen at autosomal hotspots was observed between the X and the ePAR borne on the Y chromosome both at a distal and a proximal site within the 110-kb extension. Other hallmarks of classic recombination hotspots included evidence of transmission distortion and GC-biased gene conversion. We observed good correspondence between the male DSB clusters and historical recombination activity of this region in the X chromosomes of females, as ascertained from linkage disequilibrium analysis; this suggests that this region is similarly primed for crossover in both male and female germlines, although sex-specific differences may also exist. Extensive resequencing and inference of ePAR haplotypes, placed in the framework of the Y phylogeny as ascertained by both Y microsatellites and single nucleotide polymorphisms, allowed us to estimate a minimum rate of crossover over the entire ePAR region of 6-fold greater than genome average, comparable with pedigree estimates of PAR1 activity generally. We conclude ePAR very likely contributes to the critical crossover function of PAR1.

95% of our genome is contained in 22 pairs of chromosomes shared by all humans. However, women and men differ in their sex chromosomes: while women have two X chromosomes, men have an X and a smaller, sex-determining Y chromosome. To ensure correct partition of X and Y into sperm, genetic exchange (crossover) must occur between these very different chromosomes in a short, shared region. The location of the boundary of this region was thought to have been conserved since before the divergence from old world monkeys at least 27 million years ago, but recently it has been shown that some human males carry an extended version on their Y chromosomes, thanks to the transposition of a piece of DNA from the X chromosome. Here, we asked if genetic exchange occurs in this newly extended region. To do this, we used previously published information that signposted the positions within the X chromosome segment which exhibit the hallmarks of crossover initiation. We then sought direct evidence of crossover in the sperm of men carrying the extension. This work showed that the signposts were accurate, pointing to frequent crossover in this novel shared sex-chromosomal domain.

A phylogenetic framework facilitates Y-STR variant discovery and classification via massively parallel sequencing

•

23 Y-chromosomal STRs (PPY23) reanalysed by massively parallel sequencing.

•

Phylogeny-based approach captures wide range of sequence variants in 100 samples.

•

STR variants described in phase with their flanking sequences.

•

Phylogenetic framework clarifies allele nomenclature and mutation processes.

Short tandem repeats on the male-specific region of the Y chromosome (Y-STRs) are permanently linked as haplotypes, and therefore Y-STR sequence diversity can be considered within the robust framework of a phylogeny of haplogroups defined by single nucleotide polymorphisms (SNPs). Here we use massively parallel sequencing (MPS) to analyse the 23 Y-STRs in Promega’s prototype PowerSeq™ Auto/Mito/Y System kit (containing the markers of the PowerPlex® Y23 [PPY23] System) in a set of 100 diverse Y chromosomes whose phylogenetic relationships are known from previous megabase-scale resequencing. Including allele duplications and alleles resulting from likely somatic mutation, we characterised 2311 alleles, demonstrating 99.83% concordance with capillary electrophoresis (CE) data on the same sample set. The set contains 267 distinct sequence-based alleles (an increase of 58% compared to the 169 detectable by CE), including 60 novel Y-STR variants phased with their flanking sequences which have not been reported previously to our knowledge. Variation includes 46 distinct alleles containing non-reference variants of SNPs/indels in both repeat and flanking regions, and 145 distinct alleles containing repeat pattern variants (RPV). For DYS385a,b, DYS481 and DYS390 we observed repeat count variation in short flanking segments previously considered invariable, and suggest new MPS-based structural designations based on these. We considered the observed variation in the context of the Y phylogeny: several specific haplogroup associations were observed for SNPs and indels, reflecting the low mutation rates of such variant types; however, RPVs showed less phylogenetic coherence and more recurrence, reflecting their relatively high mutation rates. In conclusion, our study reveals considerable additional diversity at the Y-STRs of the PPY23 set via MPS analysis, demonstrates high concordance with CE data, facilitates nomenclature standardisation, and places Y-STR sequence variants in their phylogenetic context.