Recombinational landscape and population genomics of Caenorhabditis elegans

Hill-Robertson interference may bias the inference of fitness effects of new mutations in highly selfing species

The accurate estimation of the distribution of fitness effects (DFE) of new mutations is critical for population genetic inference but remains a challenging task. While various methods have been developed for DFE inference using the site frequency spectrum of putatively neutral and selected sites, their applicability in species with diverse life history traits and complex demographic scenarios is not well understood. Selfing is common among eukaryotic species and can lead to decreased effective recombination rates, increasing the effects of selection at linked sites, including interference between selected alleles. We employ forward simulations to investigate the limitations of current DFE estimation approaches in the presence of selfing and other model violations, such as linkage, departures from semidominance, population structure, and uneven sampling. We find that distortions of the site frequency spectrum due to Hill-Robertson interference in highly selfing populations lead to mis-inference of the deleterious DFE of new mutations. Specifically, when inferring the distribution of selection coefficients, there is an overestimation of nearly neutral and strongly deleterious mutations and an underestimation of mildly deleterious mutations when interference between selected alleles is pervasive. In addition, the presence of cryptic population structure with low rates of migration and uneven sampling across subpopulations leads to the false inference of a deleterious DFE skewed towards effectively neutral/mildly deleterious mutations. Finally, the proportion of adaptive substitutions estimated at high rates of selfing is substantially overestimated. Our observations apply broadly to species and genomic regions with little/no recombination and where interference might be pervasive.

Regulatory logic and transposable element dynamics in nematode worm genomes

Genome sequencing has revealed a tremendous diversity of transposable elements (TEs) in eukaryotes but there is little understanding of the evolutionary processes responsible for TE diversity. Non-autonomous TEs have lost the machinery necessary for transposition and rely on closely related autonomous TEs for critical proteins. We studied two mathematical models of TE regulation, one assuming that both autonomous tranposons and their non-autonomous relatives operate under the same regulatory logic, competing for transposition resources, and one assuming that autonomous TEs self-attenuate transposition while non-autonomous transposons continually increase, parasitizing their autonomous relatives. We implemented these models in stochastic simulations and studied how TE regulatory relationships influence transposons and populations. We found that only outcrossing populations evolving with Parasitic TE regulation resulted in stable maintenance of TEs. We tested our model predictions in Caenorhabditis genomes by annotating TEs in two focal families, autonomous LINEs and their non-autonomous SINE relatives and the DNA transposon Mutator . We found broad variation in autonomous - non-autonomous relationships and rapid mutational decay in the sequences that allow non-autonomous TEs to transpose. Together, our results suggest that individual TE families evolve according to disparate regulatory rules that are relevant in the early, acute stages of TE invasion.

The chromatin-associated 53BP1 ortholog, HSR-9, regulates recombinational repair and X chromosome segregation in the Caenorhabditis elegans germ line

53BP1 plays a crucial role in regulating DNA damage repair pathway choice and checkpoint signaling in somatic cells; however, its role in meiosis has remained enigmatic. In this study, we demonstrate that the Caenorhabditis elegans ortholog of 53BP1, HSR-9, associates with chromatin in both proliferating and meiotic germ cells. Notably, HSR-9 is enriched on the X chromosome pair in pachytene oogenic germ cells. HSR-9 is also present at kinetochores during both mitotic and meiotic divisions but does not appear to be essential for monitoring microtubule-kinetochore attachments or tension. Using cytological markers of different steps in recombinational repair, we found that HSR-9 influences the processing of a subset of meiotic double-stranded breaks into COSA-1-marked crossovers. Additionally, HSR-9 plays a role in meiotic X chromosome segregation under conditions where X chromosomes fail to pair, synapse, and recombine. Together, these results highlight that chromatin-associated HSR-9 has both conserved and unique functions in the regulation of meiotic chromosome behavior.

The recombination landscape of introgression in yeast

Meiotic recombination is an evolutionary force that acts by breaking up genomic linkage, increasing the efficacy of selection. Recombination is initiated with a double-strand break which is resolved via a crossover, which involves the reciprocal exchange of genetic material between homologous chromosomes, or a non-crossover, which results in small tracts of non-reciprocal exchange of genetic material. Crossover and non-crossover rates vary between species, populations, individuals, and across the genome. In recent years, recombination rate has been associated with the distribution of ancestry derived from past interspecific hybridization (introgression) in a variety of species. We explore this interaction of recombination and introgression by sequencing spores and detecting crossovers and non-crossovers from two crosses of the yeast Saccharomyces uvarum. One cross is between strains which each contain introgression from their sister species, S. eubayanus, while the other cross has no introgression present. We find that the recombination landscape is significantly different between S. uvarum crosses, and that some of these differences can be explained by the presence of introgression in one cross. Crossovers are reduced and non-crossovers are increased in heterozygous introgression compared to syntenic regions in the cross without introgression. This translates to reduced allele shuffling within introgressed regions, and an overall reduction of shuffling on most chromosomes with introgression compared to the syntenic regions and chromosomes without introgression. Our results suggest that hybridization can significantly influence the recombination landscape, and that the reduction in allele shuffling contributes to the initial purging of introgression in the generations following a hybridization event.

Pervasive conservation of intron number and other genetic elements revealed by a chromosome-level genomic assembly of the hyper-polymorphic nematode Caenorhabditis brenneri

With within-species genetic diversity estimates that span the gambit of that seen across the entirety of animals, the Caenorhabditis genus of nematodes holds unique potential to provide insights into how population size and reproductive strategies influence gene and genome organization and evolution. Our study focuses on Caenorhabditis brenneri, currently known as one of the most genetically diverse nematodes within its genus and metazoan phyla. Here, we present a high-quality gapless genome assembly and annotation for C. brenneri, revealing a common nematode chromosome arrangement characterized by gene-dense central regions and repeat rich peripheral parts. Comparison of C. brenneri with other nematodes from the 'Elegans' group revealed conserved macrosynteny but a lack of microsynteny, characterized by frequent rearrangements and low correlation iof orthogroup sizes, indicative of high rates of gene turnover. We also assessed genome organization within corresponding syntenic blocks in selfing and outcrossing species, affirming that selfing species predominantly experience loss of both genes and intergenic DNA. Comparison of gene structures revealed strikingly small number of shared introns across species, yet consistent distributions of intron number and length, regardless of population size or reproductive mode, suggesting that their evolutionary dynamics are primarily reflective of functional constraints. Our study provides valuable insights into genome evolution and expands the nematode genome resources with the highly genetically diverse C. brenneri, facilitating research into various aspects of nematode biology and evolutionary processes.

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).

The impact of differential transposition activities of autonomous and nonautonomous hAT transposable elements on genome architecture and gene expression in Caenorhabditis inopinata

Transposable elements are DNA sequences capable of moving within genomes and significantly influence genomic evolution. The nematode Caenorhabditis inopinata exhibits a much higher transposable element copy number than its sister species, Caenorhabditis elegans. In this study, we identified a novel autonomous transposable element belonging to the hAT superfamily from a spontaneous transposable element-insertion mutant in C. inopinata and named this transposon Ci-hAT1. Further bioinformatic analyses uncovered 3 additional autonomous hAT elements-Ci-hAT2, Ci-hAT3, and Ci-hAT4-along with over 1,000 copies of 2 nonautonomous miniature inverted-repeat transposable elements, mCi-hAT1 and mCi-hAT4, likely derived from Ci-hAT1 and Ci-hAT4 through internal deletion. We tracked at least 3 sequential transpositions of Ci-hAT1 over several years. However, the transposition rates of the other 3 autonomous hAT elements were lower, suggesting varying activity levels. Notably, the distribution patterns of the 2 miniature inverted-repeat transposable element families differed significantly: mCi-hAT1 was primarily located in the chromosome arms, a pattern observed in the transposable elements of other Caenorhabditis species, whereas mCi-hAT4 was more evenly distributed across chromosomes. Additionally, interspecific transcriptome analysis indicated that C. inopinata genes with upstream or intronic these miniature inverted-repeat transposable element insertions tend to be more highly expressed than their orthologous genes in C. elegans. These findings highlight the significant role of de-silenced transposable elements in driving the evolution of genomes and transcriptomes, leading to species-specific genetic diversity.

Chromosome segregation during spermatogenesis occurs through a unique center-kinetic mechanism in holocentric moth species

Precise regulation of chromosome dynamics in the germline is essential for reproductive success across species. Yet, the mechanisms underlying meiotic chromosomal events such as homolog pairing and chromosome segregation are not fully understood in many species. Here, we employ Oligopaint DNA FISH to investigate mechanisms of meiotic homolog pairing and chromosome segregation in the holocentric pantry moth, Plodia interpunctella, and compare our findings to new and previous studies in the silkworm moth, Bombyx mori, which diverged from P. interpunctella over 100 million years ago. We find that pairing in both Bombyx and Plodia spermatogenesis is initiated at gene-rich chromosome ends. Additionally, both species form rod shaped cruciform-like bivalents at metaphase I. However, unlike the telomere-oriented chromosome segregation mechanism observed in Bombyx, Plodia can orient bivalents in multiple different ways at metaphase I. Surprisingly, in both species we find that kinetochores consistently assemble at non-telomeric loci toward the center of chromosomes regardless of where chromosome centers are located in the bivalent. Additionally, sister kinetochores do not seem to be paired in these species. Instead, four distinct kinetochores are easily observed at metaphase I. Despite this, we find clear end-on microtubule attachments and not lateral microtubule attachments co-orienting these separated kinetochores. These findings challenge the classical view of segregation where paired, poleward-facing kinetochores are required for accurate homolog separation in meiosis I. Our studies here highlight the importance of exploring fundamental processes in non-model systems, as employing novel organisms can lead to the discovery of novel biology.

Chromosome fusion and programmed DNA elimination shape karyotypes of nematodes

An increasing number of metazoans undergo programmed DNA elimination (PDE), where a significant amount of DNA is selectively lost from the somatic genome during development. In some nematodes, PDE leads to the removal and remodeling of the ends of all germline chromosomes. In several species, PDE also generates internal breaks that lead to sequence loss and increased numbers of somatic chromosomes. The biological significance of these karyotype changes associated with PDE and the origin and evolution of nematode PDE remain largely unknown. Here, we assembled the single germline chromosome of the nematode Parascaris univalens and compared the karyotypes, chromosomal gene organization, and PDE features among other nematodes. We show that PDE in Parascaris converts an XX/XY sex-determination system in the germline into an XX/XO system in the somatic cells. Comparisons of Ascaris, Parascaris, and Baylisascaris ascarid chromosomes suggest that PDE existed in the ancestor of these nematodes, and their current distinct germline karyotypes were derived from fusion events of smaller ancestral chromosomes. The DNA breaks involved in PDE resolve these fused germline chromosomes into their pre-fusion karyotypes. These karyotype changes may lead to alterations in genome architecture and gene expression in the somatic cells. Cytological and genomic analyses further suggest that satellite DNA and the heterochromatic chromosome arms are dynamic and may play a role during meiosis. Overall, our results show that chromosome fusion and PDE have been harnessed in these ascarids to sculpt their karyotypes, altering the genome organization and serving specific functions in the germline and somatic cells.

The chromatin-associated 53BP1 ortholog, HSR-9, regulates recombinational repair and X chromosome segregation in the Caenorhabditis elegans germ line

53BP1 plays a crucial role in regulating DNA damage repair pathway choice and checkpoint signaling in somatic cells; however, its role in meiosis has remained enigmatic. In this study, we demonstrate that the Caenorhabditis elegans ortholog of 53BP1, HSR-9, associates with chromatin in both proliferating and meiotic germ cells. Notably, HSR-9 is enriched on the X chromosome pair in pachytene oogenic germ cells. HSR-9 is also present at kinetochores during both mitotic and meiotic divisions but does not appear to be essential for monitoring microtubule-kinetochore attachments or tension. Using cytological markers of different steps in recombinational repair, we found that HSR-9 influences the processing of a subset of meiotic double strand breaks into COSA-1-marked crossovers. Additionally, HSR-9 plays a role in meiotic X chromosome segregation under conditions where X chromosomes fail to pair, synapse, and recombine. Together, these results highlight that chromatin-associated HSR-9 has both conserved and unique functions in the regulation of meiotic chromosome behavior.

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.

Meiotic recombination dynamics in plants with repeat-based holocentromeres shed light on the primary drivers of crossover patterning

Centromeres strongly affect (epi)genomic architecture and meiotic recombination dynamics, influencing the overall distribution and frequency of crossovers. Here we show how recombination is regulated and distributed in the holocentric plant Rhynchospora breviuscula, a species with diffused centromeres. Combining immunocytochemistry, chromatin analysis and high-throughput single-pollen sequencing, we discovered that crossover frequency is distally biased, in sharp contrast to the diffused distribution of hundreds of centromeric units and (epi)genomic features. Remarkably, we found that crossovers were abolished inside centromeric units but not in their proximity, indicating the absence of a canonical centromere effect. We further propose that telomere-led synapsis of homologues is the feature that best explains the observed recombination landscape. Our results hint at the primary influence of mechanistic features of meiotic pairing and synapsis rather than (epi)genomic features and centromere organization in determining the distally biased crossover distribution in R. breviuscula, whereas centromeres and (epi)genetic properties only affect crossover positioning locally.

Experimental evolution of hybrid populations to identify Dobzhansky-Muller incompatibility loci

Epistatic interactions between loci that reduce fitness in interspecies hybrids, Dobzhansky-Muller incompatibilities (DMIs), contribute genetically to the inviability and infertility within hybrid populations. It remains a challenge, however, to identify the loci that contribute to DMIs as causes of reproductive isolation between species. Here, we assess through forward simulation the power of evolve-and-resequence (E&R) experimental evolution of hybrid populations to map DMI loci. We document conditions under which such a mapping strategy may be most feasible and demonstrate how mapping power is sensitive to biologically relevant parameters such as one-way versus two-way incompatibility type, selection strength, recombination rate, and dominance interactions. We also assess the influence of parameters under direct control of an experimenter, including duration of experimental evolution and number of replicate populations. We conclude that an E&R strategy for mapping DMI loci, and other cases of epistasis, can be a viable option under some circumstances for study systems with short generation times like Caenorhabditis nematodes.

Patterns of Genomic Diversity in a Fig-Associated Close Relative of Caenorhabditis elegans

The evolution of reproductive mode is expected to have profound impacts on the genetic composition of populations. At the same time, ecological interactions can generate close associations among species, which can in turn generate a high degree of overlap in their spatial distributions. Caenorhabditis elegans is a hermaphroditic nematode that has enabled extensive advances in developmental genetics. Caenorhabditis inopinata, the sister species of C. elegans, is a gonochoristic nematode that thrives in figs and obligately disperses on fig wasps. Here, we describe patterns of genomic diversity in C. inopinata. We performed RAD-seq on individual worms isolated from the field across three Okinawan island populations. C. inopinata is about five times more diverse than C. elegans. Additionally, C. inopinata harbors greater differences in diversity among functional genomic regions (such as between genic and intergenic sequences) than C. elegans. Conversely, C. elegans harbors greater differences in diversity between high-recombining chromosome arms and low-recombining chromosome centers than C. inopinata. FST is low among island population pairs, and clear population structure could not be easily detected among islands, suggesting frequent migration of wasps between islands. These patterns of population differentiation appear comparable with those previously reported in its fig wasp vector. These results confirm many theoretical population genetic predictions regarding the evolution of reproductive mode and suggest C. inopinata population dynamics may be driven by wasp dispersal. This work sets the stage for future evolutionary genomic studies aimed at understanding the evolution of sex as well as the evolution of ecological interactions.

High-content phenotypic analysis of a C. elegans recombinant inbred population identifies genetic and molecular regulators of lifespan

Lifespan is influenced by complex interactions between genetic and environmental factors. Studying those factors in model organisms of a single genetic background limits their translational value for humans. Here, we mapped lifespan determinants in 85 genetically diverse C. elegans recombinant intercross advanced inbred lines (RIAILs). We assessed molecular profiles - transcriptome, proteome, and lipidome - and life-history traits, including lifespan, development, growth dynamics, and reproduction. RIAILs exhibited large variations in lifespan, which positively correlated with developmental time. Among the top candidates obtained from multi-omics data integration and QTL mapping, we validated known and novel longevity modulators, including rict-1, gfm-1 and mltn-1. We translated their relevance to humans using UK Biobank data and showed that variants in RICTOR and GFM1 are associated with an elevated risk of age-related heart disease, dementia, diabetes, kidney, and liver diseases. We organized our dataset as a resource (https://lisp-lms.shinyapps.io/RIAILs/) that allows interactive explorations for new longevity targets.

Comparative Genomics of Sex, Chromosomes, and Sex Chromosomes in Caenorhabditis elegans and Other Nematodes

The nematode phylum has evolved a remarkable diversity of reproductive modes, including the repeated emergence of asexuality and hermaphroditism across divergent clades. The species-richness and small genome size of nematodes make them ideal systems for investigating the genome-wide causes and consequences of such major transitions. The availability of functional annotations for most Caenorhabditis elegans genes further allows the linking of patterns of gene content evolution with biological processes. Such gene-centric studies were recently complemented by investigations of chromosome evolution that made use of the first chromosome-scale genome assemblies outside the Caenorhabditis genus. This review highlights recent comparative genomic studies of reproductive mode evolution addressing the hybrid origin of asexuality and the parallel gene loss following the emergence of hermaphroditism. It further summarizes ongoing efforts to characterize ancient linkage blocks called Nigon elements, which form central units of chromosome evolution. Fusions between Nigon elements have been demonstrated to impact recombination and speciation. Finally, multiple recent fusions between autosomal and the sex-linked Nigon element reveal insights into the dynamic evolution of sex chromosomes across various timescales.

Chromosome fusion and programmed DNA elimination shape karyotypes of parasitic nematodes

A growing list of metazoans undergo programmed DNA elimination (PDE), where a significant amount of DNA is selectively lost from the somatic genome during development. In some nematodes, PDE leads to the removal and remodeling of the ends of all germline chromosomes. In several species, PDE also generates internal breaks that lead to sequence loss and an increased number of somatic chromosomes. The biological significance of these karyotype changes associated with PDE and the origin and evolution of nematode PDE remain largely unknown. Here, we assembled the single germline chromosome of the horse parasite Parascaris univalens and compared the karyotypes, chromosomal gene organization, and PDE features among ascarid nematodes. We show that PDE in Parascaris converts an XX/XY sex-determination system in the germline into an XX/XO system in the somatic cells. Comparisons of Ascaris, Parascaris, and Baylisascaris ascarid chromosomes suggest that PDE existed in the ancestor of these parasites, and their current distinct germline karyotypes were derived from fusion events of smaller ancestral chromosomes. The DNA breaks involved in PDE resolve these fused germline chromosomes into their pre-fusion karyotypes, leading to alterations in genome architecture and gene expression in the somatic cells. Cytological and genomic analyses further suggest that satellite DNA and the heterochromatic chromosome arms play a dynamic role in the Parascaris germline chromosome during meiosis. Overall, our results show that chromosome fusion and PDE have been harnessed in these ascarids to sculpt their karyotypes, altering the genome organization and serving specific functions in the germline and somatic cells.

Genome-wide association and environmental suppression of the mortal germline phenotype of wild C. elegans

The animal germline lineage needs to be maintained along generations. However, some Caenorhabditis elegans wild isolates display a mortal germline phenotype, leading to sterility after several generations at 25°C. Using a genome-wide association approach, we detect a significant peak on chromosome III around 5 Mb, confirmed by introgressions. Thus, a seemingly deleterious genotype is maintained at intermediate frequency in the species. Environmental rescue is a likely explanation, and indeed associated bacteria and microsporidia suppress the phenotype of wild isolates as well as mutants in small RNA inheritance (nrde-2) and histone modifications (set-2). Escherichia coli strains of the K-12 lineage suppress the phenotype compared to B strains. By shifting a wild strain from E. coli K-12 to E. coli B, we find that memory of the suppressing condition is maintained over several generations. Thus, the mortal germline phenotype of wild C. elegans is in part revealed by laboratory conditions and may represent variation in epigenetic inheritance and environmental interactions. This study also points to the importance of non-genetic memory in the face of environmental variation.

Ancient diversity in host-parasite interaction genes in a model parasitic nematode

Host-parasite interactions exert strong selection pressures on the genomes of both host and parasite. These interactions can lead to negative frequency-dependent selection, a form of balancing selection that is hypothesised to explain the high levels of polymorphism seen in many host immune and parasite antigen loci. Here, we sequence the genomes of several individuals of Heligmosomoides bakeri, a model parasite of house mice, and Heligmosomoides polygyrus, a closely related parasite of wood mice. Although H. bakeri is commonly referred to as H. polygyrus in the literature, their genomes show levels of divergence that are consistent with at least a million years of independent evolution. The genomes of both species contain hyper-divergent haplotypes that are enriched for proteins that interact with the host immune response. Many of these haplotypes originated prior to the divergence between H. bakeri and H. polygyrus, suggesting that they have been maintained by long-term balancing selection. Together, our results suggest that the selection pressures exerted by the host immune response have played a key role in shaping patterns of genetic diversity in the genomes of parasitic nematodes.

Comprehensive detection of structural variation and transposable element differences between wild type laboratory lineages of C. elegans

Genomic structural variations (SVs) and transposable elements (TEs) can be significant contributors to genome evolution, altered gene expression, and risk of genetic diseases. Recent advancements in long-read sequencing have greatly improved the quality of de novo genome assemblies and enhanced the detection of sequence variants at the scale of hundreds or thousands of bases. Comparisons between two diverged wild isolates of Caenorhabditis elegans, the Bristol and Hawaiian strains, have been widely utilized in the analysis of small genetic variations. Genetic drift, including SVs and rearrangements of repeated sequences such as TEs, can occur over time from long-term maintenance of wild type isolates within the laboratory. To comprehensively detect both large and small structural variations as well as TEs due to genetic drift, we generated de novo genome assemblies and annotations for each strain from our lab collection using both long- and short-read sequencing and compared our assemblies and annotations with that of other lab wild type strains. Within our lab assemblies, we annotate over 3.1Mb of sequence divergence between the Bristol and Hawaiian isolates: 337,584 SNPs, 94,503 small insertion-deletions (<50bp), and 4,334 structural variations (>50bp). Further, we define the location and movement of specific DNA TEs between N2 Bristol and CB4856 Hawaiian wild type isolates. Specifically, we find the N2 Bristol genome has 20.6% more TEs from the Tc1/mariner family than the CB4856 Hawaiian genome. Moreover, we identified Zator elements as the most abundant and mobile TE family in the genome. Using specific TE sequences with unique SNPs, we also identify 38 TEs that moved intrachromosomally and 9 TEs that moved interchromosomally between the N2 Bristol and CB4856 Hawaiian genomes. By comparing the de novo genome assembly of our lab collection Bristol isolate to the VC2010 Bristol assembly, we also reveal that lab lineages display over 2 Mb of total variation: 1,162 SNPs, 1,528 indels, and 897 SVs with 95% of the variation due to SVs. Overall, our work demonstrates the unique contribution of SVs and TEs to variation and genetic drift between wild type laboratory strains assumed to be isogenic despite growing evidence of genetic drift and phenotypic variation.

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.

Effect of recombination on genetic diversity of Caenorhabditis elegans

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

Microbes are potential key players in the evolution of life histories and aging in Caenorhabditis elegans

Microbes can have profound effects on host fitness and health and the appearance of late-onset diseases. Host-microbe interactions thus represent a major environmental context for healthy aging of the host and might also mediate trade-offs between life-history traits in the evolution of host senescence. Here, we have used the nematode Caenorhabditis elegans to study how host-microbe interactions may modulate the evolution of life histories and aging. We first characterized the effects of two non-pathogenic and one pathogenic Escherichia coli strains, together with the pathogenic Serratia marcescens DB11 strain, on population growth rates and survival of C. elegans from five different genetic backgrounds. We then focused on an outbred C. elegans population, to understand if microbe-specific effects on the reproductive schedule and in traits such as developmental rate and survival were also expressed in the presence of males and standing genetic variation, which could be relevant for the evolution of C. elegans and other nematode species in nature. Our results show that host-microbe interactions have a substantial host-genotype-dependent impact on the reproductive aging and survival of the nematode host. Although both pathogenic bacteria reduced host survival in comparison with benign strains, they differed in how they affected other host traits. Host fertility and population growth rate were affected by S. marcescens DB11 only during early adulthood, whereas this occurred at later ages with the pathogenic E. coli IAI1. In both cases, these effects were largely dependent on the host genotypes. Given such microbe-specific genotypic differences in host life history, we predict that the evolution of reproductive schedules and senescence might be critically contingent on host-microbe interactions in nature.

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.

Genomic diversity landscapes in outcrossing and selfing Caenorhabditis nematodes

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

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.

Chromosomal Instability in Genome Evolution: From Cancer to Macroevolution

The integrity of the genome is crucial for the survival of all living organisms. However, genomes need to adapt to survive certain pressures, and for this purpose use several mechanisms to diversify. Chromosomal instability (CIN) is one of the main mechanisms leading to the creation of genomic heterogeneity by altering the number of chromosomes and changing their structures. In this review, we will discuss the different chromosomal patterns and changes observed in speciation, in evolutional biology as well as during tumor progression. By nature, the human genome shows an induction of diversity during gametogenesis but as well during tumorigenesis that can conclude in drastic changes such as the whole genome doubling to more discrete changes as the complex chromosomal rearrangement chromothripsis. More importantly, changes observed during speciation are strikingly similar to the genomic evolution observed during tumor progression and resistance to therapy. The different origins of CIN will be treated as the importance of double-strand breaks (DSBs) or the consequences of micronuclei. We will also explain the mechanisms behind the controlled DSBs, and recombination of homologous chromosomes observed during meiosis, to explain how errors lead to similar patterns observed during tumorigenesis. Then, we will also list several diseases associated with CIN, resulting in fertility issues, miscarriage, rare genetic diseases, and cancer. Understanding better chromosomal instability as a whole is primordial for the understanding of mechanisms leading to tumor progression.

A superior gene allele involved in abscisic acid signaling enhances drought tolerance and yield in chickpea

Identifying potential molecular tags for drought tolerance is essential for achieving higher crop productivity under drought stress. We employed an integrated genomics-assisted breeding and functional genomics strategy involving association mapping, fine mapping, map-based cloning, molecular haplotyping and transcript profiling in the introgression lines (ILs)- and near isogenic lines (NILs)-based association panel and mapping population of chickpea (Cicer arietinum). This combinatorial approach delineated a bHLH (basic helix-loop-helix) transcription factor, CabHLH10 (Cicer arietinum bHLH10) underlying a major QTL, along with its derived natural alleles/haplotypes governing yield traits under drought stress in chickpea. CabHLH10 binds to a cis-regulatory G-box promoter element to modulate the expression of RD22 (responsive to desiccation 22), a drought/abscisic acid (ABA)-responsive gene (via a trans-expression QTL), and two strong yield-enhancement photosynthetic efficiency (PE) genes. This, in turn, upregulates other downstream drought-responsive and ABA signaling genes, as well as yield-enhancing PE genes, thus increasing plant adaptation to drought with reduced yield penalty. We showed that a superior allele of CabHLH10 introgressed into the NILs improved root and shoot biomass and PE, thereby enhancing yield and productivity during drought without compromising agronomic performance. Furthermore, overexpression of CabHLH10 in chickpea and Arabidopsis (Arabidopsis thaliana) conferred enhanced drought tolerance by improving root and shoot agro-morphological traits. These findings facilitate translational genomics for crop improvement and the development of genetically tailored, climate-resilient, high-yielding chickpea cultivars.

Systematic characterization of chromodomain proteins reveals an H3K9me1/2 reader regulating aging in C. elegans

The chromatin organization modifier domain (chromodomain) is an evolutionally conserved motif across eukaryotic species. The chromodomain mainly functions as a histone methyl-lysine reader to modulate gene expression, chromatin spatial conformation and genome stability. Mutations or aberrant expression of chromodomain proteins can result in cancer and other human diseases. Here, we systematically tag chromodomain proteins with green fluorescent protein (GFP) using CRISPR/Cas9 technology in C. elegans. By combining ChIP-seq analysis and imaging, we delineate a comprehensive expression and functional map of chromodomain proteins. We then conduct a candidate-based RNAi screening and identify factors that regulate the expression and subcellular localization of the chromodomain proteins. Specifically, we reveal an H3K9me1/2 reader, CEC-5, both by in vitro biochemistry and in vivo ChIP assays. MET-2, an H3K9me1/2 writer, is required for CEC-5 association with heterochromatin. Both MET-2 and CEC-5 are required for the normal lifespan of C. elegans. Furthermore, a forward genetic screening identifies a conserved Arginine124 of CEC-5's chromodomain, which is essential for CEC-5's association with chromatin and life span regulation. Thus, our work will serve as a reference to explore chromodomain functions and regulation in C. elegans and allow potential applications in aging-related human diseases.

Chromosome fusions repatterned recombination rate and facilitated reproductive isolation during Pristionchus nematode speciation

Large-scale genome-structural evolution is common in various organisms. Recent developments in speciation genomics revealed the importance of inversions, whereas the role of other genome-structural rearrangements, including chromosome fusions, have not been well characterized. We study genomic divergence and reproductive isolation of closely related nematodes: the androdioecious (hermaphroditic) model Pristionchus pacificus and its dioecious sister species Pristionchus exspectatus. A chromosome-level genome assembly of P. exspectatus using single-molecule and Hi-C sequencing revealed a chromosome-wide rearrangement relative to P. pacificus. Strikingly, genomic characterization and cytogenetic studies including outgroup species Pristionchus occultus indicated two independent fusions involving the same chromosome, ChrIR, between these related species. Genetic linkage analysis indicated that these fusions altered the chromosome-wide pattern of recombination, resulting in large low-recombination regions that probably facilitated the coevolution between some of the ~14.8% of genes across the entire genomes. Quantitative trait locus analyses for hybrid sterility in all three sexes revealed that major quantitative trait loci mapped to the fused chromosome ChrIR. While abnormal chromosome segregations of the fused chromosome partially explain hybrid female sterility, hybrid-specific recombination that breaks linkage of genes in the low-recombination region was associated with hybrid male sterility. Thus, recent chromosome fusions repatterned recombination rate and drove reproductive isolation during Pristionchus speciation.

Heritable Cas9-induced nonhomologous recombination in C. elegans

Identification of the genetic basis of phenotypic variation within species remains challenging. In species with low recombination rates, such as Caenorhabditis elegans , genomic regions linked to a phenotype of interest by genetic mapping studies are often large, making it difficult to identify the specific genes and DNA sequence variants that underlie phenotypic differences. Here, we introduce a method that enables researchers to induce heritable targeted recombination in C. elegans with Cas9. We demonstrate that high rates of targeted nonhomologous recombination can be induced by Cas9 in a genomic region in which naturally occurring meiotic recombination events are exceedingly rare. We anticipate that Cas9-induced nonhomologous recombination (CINR) will greatly facilitate high-resolution genetic mapping in this species.

Cas9-induced nonhomologous recombination in C. elegans

Identification of the genetic basis of phenotypic variation within species remains challenging. In species with low recombination rates, such as Caenorhabditis elegans , genomic regions linked to a phenotype of interest by genetic mapping studies are often large, making it difficult to identify the specific genes and DNA sequence variants that underlie phenotypic differences. Here, we introduce a method that enables researchers to induce targeted recombination in C. elegans with Cas9. We demonstrate that high rates of targeted recombination can be induced by Cas9 in a genomic region in which naturally occurring recombination events are exceedingly rare. We anticipate that Cas9-induced nonhomologous recombination (CINR) will greatly facilitate high-resolution genetic mapping in this species.

Chromatin landscape, DSB levels, and cKU-70/80 contribute to patterning of meiotic DSB processing along chromosomes in C. elegans

Programmed DNA double-strand break (DSB) formation is essential for achieving accurate chromosome segregation during meiosis. DSB repair timing and template choice are tightly regulated. However, little is known about how DSB distribution and the choice of repair pathway are regulated along the length of chromosomes, which has direct effects on the recombination landscape and chromosome remodeling at late prophase I. Here, we use the spatiotemporal resolution of meiosis in the Caenorhabditis elegans germline along with genetic approaches to study distribution of DSB processing and its regulation. High-resolution imaging of computationally straightened chromosomes immunostained for the RAD-51 recombinase marking DSB repair sites reveals that the pattern of RAD-51 foci throughout pachytene resembles crossover distribution in wild type. Specifically, RAD-51 foci occur primarily along the gene-poor distal thirds of the chromosomes in both early and late pachytene, and on both the X and the autosomes. However, this biased off-center distribution can be abrogated by the formation of excess DSBs. Reduced condensin function, but not an increase in total physical axial length, results in a homogeneous distribution of RAD-51 foci, whereas regulation of H3K9 methylation is required for the enrichment of RAD-51 at off-center positions. Finally, the DSB recognition heterodimer cKU-70/80, but not the non-homologous end-joining canonical ligase LIG-4, contributes to the enriched off-center distribution of RAD-51 foci. Taken together, our data supports a model by which regulation of the chromatin landscape, DSB levels, and DSB detection by cKU-70/80 collaborate to promote DSB processing by homologous recombination at off-center regions of the chromosomes in C. elegans.

A role for worm cutl-24 in background- and parent-of-origin-dependent ER stress resistance

BACKGROUND

Organisms in the wild can acquire disease- and stress-resistance traits that outstrip the programs endogenous to humans. Finding the molecular basis of such natural resistance characters is a key goal of evolutionary genetics. Standard statistical-genetic methods toward this end can perform poorly in organismal systems that lack high rates of meiotic recombination, like Caenorhabditis worms.

RESULTS

Here we discovered unique ER stress resistance in a wild Kenyan C. elegans isolate, which in inter-strain crosses was passed by hermaphrodite mothers to hybrid offspring. We developed an unbiased version of the reciprocal hemizygosity test, RH-seq, to explore the genetics of this parent-of-origin-dependent phenotype. Among top-scoring gene candidates from a partial-coverage RH-seq screen, we focused on the neuronally-expressed, cuticlin-like gene cutl-24 for validation. In gene-disruption and controlled crossing experiments, we found that cutl-24 was required in Kenyan hermaphrodite mothers for ER stress tolerance in their inter-strain hybrid offspring; cutl-24 was also a contributor to the trait in purebred backgrounds.

CONCLUSIONS

These data establish the Kenyan strain allele of cutl-24 as a determinant of a natural stress-resistant state, and they set a precedent for the dissection of natural trait diversity in invertebrate animals without the need for a panel of meiotic recombinants.

Linkage mapping reveals loci that underlie differences in Caenorhabditis elegans growth

Growth rate and body size are complex traits that contribute to the fitness of organisms. The identification of loci that underlie differences in these traits provides insights into the genetic contributions to development. Leveraging Caenorhabditis elegans as a tractable metazoan model for quantitative genetics, we can identify genomic regions that underlie differences in growth. We measured postembryonic growth of the laboratory-adapted wild-type strain (N2) and a wild strain from Hawaii (CB4856) and found differences in body size. Using linkage mapping, we identified three distinct quantitative trait loci (QTL) on chromosomes IV, V, and X that are associated with variation in body growth. We further examined these growth-associated quantitative trait loci using chromosome substitution strains and near-isogenic lines and validated the chromosome X quantitative trait loci. In addition, we generated a list of candidate genes for the chromosome X quantitative trait loci. These genes could potentially contribute to differences in animal growth and should be evaluated in subsequent studies. Our work reveals the genetic architecture underlying animal growth variation and highlights the genetic complexity of growth in Caenorhabditis elegans natural populations.

Selection on modifiers of genetic architecture under migration load

Gene flow between populations adapting to differing local environmental conditions might be costly because individuals can disperse to habitats where their survival is low or because they can reproduce with locally maladapted individuals. The amount by which the mean relative population fitness is kept below one creates an opportunity for modifiers of the genetic architecture to spread due to selection. Prior work that separately considered modifiers changing dispersal, recombination rates, or altering dominance or epistasis, has typically focused on the direction of selection rather than its absolute magnitude. We here develop methods to determine the strength of selection on modifiers of the genetic architecture, including modifiers of the dispersal rate, in populations that have previously evolved local adaptation. We consider scenarios with up to five loci contributing to local adaptation and derive a new model for the deterministic spread of modifiers. We find that selection for modifiers of epistasis and dominance is stronger than selection for decreased recombination, and that selection for partial reductions in recombination are extremely weak, regardless of the number of loci contributing to local adaptation. The spread of modifiers that reduce dispersal depends on the number of loci, epistasis and extent of local adaptation in the ancestral population. We identify a novel effect, that modifiers of dominance are more strongly selected when they are unlinked to the locus that they modify. These findings help explain population differentiation and reproductive isolation and provide a benchmark to compare selection on modifiers under finite population sizes and demographic stochasticity.

When populations of a species are spread over different habitats the populations can adapt to their local conditions, provided dispersal between habitats is low enough. Natural selection allows the populations to maintain local adaptation, but dispersal and gene flow create a cost called the migration load. The migration load measures how much fitness is lost because of dispersal between different habitats, and also creates an opportunity for selection to act on the arrangement and interaction between genes that are involved in local adaptation. Modifier genes can spread in these linked populations and cause functional, local adaptation genes, to become more closely linked on a chromosome, or change the way that these genes are expressed so that the locally adapted gene copy becomes dominant. We modeled this process and found that selection on modifiers that create tighter linkage between locally adapted genes is generally weak, and modifiers that cause gene interactions are more strongly selected. Even after these gene interactions have begun to evolve, further selection for increased gene interaction is still strong. Our results show that populations are more likely to adapt to local conditions by evolving new gene interactions than by evolving tightly linked gene clusters.

A comparative transcriptomics and eQTL approach identifies SlWD40 as a tomato fruit ripening regulator

Although multiple vital genes with strong effects on the tomato (Solanum lycopersicum) ripening process have been identified via the positional cloning of ripening mutants and cloning of ripening-related transcription factors (TFs), recent studies suggest that it is unlikely that we have fully characterized the gene regulatory networks underpinning this process. Here, combining comparative transcriptomics and expression QTLs, we identified 16 candidate genes involved in tomato fruit ripening and validated them through virus-induced gene silencing analysis. To further confirm the accuracy of the approach, one potential ripening regulator, SlWD40 (WD-40 repeats), was chosen for in-depth analysis. Co-expression network analysis indicated that master regulators such as RIN (ripening inhibitor) and NOR (nonripening) as well as vital TFs including FUL1 (FRUITFUL1), SlNAC4 (NAM, ATAF1,2, and CUC2 4), and AP2a (Activating enhancer binding Protein 2 alpha) strongly co-expressed with SlWD40. Furthermore, SlWD40 overexpression and RNAi lines exhibited substantially accelerated and delayed ripening phenotypes compared with the wild type, respectively. Moreover, transcriptome analysis of these transgenics revealed that expression patterns of ethylene biosynthesis genes, phytoene synthase, pectate lyase, and branched chain amino transferase 2, in SlWD40-RNAi lines were similar to those of rin and nor fruits, which further demonstrated that SlWD40 may act as an important ripening regulator in conjunction with RIN and NOR. These results are discussed in the context of current models of ripening and in terms of the use of comparative genomics and transcriptomics as an effective route for isolating causal genes underlying differences in genotypes.

Genetic architecture and temporal analysis of Caenorhabditis briggsae hybrid developmental delay

Identifying the alleles that reduce hybrid fitness is a major goal in the study of speciation genetics. It is rare to identify systems in which hybrid incompatibilities with minor phenotypic effects are segregating in genetically diverse populations of the same biological species. Such traits do not themselves cause reproductive isolation but might initiate the process. In the nematode Caenorhabditis briggsae, a small percent of F2 generation hybrids between two natural populations suffer from developmental delay, in which adulthood is reached after approximately 33% more time than their wild-type siblings. Prior efforts to identify the genetic basis for this hybrid incompatibility assessed linkage using one or two genetic markers on chromosome III and suggested that delay is caused by a toxin-antidote element. Here, we have genotyped F2 hybrids using multiple chromosome III markers to refine the developmental delay locus. Also, to better define the developmental delay phenotype, we measured the development rate of 66 F2 hybrids and found that delay is not restricted to a particular larval developmental stage. Deviation of the developmental delay frequency from hypothetical expectations for a toxin-antidote element adds support to the assertion that the epistatic interaction is not fully penetrant. Our mapping and refinement of the delay phenotype motivates future efforts to study the genetic architecture of hybrid dysfunction between genetically distinct populations of one species by identifying the underlying loci.

Emerging mechanisms and roles of meiotic crossover repression at centromeres

Crossover events during recombination in meiosis are essential for generating genetic diversity as well as crucial to allow accurate chromosomal segregation between homologous chromosomes. Spatial control for the distribution of crossover events along the chromosomes is largely a tightly regulated process and involves many facets such as interference, repression as well as assurance, to make sure that not too many or too few crossovers are generated. Repression of crossover events at the centromeres is a highly conserved process across all species tested. Failure to inhibit such recombination events can result in chromosomal mis-segregation during meiosis resulting in aneuploid gametes that are responsible for infertility or developmental disorders such as Down's syndrome and other trisomies in humans. In the past few decades, studies to understand the molecular mechanisms behind this repression have shown the involvement of a multitude of factors ranging from the centromere-specific proteins such as the kinetochore to the flanking pericentric heterochromatin as well as DNA double-strand break repair pathways. In this chapter, we review the different mechanisms of pericentric repression mechanisms known till date as well as highlight the importance of understanding this regulation in the context of chromosomal segregation defects. We also discuss the clinical implications of dysregulation of this process, especially in human reproductive health and genetic diseases.

Evaluating the power and limitations of genome-wide association studies in Caenorhabditis elegans

Quantitative genetics in Caenorhabditis elegans seeks to identify naturally segregating genetic variants that underlie complex traits. Genome-wide association studies scan the genome for individual genetic variants that are significantly correlated with phenotypic variation in a population, or quantitative trait loci. Genome-wide association studies are a popular choice for quantitative genetic analyses because the quantitative trait loci that are discovered segregate in natural populations. Despite numerous successful mapping experiments, the empirical performance of genome-wide association study has not, to date, been formally evaluated in C. elegans. We developed an open-source genome-wide association study pipeline called NemaScan and used a simulation-based approach to provide benchmarks of mapping performance in collections of wild C. elegans strains. Simulated trait heritability and complexity determined the spectrum of quantitative trait loci detected by genome-wide association studies. Power to detect smaller-effect quantitative trait loci increased with the number of strains sampled from the C. elegans Natural Diversity Resource. Population structure was a major driver of variation in mapping performance, with populations shaped by recent selection exhibiting significantly lower false discovery rates than populations composed of more divergent strains. We also recapitulated previous genome-wide association studies of experimentally validated quantitative trait variants. Our simulation-based evaluation of performance provides the community with critical context to pursue quantitative genetic studies using the C. elegans Natural Diversity Resource to elucidate the genetic basis of complex traits in C. elegans natural populations.

Multiple Pristionchus pacificus genomes reveal distinct evolutionary dynamics between de novo candidates and duplicated genes

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.

TransposonUltimate: software for transposon classification, annotation and detection

Most genomes harbor a large number of transposons, and they play an important role in evolution and gene regulation. They are also of interest to clinicians as they are involved in several diseases, including cancer and neurodegeneration. Although several methods for transposon identification are available, they are often highly specialised towards specific tasks or classes of transposons, and they lack common standards such as a unified taxonomy scheme and output file format. We present TransposonUltimate, a powerful bundle of three modules for transposon classification, annotation, and detection of transposition events. TransposonUltimate comes as a Conda package under the GPL-3.0 licence, is well documented and it is easy to install through https://github.com/DerKevinRiehl/TransposonUltimate. We benchmark the classification module on the large TransposonDB covering 891,051 sequences to demonstrate that it outperforms the currently best existing solutions. The annotation and detection modules combine sixteen existing softwares, and we illustrate its use by annotating Caenorhabditis elegans, Rhizophagus irregularis and Oryza sativa subs. japonica genomes. Finally, we use the detection module to discover 29 554 transposition events in the genomes of 20 wild type strains of C. elegans. Databases, assemblies, annotations and further findings can be downloaded from (https://doi.org/10.5281/zenodo.5518085).

The impact of species-wide gene expression variation on Caenorhabditis elegans complex traits

Phenotypic variation in organism-level traits has been studied in Caenorhabditis elegans wild strains, but the impacts of differences in gene expression and the underlying regulatory mechanisms are largely unknown. Here, we use natural variation in gene expression to connect genetic variants to differences in organismal-level traits, including drug and toxicant responses. We perform transcriptomic analyses on 207 genetically distinct C. elegans wild strains to study natural regulatory variation of gene expression. Using this massive dataset, we perform genome-wide association mappings to investigate the genetic basis underlying gene expression variation and reveal complex genetic architectures. We find a large collection of hotspots enriched for expression quantitative trait loci across the genome. We further use mediation analysis to understand how gene expression variation could underlie organism-level phenotypic variation for a variety of complex traits. These results reveal the natural diversity in gene expression and possible regulatory mechanisms in this keystone model organism, highlighting the promise of using gene expression variation to understand how phenotypic diversity is generated.

Evolution of sexual systems, sex chromosomes and sex-linked gene transcription in flatworms and roundworms

Many species with separate male and female individuals (termed ‘gonochorism’ in animals) have sex-linked genome regions. Here, we investigate evolutionary changes when genome regions become completely sex-linked, by analyses of multiple species of flatworms (Platyhelminthes; among which schistosomes recently evolved gonochorism from ancestral hermaphroditism), and roundworms (Nematoda) which have undergone independent translocations of different autosomes. Although neither the evolution of gonochorism nor translocations fusing ancestrally autosomal regions to sex chromosomes causes inevitable loss of recombination, we document that formerly recombining regions show genomic signatures of recombination suppression in both taxa, and become strongly genetically degenerated, with a loss of most genes. Comparisons with hermaphroditic flatworm transcriptomes show masculinisation and some defeminisation in schistosome gonad gene expression. We also find evidence that evolution of sex-linkage in nematodes is accompanied by transcriptional changes and dosage compensation. Our analyses also identify sex-linked genes that could assist future research aimed at controlling some of these important parasites.

Transitions between hermaphroditic and separate sexes are relatively understudied in animals compared to pants. Here, Wang et al. reconstruct the evolution of separate sexes in the flatworms and complex changes of sex chromosomes in the roundworms.

Local adaptation and spatiotemporal patterns of genetic diversity revealed by repeated sampling of Caenorhabditis elegans across the Hawaiian Islands

The nematode Caenorhabditis elegans is among the most widely studied organisms, but relatively little is known about its natural ecology. Genetic diversity is low across much of the globe but high in the Hawaiian Islands and across the Pacific Rim. To characterize the niche and genetic diversity of C. elegans on the Hawaiian Islands and to explore how genetic diversity might be influenced by local adaptation, we repeatedly sampled nematodes over a three-year period, measured various environmental parameters at each sampling site, and whole-genome sequenced the C. elegans isolates that we identified. We found that the typical Hawaiian C. elegans niche comprises moderately moist native forests at high elevations (500-1,500 m) where ambient air temperatures are cool (15-20°C). Compared to other Caenorhabditis species found on the Hawaiian Islands (e.g., Caenorhabditis briggsae and Caenorhabditis tropicalis), we found that C. elegans were enriched in native habitats. We measured levels of genetic diversity and differentiation among Hawaiian C. elegans and found evidence of seven genetically distinct groups distributed across the islands. Then, we scanned these genomes for signatures of local adaptation and identified 18 distinct regions that overlap with hyper-divergent regions, which may be maintained by balancing selection and are enriched for genes related to environmental sensing, xenobiotic detoxification, and pathogen resistance. These results provide strong evidence of local adaptation among Hawaiian C. elegans and contribute to our understanding of the forces that shape genetic diversity on the most remote volcanic archipelago in the world.

Chromosome-level reference genomes for two strains of Caenorhabditis briggsae: an improved platform for comparative genomics

The publication of the Caenorhabditis briggsae reference genome in 2003 enabled the first comparative genomics studies between C. elegans and C. briggsae, shedding light on the evolution of genome content and structure in the Caenorhabditis genus. However, despite being widely used, the currently available C. briggsae reference genome is substantially less complete and structurally accurate than the C. elegans reference genome. Here, we used high-coverage Oxford Nanopore long-read and chromosome conformation capture data to generate chromosome-level reference genomes for two C. briggsae strains: QX1410, a new reference strain closely related to the laboratory AF16 strain, and VX34, a highly divergent strain isolated in China. We also sequenced 99 recombinant inbred lines (RILs) generated from reciprocal crosses between QX1410 and VX34 to create a recombination map and identify chromosomal domains. Additionally, we used both short- and long-read RNA sequencing (RNA-seq) data to generate high-quality gene annotations. By comparing these new reference genomes to the current reference, we reveal that hyper-divergent haplotypes cover large portions of the C. briggsae genome, similar to recent reports in C. elegans and C. tropicalis. We also show that the genomes of selfing Caenorhabditis species have undergone more rearrangement than their outcrossing relatives, which has biased previous estimates of rearrangement rate in Caenorhabditis. These new genomes provide a substantially improved platform for comparative genomics in Caenorhabditis and narrow the gap between the quality of genomic resources available for C. elegans and C. briggsae.

Nematode chromosomes

The nematode Caenorhabditis elegans has shed light on many aspects of eukaryotic biology, including genetics, development, cell biology, and genomics. A major factor in the success of C. elegans as a model organism has been the availability, since the late 1990s, of an essentially gap-free and well-annotated nuclear genome sequence, divided among 6 chromosomes. In this review, we discuss the structure, function, and biology of C. elegans chromosomes and then provide a general perspective on chromosome biology in other diverse nematode species. We highlight malleable chromosome features including centromeres, telomeres, and repetitive elements, as well as the remarkable process of programmed DNA elimination (historically described as chromatin diminution) that induces loss of portions of the genome in somatic cells of a handful of nematode species. An exciting future prospect is that nematode species may enable experimental approaches to study chromosome features and to test models of chromosome evolution. In the long term, fundamental insights regarding how speciation is integrated with chromosome biology may be revealed.

Post-insemination selection dominates pre-insemination selection in driving rapid evolution of male competitive ability

Sexual reproduction is a complex process that contributes to differences between the sexes and divergence between species. From a male’s perspective, sexual selection can optimize reproductive success by acting on the variance in mating success (pre-insemination selection) as well as the variance in fertilization success (post-insemination selection). The balance between pre- and post-insemination selection has not yet been investigated using a strong hypothesis-testing framework that directly quantifies the effects of post-insemination selection on the evolution of reproductive success. Here we use experimental evolution of a uniquely engineered genetic system that allows sperm production to be turned off and on in obligate male-female populations of Caenorhabditis elegans. We show that enhanced post-insemination competition increases the efficacy of selection and surpasses pre-insemination sexual selection in driving a polygenic response in male reproductive success. We find that after 10 selective events occurring over 30 generations post-insemination selection increased male reproductive success by an average of 5- to 7-fold. Contrary to expectation, enhanced pre-insemination competition hindered selection and slowed the rate of evolution. Furthermore, we found that post-insemination selection resulted in a strong polygenic response at the whole-genome level. Our results demonstrate that post-insemination sexual selection plays a critical role in the rapid optimization of male reproductive fitness. Therefore, explicit consideration should be given to post-insemination dynamics when considering the population effects of sexual selection.

Some of the most dramatic and diverse phenotypes observed in nature––such as head-butting in wild sheep and the elaborate tails of peacocks––are sexually dimorphic. These remarkable phenotypes are a result of sexual selection optimizing reproductive success in females and males independently. For males, total reproductive success is comprised of winning a mating event and then translating that mating event into a fertilization event. Therefore, to understand not only how male reproductive success is comprised, but also how it evolves, we must examine the interaction between pre- and post-insemination sexual selection. We combine environmentally-inducible control of sperm production within a highly reproducible factorial experimental evolution design to directly quantify the contribution of post-insemination selection to male reproductive evolution. We demonstrate that enhanced sperm competition increases the efficacy of selection and enhances the rate of male evolution. Alternatively, we show that enhanced pre-insemination competition slows the evolutionary rate. Using whole-genome approaches, we identify over 60 genes that contribute to male fertilization success. Brought together, our new approaches and results demonstrate that the unseen world of molecular interactions occurring during post-insemination are as fundamentally important as pre-mating factors.

Microbial protection favors parasite tolerance and alters host-parasite coevolutionary dynamics

Coevolution between hosts and parasites is a major driver of rapid evolutionary change¹ and diversification.^2,3 However, direct antagonistic interactions between hosts and parasites could be disrupted⁴ when host microbiota form a line of defense, a phenomenon widespread across animal and plant species.^5,6 By suppressing parasite infection, protective microbiota could reduce the need for host-based defenses and favor host support for microbiota colonization,⁶ raising the possibility that the microbiota can alter host-parasite coevolutionary patterns and processes.⁷ Here, using an experimental evolution approach, we co-passaged populations of nematode host (Caenorhabditis elegans) and parasites (Staphylococcus aureus) when hosts were colonized (or not) by protective bacteria (Enterococcus faecalis). We found that microbial protection during coevolution resulted in the evolution of host mortality tolerance—higher survival following parasite infection—and in parasites adapting to microbial defenses. Compared to unprotected host-parasite coevolution, the protected treatment was associated with reduced dominance of fluctuating selection dynamics in host populations. No differences in host recombination rate or genetic diversity were detected. Genomic divergence was observed between parasite populations coevolved in protected and unprotected hosts. These findings indicate that protective host microbiota can determine the evolution of host defense strategies and shape host-parasite coevolutionary dynamics.

•

Microbial protection resulted in the evolution of host mortality tolerance

•

Parasites adapted to counter microbial defenses within hosts

•

Protective microbes reduced fluctuating selection dynamics

•

Microbial protection did not impact host genetic diversity or recombination rates

Rapid Isolation of Wild Nematodes by Baermann Funnel

Beyond being robust experimental model organisms, Caenorhabditis elegans and its relatives are also real animals that live in nature. Studies of wild nematodes in their natural environments are valuable for understanding many aspects of biology, including the selective regimes in which distinctive genomic and phenotypic characters evolve, the genetic basis for complex trait variation, and the natural genetic diversity fundamental to all animal populations. This manuscript describes a simple and efficient method for extracting nematodes from their natural substrates, including rotting fruits, flowers, fungi, leaf litter, and soil. The Baermann funnel method, a classical nematology technique, selectively isolates active nematodes from their substrates. Because it recovers nearly all active worms from the sample, the Baermann funnel technique allows for the recovery of rare and slow-growing genotypes that co-occur with abundant and fast-growing genotypes, which might be missed in extraction methods that involve multiple generations of reproduction. The technique is also well suited to addressing metagenetic, population-genetic, and ecological questions. It captures the entire population in a sample simultaneously, allowing an unbiased view of the natural distribution of ages, sexes, and genotypes. The protocol allows for deployment at scale in the field, rapidly converting substrates into worm plates, and the authors have validated it through fieldwork on multiple continents.