Effects of Doxycycline on gene expression in Wolbachia and Brugia malayi adult female worms in vivo
© Rao et al; licensee BioMed Central Ltd. 2012
Received: 19 October 2011
Accepted: 9 February 2012
Published: 9 February 2012
Most filarial nematodes contain Wolbachia symbionts. The purpose of this study was to examine the effects of doxycycline on gene expression in Wolbachia and adult female Brugia malayi.
Brugia malayi infected gerbils were treated with doxycycline for 6-weeks. This treatment largely cleared Wolbachia and arrested worm reproduction. RNA recovered from treated and control female worms was labeled by random priming and hybridized to the Version 2- filarial microarray to obtain expression profiles.
Results and discussion
Results showed significant changes in expression for 200 Wolbachia (29% of Wolbachia genes with expression signals in untreated worms) and 546 B. malayi array elements after treatment. These elements correspond to known genes and also to novel genes with unknown biological functions. Most differentially expressed Wolbachia genes were down-regulated after treatment (98.5%). In contrast, doxycycline had a mixed effect on B. malayi gene expression with many more genes being significantly up-regulated after treatment (85% of differentially expressed genes). Genes and processes involved in reproduction (gender-regulated genes, collagen, amino acid metabolism, ribosomal processes, and cytoskeleton) were down-regulated after doxycycline while up-regulated genes and pathways suggest adaptations for survival in response to stress (energy metabolism, electron transport, anti-oxidants, nutrient transport, bacterial signaling pathways, and immune evasion).
Doxycycline reduced Wolbachia and significantly decreased bacterial gene expression. Wolbachia ribosomes are believed to be the primary biological target for doxycycline in filarial worms. B. malayi genes essential for reproduction, growth and development were also down-regulated; these changes are consistent with doxycycline effects on embryo development and reproduction. On the other hand, many B. malayi genes involved in energy production, electron-transport, metabolism, anti-oxidants, and others with unknown functions had increased expression signals after doxycycline treatment. These results suggest that female worms are able to compensate in part for the loss of Wolbachia so that they can survive, albeit without reproductive capacity. This study of doxycycline induced changes in gene expression has provided new clues regarding the symbiotic relationship between Wolbachia and B. malayi.
KeywordsDoxycycline Brugia malayi Wolbachia Filariasis Gene expression Microarray
Like many filarial nematodes, Brugia malayi contain intracellular Wolbachia bacteria (w Bm) [1–3]. These obligatory endosymbionts are essential for larval growth and development and for adult worm reproduction and survival [4–7]. They multiply by binary fission and are vertically transferred to successive generations of worms. Metabolic interdependence of endosymbionts and their hosts is common [8–10]. Two fundamental questions remain in this mutualistic relationship: what do w Bm contribute to the nematodes and what do the worms provide for the bacteria?
w Bm have been considered to be important biological targets for drug development against filarial nematodes [11, 12]. Many antibiotics have been tested in vitro and in vivo against w Bm and filarial worms. Of these, tetracycline and related drugs have been shown to be effective in vitro, in animal models, and in humans against filarial infections [6, 13–18]. Tetracyclines also have prophylactic activity, and they can cause a male-biased sex-ratio [19–21]. Doxycycline (Doxy) reduces microfilaremia and adult worm burdens by sterilizing adult worms and eventually killing them (in filarial species containing Wolbachia), either alone or in conjunction with diethylcarbamazine (DEC) or ivermectin [22–26]. In many of these studies, antibacterial effects of the drugs precede the effects on worms. In contrast, microfilariae and adult worms of filarial species that lack Wolbachia such as Acanthocheilonema viteae and Loa loa are not affected by tetracyclines [27, 28]. Therefore the effect of tetracyclines in filarial worms with w Bm is believed to be related to their effect on the bacteria and not due to direct effects on the worms. Some authors have suggested that w Bm may provide key molecules or metabolic functions that are essential for survival and reproduction in species that are w Bm dependent [29, 30]. Doxycycline primarily affects bacterial growth and survival by inhibiting the binding of amino-acyl tRNA in ribosomes which blocks protein synthesis [31, 32]. At higher concentrations, doxy has been shown to affect mitochondrial function [33, 34].
Despite their occurrence in several filarial nematodes and their emergence as a therapeutic target, little is known regarding biological interactions between w Bm and filarial worms [29, 30]. A differential display PCR-based approach has been used to look at changes in gene expression in Litomosoides sigmodontis after tetracycline treatment . A number of genes showed increased expression after doxy including a phosphate permease gene (Ls-ppe-1) which may be essential for embryo development. The genomes of w Bm and Bm provided more clues regarding potential mechanisms of interdependence [29, 36, 37]. For instance, comparative genome studies suggested that w Bm contains pathways for production of riboflavin, flavin adenine dinucleotide, heme, and nucleotides that are lacking in B. malayi. In addition, since the w Bm genome lacks complete biochemical pathways for de novo synthesis of biotin, coenzyme A, NAD, ubiquinone and folate, filarial worms may provide these and other molecules that are required for bacterial growth. w Bm also lack synthetic pathways for most amino acids , and they probably obtain amino acids from their filarial hosts .
In this study we have examined the effects of 6-week doxy treatment on w Bm and Bm, and used V2 Bm microarray to observe changes in expression profiles of many genes; a subset of these genes may play a role in survival and reproduction of female worms.
All animals were handled in strict accordance with good animal practice as defined by the Animal Welfare Act, the Guide for the Care and use of Laboratory Animals and guidelines of the Division of Comparative Medicine, Washington University School of Medicine.
Infections and parasites
Mongolian gerbils (Meriones unguiculatus) (Charles River Laboratories, Wilmington, MA) were inoculated subcutaneously (s.c.,) with 150 B. malayi L3 larvae as described previously . Adult worms at 180 days post-infection were recovered from testis, spermatic cord lymphatics, heart and lungs. Female worms were separated from males, washed in RPMI-1640, and were snap frozen in dry ice/ethanol and stored at -80 C. All worms recovered from s.c., infections used in this study were of the same age group but derived from different gerbils.
In vitro culture of female worms
Brugia malayi female worms from treated and untreated gerbils were cultured in vitro to assess microfilaria (MF) release, and the number of released MF were counted on days 2, 4, and 6 of culture . Briefly, 5 worms were incubated per well in 6-well culture plates (Costar, Cambridge, MA) in 5 ml of complete medium (RPMI-1640 supplemented with 25 mM HEPES buffer, 2 mM glutamine and 10% fetal calf serum) at 37 C in a 95% air-5% CO2 atmosphere.
Doxycycline treatment and worm recovery
Subcutaneously infected gerbils (n = 6) with circulating microfilaria at approximately 140 days post infection were treated by oral gavage with doxy (Sigma Chemical Company, St. Louis, MO) in distilled water at 100 mg/kg body weight, daily for 6 weeks. Controls were untreated. At the end of the treatment, animals were necropsied, and worms were removed from animal tissues.
Microfilaria counts and embryograms
MF in 10 μl of tail-blood of treated and untreated animals were counted in duplicate after 1, 2, 4 and 6 weeks of treatment . Brugia malayi female worms recovered after 6 weeks of doxy treatment and control worms were studied to produce embryograms as previously described . Differences between embryograms of treated and control worms were analysed by One-way ANOVA.
Female worms from treated and untreated gerbils were fixed in 100 mM phosphate with 2% paraformaldehyde and 2.5% glutaraldehyde and post-fixed in 1% osmium tetroxide in phosphate buffer. Sections were stained, dehydrated in ethanol and embedded in Eponate 812 (Ted Pella, Inc., Redding, CA). Ultrathin sections (70-90 nm) were stained with uranyl acetate and lead citrate and viewed with a JOEL 1200 transmission electron microscope .
RNA isolation, cDNA probe synthesis, and microarray analysis
Two groups of 15 female worms each from doxy treated and untreated animals (biological replicates) were crushed, and total RNA was extracted separately using TRIzol (Invitrogen Corporation, Carlsbad, CA) and treated with DNase (Ambion, Austin, TX) as per the manufacturer's instructions. RNA integrity and quantity were measured with a Bioanalyzer-model 2100 (Agilent technologies Inc., Palo Alto, CA).
The Version-2 Filarial Microarray contains 18,104 oligonucleotide elements (65mer) that uniquely represent coding sequences for 15,412 B. malayi genes or predicted genes, 804 Wolbachia genes, 1,016 Onchocerca volvulus genes, and 872 Wuchereria bancrofti genes. A recent publication provides extensive annotation for this array ; see also http://www.nematode.net).
In a preliminary experiment, we synthesized cDNA from RNA samples using oligo(dT) primers and Superscript II Reverse Transcriptase (Gibco BRL, Gaithersburg, MD). Hybridizations, scanning and data analysis were performed as previously described [40, 41]. Oligo(dT) priming has been the most widely used method for conversion of mRNA into cDNA through reverse transcription [42, 43]. However, these experiments did not provide optimal hybridization signals for bacterial transcripts which lack polyA. Therefore, later experiments employed a random priming method for probe synthesis. Briefly, cDNA was synthesized from 2 μg of RNA using a Genisphere kit (Hatfield, PA) according to the manufacturer's instructions (http://www.genisphere.com/pdf/array900mpx_protocol_v06-22-04.pdf). In order to have technical replicates, cDNA probes were dye labeled separately with Cy5 and Cy3 as per the manufacturer's instructions (PerkinElmer Life and Analytical Sciences, Waltham, MA).
Dye-labeled cDNA probes were hybridized to duplicate 65 mer oligos on the V2 filarial array. Two biological replicates and two technical replicates provided 8 hybridizations for each filarial oligomer on the array. cDNA hybridizations were done overnight at 43 C with the 3DNA Array 900 detection system. The slides were washed, and scanning and gridding were performed with a ScanArray Express HT Scanner-v3.0 (Perkin Elmer, Boston, MA). Photomultiplier tube (PMT) values were set to 68 and 59 volts for Cy3 and Cy5, respectively. An additional scan was done for each slide with the PMT set for 56 and 48 volts. The scanned fluorescence intensity values were Lowess normalized and analyzed by Genespring v7.2 (Agilent) and Partek Genomics Suite software (Partek, St. Louis, MO). An average of four arrays for each condition were considered and w Bm and Bm elements were scored as "present" if the signal was ≥ 250 or if the signal to background ratio was ≥ 2. All microarray experiments were performed in agreement with the MIAME guidelines. The microarray data sets have been deposited into GEO (accession ID; GSE34976).
Genes with fold differences in expression equal to or greater than two fold and a confidence level of 99% (P < 0.01, Student's t-test) in pair-wise comparison between treated and control worms were considered to be differentially expressed.
Gene expression analysis by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR)
Known w Bm gene sequences for 16s rRNA (w Bm9003), cell cycle gene-fts Z (w Bm0602), surface protein-wsp (AJ252061), seven heme genes (w Bm0133, w Bm 0373, w Bm0777, w Bm0728, w Bm 0001, w Bm0709 and w Bm0719), and B. malayi sequences for sheath protein-shp1 (U43568), embryonic fatty acid binding protein-FABP1 (AF178439), and actin (XM_001895760), were used to amplify the corresponding transcripts in treated and control worm RNA by qRT-PCR. These experiments were performed twice with two batches of worms. In addition, a subset of differentially expressed w Bm and Bm genes from the microarray analysis was selected for assay confirmation studies by qRT-PCR. The expression profiles for bm.03044, BMX1802, BMX9848, BMX8408, BMX9336, BMX11510, BMX4544, bm.00532, bm.01661, BMX2449 genes affected after antibiotic treatment (both in our array study and in previous reports) were also analyzed by RT-PCR using 6 wks treated vs. untreated worm cDNA. Primer pairs were designed from sequences obtained from the w Bm genome http://tools.neb.com/wolbachia/) or from GenBank with Primer Express 3.0 software (Applied Biosystems, Foster city, CA). Sequence specific primer pairs were purchased from Integrated DNA Technologies, Coralville, IA. Briefly, cDNA was synthesized using 1 μg of total RNA, random primers, and SuperScript II reverse transcriptase (Invitrogen, Carlsbad, CA). qRT-PCR reactions were performed in duplicate in 96-well optical plates in a 25 μl reaction volume containing 1 ng of cDNA with individually optimized primer concentrations and SYBR green master mix (Applied Biosystems). To ensure that there was no genomic DNA contamination in cDNA and no primer dimer artifacts, reactions containing templates generated without reverse transcriptase enzyme were included as controls. Water was also tested in each reaction as a "no template" control (NTC). Reactions were conducted with an initial step of 2 min at 50 C and 10 min at 95 C, followed by 40 cycles of 15 s at 95 C and 1 min at 60 C. A melting curve analysis was performed at the end of each cycle to assess the amplification specificity. Cycle threshold (Ct) values were obtained for the target and internal control gene. Brugia malayi histone gene (Pub_locus Bm1_49345) Ct values were used as an internal control to normalize transcript levels present in different samples. Real-time PCR efficiency (E = 10 - (1/slope)-1) was monitored by generating standard curves for Bm histone using log template concentration of cDNA from treated and untreated worms and found to be close to 100%. The relative quantitation method (2 -ΔΔCt) was used to compare expression levels of w BM and Bm genes in different samples  (http://www3.appliedbiosystems.com/cms/groups/mcb_support/documents/generaldocuments/cms_040980.pdf).
For correlation analysis, qRT-PCR results (2-ΔΔ Ct) of selected w Bm and Bm genes and their corresponding microarray results (normalized ratios) were Log10 transformed, and the results were analyzed by the nonparametric Spearman rank correlation test using GraphPad Prism v.4 software (GraphPad Software Inc, San Diego, CA). qRT-PCR results for individual genes were analysed by one-way ANOVA to assess the significance of differences between treated and control worms.
Functional classification of differentially expressed gene transcripts
Functional analysis was carried out for differentially expressed w Bm and Bm genes to identify overrepresented pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG v.43) . For each query, the top match (with P < 1.0e-10) and all the matches within 30% of the top Blast score were accepted for KEGG Ontology (KO) and pathway associations. Default parameters for InterPro Scan v13.1 were used to map against InterPro database . These mappings were used to classify the gene transcripts as biological, molecular, and cellular components. Gene Ontology mappings for w Bm genes were selected from UniProtKB database (http://www.uniprot.org/uniprot/) and from InterProScan v13.1. Statistically enriched pathways and gene ontology (GO) categories were obtained using hypergeometric analysis with P < 0.05 considered significant . For non-redundant (NR) nucleotide hits, differentially expressed genes were BLAST searched (e-value cutoff of 1.0e-10) against a NR database with no Caenorhabditis sequences. In addition, we used Wormpep v.180 with an e-value cutoff of 1.0e-05 for detecting Brugia malayi homologs of C. elegans genes having RNAi phenotypes.
Effect of doxycycline on Brugia malayi worms and Wolbachia
In vivo effects of doxycycline on Brugia malayi microfilaria (MF) and adult worms
% of pretreatment MF counts (mean ± SD)
Mean (± SD) number of worms
(% reduction in worms)
Days after treatment
121 ± 74
119 ± 30
200 ± 95
192 ± 101
167 ± 135
90 ± 18
143 ± 108
112 ± 92
29.3 ± 18 (15)
34.5 ± 10
Effects of doxycycline treatment on Wolbachia gene expression
The V2 Filarial array contains 804 w Bm genes, and hybridization signals were detected for 679 (84.5%) in control female worm cDNA (see additional file 2). Expression signals were detected for 538 (67%) of w Bm genes in doxy treated worm RNA (18% fewer than in untreated worms).
Of 679 w Bm genes with detectable expression signals in control worms, the top twenty with the highest expression signals (median values 3.08E+03 to 1.02E+04) included genes involved in protein synthesis (6), stress response (4), antioxidant (1), energy and electron transfer (1), copper acquisition and aerobic respiration (1) lipid metabolism (2), DNA binding or repair (2), host interaction and adaptation (3) and hypothetical (4) (see additional file 3). Expression signals for 140 w Bm elements were present in control worms and absent in treated worms. The top 10 genes that showed highest expression signals in control worms but 'absent' in treated worms include molecular chaperonins, transcription and translation factors, membrane and ribosomal proteins and antioxidants (see additional file 4). In contrast, signals for 98 w Bm elements were present in treated worms but absent in control worms (see additional file 4). The top 10 genes with highest expression signals here include transport or lipoprotein release (w Bm 0483), secretion (w Bm 0796), proteolysis (w Bm 0770), carbohydrate and aminoacid metabolism (w Bm 0666), cell wall biogenesis and peptidoglycan biosynthesis (w Bm 0508), nucleic acid repair (w Bm 0746), and integral membrane protein TerC (w Bm 0442) which may be involved in resistance to antibiotics.
The scatter plots displaying the mean normalized fluorescence intensity signals for all elements on the array (panel A) and for w Bm elements alone (panel B) met the criteria for gene expression in treated vs. control samples (see additional file 5).
Description of top 20 differentially expressed Wolbachia genes found in doxycycline treated Brugia malayi female worms
DNA-directed RNA polymerase sigma 32 subunit, RpoH
Ribosomal protein S14
Ribosomal protein S21
Ribosomal protein L1
Ribosomal protein S17
Ribosomal protein L16
Ribosomal protein S12
DNA replication, DNA recombination, DNA repair or SOS response
DNA polymerase III, gamma/tau subunit
Predicted EndoIII-related endonuclease
Superfamily II DNA/RNA helicase
Energy production and conversion
F0F1-type ATP synthase, beta subunit
Metabolism & Transport
Succinate dehydrogenase flavoprotein subunit, SdhA
Kef-type K+ transport system, membrane component
HesB/YadR/YfhF family protein
Type IV secretory pathway, VirB6 components
Immunogenic membrane components
Outer surface protein-wsp
KEGG data search revealed that out of 804 w Bm genes on the microarray, 499 were mapped to many pathways ranked by gene abundance, and 487 had KEGG ontology (KO) assignments (see additional file 7). Many mapped genes were associated with pathways in metabolism (221), genetic information and processing (215), environmental information and processing (46) and cellular processes (13). Out of 200 w Bm genes that were differentially expressed after doxy treatment, 141 (70%) genes were mapped to pathways, and significantly enriched pathways (P ≤ 0.05) associated with down-regulated genes include bacterial ribosome, translation factors, and the pentose phosphate pathway in carbohydrate metabolism (Table 2, see additional file 7). In contrast, two w Bm genes (w Bm 0474, w Bm 0365) with increased expression after treatment mapped to pathways involved in energy metabolism, nucleic acid replication/protein repair, respectively (see additional file 7). The third up-regulated gene (w Bm0149) is a hypothetical gene and its function is unknown. GO term assignments were concordant with KEGG results, and many biological processes and molecular functions in w Bm were down-regulated after doxy treatment (see additional file 6). InterPro scans showed a total of 366 protein domains mapped in w Bm gene transcripts, and only two of these domains related to nucleic acid binding and to an unknown function were significantly (P < 0.05) enriched in down-regulated w Bm genes (see additional file 8). The protein domains mapped to the up-regulated w Bm 0474, w Bm 0365 genes are related to the respiratory chain proteins NADH-quinone oxidoreductase and dehydrogenease subunits and DNA repair proteins.
Effect of doxycycline treatment on Brugia malayi (Bm) gene expression
Of 15,412 Bm elements on the array, 13,399 had hybridization signals in control adult female worms. Signals were detected for 14,340 elements in doxy treated worm cDNA (see additional file 2). 53 Bm elements were present in control worms and absent in treated worms. In contrast, 994 Bm elements were present in treated but absent in control worms. Of 53 elements uniquely present in control worms, elements for 15 genes have no annotation in the database and 9 were hypothetical (see additional file 9). Known genes included cytoskeletal genes (actin, myosin), the antioxidant glutathione S-transferase (GST), and ribosomal proteins (S1 and L11). KEGG analysis showed pathways associated with metabolism (biosynthesis of secondary metabolites, xenobiotic metabolism and degradation, and glycan biosynthesis), cellular processes (cell motility and communication, and immune system), and human diseases were enriched significantly (P ≤ 0.05) in control worms (see additional file 10). Of 994 genes with expression signals in treated worms but not in control worms, many were classified as hypothetical (288), ribosomal (8) and novel genes (~480) with no annotation. Known genes in this group included cuticular components, proteases (astacin and ADAMTS-like metalloproteases, cathepsin-L (cpl-1), vespid allergen antigen homolog (Val), immunogens and immunomodulators (Alt-1) (see additional file 9). KEGG pathways significantly (P ≤ 0.05) enriched in treated worms included metabolism of carbohydrates and secondary metabolites, energy metabolism (oxidative phosphorylation), and nucleotide metabolism. In addition, pathways involved in genetic information and processing (translation, folding), signaling pathways and various other cellular processes (cell motility, cell growth and communication) were enriched in treated worms (see additional file 10). Many GO terms assigned for genes present in control worms and treated worms or vice versa corresponded to pathways analysed by KEGG, and only those significantly (P < 0.05) enriched are listed in additional file 10.
Major known differentially expressed* Brugia malayi genes after doxycycline treatment
Cytochrome C oxidase
Ubiquinol-cytochrome-c oxidase activity
ATP synthase F0 subunit 6
Cathepsin L-like cysteine proteinase
Trypsin family protein
Metabolism and transport
Biopterin-dependent aromatic amino acid
Cysteine-rich, acidic integral membrane protein precursor, putative
PPE family protein, putative
Excretory/secretory protein Juv-p120
Vespid venom allergen
antigen-like protein 1
IgG and IgE immunoreactive antigen
Reproduction, development & transport
Embryonic fatty acid-binding protein
Alpha-1 collagen type IX, putative
Peptidyl-prolyl cis-trans isomerase
NADH dehydrogenase subunit 5
KEGG and InterPro analysis of differentially expressed B. malayi genes
Comparison of observed changes in Brugia malayi and Litomosoides sigmodontis gene expression after doxycycline/tetracycline treatment
We compared our results with those recently reported after i.p. tetracycline treatment of Bm and in L. sigmodontis[48, 49]. We have observed expression of more genes to be regulated after doxy treatment in Bm compared to the two previous reports (319 and 323 vs 546 respectively). However, only few genes (1% or 5/546) regulated at 6 wks doxy treatment agreed with differentially expressed genes found in Bm or L. sigmodontis at earlier time points. The expression of vespid venom allergen homolog (BMC00351) while it was down- regulated in Bm at 2 wks of treatment, found up-regulated at 5 and 6 wks in L. sigmodontis and Bm. Genes encoding hypothetical (BMX1802, BMX9848), trafficking protein (BMX 8408) and peptidyl-prolyl cis-trans isomerase (PPI-FKBP type, BMX9336) transcripts were up-regulated at 2 wks but down-regulated at 6 wks of treatment in Bm. In addition, the Bm homologs of seven-transmembrane-helix (7TM) receptor (BMX11510), zinc finger family protein (BMX4544) and 2 hypothetical proteins (bm.00532, bm.01661) found in L. sigmodontis were in agreement with regulated (up) genes at 6 wks. However, expression of nucleotide-binding protein (BMX2449) involved in signaling pathway was down in L. sigmodontis at 5 wks but up at 6 wks in Bm. Our microarray results for these genes were further confirmed by RT-PCR (data not shown).
Homologs of Caenorhabditis elegans genes with known RNAi phenotypes in doxycycline treated B. malayi
A total of 321 genes of the 994 Bm elements present in treated female worms and not in control worms have homologues in C. elegans with RNAi phenotype information. While 228 genes were wild type phenotypes, 93 had abnormal phenotypes. Most of these are involved in reproduction and development including 50 embryonic lethal (Emb), 26 growth (Gro), 35 sterility (Ste) and sterile progeny (Stp) phenotypes (see additional file 9).
272 of 546 (50%) Bm transcripts that were present in both control and treated worms but differentially expressed had homologues in C. elegans (see additional file 11). More of the genes that were up-regulated after doxy (231/272) had homologues than genes that were down-regulated (41/272). RNAi phenotype information was available for C. elegans homologues for 264 of 546 (48%) Bm genes. Of these, 207 homologues were wild type, while 57 (36 up-regulated and 21 down-regulated genes) had severe phenotypes in C. elegans including sterility, larval arrest, embryonic lethality, egg laying defect, sick, slow growth, etc. In addition, 19 down-regulated Bm genes with observed phenotypes have been reported to be female-regulated (Li et al, personal communication). These included four collagen genes (15377.m00007, 14845.m00009, 14296.m00065, 14296.m00063), one histone (13785.m00207), and one hypothetical gene (BMC01609). Many other down-regulated Bm genes are homologues of sterile phenotypes in C. elegans especially Emb, Lva, Lvl, Ste phenotypes.
Gene expression analysis by qRT-PCR
qRT-PCR studies showed significant down-regulation of w Bm 16S, fts Z and wsp expression in doxy treated female worms relative to controls. Down- regulation was more impressive for 16S and wsp genes compared to fts Z (see additional file 14). Embryo associated genes Bm-shp1 and FABP transcripts were significantly (P < 0.05) down-regulated in treated worms, while treatment did not significantly affect actin expression. Interestingly, expression of all 7 w Bm genes in the heme synthesis pathway was significantly down- regulated in treated worms by qRT-PCR (see additional file 14). Again, these results were generally consistent with those obtained by microarray. Two genes (w Bm 0709 and w Bm 0719) had no expression signal by microarray in treated worms. Expression signals for the other 5 genes were decreased in treated worms, although only one of these (w Bm 0133, 5-aminolevulinate synthase, EC 22.214.171.124) met the criterion of a 2-fold reduction with P < 0.01 (see additional file 3).
Prior studies have shown that doxy affects w Bm from filarial worms with disruption of parasite development and reproduction. Ultrastructural changes and changes in cell- surface (wsp) and cell- cycle (fts-z) DNA markers relative to controls indicate that doxy treatment reduces w Bm bacteria in female worms. The purpose of this study was to examine the effects of doxy treatment on gene expression in adult worms and in w Bm. As expected, our results confirmed previous reports regarding the effects of doxy treatment on w Bm, embryogenesis and MF production in adult filarial worms [6, 13, 15, 16, 23, 24, 50, 51]. Although there are reports that doxy kills adult worms in some in vivo systems, our results showed that Brugia malayi worms survived 6 wks of doxy treatment in gerbils. Surviving worms were not normal, and a longer post-treatment observation period might have shown that doxy shortened the lifespan of adult female worms.
We then turned our attention to gene expression. We first tested the efficiency of cDNA probes synthesized by oligo(dT) priming for hybridization and signal detection. This method detected expression signals in adult female worms for 12,791 of 18,101 elements on the array (data not shown). However, this method was not reliable for detecting w Bm genes, which lack poly-A tails. Oligo(dT) primers also have a tendency to produce truncated cDNAs through internal poly(A) priming, and this can reduce the quality of synthesized cDNA [42, 52]. Random priming of total RNA for synthesis of hybridization probes [53, 54] worked well for transcripts of both w Bm and Bm. Untreated control female worms had transcription signals for 84% of w Bm genes, with most genes showing ≥ 3.08 × 103 fluorescence units. Presumably some of these w Bm genes are critical for bacterial survival and for maintenance of the mutualistic relationship association with their filarial host. Many w Bm gene expression signals were absent or significantly reduced after doxy treatment . This finding was expected as doxy treatment kills and clears w Bm from filarial worms. It is unclear why three w Bm genes had increased expression signals after doxy treatment; w Bm0474, associated with energy metabolism and metabolism of cofactor molecules and w Bm0365 which is involved in DNA replication and repair were up-regulated. Similarly, in other studies of gene expression in tetracycline treated worms, expression of many w Bm genes was differentially expressed with 13 genes in B. malayi (at day 14) and 3 genes in L. sigmodontis (at day 36) up-regulated [48, 49]. These results suggest that some w Bm and w Bm DNA persists after at least 6 wks of doxy treatment.
We found over representation of w Bm genes in pathways such as carbohydrate metabolism, energy metabolism, aminoacid metabolism and nucleotide metabolism that were affected by doxycycline treatment. Complete pathways for de novo biosynthesis of nucleotide/nucleoside biosynthetic genes are present in the w Bm genome, and this has been postulated as a basis for Wolbachia dependence in filarial worms . Similarly, computational predictions suggested that w Bm genes (~552 genes) may be essential for certain pathways . We observed that many genes predicted to have these functions were down regulated after doxy treatment.
Two recent reports described effects of tetracycline treatment on gene expression in Bm and L. sigmodontis[48, 49]. Our study confirmed some of these observations (down-regulation of transcripts essential for transcription, translation), but there were also important differences (e.g., reduced expression of genes involved in reproduction and development, see Table 3). For example, Ghedin et al.,  reported overexpression of Bm transcripts involved in amino acid synthesis and protein translation and down-regulation of cuticle synthesis transcripts very soon after tetracycline treatment (14 days). Similarly, while many transcripts involved in transcription and translation, and in metabolic and signaling pathways were down regulated, motility, heme-binding protein, and mitochondrial transcripts encoding energy metabolism were up-regulated after treatment in Litomosoides, and a similar profile for energy genes was reported for Bm .
In addition to the timing of gene expression changes in filarial worms, the experimental design and data analysis may explain differences in microarray results in the present study from those reported in prior studies. For example, the infection method (i.p., vs. s.c.,), medication used, mode and duration of drug treatment, and methods for probe synthesis and gene expression analysis in the present study differed in many ways from those used in the other reports. It is difficult to quantify drug dosage when medications are provided when it is in drinking water, and drugs can have different effects on worms in different locations in the host. In addition, we have used a more stringent criterion for calling differential expression by using > 2 fold change with P < 0.01 value as a parameter. We think this increased stringency is important to reduce false calls for differential gene expression in large microarrays.
Our study found striking increases after treatment in expression of genes encoding proteins involved in immunoregulation, electron transfer and energy metabolism. Many immune evasion gene transcripts were highly expressed in worms (Table 3) . In addition, we found that Brugia mmc1 transcript (reported to be specific for MF stage, see , an immunogen was up-regulated in doxy treated female worms by microarray, and this result was confirmed by qRT-PCR. Activation of immune evasion functions may be needed for treated worms to survive, because the loss of w Bm might otherwise render the worms more susceptible to host immune responses. Up-regulation of mitochondrial genes involved in electron transfer and respiratory chain after treatment suggests that doxy treatment may directly affect targets other than w Bm such as worm mitochondria . This increased expression of mitochondrial transcripts after doxy treatment may partially protect filarial worms from the loss of w Bm so that they can survive. Interestingly, tetracycline has no direct effect on Wolbachia-free A. viteae worms or their mitochondrial respiratory chain transcripts . Additional research will be needed to test the hypothesis that mitochondria and w Bm may have overlapping or complimentary roles in Bm worms.
Heme is an essential factor for development, survival and fertility of many nematodes. The w Bm genome contains all but one of the genes required for heme synthesis [29, 30]. In contrast, only 1 gene in the heme biosynthetic pathway is present in the B. malayi genome (ferrochelatase, EC 126.96.36.199, Bm1_14320). This led to the hypothesis that filarial worms may depend on w Bm for heme synthesis. Recent studies that showed that inhibitors of heme synthesis by Wolbachia reduced worm viability support this hypothesis . Our qRT-PCR studies found decreased expression of all w Bm genes involved in heme synthesis in doxy-treated treated female worms. Microarray results were similar, but less impressive. Recently, Strübing et al.,  postulated that up-regulation of mitochondrial genes was essential for worms to maintain mitochondrial heme-dependent respiratory chain protein complexes for energy and homeostasis without Wolbachia. Our results support this hypothesis.
Many previous studies have reported changes in filarial worm embryogenesis and development after doxy treatment. Wolbachia is known to affect host oogenesis in insects, possibly by affecting expression of genes involved in this process [60, 61]. We observed that after the loss of w Bm, many known Bm gene transcripts essential for cytoskeleton and extracellular matrix synthesis (prolyl 4-hydroxylase, collagen/cuticulin, laminin, and caveolin) and transcripts believed to be involved in embryogenesis and reproduction (embryonic fatty acid binding protein and cathepsin-L) were down-regulated [62–66]. Although, phosphate permease family gene (Ls-ppe-1) thought to be involved in nucleotide metabolism was up-regulated after tetracycline treatment, orthologues in Bm were not up-regulated at 6 wks of treatment. However, other phosphate transport genes in Bm were down-regulated after doxy treatment. This suggests that doxy may significantly affect phosphate transport to Bm and that dysregulation of phosphate and nucleotide metabolism could impair Bm embryo development. Down-regulation of development-associated transcripts like FABP1 and shp1 after doxy by qRT-PCR further supports the array results. Similarly, Ghedin et al.  showed down-regulation of few collagen genes after tetracycline treatment in Bm. Conversely many other collagen collagen genes up-regulated in treated Litomosoides. While other studies reported that expression of troponins and several collagens was up-regulated in treated Litomosoides and Bm worms [49, 67], our results did not confirm these findings.
Zinc-finger domain genes have diverse roles in transcriptional modulation, translational control, protein degradation pathways, and growth regulation; and they are believed to have diverse, important roles in filarial biology [36, 68]. We found that doxy treatment increased expression of these genes in Bm suggesting that protein-nucleic acid and protein-protein interactions may be increased in treated worms after the depletion of w Bm. Based on NR best hits in the protein sequences, transcripts encoding at least 10 Bm genes that were differentially expressed after doxy have strong homology to mosquito (Anopheles gambiae) transcripts involved in DNA integration (14927.m00099), nucleic acid binding (15547.m00008), citric acid cycle (13478.m00071) or transport (15547.m00008, 14979.m04385, 12526.m00074). Sequence similarities to mosquito gene fragments have been reported in the Bm genome and in V2 gene expression studies ; (Li et al., unpublished observations) and in proteomic profiles of tetracycline treated worms  indicating occurrence of DNA integration of mosquitoes with Brugia nematodes.
In C. elegans, RNAi screening has demonstrated abnormal phenotypes for about 5000 genes which is about 25% of all predicted genes [69, 70]. It has been shown that about 30% of predicted genes in Bm with homologues in C. elegans have RNAi phenotypes in that species [36, 37]. Among the genes represented in the V2 filarial array, approximately 37% with homologues in C. elegans have RNAi phenotypes . The increased rate of phenotypes for filarial genes with C. elegans homologues suggest that many of these genes are necessary for core nematode pathways that are conserved across the phylum.
In summary, this study has provided detailed information on effects of doxy treatment on gene expression in w Bm and adult female Bm. As expected, treatment generally reduced w Bm gene expression. However, this was not uniform, and expression of many w Bm genes was not affected by doxy. Doxycycline had a mixed effect on Bm gene expression. Reduced expression of genes involved in reproduction, growth and development was consistent with sterility observed in doxy treated worms. The up-regulated genes and pathways in treated worms were less predictable. Up-regulation of energy production (oxidative phosphorylation, carbohydrate metabolism), electron transport, antioxidants, proteases and immunomodulatory genes may represent parasite adaptations that permit survival (albeit with sterility) when related functions normally provided by w Bm are disrupted by doxy treatment. For example, the energy and electron transport findings suggest that increased mitochondrial activity may help the worms survive the loss of w Bm. It is possible that w Bm normally function as accessory mitochondria providing extra energy in metabolically active tissues. Additional work will be needed to test this interesting hypothesis. The challenge remains to determine which of these changes in gene expression are critical adaptations and which are secondary effects or general responses to stress conditions. It would be interesting to compare doxy effects on gene expression in filarial worms to those caused by other antifilarial drugs.
List of abbreviations used
quantitative reverse transcription polymerase chain reaction
Koyoto Encyclopedia of Genes and Genomes
We would like to thank the Brugia malayi Microarray Consortium for bioinformatics work and for sharing the expense of constructing the V2 filarial array. We would like to thank Dr. W. Beatty and Ms. D. Gill, Washington University School of Medicine for electron microscopy. This study was supported in part by a grant from the Barnes-Jewish Hospital Foundation and National Institute of Health Grant (MM) and the funders have no role in study design, data collection and analysis, or in manuscript preparation.
- 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.View ArticlePubMedGoogle Scholar
- Kozek WJ: Transovarially-transmitted intracellular microorganisms in adult and larval stages of Brugia malayi. J Parasitol. 1977, 63: 992-1000. 10.2307/3279832.View ArticlePubMedGoogle Scholar
- Kozek W, Rao R: The discovery of Wolbachia in arthropods and nematodes- a historical prespective. Wolbachia: a bug's life in another bug. Edited by: AH, RU R. 2007, New York: Karger, 5:Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- Taylor MJ, Hoerauf A: Wolbachia bacteria of filarial nematodes. Parasitol Today. 1999, 15: 437-442. 10.1016/S0169-4758(99)01533-1.View ArticlePubMedGoogle Scholar
- Rao R, Weil GJ: In vitro effects of antibiotics on Brugia malayi worm survival and reproduction. J Parasitol. 2002, 88: 605-611.View ArticlePubMedGoogle Scholar
- Chirgwin SR, Coleman SU, Porthouse KH, Nowling JM, Punkosdy GA, Klei TR: Removal of Wolbachia from Brugia pahangi is closely linked to worm death and fecundity but does not result in altered lymphatic lesion formation in Mongolian gerbils (Meriones unguiculatus). Infect Immun. 2003, 71: 6986-6994. 10.1128/IAI.71.12.6986-6994.2003.PubMed CentralView ArticlePubMedGoogle Scholar
- Fenn K, Blaxter M: Wolbachia genomes: revealing the biology of parasitism and mutualism. Trends Parasitol. 2006, 22: 60-65. 10.1016/j.pt.2005.12.012.View ArticlePubMedGoogle Scholar
- Zientz E, Dandekar T, Gross R: Metabolic interdependence of obligate intracellular bacteria and their insect hosts. Microbiol Mol Biol Rev. 2004, 68: 745-770. 10.1128/MMBR.68.4.745-770.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Brownlie JC, Cass BN, Riegler M, Witsenburg JJ, Iturbe-Ormaetxe I, McGraw EA, O'Neill SL: Evidence for metabolic provisioning by a common invertebrate endosymbiont, Wolbachia pipientis, during periods of nutritional stress. PLoS Pathog. 2009, 5: e1000368-10.1371/journal.ppat.1000368.PubMed CentralView ArticlePubMedGoogle Scholar
- Taylor MJ, Bandi C, Hoerauf A: Wolbachia bacterial endosymbionts of filarial nematodes. Adv Parasitol. 2005, 60: 245-284.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Townson S, Hutton D, Siemienska J, Hollick L, Scanlon T, Tagboto SK, Taylor MJ: Antibiotics and Wolbachia in filarial nematodes: antifilarial activity of rifampicin, oxytetracycline and chloramphenicol against Onchocerca gutturosa, Onchocerca lienalis and Brugia pahangi. Ann Trop Med Parasitol. 2000, 94: 801-816. 10.1080/00034980020027988.View ArticlePubMedGoogle Scholar
- Genchi C, Sacchi L, Bandi C, Venco L: Preliminary results on the effect of tetracycline on the embryogenesis and symbiotic bacteria (Wolbachia) of Dirofilaria immitis. An update and discussion. Parassitologia. 1998, 40: 247-249.PubMedGoogle Scholar
- Chirgwin SR, Nowling JM, Coleman SU, Klei TR: Brugia pahangi and Wolbachia: the kinetics of bacteria elimination, worm viability, and host responses following tetracycline treatment. Exp Parasitol. 2003, 103: 16-26. 10.1016/S0014-4894(03)00063-8.View ArticlePubMedGoogle Scholar
- Hoerauf A, Adjei O, Buttner DW: Antibiotics for the treatment of onchocerciasis and other filarial infections. Curr Opin Investig Drugs. 2002, 3: 533-537.PubMedGoogle Scholar
- Langworthy NG, Renz A, Mackenstedt U, Henkle-Duhrsen K, de Bronsvoort MB, Tanya VN, Donnelly MJ, Trees AJ: Macrofilaricidal activity of tetracycline against the filarial nematode Onchocerca ochengi: elimination of Wolbachia precedes worm death and suggests a dependent relationship. Proc Biol Sci. 2000, 267: 1063-1069. 10.1098/rspb.2000.1110.PubMed CentralView ArticlePubMedGoogle Scholar
- Hoerauf A, Volkmann L, Nissen-Paehle K, Schmetz C, Autenrieth I, Buttner DW, Fleischer B: Targeting of Wolbachia endobacteria in Litomosoides sigmodontis: comparison of tetracyclines with chloramphenicol, macrolides and ciprofloxacin. Trop Med Int Health. 2000, 5: 275-279. 10.1046/j.1365-3156.2000.00544.x.View ArticleGoogle Scholar
- Arumugam S, Pfarr KM, Hoerauf A: Infection of the intermediate mite host with Wolbachia-depleted Litomosoides sigmodontis microfilariae: impaired L1 to L3 development and subsequent sex-ratio distortion in adult worms. Int J Parasitol. 2008, 38: 981-987. 10.1016/j.ijpara.2007.12.006.View ArticlePubMedGoogle Scholar
- Casiraghi M, McCall JW, Simoncini L, Kramer LH, Sacchi L, Genchi C, Werren JH, Bandi C: Tetracycline treatment and sex-ratio distortion: a role for Wolbachia in the moulting of filarial nematodes?. Int J Parasitol. 2002, 32: 1457-1468. 10.1016/S0020-7519(02)00158-3.View ArticlePubMedGoogle Scholar
- Bosshardt SC, McCall JW, Coleman SU, Jones KL, Petit TA, Klei TR: Prophylactic activity of tetracycline against Brugia pahangi infection in jirds (Meriones unguiculatus). J Parasitol. 1993, 79: 775-777. 10.2307/3283620.View ArticlePubMedGoogle Scholar
- Hoerauf A, Mand S, Adjei O, Fleischer B, Buttner DW: Depletion of Wolbachia endobacteria in Onchocerca volvulus by doxycycline and microfilaridermia after ivermectin treatment. Lancet. 2001, 357: 1415-1416. 10.1016/S0140-6736(00)04581-5.View ArticlePubMedGoogle Scholar
- Hoerauf A, Mand S, Fischer K, Kruppa T, Marfo-Debrekyei Y, Debrah AY, Pfarr KM, Adjei O, Buttner DW: Doxycycline as a novel strategy against bancroftian filariasis-depletion of Wolbachia endosymbionts from Wuchereria bancrofti and stop of microfilaria production. Med Microbiol Immunol. 2003, 192: 211-216. 10.1007/s00430-002-0174-6.View ArticlePubMedGoogle Scholar
- Taylor MJ, Makunde WH, McGarry HF, Turner JD, Mand S, Hoerauf A: Macrofilaricidal activity after doxycycline treatment of Wuchereria bancrofti: a double-blind, randomised placebo-controlled trial. Lancet. 2005, 365: 2116-2121. 10.1016/S0140-6736(05)66591-9.View ArticlePubMedGoogle Scholar
- Wanji S, Tendongfor N, Nji T, Esum M, Che JN, Nkwescheu A, Alassa F, Kamnang G, Enyong PA, Taylor MJ, Hoerauf A, Taylor DW: Community-directed delivery of doxycycline for the treatment of onchocerciasis in areas of co-endemicity with loiasis in Cameroon. Parasit Vectors. 2009, 2: 39-10.1186/1756-3305-2-39.PubMed CentralView ArticlePubMedGoogle Scholar
- Turner JD, Tendongfor N, Esum M, Johnston KL, Langley RS, Ford L, Faragher B, Specht S, Mand S, Hoerauf A, Enyong P, Wanji S, Taylor MJ: Macrofilaricidal activity after doxycycline only treatment of Onchocerca volvulus in an area of Loa loa co-endemicity: a randomized controlled trial. PLoS Negl Trop Dis. 2010, 4: e660-10.1371/journal.pntd.0000660.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Brouqui P, Fournier PE, Raoult D: Doxycycline and eradication of microfilaremia in patients with loiasis. Emerg Infect Dis. 2001, 7: 604-605. 10.3201/eid0703.010347.PubMed CentralView ArticlePubMedGoogle Scholar
- Foster J, Ganatra M, Kamal I, Ware J, Makarova K, Ivanova N, Bhattacharyya A, Kapatral V, Kumar S, Posfai J, Vincze T, Ingram J, Moran L, Lapidus A, Omelchenko M, Kyrpides N, Ghedin E, Wang S, Goltsman E, Joukov V, Ostrovskaya O, Tsukerman K, Mazur M, Comb D, Koonin E, Slatko B: The Wolbachia genome of Brugia malayi: endosymbiont evolution within a human pathogenic nematode. PLoS Biol. 2005, 3: e121-10.1371/journal.pbio.0030121.PubMed CentralView ArticlePubMedGoogle Scholar
- Pfarr K, Hoerauf A: The annotated genome of Wolbachia from the filarial nematode Brugia malayi: what it means for progress in antifilarial medicine. PLoS Med. 2005, 2: e110-10.1371/journal.pmed.0020110.PubMed CentralView ArticlePubMedGoogle Scholar
- Kucers A, Crowe SM, Grayson ML, Hay JF: The use of antibiotics. A clinical review of antibacterial, antifungal and antiviral drugs. 1997, Boston: Butterworth & Heinemann, Oxford, U.K.Google Scholar
- Chopra I, Roberts M: Tetracycline antibiotics: mode of action, applications, molecular biology, and epidemiology of bacterial resistance. Microbiol Mol Biol Rev. 2001, 65: 232-260. 10.1128/MMBR.65.2.232-260.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Goodman CD, Su V, McFadden GI: The effects of anti-bacterials on the malaria parasite Plasmodium falciparum. Mol Biochem Parasitol. 2007, 152: 181-191. 10.1016/j.molbiopara.2007.01.005.View ArticlePubMedGoogle Scholar
- Kiatfuengfoo R, Suthiphongchai T, Prapunwattana P, Yuthavong Y: Mitochondria as the site of action of tetracycline on Plasmodium falciparum. Mol Biochem Parasitol. 1989, 34: 109-115. 10.1016/0166-6851(89)90002-9.View ArticlePubMedGoogle Scholar
- Heider U, Blaxter M, Hoerauf A, Pfarr KM: Differential display of genes expressed in the filarial nematode Litomosoides sigmodontis reveals a putative phosphate permease up-regulated after depletion of Wolbachia endobacteria. Int J Med Microbiol. 2006, 296: 287-299. 10.1016/j.ijmm.2005.12.019.View ArticlePubMedGoogle Scholar
- Ghedin E, Wang S, Spiro D, Caler E, Zhao Q, Crabtree J, Allen JE, Delcher AL, Guiliano DB, Miranda-Saavedra D, Angiuoli SV, Creasy T, Amedeo P, Haas B, El-Sayed NM, Wortman JR, Feldblyum T, Tallon L, Schatz M, Shumway M, Koo H, Salzberg SL, Schobel S, Pertea M, Pop M, White O, Barton GJ, Carlow CK, Crawford MJ, Daub J: Draft genome of the filarial nematode parasite Brugia malayi. Science. 2007, 317: 1756-1760. 10.1126/science.1145406.PubMed CentralView ArticlePubMedGoogle Scholar
- Kumar S, Chaudhary K, Foster JM, Novelli JF, Zhang Y, Wang S, Spiro D, Ghedin E, Carlow CK: Mining predicted essential genes of Brugia malayi for nematode drug targets. PLoS ONE. 2007, 2: e1189-10.1371/journal.pone.0001189.PubMed CentralView ArticlePubMedGoogle Scholar
- Ash LR: Chronic Brugia pahangi and Brugia malayi infections in Meriones unguiculatus. J Parasitol. 1973, 59: 442-447. 10.2307/3278769.View ArticlePubMedGoogle Scholar
- Rao RU, Huang Y, Fischer K, Fischer PU, Weil GJ: Brugia malayi: Effects of nitazoxanide and tizoxanide on adult worms and microfilariae of filarial nematodes. Exp Parasitol. 2009, 121: 38-45. 10.1016/j.exppara.2008.09.020.View ArticlePubMedGoogle Scholar
- Li BW, Rush AC, Mitreva M, Yin Y, Spiro D, Ghedin E, Weil GJ: Transcriptomes and pathways associated with infectivity, survival and immunogenicity in Brugia malayi L3. BMC Genomics. 2009, 10: 267-10.1186/1471-2164-10-267.PubMed CentralView ArticlePubMedGoogle Scholar
- Li BW, Rush AC, Jiang DJ, Mitreva M, Abubucker S, Weil GJ: Gender-associated genes in filarial nematodes are important for reproduction and potential intervention targets. PLoS Negl Trop Dis. 2011, 5: e947-10.1371/journal.pntd.0000947.PubMed CentralView ArticlePubMedGoogle Scholar
- Weiss GB, Wilson GN, Steggles AW, Anderson WF: Importance of full size complementary DNA in nucleic acid hybridization. J Biol Chem. 1976, 251: 3425-3431.PubMedGoogle Scholar
- Blaxter M, Daub J, Guiliano D, Parkinson J, Whitton C: The Brugia malayi genome project: expressed sequence tags and gene discovery. Trans R Soc Trop Med Hyg. 2002, 96: 7-17. 10.1016/S0035-9203(02)90224-5.View ArticlePubMedGoogle Scholar
- Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008, 3: 1101-1108. 10.1038/nprot.2008.73.View ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M: KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010, 38: D355-360. 10.1093/nar/gkp896.PubMed CentralView ArticlePubMedGoogle Scholar
- Mulder NJ, Apweiler R, Attwood TK, Bairoch A: New developments in the InterPro database. Nucleic Acids Res. 2007, 35: D224-228. 10.1093/nar/gkl841.PubMed CentralView ArticlePubMedGoogle Scholar
- Prufer K, Muetzel B, Do HH, Weiss G, Khaitovich P, Rahm E, Paabo S, Lachmann M, Enard W: FUNC: a package for detecting significant associations between gene sets and ontological annotations. BMC Bioinformatics. 2007, 8: 41-10.1186/1471-2105-8-41.PubMed CentralView ArticlePubMedGoogle Scholar
- Ghedin E, Hailemariam T, DePasse JV, Zhang X, Oksov Y, Unnasch TR, Lustigman S: Brugia malayi gene expression in response to the targeting of the Wolbachia endosymbiont by tetracycline treatment. PLoS Negl Trop Dis. 2009, 3: e525-10.1371/journal.pntd.0000525.PubMed CentralView ArticlePubMedGoogle Scholar
- Strübing U, Lucius R, Hoerauf A, Pfarr KM: Mitochondrial genes for heme-dependent respiratory chain complexes are up-regulated after depletion of Wolbachia from filarial nematodes. Int J Parasitol. 2010, 1193-1202.Google Scholar
- Bandi C, McCall JW, Genchi C, Corona S, Venco L, Sacchi L: Effects of tetracycline on the filarial worms Brugia pahangi and Dirofilaria immitis and their bacterial endosymbionts Wolbachia. Int J Parasitol. 1999, 29: 357-364. 10.1016/S0020-7519(98)00200-8.View ArticlePubMedGoogle Scholar
- Hoerauf A, Mand S, Volkmann L, Buttner M, Marfo-Debrekyei Y, Taylor M, Adjei O, Buttner DW: Doxycycline in the treatment of human onchocerciasis: Kinetics of Wolbachia endobacteria reduction and of inhibition of embryogenesis in female Onchocerca worms. Microbes Infect. 2003, 5: 261-273. 10.1016/S1286-4579(03)00026-1.View ArticlePubMedGoogle Scholar
- Nam DK, Lee S, Zhou G, Cao X, Wang C, Clark T, Chen J, Rowley JD, Wang SM: Oligo(dT) primer generates a high frequency of truncated cDNAs through internal poly(A) priming during reverse transcription. Proc Natl Acad Sci USA. 2002, 99: 6152-6156. 10.1073/pnas.092140899.PubMed CentralView ArticlePubMedGoogle Scholar
- Richmond CS, Glasner JD, Mau R, Jin H, Blattner FR: Genome-wide expression profiling in Escherichia coli K-12. Nucleic Acids Res. 1999, 27: 3821-3835. 10.1093/nar/27.19.3821.PubMed CentralView ArticlePubMedGoogle Scholar
- Wilson M, DeRisi J, Kristensen HH, Imboden P, Rane S, Brown PO, Schoolnik GK: Exploring drug-induced alterations in gene expression in Mycobacterium tuberculosis by microarray hybridization. Proc Natl Acad Sci USA. 1999, 96: 12833-12838. 10.1073/pnas.96.22.12833.PubMed CentralView ArticlePubMedGoogle Scholar
- Rao R, Huang Y, Crosby S, Mitreva M, Yin Y, Weil GJ: Long-term doxycycline treatment affects Wolbachia and parasite gene expression in adult female Brugia malayi [abstract]. Amer J Trop Med Hyg. 2007, 77: 779-Google Scholar
- Holman AG, Davis PJ, Foster JM, Carlow CK, Kumar S: Computational prediction of essential genes in an unculturable endosymbiotic bacterium, Wolbachia of Brugia malayi. BMC Microbiol. 2009, 9: 243-10.1186/1471-2180-9-243.PubMed CentralView ArticlePubMedGoogle Scholar
- Maizels RM, Gomez-Escobar N, Gregory WF, Murray J, Zang X: Immune evasion genes from filarial nematodes. Int J Parasitol. 2001, 31: 889-898. 10.1016/S0020-7519(01)00213-2.View ArticlePubMedGoogle Scholar
- Emes R, Thompson F, Moore J, Zang X, Devaney E: Cloning and characterisation of mmc-1, a microfilarial-specific gene, from Brugia pahangi. Int J Parasitol. 2002, 32: 415-424. 10.1016/S0020-7519(02)00003-6.View ArticlePubMedGoogle Scholar
- Wu B, Novelli J, Foster J, Vaisvila R, Conway L, Ingram J, Ganatra M, Rao AU, Hamza I, Slatko B: The heme biosynthetic pathway of the obligate Wolbachia endosymbiont of Brugia malayi as a potential anti-filarial drug target. PLoS Negl Trop Dis. 2009, 3: e475-10.1371/journal.pntd.0000475.PubMed CentralView ArticlePubMedGoogle Scholar
- Dedeine F, Vavre F, Fleury F, Loppin B, Hochberg ME, Bouletreau M: Removing symbiotic Wolbachia bacteria specifically inhibits oogenesis in a parasitic wasp. Proc Natl Acad Sci USA. 2001, 98: 6247-6252. 10.1073/pnas.101304298.PubMed CentralView ArticlePubMedGoogle Scholar
- Xi Z, Gavotte L, Xie Y, Dobson SL: Genome-wide analysis of the interaction between the endosymbiotic bacterium Wolbachia and its Drosophila host. BMC Genomics. 2008, 9: 1-10.1186/1471-2164-9-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Hernandez-Bello R, Bermudez-Cruz RM, Fonseca-Linan R, Garcia-Reyna P, Le Guerhier F, Boireau P, Ortega-Pierres G: Identification, molecular characterisation and differential expression of caveolin-1 in Trichinella spiralis maturing oocytes and embryos. Int J Parasitol. 2008, 38: 191-202. 10.1016/j.ijpara.2007.07.009.View ArticlePubMedGoogle Scholar
- Michalski ML, Monsey JD, Cistola DP, Weil GJ: An embryo-associated fatty acid-binding protein in the filarial nematode Brugia malayi. Mol Biochem Parasitol. 2002, 124: 1-10. 10.1016/S0166-6851(02)00081-6.View ArticlePubMedGoogle Scholar
- Merriweather A, Guenzler V, Brenner M, Unnasch TR: Characterization and expression of enzymatically active recombinant filarial prolyl 4-hydroxylase. Mol Biochem Parasitol. 2001, 116: 185-197. 10.1016/S0166-6851(01)00317-6.View ArticlePubMedGoogle Scholar
- Britton C, Murray L: Cathepsin L protease (CPL-1) is essential for yolk processing during embryogenesis in Caenorhabditis elegans. J Cell Sci. 2004, 117: 5133-5143. 10.1242/jcs.01387.View ArticlePubMedGoogle Scholar
- Ford L, Zhang J, Liu J, Hashmi S, Fuhrman JA, Oksov Y, Lustigman S: Functional analysis of the cathepsin-like cysteine protease genes in adult Brugia malayi using RNA interference. PLoS Negl Trop Dis. 2009, 3: e377-10.1371/journal.pntd.0000377.PubMed CentralView ArticlePubMedGoogle Scholar
- Dangi A, Vedi S, Nag JK, Paithankar S, Singh MP, Kar SK, Dube A, Misra-Bhattacharya S: Tetracycline treatment targeting Wolbachia affects expression of an array of proteins in Brugia malayi parasite. Proteomics. 2009, 9: 4192-4208. 10.1002/pmic.200800324.View ArticlePubMedGoogle Scholar
- Clarke ND, Berg JM: Zinc fingers in Caenorhabditis elegans: finding families and probing pathways. Science. 1998, 282: 2018-2022.View ArticlePubMedGoogle Scholar
- Maeda I, Kohara Y, Yamamoto M, Sugimoto A: Large-scale analysis of gene function in Caenorhabditis elegans by high-throughput RNAi. Curr Biol. 2001, 11: 171-176. 10.1016/S0960-9822(01)00052-5.View ArticlePubMedGoogle Scholar
- Kamath RS, Ahringer J: Genome-wide RNAi screening in Caenorhabditis elegans. Methods. 2003, 30: 313-321. 10.1016/S1046-2023(03)00050-1.View ArticlePubMedGoogle Scholar
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.