Comparison of C. elegans and C. briggsae genome sequences reveals extensive conservation of chromosome organization and synteny

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

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

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.

A novel method of multiplex SNP genotyping assay through variable fragment length allele-specific polymerase chain reaction: Multiplex VFLASP-ARMS

Variable Fragment Length Allele-Specific Polymerase Chain Reaction (VFLASP) and Amplification Refractory Mutation System (ARMS) are reliable methods for detecting allelic variations resulting from single base changes within the genome. Due to their widespread application, allele variations caused by Single Nucleotide Polymorphisms (SNPs) can be readily detected using allele-specific primers. In the context of the current study, VFLASP was combined with ARMS method as a novel strategy to enhance the efficacy of both techniques. Clinically important base variations within SNP regions used in the study were detected by a fragment analysis method. To validate the accuracy of the developed VFLASP-ARMS method, specifically designed synthetic sequences were tested using a capillary electrophoresis system. Allele-specific primers exhibit differences solely at the 3' end based on the sequence of the SNP. Additionally, to increase the specificity of the primers, a base was intentionally added for incompatibility. Therefore, allele discrimination on fragment analysis has been made possible through the 3-6 bp differences in the amplicons. With the optimization of the system, designed synthetic sequences provided reliable and reproducible results in wild-type, heterozygous, and homozygous genotypes using the VFLASP-ARMS method. Hence, our results demonstrated that VFLASP-ARMS method, offers a novel design methodology that can be included in the content of SNP genotyping assays.

Gene duplication and evolutionary plasticity of lin-12/Notch gene function in Caenorhabditis

Gene duplication is an important substrate for the evolution of new gene functions, but the impacts of gene duplicates on their own activities and on the developmental networks in which they act are poorly understood. Here, we use a natural experiment of lin-12/Notch gene duplication within the nematode genus Caenorhabditis, combined with characterization of loss- and gain-of-function mutations, to uncover functional distinctions between the duplicate genes in 1 species (Caenorhabditis briggsae) and their single-copy ortholog in Caenorhabditis elegans. First, using improved genomic sequence and gene model characterization, we confirm that the C. briggsae genome includes 2 complete lin-12 genes, whereas most other genes encoding proteins that participate in the LIN-12 signaling pathway retain a one-to-one orthology with C. elegans. We use CRISPR-mediated genome editing to introduce alleles predicted to cause gain-of-function (gf) or loss-of-function (lf) into each C. briggsae gene and find that the gf mutations uncover functional distinctions not apparent from the lf alleles. Specifically, Cbr-lin-12.1(gf), but not Cbr-lin-12.2(gf), causes developmental defects similar to those observed in Cel-lin-12(gf). In contrast to Cel-lin-12(gf), however, the Cbr-lin-12.1(gf) alleles do not cause dominant phenotypes as compared to the wild type, and the mutant phenotype is observed only when 2 gf alleles are present. Our results demonstrate that gene duplicates can exhibit differential capacities to compensate for each other and to interfere with normal development, and uncover coincident gene duplication and evolution of developmental sensitivity to LIN-12/Notch activity.

Syntenic relationship of chromosomes in Strongyloides species and Rhabditophanes diutinus based on the chromosome-level genome assemblies

The Strongyloides clade, to which the parasitic nematode genus Strongyloides belongs, contains taxa with diverse lifestyles, ranging from free-living to obligate vertebrate parasites. Reproductive strategies are also diverse in this group of nematodes, employing not only sexual reproduction but also parthenogenesis, making it an attractive group to study genome adaptation to specific conditions. An in-depth understanding of genome evolution, however, has been hampered by fragmented genome assemblies. In this study, we generated chromosome-level genome assemblies for two Strongyloides species and the outgroup species Rhabditophanes diutinus using long-read sequencing and high-throughput chromosome conformation capture (Hi-C). Our synteny analyses revealed a clearer picture of chromosome evolution in this group, suggesting that a functional sex chromosome has been maintained throughout the group. We further investigated sex chromosome dynamics in the lifecycle of Strongyloides ratti and found that bivalent formation in oocytes appears to be important for male production in the mitotic parthenogenesis. This article is part of the Theo Murphy meeting issue 'Strongyloides: omics to worm-free populations'.

CaeNDR, the Caenorhabditis Natural Diversity Resource

Studies of model organisms have provided important insights into how natural genetic differences shape trait variation. These discoveries are driven by the growing availability of genomes and the expansive experimental toolkits afforded to researchers using these species. For example, Caenorhabditis elegans is increasingly being used to identify and measure the effects of natural genetic variants on traits using quantitative genetics. Since 2016, the C. elegans Natural Diversity Resource (CeNDR) has facilitated many of these studies by providing an archive of wild strains, genome-wide sequence and variant data for each strain, and a genome-wide association (GWA) mapping portal for the C. elegans community. Here, we present an updated platform, the Caenorhabditis Natural Diversity Resource (CaeNDR), that enables quantitative genetics and genomics studies across the three Caenorhabditis species: C. elegans, C. briggsae and C. tropicalis. The CaeNDR platform hosts several databases that are continually updated by the addition of new strains, whole-genome sequence data and annotated variants. Additionally, CaeNDR provides new interactive tools to explore natural variation and enable GWA mappings. All CaeNDR data and tools are accessible through a freely available web portal located at caendr.org.

Genomic diversity landscapes in outcrossing and selfing Caenorhabditis nematodes

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

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.

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

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

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.

Genome annotation of Caenorhabditis briggsae by TEC-RED identifies new exons, paralogs, and conserved and novel operons

The nematode Caenorhabditis briggsae is routinely used in comparative and evolutionary studies involving its well-known cousin Caenorhabditis elegans. The C. briggsae genome sequence has accelerated research by facilitating the generation of new resources, tools, and functional studies of genes. While substantial progress has been made in predicting genes and start sites, experimental evidence is still lacking in many cases. Here, we report an improved annotation of the C. briggsae genome using the trans-spliced exon coupled RNA end determination technique. In addition to identifying the 5′ ends of expressed genes, we have discovered operons and paralogs. In summary, our analysis yielded 10,243 unique 5′ end sequence tags with matches in the C. briggsae genome. Of these, 6,395 were found to represent 4,252 unique genes along with 362 paralogs and 52 previously unknown exons. These genes included 14 that are exclusively trans-spliced in C. briggsae when compared with C. elegans orthologs. A major contribution of this study is the identification of 492 high confidence operons, of which two-thirds are fully supported by tags. In addition, 2 SL1-type operons were discovered. Interestingly, comparisons with C. elegans showed that only 40% of operons are conserved. Of the remaining operons, 73 are novel, including 12 that entirely lack orthologs in C. elegans. Further analysis revealed that 4 of the 12 novel operons are conserved in Caenorhabditis nigoni. Altogether, the work described here has significantly advanced our understanding of the C. briggsae system and serves as a rich resource to aid biological studies involving this species.

The anterior Hox gene ceh-13 and elt-1/GATA activate the posterior Hox genes nob-1 and php-3 to specify posterior lineages in the C. elegans embryo

Hox transcription factors play a conserved role in specifying positional identity during animal development, with posterior Hox genes typically repressing the expression of more anterior Hox genes. Here, we dissect the regulation of the posterior Hox genes nob-1 and php-3 in the nematode C. elegans. We show that nob-1 and php-3 are co-expressed in gastrulation-stage embryos in cells that previously expressed the anterior Hox gene ceh-13. This expression is controlled by several partially redundant transcriptional enhancers. These enhancers act in a ceh-13-dependant manner, providing a striking example of an anterior Hox gene positively regulating a posterior Hox gene. Several other regulators also act positively through nob-1/php-3 enhancers, including elt-1/GATA, ceh-20/ceh-40/Pbx, unc-62/Meis, pop-1/TCF, ceh-36/Otx, and unc-30/Pitx. We identified defects in both cell position and cell division patterns in ceh-13 and nob-1;php-3 mutants, suggesting that these factors regulate lineage identity in addition to positional identity. Together, our results highlight the complexity and flexibility of Hox gene regulation and function and the ability of developmental transcription factors to regulate different targets in different stages of development.

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

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

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.

Nematode chromosomes

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

Formation and function of dauer ascarosides in the nematodes Caenorhabditis briggsae and Caenorhabditis elegans

The biosynthetic pathways and functions of ascaroside signaling molecules in the nematode Caenorhabditis elegans have been studied to better understand complex, integrative developmental decision-making. Although it is known that ascarosides play multiple roles in the development and behavior of nematode species other than C. elegans, these parallel pheromone systems have not been well-studied. Here, we show that ascarosides in the nematode Caenorhabditis briggsae are biosynthesized in the same manner as C. elegans and act to induce the alternative developmental pathway that generates the stress-resistant dauer lifestage. We show that ascr#2 is the primary component of crude dauer pheromone in C. briggsae; in contrast, C. elegans dauer pheromone relies on a combination of ascr#2, ascr#3, and several other components. We further demonstrate that Cbr-daf-22, like its C. elegans ortholog Cel-daf-22, is necessary to produce short-chain ascarosides. Moreover, Cbr-daf-22 and Cel-daf-22 mutants produce an ascaroside-independent metabolite that acts antagonistically to crude dauer pheromone and inhibits dauer formation.

Deeply conserved synteny and the evolution of metazoan chromosomes

Animal genomes show networks of deeply conserved gene linkages whose phylogenetic scope and chromosomal context remain unclear. Here, we report chromosome-scale conservation of synteny among bilaterians, cnidarians, and sponges and use comparative analysis to reconstruct ancestral chromosomes across major animal groups. Comparisons among diverse metazoans reveal the processes of chromosome evolution that produced contemporary karyotypes from their Precambrian progenitors. On the basis of these findings, we introduce a simple algebraic representation of chromosomal change and use it to establish a unified systematic framework for metazoan chromosome evolution. We find that fusion-with-mixing, a previously unappreciated mode of chromosome change, has played a central role. We find that relicts of several metazoan chromosomal units are preserved in unicellular eukaryotes. These conserved pre-metazoan linkages include the chromosomal unit that encodes the most diverse set of metazoan homeobox genes, suggesting a candidate genomic context for the early diversification of this key gene family.

Modulation of sensory perception by hydrogen peroxide enables Caenorhabditis elegans to find a niche that provides both food and protection from hydrogen peroxide

Hydrogen peroxide (H2O2) is the most common chemical threat that organisms face. Here, we show that H2O2 alters the bacterial food preference of Caenorhabditis elegans, enabling the nematodes to find a safe environment with food. H2O2 induces the nematodes to leave food patches of laboratory and microbiome bacteria when those bacterial communities have insufficient H2O2-degrading capacity. The nematode's behavior is directed by H2O2-sensing neurons that promote escape from H2O2 and by bacteria-sensing neurons that promote attraction to bacteria. However, the input for H2O2-sensing neurons is removed by bacterial H2O2-degrading enzymes and the bacteria-sensing neurons' perception of bacteria is prevented by H2O2. The resulting cross-attenuation provides a general mechanism that ensures the nematode's behavior is faithful to the lethal threat of hydrogen peroxide, increasing the nematode's chances of finding a niche that provides both food and protection from hydrogen peroxide.

Cold Survival and Its Molecular Mechanisms in a Locally Adapted Nematode Population

Since Darwin, evolutionary biologists have sought to understand the drivers and mechanisms of natural trait diversity. The field advances toward this goal with the discovery of phenotypes that vary in the wild, their relationship to ecology, and their underlying genes. Here, we established resistance to extreme low temperature in the free-living nematode Caenorhabditis briggsae as an ecological and evolutionary model system. We found that C. briggsae strains of temperate origin were strikingly more cold-resistant than those isolated from tropical localities. Transcriptional profiling revealed expression patterns unique to the resistant temperate ecotype, including dozens of genes expressed at high levels even after multiple days of cold-induced physiological slowdown. Mutational analysis validated a role in cold resistance for seven such genes. These findings highlight a candidate case of robust, genetically complex adaptation in an emerging model nematode, and shed light on the mechanisms at play.

The Exploration of Novel Regulatory Relationships Drives Haloarchaeal Operon-Like Structural Dynamics over Short Evolutionary Distances

Operons are a dominant feature of bacterial and archaeal genome organization. Numerous investigations have related aspects of operon structure to operon function, making operons exemplars for studies aimed at deciphering Nature's design principles for genomic organization at a local scale. We consider this understanding to be both fundamentally important and ultimately useful in the de novo design of increasingly complex synthetic circuits. Here we analyze the evolution of the genomic context of operon-like structures in a set of 76 sequenced and annotated species of halophilic archaea. The phylogenetic depth and breadth of this dataset allows insight into changes in operon-like structures over shorter evolutionary time scales than have been studied in previous cross-species analysis of operon evolution. Our analysis, implemented in the updated software package JContextExplorer finds that operon-like context as measured by changes in structure frequently differs from a sequence divergence model of whole-species phylogeny and that changes seem to be dominated by the exploration of novel regulatory relationships.

Signatures of adaptation to a monocot host in the plant-parasitic cyst nematode Heterodera sacchari

Interactions between plant-parasitic nematodes and their hosts are mediated by effectors, i.e. secreted proteins that manipulate the plant to the benefit of the pathogen. To understand the role of effectors in host adaptation in nematodes, we analysed the transcriptome of Heterodera sacchari, a cyst nematode parasite of rice (Oryza sativa) and sugarcane (Saccharum officinarum). A multi-gene phylogenetic analysis showed that H. sacchari and the cereal cyst nematode Heterodera avenae share a common evolutionary origin and that they evolved to parasitise monocot plants from a common dicot-parasitic ancestor. We compared the effector repertoires of H. sacchari with those of the dicot parasites Heterodera glycines and Globodera rostochiensis to understand the consequences of this transition. While, in general, effector repertoires are similar between the species, comparing effectors and non-effectors of H. sacchari and G. rostochiensis shows that effectors have accumulated more mutations than non-effectors. Although most effectors show conserved spatiotemporal expression profiles and likely function, some H. sacchari effectors are adapted to monocots. This is exemplified by the plant-peptide hormone mimics, the CLAVATA3/EMBRYO SURROUNDING REGION-like (CLE) effectors. Peptide hormones encoded by H. sacchari CLE effectors are more similar to those from rice than those from other plants, or those from other plant-parasitic nematodes. We experimentally validated the functional significance of these observations by demonstrating that CLE peptides encoded by H. sacchari induce a short root phenotype in rice, whereas those from a related dicot parasite do not. These data provide a functional example of effector evolution that co-occurred with the transition from a dicot-parasitic to a monocot-parasitic lifestyle.

Development and characterization of genomic resources for a non-model marine teleost, the red snapper (Lutjanus campechanus, Lutjanidae): Construction of a high-density linkage map, anchoring of genome contigs and comparative genomic analysis

The red snapper Lutjanus campechanus is an exploited reef fish of major economic importance in the Gulf of Mexico region. Studies of genome wide genetic variation are needed to understand the structure of wild populations and develop breeding programs for aquaculture but interpretation of these genome scans is limited by the absence of reference genome. In this work, the first draft of a reference genome was developed and characterized for the red snapper. P-454 and Illumina sequencing were conducted to produce paired-end reads that were assembled into reference contigs and scaffolds. The current assembly spans over 770 Mb, representing an estimated 69% of the red snapper genome in 67,254 scaffolds (N50 = 16,803 bp). The genome contigs were applied to map double digest Restriction-Site Associated DNA Tags and characterize Single Nucleotide Polymorphisms (SNPs) in five outbred full-sib families. The identified SNPs and 97 microsatellite loci were used to generate a high-density linkage map that includes 7,420 markers distributed across 24 linkage groups and spans 1,346.64 cM with an average inter–marker distance of 0.18 cM. Sex-specific maps revealed a 1.10:1 female to male map length ratio. A total of 4,422 genome contigs (10.5% of the assembly) were anchored to the map and used in a comparative genomic analysis of the red snapper and two model teleosts. Red snapper showed a high degree of chromosome level syntenic conservation with both medaka and spotted green puffer and a near one to one correspondence between the 24 red snapper linkage groups and corresponding medaka chromosomes was observed. This work established the first draft of a reference genome for a lutjanid fish. The obtained genomic resources will serve as a framework for the interpretation of genome scans during studies of wild populations and captive breeding programs.

Chromosome-Level Assembly of the Caenorhabditis remanei Genome Reveals Conserved Patterns of Nematode Genome Organization

The nematode Caenorhabditis elegans is one of the key model systems in biology, including possessing the first fully assembled animal genome. Whereas C. elegans is a self-reproducing hermaphrodite with fairly limited within-population variation, its relative C. remanei is an outcrossing species with much more extensive genetic variation, making it an ideal parallel model system for evolutionary genetic investigations. Here, we greatly improve on previous assemblies by generating a chromosome-level assembly of the entire C. remanei genome (124.8 Mb of total size) using long-read sequencing and chromatin conformation capture data. Like other fully assembled genomes in the genus, we find that the C. remanei genome displays a high degree of synteny with C. elegans despite multiple within-chromosome rearrangements. Both genomes have high gene density in central regions of chromosomes relative to chromosome ends and the opposite pattern for the accumulation of repetitive elements. C. elegans and C. remanei also show similar patterns of interchromosome interactions, with the central regions of chromosomes appearing to interact with one another more than the distal ends. The new C. remanei genome presented here greatly augments the use of the Caenorhabditis as a platform for comparative genomics and serves as a basis for molecular population genetics within this highly diverse species.

A Nearly Complete Genome of Ciona intestinalis Type A (C. robusta) Reveals the Contribution of Inversion to Chromosomal Evolution in the Genus Ciona

Since its initial publication in 2002, the genome of Ciona intestinalis type A (Ciona robusta), the first genome sequence of an invertebrate chordate, has provided a valuable resource for a wide range of biological studies, including developmental biology, evolutionary biology, and neuroscience. The genome assembly was updated in 2008, and it included 68% of the sequence information in 14 pairs of chromosomes. However, a more contiguous genome is required for analyses of higher order genomic structure and of chromosomal evolution. Here, we provide a new genome assembly for an inbred line of this animal, constructed with short and long sequencing reads and Hi-C data. In this latest assembly, over 95% of the 123 Mb of sequence data was included in the chromosomes. Short sequencing reads predicted a genome size of 114-120 Mb; therefore, it is likely that the current assembly contains almost the entire genome, although this estimate of genome size was smaller than previous estimates. Remapping of the Hi-C data onto the new assembly revealed a large inversion in the genome of the inbred line. Moreover, a comparison of this genome assembly with that of Ciona savignyi, a different species in the same genus, revealed many chromosomal inversions between these two Ciona species, suggesting that such inversions have occurred frequently and have contributed to chromosomal evolution of Ciona species. Thus, the present assembly greatly improves an essential resource for genome-wide studies of ascidians.

The Genome of Caenorhabditis bovis

The free-living nematode Caenorhabditis elegans is a key laboratory model for metazoan biology. C. elegans has also become a model for parasitic nematodes despite being only distantly related to most parasitic species. All of the ∼65 Caenorhabditis species currently in culture are free-living, with most having been isolated from decaying plant or fungal matter. Caenorhabditis bovis is a particularly unusual species that has been isolated several times from the inflamed ears of Zebu cattle in Eastern Africa, where it is associated with the disease bovine parasitic otitis. C. bovis is therefore of particular interest to researchers interested in the evolution of nematode parasitism. However, as C. bovis is not in laboratory culture, it remains little studied. Here, by sampling livestock markets and slaughterhouses in Western Kenya, we successfully reisolated C. bovis from the ear of adult female Zebu. We sequenced the genome of C. bovis using the Oxford Nanopore MinION platform in a nearby field laboratory and used the data to generate a chromosome-scale draft genome sequence. We exploited this draft genome sequence to reconstruct the phylogenetic relationships of C. bovis to other Caenorhabditis species and reveal the changes in genome size and content that have occurred during its evolution. We also identified expansions in several gene families that have been implicated in parasitism in other nematode species. The high-quality draft genome and our analyses thereof represent a significant advancement in our understanding of this unusual Caenorhabditis species.

Males, Outcrossing, and Sexual Selection in Caenorhabditis Nematodes

Males of Caenorhabditis elegans provide a crucial practical tool in the laboratory, but, as the rarer and more finicky sex, have not enjoyed the same depth of research attention as hermaphrodites. Males, however, have attracted the attention of evolutionary biologists who are exploiting the C. elegans system to test longstanding hypotheses about sexual selection, sexual conflict, transitions in reproductive mode, and genome evolution, as well as to make new discoveries about Caenorhabditis organismal biology. Here, we review the evolutionary concepts and data informed by study of males of C. elegans and other Caenorhabditis We give special attention to the important role of sperm cells as a mediator of inter-male competition and male-female conflict that has led to drastic trait divergence across species, despite exceptional phenotypic conservation in many other morphological features. We discuss the evolutionary forces important in the origins of reproductive mode transitions from males being common (gonochorism: females and males) to rare (androdioecy: hermaphrodites and males) and the factors that modulate male frequency in extant androdioecious populations, including the potential influence of selective interference, host-pathogen coevolution, and mutation accumulation. Further, we summarize the consequences of males being common vs rare for adaptation and for trait divergence, trait degradation, and trait dimorphism between the sexes, as well as for molecular evolution of the genome, at both micro-evolutionary and macro-evolutionary timescales. We conclude that C. elegans male biology remains underexploited and that future studies leveraging its extensive experimental resources are poised to discover novel biology and to inform profound questions about animal function and evolution.

Hybrid Assembly of the Genome of the Entomopathogenic Nematode Steinernema carpocapsae Identifies the X-Chromosome

Entomopathogenic nematodes from the genus Steinernema are lethal insect parasites that quickly kill their insect hosts with the help of their symbiotic bacteria. Steinernema carpocapsae is one of the most studied entomopathogens due to its broad lethality to diverse insect species and its effective commercial use as a biological control agent for insect pests, as well as a genetic model for studying parasitism, pathogenesis, and symbiosis. In this study, we used long-reads from the Pacific Biosciences platform and BioNano Genomics Irys system to assemble the most complete genome of the S. carpocapsae ALL strain to date, comprising 84.5 Mb in 16 scaffolds, with an N50 of 7.36 Mb. The largest scaffold, with 20.9 Mb, was identified as chromosome X based on sex-specific genome sequencing. The high level of contiguity allowed us to characterize gene density, repeat content, and GC content. RNA-seq data from 17 developmental stages, spanning from embryo to adult, were used to predict 30,957 gene models. Using this improved genome, we performed a macrosyntenic analysis to Caenorhabditis elegans and Pristionchus pacificus and found S. carpocapsae's chromosome X to be primarily orthologous to C. elegans' and P. pacificus' chromosome II and IV. We also investigated the expansion of protein families and gene expression differences between adult male and female stage nematodes. This new genome and more accurate set of annotations provide a foundation for additional comparative genomic and gene expression studies within the Steinernema clade and across the Nematoda phylum.

Genome structure predicts modular transcriptome responses to genetic and environmental conditions

Understanding the plasticity, robustness and modularity of transcriptome expression to genetic and environmental conditions is crucial to deciphering how organisms adapt in nature. To test how genome architecture influences transcriptome profiles, we quantified expression responses for distinct temperature-adapted genotypes of the nematode Caenorhabditis briggsae when exposed to chronic temperature stresses throughout development. We found that 56% of the 8,795 differentially expressed genes show genotype-specific changes in expression in response to temperature (genotype-by-environment interactions, GxE). Most genotype-specific responses occur under heat stress, indicating that cold vs. heat stress responses involve distinct genomic architectures. The 22 co-expression modules that we identified differ in their enrichment of genes with genetic vs. environmental vs. interaction effects, as well as their genomic spatial distributions, functional attributes and rates of molecular evolution at the sequence level. Genes in modules enriched for simple effects of either genotype or temperature alone tend to evolve especially rapidly, consistent with disproportionate influence of adaptation or weaker constraint on these subsets of loci. Chromosome-scale heterogeneity in nucleotide polymorphism, however, rather than the scale of individual genes predominates as the source of genetic differences among expression profiles, and natural selection regimes are largely decoupled between coding sequences and noncoding flanking sequences that contain cis-regulatory elements. These results illustrate how the form of transcriptome modularity and genome structure contribute to predictable profiles of evolutionary change.

Mutational Landscape of Spontaneous Base Substitutions and Small Indels in Experimental Caenorhabditis elegans Populations of Differing Size

Experimental investigations into the rates and fitness effects of spontaneous mutations are fundamental to our understanding of the evolutionary process. To gain insights into the molecular and fitness consequences of spontaneous mutations, we conducted a mutation accumulation (MA) experiment at varying population sizes in the nematode Caenorhabditis elegans, evolving 35 lines in parallel for 409 generations at three population sizes (N = 1, 10, and 100 individuals). Here, we focus on nuclear SNPs and small insertion/deletions (indels) under minimal influence of selection, as well as their accrual rates in larger populations under greater selection efficacy. The spontaneous rates of base substitutions and small indels are 1.84 (95% C.I. ± 0.14) × 10-9 substitutions and 6.84 (95% C.I. ± 0.97) × 10-10 changes/site/generation, respectively. Small indels exhibit a deletion bias with deletions exceeding insertions by threefold. Notably, there was no correlation between the frequency of base substitutions, nonsynonymous substitutions, or small indels with population size. These results contrast with our previous analysis of mitochondrial DNA mutations and nuclear copy-number changes in these MA lines, and suggest that nuclear base substitutions and small indels are under less stringent purifying selection compared to the former mutational classes. A transition bias was observed in exons as was a near universal base substitution bias toward A/T. Strongly context-dependent base substitutions, where 5'-Ts and 3'-As increase the frequency of A/T → T/A transversions, especially at the boundaries of A or T homopolymeric runs, manifest as higher mutation rates in (i) introns and intergenic regions relative to exons, (ii) chromosomal cores vs. arms and tips, and (iii) germline-expressed genes.

Centromere and Pericentromere Transcription: Roles and Regulation … in Sickness and in Health

The chromosomal loci known as centromeres (CEN) mediate the equal distribution of the duplicated genome between both daughter cells. Specifically, centromeres recruit a protein complex named the kinetochore, that bi-orients the replicated chromosome pairs to the mitotic or meiotic spindle structure. The paired chromosomes are then separated, and the individual chromosomes segregate in opposite direction along the regressing spindle into each daughter cell. Erroneous kinetochore assembly or activity produces aneuploid cells that contain an abnormal number of chromosomes. Aneuploidy may incite cell death, developmental defects (including genetic syndromes), and cancer (>90% of all cancer cells are aneuploid). While kinetochores and their activities have been preserved through evolution, the CEN DNA sequences have not. Hence, to be recognized as sites for kinetochore assembly, CEN display conserved structural themes. In addition, CEN nucleosomes enclose a CEN-exclusive variant of histone H3, named CENP-A, and carry distinct epigenetic labels on CENP-A and the other CEN histone proteins. Through the cell cycle, CEN are transcribed into non-coding RNAs. After subsequent processing, they become key components of the CEN chromatin by marking the CEN locus and by stably anchoring the CEN-binding kinetochore proteins. CEN transcription is tightly regulated, of low intensity, and essential for differentiation and development. Under- or overexpression of CEN transcripts, as documented for myriad cancers, provoke chromosome missegregation and aneuploidy. CEN are genetically stable and fully competent only when they are insulated from the surrounding, pericentromeric chromatin, which must be silenced. We will review CEN transcription and its contribution to faithful kinetochore function. We will further discuss how pericentromeric chromatin is silenced by RNA processing and transcriptionally repressive chromatin marks. We will report on the transcriptional misregulation of (peri)centromeres during stress, natural aging, and disease and reflect on whether their transcripts can serve as future diagnostic tools and anti-cancer targets in the clinic.

Inverted repeat structures are associated with essential and highly expressed genes on C. elegans autosome distal arms

Complementary sequences in cellular transcripts base-pair to form double-stranded RNA (dsRNA) structures. Because transposon-derived repeats often give rise to self-complementary sequences, dsRNA structures are prevalent in eukaryotic genomes, typically occurring in gene introns and untranslated regions (UTRs). However, the regulatory impact of double-stranded structures within genes is not fully understood. We used three independent methods to define loci in Caenorhabditis elegans predicted to form dsRNA and correlated these structures with patterns of gene expression, gene essentiality, and genome organization. As previously observed, dsRNA loci are enriched on distal arms of C. elegans autosomes, where genes typically show less conservation and lower overall expression. In contrast, we find that dsRNAs are associated with essential genes on autosome arms, and dsRNA-associated genes exhibit higher-than-expected expression and histone modification patterns associated with transcriptional elongation. Genes with significant repetitive sequence content are also highly expressed, and, thus, observed gene expression trends may relate either to dsRNA structures or to repeat content. Our results raise the possibility that as-yet-undescribed mechanisms promote expression of loci that produce dsRNAs, despite their well-characterized roles in gene silencing.

Development of a Simple and Cost-Effective Method Based on T7 Endonuclease Cleavage for Detection of Single Nucleotide Polymorphisms

AIMS

Single nucleotide polymorphisms (SNP) can be used as genetic markers and for risk assessment of allele-linked diseases, which can provide information for clinical diagnosis. Large-scale microarray and next-generation sequencing methods have made genome-wide SNP genotyping possible. However, in addition to their high cost, these techniques are dependent on having specialized equipment. Thus, there is a need for a simple genotyping method that can be implemented in a resource-limited environment.

METHODS

We developed a strategy for SNP genotyping based on T7 Endonuclease I cleavage and an enzyme-linked microparticle immune assay. Using this method, we genotyped two common SNP sites (rs11526468 and rs12979860). The quality of the genotyping process was validated.

RESULTS

Although a 70% false-negative rate was observed, no false-positive reactions were found. Therefore, multiple parallel repeat reactions can offset the possibility of mutation detection failure.

DISCUSSION

This method employs a duplicate reagent-dependent procedure, and therefore has the potential for integration into a portable kit for field utilization.

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.

From "the Worm" to "the Worms" and Back Again: The Evolutionary Developmental Biology of Nematodes

Since the earliest days of research on nematodes, scientists have noted the developmental and morphological variation that exists within and between species. As various cellular and developmental processes were revealed through intense focus on Caenorhabditis elegans, these comparative studies have expanded. Within the genus Caenorhabditis, they include characterization of intraspecific polymorphisms and comparisons of distinct species, all generally amenable to the same laboratory culture methods and supported by robust genomic and experimental tools. The C. elegans paradigm has also motivated studies with more distantly related nematodes and animals. Combined with improved phylogenies, this work has led to important insights about the evolution of nematode development. First, while many aspects of C. elegans development are representative of Caenorhabditis, and of terrestrial nematodes more generally, others vary in ways both obvious and cryptic. Second, the system has revealed several clear examples of developmental flexibility in achieving a particular trait. This includes developmental system drift, in which the developmental control of homologous traits has diverged in different lineages, and cases of convergent evolution. Overall, the wealth of information and experimental techniques developed in C. elegans is being leveraged to make nematodes a powerful system for evolutionary cellular and developmental biology.

Inferring synteny between genome assemblies: a systematic evaluation

BACKGROUND

Genome assemblies across all domains of life are being produced routinely. Initial analysis of a new genome usually includes annotation and comparative genomics. Synteny provides a framework in which conservation of homologous genes and gene order is identified between genomes of different species. The availability of human and mouse genomes paved the way for algorithm development in large-scale synteny mapping, which eventually became an integral part of comparative genomics. Synteny analysis is regularly performed on assembled sequences that are fragmented, neglecting the fact that most methods were developed using complete genomes. It is unknown to what extent draft assemblies lead to errors in such analysis.

RESULTS

We fragmented genome assemblies of model nematodes to various extents and conducted synteny identification and downstream analysis. We first show that synteny between species can be underestimated up to 40% and find disagreements between popular tools that infer synteny blocks. This inconsistency and further demonstration of erroneous gene ontology enrichment tests raise questions about the robustness of previous synteny analysis when gold standard genome sequences remain limited. In addition, assembly scaffolding using a reference guided approach with a closely related species may result in chimeric scaffolds with inflated assembly metrics if a true evolutionary relationship was overlooked. Annotation quality, however, has minimal effect on synteny if the assembled genome is highly contiguous.

CONCLUSIONS

Our results show that a minimum N50 of 1 Mb is required for robust downstream synteny analysis, which emphasizes the importance of gold standard genomes to the science community, and should be achieved given the current progress in sequencing technology.

Evolution of Tandem Repeat Satellite Sequences in Two Closely Related Caenorhabditis Species. Diminution of Satellites in Hermaphrodites

The availability of the genome sequence of the unisexual (male-female) Caenorhabditis nigoni offers an opportunity to compare its non-coding features with the related hermaphroditic species Caenorhabditis briggsae; to understand the evolutionary dynamics of their tandem repeat sequences (satellites), as a result of evolution from the unisexual ancestor. We take advantage of the previously developed SATFIND program to build satellite families defined by a consensus sequence. The relative number of satellites (satellites/Mb) in C. nigoni is 24.6% larger than in C. briggsae. Some satellites in C. nigoni have developed from a proto-repeat present in the ancestor species and are conserved as an isolated sequence in C. briggsae. We also identify unique satellites which occur only once and joint satellite families with a related sequence in both species. Some of these families are only found in C. nigoni, which indicates a recent appearance; they contain conserved adjacent 5′ and 3′ regions, which may favor transposition. Our results show that the number, length and turnover of satellites are restricted in the hermaphrodite C. briggsae when compared with the unisexual C. nigoni. We hypothesize that this results from differences in unequal recombination during meiotic chromosome pairing, which limits satellite turnover in hermaphrodites.

Genome Architecture and Evolution of a Unichromosomal Asexual Nematode

Asexual reproduction in animals, though rare, is the main or exclusive mode of reproduction in some long-lived lineages. The longevity of asexual clades may be correlated with the maintenance of heterozygosity by mechanisms that rearrange genomes and reduce recombination. Asexual species thus provide an opportunity to gain insight into the relationship between molecular changes, genome architecture, and cellular processes. Here we report the genome sequence of the parthenogenetic nematode Diploscapter pachys with only one chromosome pair. We show that this unichromosomal architecture is shared by a long-lived clade of asexual nematodes closely related to the genetic model organism Caenorhabditis elegans. Analysis of the genome assembly reveals that the unitary chromosome arose through fusion of six ancestral chromosomes, with extensive rearrangement among neighboring regions. Typical nematode telomeres and telomeric protection-encoding genes are lacking. Most regions show significant heterozygosity; homozygosity is largely concentrated to one region and attributed to gene conversion. Cell-biological and molecular evidence is consistent with the absence of key features of meiosis I, including synapsis and recombination. We propose that D. pachys preserves heterozygosity and produces diploid embryos without fertilization through a truncated meiosis. As a prelude to functional studies, we demonstrate that D. pachys is amenable to experimental manipulation by RNA interference.

Evolution: A Parthenogenetic Nematode Shows How Animals Become Sexless

Most animals have male and female sexes, implying that sex is ancient and beneficial; yet some have survived for millions of years without sex. The genome of the parthenogenetic nematode Diploscapter pachys gives clues as to how 'ancient asexual' animals can exist.

Trans-splicing enhances translational efficiency in C. elegans

Translational efficiency is subject to extensive regulation. However, the factors influencing such regulation are poorly understood. In Caenorhabditis elegans, 62% of genes are trans-spliced to a specific spliced leader (SL1), which replaces part of the native 5' untranslated region (5' UTR). Given the pivotal role the 5' UTR plays in the regulation of translational efficiency, we hypothesized that SL1 trans-splicing functions to regulate translational efficiency. With genome-wide analysis on Ribo-seq data, polysome profiling experiments, and CRISPR-Cas9-based genetic manipulation of trans-splicing sites, we found four lines of evidence in support of this hypothesis. First, SL1 trans-spliced genes have higher translational efficiencies than non-trans-spliced genes. Second, SL1 trans-spliced genes have higher translational efficiencies than non-trans-spliced orthologous genes in other nematode species. Third, an SL1 trans-spliced isoform has higher translational efficiency than the non-trans-spliced isoform of the same gene. Fourth, deletion of trans-splicing sites of endogenous genes leads to reduced translational efficiency. Importantly, we demonstrated that SL1 trans-splicing plays a key role in enhancing translational efficiencies of essential genes. We further discovered that SL1 trans-splicing likely enhances translational efficiency by shortening the native 5' UTRs, hence reducing the presence of upstream start codons (uAUG) and weakening mRNA secondary structures. Taken together, our study elucidates the global function of trans-splicing in enhancing translational efficiency in nematodes, paving the way for further understanding the genomic mechanisms of translational regulation.

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.

Long-read sequencing improves assembly of Trichinella genomes 10-fold, revealing substantial synteny between lineages diverged over 7 million years

Genome assemblies can form the basis of comparative analyses fostering insight into the evolutionary genetics of a parasite's pathogenicity, host-pathogen interactions, environmental constraints and invasion biology; however, the length and complexity of many parasite genomes has hampered the development of well-resolved assemblies. In order to improve Trichinella genome assemblies, the genome of the sylvatic encapsulated species Trichinella murrelli was sequenced using third-generation, long-read technology and, using syntenic comparisons, scaffolded to a reference genome assembly of Trichinella spiralis, markedly improving both. A high-quality draft assembly for T. murrelli was achieved that totalled 63·2 Mbp, half of which was condensed into 26 contigs each longer than 571 000 bp. When compared with previous assemblies for parasites in the genus, ours required 10-fold fewer contigs, which were five times longer, on average. Better assembly across repetitive regions also enabled resolution of 8 Mbp of previously indeterminate sequence. Furthermore, syntenic comparisons identified widespread scaffold misassemblies in the T. spiralis reference genome. The two new assemblies, organized for the first time into three chromosomal scaffolds, will be valuable resources for future studies linking phenotypic traits within each species to their underlying genetic bases.

The genome of Onchocerca volvulus, agent of river blindness

Human onchocerciasis is a serious neglected tropical disease caused by the filarial nematode Onchocerca volvulus that can lead to blindness and chronic disability. Control of the disease relies largely on mass administration of a single drug, and the development of new drugs and vaccines depends on a better knowledge of parasite biology. Here, we describe the chromosomes of O. volvulus and its Wolbachia endosymbiont. We provide the highest-quality sequence assembly for any parasitic nematode to date, giving a glimpse into the evolution of filarial parasite chromosomes and proteomes. This resource was used to investigate gene families with key functions that could be potentially exploited as targets for future drugs. Using metabolic reconstruction of the nematode and its endosymbiont, we identified enzymes that are likely to be essential for O. volvulus viability. In addition, we have generated a list of proteins that could be targeted by Federal-Drug-Agency-approved but repurposed drugs, providing starting points for anti-onchocerciasis drug development.

To Break or Not To Break: Sex Chromosome Hemizygosity During Meiosis in Caenorhabditis

Meiotic recombination establishes connections between homologous chromosomes to promote segregation. Hemizygous regions of sex chromosomes have no homologous chromosome to recombine with, yet must be transmitted through meiosis. An extreme case of hemizygosity exists in the genus Caenorhabditis, where males have a single X chromosome that completely lacks a homologous partner. To determine whether similar strategies have evolved to accommodate hemizygosity of the X during male meiosis in Caenorhabditis with distinct modes of sexual reproduction, we examined induction and processing of meiotic double strand breaks (DSBs) in androdioecious (hermaphrodite/male) Caenorhabditis elegans and C. briggsae, and gonochoristic (female/male) C. remanei and C. brenneri Analysis of the recombinase RAD-51 suggests more meiotic DSBs are induced in gonochoristic vs. androdioecious species. However, in late prophase in all species, chromosome pairs are restructured into bivalents around a single axis, suggesting that the holocentric nature of Caenorhabditis chromosomes dictates a single crossover per bivalent regardless of the number of DSBs induced. Interestingly, RAD-51 foci were readily observed on the X chromosome of androdioecious male germ cells, while very few were detected in gonochoristic male germ cells. As in C. elegans, the X chromosome in C. briggsae male germ cells undergoes transient pseudosynapsis and flexibility in DSB repair pathway choice. In contrast, in C. remanei and C. brenneri male germ cells, the X chromosome does not undergo pseudosynapsis and appears refractory to SPO-11-induced breaks. Together our results suggest that distinct strategies have evolved to accommodate sex chromosome hemizygosity during meiosis in closely related Caenorhabditis species.

Plasticity in the Meiotic Epigenetic Landscape of Sex Chromosomes in Caenorhabditis Species

During meiosis in the heterogametic sex in some species, sex chromosomes undergo meiotic sex chromosome inactivation (MSCI), which results in acquisition of repressive chromatin and transcriptional silencing. In Caenorhabditis elegans, MSCI is mediated by MET-2 methyltransferase deposition of histone H3 lysine 9 dimethylation. Here we examined the meiotic chromatin landscape in germ lines of four Caenorhabditis species; C. remanei and C. brenneri represent ancestral gonochorism, while C. briggsae and C. elegans are two lineages that independently evolved hermaphroditism. While MSCI is conserved across all four species, repressive chromatin modifications are distinct and do not correlate with reproductive mode. In contrast to C. elegans and C. remanei germ cells where X chromosomes are enriched for histone H3 lysine 9 dimethylation, X chromosomes in C. briggsae and C. brenneri germ cells are enriched for histone H3 lysine 9 trimethylation. Inactivation of C. briggsae MET-2 resulted in germ-line X chromosome transcription and checkpoint activation. Further, both histone H3 lysine 9 di- and trimethylation were reduced in Cbr-met-2 mutant germ lines, suggesting that in contrast to C. elegans, H3 lysine 9 di- and trimethylation are interdependent. C. briggsae H3 lysine 9 trimethylation was redistributed in the presence of asynapsed chromosomes in a sex-specific manner in the related process of meiotic silencing of unsynapsed chromatin. However, these repressive marks did not influence X chromosome replication timing. Examination of additional Caenorhabditis species revealed diverse H3 lysine 9 methylation patterns on the X, suggesting that the sex chromosome epigenome evolves rapidly.

The genomic basis of parasitism in the Strongyloides clade of nematodes

Soil-transmitted nematodes, including the Strongyloides genus, cause one of the most prevalent neglected tropical diseases. Here we compare the genomes of four Strongyloides species, including the human pathogen Strongyloides stercoralis, and their close relatives that are facultatively parasitic (Parastrongyloides trichosuri) and free-living (Rhabditophanes sp. KR3021). A significant paralogous expansion of key gene families--families encoding astacin-like and SCP/TAPS proteins--is associated with the evolution of parasitism in this clade. Exploiting the unique Strongyloides life cycle, we compare the transcriptomes of the parasitic and free-living stages and find that these same gene families are upregulated in the parasitic stages, underscoring their role in nematode parasitism.

Mitogenomics reveals high synteny and long evolutionary histories of sympatric cryptic nematode species

Species with seemingly identical morphology but with distinct genetic differences are abundant in the marine environment and frequently co-occur in the same habitat. Such cryptic species are typically delineated using a limited number of mitochondrial and/or nuclear marker genes, which do not yield information on gene order and gene content of the genomes under consideration. We used next-generation sequencing to study the composition of the mitochondrial genomes of four sympatrically distributed cryptic species of the Litoditis marina species complex (PmI, PmII, PmIII, and PmIV). The ecology, biology, and natural occurrence of these four species are well known, but the evolutionary processes behind this cryptic speciation remain largely unknown. The gene order of the mitochondrial genomes of the four species was conserved, but differences in genome length, gene length, and codon usage were observed. The atp8 gene was lacking in all four species. Phylogenetic analyses confirm that PmI and PmIV are sister species and that PmIII diverged earliest. The most recent common ancestor of the four cryptic species was estimated to have diverged 16 MYA. Synonymous mutations outnumbered nonsynonymous changes in all protein-encoding genes, with the Complex IV genes (coxI-III) experiencing the strongest purifying selection. Our mitogenomic results show that morphologically similar species can have long evolutionary histories and that PmIII has several differences in genetic makeup compared to the three other species, which may explain why it is better adapted to higher temperatures than the other species.

Evolutionary Conservation and Diversification of Puf RNA Binding Proteins and Their mRNA Targets

Reprogramming of a gene’s expression pattern by acquisition and loss of sequences recognized by specific regulatory RNA binding proteins may be a major mechanism in the evolution of biological regulatory programs. We identified that RNA targets of Puf3 orthologs have been conserved over 100–500 million years of evolution in five eukaryotic lineages. Focusing on Puf proteins and their targets across 80 fungi, we constructed a parsimonious model for their evolutionary history. This model entails extensive and coordinated changes in the Puf targets as well as changes in the number of Puf genes and alterations of RNA binding specificity including that: 1) Binding of Puf3 to more than 200 RNAs whose protein products are predominantly involved in the production and organization of mitochondrial complexes predates the origin of budding yeasts and filamentous fungi and was maintained for 500 million years, throughout the evolution of budding yeast. 2) In filamentous fungi, remarkably, more than 150 of the ancestral Puf3 targets were gained by Puf4, with one lineage maintaining both Puf3 and Puf4 as regulators and a sister lineage losing Puf3 as a regulator of these RNAs. The decrease in gene expression of these mRNAs upon deletion of Puf4 in filamentous fungi (N. crassa) in contrast to the increase upon Puf3 deletion in budding yeast (S. cerevisiae) suggests that the output of the RNA regulatory network is different with Puf4 in filamentous fungi than with Puf3 in budding yeast. 3) The coregulated Puf4 target set in filamentous fungi expanded to include mitochondrial genes involved in the tricarboxylic acid (TCA) cycle and other nuclear-encoded RNAs with mitochondrial function not bound by Puf3 in budding yeast, observations that provide additional evidence for substantial rewiring of post-transcriptional regulation. 4) Puf3 also expanded and diversified its targets in filamentous fungi, gaining interactions with the mRNAs encoding the mitochondrial electron transport chain (ETC) complex I as well as hundreds of other mRNAs with nonmitochondrial functions. The many concerted and conserved changes in the RNA targets of Puf proteins strongly support an extensive role of RNA binding proteins in coordinating gene expression, as originally proposed by Keene. Rewiring of Puf-coordinated mRNA targets and transcriptional control of the same genes occurred at different points in evolution, suggesting that there have been distinct adaptations via RNA binding proteins and transcription factors. The changes in Puf targets and in the Puf proteins indicate an integral involvement of RNA binding proteins and their RNA targets in the adaptation, reprogramming, and function of gene expression.

A map of the evolutionary history of Puf proteins and their RNA targets shows that reprogramming of global gene expression programs via adaptive mutations that affect protein-RNA interactions is an important source of biological diversity.

We set out to trace the evolutionary history of an RNA binding protein and how its interactions with targets change over evolution. Identifying this natural history is a step toward understanding the critical differences between organisms and how gene expression programs are rewired during evolution. Using bioinformatics and experimental approaches, we broadly surveyed the evolution of binding targets of a particular family of RNA binding proteins—the Puf proteins, whose protein sequences and target RNA sequences are relatively well-characterized—across 99 eukaryotic species. We found five groups of species in which targets have been conserved for at least 100 million years and then took advantage of genome sequences from a large number of fungal species to deeply investigate the conservation and changes in Puf proteins and their RNA targets. Our analyses identified multiple and extensive reconfigurations during the natural history of fungi and suggest that RNA binding proteins and their RNA targets are profoundly involved in evolutionary reprogramming of gene expression and help define distinct programs unique to each organism. Continuing to uncover the natural history of RNA binding proteins and their interactions will provide a unique window into the gene expression programs of present day species and point to new ways to engineer gene expression programs.

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.

Genome-wide SNP identification for the construction of a high-resolution genetic map of Japanese flounder (Paralichthys olivaceus): applications to QTL mapping of Vibrio anguillarum disease resistance and comparative genomic analysis

High-resolution genetic maps are essential for fine mapping of complex traits, genome assembly, and comparative genomic analysis. Single-nucleotide polymorphisms (SNPs) are the primary molecular markers used for genetic map construction. In this study, we identified 13,362 SNPs evenly distributed across the Japanese flounder (Paralichthys olivaceus) genome. Of these SNPs, 12,712 high-confidence SNPs were subjected to high-throughput genotyping and assigned to 24 consensus linkage groups (LGs). The total length of the genetic linkage map was 3,497.29 cM with an average distance of 0.47 cM between loci, thereby representing the densest genetic map currently reported for Japanese flounder. Nine positive quantitative trait loci (QTLs) forming two main clusters for Vibrio anguillarum disease resistance were detected. All QTLs could explain 5.1-8.38% of the total phenotypic variation. Synteny analysis of the QTL regions on the genome assembly revealed 12 immune-related genes, among them 4 genes strongly associated with V. anguillarum disease resistance. In addition, 246 genome assembly scaffolds with an average size of 21.79 Mb were anchored onto the LGs; these scaffolds, comprising 522.99 Mb, represented 95.78% of assembled genomic sequences. The mapped assembly scaffolds in Japanese flounder were used for genome synteny analyses against zebrafish (Danio rerio) and medaka (Oryzias latipes). Flounder and medaka were found to possess almost one-to-one synteny, whereas flounder and zebrafish exhibited a multi-syntenic correspondence. The newly developed high-resolution genetic map, which will facilitate QTL mapping, scaffold assembly, and genome synteny analysis of Japanese flounder, marks a milestone in the ongoing genome project for this species.