Pseudoautosomal region 1 length polymorphism in the human population

Evolutionary dynamics of pseudoautosomal region 1 in humans and great apes

Background

The pseudoautosomal region 1 (PAR1) is a 2.7 Mb telomeric region of human sex chromosomes. PAR1 has a crucial role in ensuring proper segregation of sex chromosomes during male meiosis, exposing it to extreme recombination and mutation processes. We investigate PAR1 evolution using population genomic datasets of extant humans, eight populations of great apes, and two archaic human genome sequences.

Results

We find that PAR1 is fast evolving and closer to evolutionary nucleotide equilibrium than autosomal telomeres. We detect a difference between substitution patterns and extant diversity in PAR1, mainly driven by the conflict between strong mutation and recombination-associated fixation bias at CpG sites. We detect excess C-to-G mutations in PAR1 of all great apes, specific to the mutagenic effect of male recombination. Despite recent evidence for Y chromosome introgression from humans into Neanderthals, we find that the Neanderthal PAR1 retained similarity to the Denisovan sequence. We find differences between substitution spectra of these archaics suggesting rapid evolution of PAR1 in recent hominin history. Frequency analysis of alleles segregating in females and males provided no evidence for recent sexual antagonism in this region. We study repeat content and double-strand break hotspot regions in PAR1 and find that they may play roles in ensuring the obligate X-Y recombination event during male meiosis.

Conclusions

Our study provides an unprecedented quantification of population genetic forces governing PAR1 biology across extant and extinct hominids. PAR1 evolutionary dynamics are predominantly governed by recombination processes with a strong impact on mutation patterns across all species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02784-x.

The genes from the pseudoautosomal region 1 (PAR1) of the mammalian sex chromosomes: Synteny, phylogeny and selection

Sex chromosomes recombine restrictly in their homologous area, the pseudoautosomal region (PAR), represented by PAR1 and PAR2, which behave like an autosome in both pairing and recombination. The PAR1, common to most of the eutherian mammals, is located at the terminus of the sex chromosomes short arm and exhibit recombination rates ~20 times higher than the autosomes. Here, we assessed the interspecific evolutionary genomic dynamics of 15 genes of the PAR1 across 41 mammalian genera (representing six orders). The strong negative selection detected in most of the assessed groups reinforces the presence of evolutionary constraints, imposed by the important function of the PAR1 genes. Indeed, mutations in these genes are associated with various diseases in humans, including stature problems (Klinefelter Syndrome), leukemia and mental diseases. Yet, a few genes exhibiting positive selection (ω-value >1) were depicted in Rodentia (ASMT and ZBED1) and Primates (CRLF2 and CSF2RA). Rodents have the smallest described PAR1, while that of simian primates/humans underwent a 3 to 5 fold size reduction. The assessment of the PAR1 genes synteny revealed differences among the mammalian species, especially in the Rodentia order where chromosomic translocations from the sex chromosomes to the autosomes were observed. Such syntenic changes may be an evidence of the rapid evolution in rodents, as previous referred in other papers, also depicted by their increased branch lengths in the phylogenetic analyses. Concluding, we suggest that genome migration is an important factor influencing the evolution of mammals and may result in changes of the selective pressures operating on the genome.

Major sex differences in allele frequencies for X chromosomal variants in both the 1000 Genomes Project and gnomAD

An unexpectedly high proportion of SNPs on the X chromosome in the 1000 Genomes Project phase 3 data were identified with significant sex differences in minor allele frequencies (sdMAF). sdMAF persisted for many of these SNPs in the recently released high coverage whole genome sequence of the 1000 Genomes Project that was aligned to GRCh38, and it was consistent between the five super-populations. Among the 245,825 common (MAF>5%) biallelic X-chromosomal SNPs in the phase 3 data presumed to be of high quality, 2,039 have genome-wide significant sdMAF (p-value <5e-8). sdMAF varied by location: non-pseudo-autosomal region (NPR) = 0.83%, pseudo-autosomal regions (PAR1) = 0.29%, PAR2 = 13.1%, and X-transposed region (XTR)/PAR3 = 0.85% of SNPs had sdMAF, and they were clustered at the NPR-PAR boundaries, among others. sdMAF at the NPR-PAR boundaries are biologically expected due to sex-linkage, but have generally been ignored in association studies. For comparison, similar analyses found only 6, 1 and 0 SNPs with significant sdMAF on chromosomes 1, 7 and 22, respectively. Similar sdMAF results for the X chromosome were obtained from the high coverage whole genome sequence data from gnomAD V 3.1.2 for both the non-Finnish European and African/African American samples. Future X chromosome analyses need to take sdMAF into account.

The human X chromosome contains over 800 genes and is the 8^th largest human chromosome. Genome-wide associations studies have generally failed to examine variants on the X chromosome for association with diseases and traits, partly due to complexities of the data analysis, and challenges with genotype imputation. We examined X chromosomal variants from the 1000 Genomes Project for sex differences in allele frequency and found that many variants showed significant differences. These variants cluster at the centromeric parts of the pseudoautosomal regions 1 and 2, as well as the putative pseudo-autosomal region 3 (also termed X-transposed region). This pattern was observed in high coverage whole genome sequence data from the same subjects that was aligned to GRCh38, suggesting that is not an artefact of low coverage sequencing or problems specific to GRCh37. In addition, we replicated this phenomenon in high coverage whole genome sequence aligned to GRCh38 from the gnomAD database in both the non-Finnish European and African/African American samples. These findings have implications for the analysis of X chromosomal variants for disease and trait associations.

Human Spermatogenesis Tolerates Massive Size Reduction of the Pseudoautosomal Region

Mammalian male meiosis requires homologous recombination between the X and Y chromosomes. In humans, such recombination occurs exclusively in the short arm pseudoautosomal region (PAR1) of 2.699 Mb in size. Although it is known that complete deletion of PAR1 causes spermatogenic arrest, no studies have addressed to what extent male meiosis tolerates PAR1 size reduction. Here, we report two families in which PAR1 partial deletions were transmitted from fathers to their offspring. Cytogenetic analyses revealed that a ∼400-kb segment at the centromeric end of PAR1, which accounts for only 14.8% of normal PAR1 and 0.26% and 0.68% of the X and Y chromosomes, respectively, is sufficient to mediate sex chromosomal recombination during spermatogenesis. These results highlight the extreme recombinogenic activity of human PAR1. Our data, in conjunction with previous findings from animal studies, indicate that the minimal size requirement of mammalian PARs to maintain male fertility is fairly small.

Heterogeneous transposable elements as silencers, enhancers and targets of meiotic recombination

During meiosis, DNA double-strand breaks are initiated by the topoisomerase-like enzyme SPO11 and are repaired by inter-sister chromatid and inter-homologue DNA repair pathways. Genome-wide maps of initiating DNA double-strand breaks and inter-homologue repair events are now available for a number of mammalian, fungal and plant species. In mammals, PRDM9 specifies the location of meiotic recombination initiation via recognition of specific DNA sequence motifs by its C2H2 zinc finger array. In fungi and plants, meiotic recombination appears to be initiated less discriminately in accessible chromatin, including at gene promoters. Generally, meiotic crossover is suppressed in highly repetitive genomic regions that are made up of transposable elements (TEs), to prevent deleterious non-allelic homologous recombination events. However, recent and older studies have revealed intriguing relationships between meiotic recombination initiation and repair, and transposable elements. For instance, gene conversion events have been detected in maize centromeric retroelements, mouse MULE-MuDR DNA transposons undergo substantial meiotic recombination initiation, Arabidopsis Helitron TEs are among the hottest of recombination initiation hotspots, and human TE sequences can modify the crossover rate at adjacent PRDM9 motifs in cis. Here, we summarize the relationship between meiotic recombination and TEs, discuss recent insights from highly divergent eukaryotes and highlight outstanding questions in the field.

Single molecule real-time (SMRT) sequencing comes of age: applications and utilities for medical diagnostics

Short read massive parallel sequencing has emerged as a standard diagnostic tool in the medical setting. However, short read technologies have inherent limitations such as GC bias, difficulties mapping to repetitive elements, trouble discriminating paralogous sequences, and difficulties in phasing alleles. Long read single molecule sequencers resolve these obstacles. Moreover, they offer higher consensus accuracies and can detect epigenetic modifications from native DNA. The first commercially available long read single molecule platform was the RS system based on PacBio's single molecule real-time (SMRT) sequencing technology, which has since evolved into their RSII and Sequel systems. Here we capsulize how SMRT sequencing is revolutionizing constitutional, reproductive, cancer, microbial and viral genetic testing.

Instability of the Pseudoautosomal Boundary in House Mice

Faithful segregation of homologous chromosomes at meiosis requires pairing and recombination. In taxa with dimorphic sex chromosomes, pairing between them in the heterogametic sex is limited to a narrow interval of residual sequence homology known as the pseudoautosomal region (PAR). Failure to form the obligate crossover in the PAR is associated with male infertility in house mice (Mus musculus) and humans. Yet despite this apparent functional constraint, the boundary and organization of the PAR is highly variable in mammals, and even between subspecies of mice. Here, we estimate the genetic map in a previously documented expansion of the PAR in the M. musculus castaneus subspecies and show that the local recombination rate is 100-fold higher than the autosomal background. We identify an independent shift in the PAR boundary in the M. musculus musculus subspecies and show that it involves a complex rearrangement, but still recombines in heterozygous males. Finally, we demonstrate pervasive copy-number variation at the PAR boundary in wild populations of M. m. domesticus, M. m. musculus, and M. m. castaneus Our results suggest that the intensity of recombination activity in the PAR, coupled with relatively weak constraints on its sequence, permit the generation and maintenance of unusual levels of polymorphism in the population of unknown functional significance.

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

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

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

Mediterranean Y-chromosome 2.0-why the Y in the Mediterranean is still relevant in the postgenomic era

CONTEXT

Due to its unique paternal inheritance, the Y-chromosome has been a highly popular marker among population geneticists for over two decades. Recently, the advent of cost-effective genome-wide methods has unlocked information-rich autosomal genomic data, paving the way to the postgenomic era. This seems to have announced the decreasing popularity of investigating Y-chromosome variation, which provides only the paternal perspective of human ancestries and is strongly influenced by genetic drift and social behaviour.

OBJECTIVE

For this special issue on population genetics of the Mediterranean, the aim was to demonstrate that the Y-chromosome still provides important insights in the postgenomic era and in a time when ancient genomes are becoming exponentially available.

METHODS

A systematic literature search on Y-chromosomal studies in the Mediterranean was performed.

RESULTS

Several applications of Y-chromosomal analysis with future opportunities are formulated and illustrated with studies on Mediterranean populations.

CONCLUSIONS

There will be no reduced interest in Y-chromosomal studies going from reconstruction of male-specific demographic events to ancient DNA applications, surname history and population-wide estimations of extra-pair paternity rates. Moreover, more initiatives are required to collect population genetic data of Y-chromosomal markers for forensic research, and to include Y-chromosomal data in GWAS investigations and studies on male infertility.

Dosage compensation in the process of inactivation/reactivation during both germ cell development and early embryogenesis in mouse

Ohno proposed that dosage compensation in mammals evolved as a two-step mechanism involving X-inactivation and X-upregulation. While X-inactivation is well characterized, it remains to further analysis whether upregulation of the single activated X chromosome in mammals occurs. We obtained RNA-seq data, including single-cell RNA-seq data, from cells undergoing inactivation/reactivation in both germ cell development and early embryogenesis stages in mouse and calculated the X: A ratio from the gene expression. Our results showed that the X: A ratio is always 1, regardless of the number of X chromosomes being transcribed for expressed genes. Furthermore, the single-cell RNA-seq data across individual cells of mouse preimplantation embryos of mixed backgrounds indicated that strain-specific SNPs could be used to distinguish transcription from maternal and paternal chromosomes and further showed that when the paternal was inactivated, the average gene dosage of the active maternal X chromosome was increased to restore the balance between the X chromosome and autosomes. In conclusion, our analysis of RNA-seq data (particularly single-cell RNA-seq) from cells undergoing the process of inactivation/reactivation provides direct evidence that the average gene dosage of the single active X chromosome is upregulated to achieve a similar level to that of two active X chromosomes and autosomes present in two copies.

Meiotic Consequences of Genetic Divergence Across the Murine Pseudoautosomal Region

The production of haploid gametes during meiosis is dependent on the homology-driven processes of pairing, synapsis, and recombination. On the mammalian heterogametic sex chromosomes, these key meiotic activities are confined to the pseudoautosomal region (PAR), a short region of near-perfect sequence homology between the X and Y chromosomes. Despite its established importance for meiosis, the PAR is rapidly evolving, raising the question of how proper X/Y segregation is buffered against the accumulation of homology-disrupting mutations. Here, I investigate the interplay of PAR evolution and function in two interfertile house mouse subspecies characterized by structurally divergent PARs, Mus musculus domesticus and M. m. castaneus. Using cytogenetic methods to visualize the sex chromosomes at meiosis, I show that intersubspecific F₁ hybrids harbor an increased frequency of pachytene spermatocytes with unsynapsed sex chromosomes. This high rate of asynapsis is due, in part, to the premature release of synaptic associations prior to completion of prophase I. Further, I show that when sex chromosomes do synapse in intersubspecific hybrids, recombination is reduced across the paired region. Together, these meiotic defects afflict ∼50% of spermatocytes from F₁ hybrids and lead to increased apoptosis in meiotically dividing cells. Despite flagrant disruption of the meiotic program, a subset of spermatocytes complete meiosis and intersubspecific F₁ males remain fertile. These findings cast light on the meiotic constraints that shape sex chromosome evolution and offer initial clues to resolve the paradox raised by the rapid evolution of this functionally significant locus.

Genetic Diversity on the Human X Chromosome Does Not Support a Strict Pseudoautosomal Boundary

Unlike the autosomes, recombination between the X chromosome and the Y chromosome is often thought to be constrained to two small pseudoautosomal regions (PARs) at the tips of each sex chromosome. PAR1 spans the first 2.7 Mb of the proximal arm of the human sex chromosomes, whereas the much smaller PAR2 encompasses the distal 320 kb of the long arm of each sex chromosome. In addition to PAR1 and PAR2, there is a human-specific X-transposed region that was duplicated from the X to the Y chromosome. The X-transposed region is often not excluded from X-specific analyses, unlike the PARs, because it is not thought to routinely recombine. Genetic diversity is expected to be higher in recombining regions than in nonrecombining regions because recombination reduces the effect of linked selection. In this study, we investigated patterns of genetic diversity in noncoding regions across the entire X chromosome of a global sample of 26 unrelated genetic females. We found that genetic diversity in PAR1 is significantly greater than in the nonrecombining regions (nonPARs). However, rather than an abrupt drop in diversity at the pseudoautosomal boundary, there is a gradual reduction in diversity from the recombining through the nonrecombining regions, suggesting that recombination between the human sex chromosomes spans across the currently defined pseudoautosomal boundary. A consequence of recombination spanning this boundary potentially includes increasing the rate of sex-linked disorders (e.g., de la Chapelle) and sex chromosome aneuploidies. In contrast, diversity in PAR2 is not significantly elevated compared to the nonPARs, suggesting that recombination is not obligatory in PAR2. Finally, diversity in the X-transposed region is higher than in the surrounding nonPARs, providing evidence that recombination may occur with some frequency between the X and Y chromosomes in the X-transposed region.

Polymerase specific error rates and profiles identified by single molecule sequencing

DNA polymerases have an innate error rate which is polymerase and DNA context specific. Historically the mutational rate and profiles have been measured using a variety of methods, each with their own technical limitations. Here we used the unique properties of single molecule sequencing to evaluate the mutational rate and profiles of six DNA polymerases at the sequence level. In addition to accurately determining mutations in double strands, single molecule sequencing also captures direction specific transversions and transitions through the analysis of heteroduplexes. Not only did the error rates vary, but also the direction specific transitions differed among polymerases.

Evolution of vertebrate sex chromosomes and dosage compensation

Differentiated sex chromosomes in mammals and other vertebrates evolved independently but in strikingly similar ways. Vertebrates with differentiated sex chromosomes share the problems of the unequal expression of the genes borne on sex chromosomes, both between the sexes and with respect to autosomes. Dosage compensation of genes on sex chromosomes is surprisingly variable - and can even be absent - in different vertebrate groups. Systems that compensate for different gene dosages include a wide range of global, regional and gene-by-gene processes that differ in their extent and their molecular mechanisms. However, many elements of these control systems are similar across distant phylogenetic divisions and show parallels to other gene silencing systems. These dosage systems cannot be identical by descent but were probably constructed from elements of ancient silencing mechanisms that are ubiquitous among vertebrates and shared throughout eukaryotes.

Microarray Analysis of Copy Number Variants on the Human Y Chromosome Reveals Novel and Frequent Duplications Overrepresented in Specific Haplogroups

BACKGROUND

The human Y chromosome is almost always excluded from genome-wide investigations of copy number variants (CNVs) due to its highly repetitive structure. This chromosome should not be forgotten, not only for its well-known relevance in male fertility, but also for its involvement in clinical phenotypes such as cancers, heart failure and sex specific effects on brain and behaviour.

RESULTS

We analysed Y chromosome data from Affymetrix 6.0 SNP arrays and found that the signal intensities for most of 8179 SNP/CN probes in the male specific region (MSY) discriminated between a male, background signals in a female and an isodicentric male containing a large deletion of the q-arm and a duplication of the p-arm of the Y chromosome. Therefore, this SNP/CN platform is suitable for identification of gain and loss of Y chromosome sequences. In a set of 1718 males, we found 25 different CNV patterns, many of which are novel. We confirmed some of these variants by PCR or qPCR. The total frequency of individuals with CNVs was 14.7%, including 9.5% with duplications, 4.5% with deletions and 0.7% exhibiting both. Hence, a novel observation is that the frequency of duplications was more than twice the frequency of deletions. Another striking result was that 10 of the 25 detected variants were significantly overrepresented in one or more haplogroups, demonstrating the importance to control for haplogroups in genome-wide investigations to avoid stratification. NO-M214(xM175) individuals presented the highest percentage (95%) of CNVs. If they were not counted, 12.4% of the rest included CNVs, and the difference between duplications (8.9%) and deletions (2.8%) was even larger.

CONCLUSIONS

Our results demonstrate that currently available genome-wide SNP platforms can be used to identify duplications and deletions in the human Y chromosome. Future association studies of the full spectrum of Y chromosome variants will demonstrate the potential involvement of gain or loss of Y chromosome sequence in different human phenotypes.

RAD mapping reveals an evolving, polymorphic and fuzzy boundary of a plant pseudoautosomal region

How loss of genetic exchanges (recombination) evolves between sex chromosomes is a long-standing question. Suppressed recombination may evolve when a sexually antagonistic (SA) polymorphism occurs in a partially sex-linked 'pseudoautosomal' region (or 'PAR'), maintaining allele frequency differences between the two sexes, and creating selection for closer linkage with the fully sex-linked region of the Y chromosome in XY systems, or the W in ZW sex chromosome systems. Most evidence consistent with the SA polymorphism hypothesis is currently indirect, and more studies of the genetics and population genetics of PAR genes are clearly needed. The sex chromosomes of the plant Silene latifolia are suitable for such studies, as they evolved recently and the loss of recombination could still be ongoing. Here, we used RAD sequencing to genetically map sequences in this plant, which has a large genome (c. 3 gigabases) and no available whole-genome sequence. We mapped 83 genes on the sex chromosomes, and comparative mapping in the related species S. vulgaris supports previous evidence for additions to an ancestral PAR and identified at least 12 PAR genes. We describe evidence that recombination rates have been reduced in meiosis of both sexes, and differences in recombination between S. latifolia families suggest ongoing recombination suppression. Large allele frequency differences between the sexes were found at several loci closely linked to the PAR boundary, and genes in different regions of the PAR showed striking sequence diversity patterns that help illuminate the evolution of the PAR.