Shivashankar H. Nagaraj

High-throughput functional annotation and data mining with the Blast2GO suite

Functional genomics technologies have been widely adopted in the biological research of both model and non-model species. An efficient functional annotation of DNA or protein sequences is a major requirement for the successful application of these approaches as functional information on gene products is often the key to the interpretation of experimental results. Therefore, there is an increasing need for bioinformatics resources which are able to cope with large amount of sequence data, produce valuable annotation results and are easily accessible to laboratories where functional genomics projects are being undertaken. We present the Blast2GO suite as an integrated and biologist-oriented solution for the high-throughput and automatic functional annotation of DNA or protein sequences based on the Gene Ontology vocabulary. The most outstanding Blast2GO features are: (i) the combination of various annotation strategies and tools controlling type and intensity of annotation, (ii) the numerous graphical features such as the interactive GO-graph visualization for gene-set function profiling or descriptive charts, (iii) the general sequence management features and (iv) high-throughput capabilities. We used the Blast2GO framework to carry out a detailed analysis of annotation behaviour through homology transfer and its impact in functional genomics research. Our aim is to offer biologists useful information to take into account when addressing the task of functionally characterizing their sequence data.

Whole–genome characterization of chemoresistant ovarian cancer

Patients with high-grade serous ovarian cancer (HGSC) have experienced little improvement in overall survival, and standard treatment has not advanced beyond platinum-based combination chemotherapy, during the past 30 years. To understand the drivers of clinical phenotypes better, here we use whole-genome sequencing of tumour and germline DNA samples from 92 patients with primary refractory, resistant, sensitive and matched acquired resistant disease. We show that gene breakage commonly inactivates the tumour suppressors RB1, NF1, RAD51B and PTEN in HGSC, and contributes to acquired chemotherapy resistance. CCNE1 amplification was common in primary resistant and refractory disease. We observed several molecular events associated with acquired resistance, including multiple independent reversions of germline BRCA1 or BRCA2 mutations in individual patients, loss of BRCA1 promoter methylation, an alteration in molecular subtype, and recurrent promoter fusion associated with overexpression of the drug efflux pump MDR1.

Proteomics Analysis of the Excretory/Secretory Component of the Blood-feeding Stage of the Hookworm, Ancylostoma caninum

Hookworms are blood-feeding intestinal parasites of mammalian hosts and are one of the major human ailments affecting ∼600 million people worldwide. These parasites form an intimate association with the host and are able to avoid vigorous immune responses in many ways including skewing of the response phenotype to promote parasite survival and longevity. The primary interface between the parasite and the host is the excretory/secretory component, a complex mixture of proteins, carbohydrates, and lipids secreted from the surface or oral openings of the parasite. The composition of this complex mixture is for the most part unknown but is likely to contain proteins important for the parasitic lifestyle and hence suitable as drug or vaccine targets. Using a strategy combining the traditional technology of one-dimensional SDS-PAGE and the newer fractionation technology of OFFGEL electrophoresis we identified 105 proteins from the excretory/secretory products of the blood-feeding stage of the dog hookworm, Ancylostoma caninum. Highly represented among the identified proteins were lectins, including three C-type lectins and three β-galactoside-specific S-type galectins, as well as a number of proteases belonging to the three major classes found in nematodes, aspartic, cysteine, and metalloproteases. Interestingly 28% of the identified proteins were homologous to activation-associated secreted proteins, a family of cysteine-rich secreted proteins belonging to the sterol carrier protein/Tpx-1/Ag5/PR-1/Sc-7 (TAPS) superfamily. Thirty-four of these proteins were identified suggesting an important role in host-parasite interactions. Other protein families identified included hyaluronidases, lysozyme-like proteins, and transthyretin-like proteins. This work identified a suite of proteins important for the parasitic lifestyle and provides new insight into the biology of hookworm infection.

Association weight matrix for the genetic dissection of puberty in beef cattle

We describe a systems biology approach for the genetic dissection of complex traits based on applying gene network theory to the results from genome-wide associations. The associations of single-nucleotide polymorphisms (SNP) that were individually associated with a primary phenotype of interest, age at puberty in our study, were explored across 22 related traits. Genomic regions were surveyed for genes harboring the selected SNP. As a result, an association weight matrix (AWM) was constructed with as many rows as genes and as many columns as traits. Each {i, j} cell value in the AWM corresponds to the z-score normalized additive effect of the ith gene (via its neighboring SNP) on the jth trait. Columnwise, the AWM recovered the genetic correlations estimated via pedigree-based restricted maximum-likelihood methods. Rowwise, a combination of hierarchical clustering, gene network, and pathway analyses identified genetic drivers that would have been missed by standard genome-wide association studies. Finally, the promoter regions of the AWM-predicted targets of three key transcription factors (TFs), estrogen-related receptor γ (ESRRG), Pal3 motif, bound by a PPAR-γ homodimer, IR3 sites (PPARG), and Prophet of Pit 1, PROP paired-like homeobox 1 (PROP1), were surveyed to identify binding sites corresponding to those TFs. Applied to our case, the AWM results recapitulate the known biology of puberty, captured experimentally validated binding sites, and identified candidate genes and gene–gene interactions for further investigation.

Regulatory impact factors

MOTIVATION Although transcription factors (TF) play a central regulatory role, their detection from expression data is limited due to their low, and often sparse, expression. In order to fill this gap, we propose a regulatory impact factor (RIF) metric to identify critical TF from gene expression data. RESULTS To substantiate the generality of RIF, we explore a set of experiments spanning a wide range of scenarios including breast cancer survival, fat, gonads and sex differentiation. We show that the strength of RIF lies in its ability to simultaneously integrate three sources of information into a single measure: (i) the change in correlation existing between the TF and the differentially expressed (DE) genes; (ii) the amount of differential expression of DE genes; and (iii) the abundance of DE genes. As a result, RIF analysis assigns an extreme score to those TF that are consistently most differentially co-expressed with the highly abundant and highly DE genes (RIF1), and to those TF with the most altered ability to predict the abundance of DE genes (RIF2). We show that RIF analysis alone recovers well-known experimentally validated TF for the processes studied. The TF identified confirm the importance of PPAR signaling in adipose development and the importance of transduction of estrogen signals in breast cancer survival and sexual differentiation. We argue that RIF has universal applicability, and advocate its use as a promising hypotheses generating tool for the systematic identification of novel TF not yet documented as critical.

A single nucleotide polymorphism-derived regulatory gene network underlying puberty in 2 tropical breeds of beef cattle

Harsh tropical environments impose serious challenges on poorly adapted species. In beef cattle, tropical adaptation in the form of temperature and disease resistance, coupled with acclimatization to seasonal and limited forage, comes at a cost to production efficiency. Prominent among these costs is delayed onset of puberty, a challenging phenotype to manipulate through traditional breeding mechanisms. Recently, system biology approaches, including gene networks, have been applied to the genetic dissection of complex phenotypes. We aimed at developing and studying gene networks underlying cattle puberty. Our starting material comprises the association results of ~50,000 SNP on 22 traits, including age at puberty, and 2 cattle breed populations: Brahman (n = 843) and Tropical Composite (n = 866). We defined age at puberty as the age at first corpus luteum (AGECL). By capturing the genes harboring mutations minimally associated (P < 0.05) to AGECL or to a set of traits related with AGECL, we derived a gene network for each breed separately and a third network for the combined data set. At the intersection of the 3 networks, we identified candidate genes and pathways that were common to both breeds. Resulting from these analyses, we identified an enrichment of genes involved in axon guidance, cell adhesion, ErbB signaling, and glutamate activity, pathways that are known to affect pulsatile release of GnRH, which is necessary for the onset of puberty. Furthermore, we employed network connectivity and centrality parameters along with a regulatory impact factor metric to identify the key transcription factors (TF) responsible for the molecular regulation of puberty. As a novel finding, we report 5 TF (HIVEP3, TOX, EYA1, NCOA2, and ZFHX4) located in the network intersecting both breeds and interacting with other TF, forming a regulatory network that harmonizes with the recent literature of puberty. Finally, we support our network predictions with evidence derived from gene expression in hypothalamic tissue of adult cows.

Transcriptional changes in the hookworm, Ancylostoma caninum, during the transition from a free-living to a parasitic larva.

Background Third-stage larvae (L3) of the canine hookworm, Ancylostoma caninum, undergo arrested development preceding transmission to a host. Many of the mRNAs up-regulated at this stage are likely to encode proteins that facilitate the transition from a free-living to a parasitic larva. The initial phase of mammalian host invasion by A. caninum L3 (herein termed “activation”) can be mimicked in vitro by culturing L3 in serum-containing medium. Methodology/Principal Findings The mRNAs differentially transcribed between activated and non-activated L3 were identified by suppression subtractive hybridisation (SSH). The analysis of these mRNAs on a custom oligonucleotide microarray printed with the SSH expressed sequence tags (ESTs) and publicly available A. caninum ESTs (non-subtracted) yielded 602 differentially expressed mRNAs, of which the most highly represented sequences encoded members of the pathogenesis-related protein (PRP) superfamily and proteases. Comparison of these A. caninum mRNAs with those of Caenorhabditis elegans larvae exiting from developmental (dauer) arrest demonstrated unexpectedly large differences in gene ontology profiles. C. elegans dauer exiting L3 up-regulated expression of mostly intracellular molecules involved in growth and development. Such mRNAs are virtually absent from activated hookworm larvae, and instead are over-represented by mRNAs encoding extracellular proteins with putative roles in host-parasite interactions. Conclusions/Significance Although this should not invalidate C. elegans dauer exit as a model for hookworm activation, it highlights the limitations of this free-living nematode as a model organism for the transition of nematode larvae from a free-living to a parasitic state.

ESTExplorer: an expressed sequence tag (EST) assembly and annotation platform

The analysis of expressed sequence tag (EST) datasets offers a rapid and cost-effective approach to elucidate the transcriptome of an organism, but requiring several computational methods for assembly and annotation. ESTExplorer is a comprehensive workflow system for EST data management and analysis. The pipeline uses a ‘distributed control approach’ in which the most appropriate bioinformatics tools are implemented over different dedicated processors. Species-specific repeat masking and conceptual translation are in-built. ESTExplorer accepts a set of ESTs in FASTA format which can be analysed using programs selected by the user. After pre-processing and assembly, the dataset is annotated at the nucleotide and protein levels, following conceptual translation. Users may optionally provide ESTExplorer with assembled contigs for annotation purposes. Functionally annotated contigs/ESTs can be analysed individually. The overall outputs are gene ontologies, protein functional identifications in terms of mapping to protein domains and metabolic pathways. ESTExplorer has been applied successfully to annotate large EST datasets from parasitic nematodes and to identify novel genes as potential targets for parasite intervention. ESTExplorer runs on a Linux cluster and is freely available for the academic community at http://estexplorer.biolinfo.org.

Gene network analyses of first service conception in Brangus heifers: Use of genome and trait associations, hypothalamic-transcriptome information, and transcription factors

Measures of heifer fertility are economically relevant traits for beef production systems and knowledge of candidate genes could be incorporated into future genomic selection strategies. Ten traits related to growth and fertility were measured in 890 Brangus heifers (3/8 Brahman × 5/8 Angus, from 67 sires). These traits were: BW and hip height adjusted to 205 and 365 d of age, postweaning ADG, yearling assessment of carcass traits (i.e., back fat thickness, intramuscular fat, and LM area), as well as heifer pregnancy and first service conception (FSC). These fertility traits were collected from controlled breeding seasons initiated with estrous synchronization and AI targeting heifers to calve by 24 mo of age. The BovineSNP50 BeadChip was used to ascertain 53,692 SNP genotypes for ∼802 heifers. Associations of genotypes and phenotypes were performed and SNP effects were estimated for each trait. Minimally associated SNP (P < 0.05) and their effects across the 10 traits formed the basis for an association weight matrix and its derived gene network related to FSC (57.3% success and heritability = 0.06 ± 0.05). These analyses yielded 1,555 important SNP, which inferred genes linked by 113,873 correlations within a network. Specifically, 1,386 SNP were nodes and the 5,132 strongest correlations (|r| ≥ 0.90) were edges. The network was filtered with genes queried from a transcriptome resource created from deep sequencing of RNA (i.e., RNA-Seq) from the hypothalamus of a prepubertal and a postpubertal Brangus heifer. The remaining hypothalamic-influenced network contained 978 genes connected by 2,560 edges or predicted gene interactions. This hypothalamic gene network was enriched with genes involved in axon guidance, which is a pathway known to influence pulsatile release of LHRH. There were 5 transcription factors with 21 or more connections: ZMAT3, STAT6, RFX4, PLAGL1, and NR6A1 for FSC. The SNP that identified these genes were intragenic and were on chromosomes 1, 5, 9, and 11. Chromosome 5 harbored both STAT6 and RFX4. The large number of interactions and genes observed with network analyses of multiple sources of genomic data (i.e., GWAS and RNA-Seq) support the concept of FSC being a polygenic trait.

Needles in the EST Haystack: Large-Scale Identification and Analysis of Excretory-Secretory (ES) Proteins in Parasitic Nematodes Using Expressed Sequence Tags (ESTs)

Background Parasitic nematodes of humans, other animals and plants continue to impose a significant public health and economic burden worldwide, due to the diseases they cause. Promising antiparasitic drug and vaccine candidates have been discovered from excreted or secreted (ES) proteins released from the parasite and exposed to the immune system of the host. Mining the entire expressed sequence tag (EST) data available from parasitic nematodes represents an approach to discover such ES targets. Methods and Findings In this study, we predicted, using EST2Secretome, a novel, high-throughput, computational workflow system, 4,710 ES proteins from 452,134 ESTs derived from 39 different species of nematodes, parasitic in animals (including humans) or plants. In total, 2,632, 786, and 1,292 ES proteins were predicted for animal-, human-, and plant-parasitic nematodes. Subsequently, we systematically analysed ES proteins using computational methods. Of these 4,710 proteins, 2,490 (52.8%) had orthologues in Caenorhabditis elegans, whereas 621 (13.8%) appeared to be novel, currently having no significant match to any molecule available in public databases. Of the C. elegans homologues, 267 had strong “loss-of-function” phenotypes by RNA interference (RNAi) in this nematode. We could functionally classify 1,948 (41.3%) sequences using the Gene Ontology (GO) terms, establish pathway associations for 573 (12.2%) sequences using Kyoto Encyclopaedia of Genes and Genomes (KEGG), and identify protein interaction partners for 1,774 (37.6%) molecules. We also mapped 758 (16.1%) proteins to protein domains including the nematode-specific protein family “transthyretin-like” and “chromadorea ALT,” considered as vaccine candidates against filariasis in humans. Conclusions We report the large-scale analysis of ES proteins inferred from EST data for a range of parasitic nematodes. This set of ES proteins provides an inventory of known and novel members of ES proteins as a foundation for studies focused on understanding the biology of parasitic nematodes and their interactions with their hosts, as well as for the development of novel drugs or vaccines for parasite intervention and control.