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

Deciphering the pH-dependent oligomerization of aspartate semialdehyde dehydrogenase from Wolbachia endosymbiont of Brugia malayi: An in vitro and in silico approaches

The enzyme aspartate semialdehyde dehydrogenase (ASDH) plays a pivotal role in the amino acid biosynthesis pathway, making it an attractive target for the development of new antimicrobial drugs due to its absence in humans. This study aims to investigate the presence of ASDH in the filarial parasite Wolbachia endosymbiont of Brugia malayi (WBm) using both in vitro and in silico approaches. The size exclusion chromatography (SEC) and Native-PAGE analysis demonstrate that WBm-ASDH undergoes pH-dependent oligomerization and dimerization. To gain a deeper understanding of this phenomenon, the modelled monomer and dimer structures were subjected to pH-dependent dynamics simulations in various conditions. The results reveal that residues Val240, Gln161, Thr159, Tyr160, and Trp316 form strong hydrogen bond contacts in the intersurface area to maintain the structure in the dimeric form. Furthermore, the binding of NADP+ induces conformational changes, leading to an open or closed conformation in the structure. Importantly, the binding of NADP+ does not disturb either the dimerization or oligomerization of the protein, a finding confirmed through both in vitro and in silico analysis. These findings shed light on the structural characteristics of WBm-ASDH and offer valuable insights for the development of new inhibitors specific to WBm, thereby contributing to the development of potential therapies for filarial parasitic infections.

Updated annotation and meta-analysis of Brugia malayi transcriptomics data reveals consistent transcriptional profiles across time and space with some study-specific differences in adult female worm transcriptional profiles

Genomics, transcriptomics, and proteomics have significantly advanced our understanding of obligately host-associated microbes, where interrogation of the biology is often limited by the complexity of the biological system and limited tools. This includes the causative agents of many neglected tropical diseases, including filarial nematodes. Therefore, numerous transcriptomics studies have been undertaken on filarial nematodes. Most of these transcriptomics studies focus on Brugia malayi, which causes lymphatic filariasis and is a laboratory model for human filarial disease. Here, we undertook a meta-analysis of the publicly available B. malayi transcriptomics data enabling the direct cross comparison of samples from almost a dozen studies. This reanalysis highlights the consistency of transcriptomics results across many different studies and experimental designs from across the globe for over a decade of research, across many different generations of a sequencing technology, library preparation protocols, and differential expression tools. Males and microfilariae across samples had similar expression profiles. However, female samples were clustered into two differential expression patterns that were significantly different from one another. Largely, we confirm previous results for all studies reanalyzed including tissue-specific gene expression and anti-Wolbachia doxycycline treatment of microfilaria. However, we did not detect previously reported differential expression upon in vitro or in vivo treatment with ivermectin, albendazole, and DEC, instead identifying a consistent lack of transcriptomic change upon exposure to these anthelminthic drugs. Updated annotation has been provided that denotes poorly supported genes including those overlapping rRNAs.

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.

Improved metagenome assemblies through selective enrichment of bacterial genomic DNA from eukaryotic host genomic DNA using ATAC-seq

Genomics can be used to study the complex relationships between hosts and their microbiota. Many bacteria cannot be cultured in the laboratory, making it difficult to obtain adequate amounts of bacterial DNA and to limit host DNA contamination for the construction of metagenome-assembled genomes (MAGs). For example, Wolbachia is a genus of exclusively obligate intracellular bacteria that live in a wide range of arthropods and some nematodes. While Wolbachia endosymbionts are frequently described as facultative reproductive parasites in arthropods, the bacteria are obligate mutualistic endosymbionts of filarial worms. Here, we achieve 50-fold enrichment of bacterial sequences using ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) with Brugia malayi nematodes, containing Wolbachia (wBm). ATAC-seq uses the Tn5 transposase to cut and attach Illumina sequencing adapters to accessible DNA lacking histones, typically thought to be open chromatin. Bacterial and mitochondrial DNA in the lysates are also cut preferentially since they lack histones, leading to the enrichment of these sequences. The benefits of this include minimal tissue input (<1 mg of tissue), a quick protocol (<4 h), low sequencing costs, less bias, correct assembly of lateral gene transfers and no prior sequence knowledge required. We assembled the wBm genome with as few as 1 million Illumina short paired-end reads with >97% coverage of the published genome, compared to only 12% coverage with the standard gDNA libraries. We found significant bacterial sequence enrichment that facilitated genome assembly in previously published ATAC-seq data sets from human cells infected with Mycobacterium tuberculosis and C. elegans contaminated with their food source, the OP50 strain of E. coli. These results demonstrate the feasibility and benefits of using ATAC-seq to easily obtain bacterial genomes to aid in symbiosis, infectious disease, and microbiome research.

Miniature spatial transcriptomics for studying parasite-endosymbiont relationships at the micro scale

Several important human infectious diseases are caused by microscale-sized parasitic nematodes like filarial worms. Filarial worms have their own spatial tissue organization; to uncover this tissue structure, we need methods that can spatially resolve these miniature specimens. Most filarial worms evolved a mutualistic association with endosymbiotic bacteria Wolbachia. However, the mechanisms underlying the dependency of filarial worms on the fitness of these bacteria remain unknown. As Wolbachia is essential for the development, reproduction, and survival of filarial worms, we spatially explored how Wolbachia interacts with the worm's reproductive system by performing a spatial characterization using Spatial Transcriptomics (ST) across a posterior region containing reproductive tissue and developing embryos of adult female Brugia malayi worms. We provide a proof-of-concept for miniature-ST to explore spatial gene expression patterns in small sample types, demonstrating the method's ability to uncover nuanced tissue region expression patterns, observe the spatial localization of key B. malayi - Wolbachia pathway genes, and co-localize the B. malayi spatial transcriptome in Wolbachia tissue regions, also under antibiotic treatment. We envision our approach will open up new avenues for the study of infectious diseases caused by micro-scale parasitic worms.

Dual RNA-seq of spleens extracted from channel catfish infected with Aeromonas veronii reveals novel insights into host-pathogen interactions

Interactions between host and pathogen are involving various dynamic changes in transcript expression and critical for understanding host immunity against infections and its associated pathogenesis. Herein, we established a model of channel catfish infected with Aeromonas veronii. The infected fish had prominent body surface bleeding, and the spleen showed hyperemia and swelling. Then, the spleen of channel catfish infected with A. veronii was analyzed by dual RNA sequencing (RNA-seq), and the transcriptome data were compared with uninfected channel catfish spleen or bacteria cultured in vitro. The transcript expression profile of pathogen-host interaction between A. veronii and channel catfish was successfully studied. During infection, the host was enriched for multiple immune-related signaling pathways, such as the Toll-like receptor signaling pathway, Cytokine-cytokine receptor interaction, and T cell receptor signaling pathway; and significantly upregulated for many innate immune-related genes, including IL-8. At the same time, we found that A. veronii mainly harmed the host spleen through hemolysin. Our current findings are of great significance in clarifying the pathogenesis of channel catfish induced by A. veronii and provide gene targets for developing preventive measures.

Ribosomal RNA-Depletion Provides an Efficient Method for Successful Dual RNA-Seq Expression Profiling of a Marine Sponge Holobiont

Investigations of host-symbiont interactions can benefit enormously from a complete and reliable holobiont gene expression profiling. The most efficient way to acquire holobiont transcriptomes is to perform RNA-Seq on both host and symbionts simultaneously. However, optimal methods for capturing both host and symbiont mRNAs are still under development, particularly when the host is a eukaryote and the symbionts are bacteria or archaea. Traditionally, poly(A)-enriched libraries have been used to capture eukaryotic mRNA, but the ability of this method to adequately capture bacterial mRNAs is unclear because of the short half-life of the bacterial transcripts. Here, we address this gap in knowledge with the aim of helping others to choose an appropriate RNA-Seq approach for analysis of animal host-bacterial symbiont transcriptomes. Specifically, we compared transcriptome bias, depth and coverage achieved by two different mRNA capture and sequencing strategies applied to the marine demosponge Amphimedon queenslandica holobiont. Annotated genomes of the sponge host and the three most abundant bacterial symbionts, which can comprise up to 95% of the adult microbiome, are available. Importantly, this allows for transcriptomes to be accurately mapped to these genomes, and thus quantitatively assessed and compared. The two strategies that we compare here are (i) poly(A) captured mRNA-Seq (Poly(A)-RNA-Seq) and (ii) ribosomal RNA depleted RNA-Seq (rRNA-depleted-RNA-Seq). For the host sponge, we find no significant difference in transcriptomes generated by the two different mRNA capture methods. However, for the symbiont transcriptomes, we confirm the expectation that the rRNA-depleted-RNA-Seq performs much better than the Poly(A)-RNA-Seq. This comparison demonstrates that RNA-Seq by ribosomal RNA depletion is an effective and reliable method to simultaneously capture gene expression in host and symbionts and thus to analyse holobiont transcriptomes.

Aspartyl Protease Inhibitors as Anti-Filarial Drugs

The current treatments for lymphatic filariasis and onchocerciasis do not effectively kill the adult parasitic nematodes, allowing these chronic and debilitating diseases to persist in millions of people. Thus, the discovery of new drugs with macrofilaricidal potential to treat these filarial diseases is critical. To facilitate this need, we first investigated the effects of three aspartyl protease inhibitors (APIs) that are FDA-approved as HIV antiretroviral drugs on the adult filarial nematode, Brugia malayi and the endosymbiotic bacteria, Wolbachia. From the three hits, nelfinavir had the best potency with an IC₅₀ value of 7.78 µM, followed by ritonavir and lopinavir with IC₅₀ values of 14.3 µM and 16.9 µM, respectively. The three APIs have a direct effect on killing adult B. malayi after 6 days of exposure in vitro and did not affect the Wolbachia titers. Sequence conservation and stage-specific gene expression analysis identified Bm8660 as the most likely primary aspartic protease target for these drug(s). Immunolocalization using antibodies raised against the Bm8660 ortholog of Onchocerca volvulus showed it is strongly expressed in female B. malayi, especially in metabolically active tissues such as lateral and dorsal/ventral chords, hypodermis, and uterus tissue. Global transcriptional response analysis using adult female B. pahangi treated with APIs identified four additional aspartic proteases differentially regulated by the three effective drugs, as well as significant enrichment of various pathways including ubiquitin mediated proteolysis, protein kinases, and MAPK/AMPK/FoxO signaling. In vitro testing against the adult gastro-intestinal nematode Trichuris muris suggested broad-spectrum potential for these APIs. This study suggests that APIs may serve as new leads to be further explored for drug discovery to treat parasitic nematode infections.

Spatial transcriptomics reveals antiparasitic targets associated with essential behaviors in the human parasite Brugia malayi

Lymphatic filariasis (LF) is a chronic debilitating neglected tropical disease (NTD) caused by mosquito-transmitted nematodes that afflicts over 60 million people. Control of LF relies on routine mass drug administration with antiparasitics that clear circulating larval parasites but are ineffective against adults. The development of effective adulticides is hampered by a poor understanding of the processes and tissues driving parasite survival in the host. The adult filariae head region contains essential tissues that control parasite feeding, sensory, secretory, and reproductive behaviors, which express promising molecular substrates for the development of antifilarial drugs, vaccines, and diagnostics. We have adapted spatial transcriptomic approaches to map gene expression patterns across these prioritized but historically intractable head tissues. Spatial and tissue-resolved data reveal distinct biases in the origins of known drug targets and secreted antigens. These data were used to identify potential new drug and vaccine targets, including putative hidden antigens expressed in the alimentary canal, and to spatially associate receptor subunits belonging to druggable families. Spatial transcriptomic approaches provide a powerful resource to aid gene function inference and seed antiparasitic discovery pipelines across helminths of relevance to human and animal health.

Dual RNA-Seq of Trunk Kidneys Extracted From Channel Catfish Infected With Yersinia ruckeri Reveals Novel Insights Into Host-Pathogen Interactions

Host-pathogen intectarions are complex, involving large dynamic changes in gene expression through the process of infection. These interactions are essential for understanding anti-infective immunity as well as pathogenesis. In this study, the host-pathogen interaction was analyzed using a model of acute infection where channel catfish were infected with Yersinia ruckeri. The infected fish showed signs of body surface hyperemia as well as hyperemia and swelling in the trunk kidney. Double RNA sequencing was performed on trunk kidneys extracted from infected channel catfish and transcriptome data was compared with data from uninfected trunk kidneys. Results revealed that the host-pathogen interaction was dynamically regulated and that the host-pathogen transcriptome fluctuated during infection. More specifically, these data revealed that the expression levels of immune genes involved in Cytokine-cytokine receptor interactions, the NF-kappa B signaling pathway, the JAK-STAT signaling pathway, Toll-like receptor signaling and other immune-related pathways were significantly upregulated. Y. ruckeri mainly promote pathogenesis through the flagellum gene fliC in channel catfish. The weighted gene co-expression network analysis (WGCNA) R package was used to reveal that the infection of catfish is closely related to metabolic pathways. This study contributes to the understanding of the host-pathogen interaction between channel catfish and Y. ruckeri, more specifically how catfish respond to infection through a transcriptional perspective and how this infection leads to enteric red mouth disease (ERM) in these fish.

Intracellular Bacterial Symbionts in Corals: Challenges and Future Directions

Corals are the main primary producers of coral reefs and build the three-dimensional reef structure that provides habitat to more than 25% of all marine eukaryotes. They harbor a complex consortium of microorganisms, including bacteria, archaea, fungi, viruses, and protists, which they rely on for their survival. The symbiosis between corals and bacteria is poorly studied, and their symbiotic relationships with intracellular bacteria are only just beginning to be acknowledged. In this review, we emphasize the importance of characterizing intracellular bacteria associated with corals and explore how successful approaches used to study such microorganisms in other systems could be adapted for research on corals. We propose a framework for the description, identification, and functional characterization of coral-associated intracellular bacterial symbionts. Finally, we highlight the possible value of intracellular bacteria in microbiome manipulation and mitigating coral bleaching.

Effects of Thelytokous Parthenogenesis-Inducing Wolbachia on the Fitness of Trichogramma dendrolimi Matsumura (Hymenoptera: Trichogrammatidae) in Superparasitised and Single-Parasitised Hosts

Thelytokous Wolbachia-infected Trichogramma species have long been considered as biological control agents against lepidopteran pests in agriculture and forestry. Wolbachia has been suggested to increase the probability of the superparasitism of Trichogramma, but the fate of infected offspring in the superparasitised host is still unknown. The present study aimed to evaluate the fitness of thelytokous Wolbachia-infected (TDW) and bisexual Wolbachia-free (TD) Trichogramma dendrolimi Matsumura (Hymenoptera: Trichogrammatidae) lines in superparasitised or single-parasitised hosts. The results showed that irrespective of whether Trichogramma wasps were developed from superparasitised or single-parasitised hosts, the TDW line was characterized by reduced fitness, including lower fecundity, shorter longevity, and smaller body size of F1 offspring, and lower emergence rate of F2 offspring than the TD line. This was not true for the survival rate and developmental time of F1 offspring. Additionally, the fitness parameters of T. dendrolimi that developed from superparasitised hosts were lower compared with that of T. dendrolimi that developed from single-parasitised hosts. Interestingly, Wolbachia-infected females had higher dispersal capacity than bisexual females when they developed from superparasitised hosts. The results indicated that Wolbachia negatively affects fitness of T. dendrolimi, but enhance dispersal capacity of T. dendrolimi females in superparasitism condition. Further studies need to be carried out to select the best line that will allow Wolbachia and their host Trichogramma to be better adapted to one another.

Best practices on the differential expression analysis of multi-species RNA-seq

Advances in transcriptome sequencing allow for simultaneous interrogation of differentially expressed genes from multiple species originating from a single RNA sample, termed dual or multi-species transcriptomics. Compared to single-species differential expression analysis, the design of multi-species differential expression experiments must account for the relative abundances of each organism of interest within the sample, often requiring enrichment methods and yielding differences in total read counts across samples. The analysis of multi-species transcriptomics datasets requires modifications to the alignment, quantification, and downstream analysis steps compared to the single-species analysis pipelines. We describe best practices for multi-species transcriptomics and differential gene expression.

The Wolbachia Symbiont: Here, There and Everywhere

Wolbachia symbionts, first observed in the 1920s, are now known to be present in about 30-70% of tested arthropod species, in about half of tested filarial nematodes (including the majority of human filarial nematodes), and some plant-parasitic nematodes. In arthropods, they are generally viewed as parasites while in nematodes they appear to be mutualists although this demarcation is not absolute. Their presence in arthropods generally leads to reproductive anomalies, while in nematodes, they are generally required for worm development and reproduction. In mosquitos, Wolbachia inhibit RNA viral infections, leading to populational reductions in human RNA virus pathogens, whereas in filarial nematodes, their requirement for worm fertility and survival has been channeled into their use as drug targets for filariasis control. While much more research on these ubiquitous symbionts is needed, they are viewed as playing significant roles in biological processes, ranging from arthropod speciation to human health.

Dual RNAseq analyses at soma and germline levels reveal evolutionary innovations in the elephantiasis-agent Brugia malayi, and adaptation of its Wolbachia endosymbionts

Brugia malayi is a human filarial nematode responsible for elephantiasis, a debilitating condition that is part of a broader spectrum of diseases called filariasis, including lymphatic filariasis and river blindness. Almost all filarial nematode species infecting humans live in mutualism with Wolbachia endosymbionts, present in somatic hypodermal tissues but also in the female germline which ensures their vertical transmission to the nematode progeny. These α-proteobacteria potentially provision their host with essential metabolites and protect the parasite against the vertebrate immune response. In the absence of Wolbachia wBm, B. malayi females become sterile, and the filarial nematode lifespan is greatly reduced. In order to better comprehend this symbiosis, we investigated the adaptation of wBm to the host nematode soma and germline, and we characterized these cellular environments to highlight their specificities. Dual RNAseq experiments were performed at the tissue-specific and ovarian developmental stage levels, reaching the resolution of the germline mitotic proliferation and meiotic differentiation stages. We found that most wBm genes, including putative effectors, are not differentially regulated between infected tissues. However, two wBm genes involved in stress responses are upregulated in the hypodermal chords compared to the germline, indicating that this somatic tissue represents a harsh environment to which wBm have adapted. A comparison of the B. malayi and C. elegans germline transcriptomes reveals a poor conservation of genes involved in the production of oocytes, with the filarial germline proliferative zone relying on a majority of genes absent from C. elegans. The first orthology map of the B. malayi genome presented here, together with tissue-specific expression enrichment analyses, indicate that the early steps of oogenesis are a developmental process involving genes specific to filarial nematodes, that likely result from evolutionary innovations supporting the filarial parasitic lifestyle.

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.

Novel insights into host-pathogen interactions of large yellow croakers ( Larimichthys crocea) and pathogenic bacterium Pseudomonas plecoglossicida using time-resolved dual RNA-seq of infected spleens

Host-pathogen interactions are highly complex, involving large dynamic changes in gene expression during infection. These interactions are fundamental to understanding anti-infection immunity of hosts, as well as the pathogenesis of pathogens. For bacterial pathogens interacting with animal hosts, time-resolved dual RNA-seq of infected tissue is difficult to perform due to low pathogen load in infected tissue. In this study, an acute infection model of Larimichthys crocea infected by Pseudomonas plecoglossicida was established. The spleens of infected fish exhibited typical symptoms, with a maximum bacterial load at two days post-injection (dpi). Time-resolved dual RNA-seq of infected spleens was successfully applied to study host-pathogen interactions between L. crocea and P. plecoglossicida. The spleens of infected L. crocea were subjected to dual RNA-seq, and transcriptome data were compared with those of noninfected spleens or in vitro cultured bacteria. Results showed that pathogen-host interactions were highly dynamically regulated, with corresponding fluctuations in host and pathogen transcriptomes during infection. The expression levels of many immunogenes involved in cytokine-cytokine receptor, Toll-like receptor signaling, and other immune-related pathways were significantly up-regulated during the infection period. Furthermore, metabolic processes and the use of oxygen in L. crocea were strongly affected by P. plecoglossicida infection. The WGCNA results showed that the metabolic process was strongly related to the entire immune process. For P. plecoglossicida, the expression levels of motility-related genes and flagellum assembly-related genes were significantly up-regulated. The results of this study may help to elucidate the interactions between L. crocea and P. plecoglossicida.

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.

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.

Spatial and morphological reorganization of endosymbiosis during metamorphosis accommodates adult metabolic requirements in a weevil

Bacterial intracellular symbiosis (endosymbiosis) is widespread in nature and impacts many biological processes. In holometabolous symbiotic insects, metamorphosis entails a complete and abrupt internal reorganization that creates a constraint for endosymbiont transmission from larvae to adults. To assess how endosymbiosis copes-and potentially evolves-throughout this major host-tissue reorganization, we used the association between the cereal weevil Sitophilus oryzae and the bacterium Sodalis pierantonius as a model system. S. pierantonius are contained inside specialized host cells, the bacteriocytes, that group into an organ, the bacteriome. Cereal weevils require metabolic inputs from their endosymbiont, particularly during adult cuticle synthesis, when endosymbiont load increases dramatically. By combining dual RNA-sequencing analyses and cell imaging, we show that the larval bacteriome dissociates at the onset of metamorphosis and releases bacteriocytes that undergo endosymbiosis-dependent transcriptomic changes affecting cell motility, cell adhesion, and cytoskeleton organization. Remarkably, bacteriocytes turn into spindle cells and migrate along the midgut epithelium, thereby conveying endosymbionts to midgut sites where future mesenteric caeca will develop. Concomitantly, endosymbiont genes encoding a type III secretion system and a flagellum apparatus are transiently up-regulated while endosymbionts infect putative stem cells and enter their nuclei. Infected cells then turn into new differentiated bacteriocytes and form multiple new bacteriomes in adults. These findings show that endosymbiosis reorganization in a holometabolous insect relies on a synchronized host-symbiont molecular and cellular "choreography" and illustrates an adaptive feature that promotes bacteriome multiplication to match increased metabolic requirements in emerging adults.

Phenotypic Response of Wolbachia pipientis in a Cell-Free Medium

Wolbachia, an obligate intracellular bacterium estimated to infect millions of arthropod species worldwide, is currently being utilized in novel control strategies to limit the transmission of Dengue and Zika viruses. A limitation for Wolbachia-based control approaches is the difficulty of transferring Wolbachia to novel hosts and the lack of tools for the genetic transformation of Wolbachia due to the inability to culture Wolbachia outside the insect host cell in an axenic media. Here, we applied extracellular Wolbachia to phenotypic microarrays to measure the metabolic response of Wolbachia in media formulations with different pH levels and supplementation with Casamino acids. Results suggested a pH of 6.5–6.8 and showed that the supplementation of 1 mg/mL casamino acids increased the survival and longevity of Wolbachia in an axenic medium. In addition, phenotypic microarrays are a useful tool to measure the phenotypic response of Wolbachia under different media conditions, as well as determine specific components that may be required for an axenic medium. This study is an initial step toward the development of a potential Wolbachia axenic culture system.

The shutting down of the insulin pathway: a developmental window for Wolbachia load and feminization

Using the isopod Armadillidium vulgare as a case study, we review the significance of the "bacterial dosage model", which connects the expression of the extended phenotype to the rise of the Wolbachia load. In isopods, the Insulin-like Androgenic Gland hormone (IAG) induces male differentiation: Wolbachia feminizes males through insulin resistance, presumably through defunct insulin receptors. This should prevent an autocrine development of the androgenic glands so that females differentiate instead: feminization should translate as IAG silencing and increased Wolbachia load in the same developmental window. In line with the autocrine model, uninfected males expressed IAG from the first larval stage on, long before the androgenic gland primordia begin to differentiate, and exponentially throughout development. In contrast in infected males, expression fully stopped at stage 4 (juvenile), when male differentiation begins. This co-occurred with the only significant rise in the Wolbachia load throughout the life-stages. Concurrently, the raw expression of the bacterial Secretion Systems co-increased, but they were not over-expressed relative to the number of bacteria. The isopod model leads to formulate the "bacterial dosage model" throughout extended phenotypes as the conjunction between bacterial load as the mode of action, timing of multiplication (pre/post-zygotic), and site of action (soma vs. germen).

Prediction pipeline for discovery of regulatory motifs associated with Brugia malayi molting

Filarial nematodes can cause debilitating diseases in humans. They have complicated life cycles involving an insect vector and mammalian hosts, and they go through a number of developmental molts. While whole genome sequences of parasitic worms are now available, very little is known about transcription factor (TF) binding sites and their cognate transcription factors that play a role in regulating development. To address this gap, we developed a novel motif prediction pipeline, Emotif Alpha, that integrates ten different motif discovery algorithms, multiple statistical tests, and a comparative analysis of conserved elements between the filarial worms Brugia malayi and Onchocerca volvulus, and the free-living nematode Caenorhabditis elegans. We identified stage-specific TF binding motifs in B. malayi, with a particular focus on those potentially involved in the L3-L4 molt, a stage important for the establishment of infection in the mammalian host. Using an in vitro molting system, we tested and validated three of these motifs demonstrating the accuracy of the motif prediction pipeline.

Diseases caused by parasitic worms such as the filariae are among the leading causes of morbidity in the developing world. Very little is known about how development is regulated in these vector-transmitted parasites. We have developed a computational method to identify motifs that correspond to transcription factor binding sites in the genome of the parasitic worm, Brugia malayi, one of the causative agents of lymphatic filariasis. Using this approach, we were able to predict stage-specific transcription factor binding sites involved in a stage of the molting process important for the establishment of the infection. We validated the role of these motifs using an in vitro molting system.

Next-generation technologies for studying host-pathogen interactions: a focus on dual transcriptomics, CRISPR/Cas9 screening and organs-on-chips

Pathogens constantly interact with their hosts and the environment, and therefore have evolved unique virulence mechanisms to target and breach host defense barriers and manipulate host immune response to establish an infection. Advances in technologies that allow genome mining, gene editing such as CRISPR/Cas9, genomic, epigenomic and transcriptomic studies such as dual RNA-seq, coupled with bioinformatics, have accelerated the field of host-pathogen interactions within a broad range of infection models. Underpinning of the molecular changes that accompany invasion of eukaryotic cells with pathogenic microorganisms at the intersection of host, pathogen and their local environment has provided a better understanding of infectious disease mechanisms and antimicrobial strategies. The recent evolution of physiologically relevant three-dimensional (3-D) tissue/organ models and microfluidic organ-on-chip devices also provided a window to a more predictive framework of infectious disease processes. These approaches combined hold the potential to highly impact discovery of novel drug targets and vaccine candidates of the future. Here, we review three of the available and emerging technologies-dual RNA-seq, CRISPR/Cas9 screening and organs-on-chips, applicable to the high throughput study and deciphering of interaction networks between pathogens and their hosts that are critical for the development of novel therapeutics.

Sex chromosome evolution in parasitic nematodes of humans

Sex determination mechanisms often differ even between related species yet the evolution of sex chromosomes remains poorly understood in all but a few model organisms. Some nematodes such as Caenorhabditis elegans have an XO sex determination system while others, such as the filarial parasite Brugia malayi, have an XY mechanism. We present a complete B. malayi genome assembly and define Nigon elements shared with C. elegans, which we then map to the genomes of other filarial species and more distantly related nematodes. We find a remarkable plasticity in sex chromosome evolution with several distinct cases of neo-X and neo-Y formation, X-added regions, and conversion of autosomes to sex chromosomes from which we propose a model of chromosome evolution across different nematode clades. The phylum Nematoda offers a new and innovative system for gaining a deeper understanding of sex chromosome evolution.

Drug Repurposing of Bromodomain Inhibitors as Potential Novel Therapeutic Leads for Lymphatic Filariasis Guided by Multispecies Transcriptomics

The current treatment regimen for lymphatic filariasis is mostly microfilaricidal. In an effort to identify new drug candidates for lymphatic filariasis, we conducted a three-way transcriptomics/systems biology study of one of the causative agents of lymphatic filariasis, Brugia malayi, its Wolbachia endosymbiont wBm, and its vector host Aedes aegypti at 16 distinct B. malayi life stages. B. malayi upregulates the expression of bromodomain-containing proteins in the adult female, embryo, and microfilaria stages. In vitro, we find that the existing cancer therapeutic JQ1(+), which is a bromodomain and extraterminal protein inhibitor, has adulticidal activity in B. malayi.

To better understand the transcriptomic interplay of organisms associated with lymphatic filariasis, we conducted multispecies transcriptome sequencing (RNA-Seq) on the filarial nematode Brugia malayi, its Wolbachia endosymbiont wBm, and its laboratory vector Aedes aegypti across the entire B. malayi life cycle. In wBm, transcription of the noncoding 6S RNA suggests that it may be a regulator of bacterial cell growth, as its transcript levels correlate with bacterial replication rates. For A. aegypti, the transcriptional response reflects the stress that B. malayi infection exerts on the mosquito with indicators of increased energy demand. In B. malayi, expression modules associated with adult female samples consistently contained an overrepresentation of genes involved in chromatin remodeling, such as the bromodomain-containing proteins. All bromodomain-containing proteins encoded by B. malayi were observed to be upregulated in the adult female, embryo, and microfilaria life stages, including 2 members of the bromodomain and extraterminal (BET) protein family. The BET inhibitor JQ1(+), originally developed as a cancer therapeutic, caused lethality of adult worms in vitro, suggesting that it may be a potential therapeutic that can be repurposed for treating lymphatic filariasis.

IMPORTANCE The current treatment regimen for lymphatic filariasis is mostly microfilaricidal. In an effort to identify new drug candidates for lymphatic filariasis, we conducted a three-way transcriptomics/systems biology study of one of the causative agents of lymphatic filariasis, Brugia malayi, its Wolbachia endosymbiont wBm, and its vector host Aedes aegypti at 16 distinct B. malayi life stages. B. malayi upregulates the expression of bromodomain-containing proteins in the adult female, embryo, and microfilaria stages. In vitro, we find that the existing cancer therapeutic JQ1(+), which is a bromodomain and extraterminal protein inhibitor, has adulticidal activity in B. malayi.

Pyruvate produced by Brugia spp. via glycolysis is essential for maintaining the mutualistic association between the parasite and its endosymbiont, Wolbachia

Human parasitic nematodes are the causative agents of lymphatic filariasis (elephantiasis) and onchocerciasis (river blindness), diseases that are endemic to more than 80 countries and that consistently rank in the top ten for the highest number of years lived with disability. These filarial nematodes have evolved an obligate mutualistic association with an intracellular bacterium, Wolbachia, a symbiont that is essential for the successful development, reproduction, and survival of adult filarial worms. Elimination of the bacteria causes adult worms to die, making Wolbachia a primary target for developing new interventional tools to combat filariases. To further explore Wolbachia as a promising indirect macrofilaricidal drug target, the essential cellular processes that define the symbiotic Wolbachia-host interactions need to be identified. Genomic analyses revealed that while filarial nematodes encode all the enzymes necessary for glycolysis, Wolbachia does not encode the genes for three glycolytic enzymes: hexokinase, 6-phosphofructokinase, and pyruvate kinase. These enzymes are necessary for converting glucose into pyruvate. Wolbachia, however, has the full complement of genes required for gluconeogenesis starting with pyruvate, and for energy metabolism via the tricarboxylic acid cycle. Therefore, we hypothesized that Wolbachia might depend on host glycolysis to maintain a mutualistic association with their parasitic host. We did conditional experiments in vitro that confirmed that glycolysis and its end-product, pyruvate, sustain this symbiotic relationship. Analysis of alternative sources of pyruvate within the worm indicated that the filarial lactate dehydrogenase could also regulate the local intracellular concentration of pyruvate in proximity to Wolbachia and thus help control bacterial growth via molecular interactions with the bacteria. Lastly, we have shown that the parasite's pyruvate kinase, the enzyme that performs the last step in glycolysis, could be a potential novel anti-filarial drug target. Establishing that glycolysis is an essential component of symbiosis in filarial worms could have a broader impact on research focused on other intracellular bacteria-host interactions where the role of glycolysis in supporting intracellular survival of bacteria has been reported.

The effect of Wolbachia on gene expression in Drosophila paulistorum and its implications for symbiont-induced host speciation

Background

The Neotropical fruit fly Drosophila paulistorum (Diptera: Drosophilidae) is a species complex in statu nascendi comprising six reproductively isolated semispecies, each harboring mutualistic Wolbachia strains. Although wild type flies of each semispecies are isolated from the others by both pre- and postmating incompatibilities, mating between semispecies and successful offspring development can be achieved once flies are treated with antibiotics to reduce Wolbachia titer. Here we use RNA-seq to study the impact of Wolbachia on D. paulistorum and investigate the hypothesis that the symbiont may play a role in host speciation. For that goal, we analyze samples of heads and abdomens of both sexes of the Amazonian, Centro American and Orinocan semispecies of D. paulistorum.

Results

We identify between 175 and 1192 differentially expressed genes associated with a variety of biological processes that respond either globally or according to tissue, sex or condition in the three semispecies. Some of the functions associated with differentially expressed genes are known to be affected by Wolbachia in other species, such as metabolism and immunity, whereas others represent putative novel phenotypes involving muscular functions, pheromone signaling, and visual perception.

Conclusions

Our results show that Wolbachia affect a large number of biological functions in D. paulistorum, particularly when present in high titer. We suggest that the significant metabolic impact of the infection on the host may cause several of the other putative and observed phenotypes. We also speculate that the observed differential expression of genes associated with chemical communication and reproduction may be associated with the emergence of pre- and postmating barriers between semispecies, which supports a role for Wolbachia in the speciation of D. paulistorum.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5816-9) contains supplementary material, which is available to authorized users.

Dirofilaria and Wolbachia in mosquitoes (Diptera: Culicidae) in central European Russia and on the Black Sea coast

Dirofilariasis is endemic in Russia, as well as in many other European countries. The aim of this study was to assess the ability of mosquitoes to transfer Dirofilaria immitis and Dirofilaria repens in regions with temperate and subtropical climates. The possible impact of the symbiotic bacterium Wolbachia on Dirofilaria transmission was also investigated. 5333 female mosquitoes were collected at 11 points in central European Russia and on the Black Sea coast during the period 2013–2017. Out of 20 mosquito species examined, 14 were infected with D. repens and 13 with D. immitis. Both species of Dirofilaria were found in different climatic regions. The total Dirofilaria spp. estimated infection rate (EIR) in the central part of Russia varied from 3.1% to 3.7% and, in the southern region, from 1.1% to 3.0%. The highest estimated infection rate was found in Anopheles messeae, the lowest in Culex pipiens. The greatest epidemiological danger was represented by Aedes aegypti, Ae. geniculatus, An. messeae and Ae. communis. Six out of 20 mosquito species were infected with Wolbachia. Pools of Aedes albopictus, Cx. pipiens and Coquillettidia richiardii were simultaneously infected with Dirofilaria and Wolbachia. After checking mosquitoes individually, it was found that there was no development of Dirofilaria to the infective larval stage in specimens infected with Wolbachia. Twenty-two Dirofilaria-infective pools were Wolbachia-free and only two mosquito pools were Wolbachia-infected. The potential for transmission of Dirofilaria in mosquito species naturally uninfected with the symbiotic bacterium Wolbachia is higher than in species infected with the bacterium.

Targeted enrichment outperforms other enrichment techniques and enables more multi-species RNA-Seq analyses

Enrichment methodologies enable the analysis of minor members in multi-species transcriptomic data. We compared the standard enrichment of bacterial and eukaryotic mRNA to a targeted enrichment using an Agilent SureSelect (AgSS) capture for Brugia malayi, Aspergillus fumigatus, and the Wolbachia endosymbiont of B. malayi (wBm). Without introducing significant systematic bias, the AgSS quantitatively enriched samples, resulting in more reads mapping to the target organism. The AgSS-enriched libraries consistently had a positive linear correlation with their unenriched counterparts (r² = 0.559–0.867). Up to a 2,242-fold enrichment of RNA from the target organism was obtained following a power law (r² = 0.90), with the greatest fold enrichment achieved in samples with the largest ratio difference between the major and minor members. While using a single total library for prokaryote and eukaryote enrichment from a single RNA sample could be beneficial for samples where RNA is limiting, we observed a decrease in reads mapping to protein coding genes and an increase in multi-mapping reads to rRNAs in AgSS enrichments from eukaryotic total RNA libraries compared to eukaryotic poly(A)-enriched libraries. Our results support a recommendation of using AgSS targeted enrichment on poly(A)-enriched libraries for eukaryotic captures, and total RNA libraries for prokaryotic captures, to increase the robustness of multi-species transcriptomic studies.

Periodic and coordinated gene expression between a diazotroph and its diatom host

In the surface ocean, light fuels photosynthetic carbon fixation of phytoplankton, playing a critical role in ecosystem processes including carbon export to the deep sea. In oligotrophic oceans, diatom-diazotroph associations (DDAs) play a keystone role in ecosystem function because diazotrophs can provide otherwise scarce biologically available nitrogen to the diatom host, fueling growth and subsequent carbon sequestration. Despite their importance, relatively little is known about the nature of these associations in situ. Here we used metatranscriptomic sequencing of surface samples from the North Pacific Subtropical Gyre (NPSG) to reconstruct patterns of gene expression for the diazotrophic symbiont Richelia and we examined how these patterns were integrated with those of the diatom host over day-night transitions. Richelia exhibited significant diel signals for genes related to photosynthesis, N₂ fixation, and resource acquisition, among other processes. N₂ fixation genes were significantly co-expressed with host nitrogen uptake and metabolism, as well as potential genes involved in carbon transport, which may underpin the exchange of nitrogen and carbon within this association. Patterns of expression suggested cell division was integrated between the host and symbiont across the diel cycle. Collectively these data suggest that symbiont-host physiological ecology is strongly interconnected in the NPSG.

Large-Scale Identification of Wolbachia pipientis Effectors

Wolbachia pipientis is an intracellular symbiont of arthropods well known for the reproductive manipulations induced in the host and, more recently, for the ability of Wolbachia to block virus replication in insect vectors. Since Wolbachia cannot yet be genetically manipulated, and due to the constraints imposed when working with an intracellular symbiont, little is known about mechanisms used by Wolbachia for host interaction. Here we employed a bioinformatics pipeline and identified 163 candidate effectors, potentially secreted by Wolbachia into the host cell. A total of 84 of these candidates were then subjected to a screen of growth defects induced in yeast upon heterologous expression which identified 14 top candidates likely secreted by Wolbachia. These predicted secreted effectors may function in concert as we find that their native expression is correlated and is highly upregulated at specific time points during Drosophila development. In addition, the evolutionary histories of some of these predicted effectors are also correlated, suggesting they may function together, or in the same pathway, during host infection. Similarly, most of these predicted effectors are limited to one or two Wolbachia strains—perhaps reflecting shared evolutionary history and strain specific functions in host manipulation. Identification of these Wolbachia candidate effectors is the first step in dissecting the mechanisms of symbiont–host interaction in this important system.

Mining Filarial Genomes for Diagnostic and Therapeutic Targets

Filarial infections of humans cause some of the most important neglected tropical diseases. The global efforts for eliminating filarial infections by mass drug administration programs may require additional tools (safe macrofilaricidal drugs, vaccines, and diagnostic biomarkers). The accurate and sensitive detection of viable parasites is essential for diagnosis and for surveillance programs. Current community-wide treatment modalities do not kill the adult filarial worms effectively; hence, there is a need to identify and develop safe macrofilaricidal drugs. High-throughput sequencing, mass spectroscopy methods and advances in computational biology have greatly accelerated the discovery process. Here, we describe post-genomic developments toward the identification of diagnostic biomarkers and drug targets for the filarial infection of humans.

Two's company: studying interspecies relationships with dual RNA-seq

Organisms do not exist isolated from each other, but constantly interact. Cells can sense the presence of interaction partners by a range of receptors and, via complex regulatory networks, specifically react by changing the expression of many of their genes. Technological advances in next-generation sequencing over the recent years now allow us to apply RNA sequencing to two species at the same time (dual RNA-seq), and thus to directly study the gene expression of two interacting species without the need to physically separate cells or RNA. In this review, we give an overview over the latest studies in interspecies interactions made possible by dual RNA-seq, ranging from pathogenic to symbiotic relationships. We summarize state-of-the-art experimental techniques, bioinformatic data analysis and data interpretation, while also highlighting potential problems and pitfalls starting from the selection of meaningful time points and number of reads to matters of rRNA depletion. A short outlook on new trends in the field of dual RNA-seq concludes this review, looking at sequencing of non-coding RNAs during host-pathogen interactions and the prediction of molecular interspecies interactions networks.