Genetic recombination is targeted towards gene promoter regions in dogs

Ancient divergence of Indian and Tibetan wolves revealed by recombination-aware phylogenomics

The grey wolf (Canis lupus) expanded its range across Holarctic regions during the late Pleistocene. Consequently, most grey wolves share recent (<100,000 years ago) maternal origins corresponding to a widespread Holarctic clade. However, two deeply divergent (200,000-700,000 years ago) mitochondrial clades are restricted, respectively, to the Indian subcontinent and the Tibetan Plateau, where remaining wolves are endangered. No genome-wide analysis had previously included wolves corresponding to the mitochondrial Indian clade or attempted to parse gene flow and phylogeny. We sequenced four Indian and two Tibetan wolves and included 31 additional canid genomes to resolve the phylogenomic history of grey wolves. Genomic analyses revealed Indian and Tibetan wolves to be distinct from each other and from broadly distributed wolf populations corresponding to the mitochondrial Holarctic clade. Despite gene flow, which was reflected disproportionately in high-recombination regions of the genome, analyses revealed Indian and Tibetan wolves to be basal to Holarctic grey wolves, in agreement with the mitochondrial phylogeny. In contrast to mitochondrial DNA, however, genomic findings suggest the possibility that the Indian wolf could be basal to the Tibetan wolf, a discordance potentially reflecting selection on the mitochondrial genome. Together, these findings imply that southern regions of Asia have been important centers for grey wolf evolution and that Indian and Tibetan wolves represent evolutionary significant units (ESUs). Further study is needed to assess whether these ESUs warrant recognition as distinct species. This question is especially urgent regarding the Indian wolf, which represents one of the world's most endangered wolf populations.

Evolution of Recombination Landscapes in Diverging Populations of Bread Wheat

Reciprocal exchanges of DNA (crossovers) that occur during meiosis are mandatory to ensure the production of fertile gametes in sexually reproducing species. They also contribute to shuffle parental alleles into new combinations thereby fueling genetic variation and evolution. However, due to biological constraints, the recombination landscape is highly heterogeneous along the genome which limits the range of allelic combinations and the adaptability of populations. An approach to better understand the constraints on the recombination process is to study how it evolved in the past. In this work, we tackled this question by constructing recombination profiles in four diverging bread wheat (Triticum aestivum L.) populations established from 371 landraces genotyped at 200,062 SNPs. We used linkage disequilibrium (LD) patterns to estimate in each population the past distribution of recombination along the genome and characterize its fine-scale heterogeneity. At the megabase scale, recombination rates derived from LD patterns were consistent with family-based estimates obtained from a population of 406 recombinant inbred lines. Among the four populations, recombination landscapes were positively correlated between each other and shared a statistically significant proportion of highly recombinant intervals. However, this comparison also highlighted that the similarity in recombination landscapes between populations was significantly decreasing with their genetic differentiation in most regions of the genome. This observation was found to be robust to SNPs ascertainment and demography and suggests a relatively rapid evolution of factors determining the fine-scale localization of recombination in bread wheat.

Recombination Rate Variation and Infrequent Sex Influence Genetic Diversity in Chlamydomonas reinhardtii

Recombination confers a major evolutionary advantage by breaking up linkage disequilibrium between harmful and beneficial mutations, thereby facilitating selection. However, in species that are only periodically sexual, such as many microbial eukaryotes, the realized rate of recombination is also affected by the frequency of sex, meaning that infrequent sex can increase the effects of selection at linked sites despite high recombination rates. Despite this, the rate of sex of most facultatively sexual species is unknown. Here, we use genomewide patterns of linkage disequilibrium to infer fine-scale recombination rate variation in the genome of the facultatively sexual green alga Chlamydomonas reinhardtii. We observe recombination rate variation of up to two orders of magnitude and find evidence of recombination hotspots across the genome. Recombination rate is highest flanking genes, consistent with trends observed in other nonmammalian organisms, though intergenic recombination rates vary by intergenic tract length. We also find a positive relationship between nucleotide diversity and physical recombination rate, suggesting a widespread influence of selection at linked sites in the genome. Finally, we use estimates of the effective rate of recombination to calculate the rate of sex that occurs in natural populations, estimating a sexual cycle roughly every 840 generations. We argue that the relatively infrequent rate of sex and large effective population size creates a population genetic environment that increases the influence of selection on linked sites across the genome.

Rat PRDM9 shapes recombination landscapes, duration of meiosis, gametogenesis, and age of fertility

BACKGROUND

Vertebrate meiotic recombination events are concentrated in regions (hotspots) that display open chromatin marks, such as trimethylation of lysines 4 and 36 of histone 3 (H3K4me3 and H3K36me3). Mouse and human PRDM9 proteins catalyze H3K4me3 and H3K36me3 and determine hotspot positions, whereas other vertebrates lacking PRDM9 recombine in regions with chromatin already opened for another function, such as gene promoters. While these other vertebrate species lacking PRDM9 remain fertile, inactivation of the mouse Prdm9 gene, which shifts the hotspots to the functional regions (including promoters), typically causes gross fertility reduction; and the reasons for these species differences are not clear.

RESULTS

We introduced Prdm9 deletions into the Rattus norvegicus genome and generated the first rat genome-wide maps of recombination-initiating double-strand break hotspots. Rat strains carrying the same wild-type Prdm9 allele shared 88% hotspots but strains with different Prdm9 alleles only 3%. After Prdm9 deletion, rat hotspots relocated to functional regions, about 40% to positions corresponding to Prdm9-independent mouse hotspots, including promoters. Despite the hotspot relocation and decreased fertility, Prdm9-deficient rats of the SHR/OlaIpcv strain produced healthy offspring. The percentage of normal pachytene spermatocytes in SHR-Prdm9 mutants was almost double than in the PWD male mouse oligospermic sterile mutants. We previously found a correlation between the crossover rate and sperm presence in mouse Prdm9 mutants. The crossover rate of SHR is more similar to sperm-carrying mutant mice, but it did not fully explain the fertility of the SHR mutants. Besides mild meiotic arrests at rat tubular stages IV (mid-pachytene) and XIV (metaphase), we also detected postmeiotic apoptosis of round spermatids. We found delayed meiosis and age-dependent fertility in both sexes of the SHR mutants.

CONCLUSIONS

We hypothesize that the relative increased fertility of rat versus mouse Prdm9 mutants could be ascribed to extended duration of meiotic prophase I. While rat PRDM9 shapes meiotic recombination landscapes, it is unnecessary for recombination. We suggest that PRDM9 has additional roles in spermatogenesis and speciation-spermatid development and reproductive age-that may help to explain male-specific hybrid sterility.

Ancient Hybridization with an Unknown Population Facilitated High-Altitude Adaptation of Canids

Genetic introgression not only provides material for adaptive evolution but also confounds our understanding of evolutionary history. This is particularly true for canids, a species complex in which genome sequencing and analysis has revealed a complex history of admixture and introgression. Here, we sequence 19 new whole genomes from high-altitude Tibetan and Himalayan wolves and dogs and combine these into a larger data set of 166 whole canid genomes. Using these data, we explore the evolutionary history and adaptation of these and other canid lineages. We find that Tibetan and Himalayan wolves are closely related to each other, and that ∼39% of their nuclear genome is derived from an as-yet-unrecognized wolf-like lineage that is deeply diverged from living Holarctic wolves and dogs. The EPAS1 haplotype, which is present at high frequencies in Tibetan dog breeds and wolves and confers an adaptive advantage to animals living at high altitudes, was probably derived from this ancient lineage. Our study underscores the complexity of canid evolution and demonstrates how admixture and introgression can shape the evolutionary trajectories of species.

R-locus for roaned coat is associated with a tandem duplication in an intronic region of USH2A in dogs and also contributes to Dalmatian spotting

Structural variations (SVs) represent a large fraction of all genetic diversity, but how this genetic diversity is translated into phenotypic and organismal diversity is unclear. Explosive diversification of dog coat color and patterns after domestication can provide a unique opportunity to explore this question; however, the major obstacle is to efficiently collect a sufficient number of individuals with known phenotypes and genotypes of hundreds of thousands of markers. Using customer-provided information about coat color and patterns of dogs tested on a commercial canine genotyping platform, we identified a genomic region on chromosome 38 that is strongly associated with a mottled coat pattern (roaning) by genome-wide association study. We identified a putative causal variant in this region, an 11-kb tandem duplication (11,131,835–11,143,237) characterized by sequence read coverage and discordant reads of whole-genome sequence data, microarray probe intensity data, and a duplication-specific PCR assay. The tandem duplication is in an intronic region of usherin gene (USH2A), which was perfectly associated with roaning but absent in non-roaned dogs. We detected strong selection signals in this region characterized by reduced nucleotide diversity (π), increased runs of homozygosity, and extended haplotype homozygosity in Wirehaired Pointing Griffons and Australian Cattle Dogs (typically roaned breeds), as well as elevated genetic difference (F_ST) between Wirehaired Pointing Griffon (roaned) and Labrador Retriever (non-roaned). Surprisingly, all Dalmatians (N = 262) carried the duplication embedded in identical or similar haplotypes with roaned dogs, indicating this region as a shared target of selection during the breed’s formation. We propose that the Dalmatian’s unique spots were a derived coat pattern by establishing a novel epistatic interaction between roaning “R-locus” on chromosome 38 and an uncharacterized modifier locus. These results highlight the utility of consumer-oriented genotype and phenotype data in the discovery of genomic regions contributing to phenotypic diversity in dogs.

Snake Recombination Landscapes Are Concentrated in Functional Regions despite PRDM9

Meiotic recombination in vertebrates is concentrated in hotspots throughout the genome. The location and stability of hotspots have been linked to the presence or absence of PRDM9, leading to two primary models for hotspot evolution derived from mammals and birds. Species with PRDM9-directed recombination have rapid turnover of hotspots concentrated in intergenic regions (i.e., mammals), whereas hotspots in species lacking PRDM9 are concentrated in functional regions and have greater stability over time (i.e., birds). Snakes possess PRDM9, yet virtually nothing is known about snake recombination. Here, we examine the recombination landscape and test hypotheses about the roles of PRDM9 in rattlesnakes. We find substantial variation in recombination rate within and among snake chromosomes, and positive correlations between recombination rate and gene density, GC content, and genetic diversity. Like mammals, snakes appear to have a functional and active PRDM9, but rather than being directed away from genes, snake hotspots are concentrated in promoters and functional regions-a pattern previously associated only with species that lack a functional PRDM9. Snakes therefore provide a unique example of recombination landscapes in which PRDM9 is functional, yet recombination hotspots are associated with functional genic regions-a combination of features that defy existing paradigms for recombination landscapes in vertebrates. Our findings also provide evidence that high recombination rates are a shared feature of vertebrate microchromosomes. Our results challenge previous assumptions about the adaptive role of PRDM9 and highlight the diversity of recombination landscape features among vertebrate lineages.

Long-read assembly of a Great Dane genome highlights the contribution of GC-rich sequence and mobile elements to canine genomes

Technological advances have allowed improvements in genome reference sequence assemblies. Here, we combined long- and short-read sequence resources to assemble the genome of a female Great Dane dog. This assembly has improved continuity compared to the existing Boxer-derived (CanFam3.1) reference genome. Annotation of the Great Dane assembly identified 22,182 protein-coding gene models and 7,049 long noncoding RNAs, including 49 protein-coding genes not present in the CanFam3.1 reference. The Great Dane assembly spans the majority of sequence gaps in the CanFam3.1 reference and illustrates that 2,151 gaps overlap the transcription start site of a predicted protein-coding gene. Moreover, a subset of the resolved gaps, which have an 80.95% median GC content, localize to transcription start sites and recombination hotspots more often than expected by chance, suggesting the stable canine recombinational landscape has shaped genome architecture. Alignment of the Great Dane and CanFam3.1 assemblies identified 16,834 deletions and 15,621 insertions, as well as 2,665 deletions and 3,493 insertions located on secondary contigs. These structural variants are dominated by retrotransposon insertion/deletion polymorphisms and include 16,221 dimorphic canine short interspersed elements (SINECs) and 1,121 dimorphic long interspersed element-1 sequences (LINE-1_Cfs). Analysis of sequences flanking the 3' end of LINE-1_Cfs (i.e., LINE-1_Cf 3'-transductions) suggests multiple retrotransposition-competent LINE-1_Cfs segregate among dog populations. Consistent with this conclusion, we demonstrate that a canine LINE-1_Cf element with intact open reading frames can retrotranspose its own RNA and that of a SINEC_Cf consensus sequence in cultured human cells, implicating ongoing retrotransposon activity as a driver of canine genetic variation.

Genomes of Pleistocene Siberian Wolves Uncover Multiple Extinct Wolf Lineages

Extant Canis lupus genetic diversity can be grouped into three phylogenetically distinct clades: Eurasian and American wolves and domestic dogs.1 Genetic studies have suggested these groups trace their origins to a wolf population that expanded during the last glacial maximum (LGM)1-3 and replaced local wolf populations.4 Moreover, ancient genomes from the Yana basin and the Taimyr peninsula provided evidence of at least one extinct wolf lineage that dwelled in Siberia during the Pleistocene.35 Previous studies have suggested that Pleistocene Siberian canids can be classified into two groups based on cranial morphology. Wolves in the first group are most similar to present-day populations, although those in the second group possess intermediate features between dogs and wolves.67 However, whether this morphological classification represents distinct genetic groups remains unknown. To investigate this question and the relationships between Pleistocene canids, present-day wolves, and dogs, we resequenced the genomes of four Pleistocene canids from Northeast Siberia dated between >50 and 14 ka old, including samples from the two morphological categories. We found these specimens cluster with the two previously sequenced Pleistocene wolves, which are genetically more similar to Eurasian wolves. Our results show that, though the four specimens represent extinct wolf lineages, they do not form a monophyletic group. Instead, each Pleistocene Siberian canid branched off the lineage that gave rise to present-day wolves and dogs. Finally, our results suggest the two previously described morphological groups could represent independent lineages similarly related to present-day wolves and dogs.

Pleistocene climate fluctuations drove demographic history of African golden wolves (Canis lupaster)

Pleistocene climate change impacted entire ecosystems throughout the world. In the northern hemisphere, the distribution of Arctic species expanded during glacial periods, while more temperate and mesic species contracted into climatic refugia, where isolation drove genetic divergence. Cycles of local cooling and warming in the Sahara region of northern Africa caused repeated contractions and expansions of savannah-like environments which connected mesic species isolated in refugia during interglacial times, possibly driving population expansions and contractions; divergence and geneflow in the associated fauna. Here, we use whole genome sequences of African golden wolves (Canis lupaster), a generalist mesopredator with a wide distribution in northern Africa to estimate their demographic history and past episodes of geneflow. We detect a correlation between divergence times and cycles of increased aridity-associated Pleistocene glacial cycles. A complex demographic history with responses to local climate change in different lineages was found, including a relict lineage north of the High Atlas Mountains of Morocco that has been isolated for more than 18,000 years, possibly a distinct ecotype.

Examination of the efficacy of small genetic panels in genomic conservation of companion animal populations

In many ways, dogs are an ideal model for the study of genetic erosion and population recovery, problems of major concern in the field of conservation genetics. Genetic diversity in many dog breeds has been declining systematically since the beginning of the 1800s, when modern breeding practices came into fashion. As such, inbreeding in domestic dog breeds is substantial and widespread and has led to an increase in recessive deleterious mutations of high effect as well as general inbreeding depression. Pedigrees can in theory be used to guide breeding decisions, though are often incomplete and do not reflect the full history of inbreeding. Small microsatellite panels are also used in some cases to choose mating pairs to produce litters with low levels of inbreeding. However, the long-term impact of such practices has not been thoroughly evaluated. Here, we use forward simulation on a model of the dog genome to examine the impact of using limited marker panels to guide pairwise mating decisions on genome-wide population-level genetic diversity. Our results suggest that in unmanaged populations, where breeding decisions are made at the pairwise-rather than population-level, such panels can lead to accelerated loss of genetic diversity at genome regions unlinked to panel markers, compared to random mating. These results demonstrate the importance of genome-wide genetic panels for managing and conserving genetic diversity in dogs and other companion animals.

Sexual dimorphism in the meiotic requirement for PRDM9: A mammalian evolutionary safeguard

In many mammals, genomic sites for recombination are determined by the histone methyltransferase PRMD9. Some mouse strains lacking PRDM9 are infertile, but instances of fertility or semifertility in the absence of PRDM9 have been reported in mice, canines, and a human female. Such findings raise the question of how the loss of PRDM9 is circumvented to maintain fertility. We show that genetic background and sex-specific modifiers can obviate the requirement for PRDM9 in mice. Specifically, the meiotic DNA damage checkpoint protein CHK2 acts as a modifier allowing female-specific fertility in the absence of PRDM9. We also report that, in the absence of PRDM9, a PRDM9-independent recombination system is compatible with female meiosis and fertility, suggesting sex-specific regulation of meiotic recombination, a finding with implications for speciation.

Double-edged sword: The evolutionary consequences of the epigenetic silencing of transposable elements

Transposable elements (TEs) are genomic parasites that selfishly replicate at the expense of host fitness. Fifty years of evolutionary studies of TEs have concentrated on the deleterious genetic effects of TEs, such as their effects on disrupting genes and regulatory sequences. However, a flurry of recent work suggests that there is another important source of TEs' harmful effects-epigenetic silencing. Host genomes typically silence TEs by the deposition of repressive epigenetic marks. While this silencing reduces the selfish replication of TEs and should benefit hosts, a picture is emerging that the epigenetic silencing of TEs triggers inadvertent spreading of repressive marks to otherwise expressed neighboring genes, ultimately jeopardizing host fitness. In this Review, we provide a long-overdue overview of the recent genome-wide evidence for the presence and prevalence of TEs' epigenetic effects, highlighting both the similarities and differences across mammals, insects, and plants. We lay out the current understanding of the functional and fitness consequences of TEs' epigenetic effects, and propose possible influences of such effects on the evolution of both hosts and TEs themselves. These unique evolutionary consequences indicate that TEs' epigenetic effect is not only a crucial component of TE biology but could also be a significant contributor to genome function and evolution.

Arctic-adapted dogs emerged at the Pleistocene-Holocene transition

Although sled dogs are one of the most specialized groups of dogs, their origin and evolution has received much less attention than many other dog groups. We applied a genomic approach to investigate their spatiotemporal emergence by sequencing the genomes of 10 modern Greenland sled dogs, an ~9500-year-old Siberian dog associated with archaeological evidence for sled technology, and an ~33,000-year-old Siberian wolf. We found noteworthy genetic similarity between the ancient dog and modern sled dogs. We detected gene flow from Pleistocene Siberian wolves, but not modern American wolves, to present-day sled dogs. The results indicate that the major ancestry of modern sled dogs traces back to Siberia, where sled dog-specific haplotypes of genes that potentially relate to Arctic adaptation were established by 9500 years ago.

Genetic differentiation and intrinsic genomic features explain variation in recombination hotspots among cocoa tree populations

Background

Recombination plays an important evolutionary role by breaking up haplotypes and shuffling genetic variation. This process impacts the ability of selection to eliminate deleterious mutations or increase the frequency of beneficial mutations in a population. To understand the role of recombination generating and maintaining haplotypic variation in a population, we can construct fine-scale recombination maps. Such maps have been used to study a variety of model organisms and proven to be informative of how selection and demographics shape species-wide variation. Here we present a fine-scale recombination map for ten populations of Theobroma cacao – a non-model, long-lived, woody crop. We use this map to elucidate the dynamics of recombination rates in distinct populations of the same species, one of which is domesticated.

Results

Mean recombination rates in range between 2.5 and 8.6 cM/Mb for most populations of T. cacao with the exception of the domesticated Criollo (525 cM/Mb) and Guianna, a more recently established population (46.5 cM/Mb). We found little overlap in the location of hotspots of recombination across populations. We also found that hotspot regions contained fewer known retroelement sequences than expected and were overrepresented near transcription start and termination sites. We find mutations in FIGL-1, a protein shown to downregulate cross-over frequency in Arabidopsis, statistically associated to higher recombination rates in domesticated Criollo.

Conclusions

We generated fine-scale recombination maps for ten populations of Theobroma cacao and used them to understand what processes are associated with population-level variation in this species. Our results provide support to the hypothesis of increased recombination rates in domesticated plants (Criollo population). We propose a testable mechanistic hypothesis for the change in recombination rate in domesticated populations in the form of mutations to a previously identified recombination-suppressing protein. Finally, we establish a number of possible correlates of recombination hotspots that help explain general patterns of recombination in this species.

Genomic regions under selection in the feralization of the dingoes

Dingoes are wild canids living in Australia, originating from domestic dogs. They have lived isolated from both the wild and the domestic ancestor, making them a unique model for studying feralization. Here, we sequence the genomes of 10 dingoes and 2 New Guinea Singing Dogs. Phylogenetic and demographic analyses show that dingoes originate from dogs in southern East Asia, which migrated via Island Southeast Asia to reach Australia around 8300 years ago, and subsequently diverged into a genetically distinct population. Selection analysis identifies 50 positively selected genes enriched in digestion and metabolism, indicating a diet change during feralization of dingoes. Thirteen of these genes have shifted allele frequencies compared to dogs but not compared to wolves. Functional assays show that an A-to-G mutation in ARHGEF7 decreases the endogenous expression, suggesting behavioral adaptations related to the transitions in environment. Our results indicate that the feralization of the dingo induced positive selection on genomic regions correlated to neurodevelopment, metabolism and reproduction, in adaptation to a wild environment.

Parallel Molecular Evolution in Pathways, Genes, and Sites in High-Elevation Hummingbirds Revealed by Comparative Transcriptomics

High-elevation organisms experience shared environmental challenges that include low oxygen availability, cold temperatures, and intense ultraviolet radiation. Consequently, repeated evolution of the same genetic mechanisms may occur across high-elevation taxa. To test this prediction, we investigated the extent to which the same biochemical pathways, genes, or sites were subject to parallel molecular evolution for 12 Andean hummingbird species (family: Trochilidae) representing several independent transitions to high elevation across the phylogeny. Across high-elevation species, we discovered parallel evolution for several pathways and genes with evidence of positive selection. In particular, positively selected genes were frequently part of cellular respiration, metabolism, or cell death pathways. To further examine the role of elevation in our analyses, we compared results for low- and high-elevation species and tested different thresholds for defining elevation categories. In analyses with different elevation thresholds, positively selected genes reflected similar functions and pathways, even though there were almost no specific genes in common. For example, EPAS1 (HIF2α), which has been implicated in high-elevation adaptation in other vertebrates, shows a signature of positive selection when high-elevation is defined broadly (>1,500 m), but not when defined narrowly (>2,500 m). Although a few biochemical pathways and genes change predictably as part of hummingbird adaptation to high-elevation conditions, independent lineages have rarely adapted via the same substitutions.

Genomic Approaches Reveal an Endemic Subpopulation of Gray Wolves in Southern China

Although gray wolves (Canis lupus) are one of the most widely distributed terrestrial mammals, their origins in China are not well understood. We sequenced six specimens from wolf skins, showing that gray wolves from Southern China (SC) derive from a single lineage, distinct from gray wolves from the Tibetan Plateau and Northern China, suggesting that SC gray wolves may form a distinct subpopulation. Of SC gray wolves, one wolf from Zhejiang carries a genetic component from a canid and had gene flow from a population related to or further diverged from wolves than the dhole. This may indicate that interspecific gene flow likely played an important role in shaping the speciation patterns and population structure in the genus Canis. Our study is the first to survey museum gray wolves' genomes from Southern China, highlighting how sequencing the paleogenome from museum specimens can help us to study extinct species.

Inference and analysis of population-specific fine-scale recombination maps across 26 diverse human populations

Fine-scale rates of meiotic recombination vary by orders of magnitude across the genome and differ between species and even populations. Studying cross-population differences has been stymied by the confounding effects of demographic history. To address this problem, we developed a demography-aware method to infer fine-scale recombination rates and applied it to 26 diverse human populations, inferring population-specific recombination maps. These maps recapitulate many aspects of the history of these populations including signatures of the trans-Atlantic slave trade and the Iberian colonization of the Americas. We also investigated modulators of the local recombination rate, finding further evidence that Polycomb group proteins and the trimethylation of H3K27 elevate recombination rates. Further differences in the recombination landscape across the genome and between populations are driven by variation in the gene that encodes the DNA binding protein PRDM9, and we quantify the weak effect of meiotic drive acting to remove its binding sites.

PRDM9 forms a trimer by interactions within the zinc finger array

PRDM9 is a trans-acting factor directing meiotic recombination to specific DNA-binding sites by its zinc finger (ZnF) array. It was suggested that PRDM9 is a multimer; however, we do not know the stoichiometry or the components inducing PRDM9 multimerization. In this work, we used in vitro binding studies and characterized with electrophoretic mobility shift assays, mass spectrometry, and fluorescence correlation spectroscopy the stoichiometry of the PRDM9 multimer of two different murine PRDM9 alleles carrying different tags and domains produced with different expression systems. Based on the migration distance of the PRDM9-DNA complex, we show that PRDM9 forms a trimer. Moreover, this stoichiometry is adapted already by the free, soluble protein with little exchange between protein monomers. The variable ZnF array of PRDM9 is sufficient for multimerization, and at least five ZnFs form already a functional trimer. Finally, we also show that only one ZnF array within the PRDM9 oligomer binds to the DNA, whereas the remaining two ZnF arrays likely maintain the trimer by ZnF-ZnF interactions.

Comparative genomics provides new insights into the remarkable adaptations of the African wild dog (Lycaon pictus)

Within the Canidae, the African wild dog (Lycaon pictus) is the most specialized with regards to cursorial adaptations (specialized for running), having only four digits on their forefeet. In addition, this species is one of the few canids considered to be an obligate meat-eater, possessing a robust dentition for taking down large prey, and displays one of the most variable coat colorations amongst mammals. Here, we used comparative genomic analysis to investigate the evolutionary history and genetic basis for adaptations associated with cursoriality, hypercanivory, and coat color variation in African wild dogs. Genome-wide scans revealed unique amino acid deletions that suggest a mode of evolutionary digit loss through expanded apoptosis in the developing first digit. African wild dog-specific signals of positive selection also uncovered a putative mechanism of molar cusp modification through changes in genes associated with the sonic hedgehog (SHH) signaling pathway, required for spatial patterning of teeth, and three genes associated with pigmentation. Divergence time analyses suggest the suite of genomic changes we identified evolved ~1.7 Mya, coinciding with the diversification of large-bodied ungulates. Our results show that comparative genomics is a powerful tool for identifying the genetic basis of evolutionary changes in Canidae.

Divergent Fine-Scale Recombination Landscapes between a Freshwater and Marine Population of Threespine Stickleback Fish

Meiotic recombination is a highly conserved process that has profound effects on genome evolution. At a fine-scale, recombination rates can vary drastically across genomes, often localized into small recombination "hotspots" with highly elevated rates, surrounded by regions with little recombination. In most species studied, the location of hotspots within genomes is highly conserved across broad evolutionary timescales. The main exception to this pattern is in mammals, where hotspot location can evolve rapidly among closely related species and even among populations within a species. Hotspot position in mammals is controlled by the gene, Prdm9, whereas in species with conserved hotspots, a functional Prdm9 is typically absent. Due to a limited number of species where recombination rates have been estimated at a fine-scale, it remains unclear whether hotspot conservation is always associated with the absence of a functional Prdm9. Threespine stickleback fish (Gasterosteus aculeatus) are an excellent model to examine the evolution of recombination over short evolutionary timescales. Using a linkage disequilibrium-based approach, we found recombination rates indeed varied at a fine-scale across the genome, with many regions organized into narrow hotspots. Hotspots had highly divergent landscapes between stickleback populations, where only ∼15% of these hotspots were shared. Our results indicate that fine-scale recombination rates may be diverging between closely related populations of threespine stickleback fish. Interestingly, we found only a weak association of a PRDM9 binding motif within hotspots, which suggests that threespine stickleback fish may possess a novel mechanism for targeting recombination hotspots at a fine-scale.

A mutation in the endonuclease domain of mouse MLH3 reveals novel roles for MutLγ during crossover formation in meiotic prophase I

During meiotic prophase I, double-strand breaks (DSBs) initiate homologous recombination leading to non-crossovers (NCOs) and crossovers (COs). In mouse, 10% of DSBs are designated to become COs, primarily through a pathway dependent on the MLH1-MLH3 heterodimer (MutLγ). Mlh3 contains an endonuclease domain that is critical for resolving COs in yeast. We generated a mouse (Mlh3DN/DN) harboring a mutation within this conserved domain that is predicted to generate a protein that is catalytically inert. Mlh3DN/DN males, like fully null Mlh3-/- males, have no spermatozoa and are infertile, yet spermatocytes have grossly normal DSBs and synapsis events in early prophase I. Unlike Mlh3-/- males, mutation of the endonuclease domain within MLH3 permits normal loading and frequency of MutLγ in pachynema. However, key DSB repair factors (RAD51) and mediators of CO pathway choice (BLM helicase) persist into pachynema in Mlh3DN/DN males, indicating a temporal delay in repair events and revealing a mechanism by which alternative DSB repair pathways may be selected. While Mlh3DN/DN spermatocytes retain only 22% of wildtype chiasmata counts, this frequency is greater than observed in Mlh3-/- males (10%), suggesting that the allele may permit partial endonuclease activity, or that other pathways can generate COs from these MutLγ-defined repair intermediates in Mlh3DN/DN males. Double mutant mice homozygous for the Mlh3DN/DN and Mus81-/- mutations show losses in chiasmata close to those observed in Mlh3-/- males, indicating that the MUS81-EME1-regulated crossover pathway can only partially account for the increased residual chiasmata in Mlh3DN/DN spermatocytes. Our data demonstrate that mouse spermatocytes bearing the MLH1-MLH3DN/DN complex display the proper loading of factors essential for CO resolution (MutSγ, CDK2, HEI10, MutLγ). Despite these functions, mice bearing the Mlh3DN/DN allele show defects in the repair of meiotic recombination intermediates and a loss of most chiasmata.

Extreme Differences in Recombination Rate between the Genomes of a Solitary and a Social Bee

Social insect genomes exhibit the highest rates of crossing over observed in plants and animals. The evolutionary causes of these extreme rates are unknown. Insight can be gained by comparing recombination rate variation across the genomes of related social and solitary insects. Here, we compare the genomic recombination landscape of the highly social honey bee, Apis mellifera, with the solitary alfalfa leafcutter bee, Megachile rotundata, by analyzing patterns of linkage disequilibrium in population-scale genome sequencing data. We infer that average recombination rates are extremely elevated in A. mellifera compared with M. rotundata. However, our results indicate that similar factors control the distribution of crossovers in the genomes of both species. Recombination rate is significantly reduced in coding regions in both species, with genes inferred to be germline methylated having particularly low rates. Genes with worker-biased patterns of expression in A. mellifera and their orthologs in M. rotundata have higher than average recombination rates in both species, suggesting that selection for higher diversity in genes involved in worker caste functions in social taxa is not the explanation for these elevated rates. Furthermore, we find no evidence that recombination has modulated the efficacy of selection among genes during bee evolution, which does not support the hypothesis that high recombination rates facilitated positive selection for new functions in social insects. Our results indicate that the evolution of sociality in insects likely entailed selection on modifiers that increased recombination rates genome wide, but that the genomic recombination landscape is determined by the same factors.

Prdm9 and Meiotic Cohesin Proteins Cooperatively Promote DNA Double-Strand Break Formation in Mammalian Spermatocytes

Meiotic recombination is required for correct segregation of chromosomes to gametes and to generate genetic diversity. In mice and humans, DNA double-strand breaks (DSBs) are initiated by SPO11 at recombination hotspots activated by PRDM9-catalyzed histone modifications on open chromatin. However, the DSB-initiating and repair proteins are associated with a linear proteinaceous scaffold called the chromosome axis, the core of which is composed of cohesin proteins. STAG3 is a stromalin subunit common to all meiosis-specific cohesin complexes. Mutations of meiotic cohesin proteins, especially STAG3, perturb both axis formation and recombination in the mouse, prompting determination of how the processes are mechanistically related. Protein interaction and genetic analyses revealed that PRDM9 interacts with STAG3 and REC8 in cooperative relationships that promote normal levels of meiotic DSBs at recombination hotspots in spermatocytes. The efficacy of the Prdm9-Stag3 genetic interaction in promoting DSB formation depends on PRDM9-mediated histone methyltransferase activity. Moreover, STAG3 deficiency has a major effect on DSB number even in the absence of PRDM9, showing that its role is not restricted to canonical PRDM9-activated hotspots. STAG3 and REC8 promote axis localization of the DSB-promoting proteins HORMAD1, IHO1, and MEI4, as well as SPO11 activity. These results establish that PRDM9 and axis-associated cohesin complexes together coordinate and facilitate meiotic recombination by recruiting key proteins for initiation of DSBs, thereby associating activated hotspots with DSB-initiating complexes on the axis.

Natural Selection and Origin of a Melanistic Allele in North American Gray Wolves

Pigmentation is often used to understand how natural selection affects genetic variation in wild populations since it can have a simple genetic basis, and can affect a variety of fitness-related traits (e.g., camouflage, thermoregulation, and sexual display). In gray wolves, the K locus, a β-defensin gene, causes black coat color via a dominantly inherited KB allele. The allele is derived from dog-wolf hybridization and is at high frequency in North American wolf populations. We designed a DNA capture array to probe the geographic origin, age, and number of introgression events of the KB allele in a panel of 331 wolves and 20 dogs. We found low diversity in KB, but not ancestral ky, wolf haplotypes consistent with a selective sweep of the black haplotype across North America. Further, North American wolf KB haplotypes are monophyletic, suggesting that a single adaptive introgression from dogs to wolves most likely occurred in the Northwest Territories or Yukon. We use a new analytical approach to date the origin of the KB allele in Yukon wolves to between 1,598 and 7,248 years ago, suggesting that introgression with early Native American dogs was the source. Using population genetic simulations, we show that the K locus is undergoing natural selection in four wolf populations. We find evidence for balancing selection, specifically in Yellowstone wolves, which could be a result of selection for enhanced immunity in response to distemper. With these data, we demonstrate how the spread of an adaptive variant may have occurred across a species' geographic range.

A common genomic code for chromatin architecture and recombination landscape

Recent findings established a link between DNA sequence composition and interphase chromatin architecture and explained the evolutionary conservation of TADs (Topologically Associated Domains) and LADs (Lamina Associated Domains) in mammals. This prompted us to analyse conformation capture and recombination rate data to study the relationship between chromatin architecture and recombination landscape of human and mouse genomes. The results reveal that: (1) low recombination domains and blocks of elevated linkage disequilibrium tend to coincide with TADs and isochores, indicating co-evolving regulatory elements and genes in insulated neighbourhoods; (2) double strand break (DSB) and recombination frequencies increase in the short loops of GC-rich TADs, whereas recombination cold spots are typical of LADs and (3) the binding and loading of proteins, which are critical for DSB and meiotic recombination (SPO11, DMC1, H3K4me3 and PRMD9) are higher in GC-rich TADs. One explanation for these observations is that the occurrence of DSB and recombination in meiotic cells are associated with compositional and epigenetic features (genomic code) that influence DNA stiffness/flexibility and appear to be similar to those guiding the chromatin architecture in the interphase nucleus of pre-leptotene cells.

Origin and recent expansion of an endogenous gammaretroviral lineage in domestic and wild canids

Background

Vertebrate genomes contain a record of retroviruses that invaded the germlines of ancestral hosts and are passed to offspring as endogenous retroviruses (ERVs). ERVs can impact host function since they contain the necessary sequences for expression within the host. Dogs are an important system for the study of disease and evolution, yet no substantiated reports of infectious retroviruses in dogs exist. Here, we utilized Illumina whole genome sequence data to assess the origin and evolution of a recently active gammaretroviral lineage in domestic and wild canids.

Results

We identified numerous recently integrated loci of a canid-specific ERV-Fc sublineage within Canis, including 58 insertions that were absent from the reference assembly. Insertions were found throughout the dog genome including within and near gene models. By comparison of orthologous occupied sites, we characterized element prevalence across 332 genomes including all nine extant canid species, revealing evolutionary patterns of ERV-Fc segregation among species as well as subpopulations.

Conclusions

Sequence analysis revealed common disruptive mutations, suggesting a predominant form of ERV-Fc spread by trans complementation of defective proviruses. ERV-Fc activity included multiple circulating variants that infected canid ancestors from the last 20 million to within 1.6 million years, with recent bursts of germline invasion in the sublineage leading to wolves and dogs.

Electronic supplementary material

The online version of this article (10.1186/s12977-019-0468-z) contains supplementary material, which is available to authorized users.

Historical Meiotic Crossover Hotspots Fueled Patterns of Evolutionary Divergence in Rice

Recombination plays an integral role in the creation of novel genetic variation in sexually reproducing species. Despite this important role, the determinants and evolution of crossover hotspots have remained poorly understood in plants. Here, we present a comparative analysis of two rice (Oryza sativa) historical recombination maps from two subspecies (indica and japonica) using 150 resequenced genomes. Fine-scale recombination rates and crossover hotspots were validated by comparison with a consensus genetic map and empirically derived crossovers, respectively. Strikingly, nearly 80% of crossover hotspots were unique to each subspecies, despite their relatively recent divergence and broad-scale correlated recombination rates. Crossover hotspots were enriched with Stowaway and P instability factor (PIF)/Harbinger transposons and overlapped accessible chromatin regions. Increased nucleotide diversity and signatures of population differentiation augmented by Stowaway and PIF/Harbinger transposons were prevalent at subspecies-specific crossover hotspots. Motifs derived from lineage-specific indica and japonica crossover hotspots were nearly identical in the two subspecies, implicating a core set of crossover motifs in rice. Finally, Stowaway and PIF/Harbinger transposons were associated with stabilized G/C bias within highly active hotspots, suggesting that hotspot activity can be fueled by de novo variation. These results provide evolutionary insight into historical crossover hotspots as potentially powerful drivers of sequence and subspecies evolution in plants.

A Molecular Portrait of De Novo Genes in Yeasts

New genes, with novel protein functions, can evolve "from scratch" out of intergenic sequences. These de novo genes can integrate the cell's genetic network and drive important phenotypic innovations. Therefore, identifying de novo genes and understanding how the transition from noncoding to coding occurs are key problems in evolutionary biology. However, identifying de novo genes is a difficult task, hampered by the presence of remote homologs, fast evolving sequences and erroneously annotated protein coding genes. To overcome these limitations, we developed a procedure that handles the usual pitfalls in de novo gene identification and predicted the emergence of 703 de novo gene candidates in 15 yeast species from 2 genera whose phylogeny spans at least 100 million years of evolution. We validated 85 candidates by proteomic data, providing new translation evidence for 25 of them through mass spectrometry experiments. We also unambiguously identified the mutations that enabled the transition from noncoding to coding for 30 Saccharomyces de novo genes. We established that de novo gene origination is a widespread phenomenon in yeasts, only a few being ultimately maintained by selection. We also found that de novo genes preferentially emerge next to divergent promoters in GC-rich intergenic regions where the probability of finding a fortuitous and transcribed ORF is the highest. Finally, we found a more than 3-fold enrichment of de novo genes at recombination hot spots, which are GC-rich and nucleosome-free regions, suggesting that meiotic recombination contributes to de novo gene emergence in yeasts.

iDog: an integrated resource for domestic dogs and wild canids

The domestic dog (Canis lupus familiaris) is indisputably one of man's best friends. It is also a fundamental model for many heritable human diseases. Here, we present iDog (http://bigd.big.ac.cn/idog), the first integrated resource dedicated to domestic dogs and wild canids. It incorporates a variety of omics data, including genome sequences assemblies for dhole and wolf, genomic variations extracted from hundreds of dog/wolf whole genomes, phenotype/disease traits curated from dog research communities and public resources, gene expression profiles derived from published RNA-Seq data, gene ontology for functional annotation, homolog gene information for multiple organisms and disease-related literature. Additionally, iDog integrates sequence alignment tools for data analyses and a genome browser for data visualization. iDog will not only benefit the global dog research community, but also provide access to a user-friendly consolidation of dog information to a large number of dog enthusiasts.

Hotspots for Initiation of Meiotic Recombination

Homologous chromosomes must pair and recombine to ensure faithful chromosome segregation during meiosis, a specialized type of cell division that occurs in sexually reproducing eukaryotes. Meiotic recombination initiates by programmed induction of DNA double-strand breaks (DSBs) by the conserved type II topoisomerase-like enzyme SPO11. A subset of meiotic DSBs are resolved as crossovers, whereby reciprocal exchange of DNA occurs between homologous chromosomes. Importantly, DSBs are non-randomly distributed along eukaryotic chromosomes, forming preferentially in permissive regions known as hotspots. In many species, including plants, DSB hotspots are located within nucleosome-depleted regions. DSB localization is governed by interconnected factors, including cis-regulatory elements, transcription factor binding, and chromatin accessibility, as well as by higher-order chromosome architecture. The spatiotemporal control of DSB formation occurs within a specialized chromosomal structure characterized by sister chromatids organized into linear arrays of chromatin loops that are anchored to a proteinaceous axis. Although SPO11 and its partner proteins required for DSB formation are bound to the axis, DSBs occur preferentially within the chromatin loops, which supports the "tethered-loop/axis model" for meiotic recombination. In this mini review, we discuss insights gained from recent efforts to define and profile DSB hotspots at high resolution in eukaryotic genomes. These advances are deepening our understanding of how meiotic recombination shapes genetic diversity and genome evolution in diverse species.

CXXC1 is not essential for normal DNA double-strand break formation and meiotic recombination in mouse

In most mammals, including mice and humans, meiotic recombination is determined by the meiosis specific histone methytransferase PRDM9, which binds to specific DNA sequences and trimethylates histone 3 at lysine-4 and lysine-36 at the adjacent nucleosomes. These actions ensure successful DNA double strand break formation and repair that occur on the proteinaceous structure forming the chromosome axis. The process of hotspot association with the axis after their activation by PRDM9 is poorly understood. Previously, we and others have identified CXXC1, an ortholog of S. cerevisiae Spp1 in mammals, as a PRDM9 interactor. In yeast, Spp1 is a histone methyl reader that links H3K4me3 sites with the recombination machinery, promoting DSB formation. Here, we investigated whether CXXC1 has a similar function in mouse meiosis. We created two Cxxc1 conditional knockout mouse models to deplete CXXC1 generally in germ cells, and before the onset of meiosis. Surprisingly, male knockout mice were fertile, and the loss of CXXC1 in spermatocytes had no effect on PRDM9 hotspot trimethylation, double strand break formation or repair. Our results demonstrate that CXXC1 is not an essential link between PRDM9-activated recombination hotspot sites and DSB machinery and that the hotspot recognition pathway in mouse is independent of CXXC1.

Meiotic recombination increases genetic diversity by ensuring novel combination of alleles passing correctly to the next generation. In most mammals, the meiotic recombination sites are determined by histone methyltransferase PRDM9. These sites are proposed to become associated with the chromosome axis with the participation of additional proteins and undergo double strand breaks, which are repaired by homologous recombination. In budding yeast, Spp1 (ortholog of CXXC1) binds to methylated H3K4 and connects these sites with the chromosome axis promoting DSB formation. However, our data suggest that even though CXXC1 interacts with PRDM9 in male germ cells, it does not play a crucial role in mouse meiotic recombination. These results indicate that, unlike in yeast, a recombination initiation pathway that includes CXXC1 could only serve as a non-essential pathway in mouse meiosis.

Interspecific Gene Flow Shaped the Evolution of the Genus Canis

The evolutionary history of the wolf-like canids of the genus Canis has been heavily debated, especially regarding the number of distinct species and their relationships at the population and species level [1-6]. We assembled a dataset of 48 resequenced genomes spanning all members of the genus Canis except the black-backed and side-striped jackals, encompassing the global diversity of seven extant canid lineages. This includes eight new genomes, including the first resequenced Ethiopian wolf (Canis simensis), one dhole (Cuon alpinus), two East African hunting dogs (Lycaon pictus), two Eurasian golden jackals (Canis aureus), and two Middle Eastern gray wolves (Canis lupus). The relationships between the Ethiopian wolf, African golden wolf, and golden jackal were resolved. We highlight the role of interspecific hybridization in the evolution of this charismatic group. Specifically, we find gene flow between the ancestors of the dhole and African hunting dog and admixture between the gray wolf, coyote (Canis latrans), golden jackal, and African golden wolf. Additionally, we report gene flow from gray and Ethiopian wolves to the African golden wolf, suggesting that the African golden wolf originated through hybridization between these species. Finally, we hypothesize that coyotes and gray wolves carry genetic material derived from a "ghost" basal canid lineage.

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.

Control of meiotic double-strand-break formation by ATM: local and global views

DNA double-strand breaks (DSBs) generated by the SPO11 protein initiate meiotic recombination, an essential process for successful chromosome segregation during gametogenesis. The activity of SPO11 is controlled by multiple factors and regulatory mechanisms, such that the number of DSBs is limited and DSBs form at distinct positions in the genome and at the right time. Loss of this control can affect genome integrity or cause meiotic arrest by mechanisms that are not fully understood. Here we focus on the DSB-responsive kinase ATM and its functions in regulating meiotic DSB numbers and distribution. We review the recently discovered roles of ATM in this context, discuss their evolutionary conservation, and examine future research perspectives.

Comparison of village dog and wolf genomes highlights the role of the neural crest in dog domestication

BACKGROUND

Domesticated from gray wolves between 10 and 40 kya in Eurasia, dogs display a vast array of phenotypes that differ from their ancestors, yet mirror other domesticated animal species, a phenomenon known as the domestication syndrome. Here, we use signatures persisting in dog genomes to identify genes and pathways possibly altered by the selective pressures of domestication.

RESULTS

Whole-genome SNP analyses of 43 globally distributed village dogs and 10 wolves differentiated signatures resulting from domestication rather than breed formation. We identified 246 candidate domestication regions containing 10.8 Mb of genome sequence and 429 genes. The regions share haplotypes with ancient dogs, suggesting that the detected signals are not the result of recent selection. Gene enrichments highlight numerous genes linked to neural crest and central nervous system development as well as neurological function. Read depth analysis suggests that copy number variation played a minor role in dog domestication.

CONCLUSIONS

Our results identify genes that act early in embryogenesis and can confer phenotypes distinguishing domesticated dogs from wolves, such as tameness, smaller jaws, floppy ears, and diminished craniofacial development as the targets of selection during domestication. These differences reflect the phenotypes of the domestication syndrome, which can be explained by alterations in the migration or activity of neural crest cells during development. We propose that initial selection during early dog domestication was for behavior, a trait influenced by genes which act in the neural crest, which secondarily gave rise to the phenotypes of modern dogs.

Coevolution between transposable elements and recombination

One of the most striking patterns of genome structure is the tight, typically negative, association between transposable elements (TEs) and meiotic recombination rates. While this is a highly recurring feature of eukaryotic genomes, the mechanisms driving correlations between TEs and recombination remain poorly understood, and distinguishing cause versus effect is challenging. Here, we review the evidence for a relation between TEs and recombination, and discuss the underlying evolutionary forces. Evidence to date suggests that overall TE densities correlate negatively with recombination, but the strength of this correlation varies across element types, and the pattern can be reversed. Results suggest that heterogeneity in the strength of selection against ectopic recombination and gene disruption can drive TE accumulation in regions of low recombination, but there is also strong evidence that the regulation of TEs can influence local recombination rates. We hypothesize that TE insertion polymorphism may be important in driving within-species variation in recombination rates in surrounding genomic regions. Furthermore, the interaction between TEs and recombination may create positive feedback, whereby TE accumulation in non-recombining regions contributes to the spread of recombination suppression. Further investigation of the coevolution between recombination and TEs has important implications for our understanding of the evolution of recombination rates and genome structure.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Sex Differences in Recombination in Sticklebacks

Recombination often differs markedly between males and females. Here we present the first analysis of sex-specific recombination in Gasterosteus sticklebacks. Using whole-genome sequencing of 15 crosses between G. aculeatus and G. nipponicus, we localized 698 crossovers with a median resolution of 2.3 kb. We also used a bioinformatic approach to infer historical sex-averaged recombination patterns for both species. Recombination is greater in females than males on all chromosomes, and overall map length is 1.64 times longer in females. The locations of crossovers differ strikingly between sexes. Crossovers cluster toward chromosome ends in males, but are distributed more evenly across chromosomes in females. Suppression of recombination near the centromeres in males causes crossovers to cluster at the ends of long arms in acrocentric chromosomes, and greatly reduces crossing over on short arms. The effect of centromeres on recombination is much weaker in females. Genomic differentiation between G. aculeatus and G. nipponicus is strongly correlated with recombination rate, and patterns of differentiation along chromosomes are strongly influenced by male-specific telomere and centromere effects. We found no evidence for fine-scale correlations between recombination and local gene content in either sex. We discuss hypotheses for the origin of sexual dimorphism in recombination and its consequences for sexually antagonistic selection and sex chromosome evolution.

Analysis of the canid Y-chromosome phylogeny using short-read sequencing data reveals the presence of distinct haplogroups among Neolithic European dogs

Background

Most genetic analyses of ancient and modern dogs have focused on variation in the autosomes or on the mitochondria. Mitochondrial DNA is more easily obtained from ancient samples than nuclear DNA and mitochondrial analyses have revealed important insights into the evolutionary history of canids. Utilizing a recently published dog Y-chromosome reference, we analyzed Y-chromosome sequence across a diverse collection of canids and determined the Y haplogroup of three ancient European dogs.

Results

We identified 1121 biallelic Y-chromosome SNVs using whole-genome sequences from 118 canids and defined variants diagnostic to distinct dog Y haplogroups. Similar to that of the mitochondria and previous more limited studies of Y diversity, we observe several deep splits in the Y-chromosome tree which may be the result of retained Y-chromosome diversity which predates dog domestication or post-domestication admixture with wolves. We find that Y-chromosomes from three ancient European dogs (4700–7000 years old) belong to distinct clades.

Conclusions

We estimate that the time to the most recent comment ancestor of dog Y haplogroups is 68–151 thousand years ago. Analysis of three Y-chromosomes from the Neolithic confirms long stranding population structure among European dogs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4749-z) contains supplementary material, which is available to authorized users.

Characterization of recombination features and the genetic basis in multiple cattle breeds

BACKGROUND

Crossover generated by meiotic recombination is a fundamental event that facilitates meiosis and sexual reproduction. Comparative studies have shown wide variation in recombination rate among species, but the characterization of recombination features between cattle breeds has not yet been performed. Cattle populations in North America count millions, and the dairy industry has genotyped millions of individuals with pedigree information that provide a unique opportunity to study breed-level variations in recombination.

RESULTS

Based on large pedigrees of Jersey, Ayrshire and Brown Swiss cattle with genotype data, we identified over 3.4 million maternal and paternal crossover events from 161,309 three-generation families. We constructed six breed- and sex-specific genome-wide recombination maps using 58,982 autosomal SNPs for two sexes in the three dairy cattle breeds. A comparative analysis of the six recombination maps revealed similar global recombination patterns between cattle breeds but with significant differences between sexes. We confirmed that male recombination map is 10% longer than the female map in all three cattle breeds, consistent with previously reported results in Holstein cattle. When comparing recombination hotspot regions between cattle breeds, we found that 30% and 10% of the hotspots were shared between breeds in males and females, respectively, with each breed exhibiting some breed-specific hotspots. Finally, our multiple-breed GWAS found that SNPs in eight loci affected recombination rate and that the PRDM9 gene associated with hotspot usage in multiple cattle breeds, indicating a shared genetic basis for recombination across dairy cattle breeds.

CONCLUSIONS

Collectively, our results generated breed- and sex-specific recombination maps for multiple cattle breeds, provided a comprehensive characterization and comparison of recombination patterns between breeds, and expanded our understanding of the breed-level variations in recombination features within an important livestock species.

Construction of PRDM9 allele-specific recombination maps in cattle using large-scale pedigree analysis and genome-wide single sperm genomics

PRDM9 contributes to hybrid sterility and species evolution. However, its role is to be confirmed in cattle, a major domesticated livestock species. We previously found an association near PRDM9 with cattle recombination features, but the causative variants are still unknown. Using millions of genotyped cattle with pedigree information, we characterized five PRDM9 alleles and generated allele-specific recombination maps. By examining allele-specific recombination patterns, we observed the impact of PRDM9 on global distribution of recombination, especially in the two ends of chromosomes. We also showed strong associations between recombination hotspot regions and functional mutations within PRDM9 zinc finger domain. More importantly, we found one allele of PRDM9 to be very different from others in both protein composition and recombination landscape, indicating the causative role of this allele on the association between PRDM9 and cattle recombination. When comparing recombination maps from sperm and pedigree data, we observed similar genome-wide recombination patterns, validating the quality of pedigree-based results. Collectively, these evidence supported PRDM9 alleles as causal variants for the reported association with cattle recombination. Our study comprehensively surveyed the bovine PRDM9 alleles, generated allele-specific recombination maps, and expanded our understanding of the role of PRDM9 on genome distribution of recombination.

PRDM9 Methyltransferase Activity Is Essential for Meiotic DNA Double-Strand Break Formation at Its Binding Sites

The programmed formation of hundreds of DNA double-strand breaks (DSBs) is essential for proper meiosis and fertility. In mice and humans, the location of these breaks is determined by the meiosis-specific protein PRDM9, through the DNA-binding specificity of its zinc-finger domain. PRDM9 also has methyltransferase activity. Here, we show that this activity is required for H3K4me3 and H3K36me3 deposition and for DSB formation at PRDM9-binding sites. By analyzing mice that express two PRDM9 variants with distinct DNA-binding specificities, we show that each variant generates its own set of H3K4me3 marks independently from the other variant. Altogether, we reveal several basic principles of PRDM9-dependent DSB site determination, in which an excess of sites are designated through PRDM9 binding and subsequent histone methylation, from which a subset is selected for DSB formation.

The PHD finger protein Spp1 has distinct functions in the Set1 and the meiotic DSB formation complexes

Histone H3K4 methylation is a feature of meiotic recombination hotspots shared by many organisms including plants and mammals. Meiotic recombination is initiated by programmed double-strand break (DSB) formation that in budding yeast takes place in gene promoters and is promoted by histone H3K4 di/trimethylation. This histone modification is recognized by Spp1, a PHD finger containing protein that belongs to the conserved histone H3K4 methyltransferase Set1 complex. During meiosis, Spp1 binds H3K4me3 and interacts with a DSB protein, Mer2, to promote DSB formation close to gene promoters. How Set1 complex- and Mer2- related functions of Spp1 are connected is not clear. Here, combining genome-wide localization analyses, biochemical approaches and the use of separation of function mutants, we show that Spp1 is present within two distinct complexes in meiotic cells, the Set1 and the Mer2 complexes. Disrupting the Spp1-Set1 interaction mildly decreases H3K4me3 levels and does not affect meiotic recombination initiation. Conversely, the Spp1-Mer2 interaction is required for normal meiotic recombination initiation, but dispensable for Set1 complex-mediated histone H3K4 methylation. Finally, we provide evidence that Spp1 preserves normal H3K4me3 levels independently of the Set1 complex. We propose a model where Spp1 works in three ways to promote recombination initiation: first by depositing histone H3K4 methylation (Set1 complex), next by “reading” and protecting histone H3K4 methylation, and finally by making the link with the chromosome axis (Mer2-Spp1 complex). This work deciphers the precise roles of Spp1 in meiotic recombination and opens perspectives to study its functions in other organisms where H3K4me3 is also present at recombination hotspots.

Meiotic recombination is a conserved pathway of sexual reproduction that is required to faithfully segregate homologous chromosomes and produce viable gametes. Recombination events between homologous chromosomes are triggered by the programmed formation of DNA breaks, which occur preferentially at places called hotspots. In many organisms, these hotspots are located close to a particular chromatin modification, the methylation of lysine 4 of histone H3 (H3K4me3). It was previously shown in the budding yeast model that one protein, Spp1, plays an important function in this process. We further explored the functional link between Spp1 and its interacting partners, and show that Spp1 shows genetically separable functions, by depositing the H3K4me3 mark on the chromatin, “reading” and protecting it, and linking it to the recombination proteins. We provide evidence that Spp1 is in distinct complexes to perform these functions. This work opens perspectives for understanding the process in other eukaryotes such as mammals, where most of the proteins involved are conserved.

Effects of Demographic History on the Detection of Recombination Hotspots from Linkage Disequilibrium

In some species, meiotic recombination is concentrated in small genomic regions. These "recombination hotspots" leave signatures in fine-scale patterns of linkage disequilibrium, raising the prospect that the genomic landscape of hotspots can be characterized from sequence variation. This approach has led to the inference that hotspots evolve rapidly in some species, but are conserved in others. Historic demographic events, such as population bottlenecks, are known to affect patterns of linkage disequilibrium across the genome, violating population genetic assumptions of this approach. Although such events are prevalent, demographic history is generally ignored when making inferences about the evolution of recombination hotspots. To determine the effect of demography on the detection of recombination hotspots, we use the coalescent to simulate haplotypes with a known recombination landscape. We measure the ability of popular linkage disequilibrium-based programs to detect hotspots across a range of demographic histories, including population bottlenecks, hidden population structure, population expansions, and population contractions. We find that demographic events have the potential to greatly reduce the power and increase the false positive rate of hotspot discovery. Neither the power nor the false positive rate of hotspot detection can be predicted without also knowing the demographic history of the sample. Our results suggest that ignoring demographic history likely overestimates the power to detect hotspots and therefore underestimates the degree of hotspot sharing between species. We suggest strategies for incorporating demographic history into population genetic inferences about recombination hotspots.

The Diverging Routes of BORIS and CTCF: An Interactomic and Phylogenomic Analysis

The CCCTC-binding factor (CTCF) is multi-functional, ubiquitously expressed, and highly conserved from Drosophila to human. It has important roles in transcriptional insulation and the formation of a high-dimensional chromatin structure. CTCF has a paralog called “Brother of Regulator of Imprinted Sites” (BORIS) or “CTCF-like” (CTCFL). It binds DNA at sites similar to those of CTCF. However, the expression profiles of the two proteins are quite different. We investigated the evolutionary trajectories of the two proteins after the duplication event using a phylogenomic and interactomic approach. We find that CTCF has 52 direct interaction partners while CTCFL only has 19. Almost all interactors already existed before the emergence of CTCF and CTCFL. The unique secondary loss of CTCF from several nematodes is paralleled by a loss of two of its interactors, the polycomb repressive complex subunit SuZ12 and the multifunctional transcription factor TYY1. In contrast to earlier studies reporting the absence of BORIS from birds, we present evidence for a multigene synteny block containing CTCFL that is conserved in mammals, reptiles, and several species of birds, indicating that not the entire lineage of birds experienced a loss of CTCFL. Within this synteny block, BORIS and its genomic neighbors seem to be partitioned into two nested chromatin loops. The high expression of SPO11, RAE1, RBM38, and PMEPA1 in male tissues suggests a possible link between CTCFL, meiotic recombination, and fertility-associated phenotypes. Using the 65,700 exomes and the 1000 genomes data, we observed a higher number of intergenic, non-synonymous, and loss-of-function mutations in CTCFL than in CTCF, suggesting a reduced strength of purifying selection, perhaps due to less functional constraint.

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.

Elevated Proportions of Deleterious Genetic Variation in Domestic Animals and Plants

A fraction of genetic variants segregating in any population are deleterious, which negatively impacts individual fitness. The domestication of animals and plants is associated with population bottlenecks and artificial selection, which are predicted to increase the proportion of deleterious variants. However, the extent to which this is a general feature of domestic species is unclear. Here, we examine the effects of domestication on the prevalence of deleterious variation using pooled whole-genome resequencing data from five domestic animal species (dog, pig, rabbit, chicken, and silkworm) and two domestic plant species (rice and soybean) compared with their wild ancestors. We find significantly reduced genetic variation and increased proportion of nonsynonymous amino acid changes in all but one of the domestic species. These differences are observable across a range of allele frequencies, both common and rare. We find proportionally more single nucleotide polymorphisms in highly conserved elements in domestic species and a tendency for domestic species to harbor a higher proportion of changes classified as damaging. Our findings most likely reflect an increased incidence of deleterious variants in domestic species, which is most likely attributable to population bottlenecks that lead to a reduction in the efficacy of selection. An exception to this pattern is displayed by European domestic pigs, which do not show traces of a strong population bottleneck and probably continued to exchange genes with wild boar populations after domestication. The results presented here indicate that an elevated proportion of deleterious variants is a common, but not ubiquitous, feature of domestic species.

The impact of recombination on human mutation load and disease

Recombination promotes genomic integrity among cells and tissues through double-strand break repair, and is critical for gamete formation and fertility through a strict regulation of the molecular mechanisms associated with proper chromosomal disjunction. In humans, congenital defects and recurrent structural abnormalities can be attributed to aberrant meiotic recombination. Moreover, mutations affecting genes involved in recombination pathways are directly linked to pathologies including infertility and cancer. Recombination is among the most prominent mechanism shaping genome variation, and is associated with not only the structuring of genomic variability, but is also tightly linked with the purging of deleterious mutations from populations. Together, these observations highlight the multiple roles of recombination in human genetics: its ability to act as a major force of evolution, its molecular potential to maintain genome repair and integrity in cell division and its mutagenic cost impacting disease evolution.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.