Molecular epidemiology, phylogeny and evolution of the filarial nematode Wuchereria bancrofti

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.

Development and validation of a long-read metabarcoding platform for the detection of filarial worm pathogens of animals and humans

BACKGROUND

Filarial worms are important vector-borne pathogens of a large range of animal hosts, including humans, and are responsible for numerous debilitating neglected tropical diseases such as, lymphatic filariasis caused by Wuchereria bancrofti and Brugia spp., as well as loiasis caused by Loa loa. Moreover, some emerging or difficult-to-eliminate filarioid pathogens are zoonotic using animals like canines as reservoir hosts, for example Dirofilaria sp. 'hongkongensis'. Diagnosis of filariasis through commonly available methods, like microscopy, can be challenging as microfilaremia may wane below the limit of detection. In contrast, conventional PCR methods are more sensitive and specific but may show limited ability to detect coinfections as well as emerging and/or novel pathogens. Use of deep-sequencing technologies obviate these challenges, providing sensitive detection of entire parasite communities, whilst also being better suited for the characterisation of rare or novel pathogens. Therefore, we developed a novel long-read metabarcoding assay for deep-sequencing the filarial nematode cytochrome c oxidase subunit I gene on Oxford Nanopore Technologies' (ONT) MinION™ sequencer. We assessed the overall performance of our assay using kappa statistics to compare it to commonly used diagnostic methods for filarial worm detection, such as conventional PCR (cPCR) with Sanger sequencing and the microscopy-based modified Knott's test (MKT).

RESULTS

We confirmed our metabarcoding assay can characterise filarial parasites from a diverse range of genera, including, Breinlia, Brugia, Cercopithifilaria, Dipetalonema, Dirofilaria, Onchocerca, Setaria, Stephanofilaria and Wuchereria. We demonstrated proof-of-concept for this assay by using blood samples from Sri Lankan dogs, whereby we identified infections with the filarioids Acanthocheilonema reconditum, Brugia sp. Sri Lanka genotype and zoonotic Dirofilaria sp. 'hongkongensis'. When compared to traditionally used diagnostics, such as the MKT and cPCR with Sanger sequencing, we identified an additional filarioid species and over 15% more mono- and coinfections.

CONCLUSIONS

Our developed metabarcoding assay may show broad applicability for the metabarcoding and diagnosis of the full spectrum of filarioids from a wide range of animal hosts, including mammals and vectors, whilst the utilisation of ONT' small and portable MinION™ means that such methods could be deployed for field use.

Advancing Strongyloides omics data: bridging the gap with Caenorhabditis elegans

Among nematodes, the free-living model organism Caenorhabditis elegans boasts the most advanced portfolio of high-quality omics data. The resources available for parasitic nematodes, including Strongyloides spp., however, are lagging behind. While C. elegans remains the most tractable nematode and has significantly advanced our understanding of many facets of nematode biology, C. elegans is not suitable as a surrogate system for the study of parasitism and it is important that we improve the omics resources available for parasitic nematode species. Here, we review the omics data available for Strongyloides spp. and compare the available resources to those for C. elegans and other parasitic nematodes. The advancements in C. elegans omics offer a blueprint for improving omics-led research in Strongyloides. We suggest areas of priority for future research that will pave the way for expansions in omics resources and technologies. This article is part of the Theo Murphy meeting issue 'Strongyloides: omics to worm-free populations'.

Lymphatic Filariasis Elimination Status: Wuchereria bancrofti Infections in Human Populations after Five Effective Rounds of Mass Drug Administration in Zambia

Lymphatic filariasis (LF), also commonly known as elephantiasis, is a neglected tropical disease (NTD) caused by filarial parasites. The disease is transmitted via a bite from infected mosquitoes. The bites of these infected mosquitoes deposit filarial parasites, Wuchereria or Brugia, whose predilection site is the lymphatic system. The damage to the lymph system causes swelling in the legs, arms, and genitalia. A mapping survey conducted between 2003 and 2011 determined LF as being endemic in Zambia in 96 out of 116 districts. Elimination of LF is known to be possible by stopping the spread of the infection through large-scale preventive chemotherapy. Therefore, mass drug administration (MDA) with diethylcarbamazine citrate (DEC) (6 mg/kg) and Albendazole (400 mg) for Zambia has been conducted and implemented in all endemic districts with five effective rounds. In order to determine whether LF prevalence has been sufficiently reduced to levels less than 2% antigenemia and less than 1% microfilaremia, a pre-transmission assessment survey (pre-TAS) was conducted. Therefore, post-MDA pre-TAS was conducted between 2021 and 2022 in 80 districts to determine the LF prevalence. We conducted a cross-sectional seroprevalence study involving 600 participants in each evaluation unit (EU) or each district. The study sites (sentinel and spot-check sites) were from districts that were the implementation units (IUs) of the LF MDA. These included 80 districts from the 9 provinces. A total of 47,235 people from sentinel and spot-check locations were tested. Of these, valid tests were 47,052, of which 27,762 (59%) were females and 19,290 (41%) were males. The survey revealed in the 79/80 endemic districts a prevalence of Wb antigens of 0.14% and 0.0% prevalence of microfilariae. All the surveyed districts had an optimum prevalence of less than 2% for antigenaemia, except for Chibombo district. The majority of participants that tested positive for Wuchereria bancrofti (Wb) Antigens (Ag) were those that had 2, 3, and 4 rounds of MDA. Surprisingly, individuals that had 1 round of MDA were not found to have circulating antigens of Wb. The study showed that all the surveyed districts, except for Chibombo, passed pre-TAS. This further implies that there is a need to conduct transmission assessment surveys (TASs) in these districts.

First Description of the Mitogenome Features of Neofoleyellides Genus (Nematoda: Onchocercidae) Isolated from a Wild Bird (Pyrrhocorax pyrrhocorax)

Simple Summary

Filarioidea, a superfamily of nematodes, presently includes 42 species divided into six genera, mainly in the family Onchocercidae, which have been reported to infect a wide range of hosts, including reptiles, birds, and mammals. Current limitations in molecular characterization methods and species identification are the main obstacles to a better understanding of the biology of Onchocercidae species, particularly in wildlife. Thus, the objective of the present study was to sequence and analyze the complete mt genome of Neofoleyellides sp. isolated from a wild bird (Pyrrhocorax pyrrhocorax) and to assess its phylogenetic position in the Onchocercidae family. The evaluated Neofoleyellides sp. mt genome was consistent with the molecular pattern of the Onchocercidae family: 36 subunits consisting of 12 PCGs, 2 rRNAs, and 22 tRNAs. Phylogenetic analyses based on the 18S rRNA gene, cox1 gene, and 12 PCGs showed consistent results, which strongly supported monophyly of the genus Neofoleyellides. These findings enriched the gene database and improved our knowledge of the molecular characteristics of the Onchocercidae family, which provide useful genetic markers to study the population genetics, molecular biology, and phylogenetics of these Onchocercidae nematodes.

Abstract

The Onchocercidae family is composed of more than 30 valid nematode species with notable zoonotic potential. Current limitations in molecular characterization methods and species identification are the main obstacles to a better understanding of the biology of Onchocercidae species, particularly in wildlife. This study describes for the first time the complete mitochondrial (mt) genome sequence of Neofoleyellides sp. isolated from a wild bird (Pyrrhocorax pyrrhocorax) and belonging to the Neofoleyellides genus (Nematoda: Onchocercidae). The mt genome of Neofoleyellides sp. (GenBank accession number: ON641583) was a typical circular DNA molecule of 13,628 bp in size with an AT content of 76.69%. The complete mt genome comprised 36 functional subunits, including 12 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 22 transfer RNA genes. The most common start codon was ATT/ATG except for nad2 with TTG, and TAA was the termination codon for all protein-coding genes (PCGs). Phylogenetic analysis of the concatenated and aligned amino acid sequences of the 12 PCGs showed that the trees generated using different methods (Bayesian inference and maximum likelihood) with different partition schemes shared similar topologies. The isolated Neofoleyellides sp. was placed in the Onchocercidae family and formed a sister branch with the genera Onchocerca and Dirofilaria. The entire mt genome of Neofoleyellides sp. presented in this study could provide useful data for studying the population genetics and phylogenetic relationships of Onchocercidae species.

Infections and glaucoma

Glaucoma is an intraocular pressure-related ophthalmic disease with multiple causes that results in an optic neuropathy and vision loss. Intraocular pressure elevation is among its strongest risk factors. While glaucoma is mostly primary in etiology, secondary glaucoma is not infrequent. Recognizing its cause is imperative, since treatment is often different depending on the pathophysiologic mechanism. Numerous clinically relevant ophthalmic infections can result in robust inflammatory responses that may result in pressure elevation or intraocular anatomic configurations that predispose to pressure elevation. Knowing the mechanisms by which these infections can lead to glaucoma is critical in treating, and we consolidate what is currently known in regards to how infectious diseases lead to glaucoma.

Genetic epidemiology of lymphatic filariasis in American Samoa after mass drug administration

Over 892 million people in 48 countries are at risk of infection by nematodes that cause lymphatic filariasis. As part of the Global Programme to Eliminate Lymphatic Filariasis, mass drug administration is distributed to communities until surveillance indicates infection rates are below target prevalence thresholds. In some countries, including American Samoa, lymphatic filariasis transmission persists despite years of mass drug administration and/or has resurged after cessation. Nothing is known about the population genetics of Wuchereria bancrofti worms in Polynesia, or whether local transmission is persisting and/or increasing due to inadequate mass drug administration coverage, expansion from residual hotspots, reintroduction from elsewhere, or a combination. We extracted DNA from microfilariae on blood slides collected during prevalence surveys in 2014 and 2016, comprising 31 pools of five microfilariae from 22 persons living in eight villages. We sequenced 1104 bp across three mitochondrial markers (ND4, COI, CYTB). We quantified parasite genetic differentiation using variant calls and estimated haplotypes using principal components analysis, F-statistics, and haplotype networks. Of the variants called, all but eight were shared across the main island of Tutuila, and three of those were from a previously described hotspot village, Fagali'i. Genotypic data did not support population genetic structure among regions or villages in 2016, although differences were observed between worms collected in Fagali'i in 2014 and those from 2016. Because estimated haplotype frequency varied between villages, these statistics suggested genetic differentiation, but were not consistent among villages. Finally, haplotype networks demonstrated American Samoan sequence clusters were related to previously published sequences from Papua New Guinea. These are, to our knowledge, the first reports of W. bancrofti genetic variation in Polynesia. The resurgent parasites circulating on the main island of American Samoa represent a single population. This study is the first step towards investigating how parasite population structure might inform strategies to manage resurgence and elimination of lymphatic filariasis.

Genomic Epidemiology in Filarial Nematodes: Transforming the Basis for Elimination Program Decisions

Onchocerciasis and lymphatic filariasis are targeted for elimination, primarily using mass drug administration at the country and community levels. Elimination of transmission is the onchocerciasis target and global elimination as a public health problem is the end point for lymphatic filariasis. Where program duration, treatment coverage, and compliance are sufficiently high, elimination is achievable for both parasites within defined geographic areas. However, transmission has re-emerged after apparent elimination in some areas, and in others has continued despite years of mass drug treatment. A critical question is whether this re-emergence and/or persistence of transmission is due to persistence of local parasites-i.e., the result of insufficient duration or drug coverage, poor parasite response to the drugs, or inadequate methods of assessment and/or criteria for determining when to stop treatment-or due to re-introduction of parasites via human or vector movement from another endemic area. We review recent genetics-based research exploring these questions in Onchocerca volvulus, the filarial nematode that causes onchocerciasis, and Wuchereria bancrofti, the major pathogen for lymphatic filariasis. We focus in particular on the combination of genomic epidemiology and genome-wide associations to delineate transmission zones and distinguish between local and introduced parasites as the source of resurgence or continuing transmission, and to identify genetic markers associated with parasite response to chemotherapy. Our ultimate goal is to assist elimination efforts by developing easy-to-use tools that incorporate genetic information about transmission and drug response for more effective mass drug distribution, surveillance strategies, and decisions on when to stop interventions to improve sustainability of elimination.

De novo Assembly of the Brugia malayi Genome Using Long Reads from a Single MinION Flowcell

Filarial nematode infections cause a substantial global disease burden. Genomic studies of filarial worms can improve our understanding of their biology and epidemiology. However, genomic information from field isolates is limited and available reference genomes are often discontinuous. Single molecule sequencing technologies can reduce the cost of genome sequencing and long reads produced from these devices can improve the contiguity and completeness of genome assemblies. In addition, these new technologies can make generation and analysis of large numbers of field isolates feasible. In this study, we assessed the performance of the Oxford Nanopore Technologies MinION for sequencing and assembling the genome of Brugia malayi, a human parasite widely used in filariasis research. Using data from a single MinION flowcell, a 90.3 Mb nuclear genome was assembled into 202 contigs with an N50 of 2.4 Mb. This assembly covered 96.9% of the well-defined B. malayi reference genome with 99.2% identity. The complete mitochondrial genome was obtained with individual reads and the nearly complete genome of the endosymbiotic bacteria Wolbachia was assembled alongside the nuclear genome. Long-read data from the MinION produced an assembly that approached the quality of a well-established reference genome using comparably fewer resources.

Potential of cell-free DNA as a screening marker for parasite infections in dog

Parasitic infections are common in stray dogs and accurate knowledge of parasite communities in dogs would provide insight into the epidemiology of parasitic diseases. In this study, we used Illumina sequencing technology to evaluate cell-free DNA (cfDNA) as a marker for screening of parasitic infections in dogs. Plasma samples from 14 stray dogs captured in Bangladesh were used in the experiments. An average of 2.3 million reads was obtained for each sample. BLASTn analysis identified 150 reads with high similarity with parasites from 19 different genera. In particular, we detected sequences of Babesia spp. in five dogs; consistent with this, a previous study using conventional PCR showed that four of these dogs were positive for B. gibsoni. Several reads with similarity to Leishmania and filarial nematodes were also identified. These findings indicate that cfDNA in blood can be a potential screening marker for identifying parasite diversity in dogs.

Mansonella ozzardi mitogenome and pseudogene characterisation provides new perspectives on filarial parasite systematics and CO-1 barcoding

Despite the broad distribution of M. ozzardi in Latin America and the Caribbean, there is still very little DNA sequence data available to study this neglected parasite’s epidemiology. Mitochondrial DNA (mtDNA) sequences, especially the cytochrome oxidase (CO1) gene’s barcoding region, have been targeted successfully for filarial diagnostics and for epidemiological, ecological and evolutionary studies. MtDNA-based studies can, however, be compromised by unrecognised mitochondrial pseudogenes, such as Numts. Here, we have used shot-gun Illumina-HiSeq sequencing to recover the first complete Mansonella genus mitogenome and to identify several mitochondrial-origin pseudogenes. Mitogenome phylogenetic analysis placed M. ozzardi in the Onchocercidae “ONC5” clade and suggested that Mansonella parasites are more closely related to Wuchereria and Brugia genera parasites than they are to Loa genus parasites. DNA sequence alignments, BLAST searches and conceptual translations have been used to compliment phylogenetic analysis showing that M. ozzardi from the Amazon and Caribbean regions are near-identical and that previously reported Peruvian M. ozzardi CO1 reference sequences are probably of pseudogene origin. In addition to adding a much-needed resource to the Mansonella genus’s molecular tool-kit and providing evidence that some M. ozzardi CO1 sequence deposits are pseudogenes, our results suggest that all Neotropical M. ozzardi parasites are closely related.

Genetic exchange in eukaryotes through horizontal transfer: connected by the mobilome

Background

All living species contain genetic information that was once shared by their common ancestor. DNA is being inherited through generations by vertical transmission (VT) from parents to offspring and from ancestor to descendant species. This process was considered the sole pathway by which biological entities exchange inheritable information. However, Horizontal Transfer (HT), the exchange of genetic information by other means than parents to offspring, was discovered in prokaryotes along with strong evidence showing that it is a very important process by which prokaryotes acquire new genes.

Main body

For some time now, it has been a scientific consensus that HT events were rare and non-relevant for evolution of eukaryotic species, but there is growing evidence supporting that HT is an important and frequent phenomenon in eukaryotes as well.

Conclusion

Here, we will discuss the latest findings regarding HT among eukaryotes, mainly HT of transposons (HTT), establishing HTT once and for all as an important phenomenon that should be taken into consideration to fully understand eukaryotes genome evolution. In addition, we will discuss the latest development methods to detect such events in a broader scale and highlight the new approaches which should be pursued by researchers to fill the knowledge gaps regarding HTT among eukaryotes.

The drug target genes show higher evolutionary conservation than non-target genes

Although evidence indicates that drug target genes share some common evolutionary features, there have been few studies analyzing evolutionary features of drug targets from an overall level. Therefore, we conducted an analysis which aimed to investigate the evolutionary characteristics of drug target genes. We compared the evolutionary conservation between human drug target genes and non-target genes by combining both the evolutionary features and network topological properties in human protein-protein interaction network. The evolution rate, conservation score and the percentage of orthologous genes of 21 species were included in our study. Meanwhile, four topological features including the average shortest path length, betweenness centrality, clustering coefficient and degree were considered for comparison analysis. Then we got four results as following: compared with non-drug target genes, 1) drug target genes had lower evolutionary rates; 2) drug target genes had higher conservation scores; 3) drug target genes had higher percentages of orthologous genes and 4) drug target genes had a tighter network structure including higher degrees, betweenness centrality, clustering coefficients and lower average shortest path lengths. These results demonstrate that drug target genes are more evolutionarily conserved than non-drug target genes. We hope that our study will provide valuable information for other researchers who are interested in evolutionary conservation of drug targets.

Ancient horizontal transfers of retrotransposons between birds and ancestors of human pathogenic nematodes

Parasite host switches may trigger disease emergence, but prehistoric host ranges are often unknowable. Lymphatic filariasis and loiasis are major human diseases caused by the insect-borne filarial nematodes Brugia, Wuchereria and Loa. Here we show that the genomes of these nematodes and seven tropical bird lineages exclusively share a novel retrotransposon, AviRTE, resulting from horizontal transfer (HT). AviRTE subfamilies exhibit 83–99% nucleotide identity between genomes, and their phylogenetic distribution, paleobiogeography and invasion times suggest that HTs involved filarial nematodes. The HTs between bird and nematode genomes took place in two pantropical waves, >25–22 million years ago (Myr ago) involving the Brugia/Wuchereria lineage and >20–17 Myr ago involving the Loa lineage. Contrary to the expectation from the mammal-dominated host range of filarial nematodes, we hypothesize that these major human pathogens may have independently evolved from bird endoparasites that formerly infected the global breadth of avian biodiversity.

Lymphatic filariasis and loiasis are diseases caused by insect-borne filarial nematodes. Here, Suh et al. identify a retrotransposon that is present in the genomes of these nematodes and seven tropical bird lineages, indicating two waves of horizontal gene transfer around 17–25 million years ago.