Iis--integrated Interactome System: A Web-based Platform For The Annotation, Analysis And Visualization Of Protein-metabolite-gene-drug Interactions By Integrating A Variety Of Data Sources And Tools.

High-throughput screening of physical, genetic and chemical-genetic interactions brings important perspectives in the Systems Biology field, as the analysis of these interactions provides new insights into protein/gene function, cellular metabolic variations and the validation of therapeutic targets and drug design. However, such analysis depends on a pipeline connecting different tools that can automatically integrate data from diverse sources and result in a more comprehensive dataset that can be properly interpreted. We describe here the Integrated Interactome System (IIS), an integrative platform with a web-based interface for the annotation, analysis and visualization of the interaction profiles of proteins/genes, metabolites and drugs of interest. IIS works in four connected modules: (i) Submission module, which receives raw data derived from Sanger sequencing (e.g. two-hybrid system); (ii) Search module, which enables the user to search for the processed reads to be assembled into contigs/singlets, or for lists of proteins/genes, metabolites and drugs of interest, and add them to the project; (iii) Annotation module, which assigns annotations from several databases for the contigs/singlets or lists of proteins/genes, generating tables with automatic annotation that can be manually curated; and (iv) Interactome module, which maps the contigs/singlets or the uploaded lists to entries in our integrated database, building networks that gather novel identified interactions, protein and metabolite expression/concentration levels, subcellular localization and computed topological metrics, GO biological processes and KEGG pathways enrichment. This module generates a XGMML file that can be imported into Cytoscape or be visualized directly on the web. We have developed IIS by the integration of diverse databases following the need of appropriate tools for a systematic analysis of physical, genetic and chemical-genetic interactions. IIS was validated with yeast two-hybrid, proteomics and metabolomics datasets, but it is also extendable to other datasets. IIS is freely available online at: http://www.lge.ibi.unicamp.br/lnbio/IIS/.

Xylella fastidiosa gene expression analysis by DNA microarrays

Xylella fastidiosa genome sequencing has generated valuable data by identifying genes acting either on metabolic pathways or in associated pathogenicity and virulence. Based on available information on these genes, new strategies for studying their expression patterns, such as microarray technology, were employed. A total of 2,600 primer pairs were synthesized and then used to generate fragments using the PCR technique. The arrays were hybridized against cDNAs labeled during reverse transcription reactions and which were obtained from bacteria grown under two different conditions (liquid XDM2 and liquid BCYE). All data were statistically analyzed to verify which genes were differentially expressed. In addition to exploring conditions for X. fastidiosa genome-wide transcriptome analysis, the present work observed the differential expression of several classes of genes (energy, protein, amino acid and nucleotide metabolism, transport, degradation of substances, toxins and hypothetical proteins, among others). The understanding of expressed genes in these two different media will be useful in comprehending the metabolic characteristics of X. fastidiosa, and in evaluating how important certain genes are for the functioning and survival of these bacteria in plants.

Validation of reference genes for gene expression studies in the honey bee, Apis mellifera, by quantitative real-time RT-PCR

For obtaining accurate and reliable gene expression results it is essential that quantitative real-time RT-PCR (qRT-PCR) data are normalized with appropriate reference genes. The current exponential increase in postgenomic studies on the honey bee, Apis mellifera, makes the standardization of qRT-PCR results an important task for ongoing community efforts. For this aim we selected four candidate reference genes (actin, ribosomal protein 49, elongation factor 1-alpha, tbp-association factor) and used three software-based approaches (geNorm, BestKeeper and NormFinder) to evaluate the suitability of these genes as endogenous controls. Their expression was examined during honey bee development, in different tissues, and after juvenile hormone exposure. Furthermore, the importance of choosing an appropriate reference gene was investigated for two developmentally regulated target genes. The results led us to consider all four candidate genes as suitable genes for normalization in A. mellifera. However, each condition evaluated in this study revealed a specific set of genes as the most appropriated ones.

Gene expression relationship between prostate cancer cells of Gleason 3, 4 and normal epithelial cells as revealed by cell type-specific transcriptomes

Background: Prostate cancer cells in primary tumors have been typed CD10(-)/CD13(-)/CD24(hi)/CD26(+)/CD38(lo)/CD44(-)/CD104(-). This CD phenotype suggests a lineage relationship between cancer cells and luminal cells. The Gleason grade of tumors is a descriptive of tumor glandular differentiation. Higher Gleason scores are associated with treatment failure. Methods: CD26(+) cancer cells were isolated from Gleason 3+3 (G3) and Gleason 4+4 (G4) tumors by cell sorting, and their gene expression or transcriptome was determined by Affymetrix DNA array analysis. Dataset analysis was used to determine gene expression similarities and differences between G3 and G4 as well as to prostate cancer cell lines and histologically normal prostate luminal cells. Results: The G3 and G4 transcriptomes were compared to those of prostatic cell types of non-cancer, which included luminal, basal, stromal fibromuscular, and endothelial. A principal components analysis of the various transcriptome datasets indicated a closer relationship between luminal and G3 than luminal and G4. Dataset comparison also showed that the cancer transcriptomes differed substantially from those of prostate cancer cell lines. Conclusions: Genes differentially expressed in cancer are potential biomarkers for cancer detection, and those differentially expressed between G3 and G4 are potential biomarkers for disease stratification given that G4 cancer is associated with poor outcomes. Differentially expressed genes likely contribute to the prostate cancer phenotype and constitute the signatures of these particular cancer cell types.

Phenotypic variation among four family members with congenital hypothyroidism caused by two distinct thyroglobulin gene mutations

Background: Thyroglobulin (Tg) is a large glycoprotein that is intimately involved in the biosynthesis of thyroxine and triiodothyronine. At least 38 mutations have been described in the Tg gene that are associated with varying degrees of hypothyroidism. We studied the Tg gene in four related subjects with congenital hypothyroidism. Summary: We found a novel compound heterozygous constellation (IVS30 + 1G>T/A2215D) in a brother and sister and one previously described related mutation (IVS30+1G>T) in their two sibling second degree cousins. The brother with the IVS30 + 1G>T/A2215D mutation and the two siblings with the IVS30+1G>T mutation had fetal or neonatal goiter and all had hypothyroidism. Conclusions: This study further confirms the association of the IVS30+1G>T mutation of the Tg gene with hypothyroidism. Computer analysis predicts that the A2215D mutation, first reported here, should cause structural instability of Tg but when present as a compound heterozygous mutation with IVS30+G>T/A its effect is unclear but is likely to be influenced by iodine intake.

Evaluating gene by sex and age interactions on cardiovascular risk factors in Brazilian families

Background: In family studies, it is important to evaluate the impact of genes and environmental factors on traits of interest. In particular, the relative influences of both genes and the environment may vary in different strata of the population of interest, such as young and old individuals, or males and females. Methods: In this paper, extensions of the variance components model are used to evaluate heterogeneity in the genetic and environmental variance components due to the effects of sex and age (the cutoff between young and old was 43 yrs). The data analyzed were from 81 Brazilian families (1,675 individuals) of the Baependi Family Heart Study. Results: The models allowing for heterogeneity of variance components by sex suggest that genetic and environmental variances are not different in males and females for diastolic blood pressure, LDL-cholesterol, and HDL-cholesterol, independent of the covariates included in the models. However, for systolic blood pressure, fasting glucose and triglycerides, the evidence for heterogeneity was dependent on the covariates in the model. For instance, in the presence of sex and age covariates, heterogeneity in the genetic variance component was suggested for fasting glucose. But, for systolic blood pressure, there was no evidence of heterogeneity in any of the two variance components. Except for the LDL-cholesterol, models allowing for heterogeneity by age provide evidence of heterogeneity in genetic variance for triglycerides and systolic and diastolic blood pressure. There was evidence of heterogeneity in environmental variance in fasting glucose and HDL-cholesterol. Conclusions: Our results suggest that heterogeneity in trait variances should not be ignored in the design and analyses of gene-finding studies involving these traits, as it may generate additional information about gene effects, and allow the investigation of more sophisticated models such as the model including sex-specific oligogenic variance components.

Random X Inactivation and Extensive Mosaicism in Human Placenta Revealed by Analysis of Allele-Specific Gene Expression along the X Chromosome

Imprinted inactivation of the paternal X chromosome in marsupials is the primordial mechanism of dosage compensation for X-linked genes between females and males in Therians. In Eutherian mammals, X chromosome inactivation (XCI) evolved into a random process in cells from the embryo proper, where either the maternal or paternal X can be inactivated. However, species like mouse and bovine maintained imprinted XCI exclusively in extraembryonic tissues. The existence of imprinted XCI in humans remains controversial, with studies based on the analyses of only one or two X-linked genes in different extraembryonic tissues. Here we readdress this issue in human term placenta by performing a robust analysis of allele-specific expression of 22 X-linked genes, including XIST, using 27 SNPs in transcribed regions. We show that XCI is random in human placenta, and that this organ is arranged in relatively large patches of cells with either maternal or paternal inactive X. In addition, this analysis indicated heterogeneous maintenance of gene silencing along the inactive X, which combined with the extensive mosaicism found in placenta, can explain the lack of agreement among previous studies. Our results illustrate the differences of XCI mechanism between humans and mice, and highlight the importance of addressing the issue of imprinted XCI in other species in order to understand the evolution of dosage compensation in placental mammals.

Activation of endogenous p53 by combined p19Arf gene transfer and nutlin-3 drug treatment modalities in the murine cell lines B16 and C6

Background: Reactivation of p53 by either gene transfer or pharmacologic approaches may compensate for loss of p19Arf or excess mdm2 expression, common events in melanoma and glioma. In our previous work, we constructed the pCLPG retroviral vector where transgene expression is controlled by p53 through a p53-responsive promoter. The use of this vector to introduce p19Arf into tumor cells that harbor p53wt should yield viral expression of p19Arf which, in turn, would activate the endogenous p53 and result in enhanced vector expression and tumor suppression. Since nutlin-3 can activate p53 by blocking its interaction with mdm2, we explored the possibility that the combination of p19Arf gene transfer and nutlin-3 drug treatment may provide an additive benefit in stimulating p53 function. Methods: B16 (mouse melanoma) and C6 (rat glioma) cell lines, which harbor p53wt, were transduced with pCLPGp19 and these were additionally treated with nutlin-3 or the DNA damaging agent, doxorubicin. Viral expression was confirmed by Western, Northern and immunofluorescence assays. p53 function was assessed by reporter gene activity provided by a p53-responsive construct. Alterations in proliferation and viability were measured by colony formation, growth curve, cell cycle and MTT assays. In an animal model, B16 cells were treated with the pCLPGp19 virus and/or drugs before subcutaneous injection in C57BL/6 mice, observation of tumor progression and histopathologic analyses. Results: Here we show that the functional activation of endogenous p53wt in B16 was particularly challenging, but accomplished when combined gene transfer and drug treatments were applied, resulting in increased transactivation by p53, marked cell cycle alteration and reduced viability in culture. In an animal model, B16 cells treated with both p19Arf and nutlin-3 yielded increased necrosis and decreased BrdU marking. In comparison, C6 cells were quite susceptible to either treatment, yet p53 was further activated by the combination of p19Arf and nutlin-3. Conclusions: To the best of our knowledge, this is the first study to apply both p19Arf and nutlin-3 for the stimulation of p53 activity. These results support the notion that a p53 responsive vector may prove to be an interesting gene transfer tool, especially when combined with p53- activating agents, for the treatment of tumors that retain wild-type p53.

Inference of gene regulatory networks from time series by Tsallis entropy

Background: The inference of gene regulatory networks (GRNs) from large-scale expression profiles is one of the most challenging problems of Systems Biology nowadays. Many techniques and models have been proposed for this task. However, it is not generally possible to recover the original topology with great accuracy, mainly due to the short time series data in face of the high complexity of the networks and the intrinsic noise of the expression measurements. In order to improve the accuracy of GRNs inference methods based on entropy (mutual information), a new criterion function is here proposed. Results: In this paper we introduce the use of generalized entropy proposed by Tsallis, for the inference of GRNs from time series expression profiles. The inference process is based on a feature selection approach and the conditional entropy is applied as criterion function. In order to assess the proposed methodology, the algorithm is applied to recover the network topology from temporal expressions generated by an artificial gene network (AGN) model as well as from the DREAM challenge. The adopted AGN is based on theoretical models of complex networks and its gene transference function is obtained from random drawing on the set of possible Boolean functions, thus creating its dynamics. On the other hand, DREAM time series data presents variation of network size and its topologies are based on real networks. The dynamics are generated by continuous differential equations with noise and perturbation. By adopting both data sources, it is possible to estimate the average quality of the inference with respect to different network topologies, transfer functions and network sizes. Conclusions: A remarkable improvement of accuracy was observed in the experimental results by reducing the number of false connections in the inferred topology by the non-Shannon entropy. The obtained best free parameter of the Tsallis entropy was on average in the range 2.5 <= q <= 3.5 (hence, subextensive entropy), which opens new perspectives for GRNs inference methods based on information theory and for investigation of the nonextensivity of such networks. The inference algorithm and criterion function proposed here were implemented and included in the DimReduction software, which is freely available at http://sourceforge.net/projects/dimreduction and http://code.google.com/p/dimreduction/.

Cuff-induced vascular intima thickening is influenced by titration of the Ace gene in mice

Lacchini S, Heimann AS, Evangelista FS, Cardoso L, Silva GJ, Krieger JE. Cuff-induced vascular intima thickening is influenced by titration of the Ace gene in mice. Physiol Genomics 37: 225-230, 2009. First published March 3, 2009; doi:10.1152/physiolgenomics.90288.2008.-We tested the hypothesis that small changes in angiotensin I-converting enzyme (ACE) expression can alter the vascular response to injury. Male mice containing one, two, three, and four copies of the Ace gene with no detectable vascular abnormality or changes in blood pressure were submitted to cuff-induced femoral artery injury. Femoral thickening was higher in 3- and 4-copy mice (42.4 +/- 4.3% and 45.7 +/- 6.5%, respectively) compared with 1- and 2-copy mice (8.3 +/- 1.3% and 8.5 +/- 0.9%, respectively). Femoral ACE levels from control and injured vessels were assessed in 1- and 3-copy Ace mice, which represent the extremes of the observed response. ACE vascular activity was higher in 3- vs. 1-copy Ace mice (2.4-fold, P < 0.05) in the control uninjured vessel. Upon injury, ACE activity significantly increased in both groups [2.41-fold and 2.14-fold (P < 0.05) for 1- and 3-copy groups, respectively] but reached higher levels in 3- vs. 1-copy Ace mice (P < 0.05). Pharmacological interventions were then used as a counterproof and to indirectly assess the role of angiotensin II (ANG II) on this response. Interestingly, ACE inhibition (enalapril) and ANG II AT(1) receptor blocker (losartan) reduced intima thickening in 3-copy mice to 1-copy mouse values (P < 0.05) while ANG II treatment significantly increased intima thickening in 1-copy mice to 3-copy mouse levels (P < 0.05). Together, these data indicate that small physiologically relevant changes in ACE, not associated with basal vascular abnormalities or blood pressure levels, do influence the magnitude of cuff-induced neointima thickening in mice.

Evaluation of the genetic variability of orchid fleck virus by single-strand conformational polymorphism analysis and nucleotide sequencing of a fragment from the nucleocapsid gene

The variability of a fragment of the nucleocapsid gene of orchid fleck virus (OFV) was investigated by single-strand conformational polymorphism (SSCP) analysis and nucleotide sequencing. Forty-eight samples of 18 genera of orchids were collected from Brazil, Costa Rica and Australia. The SSCP analysis yielded six different band patterns, and phylogenetic analysis based on the nucleotide fragment sequence obtained in this work and six available in GenBank showed two different groups, one with isolates 023Germany and So-Japan, and other with the rest of the isolates. None of the analyses showed geographic correlation among the Brazilian strains. The data obtained in this study showed a low genetic variation in this region of the genome; the d(N)/d(S) ratio of 0.251-0.405 demonstrated a negative selective pressure that maintains the stability of the analyzed fragments.

ABCB1 and ABCC1 expression in peripheral mononuclear cells is influenced by gene polymorphisms and atorvastatin treatment

This study investigated the effects of atorvastatin on ABCB1 and ABCC1 mRNA expression on peripheral blood mononuclear cells (PBMC) and their relationship with gene polymorphisms and lowering-cholesterol response. one hundred and thirty-six individuals with hypercholesterolemia were selected and treated with atorvastatin (10 mg/day/4 weeks). Blood samples were collected for serum lipids and apolipoproteins measurements and DNA and RNA extraction. ABCB1 (C3435T and G2677T/A) and ABCC1 (G2012T) gene polymorphisms were identified by polymerase chain reaction-restriction (PCR)-RFLP and mRNA expression was measured in peripheral blood mononuclear cells by singleplex real-time PCR. ABCB1 polymorphisms were associated with risk for coronary artery disease (CAD) (p < 0.05). After atorvastatin treatment, both ABCB1 and ABCC1 genes showed 50% reduction of the mRNA expression (p < 0.05). Reduction of ABCB1 expression was associated with ABCB1 G2677T/A polymorphism (p = 0.039). Basal ABCB1 mRNA in the lower quartile (<0.024) was associated with lower reduction rate of serum low-density lipoprotein (LDL) cholesterol (33.4 +/- 12.4%) and apolipoprotein B (apoB) (17.0 +/- 31.3%) when compared with the higher quartile (>0.085: LDL-c = 40.3 +/- 14.3%; apoB = 32.5 +/- 10.7%; p < 0.05). ABCB1 substrates or inhibitors did not affect the baseline expression, while ABCB1 inhibitors reversed the effects of atorvastatin on both ABCB1 and ABCC1 transporters. In conclusion, ABCB1 and ABCC1 mRNA levels in PBMC are modulated by atorvastatin and ABCB1 G2677T/A polymorphism. and ABCB1 baseline expression is related to differences in serum LDL cholesterol and apoB in response to atorvastatin. (C) 2008 Elsevier Inc. All rights reserved.

Silencing of mitochondrial alternative oxidase gene of Aspergillus fumigatus enhances reactive oxygen species production and killing of the fungus by macrophages

We previously demonstrated that conidia from Aspergillus fumigatus incubated with menadione and paraquat increases activity and expression of cyanide-insensitive alternative oxidase (AOX). Here, we employed the RNA silencing technique in A. fumigatus using the vector pALB1/aoxAf in order to down-regulate the aox gene. Positive transformants for aox gene silencing of A. fumigatus were more susceptible both to an imposed in vitro oxidative stress condition and to macrophages killing, suggesting that AOX is required for the A. fumigatus pathogenicity, mainly for the survival of the fungus conidia during host infection and resistance to reactive oxygen species generated by macrophages.

Transcription regulation of the PbGP43 gene by nitrogen in the human pathogen Paracoccidioides brasiliensis

We show indirect evidences for the possible involvement of NIT-2-like binding motifs in transcription modulation of the PbGP43 gene, which codes for an important antigen from the human fungal pathogen Paracoccidioides brasiliensis. This investigation was motivated by the finding of 23 NIT2-like sites within the proximal -2047 nucleotides of the PbGP43 5` intergenic region from the Pb339 isolate. They compose four clusters, two of them identical. We found four NIT2-containing probes that were positive in electrophoretic mobility shift assays and further analyzed them. PbGP43 could be modulated by nitrogen primary sources in Pb339, Pb3 and Pb18 isolates, as observed by reverse transcription (RT) real time-PCR. Gene reporter assays conducted in Aspergillus nidulans suggested that the minimal fragment responsible for nitrogen modulation lies within -480 bp of the PbGP43 gene. This is the first report on PbGP43 transcription modulation in response to nitrogen primary sources, which might help understand its regulation during infection. (C) 2008 Elsevier Inc. All rights reserved.

Concordance between dispersal and mitochondrial gene flow: isolation by distance in a tropical teleost, Lates calcarifer (Australian barramundi)

Patterns of population subdivision and the relationship between gene flow and geographical distance in the tropical estuarine fish Lares calcarifer (Centropomidae) were investigated using mtDNA control region sequences. Sixty-three putative haplotypes were resolved from a total of 270 individuals from nine localities within three geographical regions spanning the north Australian coastline. Despite a continuous estuarine distribution throughout the sampled range, no haplotypes were shared among regions. However, within regions, common haplotypes were often shared among localities. Both sequence-based (average Phi(ST)=0.328) and haplotype-based (average Phi(ST)=0.182) population subdivision analyses indicated strong geographical structuring. Depending on the method of calculation, geographical distance explained either 79 per cent (sequence-based) or 23 per cent (haplotype-based) of the variation in mitochondrial gene flow. Such relationships suggest that genetic differentiation of L. calcarifer has been generated via isolation-by-distance, possibly in a stepping-stone fashion. This pattern of genetic structure is concordant with expectations based on the life history of L. calcarifer and direct studies of its dispersal patterns. Mitochondrial DNA variation, although generally in agreement with patterns of allozyme variation, detected population subdivision at smaller spatial scales. Our analysis of mtDNA variation in L. calcarifer confirms that population genetic models can detect population structure of not only evolutionary significance but also of demographic significance. Further, it demonstrates the power of inferring such structure from hypervariable markers, which correspond to small effective population sizes.