Hanh Witte

The Pristionchus pacificus genome provides a unique perspective on nematode lifestyle and parasitism

Here we present a draft genome sequence of the nematode Pristionchus pacificus, a species that is associated with beetles and is used as a model system in evolutionary biology. With 169 Mb and 23,500 predicted protein-coding genes, the P. pacificus genome is larger than those of Caenorhabditis elegans and the human parasite Brugia malayi. Compared to C. elegans, the P. pacificus genome has more genes encoding cytochrome P450 enzymes, glucosyltransferases, sulfotransferases and ABC transporters, many of which were experimentally validated. The P. pacificus genome contains genes encoding cellulase and diapausin, and cellulase activity is found in P. pacificus secretions, indicating that cellulases can be found in nematodes beyond plant parasites. The relatively higher number of detoxification and degradation enzymes in P. pacificus is consistent with its necromenic lifestyle and might represent a preadaptation for parasitism. Thus, comparative genomics analysis of three ecologically distinct nematodes offers a unique opportunity to investigate the association between genome structure and lifestyle.

Multi locus analysis of Pristionchus pacificus on La Reunion Island reveals an evolutionary history shaped by multiple introductions, constrained dispersal events and rare out-crossing

Pristionchus pacificus, recently established as a model organism in evolutionary biology, is a cosmopolitan nematode that has a necromenic association with scarab beetles. The diverse array of host beetle species and habitat types occupied by P. pacificus make it a good model for investigating local adaptation to novel environments. Presence of P. pacificus on La Réunion Island, a young volcanic island with a dynamic geological history and a wide variety of ecozones, facilitates such investigation in an island biogeographic setting. Microsatellite data from 20 markers and 223 strains and mitochondrial sequence data from 272 strains reveal rich genetic diversity among La Réunion P. pacificus isolates, shaped by differentially timed introductions from diverse sources and in association with different beetle species. Distinctions between volcanic zones and between arid western and wet eastern climatic zones have likely limited westward dispersal of recently colonized lineages and maintained a genetic distinction between eastern and western clades. The highly selfing lifestyle of P. pacificus contributes to the strong fine‐scale population structure detected, with each beetle host harbouring strongly differentiated assemblages of strains. Periodic out‐crossing generates admixture between genetically diverse lineages, creating a diverse array of allelic combinations likely to increase the evolutionary potential of the species and facilitate adaptation to local environments and beetle hosts.

Microevolutionary analysis of the nematode genus Pristionchus suggests a recent evolution of redundant developmental mechanisms during vulva formation.

SUMMARY To identify the mechanisms by which molecular variation is introduced into developmental systems, microevolutionary approaches to evolutionary developmental biology have to be taken. Here, we describe the molecular and developmental characterization of laboratory strains of the nematode genus Pristionchus, which lays a foundation for a microevolutionary analysis of vulva development. We describe 13 laboratory strains of the Pristionchus genus that are derived from natural isolates from around the world. Mating experiments and ITS sequence analysis indicated that these 13 strains represent four different species: the gonochoristic species P. lheritieri and three hermaphroditic species, P. pacificus, P. maupasi, and an as yet undescribed species Pristionchus sp., respectively. P. pacificus is represented by five different strains isolated from California, Washington, Hawaii, Ontario, and Poland. Developmental differences during vulva formation are observed between strains from different species but also between strains of P. pacificus, like the strains from California and Poland. In particular, redundant developmental mechanisms present during vulva formation in P. pacificus var. California are absent in other strains. Amplified restriction fragment length polymorphism (AFLP) analyses of the P. pacificus strains revealed that the American strains are highly polymorphic. In contrast, the developmentally distinct strain from Poland is identical to the Californian strain, suggesting that the developmental differences rely on a small number of changes in developmental control genes rather than the accumulation of changes at multiple loci.

Characterization of genetic diversity in the nematode pristionchus pacificus from population-scale resequencing data

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

Natural variation in Pristionchus pacificus insect pheromone attraction involves the protein kinase EGL-4

The geographical mosaic theory of coevolution predicts that different local species interactions will shape population traits, but little is known about the molecular factors involved in mediating the specificity of these interactions. Pristionchus nematodes associate with different scarab beetles around the world, with Pristionchus pacificus isolated primarily from the oriental beetle in Japan. In particular, the constituent populations of P. pacificus represent a rare opportunity to study multiple specialized interactions and the mechanisms that influence population traits at the genetic level. We identified a component of the cGMP signaling pathway to be involved in the natural variation for sensing the insect pheromone ETDA, using targeted introgression lines, exogenous cGMP treatment, and a null egl-4 allele. Our data strongly implicate egl-4 as one of several loci involved in behavioral variation in P. pacificus populations. That EGL-4 homologs have been independently implicated for behavioral variations in other invertebrate models suggests that EGL-4 may act as a modulator for interspecies behavioral repertoires across large phylogenetic distances.

Mutation Rates and Intraspecific Divergence of the Mitochondrial Genome of Pristionchus pacificus

Evolutionary reconstruction of the natural history of an organism ultimately requires knowledge about the development, population genetics, ecology, and phylogeny of the species. Such investigations would benefit from studies of mutational processes because mutations are the source of natural variation. The nematode Pristionchus pacificus has been developed as a model organism in evolutionary biology by comparing its development with Caenorhabditis elegans. Pristionchus pacificus and related species are associated with scarab beetles, and their ecology and phylogeny are well known. More than 200 P. pacificus isolates from all over the world are available for this cosmopolitan species. We generated mutation accumulation (MA) lines in P. pacificus to study spontaneous mutation rates in the mitochondrial genome and compared mutation rate estimates with natural variation between nine representative isolates of the species. The P. pacificus mitochondrial genome is 15,955 bp in length and is typical for nematodes. Pristionchus pacificus has all known mitochondrial genes and contains an unusual suppressor transfer RNA (tRNA) for the codon UAA. This has most likely influenced the spectrum of observable mutations because 6 of 12 mutations found in the 82 MA lines analyzed are nonsense mutations that can be suppressed by the suppressor tRNA. The overall mutation rate in P. pacificus is 7.6 × 10⁻⁸ per site per generation and is less than one order of magnitude different from estimates in C. elegans and Drosophila. Using this mutation rate estimate in a comparison of the mitochondrial genome of nine P. pacificus isolates, we calculate the minimum time to the most recent common ancestor at 10⁵-10⁶ generations. The combination of mutation rate analysis with intraspecific divergence provides a powerful tool for the reconstruction of the natural history of P. pacificus, and we discuss the ecological implication of these findings.

The importance of being regular: Caenorhabditis elegans and Pristionchus pacificus defecation mutants are hypersusceptible to bacterial pathogens

Bacterial pathogens have shaped the evolution and survival of organisms throughout history, but little is known about the evolution of virulence mechanisms and the counteracting defence strategies of host species. The nematode model organisms, Caenorhabditis elegans and Pristionchus pacificus, feed on a wealth of bacteria in their natural soil environment, some of which can cause mortality. Previously, we have shown that these nematodes differ in their susceptibility to a range of human and insect pathogenic bacteria, with P. pacificus showing extreme resistance compared with C. elegans. Here, we isolated 400 strains of Bacillus from soil samples and fed their spores to both nematodes. Spores of six Bacillus strains were found to kill C. elegans but not P. pacificus. While the majority of Bacillus strains are benign to nematodes, observed pathogenicity is restricted to either the spore or the vegetative stage. We used the rapid C. elegans killer strain (Bacillus sp. 142) to conduct a screen for hypersusceptible P. pacificus mutants. Two P. pacificus mutants with severe muscle defects and an extended defecation cycle that die rapidly on Bacillus spores were isolated. These genes were identified to be homologous to C. elegans, unc-22 and unc-13. To test whether a similar relationship between defecation and bacterial pathogenesis exists in C. elegans, we used five known defecation mutants. Quantification of the defecation cycle in mutants also revealed a severe effect on survival in C. elegans. Thus, intestinal peristalsis is critical to nematode health and contributes significantly to survival when fed Gram-positive bacteria.

Chromatin remodelling and antisense-mediated up-regulation of the developmental switch gene eud-1 control predatory feeding plasticity

Phenotypic plasticity has been suggested to act through developmental switches, but little is known about associated molecular mechanisms. In the nematode Pristionchus pacificus, the sulfatase eud-1 was identified as part of a developmental switch controlling mouth-form plasticity governing a predatory versus bacteriovorous mouth-form decision. Here we show that mutations in the conserved histone-acetyltransferase Ppa-lsy-12 and the methyl-binding-protein Ppa-mbd-2 mimic the eud-1 phenotype, resulting in the absence of one mouth-form. Mutations in both genes cause histone modification defects and reduced eud-1 expression. Surprisingly, Ppa-lsy-12 mutants also result in the down-regulation of an antisense-eud-1 RNA. eud-1 and antisense-eud-1 are co-expressed and further experiments suggest that antisense-eud-1 acts through eud-1 itself. Indeed, overexpression of the antisense-eud-1 RNA increases the eud-1-sensitive mouth-form and extends eud-1 expression. In contrast, this effect is absent in eud-1 mutants indicating that antisense-eud-1 positively regulates eud-1. Thus, chromatin remodelling and antisense-mediated up-regulation of eud-1 control feeding plasticity in Pristionchus.

The Orphan Gene dauerless Regulates Dauer Development and Intraspecific Competition in Nematodes by Copy Number Variation.

Many nematodes form dauer larvae when exposed to unfavorable conditions, representing an example of phenotypic plasticity and a major survival and dispersal strategy. In Caenorhabditis elegans, the regulation of dauer induction is a model for pheromone, insulin, and steroid-hormone signaling. Recent studies in Pristionchus pacificus revealed substantial natural variation in various aspects of dauer development, i.e. pheromone production and sensing and dauer longevity and fitness. One intriguing example is a strain from Ohio, having extremely long-lived dauers associated with very high fitness and often forming the most dauers in response to other strains´ pheromones, including the reference strain from California. While such examples have been suggested to represent intraspecific competition among strains, the molecular mechanisms underlying these dauer-associated patterns are currently unknown. We generated recombinant-inbred-lines between the Californian and Ohioan strains and used quantitative-trait-loci analysis to investigate the molecular mechanism determining natural variation in dauer development. Surprisingly, we discovered that the orphan gene dauerless controls dauer formation by copy number variation. The Ohioan strain has one dauerless copy causing high dauer formation, whereas the Californian strain has two copies, resulting in strongly reduced dauer formation. Transgenic animals expressing multiple copies do not form dauers. dauerless is exclusively expressed in CAN neurons, and both CAN ablation and dauerless mutations increase dauer formation. Strikingly, dauerless underwent several duplications and acts in parallel or downstream of steroid-hormone signaling but upstream of the nuclear-hormone-receptor daf-12. We identified the novel or fast-evolving gene dauerless as inhibitor of dauer development. Our findings reveal the importance of gene duplications and copy number variations for orphan gene function and suggest daf-12 as major target for dauer regulation. We discuss the consequences of the novel vs. fast-evolving nature of orphans for the evolution of developmental networks and their role in natural variation and intraspecific competition.

Functional specificity of the nematode Hox gene mab-5

Hox genes encode evolutionarily conserved transcription factors involved in morphological specification along the anteroposterior body axis of animals. The two most striking features of Hox genes are colinearity and the strong sequence conservation. Among all animals studied so far, the nematode Caenorhabditis elegans contains one of the most divergent Hox clusters. The core cluster contains only four members, which in part deviate from the colinearity rule. In addition, orthologous and paralogous nematode Hox sequences diverged substantially. Given these nematode-specific features, we asked how these Hox proteins evolved and how they provide functional specificity. We investigated the role of MAB-5 during ray formation and established an in vivo assay using Cel-mab-5 regulatory elements to express orthologous, paralogous and chimeric cDNAs in a Cel-mab-5 mutant background. We show that the MAB-5 ortholog from Pristionchus pacificus, but not the C. elegans paralogous Hox proteins can rescue Cel-mab-5. Experiments with chimeric, truncated and mutagenized Hox proteins suggest the specificity to be conferred by the N-terminal arm and helix I, but not helix II of the homeodomain.