Search

Nematodes such as Caenorhabditis elegans and Pristionchus pacificus are extremely successful model organisms for comparative biology. Several studies have shown that phenotypic novelty but also conserved processes are controlled by taxon-restricted genes. To trace back the evolution of such new or rapidly evolving genes, a robust phylogenomic framework is indispensable. Here, we present an improved version of the genome of Parapristionchus giblindavisi which is the only known member of the sister group of Pristionchus. Relative to the previous short-read assembly, the new genome is based on long reads and displays higher levels of contiguity, completeness, and correctness. Specifically, the number of contigs dropped from over 7,303 to 735 resulting in an N50 increase from 112 to 791 kb. We made use of the new genome to revisit the evolution of multiple gene families. This revealed Pristionchus-specific expansions of several environmentally responsive gene families and a Pristionchus-specific loss of the de novo purine biosynthesis pathway. Focusing on the evolution of sulfatases and sulfotransferases, which control the mouth form plasticity in P. pacificus, reveals differences in copy number and genomic configurations between the genera Pristionchus and Parapristionchus. Altogether, this demonstrates the utility of the P. giblindavisi genome to date and polarizes lineage-specific patterns., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America.)

Because of ongoing climate change, populations of organisms are being subjected to stressful temperatures more often. This is especially problematic for ectothermic organisms, which are likely to be more sensitive to changes in temperature. Therefore, we need to know if ectotherms have adapted to environmental temperature and, if so, what are the evolutionary mechanisms behind such adaptation. Here, we use the nematode Pristionchus pacificus as a case study to investigate thermal adaptation on the Indian Ocean island of La Réunion, which experiences a range of temperatures from coast to summit. We study the evolution of high-temperature tolerance by constructing a phylogenetic tree of strains collected from many different thermal niches. We show that populations of P. pacificus at low altitudes have higher fertility at warmer temperatures. Most likely, this phenotype has arisen recently and at least twice independently, consistent with parallel evolution. We also studied low-temperature tolerance and showed that populations from high altitudes have increased their fertility at cooler temperatures. Together, these data indicate that P. pacificus strains on La Réunion are subject to divergent selection, adapting to hot and cold niches at the coast and summit of the volcano. Precisely defining these thermal niches provides essential information for models that predict the impact of future climate change on these populations., (© 2022 The Authors. Evolution published by Wiley Periodicals LLC on behalf of The Society for the Study of Evolution.)

Background: Repetitive sequences and mobile elements make up considerable fractions of individual genomes. While transposition events can be detrimental for organismal fitness, repetitive sequences form an enormous reservoir for molecular innovation. In this study, we aim to add repetitive elements to the annotation of the Pristionchus pacificus genome and assess their impact on novel gene formation., Results: Different computational approaches define up to 24% of the P. pacificus genome as repetitive sequences. While retroelements are more frequently found at the chromosome arms, DNA transposons are distributed more evenly. We found multiple DNA transposons, as well as LTR and LINE elements with abundant evidence of expression as single-exon transcripts. When testing whether transposons disproportionately contribute towards new gene formation, we found that roughly 10-20% of genes across all age classes overlap transposable elements with the strongest trend being an enrichment of low complexity regions among the oldest genes. Finally, we characterized a horizontal gene transfer of Zisupton elements into diplogastrid nematodes. These DNA transposons invaded nematodes from eukaryotic donor species and experienced a recent burst of activity in the P. pacificus lineage., Conclusions: The comprehensive annotation of repetitive elements in the P. pacificus genome builds a resource for future functional genomic analyses as well as for more detailed investigations of molecular innovations., (© 2022. The Author(s).)

The birth of new genes is a major molecular innovation driving phenotypic diversity across all domains of life. Although repurposing of existing protein-coding material by duplication is considered the main process of new gene formation, recent studies have discovered thousands of transcriptionally active sequences as a rich source of new genes. However, differential loss rates have to be assumed to reconcile the high birth rates of these incipient de novo genes with the dominance of ancient gene families in individual genomes. Here, we test this rapid turnover hypothesis in the context of the nematode model organism Pristionchus pacificus We extended the existing species-level phylogenomic framework by sequencing the genomes of six divergent P. pacificus strains. We used these data to study the evolutionary dynamics of different age classes and categories of origin at a population level. Contrasting de novo candidates with new families that arose by duplication and divergence from known genes, we find that de novo candidates are typically shorter, show less expression, and are overrepresented on the sex chromosome. Although the contribution of de novo candidates increases toward young age classes, multiple comparisons within the same age class showed significantly higher attrition in de novo candidates than in known genes. Similarly, young genes remain under weak evolutionary constraints with de novo candidates representing the fastest evolving subcategory. Altogether, this study provides empirical evidence for the rapid turnover hypothesis and highlights the importance of the evolutionary timescale when quantifying the contribution of different mechanisms toward new gene formation., (© 2022 Prabh and Rödelsperger; Published by Cold Spring Harbor Laboratory Press.)

High-resolution spatial and temporal maps of gene expression have facilitated a comprehensive understanding of animal development and evolution. In nematodes, the small body size represented a major challenge for such studies, but recent advancements have helped overcome this limitation. Here, we have implemented single worm transcriptomics (SWT) in the nematode model organism Pristionchus pacificus to provide a high-resolution map of the developmental transcriptome. We selected 38 time points from hatching of the J2 larvae to young adults to perform transcriptome analysis over 60 h of postembryonic development. A mean sequencing depth of 4.5 million read pairs allowed the detection of more than 23,135 (80%) of all genes. Nearly 3000 (10%) genes showed oscillatory expression with discrete expression levels, phases, and amplitudes. Gene age analysis revealed an overrepresentation of ancient gene classes among oscillating genes, and around one-third of them have 1:1 orthologs in C. elegans One important gene family overrepresented among oscillating genes is collagens. Several of these collagen genes are regulated by the developmental switch gene eud-1 , indicating a potential function in the regulation of mouth-form plasticity, a key developmental process in this facultative predatory nematode. Together, our analysis provides (1) an updated protocol for SWT in nematodes that is applicable to many microscopic species, (2) a 1- to 2-h high-resolution catalog of P. pacificus gene expression throughout postembryonic development, and (3) a comparative analysis of oscillatory gene expression between the two model organisms P. pacificus and C. elegans and associated evolutionary dynamics., (© 2021 Sun et al.; Published by Cold Spring Harbor Laboratory Press.)

Background: The nematode Pristionchus pacificus is an established model organism for comparative studies with Caenorhabditis elegans. Over the past years, it developed into an independent animal model organism for elucidating the genetic basis of phenotypic plasticity. Community-based curations were employed recently to improve the quality of gene annotations of P. pacificus and to more easily facilitate reverse genetic studies using candidate genes from C. elegans., Results: Here, I demonstrate that the reannotation of phylogenomic data from nine related nematode species using the community-curated P. pacificus gene set as homology data substantially improves the quality of gene annotations. Benchmarking of universal single copy orthologs (BUSCO) estimates a median completeness of 84% which corresponds to a 9% increase over previous annotations. Nevertheless, the ability to infer gene models based on homology already drops beyond the genus level reflecting the rapid evolution of nematode lineages. This also indicates that the highly curated C. elegans genome is not optimally suited for annotating non-Caenorhabditis genomes based on homology. Furthermore, comparative genomic analysis of apparently missing BUSCO genes indicates a failure of ortholog detection by the BUSCO pipeline due to the insufficient sample size and phylogenetic breadth of the underlying OrthoDB data set. As a consequence, the quality of multiple divergent nematode genomes might be underestimated., Conclusions: This study highlights the need for optimizing gene annotation protocols and it demonstrates the benefit of a high quality genome for phylogenomic data of related species.

The model organism Pristionchus pacificus and the genus Pristionchus, Kreis, 1932 have been intensively studied in the last decade with contemporary work focusing on the development, evolution, ecology, behavior, neurobiology, and genomics of this group of organisms. In particular, mechanistic studies on the development and evolution of mouth-form plasticity, predation and associated self-recognition processes enabled unique insight into life history strategies and the evolution of novelty. These studies include a comparative research agenda making use of the 39 available species of Pristionchus, all of which can be studied in living cultures. Sampling efforts revealed that Asia represents a biodiversity hotspot for Pristionchus worms. However, previous samplings have a bias towards northern and island areas, largely for logistic reasons. Here, we report on two extensive sampling trips to the Yunnan and Shaanxi provinces in Mainland China. We report the isolation of nine new Pristionchus species by morphology, morphometrics, mating experiments and genome-wide sequence analysis.

Background: Nematode model organisms such as Caenorhabditis elegans and Pristionchus pacificus are powerful systems for studying the evolution of gene function at a mechanistic level. However, the identification of P. pacificus orthologs of candidate genes known from C. elegans is complicated by the discrepancy in the quality of gene annotations, a common problem in nematode and invertebrate genomics., Results: Here, we combine comparative genomic screens for suspicious gene models with community-based curation to further improve the quality of gene annotations in P. pacificus. We extend previous curations of one-to-one orthologs to larger gene families and also orphan genes. Cross-species comparisons of protein lengths, screens for atypical domain combinations and species-specific orphan genes resulted in 4311 candidate genes that were subject to community-based curation. Corrections for 2946 gene models were implemented in a new version of the P. pacificus gene annotations. The new set of gene annotations contains 28,896 genes and has a single copy ortholog completeness level of 97.6%., Conclusions: Our work demonstrates the effectiveness of comparative genomic screens to identify suspicious gene models and the scalability of community-based approaches to improve the quality of thousands of gene models. Similar community-based approaches can help to improve the quality of gene annotations in other invertebrate species, including parasitic nematodes.

Environment shapes development through a phenomenon called developmental plasticity. Deciphering its genetic basis has potential to shed light on the origin of novel traits and adaptation to environmental change. However, molecular studies are scarce, and little is known about molecular mechanisms associated with plasticity. We investigated the gene regulatory network controlling predatory vs. non-predatory dimorphism in the nematode Pristionchus pacificus and found that it consists of genes of extremely different age classes. We isolated mutants in the conserved nuclear hormone receptor nhr-1 with previously unseen phenotypic effects. They disrupt mouth-form determination and result in animals combining features of both wild-type morphs. In contrast, mutants in another conserved nuclear hormone receptor nhr-40 display altered morph ratios, but no intermediate morphology. Despite divergent modes of control, NHR-1 and NHR-40 share transcriptional targets, which encode extracellular proteins that have no orthologs in Caenorhabditis elegans and result from lineage-specific expansions. An array of transcriptional reporters revealed co-expression of all tested targets in the same pharyngeal gland cell. Major morphological changes in this gland cell accompanied the evolution of teeth and predation, linking rapid gene turnover with morphological innovations. Thus, the origin of feeding plasticity involved novelty at the level of genes, cells and behavior., Competing Interests: The authors have declared that no competing interests exist.

Nematodes such as Caenorhabditis elegans are powerful systems to study basically all aspects of biology. Their species richness together with tremendous genetic knowledge from C. elegans facilitate the evolutionary study of biological functions using reverse genetics. However, the ability to identify orthologs of candidate genes in other species can be hampered by erroneous gene annotations. To improve gene annotation in the nematode model organism Pristionchus pacificus, we performed a genome-wide screen for C. elegans genes with potentially incorrectly annotated P. pacificus orthologs. We initiated a community-based project to manually inspect more than two thousand candidate loci and to propose new gene models based on recently generated Iso-seq and RNA-seq data. In most cases, misannotation of C. elegans orthologs was due to artificially fused gene predictions and completely missing gene models. The community-based curation raised the gene count from 25,517 to 28,036 and increased the single copy ortholog completeness level from 86% to 97%. This pilot study demonstrates how even small-scale crowdsourcing can drastically improve gene annotations. In future, similar approaches can be used for other species, gene sets, and even larger communities thus making manual annotation of large parts of the genome feasible.

Homology is a fundamental concept in comparative biology. It is extensively used at the sequence level to make phylogenetic hypotheses and functional inferences. Nonetheless, the majority of eukaryotic genomes contain large numbers of orphan genes lacking homologs in other taxa. Generally, the fraction of orphan genes is higher in genomically undersampled clades, and in the absence of closely related genomes any hypothesis about their origin and evolution remains untestable. Previously, we sequenced ten genomes with an underlying ladder-like phylogeny to establish a phylogenomic framework for studying genome evolution in diplogastrid nematodes. Here, we use this deeply sampled data set to understand the processes that generate orphan genes in our focal species Pristionchus pacificus Based on phylostratigraphic analysis and additional bioinformatic filters, we obtained 29 high-confidence candidate genes for which mechanisms of orphan origin were proposed based on manual inspection. This revealed diverse mechanisms including annotation artifacts, chimeric origin, alternative reading frame usage, and gene splitting with subsequent gain of de novo exons. In addition, we present two cases of complete de novo origination from non-coding regions, which represents one of the first reports of de novo genes in nematodes. Thus, we conclude that de novo emergence, divergence, and mixed mechanisms contribute to novel gene formation in Pristionchus nematodes., (Copyright © 2019 Prabh, Rodelsperger.)

Self-recognition is observed abundantly throughout the natural world, regulating diverse biological processes. Although ubiquitous, often little is known of the associated molecular machinery, and so far, organismal self-recognition has never been described in nematodes. We investigated the predatory nematode Pristionchus pacificus and, through interactions with its prey, revealed a self-recognition mechanism acting on the nematode surface, capable of distinguishing self-progeny from closely related strains. We identified the small peptide SELF-1, which is composed of an invariant domain and a hypervariable C terminus, as a key component of self-recognition. Modifications to the hypervariable region, including single-amino acid substitutions, are sufficient to eliminate self-recognition. Thus, the P. pacificus self-recognition system enables this nematode to avoid cannibalism while promoting the killing of competing nematodes., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Cilia are complex organelles involved in sensory perception and motility with intraflagellar transport (IFT) proteins being essential for cilia assembly and function, but little is known about cilia in an evo-devo context. For example, recent comparisons revealed conservation and divergence of IFT components in the regulation of social feeding behaviors between the nematodes Caenorhabditis elegans and Pristionchus pacificus. Here, we focus on the P. pacificus RFX transcription factor daf-19, the master regulator of ciliogenesis in C. elegans. Two CRISPR/Cas9-induced Ppa-daf-19 mutants lack ciliary structures in amphid neurons and display chemosensory defects. In contrast to IFT mutants, Ppa-daf-19 mutants do not exhibit social behavior. However, they show weak locomotive responses to shifts in oxygen concentration, suggesting partial impairment in sensing or responding to oxygen. To identify targets of Ppa-daf-19 regulation we compared the transcriptomes of Ppa-daf-19 and wild-type animals and performed a bioinformatic search for the X-box RFX binding-site across the genome. The regulatory network of Ppa-DAF-19 involves IFT genes but also many taxonomically restricted genes. We identified a conserved X-box motif as the putative binding site, which was validated for the Ppa-dyf-1 gene. Thus, Ppa-DAF-19 controls ciliogenesis, influences oxygen-induced behaviors and displays a high turnover of its regulatory network., (© 2018 Wiley Periodicals, Inc.)

The widespread identification of genes without detectable homology in related taxa is a hallmark of genome sequencing projects in animals, together with the abundance of gene duplications. Such genes have been called novel, young, taxon-restricted, or orphans, but little is known about the mechanisms accounting for their origin, age, and mode of evolution. Phylogenomic studies relying on deep and systematic taxon sampling and using the comparative method can provide insight into the evolutionary dynamics acting on novel genes. We used a phylogenomic approach for the nematode model organism Pristionchus pacificus and sequenced six additional Pristionchus and two outgroup species. This resulted in 10 genomes with a ladder-like phylogeny, sequenced in one laboratory using the same platform and analyzed by the same bioinformatic procedures. Our analysis revealed that 68%-81% of genes are assignable to orthologous gene families, the majority of which defined nine age classes with presence/absence patterns that can be explained by single evolutionary events. Contrasting different age classes, we find that older age classes are concentrated at chromosome centers, whereas novel gene families preferentially arise at the periphery, are weakly expressed, evolve rapidly, and have a high propensity of being lost. Over time, they increase in expression and become more constrained. Thus, the detailed phylogenetic resolution allowed a comprehensive characterization of the evolutionary dynamics of Pristionchus genomes indicating that distribution of age classes and their associated differences shape chromosomal divergence. This study establishes the Pristionchus system for future research on the mechanisms that drive the formation of novel genes., (© 2018 Prabh et al.; Published by Cold Spring Harbor Laboratory Press.)

Mutation and recombination are main drivers of phenotypic diversity, but the ability to create new allelic combinations is strongly dependent on the mode of reproduction. While most animals are dioecious (i.e., separated male and female sexes), in a number of evolutionary lineages females have gained the ability to self-fertilize [1, 2], with drastic consequences on effective recombination rate, genetic diversity, and the efficacy of selection [3]. In the genus Caenorhabditis, such hermaphroditic or androdioecious lineages, including C. briggsae and C. tropicalis, display a genome shrinkage relative to their dioecious sister species C. nigoni and C. brenneri, respectively [4, 5]. However, common consequences of reproductive modes on nematode genomes remain unknown, because most taxa contain single or few androdioecious species. One exception is the genus Pristionchus, with seven androdioecious species. Pristionchus worms are found in association with scarab beetles in worldwide samplings, resulting in deep taxon sampling and currently 39 culturable and available species. Here, we use phylotranscriptomics of all 39 Pristionchus species to provide a robust phylogeny based on an alignment of more than 2,000 orthologous clusters, which indicates that the seven androdioecious species represent six independent lineages. We show that gene loss is more prevalent in all hermaphroditic lineages than in dioecious relatives and that the majority of lost genes evolved recently in the Pristionchus genus. Further, we provide evidence that genes with male-biased expression are preferentially lost in hermaphroditic lineages. This supports a contribution of adaptive gene loss to shaping nematode genomes following the evolution of hermaphroditism., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Social behaviours are frequently utilised for defence and stress avoidance in nature. Both Caenorhabditis elegans and Pristionchus pacificus nematodes display social behaviours including clumping and bordering, to avoid hyperoxic stress conditions. Additionally, both species show natural variation in social behaviours with "social" and "solitary" strains. While the single solitary C. elegans N2 strain has evolved under laboratory domestication due to a gain-of-function mutation in the neuropeptide receptor gene npr-1, P. pacificus solitary strains are commonplace and likely ancestral. P. pacificus therefore provides an opportunity to further our understanding of the mechanisms regulating these complex behaviours and how they evolved within an ecologically relevant system. Using CRISPR/Cas9 engineering, we show that Ppa-npr-1 has minimal influence on social behaviours, indicating independent evolutionary pathways compared to C. elegans. Furthermore, solitary P. pacificus strains show an unexpected locomotive response to hyperoxic conditions, suggesting a novel regulatory mechanism counteracting social behaviours. By utilising both forward and reverse genetic approaches we identified 10 genes of the intraflagellar transport machinery in ciliated neurons that are essential for this inhibition. Therefore, a novel cilia-mediated environmental input adds an additional level of complexity to the regulation of hyperoxia-induced social behaviours in P. pacificus, a mechanism unknown in C. elegans.

Understanding how new species form requires investigation of evolutionary forces that cause phenotypic and genotypic changes among populations. However, the mechanisms underlying speciation vary and little is known about whether genomes diversify in the same ways in parallel at the incipient scale. We address this using the nematode, Pristionchus pacificus, which resides at an interesting point on the speciation continuum (distinct evolutionary lineages without reproductive isolation), and inhabits heterogeneous environments subject to divergent environmental pressures. Using whole genome re-sequencing of 264 strains, we estimate FST to identify outlier regions of extraordinary differentiation (∼1.725 Mb of the 172.5 Mb genome). We find evidence for shared divergent genomic regions occurring at a higher frequency than expected by chance among populations of the same evolutionary lineage. We use allele frequency spectra to find that, among lineages, 53% of divergent regions are consistent with adaptive selection, whereas 24% and 23% of such regions suggest background selection and restricted gene flow, respectively. In contrast, among populations from the same lineage, similar proportions (34-48%) of divergent regions correspond to adaptive selection and restricted gene flow, whereas 13-22% suggest background selection. Because speciation often involves phenotypic and genomic divergence, we also evaluate phenotypic variation, focusing on pH tolerance, which we find is diverging in a manner corresponding to environmental differences among populations. Taking a genome-wide association approach, we functionally validate a significant genotype-phenotype association for this trait. Our results are consistent with P. pacificus undergoing heterogeneous genotypic and phenotypic diversification related to both evolutionary and environmental processes., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Developmental plasticity, the ability of one genotype to produce distinct phenotypes in different environments, has been suggested to facilitate phenotypic diversification, and several examples in plants and animals support its macroevolutionary potential [1-8]. However, little is known about associated molecular mechanisms, because environmental effects on development are difficult to study by laboratory approaches. One promising system is the mouth dimorphism of the nematode Pristionchus pacificus [9-12]. Following an irreversible decision in larval development, these nematodes form moveable teeth that occur in either of two discrete morphs. The "eurystomatous" (Eu) form has a wide mouth and two teeth, allowing predatory feeding on other nematodes. In contrast, the alternative ("stenostomatous"; St) form has diminutive mouthparts that largely constrain its diet to microbes. The sulfatase EUD-1 was previously discovered to execute a polyphenism switch based on dosage of functional alleles [13] and confirmed a prediction of evolutionary theory about how developmental switches control plasticity [1, 3]. However, the genetic context of this single gene, and hence the molecular complexity of switch mechanisms, was previously unknown. Here we use a suppressor screen to identify factors downstream of eud-1 in mouth-form regulation. We isolated three dominant, X-linked mutants in the nuclear hormone receptor gene nhr-40 that are haploinsufficient. Both eud-1 nhr-40 double and nhr-40 single mutants are all Eu, whereas transgenic overexpression of nhr-40 does not restore the wild-type phenotype but instead results in nearly all-St lines. Thus, NHR-40 is part of a developmental switch, suggesting that switch mechanisms controlling plasticity consist of multi-component hormonal signaling systems., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Wild isolates of the nematode Caenorhabditis elegans perform social behaviours, namely clumping and bordering, to avoid hyperoxia under laboratory conditions. In contrast, the laboratory reference strain N2 has acquired a solitary behaviour in the laboratory, related to a gain-of-function variant in the neuropeptide Y-like receptor NPR-1. Here, we study the evolution and natural variation of clumping and bordering behaviours in Pristionchus pacificus nematodes in a natural context, using strains collected from 22 to 2400 metres above sea level on La Réunion Island. Through the analysis of 106 wild isolates, we show that the majority of strains display a solitary behaviour similar to C. elegans N2, whereas social behaviours are predominantly seen in strains that inhabit high-altitude locations. We show experimentally that P. pacificus social strains perform clumping and bordering to avoid hyperoxic conditions in the laboratory, suggesting that social strains may have adapted to or evolved a preference for the lower relative oxygen levels available at high altitude in nature. In contrast to C. elegans, clumping and bordering in P. pacificus do not correlate with locomotive behaviours in response to changes in oxygen conditions. Furthermore, QTL analysis indicates clumping and bordering to represent complex quantitative traits. Thus, clumping and bordering behaviours represent an example of phenotypic convergence with a different evolutionary history and distinct genetic control in both nematode species., (© 2016 The Author(s).)

The hermaphroditic nematode Pristionchus pacificus is an established model system for comparative studies with Caenorhabditis elegans in developmental biology, ecology, and population genetics. In this study, we present whole-genome sequencing data of 104 P. pacificus strains and the draft assembly of the obligate outcrossing sister species P. exspectatus. We characterize genetic diversity within P. pacificus and investigate the population genetic processes shaping this diversity. P. pacificus is 10 times more diverse than C. elegans and exhibits substantial population structure that allows us to probe its evolution on multiple timescales. Consistent with reduced effective recombination in this self-fertilizing species, we find haplotype blocks that span several megabases. Using the P. exspectatus genome as an outgroup, we polarized variation in P. pacificus and found a site frequency spectrum (SFS) that decays more rapidly than expected in neutral models. The SFS at putatively neutral sites is U shaped, which is a characteristic feature of pervasive linked selection. Based on the additional findings (i) that the majority of nonsynonymous variation is eliminated over timescales on the order of the separation between clades, (ii) that diversity is reduced in gene-rich regions, and (iii) that highly differentiated clades show very similar patterns of diversity, we conclude that purifying selection on many mutations with weak effects is a major force shaping genetic diversity in P. pacificus.

Base substitution mutations are a major source of genetic novelty and mutation accumulation line (MAL) studies revealed a nearly universal AT bias in de novo mutation spectra. While a comparison of de novo mutation spectra with the actual nucleotide composition in the genome suggests the existence of general counterbalancing mechanisms, little is known about the evolutionary and historical details of these opposing forces. Here, we correlate MAL-derived mutation spectra with patterns observed from population resequencing. Variation observed in natural populations has already been subject to evolutionary forces. Distinction between rare and common alleles, the latter of which are close to fixation and of presumably older age, can provide insight into mutational processes and their influence on genome evolution. We provide a genome-wide analysis of de novo mutations in 22 MALs of the nematode Pristionchus pacificus and compare the spectra with natural variants observed in resequencing of 104 natural isolates. MALs show an AT bias of 5.3, one of the highest values observed to date. In contrast, the AT bias in natural variants is much lower. Specifically, rare derived alleles show an AT bias of 2.4, whereas common derived alleles close to fixation show no AT bias at all. These results indicate the existence of a strong opposing force and they suggest that the GC content of the P. pacificus genome is in equilibrium. We discuss GC-biased gene conversion as a potential mechanism acting against AT-biased mutations. This study provides insight into genome evolution by combining MAL studies with natural variation.