The rise of taxon-specific epitope predictors

Computational predictors of immunogenic peptides, or epitopes, are traditionally built based on data from a broad range of pathogens without consideration for taxonomic information. While this approach may be reasonable if one aims to develop one-size-fits-all models, it may be counterproductive if the proteins for which the model is expected to generalize are known to come from a specific subset of phylogenetically related pathogens. There is mounting evidence that, for these cases, taxon-specific models can outperform generalist ones, even when trained with substantially smaller amounts of data. In this comment, we provide some perspective on the current state of taxon-specific modelling for the prediction of linear B-cell epitopes, and the challenges faced when building and deploying these predictors.

Comprehensive insights into arsenic- and iron-redox genes, their taxonomy and associated environmental drivers deciphered by a meta-analysis

In nature, arsenic (As) and iron (Fe) biotransformation are interconnected, influencing local As mobility and toxicity. While As- or Fe-metabolizing microorganisms are widely documented, knowledge concerning their cycling genes, associated with geophysicochemical data and taxonomic distribution, remains scarce. We performed a meta-analysis to explore the distribution and environmental importance of As- and Fe-redox genes (AsRGs and FeRGs) and predict their significant correlations and hosts. The most abundant and ubiquitous AsRGs and FeRGs were arsC and ccoN, respectively. The ccoN gene had the highest frequency at pH ≥ 9.1, in which dissolved Fe(II) is scarce, possibly contributing to enhanced host survival. Fe(III) oxidation genes iro and ccoN appear to be associated with As(V) detoxification in mesophilic environments. No correlation was observed between Fe(III) reduction gene omcB and arsenate reductase genes. Cytochromes with putative roles in Fe-redox reactions were identified (including yceJ and fbcH) and were significantly correlated with As(V) reduction genes under diverse geophysicochemical conditions. The taxonomies of AsRGs and FeRGs-carrying contigs revealed great diversity, among which various, such as Chlamydea (arsC) and Firmicutes (omcB), were previously undescribed. Nearly all (98.9%) of the AsRGs and FeRGs were not carried by any plasmid sequences. This meta-analysis expands our understanding of the global environmental, taxonomic and functional microbiome involved in As- and Fe-redox transformations. Moreover, these findings should help guide studies on putative in vivo functional roles of cytochromes in Fe-redox pathways.

Genome-wide analyses reveal a highly conserved Dengue virus envelope peptide which is critical for virus viability and antigenic in humans

Targeting regions of proteins that show a high degree of structural conservation has been proposed as a method of developing immunotherapies and vaccines that may bypass the wide genetic variability of RNA viruses. Despite several attempts, a vaccine that protects evenly against the four circulating Dengue virus (DV) serotypes remains elusive. To find critical conserved amino acids in dengue viruses, 120 complete genomes of each serotype were selected at random and used to calculate conservation scores for nucleotide and amino acid sequences. The identified peptide sequences were analysed for their structural conservation and localisation using crystallographic data. The longest, surface exposed, highly conserved peptide of Envelope protein was found to correspond to amino acid residues 250 to 270. Mutation of this peptide in DV1 was lethal, since no replication of the mutant virus was detected in human cells. Antibodies against this peptide were detected in DV naturally infected patients indicating its potential antigenicity. Hence, this study has identified a highly conserved, critical peptide in DV that is a target of antibodies in infected humans.

POTION: an end-to-end pipeline for positive Darwinian selection detection in genome-scale data through phylogenetic comparison of protein-coding genes

Background

Detection of genes evolving under positive Darwinian evolution in genome-scale data is nowadays a prevailing strategy in comparative genomics studies to identify genes potentially involved in adaptation processes. Despite the large number of studies aiming to detect and contextualize such gene sets, there is virtually no software available to perform this task in a general, automatic, large-scale and reliable manner. This certainly occurs due to the computational challenges involved in this task, such as the appropriate modeling of data under analysis, the computation time to perform several of the required steps when dealing with genome-scale data and the highly error-prone nature of the sequence and alignment data structures needed for genome-wide positive selection detection.

Results

We present POTION, an open source, modular and end-to-end software for genome-scale detection of positive Darwinian selection in groups of homologous coding sequences. Our software represents a key step towards genome-scale, automated detection of positive selection, from predicted coding sequences and their homology relationships to high-quality groups of positively selected genes. POTION reduces false positives through several sophisticated sequence and group filters based on numeric, phylogenetic, quality and conservation criteria to remove spurious data and through multiple hypothesis corrections, and considerably reduces computation time thanks to a parallelized design. Our software achieved a high classification performance when used to evaluate a curated dataset of Trypanosoma brucei paralogs previously surveyed for positive selection. When used to analyze predicted groups of homologous genes of 19 strains of Mycobacterium tuberculosis as a case study we demonstrated the filters implemented in POTION to remove sources of errors that commonly inflate errors in positive selection detection. A thorough literature review found no other software similar to POTION in terms of customization, scale and automation.

Conclusion

To the best of our knowledge, POTION is the first tool to allow users to construct and check hypotheses regarding the occurrence of site-based evidence of positive selection in non-curated, genome-scale data within a feasible time frame and with no human intervention after initial configuration. POTION is available at http://www.lmb.cnptia.embrapa.br/share/POTION/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1765-0) contains supplementary material, which is available to authorized users.

KOMODO: a web tool for detecting and visualizing biased distribution of groups of homologous genes in monophyletic taxa

The enrichment analysis is a standard procedure to interpret ‘omics’ experiments that generate large gene lists as outputs, such as transcriptomics and protemics. However, despite the huge success of enrichment analysis in these classes of experiments, there is a surprising lack of application of this methodology to survey other categories of large-scale biological data available. Here, we report Kegg Orthology enrichMent-Online DetectiOn (KOMODO), a web tool to systematically investigate groups of monophyletic genomes in order to detect significantly enriched groups of homologous genes in one taxon when compared with another. The results are displayed in their proper biochemical roles in a visual, explorative way, allowing users to easily formulate and investigate biological hypotheses regarding the taxonomical distribution of genomic elements. We validated KOMODO by analyzing portions of central carbon metabolism in two taxa extensively studied regarding their carbon metabolism profile (Enterobacteriaceae family and Lactobacillales order). Most enzymatic activities significantly biased were related to known key metabolic traits in these taxa, such as the distinct fates of pyruvate (the known tendency of lactate production in Lactobacillales and its complete oxidation in Enterobacteriaceae), demonstrating that KOMODO could detect biologically meaningful differences in the frequencies of shared genomic elements among taxa. KOMODO is freely available at http://komodotool.org.

Evidence for reductive genome evolution and lateral acquisition of virulence functions in two Corynebacterium pseudotuberculosis strains

Background

Corynebacterium pseudotuberculosis, a Gram-positive, facultative intracellular pathogen, is the etiologic agent of the disease known as caseous lymphadenitis (CL). CL mainly affects small ruminants, such as goats and sheep; it also causes infections in humans, though rarely. This species is distributed worldwide, but it has the most serious economic impact in Oceania, Africa and South America. Although C. pseudotuberculosis causes major health and productivity problems for livestock, little is known about the molecular basis of its pathogenicity.

Methodology and Findings

We characterized two C. pseudotuberculosis genomes (Cp1002, isolated from goats; and CpC231, isolated from sheep). Analysis of the predicted genomes showed high similarity in genomic architecture, gene content and genetic order. When C. pseudotuberculosis was compared with other Corynebacterium species, it became evident that this pathogenic species has lost numerous genes, resulting in one of the smallest genomes in the genus. Other differences that could be part of the adaptation to pathogenicity include a lower GC content, of about 52%, and a reduced gene repertoire. The C. pseudotuberculosis genome also includes seven putative pathogenicity islands, which contain several classical virulence factors, including genes for fimbrial subunits, adhesion factors, iron uptake and secreted toxins. Additionally, all of the virulence factors in the islands have characteristics that indicate horizontal transfer.

Conclusions

These particular genome characteristics of C. pseudotuberculosis, as well as its acquired virulence factors in pathogenicity islands, provide evidence of its lifestyle and of the pathogenicity pathways used by this pathogen in the infection process. All genomes cited in this study are available in the NCBI Genbank database (http://www.ncbi.nlm.nih.gov/genbank/) under accession numbers CP001809 and CP001829.

Virus-host coevolution: common patterns of nucleotide motif usage in Flaviviridae and their hosts

Virus-host biological interaction is a continuous coevolutionary process involving both host immune system and viral escape mechanisms. Flaviviridae family is composed of fast evolving RNA viruses that infects vertebrate (mammals and birds) and/or invertebrate (ticks and mosquitoes) organisms. These host groups are very distinct life forms separated by a long evolutionary time, so lineage-specific anti-viral mechanisms are likely to have evolved. Flaviviridae viruses which infect a single host lineage would be subjected to specific host-induced pressures and, therefore, selected by them. In this work we compare the genomic evolutionary patterns of Flaviviridae viruses and their hosts in an attempt to uncover coevolutionary processes inducing common features in such disparate groups. Especially, we have analyzed dinucleotide and codon usage patterns in the coding regions of vertebrate and invertebrate organisms as well as in Flaviviridae viruses which specifically infect one or both host types. The two host groups possess very distinctive dinucleotide and codon usage patterns. A pronounced CpG under-representation was found in the vertebrate group, possibly induced by the methylation-deamination process, as well as a prominent TpA decrease. The invertebrate group displayed only a TpA frequency reduction bias. Flaviviridae viruses mimicked host nucleotide motif usage in a host-specific manner. Vertebrate-infecting viruses possessed under-representation of CpG and TpA, and insect-only viruses displayed only a TpA under-representation bias. Single-host Flaviviridae members which persistently infect mammals or insect hosts (Hepacivirus and insect-only Flavivirus, respectively) were found to posses a codon usage profile more similar to that of their hosts than to related Flaviviridae. We demonstrated that vertebrates and mosquitoes genomes are under very distinct lineage-specific constraints, and Flaviviridae viruses which specifically infect these lineages appear to be subject to the same evolutionary pressures that shaped their host coding regions, evidencing the lineage-specific coevolutionary processes between the viral and host groups.

Schistosoma mansoni: Microarray analysis of gene expression induced by host sex

Schistosoma mansoni is a digenetic trematode and a human parasite responsible for high social and economic impact. Although some authors have studied the effect of host hormones on parasites, not much is known about the effects of host sex on gene expression in Schistosomes. In order to study gene transcripts associated with the host sex, we compared the gene expression profiles of both male and female unisexual adult S. mansoni parasites raised on either male or female hosts, using DNA microarrays. Our results show that host sex caused differential expression of at least 11 genes in female parasites and of 134 in male parasites. Of the differentially expressed genes in female worms, 10 were preferentially expressed in female worms from male mice, while of the 134 differentially expressed genes in male parasites, 79 (59%) were preferentially expressed in worms from female mice. Further investigation of the role of each of those genes will help understand better their importance in the pathogenesis of Schistosomiasis.

Schistosoma mansoni: Heterologous complementation of a yeast null mutant by SmRbx, a protein similar to a RING box protein involved in ubiquitination

The SCF (Skp1-Cul1-F-box) complex is one of the several E3 ligase enzymes and it catalyzes protein ubiquitination and degradation by the 26S proteasome. Rbx1 is a member of the SCF complex in humans and HRT1 is its yeast orthologue. A cDNA encoding a Schistosoma mansoni Rbx1 homolog was cloned and functionally characterized. Heterologous functional complementation in yeast showed that the worm SmRbx gene was able to complement the HRT1yeast null mutation. Gene deletion constructs for N- and C-termini truncated proteins were used to transform hrt1(-) yeast mutant strains, allowing us to observe that regions reported to be involved in the interaction with cullin1 (Cul1) were essential for SmRbx function. Yeast two-hybrid assays using SmRbx and yeast Cul1 confirmed that SmRbx, but not the mutant SmRbxDelta24N, lacking the N-terminus of the protein, was capable of interacting with Cul1. These results suggest that SmRbx protein is involved in the SCF complex formation.

Microarray analysis of gene expression induced by sexual contact in Schistosoma mansoni

Background

The parasitic trematode Schistosoma mansoni is one of the major causative agents of Schistosomiasis, a disease that affects approximately 200 million people, mostly in developing countries. Since much of the pathology is associated with eggs laid by the female worm, understanding the mechanisms involved in oogenesis and sexual maturation is an important step towards the discovery of new targets for effective drug therapy. It is known that the adult female worm only develops fully in the presence of a male worm and that the rates of oviposition and maturation of eggs are significantly increased by mating. In order to study gene transcripts associated with sexual maturation and oviposition, we compared the gene expression profiles of sexually mature and immature parasites using DNA microarrays.

Results

For each experiment, three amplified RNA microarray hybridizations and their dye swaps were analyzed. Our results show that 265 transcripts are differentially expressed in adult females and 53 in adult males when mature and immature worms are compared. Of the genes differentially expressed, 55% are expressed at higher levels in paired females while the remaining 45% are more expressed in unpaired ones and 56.6% are expressed at higher levels in paired male worms while the remaining 43.4% are more expressed in immature parasites. Real-time RT-PCR analysis validated the microarray results. Several new maturation associated transcripts were identified. Genes that were up-regulated in single-sex females were mostly related to energy generation (i.e. carbohydrate and protein metabolism, generation of precursor metabolites and energy, cellular catabolism, and organelle organization and biogenesis) while genes that were down-regulated related to RNA metabolism, reactive oxygen species metabolism, electron transport, organelle organization and biogenesis and protein biosynthesis.

Conclusion

Our results confirm previous observations related to gene expression induced by sexual maturation in female schistosome worms. They also increase the list of S. mansoni maturation associated transcripts considerably, therefore opening new and exciting avenues for the study of the conjugal biology and development of new drugs against schistosomes.

Functional complementation of a yeast knockout strain by Schistosoma mansoni Rho1 GTPase in the presence of caffeine, an agent that affects mutants defective in the protein kinase C signal transduction pathway

In a previous study, the Schistosoma mansoni Rho1 protein was able to complement Rho1 null mutant Saccharomyces cerevisiae cells at restrictive temperatures and under osmotic stress (low calcium concentration) better than the human homologue (RhoA). It is known that under osmotic stress, the S. cerevisiae Rho1 triggers two distinct pathways: activation of the membrane 1,3-beta-glucan synthase enzymatic complex and activation of the protein kinase C1 signal transduction pathway, promoting the transcription of response genes. In the present work the SmRho1 protein and its mutants smrho1E97P, smrho1L101T, and smrho1E97P, L101T were used to try to clarify the basis for the differential complementation of Rho1 knockout yeast strain by the human and S. mansoni genes. Experiments of functional complementation in the presence of caffeine and in the presence of the osmotic regulator sorbitol were conducted. SmRho1 and its mutants showed a differential complementation of the yeast cells in the presence of caffeine, since smrho1E97P and smrho1E97P, L101T mutants showed a delay in the growth when compared to the yeast complemented with the wild type SmRho1. However, in the presence of sorbitol and caffeine the wild type SmRho1 and mutants showed a similar complementation phenotype, as they allowed yeast growth in all caffeine concentrations tested.