Document Detail

A transcriptomic analysis of the adult stage of the bovine lungworm, Dictyocaulus viviparus.
Jump to Full Text
MedLine Citation:
PMID:  17784965     Owner:  NLM     Status:  MEDLINE    
BACKGROUND: Lungworms of the genus Dictyocaulus (family Dictyocaulidae) are parasitic nematodes of major economic importance. They cause pathological effects and clinical disease in various ruminant hosts, particularly in young animals. Dictyocaulus viviparus, called the bovine lungworm, is a major pathogen of cattle, with severe infections being fatal. In this study, we provide first insights into the transcriptome of the adult stage of D. viviparus through the analysis of expressed sequence tags (ESTs).
RESULTS: Using our EST analysis pipeline, we estimate that the present dataset of 4436 ESTs is derived from 2258 genes based on cluster and comparative genomic analyses of the ESTs. Of the 2258 representative ESTs, 1159 (51.3%) had homologues in the free-living nematode C. elegans, 1174 (51.9%) in parasitic nematodes, 827 (36.6%) in organisms other than nematodes, and 863 (38%) had no significant match to any sequence in the current databases. Of the C. elegans homologues, 569 had observed 'non-wildtype' RNAi phenotypes, including embryonic lethality, maternal sterility, sterility in progeny, larval arrest and slow growth. We could functionally classify 776 (35%) sequences using the Gene Ontologies (GO) and established pathway associations to 696 (31%) sequences in Kyoto Encyclopedia of Genes and Genomes (KEGG). In addition, we predicted 85 secreted proteins which could represent potential candidates for developing novel anthelmintics or vaccines.
CONCLUSION: The bioinformatic analyses of ESTs data for D. viviparus has elucidated sets of relatively conserved and potentially novel genes. The genes discovered in this study should assist research toward a better understanding of the basic molecular biology of D. viviparus, which could lead, in the longer term, to novel intervention strategies. The characterization of the D. viviparus transcriptome also provides a foundation for whole genome sequence analysis and future comparative transcriptomic analyses.
Shoba Ranganathan; Shivashankar H Nagaraj; Min Hu; Christina Strube; Thomas Schnieder; Robin B Gasser
Related Documents :
15262285 - Annotated expressed sequence tags for studies of the regulation of reproductive modes i...
18852505 - Generation of expressed sequence tags for immune gene discovery and marker development ...
21303765 - Rapid changes in gene expression: dna determinants of promoter regulation by the concen...
20081835 - Parallel, tag-directed assembly of locally derived short sequence reads.
8973355 - Genomic organization of mouse dishevelled genes.
3023365 - Molecular cloning and developmental expression of the cyclic nucleotide phosphodiestera...
Publication Detail:
Type:  Journal Article; Research Support, Non-U.S. Gov't     Date:  2007-09-05
Journal Detail:
Title:  BMC genomics     Volume:  8     ISSN:  1471-2164     ISO Abbreviation:  BMC Genomics     Publication Date:  2007  
Date Detail:
Created Date:  2007-12-12     Completed Date:  2007-12-31     Revised Date:  2013-06-06    
Medline Journal Info:
Nlm Unique ID:  100965258     Medline TA:  BMC Genomics     Country:  England    
Other Details:
Languages:  eng     Pagination:  311     Citation Subset:  IM    
Department of Chemistry and Biomolecular Sciences, Macquarie University, Sydney, New South Wales 2109, Australia.
Data Bank Information
Bank Name/Acc. No.:
GDB/EV849926;  EV849927;  EV849928;  EV849929;  EV849930;  EV849931;  EV849932;  EV849933;  EV849934;  EV849935;  EV849936;  EV849937;  EV849938;  EV849939;  EV849940;  EV849941;  EV849942;  EV849943;  EV849944;  EV849945;  EV849946;  EV849947;  EV849948;  EV849949;  EV849950;  EV849951;  EV849952;  EV849953;  EV849954;  EV849955;  EV849956;  EV849957;  EV849958;  EV849959;  EV849960;  EV849961;  EV849962;  EV849963;  EV849964;  EV849965;  EV849966;  EV849967;  EV849968;  EV849969;  EV849970;  EV849971;  EV849972;  EV849973;  EV849974;  EV849975;  EV849976;  EV849977;  EV849978;  EV849979;  EV849980;  EV849981;  EV849982;  EV849983;  EV849984;  EV849985;  EV849986;  EV849987;  EV849988;  EV849989;  EV849990;  EV849991;  EV849992;  EV849993;  EV849994;  EV849995;  EV849996;  EV849997;  EV849998;  EV849999;  EV850000;  EV850001;  EV850002;  EV850003;  EV850004;  EV850005;  EV850006;  EV850007;  EV850008;  EV850009;  EV850010;  EV850011;  EV850012;  EV850013;  EV850014;  EV850015;  EV850016;  EV850017;  EV850018;  EV850019;  EV850020;  EV850021;  EV850022;  EV850023;  EV850024;  EV850025;  EV850026;  EV850027;  EV850028;  EV850029;  EV850030;  EV850031;  EV850032;  EV850033;  EV850034;  EV850035;  EV850036;  EV850037;  EV850038;  EV850039;  EV850040;  EV850041;  EV850042;  EV850043;  EV850044;  EV850045;  EV850046;  EV850047;  EV850048;  EV850049;  EV850050;  EV850051;  EV850052;  EV850053;  EV850054;  EV850055;  EV850056;  EV850057;  EV850058;  EV850059;  EV850060;  EV850061;  EV850062;  EV850063;  EV850064;  EV850065;  EV850066;  EV850067;  EV850068;  EV850069;  EV850070;  EV850071;  EV850072;  EV850073;  EV850074;  EV850075;  EV850076;  EV850077;  EV850078;  EV850079;  EV850080;  EV850081;  EV850082;  EV850083;  EV850084;  EV850085;  EV850086;  EV850087;  EV850088;  EV850089;  EV850090;  EV850091;  EV850092;  EV850093;  EV850094;  EV850095;  EV850096;  EV850097;  EV850098;  EV850099;  EV850100;  EV850101;  EV850102;  EV850103;  EV850104;  EV850105;  EV850106;  EV850107;  EV850108;  EV850109;  EV850110;  EV850111;  EV850112;  EV850113;  EV850114;  EV850115;  EV850116;  EV850117;  EV850118;  EV850119;  EV850120;  EV850121;  EV850122;  EV850123;  EV850124;  EV850125;  EV850126;  EV850127;  EV850128;  EV850129;  EV850130;  EV850131;  EV850132;  EV850133;  EV850134;  EV850135;  EV850136;  EV850137;  EV850138;  EV850139;  EV850140;  EV850141;  EV850142;  EV850143;  EV850144;  EV850145;  EV850146;  EV850147;  EV850148;  EV850149;  EV850150;  EV850151;  EV850152;  EV850153;  EV850154;  EV850155;  EV850156;  EV850157;  EV850158;  EV850159;  EV850160;  EV850161;  EV850162;  EV850163;  EV850164;  EV850165;  EV850166;  EV850167;  EV850168;  EV850169;  EV850170;  EV850171;  EV850172;  EV850173;  EV850174;  EV850175;  EV850176;  EV850177;  EV850178;  EV850179;  EV850180;  EV850181;  EV850182;  EV850183;  EV850184;  EV850185;  EV850186;  EV850187;  EV850188;  EV850189;  EV850190;  EV850191;  EV850192;  EV850193;  EV850194;  EV850195;  EV850196;  EV850197;  EV850198;  EV850199;  EV850200;  EV850201;  EV850202;  EV850203;  EV850204;  EV850205;  EV850206;  EV850207;  EV850208;  EV850209;  EV850210;  EV850211;  EV850212;  EV850213;  EV850214;  EV850215;  EV850216;  EV850217;  EV850218;  EV850219;  EV850220;  EV850221;  EV850222;  EV850223;  EV850224;  EV850225;  EV850226;  EV850227;  EV850228;  EV850229;  EV850230;  EV850231;  EV850232;  EV850233;  EV850234;  EV850235;  EV850236;  EV850237;  EV850238;  EV850239;  EV850240;  EV850241;  EV850242;  EV850243;  EV850244;  EV850245;  EV850246;  EV850247;  EV850248;  EV850249;  EV850250;  EV850251;  EV850252;  EV850253;  EV850254;  EV850255;  EV850256;  EV850257;  EV850258;  EV850259;  EV850260;  EV850261;  EV850262;  EV850263;  EV850264;  EV850265;  EV850266;  EV850267;  EV850268;  EV850269;  EV850270;  EV850271;  EV850272;  EV850273;  EV850274;  EV850275;  EV850276;  EV850277;  EV850278;  EV850279;  EV850280;  EV850281;  EV850282;  EV850283;  EV850284;  EV850285;  EV850286;  EV850287;  EV850288;  EV850289;  EV850290;  EV850291;  EV850292;  EV850293;  EV850294;  EV850295;  EV850296;  EV850297;  EV850298;  EV850299;  EV850300;  EV850301;  EV850302;  EV850303;  EV850304;  EV850305;  EV850306;  EV850307;  EV850308;  EV850309;  EV850310;  EV850311;  EV850312;  EV850313;  EV850314;  EV850315;  EV850316;  EV850317;  EV850318;  EV850319;  EV850320;  EV850321;  EV850322;  EV850323;  EV850324;  EV850325;  EV850326;  EV850327;  EV850328;  EV850329;  EV850330;  EV850331;  EV850332;  EV850333;  EV850334;  EV850335;  EV850336;  EV850337;  EV850338;  EV850339;  EV850340;  EV850341;  EV850342;  EV850343;  EV850344;  EV850345;  EV850346;  EV850347;  EV850348;  EV850349;  EV850350;  EV850351;  EV850352;  EV850353;  EV850354;  EV850355;  EV850356;  EV850357;  EV850358;  EV850359;  EV850360;  EV850361;  EV850362;  EV850363;  EV850364;  EV850365;  EV850366;  EV850367;  EV850368;  EV850369;  EV850370;  EV850371;  EV850372;  EV850373;  EV850374;  EV850375;  EV850376;  EV850377;  EV850378;  EV850379;  EV850380;  EV850381;  EV850382;  EV850383;  EV850384;  EV850385;  EV850386;  EV850387;  EV850388;  EV850389;  EV850390;  EV850391;  EV850392;  EV850393;  EV850394;  EV850395;  EV850396;  EV850397;  EV850398;  EV850399;  EV850400;  EV850401;  EV850402;  EV850403;  EV850404;  EV850405;  EV850406;  EV850407;  EV850408;  EV850409;  EV850410;  EV850411;  EV850412;  EV850413;  EV850414;  EV850415;  EV850416;  EV850417;  EV850418;  EV850419;  EV850420;  EV850421;  EV850422;  EV850423;  EV850424;  EV850425;  EV850426;  EV850427;  EV850428;  EV850429;  EV850430;  EV850431;  EV850432;  EV850433;  EV850434;  EV850435;  EV850436;  EV850437;  EV850438;  EV850439;  EV850440;  EV850441;  EV850442;  EV850443;  EV850444;  EV850445;  EV850446;  EV850447;  EV850448;  EV850449;  EV850450;  EV850451;  EV850452;  EV850453;  EV850454;  EV850455;  EV850456;  EV850457;  EV850458;  EV850459;  EV850460;  EV850461;  EV850462;  EV850463;  EV850464;  EV850465;  EV850466;  EV850467;  EV850468;  EV850469;  EV850470;  EV850471;  EV850472;  EV850473;  EV850474;  EV850475;  EV850476;  EV850477;  EV850478;  EV850479;  EV850480;  EV850481;  EV850482;  EV850483;  EV850484;  EV850485;  EV850486;  EV850487;  EV850488;  EV850489;  EV850490;  EV850491;  EV850492;  EV850493;  EV850494;  EV850495;  EV850496;  EV850497;  EV850498;  EV850499;  EV850500;  EV850501;  EV850502;  EV850503;  EV850504;  EV850505;  EV850506;  EV850507;  EV850508;  EV850509;  EV850510;  EV850511;  EV850512;  EV850513;  EV850514;  EV850515;  EV850516;  EV850517;  EV850518;  EV850519;  EV850520;  EV850521;  EV850522;  EV850523;  EV850524;  EV850525;  EV850526;  EV850527;  EV850528;  EV850529;  EV850530;  EV850531;  EV850532;  EV850533;  EV850534;  EV850535;  EV850536;  EV850537;  EV850538;  EV850539;  EV850540;  EV850541;  EV850542;  EV850543;  EV850544;  EV850545;  EV850546;  EV850547;  EV850548;  EV850549;  EV850550;  EV850551;  EV850552;  EV850553;  EV850554;  EV850555;  EV850556;  EV850557;  EV850558;  EV850559;  EV850560;  EV850561;  EV850562;  EV850563;  EV850564;  EV850565;  EV850566;  EV850567;  EV850568;  EV850569;  EV850570;  EV850571;  EV850572;  EV850573;  EV850574;  EV850575;  EV850576;  EV850577;  EV850578;  EV850579;  EV850580;  EV850581;  EV850582;  EV850583;  EV850584;  EV850585;  EV850586;  EV850587;  EV850588;  EV850589;  EV850590;  EV850591;  EV850592;  EV850593;  EV850594;  EV850595;  EV850596;  EV850597;  EV850598;  EV850599;  EV850600;  EV850601;  EV850602;  EV850603;  EV850604;  EV850605;  EV850606;  EV850607;  EV850608;  EV850609;  EV850610;  EV850611;  EV850612;  EV850613;  EV850614;  EV850615;  EV850616;  EV850617;  EV850618;  EV850619;  EV850620;  EV850621;  EV850622;  EV850623;  EV850624;  EV850625;  EV850626;  EV850627;  EV850628;  EV850629;  EV850630;  EV850631;  EV850632;  EV850633;  EV850634;  EV850635;  EV850636;  EV850637;  EV850638;  EV850639;  EV850640;  EV850641;  EV850642;  EV850643;  EV850644;  EV850645;  EV850646;  EV850647;  EV850648;  EV850649;  EV850650;  EV850651;  EV850652;  EV850653;  EV850654;  EV850655;  EV850656;  EV850657;  EV850658;  EV850659;  EV850660;  EV850661;  EV850662;  EV850663;  EV850664;  EV850665;  EV850666;  EV850667;  EV850668;  EV850669;  EV850670;  EV850671;  EV850672;  EV850673;  EV850674;  EV850675;  EV850676;  EV850677;  EV850678;  EV850679;  EV850680;  EV850681;  EV850682;  EV850683;  EV850684;  EV850685;  EV850686;  EV850687;  EV850688;  EV850689;  EV850690;  EV850691;  EV850692;  EV850693;  EV850694;  EV850695;  EV850696;  EV850697;  EV850698;  EV850699;  EV850700;  EV850701;  EV850702;  EV850703;  EV850704;  EV850705;  EV850706;  EV850707;  EV850708;  EV850709;  EV850710;  EV850711;  EV850712;  EV850713;  EV850714;  EV850715;  EV850716;  EV850717;  EV850718;  EV850719;  EV850720;  EV850721;  EV850722;  EV850723;  EV850724;  EV850725;  EV850726;  EV850727;  EV850728;  EV850729;  EV850730;  EV850731;  EV850732;  EV850733;  EV850734;  EV850735;  EV850736;  EV850737;  EV850738;  EV850739;  EV850740;  EV850741;  EV850742;  EV850743;  EV850744;  EV850745;  EV850746;  EV850747;  EV850748;  EV850749;  EV850750;  EV850751;  EV850752;  EV850753;  EV850754;  EV850755;  EV850756;  EV850757;  EV850758;  EV850759;  EV850760;  EV850761;  EV850762;  EV850763;  EV850764;  EV850765;  EV850766;  EV850767;  EV850768;  EV850769;  EV850770;  EV850771;  EV850772;  EV850773;  EV850774;  EV850775;  EV850776;  EV850777;  EV850778;  EV850779;  EV850780;  EV850781;  EV850782;  EV850783;  EV850784;  EV850785;  EV850786;  EV850787;  EV850788;  EV850789;  EV850790;  EV850791;  EV850792;  EV850793;  EV850794;  EV850795;  EV850796;  EV850797;  EV850798;  EV850799;  EV850800;  EV850801;  EV850802;  EV850803;  EV850804;  EV850805;  EV850806;  EV850807;  EV850808;  EV850809;  EV850810;  EV850811;  EV850812;  EV850813;  EV850814;  EV850815;  EV850816;  EV850817;  EV850818;  EV850819;  EV850820;  EV850821;  EV850822;  EV850823;  EV850824;  EV850825;  EV850826;  EV850827;  EV850828;  EV850829;  EV850830;  EV850831;  EV850832;  EV850833;  EV850834;  EV850835;  EV850836;  EV850837;  EV850838;  EV850839;  EV850840;  EV850841;  EV850842;  EV850843;  EV850844;  EV850845;  EV850846;  EV850847;  EV850848;  EV850849;  EV850850;  EV850851;  EV850852;  EV850853;  EV850854;  EV850855;  EV850856;  EV850857;  EV850858;  EV850859;  EV850860;  EV850861;  EV850862;  EV850863;  EV850864;  EV850865;  EV850866;  EV850867;  EV850868;  EV850869;  EV850870;  EV850871;  EV850872;  EV850873;  EV850874;  EV850875;  EV850876;  EV850877;  EV850878;  EV850879;  EV850880;  EV850881;  EV850882;  EV850883;  EV850884;  EV850885;  EV850886;  EV850887;  EV850888;  EV850889;  EV850890;  EV850891;  EV850892;  EV850893;  EV850894;  EV850895;  EV850896;  EV850897;  EV850898;  EV850899;  EV850900;  EV850901;  EV850902;  EV850903;  EV850904;  EV850905;  EV850906;  EV850907;  EV850908;  EV850909;  EV850910;  EV850911;  EV850912;  EV850913;  EV850914;  EV850915;  EV850916;  EV850917;  EV850918;  EV850919;  EV850920;  EV850921;  EV850922;  EV850923;  EV850924;  EV850925
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Base Sequence
Cattle Diseases / parasitology*
Computational Biology / methods
Dictyocaulus / genetics*,  growth & development
Dictyocaulus Infections / parasitology*
Expressed Sequence Tags
Gene Expression Profiling / methods*
Gene Expression Regulation, Developmental
Molecular Sequence Data
Sequence Analysis, DNA
Time Factors

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-ta): BMC Genomics
ISSN: 1471-2164
Publisher: BioMed Central
Article Information
Download PDF
Copyright ? 2007 Ranganathan et al; licensee BioMed Central Ltd.
open-access: This is an Open Access article distributed under the terms of the Creative Commons Attribution License (), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Received Day: 17 Month: 4 Year: 2007
Accepted Day: 5 Month: 9 Year: 2007
collection publication date: Year: 2007
Electronic publication date: Day: 5 Month: 9 Year: 2007
Volume: 8First Page: 311 Last Page: 311
ID: 2131760
Publisher Id: 1471-2164-8-311
PubMed Id: 17784965
DOI: 10.1186/1471-2164-8-311

A transcriptomic analysis of the adult stage of the bovine lungworm, Dictyocaulus viviparus
Shoba Ranganathan12 Email:
Shivashankar H Nagaraj1 Email:
Min Hu3 Email:
Christina Strube4 Email:
Thomas Schnieder4 Email:
Robin B Gasser3 Email:
1Department of Chemistry and Biomolecular Sciences, Macquarie University, Sydney, New South Wales 2109, Australia
2Biotechnology Research Institute, Macquarie University, Sydney, New South Wales 2109, Australia
3Department of Veterinary Science, The University of Melbourne, 250 Princes Highway, Werribee, Victoria 3030, Australia
4Institute for Parasitology, University of Veterinary Medicine Hannover, Buenteweg 17, D-30559 Hannover, Germany


Parasitic nematodes of livestock cause substantial economic losses due to poor productivity, failure to thrive and deaths [1,2]. The financial losses associated with these endoparasites are estimated at billions of dollars per annum [3]. Lungworms of the genus Dictyocaulus (family Dictyocaulidae) are key parasitic nematodes which cause pathological effects and clinical disease in different ruminant hosts, particularly in young animals [4,5]. Dictyocaulus viviparus, the bovine lungworm, causes a severe and frequently fatal bronchitis (known colloquially as 'husk') which is of major importance in many countries [6]. Severe cases of dictyocaulosis lead to emphysema and pneumonia ? heavy infections can cause a mortality rate of >20% among affected cattle [1,2].

Dictyocaulus viviparus has a direct life cycle [7]. The adult stages (females and males) live in the bronchi, where the ovoviviparous females produce eggs from which first-stage larvae (L1) usually hatch rapidly whilst in the lung or the intestinal tract. The L1s are then shed in the faeces of the bovine host. Under favourable environmental conditions, L1s develop through to the infective third-stage larvae (L3s) during a period of ~4?6 days. After ingestion by the host, L3s migrate through the gut wall to the mesenteric lymph nodes, moult, and, as fourth-stage larvae (L4s), are transported to the lungs. L4s penetrate the alveoli, moult and then develop into adults. However, larval stages can remain inhibited in the lungs for up to 5 months. In cattle, the period from ingestion of L3s to reproductive maturity of the adult worms is 3?4 weeks.

While there is considerable knowledge of the morphological changes taking place during the life cycle of D. viviparus, very little is understood about the fundamental molecular and biochemical processes underlying the development and survival of this parasite and the parasite-host interplay. Insights into such processes are fundamentally important and could provide a basis for the identification of molecular targets for the rational design of nematocidal compounds, vaccines or/and for diagnosis. To date, studies of D. viviparus have been limited to individual genes and proteins. For instance, before the present study (March 2007), 323 gene sequences, 221 protein sequences and 229 research articles relating to D. viviparus were available in public databases.

Current technological advances in genomics provide exciting opportunities for exploring basic molecular biological and biochemical aspects of D. viviparus and related nematodes. For instance, expressed sequence tag (EST) data sets facilitate the prediction and categorization of key molecules, particularly those linked to development (both in pre-parasitic and parasitic stages), sexual differentiation and maturation, based on comparisons with other organisms for which sequence and functional genomic data sets are available. Also, the complete genome sequence of the free-living nematode Caenorhabditis elegans and the wealth of information on gene expression and function for this nematode [8,9] provide a means of evaluating homologues and orthologues [10], since D. viviparus cannot be maintained or propagated effectively in vitro for the functional testing of genes and gene products. Also, the potential of gene silencing techniques [11,12] provides a prospect for the functional analysis of molecules in this and other parasitic nematodes.

In the present study, we provide a first insight into the transcriptome of the adult stage of D. viviparus via EST sequencing and apply a newly established computational platform [13] for the clustering and comparative analyses of the data set against data available for a range of organisms, with an emphasis on the best characterized nematode, C. elegans. The representative ESTs from this dataset have been annotated functionally at the gene and protein levels to aid in assigning, in the main, gene ontologies, protein families and biochemical pathways. Such annotation techniques have enabled us to pin-point genes that could be considered in the development of intervention strategies. The present data provide a foundation for future investigations in areas, such as the stage-, sex- and tissue-specific gene transcription or expression, whole genome sequencing and proteomics of D. viviparus.

Results and discussion
General EST analysis

Of 5271 clones sequenced, a total of 4436 quality ESTs were obtained (Phase I, Figure 1), achieving a sequencing success of 84% (Table 1), which is consistent with previous studies [14,15]. These pre-processed ESTs ranged from 80?1164 bp, with a mean of 730 bp and a standard deviation (S.D.) of 258 bp. After clustering, the mean length of the contigs (or consensus sequences) increased to 787 (+/- 313) bp. The G+C content of the coding sequences was 43.5%, consistent with other nematodes from clade V [16,17] and slightly more than C. elegans (37%) and its congener, C. briggsae (38%) [18]. Under the assumption that the D. viviparus genome codes for ~22,000 proteins [19], the EST clusters were predicted to represent ~10?15% of the proteins encoded by this genome.

The cluster analysis of the 4436 ESTs from adult D. viviparus yielded 2258 representative ESTs (rESTs; 458 contig and 1800 singleton sequences; see Table 1), of which 1685 (74.6%) had open reading frames (ORFs). All rESTs were then subjected to analyses using ESTExplorer [13], a semi-automated bioinformatics pipeline (Fig. 1). Also, we queried rESTs (at the amino acid level) against three databases containing protein sequences from different organisms, in order to categorize the molecules from D. viviparus. Data were compared with protein sequences available for (i) C. elegans (from WORMPEP v.167 [20]), (ii) parasitic nematodes (available protein sequences and peptides from conceptually translated ESTs), and (iii) organisms other than nematodes (from NCBI non-redundant protein database) [21]. Three-way comparison of D. viviparus rESTs with homologues from C. elegans, WORMPEP and parasitic nematodes have been figuratively presented using SimiTri [22] (Figure 2). For this comparison, similarity searches of the 2258 representative sequences resulted in 1159 (51.3%) homologues to C. elegans, 1174 (51.9%) to those from other parasitic nematodes (including some strongylids), 827 (36.6%) homologues in organisms other than nematodes, and 39.5% had no significant similarity to any other organism (employing a cut-off of 1e-05) for which sequence data are currently available (Figure 2, Additional File 1). The SimiTri plot (Figure 2) shows that, given the current database contents, the sampled transcriptome from D. viviparus is equally close to available C. elegans and parasitic nematode sequences, compared with non-nematodes.

Comparative analyses with C. elegans data sets

The comparative analysis to identify homologues in C. elegans is important because D. viviparus and this free-living nematode are both considered to belong to clade V of the Nematoda [16,17], and because C. elegans also represents the best characterized nematode in many respects, particularly in terms of its genome, genetics, biology, physiology, biochemistry, as well as the localization and functions of molecules [20,23]. Specifically, the comparative analysis (at the amino acid level) of all rESTs with C. elegans data (see Additional File 1) revealed 1159 (51.3%) key, well-characterized molecules associated with various biological processes (n = 540), including development, regulation of biological processes, response to abiotic and biotic stimuli and reproduction. 'Non-wildtype' RNAi phenotypes in C. elegans (such as embryonic lethality, maternal sterility, sterility in progeny, larval arrest and slow growth) were associated with 569 (48.7%) of these 1159 molecules (Table 3). Of the 1159 C. elegans homologues, the functions for 776 (66%) rESTs could be inferred using Gene Ontologies (GO) [24], with 696 (41%) sequences being mapped to key biological pathways (including signal transduction mechanisms, antigen processing and presentation, regulation of actin cytoskeleton, ribosomal proteins and translation factors). Overall, the functional classification revealed that approximately half of the rESTs had homologues in C. elegans and parasitic nematodes, one tenth were specific to parasitic nematodes, and one third of the rESTs did not match any sequence in current databases, possibly representing novel genes.

As multiple ESTs can be derived from the same gene, it was important to predict how many unique genes were represented by the rESTs. Mapping the 2258 D. viviparus rESTs to C. elegans revealed that, of the 1159 D. viviparus rESTs with similarities to 927 C. elegans genes, the majority of these (798/1159 or 68.8%) had a one-to-one relationship to their C. elegans homologue. The remaining 369 rESTs mapped to multiple non-overlapping regions from 129 C. elegans genes (with an estimated fragmentation rate of 25%). After discounting for fragmentation, we estimated that 1694 unique D. viviparus genes were identified, with a suggested new gene discovery rate of 38.2% (1694/2258).

C. elegans (non-wild-type) RNAi phenotypes can provide some indication of the relevance and functions of orthologous genes in other nematodes, particularly in parasitic nematodes of clade V, for which the complexity of an obligate parasitic life cycle and the lack of an effective (long-term) laboratory culture system make high-throughput functional screening impractical [25]. We retrieved C. elegans RNAi data representing D. viviparus homologues. Of 1159 D. viviparus rESTs, 569 had homologues in C. elegans which could be silenced by RNAi (Additional File 1). The RNAi phenotypes (as described by Wormbase) included Adl (adult lethal), Age (ageing alteration), Bmd (body morphology defect), Dpy (dumpy), Egl (egg laying defect), Emb (embryonic lethal), Gro (slow growth), Let (larval lethal), Lvl (larval lethal), Lva (larval arrest) and Unc (uncoordinated), whereas 590 homologues had no observable RNAi phenotype in C. elegans. We also found that 23% of the most abundant (Table 2) and 22% of the transcripts predicted to represent secreted proteins of D. viviparus (Additional File 2) had C. elegans homologues with non-wildtype RNAi phenotypes.

Comparative analysis with data for other nematodes within clade V

D. viviparus is a member of clade V of the phylum Nematoda [16,17], which comprises members of the orders Strongylida, Rhabditida and Diplogasterida. We conducted a comparative analysis of all 2258 rESTs with data [17,26] for various nematodes (including Ancylostoma caninum, Ancylostoma ceylanicum, Necator americanus, Nippostrongylus brasiliensis, Haemonchus contortus, Ostertagia ostertagi, Teladorsagia circumcincta and Pristionchus pacificus) belonging to clade V (Additional File 3), to examine gene conservation within this clade. The analysis revealed 601 (26.6%) rESTs to have significant similarity to molecules from the members of clade V. These rESTs represented house-keeping (including cathepsin B-like cysteine proteases 1 and 2 and serine-threonine protein kinase) as well as nematode-specific (vitellogenin structural and major sperm protein) genes. We could assign GO terms to 461 of these 601 rEST sequences, associated with 118 different biological processes, such as embryonic development, intracellular protein transport, protein metabolism and responses to abiotic and/or biotic stimuli. Furthermore, 404 sequences could be mapped to biological pathways predicted to be associated with ribosomal proteins and the proteasome system, including 23 predicted secreted proteins mapping to cysteine proteinases, secreted protein 5 precursor and parasite pepsinogen.

Abundant transcripts in adult D. viviparus

A high level of representation in a cDNA library usually correlates with high transcript abundance in the original biological sample [27], although artefacts of library construction can result in a selection for or against representation of some transcripts. The D. viviparus clusters were ranked according to the number of contributing ESTs, and the top 30 clusters, which represented 1261 (24%) of the total number of rESTs (2258) obtained, were investigated in detail (Table 2). A number of clusters had significant alignments to known proteins, the majority of which were house-keeping genes, such as elongation factors, ribosomal proteins, aldolases, kinases, proteases and actins. Some of these genes have been identified in a number of other nematodes, including Ancylostoma caninum, Ancylostoma ceylanicum, Dirofilaria immitis, Strongyloides ratti and Meloidogyne incognita [28-32], requiring detailed characterisation to understand their functions in D. viviparus.

We analysed the most abundantly expressed transcripts from D. viviparus and found that 12 of the 30 (40%) contigs had no significant similarity to any sequence in the non-redundant protein database. As most of the nematode data are available only as ESTs and not included in the BLAST databases, we further compared these 12 contigs with sequences against Parasite Genomes using WU-Blast2 and BLASTN against the NCBI 'other-ESTs' database. We found that five sequences did not match any sequence in the databases, whereas seven entries had similarity to ESTs from other nematodes. The other 18 (60%) contigs were assigned functionality based on BLASTP against the NR database, and all of them had homologues in either non-parasitic (C. elegans and/or C. briggsae) or parasitic nematodes. A summary of these findings is provided in Table 2. Most of the homologues were found to be house-keeping or structural genes (including aldolase, histone family members, lectin, collagen and cytochrome oxidase), and four contigs were represented by molecules specific to nematodes, such as the major sperm proteins.

Comparison with cDNAs from third-stage larvae of D. viviparus

Recently, Strube et al. [33] identified and characterized 28 cDNAs differentially transcribed between experimentally induced hypobiotic and infective third-stage larvae (L3) of D. viviparus using a suppressive-subtractive hybridization (SSH) approach. We compared these 28 sequences against our dataset of 2258 rESTs, using BLASTN, to identify whether any of them were represented in the adult stage, with only one match. The sole sequence common to both datasets was L3ni 18 (accession number EG374523), which matched EST D.viviparus_42_A11 in the present dataset, being a homologue of the C. elegans hypothetical protein C05D11.10 (mitochondrial/chloroplast ribosomal S17-like protein [code KOG3447]). Thus, this EST study of adult D. viviparus represents a novel dataset for parasitic nematodes in clade V, representing the Metastrongyloidea (cf. [47]).

Functional classification of rESTs from adult D. viviparus

We annotated all 2258 rESTs systematically using a range of bioinformatic tools (details available in Figure 1). This annotation included functional classifications of rESTs using Gene Ontologies (GO) [24], pathway mapping using KEGG (Kyoto Encyclopedia of Genes and Genomes) [34], visualisation of EST data comparisons using SimiTri [22] and analyses of the D. viviparus secretome using SignalP [35], TMHMM [36] and PSORTII [37]. Results of these analyses are described in the following two sections

a. Gene Ontologies

Gene Ontology (GO) has been used widely to predict gene function and classification. GO provides a dynamic vocabulary and hierarchy that unifies descriptions of biological, cellular and molecular functions across genomes. We used BLAST2GO [38], a sequence-based tool to assign GO terms, extracting them for each BLAST hit obtained by mapping to extant annotation associations. We found that 776 (31%) of 2258 rESTs could be functionally assigned to biological processes (n = 540), cellular components (n = 328) and molecular functions (n = 457). A summary GO representation (using GO Slim) of the D. viviparus rESTs is given in Table 3.

Amongst the most common GO categories representing biological processes were: binding (GO: 0005488), catalytic activity (GO: 0003824) and structural molecule activity (GO: 0005198); development (GO: 0007275), metabolism (GO: 0008152), reproduction (GO:0000003) and growth (GO:0040007). The largest number of GO terms in cellular components was for cell part (GO:0044464), membrane-bound organelle (GO:0043227) and non-membrane-bound organelle (GO:0043228). A complete listing of GO mappings assigned for rESTs is provided in Additional File 4.

b. Pathway analysis using KEGG assignments

Biochemical functionality was predicted by mapping all 2258 rESTs to pathways, using KOBAS implemented within ESTExplorer [13], with an E-value cut-off of 1.0e-5. Enzyme commission (EC) numbers were used to appraise which sequences pertained to a specific pathway. A total of 696 (31%) sequences were mapped to 139 KEGG pathways, with 453 sequences representing metabolic enzymes characterized by unique EC numbers. The top 30 (highly represented) pathways are shown in Table 4.

Molecules involved in signal transduction mechanisms (n= 26) and glycolysis/gluconeogenesis (h= 24) had the highest representation amongst the sequences mapped to KEGG pathways. We also identified 50 predicted proteins with potential roles in host-parasite interactions, with 11 molecules predicted to be involved in antigen processing and/or presentation, six in T-cell receptor signalling pathway and seven CD molecules. Although, at this stage, the precise role of such molecules in the parasite-host interplay is unclear, they could be involved in manipulating or evading the host's immune response(s) or associated with the parasite's innate immune response. Extensive experimental work would be required to test these proposals. Furthermore, we identified families of proteins representing serine, cysteine and metallo-proteinases as well as proteinase inhibitors (e.g., cystatins). While these enzymes are inferred to mediate or modulate proteolytic functions, which, in turn, may facilitate tissue migration and other interactions with host cells, the proteinase inhibitors may protect the parasite against digestion by host or endogenous proteinases excreted/secreted in the lung of the mammalian host [39]. A complete listing of the KEGG mappings is available as supplementary data (Additional File 5).

Secretome analysis

An important starting point in the identification of potential novel drug or vaccine candidates in parasites is the prediction of molecules that are secreted or excreted in or around the host- parasite interface [40-42]. Examples of such proteins are the aspartyl protease inhibitor (API-1) [43], mi-msp-1 (similar to a venom allergen antigen AG5-like) [44] and the Ancylostoma-secreted protein (ASP) [45]. In the present data set (= 2258 rESTs), we identified 85 putatively secreted proteins representing a non-redundant catalogue of D. viviparus molecules (Additional 2). Of these, 26 (30.5%) sequences had no significant similarity to any sequence available in current databases, whereas 59 (69.4%) had homologues in nematodes, with 43 (50.5%) C. elegans and/or C. briggsae matches, and 16 (18.8%) homologues in various other parasitic nematodes, including the blood-feeding nematodes Ancylostoma ceylanicum, Necator americanus and Haemonchus contortus.

The secretome analysis (Additional File 2) revealed a number of unique features. Firstly, seven of the putative secreted protein entries were homologous to either sperm-specific family members [41] or major sperm proteins (MSP), consistent with the mass spectrometric analysis of secreted molecules from D. viviparus [41]. Secondly, acetylcholinesterase (AChE) has been identified as an important enzyme secreted by adult D. viviparus, thought to be involved in parasite survival in the host and, therefore, being a vaccine candidate [46]. We now report an entry for secreted acetylcholinesterase (D.viviparus_8_A01 of 332 amino acids, with a signal peptide length of 21) in the secretome analysis of the current EST dataset (Additional File 2), which may also represent a target for further characterization.


The present study has given us a first glimpse of the transcriptome of the adult stage of the bovine lungworm, D. viviparus, and represents a starting point for studies in a number of different fundamental and applied areas. We used a comprehensive EST analysis pipeline, ESTExplorer, for this purpose, for functional annotation at the DNA and protein levels [13]. From this single study of 5271 ESTs, we have identified 55 novel sequences (1.4%), with high confidence, with no known homologue in any other nematode or mammal for which sequence data are presently available in public databases. These molecules are particularly interesting, as they may represent genes that may be specific to parasitism or to the species. However, such molecules are very challenging to work on, as their potential functions cannot be predicted using current bioinformatic approaches. However, there is considerable scope in exploring such molecules in the future, using a combination of genomic and proteomic approaches. Insights into such molecules and/or their interaction with the bovine host could provide opportunities for developing novel intervention approaches.

From a systematic viewpoint, D. viviparus, belongs to the Metastrongyloidea ("metastrongyles" or "lungworms") based on nuclear ribosomal DNA sequence data [47], as distinct from the Trichostrongyloidea (mainly in the stomach and small intestine), Strongyloidea (mostly in the large intestine) and Ancylostomatoidea (small intestine), and thus represents, from biological, host-parasite relationship and molecular evolutionary perspectives, a very interesting species for comparative genomic analysis with nematodes from these superfamilies. Therefore, this nematode brings a number of important benefits for future investigations, particularly for genome sequencing and for subsequent comparative evolutionary analyses. Indeed, D. viviparus, among other strongylid nematodes, has recently been selected for whole genome sequencing, to be carried out at the Genome Sequencing Center of Washington University in St Louis, USA [48]. For D. viviparus, the genomic information from the present study underpins future microarray analyses, focused on exploring the transcriptional profiles among different stages (e.g., larval versus adult stages; hypobiotic stages versus those which are not arrested; free-living (L1s and L2s) versus infective versus late larval stages), sexes (female versus male) and tissues (e.g., germline versus neural versus musculature versus intestine) of the parasite. Such studies, particularly those of the molecules differentially transcribed and expressed during the transition to parasitism, the invasion of the host and hypobiosis, could provide unique insights into such key molecular developmental and reproductive processes. While there is some controversy regarding the applicability and usefulness of RNAi to some parasitic nematodes, such as the Strongylida [11,12], comparative studies of gene-silencing and transgenesis in C. elegans are considered useful for exploring the function and regulation of some relatively conserved parasite genes, provided data are interpreted with caution [10]. This is particularly the case with the continued increase in genome sequence information.

With the future availability of whole genome sequence data for D. viviparus, it will also be possible to carry out meaningful mass spectroscopic analyses of differentially expressed proteins [49,50], allowing large-scale analysis of proteins from small amounts of parasite material. Such analyses will enable the link to be made between the regulation of transcription and translation and, importantly, in the study of parasites, will allow the analysis of proteins expressed within short time frames within or external to the host animal, or within organs or micro-environments within the parasite [51,52]. Hence, the application of an integrated bioinformatic-genomic-phenomic-proteomic ("systems biology") approach, focusing on developmental processes and mechanisms, could enhance our understanding of the molecular biology of moulting, invasion of and establishment in the host, hypobiosis (arrested development), and sexual differentiation, maturation and behaviour of D. viviparus. Clearly, progress in such fundamental areas could lead to the development of exciting new ways of treating, controlling or preventing this lungworm and other parasitic nematodes, by blocking or disrupting key biological pathways in them.

Parasite material

Adults of D. viviparus (strain HannoverDv2000) were produced in helminth-free male Holstein-Friesian calves (five months of age). Four weeks after oral inoculation with 3300 L3, the calves with patent D. viviparus infection were euthanized and the worms collected from the lungs as described by Wood et al. [53]. The worms were washed extensively (five times) in large volumes (100 ml) of fresh diethyl-pyrocarbonate (DEPC)-treated saline (22?C), transferred to sterile, RNase-free screw-top cryovials? (4 ml; Roth, Karlsruhe, Germany) and frozen as 200 ?l pellets in a minimal amount of saline at -75?C or -196?C until RNA isolation.

Isolation of total RNA, cDNA synthesis and cDNA library construction

Total RNA was extracted from the adults of D. viviparus (under liquid nitrogen, employing a sterile mortar and pestle) using the TriPure isolation reagent? (Roche Molecular Biochemicals). Integrity and yields of RNA were verified and estimated, respectively, using the Bioanalyzer 2100 (Agilent). Each RNA sample (~10 ?g) was treated with 2 U of DNase I (Promega), incubated at 37?C for 30 min prior to heat denaturation of the enzyme (75?C for 5 min) and then stored at -70?C until use. The cDNA was produced from 10 ?g of total RNA from D. viviparus, using an oligo (dT) primer and Superscript II reverse transcriptase (Invitrogen) and then purified over DyeEx columns (Qiagen).

A non-directional cDNA library was constructed in the plasmid vector pGEM-T (Promega) by T-ended cloning, according to the manufacturer's protocol. Colonies were screened using blue-white selection. Clones (n = 5271) were picked randomly and patched on to grided Luria Bertani (LB) agarose plates containing 100 mg/ml ampicillin. Single-pass sequencing (using the T7 primer) was performed employing BigDye Chemistry (v3.1) in a 3730 ? l DNA analyser (Applied Biosystems).

EST analysis

The ESTs were initially analysed and annotated using ESTExplorer, an automated EST analysis platform [13,54]. In brief, the analyses comprised three phases (see Figure 1). In phase I, all ESTs were pre-processed (SeqClean, RepeatMasker), aligned/clustered using the Contig Assembly Program CAP3, employing a minimum sequence overlap length "cut-off" of 30 bases and an identity threshold of 95% (in Phase I) for the removal of flanking vector and adapter sequences, followed by assembly. Phase II of the ESTExplorer led to gene ontology inference, at the DNA-level annotation, using BLAST2GO (V 1.6.2) [38], using Gene Ontology (MySQL-DB-data release go_200609). In Phase III, rESTs were then conceptually translated into peptides (using ESTScan), which were further characterized via InterProScan (domain/motifs), and peptides mapped to respective pathways in C. elegans using KOBAS (KEGG Orthology-Based Annotation System, KEGG data release 40.0). We also retrieved KO data from KEGG and then classified enzymes from non-enzymes based on our pathway mapping results. Peptides predicted from rEST were also compared, using BLASTP, with the non-redundant protein sequence database from National Centre for Biotechnology Information (NCBI), as part of the generic ESTExplorer pipeline for systematic EST analysis and annotation.

Protein databases for 'parasitic nematodes' and 'non-nematodes' were generated in-house for similarity searches. The 'parasitic nematodes' group contains all available protein sequences for parasitic nematodes and ESTs from GenBank (17 February 2007), translated into peptide sequences, whereas the 'non-nematodes' database comprises amino acid sequences from the complete non-redundant protein database NR (17 February 2007) excluding those from nematodes. Homologues to rESTs were identified via comparisons against WormBase using BLASTX and the Parasite genome WU-BLAST2 Nematoda database (from the European Bioinformatics Institute) using BLASTN. Each EST of D. viviparus was assigned a 'statistically significant' gene homologue if the E-value from the BLAST output of the sequence alignment was <1e-05. The program SimiTri [22] was used for the comparison (at the amino acid sequence level) of D. viviparus rESTs with data in C. elegans, parasitic nematode and non-nematode protein sequence databases. SimiTri provides a two-dimensional display of relative similarity relationships among three different datasets.

From the peptides inferred from rESTs, secreted proteins were predicted using a combination of three programs, to minimize the number of false positive predictions. Firstly, SignalP 3.0 [35] was used to predict the presence of secretory signal peptides and signal anchors for each predicted rEST proteins. A signal sequence was considered present when it was predicted both by the artificial neural network and the hidden Markov model prediction approaches (SignalPNN and SignalP-HMM, available as options within SignalP). In order to exclude the erroneous prediction of putative transmembrane (TM) sequences as signal sequences, TMHMM [36], a membrane topology prediction program, was then applied. We further validated the list of secreted proteins, using extracellular localization, employing PSORT [37].

Note: The final set of quality ESTs reported in this paper 15 are available in the EMBL, GenBank and DDJB databases under accession numbers EV849926 ? EV854361.

Authors' contributions

SR, SHN and RBG conceived and designed the research plan and participated in all aspects of data collection and analysis. SHN and SR analysed, and SHN, SR, MH and RBG interpreted the data. CS and TS provided adult worms of D. viviparus. RBG and MH isolated and purified RNA for cDNA library construction and supervised and coordinated the sequencing. All authors contributed to the analyses and the writing of the manuscript. All authors read and approved the final manuscript.

Supplementary Material Additional file 1

Comparison of 2258 rESTs from Dictyocaulus viviparus with Caenorhabditis elegans data (Wormpep v 167). The table also provides corresponding RNAi phenotypes in C. elegans and information on homologues in parasitic nematodes.

Click here for additional data file (1471-2164-8-311-S1.xls)

Additional file 2

Secreted proteins predicted from rESTs from Dictyocaulus viviparus.

Click here for additional data file (1471-2164-8-311-S2.pdf)

Additional file 3

Dictyocaulus viviparus homologues in nematodes in clade V, including Ancylostoma caninum, Ancylostoma ceylanicum, Necator americanus, Nippostrongylus brasiliensis, Haemonchus contortus, Ostertagia ostertagi, Teladorsagia circumcincta and Pristionchus pacificus.

Additional file 4

Mapped metabolic pathways in adult Dictyocaulus viviparus based on Kyoto Encyclopedia of Genes and Genomes (KEGG). The pathway mapping was carried out using the program KOBAS.

Click here for additional data file (1471-2164-8-311-S4.xls)

Additional file 5

Gene Ontology mappings (using GO slim terms) for Dictyocaulus viviparus clusters generated using BLAST2GO program. Note that individual GO categories can have multiple mappings.

Click here for additional data file (1471-2164-8-311-S5.xls)


This work was supported by grants from the Australian Research Council (LP0667795 and DP0665230), Genetic Technologies Limited and Meat and Livestock Australia. Thanks to Dr Gary Cobon, Stephen Page and Joan Lloyd for continued support. SHN is grateful to Macquarie University for the award of iMURS research scholarships and an MUPGR travel grant. Funding to pay the Open Access publication charges for this article was provided by Macquarie University.

Coles GC. The future of veterinary parasitologyVet Parasitol 2001;98:31–39. [pmid: 11516577] [doi: 10.1016/S0304-4017(01)00421-6]
Panuska C. Lungworms of ruminantsVet Clin North Am Food Anim Pract 2006;22:583–593. [pmid: 17071354] [doi: 10.1016/j.cvfa.2006.06.002]
Newton SE,Munn EA. The development of vaccines against gastrointestinal nematode parasites, particularly Haemonchus contortusParasitol Today 1999;15:116–122. [pmid: 10322325] [doi: 10.1016/S0169-4758(99)01399-X]
David GP. Survey on lungworm in adult cattleVet Rec 1997;141:343–344. [pmid: 9347626]
Eysker M,van Miltenburg L. Epidemiological patterns of gastrointestinal and lung helminth infections in grazing calves in The NetherlandsVet Parasitol 1988;29:29–39. [pmid: 2972109] [doi: 10.1016/0304-4017(88)90005-2]
Corwin RM. Economics of gastrointestinal parasitism of cattleVet Parasitol 1997;72:451–457. [pmid: 9460211] [doi: 10.1016/S0304-4017(97)00110-6]
Anderson RC. Nematode Parasites of Vertebrates: Their Development and Transmission. 2001Wallingford: CABI Publishing;
Gilleard JS. The use of Caenorhabditis elegans in parasitic nematode researchParasitology 2004;128:S49–S70. [pmid: 16454899]
Gilleard JS,Woods DJ,Dow JA. Model-organism genomics in veterinary parasite drug-discoveryTrends Parasitol 2005;21:302–305. [pmid: 15923143] [doi: 10.1016/]
Britton C,Murray L. Using Caenorhabditis elegans for functional analysis of genes of parasitic nematodesInt J Parasitol 2006;36:651–659. [pmid: 16616144] [doi: 10.1016/j.ijpara.2006.02.010]
Zawadzki JL,Presidente PJ,Meeusen EN,De Veer MJ. RNAi in Haemonchus contortus: a potential method for target validationTrends Parasitol 2006;22:495–499. [pmid: 16971180] [doi: 10.1016/]
Geldhof P,Murray L,Couthier A,Gilleard JS,McLauchlan G,Knox DP,Britton C. Testing the efficacy of RNA interference in Haemonchus contortusInt J Parasitol 2006;36:801–810. [pmid: 16469321] [doi: 10.1016/j.ijpara.2005.12.004]
Nagaraj S,Deshpande N,Gasser R,Ranganathan S. ESTExplorer: an expressed sequence tag (EST) assembly and annotation platformNucleic Acids Res 2007;35:W143–W147. [pmid: 17545197] [doi: 10.1093/nar/gkm378]
Cottee PA,Nisbet AJ,Abs El-Osta YG,Webster TL,Gasser RB. Construction of gender-enriched cDNA archives for adult Oesophagostomum dentatum by suppressive-subtractive hybridization and a microarray analysis of expressed sequence tagsParasitology 2006;132:691–708. [pmid: 16426483] [doi: 10.1017/S0031182005009728]
Nisbet AJ,Gasser RB. Profiling of gender-specific gene expression for Trichostrongylus vitrinus (Nematoda: Strongylida) by microarray analysis of expressed sequence tag libraries constructed by suppressive-subtractive hybridisationInt J Parasitol 2004;34:633–643. [pmid: 15064128] [doi: 10.1016/j.ijpara.2003.12.007]
Blaxter M. Caenorhabditis elegans is a nematodeScience 1998;282:2041–2046. [pmid: 9851921] [doi: 10.1126/science.282.5396.2041]
Parkinson J,Mitreva M,Whitton C,Thomson M,Daub J,Martin J,Schmid R,Hall N,Barrell B,Waterston RH,McCarter JP,Blaxter ML. A transcriptomic analysis of the phylum NematodaNat Genet 2004;36:1259–1267. [pmid: 15543149] [doi: 10.1038/ng1472]
Stein LD,Bao Z,Blasiar D,Blumenthal T,Brent MR,Chen N,Chinwalla A,Clarke L,Clee C,Coghlan A,Coulson A,D'Eustachio P,Fitch DH,Fulton LA,Fulton RE,Griffiths-Jones S,Harris TW,Hillier LW,Kamath R,Kuwabara PE,Mardis ER,Marra MA,Miner TL,Minx P,Mullikin JC,Plumb RW,Rogers J,Schein JE,Sohrmann M,Spieth J,et al. The genome sequence of Caenorhabditis briggsae: a platform for comparative genomicsPLoS Biol 2003;1:E45. [pmid: 14624247] [doi: 10.1371/journal.pbio.0000045]
Consortium. CeSGenome sequence of the nematode C. elegans: a platform for investigating biologyScience 1998;282:2012–2018. [pmid: 9851916] [doi: 10.1126/science.282.5396.2012]
Benson DA,Karsch-Mizrachi I,Lipman DJ,Ostell J,Wheeler DL. GenBankNucleic Acids Res 2006;34:D16–D20. [pmid: 16381837] [doi: 10.1093/nar/gkj157]
Parkinson J,Blaxter M. SimiTri ? visualizing similarity relationships for groups of sequencesBioinformatics 2003;19:390–395. [pmid: 12584125] [doi: 10.1093/bioinformatics/btf870]
Ashburner M,Ball CA,Blake JA,Botstein D,Butler H,Cherry JM,Davis AP,Dolinski K,Dwight SS,Eppig JT,Harris MA,Hill DP,Issel-Tarver L,Kasarskis A,Lewis S,Matese JC,Richardson JE,Ringwald M,Rubin GM,Sherlock G. Gene ontology: tool for the unification of biology. The Gene Ontology ConsortiumNat Genet 2000;25:25–29. [pmid: 10802651] [doi: 10.1038/75556]
Geldhof P,Visser A,Clark D,Saunders G,Britton C,Gilleard J,Berriman M,Knox D. RNA interference in parasitic helminths: current situation, potential pitfalls and future prospectsParasitology 2007;134:609–619. [pmid: 17201997] [doi: 10.1017/S0031182006002071]
Mitreva M,Wendl MC,Martin J,Wylie T,Yin Y,Larson A,Parkinson J,Waterston RH,McCarter JP. Codon usage patterns in Nematoda: analysis based on over 25 million codons in thirty-two speciesGenome Biol 2006;7:R75. [pmid: 16907990] [doi: 10.1186/gb-2006-7-8-r75]
Audic S,Claverie JM. The significance of digital gene expression profilesGenome Res 1997;7:986–995. [pmid: 9331369]
Yin Y,Martin J,McCarter JP,Clifton SW,Wilson RK,Mitreva M. Identification and analysis of genes expressed in the adult filarial parasitic nematode Dirofilaria immitisInt J Parasitol 2006;36:829–839. [pmid: 16697384] [doi: 10.1016/j.ijpara.2006.03.002]
McCarter JP,Mitreva MD,Martin J,Dante M,Wylie T,Rao U,Pape D,Bowers Y,Theising B,Murphy CV,Kloek AP,Chiapelli BJ,Clifton SW,Bird DM,Waterston RH. Analysis and functional classification of transcripts from the nematode Meloidogyne incognitaGenome Biol 2003;4:R26. [pmid: 12702207] [doi: 10.1186/gb-2003-4-4-r26]
Thompson FJ,Mitreva M,Barker GL,Martin J,Waterston RH,McCarter JP,Viney ME. An expressed sequence tag analysis of the life-cycle of the parasitic nematode Strongyloides rattiMol Biochem Parasitol 2005;142:32–46. [pmid: 15907559] [doi: 10.1016/j.molbiopara.2005.03.006]
Geldhof P,Whitton C,Gregory WF,Blaxter M,Knox DP. Characterisation of the two most abundant genes in the Haemonchus contortus expressed sequence tag datasetInt J Parasitol 2005;35:513–522. [pmid: 15826643] [doi: 10.1016/j.ijpara.2005.02.009]
Mitreva M,McCarter JP,Arasu P,Hawdon J,Martin J,Dante M,Wylie T,Xu J,Stajich JE,Kapulkin W,Clifton SW,Waterston RH,Wilson RK. Investigating hookworm genomes by comparative analysis of two Ancylostoma speciesBMC Genomics 2005;6:58. [pmid: 15854223] [doi: 10.1186/1471-2164-6-58]
Strube C,Schnieder T,von Samson-Himmelstjerna G. Differential gene expression in hypobiosis-induced and non-induced third-stage larvae of the bovine lungworm Dictyocaulus viviparusInt J Parasitol 2007;37:221–231. [pmid: 17112525] [doi: 10.1016/j.ijpara.2006.09.014]
Kanehisa M,Goto S,Hattori M,Aoki-Kinoshita KF,Itoh M,Kawashima S,Katayama T,Araki M,Hirakawa M. From genomics to chemical genomics: new developments in KEGGNucleic Acids Res 2006;34:D354–D357. [pmid: 16381885] [doi: 10.1093/nar/gkj102]
Bendtsen JD,Nielsen H,von Heijne G,Brunak S. Improved prediction of signal peptides: SignalP 3.0J Mol Biol 2004;340:783–795. [pmid: 15223320] [doi: 10.1016/j.jmb.2004.05.028]
Krogh A,Larsson B,von Heijne G,Sonnhammer EL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomesJ Mol Biol 2001;305:567–580. [pmid: 11152613] [doi: 10.1006/jmbi.2000.4315]
Nakai K,Horton P. PSORT: a program for detecting sorting signals in proteins and predicting their subcellular localizationTrends Biochem Sci 1999;24:34–36. [pmid: 10087920] [doi: 10.1016/S0968-0004(98)01336-X]
Conesa A,Gotz S,Garcia-Gomez JM,Terol J,Talon M,Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics researchBioinformatics 2005;21:3674–3676. [pmid: 16081474] [doi: 10.1093/bioinformatics/bti610]
Knox DP. Proteinase inhibitors and helminth parasite infectionParasite Immunol 2007;29:57–71. [pmid: 17241394] [doi: 10.1111/j.1365-3024.2006.00913.x]
Harcus YM,Parkinson J,Fernandez C,Daub J,Selkirk ME,Blaxter ML,Maizels RM. Signal sequence analysis of expressed sequence tags from the nematode Nippostrongylus brasiliensis and the evolution of secreted proteins in parasitesGenome Biol 2004;5:R39. [pmid: 15186490] [doi: 10.1186/gb-2004-5-6-r39]
Matthews JB,Davidson AJ,Beynon RJ. The application of mass spectrometry to identify immunogenic components of excretory/secretory products from adult Dictyocaulus viviparusParasitology 2004;128:S43–S47. [pmid: 16454898]
Vanholme B,De Meutter J,Tytgat T,Van Montagu M,Coomans A,Gheysen G. Secretions of plant-parasitic nematodes: a molecular updateGene 2004;332:13–27. [pmid: 15145050] [doi: 10.1016/j.gene.2004.02.024]
Delaney A,Williamson A,Brand A,Ashcom J,Varghese G,Goud GN,Hawdon JM. Cloning and characterisation of an aspartyl protease inhibitor (API-1) from Ancylostoma hookwormsInt J Parasitol 2005;35:303–313. [pmid: 15722082] [doi: 10.1016/j.ijpara.2004.11.014]
Ding X,Shields J,Allen R,Hussey RS. A secretory cellulose-binding protein cDNA cloned from the root-knot nematode (Meloidogyne incognita)Mol Plant Microbe Interact 1998;11:952–959. [pmid: 9768512] [doi: 10.1094/MPMI.1998.11.10.952]
Zhan B,Liu Y,Badamchian M,Williamson A,Feng J,Loukas A,Hawdon JM,Hotez PJ. Molecular characterisation of the Ancylostoma-secreted protein family from the adult stage of Ancylostoma caninumInt J Parasitol 2003;33:897–907. [pmid: 12906874] [doi: 10.1016/S0020-7519(03)00111-5]
McKeand JB. Vaccine development and diagnostics of Dictyocaulus viviparusParasitology 2000;120:S17–S23. [pmid: 10874707] [doi: 10.1017/S0031182099005727]
Chilton NB,Huby-Chilton F,Gasser RB,Beveridge I. The evolutionary origins of nematodes within the order Strongylida are related to predilection sites within hostsMol Phylogenet Evol 2006;40:118–128. [pmid: 16584893] [doi: 10.1016/j.ympev.2006.01.003]
NHGRI sequencing for 10 parasitic nematode species, members of the order Strongylida
von Eggeling F,Davies H,Lomas L,Fiedler W,Junker K,Claussen U,Ernst G. Tissue-specific microdissection coupled with ProteinChip array technologies: applications in cancer researchBiotechniques 2000;29:1066–1070. [pmid: 11084869]
Yatsuda AP,Krijgsveld J,Cornelissen AW,Heck AJ,de Vries E. Comprehensive analysis of the secreted proteins of the parasite Haemonchus contortus reveals extensive sequence variation and differential immune recognitionJ Biol Chem 2003;278:16941–16951. [pmid: 12576473] [doi: 10.1074/jbc.M212453200]
Barrett J,Brophy PM. Ascaris haemoglobin: new tricks for an old proteinParasitol Today 2000;16:90–91. [pmid: 10689322] [doi: 10.1016/S0169-4758(99)01613-0]
Pernthaner A,Cole SA,Morrison L,Green R,Shaw RJ,Hein WR. Cytokine and antibody subclass responses in the intestinal lymph of sheep during repeated experimental infections with the nematode parasite Trichostrongylus colubriformisVet Immunol Immunopathol 2006;114:135–148. [pmid: 16956667] [doi: 10.1016/j.vetimm.2006.08.004]
Wood IB,Amaral NK,Bairden K,Duncan JL,Kassai T,Malone JB Jr,Pankavich JA,Reinecke RK,Slocombe O,Taylor SM,Vercruysse J. World Association for the Advancement of Veterinary Parasitology (W.A.A.V.P.) second edition of guidelines for evaluating the efficacy of anthelmintics in ruminants (bovine, ovine, caprine)Vet Parasitol 1995;58:181–213. [pmid: 7571325] [doi: 10.1016/0304-4017(95)00806-2]


[Figure ID: F1]
Figure 1 

Bioinformatics analysis of D. viviparus ESTs. ESTExplorer analysis comprising Phases I (pre-processing), II (DNA-level annotation) and III (protein-level annotation), were augmented by homologue identification from nematodes as well as parasitic nematodes, using specialized databases.

[Figure ID: F2]
Figure 2 

Comparison of D. viviparus rESTs with C. elegans, parasitic nematode and non-nematode protein sequence databases using SimiTri. The numbers at each vertex indicate the number of rESTs matching only that specific database. The numbers on the edges indicate the number of rESTs matching the two databases linked by that edge. The number within the triangle indicates the number of D. viviparus genes with matches to all three databases.

[TableWrap ID: T1] Table 1 

Preliminary analysis of the 5271 ESTs

Numbers (percentage)
Raw sequences obtained 5271
Curated sequences 4436 (84)
Clusters of multiple sequences (contigs) 458 (8.6)
Clusters of singletons 1800 (34.2)
Total 2258 (42.8)

[TableWrap ID: T3] Table 3 

Gene Ontology mappings (using GO slim terms) for Dictyocaulus viviparus clusters. Note that individual GO categories can have multiple mappings

Categories and subcategories Representation % Representation of total
a. Biological Process 540 8.53
cellular process 277 4.38
?cell communication 56 0.88
?cell recognition 1 0.02
?cell differentiation 33 0.52
response to stimulus 146 2.31
?response to biotic stimulus 11 0.17
?response to abiotic stimulus 30 0.47
?response to stress 35 0.55
?behavior 114 1.8
reproduction 180 2.84
growth 180 2.84
physiological process 424 6.7
?response to endogenous stimulus 10 0.16
?response to external stimulus 16 0.25
?response to stress 35 0.55
?death 35 0.55
?metabolism 256 4.05
?homeostasis 5 0.08
development 284 4.49
?embryonic development 228 3.6
??anatomical structure development 103 1.63
localization 54 0.85
b. Cellular component 328 5.18
cell 311 4.91
?intracellular 264 4.18
organelle 215 3.4
?vesicle 9 0.04
?membrane-bound organelle 152 2.4
??organelle envelope 1 0.02
?non-membrane-bound organelle 86 1.36
organelle part 21 0.33
?organelle lumen 18 0.28
envelope 1 0.02
protein complex 22 0.35
?unlocalized protein complex 5 0.08
cell part 280 4.42
extracellular matrix (sensu Metazoa) 10 0.16
extracellular region 17 0.27
extracellular region part 15 0.24
?extracellular space 5 0.08
c. Molecular function 457 7.22
binding 253 4
?carbohydrate binding 4 0.06
?ion binding 13 0.21
?nucleotide binding 65 1.03
?nucleic acid binding 80 1.26
?chromatin binding 2 0.03
?lipid binding 3 0.05
?protein binding 107 1.69
antioxidant activity 2 0.03
catalytic activity 214 3.38
hydrolase activity 90 1.42
transferase activity 51 0.81
enzyme regulator activity 11 0.17
electron transporter activity 10 0.16
motor activity 1 0.02
signal transducer activity 12 0.19
?receptor activity 4 0.06
structural molecule activity 214 1.69
transcription regulator activity 68 1.07
transporter activity 28 0.44
?ion transporter activity 1 0.02
?channel or pore class transporter activity 1 0.02
translation regulator activity 8 0.13

[TableWrap ID: T2] Table 2 

The most abundant transcripts in adult Dictyocaulus viviparus

No Cluster ID Percentage No of ESTs Gene ID E-value % identity Description from NR hit C. elegans Homologue Top hit in C. elegans C. elegans RNAi C. elegans gene ontology
1 DvContig 413 100 261 - - - No significant similarity (novel) No significant similarity (novel) None None None
2 DvContig 438 94.64 247 CAF29502.1 6E-64 118/126 (93%) major sperm protein [Oesophagostomum dentatum] msp-76 ? (Major Sperm Protein) ZK1251.6 fat content increased embryonic development; structural molecule activity
3 DvContig 329 31.03 81 ABA53863.1 2E-136 220/333 (66%) cathepsin B-like cysteine protease 1 [Parelaphostrongylus tenuis] cysteine protease F57F5.1 embryonic lethal (Let) locomotion abnormal (unc) larval arrest (Lva) cathepsin B activity
4 DvContig 371 27.20 71 - - - No significant similarity (novel) No significant similarity (novel) None None None
5 DvContig 366 25.67 67 CAF29502.1 4E-65 118/126 (93%) major sperm protein [Oesophagostomum dentatum] msp-76 ? (Major Sperm Protein) ZK1251.6 fat content increased embryonic development; structural molecule activity
6 DvContig 318 18.39 48 - - - No significant similarity (novel) No significant similarity (novel) None None None
7 DvContig 238 16.48 43 AAO63577.1 3E-15 59/199 (29%) secreted protein 5 precursor [Ancylostoma caninum] vap-1 ? (Venom-Allergen-like Protein) F11C7.3b No observed phenotype is found. extracellular region
8 DvContig 419 13.03 34 - - - No significant similarity (novel) No significant similarity (novel) None None None
9 DvContig 415 12.26 32 NP_500698.1 3E-26 142/237 (59%) Sperm-Specific family, class Q family member (ssq-2) [Caenorhabditis elegans] ssq-2 ? (Sperm-Specific family, class Q) T28H11.5 locomotion abnormal (unc) locomotory behavior; GPCR protein signaling pathway; integral to membrane; leukotriene receptor activity; structural molecule activity
10 DvContig 140 11.88 31 NP_579952.2 7E-29 69/86 (80%) cytochrome oxidase subunit I [Ancylostoma duodenale] cytochrome oxidase subunit I MTCE.26 No observed phenotype is found. electron transport; membrane; cytochrome C-oxidase activity
11 DvContig 385 11.88 31 NP_500520.1 3E-62 230/300 (76%) COLlagen family member (col-3) [Caenorhabditis elegans] col-34 cuticular collagen F36A4.10 organism morphology abnormal (Bmd), dumpy (Dpy), locomotion abnormal (unc) locomotory behavior; GPCR protein signaling pathway; integral to membrane; leukotriene receptor activity; structural molecule activity
12 DvContig 357 9.96 26 NP_500697.1 2E-26 70/110 (63%) Sperm-Specific family, class P family member (ssp-19) [Caenorhabditis. elegans] ssp-11 ? (Sperm-Specific family, class P) T28H11.6 No observed phenotype is found. structural molecule activity
13 DvContig 347 9.20 24 AAO63577.1 2.00E-15 59/199 (29%) secreted protein 5 precursor [Ancylostoma caninum] vap-1 ? (Venom-Allergen-like Protein) F11C7.3b No observed phenotype is found. extracellular region
14 DvContig 253 8.81 23 - - - No significant similarity (novel) No significant similarity (novel) None None None
15 DvContig 261 7.66 20 AAO63577.1 6E-16 67/198 (33%) secreted protein 5 precursor [Ancylostoma caninum] testes-specific protein like T05A10.5 No observed phenotype is found. extracellular region
16 DvContig 352 7.66 20 NP_506519.1 1E-27 82/268 (30%) Hypothetical protein CBG09313 [Caenorhabditis briggsae] Domain DUF856 T16A9.5 gene (Unnamed protein) T16A9.5 No observed phenotype is found. dicarboxylic acid transport
17 DvContig 390 7.66 20 - - - No significant similarity (novel) No significant similarity (novel) None None None
18 DvContig 258 7.28 19 - - - No significant similarity (novel) No significant similarity (novel) None None None
19 DvContig 403 7.28 19 - - - No significant similarity (novel) No significant similarity (novel) None None None
20 DvContig 178 6.90 18 - - - No significant similarity (novel) No significant similarity (novel) None None None
21 DvContig 86 5.75 15 AAO63577.1 1E-16 64/215 (29%) secreted protein 5 precursor [Ancylostoma caninum] testes-specific protein like T05A10.5 No observed phenotype is found. extracellular region
22 DvContig 346 5.75 15 - - - No significant similarity (novel) No significant similarity (novel) None None None
23 DvContig 348 5.36 14 AAP41952.1 3E-12 55/203 (27%) secreted protein ASP-2 [Necator americanus] vap-1 ? (Venom-Allergen-like Protein) F11C7.3b No observed phenotype is found. None
24 DvContig 224 4.98 13 BAA12092.1 5E-168 310/366 (84%) aldolase Ce2 [Caenorhabditis elegans] Fructose-biphosphate aldolase F01F1.12a embryonic lethal (Let), sterile progeny (Stp), egg laying abnormal (Egl), sick (Sck), locomotion abnormal (unc), (Emb) slow growth (Gro) None
25 DvContig 354 4.98 13 NP_510410.1 1E-16 49/77 (63%) HIStone family member (his-24) [Caenorhabditis elegans] his-24 histone H1 M163.3 No observed phenotype is found. nucleosome assembly; DNA binding; nucleus
26 DvContig 216 4.60 12 CAE67138.1 7E-141 233/309 (75%) casein kinase [Caenorhabditis briggsae] casein kinase C39H7.1 No observed phenotype is found. protein amino acid phosphorylation; ATP binding; protein kinase activity; protein serine/threonine kinase activity; protein-tyrosine kinase activity
27 DvContig 161 4.21 11 AAO63577.1 3E-15 64/215 (29%) secreted protein 5 precursor [Ancylostoma caninum] vap-1 ? (Venom-Allergen-like Protein) F11C7.3b No observed phenotype is found. extracellular region
28 DvContig 202 4.21 11 - - - No significant similarity (novel) No significant similarity (novel) None None None
29 DvContig 266 4.21 11 - - - No significant similarity (novel) No significant similarity (novel) None None None
30 DvContig 356 4.21 11 NP_492448.1 3E-31 76/191 (39%) C-type Lectin [Caenorhabditis elegans] clec-87 ? (C-type LECtin) C25A1.8 No observed phenotype is found. sugar binding activity
31 DvContig 154 3.83 10 NP_492457.1 0 785/852 (92%) Elongation Factor family member (eft-2) [Caenorhabditis elegans] eft-2 ? (Elongation FacTor) F25H5.4 embryonic lethal (Let), larval arrest (Lva, maternal sterile (Ste), protruding vulva (Pvl) embryonic development; translational elongation; translational termination

[TableWrap ID: T4] Table 4 

Top 30 selected metabolic pathways in adult D. viviparus mapped by Kyoto Encyclopedia of Genes and Genomes

Number KEGG Pathway rESTs sequence count Enzymes
1 Signal transduction mechanisms 26 26
2 Glycolysis/Gluconeogenesis 24 24
3 Pyruvate metabolism 19 19
4 Regulation of actin cytoskeleton 19 1
5 Ribosome 19 0
6 Protein folding and associated processing 18 14
7 Focal adhesion 17 1
8 GTP-binding proteins 15 0
9 Insulin signaling pathway 14 11
10 Carbon fixation 12 12
11 Purine metabolism 12 12
12 Antigen processing and presentation 11 6
13 Glycerophospholipid metabolism 11 11
14 MAPK signaling pathway 11 3
15 Translation factors 10 4
16 Propanoate metabolism 9 9
17 Fructose and mannose metabolism 9 9
18 Citrate cycle (TCA cycle) 9 9
19 Butanoate metabolism 9 9
20 Tight junction 9 3
21 Wnt signaling pathway 9 3
22 Oxidative phosphorylation 9 9
23 Leukocyte transendothelial migration 8 0
24 Adherens junction 8 2
25 Ion channels 8 0
26 CD molecules 7 5
27 Adipocytokine signaling pathway 7 5
28 Proteasome 7 2
29 Valine, leucine and isoleucine degradation 7 7
30 Other enzymes 37 37

Article Categories:
  • Research Article

Previous Document:  Acute graft versus host disease.
Next Document:  A systematic review of cooling for neuroprotection in neonates with hypoxic ischemic encephalopathy ...