The human Y chromosome: an evolutionary marker comes of age

Ancestry and different rates of suicide and homicide in European countries: A study with population-level data

INTRODUCTION

There are large differences in suicide rates across Europe. The current study investigated the relationship of suicide and homicide rates in different countries of Europe with ancestry as it is defined with the haplotype frequencies of Y-DNA and mtDNA.

MATERIAL AND METHODS

The mortality data were retrieved from the WHO online database. The genetic data were retrieved from http://www.eupedia.com. The statistical analysis included Forward Stepwise Multiple Linear Regression analysis and Pearson Correlation Coefficient (R).

RESULTS

In males, N and R1a Y-DNA haplotypes were positively related to both homicidal and suicidal behaviors while I1 was negatively related. The Q was positively related to the homicidal rate. Overall, 60-75% of the observed variance was explained. L, J and X mtDNA haplogroups were negatively related with suicide in females alone, with 82-85% of the observed variance described.

DISCUSSION

The current study should not be considered as a study of genetic markers but rather a study of human ancestry. Its results could mean that research on suicidality has a strong biological but locally restricted component and could be limited by the study population; generalizability of the results at an international level might not be possible. Further research with patient-level data are needed to verify whether these haplotypes could serve as biological markers to identify persons at risk to commit suicide or homicide and whether biologically-determined ancestry could serve as an intermediate grouping method or even as an endophenotype in suicide research.

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.

Copy number gain of VCX, X-linked multi-copy gene, leads to cell proliferation and apoptosis during spermatogenesis

Male factor infertility affects one-sixth of couples worldwide, and non-obstructive azoospermia (NOA) is one of the most severe forms. In recent years there has been increasing evidence to implicate the participation of X chromosome in the process of spermatogenesis. To uncover the roles of X-linked multi-copy genes in spermatogenesis, we performed systematic analysis of X-linked gene copy number variations (CNVs) and Y chromosome haplogrouping in 447 idiopathic NOA patients and 485 healthy controls. Interestingly, the frequency of individuals with abnormal level copy of Variable charge, X-linked (VCX) was significantly different between cases and controls after multiple test correction (p = 5.10 × 10⁻⁵). To discriminate the effect of gain/loss copies in these genes, we analyzed the frequency of X-linked multi-copy genes in subjects among subdivided groups. Our results demonstrated that individuals with increased copy numbers of Nuclear RNA export factor 2 (NXF2) (p = 9.21 × 10⁻⁸) and VCX (p = 1.97 × 10⁻⁴) conferred the risk of NOA. In vitro analysis demonstrated that increasing copy number of VCX could upregulate the gene expression and regulate cell proliferation and apoptosis. Our study establishes a robust association between the VCX CNVs and NOA risk.

The Geographic Origins of Ethnic Groups in the Indian Subcontinent: Exploring Ancient Footprints with Y-DNA Haplogroups

Several studies have evaluated the movements of large populations to the Indian subcontinent; however, the ancient geographic origins of smaller ethnic communities are not clear. Although historians have attempted to identify the origins of some ethnic groups, the evidence is typically anecdotal and based upon what others have written before. In this study, recent developments in DNA science were assessed to provide a contemporary perspective by analyzing the Y chromosome haplogroups of some key ethnic groups and tracing their ancient geographical origins from genetic markers on the Y-DNA haplogroup tree. A total of 2,504 Y-DNA haplotypes, representing 50 different ethnic groups in the Indian subcontinent, were analyzed. The results identified 14 different haplogroups with 14 geographic origins for these people. Moreover, every ethnic group had representation in more than one haplogroup, indicating multiple geographic origins for these communities. The results also showed that despite their varied languages and cultural differences, most ethnic groups shared some common ancestors because of admixture in the past. These findings provide new insights into the ancient geographic origins of ethnic groups in the Indian subcontinent. With about 2,000 other ethnic groups and tribes in the region, it is expected that more scientific discoveries will follow, providing insights into how, from where, and when the ancestors of these people arrived in the subcontinent to create so many different communities.

Using Y-Chromosomal Haplogroups in Genetic Association Studies and Suggested Implications

Y-chromosomal (Y-DNA) haplogroups are more widely used in population genetics than in genetic epidemiology, although associations between Y-DNA haplogroups and several traits, including cardiometabolic traits, have been reported. In apparently homogeneous populations defined by principal component analyses, there is still Y-DNA haplogroup variation which will result from population history. Therefore, hidden stratification and/or differential phenotypic effects by Y-DNA haplogroups could exist. To test this, we hypothesised that stratifying individuals according to their Y-DNA haplogroups before testing for associations between autosomal single nucleotide polymorphisms (SNPs) and phenotypes will yield difference in association. For proof of concept, we derived Y-DNA haplogroups from 6537 males from two epidemiological cohorts, Avon Longitudinal Study of Parents and Children (ALSPAC) (n = 5080; 816 Y-DNA SNPs) and the 1958 Birth Cohort (n = 1457; 1849 Y-DNA SNPs), and studied the robust associations between 32 SNPs and body mass index (BMI), including SNPs in or near Fat Mass and Obesity-associated protein (FTO) which yield the strongest effects. Overall, no association was replicated in both cohorts when Y-DNA haplogroups were considered and this suggests that, for BMI at least, there is little evidence of differences in phenotype or SNP association by Y-DNA structure. Further studies using other traits, phenome-wide association studies (PheWAS), other haplogroups and/or autosomal SNPs are required to test the generalisability and utility of this approach.

The oxytocin receptor gene polymorphism -rs237902- is associated with the severity of autism spectrum disorder: A pilot study

INTRODUCTION

Previous studies showed the association of Autism Spectrum Disorder (ASD) and oxytocin receptor (OXTR) gene. We aimed to explore the OXTR gene single nucleotide polymorphisms (SNPs) across the ASD severity categories based on DSM-5.

METHOD

The whole encoding regions of the human OXTR gene were sequenced to identify the SNPs in 100 Turkish children with ASD. Genotypes of detected SNPs were also compared with the Childhood Autism Rating Scale (CARS) scores.

RESULTS

Disease severity of the patients carrying GA and AA genotypes (GA/AA) of rs237902 were found more severe compared to those carrying GG genotype (χ² = 6.456, df = 2, p = .040). This finding was more powerful in boys (χ² = 9.288, df = 2, p = .010). Similarly, GA/AA genotypes of rs237902 were found associated with higher CARS scores in boys (U = 650.5, r = 0.24, p = .021).

CONCLUSION

Significant relationship between the ASD severity categories of DSM-5 and rs237902 was shown for the first time.

Differential distribution of Y-chromosome haplotypes in Swiss and Southern European goat breeds

The analysis of Y-chromosome variation has provided valuable clues about the paternal history of domestic animal populations. The main goal of the current work was to characterize Y-chromosome diversity in 31 goat populations from Central Eastern (Switzerland and Romania) and Southern Europe (Spain and Italy) as well as in reference populations from Africa and the Near East. Towards this end, we have genotyped seven single nucleotide polymorphisms (SNPs), mapping to the SRY, ZFY, AMELY and DDX3Y Y-linked loci, in 275 bucks from 31 populations. We have observed a low level of variability in the goat Y-chromosome, with just five haplotypes segregating in the whole set of populations. We have also found that Swiss bucks carry exclusively Y1 haplotypes (Y1A: 24%, Y1B1: 15%, Y1B2: 43% and Y1C: 18%), while in Italian and Spanish bucks Y2A is the most abundant haplotype (77%). Interestingly, in Carpathian goats from Romania the Y2A haplotype is also frequent (42%). The high Y-chromosome differentiation between Swiss and Italian/Spanish breeds might be due to the post-domestication spread of two different Near Eastern genetic stocks through the Danubian and Mediterranean corridors. Historical gene flow between Southern European and Northern African goats might have also contributed to generate such pattern of genetic differentiation.

Whole Y-chromosome sequences reveal an extremely recent origin of the most common North African paternal lineage E-M183 (M81)

E-M183 (E-M81) is the most frequent paternal lineage in North Africa and thus it must be considered to explore past historical and demographical processes. Here, by using whole Y chromosome sequences from 32 North African individuals, we have identified five new branches within E-M183. The validation of these variants in more than 200 North African samples, from which we also have information of 13 Y-STRs, has revealed a strong resemblance among E-M183 Y-STR haplotypes that pointed to a rapid expansion of this haplogroup. Moreover, for the first time, by using both SNP and STR data, we have provided updated estimates of the times-to-the-most-recent-common-ancestor (TMRCA) for E-M183, which evidenced an extremely recent origin of this haplogroup (2,000-3,000 ya). Our results also showed a lack of population structure within the E-M183 branch, which could be explained by the recent and rapid expansion of this haplogroup. In spite of a reduction in STR heterozygosity towards the West, which would point to an origin in the Near East, ancient DNA evidence together with our TMRCA estimates point to a local origin of E-M183 in NW Africa.

The genetic variation in the R1a clade among the Ashkenazi Levites' Y chromosome

Approximately 300,000 men around the globe self-identify as Ashkenazi Levites, of whom two thirds were previously shown to descend from a single male. The paucity of whole Y-chromosome sequences precluded conclusive identification of this ancestor's age, geographic origin and migration patterns. Here, we report the variation of 486 Y-chromosomes within the Ashkenazi and non-Ashkenazi Levite R1a clade, other Ashkenazi Jewish paternal lineages, as well as non-Levite Jewish and non-Jewish R1a samples. Cumulatively, the emerging profile is of a Middle Eastern ancestor, self-affiliating as Levite, and carrying the highly resolved R1a-Y2619 lineage, which was likely a minor haplogroup among the Hebrews. A star-like phylogeny, coalescing similarly to other Ashkenazi paternal lineages, ~1,743 ybp, suggests it to be one of the Ashkenazi paternal founders; to have expanded as part of the overall Ashkenazi demographic expansion, without special relation to the Levite affiliation; and to have subsequently spread to non-Ashkenazi Levites.

Linking Y-chromosomal short tandem repeat loci to human male impulsive aggression

INTRODUCTION

Men are more susceptible to impulsive behavior than women. Epidemiological studies revealed that the impulsive aggressive behavior is affected by genetic factors, and the male-specific Y chromosome plays an important role in this behavior. In this study, we investigated the association between the impulsive aggressive behavior and Y-chromosomal short tandem repeats (Y-STRs) loci.

METHODS

The collected biologic samples from 271 offenders with impulsive aggressive behavior and 492 healthy individuals without impulsive aggressive behavior were amplified by PowerPlex^RY23 PCR System and the resultant products were separated by electrophoresis and further genotyped. Then, comparisons in allele and haplotype frequencies of the selected 22 Y-STRs were made in the two groups.

RESULTS

Our results showed that there were significant differences in allele frequencies at DYS448 and DYS456 between offenders and controls (p < .05). Univariate analysis further revealed significant frequency differences for alleles 18 and 22 at DYS448 (0.18 vs 0.27, compared to the controls, p = .003, OR=0.57,95% CI=0.39-0.82; 0.03 vs 0.01, compared to the controls, p = .003, OR=7.45, 95% CI=1.57-35.35, respectively) and for allele 17 at DYS456 (0.07 vs 0.14, compared to the controls, p = .006, OR=0.48, 95% CI =0.28-0.82) between two groups. Interestingly, the frequency of haploid haplotype 22-15 on the DYS448-DYS456 (DYS448-DYS456-22-15) was significantly higher in offenders than in controls (0.033 vs 0.004, compared to the control, p = .001, OR = 8.42, 95%CI =1.81-39.24). Moreover, there were no significant differences in allele frequencies of other Y-STRs loci between two groups. Furthermore, the unconditional logistic regression analysis confirmed that alleles 18 and 22 at DYS448 and allele 17 at DYS456 are associated with male impulsive aggression. However, the DYS448-DYS456-22-15 is less related to impulsive aggression.

CONCLUSION

Our results suggest a link between Y-chromosomal allele types and male impulsive aggression.

Interaction between Y chromosome haplogroup O3* and 4-n-octylphenol exposure reduces the susceptibility to spermatogenic impairment in Han Chinese

Certain genetic background (mainly Y chromosome haplogroups, Y-hg) may modify the susceptibility of certain environmental exposure to some diseases. Compared with respective main effects of genetic background or environmental exposure, interactions between them reflect more realistic combined effects on the susceptibility to a disease. To identify the interactions on spermatogenic impairment, we performed Y chromosome haplotyping and measurement of 9 urinary phenols concentrations in 774 infertile males and 520 healthy controls in a Han Chinese population, and likelihood ratio tests were used to examine the interactions between Y-hgs and phenols. Originally, we observed that Y-hg C and Y-hg F^* might modify the susceptibility to male infertility with urinary 4-n-octylphenol (4-n-OP) level (P_inter = 0.005 and 0.019, respectively). Subsequently, based on our results, two panels were tested to identify the possible protective sub-branches of Y-hg F^* to 4-n-OP exposure, and Y-hg O3^* was uncovered to interact with 4-n-OP (P_inter = 0.019). In conclusion, while 4-n-OP shows an adverse effect on spermatogenesis, Y-hg O3^* makes individuals more adaptive to such an effect for maintaining basic reproductive capacity.

Y-STR Haplogroup Diversity in the Jat Population Reveals Several Different Ancient Origins

The Jats represent a large ethnic community that has inhabited the northwest region of India and Pakistan for several thousand years. It is estimated the community has a population of over 123 million people. Many historians and academics have asserted that the Jats are descendants of Aryans, Scythians, or other ancient people that arrived and lived in northern India at one time. Essentially, the specific origin of these people has remained a matter of contention for a long time. This study demonstrated that the origins of Jats can be clarified by identifying their Y-chromosome haplogroups and tracing their genetic markers on the Y-DNA haplogroup tree. A sample of 302 Y-chromosome haplotypes of Jats in India and Pakistan was analyzed. The results showed that the sample population had several different lines of ancestry and emerged from at least nine different geographical regions of the world. It also became evident that the Jats did not have a unique set of genes, but shared an underlying genetic unity with several other ethnic communities in the Indian subcontinent. A startling new assessment of the genetic ancient origins of these people was revealed with DNA science.

The Y chromosome: a blueprint for men's health?

The Y chromosome has long been considered a 'genetic wasteland' on a trajectory to completely disappear from the human genome. The perception of its physiological function was restricted to sex determination and spermatogenesis. These views have been challenged in recent times with the identification of multiple ubiquitously expressed Y-chromosome genes and the discovery of several unexpected associations between the Y chromosome, immune system and complex polygenic traits. The collected evidence suggests that the Y chromosome influences immune and inflammatory responses in men, translating into genetically programmed susceptibility to diseases with a strong immune component. Phylogenetic studies reveal that carriers of a common European lineage of the Y chromosome (haplogroup I) possess increased risk of coronary artery disease. This occurs amidst upregulation of inflammation and suppression of adaptive immunity in this Y lineage, as well as inferior outcomes in human immunodeficiency virus infection. From structural analysis and experimental data, the UTY (Ubiquitously Transcribed Tetratricopeptide Repeat Containing, Y-Linked) gene is emerging as a promising candidate underlying the associations between Y-chromosome variants and the immunity-driven susceptibility to complex disease. This review synthesises the recent structural, experimental and clinical insights into the human Y chromosome in the context of men's susceptibility to disease (with a particular emphasis on cardiovascular disease) and provides an overview of the paradigm shift in the perception of the Y chromosome.

High Y-chromosomal Differentiation Among Ethnic Groups of Dir and Swat Districts, Pakistan

The ethnic groups that inhabit the mountainous Dir and Swat districts of northern Pakistan are marked by high levels of cultural and phenotypic diversity. To obtain knowledge of the extent of genetic diversity in this region, we investigated Y-chromosomal diversity in five population samples representing the three main ethnic groups residing within these districts, including Gujars, Pashtuns and Kohistanis. A total of 27 Y-chromosomal short tandem repeats (Y-STRs) and 331 Y-chromosomal single nucleotide polymorphisms (Y-SNPs) were investigated. In the Y-STRs, we observed very high and significant levels of genetic differentiation in nine of the 10 pairwise between-group comparisons (R_ST 0.179-0.746), and the differences were mirrored in the Y-SNP haplogroup frequency distribution. No genetic differences were found between the two Pashtun subethnic groups Tarklanis and Yusafzais (R_ST = 0.000). Utmankhels, also considered Pashtuns culturally, were not closely related to any of the other population samples (R_ST 0.451-0.746). Thus, our findings provide examples of both associations and dissociations between cultural and genetic legacies. When analyzed within a larger continental-scale context, these five ethnic groups fall mostly outside the previously characterized Y-chromosomal gene pools of the Indo-Pakistani subcontinent. Male founder effects, coupled with culturally and topographically based constraints upon marriage and movement, are likely responsible for the high degree of genetic structure in this region.

Genetic variation and forensic characteristic analysis of 25 STRs of a novel fluorescence co-amplification system in Chinese Southern Shaanxi Han population

We analyzed the genetic polymorphisms of 15 autosomal and 10 Y-chromosomal STR loci in 214 individuals of Han population from Southern Shaanxi of China and studied the genetic relationships between Southern Shaanxi Han and other populations. We observed a total of 150 alleles at 15 autosomal STR loci with the corresponding allelic frequencies ranging from 0.0023 to 0.5210, and the combined power of discrimination and exclusion for the 15 autosomal STR loci were 0.99999999999999998866 and 0.999998491, respectively. For the 10 Y-STR loci, totally 100 different haplotypes were obtained, of which 94 were unique. The discriminatory capacity and haplotype diversity values of the 10 Y-STR loci were 0.9259 and 0.998269, respectively. The results demonstrated high genetic diversities of the 25 STR loci in the population for forensic applications. We constructed neighbor-joining tree and conducted principal component analysis based on 15 autosomal STR loci and conducted multidimensional scaling analysis and constructed neighbor-joining tree based on 10 Y-STR loci. The results of population genetic analyses based on both autosomal and Y-chromosome STRs indicated that the studied Southern Shaanxi Han population had relatively closer genetic relationship with Eastern Han population, and distant relationships with Croatian, Serbian and Moroccan populations.

Y Chromosome Uncovers the Recent Oriental Origin of Modern Stallions

The Y chromosome directly reflects male genealogies, but the extremely low Y chromosome sequence diversity in horses has prevented the reconstruction of stallion genealogies [1, 2]. Here, we resolve the first Y chromosome genealogy of modern horses by screening 1.46 Mb of the male-specific region of the Y chromosome (MSY) in 52 horses from 21 breeds. Based on highly accurate pedigree data, we estimated the de novo mutation rate of the horse MSY and showed that various modern horse Y chromosome lineages split much later than the domestication of the species. Apart from few private northern European haplotypes, all modern horse breeds clustered together in a roughly 700-year-old haplogroup that was transmitted to Europe by the import of Oriental stallions. The Oriental horse group consisted of two major subclades: the Original Arabian lineage and the Turkoman horse lineage. We show that the English Thoroughbred MSY was derived from the Turkoman lineage and that English Thoroughbred sires are largely responsible for the predominance of this haplotype in modern horses.

Bight of Benin: a Maternal Perspective of Four Beninese Populations and their Genetic Implications on the American Populations of African Ancestry

The understanding of the first movements of the ancestral populations within the African continent is still unclear, particularly in West Africa, due to several factors that have shaped the African genetic pool across time. To improve the genetic representativeness of the Beninese population and to better understand the patterns of human settlement inside West Africa and the dynamics of peopling of the Democratic Republic of Benin, we analyzed the maternal genetic variation of 193 Beninese individuals belonging to Bariba, Berba, Dendi, and Fon populations. Results support the oral traditions indicating that the western neighbouring populations have been the ancestors of the first Beninese populations, and the extant genetic structure of the Beninese populations is most likely the result of admixture between populations from neighbouring countries and native people. The present findings highlight how the Beninese populations contributed to the gene pool of the extant populations of some American populations of African ancestry. This strengthens the hypothesis that the Bight of Benin was not only an assembly point for the slave trade during the Trans-Atlantic Slave Trade but also an important slave trapping area.

Human Y-chromosome variation in the genome-sequencing era

The properties of the human Y chromosome - namely, male specificity, haploidy and escape from crossing over - make it an unusual component of the genome, and have led to its genetic variation becoming a key part of studies of human evolution, population history, genealogy, forensics and male medical genetics. Next-generation sequencing (NGS) technologies have driven recent progress in these areas. In particular, NGS has yielded direct estimates of mutation rates, and an unbiased and calibrated molecular phylogeny that has unprecedented detail. Moreover, the availability of direct-to-consumer NGS services is fuelling a rise of 'citizen scientists', whose interest in resequencing their own Y chromosomes is generating a wealth of new data.

Live birth sex ratios and father's geographic origins in Jerusalem, 1964-1976

OBJECTIVE

To examine whether ancestry influenced sex ratios of offspring in a birth cohort before parental antenatal sex selection influenced offspring sex.

METHODS

We measured the sex ratio as the percent of males according to countries of birth of paternal and maternal grandfathers in 91,459 live births from 1964 to 1976 in the Jerusalem Perinatal Study. Confidence limits (CI) were computed based on an expected sex ratio of 1.05, which is 51.4% male.

RESULTS

Of all live births recorded, 51.4% were male. Relative to Jewish ancestry (51.4% males), significantly more males (1,761) were born to Muslim ancestry (54.5, 95% CI = 52.1-56.8, P = 0.01). Among the former, sex ratios were not significantly associated with paternal or maternal age, education, or offspring's birth order. Consistent with a preference for male offspring, the sex ratio decreased despite increasing numbers of births over the 13-year period. Sex ratios were not affected by maternal or paternal origins in North Africa or Europe. However, the offspring whose paternal grandfathers were born in Western Asia included fewer males than expected (50.7, 50.1-51.3, P = 0.02), whether the father was born abroad (50.7) or in Israel (50.8). This was observed for descendents of paternal grandfathers born in Lebanon (47.6), Turkey (49.9), Yemen & Aden (50.2), Iraq (50.5), Afghanistan (50.5), Syria (50.6), and Cyprus (50.7); but not for those from India (51.5) or Iran (51.9). The West Asian group showed the strongest decline in sex ratios with increasing paternal family size.

CONCLUSIONS

A decreased sex ratio associated with ancestry in Western Asia is consistent with reduced ability to bear sons by a subset of Jewish men in the Jerusalem cohort. Lower sex ratios may be because of pregnancy stress, which may be higher in this subgroup. Alternatively, a degrading Y chromosome haplogroup or other genetic or epigenetic differences on male germ lines could affect birth ratios, such as differential exposure to an environmental agent, dietary differences, or stress. Differential stopping behaviors that favor additional pregnancies following the birth of a daughter might exacerbate these lower sex ratios.

Patterns of Inter-Chromosomal Gene Conversion on the Male-Specific Region of the Human Y Chromosome

The male-specific region of the human Y chromosome (MSY) is characterized by the lack of meiotic recombination and it has long been considered an evolutionary independent region of the human genome. In recent years, however, the idea that human MSY did not have an independent evolutionary history begun to emerge with the discovery that inter-chromosomal gene conversion (ICGC) can modulate the genetic diversity of some portions of this genomic region. Despite the study of the dynamics of this molecular mechanism in humans is still in its infancy, some peculiar features and consequences of it can be summarized. The main effect of ICGC is to increase the allelic diversity of MSY by generating a significant excess of clustered single nucleotide polymorphisms (SNPs) (defined as groups of two or more SNPs occurring in close proximity and on the same branch of the Y phylogeny). On the human MSY, 13 inter-chromosomal gene conversion hotspots (GCHs) have been identified so far, involving donor sequences mainly from the X-chromosome and, to a lesser extent, from autosomes. Most of the GCHs are evolutionary conserved and overlap with regions involved in aberrant X–Y crossing-over. This review mainly focuses on the dynamics and the current knowledge concerning the recombinational landscape of the human MSY in the form of ICGC, on how this molecular mechanism may influence the evolution of the MSY, and on how it could affect the information enclosed within a genomic region which, until recently, appeared to be an evolutionary independent unit.

Spermatogenic failure and the Y chromosome

The Y chromosome harbors a number of genes essential for testis development and function. Its highly repetitive structure predisposes this chromosome to deletion/duplication events and is responsible for Y-linked copy-number variations (CNVs) with clinical relevance. The AZF deletions remove genes with predicted spermatogenic function en block and are the most frequent known molecular causes of impaired spermatogenesis (5-10% of azoospermic and 2-5% of severe oligozoospermic men). Testing for this deletion has both diagnostic and prognostic value for testicular sperm retrieval in azoospermic men. The most dynamic region on the Yq is the AZFc region, presenting numerous NAHR hotspots leading to partial losses or gains of the AZFc genes. The gr/gr deletion (a partial AZFc deletion) negatively affects spermatogenic efficiency and it is a validated, population-dependent risk factor for oligozoospermia. In certain populations, the Y background may play a role in the phenotypic expression of partial AZFc rearrangements and similarly it may affect the predisposition to specific deletions/duplication events. Also, the Yp contains a gene array, TSPY1, with potential effect on germ cell proliferation. Despite intensive investigations during the last 20 years on the role of this sex chromosome in spermatogenesis, a number of clinical and basic questions remain to be answered. This review is aimed at providing an overview of the role of Y chromosome-linked genes, CNVs, and Y background in spermatogenesis.

Improved phylogenetic resolution for Y-chromosome Haplogroup O2a1c-002611

Y-chromosome Haplogroup O2a1c-002611 is one of the dominant lineages of East Asians and Southeast Asians. However, its internal phylogeny remains insufficiently investigated. In this study, we genotyped 89 new highly informative single nucleotide polymorphisms (SNPs) in 305 individuals with Haplogroup O2a1c-002611 identified from 2139 Han Chinese males. Two major branches were identified, O2a1c1-F18 and O2a1c2-L133.2 and the first was further divided into two main subclades, O2a1c1a-F11 and O2a1c1b-F449, accounting for 11.13% and 2.20% of Han Chinese, respectively. In Haplogroup O2a1c1a-F11, we also determined seven sublineages with quite different frequency distributions in Han Chinese ranging from 0.187% to 3.553%, implying they might have different demographic history. The reconstructed haplogroup tree for all the major clades within Haplogroup O2a1c-002611 permits better resolution of male lineages in population studies of East Asia and Southeast Asia. The dataset generated in the present study are also valuable for forensic identification and paternity tests in China.

Y-chromosomal sequences of diverse Indian populations and the ancestry of the Andamanese

We present 42 new Y-chromosomal sequences from diverse Indian tribal and non-tribal populations, including the Jarawa and Onge from the Andaman Islands, which are analysed within a calibrated Y-chromosomal phylogeny incorporating South Asian (in total 305 individuals) and worldwide (in total 1286 individuals) data from the 1000 Genomes Project. In contrast to the more ancient ancestry in the South than in the North that has been claimed, we detected very similar coalescence times within Northern and Southern non-tribal Indian populations. A closest neighbour analysis in the phylogeny showed that Indian populations have an affinity towards Southern European populations and that the time of divergence from these populations substantially predated the Indo-European migration into India, probably reflecting ancient shared ancestry rather than the Indo-European migration, which had little effect on Indian male lineages. Among the tribal populations, the Birhor (Austro-Asiatic-speaking) and Irula (Dravidian-speaking) are the nearest neighbours of South Asian non-tribal populations, with a common origin in the last few millennia. In contrast, the Riang (Tibeto-Burman-speaking) and Andamanese have their nearest neighbour lineages in East Asia. The Jarawa and Onge shared haplogroup D lineages with each other within the last ~7000 years, but had diverged from Japanese haplogroup D Y-chromosomes ~53000 years ago, most likely by a split from a shared ancestral population. This analysis suggests that Indian populations have complex ancestry which cannot be explained by a single expansion model.

A study of the Bodrogköz population in north-eastern Hungary by Y chromosomal haplotypes and haplogroups

We have determined the distribution of Y chromosomal haplotypes and haplogroups in population samples from one of the most important areas in north-eastern Hungary from many villages in the Bodrogköz. The Bodrogköz region was chosen due to its isolated nature, because this area was a moorland encircled by the Tisza, Bodrog, and Latorca Rivers and inhabitants of this part of Hungary escaped from both Tatar and Ottoman invasions, which decimated the post-Hungarian Conquest populations in many parts of the country. Furthermore, in the first half of the tenth century, this region served as the Palatial Centre and burial grounds of the Hungarian tribes. It has thus been assumed that the present population in this area is likely to be more similar to the population that lived in the Conquest period. We analysed male-specific markers, 23 Y-STRs and more than 30 Y-SNPs, that reflect the past and recent genetic history. We found that the general haplogroup distribution of the samples showed high genetic similarity to non-Bodrogköz Hungarians and neighbouring populations, despite its sheltered location and historical record. We were able to classify the Y-chromosomal haplogroups into four large groups based on STR mutation events: pre-Roman/Roman ancient lineage, Finno-Ugric speakers arriving into the Carpathian Basin, Migration period admixture, and post-Hungarian Conquest admixture. It is clear that a significantly larger database with deep haplogroup resolution, including ancient DNA data, is required to strengthen this research.

Genetic variation within the Y chromosome is not associated with histological characteristics of the atherosclerotic carotid artery or aneurysmal wall

BACKGROUND AND AIMS

Haplogroup I, a common European paternal lineage of the Y chromosome, is associated with increased risk of coronary artery disease in British men. It is unclear whether this haplogroup or any other haplogroup on the Y chromosome is associated with histological characteristics of the diseased vessel wall in other vascular manifestations of cardiovascular diseases showing a male preponderance.

METHODS

We examined Dutch men undergoing either carotid endarterectomy from the Athero-Express biobank (AE, n = 1217) or open aneurysm repair from the Aneurysm-Express biobank (AAA, n = 393). Upon resolving the Y chromosome phylogeny, each man was assigned to one of the paternal lineages based on combinations of single nucleotide polymorphisms of the male-specific region of the Y chromosome. We examined the associations between the Y chromosome and the histological characteristics of the carotid plaque and aneurysm wall, including lipid content, leukocyte infiltration and intraplaque haemorrhage, in all men.

RESULTS

A majority of men were carriers of either haplogroup I (AE: 28% AAA: 24%) or haplogroup R (AE: 59% AAA: 61%). We found no association between Y chromosomal haplogroups and histological characteristics of plaque collected from carotid arteries or tissue specimens of aneurysms. Moreover, the distribution of frequency for all Y chromosomal haplogroups in both cohorts was similar to that of a general population of Dutch men.

CONCLUSIONS

Our data show that genetic variation on the Y chromosome is not associated with histological characteristics of the plaques from carotid arteries or specimens of aneurysms in men of Dutch origin.

Phylogeography of human Y-chromosome haplogroup Q3-L275 from an academic/citizen science collaboration

Background

The Y-chromosome haplogroup Q has three major branches: Q1, Q2, and Q3. Q1 is found in both Asia and the Americas where it accounts for about 90% of indigenous Native American Y-chromosomes; Q2 is found in North and Central Asia; but little is known about the third branch, Q3, also named Q1b-L275. Here, we combined the efforts of population geneticists and genetic genealogists to use the potential of full Y-chromosome sequencing for reconstructing haplogroup Q3 phylogeography and suggest possible linkages to events in population history.

Results

We analyzed 47 fully sequenced Y-chromosomes and reconstructed the haplogroup Q3 phylogenetic tree in detail. Haplogroup Q3-L275, derived from the oldest known split within Eurasian/American haplogroup Q, most likely occurred in West or Central Asia in the Upper Paleolithic period. During the Mesolithic and Neolithic epochs, Q3 remained a minor component of the West Asian Y-chromosome pool and gave rise to five branches (Q3a to Q3e), which spread across West, Central and parts of South Asia. Around 3–4 millennia ago (Bronze Age), the Q3a branch underwent a rapid expansion, splitting into seven branches, some of which entered Europe. One of these branches, Q3a1, was acquired by a population ancestral to Ashkenazi Jews and grew within this population during the 1st millennium AD, reaching up to 5% in present day Ashkenazi.

Conclusions

This study dataset was generated by a massive Y-chromosome genotyping effort in the genetic genealogy community, and phylogeographic patterns were revealed by a collaboration of population geneticists and genetic genealogists. This positive experience of collaboration between academic and citizen science provides a model for further joint projects. Merging data and skills of academic and citizen science promises to combine, respectively, quality and quantity, generalization and specialization, and achieve a well-balanced and careful interpretation of the paternal-side history of human populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0870-2) contains supplementary material, which is available to authorized users.

Palenque de San Basilio in Colombia: genetic data support an oral history of a paternal ancestry in Congo

The Palenque, a black community in rural Colombia, have an oral history of fugitive African slaves founding a free village near Cartagena in the seventeenth century. Recently, linguists have identified some 200 words in regular use that originate in a Kikongo language, with Yombe, mainly spoken in the Congo region, being the most likely source. The non-recombining portion of the Y chromosome (NRY) and mitochondrial DNA were analysed to establish whether there was greater similarity between present-day members of the Palenque and Yombe than between the Palenque and 42 other African groups (for all individuals,n= 2799) from which forced slaves might have been taken. NRY data are consistent with the linguistic evidence that Yombe is the most likely group from which the original male settlers of Palenque came. Mitochondrial DNA data suggested substantial maternal sub-Saharan African ancestry and a strong founder effect but did not associate Palenque with any particular African group. In addition, based on cultural data including inhabitants' claims of linguistic differences, it has been hypothesized that the two districts of the village (Abajo and Arriba) have different origins, with Arriba founded by men originating in Congo and Abajo by those born in Colombia. Although significant genetic structuring distinguished the two from each other, no supporting evidence for this hypothesis was found.

Mutation Rates and Discriminating Power for 13 Rapidly-Mutating Y-STRs between Related and Unrelated Individuals

Rapidly Mutating Y-STRs (RM Y-STRs) were recently introduced in forensics in order to increase the differentiation of Y-chromosomal profiles even in case of close relatives. We estimate RM Y-STRs mutation rates and their power to discriminate between related individuals by using samples extracted from a wide set of paternal pedigrees and by comparing RM Y-STRs results with those obtained from the Y-filer set. In addition, we tested the ability of RM Y-STRs to discriminate between unrelated individuals carrying the same Y-filer haplotype, using the haplogroup R-M269 (reportedly characterised by a strong resemblance in Y-STR profiles) as a case study. Our results, despite confirming the high mutability of RM Y-STRs, show significantly lower mutation rates than reference germline ones. Consequently, their power to discriminate between related individuals, despite being higher than the one of Y-filer, does not seem to improve significantly the performance of the latter. On the contrary, when considering R-M269 unrelated individuals, RM Y-STRs reveal significant discriminatory power and retain some phylogenetic signal, allowing the correct classification of individuals for some R-M269-derived sub-lineages. These results have important implications not only for forensics, but also for molecular anthropology, suggesting that RM Y-STRs are useful tools for exploring subtle genetic variability within Y-chromosomal haplogroups.

The genetic evidence for human origin of Jivaroan shrunken heads in collections from the Polish museums

Advances in forensic identification using molecular genetics are helpful in resolving some historical mysteries. The aim of this study was to confirm the authenticity of shrunken-head artifacts exhibited by two Polish museums. Shrunken heads, known as tsantsas, were headhunting trophies of South American Indians (Jivaroan). A special preparation preserved their hair and facial appearance. However, it was quite common to offer counterfeit shrunken heads of sloths or monkeys to collectors of curiosities. We sampled small skin specimens of four shrunken-head skin from the museum collection from Warsaw and Krakow, Poland. Following genomic DNA isolation, highly polymorphic short tandem repeats were genotyped using a commercial chemistry and DNA sequencing analyzer. Haplogroups of human Y chromosome were identified. We obtained an informative genetic profile of genomic short tandem repeats from all the samples of shrunken heads. Moreover, amplification of amelogenin loci allowed for sex determination. All four studied shrunken heads were of human origin. In two ones, a shared Y-chromosome haplogroup Q characteristic for Indigenous Americans was detected. Another artifact was counterfeited because Y-chromosome haplogroup I2 was found, characteristic for the Southeastern European origin. Commercial genetic methods of identification can be applied successfully in studies on the origin and authenticity of some unusual collection items.

Y chromosome haplogroups based genome-wide association study pinpoints revelation for interactions on non-obstructive azoospermia

The Y chromosome has high genetic variability with low rates of parallel and back mutations, which make up the most informative haplotyping system. To examine whether Y chromosome haplogroups (Y-hgs) could modify the effects of autosomal variants on non-obstructive azoospermia (NOA), based on our previous genome-wide association study (GWAS), we conducted a genetic interaction analysis in GWAS subjects. Logistic regression analysis demonstrated a protective effect of Y-hg O3e^* on NOA. Then, we explored the potential interaction between Y-hg O3e^* and autosomal variants. Our results demonstrated that there was a suggestively significant interaction between Y-hg O3e^* and rs11135484 on NOA (P_inter = 9.89 × 10⁻⁵). Bioinformatic analysis revealed that genes annotated by significant single nucleotide polymorphisms (SNPs) were mainly enriched in immunological pathways. This is the first study of interactions between Y-hgs and autosomal variants on a genome-wide scale, which addresses the missing heritability in spermatogenic impairment and sheds new light on the pathogenesis of male infertility.

Molecular Genealogy of a Mongol Queen's Family and Her Possible Kinship with Genghis Khan

Members of the Mongol imperial family (designated the Golden family) are buried in a secret necropolis; therefore, none of their burial grounds have been found. In 2004, we first discovered 5 graves belonging to the Golden family in Tavan Tolgoi, Eastern Mongolia. To define the genealogy of the 5 bodies and the kinship among them, SNP and/or STR profiles of mitochondria, autosomes, and Y chromosomes were analyzed. Four of the 5 bodies were determined to carry the mitochondrial DNA haplogroup D4, while the fifth carried haplogroup CZ, indicating that this individual had no kinship with the others. Meanwhile, Y-SNP and Y-STR profiles indicate that the males examined belonged to the R1b-M343 haplogroup. Thus, their East Asian D4 or CZ matrilineal and West Eurasian R1b-M343 patrilineal origins reveal genealogical admixture between Caucasoid and Mongoloid ethnic groups, despite a Mongoloid physical appearance. In addition, Y chromosomal and autosomal STR profiles revealed that the four D4-carrying bodies bore the relationship of either mother and three sons or four full siblings with almost the same probability. Moreover, the geographical distribution of R1b-M343-carrying modern-day individuals demonstrates that descendants of Tavan Tolgoi bodies today live mainly in Western Eurasia, with a high frequency in the territories of the past Mongol khanates. Here, we propose that Genghis Khan and his family carried Y-haplogroup R1b-M343, which is prevalent in West Eurasia, rather than the Y-haplogroup C3c-M48, which is prevalent in Asia and which is widely accepted to be present in the family members of Genghis Khan. Additionally, Tavan Tolgoi bodies may have been the product of marriages between the lineage of Genghis Khan’s Borjigin clan and the lineage of either the Ongud or Hongirad clans, indicating that these individuals were members of Genghis Khan’s immediate family or his close relatives.

Punctuated bursts in human male demography inferred from 1,244 worldwide Y-chromosome sequences

We report the sequences of 1,244 human Y chromosomes randomly ascertained from 26 worldwide populations by the 1000 Genomes Project. We discovered more than 65,000 variants, including single-nucleotide variants, multiple-nucleotide variants, insertions and deletions, short tandem repeats, and copy number variants. Of these, copy number variants contribute the greatest predicted functional impact. We constructed a calibrated phylogenetic tree on the basis of binary single-nucleotide variants and projected the more complex variants onto it, estimating the number of mutations for each class. Our phylogeny shows bursts of extreme expansion in male numbers that have occurred independently among each of the five continental superpopulations examined, at times of known migrations and technological innovations.