Sex chromosome evolution in parasitic nematodes of humans

Hybridization and introgression of the mitochondrial genome between the two species Anisakis pegreffii and A. simplex (s.s.) using a wide genotyping approach: evolutionary and ecological implications

Anisakis pegreffii and A. simplex (s.s.) are the two zoonotic anisakids infecting cetaceans as well as pelagic/demersal fish and squids. In European waters, A. pegreffii prevails in the Mediterranean Sea, while A. simplex (s.s.) in the NE Atlantic Ocean. Abiotic conditions likely play a significant role in shaping their geographical distribution. The Iberian Atlantic and Alboran Sea waters are sympatric areas of the two species. A total of 429 adults and L3 stage from both sympatric and allopatric areas were studied by a wide nuclear genotyping approach (including newly and previously found diagnostic single nucleotide polymorphisms (SNPs) at nuclear DNA (nDNA) and microsatellite DNA loci) and sequenced at mitochondrial DNA (mtDNA) cox2. Admixture between the two species was detected in the sympatric areas studied by STRUCTURE Bayesian analysis; NEWHYBRIDS revealed different categories of hybridization between the two species, representing approximately 5%. A tendency for F1 female hybrids to interbreed with the parental species at the geographical distribution limits of both species was observed. This finding suggests that hybridization occurs when the two parental species significantly differ in abundance. Mitochondrial introgression of A. simplex (s.s.) in A. pegreffii from Mediterranean waters was also detected, likely as a result of past and/or paleo-introgression events. The high level of genetic differentiation between the two species and their backcrosses indicates that, despite current hybridization, reproductive isolation which maintains evolutionary boundaries between the two species, exists. Possible causes of hybridization phenomena are attempted, as well as their evolutionary and ecological implications, also considering a sea warming scenario in European waters.

Rapid chromosome evolution and acquisition of thermosensitive stochastic sex determination in nematode androdioecious hermaphrodites

The factors contributing to evolution of androdioecy, the coexistence of hermaphrodites and males such as in Caenorhabditis elegans, remains poorly known. However, nematodes exhibit androdioecy in at last 13 genera with the predatory genus Pristionchus having seven independent transitions towards androdioecy. Nonetheless, associated genomic architecture and sex determination mechanisms are largely known from Caenorhabditis. Here, studying 47 Pristionchus species, we observed repeated chromosome evolution which abolished the ancestral XX/XO sex chromosome system. Two phylogenetically unrelated androdioecious Pristionchus species have no genomic differences between sexes and mating hermaphrodites with males resulted in hermaphroditic offspring only. We demonstrate that stochastic sex determination is influenced by temperature in P. mayeri and P. entomophagus, and CRISPR engineering indicated a conserved role of the transcription factor TRA-1 in P. mayeri. Chromosome-level genome assemblies and subsequent genomic analysis of related Pristionchus species revealed stochastic sex determination to be derived from XY sex chromosome systems through sex chromosome-autosome fusions. Thus, rapid karyotype evolution, sex chromosome evolution and evolvable sex determination mechanisms are general features of this genus, and represent a dynamic background against which androdioecy has evolved recurrently. Future studies might indicate that stochastic sex determination is more common than currently appreciated.

Sex and neo-sex chromosome evolution in beetles

Beetles are the most species-rich group of animals and harbor diverse karyotypes. Most species have XY sex chromosomes, but X0 sex determination mechanisms are also common in some groups. We generated a whole-chromosome assembly of Tribolium confusum, which has a neo-sex chromosome, and utilize eleven additional beetle genomes to reconstruct karyotype evolution across Coleoptera. We identify ancestral linkage groups, termed Stevens elements, that share a conserved set of genes across beetles. While the ancestral X chromosome is maintained across beetles, we find independent additions of autosomes to the ancestral sex chromosomes. These neo-sex chromosomes evolve the stereotypical properties of sex chromosomes, including the evolution of dosage compensation and a non-random distribution of genes with sex-biased expression. Beetles thus provide a novel model to gain a better understanding of the diverse forces driving sex chromosome evolution.

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.

Chromosome fusion and programmed DNA elimination shape karyotypes of nematodes

An increasing number of metazoans undergo programmed DNA elimination (PDE), where a significant amount of DNA is selectively lost from the somatic genome during development. In some nematodes, PDE leads to the removal and remodeling of the ends of all germline chromosomes. In several species, PDE also generates internal breaks that lead to sequence loss and increased numbers of somatic chromosomes. The biological significance of these karyotype changes associated with PDE and the origin and evolution of nematode PDE remain largely unknown. Here, we assembled the single germline chromosome of the nematode Parascaris univalens and compared the karyotypes, chromosomal gene organization, and PDE features among other nematodes. We show that PDE in Parascaris converts an XX/XY sex-determination system in the germline into an XX/XO system in the somatic cells. Comparisons of Ascaris, Parascaris, and Baylisascaris ascarid chromosomes suggest that PDE existed in the ancestor of these nematodes, and their current distinct germline karyotypes were derived from fusion events of smaller ancestral chromosomes. The DNA breaks involved in PDE resolve these fused germline chromosomes into their pre-fusion karyotypes. These karyotype changes may lead to alterations in genome architecture and gene expression in the somatic cells. Cytological and genomic analyses further suggest that satellite DNA and the heterochromatic chromosome arms are dynamic and may play a role during meiosis. Overall, our results show that chromosome fusion and PDE have been harnessed in these ascarids to sculpt their karyotypes, altering the genome organization and serving specific functions in the germline and somatic cells.

Dual RNA-seq in filarial nematodes and Wolbachia endosymbionts using RNase H based ribosomal RNA depletion

Lymphatic filariasis is caused by parasitic nematodes and is a leading cause of disability worldwide. Many filarial worms contain the bacterium Wolbachia as an obligate endosymbiont. RNA sequencing is a common technique used to study their molecular relationships and to identify potential drug targets against the nematode and bacteria. Ribosomal RNA (rRNA) is the most abundant RNA species, accounting for 80-90% of the RNA in a sample. To reduce sequencing costs, it is necessary to remove ribosomal reads through poly-A enrichment or ribosomal depletion. Bacterial RNA does not contain a poly-A tail, making it difficult to sequence both the nematode and Wolbachia from the same library preparation using standard poly-A selection. Ribosomal depletion can utilize species-specific oligonucleotide probes to remove rRNA through pull-down or degradation methods. While species-specific probes are commercially available for many commonly studied model organisms, there are currently limited depletion options for filarial parasites. Here, we performed total RNA sequencing from Brugia malayi containing the Wolbachia symbiont (wBm) and designed ssDNA depletion probes against their rRNA sequences. We compared the total RNA library to poly-A enriched, Terminator 5'-Phosphate-Dependent Exonuclease treated, NEBNext Human/Bacteria rRNA depleted and our custom nematode probe depleted libraries. The custom nematode depletion library had the lowest percentage of ribosomal reads across all methods, with a 300-fold decrease in rRNA when compared to the total RNA library. The nematode depletion libraries also contained the highest percentage of Wolbachia mRNA reads, resulting in a 16-1,000-fold increase in bacterial reads compared to the other enrichment and depletion methods. Finally, we found that the Brugia malayi depletion probes can remove rRNA from the filarial worm Dirofilaria immitis and the majority of rRNA from the more distantly related free living nematode Caenorhabditis elegans. These custom filarial probes will allow for future dual RNA-seq experiments between nematodes and their bacterial symbionts from a single sequencing library.

Genome structure and population genomics of the canine heartworm Dirofilaria immitis

The heartworm, Dirofilaria immitis, is a filarial parasitic nematode responsible for significant morbidity and mortality in wild and domesticated canids. Resistance to macrocyclic lactone drug prevention represents a significant threat to parasite control and has prompted investigations to understand the genetic determinants of resistance. This study aimed to improve the genomic resources of D. immitis to enable a more precise understanding of how genetic variation is distributed within and between parasite populations worldwide, which will inform the likelihood and rate by which parasites, and in turn, resistant alleles, might spread. We have guided the scaffolding of a recently published genome assembly for D. immitis (ICBAS_JMDir_1.0) using the chromosomal-scale reference genomes of Brugia malayi and Onchocerca volvulus, resulting in an 89.5 Mb assembly composed of four autosomal- and one sex-linked chromosomal-scale scaffolds representing 99.7% of the genome. Publicly available and new whole-genome sequencing data from 32 D. immitis samples from Australia, Italy and the USA were assessed using principal component analysis, nucleotide diversity (Pi) and absolute genetic divergence (Dxy) to characterise the global genetic structure and measure within- and between-population diversity. These population genetic analyses revealed broad-scale genetic structure among globally diverse samples and differences in genetic diversity between populations; however, fine-scale subpopulation analysis was limited and biased by differences between sample types. Finally, we mapped single nucleotide polymorphisms previously associated with macrocyclic lactone resistance in the new genome assembly, revealing the physical linkage of high-priority variants on chromosome 3, and determined their frequency in the studied populations. This new chromosomal assembly for D. immitis now allows for a more precise investigation of selection on genome-wide genetic variation and will enhance our understanding of parasite transmission and the spread of genetic variants responsible for resistance to treatment.

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'.

The genome of Litomosoides sigmodontis illuminates the origins of Y chromosomes in filarial nematodes

Heteromorphic sex chromosomes are usually thought to have originated from a pair of autosomes that acquired a sex-determining locus and subsequently stopped recombining, leading to degeneration of the sex-limited chromosome. The majority of nematode species lack heteromorphic sex chromosomes and determine sex using an X-chromosome counting mechanism, with males being hemizygous for one or more X chromosomes (XX/X0). Some filarial nematode species, including important parasites of humans, have heteromorphic XX/XY karyotypes. It has been assumed that sex is determined by a Y-linked locus in these species. However, karyotypic analyses suggested that filarial Y chromosomes are derived from the unfused homologue of an autosome involved in an X-autosome fusion event. Here, we generated a chromosome-level reference genome for Litomosoides sigmodontis, a filarial nematode with the ancestral filarial karyotype and sex determination mechanism (XX/X0). By mapping the assembled chromosomes to the rhabditid nematode ancestral linkage (or Nigon) elements, we infer that the ancestral filarial X chromosome was the product of a fusion between NigonX (the ancestrally X-linked element) and NigonD (ancestrally autosomal). In the two filarial lineages with XY systems, there have been two independent X-autosome chromosome fusion events involving different autosomal Nigon elements. In both lineages, the region shared by the neo-X and neo-Y chromosomes is within the ancestrally autosomal portion of the X, confirming that the filarial Y chromosomes are derived from the unfused homologue of the autosome. Sex determination in XY filarial nematodes therefore likely continues to operate via the ancestral X-chromosome counting mechanism, rather than via a Y-linked sex-determining locus.

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.

Chromosome fusion and programmed DNA elimination shape karyotypes of parasitic nematodes

A growing list of metazoans undergo programmed DNA elimination (PDE), where a significant amount of DNA is selectively lost from the somatic genome during development. In some nematodes, PDE leads to the removal and remodeling of the ends of all germline chromosomes. In several species, PDE also generates internal breaks that lead to sequence loss and an increased number of somatic chromosomes. The biological significance of these karyotype changes associated with PDE and the origin and evolution of nematode PDE remain largely unknown. Here, we assembled the single germline chromosome of the horse parasite Parascaris univalens and compared the karyotypes, chromosomal gene organization, and PDE features among ascarid nematodes. We show that PDE in Parascaris converts an XX/XY sex-determination system in the germline into an XX/XO system in the somatic cells. Comparisons of Ascaris, Parascaris, and Baylisascaris ascarid chromosomes suggest that PDE existed in the ancestor of these parasites, and their current distinct germline karyotypes were derived from fusion events of smaller ancestral chromosomes. The DNA breaks involved in PDE resolve these fused germline chromosomes into their pre-fusion karyotypes, leading to alterations in genome architecture and gene expression in the somatic cells. Cytological and genomic analyses further suggest that satellite DNA and the heterochromatic chromosome arms play a dynamic role in the Parascaris germline chromosome during meiosis. Overall, our results show that chromosome fusion and PDE have been harnessed in these ascarids to sculpt their karyotypes, altering the genome organization and serving specific functions in the germline and somatic cells.

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.

Thousands of Pristionchus pacificus orphan genes were integrated into developmental networks that respond to diverse environmental microbiota

Adaptation of organisms to environmental change may be facilitated by the creation of new genes. New genes without homologs in other lineages are known as taxonomically-restricted orphan genes and may result from divergence or de novo formation. Previously, we have extensively characterized the evolution and origin of such orphan genes in the nematode model organism Pristionchus pacificus. Here, we employ large-scale transcriptomics to establish potential functional associations and to measure the degree of transcriptional plasticity among orphan genes. Specifically, we analyzed 24 RNA-seq samples from adult P. pacificus worms raised on 24 different monoxenic bacterial cultures. Based on coexpression analysis, we identified 28 large modules that harbor 3,727 diplogastrid-specific orphan genes and that respond dynamically to different bacteria. These coexpression modules have distinct regulatory architecture and also exhibit differential expression patterns across development suggesting a link between bacterial response networks and development. Phylostratigraphy revealed a considerably high number of family- and even species-specific orphan genes in certain coexpression modules. This suggests that new genes are not attached randomly to existing cellular networks and that integration can happen very fast. Integrative analysis of protein domains, gene expression and ortholog data facilitated the assignments of biological labels for 22 coexpression modules with one of the largest, fast-evolving module being associated with spermatogenesis. In summary, this work presents the first functional annotation for thousands of P. pacificus orphan genes and reveals insights into their integration into environmentally responsive gene networks.

Genomes of the human filarial parasites Mansonella perstans and Mansonella ozzardi

The filarial parasites Mansonella ozzardi and Mansonella perstans, causative agents of mansonellosis, infect hundreds of millions of people worldwide, yet remain among the most understudied of the human filarial pathogens. M. ozzardi is highly prevalent in Latin American countries and Caribbean Islands, while M. perstans is predominantly found in sub-Saharan Africa as well as in a few areas in South America. In addition to the differences in their geographical distribution, the two parasites are transmitted by different insect vectors, as well as exhibit differences in their responses to commonly used anthelminthic drugs. The lack of genome information has hindered investigations into the biology and evolution of Mansonella parasites and understanding the molecular basis of the clinical differences between species. In the current study, high quality genomes of two independent clinical isolates of M. perstans from Cameroon and two M. ozzardi isolates one from Brazil and one from Venezuela are reported. The genomes are approximately 76 Mb in size, encode about 10,000 genes each, and are largely complete based on BUSCO scores of about 90%, similar to other completed filarial genomes. These sequences represent the first genomes from Mansonella parasites and enabled a comparative genomic analysis of the similarities and differences between Mansonella and other filarial parasites. Horizontal DNA transfers (HDT) from mitochondria (nuMTs) as well as transfers from genomes of endosymbiotic Wolbachia bacteria (nuWTs) to the host nuclear genome were identified and analyzed. Sequence comparisons and phylogenetic analysis of known targets of anti-filarial drugs diethylcarbamazine (DEC), ivermectin and mebendazole revealed that all known target genes were present in both species, except for the DEC target encoded by gon-2 gene, which is fragmented in genome assemblies from both M. ozzardi isolates. These new reference genome sequences will provide a valuable resource for further studies on biology, symbiosis, evolution and drug discovery.

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.

Evolution and Diversity of TGF-β Pathways are Linked with Novel Developmental and Behavioral Traits

Transforming growth factor-β (TGF-β) signaling is essential for numerous biologic functions. It is a highly conserved pathway found in all metazoans including the nematode Caenorhabditis elegans, which has also been pivotal in identifying many components. Utilizing a comparative evolutionary approach, we explored TGF-β signaling in nine nematode species and revealed striking variability in TGF-β gene frequency across the lineage. Of the species analyzed, gene duplications in the DAF-7 pathway appear common with the greatest disparity observed in Pristionchus pacificus. Specifically, multiple paralogues of daf-3, daf-4 and daf-7 were detected. To investigate this additional diversity, we induced mutations in 22 TGF-β components and generated corresponding double, triple, and quadruple mutants revealing both conservation and diversification in function. Although the DBL-1 pathway regulating body morphology appears highly conserved, the DAF-7 pathway exhibits functional divergence, notably in some aspects of dauer formation. Furthermore, the formation of the phenotypically plastic mouth in P. pacificus is partially influenced through TGF-β with the strongest effect in Ppa-tag-68. This appears important for numerous processes in P. pacificus but has no known function in C. elegans. Finally, we observe behavioral differences in TGF-β mutants including in chemosensation and the establishment of the P. pacificus kin-recognition signal. Thus, TGF-β signaling in nematodes represents a stochastic genetic network capable of generating novel functions through the duplication and deletion of associated genes.

Improving helminth genome resources in the post-genomic era

Rapid advancement in high-throughput sequencing and analytical approaches has seen a steady increase in the generation of genomic resources for helminth parasites. Now, helminth genomes and their annotations are a cornerstone of numerous efforts to compare genetic and transcriptomic variation, from single cells to populations of globally distributed parasites, to genome modifications to understand gene function. Our understanding of helminths is increasingly reliant on these genomic resources, which are primarily static once published and vary widely in quality and completeness between species. This article seeks to highlight the cause and effect of this variation and argues for the continued improvement of these genomic resources - even after their publication - which is necessary to provide a more accurate and complete understanding of the biology of these important pathogens.

Multiple Pristionchus pacificus genomes reveal distinct evolutionary dynamics between de novo candidates and duplicated genes

The birth of new genes is a major molecular innovation driving phenotypic diversity across all domains of life. Although repurposing of existing protein-coding material by duplication is considered the main process of new gene formation, recent studies have discovered thousands of transcriptionally active sequences as a rich source of new genes. However, differential loss rates have to be assumed to reconcile the high birth rates of these incipient de novo genes with the dominance of ancient gene families in individual genomes. Here, we test this rapid turnover hypothesis in the context of the nematode model organism Pristionchus pacificus We extended the existing species-level phylogenomic framework by sequencing the genomes of six divergent P. pacificus strains. We used these data to study the evolutionary dynamics of different age classes and categories of origin at a population level. Contrasting de novo candidates with new families that arose by duplication and divergence from known genes, we find that de novo candidates are typically shorter, show less expression, and are overrepresented on the sex chromosome. Although the contribution of de novo candidates increases toward young age classes, multiple comparisons within the same age class showed significantly higher attrition in de novo candidates than in known genes. Similarly, young genes remain under weak evolutionary constraints with de novo candidates representing the fastest evolving subcategory. Altogether, this study provides empirical evidence for the rapid turnover hypothesis and highlights the importance of the evolutionary timescale when quantifying the contribution of different mechanisms toward new gene formation.

Evolution of sexual systems, sex chromosomes and sex-linked gene transcription in flatworms and roundworms

Many species with separate male and female individuals (termed ‘gonochorism’ in animals) have sex-linked genome regions. Here, we investigate evolutionary changes when genome regions become completely sex-linked, by analyses of multiple species of flatworms (Platyhelminthes; among which schistosomes recently evolved gonochorism from ancestral hermaphroditism), and roundworms (Nematoda) which have undergone independent translocations of different autosomes. Although neither the evolution of gonochorism nor translocations fusing ancestrally autosomal regions to sex chromosomes causes inevitable loss of recombination, we document that formerly recombining regions show genomic signatures of recombination suppression in both taxa, and become strongly genetically degenerated, with a loss of most genes. Comparisons with hermaphroditic flatworm transcriptomes show masculinisation and some defeminisation in schistosome gonad gene expression. We also find evidence that evolution of sex-linkage in nematodes is accompanied by transcriptional changes and dosage compensation. Our analyses also identify sex-linked genes that could assist future research aimed at controlling some of these important parasites.

Transitions between hermaphroditic and separate sexes are relatively understudied in animals compared to pants. Here, Wang et al. reconstruct the evolution of separate sexes in the flatworms and complex changes of sex chromosomes in the roundworms.

Recurrent chromosome reshuffling and the evolution of neo-sex chromosomes in parrots

The karyotype of most birds has remained considerably stable during more than 100 million years' evolution, except for some groups, such as parrots. The evolutionary processes and underlying genetic mechanism of chromosomal rearrangements in parrots, however, are poorly understood. Here, using chromosome-level assemblies of four parrot genomes, we uncover frequent chromosome fusions and fissions, with most of them occurring independently among lineages. The increased activities of chromosomal rearrangements in parrots are likely associated with parrot-specific loss of two genes, ALC1 and PARP3, that have known functions in the repair of double-strand breaks and maintenance of genome stability. We further find that the fusion of the ZW sex chromosomes and chromosome 11 has created a pair of neo-sex chromosomes in the ancestor of parrots, and the chromosome 25 has been further added to the sex chromosomes in monk parakeet. Together, the combination of our genomic and cytogenetic analyses characterizes the complex evolutionary history of chromosomal rearrangements and sex chromosomes in parrots.

The genome sequence of the Australian filarial nematode, Cercopithifilaria johnstoni

We present a genome assembly and annotation of an individual female Cercopithifilaria johnstoni, a parasitic filarial nematode that is transmitted by hard ticks (Ixodidae) to infect a broad range of native Australian murid and marsupial hosts. The genome sequence is 76.9 Mbp in length, and although in draft form (N50 = 99 kbp, N50[n] = 232), is largely complete based on universally conserved orthologs (BUSCOs; genome = 94.9%, protein = 96.5%) and relative to other related filarial species. These data represent the first genomic resources for the genus Cercopithifilaria, a group of parasites with a broad host range, and form the basis for comparative analysis with the human-infective parasite, Onchocerca volvulus, both of which are responsible for similar eye and skin pathologies in their respective hosts.

The genome sequence of the Australian filarial nematode, Cercopithifilaria johnstoni

We present a genome assembly and annotation of an individual female Cercopithifilaria johnstoni, a parasitic filarial nematode that is transmitted by hard ticks (Ixodidae) to infect a broad range of native Australian murid and marsupial hosts. The genome sequence is 76.9 Mbp in length, and although in draft form (N50 = 99 kbp, N50[n] = 232), is largely complete based on universally conserved orthologs (BUSCOs; genome = 94.9%, protein = 96.5%) and relative to other related filarial species. These data represent the first genomic resources for the genus Cercopithifilaria, a group of parasites with a broad host range, and form the basis for comparative analysis with the human-infective parasite, Onchocerca volvulus, both of which are responsible for similar eye and skin pathologies in their respective hosts.

X-treme loss of sequence diversity linked to neo-X chromosomes in filarial nematodes

The sequence diversity of natural and laboratory populations of Brugia pahangi and Brugia malayi was assessed with Illumina resequencing followed by mapping in order to identify single nucleotide variants and insertions/deletions. In natural and laboratory Brugia populations, there is a lack of sequence diversity on chromosome X relative to the autosomes (πX/πA = 0.2), which is lower than the expected (πX/πA = 0.75). A reduction in diversity is also observed in other filarial nematodes with neo-X chromosome fusions in the genera Onchocerca and Wuchereria, but not those without neo-X chromosome fusions in the genera Loa and Dirofilaria. In the species with neo-X chromosome fusions, chromosome X is abnormally large, containing a third of the genetic material such that a sizable portion of the genome is lacking sequence diversity. Such profound differences in genetic diversity can be consequential, having been associated with drug resistance and adaptability, with the potential to affect filarial eradication.

Analysis of meiosis in Pristionchus pacificus reveals plasticity in homolog pairing and synapsis in the nematode lineage

Meiosis is conserved across eukaryotes yet varies in the details of its execution. Here we describe a new comparative model system for molecular analysis of meiosis, the nematode Pristionchus pacificus, a distant relative of the widely studied model organism Caenorhabditis elegans. P. pacificus shares many anatomical and other features that facilitate analysis of meiosis in C. elegans. However, while C. elegans has lost the meiosis-specific recombinase Dmc1 and evolved a recombination-independent mechanism to synapse its chromosomes, P. pacificus expresses both DMC-1 and RAD-51. We find that SPO-11 and DMC-1 are required for stable homolog pairing, synapsis, and crossover formation, while RAD-51 is dispensable for these key meiotic processes. RAD-51 and DMC-1 localize sequentially to chromosomes during meiotic prophase and show nonoverlapping functions. We also present a new genetic map for P. pacificus that reveals a crossover landscape very similar to that of C. elegans, despite marked divergence in the regulation of synapsis and crossing-over between these lineages.

Fisher vs. the Worms: Extraordinary Sex Ratios in Nematodes and the Mechanisms that Produce Them

Parker, Baker, and Smith provided the first robust theory explaining why anisogamy evolves in parallel in multicellular organisms. Anisogamy sets the stage for the emergence of separate sexes, and for another phenomenon with which Parker is associated: sperm competition. In outcrossing taxa with separate sexes, Fisher proposed that the sex ratio will tend towards unity in large, randomly mating populations due to a fitness advantage that accrues in individuals of the rarer sex. This creates a vast excess of sperm over that required to fertilize all available eggs, and intense competition as a result. However, small, inbred populations can experience selection for skewed sex ratios. This is widely appreciated in haplodiploid organisms, in which females can control the sex ratio behaviorally. In this review, we discuss recent research in nematodes that has characterized the mechanisms underlying highly skewed sex ratios in fully diploid systems. These include self-fertile hermaphroditism and the adaptive elimination of sperm competition factors, facultative parthenogenesis, non-Mendelian meiotic oddities involving the sex chromosomes, and environmental sex determination. By connecting sex ratio evolution and sperm biology in surprising ways, these phenomena link two "seminal" contributions of G. A. Parker.

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.

Investigating the genetic structure of the parasites Anisakis pegreffii and A. berlandi (Nematoda: Anisakidae) in a sympatric area of the southern Pacific Ocean waters using a multilocus genotyping approach: first evidence of their interspecific hybridization

The southern Pacific Ocean, off the New Zealand coast, has been reported as one sympatric area of the two parasite species Anisakis pegreffii and A. berlandi. Here, a multilocus genotyping approach, based on a panel of eleven DNA microsatellite (SSR) loci plus the sequences analysis of the nuclear nas10 nDNA and the mitochondrial mtDNA cox2 gene loci, was applied to a total of N = 344 adults and larvae of Anisakis spp. from cetacean and fish species, respectively. Out of the newly scored SSR loci, Anisl 15 and Anisl 2 showed fixed alternative alleles between A. pegreffii and A. berlandi resulting as 100% diagnostic loci. Out of SSRs Anisl 00314 and Anisl 7 previously disclosed, two additional loci, i.e., Anisl 4 and Anisl 22, were found to be sex-linked. The Bayesian genotypes clustering approach (STRUCTURE) allowed identification, with a 100% of probability value, N = 208 specimens to the "pure parental" A. pegreffii, N = 133 to the "pure parental" A. berlandi, while one adult and two larval stages showed mixed ancestry between the two groups having, in all cases, a Q-value = 0.50. NEWHYBRIDS analysis assigned (100% of probability) those specimens to their F1 hybrid category. This represents the first evidence of contemporary hybridization between the two parasite species in a sympatric area. The pairwise F_ST values estimated at intraspecific and interspecific level, inferred from both SSR loci and mitochondrial mtDNA cox2 sequences, have also demonstrated the existence of two distinct panmictic units in this study area, corresponding respectively to A. pegreffii and A. berlandi. The results obtained support the useful application of a multilocus approach in the identification of sibling species and their hybrid categories in sympatric areas. The possible use of sex-linked SSR loci of the two species of the A. simplex (s. l.), for sex determination of their larval stages, is also suggested.

Transcriptional profiles in Strongyloides stercoralis males reveal deviations from the Caenorhabditis sex determination model

The human and canine parasitic nematode Strongyloides stercoralis utilizes an XX/XO sex determination system, with parasitic females reproducing by mitotic parthenogenesis and free-living males and females reproducing sexually. However, the genes controlling S. stercoralis sex determination and male development are unknown. We observed precocious development of rhabditiform males in permissive hosts treated with corticosteroids, suggesting that steroid hormones can regulate male development. To examine differences in transcript abundance between free-living adult males and other developmental stages, we utilized RNA-Seq. We found two clusters of S. stercoralis-specific genes encoding predicted transmembrane proteins that are only expressed in free-living males. We additionally identified homologs of several genes important for sex determination in Caenorhabditis species, including mab-3, tra-1, fem-2, and sex-1, which may have similar functions. However, we identified three paralogs of gld-1; Ss-qki-1 transcripts were highly abundant in adult males, while Ss-qki-2 and Ss-qki-3 transcripts were highly abundant in adult females. We also identified paralogs of pumilio domain-containing proteins with sex-specific transcripts. Intriguingly, her-1 appears to have been lost in several parasite lineages, and we were unable to identify homologs of tra-2 outside of Caenorhabditis species. Together, our data suggest that different mechanisms control male development in S. stercoralis and Caenorhabditis species.

The community-curated Pristionchus pacificus genome facilitates automated gene annotation improvement in related nematodes

Background

The nematode Pristionchus pacificus is an established model organism for comparative studies with Caenorhabditis elegans. Over the past years, it developed into an independent animal model organism for elucidating the genetic basis of phenotypic plasticity. Community-based curations were employed recently to improve the quality of gene annotations of P. pacificus and to more easily facilitate reverse genetic studies using candidate genes from C. elegans.

Results

Here, I demonstrate that the reannotation of phylogenomic data from nine related nematode species using the community-curated P. pacificus gene set as homology data substantially improves the quality of gene annotations. Benchmarking of universal single copy orthologs (BUSCO) estimates a median completeness of 84% which corresponds to a 9% increase over previous annotations. Nevertheless, the ability to infer gene models based on homology already drops beyond the genus level reflecting the rapid evolution of nematode lineages. This also indicates that the highly curated C. elegans genome is not optimally suited for annotating non-Caenorhabditis genomes based on homology. Furthermore, comparative genomic analysis of apparently missing BUSCO genes indicates a failure of ortholog detection by the BUSCO pipeline due to the insufficient sample size and phylogenetic breadth of the underlying OrthoDB data set. As a consequence, the quality of multiple divergent nematode genomes might be underestimated.

Conclusions

This study highlights the need for optimizing gene annotation protocols and it demonstrates the benefit of a high quality genome for phylogenomic data of related species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07529-x.

FADU: a Quantification Tool for Prokaryotic Transcriptomic Analyses

Most currently available quantification tools for transcriptomics analyses have been designed for human data sets, in which full-length transcript sequences, including the untranslated regions, are well annotated. In most prokaryotic systems, full-length transcript sequences have yet to be characterized, leading to prokaryotic transcriptomics analyses being performed based on only the coding sequences.

Quantification tools for RNA sequencing (RNA-Seq) analyses are often designed and tested using human transcriptomics data sets, in which full-length transcript sequences are well annotated. For prokaryotic transcriptomics experiments, full-length transcript sequences are seldom known, and coding sequences must instead be used for quantification steps in RNA-Seq analyses. However, operons confound accurate quantification of coding sequences since a single transcript does not necessarily equate to a single gene. Here, we introduce FADU (Feature Aggregate Depth Utility), a quantification tool designed specifically for prokaryotic RNA-Seq analyses. FADU assigns partial count values proportional to the length of the fragment overlapping the target feature. To assess the ability of FADU to quantify genes in prokaryotic transcriptomics analyses, we compared its performance to those of eXpress, featureCounts, HTSeq, kallisto, and Salmon across three paired-end read data sets of (i) Ehrlichia chaffeensis, (ii) Escherichia coli, and (iii) the Wolbachia endosymbiont wBm. Across each of the three data sets, we find that FADU can more accurately quantify operonic genes by deriving proportional counts for multigene fragments within operons. FADU is available at https://github.com/IGS/FADU.

IMPORTANCE Most currently available quantification tools for transcriptomics analyses have been designed for human data sets, in which full-length transcript sequences, including the untranslated regions, are well annotated. In most prokaryotic systems, full-length transcript sequences have yet to be characterized, leading to prokaryotic transcriptomics analyses being performed based on only the coding sequences. In contrast to eukaryotes, prokaryotes contain polycistronic transcripts, and when genes are quantified based on coding sequences instead of transcript sequences, this leads to an increased abundance of improperly assigned ambiguous multigene fragments, specifically those mapping to multiple genes in operons. Here, we describe FADU, a quantification tool for prokaryotic RNA-Seq analyses designed to assign proportional counts with the purpose of better quantifying operonic genes while minimizing the pitfalls associated with improperly assigning fragment counts from ambiguous transcripts.

Improving Transgenesis Efficiency and CRISPR-Associated Tools Through Codon Optimization and Native Intron Addition in Pristionchus Nematodes

A lack of appropriate molecular tools is one obstacle that prevents in-depth mechanistic studies in many organisms. Transgenesis, clustered regularly interspaced short palindromic repeats (CRISPR)-associated engineering, and related tools are fundamental in the modern life sciences, but their applications are still limited to a few model organisms. In the phylum Nematoda, transgenesis can only be performed in a handful of species other than Caenorhabditis elegans, and additionally, other species suffer from significantly lower transgenesis efficiencies. We hypothesized that this may in part be due to incompatibilities of transgenes in the recipient organisms. Therefore, we investigated the genomic features of 10 nematode species from three of the major clades representing all different lifestyles. We found that these species show drastically different codon usage bias and intron composition. With these findings, we used the species Pristionchus pacificus as a proof of concept for codon optimization and native intron addition. Indeed, we were able to significantly improve transgenesis efficiency, a principle that may be usable in other nematode species. In addition, with the improved transgenes, we developed a fluorescent co-injection marker in P. pacificus for the detection of CRISPR-edited individuals, which helps considerably to reduce associated time and costs.

Long-read RNA sequencing of human and animal filarial parasites improves gene models and discovers operons

Filarial parasitic nematodes (Filarioidea) cause substantial disease burden to humans and animals around the world. Recently there has been a coordinated global effort to generate, annotate, and curate genomic data from nematode species of medical and veterinary importance. This has resulted in two chromosome-level assemblies (Brugia malayi and Onchocerca volvulus) and 11 additional draft genomes from Filarioidea. These reference assemblies facilitate comparative genomics to explore basic helminth biology and prioritize new drug and vaccine targets. While the continual improvement of genome contiguity and completeness advances these goals, experimental functional annotation of genes is often hindered by poor gene models. Short-read RNA sequencing data and expressed sequence tags, in cooperation with ab initio prediction algorithms, are employed for gene prediction, but these can result in missing clade-specific genes, fragmented models, imperfect mapping of gene ends, and lack of isoform resolution. Long-read RNA sequencing can overcome these drawbacks and greatly improve gene model quality. Here, we present Iso-Seq data for B. malayi and Dirofilaria immitis, etiological agents of lymphatic filariasis and canine heartworm disease, respectively. These data cover approximately half of the known coding genomes and substantially improve gene models by extending untranslated regions, cataloging novel splice junctions from novel isoforms, and correcting mispredicted junctions. Furthermore, we validated computationally predicted operons, manually curated new operons, and merged fragmented gene models. We carried out analyses of poly(A) tails in both species, leading to the identification of non-canonical poly(A) signals. Finally, we prioritized and assessed known and putative anthelmintic targets, correcting or validating gene models for molecular cloning and target-based anthelmintic screening efforts. Overall, these data significantly improve the catalog of gene models for two important parasites, and they demonstrate how long-read RNA sequencing should be prioritized for ongoing improvement of parasitic nematode genome assemblies.

Filarial parasitic nematodes are vector-borne parasites that infect humans and animals. Brugia malayi and Dirofilaria immitis are transmitted by mosquitoes and cause human lymphatic filariasis and canine heartworm disease, respectively. Recent years have seen a dramatic increase in genomic and transcriptomic data sets and the concomitant increase in innovative strategies for drug target identification, validation, and screening. However, while the completeness of genome assemblies of filarial parasitic nematodes has seen steady improvements, the reliability of gene models has not kept pace, hindering cloning efforts. Long-read RNA sequencing technologies are uniquely able to improve gene models, but have not been widely used for the causative agents of neglected tropical diseases. Here, we report the improvement of gene models in both B. malayi and D. immitis by long-read RNA sequencing. We identified novel operons, deprecated false positive operons, identified dozens of novel genes, and described the parameters of polyadenylation. We also focused on putative anthelmintic targets, identifying novel isoforms and correcting gene models. These data substantially increase the trustworthiness of gene models in these two species and demonstrate how long-read sequencing approaches should be prioritized in the continued improvement of genome assemblies and their gene annotations.

Genomic and transcriptomic variation defines the chromosome-scale assembly of Haemonchus contortus, a model gastrointestinal worm

Haemonchus contortus is a globally distributed and economically important gastrointestinal pathogen of small ruminants and has become a key nematode model for studying anthelmintic resistance and other parasite-specific traits among a wider group of parasites including major human pathogens. Here, we report using PacBio long-read and OpGen and 10X Genomics long-molecule methods to generate a highly contiguous 283.4 Mbp chromosome-scale genome assembly including a resolved sex chromosome for the MHco3(ISE).N1 isolate. We show a remarkable pattern of conservation of chromosome content with Caenorhabditis elegans, but almost no conservation of gene order. Short and long-read transcriptome sequencing allowed us to define coordinated transcriptional regulation throughout the parasite's life cycle and refine our understanding of cis- and trans-splicing. Finally, we provide a comprehensive picture of chromosome-wide genetic diversity both within a single isolate and globally. These data provide a high-quality comparison for understanding the evolution and genomics of Caenorhabditis and other nematodes and extend the experimental tractability of this model parasitic nematode in understanding helminth biology, drug discovery and vaccine development, as well as important adaptive traits such as drug resistance.

Comprehensive Chromosome End Remodeling during Programmed DNA Elimination

Germline and somatic genomes are in general the same in a multicellular organism. However, programmed DNA elimination leads to a reduced somatic genome compared to germline cells. Previous work on the parasitic nematode Ascaris demonstrated that programmed DNA elimination encompasses high-fidelity chromosomal breaks and loss of specific genome sequences including a major tandem repeat of 120 bp and ~1,000 germline-expressed genes. However, the precise chromosomal locations of these repeats, breaks regions, and eliminated genes remained unknown. We used PacBio long-read sequencing and chromosome conformation capture (Hi-C) to obtain fully assembled chromosomes of Ascaris germline and somatic genomes, enabling a complete chromosomal view of DNA elimination. We found that all 24 germline chromosomes undergo comprehensive chromosome end remodeling with DNA breaks in their subtelomeric regions and loss of distal sequences including the telomeres at both chromosome ends. All new Ascaris somatic chromosome ends are recapped by de novo telomere healing. We provide an ultrastructural analysis of Ascaris DNA elimination and show that eliminated DNA is incorporated into double membrane-bound structures, similar to micronuclei, during telophase of a DNA elimination mitosis. These micronuclei undergo dynamic changes including loss of active histone marks and localize to the cytoplasm following daughter nuclei formation and cytokinesis where they form autophagosomes. Comparative analysis of nematode chromosomes suggests that chromosome fusions occurred, forming Ascaris sex chromosomes that become independent chromosomes following DNA elimination breaks in somatic cells. These studies provide the first chromosomal view and define novel features and functions of metazoan programmed DNA elimination.

Modeling the metabolic interplay between a parasitic worm and its bacterial endosymbiont allows the identification of novel drug targets

The filarial nematode Brugia malayi represents a leading cause of disability in the developing world, causing lymphatic filariasis in nearly 40 million people. Currently available drugs are not well-suited to mass drug administration efforts, so new treatments are urgently required. One potential vulnerability is the endosymbiotic bacteria Wolbachia-present in many filariae-which is vital to the worm. Genome scale metabolic networks have been used to study prokaryotes and protists and have proven valuable in identifying therapeutic targets, but have only been applied to multicellular eukaryotic organisms more recently. Here, we present iDC625, the first compartmentalized metabolic model of a parasitic worm. We used this model to show how metabolic pathway usage allows the worm to adapt to different environments, and predict a set of 102 reactions essential to the survival of B. malayi. We validated three of those reactions with drug tests and demonstrated novel antifilarial properties for all three compounds.