Infection by the castrating parasitic nematode Sphaerularia bombi changes gene expression in Bombus terrestris bumblebee queens.

Parasitism can result in dramatic changes in host phenotype, which are themselves underpinned by genes and their expression. Understanding how hosts respond at the molecular level to parasites can therefore reveal the molecular architecture of an altered host phenotype. The entomoparasitic nematode Sphaerularia bombi is a parasite of bumblebee (Bombus) hosts where it induces complex behavioural changes and host castration. To examine this interaction at the molecular level, we performed genome-wide transcriptional profiling using RNA-Sequencing (RNA-Seq) of S. bombi-infected Bombus terrestris queens at two critical time-points: during and just after overwintering diapause. We found that infection by S. bombi affects the transcription of genes underlying host biological processes associated with energy usage, translation, and circadian rhythm. We also found that the parasite affects the expression of immune genes, including members of the Toll signalling pathway providing evidence for a novel interaction between the parasite and the host immune response. Taken together, our results identify host biological processes and genes affected by an entomoparasitic nematode providing the first steps towards a molecular understanding of this ecologically important host-parasite interaction.

Functional insights into the infective larval stage of Anisakis simplex s.s., Anisakis pegreffii and their hybrids based on gene expression patterns.

BACKGROUND: Anisakis simplex sensu stricto and Anisakis pegreffii are sibling species of nematodes parasitic on marine mammals. Zoonotic human infection with third stage infective larvae causes anisakiasis, a debilitating and potentially fatal disease. These 2 species show evidence of hybridisation in geographical areas where they are sympatric. How the species and their hybrids differ is still poorly understood. RESULTS: Third stage larvae of Anisakis simplex s.s., Anisakis pegreffii and hybrids were sampled from Merluccius merluccius (Teleosti) hosts captured in waters of the FAO 27 geographical area. Specimens of each species and hybrids were distinguished with a diagnostic genetic marker (ITS). RNA was extracted from pools of 10 individuals of each taxon. Transcriptomes were generated using Illumina RNA-Seq, and assembled de novo. A joint assembly (here called merged transcriptome) of all 3 samples was also generated. The inferred transcript sets were functionally annotated and compared globally and also on subsets of secreted proteins and putative allergen families. While intermediary metabolism appeared to be typical for nematodes in the 3 evaluated taxa, their transcriptomes present strong levels of differential expression and enrichment, mainly of transcripts related to metabolic pathways and gene ontologies associated to energy metabolism and other pathways, with significant presence of excreted/secreted proteins, most of them allergens. The allergome of the 2 species and their hybrids has also been thoroughly studied; at least 74 different allergen families were identified in the transcriptomes. CONCLUSIONS: A. simplex s.s., A. pegreffi and their hybrids differ in gene expression patterns in the L3 stage. Strong parent-of-origin effects were observed: A. pegreffi alleles dominate in the expression patterns of hybrids albeit the latter, and A. pegreffii also display significant differences indicating that hybrids are intermediate biological entities among their parental species, and thus of outstanding interest in the study of speciation in nematodes. Analyses of differential expression based on genes coding for secreted proteins suggests that co-infections presents different repertoires of released protein to the host environment. Both species and their hybrids, share more allergen genes than previously thought and are likely to induce overlapping disease responses.

Genomics and transcriptomics across the diversity of the Nematoda.

The diversity of biology in nematodes is reflected in the diversity of their genomes. Parasitic species in particular have evolved mechanisms to invade and outwit their hosts, and these offer opportunities for the development of control measures. Genomic analyses can reveal the molecular underpinnings of phenotypes such as parasitism and thus, initiate and support research programmes that explore the manipulation of host and parasite physiologies to achieve favourable outcomes. Wide sampling across nematode diversity allows phylogenetically informed formulation of research hypotheses, identification of core features shared by all species or important evolutionary novelties present in isolated clades. Many nematode species have been investigated through the use of the expressed sequence tag approach, which samples from the transcribed genome. Gene catalogues generated in this way can be explored to reveal the patterns of expression associated with parasitism and candidates for testing as drug targets or vaccine components. Analysis environments, such as NEMBASE facilitate exploitation of these data. The development of new high-throughput DNA-sequencing technologies has facilitated transcriptomic and genomic approaches to parasite biology. Whole genome sequencing offers more complete catalogues of genes and assists a systems approach to phenotype dissection. These efforts are being coordinated through the 959 Nematode Genomes initiative.

Ultrasequencing of the meiofaunal biosphere: practice, pitfalls and promises.

Biodiversity assessment is the key to understanding the relationship between biodiversity and ecosystem functioning, but there is a well-acknowledged biodiversity identification gap related to eukaryotic meiofaunal organisms. Meiofaunal identification is confounded by the small size of taxa, morphological convergence and intraspecific variation. However, the most important restricting factor in meiofaunal ecological research is the mismatch between diversity and the number of taxonomists that are able to simultaneously identify and catalogue meiofaunal diversity. Accordingly, a molecular operational taxonomic unit (MOTU)-based approach has been advocated for en mass meiofaunal biodiversity assessment, but it has been restricted by the lack of throughput afforded by chain termination sequencing. Contemporary pyrosequencing offers a solution to this problem in the form of environmental metagenetic analyses, but this represents a novel field of biodiversity assessment. Here, we provide an overview of meiofaunal metagenetic analyses, ranging from sample preservation and DNA extraction to PCR, sequencing and the bioinformatic interrogation of multiple, independent samples using 454 Roche sequencing platforms. We report two examples of environmental metagenetic nuclear small subunit 18S (nSSU) analyses of marine and tropical rainforest habitats and provide critical appraisals of the level of putative recombinant DNA molecules (chimeras) in metagenetic data sets. Following stringent quality control measures, environmental metagenetic analyses achieve MOTU formation across the eukaryote domain of life at a fraction of the time and cost of traditional approaches. The effectiveness of Roche 454 sequencing brings substantial advantages to studies aiming to elucidate the molecular genetic richness of not only meiofaunal, but also all complex eukaryotic communities.

Duplication and divergence: the evolution of nematode globins.

In common with many other groups, nematodes express globins with unknown functions. Nematode globin-like genes can be divided into class 1 globins, similar to vertebrate myoglobins, and a wide range of additional classes. Here we show that class 1 nematode globins possess a huge amount of diversity in gene sequence and structure. There is evidence for multiple events of gene duplication, intron insertion and loss between species, and for allelic variation effecting both synonymous and non-synonymous sites within species. We have also examined gene expression patterns in class I globins from a variety of species. The results show variation in the degree of gene expression, but the tissue specificity and temporal specificity of expression may be more conserved in the phylum. Because the structure-function relationships for the binding and transport of oxygen by globins are well understood, the consequences of genetic variation causing amino acid changes are explored. The gene family shows great promise for discovering unique insights into both structure-function relationships of globins and their physiologial roles.

Comparative transcriptomic responses to chronic cadmium, fluoranthene, and atrazine exposure in Lumbricus rubellus.

Transcriptional responses of a soil-dwelling organism (the earthworm Lumbricus rubellus) to three chemicals, cadmium (Cd), fluoranthene (FA), and atrazine (AZ), were measured following chronic exposure, with the aim of identifying the nature of any shared transcriptional response. Principal component analysis indicated full or partial separation of control and exposed samples for each compound but not for the composite set of all control and exposed samples. Partial least-squares discriminant analysis allowed separation of the control and exposed samples for each chemical and also for the composite data set, suggesting a common transcriptional response to exposure. Genes identified as changing in expression level (by the least stringent test for significance) following exposure to two chemicals indicated a substantial number of common genes (> 127). The three compound overlapping gene set, however, comprised only 25 genes. We suggest that the low commonality in transcriptional response may be linked to the chronic concentrations (approximately 10% EC50) and chronic duration (28 days) used. Annotations of the three compound overlapping gene set indicated that genes from pathways most often associated with responses to environmental stress, such as heat shock, phase I and II metabolism, antioxidant defense, and cation balance, were not represented. The strongest annotation signature was for genes important in mitochondrial function and energy metabolism.

An expressed sequence tag survey of gene expression in the pond snail Lymnaea stagnalis, an intermediate vector of trematodes [corrected].

The pond snail Lymnaea stagnalis is an intermediate vector for the liver fluke Fasciola hepatica, a common parasite of ruminants and humans. Yet, despite being a disease of medical and economic importance, as well as a potentially useful comparative tool, the genetics of the relationship between Lymnaea and Fasciola has barely been investigated. As a complement to forthcoming F. hepatica expressed sequence tags (ESTs), we generated 1320 ESTs from L. stagnalis central nervous system (CNS) libraries. We estimate that these sequences derive from 771 different genes, of which 374 showed significant similarity to proteins in public databases, and 169 were similar to ESTs from the snail vector Biomphalaria glabrata. These L. stagnalis ESTs will provide insight into the function of the snail CNS, as well as the molecular components of behaviour and response to parasitism. In the future, the comparative analysis of Lymnaea/Fasciola with Biomphalaria/Schistosoma will help to understand both conserved and divergent aspects of the host-parasite relationship. The L. stagnalis ESTs will also assist gene prediction in the forthcoming B. glabrata genome sequence. The dataset is available for searching on the world-wide web at http://zeldia.cap.ed.ac.uk/mollusca.html.

Characterisation of the two most abundant genes in the Haemonchus contortus expressed sequence tag dataset.

Analysis of the Haemonchus contortus Expressed sequence tag (EST) dataset revealed that almost 10% of all ESTs (1719 ESTs) belong to a family of related genes. Close analysis of the ESTs suggests that these represent two genes (called here Hc-nim-1 and Hc-nim-2) with multiple alleles of each. These genes show significant similarity to two genes from Caenorhabditis elegans, F54D5.3 (Wormbase accession WBGene00010049, corresponding protein WP:CE28033) and F54D5.4 (WBGene00010050, WP:CE03409) of unknown function. Reverse transcriptase coupled-PCR showed that both genes are transcribed from the L4 stage onwards and are transcribed in both male and female adult worms. A partial bacterial recombinant of the Hc-NIM-1 protein was made and used to raise antiserum in rabbits which recognised a 19 kDa antigen in the water soluble protein fraction of adult worms. By immunohistochemistry, the Hc-NIM-1 protein was localised in the hypodermis of the pharyngeal region of adult worms but not posterior in the hypodermis surrounding the reproductive tract. To investigate the function of this novel protein family we conducted a RNA interference experiment for the homologuous proteins in C. elegans. No visible phenotype was detected after simultaneous RNAi treatment for both Ce-F54D5.3 and Ce-F54D5.4.

Expressed sequence tag survey of gene expression in the scab mite Psoroptes ovis--allergens, proteases and free-radical scavengers.

Psoroptes ovis, the causative agent of sheep scab, is an important ectoparasitic mite infecting sheep, goats and cattle. Infection is characterized by an extensive dermatitis, scab formation and intense itching. Initial focal lesions spread outwards, coalesce and may extend over the whole body. The host response to infestation has all the characteristics of an immediate-type hypersensitivity reaction but the mite antigens and allergens which initiate this response are almost completely undefined. Here, 507 randomly selected cDNAs derived from a mixed population of P. ovis were sequenced and the resultant nucleotide sequences subjected to Cluster analysis and Blast searches. This analysis yielded 280 clusters of which 49 had > 1 sequence with 24 showing significant Blast X homology to another protein in the databases. There were 231 sequences which appeared on one occasion and 109 of these showed significant Blast X homology to other sequences in the databases. This analysis identified homologues of 9 different types of allergens which have been characterized in other allergic conditions such as responses to house dust mites. It also identified a number of cysteine proteases which may contribute to lesion development as well as several free-radical scavenging enzymes which may protect the mite from host immune effector responses.

Comparison of biological, molecular, and morphological methods of species identification in a set of cultured panagrolaimus isolates.

We have developed a molecular barcode system that uses the small subunit ribosomal RNA (SSU) sequence to define molecular operational taxonomic units (MOTU) of soil nematodes. Here we attempt to differentiate five cultured isolates of a taxonomically difficult genus, Panagrolaimus, using morphological, molecular, and biological (breeding) criteria. The results indicated that the five culture populations belonged to two reproductively isolated species. The available morphological criteria, including scanning electron microscopy (SEM), were insufficient to differentiate among them, and all five could be classified as one morphospecies. Within-culture variation of the morphometrical data did not discern between the two biological species. Sequence data clearly separated the populations into two groups that supported the breeding results. Given this study represented only five populations of one genus, we suggest a congruence of MOTU analysis with the biological species concept. This multifaceted approach is promising for future identification of nematodes as it is simple, comparable, and transferable.

Conservation of long-range synteny and microsynteny between the genomes of two distantly related nematodes.

BACKGROUND: Comparisons between the genomes of the closely related nematodes Caenorhabditis elegans and Caenorhabditis briggsae reveal high rates of rearrangement, with a bias towards within-chromosome events. To assess whether this pattern is true of nematodes in general, we have used genome sequence to compare two nematode species that last shared a common ancestor approximately 300 million years ago: the model C. elegans and the filarial parasite Brugia malayi. RESULTS: An 83 kb region flanking the gene for Bm-mif-1 (macrophage migration inhibitory factor, a B. malayi homolog of a human cytokine) was sequenced. When compared to the complete genome of C. elegans, evidence for conservation of long-range synteny and microsynteny was found. Potential C. elegans orthologs for II of the 12 protein-coding genes predicted in the B. malayi sequence were identified. Ten of these orthologs were located on chromosome I, with eight clustered in a 2.3 Mb region. While several, relatively local, intrachromosomal rearrangements have occurred, the order, composition, and configuration of two gene clusters, each containing three genes, was conserved. Comparison of B. malayi BAC-end genome survey sequence to C. elegans also revealed a bias towards intrachromosome rearrangements. CONCLUSIONS: We suggest that intrachromosomal rearrangement is a major force driving chromosomal organization in nematodes, but is constrained by the interdigitation of functional elements of neighboring genes.

Immunological genomics of Brugia malayi: filarial genes implicated in immune evasion and protective immunity.

Filarial nematodes are metazoan parasites with genome sizes of> 100 million base pairs, probably encoding 15 000-20 000 genes. Within this considerable gene complement, it seems likely that filariae have evolved a spectrum of immune evasion products which underpin their ability to live for many years within the human host. Moreover, no suitable vaccine currently exists for human filarial diseases, and few markers have yet been established for diagnostic use. In this review, we bring together biochemical and immunological data on prominent filarial proteins with the exciting new information provided by the Filarial Genome Project's expressed sequence tag (EST) database. In this discussion, we focus on those genes with the highest immunological profile, such as inhibitors of host enzymes, cytokine homologues and stage-specific surface proteins, as well as products associated with the mosquito-borne infective larva which offer the best opportunity for an anti-filarial vaccine. These gene products provide a fascinating glimpse of the molecular repertoire which helminth parasites have evolved to manipulate and evade the mammalian immune response.

Analysis of genes expressed at the infective larval stage validates utility of Litomosoides sigmodontis as a murine model for filarial vaccine development.

We used an expressed sequence tag approach to analyze genes expressed by the infective larvae of the rodent filarial parasite Litomosoides sigmodontis. One hundred fifty two new genes were identified, including several proposed as vaccine candidates in studies with human filarial parasites. Our findings have important implications for the use of L. sigmodontis as a model for filarial infection.

The small subunit ribosomal RNA sequence of Strongyloides stercoralis.

The published small subunit rRNA (ssrRNA) gene sequences for Strongyloides ratti and Strongyloides stercoralis are remarkably divergent, particularly in the 5' 400 bases of the approximately 1700 base pair (bp) sequences. This level of divergence between species nominally in the same genus was unprecedented. We have redetermined the ssrRNA sequence of S. stercoralis and find that the published sequence is a chimaera of parasite and fungal segments. The true sequence for S. stercoralis ssrRNA is very similar to that of S. ratti.