High mutation rates have driven extensive structural polymorphism among human Y chromosomes

Why Y matters? The implication of loss of Y chromosome in blood and cancer

Hematopoietic mosaic loss of Y chromosome (mLOY) has emerged as a potential male-specific accelerator of biological aging, increasing the risk of various age-related diseases, including cancer. Importantly, mLOY is not confined to hematopoietic cells; its presence has also been observed in nonhematological cancer cells, with the impact of this presence previously unknown. Recent studies have revealed that, whether occurring in leukocytes or cancer cells, mLOY plays a role in promoting the development of an immunosuppressive tumor microenvironment. This occurs through the modulation of tumor-infiltrating immune cells, ultimately enabling cancer cells to evade the vigilant immune system. In this review, we illuminate recent progress concerning the effects of hematopoietic mLOY and cancer mLOY on cancer progression. Examining cancer progression from the perspective of LOY adds a new layer to our understanding of cancer immunity, promising insights that hold the potential to identify innovative and potent immunotherapy targets for cancer.

Primary azoospermia factor C duplication associated with spermatogenic impariment: a case-control study based on Y-chromosome haplogrouping in a Han Chinese population

BACKGROUND

Azoospermia factor C (AZFc) in the male-specific region of Y-chromosome (MSY) presents wide structure variation mainly due to frequent non-allele homologous recombination, leading to significant copy number variation of the AZFc-linked coding sequences involving in spermatogenesis. A large number of studies had been conducted to investigate the association between AZFc deletions and male infertility in certain Y chromosome genetic backgrounds, however, the influence of primary AZFc duplication on spermatogenesis remained controversial and the cause of the discrepant outcomes is unknown.

METHODS

In the present study, a total of 1,102 unrelated Han Chinese males without any detectable AZF deletions were recruited from 2014 to 2019, including 411 controls with normozoospermia and 691 patients with idiopathic spermatogenic failure. Using multiple paralog ratio tests (PRTs), the structure duplications were classified by the copy number of the AZFc-linked amplicons and genes. The Y-chromosome haplogroup (Y-hg) was categorized by genetyping of MSY-linked polymorphism loci. The association of primary AZFc duplication with spermatogenic phenotype was investigated in males with the same Y-hg.

RESULTS

Within Y-hg O3* group, the frequency of the gr/gr duplication in patients is significantly higher than that of controls (P = 1.29×10-3 , odds ratio (OR) 7.64, 95% confidence interval (CI) 1.79-32.57). Moreover, Y-hg O3* males with the gr/gr duplication presented a significantly lower sperm production compared with non-AZFc duplicated ones (sperm concentration: P = 1.46×10-3 ; total sperm count: P = 1.82 ×10-3 ). The b2/b3 duplication were identified clustered in Y-hg Cα2*, and the significant difference in the distribution was not observed between patients with spermatogenic failure and controls.

CONCLUSION

The results suggest that, in the Han Chinese population, the gr/gr duplication is a predisposing genetic factor for spermatogenic impairment in males harboring Y-hg O3* . Meanwhile, the b2/b3 duplication may be fixed on a yet-unidentified subbranch of Y-hg Cα2* without significantly deleterious effect on spermatogenesis. Our findings provide evidence that the difference in the Y-hg composition may cause the discrepancy on the association of AZFc duplication with spermatogenic failure among the studied populations.

The human Y and inactive X chromosomes similarly modulate autosomal gene expression

Somatic cells of human males and females have 45 chromosomes in common, including the "active" X chromosome. In males the 46th chromosome is a Y; in females it is an "inactive" X (Xi). Through linear modeling of autosomal gene expression in cells from individuals with zero to three Xi and zero to four Y chromosomes, we found that Xi and Y impact autosomal expression broadly and with remarkably similar effects. Studying sex chromosome structural anomalies, promoters of Xi- and Y-responsive genes, and CRISPR inhibition, we traced part of this shared effect to homologous transcription factors-ZFX and ZFY-encoded by Chr X and Y. This demonstrates sex-shared mechanisms by which Xi and Y modulate autosomal expression. Combined with earlier analyses of sex-linked gene expression, our studies show that 21% of all genes expressed in lymphoblastoid cells or fibroblasts change expression significantly in response to Xi or Y chromosomes.

Uncovering the Association Between Complete AZFc Microduplications and Spermatogenic Ability: The First Reported Series

Purpose This article aims to report the first series of men with complete AZFc microduplications and their clinical and reproductive characteristics. Methods We sampled 3000 men who presented for reproductive urology evaluation from 2012-2020, of which 104 men underwent high-resolution Y-chromosome microarray testing, and five men were identified to have complete AZFc microduplications. Medical, surgical, and reproductive histories were obtained. Semen and hormonal parameters as well as response to fertility therapies were recorded. Results Five men were identified as having complete AZFc microduplications. The mean age was 33.75 years, representing 0.2% (5/3000) of men presenting for fertility investigation, 4.8% (5/104) of men undergoing microarray testing, and 21% (5/24) of men with AZFc abnormalities. Two of the men had prior undescended testicles and one had several autoimmune processes. The mean follicle-stimulating hormone (FSH) was 5.5 IU/L, luteinizing hormone (LH) 3.6 IU/L, and testosterone 14.56 nmol/L. One man was azoospermic, one man alternated between severe oligospermia and rare non-motile sperm, one had variable parameters, with one semen analysis demonstrating azoospermia and a second demonstrating a total motile sperm count (TMSC) of 4 ×106, one man was persistently oligospermic with TMSCs ranging 3.96-12.6 ×106, and one man initially had severe oligospermia, with a mean TMSC of 1.5 ×106, which increased to 21.7 ×106 after intervention (varicocele embolization, clomiphene citrate). This last man then fathered a spontaneous pregnancy. Conclusion AZFc complete microduplications are a rare cause of spermatogenic failure but not an uncommon form of AZFc abnormality. Clinically, they represent a heterogeneous group, having a variable reproductive potential. Cases should be managed on an individual basis.

Characterizing the evolution and phenotypic impact of ampliconic Y chromosome regions

A major part of the human Y chromosome consists of palindromes with multiple copies of genes primarily expressed in testis, many of which have been claimed to affect male fertility. Here we examine copy number variation in these palindromes based on whole genome sequence data from 11,527 Icelandic men. Using a subset of 7947 men grouped into 1449 patrilineal genealogies, we infer 57 large scale de novo copy number mutations affecting palindrome 1. This corresponds to a mutation rate of 2.34 × 10-3 mutations per meiosis, which is 4.1 times larger than our phylogenetic estimate of the mutation rate (5.72 × 10-4), suggesting that de novo mutations on the Y are lost faster than expected under neutral evolution. Although simulations indicate a selection coefficient of 1.8% against non-reference copy number carriers, we do not observe differences in fertility among sequenced men associated with their copy number genotype, but we lack statistical power to detect differences resulting from weak negative selection. We also perform association testing of a diverse set of 341 traits to palindromic copy number without any significant associations. We conclude that large-scale palindrome copy number variation on the Y chromosome has little impact on human phenotype diversity.

The human Y and inactive X chromosomes similarly modulate autosomal gene expression

Somatic cells of human males and females have 45 chromosomes in common, including the "active" X chromosome. In males the 46th chromosome is a Y; in females it is an "inactive" X (Xi). Through linear modeling of autosomal gene expression in cells from individuals with zero to three Xi and zero to four Y chromosomes, we found that Xi and Y impact autosomal expression broadly and with remarkably similar effects. Studying sex-chromosome structural anomalies, promoters of Xi- and Y-responsive genes, and CRISPR inhibition, we traced part of this shared effect to homologous transcription factors - ZFX and ZFY - encoded by Chr X and Y. This demonstrates sex-shared mechanisms by which Xi and Y modulate autosomal expression. Combined with earlier analyses of sex-linked gene expression, our studies show that 21% of all genes expressed in lymphoblastoid cells or fibroblasts change expression significantly in response to Xi or Y chromosomes.

Copy number variation of horse Y chromosome genes in normal equine populations and in horses with abnormal sex development and subfertility: relationship of copy number variations with Y haplogroups

Structural rearrangements like copy number variations in the male-specific Y chromosome have been associated with male fertility phenotypes in human and mouse but have been sparsely studied in other mammalian species. Here, we designed digital droplet PCR assays for 7 horse male-specific Y chromosome multicopy genes and SRY and evaluated their absolute copy numbers in 209 normal male horses of 22 breeds, 73 XY horses with disorders of sex development and/or infertility, 5 Przewalski's horses and 2 kulans. This established baseline copy number for these genes in horses. The TSPY gene showed the highest copy number and was the most copy number variable between individuals and breeds. SRY was a single-copy gene in most horses but had 2-3 copies in some indigenous breeds. Since SRY is flanked by 2 copies of RBMY, their copy number variations were interrelated and may lead to SRY-negative XY disorders of sex development. The Przewalski's horse and kulan had 1 copy of SRY and RBMY. TSPY and ETSTY2 showed significant copy number variations between cryptorchid and normal males (P < 0.05). No significant copy number variations were observed in subfertile/infertile males. Notably, copy number of TSPY and ETSTY5 differed between successive male generations and between cloned horses, indicating germline and somatic mechanisms for copy number variations. We observed no correlation between male-specific Y chromosome gene copy number variations and male-specific Y chromosome haplotypes. We conclude that the ampliconic male-specific Y chromosome reference assembly has deficiencies and further studies with an improved male-specific Y chromosome assembly are needed to determine selective constraints over horse male-specific Y chromosome gene copy number and their relation to stallion reproduction and male biology.

Chromosomal polymorphisms have no negative effect on reproductive outcomes after IVF/ICSI-ET/FET

The present study aimed to explore whether chromosomal polymorphisms (CPs) have negative effects on reproductive outcomes of in vitro fertilization/intracytoplasmic sperm injection-embryo transfer (IVF/ICSI-ET)/frozen-thawing embryo transfer (FET)? We conducted a retrospective study consisting of 21,867 assisted reproductive technology treatment cycles, among which, fresh embryo transfer cycles accounted for 10,400, and the rest were FET cycles. According to karyotype of CPs, the former was grouped as: group 1 (male carrier, n = 425), group 2 (female carrier, n = 262), and group 3 (couple without CPs, n = 9713). Accordingly, FET cycles were divided into 3 groups: group 4 (male carrier, n = 298), group 5 (female carrier, n = 311), and group 6 (couple without CPs, n = 10,858). The embryo implantation rate (IR), clinical pregnancy rate (CPR), live birth rate (LBR), and early miscarriage rate (EMR) were compared among the groups. In fresh embryo transfer cycles after IVF/ICSI, there were no significant differences in the infertility duration, BMI, basal FSH, no. of oocyte, no. of 2PN, endometrial thickness on trigger day, serum E2, P, and LH level on trigger day (P > 0.05). The female age, no. of 2PN embryo cleavage, top-quality embryo, and no. of embryo transferred were significantly different among groups (P < 0.05). The IR was 38.8%, 36.2%, and 34.0% in groups 1, 2, and 3, respectively. The CPR was 55.1%, 52.3%, and 49.7%, respectively. The LBR was 36.9%, 37.4%, and 36.4%, respectively. The CPR and LBR showed no significant differences among groups. The IR was lower and the EMR was higher in group 3 than those of groups 1 and 2. Binary logistic regression analysis indicated that female age, no. of embryo transferred, EMT, LH, and P on the trigger day were independently factors associated with CPR. Besides, no. of embryo transferred, and EMT on trigger day were associated with LBR, while the CPs was not related with CPR and LBR after IVF/ICSI-ET. In FET cycles, the infertility duration was similar (P > 0.05), but the female age, BMI, no. of embryo transferred were significantly different among groups (P > 0.05). The IR was 24.3%, 23.6% and 22.3% in group 4, 5, and 6, receptivity. The CPR was 31.8%, 30.9%, and 30.0%, the LBR was 23.8%,26.3%, and 23.8%, while the EMR was 12.6%, 13.1%, 14.4%, respectively. The IR, CPR, EMR, and LBR showed no significant differences among groups (P > 0.05). Binary logistic regression analysis indicated that female age, infertility duration, and no. of embryo transferred were independently factors affecting CPR and LBR after FET. The CPs were not associated with CPR and LBR after FET. The results suggested that uniparental carrying of CPs have no effects on the reproductive outcomes after IVF/ICSI-ET/FET. However, it is not clear whether both parents carrying CPs would affect pregnancy outcome.

Satellite DNAs and human sex chromosome variation

Satellite DNAs are present on every chromosome in the cell and are typically enriched in repetitive, heterochromatic parts of the human genome. Sex chromosomes represent a unique genomic and epigenetic context. In this review, we first report what is known about satellite DNA biology on human X and Y chromosomes, including repeat content and organization, as well as satellite variation in typical euploid individuals. Then, we review sex chromosome aneuploidies that are among the most common types of aneuploidies in the general population, and are better tolerated than autosomal aneuploidies. This is demonstrated also by the fact that aging is associated with the loss of the X, and especially the Y chromosome. In addition, supernumerary sex chromosomes enable us to study general processes in a cell, such as analyzing heterochromatin dosage (i.e. additional Barr bodies and long heterochromatin arrays on Yq) and their downstream consequences. Finally, genomic and epigenetic organization and regulation of satellite DNA could influence chromosome stability and lead to aneuploidy. In this review, we argue that the complete annotation of satellite DNA on sex chromosomes in human, and especially in centromeric regions, will aid in explaining the prevalence and the consequences of sex chromosome aneuploidies.

Sex-limited chromosomes and non-reproductive traits

Sex chromosomes are typically viewed as having originated from a pair of autosomes, and differentiated as the sex-limited chromosome (e.g. Y) has degenerated by losing most genes through cessation of recombination. While often thought that degenerated sex-limited chromosomes primarily affect traits involved in sex determination and sex cell production, accumulating evidence suggests they also influence traits not sex-limited or directly involved in reproduction. Here, we provide an overview of the effects of sex-limited chromosomes on non-reproductive traits in XY, ZW or UV sex determination systems, and discuss evolutionary processes maintaining variation at sex-limited chromosomes and molecular mechanisms affecting non-reproductive traits.

Recurrent inversion polymorphisms in humans associate with genetic instability and genomic disorders

Unlike copy number variants (CNVs), inversions remain an underexplored genetic variation class. By integrating multiple genomic technologies, we discover 729 inversions in 41 human genomes. Approximately 85% of inversions <2 kbp form by twin-priming during L1 retrotransposition; 80% of the larger inversions are balanced and affect twice as many nucleotides as CNVs. Balanced inversions show an excess of common variants, and 72% are flanked by segmental duplications (SDs) or retrotransposons. Since flanking repeats promote non-allelic homologous recombination, we developed complementary approaches to identify recurrent inversion formation. We describe 40 recurrent inversions encompassing 0.6% of the genome, showing inversion rates up to 2.7 × 10^-4 per locus per generation. Recurrent inversions exhibit a sex-chromosomal bias and co-localize with genomic disorder critical regions. We propose that inversion recurrence results in an elevated number of heterozygous carriers and structural SD diversity, which increases mutability in the population and predisposes specific haplotypes to disease-causing CNVs.

Novel Mutations of TSPY1 Gene Associate Spermatogenic Failure Among Men

Etiology of male infertility is intriguing owing to complex genetic regulation of human spermatogenesis and ethnic variations in genetic architecture of human populations. The present study characterizes the role of Y chromosome specific spermatogenic regulator testis-specific protein Y-encoded 1 (TSPY1) gene mutation in spermatogenic failure. This case-control study includes 163 cases of spermatogenic failure and 175 age-matched fertile men as controls. We found five novel base substitutions, namely, MT162695, MN879413, MN889982, MN889983, MN719943, two deletions MN734578 and MN734579, three novel insertions MN719941, MN719942 and MN719944 through Sanger's dideoxy sequencing of TSPY1 gene reading frame. All these mutations exhibited strong association with male infertility. In silico analyses suggest prospective disruption in splice sites and alteration in different isoforms of TSPY1 transcripts and amino acid sequence in TSPY1 protein. The study provides novel evidence in favour of implication of TSPY1 gene in male fertility. The outcome sheds light to get insight into the issue of idiopathic male infertility in Bengali population.

Copy number evolution in simple and complex tandem repeats across the C57BL/6 and C57BL/10 inbred mouse lines

Simple sequence tandem repeats are among the most rapidly evolving compartments of the genome. Some repeat expansions are associated with mammalian disease or meiotic segregation distortion, yet the rates of copy number change across generations are not well known. Here, we use 14 distinct sublineages of the C57BL/6 and C57BL/10 inbred mouse strains, which have been evolving independently over about 300 generations, to estimate the rates of copy number changes in genome-wide tandem repeats. Rates of change varied across repeats and across lines. Notably, CAG, whose expansions in coding regions are associated with many neurological and genetic disorders, was highly stable in copy number, likely indicating stabilizing selection. Rates of change were positively correlated with copy number, but the direction and magnitude of changes varied across lines. Some mouse lines experienced consistent losses or gains across most simple repeats, but this did not correlate with copy number changes in complex repeats. Rates of copy number change were similar between simple repeats and the more abundant complex repeats after normalization by copy number. Finally, the Y-specific centromeric repeat had a fourfold higher rate of change than the homologous centromeric repeat on other chromosomes. Structural differences in satellite complexity, or restriction to the Y chromosome and elevated mutation rates of the male germline, may explain the higher rate of change. Overall, our work underscores the mutational fluidity of long tandem arrays of repeats, and the correlations and constraints between genome-wide tandem repeats, which suggest that turnover is not a completely neutral process.

Ampliconic Genes on the Great Ape Y Chromosomes: Rapid Evolution of Copy Number but Conservation of Expression Levels

Multicopy ampliconic gene families on the Y chromosome play an important role in spermatogenesis. Thus, studying their genetic variation in endangered great ape species is critical. We estimated the sizes (copy number) of nine Y ampliconic gene families in population samples of chimpanzee, bonobo, and orangutan with droplet digital polymerase chain reaction, combined these estimates with published data for human and gorilla, and produced genome-wide testis gene expression data for great apes. Analyzing this comprehensive data set within an evolutionary framework, we, first, found high inter- and intraspecific variation in gene family size, with larger families exhibiting higher variation as compared with smaller families, a pattern consistent with random genetic drift. Second, for four gene families, we observed significant interspecific size differences, sometimes even between sister species—chimpanzee and bonobo. Third, despite substantial variation in copy number, Y ampliconic gene families’ expression levels did not differ significantly among species, suggesting dosage regulation. Fourth, for three gene families, size was positively correlated with gene expression levels across species, suggesting that, given sufficient evolutionary time, copy number influences gene expression. Our results indicate high variability in size but conservation in gene expression levels in Y ampliconic gene families, significantly advancing our understanding of Y-chromosome evolution in great apes.

Structural and numerical Y chromosomal variations in elderly men identified through multiplex ligation-dependent probe amplification

Human Y chromosomes frequently acquire structural and numerical alterations. Known alterations include germline copy-number variations (CNVs) in the azoospermia factor (AZF) region and somatic mosaic loss of the Y chromosome (mLOY). Here, we explored Y chromosomal variations in 160 Japanese men aged 75-90 years. Multiplex ligation-dependent probe amplification (MLPA) identified ten types of AZF-linked CNVs in 77 men and mLOY of various degrees in 37. Seventeen men carried both a CNV and mLOY. MLOY levels estimated by MLPA were closely correlated with those determined by droplet digital PCR. No association was found between AZF-linked CNVs and the frequency or levels of mLOY. These results emphasize the high frequency and large inter-individual variability of AZF-linked CNVs and mLOY, and demonstrate the usefulness of MLPA in the detection of these variations. More importantly, this study provides the first evidence that AZF-linked CNVs do not increase the risk of aging-related mLOY.

Human AZFb deletions cause distinct testicular pathologies depending on their extensions in Yq11 and the Y haplogroup: new cases and review of literature

Genomic AZFb deletions in Yq11 coined “classical” (i.e. length of Y DNA deletion: 6.23 Mb) are associated with meiotic arrest (MA) of patient spermatogenesis, i.e., absence of any postmeiotic germ cells. These AZFb deletions are caused by non-allelic homologous recombination (NAHR) events between identical sequence blocks located in the proximal arm of the P5 palindrome and within P1.2, a 92 kb long sequence block located in the P1 palindrome structure of AZFc in Yq11. This large genomic Y region includes deletion of 6 protein encoding Y genes, EIFA1Y, HSFY, PRY, RBMY1, RPS4Y, SMCY. Additionally, one copy of CDY2 and XKRY located in the proximal P5 palindrome and one copy of BPY1, two copies of DAZ located in the P2 palindrome, and one copy of CDY1 located proximal to P1.2 are included within this AZFb microdeletion. It overlaps thus distally along 2.3 Mb with the proximal part of the genomic AZFc deletion. However, AZFb deletions have been also reported with distinct break sites in the proximal and/or distal AZFb breakpoint intervals on the Y chromosome of infertile men. These so called “non-classical” AZFb deletions are associated with variable testicular pathologies, including meiotic arrest, cryptozoospermia, severe oligozoospermia, or oligoasthenoteratozoospermia (OAT syndrome), respectively. This raised the question whether there are any specific length(s) of the AZFb deletion interval along Yq11 required to cause meiotic arrest of the patient’s spermatogenesis, respectively, whether there is any single AZFb Y gene deletion also able to cause this “classical” AZFb testicular pathology? Review of the literature and more cases with “classical” and “non-classical” AZFb deletions analysed in our lab since the last 20 years suggests that the composition of the genomic Y sequence in AZFb is variable in men with distinct Y haplogroups especially in the distal AZFb region overlapping with the proximal AZFc deletion interval and that its extension can be “polymorphic” in the P3 palindrome. That means this AZFb subinterval can be rearranged or deleted also on the Y chromosome of fertile men. Any AZFb deletion observed in infertile men with azoospermia should therefore be confirmed as “de novo” mutation event, i.e., not present on the Y chromosome of the patient’s father or fertile brother before it is considered as causative agent for man’s infertility. Moreover, its molecular length in Yq11 should be comparable to that of the “classical” AZFb deletion, before meiotic arrest is prognosed as the patient’s testicular pathology.

Mutation Rate Variability across Human Y-Chromosome Haplogroups

A common assumption in dating patrilineal events using Y-chromosome sequencing data is that the Y-chromosome mutation rate is invariant across haplogroups. Previous studies revealed interhaplogroup heterogeneity in phylogenetic branch length. Whether this heterogeneity is caused by interhaplogroup mutation rate variation or nongenetic confounders remains unknown. Here, we analyzed whole-genome sequences from cultured cells derived from >1,700 males. We confirmed the presence of branch length heterogeneity. We demonstrate that sex-chromosome mutations that appear within cell lines, which likely occurred somatically or in vitro (and are thus not influenced by nongenetic confounders) are informative for germline mutational processes. Using within-cell-line mutations, we computed a relative Y-chromosome somatic mutation rate, and uncovered substantial variation (up to 83.3%) in this proxy for germline mutation rate among haplogroups. This rate positively correlates with phylogenetic branch length, indicating that interhaplogroup mutation rate variation is a likely cause of branch length heterogeneity.

An atlas of CNV maps in cattle, goat and sheep

Copy number variation (CNV) is the most prevalent type of genetic structural variation that has been recognized as an important source of phenotypic variation in humans, animals and plants. However, the mechanisms underlying the evolution of CNVs and their function in natural or artificial selection remain unknown. Here, we generated CNV region (CNVR) datasets which were diverged or shared among cattle, goat, and sheep, including 886 individuals from 171 diverse populations. Using 9 environmental factors for genome-wide association study (GWAS), we identified a series of candidate CNVRs, including genes relating to immunity, tick resistance, multi-drug resistance, and muscle development. The number of CNVRs shared between species is significantly higher than expected (P<0.00001), and these CNVRs may be more persist than the single nucleotide polymorphisms (SNPs) shared between species. We also identified genomic regions under long-term balancing selection and uncovered the potential diversity of the selected CNVRs close to the important functional genes. This study provides the evidence that balancing selection might be more common in mammals than previously considered, and might play an important role in the daily activities of these ruminant species.

A Role for Chromosomal Microarray Testing in the Workup of Male Infertility

Genetic analysis is a critical component in the male infertility workup. For male infertility due to oligospermia/azoospermia, standard guidelines recommend karyotype and Y-chromosome microdeletion analyses. A karyotype is used to identify structural and numerical chromosome abnormalities, whereas Y-chromosome microdeletions are commonly evaluated by multiplex PCR analysis because of their submicroscopic size. Because these assays often require different Vacutainer tubes to be sent to different laboratories, ordering is prone to errors. In addition, this workflow limits the ability for sequential testing and a comprehensive test result. A potential solution includes performing Y-microdeletion and numerical chromosome analysis-the most common genetic causes of oligospermia/azoospermia-by chromosomal microarray (CMA) and reflexing to karyotype as both assays are often offered in the cytogenetics laboratory. Such analyses can be performed using one sodium heparin Vacutainer tube sample. To determine the effectiveness of CMA for the detection of clinically significant Y-chromosome microdeletions, 21 cases with known Y microdeletions were tested by CytoScan HD platform. CMA studies identified all known Y-chromosome microdeletions, and in 11 cases (52%) identified additional clinically important cytogenetic anomalies, including six cases of 46, XX males, one case of isodicentric Y, two cases of a dicentric Y, and three cases of terminal Yq deletions. These findings demonstrate that this testing strategy would simplify ordering and allow for an integrated interpretation of test results.

Battle of the sexes: contrasting roles of testis-specific protein Y-encoded (TSPY) and TSPX in human oncogenesis

The Y-located testis-specific protein Y-encoded (TSPY) and its X-homologue TSPX originated from the same ancestral gene, but act as a proto-oncogene and a tumor suppressor gene, respectively. TSPY has specialized in male-specific functions, while TSPX has assumed the functions of the ancestral gene. Both TSPY and TSPX harbor a conserved SET/NAP domain, but are divergent at flanking structures. Specifically, TSPX contains a C-terminal acidic domain, absent in TSPY. They possess contrasting properties, in which TSPY and TSPX, respectively, accelerate and arrest cell proliferation, stimulate and inhibit cyclin B-CDK1 phosphorylation activities, have no effect and promote proteosomal degradation of the viral HBx oncoprotein, and exacerbate and repress androgen receptor (AR) and constitutively active AR variant, such as AR-V7, gene transactivation. The inhibitory domain has been mapped to the carboxyl acidic domain in TSPX, truncation of which results in an abbreviated TSPX exerting positive actions as TSPY. Transposition of the acidic domain to the C-terminus of TSPY results in an inhibitory protein as intact TSPX. Hence, genomic mutations/aberrant splicing events could generate TSPX proteins with truncated acidic domain and oncogenic properties as those for TSPY. Further, TSPY is upregulated by AR and AR-V7 in ligand-dependent and ligand-independent manners, respectively, suggesting the existence of a positive feedback loop between a Y-located proto-oncogene and male sex hormone/receptors, thereby amplifying the respective male oncogenic actions in human cancers and diseases. TSPX counteracts such positive feedback loop. Hence, TSPY and TSPX are homologues on the sex chromosomes that function at the two extremes of the human oncogenic spectrum.

Role of genetics and epigenetics in male infertility

Male infertility is a complex condition with a strong genetic and epigenetic background. This review discusses the importance of genetic and epigenetic factors in the pathophysiology of male infertility. The interplay between thousands of genes, the epigenetic control of gene expression, and environmental and lifestyle factors, which influence genetic and epigenetic variants, determines the resulting male infertility phenotype. Currently, karyotyping, Y-chromosome microdeletion screening and CFTR gene mutation tests are routinely performed to investigate a possible genetic aetiology in patients with azoospermia and severe oligozoospermia. However, current testing is limited in its ability to identify a variety of genetic and epigenetic conditions that might be implicated in both idiopathic and unexplained infertility. Several epimutations of imprinting genes and developmental genes have been postulated to be candidate markers for male infertility. As such, development of novel diagnostic panels is essential to change the current landscape with regard to prevention, diagnosis and management. Understanding the underlying genetic mechanisms related to the pathophysiology of male infertility, and the impact of environmental exposures and lifestyle factors on gene expression might aid clinicians in developing individualised treatment strategies.

The Y-chromosome F haplogroup contributes to the development of Barrett's esophagus-associated esophageal adenocarcinoma in a white male population

Barrett's esophagus (BE) is a metaplastic condition of the distal esophagus, resulting from longstanding gastroesophageal reflux disease (GERD). BE predisposes for the highly malignant esophageal adenocarcinoma (EAC). Both BE and EAC have the highest frequencies in white males. Only a subset of patients with GERD develop BE, while <0.5% of BE will progress to EAC. Therefore, it is most likely that the development of BE and EAC is associated with underlying genetic factors. We hypothesized that in white males, Y-chromosomal haplogroups are associated with BE and EAC. To investigate this we conducted a multicenter study studying the frequencies of the Y-chromosomal haplogroups in GERD, BE, and EAC patients. We used genomic analysis by polymerase chain reaction and restriction fragment length polymorphism to determine the frequency of six Y-chromosomal haplogroups (DE, F(xJ,xK), K(xP), J, P(xR1a), and R1a) between GERD, BE, and EAC in a cohort of 1,365 white males, including 612 GERD, 753 BE patients, while 178 of the BE patients also had BE-associated EAC. Univariate logistic regression analysis was used to compare the outcomes. In this study, we found the R1a (6% vs. 9%, P = 0.04) and K (3% vs. 6%, P = 0.035) to be significantly underrepresented in BE patients as compared to GERD patients with an odds ratio (OR) of 0.63 (95% CI 0.42-0.95, P = 0.03) and of 0.56 (95% CI 0.33-0.96, P = 0.03), respectively, while the K haplogroup was protective against EAC (OR 0.30; 95% CI 0.07-0.86, P = 0.05). A significant overrepresentation of the F haplogroup was found in EAC compared to BE and GERD patients (34% vs. 27% and 23%, respectively). The F haplogroup was found to be a risk factor for EAC with an OR of 1.5 (95% CI 1.03-2.19, P = 0.03). We identified the R1a and K haplogroups as protective factors against development of BE. These haplogroups have low frequencies in white male populations. Of importance is that we could link the presence of the predominantly occurring F haplogroup in white males to EAC. It is possible that this F haplogroup is associated to genetic variants that predispose for the EAC development. In future, the haplogroups could be applied to improve stratification of BE and GERD patients with increased risk to develop BE and/or EAC.

Genetic Landscape of Nonobstructive Azoospermia and New Perspectives for the Clinic

Nonobstructive azoospermia (NOA) represents the most severe expression of male infertility, involving around 1% of the male population and 10% of infertile men. This condition is characterised by the inability of the testis to produce sperm cells, and it is considered to have an important genetic component. During the last two decades, different genetic anomalies, including microdeletions of the Y chromosome, karyotype defects, and missense mutations in genes involved in the reproductive function, have been described as the primary cause of NOA in many infertile men. However, these alterations only explain around 25% of azoospermic cases, with the remaining patients showing an idiopathic origin. Recent studies clearly suggest that the so-called idiopathic NOA has a complex aetiology with a polygenic inheritance, which may alter the spermatogenic process. Although we are far from a complete understanding of the molecular mechanisms underlying NOA, the use of the new technologies for genetic analysis has enabled a considerable increase in knowledge during the last years. In this review, we will provide a comprehensive and updated overview of the genetic basis of NOA, with a special focus on the possible application of the recent insights in clinical practice.

AZF deletions in Indian populations: original study and meta-analyses

PURPOSE

To identify the frequency of Y chromosome microdeletions in Indian populations and to quantitatively estimate the significance of association between these deletions and male infertility.

METHODS

A total of 379 infertile males (302 azoospermic and 77 oligozoospermic infertile males) and 265 normozoospermic fertile males were evaluated for Y chromosome microdeletions (YCD) using PCR amplification and gel electrophoresis. Meta-analyses were performed on AZFa (2079 cases and 1217 controls), AZFb (2212 cases and 1267 controls), AZFc (4131 cases and 2008 controls), and AZFb+c (1573 cases and 942 controls) deletions data to quantitatively estimate the significance of association between these deletions and male infertility in Indian populations.

RESULTS

The results revealed that out of 379 infertile azoospermic and oligozoospermic males, 38 (10.02%) had AZF deletions. No deletion was found in control samples. The highest percentage of deletions was observed in the AZFc region, followed by AZFa and AZFb. Qualitative analysis showed that AZF deletions were present in 0.59 to 32.62% (average 13.48%) of infertile cases in Indian populations. Meta-analysis revealed a significant association of AZFa (OR = 6.74, p value = 0.001), AZFb (OR = 4.694, p value = 0.004), AZFc (OR = 13.575, p value = 0.000), and AZFb+c (OR = 5.946, p value = 0.018) deletions with male infertility.

CONCLUSION

AZF deletions were seen in 10.02% of azoospermic and oligozoospermic cases with the highest frequency of AZFc deletions. Pooled analysis for all studies showed deletion frequency from 0.59 to 32.62% (average = 13.48%). Meta-analysis showed significant association of AZFa, AZFb, and AZFb+c deletions with male infertility. Analysis of Y chromosome microdeletions should be reckoned as an essential testing for diagnostic and therapeutic purposes.

The Role of Number of Copies, Structure, Behavior and Copy Number Variations (CNV) of the Y Chromosome in Male Infertility

The World Health Organization (WHO) defines infertility as the inability of a sexually active, non-contracepting couple to achieve spontaneous pregnancy within one year. Statistics show that the two sexes are equally at risk. Several causes may be responsible for male infertility; however, in 30–40% of cases a diagnosis of idiopathic male infertility is made in men with normal urogenital anatomy, no history of familial fertility-related diseases and a normal panel of values as for endocrine, genetic and biochemical markers. Idiopathic male infertility may be the result of gene/environment interactions, genetic and epigenetic abnormalities. Numerical and structural anomalies of the Y chromosome represent a minor yet significant proportion and are the topic discussed in this review. We searched the PubMed database and major search engines for reports about Y-linked male infertility. We present cases of Y-linked male infertility in terms of (i) anomalies of the Y chromosome structure/number; (ii) Y chromosome misbehavior in a normal genetic background; (iii) Y chromosome copy number variations (CNVs). We discuss possible explanations of male infertility caused by mutations, lower or higher number of copies of otherwise wild type, Y-linked sequences. Despite Y chromosome structural anomalies are not a major cause of male infertility, in case of negative results and of normal DNA sequencing of the ascertained genes causing infertility and mapping on this chromosome, we recommend an analysis of the karyotype integrity in all cases of idiopathic fertility impairment, with an emphasis on the structure and number of this chromosome.

Evolutionary and functional impact of common polymorphic inversions in the human genome

Inversions are one type of structural variants linked to phenotypic differences and adaptation in multiple organisms. However, there is still very little information about polymorphic inversions in the human genome due to the difficulty of their detection. Here, we develop a new high-throughput genotyping method based on probe hybridization and amplification, and we perform a complete study of 45 common human inversions of 0.1–415 kb. Most inversions promoted by homologous recombination occur recurrently in humans and great apes and they are not tagged by SNPs. Furthermore, there is an enrichment of inversions showing signatures of positive or balancing selection, diverse functional effects, such as gene disruption and gene-expression changes, or association with phenotypic traits. Therefore, our results indicate that the genome is more dynamic than previously thought and that human inversions have important functional and evolutionary consequences, making possible to determine for the first time their contribution to complex traits.

Inversions are a little-studied type of genomic variation that could contribute to phenotypic traits. Here the authors characterize 45 common polymorphic inversions in human populations and investigate their evolutionary and functional impact.

High frequencies of Non Allelic Homologous Recombination (NAHR) events at the AZF loci and male infertility risk in Indian men

Deletions in the AZoospermia Factor (AZF) regions (spermatogenesis loci) on the human Y chromosome are reported as one of the most common causes of severe testiculopathy and spermatogenic defects leading to male infertility, yet not much data is available for Indian infertile men. Therefore, we screened for AZF region deletions in 973 infertile men consisting of 771 azoospermia, 105 oligozoospermia and 97 oligoteratozoospermia cases, along with 587 fertile normozoospermic men. The deletion screening was carried out using AZF-specific markers: STSs (Sequence Tagged Sites), SNVs (Single Nucleotide Variations), PCR-RFLP (Polymerase Chain Reaction - Restriction Fragment Length Polymorphism) analysis of STS amplicons, DNA sequencing and Southern hybridization techniques. Our study revealed deletion events in a total of 29.4% of infertile Indian men. Of these, non-allelic homologous recombination (NAHR) events accounted for 25.8%, which included 3.5% AZFb deletions, 2.3% AZFbc deletions, 6.9% complete AZFc deletions, and 13.1% partial AZFc deletions. We observed 3.2% AZFa deletions and a rare long AZFabc region deletion in 0.5% azoospermic men. This study illustrates how the ethnicity, endogamy and long-time geographical isolation of Indian populations might have played a major role in the high frequencies of deletion events.

The evolutionary history of grey wolf Y chromosomes

Analyses of Y chromosome haplotypes uniquely provide a paternal picture of evolutionary histories and offer a very useful contrast to studies based on maternally inherited mitochondrial DNA (mtDNA). Here we used a bioinformatic approach based on comparison of male and female sequence coverage to identify 4.7 Mb from the grey wolf (Canis lupis) Y chromosome, probably representing most of the male-specific, nonampliconic sequence from the euchromatic part of the chromosome. We characterized this sequence and then identified ≈1,500 Y-linked single nucleotide polymorphisms in a sample of 145 resequenced male wolves, including 75 Finnish wolf genomes newly sequenced in this study, and in 24 dogs and eight other canids. We found 53 Y chromosome haplotypes, of which 26 were seen in grey wolves, that clustered in four major haplogroups. All four haplogroups were represented in samples of Finnish wolves, showing that haplogroup lineages were not partitioned on a continental scale. However, regional population structure was indicated because individual haplotypes were never shared between geographically distant areas, and genetically similar haplotypes were only found within the same geographical region. The deepest split between grey wolf haplogroups was estimated to have occurred 125,000 years ago, which is considerably older than recent estimates of the time of divergence of wolf populations. The distribution of dogs in a phylogenetic tree of Y chromosome haplotypes supports multiple domestication events, or wolf paternal introgression, starting 29,000 years ago. We also addressed the disputed origin of a recently founded population of Scandinavian wolves and observed that founding as well as most recent immigrant haplotypes were present in the neighbouring Finnish population, but not in sequenced wolves from elsewhere in the world, or in dogs.

Genetic Instability and Chromatin Remodeling in Spermatids

The near complete replacement of somatic chromatin in spermatids is, perhaps, the most striking nuclear event known to the eukaryotic domain. The process is far from being fully understood, but research has nevertheless unraveled its complexity as an expression of histone variants and post-translational modifications that must be finely orchestrated to promote the DNA topological change and compaction provided by the deposition of protamines. That this major transition may not be genetically inert came from early observations that transient DNA strand breaks were detected in situ at chromatin remodeling steps. The potential for genetic instability was later emphasized by our demonstration that a significant number of DNA double-strand breaks (DSBs) are formed and then repaired in the haploid context of spermatids. The detection of DNA breaks by 3'OH end labeling in the whole population of spermatids suggests that a reversible enzymatic process is involved, which differs from canonical apoptosis. We have set the stage for a better characterization of the genetic impact of this transition by showing that post-meiotic DNA fragmentation is conserved from human to yeast, and by providing tools for the initial mapping of the genome-wide DSB distribution in the mouse model. Hence, the molecular mechanism of post-meiotic DSB formation and repair in spermatids may prove to be a significant component of the well-known male mutation bias. Based on our recent observations and a survey of the literature, we propose that the chromatin remodeling in spermatids offers a proper context for the induction of de novo polymorphism and structural variations that can be transmitted to the next generation.

Selection Has Countered High Mutability to Preserve the Ancestral Copy Number of Y Chromosome Amplicons in Diverse Human Lineages

Amplicons-large, highly identical segmental duplications-are a prominent feature of mammalian Y chromosomes. Although they encode genes essential for fertility, these amplicons differ vastly between species, and little is known about the selective constraints acting on them. Here, we develop computational tools to detect amplicon copy number with unprecedented accuracy from high-throughput sequencing data. We find that one-sixth (16.9%) of 1,216 males from the 1000 Genomes Project have at least one deleted or duplicated amplicon. However, each amplicon's reference copy number is scrupulously maintained among divergent branches of the Y chromosome phylogeny, including the ancient branch A00, indicating that the reference copy number is ancestral to all modern human Y chromosomes. Using phylogenetic analyses and simulations, we demonstrate that this pattern of variation is incompatible with neutral evolution and instead displays hallmarks of mutation-selection balance. We also observe cases of amplicon rescue, in which deleted amplicons are restored through subsequent duplications. These results indicate that, contrary to the lack of constraint suggested by the differences between species, natural selection has suppressed amplicon copy number variation in diverse human lineages.

Sequence and Structural Diversity of Mouse Y Chromosomes

Over the 180 My since their origin, the sex chromosomes of mammals have evolved a gene repertoire highly specialized for function in the male germline. The mouse Y chromosome is unique among mammalian Y chromosomes characterized to date in that it is large, gene-rich and euchromatic. Yet, little is known about its diversity in natural populations. Here, we take advantage of published whole-genome sequencing data to survey the diversity of sequence and copy number of sex-linked genes in three subspecies of house mice. Copy number of genes on the repetitive long arm of both sex chromosomes is highly variable, but sequence diversity in nonrepetitive regions is decreased relative to expectations based on autosomes. We use simulations and theory to show that this reduction in sex-linked diversity is incompatible with neutral demographic processes alone, but is consistent with recent positive selection on genes active during spermatogenesis. Our results support the hypothesis that the mouse sex chromosomes are engaged in ongoing intragenomic conflict.

Dynamic Copy Number Evolution of X- and Y-Linked Ampliconic Genes in Human Populations

Ampliconic genes are multicopy genes often located on sex chromosomes and enriched for testis-expressed genes. Here, Lucotte et al. developed new bioinformatic approaches to investigate the ampliconic gene copy number and their coding...

Ampliconic genes are multicopy, with the majority found on sex chromosomes and enriched for testis-expressed genes. While ampliconic genes have been associated with the emergence of hybrid incompatibilities, we know little about their copy number distribution and their turnover in human populations. Here, we explore the evolution of human X- and Y-linked ampliconic genes by investigating copy number variation (CNV) and coding variation between populations using the Simons Genome Diversity Project. We develop a method to assess CNVs using the read depth on modified X and Y chromosome targets containing only one repetition of each ampliconic gene. Our results reveal extensive standing variation in copy number both within and between human populations for several ampliconic genes. For the Y chromosome, we can infer multiple independent amplifications and losses of these gene copies even within closely related Y haplogroups, that diversified < 50,000 years ago. Moreover, X- and Y-linked ampliconic genes seem to have a faster amplification dynamic than autosomal multicopy genes. Looking at expression data from another study, we also find that X- and Y-linked ampliconic genes with extensive CNV are significantly more expressed than genes with no CNV during meiotic sex chromosome inactivation (for both X and Y) and postmeiotic sex chromosome repression (for the Y chromosome only). While we cannot rule out that the XY-linked ampliconic genes are evolving neutrally, this study gives insights into the distribution of copy number within human populations and demonstrates an extremely fast turnover in copy number of these regions.

TSPY1 suppresses USP7-mediated p53 function and promotes spermatogonial proliferation

Testis-specific protein Y-linked 1 (TSPY1) is expressed predominantly in adult human spermatogonia and functions in the process of spermatogenesis; however, our understanding of the underlying mechanism is limited. Here we observed that TSPY1, as an interacting partner of TSPY-like 5 (TSPYL5), enhanced the competitive binding of TSPYL5 to ubiquitin-specific peptidase 7 (USP7) in conjunction with p53. This activity, together with its promotion of TSPYL5 expression by acting as a transcription factor, resulted in increased p53 ubiquitylation. Moreover, TSPY1 could decrease the p53 level by inducing the degradation of ubiquitinated USP7. We demonstrated that the promotion of p53 degradation by TSPY1 influenced the activity of p53 target molecules (CDK1, p21, and BAX) to expedite the G2/M phase transition and decrease cell apoptosis, accelerating cell proliferation. Taken together, the observations reveal the significance of TSPY1 as a suppressor of USP7-mediated p53 function in inhibiting p53-dependent cell proliferation arrest. By simulating TSPY1 function in Tspy1-deficient spermatogonia derived from mouse testes, we found that TSPY1 could promote spermatogonial proliferation by decreasing the Usp7-modulated p53 level. The findings suggest an additional mechanism underlying the regulation of spermatogonial p53 function, indicating the significance of TSPY1 in germline homeostasis maintenance and the potential of TSPY1 in regulating human spermatogonial proliferation via the USP7-mediated p53 signaling pathway.

High Levels of Copy Number Variation of Ampliconic Genes across Major Human Y Haplogroups

Because of its highly repetitive nature, the human male-specific Y chromosome remains understudied. It is important to investigate variation on the Y chromosome to understand its evolution and contribution to phenotypic variation, including infertility. Approximately 20% of the human Y chromosome consists of ampliconic regions which include nine multi-copy gene families. These gene families are expressed exclusively in testes and usually implicated in spermatogenesis. Here, to gain a better understanding of the role of the Y chromosome in human evolution and in determining sexually dimorphic traits, we studied ampliconic gene copy number variation in 100 males representing ten major Y haplogroups world-wide. Copy number was estimated with droplet digital PCR. In contrast to low nucleotide diversity observed on the Y in previous studies, here we show that ampliconic gene copy number diversity is very high. A total of 98 copy-number-based haplotypes were observed among 100 individuals, and haplotypes were sometimes shared by males from very different haplogroups, suggesting homoplasies. The resulting haplotypes did not cluster according to major Y haplogroups. Overall, only two gene families (RBMY and TSPY) showed significant differences in copy number among major Y haplogroups, and the haplogroup of a male could not be predicted based on his ampliconic gene copy numbers. Finally, we did not find significant correlations either between copy number variation and individual's height, or between the former and facial masculinity/femininity. Our results suggest rapid evolution of ampliconic gene copy numbers on the human Y, and we discuss its causes.

Cost-effective high-throughput single-haplotype iterative mapping and sequencing for complex genomic structures

The reference sequences of structurally complex regions can be obtained only through highly accurate clone-based approaches. We and others have successfully used single-haplotype iterative mapping and sequencing (SHIMS) 1.0 to assemble structurally complex regions across the sex chromosomes of several vertebrate species and to allow for targeted improvements to the reference sequences of human autosomes. However, SHIMS 1.0 is expensive and time consuming, requiring resources that only a genome center can provide. Here we introduce SHIMS 2.0, an improved SHIMS protocol that allows even a small laboratory to generate high-quality reference sequence from complex genomic regions. Using a streamlined and parallelized library-preparation protocol, and taking advantage of inexpensive high-throughput short-read-sequencing technologies, a small laboratory with both molecular biology and bioinformatics experience can sequence and assemble 192 large-insert bacterial artificial chromosome (BAC) or fosmid clones in 1 week. In SHIMS 2.0, in contrast to other pooling strategies, each clone is sequenced with a unique barcode, thus enabling clones containing nearly identical sequences to be multiplexed in a single sequencing run and assembled separately. Relative to SHIMS 1.0, SHIMS 2.0 decreases the required cost and time by two orders of magnitude while preserving high sequencing accuracy.

Copy number variation arising from gene conversion on the human Y chromosome

We describe the variation in copy number of a ~ 10 kb region overlapping the long intergenic noncoding RNA (lincRNA) gene, TTTY22, within the IR3 inverted repeat on the short arm of the human Y chromosome, leading to individuals with 0-3 copies of this region in the general population. Variation of this CNV is common, with 266 individuals having 0 copies, 943 (including the reference sequence) having 1, 23 having 2 copies, and two having 3 copies, and was validated by breakpoint PCR, fibre-FISH, and 10× Genomics Chromium linked-read sequencing in subsets of 1234 individuals from the 1000 Genomes Project. Mapping the changes in copy number to the phylogeny of these Y chromosomes previously established by the Project identified at least 20 mutational events, and investigation of flanking paralogous sequence variants showed that the mutations involved flanking sequences in 18 of these, and could extend over > 30 kb of DNA. While either gene conversion or double crossover between misaligned sister chromatids could formally explain the 0-2 copy events, gene conversion is the more likely mechanism, and these events include the longest non-allelic gene conversion reported thus far. Chromosomes with three copies of this CNV have arisen just once in our data set via another mechanism: duplication of 420 kb that places the third copy 230 kb proximal to the existing proximal copy. Our results establish gene conversion as a previously under-appreciated mechanism of generating copy number changes in humans and reveal the exceptionally large size of the conversion events that can occur.

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.

Analysis of 62 hybrid assembled human Y chromosomes exposes rapid structural changes and high rates of gene conversion

The human Y-chromosome does not recombine across its male-specific part and is therefore an excellent marker of human migrations. It also plays an important role in male fertility. However, its evolution is difficult to fully understand because of repetitive sequences, inverted repeats and the potentially large role of gene conversion. Here we perform an evolutionary analysis of 62 Y-chromosomes of Danish descent sequenced using a wide range of library insert sizes and high coverage, thus allowing large regions of these chromosomes to be well assembled. These include 17 father-son pairs, which we use to validate variation calling. Using a recent method that can integrate variants based on both mapping and de novo assembly, we genotype 10898 SNVs and 2903 indels (max length of 27241 bp) in our sample and show by father-son concordance and experimental validation that the non-recurrent SNP and indel variation on the Y chromosome tree is called very accurately. This includes variation called in a 0.9 Mb centromeric heterochromatic region, which is by far the most variable in the Y chromosome. Among the variation is also longer sequence-stretches not present in the reference genome but shared with the chimpanzee Y chromosome. We analyzed 2.7 Mb of large inverted repeats (palindromes) for variation patterns among the two palindrome arms and identified 603 mutation and 416 gene conversions events. We find clear evidence for GC-biased gene conversion in the palindromes (and a balancing AT mutation bias), but irrespective of this, also a strong bias towards gene conversion towards the ancestral state, suggesting that palindromic gene conversion may alleviate Muller’s ratchet. Finally, we also find a large number of large-scale gene duplications and deletions in the palindromic regions (at least 24) and find that such events can consist of complex combinations of simultaneous insertions and deletions of long stretches of the Y chromosome.

The Y chromosome is extraordinary in many respects; it is non-recombining along most of its length, it carries many testis-expressed genes that are often found in palindromes and thus in several copies, and it is generally highly repetitive with very few unique genes. Its evolutionary process is not well understood in general because short-read mapping in such complex sequence is difficult. We combine de novo assembly and mapping to investigate evolution in more than 60% of the length of 62 Y chromosomes of Danish descent. We find that Y chromosome evolution is very dynamic even among the set of closely related Y chromosomes in Denmark with many cases of complex duplications and deletions of large regions including whole genes, clear evidence of GC-biased gene conversion in the palindromes and a tendency for gene conversion to revert mutations to their ancestral state.

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.

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.

Human Y chromosome copy number variation in the next generation sequencing era and beyond

The human Y chromosome provides a fertile ground for structural rearrangements owing to its haploidy and high content of repeated sequences. The methodologies used for copy number variation (CNV) studies have developed over the years. Low-throughput techniques based on direct observation of rearrangements were developed early on, and are still used, often to complement array-based or sequencing approaches which have limited power in regions with high repeat content and specifically in the presence of long, identical repeats, such as those found in human sex chromosomes. Some specific rearrangements have been investigated for decades; because of their effects on fertility, or their outstanding evolutionary features, the interest in these has not diminished. However, following the flourishing of large-scale genomics, several studies have investigated CNVs across the whole chromosome. These studies sometimes employ data generated within large genomic projects such as the DDD study or the 1000 Genomes Project, and often survey large samples of healthy individuals without any prior selection. Novel technologies based on sequencing long molecules and combinations of technologies, promise to stimulate the study of Y-CNVs in the immediate future.

The Y chromosomes of the great apes

The great apes (orangutans, gorillas, chimpanzees, bonobos and humans) descended from a common ancestor around 13 million years ago, and since then their sex chromosomes have followed very different evolutionary paths. While great-ape X chromosomes are highly conserved, their Y chromosomes, reflecting the general lability and degeneration of this male-specific part of the genome since its early mammalian origin, have evolved rapidly both between and within species. Understanding great-ape Y chromosome structure, gene content and diversity would provide a valuable evolutionary context for the human Y, and would also illuminate sex-biased behaviours, and the effects of the evolutionary pressures exerted by different mating strategies on this male-specific part of the genome. High-quality Y-chromosome sequences are available for human and chimpanzee (and low-quality for gorilla). The chromosomes differ in size, sequence organisation and content, and while retaining a relatively stable set of ancestral single-copy genes, show considerable variation in content and copy number of ampliconic multi-copy genes. Studies of Y-chromosome diversity in other great apes are relatively undeveloped compared to those in humans, but have nevertheless provided insights into speciation, dispersal, and mating patterns. Future studies, including data from larger sample sizes of wild-born and geographically well-defined individuals, and full Y-chromosome sequences from bonobos, gorillas and orangutans, promise to further our understanding of population histories, male-biased behaviours, mutation processes, and the functions of Y-chromosomal genes.

Azoospermia Factor C Subregion of the Y Chromosome

The azoospermia factor (AZF) region on the Y chromosome consists of genes required for spermatogenesis. Among the three subregions, the AZFc subregion located at the distal portion of AZF is the driver for genetic variation in Y chromosome. The candidate gene of AZFc is known as deleted in azoospermia gene, which is studied with interest because it is involved in germ cell development and most frequently deleted genes leading to oligozoospermia and azoospermia. Recently, two partial deletions in AZFc gr/gr and b2/b3 are characterized at the molecular level which showed homologous recombination between amplicons, affecting spermatogenesis process. There are novel methods and commercially available kits for accurate screening and characterization of microdeletions. It is important to detect the AZFc microdeletions through genetic screening and counseling those infertile men who planned to avail assisted reproduction techniques such as undergoing intracytoplasmic sperm injection or in vitro fertilization.