Abundant gene conversion between arms of palindromes in human and ape Y chromosomes

The complete sequence and comparative analysis of ape sex chromosomes

Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility1. The X chromosome is vital for reproduction and cognition2. Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species.

EAA/EMQN best practice guidelines for molecular diagnosis of Y-chromosomal microdeletions: State of the art 2023

Testing for AZoospermia Factor (AZF) deletions of the Y chromosome is a key component of the diagnostic workup of azoospermic and severely oligozoospermic men. This revision of the 2013 European Academy of Andrology (EAA) and EMQN CIC (previously known as the European Molecular Genetics Quality Network) laboratory guidelines summarizes recent clinically relevant advances and provides an update on the results of the external quality assessment program jointly offered by both organizations. A basic multiplex PCR reaction followed by a deletion extension analysis remains the gold-standard methodology to detect and correctly interpret AZF deletions. Recent data have led to an update of the sY84 reverse primer sequence, as well as to a refinement of what were previously considered as interchangeable border markers for AZFa and AZFb deletion breakpoints. More specifically, sY83 and sY143 are no longer recommended for the deletion extension analysis, leaving sY1064 and sY1192, respectively, as first-choice markers. Despite the transition, currently underway in several countries, toward a diagnosis based on certified kits, it should be noted that many of these commercial products are not recommended due to an unnecessarily high number of tested markers, and none of those currently available are, to the best of our knowledge, in accordance with the new first-choice markers for the deletion extension analysis. The gr/gr partial AZFc deletion remains a population-specific risk factor for impaired sperm production and a predisposing factor for testicular germ cell tumors. Testing for this deletion type is, as before, left at the discretion of the diagnostic labs and referring clinicians. Annual participation in an external quality control program is strongly encouraged, as the 22-year experience of the EMQN/EAA scheme clearly demonstrates a steep decline in diagnostic errors and an improvement in reporting practice.

De novo transcriptome assembly of mouse male germ cells reveals novel genes, stage-specific bidirectional promoter activity, and noncoding RNA expression

In mammals, the adult testis is the tissue with the highest diversity in gene expression. Much of that diversity is attributed to germ cells, primarily meiotic spermatocytes and postmeiotic haploid spermatids. Exploiting a newly developed cell purification method, we profiled the transcriptomes of such postmitotic germ cells of mice. We used a de novo transcriptome assembly approach and identified thousands of novel expressed transcripts characterized by features distinct from those of known genes. Novel loci tend to be short in length, monoexonic, and lowly expressed. Most novel genes have arisen recently in evolutionary time and possess low coding potential. Nonetheless, we identify several novel protein-coding genes harboring open reading frames that encode proteins containing matches to conserved protein domains. Analysis of mass-spectrometry data from adult mouse testes confirms protein production from several of these novel genes. We also examine overlap between transcripts and repetitive elements. We find that although distinct families of repeats are expressed with differing temporal dynamics during spermatogenesis, we do not observe a general mode of regulation wherein repeats drive expression of nonrepetitive sequences in a cell type-specific manner. Finally, we observe many fairly long antisense transcripts originating from canonical gene promoters, pointing to pervasive bidirectional promoter activity during spermatogenesis that is distinct and more frequent compared with somatic cells.

The Complete Sequence and Comparative Analysis of Ape Sex Chromosomes

Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.

Genetic conflict and the origin of multigene families: implications for sex chromosome evolution

Sex chromosomes are havens for intragenomic conflicts. The absence of recombination between sex chromosomes creates the opportunity for the evolution of segregation distorters: selfish genetic elements that hijack different aspects of an individual's reproduction to increase their own transmission. Biased (non-Mendelian) segregation, however, often occurs at a detriment to their host's fitness, and therefore can trigger evolutionary arms races that can have major consequences for genome structure and regulation, gametogenesis, reproductive strategies and even speciation. Here, we review an emerging feature from comparative genomic and sex chromosome evolution studies suggesting that meiotic drive is pervasive: the recurrent evolution of paralogous sex-linked gene families. Sex chromosomes of several species independently acquire and co-amplify rapidly evolving gene families with spermatogenesis-related functions, consistent with a history of intragenomic conflict over transmission. We discuss Y chromosome features that might contribute to the tempo and mode of evolution of X/Y co-amplified gene families, as well as their implications for the evolution of complexity in the genome. Finally, we propose a framework that explores the conditions that might allow for recurrent bouts of fixation of drivers and suppressors, in a dosage-sensitive fashion, and therefore the co-amplification of multigene families on sex 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.

Eighty million years of rapid evolution of the primate Y chromosome

The Y chromosome usually plays a critical role in determining male sex and comprises sequence classes that have experienced unique evolutionary trajectories. Here we generated 19 new primate sex chromosome assemblies, analysed them with 10 existing assemblies and report rapid evolution of the Y chromosome across primates. The pseudoautosomal boundary has shifted at least six times during primate evolution, leading to the formation of a Simiiformes-specific evolutionary stratum and to the independent start of young strata in Catarrhini and Platyrrhini. Different primate lineages experienced different rates of gene loss and structural and chromatin change on their Y chromosomes. Selection on several Y-linked genes has contributed to the evolution of male developmental traits across the primates. Additionally, lineage-specific expansions of ampliconic regions have further increased the diversification of the structure and gene composition of the Y chromosome. Overall, our comprehensive analysis has broadened our knowledge of the evolution of the primate Y chromosome.

Dynamic alternative DNA structures in biology and disease

Repetitive elements in the human genome, once considered 'junk DNA', are now known to adopt more than a dozen alternative (that is, non-B) DNA structures, such as self-annealed hairpins, left-handed Z-DNA, three-stranded triplexes (H-DNA) or four-stranded guanine quadruplex structures (G4 DNA). These dynamic conformations can act as functional genomic elements involved in DNA replication and transcription, chromatin organization and genome stability. In addition, recent studies have revealed a role for these alternative structures in triggering error-generating DNA repair processes, thereby actively enabling genome plasticity. As a driving force for genetic variation, non-B DNA structures thus contribute to both disease aetiology and evolution.

The first gapless, reference-quality, fully annotated genome from a Southern Han Chinese individual

We used long-read DNA sequencing to assemble the genome of a Southern Han Chinese male. We organized the sequence into chromosomes and filled in gaps using the recently completed T2T-CHM13 genome as a guide, yielding a gap-free genome, Han1, containing 3,099,707,698 bases. Using the T2T-CHM13 annotation as a reference, we mapped all genes onto the Han1 genome and identified additional gene copies, generating a total of 60,708 putative genes, of which 20,003 are protein-coding. A comprehensive comparison between the genes revealed that 235 protein-coding genes were substantially different between the individuals, with frameshifts or truncations affecting the protein-coding sequence. Most of these were heterozygous variants in which one gene copy was unaffected. This represents the first gene-level comparison between two finished, annotated individual human genomes.

Haplotypic polymorphisms and forensic applications in Chinese Hunan Han population based on a series of Y-STR loci: a perspective of paternal inheritance

BACKGROUND

Hunan, a multinational province in China, possesses more than fifty ethnic groups, such as the Han, Yao, Tujia, Miao, and so on.

AIMS

To evaluate the forensic efficiency of the novel panel and investigate the genetic relationships between Hunan Han population and 12 other reference populations from China.

SUBJECTS AND METHODS

Haplotypic data of 153 unrelated males of Hunan Han population were investigated using the AGCU Y SUPP Plus amplification system containing 27 Y-chromosomal short tandem repeat (STR) loci. Forensic parameters were calculated to evaluate the application efficiency of this panel in Hunan Han population.

RESULTS

Haplotype diversity, discrimination capability, and match probability values were 0.9999999977, 1.0000, and 0.0065, respectively. Pairwise fixation index values demonstrated that the minimal genetic differentiation (0.0073) was found between Hunan Han population and Hunan Yao group, while the maximal genetic differentiation (0.0651) was observed between Hunan Han and Guangxi Yao group from the perspective of the patrilineal DNA analysis.

CONCLUSIONS

The haplotype distributions of 27 Y-STR loci in Hunan Han population exhibited remarkable polymorphisms. Moreover, this panel has potential advantages for the forensic applications regarding family investigations, paternity testing of the paternal line, and population genetics.

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.

Developmental validation of a 381 Y-chromosome SNP panel for haplogroup analysis in the Chinese populations

Y-chromosome single nucleotide polymorphism (Y-SNP) shows great variation in geographical distribution and population heterogeneity and can be used to map population genetics around the world. Massive parallel sequencing (MPS) methodology enables high-resolution Y-SNP haplogrouping for a certain male and is widely used in forensic genetics and evolutionary studies. In this present study, we used MPS to develop a customized 381 Y-SNP panel (SifaMPS 381 Y-SNP panel) to investigate the basic structure and subbranches of the haplogroup tree of the Chinese populations. The SifaMPS 381 Y-SNP panel covers all the Y-SNPs from our previously designed 183 Y-SNP panel and additional SNPs under the predominant haplogroups in the Chinese populations based on certain criteria. We also evaluated the sequencing matrix, concordance, sensitivity, repeatability of this panel and the ability to analyze mixed and case-type samples based on the Illumina MiSeq System. The results demonstrated that the novel MPS Y-SNP panel possessed good sequencing performance and generated accurate Y-SNP genotyping results. Although the recommended DNA input was greater than 1.25 ng, we observed that a lower DNA amount could still be used to analyze haplogroups correctly. In addition, this panel could handle mixed samples and common case-type samples and had higher resolution among Chinese Han males than previously reported. In conclusion, the SifaMPS 381 Y-SNP panel showed an overall good performance and offers a better choice for Y-SNP haplogrouping of the Chinese population, thereby facilitating paternal lineage classification, familial searching and other forensic applications.

Chromosome-level genome and recombination map of the male buffalo

BACKGROUND

The swamp buffalo (Bubalus bubalis carabanesis) is an economically important livestock supplying milk, meat, leather, and draft power. Several female buffalo genomes have been available, but the lack of high-quality male genomes hinders studies on chromosome evolution, especially Y, as well as meiotic recombination.

RESULTS

Here, a chromosome-level genome with a contig N50 of 72.2 Mb and a fine-scale recombination map of male buffalo were reported. We found that transposable elements (TEs) and structural variants (SVs) may contribute to buffalo evolution by influencing adjacent gene expression. We further found that the pseudoautosomal region (PAR) of the Y chromosome is subject to stronger purification selection. The meiotic recombination map showed that there were 2 obvious recombination hotspots on chromosome 8, and the genes around them were mainly related to tooth development, which may have helped to enhance the adaption of buffalo to inferior feed. Among several genomic features, TE density has the strongest correlation with recombination rates. Moreover, the TE subfamily, SINE/tRNA, is likely to play a role in driving recombination into SVs.

CONCLUSIONS

The male genome and sperm sequencing will facilitate the understanding of the buffalo genomic evolution and functional research.

Chromosome-scale assemblies of the male and female Populus euphratica genomes reveal the molecular basis of sex determination and sexual dimorphism

Reference-quality genomes of both sexes are essential for studying sex determination and sex-chromosome evolution, as their gene contents and expression profiles differ. Here, we present independent chromosome-level genome assemblies for the female (XX) and male (XY) genomes of desert poplar (Populus euphratica), resolving a 22.7-Mb X and 24.8-Mb Y chromosome. We also identified a relatively complete 761-kb sex-linked region (SLR) in the peritelomeric region on chromosome 14 (Y). Within the SLR, recombination around the partial repeats for the feminizing factor ARR17 (ARABIDOPSIS RESPONSE REGULATOR 17) was potentially suppressed by flanking palindromic arms and the dense accumulation of retrotransposons. The inverted small segments S1 and S2 of ARR17 exhibited relaxed selective pressure and triggered sex determination by generating 24-nt small interfering RNAs that induce male-specific hyper-methylation at the promoter of the autosomal targeted ARR17. We also detected two male-specific fusion genes encoding proteins with NB-ARC domains at the breakpoint region of an inversion in the SLR that may be responsible for the observed sexual dimorphism in immune responses. Our results show that the SLR appears to follow proposed evolutionary dynamics for sex chromosomes and advance our understanding of sex determination and the evolution of sex chromosomes in Populus.

Reference-quality genomes of both sexes of the dioecious tree species, Populus euphratica, provide further insight into the evolution of Populus sex chromosomes and highlight male-specific fusion genes that may contribute to sexual dimorphism.

The rattlesnake W chromosome: A GC-rich retroelement refugium with retained gene function across ancient evolutionary strata

Sex chromosomes diverge after the establishment of recombination suppression, resulting in differential sex-linkage of genes involved in genetic sex determination and dimorphic traits. This process produces systems of male or female heterogamety wherein the Y and W chromosomes are only present in one sex and are often highly degenerated. Sex-limited Y and W chromosomes contain valuable information about the evolutionary transition from autosomes to sex chromosomes, yet detailed characterizations of the structure, composition, and gene content of sex-limited chromosomes are lacking for many species. In this study, we characterize the female-specific W chromosome of the prairie rattlesnake (Crotalus viridis) and evaluate how recombination suppression and other processes have shaped sex chromosome evolution in ZW snakes. Our analyses indicate that the rattlesnake W chromosome is over 80% repetitive and that an abundance of GC-rich mdg4 elements has driven an overall high degree of GC-richness despite a lack of recombination. The W chromosome is also highly enriched for repeat sequences derived from endogenous retroviruses and likely acts as a “refugium” for these and other retroelements. We annotated 219 putatively functional W-linked genes across at least two evolutionary strata identified based on estimates of sequence divergence between Z and W gametologs. The youngest of these strata is relatively gene-rich, however gene expression across strata suggests retained gene function amidst a greater degree of degeneration following ancient recombination suppression. Functional annotation of W-linked genes indicates a specialization of the W chromosome for reproductive and developmental function since recombination suppression from the Z chromosome.

Repetitive DNA Sequences in the Human Y Chromosome and Male Infertility

The male-specific Y chromosome, which is well known for its diverse and complex repetitive sequences, has different sizes, genome structures, contents and evolutionary trajectories from other chromosomes and is of great significance for testis development and function. The large number of repetitive sequences and palindrome structure of the Y chromosome play an important role in maintaining the stability of male sex determining genes, although they can also cause non-allelic homologous recombination within the chromosome. Deletion of certain Y chromosome sequences will lead to spermatogenesis disorders and male infertility. And Y chromosome genes are also involved in the occurrence of reproductive system cancers and can increase the susceptibility of other tumors. In addition, the Y chromosome has very special value in the personal identification and parentage testing of male-related cases in forensic medicine because of its unique paternal genetic characteristics. In view of the extremely high frequency and complexity of gene rearrangements and the limitations of sequencing technology, the analysis of Y chromosome sequences and the study of Y-gene function still have many unsolved problems. This article will introduce the structure and repetitive sequence of the Y chromosome, summarize the correlation between Y chromosome various sequence deletions and male infertility for understanding the repetitive sequence of Y chromosome more systematically, in order to provide research motivation for further explore of the molecules mechanism of Y-deletion and male infertility and theoretical foundations for the transformation of basic research into applications in clinical medicine and forensic medicine.

Copy number variation of ZNF280BY across eight sheep breeds and its association with testicular size of Hu sheep

ZNF280BY, a bovid-specific Y chromosome gene, was firstly found to be highly expressed in bovine testis, indicating it may play important roles in testicular development and male fertility. In this study, we firstly cloned the full-length cDNA of ovine ZNF280BY containing 1993 bp, and with a 1632 bp open reading frame. ZNF280BY was predominantly expressed in the testis, and its expression level was significantly higher in large testis than in small testis in Hu sheep at 6 months of age. In addition, the expression level of ZNF280BY significantly increased during testicular development, showing the highest expression level at 12 months of age. ZNF280BY showed copy number variation (CNV) in 723 rams from eight sheep breeds, ranging from 17 to 514 copies, with a median copy number of 188. Pearson correlation analysis showed that the CNV of ZNF280BY was negatively correlated with testis size in Hu sheep. Furthermore, its mRNA expression level in testis had no significant correlation with the CNV but was significantly correlated with testis size. This study concluded that the expression of ZNF280BY was closely related to testicular development, and the CNV of ZNF280BY could be used as an important genetic marker to evaluate the ram reproductive capacity at an early stage in Hu sheep.

Sex chromosome systems in Neotropical Primates: What have we learnt so far from cytogenetics and genomics?

Neotropical Primates (Platyrrhini) show great diversity in their life histories, ecology, behaviour and genetics. This diversity extends to their chromosome complements, both to autosomes and to sex chromosomes. In this contribution, we will review what is currently known about sex chromosomes in this group, both from cytogenetic and from genomic evidence. The X and Y chromosomes in Neotropical Primates, also known as New World Monkeys, have striking structural differences compared with Old World Monkeys when Catarrhini sex chromosomes are considered. The XY bivalent displays a different meiotic behaviour in prophase I, and their Y chromosome shows extensive genomic differences. Even though the most widespread sex chromosome system is the XX/XY and thus considered the ancestral one for Platyrrhini, modifications of this sexual system are observed within this group. Multiple sex chromosome systems originated from Y-autosome translocations were described in several genera (Aotus, Callimico and Alouatta). In the howler monkeys, genus Alouatta, an independent origin of the sexual systems in South American and Mesoamerican species was postulated. All the above-mentioned evidence suggests that the Y chromosome of Platyrrhini has a different evolutionary history compared with the Catarrhini Y. There is still much to understand regarding their sex chromosome systems.

Representing sex chromosomes in genome assemblies

Sex chromosomes have evolved hundreds of independent times across eukaryotes. As genome sequencing, assembly, and scaffolding techniques rapidly improve, it is now feasible to build fully phased sex chromosome assemblies. Despite technological advances enabling phased assembly of whole chromosomes, there are currently no standards for representing sex chromosomes when publicly releasing a genome. Furthermore, most computational analysis tools are unable to efficiently investigate their unique biology relative to autosomes. We discuss a diversity of sex chromosome systems and consider the challenges of representing sex chromosome pairs in genome assemblies. By addressing these issues now as technologies for full phasing of chromosomal assemblies are maturing, we can collectively ensure that future genome analysis toolkits can be broadly applied to all eukaryotes with diverse types of sex chromosome systems. Here we provide best practice guidelines for presenting a genome assembly that contains sex chromosomes. These guidelines can also be applied to other non-recombining genomic regions, such as S-loci in plants and mating-type loci in fungi and algae.

Evolution of the sex-determining region in Ginkgo biloba

Sex chromosomes or sex-determining regions (SDR) have been discovered in many dioecious plant species, including the iconic ‘living fossil' Ginkgo biloba, though the location and size of the SDR in G. biloba remain contradictory. Here we resolve these controversies and analyse the evolution of the SDR in this species. Based on transcriptome sequencing data from four genetic crosses we reconstruct male- and female-specific genetic maps and locate the SDR to the middle of chromosome 2. Integration of the genetic maps with the genome sequence reveals that recombination in and around the SDR is suppressed in a region of about 50 Mb in both males and females. However, occasional recombination does occur except a small, less than 5 Mb long region that does not recombine in males. Based on synonymous divergence between homologous X- and Y-linked genes in this region, we infer that the Ginkgo SDR is fairly old—at least of Cretaceous origin. The analysis of substitution rates and gene expression reveals only slight Y-degeneration. These results are consistent with findings in other dioecious plants with homomorphic sex chromosomes, where the SDR is typically small and evolves in a region with pre-existing reduced recombination, surrounded by long actively recombining pseudoautosomal regions.

This article is part of the theme issue ‘Sex determination and sex chromosome evolution in land plants’.

State of the art on the physical mapping of the Y-chromosome in the <i>Bovidae</i> and comparison with other species

The next generation sequencing has significantly contributed to clarify the genome structure of many species of zootechnical interest. However, to date, some portions of the genome, especially those linked to a heterogametic nature such as the Y chromosome, are difficult to assemble and many gaps are still present. It is well known that the fluorescence in situ hybridization (FISH) is an excellent tool for identifying genes unequivocably mapped on chromosomes. Therefore, FISH can contribute to the localization of unplaced genome sequences, as well as to correct assembly errors generated by comparative bioinformatics. To this end, it is necessary to have starting points; therefore, in this study, we reviewed the physically mapped genes on the Y chromosome of cattle, buffalo, sheep, goats, pigs, horses and alpacas. A total of 208 loci were currently mapped by FISH. 89 were located in the male-specific region of the Y chromosome (MSY) and 119 were identified in the pseudoautosomal region (PAR). The loci reported in MSY and PAR were respectively: 18 and 25 in Bos taurus, 5 and 7 in Bubalus bubalis, 5 and 24 in Ovis aries, 5 and 19 in Capra hircus, 10 and 16 in Sus scrofa, 46 and 18 in Equus caballus. While in Vicugna pacos only 10 loci are reported in the PAR region. The correct knowledge and assembly of all genome sequences, including those of genes mapped on the Y chromosome, will help to elucidate their biological processes, as well as to discover and exploit potentially epistasis effects useful for selection breeding programs.

Decoding sex: Elucidating sex determination and how high-quality genome assemblies are untangling the evolutionary dynamics of sex chromosomes

Sexual reproduction is a diverse and widespread process. In gonochoristic species, the differentiation of sexes occurs through diverse mechanisms, influenced by environmental and genetic factors. In most vertebrates, a master-switch gene is responsible for triggering a sex determination network. However, only a few genes have acquired master-switch functions, and this process is associated with the evolution of sex-chromosomes, which have a significant influence in evolution. Additionally, their highly repetitive regions impose challenges for high-quality sequencing, even using high-throughput, state-of-the-art techniques. Here, we review the mechanisms involved in sex determination and their role in the evolution of species, particularly vertebrates, focusing on sex chromosomes and the challenges involved in sequencing these genomic elements. We also address the improvements provided by the growth of sequencing projects, by generating a massive number of near-gapless, telomere-to-telomere, chromosome-level, phased assemblies, increasing the number and quality of sex-chromosome sequences available for further studies.

Unique structure and positive selection promote the rapid divergence of Drosophila Y chromosomes

Y chromosomes across diverse species convergently evolve a gene-poor, heterochromatic organization enriched for duplicated genes, LTR retrotransposons, and satellite DNA. Sexual antagonism and a loss of recombination play major roles in the degeneration of young Y chromosomes. However, the processes shaping the evolution of mature, already degenerated Y chromosomes are less well-understood. Because Y chromosomes evolve rapidly, comparisons between closely related species are particularly useful. We generated de novo long-read assemblies complemented with cytological validation to reveal Y chromosome organization in three closely related species of the Drosophila simulans complex, which diverged only 250,000 years ago and share >98% sequence identity. We find these Y chromosomes are divergent in their organization and repetitive DNA composition and discover new Y-linked gene families whose evolution is driven by both positive selection and gene conversion. These Y chromosomes are also enriched for large deletions, suggesting that the repair of double-strand breaks on Y chromosomes may be biased toward microhomology-mediated end joining over canonical non-homologous end-joining. We propose that this repair mechanism contributes to the convergent evolution of Y chromosome organization across organisms.

Protein innovation through template switching in the Saccharomyces cerevisiae lineage

DNA polymerase template switching between short, non-identical inverted repeats (IRs) is a genetic mechanism that leads to the homogenization of IR arms and to IR spacer inversion, which cause multinucleotide mutations (MNMs). It is unknown if and how template switching affects gene evolution. In this study, we performed a phylogenetic analysis to determine the effect of template switching between IR arms on coding DNA of Saccharomyces cerevisiae. To achieve this, perfect IRs that co-occurred with MNMs between a strain and its parental node were identified in S. cerevisiae strains. We determined that template switching introduced MNMs into 39 protein-coding genes through S. cerevisiae evolution, resulting in both arm homogenization and inversion of the IR spacer. These events in turn resulted in nonsynonymous substitutions and up to five neighboring amino acid replacements in a single gene. The study demonstrates that template switching is a powerful generator of multiple substitutions within codons. Additionally, some template switching events occurred more than once during S. cerevisiae evolution. Our findings suggest that template switching constitutes a general mutagenic mechanism that results in both nonsynonymous substitutions and parallel evolution, which are traditionally considered as evidence for positive selection, without the need for adaptive explanations.

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.

Muller's ratchet of the Y chromosome with gene conversion

Muller's ratchet is a process in which deleterious mutations are fixed irreversibly in the absence of recombination. The degeneration of the Y chromosome, and the gradual loss of its genes, can be explained by Muller's ratchet. However, most theories consider single-copy genes, and may not be applicable to Y chromosomes, which have a number of duplicated genes in many species, which are probably undergoing concerted evolution by gene conversion. We developed a model of Muller's ratchet to explore the evolution of the Y chromosome. The model assumes a non-recombining chromosome with both single-copy and duplicated genes. We used analytical and simulation approaches to obtain the rate of gene loss in this model, with special attention to the role of gene conversion. Homogenization by gene conversion makes both duplicated copies either mutated or intact. The former promotes the ratchet, and the latter retards, and we ask which of these counteracting forces dominates under which conditions. We found that the effect of gene conversion is complex, and depends upon the fitness effect of gene duplication. When duplication has no effect on fitness, gene conversion accelerates the ratchet of both single-copy and duplicated genes. If duplication has an additive fitness effect, the ratchet of single-copy genes is accelerated by gene duplication, regardless of the gene conversion rate, whereas gene conversion slows the degeneration of duplicated genes. Our results suggest that the evolution of the Y chromosome involves several parameters, including the fitness effect of gene duplication by increasing dosage and gene conversion rate.

Genetics of Male Infertility - Present and Future: A Narrative Review

Infertility affects 8%–12% of couples worldwide with a male factor contributing to nearly 50% of couples either as a primary or contributing cause. Several genetic factors that include single-gene and multiple-gene defects associated with male infertility were reported in the past two decades. However, the etiology remains ambiguous in a majority of infertile men (~40%). The objective of this narrative review is to provide an update on the genetic factors associated with idiopathic male infertility and male reproductive system abnormalities identified in the last two decades. We performed a thorough literature search in online databases from January 2000 to July 2021. We observed a total of 13 genes associated with nonobstructive azoospermia due to maturation/meiotic arrest. Several studies that reported novel genes associated with multiple morphological abnormalities of the sperm flagella are also discussed in this review. ADGRG2, PANK2, SCNN1B, and CA12 genes are observed in non-CFTR-related vas aplasia. The genomic analysis should be quickly implemented in clinical practice as the detection of gene abnormalities in different male infertility phenotypes will facilitate genetic counseling.

GC-biased gene conversion in X-chromosome palindromes conserved in human, chimpanzee, and rhesus macaque

Gene conversion is GC-biased across a wide range of taxa. Large palindromes on mammalian sex chromosomes undergo frequent gene conversion that maintains arm-to-arm sequence identity greater than 99%, which may increase their susceptibility to the effects of GC-biased gene conversion. Here, we demonstrate a striking history of GC-biased gene conversion in 12 palindromes conserved on the X chromosomes of human, chimpanzee, and rhesus macaque. Primate X-chromosome palindrome arms have significantly higher GC content than flanking single-copy sequences. Nucleotide replacements that occurred in human and chimpanzee palindrome arms over the past 7 million years are one-and-a-half times as GC-rich as the ancestral bases they replaced. Using simulations, we show that our observed pattern of nucleotide replacements is consistent with GC-biased gene conversion with a magnitude of 70%, similar to previously reported values based on analyses of human meioses. However, GC-biased gene conversion since the divergence of human and rhesus macaque explains only a fraction of the observed difference in GC content between palindrome arms and flanking sequence, suggesting that palindromes are older than 29 million years and/or had elevated GC content at the time of their formation. This work supports a greater than 2:1 preference for GC bases over AT bases during gene conversion and demonstrates that the evolution and composition of mammalian sex chromosome palindromes is strongly influenced by GC-biased gene conversion.

Y chromosome functions in mammalian spermatogenesis

The mammalian Y chromosome is critical for male sex determination and spermatogenesis. However, linking each Y gene to specific aspects of male reproduction has been challenging. As the Y chromosome is notoriously hard to sequence and target, functional studies have mostly relied on transgene-rescue approaches using mouse models with large multi-gene deletions. These experimental limitations have oriented the field toward the search for a minimum set of Y genes necessary for male reproduction. Here, considering Y-chromosome evolutionary history and decades of discoveries, we review the current state of research on its function in spermatogenesis and reassess the view that many Y genes are disposable for male reproduction.

Optical genome mapping enables constitutional chromosomal aberration detection

Chromosomal aberrations including structural variations (SVs) are a major cause of human genetic diseases. Their detection in clinical routine still relies on standard cytogenetics. Drawbacks of these tests are a very low resolution (karyotyping) and the inability to detect balanced SVs or indicate the genomic localization and orientation of duplicated segments or insertions (copy number variant [CNV] microarrays). Here, we investigated the ability of optical genome mapping (OGM) to detect known constitutional chromosomal aberrations. Ultra-high-molecular-weight DNA was isolated from 85 blood or cultured cells and processed via OGM. A de novo genome assembly was performed followed by structural variant and CNV calling and annotation, and results were compared to known aberrations from standard-of-care tests (karyotype, FISH, and/or CNV microarray). In total, we analyzed 99 chromosomal aberrations, including seven aneuploidies, 19 deletions, 20 duplications, 34 translocations, six inversions, two insertions, six isochromosomes, one ring chromosome, and four complex rearrangements. Several of these variants encompass complex regions of the human genome involved in repeat-mediated microdeletion/microduplication syndromes. High-resolution OGM reached 100% concordance compared to standard assays for all aberrations with non-centromeric breakpoints. This proof-of-principle study demonstrates the ability of OGM to detect nearly all types of chromosomal aberrations. We also suggest suited filtering strategies to prioritize clinically relevant aberrations and discuss future improvements. These results highlight the potential for OGM to provide a cost-effective and easy-to-use alternative that would allow comprehensive detection of chromosomal aberrations and structural variants, which could give rise to an era of "next-generation cytogenetics."

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.

Complex Evolutionary History of the Y Chromosome in Flies of the Drosophila obscura Species Group

The Drosophila obscura species group shows dramatic variation in karyotype, including transitions among sex chromosomes. Members of the affinis and pseudoobscura subgroups contain a neo-X chromosome (a fusion of the X with an autosome), and ancestral Y genes have become autosomal in species harboring the neo-X. Detailed analysis of species in the pseudoobscura subgroup revealed that ancestral Y genes became autosomal through a translocation to the small dot chromosome. Here, we show that the Y-dot translocation is restricted to the pseudoobscura subgroup, and translocation of ancestral Y genes in the affinis subgroup likely followed a different route. We find that most ancestral Y genes have translocated to unique autosomal or X-linked locations in different taxa of the affinis subgroup, and we propose a dynamic model of sex chromosome formation and turnover in the obscura species group. Our results suggest that Y genes can find unique paths to escape unfavorable genomic environments that form after sex chromosome–autosome fusions.

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.

Panorama of stepwise involvement of the IgH 3' regulatory region in murine B cells

Among the multiple events leading to immunoglobulin (Ig) expression in B cells, stepwise activation of the Ig heavy chain locus (IgH) is of critical importance. Transcription regulation of the complex IgH locus has always been an interesting viewpoint to unravel the multiple and complex events required for IgH expression. First, regulatory germline transcripts (GLT) assist DNA remodeling events such as VDJ recombination, class switch recombination (CSR) and somatic hypermutation (SHM). Second, productive spliced transcripts restrict heavy chain protein expression associated either with the surface receptor of developing B cells or secreted in large amounts in plasma cells. One main transcriptional regulator for IgH lies at its 3' extremity and includes both a set of enhancers grouped in a large 3' regulatory region (3'RR) and a cluster of 3'CTCF-binding elements (3'CBEs). In this focused review, we will preferentially refer to evidence reported for the murine endogenous IgH locus, whether it is wt or carries deletions or insertions within the IgH 3' boundary and associated regulatory region.

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.

Were eukaryotes made by sex?: Sex might have been vital for merging endosymbiont and host genomes giving rise to eukaryotes

I hypothesize that the appearance of sex facilitated the merging of the endosymbiont and host genomes during early eukaryote evolution. Eukaryotes were formed by symbiosis between a bacterium that entered an archaeon, eventually giving rise to mitochondria. This entry was followed by the gradual transfer of most bacterial endosymbiont genes into the archaeal host genome. I argue that the merging of the mitochondrial genes into the host genome was vital for the evolution of genuine eukaryotes. At the time this process commenced it was unprecedented and required a novel mechanism. I suggest that this mechanism was meiotic sex, and that its appearance might have been THE crucial step that enabled the evolution of proper eukaryotes from early endosymbiont containing proto-eukaryotes. Sex might continue to be essential today for keeping genome insertions in check. Also see the video abstract here: https://youtu.be/aVMvWMpomac.

Large X-Linked Palindromes Undergo Arm-to-Arm Gene Conversion across Mus Lineages

Large (>10 kb), nearly identical (>99% nucleotide identity), palindromic sequences are enriched on mammalian sex chromosomes. Primate Y-palindromes undergo high rates of arm-to-arm gene conversion, a proposed mechanism for maintaining their sequence integrity in the absence of X-Y recombination. It is unclear whether X-palindromes, which can freely recombine in females, undergo arm-to-arm gene conversion and, if so, at what rate. We generated high-quality sequence assemblies of Mus molossinus and M. spretus X-palindromic regions and compared them with orthologous M. musculus X-palindromes. Our evolutionary sequence comparisons find evidence of X-palindrome arm-to-arm gene conversion at rates comparable to autosomal allelic gene conversion rates in mice. Mus X-palindromes also carry more derived than ancestral variants between species, suggesting that their sequence is rapidly diverging. We speculate that in addition to maintaining genes' sequence integrity via sequence homogenization, palindrome arm-to-arm gene conversion may also facilitate rapid sequence divergence.

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.

Multi-Copy Gene Family Evolution on the Avian W Chromosome

The sex chromosomes often follow unusual evolutionary trajectories. In particular, the sex-limited chromosomes frequently exhibit a small but unusual gene content in numerous species, where many genes have undergone massive gene amplification. The reasons for this remain elusive with a number of recent studies implicating meiotic drive, sperm competition, genetic drift, and gene conversion in the expansion of gene families. However, our understanding is primarily based on Y chromosome studies as few studies have systematically tested for copy number variation on W chromosomes. Here, we conduct a comprehensive investigation into the abundance, variability, and evolution of ampliconic genes on the avian W. First, we quantified gene copy number and variability across the duck W chromosome. We find a limited number of gene families as well as conservation in W-linked gene copy number across duck breeds, indicating that gene amplification may not be such a general feature of sex chromosome evolution as Y studies would initially suggest. Next, we investigated the evolution of HINTW, a prominent ampliconic gene family hypothesized to play a role in female reproduction and oogenesis. In particular, we investigated the factors driving the expansion of HINTW using contrasts between modern chicken and duck breeds selected for different female-specific selection regimes and their wild ancestors. Although we find the potential for selection related to fecundity in explaining small-scale gene amplification of HINTW in the chicken, purifying selection seems to be the dominant mode of evolution in the duck. Together, this challenges the assumption that HINTW is key for female fecundity across the avian phylogeny.

Genomic Features of Parthenogenetic Animals

Evolution without sex is predicted to impact genomes in numerous ways. Case studies of individual parthenogenetic animals have reported peculiar genomic features that were suggested to be caused by their mode of reproduction, including high heterozygosity, a high abundance of horizontally acquired genes, a low transposable element load, or the presence of palindromes. We systematically characterized these genomic features in published genomes of 26 parthenogenetic animals representing at least 18 independent transitions to asexuality. Surprisingly, not a single feature was systematically replicated across a majority of these transitions, suggesting that previously reported patterns were lineage-specific rather than illustrating the general consequences of parthenogenesis. We found that only parthenogens of hybrid origin were characterized by high heterozygosity levels. Parthenogens that were not of hybrid origin appeared to be largely homozygous, independent of the cellular mechanism underlying parthenogenesis. Overall, despite the importance of recombination rate variation for the evolution of sexual animal genomes, the genome-wide absence of recombination does not appear to have had the dramatic effects which are expected from classical theoretical models. The reasons for this are probably a combination of lineage-specific patterns, the impact of the origin of parthenogenesis, and a survivorship bias of parthenogenetic lineages.

Y chromosome copy number variation and its effects on fertility and other health factors: a review

The Y chromosome is essential for testis development and spermatogenesis. It is a chromosome with the lowest gene density owing to its medium size but paucity of coding genes. The Y chromosome is unique in that the majority of its structure is highly repetitive sequences, with the majority of these limited genes occurring in 9 amplionic sequences throughout the chromosome. The repetitive nature has its benefits as it can be protective against gene loss over many generations, but it can also predispose the Y chromosome to having wide variations of the number of gene copies present in these repeated sequences. This is known as copy number variation. Copy number variation is not unique to the Y chromosome but copy number variation is a well-known cause of male infertility and having effects on spermatogenesis. This is most commonly seen as deletions of the AZF sequences on the Y chromosome. However, there are other implications for copy number variation beyond just the AZF deletions that can affect spermatogenesis and potentially have other health implications. Copy number variations of TSPY1, DAZ, CDY1, RBMY1, the DYZ1 array, along with minor deletions of gr/gr, b1/b3, and b2/b3 have all be implicated in affecting spermatogenesis. UTY copy number variations have been implicated in risk for cardiovascular disease, and other deletions within gr/gr and the AZF sequences have been implicated in cancer and neuropsychiatric diseases. This review sets out to describe the Y chromosome and unique susceptibility to copy number variation and then to examine how this growing body of research impacts spermatogenesis and other health factors.

IUPACpal: efficient identification of inverted repeats in IUPAC-encoded DNA sequences

BACKGROUND

An inverted repeat is a DNA sequence followed downstream by its reverse complement, potentially with a gap in the centre. Inverted repeats are found in both prokaryotic and eukaryotic genomes and they have been linked with countless possible functions. Many international consortia provide a comprehensive description of common genetic variation making alternative sequence representations, such as IUPAC encoding, necessary for leveraging the full potential of such broad variation datasets.

RESULTS

We present IUPACPAL, an exact tool for efficient identification of inverted repeats in IUPAC-encoded DNA sequences allowing also for potential mismatches and gaps in the inverted repeats.

CONCLUSION

Within the parameters that were tested, our experimental results show that IUPACPAL compares favourably to a similar application packaged with EMBOSS. We show that IUPACPAL identifies many previously unidentified inverted repeats when compared with EMBOSS, and that this is also performed with orders of magnitude improved speed.

A new duck genome reveals conserved and convergently evolved chromosome architectures of birds and mammals

BACKGROUND

Ducks have a typical avian karyotype that consists of macro- and microchromosomes, but a pair of much less differentiated ZW sex chromosomes compared to chickens. To elucidate the evolution of chromosome architectures between ducks and chickens, and between birds and mammals, we produced a nearly complete chromosomal assembly of a female Pekin duck by combining long-read sequencing and multiplatform scaffolding techniques.

RESULTS

A major improvement of genome assembly and annotation quality resulted from the successful resolution of lineage-specific propagated repeats that fragmented the previous Illumina-based assembly. We found that the duck topologically associated domains (TAD) are demarcated by putative binding sites of the insulator protein CTCF, housekeeping genes, or transitions of active/inactive chromatin compartments, indicating conserved mechanisms of spatial chromosome folding with mammals. There are extensive overlaps of TAD boundaries between duck and chicken, and also between the TAD boundaries and chromosome inversion breakpoints. This suggests strong natural selection pressure on maintaining regulatory domain integrity, or vulnerability of TAD boundaries to DNA double-strand breaks. The duck W chromosome retains 2.5-fold more genes relative to chicken. Similar to the independently evolved human Y chromosome, the duck W evolved massive dispersed palindromic structures, and a pattern of sequence divergence with the Z chromosome that reflects stepwise suppression of homologous recombination.

CONCLUSIONS

Our results provide novel insights into the conserved and convergently evolved chromosome features of birds and mammals, and also importantly add to the genomic resources for poultry studies.

Population data and genetic diversity analysis of 17 Y-STR loci in Saudi population

Background

The Y chromosome polymorphism has been widely studied for human migrations, population genetics, forensic applications, and paternity analysis. However, studies regarding genetic lineage and population genetic structure of the Y chromosome in different regions of Saudi Arabia are limited.

Aim

This study aimed to analyze the distribution of Y chromosome haplotypes in a sample of 125 native Saudi males from different geographic regions of Saudi Arabia and compare to previously published Y chromosome haplotype data from Saudi Arabia and some neighboring Arab populations.

Materials and methods

Buccal swabs were collected from 125 healthy unrelated native Saudi males from different geographic regions of Saudi Arabia. Genomic DNA was extracted by Chelex®100; 17 Y-STR loci were amplified using the AmpFℓlSTR Yfiler PCR amplification kit and detected on the 3130 Genetic AnalyzerTM. Allele frequency and gene diversity were calculated with online tool STRAF. The Saudi population data were compared with the neighboring populations using pairwise genetic distances and associated probability values were calculated using the Y Chromosome Haplotype Reference Database Website (YHRD) software.

Results and conclusion

One hundred six YSTR haplotypes and 102 YSTR alleles (excluding 4 null alleles) were identified having a discrimination capacity (DC) of 85.8%. The highest haplotype diversity (HD) and gene diversity (GD) were observed at the loci DYS 458 (0.817) and DYS385b (0.807), respectively. According to our results, the Iraqi and Qena (Egypt) populations appeared to have closer relatedness to the Saudi population as compared with Yemen. The UAE and Kuwait populations showed the same degree of relatedness to the Saudi population followed by Bahrain. On the contrary, the Adnanit and Qahtanit populations of Jordan demonstrated low genetic distance from the Saudi population. In short, studying a population sample of pure Saudi ethnicity enabled us to identify a unique set of haplotypes which may help in establishing genetic relatedness between Saudi and the neighboring Arab populations. The present paper, therefore, highlights the importance of ensuring ethnic originality of the study sample while conducting population genetics studies.

Sequence analysis in Bos taurus reveals pervasiveness of X-Y arms races in mammalian lineages

Studies of Y Chromosome evolution have focused primarily on gene decay, a consequence of suppression of crossing-over with the X Chromosome. Here, we provide evidence that suppression of X-Y crossing-over unleashed a second dynamic: selfish X-Y arms races that reshaped the sex chromosomes in mammals as different as cattle, mice, and men. Using super-resolution sequencing, we explore the Y Chromosome of Bos taurus (bull) and find it to be dominated by massive, lineage-specific amplification of testis-expressed gene families, making it the most gene-dense Y Chromosome sequenced to date. As in mice, an X-linked homolog of a bull Y-amplified gene has become testis-specific and amplified. This evolutionary convergence implies that lineage-specific X-Y coevolution through gene amplification, and the selfish forces underlying this phenomenon, were dominatingly powerful among diverse mammalian lineages. Together with Y gene decay, X-Y arms races molded mammalian sex chromosomes and influenced the course of mammalian evolution.

Comparative analysis reveals the expansion of mitochondrial DNA control region containing unusually high G-C tandem repeat arrays in Nasonia vitripennis

Insect mitochondrial DNA (mtDNA) ranges from 14 to 19 kbp, and the size difference is attributed to the AT-rich control region. Jewel wasps have a parasitoid lifestyle, which may affect mitochondria function and evolution. We sequenced, assembled, and annotated mitochondrial genomes in Nasonia and outgroup species. Gene composition and order are conserved within Nasonia, but they differ from other parasitoids by two large inversion events that were not reported before. We observed a much higher substitution rate relative to the nuclear genome and mitochondrial introgression between N. giraulti and N. oneida, which is consistent with previous studies. Most strikingly, N. vitripennis mtDNA has an extremely long control region (7665 bp), containing twenty-nine 217 bp tandem repeats and can fold into a super-cruciform structure. In contrast to tandem repeats commonly found in other mitochondria, these high-copy repeats are highly conserved (98.7% sequence identity), much longer in length (approximately 8 Kb), extremely GC-rich (50.7%), and CpG-rich (percent CpG 19.4% vs. 1.1% in coding region), resulting in a 23 kbp mtDNA beyond the typical size range in insects. These N. vitripennis-specific mitochondrial repeats are not related to any known sequences in insect mitochondria. Their evolutionary origin and functional consequences warrant further investigations.

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.

Sequence Transpositions Restore Genes on the Highly Degenerated W Chromosomes of Songbirds

The female-specific W chromosomes of most Neognathae birds are highly degenerated and gene-poor. Previous studies have demonstrated that the gene repertoires of the Neognathae bird W chromosomes, despite being in small numbers, are conserved across bird species, likely due to purifying selection maintaining the regulatory and dosage-sensitive genes. Here we report the discovery of DNA-based sequence duplications from the Z to the W chromosome in birds-of-paradise (Paradisaeidae, Passeriformes), through sequence transposition. The original transposition involved nine genes, but only two of them (ANXA1 and ALDH1A1) survived on the W chromosomes. Both ANXA1 and ALDH1A1 are predicted to be dosage-sensitive, and the expression of ANXA1 is restricted to ovaries in all the investigated birds. These analyses suggest the newly transposed gene onto the W chromosomes can be favored for their role in restoring dosage imbalance or through female-specific selection. After examining seven additional songbird genomes, we further identified five other transposed genes on the W chromosomes of Darwin's finches and one in the great tit, expanding the observation of the Z-to-W transpositions to a larger range of bird species, but not all transposed genes exhibit dosage-sensitivity or ovary-biased expression We demonstrate a new mechanism by which the highly degenerated W chromosomes of songbirds can acquire genes from the homologous Z chromosomes, but further functional investigations are needed to validate the evolutionary forces underlying the transpositions.

Dynamic evolution of great ape Y chromosomes

The male-specific Y chromosome harbors genes important for sperm production. Because Y is repetitive, its DNA sequence was deciphered for only a few species, and its evolution remains elusive. Here we compared the Y chromosomes of great apes (human, chimpanzee, bonobo, gorilla, and orangutan) and found that many of their repetitive sequences and multicopy genes were likely already present in their common ancestor. Y repeats had increased intrachromosomal contacts, which might facilitate preservation of genes and gene regulatory elements. Chimpanzee and bonobo, experiencing high sperm competition, underwent many DNA changes and gene losses on the Y. Our research is significant for understanding the role of the Y chromosome in reproduction of nonhuman great apes, all of which are endangered.

The mammalian male-specific Y chromosome plays a critical role in sex determination and male fertility. However, because of its repetitive and haploid nature, it is frequently absent from genome assemblies and remains enigmatic. The Y chromosomes of great apes represent a particular puzzle: their gene content is more similar between human and gorilla than between human and chimpanzee, even though human and chimpanzee share a more recent common ancestor. To solve this puzzle, here we constructed a dataset including Ys from all extant great ape genera. We generated assemblies of bonobo and orangutan Ys from short and long sequencing reads and aligned them with the publicly available human, chimpanzee, and gorilla Y assemblies. Analyzing this dataset, we found that the genus Pan, which includes chimpanzee and bonobo, experienced accelerated substitution rates. Pan also exhibited elevated gene death rates. These observations are consistent with high levels of sperm competition in Pan. Furthermore, we inferred that the great ape common ancestor already possessed multicopy sequences homologous to most human and chimpanzee palindromes. Nonetheless, each species also acquired distinct ampliconic sequences. We also detected increased chromatin contacts between and within palindromes (from Hi-C data), likely facilitating gene conversion and structural rearrangements. Our results highlight the dynamic mode of Y chromosome evolution and open avenues for studies of male-specific dispersal in endangered great ape species.