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

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.

Assembly of 43 human Y chromosomes reveals extensive complexity and variation

The prevalence of highly repetitive sequences within the human Y chromosome has prevented its complete assembly to date1 and led to its systematic omission from genomic analyses. Here we present de novo assemblies of 43 Y chromosomes spanning 182,900 years of human evolution and report considerable diversity in size and structure. Half of the male-specific euchromatic region is subject to large inversions with a greater than twofold higher recurrence rate compared with all other chromosomes2. Ampliconic sequences associated with these inversions show differing mutation rates that are sequence context dependent, and some ampliconic genes exhibit evidence for concerted evolution with the acquisition and purging of lineage-specific pseudogenes. The largest heterochromatic region in the human genome, Yq12, is composed of alternating repeat arrays that show extensive variation in the number, size and distribution, but retain a 1:1 copy-number ratio. Finally, our data suggest that the boundary between the recombining pseudoautosomal region 1 and the non-recombining portions of the X and Y chromosomes lies 500 kb away from the currently established1 boundary. The availability of fully sequence-resolved Y chromosomes from multiple individuals provides a unique opportunity for identifying new associations of traits with specific Y-chromosomal variants and garnering insights into the evolution and function of complex regions of the human genome.

The complete sequence of a human Y chromosome

The human Y chromosome has been notoriously difficult to sequence and assemble because of its complex repeat structure that includes long palindromes, tandem repeats and segmental duplications1-3. As a result, more than half of the Y chromosome is missing from the GRCh38 reference sequence and it remains the last human chromosome to be finished4,5. Here, the Telomere-to-Telomere (T2T) consortium presents the complete 62,460,029-base-pair sequence of a human Y chromosome from the HG002 genome (T2T-Y) that corrects multiple errors in GRCh38-Y and adds over 30 million base pairs of sequence to the reference, showing the complete ampliconic structures of gene families TSPY, DAZ and RBMY; 41 additional protein-coding genes, mostly from the TSPY family; and an alternating pattern of human satellite 1 and 3 blocks in the heterochromatic Yq12 region. We have combined T2T-Y with a previous assembly of the CHM13 genome4 and mapped available population variation, clinical variants and functional genomics data to produce a complete and comprehensive reference sequence for all 24 human chromosomes.

High incidence of AZF duplications in clan-structured Iranian populations detected through Y chromosome sequencing read depth analysis

The ampliconic region of the human Y chromosome consists of large duplicated sequences that can undergo non-allelic homologous recombination (NAHR), resulting in structural rearrangements that may cause infertility, especially when they occur in the azoospermia factor b/c (AZFb/c) region. Although AZF duplications have long been neglected due to the technical limitations of STS-based studies that focused mainly on deletions, recent next generation sequencing (NGS) technologies provided evidence for their importance in fertility. In this study, a NGS read depth approach was used to detect AZFb/c rearrangements in 87 Iranians from different ethnic groups. The duplication frequency in Iran proved to be twice as high as in the "1000 Genomes" dataset. Interestingly, most duplications were found in patrilineal ethnic groups, possibly as a consequence of their lower male effective population size which can counteract negative selection. Moreover, we found a large 8.0 Mb duplication, resulting in a fourfold increase in the copy number of AZFc genes, which to our knowledge is the largest duplication ever reported in this region. Overall, our results suggest that it is important to consider not only AZF deletions but also duplications to investigate the causes of male infertility, especially in patrilineal clan-based populations.

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.

Genomic instability in individuals with sex determination defects and germ cell cancer

The ability to transmit genetic information through generations depends on the preservation of genome integrity. Genetic abnormalities affect cell differentiation, causing tissue specification defects and cancer. We addressed genomic instability in individuals with Differences of Sex Development (DSD), characterized by gonadal dysgenesis, infertility, high susceptibility for different types of cancer, especially Germ Cell Tumors (GCT), and in men with testicular GCTs. Whole proteome analysis of leukocytes, supported by specific gene expression assessment, and dysgenic gonads characterization, uncovered DNA damage phenotypes with altered innate immune response and autophagy. Further examination of DNA damage response revealed a reliance on deltaTP53, which was compromised by mutations in the transactivation domain in DSD-individuals with GCT. Accordingly, drug-induced rescue of DNA damage was achieved by autophagy inhibition but not by stabilization of TP53 in DSD-individuals' blood in vitro. This study elucidates possibilities for prophylactic treatments of DSD-individuals, as well as new diagnostic approaches of GCT.

Disclosing complex mutational dynamics at a Y chromosome palindrome evolving through intra- and inter-chromosomal gene conversion

The human MSY ampliconic region is mainly composed of large duplicated sequences that are organized in eight palindromes (termed P1-P8), and may undergo arm-to-arm gene conversion. Although the importance of these elements is widely recognized, their evolutionary dynamics are still nuanced. Here, we focused on the P8 palindrome, which shows a complex evolutionary history, being involved in intra- and inter-chromosomal gene conversion. To disclose its evolutionary complexity, we performed a high-depth (50×) targeted next-generation sequencing of this element in 157 subjects belonging to the most divergent lineages of the Y chromosome tree. We found a total of 72 polymorphic paralogous sequence variants that have been exploited to identify 41 Y-Y gene conversion events that occurred during recent human history. Through our analysis, we were able to categorize P8 arms into three portions, whose molecular diversity was modelled by different evolutionary forces. Notably, the outer region of the palindrome is not involved in any gene conversion event and evolves exclusively through the action of mutational pressure. The inner region is affected by Y-Y gene conversion occurring at a rate of 1.52 × 10-5 conversions/base/year, with no bias towards the retention of the ancestral state of the sequence. In this portion, GC-biased gene conversion is counterbalanced by a mutational bias towards AT bases. Finally, the middle region of the arms, in addition to intra-chromosomal gene conversion, is involved in X-to-Y gene conversion (at a rate of 6.013 × 10-8 conversions/base/year) thus being a major force in the evolution of the VCY/VCX gene family.

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.

New insights into the evolution of human Y chromosome palindromes through mutation and gene conversion

About one-quarter of the euchromatic portion of the male-specific region of the human Y chromosome consists of large duplicated sequences that are organized in eight palindromes (termed P1-P8), which undergo arm-to arm gene conversion, a proposed mechanism for maintaining their sequence integrity. Although the relevance of gene conversion in the evolution of palindromic sequences has been profoundly recognized, the dynamic of this mechanism is still nuanced. To shed light into the evolution of these genomic elements, we performed a high-depth (50×) targeted next-generation sequencing of the palindrome P6 in 157 subjects belonging to the most divergent evolutionary lineages of the Y chromosome. We found 118 new paralogous sequence variants, which were placed into the context of a robust Y chromosome phylogeny based on 7240 SNPs of the X-degenerate region. We mapped along the phylogeny 80 gene conversion events that shaped the diversity of P6 arms during recent human history. In contrast to previous studies, we demonstrated that arm-to-arm gene conversion, which occurs at a rate of 6.01 × 10 -6 conversions/base/year, is not biased toward the retention of the ancestral state of sequences. We also found a significantly lower mutation rate of the arms (6.18 × 10-10 mutations/base/year) compared with the spacer (9.16 × 10-10 mutations/base/year), a finding that may explain the observed higher inter-species conservation of arms, without invoking any bias of conversion. Finally, by formally testing the mutation/conversion balance in P6, we found that the arms of this palindrome reached a steady-state equilibrium between mutation and gene conversion.

Large palindromes on the primate X Chromosome are preserved by natural selection

Mammalian sex chromosomes carry large palindromes that harbor protein-coding gene families with testis-biased expression. However, there are few known examples of sex-chromosome palindromes conserved between species. We identified 26 palindromes on the human X Chromosome, constituting more than 2% of its sequence, and characterized orthologous palindromes in the chimpanzee and the rhesus macaque using a clone-based sequencing approach that incorporates full-length nanopore reads. Many of these palindromes are missing or misassembled in the current reference assemblies of these species' genomes. We find that 12 human X palindromes have been conserved for at least 25 million years, with orthologs in both chimpanzee and rhesus macaque. Insertions and deletions between species are significantly depleted within the X palindromes' protein-coding genes compared to their noncoding sequence, demonstrating that natural selection has preserved these gene families. The spacers that separate the left and right arms of palindromes are a site of localized structural instability, with seven of 12 conserved palindromes showing no spacer orthology between human and rhesus macaque. Analysis of the 1000 Genomes Project data set revealed that human X-palindrome spacers are enriched for deletions relative to arms and flanking sequence, including a common spacer deletion that affects 13% of human X Chromosomes. This work reveals an abundance of conserved palindromes on primate X Chromosomes and suggests that protein-coding gene families in palindromes (most of which remain poorly characterized) promote X-palindrome survival in the face of ongoing structural instability.

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.

Improved Models of Coalescence Ages of Y-DNA Haplogroups

Databases of commercial DNA-testing companies now contain more customers with sequenced DNA than any completed academic study, leading to growing interest from academic and forensic entities. An important result for both these entities and the test takers themselves is how closely two individuals are related in time, as calculated through one or more molecular clocks. For Y-DNA, existing interpretations of these clocks are insufficiently accurate to usefully measure relatedness in historic times. In this article, I update the methods used to calculate coalescence ages (times to most-recent common ancestor, or TMRCAs) using a new, probabilistic statistical model that includes Y-SNP, Y-STR and ancilliary historical data, and provide examples of its use.

Understanding the Early Evolutionary Stages of a Tandem Drosophilamelanogaster-Specific Gene Family: A Structural and Functional Population Study

Gene families underlie genetic innovation and phenotypic diversification. However, our understanding of the early genomic and functional evolution of tandemly arranged gene families remains incomplete as paralog sequence similarity hinders their accurate characterization. The Drosophila melanogaster-specific gene family Sdic is tandemly repeated and impacts sperm competition. We scrutinized Sdic in 20 geographically diverse populations using reference-quality genome assemblies, read-depth methodologies, and qPCR, finding that ∼90% of the individuals harbor 3-7 copies as well as evidence of population differentiation. In strains with reliable gene annotations, copy number variation (CNV) and differential transposable element insertions distinguish one structurally distinct version of the Sdic region per strain. All 31 annotated copies featured protein-coding potential and, based on the protein variant encoded, were categorized into 13 paratypes differing in their 3' ends, with 3-5 paratypes coexisting in any strain examined. Despite widespread gene conversion, the only copy present in all strains has functionally diverged at both coding and regulatory levels under positive selection. Contrary to artificial tandem duplications of the Sdic region that resulted in increased male expression, CNV in cosmopolitan strains did not correlate with expression levels, likely as a result of differential genome modifier composition. Duplicating the region did not enhance sperm competitiveness, suggesting a fitness cost at high expression levels or a plateau effect. Beyond facilitating a minimally optimal expression level, Sdic CNV acts as a catalyst of protein and regulatory diversity, showcasing a possible evolutionary path recently formed tandem multigene families can follow toward long-term consolidation in eukaryotic genomes.

Sequence Read Depth Analysis of a Monophyletic Cluster of Y Chromosomes Characterized by Structural Rearrangements in the AZFc Region Resulting in DYS448 Deletion and DYF387S1 Duplication

The azoospermia factor c region (AZFc), located in the long arm of the human Y chromosome, is frequently involved in chromosome rearrangements, mainly due to non-allelic homologous recombination events that occur between the nearly identical sequences (amplicon) that comprises it. These rearrangements may have major phenotypic effects like spermatogenic failure or other pathologies linked to male infertility. Moreover, they may also be relevant in forensic genetics, since some of the Y chromosome short tandem repeats (Y-STRs) commonly used in forensic analysis are located in amplicons or in inter-amplicon sequences of the AZFc. In a previous study, we identified four phylogenetically related samples with a null allele at DYS448 and a tetrallelic pattern at DYF387S1, two Y-STRs located in the AZFc. Through NGS read depth analysis, we found that the unusual Y-STR pattern may be due to a 1.6 Mb deletion arising concurrently or after a 3.5 Mb duplication event. The observed large genomic rearrangement results in copy number reduction for the RBMY gene family as well as duplication of other AZFc genes. Based on the diversity of 16 additional Y-STRs, we estimated that the duplication/deletion event occurred at least twenty generations ago, suggesting that it has not been affected by negative selection.

A common 1.6 mb Y-chromosomal inversion predisposes to subsequent deletions and severe spermatogenic failure in humans

Male infertility is a prevalent condition, affecting 5–10% of men. So far, few genetic factors have been described as contributors to spermatogenic failure. Here, we report the first re-sequencing study of the Y-chromosomal Azoospermia Factor c (AZFc) region, combined with gene dosage analysis of the multicopy DAZ, BPY2, and CDYgenes and Y-haplogroup determination. In analysing 2324 Estonian men, we uncovered a novel structural variant as a high-penetrance risk factor for male infertility. The Y lineage R1a1-M458, reported at >20% frequency in several European populations, carries a fixed ~1.6 Mb r2/r3 inversion, destabilizing the AZFc region and predisposing to large recurrent microdeletions. Such complex rearrangements were significantly enriched among severe oligozoospermia cases. The carrier vs non-carrier risk for spermatogenic failure was increased 8.6-fold (p=6.0×10⁻⁴). This finding contributes to improved molecular diagnostics and clinical management of infertility. Carrier identification at young age will facilitate timely counselling and reproductive decision-making.

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.

Suppression of Non-Random Fertilization by MHC Class I Antigens

Hermaphroditic invertebrates and plants have a self-recognition system on the cell surface of sperm and eggs, which prevents their self-fusion and enhances non-self-fusion, thereby contributing to genetic variation. However, the system of sperm–egg recognition in mammals is under debate. To address this issue, we explored the role of major histocompatibility complex class I (MHC class I, also known as histocompatibility 2-K^b or H2-K^b and H2-D^b in mice) antigens by analyzing H2-K^b-/-H2-D^b-/-β2-microglobulin (β2M)^-/- triple-knockout (T-KO) male mice with full fertility. T-KO sperm exhibited an increased sperm number in the perivitelline space of wild-type (WT) eggs in vitro. Moreover, T-KO sperm showed multiple fusion with zona pellucida (ZP)-free WT eggs, implying that the ability of polyspermy block for sperm from T-KO males was weakened in WT eggs. When T-KO male mice were intercrossed with WT female mice, the percentage of females in progeny increased. We speculate that WT eggs prefer fusion with T-KO sperm, more specifically X-chromosome-bearing sperm (X sperm), suggesting the presence of preferential (non-random) fertilization in mammals, including humans.

The Y Chromosome: A Complex Locus for Genetic Analyses of Complex Human Traits

The Human Y chromosome (ChrY) has been demonstrated to be a powerful tool for phylogenetics, population genetics, genetic genealogy and forensics. However, the importance of ChrY genetic variation in relation to human complex traits is less clear. In this review, we summarise existing evidence about the inherent complexities of ChrY variation and their use in association studies of human complex traits. We present and discuss the specific particularities of ChrY genetic variation, including Y chromosomal haplogroups, that need to be considered in the design and interpretation of genetic epidemiological studies involving ChrY.

Extinction of chromosomes due to specialization is a universal occurrence

The human X and Y chromosomes evolved from a pair of autosomes approximately 180 million years ago. Despite their shared evolutionary origin, extensive genetic decay has resulted in the human Y chromosome losing 97% of its ancestral genes while gene content and order remain highly conserved on the X chromosome. Five 'stratification' events, most likely inversions, reduced the Y chromosome's ability to recombine with the X chromosome across the majority of its length and subjected its genes to the erosive forces associated with reduced recombination. The remaining functional genes are ubiquitously expressed, functionally coherent, dosage-sensitive genes, or have evolved male-specific functionality. It is unknown, however, whether functional specialization is a degenerative phenomenon unique to sex chromosomes, or if it conveys a potential selective advantage aside from sexual antagonism. We examined the evolution of mammalian orthologs to determine if the selective forces that led to the degeneration of the Y chromosome are unique in the genome. The results of our study suggest these forces are not exclusive to the Y chromosome, and chromosomal degeneration may have occurred throughout our evolutionary history. The reduction of recombination could additionally result in rapid fixation through isolation of specialized functions resulting in a cost-benefit relationship during times of intense selective pressure.

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.

Birth, expansion, and death of VCY-containing palindromes on the human Y chromosome

BACKGROUND

Large palindromes (inverted repeats) make up substantial proportions of mammalian sex chromosomes, often contain genes, and have high rates of structural variation arising via ectopic recombination. As a result, they underlie many genomic disorders. Maintenance of the palindromic structure by gene conversion between the arms has been documented, but over longer time periods, palindromes are remarkably labile. Mechanisms of origin and loss of palindromes have, however, received little attention.

RESULTS

Here, we use fiber-FISH, 10x Genomics Linked-Read sequencing, and breakpoint PCR sequencing to characterize the structural variation of the P8 palindrome on the human Y chromosome, which contains two copies of the VCY (Variable Charge Y) gene. We find a deletion of almost an entire arm of the palindrome, leading to death of the palindrome, a size increase by recruitment of adjacent sequence, and other complex changes including the formation of an entire new palindrome nearby. Together, these changes are found in ~ 1% of men, and we can assign likely molecular mechanisms to these mutational events. As a result, healthy men can have 1-4 copies of VCY.

CONCLUSIONS

Gross changes, especially duplications, in palindrome structure can be relatively frequent and facilitate the evolution of sex chromosomes in humans, and potentially also in other mammalian species.

Dosage regulation, and variation in gene expression and copy number of human Y chromosome ampliconic genes

The Y chromosome harbors nine multi-copy ampliconic gene families expressed exclusively in testis. The gene copies within each family are >99% identical to each other, which poses a major challenge in evaluating their copy number. Recent studies demonstrated high variation in Y ampliconic gene copy number among humans. However, how this variation affects expression levels in human testis remains understudied. Here we developed a novel computational tool Ampliconic Copy Number Estimator (AmpliCoNE) that utilizes read sequencing depth information to estimate Y ampliconic gene copy number per family. We applied this tool to whole-genome sequencing data of 149 men with matched testis expression data whose samples are part of the Genotype-Tissue Expression (GTEx) project. We found that the Y ampliconic gene families with low copy number in humans were deleted or pseudogenized in non-human great apes, suggesting relaxation of functional constraints. Among the Y ampliconic gene families, higher copy number leads to higher expression. Within the Y ampliconic gene families, copy number does not influence gene expression, rather a high tolerance for variation in gene expression was observed in testis of presumably healthy men. No differences in gene expression levels were found among major Y haplogroups. Age positively correlated with expression levels of the HSFY and PRY gene families in the African subhaplogroup E1b, but not in the European subhaplogroups R1b and I1. We also found that expression of five Y ampliconic gene families is coordinated with that of their non-Y (i.e. X or autosomal) homologs. Indeed, five ampliconic gene families had consistently lower expression levels when compared to their non-Y homologs suggesting dosage regulation, while the HSFY family had higher expression levels than its X homolog and thus lacked dosage regulation.

Evolutionary and functional analysis of RBMY1 gene copy number variation on the human Y chromosome

Human RBMY1 genes are located in four variable-sized clusters on the Y chromosome, expressed in male germ cells and possibly associated with sperm motility. We have re-investigated the mutational background and evolutionary history of the RBMY1 copy number distribution in worldwide samples and its relevance to sperm parameters in an Estonian cohort of idiopathic male factor infertility subjects. We estimated approximate RBMY1 copy numbers in 1218 1000 Genomes Project phase 3 males from sequencing read-depth, then chose 14 for valid ation by multicolour fibre-FISH. These fibre-FISH samples provided accurate calibration standards for the entire panel and led to detailed insights into population variation and mutational mechanisms. RBMY1 copy number worldwide ranged from 3 to 13 with a mode of 8. The two larger proximal clusters were the most variable, and additional duplications, deletions and inversions were detected. Placing the copy number estimates onto the published Y-SNP-based phylogeny of the same samples suggested a minimum of 562 mutational changes, translating to a mutation rate of 2.20 × 10-3 (95% CI 1.94 × 10-3 to 2.48 × 10-3) per father-to-son Y-transmission, higher than many short tandem repeat (Y-STRs), and showed no evidence for selection for increased or decreased copy number, but possible copy number stabilizing selection. An analysis of RBMY1 copy numbers among 376 infertility subjects failed to replicate a previously reported association with sperm motility and showed no significant effect on sperm count and concentration, serum follicle stimulating hormone (FSH), luteinizing hormone (LH) and testosterone levels or testicular and semen volume. These results provide the first in-depth insights into the structural rearrangements underlying RBMY1 copy number variation across diverse human lineages.

gr/gr deletion predisposes to testicular germ cell tumour independently from altered spermatogenesis: results from the largest European study

The association between impaired spermatogenesis and TGCT has stimulated research on shared genetic factors. Y chromosome-linked partial AZFc deletions predispose to oligozoospermia and were also studied in TGCT patients with controversial results. In the largest study reporting the association between gr/gr deletion and TGCT, sperm parameters were unknown. Hence, it remains to be established whether this genetic defect truly represents a common genetic link between TGCT and impaired sperm production. Our aim was to explore the role of the following Y chromosome-linked factors in the predisposition to TGCT: (i) gr/gr deletion in subjects with known sperm parameters; (ii) other partial AZFc deletions and, for the first time, the role of partial AZFc duplications; (iii) DAZ gene dosage variation. 497 TGCT patients and 2030 controls from two Mediterranean populations with full semen/andrological characterization were analyzed through a series of molecular genetic techniques. Our most interesting finding concerns the gr/gr deletion and DAZ gene dosage variation (i.e., DAZ copy number is different from the reference sequence), both conferring TGCT susceptibility. In particular, the highest risk was observed when normozoospermic TGCT and normozoospermic controls were compared (OR = 3.7; 95% CI = 1.5-9.1; p = 0.006 for gr/gr deletion and OR = 1.8; 95% CI = 1.1-3.0; p = 0.013 for DAZ gene dosage alteration). We report in the largest European study population the predisposing effect of gr/gr deletion to TGCT as an independent risk factor from impaired spermatogenesis. Our finding implies regular tumour screening/follow-up in male family members of TGCT patients with gr/gr deletion and in infertile gr/gr deletion carriers.