An Improved Statistic for Detecting Over-Represented Gene Ontology Annotations in Gene Sets

Identification of the Schistosoma mansoni TNF-alpha receptor gene and the effect of human TNF-alpha on the parasite gene expression profile

Background

Schistosoma mansoni is the major causative agent of schistosomiasis. The parasite takes advantage of host signals to complete its development in the human body. Tumor necrosis factor-alpha (TNF-α) is a human cytokine involved in skin inflammatory responses, and although its effect on the adult parasite's metabolism and egg-laying process has been previously described, a comprehensive assessment of the TNF-α pathway and its downstream molecular effects is lacking.

Methodology/Principal Findings

In the present work we describe a possible TNF-α receptor (TNFR) homolog gene in S. mansoni (SmTNFR). SmTNFR encodes a complete receptor sequence composed of 599 amino acids, and contains four cysteine-rich domains as described for TNFR members. Real-time RT-PCR experiments revealed that SmTNFR highest expression level is in cercariae, 3.5 (±0.7) times higher than in adult worms. Downstream members of the known human TNF-α pathway were identified by an in silico analysis, revealing a possible TNF-α signaling pathway in the parasite. In order to simulate parasite's exposure to human cytokine during penetration of the skin, schistosomula were exposed to human TNF-α just 3 h after cercariae-to-schistosomula in vitro transformation, and large-scale gene expression measurements were performed with microarrays. A total of 548 genes with significantly altered expression were detected, when compared to control parasites. In addition, treatment of adult worms with TNF-α caused a significantly altered expression of 1857 genes. Interestingly, the set of genes altered in adults is different from that of schistosomula, with 58 genes in common, representing 3% of altered genes in adults and 11% in 3 h-old early schistosomula.

Conclusions/Significance

We describe the possible molecular elements and targets involved in human TNF-α effect on S. mansoni, highlighting the mechanism by which recently transformed schistosomula may sense and respond to this host mediator at the site of cercarial penetration into the skin.

Schistosoma mansoni is the major causative agent of schistosomiasis in the Americas. This parasite takes advantage of host signaling molecules such as cytokines and hormones to complete its development inside the host. Tumor necrosis factor-alpha (TNF-α) is one of the most important host cytokines involved in the inflammatory response. When cercariae, the infective stage, penetrates the human skin the release of TNF-α is started. In this work the authors describe the complete sequence of a possible TNF-α receptor in S. mansoni and detect that the receptor is most highly expressed in cercariae among all life cycle stages. Aiming to mimic the situation at the site of skin penetration, cercariae were mechanically transformed in vitro into schistosomula and exposed to human TNF-α. Exposure of early-developing schistosomula to the human hormone caused a large-scale change in the expression of parasite genes. Exposure of adult worms to human TNF-α caused gene expression changes as well, and the set of parasite altered genes in the adult parasite was different from that of schistosomula. This work increases the number of known signaling pathways of the parasite, and opens new perspectives into understanding the molecular components of TNF-α response as well as into possibly interfering with parasite–host interaction.

The Aedes aegypti larval transcriptome: a comparative perspective with emphasis on trypsins and the domain structure of peritrophins

The genome sequence of Aedes aegypti was recently reported. A significant amount of Expressed Sequence Tags (ESTs) were sequenced to aid in the gene prediction process. In the present work we describe an integrated analysis of the genomic and EST data, focusing on genes with preferential expression in larvae (LG), adults (AG) and in both stages (SG). A total of 913 genes (5.4% of the transcript complement) are LG, including ion transporters and cuticle proteins that are important for ion homeostasis and defense. From a starting set of 245 genes encoding the trypsin domain, we identified 66 putative LG, AG, and SG trypsins by manual curation. Phylogenetic analyses showed that AG trypsins are divergent from their larval counterparts (LG), grouping with blood-induced trypsins from Anopheles gambiae and Simulium vittatum. These results support the hypothesis that blood-feeding arose only once, in the ancestral Culicomorpha. Peritrophins are proteins that interlock chitin fibrils to form the peritrophic membrane (PM) that compartmentalizes the food in the midgut. These proteins are recognized by having chitin-binding domains with 6 conserved Cys and may also present mucin-like domains (regions expected to be highly O-glycosylated). PM may be formed by a ring of cells (type 2, seen in Ae. aegypti larvae and Drosophila melanogaster) or by most midgut cells (type 1, found in Ae. aegypti adult and Tribolium castaneum). LG and D. melanogaster peritrophins have more complex domain structures than AG and T. castaneum peritrophins. Furthermore, mucin-like domains of peritrophins from T. castaneum (feeding on rough food) are lengthier than those of adult Ae. aegypti (blood-feeding). This suggests, for the first time, that type 1 and type 2 PM may have variable molecular architectures determined by different peritrophins and/or ancillary proteins, which may be partly modulated by diet.

The landscape of human proteins interacting with viruses and other pathogens

Infectious diseases result in millions of deaths each year. Mechanisms of infection have been studied in detail for many pathogens. However, many questions are relatively unexplored. What are the properties of human proteins that interact with pathogens? Do pathogens interact with certain functional classes of human proteins? Which infection mechanisms and pathways are commonly triggered by multiple pathogens? In this paper, to our knowledge, we provide the first study of the landscape of human proteins interacting with pathogens. We integrate human-pathogen protein-protein interactions (PPIs) for 190 pathogen strains from seven public databases. Nearly all of the 10,477 human-pathogen PPIs are for viral systems (98.3%), with the majority belonging to the human-HIV system (77.9%). We find that both viral and bacterial pathogens tend to interact with hubs (proteins with many interacting partners) and bottlenecks (proteins that are central to many paths in the network) in the human PPI network. We construct separate sets of human proteins interacting with bacterial pathogens, viral pathogens, and those interacting with multiple bacteria and with multiple viruses. Gene Ontology functions enriched in these sets reveal a number of processes, such as cell cycle regulation, nuclear transport, and immune response that participate in interactions with different pathogens. Our results provide the first global view of strategies used by pathogens to subvert human cellular processes and infect human cells. Supplementary data accompanying this paper is available at http://staff.vbi.vt.edu/dyermd/publications/dyer2008a.html.

Use of a 44k oligoarray to explore the transcriptome of Schistosoma mansoni adult worms

Recent advances in the study of Schistosoma mansoni genome and transcriptome have led to a better description of the S. mansoni gene complement. In this work, we report the design and use of a new S. mansoni 60-mer oligonucleotide microarray platform with approximately 44,000 probes, based on all publicly available cDNA sequence data for S. mansoni and Schistosoma japonicum. The large number of probes combined with the extensive sequence annotation available allowed a comprehensive approach, where most of the S. mansoni transcriptome is represented. Hybridization with adult worm RNA pointed to a set of genes transcriptionally active in this stage of the parasite's life cycle. Interestingly, a large proportion (43%) of genes for which transcription was detected in adults is comprised of "no match" genes, i.e. S. mansoni genes with unknown function and no identifiable orthologs in GenBank. Moreover, detection of bi-directional transcription for 7% of the active "no match" genes in adults leads us to hypothesize a widespread production of antisense RNA in S. mansoni, with possible regulatory roles.

Similarity of synonymous substitution rates across mammalian genomes

Given that a gene has a high (or low) synonymous substitution rate in one mammalian species, will it also have a high (or low) synonymous substitution rate in another mammalian species? Such similarities in the rate of synonymous substitution can reveal both selective pressures and neutral processes acting on mammalian gene sequences; however, the existence of such an effect has been a matter of disagreement. We resolve whether such synonymous substitution rate similarities exist using 7462 ortholog triplets aligned across rat, mouse, and human, a dataset two orders of magnitude larger than previous studies. We find that a gene's synonymous substitution rate in the rat-mouse branch of the phylogeny is correlated with its rate in the branch connecting human and the rat-mouse ancestor. We confirm this for several different measures of synonymous substitution rate, including corrections for base composition and CpG dinucleotides, and we verify the results in the larger mouse-human-rat-dog phylogeny. This similarity of rates is most apparent for genes in which synonymous sites are well conserved across species, suggesting that a significant component of the effect is due to purifying selection. We observe rate correlations at a resolution as fine as a few hundred kilobases, and the genes with the most similar synonymous substitution rates are enriched for regulatory functions. Genes with above-average substitution rates also exhibit significant, though somewhat weaker, rate correlations, suggesting that some neutral processes may have persisted in the phylogeny as well.