Skip to main content

Advertisement

We’d like to understand how you use our websites in order to improve them. Register your interest.
Download PDF

Transcriptomes and pathways associated with infectivity, survival and immunogenicity in Brugia malayi L3

BMC Genomics volume 10, Article number: 267 (2009) Cite this article

Abstract

Background

Filarial nematode parasites cause serious diseases such as elephantiasis and river blindness in humans, and heartworm infections in dogs. Third stage filarial larvae (L3) are a critical stage in the life cycle of filarial parasites, because this is the stage that is transmitted by arthropod vectors to initiate infections in mammals. Improved understanding of molecular mechanisms associated with this transition may provide important leads for development of new therapies and vaccines to prevent filarial infections. This study explores changes in gene expression associated with the transition of Brugia malayi third stage larvae (BmL3) from mosquitoes into mammalian hosts and how these changes are affected by radiation. Radiation effects are especially interesting because irradiated L3 induce partial immunity to filarial infections. The underlying molecular mechanisms responsible for the efficacy of such vaccines are unkown.

Results

Expression profiles were obtained using a new filarial microarray with 18, 104 64-mer elements. 771 genes were identified as differentially expressed in two-way comparative analyses of the three L3 types. 353 genes were up-regulated in mosquito L3 (L3i) relative to cultured L3 (L3c). These genes are important for establishment of filarial infections in mammalian hosts. Other genes were up-regulated in L3c relative to L3i (234) or irradiated L3 (L3ir) (22). These culture-induced transcripts include key molecules required for growth and development. 165 genes were up-regulated in L3ir relative to L3c; these genes encode highly immunogenic proteins and proteins involved in radiation repair. L3ir and L3i have similar transcription profiles for genes that encode highly immunogenic proteins, antioxidants and cuticle components.

Conclusion

Changes in gene expression that normally occur during culture under conditions that support L3 development and molting are prevented or delayed by radiation. This may explain the enhanced immunogenicity of L3ir. Gene Ontology and KEGG analyses revealed altered pathways between L3 types. Energy and "immune pathways" are up-regulated and may be needed for L3i invasion and survival, while growth and development are priorities for L3c. This study has improved our understanding of molecules involved in parasite invasion and immune evasion, potential targets of protective immunity, and molecules required for parasite growth and development.

Background

Filarial nematodes cause a number of serious diseases in humans (e.g., lymphatic filariasis and onchocerciasis) and animals (e.g., canine heartworm disease). Filarial parasites have complex life cycles that share the following general features: Arthropod vectors (ticks, mites, mosquitoes, black flies) ingest microfilariae (Mf) that circulate in the blood or live in the skin of definitive vertebrate hosts. Ingested Mf penetrate the arthropod midgut and develop over a period of 10 to 14 days to become third-stage infective larvae (L3). Filarial L3, developmentally arrested in arthropod vectors, resume development after they are transferred to a permissive vertebrate host. The L3 is the critical life stage for infection, as it is the transitional stage between the vector and the mammalian host. Molecular mechanisms associated with this transition are poorly understood; they may provide important clues for new therapies and vaccines to prevent filarial infections.

Prior studies have used different methods to investigate adaptations that occur in filarial L3 after they enter mammalian hosts (or in response to culture conditions that mimic the mammalian environment). Early studies used metabolic labeling and surface labeling techniques with SDS-PAGE to detect altered protein patterns [1–3]. Molecular biological approaches have included differential hybridization [2] and analysis of expressed sequence tags (ESTs) from cDNA libraries from pre- and post-infective L3 [4, 5]. These studies provided interesting, preliminary information regarding adaptations associated with the transition of L3 from arthropod to mammalian hosts.

In addition to being the infective stage for vertebrate hosts, L3 are also important targets of protective immunity [6–8]. A number of studies have shown that irradiated L3 (L3ir) are more effective for inducing immunity to infection than L3 freshly isolated from vectors (L3i) [9–11]. Potential targets of such protective immunity have been identified by screening crude antigen preparations and expression libraries with sera from animals vaccinated with irradiated L3s [11, 12]. We have previously used quantitative real-time PCR (qRT-PCR) analysis to identify a limited number of genes that were differentially expressed by cultured or irradiated L3 [5]. The present study builds on this prior work to more broadly describe changes in L3 gene expression associated with the transition from arthropod to mammalian host conditions (presumably related to invasion, growth, and immune evasion) and effects of radiation on gene expression (which may provide insight into targets of protective immunity). While prior studies focused on ESTs from cDNA libraries, the current study benefited from recent advances in filarial genomics [13] and microarray technology to provide more complete expression profiles. We have also compared changes in filarial L3 gene expression associated with the vector-vertebrate transition with changes associated with the soil-vertebrate transition in Ancylostoma caninum (canine hookworm) and with dauer exit in Caenorhabditis elegans.

Results and discussion

General properties of the B. malayi Version 2 array (BmV2 array)

The oligonucleotide elements on the BmV2 array represent ~86% of currently annotated B. malayi gene models defined by Ghedin et al [13]; 11,975 of 15,412 B. malayi sequences correspond to 9,798 unique annotated gene models, and the remaining 3,437 sequences that do not match annotated gene models represent genes that have not yet been sequenced or annotated in the Brugia genomic DNA. The BmV2array elements included all genes in the Wolbachia genome and ESTs from O. volvulus and W. bancrofti with little homology to known B. malayi sequences. A homology search of the 18,104 sequences from which the oligos were derived vs. the NCBI non-redundancy (NR) database revealed that ~69% (12, 423/18,104) had significant homology to publicly available known or unknown proteins from other species (e-05) [See additional file 1]. Of 15,412 Brugia sequences, 9,824 sequences had best matches to 8,038 unique C. elegans genes with RNAi information; 37% of these genes had RNAi phenotypes [See additional file 2].

Functional assignments for elements on the microarray were performed using InterPro, KEGG, and GO analysis. The complete results are reported in supplemental materials [protein domain matches in additional file 3; pathways in additional file 4]. The top 20 most abundant protein domain matches in the oligoarray are presented in additional file 5. Based on InterPro protein domain matches, 48% of the products of genes on the array mapped to one or more gene ontology terms [14, 15]. These GO associations are available online through the Amigo viewer http://nematode.net/cgi-bin/amigo/go_brugia_malayi_microarray_v2/go.cgi?session_id=94439521169750740719

Overview of differences in gene expression by L3 type

A complete list of elements for the array with oligonucleotide sequences and hybridization data are available online through the link: http://nematode.net/Microarray/index.php. The expression profile analysis included elements with hybridization signals above background for at least one condition in pair-wise comparisons (L3i versus L3c and L3ir versus L3c). A total of 774 genes exhibited significant differential expression in pair-wise comparisons (≥ 2 fold differences with P-values < 0.01) [See additional file 6]. 353 of these genes (~46%) were up-regulated in L3i relative to L3c (the "L3i gene set"); 266 (~33%) were up-regulated in L3c (the "L3c gene set", with 244 genes up relative to L3i and 22 up relative to L3ir); 165 (~21.3%) were L3ir-up-regulated relative to L3c (the "L3ir gene set') (Figure 1). The L3i gene set contains genes required for L3 survival in mosquitoes and for infection of mammalian hosts. The L3c gene set represents genes induced by culture conditions that mimic the mammalian environment, and the L3ir gene set represents genes induced by irradiation. Interestingly, some genes exhibited up-regulated expression in more than one condition (Figure 1).

Figure 1
figure1

Venn diagram show distribution of up-regulated genes in three L3 type. For example, out of 353 L3i up-regulated genes, 75 genes are up-regulated in both L3i and L3ir relative to L3c.

Full size image

We were surprised to see that the L3i list was the largest of the three L3 gene sets. However, this was consistent with the finding that more Ancylostoma canium L3 (AcL3) genes had reduced transcription after serum stimulation (equivalent to our L3c) relative to non-serum stimulated AcL3 (equivalent to our L3i) [16, 17]. This result suggests that the L3i is not simply a quiescent or resting stage. L3i are actively preparing for their transition into the mammalian host by expressing proteins that are needed for processes like invasion and immune evasion and proteins for survival (stress resistance and antioxidants).

As expected, irradiation repressed gene expression globally; the L3ir gene set was the smallest of the three L3 gene sets. Of 165 L3ir-up-regulated genes, 75 (45%) genes were also up-regulated in L3i relative to L3c (Figure 1). Most of these genes (54/75) were expressed at higher levels in L3i than in L3ir, and only 4 (intermediate filament, collagen family, and Gly-rich proteins) had higher expression in L3ir; 17 were expressed at similar levels in both conditions [See additional file 6].

We previously reported partial conservation of L3i expression patterns in L3ir based on an EST analysis with qRT-PCR confirmation data [5]. The current analysis has confirmed and extended this observation. Several collagen genes and other so-called immunogenic genes (Gln-rich, spx-1), were over-expressed in L3ir relative to L3c [See additional file 6]. Persistent expression of important targets of protective immunity may explain why L3ir are better inducers of anti-L3 immunity than L3i, which are likely to change gene expression to resemble L3c with reduced expression of these proteins shortly after infection.

The taxonomic distribution of proteins encoded by the up-regulated gene models in each condition are shown in Table 1. More proteins encoded by B. malayi gene models in the L3i and L3ir gene sets are nematode-specific as defined by Ghedin et al [13] compared to those in the L3c gene set. Nematode-specific gene products are attractive candidate targets for new drugs and vaccines because they lack homology to mammalian proteins.

Table 1 Taxonomic distribution of known proteins encoded by up-regulated gene models in the L3 gene set
Full size table

RNAi information for C. elegans homologues of genes associated with different BmL3 types are listed in additional file 7. RNAi resulted in phenotypes in 10%–14% of all C. elegans predicted genes [18–20] and in 27% of genes with evidence of expression [21]. In contrast, a much higher percentage of C. elegans genes with strong homology to genes associated with different BmL3 types had RNAi phenotypes (Table 2). C. elegans genes homologous to genes up-regulated in BmL3c were more likely to have RNAi phenotypes (60%) than those with homologues upregulated in L3i or L3ir. This was consistent with results reported for Strongyloides stercoralis [22]. The most common RNAi phenotypes observed for C. elegans genes with homology to genes in the BmL3 gene sets were developmental delay (slow growth, larval arrest and lethal), embryonic lethal, and sterile [See additional file 7].

Table 2 RNAi phenotypes for differentially expressed filarial genes with close C. elegans homolog
Full size table

Distinct functions of known genes associated with different L3 gene sets

The L3c gene set had a higher percentage of novel genes (42%) than the L3ir (27%) and L3i (28%) gene sets (Table 3). This is consistent with results from our study of ESTs from BmL3 cDNA libraries [5] and a recent study of transcriptional changes in AcL3 [17]. These novel genes may be filaria-specific, because they have no similarity to known genes in public databases. Additional work will be needed to determine the function of these genes. In the meantime, we can learn a lot from analysis of differentially expressed known genes. We grouped the most abundant up-regulated known genes in each L3 gene set by biological function (Table 4, 5, 6), and this information is discussed below.

Table 3 Distribution of BLAST matches in the L3 gene set
Full size table
Table 4 Major known genes by functional class in the L3i gene set
Full size table
Table 5 Major known genes by functional class in the L3ir gene set
Full size table
Table 6 Major known genes by functional class in the L3c gene set
Full size table

The L3i gene set contains genes associated with immune evasion, stress resistance and infectivity

Several classes of known genes were highly expressed in the L3i gene set as shown in Table 4. Immunomodulatory genes are important for invasion and establishment of infection in mammalian hosts [23], and prior studies have suggested that L3i modulate the innate immune system [6]. Microarray experiments showed that live BmL3 induced functional alterations of epidermal Langerhans cells (LC) and diminished their ability to present antigens to T-cells. Thus, L3 may actively suppress host immune responses to the parasite immediately after invasion [24]. Our findings suggest a molecular basis for such immune evasion, because a number of L3i up-regulated genes encode immune-modulators. These molecules may suppress local immunity to promote parasite survival just after invasion. For example, filarial ALT proteins encoded by the Bm-alt gene family modulate cytokine-induced signaling in immune cells and may render incoming L3 more resistant to IFN-γ induced killing by macrophages, as shown in a transgenic L. mexicana system [25]. Filarial cystatins are well-described pathogenicity factors that down-regulate T-cell proliferation and induce anti-inflammatory cytokine responses [26, 27]. The cystatin protein encoded by Bm-cpi-2 interferes with antigen presentation by inhibiting multiple cysteine protease activities including the key enzyme asparaginyl endopeptidase [28]. Bm-spn-2, a member of the B. malayi serpin family, has been reported to impair granulocyte function by inhibiting two neutrophil enzymes, cathepsin G and neutral elastase [29].

The L3i gene list includes a broad spectrum of anti-oxidant and detoxification gene families such as superoxide dismutase (Bm-sods) [30], thioredoxin peroxidase (Bm-tpx-2) [31] glutathione S-transferase (Bm-gst) [32, 33], and aldo/keto reductase family proteins [34]. Since non-feeding filarial L3 are protected by a resistant cuticle, increased transcriptional activity of these genes may be needed for detoxifying endogenous compounds such as metabolites of fatty acids [35]. Fatty acids appear to be an important energy source in L3i as shown below. These enzymes also protect L3 from reactive oxygen species (ROS) generated in energy production process (oxidative phosphorylation) [36]. In addition, L3i are preparing to cope with oxidative stress upon host entry. Prior studies have shown that filarial larvae (L1 and L3) are more susceptible to killing than adult worms via antibody-dependent cellular cytotoxicity in vitro. This killing is believed to depend on reactive oxygen and nitrogen intermediates produced by activated myeloid cells [37, 38] L3 must be able to resist oxidative stress to survive in mammalian hosts.

The L3i gene set also includes several proteases that may be important for invasion and molting. For example, filarial cathepsin L-like proteases (Bm-CPL-1,-4 and -5) are believed to be important for larval migration and molting [39–41]. L3i express a protease (BMC12240), which is analogous to a matrix metalloprotease in the canine hookworm (Ac-mtp-1) that is highly expressed in non-activated AcL3 and involved in skin penetration [17, 42]. Thus, our results support the concept that proteases involved in tissue migration are over-expressed in helminth stages that infect mammalian hosts [43].

The L3i gene set also includes genes that are believed to be involved in pathogenesis, parasitism and stress resistance. For example, the highly expressed Bm-val-1 encodes a venom allergen-like protein that is related to genes in A. caninum that encode Ancylostoma secreted-proteins (ASPs) [44]. ASPs belong to the pathogenesis-related proteins superfamily (PRPs) [45]; they are associated with the transition to parasitism in hookworms and highly expressed in non-activated AcL3 [16, 17, 46]. It is interesting that the O. volvolus homologue of Bm-val-1 (Ov-ASP-1) has been reported to induce anti-L3 immunity in an animal model [47].

The L3i gene set also includes a unique small heat shock protein (sHSP) of B. malayi that was represented by two ESTs on the array (BMC02059 and TC3642). This gene was highly expressed in L3i relative to L3c, and its over-expression was confirmed by qRT-PCR (data not shown). This sHSP is developmentally regulated and inducible in B. malayi [48]. Small heat shock proteins are among the most highly expressed dauer-specific transcripts in C. elegans, and they have a demonstrated role in longevity and stress resistance [49–51].

The last two categories of proteins highly expressed in L3i were surface or cuticular proteins and highly immunogenic proteins. Genes encoding cuticular collagens accounted for 19% of L3i up-regulated known genes (47/249). Their dominant presence is understandable, as the nematode cuticle protects the worms from environmental stress. Cuticular collagens are also needed for the development of the next stage of the parasite (L4) after invasion. C. elegans homologs of these genes have RNAi phenotypes of developmental delay such as larval arrest (Lva), larval lethal (Lvl) (BMC0572), and slow growth (Gro) (BMC12441) [See additional file 6]. Collagen genes are over-expressed just before molting in C. elegans [51]. The L3i gene set also contains a highly expressed ground-like family member (grl-4) [52] (4083.m00057). Grl-4 is believed to function during the dauer-to-L4 molt in C. elegans [51]. It would be interesting to test whether this gene is also required for molting in filarial worms. The highly immunogenic proteins, such as Bm-VAL, intermediate filament, and microfilaria surface-associated proteins may be targets of protective immunity [28, 53]. Filarial SPX-1 and Gln-rich protein have been used as diagnostic antigens [54, 55]. These proteins have also been suggested as targets of protective immunity to Brugia and Ascaris [56, 57]. Brugia larval allergen (gene 14323.m00082) was differentially recognized by IgE antibody from putatively immune endemic normal humans; this protein induces specific IgE antibody responses in animals (authors' unpublished data).

The L3ir gene set contains irradiation responsive (IR) genes and highly immunogenic genes that are also up-regulated in L3i

The L3ir gene set includes genes involved in DNA repair and cuticle regeneration after radiation. Genes reported to be IR in other systems were also up-regulated in L3ir as shown in Table 5. These included cytidine deaminase, uridine phosphorylase, aquaporin and finger proteins [58–60]. Notably, various collagen gene families were also over-represented, and these accounted for more than 30% of the total and 50% of L3ir known genes, respectively. Collagen is a principle component of the nematode cuticle that is secreted by the underlying hypodermis, and collagen synthesis is often increased after radiation injury in animals [61, 62]. The nematode cuticle protects the parasite from environmental insults, and it may be an important target of immune responses [63, 64]. The collagen gene Ovcol-1 of O. volvulus, a homologue of the B. malayi collagen gene Bmcol-2 (BMC12441, TC2798) up-regulated in L3ir, is preferentially recognized by immunoglobulin G3 (IgG3) from putatively immune individuals [65]. Cuticular collagens of nematodes are very different from vertebrate collagens, and they are potential targets of protective immune responses [66]. In addition to collagens, the L3ir gene set contains highly immunogenic secretory and membrane proteins; most of these are also present in the L3i gene set. It is interesting that these genes were more highly expressed in L3ir (which induce protective immunity) than in L3c (which do not). It is possible that the enhanced immunogenicity of L3ir is partially due to their persistent high level expression of L3i proteins that are normally down-regulated shortly after L3 enter the mammalian host. This could prolong and enhance priming of the host immune system to critical protective antigens.

The L3c gene set contains genes that are involved in growth and molting

In contrast to L3i and L3ir, most of the known L3c up-regulated genes were involved in protein expression, reflecting increased protein synthesis during culture (Table 6). Genes encoding various ribosomal proteins accounted for 25% of known L3c genes (38/154). It is interesting that some of these genes have C. elegans homologues involved in molting as shown in Table 7[67]. C. elegans homologues of those genes have severe embryonic lethal and developmental delay phenotypes such as sterile and larval arrest (13787.m00042 and BMC00239) [See additional file 7]. Our data support the idea that molting requires a burst of biosynthetic activity, presumably to make components for the new cuticle and for growth [67, 68]. Other over-represented categories were genes involved in cell cycle and phagocytosis and genes involved in fiber growth. Many of the L3c up-regulated genes are homologues of genes that are also up-regulated in activated (post infection) hookworm L3. These include genes that encode heat shock proteins, tubulin, and rRNA [16, 17]. In addition, two oligos matched to small heat shock proteins (sHSPs) (BMC03771 and 15334.m00009) were up-regulated, and this is consistent with previous report [2].

Table 7 L3c genes for protein synthesis with C. elegans homologs identified in a molting screen*
Full size table

Overview of pathways and GO analysis

KEGG and InterPro analyses were performed to functionally classify and categorize genes up-regulated in each L3 type. The complete KEGG pathways identified for each gene set are listed in additional file 8. Significant protein domain matches for each gene set (P < 0.05) are shown in additional file 9. The GO associations for whole L3 gene set are available online through the Amigo viewer http://nematode.net/cgi-bin/amigo/go_brugia_malayi_L3/go.cgi?search_constraint=terms&action=replace_tree&session_id=91575331171470308620. Differences in pathways and GO terms inferred by transcriptional expression profiling were observed among different L3 types. It is important to note that the RNA levels do not correlate with protein level in all cases; therefore, such differences are needed to be confirmed using appropriate methods.

Altered pathways and significant GO terms in L3i gene set

KEGG pathways that were significantly enriched in the L3i gene set are shown in Table 8. An overview of the main up-regulated pathways for energy metabolism in L3i is presented in Figure 2. Developmentally arrested and non-feeding infective L3 were expected to have low energy production. We were surprised to see that the most enriched pathway mapped in the L3i gene set was energy metabolism (oxidative phosphorylation (ko00190)). This suggests that energy production and consumption might be high in L3i. High energy consumption might be needed in part to detoxify endogenous metabolic products or host proteins. For example, the detoxification reactions catalyzed by GSTs and dehydrogenases/reductases are metabolically costly [35]. Energy is also required for invasion as filarial L3 actively migrate immediately after entering the host through subcutaneous tissues en route to lymphatic vessels [69]. However, additional studies will be needed to test these hypotheses. The microarray data show that BmL3i probably use oxidative phosphorylation and the TCA cycle to produce and store energy through aerobic cellular respiration in the form of ATP as shown in Figure 2. All of the enzymes needed for electron transport in the oxidative phosphorylation pathway are up-regulated in L3i. This explains the association of the significant GO term "mitochondria" with these enzymes [See additional file 10]. This also explains high expression of enzymes to detoxify ROS such as superoxide dismutase and catalase (Table 4) in L3i, since complexes I and III of the respiratory chain generate these ROS [36]. NADH and succinate needed for ATP production can be generated in the citric acid cycle (TCA, Ko00020). Genes encoding three TCA enzymes (malate dehydrogenase, succinate dehydrogenase, and isocitrate dehydrogenase) were over-expressed by a factor of 4 in L3i relative to L3c. Acetyl-CoA, the first molecule entering the TCA, is produced through fatty acid β-oxidation; acetyl-CoA C-acyltransferase (K00632, E2.3.1.16, fadA), the enzyme catalyzing the final step of β-oxidation, was highly expressed in L3i [See additional file 8]. This is consistent with the finding in dauer stage of C. elegans [51, 70]. This suggests that the metabolism of non-feeding BmL3 is adapted to utilize internal energy reserves, predominantly fatty acids. However, there may be the alternate methods to produce Acetyl-CoA in BmL3. For example, Acetyl-CoA might be produced through the pyruvate dehydrogenase complex, as genes encoding two enzymes in the complex were up-regulated; Acetyl-CoA also could be formed by acetate and coenzyme A, as acetyl-coenzyme A synthetase (k01895, E 6.2.1.1) was also over-expressed in L3i. We also found up-regulation of genes encoding enzymes for pyruvate metabolism (ko00620) and gluconeogenesis/glycolysis (ko00010) in L3i including phosphoenolpyruvate carboxykinase (PEPCK) and triose-phosphate isomerase [See additional file 11]. This is also consistent with results reported for dauer larvae in C. elegans [[35, 49], and [51]].

Table 8 Significant KEGG pathways in the L3i gene set
Full size table
Figure 2
figure2

Overview of the pathways involved in energy generation in BmL3i. The red arrows and numbers indicate increased levels of enzymatic activity as inferred by fold changes in gene expression. ETC: electron transport chain; Q: electron carrier. The numbers in parentheses are fold changes: 1, isocitrate dehydrogenase-EC 1.1.1.42: 14990.m07912 (3 fold); 2, succinate-ubiquinone oxidoreductase – EC 1.3.5.1: BMC09915 (3 fold)and WB-contig_1179 (4.3); 3, malate dehydrogenase – EC 1.1.1.37: BMC05093 (3 fold), BMC11994 (4 fold); 4, NADH dehydrogenase (ubiquinone): – EC1.6.5.3: there are four elements encoding the enzyme – BMC07652 (3.3 fold), TC3428 (2.9 fold), BMC04790 (4 fold) and 14062.m00071(2.6 fold); 5, Ubiquinol – cytochrome-c reductase or Cytochrome bc1 complex – EC 1.10.2.2: WB-contig_1179 (4.3); 6, Cytochrome c oxidase-EC 1.9.3.1: BMC06398 (2.7); 7, ATP synthase-EC 3.6.3.14: 14980.m02712 (3.4), 14982.m02237 (3.3) and WB-contig_234 (2); 8, Pyruvate dehydrogenase complexes-EC 1.2.4.1: BMC08144 (4.3) and WB-contig_1351 (5.4); 9, Acetyl-CoA C-acyltransferase-EC 2.3.1.16: TC3583 (6.2); 10, Acetyl-coenzyme A synthetase-EC 6.2.1.1: BMC03124 (2.6).

Full size image

GO analysis indicated that genes up-regulated in L3i were significantly involved in four biological processes (namely generation of precursor metabolites and energy GO:0006091; establishment of localization GO:0051234; proteolysis GO:0006508 and coenzyme metabolic process GO:0006732), and associated with three molecular functions (structural constituent of cuticle GO: 0042302; catalytic activity GO: 0003824; and enzyme regulator activity GO: 0030234) [See additional file 10]. This is consistent with the results of the KEGG analysis. BmL3 generate metabolites and energy by oxidative phosphorylation (GO:0006119). Coordinately, genes involved in the molecular function "electron carrier activity" (GO:0009055) are over-expressed in L3i. Most of the genes involved in the molecular function "structural constituent of cuticle" (GO: 0042302) encoded collagens (cuticular collagen 34 and 14 and 2). Interestingly, these collagens were also mapped into the biological process "phosphate transport" (GO: 0006817), and some of these genes were nematode – (BMC09088) or filarial-specific (BMC06834, TC2884). C. elegans homologs of up-regulated collagens in L3i had RNAi phenotypes such as embryonic lethal (Emb, BMC05026; BMC12506), developmental delay, and larval lethal (Lvl) (BMC05172) [See additional file 7]. These B. malayi collagens merit further study as potential vaccine targets.

Several enzyme activities such as oxidoreductase, hydrolase, and kinases were obviously increased in L3i. These molecular functions were well related to parasite survival (oxidoreductase), parasite invasion and development (hydrolase), and energy production (kinase). The increased enzyme regulator activity (GO:0030234) was mainly for enzyme inhibitor activity (GO: 0004857), especially increased activities of serine-type endopeptidase inhibitor (serpin) and a cysteine protease inhibitor (cystatin). In addition, the products of 14965.m00429 and BMC02136, up-regulated in L3i, were mitochondrial ATPase inhibitors (ATP family protein) associated with down-regulation of nucleotide metabolic processes (GO:0045980, biological process). These molecules might contribute to developmental arrest of L3i in mosquitoes through negative regulation of nucleotide metabolism. GO analysis identified the mitochondrion (GO:0005739) and extracellular region (GO:0005576) as significant cellular localizations for L3i up-regulated genes.

Altered pathways and GO terms in L3ir and L3c

The same approach was followed to determine functional assignments for the L3ir and L3c gene sets. To our surprise, L3ir may rely heavily on fatty acids for energy during in vitro culture as shown in Table 9. Two genes encoding enzymes involved in fatty acid β-oxidation, namely acetyl-CoA C-acyltransferase (TC3583, 3 fold) and carnitine O-palmitoyltransferase (BMC04165, 2.3 fold; Pub_Locus Bm1_02015, K00636, EC2.3.1.21), were over- expressed in L3ir relative to L3c [See additional file 6]. In addition, genes encoding cellular membrane proteins (BMC02613, 3 fold; 13888.m00008, 4 fold), a family of transmembrane channels called aquaporins (AQPs), were also over-expressed in L3ir relative to L3c; these genes mapped into the "pores ion channels" KEGG pathway (K02440) that facilitates glycerol uptake [71]. Pathways involved in amino acid degradation were also significantly enriched in L3ir. The amino acids include both ketogenic (leucine, lysine and isoleucine) and glucogenic (histidine and cysteine) amino acids. The gene encoding dihydrolipoamide branched chain transacylase (DBT, k00680, E E2.3.1.-) represented by two oligos elements on the array (BMC03018; 14977.m05029, Pub_Locus-Bm1_43910) was also over-expressed in L3ir relative to L3c (3 fold). DBT is involved in the metabolism of amino acids (ko00280, ko00340, ko00310) other than cysteine (ko00272) and also involved in pathways "xenobiotics biodegradation and metabolism" (ko00642, ko00624 and ko00632) and "biosynthesis of secondary metabolites" (ko00903 and ko00960). It is unlikely that all of these pathways are altered in L3ir as a result of up-regulation of DBT. However, there are good reasons to believe that some of the pathways linked to DBT are altered in L3ir. For example, the biosynthesis of many classes of secondary metabolites in plants can be induced by stress or environmental changes [72–74]. In addition, exposure to UV radiation induces the biosynthesis of UV-absorbing compounds such as secondary metabolites (flavonoids) in plants [75]. Cysteine has a crucial role in inducible, endogenous detoxication mechanisms in the body [76]. Taken together, stress responses or IR seem to be appropriately altered in L3ir.

Table 9 Significant KEGG pathways identified in the L3ir gene set
Full size table

The statistically enriched gene ontology list for L3ir genes shared some features with the L3i list [See additional file 12]. However, most of the L3ir up-regulated genes were located in cytoplasm (GO:0005737) in contrast to extracellular region (GO:0005576) for L3i genes. L3ir genes had fewer significant GO terms than L3i or L3c. GO terms related to catalytic activity and enzyme regulatory activity were not enriched in the L3ir gene set.

The only KEGG pathway mapped in the L3c gene set was "ribosome for protein translation". In contrast, L3c had a long list of GO terms [See additional file 13]. However, most of the up-regulated genes were involved in biological processes of protein synthesis and growth such as translation (GO:0006412) and cell development (GO:0048468). Coordinately, the molecular functions related to these biological events were also enriched such as "structural constituent of ribosome" (GO: 0003735), "binding activities of nucleic acid binding" (GO:0003676), and "nucleotide binding" (GO:0000166)). In contrast to L3i, most genes up-regulated in L3c appear to be intracellular (GO:0005622) in various compartments such as organelle (GO:0044446) and microtubule cytoskeleton (GO:0015630).

Additional comparisons of BmL3 gene expression results with those reported for C. elegans and hookworm larvae

We compared gene expression data from BmL3 with data from infective larvae of A. caninum (Ac) and L3 dauer larvae of C. elegans. Datu et al identified 602 mRNAs differentially transcribed between non-activated- and activated-Ac L3 by serum [17]. We compared transcription profiles of the non-activated AcL3 and mosquito derived BmL3 (L3i). The genes coding cytochrome C oxidase, genes encoding PRPs family member ASP, and alpha-crystallin domain containing proteins (small heat shock proteins except hsp-12.6) were up-regulated in both non-activated AcL3 and BmL3i. However, genes encoding catalytic proteases such as cysteine, serine and metalloprotease were up-regulated in BmL3i relative to BmL3c, which contrasts with finding that these genes were up-regulated in activated AcL3 relative to non-activated AcL3. In addition, the lack of similarity between activated AcL3 and BmL3c was observed in gene-by-gene comparisons. For example, PRPs and cysteine protease were up-regulated in activated AcL3 but not in BmL3c. Also BmL3i actively expressed genes involved in invasion and immune evasion, while genes that encode such PRPs were activation-associated in AcL3. The biological significance of these differences remains to be determined. However, it is important to note that AcL3 are free living organisms while BmL3 are parasitic in mosquitoes.

SAGE and microarray studies have identified gene expression profiles associated with dauer arrest and exit in C. elegans [49, 51]. A total of 1,984 genes identified as dauer exit-specific genes had been further divided into five subclasses: dauer-enriched, transient, early climbing, and late. We compared BmL3 expression data to the dauer exit gene lists, and concordant functional classes are shown in Table 10. Genes involved in stress resistance and lipid metabolism was up-regulated in BmL3i and also in the dauer-enriched subclass. Other BmL3i up-regulated genes were identified in other subclasses especially the "late climbing group". Of course, there were also important differences between BmL3i and dauer exit gene lists. For example, the dauer-enriched gene encoding cytochrome P450 linked to dauer formation [77, 78] was not up-regulated in B. malayi L3i despite the presence of eight elements on the array annotated as cytochrome P450. Furthermore, Ce-hsp-12.6 (designated DAF-16) encodes a Forkhead transcription factor and is a member of the transient gene class [79]. Daf-16 promotes dauer formation and functions downstream of the daf-2/insulin-like receptor signal transduction pathway that regulates dauer development and longevity. There are 6 ESTs annotated as B. malayi homologs of small heat shock 12.6 on the array, but none was up-regulated in L3i. Interestingly, the transcription patterns for these two genes are similar in AcL3 and C. elegans dauer larvae [16, 17]. This suggests that non-activated AcL3 may be more similar to C. elegans dauer than BmL3. Again, this may be related to the fact that both AcL3 and C. elegans dauer are free-living nematode larvae while BmL3i are parasitic in mosquitoes.

Table 10 Overlap of functional classes in L3i genes and C. elegans dauer exit-specific genes*
Full size table

Consistency with prior studies and confirmation by real-time qRT-PCR

The fact that the findings from this study are consistent with previous reports tends to validate our gene expression profiles. For example, prior studies have shown those BmL3i genes such as cathepsin L-like proteases-5 (Bm-cpl-5), serine protease inhibitor-1 (Bm-spn-1), and several Bm-alt family genes are highly expressed in BmL3i [41, 55, 80]. Irradiation responsive genes in the BmL3ir gene set such as cytidine deaminase, uridine phosphorylase and collagens are also up-regulated after irradiation in other organisms [58–60, 62]. In addition, expression profiles were confirmed by real-time qRT-PCR for selected gene candidates as shown in additional file 14.

Conclusion

In summary, our data show that the transition of BmL3 from the mosquito host to culture conditions that mimic the mammalian host is accompanied by important changes in gene expression associated with several biochemical pathways and functions. Major differences involve genes that encode proteins with various functions such as invasion and immune evasion in L3i and growth and development in L3c. Irradiation significantly altered these changes with persistent expression of many L3i genes that are normally down-regulated in culture and increased expression of irradiation response genes. These alterations in gene expression may contribute to the efficacy of L3ir vaccines. Our data support the findings of prior studies that suggested that L3i enriched genes such as Bm-CPI-2 and Bm-ALT-2 may be involved in establishment of infection and immune evasion. In addition, our analysis greatly extends the number of candidate genes involved in such roles. L3i genes encode proteins newly identified in the secretome of Brugia adult worms [81] such as transthyretin-like family proteins (TLP) (TC3177), triose phosphate isomerase (TPI) (BMC11888) and galectin (BMC04762 and BMC12477), each of which fully warrant further investigation. Many of the genes up-regulated in the different L3 types have C. elegans homologues with severe RNAi phenotypes [See additional file 7] and some have been previously suggested as potential drug targets [See additional file 15] [82]. Altered pathways associated with L3 types provide clues regarding the molecular basis of important aspects of BmL3 biology such as invasion and evasion (for BmL3i), and growth and development (for BmL3c).

Filarial genomics is a work in progress. The B. malayi genome contains all of the puzzle pieces. This study represents a step in the long process of assembling these pieces to provide insight into the biology of filarial nematodes.

Methods

Parasite materials

BmL3 were isolated by dissection from Aedes aegypti mosquitoes that had been fed 14 days earlier on microfilaremic jirds as previously described [83]. Infective larvae were washed three times with RPMI-1640 (GIBCO-BRL, Grand Island, NY and immediately snap frozen on dry ice and stored at -80°C as mosquito derived infective larvae (L3i) as previously described [5]. Parasite cultures were performed as previously described [84]. Briefly, 2000–2500 L3i were incubated in 5 ml complete NI-medium [NCTC 135 (Sigma Chemical Co., St Louis, MO) and IMDM (GIBCO-BRL, Life Technologies, Grand Island, NY) (1:1) containing 2 mM glutamine, 200 U/ml streptomycin, 200 ug/ml penicillin, 0.5 ug/ml amphotericin B, and 10% fetal calf serum in 15 ml conical polystyrene tubes (Costar, Cambridge, Mass) at 37°C in a 95% air, 5% CO2 atmosphere for two days. Cultured worms (L3c) were recovered after 2 days, frozen immediately, and stored at -80°C. For irradiated L3 (L3ir), L3i were placed in 2 ml of culture medium in a 5 ml polystyrene tube and exposed to a laboratory certified 137 Cesium source (J.L. Shepherd and Associates, San Fernando, CA) at a dose of 75 krad and then cultured for 2 days as described above for L3c.

RNA isolation and probe preparation

Total RNA from L3i, L3c and L3ir parasites was prepared using TRIzol (GibcoBRL, Life Technologies) as previously described [5]. cDNA was synthesized from 5–7 ug total RNA from each sample using 3DNA capture sequence primers (3DNA Array 350 Detection system, Genisphere, Hatfield, PA) and SuperScript II Reverse Transcriptase (Gibco BRL, Gaithersburg, MD) for each probe according to standard protocols. cDNA was concentrated by Microcon YM-100 filter (Millipore) and either used immediately or stored at -80°C. cDNA was synthesized from two different RNA samples from each L3 type. A two-step protocol was used for hybridization (3DNA Array 350 Detection system, Genisphere, Hatfield, PA) as previously described [85]. Each experiment consisted of a pair-wise competitive hybridization of cDNA samples (L3i/L3c and L3ir/L3c) with reciprocal dye-flip replicates. Since biological replicates and dye-flip replicates were tested, a total of four DNA microarrays were used for each comparison of two types of cDNA.

Microarray fabrication

The BmV2array contains 18,104 elements derived from B. malayi (15,412), Onchocerca volvulus (O. volvulus, 1,016), Wuchereria bancrofti (W. bancrofti, 872) and Wolbachia (wBm, 804 genomic elements). The Brugia elements represent all B. malayi non-redundant gene models annotated in the course of the B. malayi genome project http://www.tigr.org/tdb/e2k1/bma1/bma1.shtml and ESTs from the TIGR Gene Indices (B. malayi v. 5.0 http://compbio.dfci.harvard.edu/tgi/tgipage.html) that did not match B. malayi gene models (which may represent un-annotated genes). Other elements represented O. volvulus and W. bancrofti ESTs that differ significantly from known B. malayi genes. The array also displays oligos designed on each wBm gene model. ArrayOligoselector http://arrayoligosel.sourceforge.net/[86] was used according to default parameters to choose unique 65 mers for each gene model and EST. The oligos were screened to avoid cross-hybridization to non-redundant human, bacterial and viral data nucleotide databases http://www.tigr.org/software/grid.shtml, and synthesized by standard methods by Illumina (San Diego, CA). The oligonucleotides (50 nM in 3× SSC with 0.75 M betaine) were printed in duplicate on MWG Epoxy slides (MWG Bioteche Inc, High Point, NC) by a locally constructed linear servo-arrayer (after the DeRisi model, http://derisilab.ucsf.edu/).

Data processing and analysis

Slide scanning and image analysis were performed as described previously (85). Briefly, slides were scanned immediately after hybridization on a ScanArray Express HT Scanner (Perkin Elmer, Boston, MA) to detect Cy3 and Cy5 fluorescence at 543 and 633 nm, respectively. The scanner produces green Cy3 and red Cy5 16-bit TIFF image files and extracts intensities from the scanner images for both dyes. The resultant values were background subtracted and Lowess [86, 87] normalized by using GeneSpring version 6.1 software (Silicon Genetics, Redwood City, CA). Twenty percent of the data were used to calculate the Lowess fit at each point. Genes with fold differences ≥ 2 and P < 0.01 based on pair-wise comparisons (L3i/L3c, L3c/L3ir) were considered to be differentially expressed. The best P-value obtained with high PMT scans and its accompanying arithmetic ratio was reported.

Annotation and Functional Assignments

A similarity search was performed using WU-TBLASTX https://blast.wustl.edu[88] with the 18,104 elements on the chip as queries versus the NCBI NR database (downloaded on Nov 14, 2006). Homologies were reported with a probability of 1E-05 or better. To identify cases where Brugia homologs in C. elegans have been surveyed for knockout phenotypes using RNA interference (RNAi) [18, 19, 21], Wormpep BLAST matches were cross referenced to a list of 19, 793 C. elegans genes with available RNAi information (Version 156). For each Brugia sequence, only the highest-scoring C. elegans match was considered.

Default parameters for InterProScan v13.1 were used to search against the InterPro database [89]. Gene ontology (GO) was assigned and displayed graphically by AmiGO with default parameters and ontology data released on March 15, 2007 [15]. For each GO term, its enrichment in a subset group (such as the L3i group) was measured over background using a hypergeometric analysis with the P-value cutoff e-05. Less informative ontology terms, including those at level 4 or higher for Biological Process or Molecular Function, and those at level 2 or higher for Cellular Component, were removed from the enrichment list.

An E-value cut-off of 1.0e-10 reported by WU-BLASTP against Genes Database Release 39.0 from Kyoto Encyclopedia of Genes and Genomes was used for pathway mapping; the top match and all the matches within a range of 30% of the top BLAST score that met the cut-off were accepted as valid KEGG associations [90, 91]. A hypergeometric analysis (measuring the coverage of KEGG Orthology or KO gene groups for each KEGG pathway compared to the complete gene set on the V2 chip) was implemented to identify enriched pathways for each gene set [92].

Validation of transcription by real-time qRT-PCR

Reverse transcription real-time PCR was used for the validation of microarray data as preciously described [85]. The gene coding for histone (BMC12346) was verified and used as endogenous control as previously reported. The sequences of all of the primers used in the real-time PCR are listed in additional file 16.

Abbreviations

L3i:

infective L3

L3ir:

irradiated and cultured L3

L3c:

cultured L3

KEGG:

Kyoto Encyclopedia of Genes and Genomes

qRT-PCR:

quantitative reverse transcription polymerase chain reaction

RNAi:

RNA interference.

References

  1. 1.

    Bianco AE, Robertson BD, Kuo YM, Townson S, Ham PJ: Developmentally regulated expression and secretion of a polymorphic antigen by Onchocerca infective-stage larvae. Mol Biochem Parasitol. 1990, 39: 203-211. 10.1016/0166-6851(90)90059-U.

  2. 2.

    Devaney E, Martin SA, Thompson FJ: Stage-specific gene expression in lymphatic filarial nematodes. Parasitol Today. 1996, 12: 418-424. 10.1016/0169-4758(96)10065-X.

  3. 3.

    Proudfoot L, Kusel JR, Smith HV, Harnett W, Worms MJ, Kennedy MW: Rapid changes in the surface of parasitic nematodes during transition from pre- to post-parasitic forms. Parasitology. 1993, 107 (Pt 1): 107-117. 10.1017/S0031182000079464.

  4. 4.

    Blaxter ML, Raghavan N, Ghosh I, Guiliano D, Lu W, Williams SA, Slatko B, Scott AL: Genes expressed in Brugia malayi infective third stage larvae. Mol Biochem Parasitol. 1996, 77: 77-93. 10.1016/0166-6851(96)02571-6.

  5. 5.

    Li BW, Rush AC, Weil GJ, McCarter JP, Mitreva M: Brugia malayi: effects of radiation and culture on gene expression in infective larvae. Mol Biochem Parasitol. 2006, 149: 201-207. 10.1016/j.molbiopara.2006.05.013.

  6. 6.

    Devaney E, Osborne J: The third-stage larva (L3) of Brugia: its role in immune modulation and protective immunity. Microbes & Infections. 2000, 1363-1371. 10.1016/S1286-4579(00)01290-9. 11

  7. 7.

    Eisenbeiss WF, Apfel H, Meyer TF: Protective immunity linked with a distinct developmental stage of a filarial parasite. J Immunol. 1994, 152: 735-742.

  8. 8.

    Gnanasekar M, Rao KV, He YX, Mishra PK, Nutman TB, Kaliraj P, Ramaswamy K: Novel phage display-based subtractive screening to identify vaccine candidates of Brugia malayi. Infect Immun. 2004, 72: 4707-4715. 10.1128/IAI.72.8.4707-4715.2004.

  9. 9.

    Oothuman P, Denham D, McGreevy P: Successful vaccination of cats against Brugia pahangi with larvae attenuated by irradiation with 10 kRas Cobalt 60. Parasite Immunol. 1979, 1: 209-216. 10.1111/j.1365-3024.1979.tb00707.x.

  10. 10.

    Yates J, Hagashi G: : Vaccination of jirds with 60 Cobalt-attenuated infective stage larvae protects against homologous challenge. Brugia malayi. 1985, 34: 1132-1137.

  11. 11.

    Weil GJ, Li BW, Liftis F, Chandrashekar R: Brugia malayi: antibody responses to larval antigens in infected and immunized jirds. Exp Parasitol. 1992, 74: 315-323. 10.1016/0014-4894(92)90155-4.

  12. 12.

    Li BW, Chandrashekar R, Alvarez RM, Liftis F, Weil GJ: Identification of paramyosin as a potential protective antigen against Brugia malayi infection in jirds. Mol Biochem Parasitol. 1991, 49: 315-323. 10.1016/0166-6851(91)90075-H.

  13. 13.

    Ghedin E, Wang S, Spiro D, Caler E, Zhao Q, Crabtree J, Allen JE, Delcher AL, Guiliano DB, Miranda-Saavedra D, et al: Draft genome of the filarial nematode parasite Brugia malayi. Science. 2007, 317: 1756-1760. 10.1126/science.1145406.

  14. 14.

    Camon E, Magrane M, Barrell D, Binns D, Fleischmann W, Kersey P, Mulder N, Oinn T, Maslen J, Cox A, Apweiler R: The Gene Ontology Annotation (GOA) project: implementation of GO in SWISS-PROT, TrEMBL, and InterPro. Genome Res. 2003, 13: 662-672. 10.1101/gr.461403.

  15. 15.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.

  16. 16.

    Moser JM, Freitas T, Arasu P, Gibson G: Gene expression profiles associated with the transition to parasitism in Ancylostoma caninum larvae. Mol Biochem Parasitol. 2005, 143: 39-48. 10.1016/j.molbiopara.2005.04.012.

  17. 17.

    Datu BJ, Gasser RB, Nagaraj SH, Ong EK, O'Donoghue P, McInnes R, Ranganathan S, Loukas A: Transcriptional Changes in the Hookworm, Ancylostoma caninum, during the transition from a free-living to a parasitic larva. PLoS Negl Trop Dis. 2008, 2: e130-10.1371/journal.pntd.0000130.

  18. 18.

    Kamath RS, Fraser AG, Dong Y, Poulin G, Durbin R, Gotta M, Kanapin A, Le Bot N, Moreno S, Sohrmann M, et al: Systematic functional analysis of the Caenorhabditis elegans genome using RNAi. Nature. 2003, 421: 231-237. 10.1038/nature01278.

  19. 19.

    Fraser AG, Kamath RS, Zipperlen P, Martinez-Campos M, Sohrmann M, Ahringer J: Functional genomic analysis of C. elegans chromosome I by systematic RNA interference. Nature. 2000, 408: 325-330. 10.1038/35042517.

  20. 20.

    Simmer F, Moorman C, Linden van der AM, Kuijk E, Berghe van den PV, Kamath RS, Fraser AG, Ahringer J, Plasterk RH: Genome-wide RNAi of C. elegans using the hypersensitive rrf-3 strain reveals novel gene functions. PLoS Biol. 2003, 1: E12-10.1371/journal.pbio.0000012.

  21. 21.

    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.

  22. 22.

    Mitreva M, McCarter JP, Martin J, Dante M, Wylie T, Chiapelli B, Pape D, Clifton SW, Nutman TB, Waterston RH: Comparative genomics of gene expression in the parasitic and free-living nematodes Strongyloides stercoralis and Caenorhabditis elegans. Genome Res. 2004, 14: 209-220. 10.1101/gr.1524804.

  23. 23.

    Maizels RM, Balic A, Gomez-Escobar N, Nair M, Taylor MD, Allen JE: Helminth parasites – masters of regulation. Immunol Rev. 2004, 201: 89-116. 10.1111/j.0105-2896.2004.00191.x.

  24. 24.

    Semnani RT, Law M, Kubofcik J, Nutman TB: Filaria-induced immune evasion: suppression by the infective stage of Brugia malayi at the earliest host-parasite interface. J Immunol. 2004, 172: 6229-6238.

  25. 25.

    Gomez-Escobar N, Bennett C, Prieto-Lafuente L, Aebischer T, Blackburn CC, Maizels RM: Heterologous expression of the filarial nematode alt gene products reveals their potential to inhibit immune function. BMC Biol. 2005, 3: 8-10.1186/1741-7007-3-8.

  26. 26.

    Hartmann S, Kyewski B, Sonnenburg B, Lucius R: A filarial cysteine protease inhibitor down-regulates T cell proliferation and enhances interleukin-10 production. Eur J Immunol. 1997, 27: 2253-2260. 10.1002/eji.1830270920.

  27. 27.

    Schierack P, Lucius R, Sonnenburg B, Schilling K, Hartmann S: Parasite-specific immunomodulatory functions of filarial cystatin. Infect Immun. 2003, 71: 2422-2429. 10.1128/IAI.71.5.2422-2429.2003.

  28. 28.

    Manoury B, Gregory WF, Maizels RM, Watts C: Bm-CPI-2, a cystatin homolog secreted by the filarial parasite Brugia malayi, inhibits class II MHC-restricted antigen processing. Curr Biol. 2001, 11: 447-451. 10.1016/S0960-9822(01)00118-X.

  29. 29.

    Zang X, Maizels RM: Serine proteinase inhibitors from nematodes and the arms race between host and pathogen. Trends Biochem Sci. 2001, 26: 191-197. 10.1016/S0968-0004(00)01761-8.

  30. 30.

    Ou X, Tang L, McCrossan M, Henkle-Duhrsen K, Selkirk ME: Brugia malayi: localisation and differential expression of extracellular and cytoplasmic CuZn superoxide dismutases in adults and microfilariae. Exp Parasitol. 1995, 80: 515-529. 10.1006/expr.1995.1064.

  31. 31.

    Ghosh I, Eisinger SW, Raghavan N, Scott AL: Thioredoxin peroxidases from Brugia malayi. Mol Biochem Parasitol. 1998, 91: 207-220. 10.1016/S0166-6851(97)00213-2.

  32. 32.

    Tracy JW, Catto BA, Webster LT: Reductive metabolism of niridazole by adult Schistosoma mansoni. Correlation with covalent drug binding to parasite macromolecules. Mol Pharmacol. 1983, 24: 291-299.

  33. 33.

    Gupta S, Rathaur S: Filarial glutathione S-transferase: its induction by xenobiotics and potential as drug target. Acta Biochim Pol. 2005, 52: 493-500.

  34. 34.

    Hyndman D, Bauman DR, Heredia VV, Penning TM: The aldo-keto reductase superfamily homepage. Chem Biol Interact. 2003, 143–144: 621-631. 10.1016/S0009-2797(02)00193-X.

  35. 35.

    Burnell AM, Houthoofd K, O'Hanlon K, Vanfleteren JR: Alternate metabolism during the dauer stage of the nematode Caenorhabditis elegans. Exp Gerontol. 2005, 40: 850-856. 10.1016/j.exger.2005.09.006.

  36. 36.

    Raha S, Robinson BH: Mitochondria, oxygen free radicals, disease and ageing. Trends Biochem Sci. 2000, 25: 502-508. 10.1016/S0968-0004(00)01674-1.

  37. 37.

    Selkirk ME, Smith VP, Thomas GR, Gounaris K: Resistance of filarial nematode parasites to oxidative stress. Int J Parasitol. 1998, 28: 1315-1332. 10.1016/S0020-7519(98)00107-6.

  38. 38.

    Ou X, Thomas GR, Chacon MR, Tang L, Selkirk ME: Brugia malayi: differential susceptibility to and metabolism of hydrogen peroxide in adults and microfilariae. Exp Parasitol. 1995, 80: 530-540. 10.1006/expr.1995.1065.

  39. 39.

    Richer JK, Sakanari JA, Frank GR, Grieve RB: Dirofilaria immitis: proteases produced by third- and fourth-stage larvae. Exp Parasitol. 1992, 75: 213-222. 10.1016/0014-4894(92)90181-9.

  40. 40.

    Lustigman S, Zhang J, Liu J, Oksov Y, Hashmi S: RNA interference targeting cathepsin L and Z-like cysteine proteases of Onchocerca volvulus confirmed their essential function during L3 molting. Mol Biochem Parasitol. 2004, 138: 165-170. 10.1016/j.molbiopara.2004.08.003.

  41. 41.

    Guiliano DB, Hong X, McKerrow JH, Blaxter ML, Oksov Y, Liu J, Ghedin E, Lustigman S: A gene family of cathepsin L-like proteases of filarial nematodes are associated with larval molting and cuticle and eggshell remodeling. Mol Biochem Parasitol. 2004, 136: 227-242. 10.1016/j.molbiopara.2004.03.015.

  42. 42.

    Williamson AL, Lustigman S, Oksov Y, Deumic V, Plieskatt J, Mendez S, Zhan B, Bottazzi ME, Hotez PJ, Loukas A: Ancylostoma caninum MTP-1, an astacin-like metalloprotease secreted by infective hookworm larvae, is involved in tissue migration. Infect Immun. 2006, 74: 961-967. 10.1128/IAI.74.2.961-967.2006.

  43. 43.

    McKerrow JH, Salter J: Invasion of skin by Schistosoma cercariae. Trends Parasitol. 2002, 18: 193-195. 10.1016/S1471-4922(02)02309-7.

  44. 44.

    Murray J, Gregory WF, Gomez-Escobar N, Atmadja AK, Maizels RM: Expression and immune recognition of Brugia malayi VAL-1, a homologue of vespid venom allergens and Ancylostoma secreted proteins. Mol Biochem Parasitol. 2001, 118: 89-96. 10.1016/S0166-6851(01)00374-7.

  45. 45.

    Henriksen A, King TP, Mirza O, Monsalve RI, Meno K, Ipsen H, Larsen JN, Gajhede M, Spangfort MD: Major venom allergen of yellow jackets, Ves v 5: structural characterization of a pathogenesis-related protein superfamily. Proteins. 2001, 45: 438-448. 10.1002/prot.1160.

  46. 46.

    Hawdon JM, Jones BF, Hoffman DR, Hotez PJ: Cloning and characterization of Ancylostoma-secreted protein. A novel protein associated with the transition to parasitism by infective hookworm larvae. J Biol Chem. 1996, 271: 6672-6678. 10.1074/jbc.271.12.6672.

  47. 47.

    MacDonald AJ, Tawe W, Leon O, Cao L, Liu J, Oksov Y, Abraham D, Lustigman S: Ov-ASP-1, the Onchocerca volvulus homologue of the activation associated secreted protein family is immunostimulatory and can induce protective anti-larval immunity. Parasite Immunol. 2004, 26: 53-62. 10.1111/j.0141-9838.2004.00685.x.

  48. 48.

    Raghavan N, Ghosh I, Eisinger WS, Pastrana D, Scott AL: Developmentally regulated expression of a unique small heat shock protein in Brugia malayi. Mol Biochem Parasitol. 1999, 104: 233-246. 10.1016/S0166-6851(99)00150-4.

  49. 49.

    Jones SJ, Riddle DL, Pouzyrev AT, Velculescu VE, Hillier L, Eddy SR, Stricklin SL, Baillie DL, Waterston R, Marra MA: Changes in gene expression associated with developmental arrest and longevity in Caenorhabditis elegans. Genome Res. 2001, 11: 1346-1352. 10.1101/gr.184401.

  50. 50.

    McElwee JJ, Schuster E, Blanc E, Thomas JH, Gems D: Shared transcriptional signature in Caenorhabditis elegans Dauer larvae and long-lived daf-2 mutants implicates detoxification system in longevity assurance. J Biol Chem. 2004, 279: 44533-44543. 10.1074/jbc.M406207200.

  51. 51.

    Wang J, Kim SK: Global analysis of dauer gene expression in Caenorhabditis elegans. Development. 2003, 130: 1621-1634. 10.1242/dev.00363.

  52. 52.

    Aspock G, Kagoshima H, Niklaus G, Burglin TR: Caenorhabditis elegans has scores of hedgehog-related genes: sequence and expression analysis. Genome Res. 1999, 9: 909-923. 10.1101/gr.9.10.909.

  53. 53.

    Lustigman S, Brotman B, Johnson EH, Smith AB, Huima T, Prince AM: Identification and characterization of an Onchocerca volvulus cDNA clone encoding a microfilarial surface-associated antigen. Mol Biochem Parasitol. 1992, 50: 79-93. 10.1016/0166-6851(92)90246-G.

  54. 54.

    Chandrashekar R, Curtis KC, Ramzy RM, Liftis F, Li BW, Weil GJ: Molecular cloning of Brugia malayi antigens for diagnosis of lymphatic filariasis. Mol Biochem Parasitol. 1994, 64: 261-271. 10.1016/0166-6851(94)00035-2.

  55. 55.

    Joseph GT, McCarthy JS, Huima T, Mair KF, Kass PH, Boussinesq M, Goodrick L, Bradley JE, Lustigman S: Onchocera volvulus: characterization of a highly immunogenic Gln-rich protein. Mol Biochem Parasitol. 1997, 90: 55-68. 10.1016/S0166-6851(97)00132-1.

  56. 56.

    Garraud O, Nkenfou C, Bradley JE, Perler FB, Nutman TB: Identification of recombinant filarial proteins capable of inducing polyclonal and antigen-specific IgE and IgG4 antibodies. J Immunol. 1995, 155: 1316-1325.

  57. 57.

    Tsuji N, Suzuki K, Kasuga-Aoki H, Isobe T, Arakawa T, Matsumoto Y: Mice intranasally immunized with a recombinant 16-kilodalton antigen from roundworm Ascaris parasites are protected against larval migration of Ascaris suum. Infect Immun. 2003, 71: 5314-5323. 10.1128/IAI.71.9.5314-5323.2003.

  58. 58.

    Rieger KE, Chu G: Portrait of transcriptional responses to ultraviolet and ionizing radiation in human cells. Nucleic Acids Res. 2004, 32: 4786-4803. 10.1093/nar/gkh783.

  59. 59.

    Khodarev NN, Park JO, Yu J, Gupta N, Nodzenski E, Roizman B, Weichselbaum RR: Dose-dependent and independent temporal patterns of gene responses to ionizing radiation in normal and tumor cells and tumor xenografts. Proc Natl Acad Sci USA. 2001, 98: 12665-12670. 10.1073/pnas.211443698.

  60. 60.

    Otomo T, Hishii M, Arai H, Sato K, Sasai K: Microarray analysis of temporal gene responses to ionizing radiation in two glioblastoma cell lines: up-regulation of DNA repair genes. J Radiat Res (Tokyo). 2004, 45: 53-60. 10.1269/jrr.45.53.

  61. 61.

    Remy J, Wegrowski J, Crechet F, Martin M, Daburon F: Long-term overproduction of collagen in radiation-induced fibrosis. Radiat Res. 1991, 125: 14-19. 10.2307/3577976.

  62. 62.

    Followill DS, Travis EL: Differential expression of collagen types I and III in consequential and primary fibrosis in irradiated mouse colon. Radiat Res. 1995, 144: 318-328. 10.2307/3578952.

  63. 63.

    Selkirk M, L KC, Partono F, Sayers G, Maizels G: Identification, synthesis and immunogenicity of cuticular collagens from the filarial nematodes Brugia malagi and Brugia pahang. Mol Biochem Parasitol. 1989, 32: 230-246. 10.1016/0166-6851(89)90073-X.

  64. 64.

    Petralanda I, Piessens WF: Onchocerca volvulus, O. gutturosa and, Brugia malayi, and Dirofilaria immitis: a comparative study of the immunochemical properties of cuticular proteins from filarial parasites. Exp Parasitol. 1991, 72: 164-173. 10.1016/0014-4894(91)90134-I.

  65. 65.

    Stewart GR, Zhu Y, Parredes W, Tree TI, Guderian R, Bradley JE: The novel cuticular collagen Ovcol-1 of Onchocerca volvulus is preferentially recognized by immunoglobulin G3 from putatively immune individuals. Infect Immun. 1997, 65: 164-170.

  66. 66.

    Fetterer RH, Hill DE, Urban JF: The cuticular biology in developmental stages of Ascaris suum. Acta Trop. 1990, 47: 289-295. 10.1016/0001-706X(90)90030-4.

  67. 67.

    Frand AR, Russel S, Ruvkun G: Functional genomic analysis of C. elegans molting. PLoS Biol. 2005, 3: e312-10.1371/journal.pbio.0030312.

  68. 68.

    Singh JES: Some observations on molting in Caenorhabditis elegans. Nematologica. 1978, 24: 63-71.

  69. 69.

    Chirgwin SR, Coleman SU, Klei TR: Brugia pahangi: In vivo tissue migration of early L3 alters gene expression. Exp Parasitol. 2008, 118: 89-95. 10.1016/j.exppara.2007.06.007.

  70. 70.

    O'riordan V, Burnell A: Intermediary metabolism in the dauer larva of the nematode Caenorhabditis elegans – The Glyoxylate cycle and fatty-acidoxidation. Comp Biochem Physiol. 1990, 95: 125-130.

  71. 71.

    Grayson P, Tajkhorshid E, Schulten K: Mechanisms of selectivity in channels and enzymes studied with interactive molecular dynamics. Biophys J. 2003, 85: 36-48. 10.1016/S0006-3495(03)74452-X.

  72. 72.

    Fits Van der L, Memelink J: ORCA3, a jasmonate-responsive transcriptional regulator of plant primary and secondary metabolism. Science. 2000, 289: 295-297. 10.1126/science.289.5477.295.

  73. 73.

    Kutchan TM: Alkaloid Biosynthesis[mdash]The Basis for Metabolic Engineering of Medicinal Plants. Plant Cell. 1995, 7: 1059-1070. 10.1105/tpc.7.7.1059.

  74. 74.

    Blechert S, Brodschelm W, Holder S, Kammerer L, Kutchan TM, Mueller MJ, Xia ZQ, Zenk MH: The octadecanoic pathway: signal molecules for the regulation of secondary pathways. Proc Natl Acad Sci USA. 1995, 92: 4099-4105. 10.1073/pnas.92.10.4099.

  75. 75.

    Lavola A: Accumulation of flavonoids and related compounds in birch induced by UV-B irradiance. Tree Physiol. 1998, 18: 53-58.

  76. 76.

    Quig D: Cysteine metabolism and metal toxicity. Altern Med Rev. 1998, 3: 262-270.

  77. 77.

    Gerisch B, Antebi A: Hormonal signals produced by DAF-9/cytochrome P450 regulate C. elegans dauer diapause in response to environmental cues. Development. 2004, 131: 1765-1776. 10.1242/dev.01068.

  78. 78.

    Jia K, Albert PS, Riddle DL: DAF-9, a cytochrome P450 regulating C. elegans larval development and adult longevity. Development. 2002, 129: 221-231.

  79. 79.

    Ogg S, Paradis S, Gottlieb S, Patterson GI, Lee L, Tissenbaum HA, Ruvkun G: The Fork head transcription factor DAF-16 transduces insulin-like metabolic and longevity signals in C. elegans. Nature. 1997, 389: 994-999. 10.1038/40194.

  80. 80.

    Gregory WF, Blaxter ML, Maizels RM: Differentially expressed, abundant trans-spliced cDNAs from larval Brugia malayi. Mol Biochem Parasitol. 1997, 87: 85-95. 10.1016/S0166-6851(97)00050-9.

  81. 81.

    Hewitson JP, Harcus YM, Curwen RS, Dowle AA, Atmadja AK, Ashton PD, Wilson A, Maizels RM: The secretome of the filarial parasite, Brugia malayi: proteomic profile of adult excretory-secretory products. Mol Biochem Parasitol. 2008, 160: 8-21. 10.1016/j.molbiopara.2008.02.007.

  82. 82.

    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.

  83. 83.

    MaCall J: The role of arthropods in the development of animal models for filaria research. Journal of Georgia Entomology Society. 1981, 16: 283-293.

  84. 84.

    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.

  85. 85.

    Li BW, Rush AC, Crosby SD, Warren WC, Williams SA, Mitreva M, Weil GJ: Profiling of gender-regulated gene transcripts in the filarial nematode Brugia malayi by cDNA oligonucleotide array analysis. Mol Biochem Parasitol. 2005, 143: 49-57. 10.1016/j.molbiopara.2005.05.005.

  86. 86.

    Bozdech Z, Zhu J, Joachimiak M, Cohen F, Pulliam BJD: Expression profiling of the schizont and trophozoite stages of Plasmodium falciparum with a long-oligonucleotide microarray. Genome Biol. 2003, 4: R9-10.1186/gb-2003-4-2-r9.

  87. 87.

    Clevel W, Devlin S: Locally-weighted regression: an approach to regression analysis by local fitting. J American Statistic. J American Statistic Association. 1988, 83: 596-610. 10.2307/2289282.

  88. 88.

    Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.

  89. 89.

    Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N: InterProScan: protein domains identifier. Nucleic Acids Res. 2005, 33: W116-120. 10.1093/nar/gki442.

  90. 90.

    Bono H, Ogata H, Goto S, Kanehisa M: Reconstruction of amino acid biosynthesis pathways from the complete genome sequence. Genome Res. 1998, 8: 203-210.

  91. 91.

    Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28: 27-30. 10.1093/nar/28.1.27.

  92. 92.

    Mao X, Cai T, Olyarchuk J, Wei L: Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary. Bioinformatics. 2005, 21: 3787-3793. 10.1093/bioinformatics/bti430.

Download references

Acknowledgements

The Brugia malayi Microarray Consortium worked together to plan and finance construction of the array. The authors wish to thank Michael Heinz for his excellent work in performing the hybridizations to the microarrays.

Author information

Affiliations

  1. Department of internal medicine, Washington University School of Medicine, St. Louis, MO, 63110, USA
    • Ben-Wen Li
    • , Amy C Rush
    •  & Gary J Weil
  2. The Genome Sequencing Center, Department of Genetics, Washington University School of Medicine, St. Louis, MO, 63110, USA
    • Makedonka Mitreva
    •  & Yong Yin
  3. J. Craig Venter Institute, 9704 Medical Center Drive, Rockville, MD, 20850, USA
    • David Spiro
  4. Divisions of Infectious Diseases, University of Pittsburgh School of Medicine, Pittsburgh, PA, 15261, USA
    • Elodie Ghedin
Authors
  1. Ben-Wen Li
    View author publications
    You can also search for this author in
  2. Amy C Rush
    View author publications
    You can also search for this author in
  3. Makedonka Mitreva
    View author publications
    You can also search for this author in
  4. Yong Yin
    View author publications
    You can also search for this author in
  5. David Spiro
    View author publications
    You can also search for this author in
  6. Elodie Ghedin
    View author publications
    You can also search for this author in
  7. Gary J Weil
    View author publications
    You can also search for this author in

Corresponding author

Correspondence to Ben-Wen Li.

Additional information

Authors' contributions

BWL: experimental design, interpretation of results, data analysis, drafted manuscript. ACR: RNA extraction, qRT-PCR, data mining. MM: genome mining, manuscript preparation. YY: genome mining, statistical analysis, manuscript preparation. DS: oligo calling, genome mining. EG: oligo calling, genome mining, manuscript preparation. GJW: experimental design, interpretation of results, manuscript preparation. All authors have read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Summary of BLAST output of BmV2array elements vs. NCBI_NR. The data provided represent a similarity search of BmV2 array sequences against NCBI non-redundant database (downloaded Nov 14, 2006). (XLS 3 MB)

Additional file 2: Complete RNAi information of C. elegans homologues of Brugia genes on BmV2array. The C. elegans homologues of Brugia genes on the BmV2array were found to have various RNAi phenotypes. (XLS 942 KB)

Additional file 3: Summary of protein domain matches for BmV2array elements. The protein domain matches of BmV2array against InterProScan v13.1 database. (XLS 826 KB)

Additional file 4: Summary of KEGG analysis of BmV2array. The top match of pathway mapping for BmV2array sequences. (XLS 143 KB)

Additional file 5: Top 20 protein domain matches for BmV2array elements. The data provided represent top 20 protein domains matched in BmV2array sequences. (DOCX 13 KB)

Additional file 6: Up-regulated genes in the L3 type. The up-regulated genes in pair-wise comparison of L3i vs L3c and L3ir vs L3c. (XLS 196 KB)

Additional file 7: RNAi phenotypes of C. elegans homologues in L3 the gene sets. The C. elegans homologues of Brugia genes up-regulated in L3 (L3i, L3ir and L3c) were found to have various RNAi phenotypes. (XLS 58 KB)

Additional file 8: Complete KEGG analysis in the L3 gene set. The KEGG pathway mapped in the L3 gene set. (XLS 117 KB)

Additional file 9: Significant protein domains in each L3 gene set. The most significant protein domains matched in each L3 gene set. (XLS 22 KB)

Additional file 10: Significant GO terms in the L3i gene set. The statistically significant Gene Ontology terms in the L3i gene set. (XLS 24 KB)

Additional file 11: L3i genes involved in gluconeogenesis. The genes up-regulated in L3i were involved in gluconeogenesis. (XLS 20 KB)

Additional file 12: Significant GO terms in the L3ir gene set. The statistically significant Gene Ontology terms were mapped in the L3ir gene set. (XLS 18 KB)

Additional file 13: Significant GO terms in the L3c gene set. The statistically significant Gene Ontology terms in the L3c gene set. (XLS 22 KB)

Additional file 14: Up-regulated genes confirmed by qRT-PCR. Genes up-regulated by microarray were confirmed by real-time qRT-PCR. (XLS 22 KB)

Additional file 15: Potential drug targets in each L3 gene set. The potential drug targets were identified in each L3 gene set. (XLS 52 KB)

Additional file 16: Primers of selected genes confirmed by real-time qRT-PCR. The primers used in real-time qRT-PCR for the confirmation. (XLS 24 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

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.

Reprints and Permissions

About this article

Cite this article

Li, B., Rush, A.C., Mitreva, M. et al. Transcriptomes and pathways associated with infectivity, survival and immunogenicity in Brugia malayi L3. BMC Genomics 10, 267 (2009). https://doi.org/10.1186/1471-2164-10-267

Download citation

Keywords

  • Gene Ontology
  • Immune Evasion
  • Mammalian Host
  • Small Heat Shock Protein
  • Filarial Nematode

BMC Genomics

ISSN: 1471-2164

Contact us