A toolkit for rapid gene mapping in the nematode Caenorhabditis briggsae

Natural variation in infection specificity of Caenorhabditis briggsae isolates by two RNA viruses

Antagonistic relationships such as host-virus interactions potentially lead to rapid evolution and specificity in interactions. The Orsay virus is so far the only horizontal virus naturally infecting the nematode C. elegans. In contrast, several related RNA viruses infect its congener C. briggsae, including Santeuil (SANTV) and Le Blanc (LEBV) viruses. Here we focus on the host's intraspecific variation in sensitivity to these two intestinal viruses. Many temperate-origin C. briggsae strains, including JU1264 and JU1498, are sensitive to both, while many tropical strains, such as AF16, are resistant to both. Interestingly, some C. briggsae strains exhibit a specific resistance, such as the HK104 strain, specifically resistant to LEBV. The viral sensitivity pattern matches the strains' geographic and genomic relationships. The heavily infected strains mount a seemingly normal small RNA response that is insufficient to suppress viral infection, while the resistant strains show no small RNA response, suggesting an early block in viral entry or replication. We use a genetic approach from the host side to map genomic regions participating in viral resistance polymorphisms. Using Advanced Intercrossed Recombinant Inbred Lines (RILs) between virus-resistant AF16 and SANTV-sensitive HK104, we detect Quantitative Trait Loci (QTLs) on chromosomes IV and III. Building RILs between virus-sensitive JU1498 and LEBV-resistant HK104 followed by bulk segregant analysis, we identify a chromosome II QTL. In both cases, further introgressions of the regions confirmed the QTLs. This diversity provides an avenue for studying virus entry, replication, and exit mechanisms, as well as host-virus specificity and the host response to a specific virus infection.

Virus-like transposons cross the species barrier and drive the evolution of genetic incompatibilities

Horizontal gene transfer, the movement of genetic material between species, has been reported across all major eukaryotic lineages. However, the underlying mechanisms of transfer and their impact on genome evolution are still poorly understood. While studying the evolutionary origin of a selfish element in the nematode Caenorhabditis briggsae, we discovered that Mavericks, ancient virus-like transposons related to giant viruses and virophages, are one of the long-sought vectors of horizontal gene transfer. We found that Mavericks gained a novel herpesvirus-like fusogen in nematodes, leading to the widespread exchange of cargo genes between extremely divergent species, bypassing sexual and genetic barriers spanning hundreds of millions of years. Our results show how the union between viruses and transposons causes horizontal gene transfer and ultimately genetic incompatibilities in natural populations.

Characterization of two new C. briggsae multivulva genes

The nematode C. briggsae is an excellent genetic model for comparative and evolutionary studies involving its well-known cousin C. elegans . The vulval system in these two species has been used extensively to investigate genes and pathways involved in cell proliferation and cell differentiation. Here we report initial characterization of two C. briggsae multivulva (Muv) mutants, Cbr-lin(bh1) and Cbr-lin(bh3) .

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.

Fast genetic mapping using insertion-deletion polymorphisms in Caenorhabditis elegans

Genetic mapping is used in forward genetics to narrow the list of candidate mutations and genes corresponding to the mutant phenotype of interest. Even with modern advances in biology such as efficient identification of candidate mutations by whole-genome sequencing, mapping remains critical in pinpointing the responsible mutation. Here we describe a simple, fast, and affordable mapping toolkit that is particularly suitable for mapping in Caenorhabditis elegans. This mapping method uses insertion-deletion polymorphisms or indels that could be easily detected instead of single nucleotide polymorphisms in commonly used Hawaiian CB4856 mapping strain. The materials and methods were optimized so that mapping could be performed using tiny amount of genetic material without growing many large populations of mutants for DNA purification. We performed mapping of previously known and unknown mutations to show strengths and weaknesses of this method and to present examples of completed mapping. For situations where Hawaiian CB4856 is unsuitable, we provide an annotated list of indels as a basis for fast and easy mapping using other wild isolates. Finally, we provide rationale for using this mapping method over other alternatives as a part of a comprehensive strategy also involving whole-genome sequencing and other methods.

A broad mutational target explains a fast rate of phenotypic evolution

The rapid evolution of a trait in a clade of organisms can be explained by the sustained action of natural selection or by a high mutational variance, that is the propensity to change under spontaneous mutation. The causes for a high mutational variance are still elusive. In some cases, fast evolution depends on the high mutation rate of one or few loci with short tandem repeats. Here, we report on the fastest evolving cell fate among vulva precursor cells in Caenorhabditis nematodes, that of P3.p. We identify and validate causal mutations underlying P3.p's high mutational variance. We find that these positions do not present any characteristics of a high mutation rate, are scattered across the genome and the corresponding genes belong to distinct biological pathways. Our data indicate that a broad mutational target size is the cause of the high mutational variance and of the corresponding fast phenotypic evolutionary rate.

Heritable characteristics or traits of a group of organisms, for example the large brain size of primates or the hooves of a horse, are determined by genes, the environment, and by the interactions between them. Traits can change over time and generations when enough mutations in these genes have spread in a species to result in visible differences.

However, some traits, such as the large brain of primates, evolve faster than others, but why this is the case has been unclear. It could be that a few specific genes important for that trait in question mutate at a high rate, or, that many genes affect the trait, creating a lot of variation for natural selection to choose from.

Here, Besnard, Picao-Osorio et al. studied the roundworm Caenorhabditis elegans to better understand the causes underlying the different rates of trait evolution. These worms have a short life cycle and evolve quickly over many generations, making them an ideal candidate for studying mutation rates in different traits.

Previous studies have shown that one of C. elegans’ six cells of the reproductive system evolves faster than the others. To investigate this further, Besnard, Picao-Osorio et al. analysed the genetic mutations driving change in this cell in 250 worm generations. The results showed that five mutations in five different genes – all responsible for different processes in the cells – were behind the supercharged evolution of this particular cell. This suggests that fast evolution results from natural selection acting upon a collection of genes, rather than one gene, and that many genes and pathways shape this trait.

In conclusion, these results demonstrate that how traits are coded at the molecular level, in one gene or many, can influence the rate at which they evolve.

Genetically Distinct Behavioral Modules Underlie Natural Variation in Thermal Performance Curves

Thermal reaction norms pervade organismal traits as stereotyped responses to temperature, a fundamental environmental input into sensory and physiological systems. Locomotory behavior represents an especially plastic read-out of animal response, with its dynamic dependence on environmental stimuli presenting a challenge for analysis and for understanding the genomic architecture of heritable variation. Here we characterize behavioral reaction norms as thermal performance curves for the nematode Caenorhabditis briggsae, using a collection of 23 wild isolate genotypes and 153 recombinant inbred lines to quantify the extent of genetic and plastic variation in locomotory behavior to temperature changes. By reducing the dimensionality of the multivariate phenotypic response with a function-valued trait framework, we identified genetically distinct behavioral modules that contribute to the heritable variation in the emergent overall behavioral thermal performance curve. Quantitative trait locus mapping isolated regions on Chromosome II associated with locomotory activity at benign temperatures and Chromosome V loci related to distinct aspects of sensitivity to high temperatures, with each quantitative trait locus explaining up to 28% of trait variation. These findings highlight how behavioral responses to environmental inputs as thermal reaction norms can evolve through independent changes to genetically distinct modular components of such complex phenotypes.

The Genetic Architecture of Intra-Species Hybrid Mito-Nuclear Epistasis

Genetic variants that are neutral within, but deleterious between, populations (Dobzhansky-Muller Incompatibilities) are thought to initiate hybrid dysfunction and then to accumulate and complete the speciation process. To identify the types of genetic differences that might initiate speciation, it is useful to study inter-population (intra-species) hybrids that exhibit reduced fitness. In Caenorhabditis briggsae, a close relative of the nematode C. elegans, such minor genetic incompatibilities have been identified. One incompatibility between the mitochondrial and nuclear genomes reduces the fitness of some hybrids. To understand the nuclear genetic architecture of this epistatic interaction, we constructed two sets of recombinant inbred lines by hybridizing two genetically diverse wild populations. In such lines, selection is able to eliminate deleterious combinations of alleles derived from the two parental populations. The genotypes of surviving hybrid lines thus reveal favorable allele combinations at loci experiencing selection. Our genotype data from the resulting lines are consistent with the interpretation that the X alleles participate in epistatic interactions with autosomes and the mitochondrial genome. We evaluate this possibility given predictions that mitochondria-X epistasis should be more prevalent than between mitochondria and autosomes. Our empirical identification of inter-genomic linkage disequilibrium supports the body of literature indicating that the accumulation of mito-nuclear genetic incompatibilities might initiate the speciation process through the generation of less-fit inter-population hybrids.

Toward Universal Forward Genetics: Using a Draft Genome Sequence of the Nematode Oscheius tipulae To Identify Mutations Affecting Vulva Development

Mapping-by-sequencing has become a standard method to map and identify phenotype-causing mutations in model species. Here, we show that a fragmented draft assembly is sufficient to perform mapping-by-sequencing in nonmodel species. We generated a draft assembly and annotation of the genome of the free-living nematode Oscheius tipulae, a distant relative of the model Caenorhabditis elegans We used this draft to identify the likely causative mutations at the O. tipulae cov-3 locus, which affect vulval development. The cov-3 locus encodes the O. tipulae ortholog of C. elegans mig-13, and we further show that Cel-mig-13 mutants also have an unsuspected vulval-development phenotype. In a virtuous circle, we were able to use the linkage information collected during mutant mapping to improve the genome assembly. These results showcase the promise of genome-enabled forward genetics in nonmodel species.

Revisiting Suppression of Interspecies Hybrid Male Lethality in Caenorhabditis Nematodes

Within the nematode genus Caenorhabditis, Caenorhabditis briggsae and C. nigoni are among the most closely related species known. They differ in sexual mode, with C. nigoni retaining the ancestral XO male-XX female outcrossing system, while C. briggsae recently evolved self-fertility and an XX-biased sex ratio. Wild-type C. briggsae and C. nigoni can produce fertile hybrid XX female progeny, but XO progeny are either 100% inviable (when C. briggsae is the mother) or viable but sterile (when C. nigoni is the mother). A recent study provided evidence suggesting that loss of the Cbr-him-8 meiotic regulator in C. briggsae hermaphrodites allowed them to produce viable and fertile hybrid XO male progeny when mated to C. nigoni Because such males would be useful for a variety of genetic experiments, we sought to verify this result. Preliminary crosses with wild-type C. briggsae hermaphrodites occasionally produced fertile males, but they could not be confirmed to be interspecies hybrids. Using an RNA interference (RNAi) protocol that eliminates any possibility of self-progeny in Cbr-him-8 hermaphrodites, we found sterile males bearing the C. nigoni X chromosome, but no fertile males bearing the C. briggsae X, as in wild-type crosses. Our results suggest that the apparent rescue of XO hybrid viability and fertility is due to incomplete purging of self-sperm prior to mating.

Differing roles for sur-2/MED23 in C. elegans and C. briggsae vulval development

Normal vulval development in the nematode Caenorhabditis briggsae is identical to that in the related Caenorhabditis elegans. However, several experiments suggest that there are differences between the two species with respect to the contribution of EGF/Ras signaling. To investigate these differences genetically, we have characterized a C. briggsae mutant strain that phenocopies the effect observed when C. briggsae animals are treated with U0126, an inhibitor of the EGF pathway component MEK. We identify that the gene affected in the mutant strain is Cbr-sur-2, which encodes a MED23 mediator complex protein that acts downstream of EGF signaling in C. elegans and other organisms, such as mammals. When Cbr-sur-2 and Cel-sur-2 mutants are compared, we find that the production of additional vulval cells from P5.p and P7.p in C. elegans is dependent on proper development of P6.p, while C. briggsae does not have a similar requirement. Combined chemical and genetic interference with the EGF pathway completely eliminates vulval development in C. elegans but not in C. briggsae. Our results provide genetic evidence for the differing requirements for EGF signaling in the two species.

Paternal Mitochondrial Transmission in Intra-Species Caenorhabditis briggsae Hybrids

To study mitochondrial–nuclear genetic interactions in the nematode Caenorhabditis briggsae, our three laboratories independently created 38 intra-species cytoplasmic–nuclear hybrid (cybrid) lines. Although the cross design combines maternal mitotypes with paternal nuclear genotypes, eight lines (21%) unexpectedly contained paternal mitotypes. All eight share in common ancestry of one of two genetically related strains. This unexpected parallel observation of paternal mitochondrial transmission, undesirable given our intent of creating cybrids, provides a serendipitous experimental model and framework to study the molecular and evolutionary basis of uniparental mitochondrial inheritance.

Genetics: modes of reproduction and genetic analysis

Classical and reverse genetics remain invaluable tools for the scientific investigation of model organisms. Genetic analysis of endoparasites is generally difficult because the sexual adults required for crossing and other manipulations are usually hidden within their host. Strongyloides spp. and Parastrongyloides spp. are notable exceptions to this and their free-living adults offer unique opportunities to manipulate these parasites experimentally. Here I review the modes of inheritance in the two generations of Strongyloides/Parastrongyloides and I discuss the opportunities and the limitations of the currently available methodology for the genetic analysis of these two genera.

Mitochondrial-Nuclear Epistasis Impacts Fitness and Mitochondrial Physiology of Interpopulation Caenorhabditis briggsae Hybrids

In order to identify the earliest genetic changes that precipitate species formation, it is useful to study genetic incompatibilities that cause only mild dysfunction when incompatible alleles are combined in an interpopulation hybrid. Such hybridization within the nematode species Caenorhabditis briggsae has been suggested to result in selection against certain combinations of nuclear and mitochondrial alleles, raising the possibility that mitochondrial-nuclear (mitonuclear) epistasis reduces hybrid fitness. To test this hypothesis, cytoplasmic-nuclear hybrids (cybrids) were created to purposefully disrupt any epistatic interactions. Experimental analysis of the cybrids suggests that mitonuclear discord can result in decreased fecundity, increased lipid content, and increased mitochondrial reactive oxygen species levels. Many of these effects were asymmetric with respect to cross direction, as expected if cytoplasmic-nuclear Dobzhansky-Muller incompatibilities exist. One such effect is consistent with the interpretation that disrupting coevolved mitochondrial and nuclear loci impacts mitochondrial function and organismal fitness. These findings enhance efforts to study the genesis, identity, and maintenance of genetic incompatibilities that precipitate the speciation process.

Full-genome evolutionary histories of selfing, splitting, and selection in Caenorhabditis

The nematode Caenorhabditis briggsae is a model for comparative developmental evolution with C. elegans. Worldwide collections of C. briggsae have implicated an intriguing history of divergence among genetic groups separated by latitude, or by restricted geography, that is being exploited to dissect the genetic basis to adaptive evolution and reproductive incompatibility; yet, the genomic scope and timing of population divergence is unclear. We performed high-coverage whole-genome sequencing of 37 wild isolates of the nematode C. briggsae and applied a pairwise sequentially Markovian coalescent (PSMC) model to 703 combinations of genomic haplotypes to draw inferences about population history, the genomic scope of natural selection, and to compare with 40 wild isolates of C. elegans. We estimate that a diaspora of at least six distinct C. briggsae lineages separated from one another approximately 200,000 generations ago, including the “Temperate” and “Tropical” phylogeographic groups that dominate most samples worldwide. Moreover, an ancient population split in its history approximately 2 million generations ago, coupled with only rare gene flow among lineage groups, validates this system as a model for incipient speciation. Low versus high recombination regions of the genome give distinct signatures of population size change through time, indicative of widespread effects of selection on highly linked portions of the genome owing to extreme inbreeding by self-fertilization. Analysis of functional mutations indicates that genomic context, owing to selection that acts on long linkage blocks, is a more important driver of population variation than are the functional attributes of the individually encoded genes.

A Genome-wide hybrid incompatibility landscape between Caenorhabditis briggsae and C. nigoni

Systematic characterization of ẖybrid incompatibility (HI) between related species remains the key to understanding speciation. The genetic basis of HI has been intensively studied in Drosophila species, but remains largely unknown in other species, including nematodes, which is mainly due to the lack of a sister species with which C. elegans can mate and produce viable progeny. The recent discovery of a C. briggsae sister species, C. nigoni, has opened up the possibility of dissecting the genetic basis of HI in nematode species. However, the paucity of dominant and visible marker prevents the efficient mapping of HI loci between the two species. To elucidate the genetic basis of speciation in nematode species, we first generated 96 chromosomally integrated GFP markers in the C. briggsae genome and mapped them into the defined locations by PCR and Next-Generation Sequencing (NGS). Aided by the marker, we backcrossed the GFP-associated C. briggsae genomic fragments into C. nigoni for at least 15 generations and produced 111 independent introgressions. The introgression fragments cover most of the C. briggsae genome. We finally dissected the patterns of HI by scoring the embryonic lethality, larval arrest, sex ratio and male sterility for each introgression line, through which we identified pervasive HI loci and produced a genome-wide landscape of HI between the two nematode species, the first of its type for any non-Drosophila species. The HI data not only provided insights into the genetic basis of speciation, but also established a framework for the possible cloning of HI loci between the two nematode species. Furthermore, the data on hybrids confirmed Haldane’s rule and suggested the presence of a large X effect in terms of fertility between the two species. Importantly, this work opens a new avenue for studying speciation genetics between nematode species and allows parallel comparison of the HI with that in Drosophila and other species.

Hybrid incompatibility (HI) has been intensively studied among Drosophila species, but remains largely unknown in other species. Model organism is a species of choice for the HI study because these species provide sophisticated molecular and genetic tools for illustrating mechanism underlying a given HI. C. elegans as a model organism contributed little to the field due to the lack of a sister species with which it can mate and produce viable progeny. The recent discovery of a C. briggsae sister species, C. nigoni, has made it possible to study the HI between nematode species. However, the paucity of molecular and genetic tools in both species inhibits their use in such studies. To empower the use of this species pair in HI studies, we first created a collection of fluorescent markers over the C. briggsae genome to facilitate the directional introduction of the marker-associated C. briggsae genomic fragments into the C. nigoni background. We next mapped the marker insertion sites and introduced the markers into C. nigoni by repeated crossings. Finally, we generated a genome-wide HI landscape between the two species by scoring the HI phenotypes of their hybrid progeny. The study not only provides an invaluable resource for the molecular cloning of HI loci between C. briggsae and C. nigoni, but also permits comparative analysis of speciation genetics between nematode and other species.

Mutations in Caenorhabditis briggsae identify new genes important for limiting the response to EGF signaling during vulval development

Studies of vulval development in the nematode C. elegans have identified many genes that are involved in cell division and differentiation processes. Some of these encode components of conserved signal transduction pathways mediated by EGF, Notch, and Wnt. To understand how developmental mechanisms change during evolution, we are doing a comparative analysis of vulva formation in C. briggsae, a species that is closely related to C. elegans. Here, we report 14 mutations in 7 Multivulva (Muv) genes in C. briggsae that inhibit inappropriate division of vulval precursors. We have developed a new efficient and cost-effective gene mapping method to localize Muv mutations to small genetic intervals on chromosomes, thus facilitating cloning and functional studies. We demonstrate the utility of our method by determining molecular identities of three of the Muv genes that include orthologs of Cel-lin-1 (ETS) and Cel-lin-31 (Winged-Helix) of the EGF-Ras pathway and Cel-pry-1 (Axin), of the Wnt pathway. The remaining four genes reside in regions that lack orthologs of known C. elegans Muv genes. Inhibitor studies demonstrate that the Muv phenotype of all four new genes is dependent on the activity of the EGF pathway kinase, MEK. One of these, Cbr-lin(gu167), shows modest increase in the expression of Cbr-lin-3/EGF compared to wild type. These results argue that while Cbr-lin(gu167) may act upstream of Cbr-lin-3/EGF, the other three genes influence the EGF pathway downstream or in parallel to Cbr-lin-3. Overall, our findings demonstrate that the genetic program underlying a conserved developmental process includes both conserved and divergent functional contributions.

Fundamentals of Comparative Genome Analysis in Caenorhabditis Nematodes

The genome of the nematode Caenorhabditis elegans was the first of any animal to be sequenced completely, and it remains the "gold standard" for completeness and annotations. Even before the C. elegans genome was completed, however, biologists began examining the generality of its features in the genomes of other Caenorhabditis species. With many such genomes now sequenced and available via WormBase, C. elegans researchers are often confronted with how to interpret comparative genomic data. In this article, we present practical approaches to addressing several common issues, including possible sources of error in homology annotations, the often complex relationships between sequence similarity, orthology, paralogy, and gene family evolution, the impact of sexual mode on genome assemblies and content, and the determination and use of synteny as a tool.

Comparative functional characterization of the CSR-1 22G-RNA pathway in Caenorhabditis nematodes

As a champion of small RNA research for two decades, Caenorhabditis elegans has revealed the essential Argonaute CSR-1 to play key nuclear roles in modulating chromatin, chromosome segregation and germline gene expression via 22G-small RNAs. Despite CSR-1 being preserved among diverse nematodes, the conservation and divergence in function of the targets of small RNA pathways remains poorly resolved. Here we apply comparative functional genomic analysis between C. elegans and Caenorhabditis briggsae to characterize the CSR-1 pathway, its targets and their evolution. C. briggsae CSR-1-associated small RNAs that we identified by immunoprecipitation-small RNA sequencing overlap with 22G-RNAs depleted in cbr-csr-1 RNAi-treated worms. By comparing 22G-RNAs and target genes between species, we defined a set of CSR-1 target genes with conserved germline expression, enrichment in operons and more slowly evolving coding sequences than other genes, along with a small group of evolutionarily labile targets. We demonstrate that the association of CSR-1 with chromatin is preserved, and show that depletion of cbr-csr-1 leads to chromosome segregation defects and embryonic lethality. This first comparative characterization of a small RNA pathway in Caenorhabditis establishes a conserved nuclear role for CSR-1 and highlights its key role in germline gene regulation across multiple animal species.

Evidence of hermaphroditism and sex ratio distortion in the fungal feeding nematode Bursaphelenchus okinawaensis

Nematodes have many different reproductive strategies along with their divergent life histories; the ability of hermaphrodite to self- and cross-fertilize is useful for genetic manipulation. Here, we demonstrate the hermaphroditism of the fungal feeding nematode Bursaphelenchus okinawaensis, which was formerly described as a parthenogenetic nematode, and we show its other unique sexual characteristics. To determine that it is hermaphroditic, we performed the following experiments: observation of the pronuclear and chromosome behavior during oogenesis and early embryogenesis; observation of spermatogenesis during the fourth larval stage; investigation of sperm utilization; and investigation of phenotypic segregation after cross-mating using a chemically induced visible mutant. We then investigated the mating preferences and spermatid size difference between males and hermaphrodites. B. okinawaensis males successfully mated only with sperm-depleted old hermaphrodites, and the spermatid sizes of males were almost the same as those of hermaphrodites. Moreover, the sex ratio of cross-fertilized progeny was highly skewed toward hermaphrodites. B. okinawaensis is phylogenetically distant from established model nematodes such as C. elegans and is more closely related to some economically relevant parasitic nematodes. This newly discovered hermaphroditic nematode has great potential for evolutionary and parasitological research.

Sex-biased gene expression and evolution of the x chromosome in nematodes

Studies of X chromosome evolution in various organisms have indicated that sex-biased genes are nonrandomly distributed between the X and autosomes. Here, to extend these studies to nematodes, we annotated and analyzed X chromosome gene content in four Caenorhabditis species and in Pristionchus pacificus. Our gene expression analyses comparing young adult male and female mRNA-seq data indicate that, in general, nematode X chromosomes are enriched for genes with high female-biased expression and depleted of genes with high male-biased expression. Genes with low sex-biased expression do not show the same trend of X chromosome enrichment and depletion. Combined with the observation that highly sex-biased genes are primarily expressed in the gonad, differential distribution of sex-biased genes reflects differences in evolutionary pressures linked to tissue-specific regulation of X chromosome transcription. Our data also indicate that X dosage imbalance between males (XO) and females (XX) is influential in shaping both expression and gene content of the X chromosome. Predicted upregulation of the single male X to match autosomal transcription (Ohno's hypothesis) is supported by our observation that overall transcript levels from the X and autosomes are similar for highly expressed genes. However, comparison of differentially located one-to-one orthologs between C. elegans and P. pacificus indicates lower expression of X-linked orthologs, arguing against X upregulation. These contradicting observations may be reconciled if X upregulation is not a global mechanism but instead acts locally on a subset of tissues and X-linked genes that are dosage sensitive.

Rapid creation of forward-genetics tools for C. briggsae using TALENs: lessons for nonmodel organisms

Although evolutionary studies of gene function often rely on RNA interference, the ideal approach would use reverse genetics to create null mutations for cross-species comparisons and forward genetics to identify novel genes in each species. We have used transcription activator-like effector nucleases (TALENs) to facilitate both approaches in Caenorhabditis nematodes. First, by combining golden gate cloning and TALEN technology, we can induce frameshifting mutations in any gene. Second, by combining this approach with bioinformatics we can predict and create the resources needed for forward genetic analysis in species like Caenorhabditis briggsae. Although developing genetic model organisms used to require years to isolate marker mutations, balancers, and tools, with TALENs, these reagents can now be produced in months. Furthermore, the analysis of nonsense mutants in related model organisms allows a directed approach for making these markers and tools. When used together, these methods could simplify the adaptation of other organisms for forward and reverse genetics.

Evolutionary change within a bipotential switch shaped the sperm/oocyte decision in hermaphroditic nematodes

A subset of transcription factors like Gli2 and Oct1 are bipotential--they can activate or repress the same target, in response to changing signals from upstream genes. Some previous studies implied that the sex-determination protein TRA-1 might also be bipotential; here we confirm this hypothesis by identifying a co-factor, and use it to explore how the structure of a bipotential switch changes during evolution. First, null mutants reveal that C. briggsae TRR-1 is required for spermatogenesis, RNA interference implies that it works as part of the Tip60 Histone Acetyl Transferase complex, and RT-PCR data show that it promotes the expression of Cbr-fog-3, a gene needed for spermatogenesis. Second, epistasis tests reveal that TRR-1 works through TRA-1, both to activate Cbr-fog-3 and to control the sperm/oocyte decision. Since previous studies showed that TRA-1 can repress fog-3 as well, these observations demonstrate that it is bipotential. Third, TRR-1 also regulates the development of the male tail. Since Cbr-tra-2 Cbr-trr-1 double mutants resemble Cbr-tra-1 null mutants, these two regulatory branches control all tra-1 activity. Fourth, striking differences in the relationship between these two branches of the switch have arisen during recent evolution. C. briggsae trr-1 null mutants prevent hermaphrodite spermatogenesis, but not Cbr-fem null mutants, which disrupt the other half of the switch. On the other hand, C. elegans fem null mutants prevent spermatogenesis, but not Cel-trr-1 mutants. However, synthetic interactions confirm that both halves of the switch exist in each species. Thus, the relationship between the two halves of a bipotential switch can shift rapidly during evolution, so that the same phenotype is produce by alternative, complementary mechanisms.

Genetic control of vulval development in Caenorhabditis briggsae

The nematode Caenorhabditis briggsae is an excellent model organism for the comparative analysis of gene function and developmental mechanisms. To study the evolutionary conservation and divergence of genetic pathways mediating vulva formation, we screened for mutations in C. briggsae that cause the egg-laying defective (Egl) phenotype. Here, we report the characterization of 13 genes, including three that are orthologs of Caenorhabditis elegans unc-84 (SUN domain), lin-39 (Dfd/Scr-related homeobox), and lin-11 (LIM homeobox). Based on the morphology and cell fate changes, the mutants were placed into four different categories. Class 1 animals have normal-looking vulva and vulva-uterine connections, indicating defects in other components of the egg-laying system. Class 2 animals frequently lack some or all of the vulval precursor cells (VPCs) due to defects in the migration of P-cell nuclei into the ventral hypodermal region. Class 3 animals show inappropriate fusion of VPCs to the hypodermal syncytium, leading to a reduced number of vulval progeny. Finally, class 4 animals exhibit abnormal vulval invagination and morphology. Interestingly, we did not find mutations that affect VPC induction and fates. Our work is the first study involving the characterization of genes in C. briggsae vulva formation, and it offers a basis for future investigations of these genes in C. elegans.

A toolkit for bulk PCR-based marker design from next-generation sequence data: application for development of a framework linkage map in bulb onion (Allium cepa L.)

Background

Although modern sequencing technologies permit the ready detection of numerous DNA sequence variants in any organisms, converting such information to PCR-based genetic markers is hampered by a lack of simple, scalable tools. Onion is an example of an under-researched crop with a complex, heterozygous genome where genome-based research has previously been hindered by limited sequence resources and genetic markers.

Results

We report the development of generic tools for large-scale web-based PCR-based marker design in the Galaxy bioinformatics framework, and their application for development of next-generation genetics resources in a wide cross of bulb onion (Allium cepa L.). Transcriptome sequence resources were developed for the homozygous doubled-haploid bulb onion line ‘CUDH2150’ and the genetically distant Indian landrace ‘Nasik Red’, using 454™ sequencing of normalised cDNA libraries of leaf and shoot. Read mapping of ‘Nasik Red’ reads onto ‘CUDH2150’ assemblies revealed 16836 indel and SNP polymorphisms that were mined for portable PCR-based marker development. Tools for detection of restriction polymorphisms and primer set design were developed in BioPython and adapted for use in the Galaxy workflow environment, enabling large-scale and targeted assay design. Using PCR-based markers designed with these tools, a framework genetic linkage map of over 800cM spanning all chromosomes was developed in a subset of 93 F₂ progeny from a very large F₂ family developed from the ‘Nasik Red’ x ‘CUDH2150’ inter-cross. The utility of tools and genetic resources developed was tested by designing markers to transcription factor-like polymorphic sequences. Bin mapping these markers using a subset of 10 progeny confirmed the ability to place markers within 10 cM bins, enabling increased efficiency in marker assignment and targeted map refinement. The major genetic loci conditioning red bulb colour (R) and fructan content (Frc) were located on this map by QTL analysis.

Conclusions

The generic tools developed for the Galaxy environment enable rapid development of sets of PCR assays targeting sequence variants identified from Illumina and 454 sequence data. They enable non-specialist users to validate and exploit large volumes of next-generation sequence data using basic equipment.

Temperature-dependent behaviours are genetically variable in the nematode Caenorhabditis briggsae

Temperature-dependent behaviours in Caenorhabditis elegans, such as thermotaxis and isothermal tracking, are complex behavioural responses that integrate sensation, foraging and learning, and have driven investigations to discover many essential genetic and neural pathways. The ease of manipulation of the Caenorhabditis model system also has encouraged its application to comparative analyses of phenotypic evolution, particularly contrasts of the classic model C. elegans with C. briggsae. And yet few studies have investigated natural genetic variation in behaviour in any nematode. Here we measure thermotaxis and isothermal tracking behaviour in genetically distinct strains of C. briggsae, further motivated by the latitudinal differentiation in C. briggsae that is associated with temperature-dependent fitness differences in this species. We demonstrate that C. briggsae performs thermotaxis and isothermal tracking largely similar to that of C. elegans, with a tendency to prefer its rearing temperature. Comparisons of these behaviours among strains reveal substantial heritable natural variation within each species that corresponds to three general patterns of behavioural response. However, intraspecific genetic differences in thermal behaviour often exceed interspecific differences. These patterns of temperature-dependent behaviour motivate further development of C. briggsae as a model system for dissecting the genetic underpinnings of complex behavioural traits.

A method for rapid and simultaneous mapping of genetic loci and introgression sizes in nematode species

Caenorhabditis briggsae is emerging as an attractive model organism not only in studying comparative biology against C. elegans, but also in developing novel experimentation avenues. In particular, recent identification of a new Caenorhabditis species, C. sp.9 with which it can mate and produce viable progeny provides an opportunity for studying the genetics of hybrid incompatibilities (HI) between the two. Mapping of a specific HI locus demands repeated backcrossing to get hold of the specific genomic region underlying an observed phenotype. To facilitate mapping of HI loci between C. briggsae and C. sp.9, an efficient mapping method and a genetic map ideally consisting of dominant markers are required for systematic introgression of genomic fragments between the two species. We developed a fast and cost-effective method for high throughput mapping of dominant loci with resolution up to 1 million bps in C. briggsae. The method takes advantage of the introgression between C. briggsae and C. sp.9 followed by PCR genotyping using C. briggsae specific primers. Importantly, the mapping results can not only serve as an effective way for estimating the chromosomal position of a genetic locus in C. briggsae, but also provides size information for the introgression fragment in an otherwise C. sp.9 background. In addition, it also helps generate introgression line as a side-product that is invaluable for the subsequent mapping of HI loci. The method will greatly facilitate the construction of a genetic map consisting of dominant markers and pave the way for systematic isolation of HI loci between C. briggsae and C. sp.9 which has so far not been attempted between nematode species. The method is designed for mapping of a dominant allele, but can be easily adapted for mapping of any other type of alleles in any other species if introgression between a sister species pair is feasible.

In vivo quantification reveals extensive natural variation in mitochondrial form and function in Caenorhabditis briggsae

We have analyzed natural variation in mitochondrial form and function among a set of Caenorhabditis briggsae isolates known to harbor mitochondrial DNA structural variation in the form of a heteroplasmic nad5 gene deletion (nad5Δ) that correlates negatively with organismal fitness. We performed in vivo quantification of 24 mitochondrial phenotypes including reactive oxygen species level, membrane potential, and aspects of organelle morphology, and observed significant among-isolate variation in 18 traits. Although several mitochondrial phenotypes were non-linearly associated with nad5Δ levels, most of the among-isolate phenotypic variation could be accounted for by phylogeographic clade membership. In particular, isolate-specific mitochondrial membrane potential was an excellent predictor of clade membership. We interpret this result in light of recent evidence for local adaptation to temperature in C. briggsae. Analysis of mitochondrial-nuclear hybrid strains provided support for both mtDNA and nuclear genetic variation as drivers of natural mitochondrial phenotype variation. This study demonstrates that multicellular eukaryotic species are capable of extensive natural variation in organellar phenotypes and highlights the potential of integrating evolutionary and cell biology perspectives.

Variation in base-substitution mutation in experimental and natural lineages of Caenorhabditis nematodes

Variation among lineages in the mutation process has the potential to impact diverse biological processes ranging from susceptibilities to genetic disease to the mode and tempo of molecular evolution. The combination of high-throughput DNA sequencing (HTS) with mutation-accumulation (MA) experiments has provided a powerful approach to genome-wide mutation analysis, though insights into mutational variation have been limited by the vast evolutionary distances among the few species analyzed. We performed a HTS analysis of MA lines derived from four Caenorhabditis nematode natural genotypes: C. elegans N2 and PB306 and C. briggsae HK104 and PB800. Total mutation rates did not differ among the four sets of MA lines. A mutational bias toward G:C→A:T transitions and G:C→T:A transversions was observed in all four sets of MA lines. Chromosome-specific rates were mostly stable, though there was some evidence for a slightly elevated X chromosome mutation rate in PB306. Rates were homogeneous among functional coding sequence types and across autosomal cores, arms, and tips. Mutation spectra were similar among the four MA line sets but differed significantly when compared with patterns of natural base-substitution polymorphism for 13/14 comparisons performed. Our findings show that base-substitution mutation processes in these closely related animal lineages are mostly stable but differ from natural polymorphism patterns in these two species.

Caenorhabditis briggsae recombinant inbred line genotypes reveal inter-strain incompatibility and the evolution of recombination

The nematode Caenorhabditis briggsae is an emerging model organism that allows evolutionary comparisons with C. elegans and exploration of its own unique biological attributes. To produce a high-resolution C. briggsae recombination map, recombinant inbred lines were generated from reciprocal crosses between two strains and genotyped at over 1,000 loci. A second set of recombinant inbred lines involving a third strain was also genotyped at lower resolution. The resulting recombination maps exhibit discrete domains of high and low recombination, as in C. elegans, indicating these are a general feature of Caenorhabditis species. The proportion of a chromosome's physical size occupied by the central, low-recombination domain is highly correlated between species. However, the C. briggsae intra-species comparison reveals striking variation in the distribution of recombination between domains. Hybrid lines made with the more divergent pair of strains also exhibit pervasive marker transmission ratio distortion, evidence of selection acting on hybrid genotypes. The strongest effect, on chromosome III, is explained by a developmental delay phenotype exhibited by some hybrid F2 animals. In addition, on chromosomes IV and V, cross direction-specific biases towards one parental genotype suggest the existence of cytonuclear epistatic interactions. These interactions are discussed in relation to surprising mitochondrial genome polymorphism in C. briggsae, evidence that the two strains diverged in allopatry, the potential for local adaptation, and the evolution of Dobzhansky-Muller incompatibilities. The genetic and genomic resources resulting from this work will support future efforts to understand inter-strain divergence as well as facilitate studies of gene function, natural variation, and the evolution of recombination in Caenorhabditis nematodes.

Comparative genetics and genomics of nematodes: genome structure, development, and lifestyle

Nematodes are found in virtually all habitats on earth. Many of them are parasites of plants and animals, including humans. The free-living nematode, Caenorhabditis elegans, is one of the genetically best-studied model organisms and was the first metazoan whose genome was fully sequenced. In recent years, the draft genome sequences of another six nematodes representing four of the five major clades of nematodes were published. Compared to mammalian genomes, all these genomes are very small. Nevertheless, they contain almost the same number of genes as the human genome. Nematodes are therefore a very attractive system for comparative genetic and genomic studies, with C. elegans as an excellent baseline. Here, we review the efforts that were made to extend genetic analysis to nematodes other than C. elegans, and we compare the seven available nematode genomes. One of the most striking findings is the unexpectedly high incidence of gene acquisition through horizontal gene transfer (HGT).

Insights into species divergence and the evolution of hermaphroditism from fertile interspecies hybrids of Caenorhabditis nematodes

The architecture of both phenotypic variation and reproductive isolation are important problems in evolutionary genetics. The nematode genus Caenorhabditis includes both gonochoristic (male/female) and androdioecious (male/hermaprodite) species. However, the natural genetic variants distinguishing reproductive mode remain unknown, and nothing is known about the genetic basis of postzygotic isolation in the genus. Here we describe the hybrid genetics of the first Caenorhabditis species pair capable of producing fertile hybrid progeny, the gonochoristic Caenorhabditis sp. 9 and the androdioecious C. briggsae. Though many interspecies F(1) arrest during embryogenesis, a viable subset develops into fertile females and sterile males. Reciprocal parental crosses reveal asymmetry in male-specific viability, female fertility, and backcross viability. Selfing and spermatogenesis are extremely rare in XX F(1), and almost all hybrid self-progeny are inviable. Consistent with this, F(1) females do not express male-specific molecular germline markers. We also investigated three approaches to producing hybrid hermaphrodites. A dominant mutagenesis screen for self-fertile F(1) hybrids was unsuccessful. Polyploid F(1) hybrids with increased C. briggsae genomic material did show elevated rates of selfing, but selfed progeny were mostly inviable. Finally, the use of backcrosses to render the hybrid genome partial homozygous for C. briggsae alleles did not increase the incidence of selfing or spermatogenesis relative to the F(1) generation. These hybrid animals were genotyped at 23 loci, and significant segregation distortion (biased against C. briggsae) was detected at 13 loci. This, combined with an absence of productive hybrid selfing, prevents formulation of simple hypotheses about the genetic architecture of hermaphroditism. In the near future, this hybrid system will likely be fruitful for understanding the genetics of reproductive isolation in Caenorhabditis.

Conserved mechanism of Wnt signaling function in the specification of vulval precursor fates in C. elegans and C. briggsae

The C. elegans hermaphrodite vulva serves as a paradigm for understanding how signaling pathways control organ formation. Previous studies have shown that Wnt signaling plays important roles in vulval development. To understand the function and evolution of Wnt signaling in Caenorhabditis nematodes we focused on C. briggsae, a species that is substantially divergent from C. elegans in terms of the evolutionary time scale yet shares almost identical morphology. We isolated mutants in C. briggsae that display multiple pseudo-vulvae resulting from ectopic VPC induction. We cloned one of these loci and found that it encodes an Axin homolog, Cbr-PRY-1. Our genetic studies revealed that Cbr-pry-1 functions upstream of the canonical Wnt pathway components Cbr-bar-1 (beta-catenin) and Cbr-pop-1(tcf/lef) as well as the Hox target Cbr-lin-39 (Dfd/Scr). We further characterized the pry-1 vulval phenotype in C. briggsae and C. elegans using 8 cell fate markers, cell ablation, and genetic interaction approaches. Our results show that ectopically induced VPCs in pry-1 mutants adopt 2° fates independently of the gonad-derived inductive and LIN-12/Notch-mediated lateral signaling pathways. We also found that Cbr-pry-1 mutants frequently show a failure of P7.p induction. A similar, albeit low penetrant, defect is also observed in C. elegans pry-1 mutants. The genetic analysis of the P7.p induction defect revealed that it was caused by altered regulation of lin-12 and its transcriptional target lip-1 (MAP kinase phosphatase). Thus, our results provide evidence for LIN-12/Notch-dependent and independent roles of Wnt signaling in promoting 2 degrees VPC fates in both nematode species.