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

Highly efficient transgenesis with miniMos in Caenorhabditis briggsae

Caenorhabditis briggsae as a companion species for Caenorhabditis elegans has played an increasingly important role in study of evolution of development and genome and gene regulation. Aided by the isolation of its sister spices, it has recently been established as a model for speciation study. To take full advantage of the species for comparative study, an effective transgenesis method especially those with single-copy insertion is important for functional comparison. Here, we improved a transposon-based transgenesis methodology that had been originally developed in C. elegans but worked marginally in C. briggsae. By incorporation of a heat shock step, the transgenesis efficiency in C. briggsae with a single-copy insertion is comparable to that in C. elegans. We used the method to generate 54 independent insertions mostly consisting of a mCherry tag over the C. briggsae genome. We demonstrated the use of the tags in identifying interacting loci responsible for hybrid male sterility between C. briggsae and Caenorhabditis nigoni when combined with the GFP tags we generated previously. Finally, we demonstrated that C. briggsae tolerates the C. elegans toxin, PEEL-1, but not SUP-35, making the latter a potential negative selection marker against extrachromosomal array.

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.

Genomic architecture of 5S rDNA cluster and its variations within and between species

BACKGROUND

Ribosomal DNAs (rDNAs) are arranged in purely tandem repeats, preventing them from being reliably assembled onto chromosomes during generation of genome assembly. The uncertainty of rDNA genomic structure presents a significant barrier for studying their function and evolution.

RESULTS

Here we generate ultra-long Oxford Nanopore Technologies (ONT) and short NGS reads to delineate the architecture and variation of the 5S rDNA cluster in the different strains of C. elegans and C. briggsae. We classify the individual rDNA's repeating units into 25 types based on the unique sequence variations in each unit of C. elegans (N2). We next perform assembly of the cluster by taking advantage of the long reads that carry these units, which led to an assembly of 5S rDNA cluster consisting of up to 167 consecutive 5S rDNA units in the N2 strain. The ordering and copy number of various rDNA units are consistent with the separation time between strains. Surprisingly, we observed a drastically reduced level of variation in the unit composition in the 5S rDNA cluster in the C. elegans CB4856 and C. briggsae AF16 strains than in the C. elegans N2 strain, suggesting that N2, a widely used reference strain, is likely to be defective in maintaining the 5S rDNA cluster stability compared with other wild isolates of C. elegans or C. briggsae.

CONCLUSIONS

The results demonstrate that Nanopore DNA sequencing reads are capable of generating assembly of highly repetitive sequences, and rDNA units are highly dynamic both within and between population(s) of the same species in terms of sequence and copy number. The detailed structure and variation of the 5S rDNA units within the rDNA cluster pave the way for functional and evolutionary studies.

Specific Interactions Between Autosome and X Chromosomes Cause Hybrid Male Sterility in Caenorhabditis Species

Hybrid male progeny from interspecies crosses are more prone to sterility or inviability than hybrid female progeny, and the male sterility and inviability often demonstrate parent-of-origin asymmetry. However, the underlying genetic mechanism of asymmetric sterility or inviability remains elusive. We previously established a genome-wide hybrid incompatibility (HI) landscape between Caenorhabditis briggsae and C. nigoni by phenotyping a large collection of C. nigoni strains each carrying a C. briggsae introgression. In this study, we systematically dissect the genetic mechanism of asymmetric sterility and inviability in both hybrid male and female progeny between the two species. Specifically, we performed reciprocal crosses between C . briggsae and different C. nigoni strains that each carry a GFP-labeled C. briggsae genomic fragment referred to as introgression, and scored the HI phenotypes in the F1 progeny. The aggregated introgressions cover 94.6% of the C. briggsae genome, including 100% of the X chromosome. Surprisingly, we observed that two C. briggsae X fragments that produce C. nigoni male sterility as an introgression rescued hybrid F1 sterility in males fathered by C. briggsae Subsequent backcrossing analyses indicated that a specific interaction between the X-linked interaction and one autosome introgression is required to rescue the hybrid male sterility. In addition, we identified another two C. briggsae genomic intervals on chromosomes II and IV that can rescue the inviability, but not the sterility, of hybrid F1 males fathered by C. nigoni, suggesting the involvement of differential epistatic interactions in the asymmetric hybrid male fertility and inviability. Importantly, backcrossing of the rescued sterile males with C. nigoni led to the isolation of a 1.1-Mb genomic interval that specifically interacts with an X-linked introgression, which is essential for hybrid male fertility. We further identified three C. briggsae genomic intervals on chromosome I, II, and III that produced inviability in all F1 progeny, dependent on or independent of the parent-of-origin. Taken together, we identified multiple independent interacting loci that are responsible for asymmetric hybrid male and female sterility, and inviability, which lays a foundation for their molecular characterization.

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.

Genomic basis of recombination suppression in the hybrid between Caenorhabditis briggsae and C. nigoni

DNA recombination is required for effective segregation and diversification of genomes and for the successful completion of meiosis. Recent studies in various species hybrids have demonstrated a genetic link between DNA recombination and speciation. Consistent with this, we observed a striking suppression of recombination in the hybrids between two nematodes, the hermaphroditic Caenorhabditis briggsae and the gonochoristic C. nigoni. To unravel the molecular basis underlying the recombination suppression in their hybrids, we generated a C. nigoni genome with chromosome-level contiguity and produced an improved C. briggsae genome with resolved gaps up to 2.8 Mb. The genome alignment reveals not only high sequence divergences but also pervasive intra- and inter-chromosomal sequence re-arrangements between the two species, which are plausible culprits for the observed suppression. Comparison of recombination boundary sequences suggests that recombination in the hybrid requires extensive sequence homology, which is rarely seen between the two genomes. The new genomes and genomic libraries form invaluable resources for studying genome evolution, hybrid incompatibilities and sex evolution for this pair of model species.

Rapid genome shrinkage in a self-fertile nematode reveals sperm competition proteins

To reveal impacts of sexual mode on genome content, we compared chromosome-scale assemblies of the outcrossing nematode Caenorhabditis nigoni to its self-fertile sibling species, C. briggsaeC. nigoni's genome resembles that of outcrossing relatives but encodes 31% more protein-coding genes than C. briggsaeC. nigoni genes lacking C. briggsae orthologs were disproportionately small and male-biased in expression. These include the male secreted short (mss) gene family, which encodes sperm surface glycoproteins conserved only in outcrossing species. Sperm from mss-null males of outcrossing C. remanei failed to compete with wild-type sperm, despite normal fertility in noncompetitive mating. Restoring mss to C. briggsae males was sufficient to enhance sperm competitiveness. Thus, sex has a pervasive influence on genome content that can be used to identify sperm competition factors.

Comparative proteome analysis between C . briggsae embryos and larvae reveals a role of chromatin modification proteins in embryonic cell division

Caenorhabditis briggsae has emerged as a model for comparative biology against model organism C. elegans. Most of its cell fate specifications are completed during embryogenesis whereas its cell growth is achieved mainly in larval stages. The molecular mechanism underlying the drastic developmental changes is poorly understood. To gain insights into the molecular changes between the two stages, we compared the proteomes between the two stages using iTRAQ. We identified a total of 2,791 proteins in the C. briggsae embryos and larvae, 247 of which undergo up- or down-regulation between the two stages. The proteins that are upregulated in the larval stages are enriched in the Gene Ontology categories of energy production, protein translation, and cytoskeleton; whereas those upregulated in the embryonic stage are enriched in the categories of chromatin dynamics and posttranslational modification, suggesting a more active chromatin modification in the embryos than in the larva. Perturbation of a subset of chromatin modifiers followed by cell lineage analysis suggests their roles in controlling cell division pace. Taken together, we demonstrate a general molecular switch from chromatin modification to metabolism during the transition from C. briggsae embryonic to its larval stages using iTRAQ approach. The switch might be conserved across metazoans.

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.

Specific down-regulation of spermatogenesis genes targeted by 22G RNAs in hybrid sterile males associated with an X-Chromosome introgression

Hybrid incompatibility (HI) prevents gene flow between species, thus lying at the heart of speciation genetics. One of the most common HIs is male sterility. Two superficially contradictory observations exist for hybrid male sterility. First, an introgression on the X Chromosome is more likely to produce male sterility than on autosome (so-called large-X theory); second, spermatogenesis genes are enriched on the autosomes but depleted on the X Chromosome (demasculinization of X Chromosome). Analysis of gene expression in Drosophila hybrids suggests a genetic interaction between the X Chromosome and autosomes that is essential for male fertility. However, the prevalence of such an interaction and its underlying mechanism remain largely unknown. Here we examine the interaction in nematode species by contrasting the expression of both coding genes and transposable elements (TEs) between hybrid sterile males and its parental nematode males. We use two lines of hybrid sterile males, each carrying an independent introgression fragment from Caenorhabditis briggsae X Chromosome in an otherwise Caenorhabditis nigoni background, which demonstrate similar defects in spermatogenesis. We observe a similar pattern of down-regulated genes that are specific for spermatogenesis between the two hybrids. Importantly, the down-regulated genes caused by the X Chromosome introgressions show a significant enrichment on the autosomes, supporting an epistatic interaction between the X Chromosome and autosomes. We investigate the underlying mechanism of the interaction by measuring small RNAs and find that a subset of 22G RNAs specifically targeting the down-regulated spermatogenesis genes is significantly up-regulated in hybrids, suggesting that perturbation of small RNA-mediated regulation may contribute to the X-autosome interaction.

Suppression of F1 Male-Specific Lethality in Caenorhabditis Hybrids by cbr-him-8

Haldane's Rule and Darwin's Corollary to Haldane's Rule are the observations that heterogametic F1 hybrids are frequently less fit than their homogametic siblings, and that asymmetric results are often obtained from reciprocal hybrid crosses. In Caenorhabditis, Haldane's Rule and Darwin's Corollary have been observed in several hybrid crosses, including crosses of Caenorhabditis briggsae and C. nigoni. Fertile F1 females are obtained from reciprocal crosses. However, F1 males obtained from C. nigoni mothers are sterile and F1 males obtained from C. briggsae die during embryogenesis. We have identified cbr-him-8 as a recessive maternal-effect suppressor of F1 hybrid male-specific lethality in this combination of species. This result implicates epigenetic meiotic silencing in the suppression of F1 male-specific lethality. It is also shown that F1 males bearing a C. briggsae X chromosome are fertile. When crossed to C. briggsae hermaphrodites or F1 females derived from C. briggsae hermaphrodites, viable F2 and backcross (B2) progeny were obtained. Sibling males that possessed a C. nigoni X chromosome were sterile. Therefore, the sterility of F1 males bearing a C. nigoni X chromosome must result from dysgenic interactions between the X chromosome of C. nigoni and the autosomes of C. briggsae. The fertility of F1 males bearing a C. briggsae X chromosome provides an opportunity to identify C. nigoni loci that prevent spermatogenesis, and hence hermaphroditic reproduction, in diplo-X hybrids.

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.

Illumina Synthetic Long Read Sequencing Allows Recovery of Missing Sequences even in the "Finished" C. elegans Genome

Most next-generation sequencing platforms permit acquisition of high-throughput DNA sequences, but the relatively short read length limits their use in genome assembly or finishing. Illumina has recently released a technology called Synthetic Long-Read Sequencing that can produce reads of unusual length, i.e., predominately around 10 Kb. However, a systematic assessment of their use in genome finishing and assembly is still lacking. We evaluate the promise and deficiency of the long reads in these aspects using isogenic C. elegans genome with no gap. First, the reads are highly accurate and capable of recovering most types of repetitive sequences. However, the presence of tandem repetitive sequences prevents pre-assembly of long reads in the relevant genomic region. Second, the reads are able to reliably detect missing but not extra sequences in the C. elegans genome. Third, the reads of smaller size are more capable of recovering repetitive sequences than those of bigger size. Fourth, at least 40 Kbp missing genomic sequences are recovered in the C. elegans genome using the long reads. Finally, an N50 contig size of at least 86 Kbp can be achieved with 24 × reads but with substantial mis-assembly errors, highlighting a need for novel assembly algorithm for the long reads.

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.

Gonad morphogenesis defects drive hybrid male sterility in asymmetric hybrid breakdown of Caenorhabditis nematodes

Determining the causes and evolution of reproductive barriers to gene flow between populations, speciation, is the key to understanding the origin of diversity in nature. Many species manifest hybrid breakdown when they intercross, characterized by increasingly exacerbated problems in later generations of hybrids. Recently, Caenorhabditis nematodes have emerged as a genetic model for studying speciation, and here we investigate the nature and causes of hybrid breakdown between Caenorhabditis remanei and C. latens. We quantify partial F1 hybrid inviability and extensive F2 hybrid inviability; the ~75% F2 embryonic arrest occurs primarily during gastrulation or embryonic elongation. Moreover, F1 hybrid males exhibit Haldane's rule asymmetrically for both sterility and inviability, being strongest when C. remanei serves as maternal parent. We show that the mechanism by which sterile hybrid males are incapable of transferring sperm or a copulatory plug involves defective gonad morphogenesis, which we hypothesize results from linker cell defects in migration and/or cell death during development. This first documented case of partial hybrid male sterility in Caenorhabditis follows expectations of Darwin's corollary to Haldane's rule for asymmetric male fitness, providing a powerful foundation for molecular dissection of intrinsic reproductive barriers and divergence of genetic pathways controlling organ morphogenesis.

Worm variation made accessible: Take your shopping cart to store, link, and investigate!

In Caenorhabditis elegans, the recent advances in high-throughput quantitative analyses of natural genetic and phenotypic variation have led to a wealth of data on genotype phenotype relations. This data has resulted in the discovery of genes with major allelic effects and insights in the effect of natural genetic variation on a whole range of complex traits as well as how this variation is distributed across the genome. Regardless of the advances presented in specific studies, the majority of the data generated in these studies had yet to be made easily accessible, allowing for meta-analysis. Not only data in figures or tables but meta-data should be accessible for further investigation and comparison between studies. A platform was created where all the data, phenotypic measurements, genotypes, and mappings can be stored, compared, and new linkages within and between published studies can be discovered. WormQTL focuses on quantitative genetics in Caenorhabditis and other nematode species, whereas WormQTL^HD quantitatively links gene expression quantitative trait loci (eQTL) in C. elegans to gene–disease associations in humans.

Mainstreaming Caenorhabditis elegans in experimental evolution

Experimental evolution provides a powerful manipulative tool for probing evolutionary process and mechanism. As this approach to hypothesis testing has taken purchase in biology, so too has the number of experimental systems that use it, each with its own unique strengths and weaknesses. The depth of biological knowledge about Caenorhabditis nematodes, combined with their laboratory tractability, positions them well for exploiting experimental evolution in animal systems to understand deep questions in evolution and ecology, as well as in molecular genetics and systems biology. To date, Caenorhabditis elegans and related species have proved themselves in experimental evolution studies of the process of mutation, host-pathogen coevolution, mating system evolution and life-history theory. Yet these organisms are not broadly recognized for their utility for evolution experiments and remain underexploited. Here, we outline this experimental evolution work undertaken so far in Caenorhabditis, detail simple methodological tricks that can be exploited and identify research areas that are ripe for future discovery.