In silico gene function prediction using ontology-based pattern identification

MOTIVATION: With the emergence of genome-wide expression profiling data sets, the guilt by association (GBA) principle has been a cornerstone for deriving gene functional interpretations in silico. Given the limited success of traditional methods for producing clusters of genes with great amounts of functional similarity, new data-mining algorithms are required to fully exploit the potential of high-throughput genomic approaches. RESULTS: Ontology-based pattern identification (OPI) is a novel data-mining algorithm that systematically identifies expression patterns that best represent existing knowledge of gene function. Instead of relying on a universal threshold of expression similarity to define functionally related groups of genes, OPI finds the optimal analysis settings that yield gene expression patterns and gene lists that best predict gene function using the principle of GBA. We applied OPI to a publicly available gene expression data set on the life cycle of the malarial parasite Plasmodium falciparum and systematically annotated genes for 320 functional categories based on current Gene Ontology annotations. An ontology-based hierarchical tree of the 320 categories provided a systems-wide biological view of this important malarial parasite.     

A Bovine Lymphosarcoma Cell Line Infected with Theileria annulata Exhibits an Irreversible Reconfiguration of Host Cell Gene Expression

Theileria annulata, an intracellular parasite of bovine lymphoid cells, induces substantial phenotypic alterations to its host cell including continuous proliferation, cytoskeletal changes and resistance to apoptosis. While parasite induced modulation of host cell signal transduction pathways and NFκB activation are established, there remains considerable speculation on the complexities of the parasite directed control mechanisms that govern these radical changes to the host cell. Our objectives in this study were to provide a comprehensive analysis of the global changes to host cell gene expression with emphasis on those that result from direct intervention by the parasite. By using comparative microarray analysis of an uninfected bovine cell line and its Theileria infected counterpart, in conjunction with use of the specific parasitacidal agent, buparvaquone, we have identified a large number of host cell gene expression changes that result from parasite infection. Our results indicate that the viable parasite can irreversibly modify the transformed phenotype of a bovine cell line. Fifty percent of genes with altered expression failed to show a reversible response to parasite death, a possible contributing factor to initiation of host cell apoptosis. The genes that did show an early predicted response to loss of parasite viability highlighted a sub-group of genes that are likely to be under direct control by parasite infection. Network and pathway analysis demonstrated that this sub-group is significantly enriched for genes involved in regulation of chromatin modification and gene expression. The results provide evidence that the Theileria parasite has the regulatory capacity to generate widespread change to host cell gene expression in a complex and largely irreversible manner.     

Compounds of the upper gastrointestinal tract induce rapid and efficient excystation of Entamoeba invadens

The infective stage of Entamoeba parasites is an encysted form. This stage can be readily generated in vitro, which has allowed identification of stimuli that trigger the differentiation of the parasite trophozoite stage into the cyst stage. Studies of the second differentiation event, emergence of the parasite from the cyst upon infection of a host, have been hampered by the lack of an efficient means to excyst the parasite and complete the life cycle in vitro. We have determined that a combination of exposures to water, bicarbonate and bile induces rapid excystment of Entamoeba invadens cysts. The high efficiency of this method has allowed the visualization of the dynamics of the process by electron and confocal microscopy, and should permit the analysis of stage-specific gene expression and high-throughput screening of inhibitory compounds.     

Specific Gene Expression Responses to Parasite Genotypes Reveal Redundancy of Innate Immunity in Vertebrates

Vertebrate innate immunity is the first line of defense against an invading pathogen and has long been assumed to be largely unspecific with respect to parasite/pathogen species. However, recent phenotypic evidence suggests that immunogenetic variation, i.e. allelic variability in genes associated with the immune system, results in host-parasite genotype-by-genotype interactions and thus specific innate immune responses. Immunogenetic variation is common in all vertebrate taxa and this reflects an effective immunological function in complex environments. However, the underlying variability in host gene expression patterns as response of innate immunity to within-species genetic diversity of macroparasites in vertebrates is unknown. We hypothesized that intra-specific variation among parasite genotypes must be reflected in host gene expression patterns. Here we used high-throughput RNA-sequencing to examine the effect of parasite genotypes on gene expression patterns of a vertebrate host, the three-spined stickleback (Gasterosteus aculeatus). By infecting naïve fish with distinct trematode genotypes of the species Diplostomum pseudospathaceum we show that gene activity of innate immunity in three-spined sticklebacks depended on the identity of an infecting macroparasite genotype. In addition to a suite of genes indicative for a general response against the trematode we also find parasite-strain specific gene expression, in particular in the complement system genes, despite similar infection rates of single clone treatments. The observed discrepancy between infection rates and gene expression indicates the presence of alternative pathways which execute similar functions. This suggests that the innate immune system can induce redundant responses specific to parasite genotypes.     

Structural RNAs of known and unknown function identified in malaria parasites by comparative genomics and RNA analysis

As the genomes of more eukaryotic pathogens are sequenced, understanding how molecular differences between parasite and host might be exploited to provide new therapies has become a major focus. Central to cell function are RNA-containing complexes involved in gene expression, such as the ribosome, the spliceosome, snoRNAs, RNase P, and telomerase, among others. In this article we identify by comparative genomics and validate by RNA analysis numerous previously unknown structural RNAs encoded by the Plasmodium falciparum genome, including the telomerase RNA, U3, 31 snoRNAs, as well as previously predicted spliceosomal snRNAs, SRP RNA, MRP RNA, and RNAse P RNA. Furthermore, we identify six new RNA coding genes of unknown function. To investigate the relationships of the RNA coding genes to other genomic features in related parasites, we developed a genome browser for P. falciparum (http://areslab.ucsc.edu/cgi-bin/hgGateway). Additional experiments provide evidence supporting the prediction that snoRNAs guide methylation of a specific position on U4 snRNA, as well as predicting an snRNA promoter element particular to Plasmodium sp. These findings should allow detailed structural comparisons between the RNA components of the gene expression machinery of the parasite and its vertebrate hosts.     

Haemonchus contortus: molecular cloning, sequencing, and expression analysis of the gene coding for the small subunit of ribonucleotide reductase

Gastro-intestinal (GI) parasites are of great agricultural importance, annually costing the livestock industry vast amounts in resources to control parasitism. One such GI parasite, Haemonchus contortus, is principally pathogenic to sheep; with the parasite's blood-feeding behaviour causing effects ranging from mild anaemia to mortality in young animals. Current means of control, which are dependent on repeated treatment with anthelmintic chemicals, have led to increasing drug resistance. Together with the growing concern over residual chemicals in the environment and food chain, this has led to attempts to better understand the biology of the parasite with the aim to develop alternate or supplementary means of control, including the development of molecular vaccines. As a first step towards the understanding of the synthesis of deoxyribonucleotides in H. contortus, and its potential application to therapeutic control of this economically important parasite, we report the cloning, sequencing, and mRNA expression analysis of the ribonucleotide reductase R2 gene.     

Plasmodium yoelii: semiquantitative analyses of circumsporozoite protein gene expression during the sporogonic development of P. y. yoelii and P. y. nigeriensis in the mosquito vector Anopheles stephensi

Malaria infection in the mosquito vector can be modulated by the vertebrate host, mosquito factors, and interactions between different parasite populations. Modulation of parasite development can be assessed through the study of gene expression. The present report describes a specific, sensitive, and nonradioactive method that permits assessment of parasite load and quantification of circumsporozoite protein gene expression during the sporogonic stages of Plasmodium yoelii yoelii and P. y. nigeriensis. A decrease in parasite load was observed when comparing DNA of oocysts on day 7 postinfection with that of oocysts and sporozoites on day 19. On day 7, parasites (oocysts) showed a marked increase of circumsporozoite protein expression when compared with that (sporozoites and oocysts) on day 19. The method developed in this work can be a valuable tool to understand parasite interaction mechanisms that are involved in mosquito malaria infections.     

Schistosoma mansoni: the dicer gene and its expression

RNA interference (RNAi) is a gene silencing mechanism that plays an important role in regulating gene expression in many eukaryotes and has become a valuable molecular tool for analyzing gene function. Multi-domain nucleases called Dicer proteins play pivotal roles in RNAi. In this paper, we characterize the structure and expression of the Dicer gene from the platyhelminth parasite Schistosoma mansoni. The gene (SmDicer) is over 54kb long and comprises 30 exons that potentially encode a 2641 amino acid protein. This is the largest Dicer protein yet described. SmDicer contains all domains that are characteristic of metazoan dicers including an amino terminal helicase domain, DUF283, a PAZ domain, two RNAse III domains and an RNA binding domain. An examination of the available S. mansoni genome sequence suggests that the Dicer gene described here is the only Dicer gene in the parasite genome. SmDicer is expressed throughout schistosome development suggesting that RNAi technologies might be employed in deciphering gene function in all life stages of this parasite.     

Patterning of Caenorhabditis elegans posterior structures by the Abdominal-B homolog, egl-5

The Caenorhabditis elegans body axis, like that of other animals, is patterned by the action of Hox genes. In order to examine the function of one C. elegans Hox gene in depth, we determined the postembryonic expression pattern of egl-5, the C. elegans member of the Abdominal-B Hox gene paralog group, by means of whole-mount staining with a polyclonal antibody. A major site of egl-5 expression and function is in the epithelium joining the posterior digestive tract with the external epidermis. Patterning this region and its derived structures is a conserved function of Abd-B paralog group genes in other animals. Cells that initiate egl-5 expression during embryogenesis are clustered around the presumptive anus. Expression is initiated postembryonically in four additional mesodermal and ectodermal cell lineages or tissues. Once initiated in a lineage, egl-5 expression continues throughout development, suggesting that the action of egl-5 can be regarded as defining a positional cell identity. A variety of cross-regulatory interactions between egl-5 and the next more anterior Hox gene, mab-5, help define the expression domains of their respective gene products. In its expression in a localized body region, function as a marker of positional cell identity, and interactions with another Hox gene, egl-5 resembles Hox genes of other animals. This suggests that C. elegans, in spite of its small cell number and reproducible cell lineages, may not differ greatly from other animals in the way it employs Hox genes for regional specification during development.