Tissue-specific transcriptomics and proteomics of a filarial nematode and its Wolbachia endosymbiont

Transcriptome profiling of male and female Ascaris lumbricoides reproductive tissues

BACKGROUND

Ascaris lumbricoides causes human ascariasis, the most prevalent helminth disease, infecting approximately 1 billion individuals globally. In 2019 the global disease burden was estimated to be 754,000 DALYs and resulted in 2090 deaths. In the absence of a vaccination strategy, treatment of ascariasis has relied on anthelminthic chemotherapy, but drug resistance is a concern. The propensity for reinfection is also a major challenge to disease control; female worms lay up to 200,000 eggs daily, which contaminate surrounding environments and remain viable for years, resulting in high transmission rates. Understanding the molecular mechanisms of reproductive processes, including control of egg production, spermatogenesis, oogenesis and embryogenesis, will drive the development of new drugs and/or vaccine targets for future ascariasis control.

METHODS

Transcriptome profiles of discrete reproductive and somatic tissue samples were generated from adult male and female worms using Illumina HiSeq with 2 × 150 bp paired-end sequencing. Male tissues included: testis germinal zone, testis part of vas deferens, seminal vesicle and somatic tissue. Female tissues included: ovary germinal zone, ovary part of the oviduct, uterus and somatic tissue. Differentially expressed genes (DEGs) were identified from the fragments per kilobases per million reads (FPKM) profiles. Hierarchical analysis was performed to identify tissue-specific genes. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were employed to identify significant terms and pathways for the DEGs.

RESULTS

DEGs involved in protein phosphorylation and adhesion molecules were indicated to play a crucial role in spermatogenesis and fertilization, respectively. Those genes associated with the G-protein-coupled receptor (GPCR) signaling pathway and small GTPase-mediated signal transduction pathway play an essential role in cytoskeleton organization during oogenesis. Additionally, DEGs associated with the SMA genes and TGF-β signaling pathway are crucial in adult female embryogenesis. Some genes associated with particular biological processes and pathways that were identified in this study have been linked to defects in germline development, embryogenesis and reproductive behavior. In the enriched KEGG pathway analysis, Hippo signaling, oxytocin signaling and tight junction pathways were identified to play a role in Ascaris male and female reproductive systems.

CONCLUSIONS

This study has provided comprehensive transcriptome profiles of discrete A. lumbricoides reproductive tissue samples, revealing the molecular basis of these functionally important tissues. The data generated from this study will provide fundamental knowledge on the reproductive biology of Ascaris and will inform future target identification for anti-ascariasis drugs and/or vaccines.

Identification and profiling of Trichinella spiralis circulating antigens and proteins in sera of mice with trichinellosis

Trichinellosis is a zoonotic disease caused by the ingestion of the Trichinella nematode. With a worldwide incidence of approximately 10,000 cases per year, Trichinella spiralis is responsible for most human infections. There are no specific signs or symptoms of this parasitic infection. Muscle biopsy is the gold diagnostic standard for trichinellosis, but the technique is invasive and unable to detect the early stage of infection. Although immunodiagnostics are also available, antibody detection usually occurs after 3 weeks and prolonged up to 19 years after the acute phase. Therefore, additional diagnostic biomarkers must be identified to improve trichinellosis diagnosis. This study aimed to measure concentration changes in mouse serum proteins prior to T. spiralis infection and 2, 4 and 8 weeks after infection, and to identify T. spiralis circulating proteins and antigens using mass spectrometry-based proteomics. Mouse muscle-related proteins including inter-alpha-trypsin inhibitor heavy chain H2, a protein involved in the response to muscle tissue damage, were up-regulated in mouse sera during the T. spiralis larvae invasion. Additionally, 33 circulatory parasite proteins were identified in infected mouse sera. Notably, T. spiralis long-chain fatty acid transport protein 1 could be detected in the early stage of infection and peroxidasin-like protein was identified 2, 4 and 8 weeks after infection. Seventeen T. spiralis circulating antigens were detected in mouse immune complexes, with PX domain protein being found 2, 4 and 8 weeks after infection. Because peroxidasin-like protein and PX domain protein were detected at all post-infection time points, sequence alignments of these proteins were performed, which showed they are conserved among Trichinella spp. and have less similarity to the human and murine sequences. Integrative analysis of T. spiralis biomarkers throughout the course of infection may reveal additional diagnostic targets to improve early diagnosis of trichinellosis.

Heartworm disease - Overview, intervention, and industry perspective

Dirofilaria immitis, also known as heartworm, is a major parasitic threat for dogs and cats around the world. Because of its impact on the health and welfare of companion animals, heartworm disease is of huge veterinary and economic importance especially in North America, Europe, Asia and Australia. Within the animal health market many different heartworm preventive products are available, all of which contain active components of the same drug class, the macrocyclic lactones. In addition to compliance issues, such as under-dosing or irregular treatment intervals, the occurrence of drug-resistant heartworms within the populations in the Mississippi River areas adds to the failure of preventive treatments. The objective of this review is to provide an overview of the disease, summarize the current disease control measures and highlight potential new avenues and best practices for treatment and prevention.

Immunoproteomic Analysis of Dirofilaria repens Microfilariae and Adult Parasite Stages

Dirofilariarepens is a parasitic nematode causing a vector-borne zoonotic infection (dirofilariosis), considered an emerging problem in human and veterinary medicine. Currently, diagnosis is based on the detection of the adult parasite and microfilariae in the host tissues. However, the efficacy of tests relying on microfilariae detection is limited by microfilariae periodic occurrence. Therefore, a new reliable and affordable serological diagnostic method is needed. Better characteristic of the parasite biology and its interaction with host immune system should help to achieve this goal. This study analyzes adult and microfilariae proteomes, and the use of one-dimensional electrophoresis (1-DE) and two-dimensional electrophoresis (2-DE) proteomics, immunoproteomics, and LC-MS/MS mass spectrometry allowed us to identify 316 potentially immunogenic proteins (75 belong to adult stage, 183 to microfilariae, and 58 are common for both). Classified by their ontology, the proteins showed important similarities and differences between both parasite stages. The most frequently identified proteins are structural, metabolic, and heat shock proteins. Additionally, real-time PCR analysis of some immunogenic targets revealed significant differences between microfilariae and adult life stages. We indicated molecules involved in parasite-host interactions and discussed their importance in parasite biology, which may help to reveal potential diagnostic antigens or select drug and vaccine targets.

Diminutive, degraded but dissimilar: Wolbachia genomes from filarial nematodes do not conform to a single paradigm

Wolbachia are alpha-proteobacteria symbionts infecting a large range of arthropod species and two different families of nematodes. Interestingly, these endosymbionts are able to induce diverse phenotypes in their hosts: they are reproductive parasites within many arthropods, nutritional mutualists within some insects and obligate mutualists within their filarial nematode hosts. Defining Wolbachia 'species' is controversial and so they are commonly classified into 17 different phylogenetic lineages, termed supergroups, named A-F, H-Q and S. However, available genomic data remain limited and not representative of the full Wolbachia diversity; indeed, of the 24 complete genomes and 55 draft genomes of Wolbachia available to date, 84 % belong to supergroups A and B, exclusively composed of Wolbachia from arthropods. For the current study, we took advantage of a recently developed DNA-enrichment method to produce four complete genomes and two draft genomes of Wolbachia from filarial nematodes. Two complete genomes, wCtub and wDcau, are the smallest Wolbachia genomes sequenced to date (863 988 bp and 863 427 bp, respectively), as well as the first genomes representing supergroup J. These genomes confirm the validity of this supergroup, a controversial clade due to weaknesses of the multilocus sequence typing approach. We also produced the first draft Wolbachia genome from a supergroup F filarial nematode representative (wMhie), two genomes from supergroup D (wLsig and wLbra) and the complete genome of wDimm from supergroup C. Our new data confirm the paradigm of smaller Wolbachia genomes from filarial nematodes containing low levels of transposable elements and the absence of intact bacteriophage sequences, unlike many Wolbachia from arthropods, where both are more abundant. However, we observe differences among the Wolbachia genomes from filarial nematodes: no global co-evolutionary pattern, strong synteny between supergroup C and supergroup J Wolbachia, and more transposable elements observed in supergroup D Wolbachia compared to the other supergroups. Metabolic pathway analysis indicates several highly conserved pathways (haem and nucleotide biosynthesis, for example) as opposed to more variable pathways, such as vitamin B biosynthesis, which might be specific to certain host-symbiont associations. Overall, there appears to be no single Wolbachia-filarial nematode pattern of co-evolution or symbiotic relationship.

A Meta-Analysis of Wolbachia Transcriptomics Reveals a Stage-Specific Wolbachia Transcriptional Response Shared Across Different Hosts

Wolbachia is a genus containing obligate, intracellular endosymbionts with arthropod and nematode hosts. Numerous studies have identified differentially expressed transcripts in Wolbachia endosymbionts that potentially inform the biological interplay between these endosymbionts and their hosts, albeit with discordant results. Here, we re-analyze previously published Wolbachia RNA-Seq transcriptomics data sets using a single workflow consisting of the most up-to-date algorithms and techniques, with the aim of identifying trends or patterns in the pan-Wolbachia transcriptional response. We find that data from one of the early studies in filarial nematodes did not allow for robust conclusions about Wolbachia differential expression with these methods, suggesting the original interpretations should be reconsidered. Across datasets analyzed with this unified workflow, there is a general lack of global gene regulation with the exception of a weak transcriptional response resulting in the upregulation of ribosomal proteins in early larval stages. This weak response is observed across diverse Wolbachia strains from both nematode and insect hosts suggesting a potential pan-Wolbachia transcriptional response during host development that diverged more than 700 million years ago.

Elucidating the molecular and developmental biology of parasitic nematodes: Moving to a multiomics paradigm

In the past two decades, significant progress has been made in the sequencing, assembly, annotation and analyses of genomes and transcriptomes of parasitic worms of socioeconomic importance. This progress has somewhat improved our knowledge and understanding of these pathogens at the molecular level. However, compared with the free-living nematode Caenorhabditis elegans, the areas of functional genomics, transcriptomics, proteomics and metabolomics of parasitic nematodes are still in their infancy, and there are major gaps in our knowledge and understanding of the molecular biology of parasitic nematodes. The information on signalling molecules, molecular pathways and microRNAs (miRNAs) that are known to be involved in developmental processes in C. elegans and the availability of some molecular resources (draft genomes, transcriptomes and some proteomes) for selected parasitic nematodes provide a basis to start exploring the developmental biology of parasitic nematodes. Indeed, some studies have identified molecules and pathways that might associate with developmental processes in related, parasitic nematodes, such as Haemonchus contortus (barber's pole worm). However, detailed information is often scant and 'omics resources are limited, preventing a proper integration of 'omic data sets and comprehensive analyses. Moreover, little is known about the functional roles of pheromones, hormones, signalling pathways and post-transcriptional/post-translational regulations in the development of key parasitic nematodes throughout their entire life cycles. Although C. elegans is an excellent model to assist molecular studies of parasitic nematodes, its use is limited when it comes to explorations of processes that are specific to parasitism within host animals. A deep understanding of parasitic nematodes, such as H. contortus, requires substantially enhanced resources and the use of integrative 'omics approaches for analyses. The improved genome and well-established in vitro larval culture system for H. contortus provide unprecedented opportunities for comprehensive studies of the transcriptomes (mRNA and miRNA), proteomes (somatic, excretory/secretory and phosphorylated proteins) and lipidomes (e.g., polar and neutral lipids) of this nematode. Such resources should enable in-depth explorations of its developmental biology at a level, not previously possible. The main aims of this review are (i) to provide a background on the development of nematodes, with a particular emphasis on the molecular aspects involved in the dauer formation and exit in C. elegans; (ii) to critically appraise the current state of knowledge of the developmental biology of parasitic nematodes and identify key knowledge gaps; (iii) to cover salient aspects of H. contortus, with a focus on the recent advances in genomics, transcriptomics, proteomics and lipidomics as well as in vitro culturing systems; (iv) to review recent advances in our knowledge and understanding of the molecular and developmental biology of H. contortus using an integrative multiomics approach, and discuss the implications of this approach for detailed explorations of signalling molecules, molecular processes and pathways likely associated with nematode development, adaptation and parasitism, and for the identification of novel intervention targets against these pathogens. Clearly, the multiomics approach established recently is readily applicable to exploring a wide range of interesting and socioeconomically significant parasitic worms (including also trematodes and cestodes) at the molecular level, and to elucidate host-parasite interactions and disease processes.

Comparative Genomics Guides Elucidation of Vitamin B12 Biosynthesis in Novel Human-Associated Akkermansia Strains

There is significant interest in the therapeutic and probiotic potential of the common gut bacterium Akkermansia muciniphila. However, knowledge of both the genomic and physiological diversity of this bacterial lineage is limited. Using a combination of genomic, molecular biological, and traditional microbiological approaches, we identified at least four species-level phylogroups with differing functional potentials that affect how these bacteria interact with both their human host and other members of the human gut microbiome. Specifically, we identified and isolated Akkermansia strains that were able to synthesize vitamin B₁₂. The ability to synthesize this important cofactor broadens the physiological capabilities of human-associated Akkermansia strains, fundamentally altering our understanding of how this important bacterial lineage may affect human health.

Akkermansia muciniphila is a mucin-degrading bacterium found in the gut of most humans and is considered a “next-generation probiotic.” However, knowledge of the genomic and physiological diversity of human-associated Akkermansia sp. strains is limited. Here, we reconstructed 35 metagenome-assembled genomes and combined them with 40 publicly available genomes for comparative genomic analysis. We identified at least four species-level phylogroups (AmI to AmIV), with distinct functional potentials. Most notably, we identified genes for cobalamin (vitamin B₁₂) biosynthesis within the AmII and AmIII phylogroups. To verify these predictions, 10 Akkermansia strains were isolated from adults and screened for vitamin B₁₂ biosynthesis genes via PCR. Two AmII strains were positive for the presence of cobalamin biosynthesis genes, while all 9 AmI strains tested were negative. To demonstrate vitamin B₁₂ biosynthesis, we measured the production of acetate, succinate, and propionate in the presence and absence of vitamin supplementation in representative strains of the AmI and AmII phylogroups, since cobalamin is an essential cofactor in propionate metabolism. Results showed that the AmII strain produced acetate and propionate in the absence of supplementation, which is indicative of vitamin B₁₂ biosynthesis. In contrast, acetate and succinate were the main fermentation products for the AmI strains when vitamin B₁₂ was not supplied in the culture medium. Lastly, two bioassays were used to confirm vitamin B₁₂ production by the AmII phylogroup. This novel physiological trait of human-associated Akkermansia strains may affect how these bacteria interact with the human host and other members of the human gut microbiome.

IMPORTANCE There is significant interest in the therapeutic and probiotic potential of the common gut bacterium Akkermansia muciniphila. However, knowledge of both the genomic and physiological diversity of this bacterial lineage is limited. Using a combination of genomic, molecular biological, and traditional microbiological approaches, we identified at least four species-level phylogroups with differing functional potentials that affect how these bacteria interact with both their human host and other members of the human gut microbiome. Specifically, we identified and isolated Akkermansia strains that were able to synthesize vitamin B₁₂. The ability to synthesize this important cofactor broadens the physiological capabilities of human-associated Akkermansia strains, fundamentally altering our understanding of how this important bacterial lineage may affect human health.

Defining the target and the effect of imatinib on the filarial c-Abl homologue

BACKGROUND

Previously we demonstrated the micro- and macrofilaricidal properties of imatinib in vitro. Here we use electron and multiphoton microscopy to define the target of imatinib in the adult and microfilarial stages of Brugia malayi and assess the effects of pharmacologically relevant levels of imatinib on the adult parasites.

METHODS

After fixation of adult B. malayi males and females, sections were stained with polyclonal rabbit anti-c-Abl antibody (or isotype control) and imaged with multiphoton fluorescent microscopy. Microfilariae were fixed and labeled with rabbit anti-c-Abl IgG primary antibody followed by anti-rabbit gold conjugated secondary antibody and imaged using transmission electron microscopy (TEM; immunoEM). In addition, adult B. malayi males and females were exposed to 0 or 10μM of imatinib for 7 days following which they were prepared for transmission electron microscopy (TEM) to assess the drug's effect on filarial ultrastructure.

RESULTS

Fluorescent localization of anti-c-Abl antibody demonstrated widespread uptake in the adult filariae, but the most intense signal was seen in the reproductive organs, muscle, and intestine of both male and female worms. Fluorescence was significantly more intense in the early microfilarial stage (i.e. early morula) compared with later development stages (i.e. pretzel). Anti-c-Abl antibody in the microfilariae localized to the nuclei. Based on TEM assessment following imatinib exposure, imatinib appeared to be detrimental to embryogenesis in the adult female B. malayi.

CONCLUSIONS

At pharmacologically achievable concentrations of imatinib, embryogenesis is impaired and possibly halted in adult filariae. Imatinib is likely a slow microfilaricide due to interference in intra-nuclear processes, which are slowly detrimental to the parasite and not immediately lethal, and thus may be used to lower the levels of L. loa microfilariae before they are treated within the context of conventional mass drug administration.

Defining Brugia malayi and Wolbachia symbiosis by stage-specific dual RNA-seq

Background

Filarial nematodes currently infect up to 54 million people worldwide, with millions more at risk for infection, representing the leading cause of disability in the developing world. Brugia malayi is one of the causative agents of lymphatic filariasis and remains the only human filarial parasite that can be maintained in small laboratory animals. Many filarial nematode species, including B. malayi, carry an obligate endosymbiont, the alpha-proteobacteria Wolbachia, which can be eliminated through antibiotic treatment. Elimination of the endosymbiont interferes with development, reproduction, and survival of the worms within the mamalian host, a clear indicator that the Wolbachia are crucial for survival of the parasite. Little is understood about the mechanism underlying this symbiosis.

Methodology/ Principle findings

To better understand the molecular interplay between these two organisms we profiled the transcriptomes of B. malayi and Wolbachia by dual RNA-seq across the life cycle of the parasite. This helped identify functional pathways involved in this essential symbiotic relationship provided by the co-expression of nematode and bacterial genes. We have identified significant stage-specific and gender-specific differential expression in Wolbachia during the nematode’s development. For example, during female worm development we find that Wolbachia upregulate genes involved in ATP production and purine biosynthesis, as well as genes involved in the oxidative stress response.

Conclusions/ Significance

This global transcriptional analysis has highlighted specific pathways to which both Wolbachia and B. malayi contribute concurrently over the life cycle of the parasite, paving the way for the development of novel intervention strategies.

Filarial nematodes currently infect millions of people worldwide and represent a leading cause of disability. Currently available medications are insufficient in reaching elimination of these parasites. Many filarial nematodes, including Brugia malayi, have an Achilles heel of sorts—that is their obligate symbiotic relationship with the bacteria Wolbachia. While it is known that the nematode and the bacteria are co-dependent, the molecular basis of this relationship remains poorly understood. Using deep sequencing, we profiled the transcriptomes of B. malayi and Wolbachia across the life cycle of the parasite to determine the functional pathways necessary for parasite survival provided by the co-expression of nematode and bacterial genes. Defining the mechanisms of endosymbiosis between these two organisms will allow for the exploitation of this relationship for the development of new intervention strategies.

Removing the needle from the haystack: Enrichment of Wolbachia endosymbiont transcripts from host nematode RNA by Cappable-seq™

Efficient transcriptomic sequencing of microbial mRNA derived from host-microbe associations is often compromised by the much lower relative abundance of microbial RNA in the mixed total RNA sample. One solution to this problem is to perform extensive sequencing until an acceptable level of transcriptome coverage is obtained. More cost-effective methods include use of prokaryotic and/or eukaryotic rRNA depletion strategies, sometimes in conjunction with depletion of polyadenylated eukaryotic mRNA. Here, we report use of Cappable-seq™ to specifically enrich, in a single step, Wolbachia endobacterial mRNA transcripts from total RNA prepared from the parasitic filarial nematode, Brugia malayi. The obligate Wolbachia endosymbiont is a proven drug target for many human filarial infections, yet the precise nature of its symbiosis with the nematode host is poorly understood. Insightful analysis of the expression levels of Wolbachia genes predicted to underpin the mutualistic association and of known drug target genes at different life cycle stages or in response to drug treatments is typically challenged by low transcriptomic coverage. Cappable-seq resulted in up to ~ 5-fold increase in the number of reads mapping to Wolbachia. On average, coverage of Wolbachia transcripts from B. malayi microfilariae was enriched ~40-fold by Cappable-seq. Additionally, this method has an additional benefit of selectively removing abundant prokaryotic ribosomal RNAs.The deeper microbial transcriptome sequencing afforded by Cappable-seq facilitates more detailed characterization of gene expression levels of pathogens and symbionts present in animal tissues.

Lessons from the genomes and transcriptomes of filarial nematodes

Human filarial infections are a leading cause of morbidity in the developing world. While a small arsenal of drugs exists to treat these infections, there remains a tremendous need for the development of additional interventions. Recent genome sequences and transcriptome analyses of filarial nematodes have provided novel biological insight and allowed for the prediction of novel drug targets as well as potential vaccine candidates. In this review, we discuss the currently available data, insights gained into the metabolism of these organisms, and how the filaria field can move forward by leveraging these data.

Perusal of parasitic nematode 'omics in the post-genomic era

The advent of high-throughput, next-generation sequencing methods combined with advances in computational biology and bioinformatics have greatly accelerated discovery within biomedical research. This "post-genomics" era has ushered in powerful approaches allowing one to quantify RNA transcript and protein abundance for every gene in the genome - often for multiple conditions. Herein, we chronicle how the post-genomics era has advanced our overall understanding of parasitic nematodes through transcriptomics and proteomics and highlight some of the important advances made in each major nematode clade. We primarily focus on organisms relevant to human health, given that nematode infections significantly impact disability-adjusted life years (DALY) scores within the developing world, but we also discuss organisms of veterinary importance as well as those used as laboratory models. As such, we envision that this review will serve as a comprehensive resource for those seeking a better understanding of basic parasitic nematode biology as well as those interested in targets for vaccination and pharmacological intervention.

Heme acquisition in the parasitic filarial nematode Brugia malayi

Nematodes lack a heme biosynthetic pathway and must acquire heme from exogenous sources. Given the indispensable role of heme, this auxotrophy may be exploited to develop drugs that interfere with heme uptake in parasites. Although multiple heme-responsive genes (HRGs) have been characterized within the free-living nematode Caenorhabditis elegans, we have undertaken the first study of heme transport in Brugia malayi, a causative agent of lymphatic filariasis. Through functional assays in yeast, as well as heme analog, RNAi, and transcriptomic experiments, we have shown that the heme transporter B. malayi HRG-1 (BmHRG-1) is indeed functional in B. malayi In addition, BmHRG-1 localizes both to the endocytic compartments and cell membrane when expressed in yeast cells. Transcriptomic sequencing revealed that BmHRG-1, BmHRG-2, and BmMRP-5 (all orthologs of HRGs in C. elegans) are down-regulated in heme-treated B. malayi, as compared to non-heme-treated control worms. Likely because of short gene lengths, multiple exons, other HRGs in B. malayi (BmHRG-3-6) remain unidentified. Although the precise mechanisms of heme homeostasis in a nematode with the ability to acquire heme remains unknown, this study clearly demonstrates that the filarial nematode B. malayi is capable of transporting exogenous heme.-Luck, A. N., Yuan, X., Voronin, D., Slatko, B. E., Hamza, I., Foster, J. M. Heme acquisition in the parasitic filarial nematode Brugia malayi.