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

Pervasive conservation of intron number and other genetic elements revealed by a chromosome-level genomic assembly of the hyper-polymorphic nematode Caenorhabditis brenneri

With within-species genetic diversity estimates that span the gambit of that seen across the entirety of animals, the Caenorhabditis genus of nematodes holds unique potential to provide insights into how population size and reproductive strategies influence gene and genome organization and evolution. Our study focuses on Caenorhabditis brenneri, currently known as one of the most genetically diverse nematodes within its genus and metazoan phyla. Here, we present a high-quality gapless genome assembly and annotation for C. brenneri, revealing a common nematode chromosome arrangement characterized by gene-dense central regions and repeat rich peripheral parts. Comparison of C. brenneri with other nematodes from the 'Elegans' group revealed conserved macrosynteny but a lack of microsynteny, characterized by frequent rearrangements and low correlation iof orthogroup sizes, indicative of high rates of gene turnover. We also assessed genome organization within corresponding syntenic blocks in selfing and outcrossing species, affirming that selfing species predominantly experience loss of both genes and intergenic DNA. Comparison of gene structures revealed strikingly small number of shared introns across species, yet consistent distributions of intron number and length, regardless of population size or reproductive mode, suggesting that their evolutionary dynamics are primarily reflective of functional constraints. Our study provides valuable insights into genome evolution and expands the nematode genome resources with the highly genetically diverse C. brenneri, facilitating research into various aspects of nematode biology and evolutionary processes.

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.

SMC-mediated dosage compensation in C. elegans evolved in the presence of an ancestral nematode mechanism

Mechanisms of X chromosome dosage compensation have been studied extensively in a few model species representing clades of shared sex chromosome ancestry. However, the diversity within each clade as a function of sex chromosome evolution is largely unknown. Here, we anchor ourselves to the nematode Caenorhabditis elegans, for which a well-studied mechanism of dosage compensation occurs through a specialized structural maintenance of chromosomes (SMC) complex, and explore the diversity of dosage compensation in the surrounding phylogeny of nematodes. Through phylogenetic analysis of the C. elegans dosage compensation complex and a survey of its epigenetic signatures, including X-specific topologically associating domains (TADs) and X-enrichment of H4K20me1, we found that the condensin-mediated mechanism evolved recently in the lineage leading to Caenorhabditis through an SMC-4 duplication. Intriguingly, an independent duplication of SMC-4 and the presence of X-specific TADs in Pristionchus pacificus suggest that condensin-mediated dosage compensation arose more than once. mRNA-seq analyses of gene expression in several nematode species indicate that dosage compensation itself is ancestral, as expected from the ancient XO sex determination system. Indicative of the ancestral mechanism, H4K20me1 is enriched on the X chromosomes in Oscheius tipulae, which does not contain X-specific TADs or SMC-4 paralogs. Together, our results indicate that the dosage compensation system in C. elegans is surprisingly new, and condensin may have been co-opted repeatedly in nematodes, suggesting that the process of evolving a chromosome-wide gene regulatory mechanism for dosage compensation is constrained.

Significance statement

X chromosome dosage compensation mechanisms evolved in response to Y chromosome degeneration during sex chromosome evolution. However, establishment of dosage compensation is not an endpoint. As sex chromosomes change, dosage compensation strategies may have also changed. In this study, we performed phylogenetic and epigenomic analyses surrounding Caenorhabditis elegans and found that the condensin-mediated dosage compensation mechanism in C. elegans is surprisingly new, and has evolved in the presence of an ancestral mechanism. Intriguingly, condensin-based dosage compensation may have evolved more than once in the nematode lineage, the other time in Pristionchus. Together, our work highlights a previously unappreciated diversity of dosage compensation mechanisms within a clade, and suggests constraints in evolving new mechanisms in the presence of an existing one.

Experimental evolution of hybrid populations to identify Dobzhansky-Muller incompatibility loci

Epistatic interactions between loci that reduce fitness in interspecies hybrids, Dobzhansky-Muller incompatibilities (DMIs), contribute genetically to the inviability and infertility within hybrid populations. It remains a challenge, however, to identify the loci that contribute to DMIs as causes of reproductive isolation between species. Here, we assess through forward simulation the power of evolve-and-resequence (E&R) experimental evolution of hybrid populations to map DMI loci. We document conditions under which such a mapping strategy may be most feasible and demonstrate how mapping power is sensitive to biologically relevant parameters such as one-way versus two-way incompatibility type, selection strength, recombination rate, and dominance interactions. We also assess the influence of parameters under direct control of an experimenter, including duration of experimental evolution and number of replicate populations. We conclude that an E&R strategy for mapping DMI loci, and other cases of epistasis, can be a viable option under some circumstances for study systems with short generation times like Caenorhabditis nematodes.

Comparative Genomics of Sex, Chromosomes, and Sex Chromosomes in Caenorhabditis elegans and Other Nematodes

The nematode phylum has evolved a remarkable diversity of reproductive modes, including the repeated emergence of asexuality and hermaphroditism across divergent clades. The species-richness and small genome size of nematodes make them ideal systems for investigating the genome-wide causes and consequences of such major transitions. The availability of functional annotations for most Caenorhabditis elegans genes further allows the linking of patterns of gene content evolution with biological processes. Such gene-centric studies were recently complemented by investigations of chromosome evolution that made use of the first chromosome-scale genome assemblies outside the Caenorhabditis genus. This review highlights recent comparative genomic studies of reproductive mode evolution addressing the hybrid origin of asexuality and the parallel gene loss following the emergence of hermaphroditism. It further summarizes ongoing efforts to characterize ancient linkage blocks called Nigon elements, which form central units of chromosome evolution. Fusions between Nigon elements have been demonstrated to impact recombination and speciation. Finally, multiple recent fusions between autosomal and the sex-linked Nigon element reveal insights into the dynamic evolution of sex chromosomes across various timescales.

Ancient diversity in host-parasite interaction genes in a model parasitic nematode

Host-parasite interactions exert strong selection pressures on the genomes of both host and parasite. These interactions can lead to negative frequency-dependent selection, a form of balancing selection that is hypothesised to explain the high levels of polymorphism seen in many host immune and parasite antigen loci. Here, we sequence the genomes of several individuals of Heligmosomoides bakeri, a model parasite of house mice, and Heligmosomoides polygyrus, a closely related parasite of wood mice. Although H. bakeri is commonly referred to as H. polygyrus in the literature, their genomes show levels of divergence that are consistent with at least a million years of independent evolution. The genomes of both species contain hyper-divergent haplotypes that are enriched for proteins that interact with the host immune response. Many of these haplotypes originated prior to the divergence between H. bakeri and H. polygyrus, suggesting that they have been maintained by long-term balancing selection. Together, our results suggest that the selection pressures exerted by the host immune response have played a key role in shaping patterns of genetic diversity in the genomes of parasitic nematodes.

Novel and improved Caenorhabditis briggsae gene models generated by community curation

BACKGROUND

The nematode Caenorhabditis briggsae has been used as a model in comparative genomics studies with Caenorhabditis elegans because of their striking morphological and behavioral similarities. However, the potential of C. briggsae for comparative studies is limited by the quality of its genome resources. The genome resources for the C. briggsae laboratory strain AF16 have not been developed to the same extent as C. elegans. The recent publication of a new chromosome-level reference genome for QX1410, a C. briggsae wild strain closely related to AF16, has provided the first step to bridge the gap between C. elegans and C. briggsae genome resources. Currently, the QX1410 gene models consist of software-derived gene predictions that contain numerous errors in their structure and coding sequences. In this study, a team of researchers manually inspected over 21,000 gene models and underlying transcriptomic data to repair software-derived errors.

RESULTS

We designed a detailed workflow to train a team of nine students to manually curate gene models using RNA read alignments. We manually inspected the gene models, proposed corrections to the coding sequences of over 8,000 genes, and modeled thousands of putative isoforms and untranslated regions. We exploited the conservation of protein sequence length between C. briggsae and C. elegans to quantify the improvement in protein-coding gene model quality and showed that manual curation led to substantial improvements in the protein sequence length accuracy of QX1410 genes. Additionally, collinear alignment analysis between the QX1410 and AF16 genomes revealed over 1,800 genes affected by spurious duplications and inversions in the AF16 genome that are now resolved in the QX1410 genome.

CONCLUSIONS

Community-based, manual curation using transcriptome data is an effective approach to improve the quality of software-derived protein-coding genes. The detailed protocols provided in this work can be useful for future large-scale manual curation projects in other species. Our manual curation efforts have brought the QX1410 gene models to a comparable level of quality as the extensively curated AF16 gene models. The improved genome resources for C. briggsae provide reliable tools for the study of Caenorhabditis biology and other related nematodes.

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.

The CERV protein of Cer1, a C. elegans LTR retrotransposon, is required for nuclear export of viral genomic RNA and can form giant nuclear rods

Retroviruses and closely related LTR retrotransposons export full-length, unspliced genomic RNA (gRNA) for packaging into virions and to serve as the mRNA encoding GAG and POL polyproteins. Because gRNA often includes splice acceptor and donor sequences used to splice viral mRNAs, retroelements must overcome host mechanisms that retain intron-containing RNAs in the nucleus. Here we examine gRNA expression in Cer1, an LTR retrotransposon in C. elegans which somehow avoids silencing and is highly expressed in germ cells. Newly exported Cer1 gRNA associates rapidly with the Cer1 GAG protein, which has structural similarity with retroviral GAG proteins. gRNA export requires CERV (C. elegans regulator of viral expression), a novel protein encoded by a spliced Cer1 mRNA. CERV phosphorylation at S214 is essential for gRNA export, and phosphorylated CERV colocalizes with nuclear gRNA at presumptive sites of transcription. By electron microscopy, tagged CERV proteins surround clusters of distinct, linear fibrils that likely represent gRNA molecules. Single fibrils, or groups of aligned fibrils, also localize near nuclear pores. During the C. elegans self-fertile period, when hermaphrodites fertilize oocytes with their own sperm, CERV concentrates in two nuclear foci that are coincident with gRNA. However, as hermaphrodites cease self-fertilization, and can only produce cross-progeny, CERV undergoes a remarkable transition to form giant nuclear rods or cylinders that can be up to 5 microns in length. We propose a novel mechanism of rod formation, in which stage-specific changes in the nucleolus induce CERV to localize to the nucleolar periphery in flattened streaks of protein and gRNA; these streaks then roll up into cylinders. The rods are a widespread feature of Cer1 in wild strains of C. elegans, but their function is not known and might be limited to cross-progeny. We speculate that the adaptive strategy Cer1 uses for the identical self-progeny of a host hermaphrodite might differ for heterozygous cross-progeny sired by males. For example, mating introduces male chromosomes which can have different, or no, Cer1 elements.

Novel and improved Caenorhabditis briggsae gene models generated by community curation

Background

The nematode Caenorhabditis briggsae has been used as a model for genomics studies compared to Caenorhabditis elegans because of its striking morphological and behavioral similarities. These studies yielded numerous findings that have expanded our understanding of nematode development and evolution. However, the potential of C. briggsae to study nematode biology is limited by the quality of its genome resources. The reference genome and gene models for the C. briggsae laboratory strain AF16 have not been developed to the same extent as C. elegans . The recent publication of a new chromosome-level reference genome for QX1410, a C. briggsae wild strain closely related to AF16, has provided the first step to bridge the gap between C. elegans and C. briggsae genome resources. Currently, the QX1410 gene models consist of protein-coding gene predictions generated from short- and long-read transcriptomic data. Because of the limitations of gene prediction software, the existing gene models for QX1410 contain numerous errors in their structure and coding sequences. In this study, a team of researchers manually inspected over 21,000 software-derived gene models and underlying transcriptomic data to improve the protein-coding gene models of the C. briggsae QX1410 genome.

Results

We designed a detailed workflow to train a team of nine students to manually curate genes using RNA read alignments and predicted gene models. We manually inspected the gene models using the genome annotation editor, Apollo, and proposed corrections to the coding sequences of over 8,000 genes. Additionally, we modeled thousands of putative isoforms and untranslated regions. We exploited the conservation of protein sequence length between C. briggsae and C. elegans to quantify the improvement in protein-coding gene model quality before and after curation. Manual curation led to a substantial improvement in the protein sequence length accuracy of QX1410 genes. We also compared the curated QX1410 gene models against the existing AF16 gene models. The manual curation efforts yielded QX1410 gene models that are similar in quality to the extensively curated AF16 gene models in terms of protein-length accuracy and biological completeness scores. Collinear alignment analysis between the QX1410 and AF16 genomes revealed over 1,800 genes affected by spurious duplications and inversions in the AF16 genome that are now resolved in the QX1410 genome.

Conclusions

Community-based, manual curation using transcriptome data is an effective approach to improve the quality of software-derived protein-coding genes. Comparative genomic analysis using a related species with high-quality reference genome(s) and gene models can be used to quantify improvements in gene model quality in a newly sequenced genome. The detailed protocols provided in this work can be useful for future large-scale manual curation projects in other species. The chromosome-level reference genome for the C. briggsae strain QX1410 far surpasses the quality of the genome of the laboratory strain AF16, and our manual curation efforts have brought the QX1410 gene models to a comparable level of quality to the previous reference, AF16. The improved genome resources for C. briggsae provide reliable tools for the study of Caenorhabditis biology and other related nematodes.

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.

The compact genome of Caenorhabditis niphades n. sp., isolated from a wood-boring weevil, Niphades variegatus

BACKGROUND

The first metazoan genome sequenced, that of Caenorhabditis elegans, has motivated animal genome evolution studies. To date > 50 species from the genus Caenorhabditis have been sequenced, allowing research on genome variation.

RESULTS

In the present study, we describe a new gonochoristic species, Caenorhabditis niphades n. sp., previously referred as C. sp. 36, isolated from adult weevils (Niphades variegatus), with whom they appear to be tightly associated during its life cycle. Along with a species description, we sequenced the genome of C. niphades n. sp. and produced a chromosome-level assembly. A genome comparison highlighted that C. niphades n. sp. has the smallest genome (59 Mbp) so far sequenced in the Elegans supergroup, despite being closely related to a species with an exceptionally large genome, C. japonica.

CONCLUSIONS

The compact genome of C. niphades n. sp. can serve as a key resource for comparative evolutionary studies of genome and gene number expansions in Caenorhabditis species.

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.

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.

Slow Recovery from Inbreeding Depression Generated by the Complex Genetic Architecture of Segregating Deleterious Mutations

The deleterious effects of inbreeding have been of extreme importance to evolutionary biology, but it has been difficult to characterize the complex interactions between genetic constraints and selection that lead to fitness loss and recovery after inbreeding. Haploid organisms and selfing organisms like the nematode Caenorhabditis elegans are capable of rapid recovery from the fixation of novel deleterious mutation; however, the potential for recovery and genomic consequences of inbreeding in diploid, outcrossing organisms are not well understood. We sought to answer two questions: 1) Can a diploid, outcrossing population recover from inbreeding via standing genetic variation and new mutation? and 2) How does allelic diversity change during recovery? We inbred C. remanei, an outcrossing relative of C. elegans, through brother-sister mating for 30 generations followed by recovery at large population size. Inbreeding reduced fitness but, surprisingly, recovery from inbreeding at large populations sizes generated only very moderate fitness recovery after 300 generations. We found that 65% of ancestral single nucleotide polymorphisms (SNPs) were fixed in the inbred population, far fewer than the theoretical expectation of ∼99%. Under recovery, 36 SNPs across 30 genes involved in alimentary, muscular, nervous, and reproductive systems changed reproducibly across replicates, indicating that strong selection for fitness recovery does exist. Our results indicate that recovery from inbreeding depression via standing genetic variation and mutation is likely to be constrained by the large number of segregating deleterious variants present in natural populations, limiting the capacity for recovery of small populations.

Natural genetic variation as a tool for discovery in Caenorhabditis nematodes

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

Optimizing experimental design for genome sequencing and assembly with Oxford Nanopore Technologies

High quality reference genome sequences are the core of modern genomics. Oxford Nanopore Technologies (ONT) produces inexpensive DNA sequences, but has high error rates, which make sequence assembly and analysis difficult as genome size and complexity increases. Robust experimental design is necessary for ONT genome sequencing and assembly, but few studies have addressed eukaryotic organisms. Here, we present novel results using simulated and empirical ONT and DNA libraries to identify best practices for sequencing and assembly for several model species. We find that the unique error structure of ONT libraries causes errors to accumulate and assembly statistics plateau as sequence depth increases. High-quality assembled eukaryotic sequences require high-molecular-weight DNA extractions that increase sequence read length, and computational protocols that reduce error through pre-assembly correction and read selection. Our quantitative results will be helpful for researchers seeking guidance for de novo assembly projects.

A telomere-to-telomere assembly of Oscheius tipulae and the evolution of rhabditid nematode chromosomes

Eukaryotic chromosomes have phylogenetic persistence. In many taxa, each chromosome has a single functional centromere with essential roles in spindle attachment and segregation. Fusion and fission can generate chromosomes with no or multiple centromeres, leading to genome instability. Groups with holocentric chromosomes (where centromeric function is distributed along each chromosome) might be expected to show karyotypic instability. This is generally not the case, and in Caenorhabditis elegans, it has been proposed that the role of maintenance of a stable karyotype has been transferred to the meiotic pairing centers, which are found at one end of each chromosome. Here, we explore the phylogenetic stability of nematode chromosomes using a new telomere-to-telomere assembly of the rhabditine nematode Oscheius tipulae generated from nanopore long reads. The 60-Mb O. tipulae genome is resolved into six chromosomal molecules. We find the evidence of specific chromatin diminution at all telomeres. Comparing this chromosomal O. tipulae assembly with chromosomal assemblies of diverse rhabditid nematodes, we identify seven ancestral chromosomal elements (Nigon elements) and present a model for the evolution of nematode chromosomes through rearrangement and fusion of these elements. We identify frequent fusion events involving NigonX, the element associated with the rhabditid X chromosome, and thus sex chromosome-associated gene sets differ markedly between species. Despite the karyotypic stability, gene order within chromosomes defined by Nigon elements is not conserved. Our model for nematode chromosome evolution provides a platform for investigation of the tensions between local genome rearrangement and karyotypic evolution in generating extant genome architectures.

Chromosome-Scale Genome Assemblies of Aphids Reveal Extensively Rearranged Autosomes and Long-Term Conservation of the X Chromosome

Chromosome rearrangements are arguably the most dramatic type of mutations, often leading to rapid evolution and speciation. However, chromosome dynamics have only been studied at the sequence level in a small number of model systems. In insects, Diptera and Lepidoptera have conserved genome structure at the scale of whole chromosomes or chromosome arms. Whether this reflects the diversity of insect genome evolution is questionable given that many species exhibit rapid karyotype evolution. Here, we investigate chromosome evolution in aphids-an important group of hemipteran plant pests-using newly generated chromosome-scale genome assemblies of the green peach aphid (Myzus persicae) and the pea aphid (Acyrthosiphon pisum), and a previously published assembly of the corn-leaf aphid (Rhopalosiphum maidis). We find that aphid autosomes have undergone dramatic reorganization over the last 30 My, to the extent that chromosome homology cannot be determined between aphids from the tribes Macrosiphini (Myzus persicae and Acyrthosiphon pisum) and Aphidini (Rhopalosiphum maidis). In contrast, gene content of the aphid sex (X) chromosome remained unchanged despite rapid sequence evolution, low gene expression, and high transposable element load. To test whether rapid evolution of genome structure is a hallmark of Hemiptera, we compared our aphid assemblies with chromosome-scale assemblies of two blood-feeding Hemiptera (Rhodnius prolixus and Triatoma rubrofasciata). Despite being more diverged, the blood-feeding hemipterans have conserved synteny. The exceptional rate of structural evolution of aphid autosomes renders them an important emerging model system for studying the role of large-scale genome rearrangements in evolution.

Degradation of the Repetitive Genomic Landscape in a Close Relative of Caenorhabditis elegans

The abundance, diversity, and genomic distribution of repetitive elements is highly variable among species. These patterns are thought to be driven in part by reproductive mode and the interaction of selection and recombination, and recombination rates typically vary by chromosomal position. In the nematode Caenorhabditis elegans, repetitive elements are enriched at chromosome arms and depleted on centers, and this mirrors the chromosomal distributions of other genomic features such as recombination rate. How conserved is this genomic landscape of repeats, and what evolutionary forces maintain it? To address this, we compared the genomic organization of repetitive elements across five Caenorhabditis species with chromosome-level assemblies. As previously reported, repeat content is enriched on chromosome arms in most Caenorhabditis species, and no obvious patterns of repeat content associated with reproductive mode were observed. However, the fig-associated C. inopinata has experienced repetitive element expansion and reveals no association of global repeat density with chromosome position. Patterns of repeat superfamily specific distributions reveal this global pattern is driven largely by a few repeat superfamilies that in C. inopinata have expanded in number and have weak associations with chromosome position. Additionally, 15% of predicted protein-coding genes in C. inopinata align to transposon-related proteins. When these are excluded, C. inopinata has no enrichment of genes in chromosome centers, in contrast to its close relatives who all have such clusters. Forward evolutionary simulations reveal that chromosomal heterogeneity in recombination rate alone can generate structured repetitive genomic landscapes when insertions are weakly deleterious, whereas chromosomal heterogeneity in the fitness effects of transposon insertion can promote such landscapes across a variety of evolutionary scenarios. Thus, patterns of gene density along chromosomes likely contribute to global repetitive landscapes in this group, although other historical or genomic factors are needed to explain the idiosyncrasy of genomic organization of various transposable element taxa within C. inopinata. Taken together, these results highlight the power of comparative genomics and evolutionary simulations in testing hypotheses regarding the causes of genome organization.