Opt-in or out? Public perspectives on forensic DNA kinship investigations within the Dutch-speaking community

Forensic DNA kinship investigation involves analyzing genetic relationships between individuals to offer new leads for solving (cold) cases. Familial DNA matching has become a valuable asset in criminal case investigations, especially when traditional DNA methods hit dead ends. However, concerns surrounding ethical and privacy implications raised questions about its implementation and acceptance among the general public. The present study investigated the public perspectives regarding forensic DNA kinship investigations among 1710 Dutch-speaking Belgians using an online cross-sectional survey. The questionnaire consisted of three categories, including personal information, DNA knowledge, and their opinion on several familial DNA searching and investigative genetic genealogy related questions. The participants' average DNA knowledge score was 71 %, indicating a relatively high level of understanding of DNA-related concepts. Remarkably, the study revealed that 92 % of the participants expressed willingness to cooperate as a volunteer in a forensic DNA kinship investigation, irrespective of their scientific background or educational level. Key factors influencing participation included assurance of painless sampling and robust privacy safeguards. Participants lacking familiarity with DNA hesitated more towards participating in forensic DNA analysis, referring to "the fear of the unknown". Despite ethical and privacy concerns, the highly positive attitude towards forensic DNA analysis reflects a level of empathy and willingness to contribute to the pursuit of justice. Nearly all participants (95 %) agreed to use online DNA databases for resolving violent crimes with forensic genetic genealogy, but half emphasized the need for prior informed consent, referring to the current "opt-in" system. The results underscore the need for stringent regulations and ethical oversight to ensure the responsible use of genetic data while striking a balance between public safety and the protection of individuals' privacy rights. These findings add to the growing body of evidence regarding the potential benefits of forensic DNA kinship matching as a tool in criminal investigations, suggesting its potential future utilization and legalization.

Age prediction using DNA methylation of Y-chromosomal CpGs in semen samples

In cases of sexual assault, the evidence often exists as a mixture of female and male body fluids, and in many cases, contains a higher proportion of female body fluids than males. In these cases, Y-STR, rather than autosomal STRs, can provide useful information. It becomes very difficult to identify the true suspect if there is no match among known suspects or if a match exists for two or more suspects, e.g. two suspects from the same paternal lineage. However, age prediction using the DNA methylation of Y-chromosomal CpGs can help narrow the search for unknown suspects and discriminate between older and younger suspects. Therefore, the DNA methylation profiles of semen samples from 56 healthy Korean males were generated using Illumina's Infinium MethylationEPIC BeadChip Array. Among the ten identified age-associated CpG markers located in the Y-chromosome, nine were used to construct age prediction models. The identified markers were further investigated in the MPS analysis of 147 semen samples, and the multiplex assay was validated with the reliability, reproducibility and sensitivity tests. Several age prediction models were constructed using the MPS data with the multiple linear regression, stepwise linear regression, ridge linear regression, lasso regression, elastic net linear regression and support vector machine analyses, and all showed MAEs of 5 to 7 years in the test set samples. Six single-source female samples were also subjected to MPS analysis but showed very low coverage that could not affect the analysis of the mixed samples. Therefore, the age prediction models of the present study are expected to provide useful investigative leads, especially in mixed male and female samples from sexual assault cases.

Genetic Analysis of Mingrelians Reveals Long-Term Continuity of Populations in Western Georgia (Caucasus)

To elucidate the population history of the Caucasus, we conducted a survey of genetic diversity in Samegrelo (Mingrelia), western Georgia. We collected DNA samples and genealogical information from 485 individuals residing in 30 different locations, the vast majority of whom being Mingrelian speaking. From these DNA samples, we generated mitochondrial DNA (mtDNA) control region sequences for all 485 participants (female and male), Y-short tandem repeat haplotypes for the 372 male participants, and analyzed all samples at nearly 590,000 autosomal single nucleotide polymorphisms (SNPs) plus around 33,000 on the sex chromosomes, with 27,000 SNP removed for missingness, using the GenoChip 2.0+ microarray. The resulting data were compared with those from populations from Anatolia, the Caucasus, the Near East, and Europe. Overall, Mingrelians exhibited considerable mtDNA haplogroup diversity, having high frequencies of common West Eurasian haplogroups (H, HV, I, J, K, N1, R1, R2, T, U, and W. X2) and low frequencies of East Eurasian haplogroups (A, C, D, F, and G). From a Y-chromosome standpoint, Mingrelians possessed a variety of haplogroups, including E1b1b, G2a, I2, J1, J2, L, Q, R1a, and R1b. Analysis of autosomal SNP data further revealed that Mingrelians are genetically homogeneous and cluster with other modern-day South Caucasus populations. When compared with ancient DNA samples from Bronze Age archaeological contexts in the broader region, these data indicate that the Mingrelian gene pool began taking its current form at least by this period, probably in conjunction with the formation of a distinct linguistic community.

Developmental validation of a high-resolution panel genotyping 639 Y-chromosome SNP and InDel markers and its evolutionary features in Chinese populations

Uniparental-inherited haploid genetic marker of Y-chromosome single nucleotide polymorphisms (Y-SNP) have the power to provide a deep understanding of the human evolutionary past, forensic pedigree, and bio-geographical ancestry information. Several international cross-continental or regional Y-panels instead of Y-whole sequencing have recently been developed to promote Y-tools in forensic practice. However, panels based on next-generation sequencing (NGS) explicitly developed for Chinese populations are insufficient to represent the Chinese Y-chromosome genetic diversity and complex population structures, especially for Chinese-predominant haplogroup O. We developed and validated a 639-plex panel including 633 Y-SNPs and 6 Y-Insertion/deletions, which covered 573 Y haplogroups on the Y-DNA haplogroup tree. In this panel, subgroups from haplogroup O accounted for 64.4% of total inferable haplogroups. We reported the sequencing metrics of 354 libraries sequenced with this panel, with the average sequencing depth among 226 individuals being 3,741×. We illuminated the high level of concordance, accuracy, reproducibility, and specificity of the 639-plex panel and found that 610 loci were genotyped with as little as 0.03 ng of genomic DNA in the sensitivity test. 94.05% of the 639 loci were detectable in male-female mixed DNA samples with a mix ratio of 1:500. Nearly all of the loci were genotyped correctly when no more than 25 ng/μL tannic acid, 20 ng/μL humic acid, or 37.5 μM hematin was added to the amplification mixture. More than 80% of genotypes were obtained from degraded DNA samples with a degradation index of 11.76. Individuals from the same pedigree shared identical genotypes in 11 male pedigrees. Finally, we presented the complex evolutionary history of 183 northern Chinese Hans and six other Chinese populations, and found multiple founding lineages that contributed to the northern Han Chinese gene pool. The 639-plex panel proved an efficient tool for Chinese paternal studies and forensic applications.

Comprehensive analyses of genetic diversities and population structure of the Guizhou Dong group based on 44 Y-markers

Background

The non-recombining region of the human Y chromosome (NRY) is a strictly paternally inherited genetic marker and the best material to trace the paternal lineages of populations. Y chromosomal short tandem repeat (Y-STR) is characterized by high polymorphism and paternal inheritance pattern, so it has been widely used in forensic medicine and population genetic research. This study aims to understand the genetic distribution of Y-STRs in the Guizhou Dong population, provide reference data for forensic application, and explore the phylogenetic relationships between the Guizhou Dong population and other comparison populations.

Methods

Based on the allele profile of 44 Y-markers in the Guizhou Dong group, we estimate their allele frequencies and haplotype frequencies. In addition, we also compare the forensic application efficiency of different Y-STR sets in the Guizhou Dong group. Finally, genetic relationships among Guizhou Dong and other reference populations are dissected by the multi-dimensional scaling and the phylogenetic tree.

Results

A total of 393 alleles are observed in 312 Guizhou Dong individuals for these Y-markers, with allele frequencies ranging from 0.0032 to 0.9679. The haplotype diversity and discriminatory capacity for these Y-markers in the Guizhou Dong population are 0.99984 and 0.97440, respectively. The population genetic analyses of the Guizhou Dong group and other reference populations show that the Guizhou Dong group has the closest genetic relationship with the Hunan Dong population, and followed by the Guizhou Tujia population.

Conclusions

In conclusion, these 44 Y-markers can be used as an effective tool for male differentiation in the Guizhou Dong group. The haplotype data in this study not only enrich the Y-STR data of different ethnic groups in China, but also have important significance for population genetics and forensic research.

Origins of East Caucasus Gene Pool: Contributions of Autochthonous Bronze Age Populations and Migrations from West Asia Estimated from Y-Chromosome Data

The gene pool of the East Caucasus, encompassing modern-day Azerbaijan and Dagestan populations, was studied alongside adjacent populations using 83 Y-chromosome SNP markers. The analysis of genetic distances among 18 populations (N = 2216) representing Nakh-Dagestani, Altaic, and Indo-European language families revealed the presence of three components (Steppe, Iranian, and Dagestani) that emerged in different historical periods. The Steppe component occurs only in Karanogais, indicating a recent medieval migration of Turkic-speaking nomads from the Eurasian steppe. The Iranian component is observed in Azerbaijanis, Dagestani Tabasarans, and all Iranian-speaking peoples of the Caucasus. The Dagestani component predominates in Dagestani-speaking populations, except for Tabasarans, and in Turkic-speaking Kumyks. Each component is associated with distinct Y-chromosome haplogroup complexes: the Steppe includes C-M217, N-LLY22g, R1b-M73, and R1a-M198; the Iranian includes J2-M172(×M67, M12) and R1b-M269; the Dagestani includes J1-Y3495 lineages. We propose J1-Y3495 haplogroup's most common lineage originated in an autochthonous ancestral population in central Dagestan and splits up ~6 kya into J1-ZS3114 (Dargins, Laks, Lezgi-speaking populations) and J1-CTS1460 (Avar-Andi-Tsez linguistic group). Based on the archeological finds and DNA data, the analysis of J1-Y3495 phylogeography suggests the growth of the population in the territory of modern-day Dagestan that started in the Bronze Age, its further dispersal, and the microevolution of the diverged population.

Relationship of the gene pool of the Khants with the peoples of Western Siberia, Cis-Urals and the Altai-Sayan Region according to the data on the polymorphism of autosomic locus and the Y-chromosome

Khanty are indigenous Siberian people living on the territory of Western Siberia, mainly on the territory of the Khanty-Mansiysk and Yamalo-Nenets Autonomous Okrugs. The present study is aimed at a comprehensive analysis of the structure of the Khanty gene pool and their comparison with other populations of the indigenous population of Southern and Western Siberia. To address the issues of genetic proximity of the Khanty with other indigenous peoples, we performed genotyping of a wide genomic set of autosomal markers using high-density biochips, as well as an expanded set of SNP and STR markers of the Y-chromosome in various ethnic groups: Khakas, Tuvans, Southern Altaians, Siberian Tatars, Chulyms (Turkic language family) and Kets (Yeniseian language family). The structure of the gene pool of the Khanty and other West Siberian and South Siberian populations was studied using a genome-wide panel of autosomal single nucleotide polymorphic markers and Y-chromosome markers. The results of the analysis of autosomal SNPs frequencies by various methods, the similarities in the composition of the Y-chromosome haplogroups and YSTR haplotypes indicate that the Khanty gene pool is quite specific. When analyzing autosomal SNPs, the Ugrian genetic component completely dominates in both samples (up to 99-100 %). The samples of the Khanty showed the maximum match in IBD blocks with each other, with a sample of the Kets, Chulyms, Tuvans, Tomsk Tatars, Khakas, Kachins, and Southern Altaians. The degree of coincidence of IBD blocks between the Khanty, Kets, and Tomsk Tatars is consistent with the results of the distribution of allele frequencies and common genetic components in these populations. According to the composition of the Y-chromosome haplogroups, the two samples of the Khanty differ significantly from each other. A detailed phylogenetic analysis of various Y-chromosome haplogroups made it possible to describe and clarify the differences in the phylogeny and structure of individual ethnospecific sublines, to determine their relationship, traces of population expansion in the Khanty gene pool. Variants of different haplogroups of the Y-chromosome in the Khanty, Khakas and Tuvans go back to their common ancestral lines. The results of a comparative analysis of male samples indicate a close genetic relationship between the Khanty and Nenets, Komi, Udmurts and Kets. The specificity of haplotypes, the discovery of various terminal SNPs confirms that the Khanty did not come into contact with other ethnic groups for a long time, except for the Nenets, which included many Khanty clans.

Structure and origin of Tuvan gene pool according to autosome SNP and Y-chromosome haplogroups

Tuvans are one of the most compactly living peoples of Southern Siberia, settled mainly in the territory of Tuva. The gene pool of the Tuvans is quite isolated, due to endogamy and a very low frequency of interethnic marriages. The structure of the gene pool of the Tuvans and other Siberian populations was studied using a genome-wide panel of autosomal single nucleotide polymorphic markers and Y-chromosome markers. The results of the analysis of the frequencies of autosomal SNPs by various methods, the similarities in the composition of the Y-chromosome haplogroups and YSTR haplotypes show that the gene pool of the Tuvans is very heterogeneous in terms of the composition of genetic components. It includes the ancient autochthonous Yeniseian component, which dominates among the Chulym Turks and Kets, the East Siberian component, which prevails among the Yakuts and Evenks, and the Far Eastern component, the frequency of which is maximum among the Nivkhs and Udeges. Analysis of the composition of IBD-blocks on autosomes shows the maximum genetic relationship of the Tuvans with the Southern Altaians, Khakas and Shors, who were formed during the settlement of the Turkic groups of populations on the territory of the Altai-Sayan region. A very diverse composition of the Tuvan gene pool is shown for various sublines of Y-chromosomal haplogroups, most of which show strong ethnic specificity. Phylogenetic analysis of individual Y-chromosome haplogroups demonstrates the maximum proximity of the gene pool of the Tuvans with the Altaians, Khakas and Shors. Differences in frequencies of Y-chromosome haplogroups between the Todzhans and Tuvans and a change in the frequencies of haplogroups from south to north associated with the East Asian component were found. The majority of the most frequent Y-chromosome haplogroups in the Tuvans demonstrate the founder effect, the formation age of which is fully consistent with the data on their ethnogenesis.

Weight of evidence of Y-STR matches computed with the discrete Laplace method: Impact of adding a suspect's profile to a reference database

The discrete Laplace method is recommended by multiple parties (including the International Society for Forensic Genetics, ISFG) to estimate the weight of evidence in criminal cases when a suspect's Y-STR profile matches the crime scene Y-STR profile. Unfortunately, modelling the distribution of Y-STR profiles in the population reference database is time-consuming and requires expert knowledge. When the suspect's Y-STR profile is added to the database, as would be the protocol in many cases, the parameters of the discrete Laplace model must be re-estimated. We found that the likelihood ratios with and without adding the suspect's Y-STR profile were almost identical with 1,000 or more Y-STR profiles in the database for Y-STR profiles with 8, 12, and 17 loci. Thus, likelihood ratio calculations can be performed in seconds if an established discrete Laplace model based on at least 1,000 Y-STR profiles is used. A match in a population reference database with 17 Y-STR loci from at least 1,000 male individuals results in a likelihood ratio above 10,000 in approximately 94% of the cases, and above 100,000 in approximately 82% of the cases. We offer free software accessible without restrictions to estimate a discrete Laplace model using a Y-STR reference database and subsequently to calculate likelihood ratios.

High Coverage Mitogenomes and Y-Chromosomal Typing Reveal Ancient Lineages in the Modern-Day Székely Population in Romania

Here we present 115 whole mitogenomes and 92 Y-chromosomal Short Tandem Repeat (STR) and Single Nucleotide Polymorphism (SNP) profiles from a Hungarian ethnic group, the Székelys (in Romanian: Secuii, in German: Sekler), living in southeast Transylvania (Romania). The Székelys can be traced back to the 12th century in the region, and numerous scientific theories exist as to their origin. We carefully selected sample providers that had local ancestors inhabiting small villages in the area of Odorheiu Secuiesc/Székelyudvarhely in Romania. The results of our research and the reported data signify a qualitative leap compared to previous studies since it presents the first complete mitochondrial DNA sequences and Y-chromosomal profiles of 23 STRs from the region. We evaluated the results with population genetic and phylogenetic methods in the context of the modern and ancient populations that are either geographically or historically related to the Székelys. Our results demonstrate a predominantly local uniparental make-up of the population that also indicates limited admixture with neighboring populations. Phylogenetic analyses confirmed the presumed eastern origin of certain maternal (A, C, D) and paternal (Q, R1a) lineages, and, in some cases, they could also be linked to ancient DNA data from the Migration Period (5th-9th centuries AD) and Hungarian Conquest Period (10th century AD) populations.

Reconstruction of the Major Maternal and Paternal Lineages in the Feral New Zealand Kaimanawa Horses

New Zealand has the fourth largest feral horse population in the world. The Kaimanawas (KHs) are feral horses descended from various domestic horse breeds released into the Kaimanawa ranges in the 19th and 20th centuries. Over time, the population size has fluctuated dramatically due to hunting, large-scale farming and forestry. Currently, the herd is managed by an annual round-up, limiting the number to 300 individuals to protect the native ecosystem. Here, we genotyped 96 KHs for uniparental markers (mitochondrial DNA, Y-chromosome) and assessed their genetic similarity with respect to other domestic horses. We show that at least six maternal and six paternal lineages contributed unequally to the KH gene pool, and today's KH population possibly represents two sub-populations. Our results indicate that three horse breeds, namely Welsh ponies, Thoroughbreds and Arabian horses had a major influence in the genetic-makeup of the extant KH population. We show that mitochondrial genetic diversity in KHs (π = 0.00687 ± 0.00355) is closer to that of the Sable Island horses (π = 0.0034 ± 0.00301), and less than other feral horse populations around the world. Our current findings, combined with ongoing genomic research, will provide insight into the population-specific genetic variation and inbreeding among KHs. This will largely advance equine research and improve the management of future breeding programs of these treasured New Zealand horse.

Clinical, cytogenomic, and molecular characterization of isodicentric Y-chromosome and prediction of testicular sperm retrieval outcomes in azoospermic and severe oligozoospermic infertile men

PURPOSE

Sex chromosome abnormalities are associated with male infertility. The aim of this study was to characterize the clinical, cytogenetic, and molecular findings of 12 infertile men with isodicentric Y-chromosome [idic(Y)] abnormalities diagnosed over a period of 13 years.

MATERIALS AND METHODS

Chromosomal analyses of peripheral blood samples were done using standard procedures. Fluorescence in situ hybridization (FISH) analysis was performed on metaphase spreads of the patients. Multiplex polymerase chain reaction (PCR) using several sequence-tagged site (STS) primer sets within the long arm of Y-chromosome was used to detect AZF deletions.The breakpoints and copy number variations (CNV) were identified by array comparative genomic hybridization analysis (aCGH) analysis.The short-stature homeobox (SHOX) gene deletions were verified using multiplex ligation-dependent probe amplification (MLPA) analysis.

RESULTS

Twelve infertile men were diagnosed cytogenetically with idic(Y). The karyotypes of two of the patients were non-mosaic, and the remaining karyotypes showed various degrees of mosaicism. SHOX gene deletion was found in two of the four patients with short stature, and the remaining two patients had shown a 45,X dominant cell line (33.3%). The most common breakpoints for idic(Yq) and idic(Yp) were found to be in Yq11.222 and Yp11.32, respectively. Semen analysis of ten patients (83.3%) demonstrated azoospermia, and the remaining two patients (16.7%) showed severe oligoasthenoteratozoospermia (OAT). In total, 33% (4/12) of idic(Y) patients with or without microsurgical testicular sperm extraction (microTESE) had sperm retrieval.

CONCLUSIONS

Twelve patients with idic(Y) and different breakpoints of Y-chromosome were characterized using multiple detection strategies. Sperm retrieval outcomes of patients either with idic(Yp) or idic(Yq) showed the possibility to find sperm by microTESE.

Ancient Components and Recent Expansion in the Eurasian Heartland: Insights into the Revised Phylogeny of Y-Chromosomes from Central Asia

In the past two decades, studies of Y chromosomal single nucleotide polymorphisms (Y-SNPs) and short tandem repeats (Y-STRs) have shed light on the demographic history of Central Asia, the heartland of Eurasia. However, complex patterns of migration and admixture have complicated population genetic studies in Central Asia. Here, we sequenced and analyzed the Y-chromosomes of 187 male individuals from Kazakh, Kyrgyz, Uzbek, Karakalpak, Hazara, Karluk, Tajik, Uyghur, Dungan, and Turkmen populations. High diversity and admixture from peripheral areas of Eurasia were observed among the paternal gene pool of these populations. This general pattern can be largely attributed to the activities of ancient people in four periods, including the Neolithic farmers, Indo-Europeans, Turks, and Mongols. Most importantly, we detected the consistent expansion of many minor lineages over the past thousand years, which may correspond directly to the formation of modern populations in these regions. The newly discovered sub-lineages and variants provide a basis for further studies of the contributions of minor lineages to the formation of modern populations in Central Asia.

Assessing the utility of quantitative and qualitative metrics in the DNA quantification process of skeletal remains for autosomal and Y-chromosome STR amplification purposes

In historical cases, ancient DNA investigations and missing persons identification, teeth or bone samples are often the only and almost always the best biological material available for DNA typing. On the other hand, DNA obtained from bone material may be characterized by a high degradation index (DI) or its low content, or DNA tests cannot be repeated due to bone piece size limitation. That is often the effect of the environment in which the material was placed and the time during which exposure to unfavorable environmental factors took place. Therefore, it is very important to use appropriate procedures related to STR analysis. For our study, we selected 80 challenging bone samples. The amount of DNA was compared in qPCR using Quantifiler™ Trio DNA Quantification Kit and Investigator® Quantiplex® Pro RGQ. All qPCR results were confirmed by PCR-CE. The results of DNA concentrations and the assigned degradation index (DI) differed significantly within analyzed samples (~10%). Additionally, the Y-chromosome DI also differed from the autosomal DI in the samples. The difference in degradation indexes could explain the lower Y-chromosome amplification success rate compared to autosomal e.g. during human identification process. The results indicate that performing two DNA quantifications with the use of two different kits (primers sets) allows for a much more precise evaluation of the DNA quality and quantity in the isolate. We suggest that at least one of two suggested DNA concentration measurements should be based on an additional determination of the Y chromosome degradation index. Altogether, it allows for rational isolate management, especially when the volume is limited and the sample is unique.

Geographical contrasts of Y-chromosomal haplogroups from wild and domestic goats reveal ancient migrations and recent introgressions

By their paternal transmission, Y-chromosomal haplotypes are sensitive markers of population history and male-mediated introgression. Previous studies identified biallelic single-nucleotide variants in the SRY, ZFY, DDX3Y genes, which in domestic goats identified four major Y-chromosomal haplotypes Y1A, Y1B, Y2A and Y2B with a marked geographic partitioning. Here, we extracted goat Y-chromosomal variants from whole-genome sequences of 386 domestic goats (75 breeds) and 7 wild goat species, which were generated by the VarGoats goat genome project. Phylogenetic analyses indicated domestic haplogroups corresponding to Y1B, Y2A and Y2B, respectively, whereas Y1A is split into Y1AA and Y1AB. All five haplogroups were detected in 26 ancient DNA samples from southeast Europe or Asia. Haplotypes from present-day bezoars are not shared with domestic goats and are attached to deep nodes of the trees and networks. Haplogroup distributions for 186 domestic breeds indicate ancient paternal population bottlenecks and expansions during the migrations into northern Europe, eastern and southern Asia and Africa south of the Sahara. In addition, sharing of haplogroups indicates male-mediated introgressions, most notably an early gene flow from Asian goats into Madagascar and the crossbreeding that in the 19^th century resulted in the popular Boer and Anglo-Nubian breeds. More recent introgressions are those from European goats into the native Korean goat population and from Boer goat into Uganda, Kenya, Tanzania, Malawi and Zimbabwe. This study illustrates the power of the Y-chromosomal variants for reconstructing the history of domestic species with a wide geographic range.

The development of an indel panel for microchimerism detection

OBJECTIVES

The aim of the study was to create a simple assay for microchimerism detection independent of sex and without HLA genotyping.

METHODS

The method is based on detection of insertion or deletions utilizing a multiplex PCR followed by fragment analysis by capillary electrophoresis, and probe-based qPCR assays. A total of 192 samples, taken either before pregnancy, during 1st trimester, or either during 2nd trimester or at miscarriage, obtained from a cohort of 97 female patients with either primary or secondary recurrent pregnancy loss, were screened for fetal microchimerism by the indel panel as well as an existing assay based on detection of the Y-chromosome marker; DYS14.

RESULTS

The overall prevalence of DYS14 positive samples was 29% (55/192) whereas 32% (61/192) tested positive by the indel method. There was an overall agreement of 64% (122/192) between the results obtained by the two methods. A Fisher's Exact test showed no statistic significant difference in the prevalence of microchimerism detected by the two methods at any of the three times of sampling. The distribution of the number of positive wells detected by both methods were compared by a Mann-Whitney U test, which showed no statistically significant difference at any of the three times of sampling.

CONCLUSION

The data indicates that microchimerism can be detected efficiently by the indel method. This makes it possible to detect both female and male cells without the need of HLA-genotyping. Furthermore, the indel method has potential to be implemented as a routine analysis. This will remove the sex bias in future explorations of the role microchimerism plays in health and disease.

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.

The surname structure of Trentino (Italy) and its relationship with dialects and genes

BACKGROUND

Thanks to the availability of rich surname, linguistic and genetic information, together with its geographic and cultural complexity, Trentino (North-Eastern Italy) is an ideal place to test the relationships between genetic and cultural traits.

AIM

We provide a comprehensive study of population structures based on surname and dialect variability and evaluate their relationships with genetic diversity in Trentino.

SUBJECTS AND METHODS

Surname data were collected for 363 parishes, linguistic data for 57 dialects and genetic data for different sets of molecular markers (Y-chromosome, mtDNA, autosomal) in 10 populations. Analyses relied on different multivariate methods and correlation tests.

RESULTS

Besides the expected isolation-by-distance-like patterns (with few local exceptions, likely related to sociocultural instances), we detected a significant and geography-independent association between dialects and surnames. As for molecular markers, only Y-chromosomal STRs seem to be associated with the dialects, although no significant result was obtained. No evidence for correlation between molecular markers and surnames was observed.

CONCLUSION

Surnames act as cultural markers as do other words, although in this context they cannot be used as reliable proxies for genetic variability at a local scale.

YMrCA: Improving Y-chromosomal ancestor time estimation for DNA kinship research

The Y-chromosome is a valuable kinship indicator in family history and forensic research. To reconstruct genealogies, the time to the most recent common ancestor (tMRCA) between paternal relatives can be estimated through Y-STR analysis. Existing models are the stepwise mutation model (SMM, only one-step Y-STR changes) and the infinite allele model (IAM, new allele per Y-STR change). In this study, these mutation models and all existing tMRCA calculators were validated through a genetic-genealogy database containing 1,120 biologically related genealogical pairs confirmed by 46 Y-STRs with known tMRCA (18,109 generations). Consistent under- and overestimation and broad confidence intervals were observed, leading to dubious tMRCA estimates. This is because they do not include individual mutation rates or multi-step changes and ignore hidden multiple, back, or parallel modifications. To improve tMRCA estimation, we developed a user-friendly calculator, the "YMrCA", including all previously mentioned mutation characteristics. After extensive validation, we observed that the YMrCA calculator demonstrated a promising performance. The YMrCA yields a significantly higher tMRCA success rate (96%; +20%) and a lower tMRCA error (7; -3) compared to the mutation models and all online tMRCA calculators. Therefore, YMrCA offers the next step towards more objective tMRCA estimation for DNA kinship research.

Uniparental genetic markers in Native Americans: A summary of all available data from ancient and contemporary populations

OBJECTIVES

The aim of this study was to create a comprehensive summary of available mtDNA and Y-chromosome data for Native Americans from North, Central, and South America, including both modern and ancient DNA. To illustrate the usefulness of this dataset we present a broad picture of the genetic variation for both markers across the Americas.

METHODS

We searched PubMed, ResearchGate, Google Scholar for studies about mtDNA or Y-chromosome variation in Native American populations, including geographic, linguistic, ecological (ecoregion), archeological and chronological information. We used AMOVA to estimate the genetic structure associated with language and ecoregion grouping and Mantel tests to evaluate the correlation between genetic and geographic distances.

RESULTS

Genetic data were obtained from 321 primary sources, including 22,569 individuals from 298 contemporary populations, and 3628 individuals from 202 archeological populations. MtDNA lineages of probable non-Amerindian origin were rare, in contrast with Y-chromosome lineages. Mantel tests showed a statistically significant correlation for the whole continent considering mtDNA but not the Y-chromosome. Genetic structure between groups was always stronger for mtDNA than for the Y-chromosome.

CONCLUSIONS

This study summarizes decades of research conducted in Native American populations for both mtDNA and the Y-chromosome. Continental or sub-continental patterns of variation reveal that most of the genetic variation occurs within populations rather than among linguistic or ecoregional groups, and that isolation by distance is barely detectable in most population sets. The genetic structure among groups was always larger for mtDNA than for the Y-chromosome, suggesting between-sex differences in gene flow.

Analysis of pedigree data and whole-genome sequences in 12 cattle breeds reveals extremely low within-breed Y-chromosome diversity

In this article, we analyzed pedigree information on males from 12 bovine breeds born in France between 2015 and 2019. We report an overall small number of paternal lineages with, for example, a minimal number of ancestors accounting for 95% of the Y‐chromosome pool of their breed ranging from only 2 to 15 individuals. Then, we mined whole‐genome sequence data from 811 sires (2 ≤ n ≤ 510 per breed) and built a median‐joining network using 1411 SNPs. Most branches were breed‐specific and in agreement with the geographic and genetic relatedness of these populations. The within‐breed haplotype diversity was lower than expected based on genealogical information, which supports the existence of major male founder effects predating pedigree recording. In addition, we observed de novo mutation events among the descendants of the same ancestors, which are of interest to define paternal sub‐lineages. Our results pave the way to future studies on the estimation of the effects of Y‐chromosome haplotypes on male reproductive performances and on the conservation of Y‐chromosome diversity.

Y-chromosome variability and genetic history of Commons from Northern Italy

Objectives

Genetic drift and admixture are driving forces in human evolution, but their concerted impact to population evolution in historical times and at a micro‐geographic scale is poorly assessed. In this study we test a demographic model encompassing both admixture and drift to the case of social‐cultural isolates such as the so‐called “Commons.”

Materials and methods

Commons are peculiar institutions of medieval origins whose key feature is the tight relationship between population and territory, mediated by the collective property of shared resources. Here, we analyze the Y‐chromosomal genetic structure of four Commons (for a total of 366 samples) from the Central and Eastern Padana plain in Northern Italy.

Results

Our results reveal that all these groups exhibit patterns of significant diversity reduction, peripheral/outlier position within the Italian/European genetic space and high frequency of Common‐specific haplogroups. By explicitly testing different drift‐admixture models, we show that a drift‐only model is more probable for Central Padana Commons, while additional admixture (~20%) from external population around the same time of their foundation cannot be excluded for the Eastern ones.

Discussion

Building on these results, we suggest central Middle Ages as the most probable age of foundation for three of the considered Commons, the remaining one pointing to late antiquity. We conclude that an admixture‐drift model is particularly useful for interpreting the genetic structure and recent demographic history of small‐scale populations in which social‐cultural features play a significant role.

Male-specific age estimation based on Y-chromosomal DNA methylation

Although DNA methylation variation of autosomal CpGs provides robust age predictive biomarkers, no male-specific age predictor exists based on Y-CpGs yet. Since sex chromosomes play an important role in aging, a Y-chromosome-based age predictor would allow studying male-specific aging effects and would also be useful in forensics. Here, we used blood-based DNA methylation microarray data of 1,057 males from six cohorts aged 15-87 and identified 75 Y-CpGs with an interquartile range of ≥0.1. Of these, 22 and six were significantly hyper- and hypomethylated with age (p(cor)<0.05, Bonferroni), respectively. Amongst several machine learning algorithms, a model based on support vector machines with radial kernel performed best in male-specific age prediction. We achieved a mean absolute deviation (MAD) between true and predicted age of 7.54 years (cor=0.81, validation) when using all 75 Y-CpGs, and a MAD of 8.46 years (cor=0.73, validation) based on the most predictive 19 Y-CpGs. The accuracies of both age predictors did not worsen with increased age, in contrast to autosomal CpG-based age predictors that are known to predict age with reduced accuracy in the elderly. Overall, we introduce the first-of-its-kind male-specific epigenetic age predictor for future applications in aging research and forensics.

Genetic variation of Nigerian cattle inferred from maternal and paternal genetic markers

The African cattle provide unique genetic resources shaped up by both diverse tropical environmental conditions and human activities, the assessment of their genetic diversity will shade light on the mechanism of their remarkable adaptive capacities. We therefore analyzed the genetic diversity of cattle samples from Nigeria using both maternal and paternal DNA markers. Nigerian cattle can be assigned to 80 haplotypes based on the mitochondrial DNA (mtDNA) D-loop sequences and haplotype diversity was 0.985 + 0.005. The network showed two major matrilineal clustering: the dominant cluster constituting the Nigerian cattle together with other African cattle while the other clustered Eurasian cattle. Paternal analysis indicates only zebu haplogroup in Nigerian cattle with high genetic diversity 1.000 ± 0.016 compared to other cattle. There was no signal of maternal genetic structure in Nigerian cattle population, which may suggest an extensive genetic intermixing within the country. The absence of Bos indicus maternal signal in Nigerian cattle is attributable to vulnerability bottleneck of mtDNA lineages and concordance with the view of male zebu genetic introgression in African cattle. Our study shades light on the current genetic diversity in Nigerian cattle and population history in West Africa.

A review of studies examining the association between genetic biomarkers (short tandem repeats and single-nucleotide polymorphisms) and risk of prostate cancer: the need for valid predictive biomarkers

Prostate cancer (PCa) is a challenging polygenic disease because the genes that cause PCa remain largely elusive and are affected by several causal factors. Consequently, research continuously strives to identify a genetic marker which could be used as an indicator to predict the most vulnerable (i.e., predisposed) segments of the population to the disease or for the gene which may be directly responsible for PCa. To enhance the genetic etiology of PCa, this research sought to discover the key studies conducted in this field using data from the main journal publication search engines, as it was hoped that this could shed light on the main research findings from these studies, which in turn could assist in determining these genes or markers. From the research highlighted, the studies primarily used two kinds of markers: short tandem repeats or single-nucleotide polymorphisms. These markers were found to be quite prevalent in all the chromosomes within the research carried out. It also became apparent that the studies differed in both quantity and quality, as well as being conducted in a variety of societies. Links were also determined between the degree and strength of the relationship between these markers and the occurrence of the disease. From the studies identified, most recommended a larger and more diverse survey for the parameters which had not been studied before, as well as an increase in the size of the community (i.e., the population) being studied. This is an indication that work in this field is far from complete, and thus, current research remains committed toward finding genetic markers that can be used clinically for the diagnosis and screening of patients with PCa.

Assessing the Diversity and Population Substructure of Sarda Breed Bucks by Using Mtdna and Y-Chromosome Markers

A sample of 146 Sarda bucks from eight subregions of Sardinia, Italy (Nuorese, Barbagia, Baronia, Ogliastra, Sarrabus, Guspinese, Iglesiente, Sulcis) were characterized for Y-chromosome and mtDNA markers to assess the levels of population substructure. Five polymorphic loci (SRY, AMELY, ZFY, and DDX3Y) on the Y-chromosome were genotyped. The control region of mtDNA was sequenced as a source of complementary information. Analysis of Y-chromosome data revealed the segregation of 5 haplotypes: Y1A (66.43%), Y2 (28.57%), Y1C (3.57%), Y1B1 (0.71%), and Y1B2 (0.71%). High levels of Y-chromosome diversity were observed in populations from Southwest Sardinia. The FST values based on Y-chromosome and mtDNA data were low, although a paternal genetic differentiation was observed when comparing the Nuorese and Barbagia populations (Central Sardinia) with the Sulcis, Iglesiente, and Sarrabus populations (Southern Sardinia). AMOVA analysis supported the lack of population substructure. These results suggest the occurrence of a historical and extensive gene flow between Sarda goat populations from different locations of Sardinia, despite the fact that this island is covered by several large mountain ranges. Introgression with foreign caprine breeds in order to improve milk production might have also contributed to avoiding the genetic differentiation amongst Sarda populations.

fascicularis, and their Indochinese hybrids

Species in the genus Macaca typically live in multimale-multifemale social groups with male macaques exhibiting some of the largest testis: body weight ratios among primates. Males are believed to experience intense levels of sperm competition. Several spermatogenesis genes are located on the Y-chromosome and, interestingly, occasional hybridization between two species has led to the introgression of the rhesus macaque (Macaca mulatta) Y-chromosome deep into the range of the long-tailed macaque (M. fascicularis). These observations have led to the prediction that the successful introgression of the rhesus Y-haplotype is due to functional differences in spermatogenesis genes compared to those of the native long-tailed Y-haplotype. We examine here four Y-chromosomal loci-RBMY, XKRY, and two nearly identical copies of CDY-and their corresponding protein sequences. The genes were surveyed in representative animals from north of, south of, and within the rhesus x long-tailed introgression zone. Our results show a series of non-synonymous amino acid substitutions present between the two Y-haplotypes. Protein structure modeling via I-TASSER revealed different folding patterns between the two species' Y-proteins, and functional predictions via TreeSAAP further reveal physicochemical differences as a result of non-synonymous substitutions. These differences inform our understanding of the evolution of primate Y-proteins involved in spermatogenesis and, in turn, have biomedical implications for human male fertility.

Y-Chromosome Genetic Analysis of Modern Polish Population

The study presents a full analysis of the Y-chromosome variability of the modern male Polish population. It is the first study of the Polish population to be conducted with such a large set of data (2,705 individuals), which includes genetic information from inhabitants of all voivodeships, i.e., the first administrative level, in the country and the vast majority of its counties, i.e., the second level. In addition, the available data were divided into clusters corresponding to more natural geographic regions. Genetic analysis included the estimation of F_ST distances, the visualization with the use of multidimensional scaling plots and analysis of molecular variance. Y-chromosome binary haplogroups were classified and visualized with the use of interpolation maps. Results showed that the level of differentiation within Polish population is quite low, but some differences were indicated. It was confirmed that the Polish population is characterized by a high degree of homogeneity, with only slight genetic differences being observed at the regional level. The use of regional clustering as an alternative to counties and voivodeships provided a more detailed view of the genetic structure of the population. Those regional differences identified in the present study highlighted the need for additional division of the population by cultural and ethnic criteria in such studies rather than just by geographical or administrative regionalization.

The medieval Mongolian roots of Y-chromosomal lineages from South Kazakhstan

Background

The majority of the Kazakhs from South Kazakhstan belongs to the 12 clans of the Senior Zhuz. According to traditional genealogy, nine of these clans have a common ancestor and constitute the Uissun tribe. There are three main hypotheses of the clans’ origin, namely, origin from early Wusuns, from Niru’un Mongols, or from Darligin Mongols. We genotyped 490 samples of South Kazakhs by 35 Y-chromosomal SNPs (single nucleotide polymorphism) and 17 STRs (short tandem repeat). Additionally, 133 samples from citizen science projects were included into the study.

Results

We found that three Uissun clans have unique Y-chromosomal profiles, but the remaining six Uissun clans and one non-Uissun clan share a common paternal gene pool. They share a high frequency (> 40%) of the C2*-ST haplogroup (marked by the SNP F3796), which is associated with the early Niru’un Mongols. Phylogenetic analysis of this haplogroup carried out on 743 individuals from 25 populations of Eurasia has revealed a set of haplotype clusters, three of which contain the Uissun haplotypes. The demographic expansion of these clusters dates back to the 13-fourteenth century, coinciding with the time of the Uissun’s ancestor Maiky-biy known from historical sources. In addition, it coincides with the expansion period of the Mongol Empire in the Late Middle Ages. A comparison of the results with published aDNA (ancient deoxyribonucleic acid) data and modern Y haplogroups frequencies suggest an origin of Uissuns from Niru’un Mongols rather than from Wusuns or Darligin Mongols.

Conclusions

The Y-chromosomal variation in South Kazakh clans indicates their common origin in 13th–14th centuries AD, in agreement with the traditional genealogy. Though genetically there were at least three ancestral lineages instead of the traditional single ancestor. The majority of the Y-chromosomal lineages of South Kazakhstan was brought by the migration of the population related to the medieval Niru’un Mongols.

Y-DNA genetic evidence reveals several different ancient origins in the Brahmin population

The ancient geographical origins of Brahmins-a prominent ethnic group in the Indian subcontinent-have remained controversial for a long time. This study employed the AMOVA (analysis of molecular variance) test to evaluate genetic affinities of this group with thirty populations of Central Asia and Europe. A domestic comparison was performed with fifty non-Brahmin groups in India. The results showed that Brahmins had genetic affinities with several foreign populations and also shared their genetic heritage with several domestic non-Brahmin groups. The study identified the deep ancient origins of Brahmins by tracing their Y-chromosome haplogroups and genetic markers on the Y-DNA phylogenetic tree. It was confirmed that the progenitors of this group emerged from at least 12 different geographic regions of the world. The study concluded that about 83% of the Brahmins in the dataset belonged to four major haplogroups, of which two emerged from Central Asia, one from the Fertile Crescent, and one was of an indigenous Indian origin.

The H-Y Antigen in Embryonic Stem Cells Causes Rejection in Syngeneic Female Recipients

Pluripotent stem cells are promising candidates for cell-based regenerative therapies. To avoid rejection of transplanted cells, several approaches are being pursued to reduce immunogenicity of the cells or modulate the recipient's immune response. These include gene editing to reduce the antigenicity of cell products, immunosuppression of the host, or using major histocompatibility complex-matched cells from cell banks. In this context, we have investigated the antigenicity of H-Y antigens, a class of minor histocompatibility antigens encoded by the Y chromosome, to assess whether the gender of the donor affects the cell's antigenicity. In a murine transplant model, we show that the H-Y antigen in undifferentiated embryonic stem cells (ESCs), as well as ESC-derived endothelial cells, provokes T- and B cell responses in female recipients.

Structure and Sequence of the Sex Determining Locus in Two Wild Populations of Nile Tilapia

In domesticated strains of the Nile tilapia, phenotypic sex has been linked to genetic variants on linkage groups 1, 20 and 23. This diversity of sex-loci might reflect a naturally polymorphic sex determination system in Nile tilapia, or it might be an artefact arising from the process of domestication. Here, we searched for sex-determiners in wild populations from Kpandu, Lake Volta (Ghana-West Africa), and from Lake Koka (Ethiopia-East Africa) that have not been subjected to any genetic manipulation. We analysed lab-reared families using double-digest Restriction Associated DNA sequencing (ddRAD) and analysed wild-caught males and females with pooled whole-genome sequencing (WGS). Strong sex-linked signals were found on LG23 in both populations, and sex-linked signals with LG3 were observed in Kpandu samples. WGS uncovered blocks of high sequence coverage, suggesting the presence of B chromosomes. We confirmed the existence of a tandem amh duplication in LG23 in both populations and determined its breakpoints between the oaz1 and dot1l genes. We found two common deletions of ~5 kb in males and confirmed the presence of both amhY and amh∆Y genes. Males from Lake Koka lack both the previously reported 234 bp deletion and the 5 bp frameshift-insertion that creates a premature stop codon in amh∆Y.

Paternal Origins and Migratory Episodes of Domestic Sheep

The domestication and subsequent global dispersal of livestock are crucial events in human history, but the migratory episodes during the history of livestock remain poorly documented [1-3]. Here, we first developed a set of 493 novel ovine SNPs of the male-specific region of Y chromosome (MSY) by genome mapping. We then conducted a comprehensive genomic analysis of Y chromosome, mitochondrial DNA, and whole-genome sequence variations in a large number of 595 rams representing 118 domestic populations across the world. We detected four different paternal lineages of domestic sheep and resolved, at the global level, their paternal origins and differentiation. In Northern European breeds, several of which have retained primitive traits (e.g., a small body size and short or thin tails), and fat-tailed sheep, we found an overrepresentation of MSY lineages y-HC and y-HB, respectively. Using an approximate Bayesian computation approach, we reconstruct the demographic expansions associated with the segregation of primitive and fat-tailed phenotypes. These results together with archaeological evidence and historical data suggested the first expansion of early domestic hair sheep and the later expansion of fat-tailed sheep occurred ∼11,800-9,000 years BP and ∼5,300-1,700 years BP, respectively. These findings provide important insights into the history of migration and pastoralism of sheep across the Old World, which was associated with different breeding goals during the Neolithic agricultural revolution.

A genome database for a Japanese population of the larvacean Oikopleura dioica

The larvacean Oikopleura dioica is a planktonic chordate and is a tunicate that belongs to the closest relatives to vertebrates. Its simple and transparent body, invariant embryonic cell lineages, and short life cycle of 5 days make it a promising model organism for the study of developmental biology. The genome browser OikoBase was established in 2013 using Norwegian O. dioica. However, genome information for other populations is not available, even though many researchers have studied local populations. In the present study, we sequenced using Illumina and PacBio RSII technologies the genome of O. dioica from a southwestern Japanese population that was cultured in our laboratory for 3 years. The genome of Japanese O. dioica was assembled into 576 scaffold sequences with a total length and N50 length of 56.6 and 1.5 Mb, respectively. A total of 18,743 gene models (transcript models) were predicted in the genome assembly, named OSKA2016. In addition, 19,277 non-redundant transcripts were assembled using RNA-seq data. The OSKA2016 has global sequence similarity of only 86.5% when compared with the OikoBase, highlighting the sequence difference between the two far distant O. dioica populations on the globe. The genome assembly, transcript assembly, and transcript models were incorporated into ANISEED (https://www.aniseed.cnrs.fr/) for genome browsing and BLAST searches. Mapping of reads obtained from male- or female-specific genome libraries yielded male-specific scaffolds in the OSKA2016 and revealed that over 2.6 Mb of sequence were included in the male-specific Y-region. The genome and transcriptome resources from two distinct populations will be useful datasets for developmental biology, evolutionary biology, and molecular ecology using this model organism.

Population Genetics in Malaysia and Japanese Populations Using Power Plex Y23 System

Population flow between Southeast Asian countries and Japan continues to gather pace. Accordingly, the number of foreigners involved in incidents in Japan has markedly increased, which means that forensic dentistry is now increasingly being faced with the need to analyze DNA from persons of non-Japanese extraction. The DNA test currently used for personal identification mainly utilizes short tandem repeats (STRs) on autosomal chromosomes and the Y-chromosome. This test was developed for use in personal identification, not for distinguishing among races; nonetheless, the standard method for personal identification is often used because the procedure has been established. To determine the degree to which racial differences can be distinguished by standard DNA analysis, 23 STRs located on the Y chromosome were investigated in 218 Malay and 426 Japanese males. The frequencies of each STR were calculated in the two populations. The difference in the power of discrimination between the Malay and Japanese populations ranged from a minimum of 0.01 to a maximum of 0.27; the difference in polymorphic information content ranged from 0.01 (minimum) to 0.23 (maximum). No major differences were noted in the polymorphisms in these two Mongoloid populations, but the distributions of the 17 STRs differed significantly. Short tandem repeat types demonstrating a likelihood of racial differences were identified in 14 of the STRs. Race-specific STR types were identified in 10 STRs. These results suggest that the likelihood of Malay or Japanese genetic background can be judged based on Y-chromosome STR test results.

Ancient Patrilineal Lines and Relatively High ECAY Diversity Preserved in Indigenous Horses Revealed With Novel Y-Chromosome Markers

Extremely low nucleotide diversity of modern horse Y-chromosome has been reported, and only poor phylogenetic resolution could be resulted from limited Y-chromosome markers. In this study, three types of horse Y-chromosome markers, including Single-nucleotide polymorphisms (SNPs), copy number variants (CNVs), and allele-specific CNVs, were developed by screening more than 300 male horses from 23 indigenous Chinese horse populations and 4 imported horse breeds. Fourteen segregating sites including a novel SNP in the AMELY gene were found in approximately 53 kb of male-specific Y-chromosome sequences. CNVs were detected at 11 of 14 sites, while allele-specific CNVs at 6 polymorphic sites in repeated fragments were also determined. The phylogenetic analyses with the SNPs identified in this study and previously published 51 SNPs obtained mainly from European horses showed that indigenous Chinese horses exhibit much deeper divergence than European and Middle Eastern horses, while individuals of Chinese horses with the C allele of the AMELY gene constituted the most ancient group. Via SNPs, CNVs, and allele-specific CNVs, much higher diversity of paternal lines can be detected than those identified with merely SNPs. Our results indicated that there are ancient paternal horse lines preserved in southwestern China, which sheds new light on the domestication and immigration of horses, and suggest that the priorities of the conservation should be given to the ancient and rare paternal lines. These three marker types provided finer phylogenetic resolution of horse patrilineal lines, and the strategies used in the present study also provide valuable reference for the genetic studies of other mammalian patrilineages.

Paternal heritage of the Han Chinese in Henan province (Central China): high diversity and evidence of in situ Neolithic expansions

Background: Due to their long history, complex admixture processes and large population sizes, more research is required to explore the fine genetic structure of Han populations from different geographic locations of China.Aim: To characterise the paternal genetic structure of the Han Chinese in Henan province, which was once the central living region of the ancient Huaxia population, the precursors of the Han Chinese.Subjects and methods: We sequenced Y chromosomes of 60 males from Zhengzhou, Henan Province, and reconstructed a phylogenetic tree for these samples with age estimation.Results: We observed high diversity of paternal lineages in our collection. We found that the in situ Neolithic expansion of the "Major lineages" contributed to a large portion of the paternal gene pool of the Han population in Henan Province. We also detected a large number of "Minor lineages" that diverged in the Palaeolithic Age.Conclusion: We suggest that the high genetic diversity in the paternal gene pool of modern Han populations is mainly attributed to the reservation of a larger number of lineages that diverged in the Palaeolithic Age, while the recent expansion of limited lineages contributed to the majority of the gene pool of modern Han populations. We propose that such a structure is a basal characteristic for the genetic structure of modern Han populations.

Genetic analysis of 27 Y-chromosomal STR loci in a Zimbabwean Shona ethnic group

Buccal swabs from 200 unrelated Zimbabwean males were collected from voluntary participants located in Harare province. The 5-dye SureID® 27Y Human STR Identification Kit was used to perform multiplex polymerase chain reactions (PCR) and generate Y-chromosomal DNA profiles. This kit targets markers DYS456, DYS576, DYS570, DYS481, DYF387S1, DYS627, DYS393, DYS391, DYS390, DYS635, DYS449, DYS533, DYS438, DYS389I, DYS448, DYS389II, DYS19, GATA_H4, DYS518, DYS458, DYS460, DYS437, DYS439, DYS392, and DYS385, similar to the Yfiler® Plus Amplification Kit. A total of 161 haplotypes were generated with the PowerPlex® Y system, whereas 159 complete haplotypes were generated for the Yfiler® Plus system. Haplotype Discrimination Capacity (DC) with the Yfiler® Plus system was determined to be 0.9686, while the Genetic Diversity (GD) of the targeted loci ranged from 0.03748 at DYS392 to 0.867239 at DYS449. One haplotype contained the triallelic pattern 37, 38, and 39 at DYS387S1. In addition, marker DYS387S1 and marker DYS385 had 13 counts of microvariant alleles overall, while 9 null allele counts were noted at marker DYS448. Genetic distances between our population data and 22 other data sets from African countries and people of African descent were estimated and results showed significant genetic variation.

Identification of Y-chromosome scaffolds of the Queensland fruit fly reveals a duplicated gyf gene paralogue common to many Bactrocera pest species

Bactrocera tryoni (Queensland fruit fly) are polyphagous horticultural pests of eastern Australia. Heterogametic males contain a sex-determining Y-chromosome thought to be gene poor and repetitive. Here, we report 39 Y-chromosome scaffolds (~700 kb) from B. tryoni identified using genotype-by-sequencing data and whole-genome resequencing. Male diagnostic PCR assays validated eight Y-scaffolds, and one (Btry4096) contained a novel gene with five exons that encode a predicted 575 amino acid protein. The Y-gene, referred to as typo-gyf, is a truncated Y-chromosome paralogue of X-chromosome gene gyf (1773 aa). The Y-chromosome contained ~41 copies of typo-gyf, and expression occurred in male flies and embryos. Analysis of 13 tephritid transcriptomes confirmed typo-gyf expression in six additional Bactrocera species, including Bactrocera latifrons, Bactrocera dorsalis and Bactrocera zonata. Molecular dating estimated typo-gyf evolved within the past 8.02 million years (95% highest posterior density 10.56-5.52 million years), after the split with Bactrocera oleae. Phylogenetic analysis also highlighted complex evolutionary histories among several Bactrocera species, as discordant nuclear (116 genes) and mitochondrial (13 genes) topologies were observed. B. tryoni Y-sequences may provide useful sites for future transgene insertions, and typo-gyf could act as a Y-chromosome diagnostic marker for many Bactrocera species, although its function is unknown.

Ysurnames? The patrilineal Y-chromosome and surname correlation for DNA kinship research

The Y-chromosome is a widely studied and useful small part of the genome providing different applications for interdisciplinary research. In many (Western) societies, the Y-chromosome and surnames are paternally co-inherited, suggesting a corresponding Y-haplotype for every namesake. While it has already been observed that this correlation may be disrupted by a false-paternity event, adoption, anonymous sperm donor or the co-founding of surnames, extensive information on the strength of the surname match frequency (SMF) with the Y-chromosome remains rather unknown. For the first time in Belgium and the Netherlands, we were able to study this correlation using 2,401 males genotyped for 46 Y-STRs and 183 Y-SNPs. The SMF was observed to be dependent on the number of Y-STRs analyzed, their mutation rates and the number of Y-STR differences allowed for a kinship. For a perfect match, the Yfiler® Plus and our in-house YForGen kit gave a similar high SMF of 98%, but for non-perfect matches, the latter could overall be identified as the best kit. The SMF generally increased due to less mismatches when encountering [1] deep Y-subhaplogroups, [2] less frequently occurring surnames, and [3] small geographical distances between relatives. This novel information enabled the design of a surname prediction model based on genetic and geographical distances of a kinship. The prediction model has an area under the curve (AUC) of 0.9 and is therefore useable for DNA kinship priority listing in estimation applications like forensic familial searching.

Profiling of 17 Y-STR loci in Mazandaran and Gilan provinces of Iran

Background/aim

The Y-chromosome mainly consists of heterochromatin regions that have a father-to-son inheritance. Short tandem repeat polymorphic (STRP) markers distributed all over the chromosome provide the opportunity for investigations in forensic medicine and ancestral lineage studies. Due to the existence of wide varieties of geographical and ethnic groups in Iran, studying Y-STRP markers is necessary for further applications. Here we investigated the provinces of Mazandaran and Gilan for the first time.

Materials and methods

Samples included 119 and 90 unrelated males from Mazandaran and Gilan, respectively. Using a PCR amplification kit, 17 Y-STRP markers were amplified and genotyping was conducted by capillary electrophoresis. Allele frequency, haplotype diversity (HD), and haplotype discrimination capacity (DC) were calculated. The populations were compared together and to neighboring countries including Afghanistan and Azerbaijan by F_ST index.

Results

A total of 204 unique haplotypes were observed. No uniqueness was observed between the two provinces. HD was 0.9993 and 0.9998 in Mazandaran and Gilan, respectively. DC was 0.9666 and 0.9888 for Mazandaran and Gilan, respectively. DYS385b and DYS391 had the most and least polymorphic content in both provinces, respectively. There was not a significant difference between these two provinces (F_ST = 0.0006 and P = 0.00) and neighboring countries.

Conclusion

The results highlight the effectiveness of these Y-STRP markers for male discrimination in the north of Iran. Using additional markers along with extended sample size would provide a better opportunity for removing matched haplotypes and introducing the best polymorphic markers in this specific population.

Novel Y-Chromosome Long Non-Coding RNAs Expressed in Human Male CNS During Early Development

Global microarray gene expression analyses previously demonstrated differences in female and male embryos during neurodevelopment. In particular, before sexual maturation of the gonads, the differences seem to concentrate on the expression of genes encoded on the X- and Y-chromosomes. To investigate genome-wide differences in expression during this early developmental window, we combined high-resolution RNA sequencing with qPCR to analyze brain samples from human embryos during the first trimester of development. Our analysis was tailored for maximum sensitivity to discover Y-chromosome gene expression, but at the same time, it was underpowered to detect X-inactivation escapees. Using this approach, we found that 5 out of 13 expressed gametolog pairs showed unbalanced gene dosage, and as a consequence, a male-biased expression. In addition, we found six novel non-annotated long non-coding RNAs on the Y-chromosome with conserved expression patterns in newborn chimpanzee. The tissue specific and time-restricted expression of these long non-coding RNAs strongly suggests important functions during central nervous system development in human males.

Phylogeography of Y-chromosome haplogroup Q1a1a-M120, a paternal lineage connecting populations in Siberia and East Asia

Background: Previous studies have suggested that the human Y-chromosome haplogroup Q1a1a-M120, a widespread paternal lineage in East Asian populations, originated in South Siberia. However, much uncertainty remains regarding the origin, diversification, and expansion of this paternal lineage.Aim: To explore the origin and diffusion of paternal Q-M120 lineages in East Asia.Subjects and methods: The authors generated 26 new Y chromosome sequences of Q-M120 males and co-analysed 45 Y chromosome sequences of this haplogroup. A highly-revised phylogenetic tree of haplogroup Q-M120 with age estimates was reconstructed. Additionally, a comprehensive phylogeographic analysis of this lineage was performed including 15,007 samples from 440 populations in eastern Eurasia.Results: An ancient connection of this lineage with populations in Siberia was revealed. However, this paternal lineage experienced an in-situ expansion between 5000 and 3000 years ago in northwestern China. Ancient populations with high frequencies of Q-M120 were involved in the formation of ancient Huaxia populations before 2000 years ago; this haplogroup eventually became one of the founding paternal lineages of modern Han populations.Conclusion: This study provides a clear pattern of the origin and diffusion process of haplogroup Q1a1a-M120, as well as the role of this paternal lineage during the formation of ancient Huaxia populations and modern Han populations.

Genes located in Y-chromosomal regions important for male fertility show altered transcript levels in cryptorchidism and respond to curative hormone treatment

Background

Undescended (cryptorchid) testes in patients with defective mini-puberty and low testosterone levels contain gonocytes that fail to differentiate normally, which impairs the development of Ad spermatogonia and ultimately leads to adult infertility. Treatment with the gonadotropin-releasing hormone agonist GnRHa increases luteinizing hormone and testosterone and rescues fertility in the majority of pathological cryptorchid testes. Several Y-chromosomal genes in the male-specific Y region (MSY) are essential for spermatogenesis, testis development and function, and are associated with azoospermia, infertility and cryptorchidism. In this study, we analyzed the expression of MSY genes in testes with Ad spermatogonia (low infertility risk patients) as compared to testes lacking Ad spermatogonia (high infertility risk) before and after curative GnRHa treatment, and in correlation to their location on the Y-chromosome.

Results

Twenty genes that are up- or down-regulated in the Ad- group are in the X-degenerate or the ampliconic region, respectively. GnRHa treatment increases mRNA levels of 14 genes in the ampliconic region and decreases mRNA levels of 10 genes in the X-degenerate region.

Conclusion

Our findings implicate Y-chromosomal genes, including USP9Y, UTY, TXLNGY, RBMY1B, RBMY1E, RBMY1J and TSPY4, some of which are known to be important for spermatogenesis, in the curative hormonal treatment of cryptorchidism-induced infertility.

The Red Fox Y-Chromosome in Comparative Context

While the number of mammalian genome assemblies has proliferated, Y-chromosome assemblies have lagged behind. This discrepancy is caused by biological features of the Y-chromosome, such as its high repeat content, that present challenges to assembly with short-read, next-generation sequencing technologies. Partial Y-chromosome assemblies have been developed for the cat (Feliscatus), dog (Canislupusfamiliaris), and grey wolf (Canislupuslupus), providing the opportunity to examine the red fox (Vulpesvulpes) Y-chromosome in the context of closely related species. Here we present a data-driven approach to identifying Y-chromosome sequence among the scaffolds that comprise the short-read assembled red fox genome. First, scaffolds containing genes found on the Y-chromosomes of cats, dogs, and wolves were identified. Next, analysis of the resequenced genomes of 15 male and 15 female foxes revealed scaffolds containing male-specific k-mers and patterns of inter-sex copy number variation consistent with the heterogametic chromosome. Analyzing variation across these two metrics revealed 171 scaffolds containing 3.37 Mbp of putative Y-chromosome sequence. The gene content of these scaffolds is consistent overall with that of the Y-chromosome in other carnivore species, though the red fox Y-chromosome carries more copies of BCORY2 and UBE1Y than has been reported in related species and fewer copies of SRY than in other canids. The assignment of these scaffolds to the Y-chromosome serves to further characterize the content of the red fox draft genome while providing resources for future analyses of canid Y-chromosome evolution.

Analysis of biogeographic ancestry reveals complex genetic histories for indigenous communities of St. Vincent and Trinidad

OBJECTIVES

From a genetic perspective, relatively little is known about how mass emigrations of African, European, and Asian peoples beginning in the 16th century affected Indigenous Caribbean populations. Therefore, we explored the impact of serial colonization on the genetic variation of the first Caribbean islanders.

MATERIALS AND METHODS

Sixty-four members of St. Vincent's Garifuna Community and 36 members of Trinidad's Santa Rosa First People's Community (FPC) of Arima were characterized for mitochondrial DNA and Y-chromosome diversity via direct sequencing and targeted SNP and STR genotyping. A subset of 32 Garifuna and 18 FPC participants were genotyped using the GenoChip 2.0 microarray. The resulting data were used to examine genetic diversity, admixture, and sex biased gene flow in the study communities.

RESULTS

The Garifuna were most genetically comparable to African descendant populations, whereas the FPC were more similar to admixed American groups. Both communities also exhibited moderate frequencies of Indigenous American matrilines and patrilines. Autosomal SNP analysis indicated modest Indigenous American ancestry in these populations, while both showed varying degrees of African, European, South Asian, and East Asian ancestry, with patterns of sex-biased gene flow differing between the island communities.

DISCUSSION

These patterns of genetic variation are consistent with historical records of migration, forced, or voluntary, and suggest that different migration events shaped the genetic make-up of each island community. This genomic study is the highest resolution analysis yet conducted with these communities, and provides a fuller understanding of the complex bio-histories of Indigenous Caribbean peoples in the Lesser Antilles.

Paternal lineage of the Berbers from Aurès in Algeria: estimate of their genetic variation

Background: Aurès is a vast territory in the east of Algeria, characterised by its traditional Berber settlement which has preserved its language and its rich history; its name goes back to antiquity and before the Roman conquest it was part of the territory of ancient Numidia. The Chaoui people in this region are one of Algeria's largest Berber groups. Aim: The aims were to investigate the level of genetic diversity of the Berbers of Aurès through the analysis of the paternal gene pool and to estimate the percentage of genetic variation among different geographical regions and linguistic groups from Algeria. Subjects and methods: Twenty-three Y-STRs were genotyped in a sample of 218 unrelated males of the Berbers of Aurès. Algorithms were used to estimate the Y-chromosome haplogroups. Genetic distance, non-metric MDS and AMOVA were used to analyse the genetic relationships between sample groups. Results: The paternal lineage of this sample of the Aurès region did not exhibit strong signals of differentiation with other samples from North-central, Northwest, and South Algeria. However, significant differences were found within this sample, demonstrating a high degree of heterogeneity. Conclusion: The results demonstrate that Aurès people are isolated and closed, but nevertheless have quite different genetic profiles.

The Y-linked proto-oncogene TSPY contributes to poor prognosis of the male hepatocellular carcinoma patients by promoting the pro-oncogenic and suppressing the anti-oncogenic gene expression

Background

Liver cancer is one of the major causes of cancer death worldwide, with significantly higher incidence and mortality among the male patients. Although sex hormones and their receptors could contribute to such sex differences, the story is incomplete. Genes on the male-specific region of the Y chromosome could play a role(s) in this cancer. TSPY is the putative gene for the gonadoblastoma locus on the Y chromosome (GBY) that is ectopically expressed in a subset of male hepatocellular carcinomas (HCCs). Although various studies showed that TSPY expression is associated with poor prognosis in the patients and its overexpression promotes cell proliferation of various cancer cell lines, it remains unclear how TSPY contributes to the clinical outcomes of the HCC patients. Identifying the downstream genes and pathways of TSPY actions would provide novel insights on its contribution(s) to male predominance in this deadly cancer.

Results

To determine the effects of TSPY on HCC, a TSPY transgene was introduced to the HCC cell line, HuH-7, and studied with RNA-Seq transcriptome analysis. The results showed that TSPY upregulates various genes associated with cell-cycle and cell-viability, and suppresses cell-death related genes. To correlate the experimental observations with those of clinical specimens, transcriptomes of male HCCs with high TSPY expression were analyzed with reference to those with silent TSPY expression from the Cancer Genome Atlas (TCGA). The comparative analysis identified 49 genes, which showed parallel expression patterns between HuH-7 cells overexpressing TSPY and clinical specimens with high TSPY expression. Among these 49 genes, 16 likely downstream genes could be associated with survival rates in HCC patients. The major upregulated targets were cell-cycle related genes and growth factor receptor genes, including CDC25B and HMMR, whose expression levels are negatively correlated with the patient survival rates. In contrast, PPARGC1A, SLC25A25 and SOCS2 were downregulated with TSPY expression, and possess favorable prognoses for HCC patients.

Conclusion

We demonstrate that TSPY could exacerbate the oncogenesis of HCC by differentially upregulate the expression of pro-oncogenic genes and downregulate those of anti-oncogenic genes in male HCC patients, thereby contributing to the male predominance in this deadly cancer.

Electronic supplementary material

The online version of this article (10.1186/s13578-019-0287-x) contains supplementary material, which is available to authorized users.

Site-specific transgenesis of the Drosophila melanogaster Y-chromosome using CRISPR/Cas9

Despite the importance of Y-chromosomes in evolution and sex determination, their heterochromatic, repeat-rich nature makes them difficult to sequence (due, in part, to ambiguities in sequence alignment and assembly) and to genetically manipulate. Therefore, they generally remain poorly understood. For example, the Drosophila melanogaster Y-chromosome, one of the most extensively studied Y-chromosomes, is widely heterochromatic and composed mainly of highly repetitive sequences, with only a handful of expressed genes scattered throughout its length. Efforts to insert transgenes on this chromosome have thus far relied on either random insertion of transposons (sometimes harbouring 'landing sites' for subsequent integrations) with limited success or on chromosomal translocations, thereby limiting the types of Y-chromosome-related questions that could be explored. Here, we describe a versatile approach to site-specifically insert transgenes on the Y-chromosome in D. melanogaster via CRISPR/Cas9-mediated homology-directed repair. We demonstrate the ability to insert, and detect expression from, fluorescently marked transgenes at two specific locations on the Y-chromosome, and we utilize these marked Y-chromosomes to detect and quantify rare chromosomal nondisjunction effects. Finally, we discuss how this Y-docking technique could be adapted to other insects to aid in the development of genetic control technologies for the management of insect disease vectors and pests.