Selection on a Subunit of the NURF Chromatin Remodeler Modifies Life History Traits in a Domesticated Strain of Caenorhabditis elegans

A toxin-antidote selfish element increases fitness of its host

Selfish genetic elements can promote their transmission at the expense of individual survival, creating conflict between the element and the rest of the genome. Recently, a large number of toxin-antidote (TA) post-segregation distorters have been identified in non-obligate outcrossing nematodes. Their origin and the evolutionary forces that keep them at intermediate population frequencies are poorly understood. Here, we study a TA element in Caenorhabditis elegans called zeel-1;peel-1. Two major haplotypes of this locus, with and without the selfish element, segregate in C. elegans. We evaluate the fitness consequences of the zeel-1;peel-1 element outside of its role in gene drive in non-outcrossing animals and demonstrate that loss of the toxin peel-1 decreased fitness of hermaphrodites and resulted in reductions in fecundity and body size. These findings suggest a biological role for peel-1 beyond toxin lethality. This work demonstrates that a TA element can provide a fitness benefit to its hosts either during their initial evolution or by being co-opted by the animals following their selfish spread. These findings guide our understanding on how TA elements can remain in a population where gene drive is minimized, helping resolve the mystery of prevalent TA elements in selfing animals.

Chromatin: the old and young of it

Aging affects nearly all aspects of our cells, from our DNA to our proteins to how our cells handle stress and communicate with each other. Age-related chromatin changes are of particular interest because chromatin can dynamically respond to the cellular and organismal environment, and many modifications at chromatin are reversible. Changes at chromatin occur during aging, and evidence from model organisms suggests that chromatin factors could play a role in modulating the aging process itself, as altering proteins that work at chromatin often affect the lifespan of yeast, worms, flies, and mice. The field of chromatin and aging is rapidly expanding, and high-resolution genomics tools make it possible to survey the chromatin environment or track chromatin factors implicated in longevity with precision that was not previously possible. In this review, we discuss the state of chromatin and aging research. We include examples from yeast, Drosophila, mice, and humans, but we particularly focus on the commonly used aging model, the worm Caenorhabditis elegans, in which there are many examples of chromatin factors that modulate longevity. We include evidence of both age-related changes to chromatin and evidence of specific chromatin factors linked to longevity in core histones, nuclear architecture, chromatin remodeling, and histone modifications.

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.

Natural genetic variation as a tool for discovery in Caenorhabditis nematodes

Over the last 20 years, studies of Caenorhabditis elegans natural diversity have demonstrated the power of quantitative genetic approaches to reveal the evolutionary, ecological, and genetic factors that shape traits. These studies complement the use of the laboratory-adapted strain N2 and enable additional discoveries not possible using only one genetic background. In this chapter, we describe how to perform quantitative genetic studies in Caenorhabditis, with an emphasis on C. elegans. These approaches use correlations between genotype and phenotype across populations of genetically diverse individuals to discover the genetic causes of phenotypic variation. We present methods that use linkage, near-isogenic lines, association, and bulk-segregant mapping, and we describe the advantages and disadvantages of each approach. The power of C. elegans quantitative genetic mapping is best shown in the ability to connect phenotypic differences to specific genes and variants. We will present methods to narrow genomic regions to candidate genes and then tests to identify the gene or variant involved in a quantitative trait. The same features that make C. elegans a preeminent experimental model animal contribute to its exceptional value as a tool to understand natural phenotypic variation.

The genetics of gene expression in a Caenorhabditis elegans multiparental recombinant inbred line population

Studying genetic variation of gene expression provides a powerful way to unravel the molecular components underlying complex traits. Expression quantitative trait locus (eQTL) studies have been performed in several different model species, yet most of these linkage studies have been based on the genetic segregation of two parental alleles. Recently, we developed a multiparental segregating population of 200 recombinant inbred lines (mpRILs) derived from four wild isolates (JU1511, JU1926, JU1931, and JU1941) in the nematode Caenorhabditis elegans. We used RNA-seq to investigate how multiple alleles affect gene expression in these mpRILs. We found 1789 genes differentially expressed between the parental lines. Transgression, expression beyond any of the parental lines in the mpRILs, was found for 7896 genes. For expression QTL mapping almost 9000 SNPs were available. By combining these SNPs and the RNA-seq profiles of the mpRILs, we detected almost 6800 eQTLs. Most trans-eQTLs (63%) co-locate in six newly identified trans-bands. The trans-eQTLs found in previous two-parental allele eQTL experiments and this study showed some overlap (17.5–46.8%), highlighting on the one hand that a large group of genes is affected by polymorphic regulators across populations and conditions, on the other hand, it shows that the mpRIL population allows identification of novel gene expression regulatory loci. Taken together, the analysis of our mpRIL population provides a more refined insight into C. elegans complex trait genetics and eQTLs in general, as well as a starting point to further test and develop advanced statistical models for detection of multiallelic eQTLs and systems genetics studying the genotype–phenotype relationship.

From QTL to gene: C. elegans facilitates discoveries of the genetic mechanisms underlying natural variation

Although many studies have examined quantitative trait variation across many species, only a small number of genes and thereby molecular mechanisms have been discovered. Without these data, we can only speculate about evolutionary processes that underlie trait variation. Here, we review how quantitative and molecular genetics in the nematode Caenorhabditis elegans led to the discovery and validation of 37 quantitative trait genes over the past 15 years. Using these data, we can start to make inferences about evolution from these quantitative trait genes, including the roles that coding versus noncoding variation, gene family expansion, common versus rare variants, pleiotropy, and epistasis play in trait variation across this species.

Stressful development: integrating endoderm development, stress, and longevity

In addition to performing digestion and nutrient absorption, the intestine serves as one of the first barriers to the external environment, crucial for protecting the host from environmental toxins, pathogenic invaders, and other stress inducers. The gene regulatory network (GRN) governing embryonic development of the endoderm and subsequent differentiation and maintenance of the intestine has been well-documented in C. elegans. A key regulatory input that initiates activation of the embryonic GRN for endoderm and mesoderm in this animal is the maternally provided SKN-1 transcription factor, an ortholog of the vertebrate Nrf1 and 2, which, like C. elegans SKN-1, perform conserved regulatory roles in mediating a variety of stress responses across metazoan phylogeny. Other key regulatory factors in early gut development also participate in stress response as well as in innate immunity and aging and longevity. In this review, we discuss the intersection between genetic nodes that mediate endoderm/intestine differentiation and regulation of stress and homeostasis. We also consider how direct signaling from the intestine to the germline, in some cases involving SKN-1, facilitates heritable epigenetic changes, allowing transmission of adaptive stress responses across multiple generations. These connections between regulation of endoderm/intestine development and stress response mechanisms suggest that varying selective pressure exerted on the stress response pathways may influence the architecture of the endoderm GRN, thereby leading to genetic and epigenetic variation in early embryonic GRN regulatory events.

Starvation Responses Throughout the Caenorhabditiselegans Life Cycle

Caenorhabditis elegans survives on ephemeral food sources in the wild, and the species has a variety of adaptive responses to starvation. These features of its life history make the worm a powerful model for studying developmental, behavioral, and metabolic starvation responses. Starvation resistance is fundamental to life in the wild, and it is relevant to aging and common diseases such as cancer and diabetes. Worms respond to acute starvation at different times in the life cycle by arresting development and altering gene expression and metabolism. They also anticipate starvation during early larval development, engaging an alternative developmental program resulting in dauer diapause. By arresting development, these responses postpone growth and reproduction until feeding resumes. A common set of signaling pathways mediates systemic regulation of development in each context but with important distinctions. Several aspects of behavior, including feeding, foraging, taxis, egg laying, sleep, and associative learning, are also affected by starvation. A variety of conserved signaling, gene regulatory, and metabolic mechanisms support adaptation to starvation. Early life starvation can have persistent effects on adults and their descendants. With its short generation time, C. elegans is an ideal model for studying maternal provisioning, transgenerational epigenetic inheritance, and developmental origins of adult health and disease in humans. This review provides a comprehensive overview of starvation responses throughout the C. elegans life cycle.

Dissecting the eQTL Micro-Architecture in Caenorhabditis elegans

The study of expression quantitative trait loci (eQTL) using natural variation in inbred populations has yielded detailed information about the transcriptional regulation of complex traits. Studies on eQTL using recombinant inbred lines (RILs) led to insights on cis and trans regulatory loci of transcript abundance. However, determining the underlying causal polymorphic genes or variants is difficult, but ultimately essential for the understanding of regulatory networks of complex traits. This requires insight into whether associated loci are single eQTL or a combination of closely linked eQTL, and how this QTL micro-architecture depends on the environment. We addressed these questions by testing for independent replication of previously mapped eQTL in Caenorhabditis elegans using new data from introgression lines (ILs). Both populations indicate that the overall heritability of gene expression, number, and position of eQTL differed among environments. Across environments we were able to replicate 70% of the cis- and 40% of the trans-eQTL using the ILs. Testing eight different simulation models, we suggest that additive effects explain up to 60-93% of RIL/IL heritability for all three environments. Closely linked eQTL explained up to 40% of RIL/IL heritability in the control environment whereas only 7% in the heat-stress and recovery environments. In conclusion, we show that reproducibility of eQTL was higher for cis vs. trans eQTL and that the environment affects the eQTL micro-architecture.

The Ancestral Caenorhabditis elegans Cuticle Suppresses rol-1

Genetic background commonly modifies the effects of mutations. We discovered that worms mutant for the canonical rol-1 gene, identified by Brenner in 1974, do not roll in the genetic background of the wild strain CB4856. Using linkage mapping, association analysis and gene editing, we determined that N2 carries an insertion in the collagen gene col-182 that acts as a recessive enhancer of rol-1 rolling. From population and comparative genomics, we infer the insertion is derived in N2 and related laboratory lines, likely arising during the domestication of Caenorhabditis elegans, and breaking a conserved protein. The ancestral version of col-182 also modifies the phenotypes of four other classical cuticle mutant alleles, and the effects of natural genetic variation on worm shape and locomotion. These results underscore the importance of genetic background and the serendipity of Brenner's choice of strain.

A spontaneous complex structural variant in rcan-1 increases exploratory behavior and laboratory fitness of Caenorhabditis elegans

Over long evolutionary timescales, major changes to the copy number, function, and genomic organization of genes occur, however, our understanding of the individual mutational events responsible for these changes is lacking. In this report, we study the genetic basis of adaptation of two strains of C. elegans to laboratory food sources using competition experiments on a panel of 89 recombinant inbred lines (RIL). Unexpectedly, we identified a single RIL with higher relative fitness than either of the parental strains. This strain also displayed a novel behavioral phenotype, resulting in higher propensity to explore bacterial lawns. Using bulk-segregant analysis and short-read resequencing of this RIL, we mapped the change in exploration behavior to a spontaneous, complex rearrangement of the rcan-1 gene that occurred during construction of the RIL panel. We resolved this rearrangement into five unique tandem inversion/duplications using Oxford Nanopore long-read sequencing. rcan-1 encodes an ortholog to human RCAN1/DSCR1 calcipressin gene, which has been implicated as a causal gene for Down syndrome. The genomic rearrangement in rcan-1 creates two complete and two truncated versions of the rcan-1 coding region, with a variety of modified 5’ and 3’ non-coding regions. While most copy-number variations (CNVs) are thought to act by increasing expression of duplicated genes, these changes to rcan-1 ultimately result in the reduction of its whole-body expression due to changes in the upstream regions. By backcrossing this rearrangement into a common genetic background to create a near isogenic line (NIL), we demonstrate that both the competitive advantage and exploration behavioral changes are linked to this complex genetic variant. This NIL strain does not phenocopy a strain containing an rcan-1 loss-of-function allele, which suggests that the residual expression of rcan-1 is necessary for its fitness effects. Our results demonstrate how colonization of new environments, such as those encountered in the laboratory, can create evolutionary pressure to modify gene function. This evolutionary mismatch can be resolved by an unexpectedly complex genetic change that simultaneously duplicates and diversifies a gene into two uniquely regulated genes. Our work shows how complex rearrangements can act to modify gene expression in ways besides increased gene dosage.

Evolution acts on genetic variants that modify phenotypes that increase the likelihood of staying alive and passing on these genetic changes to subsequent generations (i.e. fitness). There is general interest in understanding the types of genetic variants that can increase fitness in specific environments. One route that fitness can be increased is through changes in behavior, such as finding new food sources. Here, we identify a spontaneous genetic change that increases exploration behavior and fitness of animals in laboratory environments. Interestingly, this genetic change is not a simple genetic change that deletes or changes the sequence of a protein product, but rather a complex structural variant that simultaneously duplicates the rcan-1 gene and also modifies its expression in a number of tissues. Our work demonstrates how a complex structural change can duplicate a gene, modify the DNA control regions that determine its cellular sites of action, and confer a fitness advantage that could lead to its spread in a population.

WormQTL2: an interactive platform for systems genetics in Caenorhabditis elegans

Quantitative genetics provides the tools for linking polymorphic loci to trait variation. Linkage analysis of gene expression is an established and widely applied method, leading to the identification of expression quantitative trait loci (eQTLs). (e)QTL detection facilitates the identification and understanding of the underlying molecular components and pathways, yet (e)QTL data access and mining often is a bottleneck. Here, we present WormQTL2, a database and platform for comparative investigations and meta-analyses of published (e)QTL data sets in the model nematode worm C. elegans. WormQTL2 integrates six eQTL studies spanning 11 conditions as well as over 1000 traits from 32 studies and allows experimental results to be compared, reused and extended upon to guide further experiments and conduct systems-genetic analyses. For example, one can easily screen a locus for specific cis-eQTLs that could be linked to variation in other traits, detect gene-by-environment interactions by comparing eQTLs under different conditions, or find correlations between QTL profiles of classical traits and gene expression. WormQTL2 makes data on natural variation in C. elegans and the identified QTLs interactively accessible, allowing studies beyond the original publications. Database URL: www.bioinformatics.nl/WormQTL2/.

A fln-2 mutation affects lethal pathology and lifespan in C. elegans

Differences in genetic background in model organisms can have complex effects on phenotypes of interest. We previously reported a difference in hermaphrodite lifespan between two wild-type lines widely used by C. elegans researchers (N2 hermaphrodite and male stocks). Here, using pathology-based approaches and genome sequencing, we identify the cause of this difference as a nonsense mutation in the filamin gene fln-2 in the male stock, which reduces early mortality caused by pharyngeal infection. We show how fln-2 variation explains previous discrepancies involving effects of sir-2.1 (sirtuin deacetylase) on ageing, and show that in a fln-2(+) background, sir-2.1 over-expression causes an FUDR (DNA synthesis inhibitor)-dependent reduction in pharyngeal infection and increase in lifespan. In addition we show how fln-2 variation confounds effects on lifespan of daf-2 (insulin/IGF-1 signalling), daf-12 (steroid hormone signalling), and eat-2 (putative dietary restriction). These findings underscore the importance of identifying and controlling genetic background variation.

Natural Variation and Genetic Determinants of Caenorhabditis elegans Sperm Size

The diversity in sperm shape and size represents a powerful paradigm to understand how selection drives the evolutionary diversification of cell morphology. Experimental work on the sperm biology of the male-hermaphrodite nematode Caenorhabditis elegans has elucidated diverse factors important for sperm fertilization success, including the competitive superiority of larger sperm. Yet despite extensive research, the molecular mechanisms regulating C. elegans sperm size and the genetic basis underlying natural variation in sperm size remain unknown. To address these questions, we quantified male sperm size variation of a worldwide panel of 97 genetically distinct C. elegans strains, allowing us to uncover significant genetic variation in male sperm size. Aiming to characterize the molecular genetic basis of C. elegans male sperm size variation using a genome-wide association study, we did not detect any significant quantitative trait loci. We therefore focused on the genetic analysis of pronounced sperm size differences observed between recently diverged laboratory strains (N2 vs. LSJ1/2). Using mutants and quantitative complementation tests, we demonstrate that variation in the gene nurf-1 underlies the evolution of small sperm in the LSJ lineage. Given the previous discovery that this same nurf-1 variation was central for hermaphrodite laboratory adaptation, the evolution of reduced male sperm size in LSJ strains likely reflects a pleiotropic consequence. Together, our results provide a comprehensive quantification of natural variation in C. elegans sperm size and first insights into the genetic determinants of Caenorhabditis sperm size, pointing at an involvement of the NURF chromatin remodeling complex.

Evolution of Yin and Yang isoforms of a chromatin remodeling subunit precedes the creation of two genes

Genes can encode multiple isoforms, broadening their functions and providing a molecular substrate to evolve phenotypic diversity. Evolution of isoform function is a potential route to adapt to new environments. Here we show that de novo, beneficial alleles in the nurf-1 gene became fixed in two laboratory lineages of C. elegans after isolation from the wild in 1951, before methods of cryopreservation were developed. nurf-1 encodes an ortholog of BPTF, a large (>300 kD) multidomain subunit of the NURF chromatin remodeling complex. Using CRISPR-Cas9 genome editing and transgenic rescue, we demonstrate that in C. elegans, nurf-1 has split into two, largely non-overlapping isoforms (NURF-1.D and NURF-1.B, which we call Yin and Yang, respectively) that share only two of 26 exons. Both isoforms are essential for normal gametogenesis but have opposite effects on male/female gamete differentiation. Reproduction in hermaphrodites, which involves production of both sperm and oocytes, requires a balance of these opposing Yin and Yang isoforms. Transgenic rescue and genetic position of the fixed mutations suggest that different isoforms are modified in each laboratory strain. In a related clade of Caenorhabditis nematodes, the shared exons have duplicated, resulting in the split of the Yin and Yang isoforms into separate genes, each containing approximately 200 amino acids of duplicated sequence that has undergone accelerated protein evolution following the duplication. Associated with this duplication event is the loss of two additional nurf-1 transcripts, including the long-form transcript and a newly identified, highly expressed transcript encoded by the duplicated exons. We propose these lost transcripts are non-functional side products necessary to transcribe the Yin and Yang transcripts in the same cells. Our work demonstrates how gene sharing, through the production of multiple isoforms, can precede the creation of new, independent genes.

Natural variation in C. elegans arsenic toxicity is explained by differences in branched chain amino acid metabolism

We find that variation in the dbt-1 gene underlies natural differences in Caenorhabditis elegans responses to the toxin arsenic. This gene encodes the E2 subunit of the branched-chain α-keto acid dehydrogenase (BCKDH) complex, a core component of branched-chain amino acid (BCAA) metabolism. We causally linked a non-synonymous variant in the conserved lipoyl domain of DBT-1 to differential arsenic responses. Using targeted metabolomics and chemical supplementation, we demonstrate that differences in responses to arsenic are caused by variation in iso-branched chain fatty acids. Additionally, we show that levels of branched chain fatty acids in human cells are perturbed by arsenic treatment. This finding has broad implications for arsenic toxicity and for arsenic-focused chemotherapeutics across human populations. Our study implicates the BCKDH complex and BCAA metabolism in arsenic responses, demonstrating the power of C. elegans natural genetic diversity to identify novel mechanisms by which environmental toxins affect organismal physiology.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

A multi-parent recombinant inbred line population of C. elegans allows identification of novel QTLs for complex life history traits

BACKGROUND

The nematode Caenorhabditis elegans has been extensively used to explore the relationships between complex traits, genotypes, and environments. Complex traits can vary across different genotypes of a species, and the genetic regulators of trait variation can be mapped on the genome using quantitative trait locus (QTL) analysis of recombinant inbred lines (RILs) derived from genetically and phenotypically divergent parents. Most RILs have been derived from crossing two parents from globally distant locations. However, the genetic diversity between local C. elegans populations can be as diverse as between global populations and could thus provide means of identifying genetic variation associated with complex traits relevant on a broader scale.

RESULTS

To investigate the effect of local genetic variation on heritable traits, we developed a new RIL population derived from 4 parental wild isolates collected from 2 closely located sites in France: Orsay and Santeuil. We crossed these 4 genetically diverse parental isolates to generate a population of 200 multi-parental RILs and used RNA-seq to obtain sequence polymorphisms identifying almost 9000 SNPs variable between the 4 genotypes with an average spacing of 11 kb, doubling the mapping resolution relative to currently available RIL panels for many loci. The SNPs were used to construct a genetic map to facilitate QTL analysis. We measured life history traits such as lifespan, stress resistance, developmental speed, and population growth in different environments, and found substantial variation for most traits. We detected multiple QTLs for most traits, including novel QTLs not found in previous QTL analysis, including those for lifespan and pathogen responses. This shows that recombining genetic variation across C. elegans populations that are in geographical close proximity provides ample variation for QTL mapping.

CONCLUSION

Taken together, we show that using more parents than the classical two parental genotypes to construct a RIL population facilitates the detection of QTLs and that the use of wild isolates facilitates the detection of QTLs. The use of multi-parent RIL populations can further enhance our understanding of local adaptation and life history trade-offs.

Analysis of Epistasis in Natural Traits Using Model Organisms

The ability to detect and understand epistasis in natural populations is important for understanding how biological traits are influenced by genetic variation. However, identification and characterization of epistasis in natural populations remains difficult due to statistical issues that arise as a result of multiple comparisons, and the fact that most genetic variants segregate at low allele frequencies. In this review, we discuss how model organisms may be used to manipulate genotypic combinations to power the detection of epistasis as well as test interactions between specific genes. Findings from a number of species indicate that statistical epistasis is pervasive between natural genetic variants. However, the properties of experimental systems that enable analysis of epistasis also constrain extrapolation of these results back into natural populations.

Extreme allelic heterogeneity at a Caenorhabditis elegans beta-tubulin locus explains natural resistance to benzimidazoles

Benzimidazoles (BZ) are essential components of the limited chemotherapeutic arsenal available to control the global burden of parasitic nematodes. The emerging threat of BZ resistance among multiple nematode species necessitates the development of novel strategies to identify genetic and molecular mechanisms underlying this resistance. All detection of parasitic helminth resistance to BZ is focused on the genotyping of three variant sites in the orthologs of the β-tubulin gene found to confer resistance in the free-living nematode Caenorhabditis elegans. Because of the limitations of laboratory and field experiments in parasitic nematodes, it is difficult to look beyond these three sites to identify additional mechanisms that might contribute to BZ resistance in the field. Here, we took an unbiased genome-wide mapping approach in the free-living nematode species C. elegans to identify the genetic underpinnings of natural resistance to the commonly used BZ, albendazole (ABZ). We found a wide range of natural variation in ABZ resistance in natural C. elegans populations. In agreement with known mechanisms of BZ resistance in parasites, we found that a majority of the variation in ABZ resistance among wild C. elegans strains is caused by variation in the β-tubulin gene ben-1. This result shows empirically that resistance to ABZ naturally exists and segregates within the C. elegans population, suggesting that selection in natural niches could enrich for resistant alleles. We identified 25 distinct ben-1 alleles that are segregating at low frequencies within the C. elegans population, including many novel molecular variants. Population genetic analyses indicate that ben-1 variation arose multiple times during the evolutionary history of C. elegans and provide evidence that these alleles likely occurred recently because of local selective pressures. Additionally, we find purifying selection at all five β-tubulin genes, despite predicted loss-of-function variants in ben-1, indicating that BZ resistance in natural niches is a stronger selective pressure than loss of one β-tubulin gene. Furthermore, we used genome-editing to show that the most common parasitic nematode β-tubulin allele that confers BZ resistance, F200Y, confers resistance in C. elegans. Importantly, we identified a novel genomic region that is correlated with ABZ resistance in the C. elegans population but independent of ben-1 and the other β-tubulin loci, suggesting that there are multiple mechanisms underlying BZ resistance. Taken together, our results establish a population-level resource of nematode natural diversity as an important model for the study of mechanisms that give rise to BZ resistance.

Changes to social feeding behaviors are not sufficient for fitness gains of the Caenorhabditis elegans N2 reference strain

The standard reference Caenorhabditis elegans strain, N2, has evolved marked behavioral changes in social feeding behavior since its isolation from the wild. We show that the causal, laboratory-derived mutations in two genes, npr-1 and glb-5, confer large fitness advantages in standard laboratory conditions. Using environmental manipulations that suppress social/solitary behavior differences, we show the fitness advantages of the derived alleles remained unchanged, suggesting selection on these alleles acted through pleiotropic traits. Transcriptomics, developmental timing, and food consumption assays showed that N2 animals mature faster, produce more sperm, and consume more food than a strain containing ancestral alleles of these genes regardless of behavioral strategies. Our data suggest that the pleiotropic effects of glb-5 and npr-1 are a consequence of changes to O₂ -sensing neurons that regulate both aerotaxis and energy homeostasis. Our results demonstrate how pleiotropy can lead to profound behavioral changes in a popular laboratory model.

Why do humans walk on two feet? And what makes us smarter than our ape ancestors? The answers to these questions, and countless others about the particular traits of any number of species, is often said to be natural selection – a process where genes that ensure the survival of a species are favored of others. But it is not always the answer. Other evolutionary forces, such as random changes to the frequency of certain gene variants, restrictions on the development of a certain trait and pleiotropy (where one gene influences other, seemingly unrelated traits) can also cause differences between species.

Designing experiments to test whether a trait difference is due to natural selection or other factors is notoriously difficult. However, the humble nematode worm, Caenorhabditis elegans, has proven to be particularly useful in this respect. One subtype or strain of C. elegans with certain changes to its genes is used internationally as a ‘reference strain’, to ensure results between labs are comparable. This strain, N2, has been bred in the laboratory for hundreds of generations, isolated from its wild counterparts. N2 shows several differences in behavior from the wildtype, including its feeding habits. Wild C. elegans tend to feed together socially, whereas N2 prefers to feed alone.

In 1998 and 2009, researchers – including some involved in the current study – have identified the genetic modifications responsible for this change in behavior. Now, Zhao et al. set out to determine whether this was due to natural selection, and if so, was there a benefit to solitary feeding in laboratory conditions that was driving this genetic change? Zhao et al. found that the genetic changes in the N2 strain gave the worms a considerable advantage in the artificial environment. However, experiments to modify the conditions the animals grew in revealed that the solitary feeding habits were not necessary for the fitness advantage. In other words, the changes in feeding habits were a symptom of the genetic changes that gave N2 a selective advantage, but they were not the cause. In other words, the changes in feeding behavior were not a result of natural selection, but rather of pleiotropy.

The findings highlight that not every change in a trait is down to natural selection and must therefore be put to the test. With declining costs of DNA sequencing, researchers can now easily identify genes and regions of DNA that are likely to be under selection. However, they must be careful before leaping to the conclusion that behavioral differences linked to genetic changes are adaptive. In addition, the findings show that the laboratories relying on N2 as a model organism should be aware that the strain has evolved fundamental differences in its brain connections compared with the wildtype.

The Chromatin Remodeler BPTF Activates a Stemness Gene-Expression Program Essential for the Maintenance of Adult Hematopoietic Stem Cells

Self-renewal and differentiation of adult stem cells are tightly regulated partly through configuration of chromatin structure by chromatin remodelers. Using knockout mice, we here demonstrate that bromodomain PHD finger transcription factor (BPTF), a component of the nucleosome remodeling factor (NURF) chromatin-remodeling complex, is essential for maintaining the population size of hematopoietic stem/progenitor cells (HSPCs), including long-term hematopoietic stem cells (HSCs). Bptf-deficient HSCs are defective in reconstituted hematopoiesis, and hematopoietic-specific knockout of Bptf caused profound defects including bone marrow failure and anemia. Genome-wide transcriptome profiling revealed that BPTF loss caused downregulation of HSC-specific gene-expression programs, which contain several master transcription factors (Meis1, Pbx1, Mn1, and Lmo2) required for HSC maintenance and self-renewal. Furthermore, we show that BPTF potentiates the chromatin accessibility of key HSC “stemness” genes. These results demonstrate an essential requirement of the chromatin remodeler BPTF and NURF for activation of “stemness” gene-expression programs and proper function of adult HSCs.

•

The chromatin remodeler gene Bptf shows preferential expression in primitive HSPCs

•

Bptf ensures maintenance and functionality of HSCs in a cell-autonomous manner

•

Bptf is required for the activation of HSC “stemness” genes including master TFs

•

Hematopoietic-specific loss of Bptf results in bone marrow failure and anemia

Natural Variation in the Distribution and Abundance of Transposable Elements Across the Caenorhabditis elegans Species

Transposons are mobile DNA elements that generate both adaptive and deleterious phenotypic variation thereby driving genome evolution. For these reasons, genomes have mechanisms to regulate transposable element (TE) activity. Approximately 12–16% of the Caenorhabditis elegans genome is composed of TEs, of which the majority are likely inactive. However, most studies of TE activity have been conducted in the laboratory strain N2, which limits our knowledge of the effects of these mobile elements across natural populations. We analyzed the distribution and abundance of TEs in 208 wild C. elegans strains to better understand how transposons contribute to variation in natural populations. We identified 3,397 TEs as compared with the reference strain, of which 2,771 are novel insertions and 241 are TEs that have been excised in at least one wild strain. Likely because of their hypothesized deleterious effects, we find that TEs are found at low allele frequencies throughout the population, and we predict functional effects of TE insertions. The abundances of TEs reflect their activities, and these data allowed us to perform both genome-wide association mappings and rare variant correlations to reveal several candidate genes that impact TE regulation, including small regulatory piwi-interacting RNAs and chromatin factors. Because TE variation in natural populations could underlie phenotypic variation for organismal and behavioral traits, the transposons that we identified and their regulatory mechanisms can be used in future studies to explore the genomics of complex traits and evolutionary changes.

Experimental Evolution with Caenorhabditis Nematodes

The hermaphroditic nematode Caenorhabditis elegans has been one of the primary model systems in biology since the 1970s, but only within the last two decades has this nematode also become a useful model for experimental evolution. Here, we outline the goals and major foci of experimental evolution with C. elegans and related species, such as C. briggsae and C. remanei, by discussing the principles of experimental design, and highlighting the strengths and limitations of Caenorhabditis as model systems. We then review three exemplars of Caenorhabditis experimental evolution studies, underlining representative evolution experiments that have addressed the: (1) maintenance of genetic variation; (2) role of natural selection during transitions from outcrossing to selfing, as well as the maintenance of mixed breeding modes during evolution; and (3) evolution of phenotypic plasticity and its role in adaptation to variable environments, including host-pathogen coevolution. We conclude by suggesting some future directions for which experimental evolution with Caenorhabditis would be particularly informative.

Microfluidic platform with spatiotemporally controlled micro-environment for studying long-term C. elegans developmental arrests

Animals' long-term survival is dependent on their ability to sense, filter and respond to their environment at multiple timescales. For example, during development, animals integrate environmental information, which can then modulate adult behavior and developmental trajectory. The neural and molecular mechanisms that underlie these changes are poorly understood. C. elegans is a powerful model organism to study such mechanisms; however, conventional plate-based culturing techniques are limited in their ability to consistently control and modulate an animal's environmental conditions. To address this need, we developed a microfluidics-based experimental platform capable of long-term culture of populations of developing C. elegans covering the L1 larval stage to adulthood, while achieving spatial consistency and temporal control of their environment. To prevent bacterial accumulation and maintain optimal flow characteristics and nutrient consistency over the operational period of over one hundred and fifty hours, several features of the microfluidic system and the peripheral equipment were optimized. By manipulating food and pheromone exposure over several days, we were able to demonstrate environmental-dependent changes to growth rate and entry to dauer, an alternative developmental state. We envision this system to be useful in studying the mechanisms underlying long timescale changes to behavior and development in response to environmental changes.

Modeling of a negative feedback mechanism explains antagonistic pleiotropy in reproduction in domesticated Caenorhabditis elegans strains

Most biological traits and common diseases have a strong but complex genetic basis, controlled by large numbers of genetic variants with small contributions to a trait or disease risk. The effect-size of most genetic variants is not absolute and is instead dependent upon multiple factors such as the age and genetic background of an organism. In order to understand the mechanistic basis of these changes, we characterized heritable trait differences between two domesticated strains of C. elegans. We previously identified a major effect locus, caused in part by a mutation in a component of the NURF chromatin remodeling complex, that regulates reproductive output in an age-dependent manner. The effect-size of this locus changes from positive to negative over the course of an animal’s reproductive lifespan. Here, we use a previously published macroscale model of the egg-laying rate in C. elegans to show that time-dependent effect-size is explained by an unequal use of sperm combined with negative feedback between sperm and ovulation rate. We validate key predictions of this model with controlled mating experiments and quantification of oogenesis and sperm use. Incorporation of this model into QTL mapping allows us to identify and partition new QTLs into specific aspects of the egg-laying process. Finally, we show how epistasis between two genetic variants is predicted by this modeling as a consequence of the unequal use of sperm. This work demonstrates how modeling of multicellular communication systems can improve our ability to predict and understand the role of genetic variation on a complex phenotype. Negative autoregulatory feedback loops, common in transcriptional regulation, could play an important role in modifying genetic architecture in other traits.

Complex traits are influenced by the individual effects of genetic variants in addition to the interactions of the variants with the environment, age, and each other. While complex genetic architectures are ubiquitous in natural traits, little is known about the causal mechanisms that create their complex genetic architectures. Here we identify an example of age-dependent genetic architecture controlling the rate and timing of reproduction in the hermaphroditic nematode C. elegans. We use computational modeling to demonstrate how age-dependent genetic architecture can arise as a consequence of two factors: hormonal feedback on oocytes mediated by major sperm protein (MSP) released by sperm stored in the spermatheca and life history differences in sperm use caused by genetic variants. Our work also suggests how antagonistic pleiotropy can emerge from multicellular feedback systems.