Crossovers are associated with mutation and biased gene conversion at recombination hotspots

Population size interacts with reproductive longevity to shape the germline mutation rate

Mutation rates vary across the tree of life by many orders of magnitude, with fewer mutations occurring each generation in species that reproduce quickly and maintain large effective population sizes. A compelling explanation is that large effective population sizes facilitate selection against weakly deleterious "mutator alleles" such as variants that modulate cell division or interfere with the molecular efficacy of DNA repair. However, while the fidelity of a single cell division largely determines microorganisms' mutation rates, the relationship of the mutation rate to the molecular determinants of DNA damage and repair is more complex in multicellular species with long generation times. Since long generations leave more time for mutations to accrue each generation, we posit that a long generation time likely amplifies the fitness consequences of any damage agent or DNA repair defect that creates extra mutations in the spermatogonia or oocytes. This leads to the counterintuitive prediction that the species with the highest germline mutation rates per generation are also the species with most effective mechanisms for avoiding and repairing mutations in their reproductive cells. Consistent with this, we show that mutation rates in the reproductive cells are inversely correlated with generation time; in contrast, the number of germline mutations that occur during prepuberty development trends weakly upward as generation time increases. Our results parallel recent findings that the longest-lived species have the lowest mutation rates in adult somatic tissues, potentially due to selection to keep the lifetime mutation load below a harmful threshold.

Independent genetic basis of meiotic crossover positioning and interference in domestic pigs

Meiotic crossover patterning shows huge variation within and between chromosomes, individuals, and species, yet the molecular and evolutionary causes and consequences of this variation remain poorly understood. A key step is to understand the genetic architecture of the crossover rate, positioning, and interference to determine if these factors are governed by common or distinct genetic processes. Here, we investigate individual variation in autosomal crossover count, crossover position (measured as both intra-chromosomal shuffling and distance to telomere), and crossover interference in a large breeding population of domestic pigs (N = 82,474 gametes). We show that all traits are heritable in females at the gamete (h2 = 0.07-0.11) and individual mean levels (h2 = 0.08-0.41). In females, crossover count, and interference are strongly associated with RNF212, but crossover positioning is associated with SYCP2, MEI4, and PRDM9. Our results show that crossover positioning and rate/interference are driven by distinct genetic processes in female pigs and have the capacity to evolve independently.

Estimating Gene Conversion Tract Length and Rate From PacBio HiFi Data

Gene conversions are broadly defined as the transfer of genetic material from a "donor" to an "acceptor" sequence and can happen both in meiosis and mitosis. They are a subset of noncrossover (NCO) events and, like crossover (CO) events, gene conversion can generate new combinations of alleles and counteract mutation load by reverting germline mutations through GC-biased gene conversion. Estimating gene conversion rate and the distribution of gene conversion tract lengths remains challenging. We present a new method for estimating tract length, rate, and detection probability of NCO events directly in HiFi PacBio long read data. The method can be used to make inference from sequencing of gametes from a single individual. The method is unbiased even under low single nucleotide variant (SNV) densities and does not necessitate any demographic or evolutionary assumptions. We test the accuracy and robustness of our method using simulated datasets where we vary length of tracts, number of tracts, the genomic SNV density, and levels of correlation between SNV density and NCO event position. Our simulations show that under low SNV densities, like those found in humans, only a minute fraction (∼2%) of NCO events are expected to become visible as gene conversions by moving at least 1 SNV. We finally illustrate our method by applying it to PacBio sequencing data from human sperm.

The recombination landscape of the barn owl, from families to populations

Homologous recombination is a meiotic process that generates diversity along the genome and interacts with all evolutionary forces. Despite its importance, studies of recombination landscapes are lacking due to methodological limitations and limited data. Frequently used approaches include linkage mapping based on familial data that provides sex-specific broad-scale estimates of realized recombination and inferences based on population linkage disequilibrium that reveal a more fine-scale resolution of the recombination landscape, albeit dependent on the effective population size and the selective forces acting on the population. In this study, we use a combination of these 2 methods to elucidate the recombination landscape for the Afro-European barn owl (Tyto alba). We find subtle differences in crossover placement between sexes that lead to differential effective shuffling of alleles. Linkage disequilibrium-based estimates of recombination are concordant with family-based estimates and identify large variation in recombination rates within and among linkage groups. Larger chromosomes show variation in recombination rates, while smaller chromosomes have a universally high rate that shapes the diversity landscape. We find that recombination rates are correlated with gene content, genetic diversity, and GC content. We find no conclusive differences in the recombination landscapes between populations. Overall, this comprehensive analysis enhances our understanding of recombination dynamics, genomic architecture, and sex-specific variation in the barn owl, contributing valuable insights to the broader field of avian genomics.

Patterns of chromosome evolution in ruminants

Studying when and where gross genomic rearrangements occurred during evolution is key to understanding changes in genome structure with functional consequences that might eventually lead to speciation. Here we identified chromosome rearrangements in ruminants, a clade characterized by large chromosome differences. Using 26 genome assemblies, we reconstructed five ancestral karyotypes and classified the rearrangement events occurring in each lineage. With these reconstructions, we then identified evolutionary breakpoints regions (EBRs) and synteny fragments. Ruminant karyotype evolution is characterized by inversions, while interchromosomal rearrangements occurred preferentially in the oldest ancestor of ruminants. We found that EBRs are depleted of protein coding genes, including housekeeping genes. Similarly, EBRs are not enriched in high GC regions, suggesting that meiotic double strand breaks might not be their origin. Overall, our results characterize at fine detail the location of chromosome rearrangements in ruminant evolution and provide new insights into the formation of EBRs.

Population size interacts with reproductive longevity to shape the germline mutation rate

Mutation rates vary across the tree of life by many orders of magnitude, with lower mutation rates in species that reproduce quickly and maintain large effective population sizes. A compelling explanation for this trend is that large effective population sizes facilitate selection against weakly deleterious "mutator alleles" such as variants that interfere with the molecular efficacy of DNA repair. However, in multicellular organisms, the relationship of the mutation rate to DNA repair efficacy is complicated by variation in reproductive age. Long generation times leave more time for mutations to accrue each generation, and late reproduction likely amplifies the fitness consequences of any DNA repair defect that creates extra mutations in the sperm or eggs. Here, we present theoretical and empirical evidence that a long generation time amplifies the strength of selection for low mutation rates in the spermatocytes and oocytes. This leads to the counterintuitive prediction that the species with the highest germline mutation rates per generation are also the species with most effective mechanisms for DNA proofreading and repair in their germ cells. In contrast, species with different generation times accumulate similar mutation loads during embryonic development. Our results parallel recent findings that the longest-lived species have the lowest mutation rates in adult somatic tissues, potentially due to selection to keep the lifetime mutation load below a harmful threshold.

Male-biased recombination at chromosome ends in a songbird revealed by precisely mapping crossover positions

Recombination plays a crucial role in evolution by generating novel haplotypes and disrupting linkage between genes, thereby enhancing the efficiency of selection. Here, we analyze the genomes of 12 great reed warblers (Acrocephalus arundinaceus) in a 3-generation pedigree to identify precise crossover positions along the chromosomes. We located more than 200 crossovers and found that these were highly concentrated toward the telomeric ends of the chromosomes. Apart from this major pattern in the recombination landscape, we found significantly higher frequencies of crossovers in genic compared with intergenic regions, and in exons compared with introns. Moreover, while the number of recombination events was similar between the sexes, the crossovers were located significantly closer to the ends of paternal compared with maternal chromosomes. In conclusion, our study of the great reed warbler revealed substantial variation in crossover frequencies within chromosomes, with a distinct bias toward the sub-telomeric regions, particularly on the paternal side. These findings emphasize the importance of thoroughly screening the entire length of chromosomes to characterize the recombination landscape and uncover potential sex-biases in recombination.

Recurrent selection shapes the genomic landscape of differentiation between a pair of host-specialized haplodiploids that diverged with gene flow

Understanding the genetics of adaptation and speciation is critical for a complete picture of how biodiversity is generated and maintained. Heterogeneous genomic differentiation between diverging taxa is commonly documented, with genomic regions of high differentiation interpreted as resulting from differential gene flow, linked selection and reduced recombination rates. Disentangling the roles of each of these non-exclusive processes in shaping genome-wide patterns of divergence is challenging but will enhance our knowledge of the repeatability of genomic landscapes across taxa. Here, we combine whole-genome resequencing and genome feature data to investigate the processes shaping the genomic landscape of differentiation for a sister-species pair of haplodiploid pine sawflies, Neodiprion lecontei and Neodiprion pinetum. We find genome-wide correlations between genome features and summary statistics are consistent with pervasive linked selection, with patterns of diversity and divergence more consistently predicted by exon density and recombination rate than the neutral mutation rate (approximated by dS). We also find that both global and local patterns of FST, dXY and π provide strong support for recurrent selection as the primary selective process shaping variation across pine sawfly genomes, with some contribution from balancing selection and lineage-specific linked selection. Because inheritance patterns for haplodiploid genomes are analogous to those of sex chromosomes, we hypothesize that haplodiploids may be especially prone to recurrent selection, even if gene flow occurred throughout divergence. Overall, our study helps fill an important taxonomic gap in the genomic landscape literature and contributes to our understanding of the processes that shape genome-wide patterns of genetic variation.

Causes and consequences of a complex recombinational landscape in the ant Cardiocondyla obscurior

Eusocial Hymenoptera have the highest recombination rates among all multicellular animals studied so far, but it is unclear why this is and how this affects the biology of individual species. A high-resolution linkage map for the ant Cardiocondyla obscurior corroborates genome-wide high recombination rates reported for ants (8.1 cM/Mb). However, recombination is locally suppressed in regions that are enriched with TEs, that have strong haplotype divergence, or that show signatures of epistatic selection in C. obscurior The results do not support the hypotheses that high recombination rates are linked to phenotypic plasticity or to modulating selection efficiency. Instead, genetic diversity and the frequency of structural variants correlate positively with local recombination rates, potentially compensating for the low levels of genetic variation expected in haplodiploid social Hymenoptera with low effective population size. Ultimately, the data show that recombination contributes to within-population polymorphism and to the divergence of the lineages within C. obscurior.

High prevalence of PRDM9-independent recombination hotspots in placental mammals

In many mammals, recombination events are concentrated in hotspots directed by a sequence-specific DNA-binding protein named PRDM9. Intriguingly, PRDM9 has been lost several times in vertebrates, and notably among mammals, it has been pseudogenized in the ancestor of canids. In the absence of PRDM9, recombination hotspots tend to occur in promoter-like features such as CpG islands. It has thus been proposed that one role of PRDM9 could be to direct recombination away from PRDM9-independent hotspots. However, the ability of PRDM9 to direct recombination hotspots has been assessed in only a handful of species, and a clear picture of how much recombination occurs outside of PRDM9-directed hotspots in mammals is still lacking. In this study, we derived an estimator of past recombination activity based on signatures of GC-biased gene conversion in substitution patterns. We quantified recombination activity in PRDM9-independent hotspots in 52 species of boreoeutherian mammals. We observe a wide range of recombination rates at these loci: several species (such as mice, humans, some felids, or cetaceans) show a deficit of recombination, while a majority of mammals display a clear peak of recombination. Our results demonstrate that PRDM9-directed and PRDM9-independent hotspots can coexist in mammals and that their coexistence appears to be the rule rather than the exception. Additionally, we show that the location of PRDM9-independent hotspots is relatively more stable than that of PRDM9-directed hotspots, but that PRDM9-independent hotspots nevertheless evolve slowly in concert with DNA hypomethylation.

Mutation and evolution: Conceptual possibilities

Although random mutation is central to models of evolutionary change, a lack of clarity remains regarding the conceptual possibilities for thinking about the nature and role of mutation in evolution. We distinguish several claims at the intersection of mutation, evolution, and directionality and then characterize a previously unrecognized category: complex conditioned mutation. Empirical evidence in support of this category suggests that the historically famous fluctuation test should be revisited, and new experiments should be undertaken with emerging experimental techniques to facilitate detecting mutation rates within specific loci at an ultra-high, individual base pair resolution.

Exploring FGFR3 Mutations in the Male Germline: Implications for Clonal Germline Expansions and Paternal Age-Related Dysplasias

Delayed fatherhood results in a higher risk of inheriting a new germline mutation that might result in a congenital disorder in the offspring. In particular, some FGFR3 mutations increase in frequency with age, but there are still a large number of uncharacterized FGFR3 mutations that could be expanding in the male germline with potentially early- or late-onset effects in the offspring. Here, we used digital polymerase chain reaction to assess the frequency and spatial distribution of 10 different FGFR3 missense substitutions in the sexually mature male germline. Our functional assessment of the receptor signaling of the variants with biophysical methods showed that 9 of these variants resulted in a higher activation of the receptor´s downstream signaling, resulting in 2 different expansion behaviors. Variants that form larger subclonal expansions in a dissected postmortem testis also showed a positive correlation of the substitution frequency with the sperm donor's age, and a high and ligand-independent FGFR3 activation. In contrast, variants that measured high FGFR3 signaling and elevated substitution frequencies independent of the donor's age did not result in measurable subclonal expansions in the testis. This suggests that promiscuous signal activation might also result in an accumulation of mutations before the sexual maturation of the male gonad with clones staying relatively constant in size throughout time. Collectively, these results provide novel insights into our understanding of the mutagenesis of driver mutations and their resulting mosaicism in the male germline with important consequences for the transmission and recurrence of associated disorders.

Exploring the Micro-Mosaic Landscape of FGFR3 Mutations in the Ageing Male Germline and Their Potential Implications in Meiotic Differentiation

Advanced paternal age increases the risk of transmitting de novo germline mutations, particularly missense mutations activating the receptor tyrosine kinase (RTK) signalling pathway, as exemplified by the FGFR3 mutation, which is linked to achondroplasia (ACH). This risk is attributed to the expansion of spermatogonial stem cells carrying the mutation, forming sub-clonal clusters in the ageing testis, thereby increasing the frequency of mutant sperm and the number of affected offspring from older fathers. While prior studies proposed a correlation between sub-clonal cluster expansion in the testis and elevated mutant sperm production in older donors, limited data exist on the universality of this phenomenon. Our study addresses this gap by examining the testis-expansion patterns, as well as the increases in mutations in sperm for two FGFR3 variants-c.1138G>A (p.G380R) and c.1948A>G (p.K650E)-which are associated with ACH or thanatophoric dysplasia (TDII), respectively. Unlike the ACH mutation, which showed sub-clonal expansion events in an aged testis and a significant increase in mutant sperm with the donor's age, as also reported in other studies, the TDII mutation showed focal mutation pockets in the testis but exhibited reduced transmission into sperm and no significant age-related increase. The mechanism behind this divergence remains unclear, suggesting potential pleiotropic effects of aberrant RTK signalling in the male germline, possibly hindering differentiation requiring meiosis. This study provides further insights into the transmission risks of micro-mosaics associated with advanced paternal age in the male germline.

Genomic Architecture Predicts Tree Topology, Population Structuring, and Demographic History in Amazonian Birds

Geographic barriers are frequently invoked to explain genetic structuring across the landscape. However, inferences on the spatial and temporal origins of population variation have been largely limited to evolutionary neutral models, ignoring the potential role of natural selection and intrinsic genomic processes known as genomic architecture in producing heterogeneity in differentiation across the genome. To test how variation in genomic characteristics (e.g. recombination rate) impacts our ability to reconstruct general patterns of differentiation between species that cooccur across geographic barriers, we sequenced the whole genomes of multiple bird populations that are distributed across rivers in southeastern Amazonia. We found that phylogenetic relationships within species and demographic parameters varied across the genome in predictable ways. Genetic diversity was positively associated with recombination rate and negatively associated with species tree support. Gene flow was less pervasive in genomic regions of low recombination, making these windows more likely to retain patterns of population structuring that matched the species tree. We further found that approximately a third of the genome showed evidence of selective sweeps and linked selection, skewing genome-wide estimates of effective population sizes and gene flow between populations toward lower values. In sum, we showed that the effects of intrinsic genomic characteristics and selection can be disentangled from neutral processes to elucidate spatial patterns of population differentiation.

Meiotic DNA breaks drive multifaceted mutagenesis in the human germ line

Meiotic recombination commences with hundreds of programmed DNA breaks; however, the degree to which they are accurately repaired remains poorly understood. We report that meiotic break repair is eightfold more mutagenic for single-base substitutions than was previously understood, leading to de novo mutation in one in four sperm and one in 12 eggs. Its impact on indels and structural variants is even higher, with 100- to 1300-fold increases in rates per break. We uncovered new mutational signatures and footprints relative to break sites, which implicate unexpected biochemical processes and error-prone DNA repair mechanisms, including translesion synthesis and end joining in meiotic break repair. We provide evidence that these mechanisms drive mutagenesis in human germ lines and lead to disruption of hundreds of genes genome wide.

Homologous Recombination and Repair Functions Required for Mutagenicity during Yeast Meiosis

Several meiotic events reshape the genome prior to its transfer (via gametes) to the next generation. The occurrence of new meiotic mutations is tightly linked to homologous recombination (HR) and firmly depends on Spo11-induced DNA breaks. To gain insight into the molecular mechanisms governing mutagenicity during meiosis, we examined the timing of mutation and recombination events in cells deficient in various DNA HR-repair genes, which represent distinct functions along the meiotic recombination process. Despite sequence similarities and overlapping activities of the two DNA translocases, Rad54 and Tid1, we observed essential differences in their roles in meiotic mutation occurrence: in the absence of Rad54, meiotic mutagenicity was elevated 8-fold compared to the wild type (WT), while in the tid1Δ mutant, there were few meiotic mutations, nine percent compared to the WT. We propose that the presence of Rad54 channels recombinational repair to a less mutagenic pathway, whereas repair assisted by Tid1 is more mutagenic. A 3.5-fold increase in mutation level was observed in dmc1∆ cells, suggesting that single-stranded DNA (ssDNA) may be a potential source for mutagenicity during meiosis. Taken together, we suggest that the introduction of de novo mutations also contributes to the diversification role of meiotic recombination. These rare meiotic mutations revise genomic sequences and may contribute to long-term evolutionary changes.

Sexual dimorphic regulation of recombination by the synaptonemal complex in C. elegans

In sexually reproducing organisms, germ cells faithfully transmit the genome to the next generation by forming haploid gametes, such as eggs and sperm. Although most meiotic proteins are conserved between eggs and sperm, many aspects of meiosis are sexually dimorphic, including the regulation of recombination. The synaptonemal complex (SC), a large ladder-like structure that forms between homologous chromosomes, is essential for regulating meiotic chromosome organization and promoting recombination. To assess whether sex-specific differences in the SC underpin sexually dimorphic aspects of meiosis, we examined Caenorhabditis elegans SC central region proteins (known as SYP proteins) in oogenesis and spermatogenesis and uncovered sex-specific roles for the SYPs in regulating meiotic recombination. We find that SC composition, specifically SYP-2, SYP-3, SYP-5, and SYP-6, is regulated by sex-specific mechanisms throughout meiotic prophase I. During pachytene, both oocytes and spermatocytes differentially regulate the stability of SYP-2 and SYP-3 within an assembled SC. Further, we uncover that the relative amount of SYP-2 and SYP-3 within the SC is independently regulated in both a sex-specific and a recombination-dependent manner. Specifically, we find that SYP-2 regulates the early steps of recombination in both sexes, while SYP-3 controls the timing and positioning of crossover recombination events across the genomic landscape in only oocytes. Finally, we find that SYP-2 and SYP-3 dosage can influence the composition of the other SYPs in the SC via sex-specific mechanisms during pachytene. Taken together, we demonstrate dosage-dependent regulation of individual SC components with sex-specific functions in recombination. These sexual dimorphic features of the SC provide insights into how spermatogenesis and oogenesis adapted similar chromosome structures to differentially regulate and execute recombination.

Effect of recombination on genetic diversity of Caenorhabditis elegans

Greater molecular divergence and genetic diversity are present in regions of high recombination in many species. Studies describing the correlation between variant abundance and recombination rate have long focused on recombination in the context of linked selection models, whereby interference between linked sites under positive or negative selection reduces genetic diversity in regions of low recombination. Here, we show that indels, especially those of intermediate sizes, are enriched relative to single nucleotide polymorphisms in regions of high recombination in C. elegans. To explain this phenomenon, we reintroduce an alternative model that emphasizes the mutagenic effect of recombination. To extend the analysis, we examine the variants with a phylogenetic context and discuss how different models could be examined together. The number of variants generated by recombination in natural populations could be substantial including possibly the majority of some indel subtypes. Our work highlights the potential importance of a mutagenic effect of recombination, which could have a significant role in the shaping of natural genetic diversity.

Genomic diversity landscapes in outcrossing and selfing Caenorhabditis nematodes

Caenorhabditis nematodes form an excellent model for studying how the mode of reproduction affects genetic diversity, as some species reproduce via outcrossing whereas others can self-fertilize. Currently, chromosome-level patterns of diversity and recombination are only available for self-reproducing Caenorhabditis, making the generality of genomic patterns across the genus unclear given the profound potential influence of reproductive mode. Here we present a whole-genome diversity landscape, coupled with a new genetic map, for the outcrossing nematode C. remanei. We demonstrate that the genomic distribution of recombination in C. remanei, like the model nematode C. elegans, shows high recombination rates on chromosome arms and low rates toward the central regions. Patterns of genetic variation across the genome are also similar between these species, but differ dramatically in scale, being tenfold greater for C. remanei. Historical reconstructions of variation in effective population size over the past million generations echo this difference in polymorphism. Evolutionary simulations demonstrate how selection, recombination, mutation, and selfing shape variation along the genome, and that multiple drivers can produce patterns similar to those observed in natural populations. The results illustrate how genome organization and selection play a crucial role in shaping the genomic pattern of diversity whereas demographic processes scale the level of diversity across the genome as a whole.

Modeling Recombination Rate as a Quantitative Trait Reveals New Insight into Selection in Humans

Meiotic recombination is both a fundamental biological process required for proper chromosomal segregation during meiosis and an important genomic parameter that shapes major features of the genomic landscape. However, despite the central importance of this phenotype, we lack a clear understanding of the selective pressures that shape its variation in natural populations, including humans. While there is strong evidence of fitness costs of low rates of recombination, the possible fitness costs of high rates of recombination are less defined. To determine whether a single lower fitness bound can explain the variation in recombination rates observed in human populations, we simulated the evolution of recombination rates as a sexually dimorphic quantitative trait. Under each scenario, we statistically compared the resulting trait distribution with the observed distribution of recombination rates from a published study of the Icelandic population. To capture the genetic architecture of recombination rates in humans, we modeled it as a moderately complex trait with modest heritability. For our fitness function, we implemented a hyperbolic tangent curve with several flexible parameters to capture a wide range of existing hypotheses. We found that costs of low rates of recombination alone are likely insufficient to explain the current variation in recombination rates in both males and females, supporting the existence of fitness costs of high rates of recombination in humans. With simulations using both upper and lower fitness boundaries, we describe a parameter space for the costs of high recombination rates that produces results consistent with empirical observations.

The roles of sexual selection and sexual conflict in shaping patterns of genome and transcriptome variation

Sexual dimorphism is one of the most prevalent, and often the most extreme, examples of phenotypic variation within species, and arises primarily from genomic variation that is shared between females and males. Many sexual dimorphisms arise through sex differences in gene expression, and sex-biased expression is one way that a single, shared genome can generate multiple, distinct phenotypes. Although many sexual dimorphisms are expected to result from sexual selection, and many studies have invoked the possible role of sexual selection to explain sex-specific traits, the role of sexual selection in the evolution of sexually dimorphic gene expression remains difficult to differentiate from other forms of sex-specific selection. In this Review, we propose a holistic framework for the study of sex-specific selection and transcriptome evolution. We advocate for a comparative approach, across tissues, developmental stages and species, which incorporates an understanding of the molecular mechanisms, including genomic variation and structure, governing gene expression. Such an approach is expected to yield substantial insights into the evolution of genetic variation and have important applications in a variety of fields, including ecology, evolution and behaviour.

Pervasive incomplete lineage sorting illuminates speciation and selection in primates

Incomplete lineage sorting (ILS) causes the phylogeny of some parts of the genome to differ from the species tree. In this work, we investigate the frequencies and determinants of ILS in 29 major ancestral nodes across the entire primate phylogeny. We find up to 64% of the genome affected by ILS at individual nodes. We exploit ILS to reconstruct speciation times and ancestral population sizes. Estimated speciation times are much more recent than genomic divergence times and are in good agreement with the fossil record. We show extensive variation of ILS along the genome, mainly driven by recombination but also by the distance to genes, highlighting a major impact of selection on variation along the genome. In many nodes, ILS is reduced more on the X chromosome compared with autosomes than expected under neutrality, which suggests higher impacts of natural selection on the X chromosome. Finally, we show an excess of ILS in genes with immune functions and a deficit of ILS in housekeeping genes. The extensive ILS in primates discovered in this study provides insights into the speciation times, ancestral population sizes, and patterns of natural selection that shape primate evolution.

Repeat-Induced Point Mutation and Gene Conversion Coinciding with Heterochromatin Shape the Genome of a Plant-Pathogenic Fungus

Meiosis is associated with genetic changes in the genome-via recombination, gene conversion, and mutations. The occurrence of gene conversion and mutations during meiosis may further be influenced by the chromatin conformation, similar to the effect of the chromatin conformation on the mitotic mutation rate. To date, however, the exact distribution and type of meiosis-associated changes and the role of the chromatin conformation in this context are largely unexplored. Here, we determine recombination, gene conversion, and de novo mutations using whole-genome sequencing of all meiotic products of 23 individual meioses in Zymoseptoria tritici, an important pathogen of wheat. We confirm a high genome-wide recombination rate of 65 centimorgan (cM)/Mb and see higher recombination rates on the accessory compared to core chromosomes. A substantial fraction of 0.16% of all polymorphic markers was affected by gene conversions, showing a weak GC-bias and occurring at higher frequency in regions of constitutive heterochromatin, indicated by the histone modification H3K9me3. The de novo mutation rate associated with meiosis was approximately three orders of magnitude higher than the corresponding mitotic mutation rate. Importantly, repeat-induced point mutation (RIP), a fungal defense mechanism against duplicated sequences, is active in Z. tritici and responsible for the majority of these de novo meiotic mutations. Our results indicate that the genetic changes associated with meiosis are a major source of variability in the genome of an important plant pathogen and shape its evolutionary trajectory. IMPORTANCE The impact of meiosis on the genome composition via gene conversion and mutations is mostly poorly understood, in particular, for non-model species. Here, we sequenced all four meiotic products for 23 individual meioses and determined the genetic changes caused by meiosis for the important fungal wheat pathogen Zymoseptoria tritici. We found a high rate of gene conversions and an effect of the chromatin conformation on gene conversion rates. Higher conversion rates were found in regions enriched with the H3K9me3-a mark for constitutive heterochromatin. Most importantly, meiosis was associated with a much higher frequency of de novo mutations than mitosis; 78% of the meiotic mutations were caused by repeat-induced point mutations-a fungal defense mechanism against duplicated sequences. In conclusion, the genetic changes associated with meiosis are therefore a major factor shaping the genome of this fungal pathogen.

High-Resolution Genomic Profiling of Liver Cancer Links Etiology With Mutation and Epigenetic Signatures

BACKGROUND & AIMS

Hepatocellular carcinoma (HCC) is a model of a diverse spectrum of cancers because it is induced by well-known etiologies, mainly hepatitis C virus (HCV) and hepatitis B virus. Here, we aimed to identify HCV-specific mutational signatures and explored the link between the HCV-related regional variation in mutations rates and HCV-induced alterations in genome-wide chromatin organization.

METHODS

To identify an HCV-specific mutational signature in HCC, we performed high-resolution targeted sequencing to detect passenger mutations on 64 HCC samples from 3 etiology groups: hepatitis B virus, HCV, or other. To explore the link between the genomic signature and genome-wide chromatin organization we performed chromatin immunoprecipitation sequencing for the transcriptionally permissive H3K4Me3, H3K9Ac, and suppressive H3K9Me3 modifications after HCV infection.

RESULTS

Regional variation in mutation rate analysis showed significant etiology-dependent regional mutation rates in 12 genes: LRP2, KRT84, TMEM132B, DOCK2, DMD, INADL, JAK2, DNAH6, MTMR9, ATM, SLX4, and ARSD. We found an enrichment of C->T transversion mutations in the HCV-associated HCC cases. Furthermore, these cases showed regional variation in mutation rates associated with genomic intervals in which HCV infection dictated epigenetic alterations. This signature may be related to the HCV-induced decreased expression of genes encoding key enzymes in the base excision repair pathway.

CONCLUSIONS

We identified novel distinct HCV etiology-dependent mutation signatures in HCC associated with HCV-induced alterations in histone modification. This study presents a link between cancer-causing mutagenesis and the increased predisposition to liver cancer in chronic HCV-infected individuals, and unveils novel etiology-specific mechanisms leading to HCC and cancer in general.

Crossover patterning in plants

KEY MESSAGE

Chromatin state, and dynamic loading of pro-crossover protein HEI10 at recombination intermediates shape meiotic chromosome patterning in plants. Meiosis is the basis of sexual reproduction, and its basic progression is conserved across eukaryote kingdoms. A key feature of meiosis is the formation of crossovers which result in the reciprocal exchange of segments of maternal and paternal chromosomes. This exchange generates chromosomes with new combinations of alleles, increasing the efficiency of both natural and artificial selection. Crossovers also form a physical link between homologous chromosomes at metaphase I which is critical for accurate chromosome segregation and fertility. The patterning of crossovers along the length of chromosomes is a highly regulated process, and our current understanding of its regulation forms the focus of this review. At the global scale, crossover patterning in plants is largely governed by the classically observed phenomena of crossover interference, crossover homeostasis and the obligatory crossover which regulate the total number of crossovers and their relative spacing. The molecular actors behind these phenomena have long remained obscure, but recent studies in plants implicate HEI10 and ZYP1 as key players in their coordination. In addition to these broad forces, a wealth of recent studies has highlighted how genomic and epigenomic features shape crossover formation at both chromosomal and local scales, revealing that crossovers are primarily located in open chromatin associated with gene promoters and terminators with low nucleosome occupancy.

Why does the X chromosome lag behind autosomes in GWAS findings?

The X-chromosome is among the largest human chromosomes. It differs from autosomes by a number of important features including hemizygosity in males, an almost complete inactivation of one copy in females, and unique patterns of recombination. We used data from the Catalog of Published Genome Wide Association Studies to compare densities of the GWAS-detected SNPs on the X-chromosome and autosomes. The density of GWAS-detected SNPs on the X-chromosome is 6-fold lower compared to the density of the GWAS-detected SNPs on autosomes. Differences between the X-chromosome and autosomes cannot be explained by differences in the overall SNP density, lower X-chromosome coverage by genotyping platforms or low call rate of X-chromosomal SNPs. Similar differences in the density of GWAS-detected SNPs were found in female-only GWASs (e.g. ovarian cancer GWASs). We hypothesized that the lower density of GWAS-detected SNPs on the X-chromosome compared to autosomes is not a result of a methodological bias, e.g. differences in coverage or call rates, but has a real underlying biological reason-a lower density of functional SNPs on the X-chromosome versus autosomes. This hypothesis is supported by the observation that (i) the overall SNP density of X-chromosome is lower compared to the SNP density on autosomes and that (ii) the density of genic SNPs on the X-chromosome is lower compared to autosomes while densities of intergenic SNPs are similar.

Variation and genetic control of individual recombination rates in Norwegian Red dairy cattle

Meiotic recombination is an important evolutionary mechanism that breaks up linkages between loci and creates novel haplotypes for selection to act upon. Understanding the genetic control of variation in recombination rates is therefore of great interest in both natural and domestic breeding populations. In this study, we used pedigree information and medium-density (∼50K) genotyped data in a large cattle (Bos taurus) breeding population in Norway (Norwegian Red cattle) to investigate recombination rate variation between sexes and individual animals. Sex-specific linkage mapping showed higher rates in males than in females (total genetic length of autosomes = 2,492.9 cM in males and 2,308.9 cM in females). However, distribution of recombination along the genome showed little variation between males and females compared with that in other species. The heritability of autosomal crossover count was low but significant in both sexes (h² = 0.04 and 0.09 in males and females, respectively). We identified 2 loci associated with variation in individual crossover counts in female, one close to the candidate gene CEP55 and one close to both MLH3 and NEK9. All 3 genes have been associated with recombination rates in other cattle breeds. Our study contributes to the understanding of how recombination rates are controlled and how they may vary between closely related breeds as well as between species.

Associations between polymorphisms in the myostatin gene with calving difficulty and carcass merit in cattle

A fully functional myostatin gene inhibits muscle fiber growth. The objective of the present study was to quantify the association between 21 known myostatin mutations with both calving and carcass traits in 12 cattle breeds. The myostatin genotypes of 32,770 dam-progeny combinations were used in the association analysis of calving dystocia, with the genotypes of 129,803 animals used in the mixed model association analyses of carcass weight, conformation, and fat score. The mixed model included additive genetic, maternal, and permanent environmental effects where appropriate. The mutant genotypes of nt821, Q204X, and F94L were all associated (P < 0.01) with more calving difficulty when present in either the dam or the progeny. The nt821 deletion had the greatest association with calving difficulty when the homozygous deletion was present in either the calf (0.37 points greater calving difficulty score relative to calves carrying no copies of the deletion based on a one to four scale) or the dam (1.30 points greater calving difficulty score relative to dams carrying no copies of the deletion), although the association between the calf's nt821 genotype and calving difficulty differed depending on the nt821 genotype of the dam. With the exception of nt748_78, nt414, and nt374_51, all other seven segregating myostatin variants were associated (range of allele substitution effect size relative to animals with no copies of the mutant allele) with carcass weight (2.36 kg lighter to 15.56 kg heavier), all 10 segregating variants with conformation (0.15 units less conformed to 2.24 units more conformed assessed on a scale of 1 to 15), and all segregating variants other than E226X with carcass fat (0.23 units less carcass fat cover to 3.85 units more carcass fat cover assessed on a scale of 1 to 15). Of these, the F94L, Q204X, and nt821 mutations generally had the greatest association with all three carcass traits, giving rise to heavier and more conformed carcasses. Despite the antagonistic genetic relationship between calving difficulty and carcass traits, the nt374_51, F94L, and E226X mutations were all associated with improved carcass merit while having minimal expected consequences on calving difficulty. Thus, animals carrying these mutation(s) may have favorable genetic merit for calving difficulty and carcass merit. Furthermore, depending on the dam genotype, a bull with two copies of the nt821 mutation can produce progeny with improved carcass merit while minimizing calving problems.

Independent Evolution of Sex Chromosomes and Male Pregnancy-Related Genes in Two Seahorse Species

Unlike birds and mammals, many teleosts have homomorphic sex chromosomes, and changes in the chromosome carrying the sex-determining locus, termed "turnovers", are common. Recent turnovers allow studies of several interesting questions. One question is whether the new sex-determining regions evolve to become completely non-recombining, and if so, how and why. Another is whether (as predicted) evolutionary changes that benefit one sex accumulate in the newly sex-linked region. To study these questions, we analyzed the genome sequences of two seahorse species of the Syngnathidae, a fish group in which many species evolved a unique structure, the male brood pouch. We find that both seahorse species have XY sex chromosome systems, but their sex chromosome pairs are not homologs, implying that at least one turnover event has occurred. The Y-linked regions occupy 63.9% and 95.1% of the entire sex chromosome of the two species and do not exhibit extensive sequence divergence with their X-linked homologs. We find evidence for occasional recombination between the extant sex chromosomes that may account for their homomorphism. We argue that these Y-linked regions did not evolve by recombination suppression after the turnover, but by the ancestral nature of the low crossover rates in these chromosome regions. With such an ancestral crossover landscape, a turnover can instantly create an extensive Y-linked region. Finally, we test for adaptive evolution of male pouch-related genes after they became Y-linked in the seahorse.

Rapid Evolution of the Fine-scale Recombination Landscape in Wild House Mouse (Mus musculus) Populations

Meiotic recombination is an important evolutionary force and an essential meiotic process. In many species, recombination events concentrate into hotspots defined by the site-specific binding of PRMD9. Rapid evolution of Prdm9's zinc finger DNA-binding array leads to remarkably abrupt shifts in the genomic distribution of hotspots between species, but the question of how Prdm9 allelic variation shapes the landscape of recombination between populations remains less well understood. Wild house mice (Mus musculus) harbor exceptional Prdm9 diversity, with >150 alleles identified to date, and pose a particularly powerful system for addressing this open question. We employed a coalescent-based approach to construct broad- and fine-scale sex-averaged recombination maps from contemporary patterns of linkage disequilibrium in nine geographically isolated wild house mouse populations, including multiple populations from each of three subspecies. Comparing maps between wild mouse populations and subspecies reveals several themes. First, we report weak fine- and broad-scale recombination map conservation across subspecies and populations, with genetic divergence offering no clear prediction for recombination map divergence. Second, most hotspots are unique to one population, an outcome consistent with minimal sharing of Prdm9 alleles between surveyed populations. Finally, by contrasting aggregate hotspot activity on the X versus autosomes, we uncover evidence for population-specific differences in the degree and direction of sex dimorphism for recombination. Overall, our findings illuminate the variability of both the broad- and fine-scale recombination landscape in M. musculus and underscore the functional impact of Prdm9 allelic variation in wild mouse populations.

Adaptive sequence divergence forged new neurodevelopmental enhancers in humans

Searches for the genetic underpinnings of uniquely human traits have focused on human-specific divergence in conserved genomic regions, which reflects adaptive modifications of existing functional elements. However, the study of conserved regions excludes functional elements that descended from previously neutral regions. Here, we demonstrate that the fastest-evolved regions of the human genome, which we term "human ancestor quickly evolved regions" (HAQERs), rapidly diverged in an episodic burst of directional positive selection prior to the human-Neanderthal split, before transitioning to constraint within hominins. HAQERs are enriched for bivalent chromatin states, particularly in gastrointestinal and neurodevelopmental tissues, and genetic variants linked to neurodevelopmental disease. We developed a multiplex, single-cell in vivo enhancer assay to discover that rapid sequence divergence in HAQERs generated hominin-unique enhancers in the developing cerebral cortex. We propose that a lack of pleiotropic constraints and elevated mutation rates poised HAQERs for rapid adaptation and subsequent susceptibility to disease.

The Genomic Landscapes of Desert Birds Form over Multiple Time Scales

Spatial models show that genetic differentiation between populations can be explained by factors ranging from geographic distance to environmental resistance across the landscape. However, genomes exhibit a landscape of differentiation, indicating that multiple processes may mediate divergence in different portions of the genome. We tested this idea by comparing alternative geographic predctors of differentiation in ten bird species that co-occur in Sonoran and Chihuahuan Deserts of North America. Using population-level genomic data, we described the genomic landscapes across species and modeled conditions that represented historical and contemporary mechanisms. The characteristics of genomic landscapes differed across species, influenced by varying levels of population structuring and admixture between deserts, and the best-fit models contrasted between the whole genome and partitions along the genome. Both historical and contemporary mechanisms were important in explaining genetic distance, but particularly past and current environments, suggesting that genomic evolution was modulated by climate and habitat There were also different best-ftit models across genomic partitions of the data, indicating that these regions capture different evolutionary histories. These results show that the genomic landscape of differentiation can be associated with alternative geographic factors operating on different portions of the genome, which reflect how heterogeneous patterns of genetic differentiation can evolve across species and genomes.

Some complexities in interpreting apparent effects of hitchhiking: A commentary on Gompert et al. (2022)

We write to address recent claims by regarding the potentially important and underappreciated phenomena of "indirect selection," the observation that neutral regions may be affected by natural selection. We argue both that this phenomenon-generally known as genetic hitchhiking-is neither new nor poorly studied, and that the patterns described by the authors have multiple alternative explanations.

Natural and artificial sources of genetic variation used in crop breeding: A baseline comparator for genome editing

Traditional breeding has successfully selected beneficial traits for food, feed, and fibre crops over the last several thousand years. The last century has seen significant technological advancements particularly in marker assisted selection and the generation of induced genetic variation, including over the last few decades, through mutation breeding, genetic modification, and genome editing. While regulatory frameworks for traditional varietal development and for genetic modification with transgenes are broadly established, those for genome editing are lacking or are still evolving in many regions. In particular, the lack of "foreign" recombinant DNA in genome edited plants and that the resulting SNPs or INDELs are indistinguishable from those seen in traditional breeding has challenged development of new legislation. Where products of genome editing and other novel breeding technologies possess no transgenes and could have been generated via traditional methods, we argue that it is logical and proportionate to apply equivalent legislative oversight that already exists for traditional breeding and novel foods. This review analyses the types and the scale of spontaneous and induced genetic variation that can be selected during traditional plant breeding activities. It provides a base line from which to judge whether genetic changes brought about by techniques of genome editing or other reverse genetic methods are indeed comparable to those routinely found using traditional methods of plant breeding.

Enabling Inference for Context-Dependent Models of Mutation by Bounding the Propagation of Dependency

Although the rates at which positions in the genome mutate are known to depend not only on the nucleotide to be mutated, but also on neighboring nucleotides, it remains challenging to do phylogenetic inference using models of context-dependent mutation. In these models, the effects of one mutation may in principle propagate to faraway locations, making it difficult to compute exact likelihoods. This article shows how to use bounds on the propagation of dependency to compute likelihoods of mutation of a given segment of genome by marginalizing over sufficiently long flanking sequence. This can be used for maximum likelihood or Bayesian inference. Protocols examining residuals and iterative model refinement are also discussed. Tools for efficiently working with these models are provided in an R package, which could be used in other applications. The method is used to examine context dependence of mutations since the common ancestor of humans and chimpanzee.

How Can We Resolve Lewontin's Paradox?

We discuss the genetic, demographic, and selective forces that are likely to be at play in restricting observed levels of DNA sequence variation in natural populations to a much smaller range of values than would be expected from the distribution of census population sizes alone-Lewontin's Paradox. While several processes that have previously been strongly emphasized must be involved, including the effects of direct selection and genetic hitchhiking, it seems unlikely that they are sufficient to explain this observation without contributions from other factors. We highlight a potentially important role for the less-appreciated contribution of population size change; specifically, the likelihood that many species and populations may be quite far from reaching the relatively high equilibrium diversity values that would be expected given their current census sizes.

The megabase-scale crossover landscape is largely independent of sequence divergence

Meiotic recombination frequency varies along chromosomes and strongly correlates with sequence divergence. However, the causal relationship between recombination landscapes and polymorphisms is unclear. Here, we characterize the genome-wide recombination landscape in the quasi-absence of polymorphisms, using Arabidopsis thaliana homozygous inbred lines in which a few hundred genetic markers were introduced through mutagenesis. We find that megabase-scale recombination landscapes in inbred lines are strikingly similar to the recombination landscapes in hybrids, with the notable exception of heterozygous large rearrangements where recombination is prevented locally. In addition, the megabase-scale recombination landscape can be largely explained by chromatin features. Our results show that polymorphisms are not a major determinant of the shape of the megabase-scale recombination landscape but rather favour alternative models in which recombination and chromatin shape sequence divergence across the genome.

The frequency of recombination varies along chromosomes and highly correlates with sequence divergence. Here, the authors show that polymorphisms are not a major determinant of the megabase-scale recombination landscape in Arabidopsis, which is rather determined by chromatin accessibility and DNA methylation.

On the prospect of achieving accurate joint estimation of selection with population history

As both natural selection and population history can affect genome-wide patterns of variation, disentangling the contributions of each has remained as a major challenge in population genetics. We here discuss historical and recent progress towards this goal-highlighting theoretical and computational challenges that remain to be addressed, as well as inherent difficulties in dealing with model complexity and model violations-and offer thoughts on potentially fruitful next steps.

Extreme variation in recombination rate and genetic diversity along the Sylvioidea neo-sex chromosome

Recombination strongly impacts sequence evolution by affecting the extent of linkage and the efficiency of selection. Here, we study recombination over the Z chromosome in great reed warblers (Acrocephalus arundinaceus) using pedigree-based linkage mapping. This species has extended Z and W chromosomes ("neo-sex chromosomes") formed by a fusion between a part of chromosome 4A and the ancestral sex chromosomes, which provides a unique opportunity to assess recombination and sequence evolution in sex-linked regions of different ages. We assembled an 87.54 Mbp and 90.19 cM large Z with a small pseudoautosomal region (0.89 Mbp) at one end and the fused Chr4A-part at the other end of the chromosome. A prominent feature in our data was an extreme variation in male recombination rate along Z with high values at both chromosome ends, but an apparent lack of recombination over a substantial central section, covering 78% of the chromosome. The nonrecombining region showed a drastic loss of genetic diversity and accumulation of repeats compared to the recombining parts. Thus, our data emphasize a key role of recombination in affecting local levels of polymorphism. Nonetheless, the evolutionary rate of genes (dN/dS) did not differ between high and low recombining regions, suggesting that the efficiency of selection on protein-coding sequences can be maintained also at very low levels of recombination. Finally, the Chr4A-derived part showed a similar recombination rate as the part of the ancestral Z that did recombine, but its sequence characteristics reflected both its previous autosomal, and current Z-linked, recombination patterns.

Chromatin accessibility shapes meiotic recombination in mouse primordial germ cells through assisting double-strand breaks and loop formation

Meiotic recombination is a driver of evolution, and aberrant recombination is a major contributor to aneuploidy in mammals. Mechanism of recombination remains elusive yet. Here, we present a computational analysis to explore recombination-related dynamics of chromatin accessibility in mouse primordial germ cells (PGCs). Our data reveals that: (1) recombination hotspots which get accessible at meiosis-specific DNase I-hypersensitive sites (DHSs) only when PGCs enter meiosis are located preferentially in intronic and distal intergenic regions; (2) stable DHSs maintained stably across PGC differentiation are enriched by CTCF motifs and CTCF binding and mediate chromatin loop formation; (3) compared with the specific DHSs aroused at meiotic stage, stable DHSs are largely encoded in DNA sequence and also enriched by epigenetic marks; (4) PRDM9 is likely to target nucleosome-occupied hotspot regions and remodels local chromatin structure to make them accessible for recombination machinery; and (5) cells undergoing meiotic recombination are deficient in TAD structure and chromatin loop arrays are organized regularly along the axis formed between homologous chromosomes. Taken together, by analyzing DHS-related DNA features, epigenetic marks and 3D genome structure, we revealed some specific roles of chromatin accessibility in recombination, which would expand our understanding of recombination mechanism.

Admixture of evolutionary rates across a butterfly hybrid zone

Hybridization is a major evolutionary force that can erode genetic differentiation between species, whereas reproductive isolation maintains such differentiation. In studying a hybrid zone between the swallowtail butterflies Papilio syfanius and Papilio maackii (Lepidoptera: Papilionidae), we made the unexpected discovery that genomic substitution rates are unequal between the parental species. This phenomenon creates a novel process in hybridization, where genomic regions most affected by gene flow evolve at similar rates between species, while genomic regions with strong reproductive isolation evolve at species-specific rates. Thus, hybridization mixes evolutionary rates in a way similar to its effect on genetic ancestry. Using coalescent theory, we show that the rate-mixing process provides distinct information about levels of gene flow across different parts of genomes, and the degree of rate-mixing can be predicted quantitatively from relative sequence divergence ([Formula: see text]) between the hybridizing species at equilibrium. Overall, we demonstrate that reproductive isolation maintains not only genomic differentiation, but also the rate at which differentiation accumulates. Thus, asymmetric rates of evolution provide an additional signature of loci involved in reproductive isolation.

The Evolution of Widespread Recombination Suppression on the Dwarf Hamster (Phodopus) X Chromosome

The X chromosome of therian mammals shows strong conservation among distantly related species, limiting insights into the distinct selective processes that have shaped sex chromosome evolution. We constructed a chromosome-scale de novo genome assembly for the Siberian dwarf hamster (Phodopus sungorus), a species reported to show extensive recombination suppression across an entire arm of the X chromosome. Combining a physical genome assembly based on shotgun and long-range proximity ligation sequencing with a dense genetic map, we detected widespread suppression of female recombination across ∼65% of the Phodopus X chromosome. This region of suppressed recombination likely corresponds to the Xp arm, which has previously been shown to be highly heterochromatic. Using additional sequencing data from two closely related species (P. campbelli and P. roborovskii), we show that recombination suppression on Xp appears to be independent of major structural rearrangements. The suppressed Xp arm was enriched for several transposable element families and de-enriched for genes primarily expressed in placenta, but otherwise showed similar gene densities, expression patterns, and rates of molecular evolution when compared to the recombinant Xq arm. Phodopus Xp gene content and order was also broadly conserved relative to the more distantly related rat X chromosome. These data suggest that widespread suppression of recombination has likely evolved through the transient induction of facultative heterochromatin on the Phodopus Xp arm without major changes in chromosome structure or genetic content. Thus, substantial changes in the recombination landscape have so far had relatively subtle influences on patterns of X-linked molecular evolution in these species.

Recombination rates in pigs differ between breeds, sexes and individuals, and are associated with the RNF212, SYCP2, PRDM7, MEI1 and MSH4 loci

Background

Recombination is a fundamental part of mammalian meiosis that leads to the exchange of large segments of DNA between homologous chromosomes and is therefore an important driver of genetic diversity in populations. In breeding populations, understanding recombination is of particular interest because it can break up unfavourable linkage phases between alleles and produce novel combinations of alleles that could be exploited in selection. In this study, we used dense single nucleotide polymorphism (SNP) genotype data and pedigree information to analyse individual and sex-specific variation and genetic architecture of recombination rates within and between five commercially selected pig breeds.

Results

In agreement with previous studies, recombination rates were higher in females than in males for all breeds and for all chromosomes, except 1 and 13, for which male rates were slightly higher. Total recombination rate differed between breeds but the pattern of recombination along the chromosomes was well conserved across breeds for the same sex. The autosomal linkage maps spanned a total length of 1731 to 1887 cM for males and of 2231 to 2515 cM for females. Estimates of heritability for individual autosomal crossover count ranged from 0.04 to 0.07 for males and from 0.08 to 0.11 for females. Fourteen genomic regions were found to be associated with individual autosomal crossover count. Of these, four were close to or within candidate genes that have previously been associated with individual recombination rates in pigs and other mammals, namely RNF212, SYCP2 and MSH4. Two of the identified regions included the PRDM7 and MEI1 genes, which are known to be involved in meiosis but have not been previously associated with variation in individual recombination rates.

Conclusions

This study shows that genetic variation in autosomal recombination rate persists in domesticated species under strong selection, with differences between closely-related breeds and marked differences between the sexes. Our findings support results from other studies, i.e., that individual crossover counts are associated with the RNF212, SYCP2 and MSH4 genes in pig. In addition, we have found two novel candidate genes associated with the trait, namely PRDM7 and MEI1.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00723-9.

Familial long-read sequencing increases yield of de novo mutations

Studies of de novo mutation (DNM) have typically excluded some of the most repetitive and complex regions of the genome because these regions cannot be unambiguously mapped with short-read sequencing data. To better understand the genome-wide pattern of DNM, we generated long-read sequence data from an autism parent-child quad with an affected female where no pathogenic variant had been discovered in short-read Illumina sequence data. We deeply sequenced all four individuals by using three sequencing platforms (Illumina, Oxford Nanopore, and Pacific Biosciences) and three complementary technologies (Strand-seq, optical mapping, and 10X Genomics). Using long-read sequencing, we initially discovered and validated 171 DNMs across two children-a 20% increase in the number of de novo single-nucleotide variants (SNVs) and indels when compared to short-read callsets. The number of DNMs further increased by 5% when considering a more complete human reference (T2T-CHM13) because of the recovery of events in regions absent from GRCh38 (e.g., three DNMs in heterochromatic satellites). In total, we validated 195 de novo germline mutations and 23 potential post-zygotic mosaic mutations across both children; the overall true substitution rate based on this integrated callset is at least 1.41 × 10-8 substitutions per nucleotide per generation. We also identified six de novo insertions and deletions in tandem repeats, two of which represent structural variants. We demonstrate that long-read sequencing and assembly, especially when combined with a more complete reference genome, increases the number of DNMs by >25% compared to previous studies, providing a more complete catalog of DNM compared to short-read data alone.

Chromosome size affects sequence divergence between species through the interplay of recombination and selection

The structure of the genome shapes the distribution of genetic diversity and sequence divergence. To investigate how the relationship between chromosome size and recombination rate affects sequence divergence between species, we combined empirical analyses and evolutionary simulations. We estimated pairwise sequence divergence among 15 species from three different mammalian clades-Peromyscus rodents, Mus mice, and great apes-from chromosome-level genome assemblies. We found a strong significant negative correlation between chromosome size and sequence divergence in all species comparisons within the Peromyscus and great apes clades but not the Mus clade, suggesting that the dramatic chromosomal rearrangements among Mus species may have masked the ancestral genomic landscape of divergence in many comparisons. Our evolutionary simulations showed that the main factor determining differences in divergence among chromosomes of different sizes is the interplay of recombination rate and selection, with greater variation in larger populations than in smaller ones. In ancestral populations, shorter chromosomes harbor greater nucleotide diversity. As ancestral populations diverge, diversity present at the onset of the split contributes to greater sequence divergence in shorter chromosomes among daughter species. The combination of empirical data and evolutionary simulations revealed that chromosomal rearrangements, demography, and divergence times may also affect the relationship between chromosome size and divergence, thus deepening our understanding of the role of genome structure in the evolution of species divergence.

Supergene origin and maintenance in Atlantic cod

Supergenes are sets of genes that are inherited as a single marker and encode complex phenotypes through their joint action. They are identified in an increasing number of organisms, yet their origins and evolution remain enigmatic. In Atlantic cod, four megabase-scale supergenes have been identified and linked to migratory lifestyle and environmental adaptations. Here we investigate the origin and maintenance of these four supergenes through analysis of whole-genome-sequencing data, including a new long-read-based genome assembly for a non-migratory Atlantic cod individual. We corroborate the finding that chromosomal inversions underlie all four supergenes, and we show that they originated at different times between 0.40 and 1.66 million years ago. We reveal gene flux between supergene haplotypes where migratory and stationary Atlantic cod co-occur and conclude that this gene flux is driven by gene conversion, on the basis of an increase in GC content in exchanged sites. Additionally, we find evidence for double crossover between supergene haplotypes, leading to the exchange of an ~275 kilobase fragment with genes potentially involved in adaptation to low salinity in the Baltic Sea. Our results suggest that supergenes can be maintained over long timescales in the same way as hybridizing species, through the selective purging of introduced genetic variation.

Atlantic cod carries four supergenes linked to migratory lifestyle and environmental adaptations. Using whole-genome sequencing, the authors show that the genome inversions that underlie the supergenes originated at different times and show gene flux between supergene haplotypes.

Quantitative modelling of fine-scale variations in the Arabidopsis thaliana crossover landscape

In, essentially, all species where meiotic crossovers (COs) have been studied, they occur preferentially in open chromatin, typically near gene promoters and to a lesser extent, at the end of genes. Here, in the case of Arabidopsis thaliana, we unveil further trends arising when one considers contextual information, namely summarised epigenetic status, gene or intergenic region size, and degree of divergence between homologs. For instance, we find that intergenic recombination rate is reduced if those regions are less than 1.5 kb in size. Furthermore, we propose that the presence of single nucleotide polymorphisms enhances the rate of CO formation compared to when homologous sequences are identical, in agreement with previous works comparing rates in adjacent homozygous and heterozygous blocks. Lastly, by integrating these different effects, we produce a quantitative and predictive model of the recombination landscape that reproduces much of the experimental variation.

Genome Fractionation and Loss of Heterozygosity in Hybrids and Polyploids: Mechanisms, Consequences for Selection, and Link to Gene Function

Hybridization and genome duplication have played crucial roles in the evolution of many animal and plant taxa. The subgenomes of parental species undergo considerable changes in hybrids and polyploids, which often selectively eliminate segments of one subgenome. However, the mechanisms underlying these changes are not well understood, particularly when the hybridization is linked with asexual reproduction that opens up unexpected evolutionary pathways. To elucidate this problem, we compared published cytogenetic and RNAseq data with exome sequences of asexual diploid and polyploid hybrids between three fish species; Cobitis elongatoides, C. taenia, and C. tanaitica. Clonal genomes remained generally static at chromosome-scale levels but their heterozygosity gradually deteriorated at the level of individual genes owing to allelic deletions and conversions. Interestingly, the impact of both processes varies among animals and genomic regions depending on ploidy level and the properties of affected genes. Namely, polyploids were more tolerant to deletions than diploid asexuals where conversions prevailed, and genomic restructuring events accumulated preferentially in genes characterized by high transcription levels and GC-content, strong purifying selection and specific functions like interacting with intracellular membranes. Although hybrids were phenotypically more similar to C. taenia, we found that they preferentially retained C. elongatoides alleles. This demonstrates that favored subgenome is not necessarily the transcriptionally dominant one. This study demonstrated that subgenomes in asexual hybrids and polyploids evolve under a complex interplay of selection and several molecular mechanisms whose efficiency depends on the organism's ploidy level, as well as functional properties and parental ancestry of the genomic region.