Characterization of Prdm9 in equids and sterility in mules

Fertility of hybrids of dromedary and Bactrian camels: A possible role of conserved architecture of zinc finger domain of recombination regulator PRDM9

Recombination regulator, PRDM9, has been regarded as the most rapidly evolving gene in the genomes of many metazoans, in addition to being acknowledged as the sole speciation gene in vertebrates. It has become the focus of many scientific investigations because of exceptional numerical and sequence variability in its zinc finger (ZF) domain within and across species that contributes to reproductive isolation between species. This study is the maiden attempt to explore the architecture of PRDM9 ZF domain in two Camelid species (Camelus dromedarius and Camelus bactrianus). Sequence analysis revealed highly conserved domain architecture with presence of 3 and 4 ZFs in dromedary and Bactrian camels, respectively. Typical evolutionary features of PRDM9 ZF domain i.e. concerted evolution and positive selection were invariably absent in both the one-humped dromedary and the two-humped Bactrian camels. Fertility of hybrids of dromedary and Bactrian camels, despite being taxonomically distinct species can be attributed to the lack of sequence variability in PRDM9 in these species. Phylogenetic analysis underpinned clear demarcation of camels from other livestock species. The results of the present study defy what has been learnt so far about PRDM9 and add to the enigma surrounding the most intriguing gene in the genome.

Expression Analysis of Molecular Chaperones Hsp70 and Hsp90 on Development and Metabolism of Different Organs and Testis in Cattle (Cattle-yak and Yak)

Hsp70 and Hsp90 play an important role in testis development and spermatogenesis regulation, but the exact connection between Hsp70 and Hsp90 and metabolic stress in cattle is unclear. Here, we focused on the male cattle−yak and yak, investigated the expression and localization of Hsp70 and Hsp90 in their tissues, and explored the influence of these factors on development and metabolism. In our study, a total of 54 cattle (24 cattle−yaks and 30 yaks; aged 1 day to 10 years) were examined. The Hsp90 mRNA of the cattle−yak was first cloned and compared with that of the yak, and variation in the amino acid sequence was found, which led to differences in protein spatial structure. Using real-time quantitative PCR (RT-qPCR) and Western blot (WB) techniques, we investigated whether the expression of Hsp70 and Hsp90 mRNA and protein are different in the cattle−yak and yak. We found a disparity in Hsp70 and Hsp90 mRNA and protein expression in different non-reproductive organs and in testicular tissues at different stages of development, while high expression was observed in the testes of both juveniles and adults. Moreover, it was intriguing to observe that Hsp70 expression was significantly high in the yak, whereas Hsp90 was high in the cattle−yak (p < 0.01). We also examined the location of Hsp70 and Hsp90 in the testis by immunohistochemical (IHC) and immunofluorescence (IF) techniques, and the results showed that Hsp70 and Hsp90 were positive in the epithelial cells, spermatogenic cells, and mesenchymal cells. In summary, our study proved that Hsp70 and Hsp90 expressions were different in different tissues (kidney, heart, cerebellum, liver, lung, spleen, and testis), and Hsp90 expression was high in the testis of the cattle−yak, suggesting that dysplasia of the cattle−yak may correlate with an over-metabolism of Hsp90.

Evolution of the recombination regulator PRDM9 in minke whales

Background

PRDM9 is a key regulator of meiotic recombination in most metazoans, responsible for reshuffling parental genomes. During meiosis, the PRDM9 protein recognizes and binds specific target motifs via its array of C₂H₂ zinc-fingers encoded by a rapidly evolving minisatellite. The gene coding for PRDM9 is the only speciation gene identified in vertebrates to date and shows high variation, particularly in the DNA-recognizing positions of the zinc-finger array, within and between species. Across all vertebrate genomes studied for PRDM9 evolution, only one genome lacks variability between repeat types – that of the North Pacific minke whale. This study aims to understand the evolution and diversity of Prdm9 in minke whales, which display the most unusual genome reference allele of Prdm9 so far discovered in mammals.

Results

Minke whales possess all the features characteristic of PRDM9-directed recombination, including complete KRAB, SSXRD and SET domains and a rapidly evolving array of C₂H₂-type-Zincfingers (ZnF) with evidence of rapid evolution, particularly at DNA-recognizing positions that evolve under positive diversifying selection. Seventeen novel PRDM9 variants were identified within the Antarctic minke whale species, plus a single distinct PRDM9 variant in Common minke whales – shared across North Atlantic and North Pacific minke whale subspecies boundaries.

Conclusion

The PRDM9 ZnF array evolves rapidly, in minke whales, with at least one DNA-recognizing position under positive selection. Extensive PRDM9 diversity is observed, particularly in the Antarctic in minke whales. Common minke whales shared a specific Prdm9 allele across subspecies boundaries, suggesting incomplete speciation by the mechanisms associated with PRDM9 hybrid sterility.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08305-1.

Genetics, Evolution, and Physiology of Donkeys and Mules

The genus Equus is made up of donkeys, horses, and zebras. Despite significant variation in chromosome number across these species, interspecies breeding results in healthy, although infertile, hybrid offspring. Most notable among these are the horse-donkey hybrids, the mule and hinny. Donkeys presently are used for everything from companion animals to beasts of burden. Although closely related from an evolutionary standpoint, differences in anatomy and physiology preclude the assumption that they can be treated identically to the domestic horse. Veterinarians should be aware of these differences and adjust their practice accordingly.

Unravelling the hybrid vigor in domestic equids: the effect of hybridization on bone shape variation and covariation

BACKGROUND

Hybridization has been widely practiced in plant and animal breeding as a means to enhance the quality and fitness of the organisms. In domestic equids, this hybrid vigor takes the form of improved physical and physiological characteristics, notably for strength or endurance. Because the offspring of horse and donkey is generally sterile, this widely recognized vigor is expressed in the first generation (F1). However, in the absence of recombination between the two parental genomes, F1 hybrids can be expected to be phenotypically intermediate between their parents which could potentially restrict the possibilities of an increase in overall fitness. In this study, we examine the morphology of the main limb bones of domestic horses, donkeys and their hybrids to investigate the phenotypic impact of hybridization on the locomotor system. We explore bone shape variation and covariation to gain insights into the morphological and functional expressions of the hybrid vigor commonly described in domestic equids.

RESULTS

Our data reveal the occurrence of transgressive effects on several bones in the F1 generation. The patterns of morphological integration further demonstrate that the developmental processes producing covariation are not disrupted by hybridization, contrary to functional ones.

CONCLUSIONS

These results suggest that an increase in overall fitness could be related to more flexibility in shape change in hybrids, except for the main forelimb long bones of which the morphology is strongly driven by muscle interactions. More broadly, this study illustrates the interest of investigating not only bone shape variation but also underlying processes, in order to contribute to better understanding how developmental and functional mechanisms are affected by hybridization.

Exploration of fine-scale recombination rate variation in the domestic horse

Total genetic map length and local recombination landscapes typically vary within and across populations. As a first step to understanding the recombination landscape in the domestic horse, we calculated population recombination rates and identified likely recombination hotspots using approximately 1.8 million SNP genotypes for 485 horses from 32 distinct breeds. The resulting breed-averaged recombination map spans 2.36 Gb and accounts for 2939.07 cM. Recombination hotspots occur once per 23.8 Mb on average and account for ∼9% of the physical map length. Regions with elevated recombination rates in the entire cohort were enriched for genes in pathways involving interaction with the environment: immune system processes (specifically, MHC class I and class II genes), responses to stimuli, and serotonin receptor pathways. We found significant correlations between differences in local recombination rates and population differentiation quantified by F _ST Analysis of breed-specific maps revealed thousands of hotspot regions unique to particular breeds, as well as unique "coldspots," regions where a particular breed showed below-average recombination, whereas all other breeds had evidence of a hotspot. Finally, we identified relative enrichment (P = 5.88 × 10^-27) for the in silico-predicted recognition motif for equine PR/SET domain 9 (PRDM9) in recombination hotspots. These results indicate that selective pressures and PRDM9 function contribute to variation in recombination rates across the domestic horse genome.

Regulation by Hsp27/P53 in testis development and sperm apoptosis of male cattle (cattle-yak and yak)

Heat shock protein 27 (Hsp27)/protein 53 (P53) plays an important role in testis development and spermatozoa regulation, but the relationship between Hsp27/P53 and infertility in cattle is unclear. Here, we focus on male cattle-yak and yak to investigate the expression and localization of Hsp27/P53 in testis tissues and to explore the influence of Hsp27/P53 on infertility. In our study, a total of 54 cattle (24 cattle-yak and 30 yak) were examined. The Hsp27 and P53 messenger RNA (mRNA) of cattle-yak were cloned, and amino acid variations in Hsp27 and P53 were found; the variations led to differences in the protein spatial structure compared with yak. We used real-time quantitative polymerase chain reaction and western blot to investigate whether the expression of Hsp27/P53 mRNA and protein was different in cattle-yak and yak. We found that the expression levels of Hsp27/P53 mRNA and protein were different in the testis developmental stages and the highest expression was observed in testicles during adulthood. Moreover, the Hsp27 expression was significantly higher in yak, whereas P53 expression was higher in cattle-yak (p < 0.01). On this basis, we detected the location of Hsp27/P53 in the testis by immunohistochemistry and immunofluorescence. The results demonstrated that Hsp27 was located in spermatogenic cells at different developmental stages and mesenchymal cells of the yak testicles. However, P53 was located in the primary spermatocyte and interstitial cells of the cattle-yak testicles. In summary, our study proved that the expression of Hsp27/P53 differed across the testis developmental stages and the expression of P53 was higher in the testis of cattle-yak, which suggested that the infertility of cattle-yak may be caused by the upregulation of P53.

PRDM9, a driver of the genetic map

During meiosis, maternal and paternal chromosomes undergo exchanges by homologous recombination. This is essential for fertility and contributes to genome evolution. In many eukaryotes, sites of meiotic recombination, also called hotspots, are regions of accessible chromatin, but in many vertebrates, their location follows a distinct pattern and is specified by PR domain-containing protein 9 (PRDM9). The specification of meiotic recombination hotspots is achieved by the different activities of PRDM9: DNA binding, histone methyltransferase, and interaction with other proteins. Remarkably, PRDM9 activity leads to the erosion of its own binding sites and the rapid evolution of its DNA-binding domain. PRDM9 may also contribute to reproductive isolation, as it is involved in hybrid sterility potentially due to a reduction of its activity in specific heterozygous contexts.

Horse Y chromosome assembly displays unique evolutionary features and putative stallion fertility genes

Dynamic evolutionary processes and complex structure make the Y chromosome among the most diverse and least understood regions in mammalian genomes. Here, we present an annotated assembly of the male specific region of the horse Y chromosome (eMSY), representing the first comprehensive Y assembly in odd-toed ungulates. The eMSY comprises single-copy, equine specific multi-copy, PAR transposed, and novel ampliconic sequence classes. The eMSY gene density approaches that of autosomes with the highest number of retained X-Y gametologs recorded in eutherians, in addition to novel Y-born and transposed genes. Horse, donkey and mule testis RNAseq reveals several candidate genes for stallion fertility. A novel testis-expressed XY ampliconic sequence class, ETSTY7, is shared with the parasite Parascaris genome, providing evidence for eukaryotic horizontal transfer and inter-chromosomal mobility. Our study highlights the dynamic nature of the Y and provides a reference sequence for improved understanding of equine male development and fertility.

The consequences of sequence erosion in the evolution of recombination hotspots

Meiosis is initiated by a double-strand break (DSB) introduced in the DNA by a highly controlled process that is repaired by recombination. In many organisms, recombination occurs at specific and narrow regions of the genome, known as recombination hotspots, which overlap with regions enriched for DSBs. In recent years, it has been demonstrated that conversions and mutations resulting from the repair of DSBs lead to a rapid sequence evolution at recombination hotspots eroding target sites for DSBs. We still do not fully understand the effect of this erosion in the recombination activity, but evidence has shown that the binding of trans-acting factors like PRDM9 is affected. PRDM9 is a meiosis-specific, multi-domain protein that recognizes DNA target motifs by its zinc finger domain and directs DSBs to these target sites. Here we discuss the changes in affinity of PRDM9 to eroded recognition sequences, and explain how these changes in affinity of PRDM9 can affect recombination, leading sometimes to sterility in the context of hybrid crosses. We also present experimental data showing that DNA methylation reduces PRDM9 binding in vitro. Finally, we discuss PRDM9-independent hotspots, posing the question how these hotspots evolve and change with sequence erosion.

This article is part of the themed issue ‘Evolutionary causes and consequences of recombination rate variation in sexual organisms’.

Mammalian Comparative Genomics Reveals Genetic and Epigenetic Features Associated with Genome Reshuffling in Rodentia

Understanding how mammalian genomes have been reshuffled through structural changes is fundamental to the dynamics of its composition, evolutionary relationships between species and, in the long run, speciation. In this work, we reveal the evolutionary genomic landscape in Rodentia, the most diverse and speciose mammalian order, by whole-genome comparisons of six rodent species and six representative outgroup mammalian species. The reconstruction of the evolutionary breakpoint regions across rodent phylogeny shows an increased rate of genome reshuffling that is approximately two orders of magnitude greater than in other mammalian species here considered. We identified novel lineage and clade-specific breakpoint regions within Rodentia and analyzed their gene content, recombination rates and their relationship with constitutive lamina genomic associated domains, DNase I hypersensitivity sites and chromatin modifications. We detected an accumulation of protein-coding genes in evolutionary breakpoint regions, especially genes implicated in reproduction and pheromone detection and mating. Moreover, we found an association of the evolutionary breakpoint regions with active chromatin state landscapes, most probably related to gene enrichment. Our results have two important implications for understanding the mechanisms that govern and constrain mammalian genome evolution. The first is that the presence of genes related to species-specific phenotypes in evolutionary breakpoint regions reinforces the adaptive value of genome reshuffling. Second, that chromatin conformation, an aspect that has been often overlooked in comparative genomic studies, might play a role in modeling the genomic distribution of evolutionary breakpoints.

PRDM9 and Its Role in Genetic Recombination

PRDM9 is a zinc finger protein that binds DNA at specific locations in the genome where it trimethylates histone H3 at lysines 4 and 36 at surrounding nucleosomes. During meiosis in many species, including humans and mice where PRDM9 has been most intensely studied, these actions determine the location of recombination hotspots, where genetic recombination occurs. In addition, PRDM9 facilitates the association of hotspots with the chromosome axis, the site of the programmed DNA double-strand breaks (DSBs) that give rise to genetic exchange between chromosomes. In the absence of PRDM9 DSBs are not properly repaired. Collectively, these actions determine patterns of genetic linkage and the possibilities for chromosome reorganization over successive generations.

A new genus of horse from Pleistocene North America

The extinct ‘New World stilt-legged’, or NWSL, equids constitute a perplexing group of Pleistocene horses endemic to North America. Their slender distal limb bones resemble those of Asiatic asses, such as the Persian onager. Previous palaeogenetic studies, however, have suggested a closer relationship to caballine horses than to Asiatic asses. Here, we report complete mitochondrial and partial nuclear genomes from NWSL equids from across their geographic range. Although multiple NWSL equid species have been named, our palaeogenomic and morphometric analyses support the idea that there was only a single species of middle to late Pleistocene NWSL equid, and demonstrate that it falls outside of crown group Equus. We therefore propose a new genus, Haringtonhippus, for the sole species H. francisci. Our combined genomic and phenomic approach to resolving the systematics of extinct megafauna will allow for an improved understanding of the full extent of the terminal Pleistocene extinction event.

The horse family – which also includes zebras, donkeys and asses – is often featured on the pages of textbooks about evolution. All living horses belong to a group, or genus, called Equus. The fossil record shows how the ancestors of these animals evolved from dog-sized, three-toed browsers to larger, one-toed grazers. This process took around 55 million years, and many members of the horse family tree went extinct along the way.

Nevertheless, the details of the horse family tree over the past 2.5 million years remain poorly understood. In North America, horses from this period – which is referred to as the Pleistocene – have been classed into two major groups: stout-legged horses and stilt-legged horses. Both groups became extinct near the end of the Pleistocene in North America, and it was not clear how they relate to one another. Based on their anatomy, many scientists suggested that stilt-legged horses were most closely related to modern-day asses living in Asia. Yet, other studies using ancient DNA placed the stilt-legged horses closer to the stout-legged horses.

Heintzman et al. set out to resolve where the stilt-legged horses sit within the horse family tree by examining more ancient DNA than the previous studies. The analyses showed that the stilt-legged horses were much more distinct than previously thought. In fact, contrary to all previous findings, these animals actually belonged outside of the genus Equus. Heintzman et al. named the new genus for the stilt-legged horses Haringtonhippus, and showed that all stilt-legged horses belonged to a single species within this genus, Haringtonhippus francisci.

Together these new findings provide a benchmark for reclassifying problematic fossil groups across the tree of life. A similar approach could be used to resolve the relationships in other problematic groups of Pleistocene animals, such as mammoths and bison. This would give scientists a more nuanced understanding of evolution and extinction during this period.

Detecting taxonomic and phylogenetic signals in equid cheek teeth: towards new palaeontological and archaeological proxies

The Plio-Pleistocene evolution of Equus and the subsequent domestication of horses and donkeys remains poorly understood, due to the lack of phenotypic markers capable of tracing this evolutionary process in the palaeontological/archaeological record. Using images from 345 specimens, encompassing 15 extant taxa of equids, we quantified the occlusal enamel folding pattern in four mandibular cheek teeth with a single geometric morphometric protocol. We initially investigated the protocol accuracy by assigning each tooth to its correct anatomical position and taxonomic group. We then contrasted the phylogenetic signal present in each tooth shape with an exome-wide phylogeny from 10 extant equine species. We estimated the strength of the phylogenetic signal using a Brownian motion model of evolution with multivariate K statistic, and mapped the dental shape along the molecular phylogeny using an approach based on squared-change parsimony. We found clear evidence for the relevance of dental phenotypes to accurately discriminate all modern members of the genus Equus and capture their phylogenetic relationships. These results are valuable for both palaeontologists and zooarchaeologists exploring the spatial and temporal dynamics of the evolutionary history of the horse family, up to the latest domestication trajectories of horses and donkeys.

Ruminant-specific multiple duplication events of PRDM9 before speciation

Background

Understanding the genetic and evolutionary mechanisms of speciation genes in sexually reproducing organisms would provide important insights into mammalian reproduction and fitness. PRDM9, a widely known speciation gene, has recently gained attention for its important role in meiotic recombination and hybrid incompatibility. Despite the fact that PRDM9 is a key regulator of recombination and plays a dominant role in hybrid incompatibility, little is known about the underlying genetic and evolutionary mechanisms that generated multiple copies of PRDM9 in many metazoan lineages.

Results

The present study reports (1) evidence of ruminant-specific multiple gene duplication events, which likely have had occurred after the ancestral ruminant population diverged from its most recent common ancestor and before the ruminant speciation events, (2) presence of three copies of PRDM9, one copy (lineages I) in chromosome 1 (chr1) and two copies (lineages II & III) in chromosome X (chrX), thus indicating the possibility of ancient inter- and intra-chromosomal unequal crossing over and gene conversion events, (3) while lineages I and II are characterized by the presence of variable tandemly repeated C2H2 zinc finger (ZF) arrays, lineage III lost these arrays, and (4) C2H2 ZFs of lineages I and II, particularly the amino acid residues located at positions −1, 3, and 6 have evolved under strong positive selection.

Conclusions

Our results demonstrated two gene duplication events of PRDM9 in ruminants: an inter-chromosomal duplication that occurred between chr1 and chrX, and an intra-chromosomal X-linked duplication, which resulted in two additional copies of PRDM9 in ruminants. The observation of such duplication between chrX and chr1 is rare and may possibly have happened due to unequal crossing-over millions of years ago when sex chromosomes were independently derived from a pair of ancestral autosomes. Two copies (lineages I & II) are characterized by the presence of variable sized tandem-repeated C2H2 ZFs and evolved under strong positive selection and concerted evolution, supporting the notion of well-established Red Queen hypothesis. Collectively, gene duplication, concerted evolution, and positive selection are the likely driving forces for the expansion of ruminant PRDM9 sub-family.

Electronic supplementary material

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

Evidence of positive selection and concerted evolution in the rapidly evolving PRDM9 zinc finger domain in goats and sheep

Meiotic recombination contributes to augmentation of genetic diversity, exclusion of deleterious alleles and proper segregation of chromatids. PRDM9 has been identified as the gene responsible for specifying the location of recombination hotspots during meiosis and is also the only known vertebrate gene associated with reproductive isolation between species. PRDM9 encodes a protein with a highly variable zinc finger (ZF) domain that varies between as well as within species. In the present study, the ZF domain of PRDM9 on chromosome 1 was characterized for the first time in 15 goat breeds and 25 sheep breeds of India. A remarkable variation in the number and sequence of ZF domains was observed. The number of ZF repeats in the ZF array varied from eight to 12 yielding five homozygous and 10 heterozygous genotypes. The number of different ZF domains was 84 and 52 producing 36 and 26 unique alleles in goats and sheep respectively. The posterior mean of dN/dS or omega values were calculated using the codeml tool of pamlx to identify amino acids that are evolving positively in goats and sheep, as positions -1, +3 and +6 in the ZF domain have been reported to experience strong positive selection across different lineages. Our study identified sites -5, -1, +3, +4 and +6 to be experiencing positive selection. Small ruminant zinc fingers were also found to be evolving under concerted evolution. Our results demonstrate the existence of a vast diversity of PRDM9 in goats and sheep, which is in concert with reports in many metazoans.

The Meiotic Recombination Activator PRDM9 Trimethylates Both H3K36 and H3K4 at Recombination Hotspots In Vivo

In many mammals, including humans and mice, the zinc finger histone methyltransferase PRDM9 performs the first step in meiotic recombination by specifying the locations of hotspots, the sites of genetic recombination. PRDM9 binds to DNA at hotspots through its zinc finger domain and activates recombination by trimethylating histone H3K4 on adjacent nucleosomes through its PR/SET domain. Recently, the isolated PR/SET domain of PRDM9 was shown capable of also trimethylating H3K36 in vitro, raising the question of whether this reaction occurs in vivo during meiosis, and if so, what its function might be. Here, we show that full-length PRDM9 does trimethylate H3K36 in vivo in mouse spermatocytes. Levels of H3K4me3 and H3K36me3 are highly correlated at hotspots, but mutually exclusive elsewhere. In vitro, we find that although PRDM9 trimethylates H3K36 much more slowly than it does H3K4, PRDM9 is capable of placing both marks on the same histone molecules. In accord with these results, we also show that PRDM9 can trimethylate both K4 and K36 on the same nucleosomes in vivo, but the ratio of K4me3/K36me3 is much higher for the pair of nucleosomes adjacent to the PRDM9 binding site compared to the next pair further away. Importantly, H3K4me3/H3K36me3-double-positive nucleosomes occur only in regions of recombination: hotspots and the pseudoautosomal (PAR) region of the sex chromosomes. These double-positive nucleosomes are dramatically reduced when PRDM9 is absent, showing that this signature is PRDM9-dependent at hotspots; the residual double-positive nucleosomes most likely come from the PRDM9-independent PAR. These results, together with the fact that PRDM9 is the only known mammalian histone methyltransferase with both H3K4 and H3K36 trimethylation activity, suggest that trimethylation of H3K36 plays an important role in the recombination process. Given the known requirement of H3K36me3 for double strand break repair by homologous recombination in somatic cells, we suggest that it may play the same role in meiosis.

Genetic recombination is the meiotic process by which novel combinations of alleles are passed on to the next generation. This process accelerates evolution by creating genetic diversity, and is also essential for successful meiosis. In mammals, the enzyme PRDM9 initiates recombination and determines the subset of sites within the genome—called recombination hotspots—that recombine. PRDM9 does this by binding DNA at hotspots and placing epigenetic marks. Previously, PRDM9 was only known to place the H3K4me3 mark at hotspots in living cells. Here, we show that PRDM9 places the H3K36me3 mark at hotspots, and that H3K4me3 and H3K36me3 coincide only at regions of recombination in germ cells. We prove that this coincidence is driven almost entirely by PRDM9; there is dramatically less coincidence between H3K4me3 and H3K36me3 when PRDM9 is absent. These results reveal a new enzymatic function for PRDM9, and a new epigenetic signature at hotspots that may restrict recombination to these sites. Since aberrant recombination can cause aneuploidy resulting in fetal loss, and abnormal genome rearrangements that underlie many congenital syndromes and some cancers, a thorough understanding of this fundamental process has potentially far-reaching implications.

Zinc Finger Domain of the PRDM9 Gene on Chromosome 1 Exhibits High Diversity in Ruminants but Its Paralog PRDM7 Contains Multiple Disruptive Mutations

PRDM9 is the sole hybrid sterility gene identified so far in vertebrates. PRDM9 gene encodes a protein with an immensely variable zinc-finger (ZF) domain that determines the site of meiotic recombination hotspots genome-wide. In this study, the terminal ZF domain of PRDM9 on bovine chromosome 1 and its paralog on chromosome 22 were characterized in 225 samples from five ruminant species (cattle, yak, mithun, sheep and goat). We found extraordinary variation in the number of PRDM9 zinc fingers (6 to 12). We sequenced PRDM9 ZF encoding region from 15 individuals (carrying the same ZF number in both copies) and found 43 different ZF domain sequences. Ruminant zinc fingers of PRDM9 were found to be diversifying under positive selection and concerted evolution, specifically at positions involved in defining their DNA-binding specificity, consistent with the reports from other vertebrates such as mice, humans, equids and chimpanzees. ZF-encoding regions of the PRDM7, a paralog of PRDM9 on bovine chromosome 22 and on unknown chromosomes in other studied species were found to contain 84 base repeat units as in PRDM9, but there were multiple disruptive mutations after the first repeat unit. The diversity of the ZFs suggests that PRDM9 may activate recombination hotspots that are largely unique to each ruminant species.

Multimer Formation Explains Allelic Suppression of PRDM9 Recombination Hotspots

Genetic recombination during meiosis functions to increase genetic diversity, promotes elimination of deleterious alleles, and helps assure proper segregation of chromatids. Mammalian recombination events are concentrated at specialized sites, termed hotspots, whose locations are determined by PRDM9, a zinc finger DNA-binding histone methyltransferase. Prdm9 is highly polymorphic with most alleles activating their own set of hotspots. In populations exhibiting high frequencies of heterozygosity, questions remain about the influences different alleles have in heterozygous individuals where the two variant forms of PRDM9 typically do not activate equivalent populations of hotspots. We now find that, in addition to activating its own hotspots, the presence of one Prdm9 allele can modify the activity of hotspots activated by the other allele. PRDM9 function is also dosage sensitive; Prdm9^+/- heterozygous null mice have reduced numbers and less active hotspots and increased numbers of aberrant germ cells. In mice carrying two Prdm9 alleles, there is allelic competition; the stronger Prdm9 allele can partially or entirely suppress chromatin modification and recombination at hotspots of the weaker allele. In cell cultures, PRDM9 protein variants form functional heteromeric complexes which can bind hotspots sequences. When a heteromeric complex binds at a hotspot of one PRDM9 variant, the other PRDM9 variant, which would otherwise not bind, can still methylate hotspot nucleosomes. We propose that in heterozygous individuals the underlying molecular mechanism of allelic suppression results from formation of PRDM9 heteromers, where the DNA binding activity of one protein variant dominantly directs recombination initiation towards its own hotspots, effectively titrating down recombination by the other protein variant. In natural populations with many heterozygous individuals, allelic competition will influence the recombination landscape.

During formation of sperm and eggs chromosomes exchange DNA in a process known as recombination, creating new combinations responsible for much of the enormous diversity in populations. In some mammals, including humans, the locations of recombination are chosen by a DNA-binding protein named PRDM9. Importantly, there are tens to hundreds of different variations of the Prdm9 gene (termed alleles), many of which are predicted to bind a unique DNA sequence. This high frequency of variation results in many individuals having two different copies of Prdm9, and several lines of evidence indicate that alleles compete to initiate recombination. In seeking to understand the mechanism of this competition we found that Prdm9 activity is sensitive to the number of gene copies present, suggesting that availability of this protein is a limiting factor during recombination. Moreover, we found that variant forms of PRDM9 protein can physically interact suggesting that when this happens one variant can influence which hotspots will become activated. Genetic crosses in mice support these observations; the presence of a dominant Prdm9 allele can completely suppress recombination at some locations. We conclude that allele-dominance of PRDM9 is a consequence of protein-protein interaction and competition for DNA binding in a limited pool of molecules, thus shaping the recombination landscape in natural populations.

Affinity-seq detects genome-wide PRDM9 binding sites and reveals the impact of prior chromatin modifications on mammalian recombination hotspot usage

Background

Genetic recombination plays an important role in evolution, facilitating the creation of new, favorable combinations of alleles and the removal of deleterious mutations by unlinking them from surrounding sequences. In most mammals, the placement of genetic crossovers is determined by the binding of PRDM9, a highly polymorphic protein with a long zinc finger array, to its cognate binding sites. It is one of over 800 genes encoding proteins with zinc finger domains in the human genome.

Results

We report a novel technique, Affinity-seq, that for the first time identifies both the genome-wide binding sites of DNA-binding proteins and quantitates their relative affinities. We have applied this in vitro technique to PRDM9, the zinc-finger protein that activates genetic recombination, obtaining new information on the regulation of hotspots, whose locations and activities determine the recombination landscape. We identified 31,770 binding sites in the mouse genome for the PRDM9^Dom2 variant. Comparing these results with hotspot usage in vivo, we find that less than half of potential PRDM9 binding sites are utilized in vivo. We show that hotspot usage is increased in actively transcribed genes and decreased in genomic regions containing H3K9me2/3 histone marks or bound to the nuclear lamina.

Conclusions

These results show that a major factor determining whether a binding site will become an active hotspot and what its activity will be are constraints imposed by prior chromatin modifications on the ability of PRDM9 to bind to DNA in vivo. These constraints lead to the presence of long genomic regions depleted of recombination.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-015-0024-6) contains supplementary material, which is available to authorized users.

PRDM9 drives evolutionary erosion of hotspots in Mus musculus through haplotype-specific initiation of meiotic recombination

Meiotic recombination generates new genetic variation and assures the proper segregation of chromosomes in gametes. PRDM9, a zinc finger protein with histone methyltransferase activity, initiates meiotic recombination by binding DNA at recombination hotspots and directing the position of DNA double-strand breaks (DSB). The DSB repair mechanism suggests that hotspots should eventually self-destruct, yet genome-wide recombination levels remain constant, a conundrum known as the hotspot paradox. To test if PRDM9 drives this evolutionary erosion, we measured activity of the Prdm9^Cst allele in two Mus musculus subspecies, M.m. castaneus, in which Prdm9^Cst arose, and M.m. domesticus, into which Prdm9^Cst was introduced experimentally. Comparing these two strains, we find that haplotype differences at hotspots lead to qualitative and quantitative changes in PRDM9 binding and activity. Using Mus spretus as an outlier, we found most variants affecting PRDM9^Cst binding arose and were fixed in M.m. castaneus, suppressing hotspot activity. Furthermore, M.m. castaneus×M.m. domesticus F1 hybrids exhibit novel hotspots, with large haplotype biases in both PRDM9 binding and chromatin modification. These novel hotspots represent sites of historic evolutionary erosion that become activated in hybrids due to crosstalk between one parent's Prdm9 allele and the opposite parent's chromosome. Together these data support a model where haplotype-specific PRDM9 binding directs biased gene conversion at hotspots, ultimately leading to hotspot erosion.

Sexually reproducing creatures need to produce germ cells, notably sperm and egg, and do so using a specialized cell division, termed meiosis. A hallmark of meiosis is the process of recombination, in which pieces of maternal and paternal genetic material are exchanged, creating new combinations that are inherited by their progeny. Recombination increases diversity in subsequent generations, facilitating evolution. However, if recombination goes awry, it can lead to genetic disorders and spontaneous abortions in humans. In most mammals the sites of recombination are directed by the enzyme PRDM9 to specific regions on DNA, termed hotspots. For several decades it has been speculated that the process of recombination should lead to the eventual inactivation of hotspots, resulting in the loss of ability to reproduce. The discovery of PRDM9 provided a potential solution to this dilemma when the appearance of new PRDM9 alleles with altered DNA binding specificity would immediately create a new set of hotspots. We have now used the power of mouse genetics and large scale measurements of PRDM9 location and activity to show that this cycle of hotspot loss and recovery does indeed occur over the course of hundreds of thousands of years, and is directed by PRDM9.

Speciation with gene flow in equids despite extensive chromosomal plasticity

Horses, asses, and zebras belong to a single genus, Equus, which emerged 4.0-4.5 Mya. Although the equine fossil record represents a textbook example of evolution, the succession of events that gave rise to the diversity of species existing today remains unclear. Here we present six genomes from each living species of asses and zebras. This completes the set of genomes available for all extant species in the genus, which was hitherto represented only by the horse and the domestic donkey. In addition, we used a museum specimen to characterize the genome of the quagga zebra, which was driven to extinction in the early 1900s. We scan the genomes for lineage-specific adaptations and identify 48 genes that have evolved under positive selection and are involved in olfaction, immune response, development, locomotion, and behavior. Our extensive genome dataset reveals a highly dynamic demographic history with synchronous expansions and collapses on different continents during the last 400 ky after major climatic events. We show that the earliest speciation occurred with gene flow in Northern America, and that the ancestor of present-day asses and zebras dispersed into the Old World 2.1-3.4 Mya. Strikingly, we also find evidence for gene flow involving three contemporary equine species despite chromosomal numbers varying from 16 pairs to 31 pairs. These findings challenge the claim that the accumulation of chromosomal rearrangements drive complete reproductive isolation, and promote equids as a fundamental model for understanding the interplay between chromosomal structure, gene flow, and, ultimately, speciation.

Diversity of Prdm9 zinc finger array in wild mice unravels new facets of the evolutionary turnover of this coding minisatellite

In humans and mice, meiotic recombination events cluster into narrow hotspots whose genomic positions are defined by the PRDM9 protein via its DNA binding domain constituted of an array of zinc fingers (ZnFs). High polymorphism and rapid divergence of the Prdm9 gene ZnF domain appear to involve positive selection at DNA-recognition amino-acid positions, but the nature of the underlying evolutionary pressures remains a puzzle. Here we explore the variability of the Prdm9 ZnF array in wild mice, and uncovered a high allelic diversity of both ZnF copy number and identity with the caracterization of 113 alleles. We analyze features of the diversity of ZnF identity which is mostly due to non-synonymous changes at codons -1, 3 and 6 of each ZnF, corresponding to amino-acids involved in DNA binding. Using methods adapted to the minisatellite structure of the ZnF array, we infer a phylogenetic tree of these alleles. We find the sister species Mus spicilegus and M. macedonicus as well as the three house mouse (Mus musculus) subspecies to be polyphyletic. However some sublineages have expanded independently in Mus musculus musculus and M. m. domesticus, the latter further showing phylogeographic substructure. Compared to random genomic regions and non-coding minisatellites, none of these patterns appears exceptional. In silico prediction of DNA binding sites for each allele, overlap of their alignments to the genome and relative coverage of the different families of interspersed repeated elements suggest a large diversity between PRDM9 variants with a potential for highly divergent distributions of recombination events in the genome with little correlation to evolutionary distance. By compiling PRDM9 ZnF protein sequences in Primates, Muridae and Equids, we find different diversity patterns among the three amino-acids most critical for the DNA-recognition function, suggesting different diversification timescales.

Meiotic recombination in mammals: localization and regulation

During meiosis, a programmed induction of DNA double-strand breaks (DSBs) leads to the exchange of genetic material between homologous chromosomes. These exchanges increase genome diversity and are essential for proper chromosome segregation at the first meiotic division. Recent findings have highlighted an unexpected molecular control of the distribution of meiotic DSBs in mammals by a rapidly evolving gene, PR domain-containing 9 (PRDM9), and genome-wide analyses have facilitated the characterization of meiotic DSB sites at unprecedented resolution. In addition, the identification of new players in DSB repair processes has allowed the delineation of recombination pathways that have two major outcomes, crossovers and non-crossovers, which have distinct mechanistic roles and consequences for genome evolution.