Variation in recombination frequency and distribution across eukaryotes: patterns and processes

On the estimation of genome-average recombination rates

The rate at which recombination events occur in a population is an indicator of its effective population size and the organism's reproduction mode. It determines the extent of linkage disequilibrium along the genome and, thereby, the efficacy of both purifying and positive selection. The population recombination rate can be inferred using models of genome evolution in populations. Classic methods based on the patterns of linkage disequilibrium provide the most accurate estimates, providing large sample sizes are used and the demography of the population is properly accounted for. Here, the capacity of approaches based on the sequentially Markov coalescent (SMC) to infer the genome-average recombination rate from as little as a single diploid genome is examined. SMC approaches provide highly accurate estimates even in the presence of changing population sizes, providing that (1) within genome heterogeneity is accounted for and (2) classic maximum-likelihood optimization algorithms are employed to fit the model. SMC-based estimates proved sensitive to gene conversion, leading to an overestimation of the recombination rate if conversion events are frequent. Conversely, methods based on the correlation of heterozygosity succeed in disentangling the rate of crossing over from that of gene conversion events, but only when the population size is constant and the recombination landscape homogeneous. These results call for a convergence of these two methods to obtain accurate and comparable estimates of recombination rates between populations.

Fine-scale contemporary recombination variation and its fitness consequences in adaptively diverging stickleback fish

Despite deep evolutionary conservation, recombination rates vary greatly across the genome and among individuals, sexes and populations. Yet the impact of this variation on adaptively diverging populations is not well understood. Here we characterized fine-scale recombination landscapes in an adaptively divergent pair of marine and freshwater populations of threespine stickleback from River Tyne, Scotland. Through whole-genome sequencing of large nuclear families, we identified the genomic locations of almost 50,000 crossovers and built recombination maps for marine, freshwater and hybrid individuals at a resolution of 3.8 kb. We used these maps to quantify the factors driving variation in recombination rates. We found strong heterochiasmy between sexes but also differences in recombination rates among ecotypes. Hybrids showed evidence of significant recombination suppression in overall map length and in individual loci. Recombination rates were lower not only within individual marine-freshwater-adaptive loci, but also between loci on the same chromosome, suggesting selection on linked gene 'cassettes'. Through temporal sampling along a natural hybrid zone, we found that recombinants showed traits associated with reduced fitness. Our results support predictions that divergence in cis-acting recombination modifiers, whose functions are disrupted in hybrids, may play an important role in maintaining differences among adaptively diverging populations.

Temporal challenges in detecting balancing selection from population genomic data

The role of balancing selection in maintaining genetic variation remains an open question in population genetics. Recent years have seen numerous studies identifying candidate loci potentially experiencing balancing selection, most predominantly in human populations. There are however numerous alternative evolutionary processes that may leave similar patterns of variation, thereby potentially confounding inference, and the expected signatures of balancing selection additionally change in a temporal fashion. Here we use forward-in-time simulations to quantify expected statistical power to detect balancing selection using both site frequency spectrum- and linkage disequilibrium-based methods under a variety of evolutionarily realistic null models. We find that whilst site frequency spectrum-based methods have little power immediately after a balanced mutation begins segregating, power increases with time since the introduction of the balanced allele. Conversely, linkage disequilibrium-based methods have considerable power whilst the allele is young, and power dissipates rapidly as the time since introduction increases. Taken together, this suggests that site frequency spectrum-based methods are most effective at detecting long-term balancing selection (>25N generations since the introduction of the balanced allele) whilst linkage disequilibrium-based methods are effective over much shorter timescales (<1N generations), thereby leaving a large time frame over which current methods have little power to detect the action of balancing selection. Finally, we investigate the extent to which alternative evolutionary processes may mimic these patterns, and demonstrate the need for caution in attempting to distinguish the signatures of balancing selection from those of both neutral processes (e.g. population structure and admixture) as well as of alternative selective processes (e.g. partial selective sweeps).

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.

Genetic background affects the strength of crossover interference in house mice

Meiotic recombination is required for faithful chromosome segregation in most sexually reproducing organisms and shapes the distribution of genetic variation in populations. Both the overall rate and the spatial distribution of crossovers vary within and between species. Adjacent crossovers on the same chromosome tend to be spaced more evenly than expected at random, a phenomenon known as crossover interference. Although interference has been observed in many taxa, the factors that influence the strength of interference are not well understood. We used house mice (Mus musculus), a well-established model system for understanding recombination, to study the effects of genetics and age on recombination rate and interference in the male germline. We analyzed crossover positions in 503 progeny from reciprocal F1 hybrids between inbred strains representing the three major subspecies of house mice. Consistent with previous studies, autosomal alleles from M. m. musculus tend to increase recombination rate, while inheriting a M. m. musculus X chromosome decreases recombination rate. Old males transmit an average of 0.6 more crossovers per meiosis (5.0%) than young males, though the effect varies across genetic backgrounds. We show that the strength of crossover interference depends on genotype, providing a rare demonstration that interference evolves over short timescales. Differences between reciprocal F1s suggest that X-linked factors modulate the strength of interference. Our findings motivate additional comparisons of interference among recently diverged species and further examination of the role of paternal age in determining the number and positioning of crossovers.

Limited genomic signatures of population collapse in the critically endangered black abalone (Haliotis cracherodii)

The black abalone, Haliotis cracherodii, is a large, long-lived marine mollusc that inhabits rocky intertidal habitats along the coast of California and Mexico. In 1985, populations were impacted by a bacterial disease known as withering syndrome (WS) that wiped out >90% of individuals, leading to the closure of all U.S. black abalone fisheries since 1993. Current conservation strategies include restoring diminished populations by translocating healthy individuals. However, population collapse on this scale may have dramatically lowered genetic diversity and strengthened geographic differentiation, making translocation-based recovery contentious. Additionally, the current prevalence of WS remains unknown. To address these uncertainties, we sequenced and analysed the genomes of 133 black abalone individuals from across their present range. We observed no spatial genetic structure among black abalone, with the exception of a single chromosomal inversion that increases in frequency with latitude. Outside the inversion, genetic differentiation between sites is minimal and does not scale with either geographic distance or environmental dissimilarity. Genetic diversity appears uniformly high across the range. Demographic inference does indicate a severe population bottleneck beginning just 15 generations in the past, but this decline is short lived, with present-day size far exceeding the pre-bottleneck status quo. Finally, we find the bacterial agent of WS is equally present across the sampled range, but only in 10% of individuals. The lack of population genetic structure, uniform diversity and prevalence of WS bacteria indicates that translocation could be a valid and low-risk means of population restoration for black abalone species' recovery.

Blessing or curse: how the epigenetic resolution of host-transposable element conflicts shapes their evolutionary dynamics

Transposable elements (TEs) are selfish genetic elements whose antagonistic interactions with hosts represent a common genetic conflict in eukaryotes. To resolve this conflict, hosts have widely adopted epigenetic silencing that deposits repressive marks at TEs. However, this mechanism is imperfect and fails to fully halt TE replication. Furthermore, TE epigenetic silencing can inadvertently spread repressive marks to adjacent functional sequences, a phenomenon considered a 'curse' of this conflict resolution. Here, we used forward simulations to explore how TE epigenetic silencing and its harmful side effects shape the evolutionary dynamics of TEs and their hosts. Our findings reveal that epigenetic silencing allows TEs and their hosts to stably coexist under a wide range of conditions, because the underlying molecular mechanisms give rise to copy-number dependency of the strength of TE silencing. Interestingly, contrary to intuitive expectations that TE epigenetic silencing should evolve to be as strong as possible, we found a selective benefit for modifier alleles that weaken TE silencing under biologically feasible conditions. These results reveal that the dual nature of TE epigenetic silencing, with both positive and negative effects, complicates its evolutionary trajectory and makes it challenging to determine whether TE epigenetic silencing is a 'blessing' or a 'curse'.

Counting the genetic ancestors from source populations in members of an admixed population

In a genetically admixed population, admixed individuals possess genealogical and genetic ancestry from multiple source groups. Under a mechanistic model of admixture, we study the number of distinct ancestors from the source populations that the admixture represents. Combining a mechanistic admixture model with a recombination model that describes the probability that a genealogical ancestor is a genetic ancestor, for a member of a genetically admixed population, we count genetic ancestors from the source populations-those genealogical ancestors from the source populations who contribute to the genome of the modern admixed individual. We compare patterns in the numbers of genealogical and genetic ancestors across the generations. To illustrate the enumeration of genetic ancestors from source populations in an admixed group, we apply the model to the African-American population, extending recent results on the numbers of African and European genealogical ancestors that contribute to the pedigree of an African-American chosen at random, so that we also evaluate the numbers of African and European genetic ancestors who contribute to random African-American genomes. The model suggests that the autosomal genome of a random African-American born in the interval 1960-1965 contains genetic contributions from a mean of 162 African (standard deviation 47, interquartile range 127-192) and 32 European ancestors (standard deviation 14, interquartile range 21-43). The enumeration of genetic ancestors can potentially be performed in other diploid species in which admixture and recombination models can be specified.

The pseudoentropy of allele frequency trajectories, the persistence of variation, and the effective population size

To concisely describe how genetic variation, at individual loci or across whole genomes, changes over time, and to follow transitory allelic changes, we introduce a quantity related to entropy, that we term pseudoentropy. This quantity emerges in a diffusion analysis of the mean time a mutation segregates in a population. For a neutral locus with an arbitrary number of alleles, the mean time of segregation is generally proportional to the pseudoentropy of initial allele frequencies. After the initial time point, pseudoentropy generally decreases, but other behaviours are possible, depending on the genetic diversity and selective forces present. For a biallelic locus, pseudoentropy and entropy coincide, but they are distinct quantities with more than two alleles. Thus for populations with multiple biallelic loci, the language of entropy suffices. Then entropy, combined across loci, serves as a concise description of genetic variation. We used individual based simulations to explore how this entropy behaves under different evolutionary scenarios. In agreement with predictions, the entropy associated with unlinked neutral loci decreases over time. However, deviations from free recombination and neutrality have clear and informative effects on the entropy's behaviour over time. Analysis of publicly available data of a natural D. melanogaster population, that had been sampled over seven years, using a sliding-window approach, yielded considerable variation in entropy trajectories of different genomic regions. These mostly follow a pattern that suggests a substantial effective population size and a limited effect of positive selection on genome-wide diversity over short time scales.

Metazoa-level USCOs as markers in species delimitation and classification

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

rec-1 loss of function increases recombination in the central gene clusters at the expense of autosomal pairing centers

Meiotic control of crossover (CO) number and position is critical for homologous chromosome segregation and organismal fertility, recombination of parental genotypes, and the generation of novel genetic combinations. We here characterize the recombination rate landscape of a rec-1 loss of function modifier of CO position in Caenorhabditis elegans, one of the first ever modifiers discovered. By averaging CO position across hermaphrodite and male meioses and by genotyping 203 single-nucleotide variants covering about 95% of the genome, we find that the characteristic chromosomal arm-center recombination rate domain structure is lost in the loss of function rec-1 mutant. The rec-1 loss of function mutant smooths the recombination rate landscape but is insufficient to eliminate the nonuniform position of CO. Lower recombination rates in the rec-1 mutant are particularly found in the autosomal arm domains containing the pairing centers. We further find that the rec-1 mutant is of little consequence for organismal fertility and egg viability and thus for rates of autosomal nondisjunction. It nonetheless increases X chromosome nondisjunction rates and thus male appearance. Our findings question the maintenance of recombination rate heritability and genetic diversity among C. elegans natural populations, and they further suggest that manipulating genetic modifiers of CO position will help find quantitative trait loci located in low-recombining genomic regions normally refractory to discovery.

Absence of long-term balancing selection on variation in EuMYB3, an R2R3-MYB gene responsible for the anther-color polymorphism in Erythronium umbilicatum

Balancing selection has been shown to be common in plants for several different types of traits, such as self-incompatibility and heterostyly. Generally, for these traits balancing selection is generated by interactions among individuals or between individuals and other species (e.g., pathogens or pollinators). However, there are phenotypic polymorphisms in plants that do not obviously involve types of interactions that generate balancing selection. Little is known about the extent to which balancing selection also acts to preserve these polymorphisms. Here we ask whether balancing selection preserves an anther-color polymorphism in Erythronium umbilicatum (Liliaceae). We identified a major gene underlying this polymorphism. We then attempted to detect signatures of balancing selection on that gene by developing a new coalescence test for balancing selection. We found that variation in anther color is in large part caused by variation in a paralog of EuMYB3, an anthocyanin-regulating R2R3-MYB transcription factor. However, we found little evidence for balancing selection having acted historically on EuMYB3. Our results thus suggest that plant polymorphisms, especially those not involved in interactions that are likely to generate negative frequency-dependent selection, may reflect a transient state in which one morph will eventually be fixed by either genetic drift or directional selection. Our results also suggest that regulation of the anthocyanin pathway is more evolutionarily labile than is generally believed.

Turnovers of Sex-Determining Mutation in the Golden Pompano and Related Species Provide Insights into Microevolution of Undifferentiated Sex Chromosome

The suppression of recombination is considered a hallmark of sex chromosome evolution. However, previous research has identified undifferentiated sex chromosomes and sex determination by single SNP in the greater amberjack (Seriola dumerili). We observed the same phenomena in the golden pompano (Trachinotus ovatus) of the same family Carangidae and discovered a different sex-determining SNP within the same gene Hsd17b1. We propose an evolutionary model elucidating the turnover of sex-determining mutations by highlighting the contrasting dynamics between purifying selection, responsible for maintaining W-linked Hsd17b1, and neutral evolution, which drives Z-linked Hsd17b1. Additionally, sporadic loss-of-function mutations in W-linked Hsd17b1 contribute to the conversion of W chromosomes into Z chromosomes. This model was directly supported by simulations, closely related species, and indirectly by zebrafish mutants. These findings shed new light on the early stages of sex chromosome evolution.

The genome-wide meiotic recombination landscape in ciliates and its implications for crossover regulation and genome evolution

Meiotic recombination is essential for sexual reproduction and its regulation has been extensively studied in many taxa. However, genome-wide recombination landscape has not been reported in ciliates and it remains unknown how it is affected by the unique features of ciliates: the synaptonemal complex (SC)-independent meiosis and the nuclear dimorphism. Here, we show the recombination landscape in the model ciliate Tetrahymena thermophila by analyzing single-nucleotide polymorphism datasets from 38 hybrid progeny. We detect 1021 crossover (CO) events (35.8 per meiosis), corresponding to an overall CO rate of 9.9 cM/Mb. However, gene conversion by non-crossover is rare (1.03 per meiosis) and not biased towards G or C alleles. Consistent with the reported roles of SC in CO interference, we find no obvious sign of CO interference. CO tends to occur within germ-soma common genomic regions and many of the 44 identified CO hotspots localize at the centromeric or subtelomeric regions. Gene ontology analyses show that CO hotspots are strongly associated with genes responding to environmental changes. We discuss these results with respect to how nuclear dimorphism has potentially driven the formation of the observed recombination landscape to facilitate environmental adaptation and the sharing of machinery among meiotic and somatic recombination.

The Promise of Inferring the Past Using the Ancestral Recombination Graph

The ancestral recombination graph (ARG) is a structure that represents the history of coalescent and recombination events connecting a set of sequences (Hudson RR. In: Futuyma D, Antonovics J, editors. Gene genealogies and the coalescent process. In: Oxford Surveys in Evolutionary Biology; 1991. p. 1 to 44.). The full ARG can be represented as a set of genealogical trees at every locus in the genome, annotated with recombination events that change the topology of the trees between adjacent loci and the mutations that occurred along the branches of those trees (Griffiths RC, Marjoram P. An ancestral recombination graph. In: Donnelly P, Tavare S, editors. Progress in population genetics and human evolution. Springer; 1997. p. 257 to 270.). Valuable insights can be gained into past evolutionary processes, such as demographic events or the influence of natural selection, by studying the ARG. It is regarded as the "holy grail" of population genetics (Hubisz M, Siepel A. Inference of ancestral recombination graphs using ARGweaver. In: Dutheil JY, editors. Statistical population genomics. New York, NY: Springer US; 2020. p. 231-266.) since it encodes the processes that generate all patterns of allelic and haplotypic variation from which all commonly used summary statistics in population genetic research (e.g. heterozygosity and linkage disequilibrium) can be derived. Many previous evolutionary inferences relied on summary statistics extracted from the genotype matrix. Evolutionary inferences using the ARG represent a significant advancement as the ARG is a representation of the evolutionary history of a sample that shows the past history of recombination, coalescence, and mutation events across a particular sequence. This representation in theory contains as much information, if not more, than the combination of all independent summary statistics that could be derived from the genotype matrix. Consistent with this idea, some of the first ARG-based analyses have proven to be more powerful than summary statistic-based analyses (Speidel L, Forest M, Shi S, Myers SR. A method for genome-wide genealogy estimation for thousands of samples. Nat Genet. 2019:51(9):1321 to 1329.; Stern AJ, Wilton PR, Nielsen R. An approximate full-likelihood method for inferring selection and allele frequency trajectories from DNA sequence data. PLoS Genet. 2019:15(9):e1008384.; Hubisz MJ, Williams AL, Siepel A. Mapping gene flow between ancient hominins through demography-aware inference of the ancestral recombination graph. PLoS Genet. 2020:16(8):e1008895.; Fan C, Mancuso N, Chiang CWK. A genealogical estimate of genetic relationships. Am J Hum Genet. 2022:109(5):812-824.; Fan C, Cahoon JL, Dinh BL, Ortega-Del Vecchyo D, Huber C, Edge MD, Mancuso N, Chiang CWK. A likelihood-based framework for demographic inference from genealogical trees. bioRxiv. 2023.10.10.561787. 2023.; Hejase HA, Mo Z, Campagna L, Siepel A. A deep-learning approach for inference of selective sweeps from the ancestral recombination graph. Mol Biol Evol. 2022:39(1):msab332.; Link V, Schraiber JG, Fan C, Dinh B, Mancuso N, Chiang CWK, Edge MD. Tree-based QTL mapping with expected local genetic relatedness matrices. bioRxiv. 2023.04.07.536093. 2023.; Zhang BC, Biddanda A, Gunnarsson ÁF, Cooper F, Palamara PF. Biobank-scale inference of ancestral recombination graphs enables genealogical analysis of complex traits. Nat Genet. 2023:55(5):768-776.). As such, there has been significant interest in the field to investigate 2 main problems related to the ARG: (i) How can we estimate the ARG based on genomic data, and (ii) how can we extract information of past evolutionary processes from the ARG? In this perspective, we highlight 3 topics that pertain to these main issues: The development of computational innovations that enable the estimation of the ARG; remaining challenges in estimating the ARG; and methodological advances for deducing evolutionary forces and mechanisms using the ARG. This perspective serves to introduce the readers to the types of questions that can be explored using the ARG and to highlight some of the most pressing issues that must be addressed in order to make ARG-based inference an indispensable tool for evolutionary research.

Anthropogenic disturbance driving population decline of a dominant tree in East Asia evergreen broadleaved forests over the last 11,000 years

Current biodiversity loss is generally considered to have been caused by anthropogenic disturbance, but it is unclear when anthropogenic activities began to affect biodiversity loss. One hypothesis suggests it began with the Industrial Revolution, whereas others propose that anthropogenic disturbance has been associated with biodiversity decline since the early Holocene. To test these hypotheses, we examined the unique vegetation of evergreen broadleaved forests (EBLFs) in East Asia, where humans have affected landscapes since the early Holocene. We adopted a genomic approach to infer the demographic history of a dominant tree (Litsea elongata) of EBLFs. We used Holocene temperature and anthropogenic disturbance factors to calculate the correlation between these variables and the historical effective population size of L. elongata with Spearman statistics and integrated the maximum-entropy niche model to determine the impact of climate change and anthropogenic disturbance on fluctuation in its effective population size. We identified 9 well-defined geographic clades for the populations of L. elongata. Based on the estimated historical population sizes of these clades, all the populations contracted, indicating persistent population decline over the last 11,000 years. Demographic history of L. elongata and human population change, change in cropland use, and change in irrigated rice area were significantly negatively correlated, whereas climate change in the Holocene was not correlated with demographic history. Our results support the early human impact hypothesis and provide comprehensive evidence that early anthropogenic disturbance may contribute to the current biodiversity crisis in East Asia.

Divergence and gene flow history at two large chromosomal inversions underlying ecotype differentiation in the long-snouted seahorse

Chromosomal inversions can play an important role in divergence and reproductive isolation by building and maintaining distinct allelic combinations between evolutionary lineages. Alternatively, they can take the form of balanced polymorphisms that segregate within populations until one arrangement becomes fixed. Many questions remain about how inversion polymorphisms arise, how they are maintained over the long term, and ultimately, whether and how they contribute to speciation. The long-snouted seahorse (Hippocampus guttulatus) is genetically subdivided into geographic lineages and marine-lagoon ecotypes, with shared structural variation underlying lineage and ecotype divergence. Here, we aim to characterize structural variants and to reconstruct their history and suspected role in ecotype formation. We generated a near chromosome-level genome assembly and described genome-wide patterns of diversity and divergence through the analysis of 112 whole-genome sequences from Atlantic, Mediterranean, and Black Sea populations. By also analysing linked-read sequencing data, we found evidence for two chromosomal inversions that were several megabases in length and showed contrasting allele frequency patterns between lineages and ecotypes across the species range. We reveal that these inversions represent ancient intraspecific polymorphisms, one likely being maintained by divergent selection and the other by pseudo-overdominance. A possible selective coupling between the two inversions was further supported by the absence of specific haplotype combinations and a putative functional interaction between the two inversions in reproduction. Lastly, we detected gene flux eroding divergence between inverted alleles at varying levels for the two inversions, with a likely impact on their dynamics and contribution to divergence and speciation.

Multiple reorganizations of the lateral elements of the synaptonemal complex facilitate homolog segregation in Bombyx mori oocytes

Although Lepidopteran females build a synaptonemal complex (SC) in pachytene, homologs do not crossover, necessitating an alternative method of homolog conjunction. In Bombyx mori oocytes, the SC breaks down at the end of pachytene, and homolog associations are maintained by a large oocyte-specific structure, which we call the bivalent bridge (BB), connecting paired homologs. The BB is derived from at least some components of the SC lateral elements (LEs). It contains the HORMAD protein HOP1 and the LE protein SYCP2 and is formed by the fusion of the two LE derivatives. As diplotene progresses, the BB increases in width and acquires a layered structure with a thick band of HOP1 separating two layers of SYCP2. The HOP1 interacting protein, PCH2, joins the BB in mid-diplotene, and by late-diplotene, it lies in the middle of the HOP1 filament. This structure is maintained through metaphase I. SYCP2 and PCH2 are lost at anaphase I, and the BB no longer connects the separating homologs. However, a key component of the BB, HOP1, remains at the metaphase I plate. These changes in organization of the BB occur simultaneously with the movement of the kinetochore protein, DSN1, from within the BB at mid-diplotene to the edge of the homologs facing the poles by metaphase I. We view these data in context of models in which SC components and regulators can be repurposed to achieve different functions, a fascinating example of evolution achieving homolog conjunction in an alternative way with recycling of SC proteins.

Novel Insights into the Landscape of Crossover and Noncrossover Events in Rhesus Macaques (Macaca mulatta)

Meiotic recombination landscapes differ greatly between distantly and closely related taxa, populations, individuals, sexes, and even within genomes; however, the factors driving this variation are yet to be well elucidated. Here, we directly estimate contemporary crossover rates and, for the first time, noncrossover rates in rhesus macaques (Macaca mulatta) from four three-generation pedigrees comprising 32 individuals. We further compare these results with historical, demography-aware, linkage disequilibrium-based recombination rate estimates. From paternal meioses in the pedigrees, 165 crossover events with a median resolution of 22.3 kb were observed, corresponding to a male autosomal map length of 2,357 cM-approximately 15% longer than an existing linkage map based on human microsatellite loci. In addition, 85 noncrossover events with a mean tract length of 155 bp were identified-similar to the tract lengths observed in the only other two primates in which noncrossovers have been studied to date, humans and baboons. Consistent with observations in other placental mammals with PRDM9-directed recombination, crossover (and to a lesser extent noncrossover) events in rhesus macaques clustered in intergenic regions and toward the chromosomal ends in males-a pattern in broad agreement with the historical, sex-averaged recombination rate estimates-and evidence of GC-biased gene conversion was observed at noncrossover sites.

Paramecium Genetics, Genomics, and Evolution

The ciliate genus Paramecium served as one of the first model systems in microbial eukaryotic genetics, contributing much to the early understanding of phenomena as diverse as genome rearrangement, cryptic speciation, cytoplasmic inheritance, and endosymbiosis, as well as more recently to the evolution of mating types, introns, and roles of small RNAs in DNA processing. Substantial progress has recently been made in the area of comparative and population genomics. Paramecium species combine some of the lowest known mutation rates with some of the largest known effective populations, along with likely very high recombination rates, thereby harboring a population-genetic environment that promotes an exceptionally efficient capacity for selection. As a consequence, the genomes are extraordinarily streamlined, with very small intergenic regions combined with small numbers of tiny introns. The subject of the bulk of Paramecium research, the ancient Paramecium aurelia species complex, is descended from two whole-genome duplication events that retain high degrees of synteny, thereby providing an exceptional platform for studying the fates of duplicate genes. Despite having a common ancestor dating to several hundred million years ago, the known descendant species are morphologically indistinguishable, raising significant questions about the common view that gene duplications lead to the origins of evolutionary novelties.

RecView: an interactive R application for locating recombination positions using pedigree data

BACKGROUND

Recombination reshuffles alleles at linked loci, allowing genes to evolve independently and consequently enhancing the efficiency of selection. This makes quantifying recombination along chromosomes an important goal for understanding how selection and drift are acting on genes and chromosomes.

RESULTS

We present RecView, an interactive R application and its homonymous R package, to facilitate locating recombination positions along chromosomes or scaffolds using whole-genome genotype data of a three-generation pedigree. RecView analyses and plots the grandparent-of-origin of all informative alleles along each chromosome of the offspring in the pedigree, and infers recombination positions with either of two built-in algorithms: one based on change in the proportion of the alleles with specific grandparent-of-origin, and one on the degree of continuity of alleles with the same grandparent-of-origin. RecView handles multiple offspring and chromosomes simultaneously, and all putative recombination positions are reported in base pairs together with an estimated precision based on the local density of informative alleles. We demonstrate RecView using genotype data of a passerine bird with an available reference genome, the great reed warbler (Acrocephalus arundinaceus), and show that recombination events can be located to specific positions.

CONCLUSIONS

RecView is an easy-to-use and highly effective application for locating recombination positions with high precision. RecView is available on GitHub ( https://github.com/HKyleZhang/RecView.git ).

Genetic architecture of individual meiotic crossover rate and distribution in Atlantic Salmon

Meiotic recombination through chromosomal crossovers ensures proper segregation of homologous chromosomes during meiosis, while also breaking down linkage disequilibrium and shuffling alleles at loci located on the same chromosome. Rates of recombination can vary between species, but also between and within individuals, sex and chromosomes within species. Indeed, the Atlantic salmon genome is known to have clear sex differences in recombination with female biased heterochiasmy and markedly different landscapes of crossovers between males and females. In male meiosis, crossovers occur strictly in the telomeric regions, whereas in female meiosis crossovers tend to occur closer to the centromeres. However, little is known about the genetic control of these patterns and how this differs at the individual level. Here, we investigate genetic variation in individual measures of recombination in > 5000 large full-sib families of a Norwegian Atlantic salmon breeding population with high-density SNP genotypes. We show that females had 1.6 × higher crossover counts (CC) than males, with autosomal linkage maps spanning a total of 2174 cM in females and 1483 cM in males. However, because of the extreme telomeric bias of male crossovers, female recombination is much more important for generation of new haplotypes with 8 × higher intra-chromosomal genetic shuffling than males. CC was heritable in females (h2 = 0.11) and males (h2 = 0.10), and shuffling was also heritable in both sex but with a lower heritability in females (h2 = 0.06) than in males (h2 = 0.11). Inter-sex genetic correlations for both traits were close to zero, suggesting that rates and distribution of crossovers are genetically distinct traits in males and females, and that there is a potential for independent genetic change in both sexes in the Atlantic Salmon. Together, these findings give novel insights into the genetic architecture of recombination in salmonids and contribute to a better understanding of how rates and distribution of recombination may evolve in eukaryotes more broadly.

The effects of mutation and recombination rate heterogeneity on the inference of demography and the distribution of fitness effects

Disentangling the effects of demography and selection has remained a focal point of population genetic analysis. Knowledge about mutation and recombination is essential in this endeavour; however, despite clear evidence that both mutation and recombination rates vary across genomes, it is common practice to model both rates as fixed. In this study, we quantify how this unaccounted for rate heterogeneity may impact inference using common approaches for inferring selection (DFE-alpha, Grapes, and polyDFE) and/or demography (fastsimcoal2 and δaδi). We demonstrate that, if not properly modelled, this heterogeneity can increase uncertainty in the estimation of demographic and selective parameters and in some scenarios may result in mis-leading inference. These results highlight the importance of quantifying the fundamental evolutionary parameters of mutation and recombination prior to utilizing population genomic data to quantify the effects of genetic drift (i.e., as modulated by demographic history) and selection; or, at the least, that the effects of uncertainty in these parameters can and should be directly modelled in downstream inference.

The chromatin determinants and Ph1 gene effect at wheat sites with contrasting recombination frequency

INTRODUCTION

Meiotic recombination is one of the most important processes of evolution and adaptation to environmental conditions. Even though there is substantial knowledge about proteins involved in the process, targeting specific DNA loci by the recombination machinery is not well understood.

OBJECTIVES

This study aims to investigate a wheat recombination hotspot (H1) in comparison with a "regular" recombination site (Rec7) on the sequence and epigenetic level in conditions with functional and non-functional Ph1 locus.

METHODS

The DNA sequence, methylation pattern, and recombination frequency were analyzed for the H1 and Rec7 in three mapping populations derived by crossing introgressive wheat line 8.1 with cv. Chinese Spring (with Ph1 and ph1 alleles) and cv. Tähti.

RESULTS

The H1 and Rec7 loci are 1.586 kb and 2.538 kb long, respectively. High-density mapping allowed to delimit the Rec7 and H1 to 19 and 574 bp and 593 and 571 bp CO sites, respectively. A new method (ddPing) allowed screening recombination frequency in almost 66 thousand gametes. The screening revealed a 5.94-fold higher recombination frequency at the H1 compared to the Rec7. The H1 was also found out of the Ph1 control, similarly as gamete distortion. The recombination was strongly affected by larger genomic rearrangements but not by the SNP proximity. Moreover, chromatin markers for open chromatin and DNA hypomethylation were found associated with crossover occurrence except for the CHH methylation.

CONCLUSION

Our results, for the first time, allowed study of wheat recombination directly on sequence, shed new light on chromatin landmarks associated with particular recombination sites, and deepened knowledge about role of the Ph1 locus in control of wheat recombination processes. The results are suggesting more than one recombination control pathway. Understanding this phenomenon may become a base for more efficient wheat genome manipulation, gene pool enrichment, breeding, and study processes of recombination itself.

Evaluating power to detect recurrent selective sweeps under increasingly realistic evolutionary null models

The detection of selective sweeps from population genomic data often relies on the premise that the beneficial mutations in question have fixed very near the sampling time. As it has been previously shown that the power to detect a selective sweep is strongly dependent on the time since fixation as well as the strength of selection, it is naturally the case that strong, recent sweeps leave the strongest signatures. However, the biological reality is that beneficial mutations enter populations at a rate, one that partially determines the mean wait time between sweep events and hence their age distribution. An important question thus remains about the power to detect recurrent selective sweeps when they are modeled by a realistic mutation rate and as part of a realistic distribution of fitness effects, as opposed to a single, recent, isolated event on a purely neutral background as is more commonly modeled. Here we use forward-in-time simulations to study the performance of commonly used sweep statistics, within the context of more realistic evolutionary baseline models incorporating purifying and background selection, population size change, and mutation and recombination rate heterogeneity. Results demonstrate the important interplay of these processes, necessitating caution when interpreting selection scans; specifically, false-positive rates are in excess of true-positive across much of the evaluated parameter space, and selective sweeps are often undetectable unless the strength of selection is exceptionally strong.

Measuring the Impact of Genetic Heterogeneity and Chromosomal Inversions on the Efficacy of CRISPR-Cas9 Gene Drives in Different Strains of Anopheles gambiae

The human malaria vector Anopheles gambiae is becoming increasingly resistant to insecticides, spurring the development of genetic control strategies. CRISPR-Cas9 gene drives can modify a population by creating double-stranded breaks at highly specific targets, triggering copying of the gene drive into the cut site ("homing"), ensuring its inheritance. The DNA repair mechanism responsible requires homology between the donor and recipient chromosomes, presenting challenges for the invasion of laboratory-developed gene drives into wild populations of target species An. gambiae species complex, which show high levels of genome variation. Two gene drives (vas2-5958 and zpg-7280) were introduced into three An. gambiae strains collected across Africa with 5.3-6.6% variation around the target sites, and the effect of this variation on homing was measured. Gene drive homing across different karyotypes of the 2La chromosomal inversion was also assessed. No decrease in gene drive homing was seen despite target site heterology, demonstrating the applicability of gene drives to wild populations.

Meiotic recombination is confirmed to be unusually high in the fission yeast Schizosaccharomyces pombe

In most eukaryotes, meiotic crossovers (COs) are limited to 1-3 per chromosome, and are prevented from occurring close to one another by CO interference. The fission yeast Schizosaccharomyces pombe, an exception to these general rules, was reported to have the highest CO number per chromosome and no or weak interference. However, global CO frequency was indirectly estimated, calling for confirmation. Here, we used an innovative strategy to determine COs genome-wide in S. pombe. We confirmed weak CO interference, acting at physical distances compatible with the patterning of recombination precursors. We revealed a slight co-variation in CO number between chromosomes, suggesting that a limiting pro-CO factor varies between meiocytes. CO number per chromosome varies proportionally with chromosome size, with the three chromosomes having, on average, 15.9, 12.5, and 7.0 COs, respectively. This reinforces S. pombe's status as the eukaryote with the highest CO number per chromosome described to date.

FREEDA: An automated computational pipeline guides experimental testing of protein innovation

Cell biologists typically focus on conserved regions of a protein, overlooking innovations that can shape its function over evolutionary time. Computational analyses can reveal potential innovations by detecting statistical signatures of positive selection that lead to rapid accumulation of beneficial mutations. However, these approaches are not easily accessible to non-specialists, limiting their use in cell biology. Here, we present an automated computational pipeline FREEDA that provides a simple graphical user interface requiring only a gene name; integrates widely used molecular evolution tools to detect positive selection in rodents, primates, carnivores, birds, and flies; and maps results onto protein structures predicted by AlphaFold. Applying FREEDA to >100 centromere proteins, we find statistical evidence of positive selection within loops and turns of ancient domains, suggesting innovation of essential functions. As a proof-of-principle experiment, we show innovation in centromere binding of mouse CENP-O. Overall, we provide an accessible computational tool to guide cell biology research and apply it to experimentally demonstrate functional innovation.

The human fungal pathogen Aspergillus fumigatus can produce the highest known number of meiotic crossovers

Sexual reproduction involving meiosis is essential in most eukaryotes. This produces offspring with novel genotypes, both by segregation of parental chromosomes as well as crossovers between homologous chromosomes. A sexual cycle for the opportunistic human pathogenic fungus Aspergillus fumigatus is known, but the genetic consequences of meiosis have remained unknown. Among other Aspergilli, it is known that A. flavus has a moderately high recombination rate with an average of 4.2 crossovers per chromosome pair, whereas A. nidulans has in contrast a higher rate with 9.3 crossovers per chromosome pair. Here, we show in a cross between A. fumigatus strains that they produce an average of 29.9 crossovers per chromosome pair and large variation in total map length across additional strain crosses. This rate of crossovers per chromosome is more than twice that seen for any known organism, which we discuss in relation to other genetic model systems. We validate this high rate of crossovers through mapping of resistance to the laboratory antifungal acriflavine by using standing variation in an undescribed ABC efflux transporter. We then demonstrate that this rate of crossovers is sufficient to produce one of the common multidrug resistant haplotypes found in the cyp51A gene (TR34/L98H) in crosses among parents harboring either of 2 nearby genetic variants, possibly explaining the early spread of such haplotypes. Our results suggest that genomic studies in this species should reassess common assumptions about linkage between genetic regions. The finding of an unparalleled crossover rate in A. fumigatus provides opportunities to understand why these rates are not generally higher in other eukaryotes.

The molecular mechanisms of recombinant chromosome 18 with parental pericentric inversions and a review of the literature

Chromosomal rearrangements mostly result from non-allelic homologous recombination mediated by low-copy repeats (LCRs) or segmental duplications (SDs). Recent studies on recombinant chromosome 18 (rec (18)) have focused on diagnoses and clinical phenotypes. We diagnosed two cases of prenatal rec (18) and identified precise breakpoint intervals using karyotype and chromosomal microarray analyses. We analyzed the distribution characteristics of breakpoint repetitive elements to infer rearrangement mechanisms and reviewed relevant literature to identify genetic trends. Among the 12 families with 25 pregnancies analyzed, 68% rec (18), 24% spontaneous abortions, and 8% normal births were reported. In the 17 rec (18) cases, 65% presented maternal origin and 35% were paternal. Short-arm breakpoints at p11.31 were reported in 10 cases, whereas the long-arm breakpoints were located at q21.3 (6 cases) and q12 (4 cases). Breakpoints of pericentric inversions on chromosome 18 are concentrated in p11.31, q21.3, and q12 regions. Rearrangements at 18p11.31 are non-recurrent events. ALUs, LINE1s, and MIRs were enriched at the breakpoint regions (1.85 to 3.42-fold enrichment over the entire chromosome 18), while SDs and LCRs were absent. ALU subfamilies had sequence identities of 85.94% and 83.01% between two pair breakpoints. Small repetitive elements may mediate recombination-coupled DNA repair processes, facilitating rearrangements on chromosome 18. Maternal inversion carriers are more prone to abnormal recombination in prenatal families with rec (18). Recombinant chromosomes may present preferential segregation during gamete formation.

Nascent evolution of recombination rate differences as a consequence of chromosomal rearrangements

Reshuffling of genetic variation occurs both by independent assortment of chromosomes and by homologous recombination. Such reshuffling can generate novel allele combinations and break linkage between advantageous and deleterious variants which increases both the potential and the efficacy of natural selection. Here we used high-density linkage maps to characterize global and regional recombination rate variation in two populations of the wood white butterfly (Leptidea sinapis) that differ considerably in their karyotype as a consequence of at least 27 chromosome fissions and fusions. The recombination data were compared to estimates of genetic diversity and measures of selection to assess the relationship between chromosomal rearrangements, crossing over, maintenance of genetic diversity and adaptation. Our data show that the recombination rate is influenced by both chromosome size and number, but that the difference in the number of crossovers between karyotypes is reduced as a consequence of a higher frequency of double crossovers in larger chromosomes. As expected from effects of selection on linked sites, we observed an overall positive association between recombination rate and genetic diversity in both populations. Our results also revealed a significant effect of chromosomal rearrangements on the rate of intergenic diversity change between populations, but limited effects on polymorphisms in coding sequence. We conclude that chromosomal rearrangements can have considerable effects on the recombination landscape and consequently influence both maintenance of genetic diversity and efficiency of selection in natural populations.

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.

Evaluating power to detect recurrent selective sweeps under increasingly realistic evolutionary null models

The detection of selective sweeps from population genomic data often relies on the premise that the beneficial mutations in question have fixed very near the sampling time. As it has been previously shown that the power to detect a selective sweep is strongly dependent on the time since fixation as well as the strength of selection, it is naturally the case that strong, recent sweeps leave the strongest signatures. However, the biological reality is that beneficial mutations enter populations at a rate, one that partially determines the mean wait time between sweep events and hence their age distribution. An important question thus remains about the power to detect recurrent selective sweeps when they are modelled by a realistic mutation rate and as part of a realistic distribution of fitness effects (DFE), as opposed to a single, recent, isolated event on a purely neutral background as is more commonly modelled. Here we use forward-in-time simulations to study the performance of commonly used sweep statistics, within the context of more realistic evolutionary baseline models incorporating purifying and background selection, population size change, and mutation and recombination rate heterogeneity. Results demonstrate the important interplay of these processes, necessitating caution when interpreting selection scans; specifically, false positive rates are in excess of true positive across much of the evaluated parameter space, and selective sweeps are often undetectable unless the strength of selection is exceptionally strong.

Teaser Text

Outlier-based genomic scans have proven a popular approach for identifying loci that have potentially experienced recent positive selection. However, it has previously been shown that an evolutionarily appropriate baseline model that incorporates non-equilibrium population histories, purifying and background selection, and variation in mutation and recombination rates is necessary to reduce often extreme false positive rates when performing genomic scans. Here we evaluate the power to detect recurrent selective sweeps using common SFS-based and haplotype-based methods under these increasingly realistic models. We find that while these appropriate evolutionary baselines are essential to reduce false positive rates, the power to accurately detect recurrent selective sweeps is generally low across much of the biologically relevant parameter space.

Patterns, mechanisms, and consequences of homoeologous exchange in allopolyploid angiosperms: a genomic and epigenomic perspective

Allopolyploids result from hybridization between different evolutionary lineages coupled with genome doubling. Homoeologous chromosomes (chromosomes with common shared ancestry) may undergo recombination immediately after allopolyploid formation and continue over successive generations. The outcome of this meiotic pairing behavior is dynamic and complex. Homoeologous exchanges (HEs) may lead to the formation of unbalanced gametes, reduced fertility, and selective disadvantage. By contrast, HEs could act as sources of novel evolutionary substrates, shifting the relative dosage of parental gene copies, generating novel phenotypic diversity, and helping the establishment of neo-allopolyploids. However, HE patterns vary among lineages, across generations, and even within individual genomes and chromosomes. The causes and consequences of this variation are not fully understood, though interest in this evolutionary phenomenon has increased in the last decade. Recent technological advances show promise in uncovering the mechanistic basis of HEs. Here, we describe recent observations of the common patterns among allopolyploid angiosperm lineages, underlying genomic and epigenomic features, and consequences of HEs. We identify critical research gaps and discuss future directions with far-reaching implications in understanding allopolyploid evolution and applying them to the development of important phenotypic traits of polyploid crops.

Divergent patterns of meiotic double strand breaks and synapsis initiation dynamics suggest an evolutionary shift in the meiosis program between American and Australian marsupials

In eutherian mammals, hundreds of programmed DNA double-strand breaks (DSBs) are generated at the onset of meiosis. The DNA damage response is then triggered. Although the dynamics of this response is well studied in eutherian mammals, recent findings have revealed different patterns of DNA damage signaling and repair in marsupial mammals. To better characterize these differences, here we analyzed synapsis and the chromosomal distribution of meiotic DSBs markers in three different marsupial species (Thylamys elegans, Dromiciops gliorides, and Macropus eugenii) that represent South American and Australian Orders. Our results revealed inter-specific differences in the chromosomal distribution of DNA damage and repair proteins, which were associated with differing synapsis patterns. In the American species T. elegans and D. gliroides, chromosomal ends were conspicuously polarized in a bouquet configuration and synapsis progressed exclusively from the telomeres towards interstitial regions. This was accompanied by sparse H2AX phosphorylation, mainly accumulating at chromosomal ends. Accordingly, RAD51 and RPA were mainly localized at chromosomal ends throughout prophase I in both American marsupials, likely resulting in reduced recombination rates at interstitial positions. In sharp contrast, synapsis initiated at both interstitial and distal chromosomal regions in the Australian representative M. eugenii, the bouquet polarization was incomplete and ephemeral, γH2AX had a broad nuclear distribution, and RAD51 and RPA foci displayed an even chromosomal distribution. Given the basal evolutionary position of T. elegans, it is likely that the meiotic features reported in this species represent an ancestral pattern in marsupials and that a shift in the meiotic program occurred after the split of D. gliroides and the Australian marsupial clade. Our results open intriguing questions about the regulation and homeostasis of meiotic DSBs in marsupials. The low recombination rates observed at the interstitial chromosomal regions in American marsupials can result in the formation of large linkage groups, thus having an impact in the evolution of their genomes.

Structure-Based GC Investigation Sheds New Light on ITS2 Evolution in Corydalis Species

Guanine and cytosine (GC) content is a fundamental component of genetic diversity and essential for phylogenetic analyses. However, the GC content of the ribosomal internal transcribed spacer 2 (ITS2) remains unknown, despite the fact that ITS2 is a widely used phylogenetic marker. Here, the ITS2 was high-throughput sequenced from 29 Corydalis species, and their GC contents were comparatively investigated in the context of ITS2's characteristic secondary structure and concerted evolution. Our results showed that the GC contents of ITS2 were 131% higher than those of their adjacent 5.8S regions, suggesting that ITS2 underwent GC-biased evolution. These GCs were distributed in a heterogeneous manner in the ITS2 secondary structure, with the paired regions being 130% larger than the unpaired regions, indicating that GC is chosen for thermodynamic stability. In addition, species with homogeneous ITS2 sequences were always GC-rich, supporting GC-biased gene conversion (gBGC), which occurred with ITS2's concerted evolution. The RNA substitution model inferred also showed a GC preference among base pair transformations, which again supports gBGC. Overall, structurally based GC investigation reveals that ITS2 evolves under structural stability and gBGC selection, significantly increasing its GC content.

The coordinated regulation of early meiotic stages is dominated by non-coding RNAs and stage-specific transcription in wheat

Reproductive success hinges on precisely coordinated meiosis, yet our understanding of how structural rearrangements of chromatin and phase transitions during meiosis are transcriptionally regulated is limited. In crop plants, detailed analysis of the meiotic transcriptome could identify regulatory genes and epigenetic regulators that can be targeted to increase recombination rates and broaden genetic variation, as well as provide a resource for comparison among eukaryotes of different taxa to answer outstanding questions about meiosis. We conducted a meiotic stage-specific analysis of messenger RNA (mRNA), small non-coding RNA (sncRNA), and long intervening/intergenic non-coding RNA (lincRNA) in wheat (Triticum aestivum L.) and revealed novel mechanisms of meiotic transcriptional regulation and meiosis-specific transcripts. Amidst general repression of mRNA expression, significant enrichment of ncRNAs was identified during prophase I relative to vegetative cells. The core meiotic transcriptome was comprised of 9309 meiosis-specific transcripts, 48 134 previously unannotated meiotic transcripts, and many known and novel ncRNAs differentially expressed at specific stages. The abundant meiotic sncRNAs controlled the reprogramming of central metabolic pathways by targeting genes involved in photosynthesis, glycolysis, hormone biosynthesis, and cellular homeostasis, and lincRNAs enhanced the expression of nearby genes. Alternative splicing was not evident in this polyploid species, but isoforms were switched at phase transitions. The novel, stage-specific regulatory controls uncovered here challenge the conventional understanding of this crucial biological process and provide a new resource of requisite knowledge for those aiming to directly modulate meiosis to improve crop plants. The wheat meiosis transcriptome dataset can be queried for genes of interest using an eFP browser located at https://bar.utoronto.ca/efp_wheat/cgi-bin/efpWeb.cgi?dataSource=Wheat_Meiosis.

Genomic regions underlying the species-specific mating songs of green lacewings

Rapid species radiations provide insight into the process of speciation and diversification. The radiation of Chrysoperla carnea-group lacewings seems to be driven, at least in part, by their species-specific pre-mating vibrational duets. We associated genetic markers from across the genome with courtship song period in the offspring of a laboratory cross between Chrysoperla plorabunda and Chrysoperla adamsi, two species primarily differentiated by their mating songs. Two genomic regions were strongly associated with the song period phenotype. Large regions of chromosomes one and two were associated with song phenotype, as fewer recombination events occurred on these chromosomes relative to the other autosomes. Candidate genes were identified by functional annotation of proteins from the C. carnea reference genome. The majority of genes that are associated with vibrational courtship signals in other insects were found within QTL for lacewing song phenotype. Together these findings suggest that decreased recombination may be acting to keep together loci important to reproductive isolation between these species. Using wild-caught individuals from both species, we identified signals of genomic divergence across the genome. We identified several candidate genes both in song-associated regions and near divergence outliers including nonA, fruitless, paralytic, period, and doublesex. Together these findings bring us one step closer to identifying the genomic basis of a mating song trait critical to the maintenance of species boundaries in green lacewings.

Variation in mutation, recombination, and transposition rates in Drosophila melanogaster and Drosophila simulans

The rates of mutation, recombination, and transposition are core parameters in models of evolution. They impact genetic diversity, responses to ongoing selection, and levels of genetic load. However, even for key evolutionary model species such as Drosophila melanogaster and Drosophila simulans, few estimates of these parameters are available, and we have little idea of how rates vary between individuals, sexes, or populations. Knowledge of this variation is fundamental for parameterizing models of genome evolution. Here, we provide direct estimates of mutation, recombination, and transposition rates and their variation in a West African and a European population of D. melanogaster and a European population of D. simulans Across 89 flies, we observe 58 single-nucleotide mutations, 286 crossovers, and 89 transposable element (TE) insertions. Compared to the European D. melanogaster, we find the West African population has a lower mutation rate (1.67 × 10-9 site-1 gen-1 vs. 4.86 × 10-9 site-1 gen-1) and a lower transposition rate (8.99 × 10-5 copy-1 gen-1 vs. 23.36 × 10-5 copy-1 gen-1), but a higher recombination rate (3.44 cM/Mb vs. 2.06 cM/Mb). The European D. simulans population has a similar mutation rate to European D. melanogaster, but a significantly higher recombination rate and a lower, but not significantly different, transposition rate. Overall, we find paternal-derived mutations are more frequent than maternal ones in both species. Our study quantifies the variation in rates of mutation, recombination, and transposition among different populations and sexes, and our direct estimates of these parameters in D. melanogaster and D. simulans will benefit future studies in population and evolutionary genetics.

A genetic linkage map and improved genome assembly of the termite symbiont Termitomyces cryptogamus

BACKGROUND

The termite-fungus symbiosis is an ancient stable mutualism of two partners that reproduce and disperse independently. With the founding of each termite colony the symbiotic association must be re-established with a new fungus partner. Complementarity in the ability to break down plant substrate may help to stabilize this symbiosis despite horizontal symbiont transmission. An alternative, non-exclusive, hypothesis is that a reduced rate of evolution may contribute to stabilize the symbiosis, the so-called Red King Effect.

METHODS

To explore this concept, we produced the first linkage map of a species of Termitomyces, using genotyping by sequencing (GBS) of 88 homokaryotic offspring. We constructed a highly contiguous genome assembly using PacBio data and a de-novo evidence-based annotation. This improved genome assembly and linkage map allowed for examination of the recombination landscape and its potential effect on the mutualistic lifestyle.

RESULTS

Our linkage map resulted in a genome-wide recombination rate of 22 cM/Mb, lower than that of other related fungi. However, the total map length of 1370 cM was similar to that of other related fungi.

CONCLUSIONS

The apparently decreased rate of recombination is primarily due to genome expansion of islands of gene-poor repetitive sequences. This study highlights the importance of inclusion of genomic context in cross-species comparisons of recombination rate.

Predicting recombination frequency from map distance

Map distance is one of the key measures in genetics and indicates the expected number of crossovers between two loci. Map distance is estimated from the observed recombination frequency using mapping functions, the most widely used of those, Haldane and Kosambi, being developed at the time when the number of markers was low and unobserved crossovers had a substantial effect on the recombination fractions. In contemporary high-density marker data, the probability of multiple crossovers between adjacent loci is negligible and different mapping functions yield the same result, that is, the recombination frequency between adjacent loci is equal to the map distance in Morgans. However, high-density linkage maps contain an interpretation problem: the map distance over a long interval is additive and its association with recombination frequency is not defined. Here, we demonstrate with high-density linkage maps from humans and stickleback fishes that the inverses of Haldane's and Kosambi's mapping functions systematically underpredict recombination frequencies from map distance. To remedy this, we formulate a piecewise function that yields more accurate predictions of recombination frequency from map distance. Our results demonstrate that the association between map distance and recombination frequency is context-dependent and without a universal solution.

The roles of recombination and selection in shaping genomic divergence in an incipient ecological species complex

Speciation genomic studies have revealed that genomes of diverging lineages are shaped jointly by the actions of gene flow and selection. These evolutionary forces acting in concert with processes such as recombination and genome features such as gene density shape a mosaic landscape of divergence. We investigated the roles of recombination and gene density in shaping the patterns of differentiation and divergence between the cyclically parthenogenetic ecological sister-taxa, Daphnia pulicaria and Daphnia pulex. First, we assembled a phased chromosome-scale genome assembly using trio-binning for D. pulicaria and constructed a genetic map using an F2-intercross panel to understand sex-specific recombination rate heterogeneity. Finally, we used a ddRADseq data set with broad geographic sampling of D. pulicaria, D. pulex, and their hybrids to understand the patterns of genome-scale divergence and demographic parameters. Our study provides the first sex-specific estimates of recombination rates for a cyclical parthenogen, and unlike other eukaryotic species, we observed male-biased heterochiasmy in D. pulicaria, which may be related to this somewhat unique breeding mode. Additionally, regions of high gene density and recombination are generally more divergent than regions of suppressed recombination. Outlier analysis indicated that divergent genomic regions are probably driven by selection on D. pulicaria, the derived lineage colonizing a novel lake habitat. Together, our study supports a scenario of selection acting on genes related to local adaptation shaping genome-wide patterns of differentiation despite high local recombination rates in this species complex. Finally, we discuss the limitations of our data in light of demographic uncertainty.

FREEDA: an automated computational pipeline guides experimental testing of protein innovation by detecting positive selection

Cell biologists typically focus on conserved regions of a protein, overlooking innovations that can shape its function over evolutionary time. Computational analyses can reveal potential innovations by detecting statistical signatures of positive selection that leads to rapid accumulation of beneficial mutations. However, these approaches are not easily accessible to non-specialists, limiting their use in cell biology. Here, we present an automated computational pipeline FREEDA (Finder of Rapidly Evolving Exons in De novo Assemblies) that provides a simple graphical user interface requiring only a gene name, integrates widely used molecular evolution tools to detect positive selection, and maps results onto protein structures predicted by AlphaFold. Applying FREEDA to >100 mouse centromere proteins, we find evidence of positive selection in intrinsically disordered regions of ancient domains, suggesting innovation of essential functions. As a proof-of-principle experiment, we show innovation in centromere binding of CENP-O. Overall, we provide an accessible computational tool to guide cell biology research and apply it to experimentally demonstrate functional innovation.

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.

Stochastic nuclear organization and host-dependent allele contribution in Rhizophagus irregularis

BACKGROUND

Arbuscular mycorrhizal (AM) fungi are arguably the most important symbionts of plants, offering a range of benefits to their hosts. However, the provisioning of these benefits does not appear to be uniform among AM fungal individuals, with genetic variation between fungal symbionts having a substantial impact on plant performance. Interestingly, genetic variation has also been reported within fungal individuals, which contain millions of haploid nuclei sharing a common cytoplasm. In the model AM fungus, Rhizophagus irregularis, several isolates have been reported to be dikaryotes, containing two genetically distinct types of nuclei recognized based on their mating-type (MAT) locus identity. However, their extremely coenocytic nature and lack of a known single nucleus stage has raised questions on the origin, distribution and dynamics of this genetic variation.

RESULTS

Here we performed DNA and RNA sequencing at the mycelial individual, single spore and single nucleus levels to gain insight into the dynamic genetic make-up of the dikaryote-like R. irregularis C3 isolate and the effect of different host plants on its genetic variation. Our analyses reveal that parallel spore and root culture batches can have widely variable ratios of two main genotypes in C3. Additionally, numerous polymorphisms were found with frequencies that deviated significantly from the general genotype ratio, indicating a diverse population of slightly different nucleotypes. Changing host plants did not show consistent host effects on nucleotype ratio's after multiple rounds of subculturing. Instead, we found a major effect of host plant-identity on allele-specific expression in C3.

CONCLUSION

Our analyses indicate a highly dynamic/variable genetic organization in different isolates of R. irregularis. Seemingly random fluctuations in nucleotype ratio's upon spore formation, recombination events, high variability of non-tandemly repeated rDNA sequences and host-dependent allele expression all add levels of variation that may contribute to the evolutionary success of these widespread symbionts.

Lessons from the meiotic recombination landscape of the ZMM deficient budding yeast Lachancea waltii

Meiotic recombination is a driving force for genome evolution, deeply characterized in a few model species, notably in the budding yeast Saccharomyces cerevisiae. Interestingly, Zip2, Zip3, Zip4, Spo16, Msh4, and Msh5, members of the so-called ZMM pathway that implements the interfering meiotic crossover pathway in S. cerevisiae, have been lost in Lachancea yeast species after the divergence of Lachancea kluyveri from the rest of the clade. In this context, after investigating meiosis in L. kluyveri, we determined the meiotic recombination landscape of Lachancea waltii. Attempts to generate diploid strains with fully hybrid genomes invariably resulted in strains with frequent whole-chromosome aneuploidy and multiple extended regions of loss of heterozygosity (LOH), which mechanistic origin is so far unclear. Despite the lack of multiple ZMM pro-crossover factors in L. waltii, numbers of crossovers and noncrossovers per meiosis were higher than in L. kluyveri but lower than in S. cerevisiae, for comparable genome sizes. Similar to L. kluyveri but opposite to S. cerevisiae, L. waltii exhibits an elevated frequency of zero-crossover bivalents. Lengths of gene conversion tracts for both crossovers and non-crossovers in L. waltii were comparable to those observed in S. cerevisiae and shorter than in L. kluyveri despite the lack of Mlh2, a factor limiting conversion tract size in S. cerevisiae. L. waltii recombination hotspots were not shared with either S. cerevisiae or L. kluyveri, showing that meiotic recombination hotspots can evolve at a rather limited evolutionary scale within budding yeasts. Finally, L. waltii crossover interference was reduced relative to S. cerevisiae, with interference being detected only in the 25 kb distance range. Detection of positive inference only at short distance scales in the absence of multiple ZMM factors required for interference-sensitive crossovers in other systems likely reflects interference between early recombination precursors such as DSBs.

The effect of DNA polymorphisms and natural variation on crossover hotspot activity in Arabidopsis hybrids

In hybrid organisms, genetically divergent homologous chromosomes pair and recombine during meiosis; however, the effect of specific types of polymorphisms on crossover is poorly understood. Here, to analyze this in Arabidopsis, we develop the seed-typing method that enables the massively parallel fine-mapping of crossovers by sequencing. We show that structural variants, observed in one of the generated intervals, do not change crossover frequency unless they are located directly within crossover hotspots. Both natural and Cas9-induced deletions result in lower hotspot activity but are not compensated by increases in immediately adjacent hotspots. To examine the effect of single nucleotide polymorphisms on crossover formation, we analyze hotspot activity in mismatch detection-deficient msh2 mutants. Surprisingly, polymorphic hotspots show reduced activity in msh2. In lines where only the hotspot-containing interval is heterozygous, crossover numbers increase above those in the inbred (homozygous). We conclude that MSH2 shapes crossover distribution by stimulating hotspot activity at polymorphic regions.

Chromosome size matters: genome evolution in the cyperid clade

BACKGROUND AND AIMS

While variation in genome size and chromosome numbers and their consequences are often investigated in plants, the biological relevance of variation in chromosome size remains poorly known. Here, we examine genome and mean chromosome size in the cyperid clade (families Cyperaceae, Juncaceae and Thurniaceae), which is the largest vascular plant lineage with predominantly holocentric chromosomes.

METHODS

We measured genome size in 436 species of cyperids using flow cytometry, and augment these data with previously published datasets. We then separately compared genome and mean chromosome sizes (2C/2n) amongst the major lineages of cyperids and analysed how these two genomic traits are associated with various environmental factors using phylogenetically informed methods.

KEY RESULTS

We show that cyperids have the smallest mean chromosome sizes recorded in seed plants, with a large divergence between the smallest and largest values. We found that cyperid species with smaller chromosomes have larger geographical distributions and that there is a strong inverse association between mean chromosome size and number across this lineage.

CONCLUSIONS

The distinct patterns in genome size and mean chromosome size across the cyperids might be explained by holokinetic drive. The numerous small chromosomes might function to increase genetic diversity in this lineage where crossovers are limited during meiosis.

Identification of sex chromosomes and primary sex ratio in the small hive beetle, a worldwide parasite of honey bees

BACKGROUND

The small hive beetle (SHB), Aethina tumida, has emerged as a worldwide threat to honey bees in the past two decades. These beetles harvest nest resources, feed on larval bees, and ultimately spoil nest resources with gelatinous slime together with the fungal symbiont Kodamaea ohmeri.

RESULTS

Here, we present the first chromosome-level genome assembly for the SHB. With a 99.1% representation of conserved (BUSCO) arthropod genes, this resource enables the study of chemosensory, digestive, and detoxification traits critical for SHB success and possible control. We use this annotated assembly to characterize features of SHB sex chromosomes and a female-skewed primary sex ratio. We also found chromosome fusion and a lower recombination rate in sex chromosomes than in autosomes.

CONCLUSIONS

Genome-enabled insights will clarify the traits that allowed this beetle to exploit hive resources successfully and will be critical for determining the causes of observed sex ratio asymmetries.

High-Density Genetic Linkage Map of the Southern Blue-ringed Octopus (Octopodidae: Hapalochlaena maculosa)

Genetic linkage maps provide a useful resource for non-model genomes and can aid in genome reassembly to form more contiguous pseudo-chromosomes. We present the first linkage map of any cephalopod, H. maculosa, composed of 47 linkage groups (LG). A total of 2166 single nucleotide polymorphisms and 2455 presence–absence variant loci were utilised by Lep-Map3 in linkage map construction. The map length spans 2016.62 cM with an average marker distance of 0.85 cM. Integration of the recent H. maculosa genome allowed 1151 scaffolds comprising 34% of the total genomic sequence to be orientated and/or placed using 1278 markers across all 47 LG. The linkage map generated provides a new perspective on HOX gene distribution in octopods. In the H. maculosa linkage map three (SCR, LOX4 and POST1) of six identified HOX genes (HOX1/LAB, SCR, LOX2, LOX4, LOX5, POST1) were located within the same LG (LG 9). The generation of a linkage map for H. maculosa has provided a valuable resource for understanding the evolution of cephalopod genomes and will provide a base for future work.