- Open Access
Localization of Wolbachia-like gene transcripts and peptides in adult Onchocerca flexuosa worms indicates tissue specific expression
Parasites & Vectors volume 6, Article number: 2 (2013)
Most filarial species in the genus Onchocerca depend on Wolbachia endobacteria to successfully carry out their life cycle. O. flexuosa is a Wolbachia-free species, but its genome contains Wolbachia-like sequences presumably obtained from Wolbachia via horizontal gene transfer. Proteogenomic studies have shown that many of these Wolbachia-like sequences are expressed in adult worms.
Six Wolbachia-like sequences in O. flexuosa were chosen for further study based on their sequence conservation with Wolbachia genes, length of predicted open reading frames, and expression at the RNA and/or protein levels. In situ hybridization and immunohistochemical labeling were used to localize Wolbachia-like transcripts and peptides in adult worm tissues.
RNA probes representing three of the six target sequences produced hybridization signals in worm tissues. These probes bound to transcripts in the intestine and lateral chords of both sexes, in the hypodermis, median chords and uteri in females, and in sperm precursor cells in males. Antibodies raised to three peptides corresponding to these transcripts bound to specific bands in a soluble extract of adult O. flexuosa by Western blot that were not labeled by control antibodies in pre-immune serum. Two of the three antibodies produced labeling patterns in adult worm sections that were similar to those of the RNA probes, while the third produced a different pattern.
A subset of the Wolbachia-like sequences present in the genome of the Wolbachia-free filarial species O. flexuosa are transcribed in tissues where Wolbachia reside in infected filarial species. Some of the peptides and/or proteins derived from these transcripts appear to be concentrated in the same tissues while others may be exported to other regions of the worm. These results suggest that horizontally transferred Wolbachia genes and gene products may replicate important Wolbachia functions in uninfected filarial worms.
Filarial nematodes comprise a superfamily of parasitic worms that infect a wide array of hosts, including humans. Three filarial species (Onchocerca volvulus, Wuchereria bancrofti, and Brugia malayi) are important human pathogens in the tropics, and many other filarial species infect wild or domestic animals throughout the world. One feature shared by many filarial pathogens is their association with Wolbachia endobacteria. Prior studies have shown that Wolbachia play important roles in worm growth, molting, reproduction, and pathogenesis [1–4].
Filarial nematodes have been divided into eight subfamilies based on classical parasitological criteria. Two of these subfamilies, the Onchocercinae and Dirofilariinae, appear to be dominated by Wolbachia-dependent species [5–7]. The abundance of Wolbachia-dependent species relative to Wolbachia-free species and the congruencies between the phylogenetic lineages of Wolbachia and their filarial hosts suggest that Wolbachia entered the filarial lineage prior to the differentiation of the Onchocercinae and Dirofilariinae [5–7]. Recent studies have shown that two Wolbachia-free filarial species, Onchocerca flexuosa (subfamily: Onchocercinae) and Acanthocheilonema viteae (subfamily: Dirofilariinae), contain Wolbachia-like DNA sequences in their nuclear genomes, indicating that these species may have been associated with Wolbachia in the ancient past . Since O. flexuosa and A. viteae are relatively distantly related [6, 9], we expect that this trend will prove consistent among other Wolbachia-free onchocercids and dirofilariids.
Our previous surveys of the genome and transcriptome of O. flexuosa identified sequence fragments with similarity to 178 different Wolbachia genes [8, 10]. qRT-PCR reactions and partial transcriptome sequencing indicated that many of these Wolbachia-like sequences are transcribed, and a mass spectrometry study with follow-up immunohistology and Western blot studies indicated that at least two Wolbachia-like sequences were translated into Wolbachia-like peptides . Wolbachia are known to be proficient at transferring genetic material to host cells and many Wolbachia-like sequences are present in the genomes of Wolbachia-dependent filarial worms; however, the Wolbachia-like sequences found in the genomes of dependent species like B. malayi are thought to be degenerate and non-functional . Retention and expression of Wolbachia-like sequences in O. flexuosa suggests that they may have essential roles in the biology of Wolbachia-free filarial worms.
Prior work has shown that Wolbachia endobacteria are restricted to specific tissues in filarial nematodes [7, 12–16]. In early development, the vertically transmitted bacteria that were present in the fertilized oocyte segregate to cells that give rise to the hypodermal lateral chords . Wolbachia from the lateral chords then invade the ovaries and testis prior to sexual maturation [12, 17]. This pattern of localization (i.e., lateral chords and reproductive tissues) may be critical to Wolbachia’s function as a mutualist and to its role in development and reproduction. Therefore, localization studies were performed to gain insight into the possible function(s) of transferred and retained Wolbachia-like sequences in the Wolbachia- free filarial species O. flexuosa.
The present study used in situ hybridization to localize expression of Wolbachia-like transcripts in adult O. flexuosa and immunohistochemical staining to localize peptides corresponding to these sequences. Thus far, all of the RNA probes that produce a signal in adult tissues stain the intestine and lateral chords in both sexes, the hypodermis, median chords, and uteri of females, and developing sperm in males. Two of the three Wolbachia-like peptides were identified in the same tissues and in developing embryos within the females. The third Wolbachia-like peptide was found in similar tissues and in the basal areas of somatic muscles.
O. flexuosa nodules were cut from the skins of freshly-shot European red deer (Cervus elaphus) following community hunts in northern Germany (Schleswig-Holstein) as previously described . Adult B. malayi were obtained from the Filariasis Research Reagent Resource Center . Worms to be used for RNA or DNA isolation were stored in TRIzol (Invitrogen, Carlsbad, CA, USA) or 1x phosphate buffered saline, respectively, at -80°C. O. flexuosa nodules intended for histological examination were fixed in 4% buffered formalin, embedded in paraffin and sectioned at 5μm thickness according to standard histological procedure.
Nucleic acid isolation
DNA was isolated using an E.Z.N.A Tissue DNA Kit (Omega Bio-Tek, Norcross, GA, USA) according to the manufacturer’s suggested protocol. Total RNA was isolated by homogenization in TRIzol (Invitrogen) and organic extraction with 1-bromo-3-chloropropane and purified using an RNeasy Mini Kit with the optional on-column DNase digest (Qiagen, Valencia, CA, USA). A second DNase treatment was performed using the Turbo DNA-free kit (Invitrogen). cDNA was synthesized using qScript cDNA SuperMix (Quanta, Gaithersburg, MD, USA) and tested for DNA contamination by PCR with intron-spanning primer sets as previously described .
The O. flexuosa adult transcriptome was sequenced using Roche/454 Titanium technology and assembled using the Newber v2.5 assembler as described previously . Sequences are available from the Genbank transcriptome shotgun assembly database (BioProject number 62565, accession numbers JI459010-JI484230) and from Nematode.net . Wolbachia-like sequences were identified from O. flexuosa transcripts by BLASTX search against the non-redundant protein database . Measurements of sequence conservation were based on the e-values of hits obtained in this blast search. The relative expression of each O. flexuosa isogroup (unique sequence locus) was estimated based on the number of reads included in each assembled isogroup normalized to the total length of the isogroup (the summed length of all contigs included in the isogroup) (Additional file 1: Table S1). Putative Wolbachia-like peptides were identified by blasting six-frame translations of O. flexuosa transcripts against the non-redundant protein database (>50% sequence identity shared with Wolbachia protein, bit score >35).
RNA probe construction and in situ hybridization
Probes were constructed as previously described . Briefly, 300-600 bp target sequences were PCR amplified from cDNA (see Additional file 2: Table S2 for primers) and cloned using the TOPO TA Cloning Kit Dual Promoter (Invitrogen). Plasmids were purified with the QIAprep Spin Miniprep Kit (Qiagen) and linearized by restriction digest with either EcoRV or BamHI. Following purification of linearized plasmid (QIAquick PCR Purification Kit, Qiagen), biotinlyated antisense RNA probes and sense controls were transcribed from the plasmid template using the MEGAscript SP6 and T7 in vitro transcription kits (Invitrogen) with biotinylated NTPs (Roche Diagnostics, Indianapolis, IN, USA). The final biotinlyated RNA probes were purified and concentrated by ethanol precipitation and stored in 1 × tris-EDTA buffer, pH 8.0.
Slides from various nodules were tested with control probes to determine the quality of RNA preservation. This step is necessary due to the short lifespan of O. flexuosa, as studies have indicated that 20% of nodules collected (even from very young deer) contain dead worms . In situ hybridizations showed inconsistent results in older, mf producing worms. This may have been due to reduced penetration of the fixative through larger, tougher nodules or due to reduced transcription rates in the older worms. Slides prepared from small, soft nodules containing young adult worms showed better RNA integrity. Therefore, we chose to focus on young worms in our in situ hybridization studies.
In situ hybridizations using formalin-fixed, paraffin-embedded specimens were carried out as previously described . Briefly, slides were hybridized with 1μg/mL of RNA probe in hybridization buffer (50% formamide, 5x SSC, 0.3mg/mL yeast tRNA, 100μg/mL heparin, 1x Denhardt’s reagent, 0.1% CHAPS and 5mM EDTA) overnight at 60°C. Stringency washes (60°C for 30m) were carried out using reagents from the “In situ Hybridization Detection System” (Dako, Carpinteria, CA, USA). The same kit was used for colorimetric detection according to the manufacturer’s suggested protocol. Slides were viewed using an Olympus-BX40 microscope (Olympus, Tokyo, Japan) and photographed with an Infinity2 digital microscope camera using Infinity Capture software (Lumenera, Ottowa, Ontario, Canada). For fluorescent detection, washed sections were incubated with 5μg/mL streptavidin Alexa Fluor 488 (Invitrogen) in 1x phosphate buffered saline with 0.1% bovine serum albumin in the dark, at room temperature. After 20 minutes, 5μg/mL wheat germ aggutinin Alexa Fluor 633 was applied to the slide to highlight cell membranes, and the incubation was allowed to proceed for 10 more minutes. Finally, sections were rinsed in 1x tris buffered saline and cover slips were mounted with ProLong Gold Antifade Reagent with DAPI (Invitrogen). Fluorescent labeling was viewed with a Zeiss Axioskop 2 Mot Plus fluorescence microscope equipped with an Axiocam MRm monochrome camera, and images were captured using Axiovision 4.6 software (Carl Zeiss Inc., Thornwood, NY, USA).
Anti-peptide antibodies were produced and purified by LifeTein LLC (South Plainfield, NJ, USA). Targeted portions of Wolbachia-like peptides were selected based on predicted chemical and immunogenic properties (Additional file 3: Table S3). Target peptides were synthesized and coupled to a keyhole limpit hemocyanin (KLH) carrier protein. Rabbits were bled to collect pre-immune serum and then immunized with the peptide/KLH conjugates. Antibody production was allowed to proceed for 12-15 weeks with 3 booster immunizations. Following the terminal bleed, polyclonal antibodies were affinity purified from serum using the target peptide (without the KLH carrier) and tested by ELISA prior to use. Polyclonal antibodies against the KLH carrier protein were produced and purified in the same manner. Total IgG was purified from rabbit pre-immune sera using the Protein A Agarose Kit (KPL, Gaithersburg, MD, USA).
Binding of antibodies to Wolbachia-like peptides present in O. flexuosa antigen extract by Western blot
Western blots were performed as previously described . Briefly, O. flexuosa total worm homogenate separated by SDS-PAGE gel electrophoresis and transferred to a nitrocellulose membrane. Blot strips were probed with peptide antibodies (2.5μg/mL for antibodies against the peptides from isotig12596 and isotig21532, 5μg/mL for antibodies against the HlyD peptide) and purified IgG from the corresponding rabbits’ pre-immune sera (5μg/mL from rabbits used to produce antibodies against the peptides from isotig12596 and isotig21532, 10μg/mL from the rabbit used to produce antibodies against the HlyD peptide) in 1x phosphate buffered saline with 0.5% Tween (PBS/T) overnight at 4°C and washed with PBS/T at room temperature. Anti-rabbit IgG(Fc) AP conjugate (Promega, Sunnyvale, CA, USA) was diluted 1:3,500 in PBS/T and applied to the strips for 1h at 37°C. Strips were again washed with PBS/T and developed using NBT/BCIP substrate (Promega).
Various dilutions of primary antibodies were tested in order to optimize signal/background. Antibodies against peptides from isotig12596, isotig21532, HlyD and KLH (negative control) were used at 9.7μg/mL, 5.7μg/mL, 7.4μg/mL, and 1.4 μg/mL, respectively. Visualization was mostly performed with alkaline phosphatase anti-alkaline phosphatase reagents (Dako) as previously described [12, 22], and slides were viewed using an Olympus-BX40 microscope (Olympus) and photographed with an Infinity2 digital microscope camera using Infinity Capture software (Lumenera). FITC labeled goat anti-mouse IgG (1:300; Sigma, St. Louis, MO, USA) was used as secondary antibodies for immunofluorescent labeling, and DAPI and wheat germ agglutinin Alexa Fluor 633 conjugate were used to visualize DNA and cell membranes, respectively (Invitrogen). Fluorescent labeling was viewed with a Zeiss Axioskop 2 Mot Plus fluorescence microscope equipped with an Axiocam MRm monochrome camera, and images were captured using Axiovision 4.6 software (Carl Zeiss Inc.).
Results and discussion
Genomic and transcriptomic surveys have shown that at least 178 different Wolbachia genes are represented by sequence fragments in the genome of O. flexuosa[8, 10]. In our previous mass specrometry analysis of O. flexuosa adult worm lysate, three unique peptides mapped to two Wolbachia proteins included in our comparative database . One peptide mapped to a lipoprotein releasing system trans-membrane protein (LolC) present in several Wolbachia strains. The localization of this peptide and its corresponding RNA were previously described . The other two peptides mapped to an HlyD family secretion protein from the Wolbachia endosymbiont of Culex quinquefasciatus (Figure 1). Like the LolC protein, HlyD is a transmembrane component of a bacterial ABC transport system . Of the two peptides that mapped to HlyD, we chose to focus on the one predicted to have the most favorable chemical and immunogenic qualities for further analysis despite the fact that no corresponding sequence has been identified in among our (as yet incomplete) set of O. flexuosa genomic and transcriptomic sequences.
None of the Wolbachia-like peptides predicted from the O. flexuosa transcriptome were identified in our previous proteomic study . While this could certainly be due to lack of expression, it is quite likely that the Wolbachia-like peptides are difficult to detect in total worm lysate due to their low abundance compared to other O. flexuosa proteins. Indeed, Wolbachia proteins have proven difficult to detect using this method even in Wolbachia-dependent species . Therefore, a targeted approach employing specific anti-peptide antibodies was used to search for evidence of protein-level expression among our Wolbachia-like transcript sequences. We decided on three potential indicators of importance to determine which of the Wolbachia-like sequences warranted further examination: sequence conservation with Wolbachia, putative size of the open reading frame, and expression level (Table 1). Purifying selection should protect sequences that are vital for the parasite. Therefore, we focused on sequences that were most similar to the presumed donor sequence from Wolbachia (or the closest proxy for the donor sequence included in public databases) and sequences with long open reading frames (i.e., without numerous stop codons and frameshift mutations). Though no exact criteria for biologically relevant expression levels are defined, it is generally assumed that abundantly expressed genes are important. Some of the isogroups containing Wolbachia-like sequences are expressed at relatively low (presumably background) levels, while others are expressed to a much higher degree (Additional file 1: Table S1). The top ten most abundantly expressed Wolbachia-like sequences all had expression levels higher than the average and median expression seen among all O. flexuosa isogroups. Ultimately, five Wolbachia-like sequences fitting our criteria were selected for analysis (Table 2). These include sequences similar to GMP synthase, aminopeptidase P, succinyl-diaminopimelate desuccinylase, citrate synthase, and a Wolbachia hypothetical protein. According to the COG and KEGG databases [25–28], GMP synthase is involved in purine biosynthesis, aminopeptidase P is involved in amino acid transport and metabolism, succinly-diaminopimelate desuccinylase is involved in lysine biosynthesis, and citrate synthase is involved in the Krebs cycle. The function of the Wolbachia hypothetical protein is, of course, unknown.
Foster et al.  reported that the Wolbachia endosymbionts of B. malayi are capable of de novo nucleotide synthesis, and they suggested that the bacteria may play a role in metabolic provisioning. A recent analysis of the genome of the Wolbachia endosymbiont of Onchocerca ochengi further supported this hypothesis, suggesting that Wolbachia may supplement the host nucleotide pool . It is interesting that enzymes related to energy metabolism and purine synthesis were identified using our selection algorithm, which did not include potential enzymatic function as a selection criterion. Future work will be required to determine whether these sequences encode enzymes that are functional in O. flexuosa.
In situ hybridization
RNA probes corresponding to the five sequences of interest selected from the O. flexuosa transcriptome and Wolbachia HlyD were constructed and used in in situ hybridization studies. Because no genomic or transcriptomic sequence related to Wolbachia HlyD has been found in O. flexuosa, the HlyD probe was amplified from the Wolbachia endosymbiont of B. malayi. Probes against isotig12596 and Wolbachia HlyD produced strong signals in O. flexuosa tissues (Figure 2 and Additional file 4: Figure S1, respectively). The probe against O. flexuosa isotig21532 produced a much weaker signal (Additional file 5: Figure S2), and the probes against isotig14485, isotig12363 and contig10287 produced no signal at all. This was surprising, because expression levels for these transcripts were not estimated to be significantly lower than that of isotig12596 (Additional file 1: Table S1). It could be that these three sequences are expressed under conditions that are not represented in the particular worm specimens we examined or that some Wolbachia-like sequences are expressed at a low level throughout the body rather than concentrating in specific tissues, thus making them difficult to detect using this method. It is also possible that technical issues with the RNA probes led to this result.
To date, all of the RNA probes against Wolbachia-like sequences that have successfully labeled O. flexuosa tissues produced very similar labeling patterns (see [8, 10], Figure 2, Additional file 4: Figure S1, and Additional file 5: Figure S2). The results of in situ hybridizations with the probe against isotig12596 (putative aminopeptidase P) are shown as an example (Figure 2). Intense labeling was seen in the lateral chords and intestine of both sexes (Figure 2B, E, and F), in the hypodermis, in the empty uteri of young females (Figure 2B), and in the testis and developing sperm of males (Figure 2D-F). Mature sperm in the vas deferens were not labeled, and the signal in the lateral chords decreased towards the tail end of the male worm where the mature sperm reside (Figure 2G).
The intestine, lateral chords and the reproductive organs are among the most metabolically active tissues in the adult worm. Since these tissues could be responsible for the production of many gene products, we tested probes unrelated to the Wolbachia-like sequences to ensure that the detected pattern is specific to this subset of sequences. As expected, unrelated probes produced dissimilar labeling patterns, as was shown for a probe targeting transcripts derived from the major sperm protein gene (Figure 3) .
A shared expression pattern can be taken as evidence that a group of genes is under the control of either the same or functionally similar promoters or other regulatory elements. Approximately 15% of the genes encoded by B. malayi are organized into operons [31, 32]. Since the complete genome of O. flexuosa has not been sequenced, we cannot rule out the possibility that some of the Wolbachia-like sequences are organized in this way, giving rise to a characteristic pattern of expression. However, clusters of Wolbachia-like sequences were not identified in genomic surveys of O. flexuosa, so it is unlikely that they are all organized into operons.
Wolbachia are commonly detected in the lateral chords, embryos, developing (but not mature) sperm, and sometimes in the intestine of infected filarial nematodes [7, 12]. Interestingly, the expression pattern of Wolbachia-like sequences in O. flexuosa mimics the distribution of Wolbachia in infected species. This may indicate that some of the same tissues (e.g., lateral chords, reproductive organs, intestine etc.) that harbor Wolbachia in infected species produce Wolbachia-related products in Wolbachia-independent species.
Western blot results
Polyclonal antibodies were raised against a peptide from HlyD (Figure 1) and translated peptides from isotig12596, and isotig21532 (see Table 2). These were used to probe O. flexuosa lysate by Western blot (Figure 4). Antibodies to the predicted peptide from isotig12596 specifically bound to a protein band at approximately 50 kDa, while antibodies against the predicted peptide from isotigs21532 and the Wolbachia HlyD peptide each bound to two protein bands (molecular weights of approximately 120 and 260 kDa for the antibodies to the predicted peptide from isotig21532, and approximately 48 and 18 kDa for the antibodies to the HlyD peptide). Though the detection of multiple bands could be a sign of cross reactivity, it could also reflect multiple splice isoforms, post-translational modification, proteolytic processing, or degradation. IgG from the same rabbits’ pre-immune sera did not bind to bands recognized by their immune sera.
It is not possible to estimate the masses of the O. flexuosa proteins based on our transcriptome data due to heavy fragmentation. The masses of the top blastx matches to isotigs12596 (aminopeptidase P from the Wolbachia endosymbiont of Muscidifurax uniraptor) and 21532 (a hypothetical protein from the Wolbachia endosymbiont of Onchocerca volvulus) and the HlyD protein are estimated to be approximately 58 kDa, 83 kDa, and 56 kDa, respectively, based on their reported amino acid sequences. The band detected by the antibodies to the predicted peptide from isotig12596 was similar in size to the corresponding Wolbachia protein (50 kDa vs. 58 kDa), as was one of the bands detected by the antibodies to the HlyD peptide (48 kDa vs. 56 kDa). The bands detected by the antibodies to the predicted peptide from isotig21532 were much larger than the putative Wolbachia homolog (120 and 260 kDa vs. 83 kDa), perhaps due to the incorporation of extra domains from the O. flexuosa genome. Of course, the full amino acid sequences of the O. flexuosa proteins detected by these antibodies are not known at this time.
Antibodies against the KLH carrier protein used for antibody production produced no tissue-specific labeling (Figure 5). In contrast, antibodies raised against the HlyD peptide and against the predicted Wolbachia-like peptide from isotig21532 both produced tissue-level labeling patterns similar to those seen in the in situ hybridization experiments. Antibodies against the HlyD peptide labeled the intestine, hypodermis, lateral chords, and median chords of adult male and female worms (Figure 6A-D and F). They also labeled tissues corresponding to the early stages of sperm development (Figure 6A, B), but did not label mature spermatozoa (Figure 6C). In microfilariae producing females, intense labeling was observed in oocytes (Figure 6E), early embryos (Figure 6F), and regions of microfilariae (Figure 6H, I). Less intense labeling was noted in morula (Figure 6F) and pretzel stages (Figure 6G). The labeling pattern produced by antibodies against the predicted peptide from isotig21532, which shows sequence similarity to a Wolbachia hypothetical protein from O. volvulus, was nearly identical except that nuclei were more intensely labeled in the stained tissues (Figure 7). With the exception of the male hypodermis, which was not labeled by the RNA probe, it appears likely that these two peptides are mostly retained in the tissues where the corresponding RNAs were synthesized. The similarities between in situ hybridization and immunolabeling results make it seem less likely that the multiple bands seen in Western blots probed with these antibodies reflect non-specific labeling or cross-reactivity.
Antibodies against the predicted peptide from isotig12596, a peptide similar to Wolbachia aminopeptidase P, produced a labeling pattern somewhat different from the other two antibodies we examined. In this case, strong labeling was seen in the hypodermis in both sexes, but was not as distinct in the lateral chords (Figure 8A, D, E, G). Some labeling was also seen in the basal layers of the muscle in males (Figure 8A, D). We reported similar results for a putative Wolbachia-like LolC peptide; while the RNA localization pattern was identical to that seen in Figure 2, immunohistology suggested that the corresponding peptide was restricted to somatic muscles . Again, antibodies against the predicted peptide from isotigs 12596 labeled spermatogonia (Figure 8A, B), but no labeling was detected in spermatocytes (Figure 8A, C), spermatids, or mature spermatozoa (Figure 8D). This could indicate that the corresponding RNA detected in spermatocytes is not translated into protein or that the protein is quickly degraded. The antibody binding patterns seen in female worms and developing offspring were similar to those reported for the other two antibodies (i.e., labeled ovaries, oocytes, early embryos and microfilariae, but not morulae) except that labeling was apparent in the uterine wall of mature females rather than that of immature females (Figure 8E, G, H, I).
It is not surprising that Wolbachia-like proteins were detected in tissues or sub-cellular compartments where the corresponding mRNA was not found (e.g., somatic muscles, uterus of mature females, nucleus vs. cytoplasm, etc). Many proteins are produced in a specific tissue (or subset of tissues) but have important functional roles in other regions of the body. In fact, Wolbachia gene products are believed to be exported from the endobacteria to interact with host cells, and some of them have been detected among the excretory/secretory products of Wolbachia-dependent worms .
Different patterns of localization observed for different putative Wolbachia-like peptides could indicate that the proteins containing these peptides are responsible for different functions in the worms. Darby et al.  proposed a dual role for Wolbachia in filarial biology based on their recent genomic analysis of the Wolbachia endosymbiont of the cattle parasite O. ochengi; they speculated that Wolbachia may engage in metabolic provisioning and also serve as a diversion for the vertebrate host immune system. If we hypothesize that Wolbachia-like sequences in O. flexuosa function as a substitute for the endosymbiont, the localization patterns of proteins with Wolbachia-like sequences may provide clues as to which of these two tasks they are most likely to perform. For example, proteins present in somatic muscles and reproductive tissues may be involved in worm metabolism, perhaps granting O. flexuosa novel biochemical capabilities compared to other, Wolbachia-dependent filarial species (e.g., de novo heme, nucleotide or riboflavin synthesis [29, 30]). Likewise, proteins with transmembrane domains or secretion peptides that are present in the hypodermis or lateral chords may be released from the worm to distract the host’s immune system by promoting an ineffective Th1 type immune response . Further work will be needed to define the roles of these proteins in the biology of Wolbachia-free filarial parasites like O. flexuosa.
The results of this study further support the notion that Wolbachia- like sequences are expressed in a Wolbachia-free filarial parasite at both the RNA and protein levels; this expression is highly regulated with regard to tissue and parasite stage. Thus far, all of the Wolbachia-like RNA probes that have successfully labeled O. flexuosa tissues have produced the same labeling pattern: intestine and lateral chords of both sexes, the hypodermis, median chords and uteri of females, and sperm precursors in the male testis. Polyclonal antibodies against predicted Wolbachia-like peptides bound to specific bands in Western blots performed with O. flexuosa adult worm lysate and to specific tissues in fixed and sectioned parasite specimens. Wolbachia-like peptides were sometimes found in the same tissues where the corresponding RNAs were produced and sometimes in other locations. The localization of these proteins within the worm may provide clues regarding their functions by suggesting a role in worm metabolism or in host immune modulation.
Taylor MJ, Bandi C, Hoerauf A: Wolbachia bacterial endosymbionts of filarial nematodes. Adv Parasitol. 2005, 60: 245-284.
Hoerauf A, Nissen-Pahle K, Schmetz C, Henkle-Duhrsen K, Blaxter ML, Buttner DW, Gallin MY, Al-Qaoud KM, Lucius R, Fleischer B: Tetracycline therapy targets intracellular bacteria in the filarial nematode Litomosoides sigmodontis and results in filarial infertility. J Clin Invest. 1999, 103: 11-18. 10.1172/JCI4768.
Slatko BE, Taylor MJ, Foster JM: The Wolbachia endosymbiont as an anti-filarial nematode target. Symbiosis. 2010, 51: 55-65. 10.1007/s13199-010-0067-1.
Rao RU, Moussa H, Weil GJ: Brugia malayi: effects of antibacterial agents on larval viability and development in vitro. Exp Parasitol. 2002, 101: 77-81. 10.1016/S0014-4894(02)00019-X.
Bandi C, Anderson TJ, Genchi C, Blaxter ML: Phylogeny of Wolbachia in filarial nematodes. Proc Biol Sci. 1998, 265: 2407-2413. 10.1098/rspb.1998.0591.
Casiraghi M, Bain O, Guerrero R, Martin C, Pocacqua V, Gardner SL, Franceschi A, Bandi C: Mapping the presence of Wolbachia pipientis on the phylogeny of filarial nematodes: evidence for symbiont loss during evolution. Int J Parasitol. 2004, 34: 191-203. 10.1016/j.ijpara.2003.10.004.
Ferri E, Bain O, Barbuto M, Martin C, Lo N, Uni S, Landmann F, Baccei SG, Guerrero R, de Souza Lima S: New Insights into the Evolution of Wolbachia Infections in Filarial Nematodes Inferred from a Large Range of Screened Species. PLoS One. 2011, 6: e20843-10.1371/journal.pone.0020843.
McNulty SN, Foster JM, Mitreva M, Dunning Hotopp JC, Martin J, Fischer K, Wu B, Davis PJ, Kumar S, Brattig NW: Endosymbiont DNA in endobacteria-free filarial nematodes indicates ancient horizontal genetic transfer. PLoS One. 2010, 5: e11029-10.1371/journal.pone.0011029.
McNulty SN, Mullin AS, Vaughan JA, Tkach VV, Weil GJ, Fischer PU: Comparing the Mitochondrial Genomes of Wolbachia-Dependent and Independent Filarial Nematode Species. BMC Genomics. 2012, 13: 145-10.1186/1471-2164-13-145.
McNulty SN, Abubucker S, Simon GM, Mitreva M, McNulty NP, Fischer K, Curtis KC, Brattig NW, Weil GJ, Fischer PU: Transcriptomic and Proteomic Analyses of a Wolbachia-Free Filarial Parasite Provide Evidence of Trans-Kingdom Horizontal Gene Transfer. PLoS One. 2012, 7: e45777-10.1371/journal.pone.0045777.
Dunning Hotopp JC, Clark ME, Oliveira DC, Foster JM, Fischer P, Munoz Torres MC, Giebel JD, Kumar N, Ishmael N, Wang S: Widespread lateral gene transfer from intracellular bacteria to multicellular eukaryotes. Science. 2007, 317: 1753-1756. 10.1126/science.1142490.
Fischer K, Beatty WL, Jiang D, Weil GJ, Fischer PU: Tissue and Stage-Specific Distribution of Wolbachia in Brugia malayi. PLoS Negl Trop Dis. 2011, 5: e1174-10.1371/journal.pntd.0001174.
Landmann F, Foster JM, Slatko B, Sullivan W: Asymmetric Wolbachia segregation during early Brugia malayi embryogenesis determines its distribution in adult host tissues. PLoS Negl Trop Dis. 2010, 4: e758-10.1371/journal.pntd.0000758.
Kozek WJ: Transovarially-transmitted intracellular microorganisms in adult and larval stages of Brugia malayi. J Parasitol. 1977, 63: 992-1000. 10.2307/3279832.
Kozek WJ, Marroquin HF: Intracytoplasmic bacteria in Onchocerca volvulus. AmJTrop Med Hyg. 1977, 26: 663-678.
McLaren DJ, Worms MJ, Laurence BR, Simpson MG: Micro-organisms in filarial larvae (Nematoda). Trans R Soc Trop Med Hyg. 1975, 69: 509-514. 10.1016/0035-9203(75)90110-8.
Landmann F, Bain O, Martin C, Uni S, Taylor MJ, Sullivan W: Both asymmetric mitotic segregation and cell-to-cell invasion are required for stable germline transmission of Wolbachia in filarial nematodes. Biology Open. 2012, 1: 536-547. 10.1242/bio.2012737.
Plenge-Bonig A, Kromer M, Buttner DW: Light and electron microscopy studies on Onchocerca jakutensis and O. flexuosa of red deer show different host-parasite interactions. Parasitol Res. 1995, 81: 66-73. 10.1007/BF00932419.
Michalski ML, Griffiths KG, Williams SA, Kaplan RM, Moorhead AR: The NIH-NIAID Filariasis Research Reagent Resource Center. PLoS Negl Trop Dis. 2011, 5: e1261-10.1371/journal.pntd.0001261.
Martin J, Abubucker S, Heizer E, Taylor CM, Mitreva M: Nematode.net update 2011: addition of data sets and tools featuring next-generation sequencing data. Nucleic Acids Res. 2012, 40: D720-D728. 10.1093/nar/gkr1194.
Jiang D, Li BW, Fischer PU, Weil GJ: Localization of gender-regulated gene expression in the filarial nematode Brugia malayi. Int J Parasitol. 2008, 38: 503-512. 10.1016/j.ijpara.2007.09.010.
Buttner DW, Wanji S, Bazzocchi C, Bain O, Fischer P: Obligatory symbiotic Wolbachia endobacteria are absent from Loa loa. Filaria J. 2003, 2: 10-10.1186/1475-2883-2-10.
Fath MJ, Kolter R: ABC transporters: bacterial exporters. Microbiol Rev. 1993, 57: 995-1017.
Bennuru S, Meng Z, Ribeiro JM, Semnani RT, Ghedin E, Chan K, Lucas DA, Veenstra TD, Nutman TB: Stage-specific proteomic expression patterns of the human filarial parasite Brugia malayi and its endosymbiont Wolbachia. Proc Natl Acad Sci USA. 2011, 108: 9649-9654. 10.1073/pnas.1011481108.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.
Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M: KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012, 40: D109-D114. 10.1093/nar/gkr988.
Tatusov RL, Fedorova ND, Jackson JD, Jacobs AR, Kiryutin B, Koonin EV, Krylov DM, Mazumder R, Mekhedov SL, Nikolskaya AN: The COG database: an updated version includes eukaryotes. BMC Bioinforma. 2003, 4: 41-10.1186/1471-2105-4-41.
Tatusov RL, Koonin EV, Lipman DJ: A genomic perspective on protein families. Science. 1997, 278: 631-637. 10.1126/science.278.5338.631.
Foster J, Ganatra M, Kamal I, Ware J, Makarova K, Ivanova N, Bhattacharyya A, Kapatral V, Kumar S, Posfai J: The Wolbachia genome of Brugia malayi: endosymbiont evolution within a human pathogenic nematode. PLoS Biol. 2005, 3: e121-10.1371/journal.pbio.0030121.
Darby AC, Armstrong SD, Bah GS, Kaur G, Hughes MA, Kay SM, Koldkjaer P, Radford AD, Blaxter ML, Tanya VN: Analysis of gene expression from the Wolbachia genome of a filarial nematode supports both metabolic and defensive roles within the symbiosis. Genome Res. 2012, 22 (12): 2467-2477. 10.1101/gr.138420.112.
Ghedin E, Wang S, Spiro D, Caler E, Zhao Q, Crabtree J, Allen JE, Delcher AL, Guiliano DB, Miranda-Saavedra D: Draft genome of the filarial nematode parasite Brugia malayi. Science. 2007, 317: 1756-1760. 10.1126/science.1145406.
Scott AL, Ghedin E: The genome of Brugia malayi - all worms are not created equal. Parasitol Int. 2009, 58: 6-11. 10.1016/j.parint.2008.09.003.
Bennuru S, Semnani R, Meng Z, Ribeiro JM, Veenstra TD, Nutman TB: Brugia malayi excreted/secreted proteins at the host/parasite interface: stage- and gender-specific proteomic profiling. PLoS Negl Trop Dis. 2009, 3: e410-10.1371/journal.pntd.0000410.
The authors would like to thank Katherine Mann for technical assistance and C. Dohr and E. Bach for organizing the deer hunt and granting us access to the hunted deer. This work was supported by a grant from the Barnes Jewish Hospital Foundation.
The authors declare that they have no competing interests.
SM analyzed the O. flexuosa transcriptome and proteome, performed bioinformatic analyses and BLAST searches, designed and constructed RNA probes, assisted in the selection of epitopes for antibody production, interpreted results, and drafted the manuscript. KF fixed, embedded and sectioned O. flexuosa nodules and performed in situ and immunohistochemical labeling studies. KC produced the Western blots. NWB identified, collected and supervised the preservation and preparation of O. flexuosa material and edited the manuscript. GW and PF designed and coordinated the study and helped to interpret the results and draft the manuscript. All authors read and approved the final version of the manuscript.
Electronic supplementary material
Additional file 1: Table S1. Relative expression of isogroups derived from the O. flexuosa transcriptome assembly, reported as sequence reads per kilobase of length. (XLSX 812 KB)
Additional file 2: Table S2. Primers used in the construction of in situ hybridization probes. Wolbachia-like sequences identified from the transcriptome of O. flexuosa were amplified from O. flexuosa cDNA. A homolog of the HlyD gene from the Wolbachia endosymbiont of Culex quinquefasciatus was identified in the Wolbachia endosymbiont of B. malayi. This sequence was amplified from B. malayi genomic DNA which also contains DNA from the Wolbachia endosymbiont. (DOCX 41 kb) (DOCX 42 KB)
Additional file 3: Table S3. Target epitopes of anti-peptide antibodies. Epitopes predicted to exhibit favorable immunogenic properties were chosen from regions of O. flexuosa peptide translations with sequence similarity to Wolbachia. A favorable epitope was also chosen from the HlyD gene from the Wolbachia endosymbiont of Culex quinquefasciatus (wCq). The underlined portion of this epitope matches one of the two peptides identified in our proteomic analysis of O. flexuosa worm lysate that mapped to HlyD from wCq. (XLSX 39 kb) (XLSX 39 KB)
Additional file 4: Figure S1. Localization of a putative Wolbachia-like HlyD transcript in adult O. flexuosa. The sense RNA probe (negative control) produced no signal in tissues in either sex (A, C). In female worms, the antisense RNA probe against Wolbachia HlyD labeled the hypodermis, lateral chords, median chords, intestine and uterus of young females (B). In male worms, labeling was seen in the lateral and median chords and in spermatocytes within the testis (D) but not in mature spermatozoa in the vas deferens (E). Abbreviations: cu, cuticle; hy, hypodermis; m, muscle; lc, lateral chords; mc, median chords; i, intestine; ut, uterus; t, testis; vd, vas deferens. Scale bar = equals 50μm. (TIFF 9072 kb) (TIFF 9 MB)
Additional file 5: Figure S2. Localization of a transcript with sequence similarity to a Wolbachia hypothetical protein in adult O. flexuosa. The sense RNA probe (negative control) produced no signal in either sex (A, C). In female worms, the antisense RNA probe against isotigs21532, a sequence similar to that of a hypothetical protein from the Wolbachia endosymbiont of O. volvulus, showed light labeling of the hypodermis, lateral chords, median chords, intestine and uterus (B). In male worms, labeling was seen in the lateral and median chords and in the germinal zone of the testis (D). Abbreviations: cu, cuticle; hy, hypodermis; m, muscle; lc, lateral chords; mc, median chords; i, intestine; ut, uterus; t, testis. Scale bar = equals 50μm. (TIFF 8374 kb) (TIFF 8 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
McNulty, S.N., Fischer, K., Curtis, K.C. et al. Localization of Wolbachia-like gene transcripts and peptides in adult Onchocerca flexuosa worms indicates tissue specific expression. Parasites Vectors 6, 2 (2013). https://doi.org/10.1186/1756-3305-6-2
- horizontal gene transfer
- in situ hybridization