Predicting cancer involvement of genes from heterogeneous data

Background: Systematic approaches for identifying proteins involved in different types of cancer are needed. Experimental techniques such as microarrays are being used to characterize cancer, but validating their results can be a laborious task. Computational approaches are used to prioritize between genes putatively involved in cancer, usually based on further analyzing experimental data. Results: We implemented a systematic method using the PIANA software that predicts cancer involvement of genes by integrating heterogeneous datasets. Specifically, we produced lists of genes likely to be involved in cancer by relying on: (i) protein-protein interactions; (ii) differential expression data; and (iii) structural and functional properties of cancer genes. The integrative approach that combines multiple sources of data obtained positive predictive values ranging from 23% (on a list of 811 genes) to 73% (on a list of 22 genes), outperforming the use of any of the data sources alone. We analyze a list of 20 cancer gene predictions, finding that most of them have been recently linked to cancer in literature. Conclusion: Our approach to identifying and prioritizing candidate cancer genes can be used to produce lists of genes likely to be involved in cancer. Our results suggest that differential expression studies yielding high numbers of candidate cancer genes can be filtered using protein interaction networks.

OSIRISv1.2: a named entity recognition system for sequence variants of genes in biomedical literature

Background: Single Nucleotide Polymorphisms, among other type of sequence variants, constitute key elements in genetic epidemiology and pharmacogenomics. While sequence data about genetic variation is found at databases such as dbSNP, clues about the functional and phenotypic consequences of the variations are generally found in biomedical literature. The identification of the relevant documents and the extraction of the information from them are hampered by the large size of literature databases and the lack of widely accepted standard notation for biomedical entities. Thus, automatic systems for the identification of citations of allelic variants of genes in biomedical texts are required. Results: Our group has previously reported the development of OSIRIS, a system aimed at the retrieval of literature about allelic variants of genes http://ibi.imim.es/osirisform.html. Here we describe the development of a new version of OSIRIS (OSIRISv1.2, http://ibi.imim.es/OSIRISv1.2.html webcite) which incorporates a new entity recognition module and is built on top of a local mirror of the MEDLINE collection and HgenetInfoDB: a database that collects data on human gene sequence variations. The new entity recognition module is based on a pattern-based search algorithm for the identification of variation terms in the texts and their mapping to dbSNP identifiers. The performance of OSIRISv1.2 was evaluated on a manually annotated corpus, resulting in 99% precision, 82% recall, and an F-score of 0.89. As an example, the application of the system for collecting literature citations for the allelic variants of genes related to the diseases intracranial aneurysm and breast cancer is presented. Conclusion: OSIRISv1.2 can be used to link literature references to dbSNP database entries with high accuracy, and therefore is suitable for collecting current knowledge on gene sequence variations and supporting the functional annotation of variation databases. The application of OSIRISv1.2 in combination with controlled vocabularies like MeSH provides a way to identify associations of biomedical interest, such as those that relate SNPs with diseases.

Tumor necrosis factor and transforming growth factor β regulate clock genes by controlling the expression of the cold inducible RNA-binding protein (CIRBP).

The circadian clock drives the rhythmic expression of a broad array of genes that orchestrate metabolism, sleep wake behavior, and the immune response. Clock genes are transcriptional regulators engaged in the generation of circadian rhythms. The cold inducible RNA-binding protein (CIRBP) guarantees high amplitude expression of clock. The cytokines TNF and TGFβ impair the expression of clock genes, namely the period genes and the proline- and acidic amino acid-rich basic leucine zipper (PAR-bZip) clock-controlled genes. Here, we show that TNF and TGFβ impair the expression of Cirbp in fibroblasts and neuronal cells. IL-1β, IL-6, IFNα, and IFNγ do not exert such effects. Depletion of Cirbp is found to increase the susceptibility of cells to the TNF-mediated inhibition of high amplitude expression of clock genes and modulates the TNF-induced cytokine response. Our findings reveal a new mechanism of cytokine-regulated expression of clock genes.

Genetik der essentiellen arteriellen Hypertonie [Essential hypertension - which genes can be held responsible for it?].

Hypertension is a common heritable cardiovascular risk factor. Some rare monogenic forms of hypertension have been described, but the majority of patients suffer from essential hypertension, for whom the underlying genetic mechanisms are not clear. Essential hypertension is a complex trait, involving multiple genes and environmental factors. Recently, progress in the identification of common genetic variants associated with essential hypertension has been made due to large-scale international collaborative projects. In this article we review the new research methods used as well as selected recent findings in this field.

Repeated exposure to Lutzomyia intermedia sand fly saliva induces local expression of interferon-inducible genes both at the site of injection in mice and in human blood.

During a blood meal, Lutzomyia intermedia sand flies transmit Leishmania braziliensis, a parasite causing tegumentary leishmaniasis. In experimental leishmaniasis, pre-exposure to saliva of most blood-feeding sand flies results in parasite establishment in absence of any skin damages in mice challenged with dermotropic Leishmania species together with saliva. In contrast, pre-immunization with Lu. intermedia salivary gland sonicate (SGS) results in enhanced skin inflammatory exacerbation upon co-inoculation of Lu. intermedia SGS and L. braziliensis. These data highlight potential unique features of both L. braziliensis and Lu. intermedia. In this study, we investigated the genes modulated by Lu. intermedia SGS immunization to understand their potential impact on the subsequent cutaneous immune response following inoculation of both SGS and L. braziliensis. The cellular recruitment and global gene expression profile was analyzed in mice repeatedly inoculated or not with Lu. intermedia. Microarray gene analysis revealed the upregulation of a distinct set of IFN-inducible genes, an immune signature not seen to the same extent in control animals. Of note this INF-inducible gene set was not induced in SGS pre-immunized mice subsequently co-inoculated with SGS and L. braziliensis. These data suggest the parasite prevented the upregulation of this Lu. intermedia saliva-related immune signature. The presence of these IFN-inducible genes was further analyzed in peripheral blood mononuclear cells (PBMCs) sampled from uninfected human individuals living in a L. braziliensis-endemic region of Brazil thus regularly exposed to Lu. intermedia bites. PBMCs were cultured in presence or absence of Lu. intermedia SGS. Using qRT-PCR we established that the IFN-inducible genes induced in the skin of SGS pre-immunized mice, were also upregulated by SGS in PBMCs from human individuals regularly exposed to Lu. intermedia bites, but not in PBMCs of control subjects. These data demonstrate that repeated exposure to Lu. intermedia SGS induces the expression of potentially host-protective IFN-inducible genes.

Association of candidate genes with phenotypic traits relevant to anorexia nervosa.

This analysis is a follow-up to an earlier investigation of 182 genes selected as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN). As those initial case-control results revealed no statistically significant differences in single nucleotide polymorphisms, herein, we investigate alternative phenotypes associated with AN. In 1762 females, using regression analyses, we examined the following: (i) lowest illness-related attained body mass index; (ii) age at menarche; (iii) drive for thinness; (iv) body dissatisfaction; (v) trait anxiety; (vi) concern over mistakes; and (vii) the anticipatory worry and pessimism versus uninhibited optimism subscale of the harm avoidance scale. After controlling for multiple comparisons, no statistically significant results emerged. Although results must be viewed in the context of limitations of statistical power, the approach illustrates a means of potentially identifying genetic variants conferring susceptibility to AN because less complex phenotypes associated with AN are more proximal to the genotype and may be influenced by fewer genes. Copyright © 2011 John Wiley & Sons, Ltd and Eating Disorders Association.

Identifying protein-coding genes in genomic sequences.

The vast majority of the biology of a newly sequenced genome is inferred from the set of encoded proteins. Predicting this set is therefore invariably the first step after the completion of the genome DNA sequence. Here we review the main computational pipelines used to generate the human reference protein-coding gene sets.

Ferripyochelin uptake genes are involved in pyochelin-mediated signalling in Pseudomonas aeruginosa.

In response to iron starvation, Pseudomonas aeruginosa produces the siderophore pyochelin. When secreted to the extracellular environment, pyochelin chelates iron and transports it to the bacterial cytoplasm via its specific outer-membrane receptor FptA and the inner-membrane permease FptX. Exogenously added pyochelin also acts as a signal which induces the expression of the pyochelin biosynthesis and uptake genes by activating PchR, a cytoplasmic regulatory protein of the AraC/XylS family. The importance of ferripyochelin uptake genes in this regulation was evaluated. The fptA and fptX genes were shown to be part of the fptABCX ferripyochelin transport operon, which is conserved in Burkholderia sp. and Rhodospirillum rubrum. The fptB and fptC genes were found to be dispensable for utilization of pyochelin as an iron source, for signalling and for pyochelin production. By contrast, mutations in fptA and fptX not only interfered with pyochelin utilization, but also affected signalling and diminished siderophore production. It is concluded from this that pyochelin-mediated signalling operates to a large extent via the ferripyochelin transport system.

Association study of 182 candidate genes in anorexia nervosa.

We performed association studies with 5,151 SNPs that were judged as likely candidate genetic variations conferring susceptibility to anorexia nervosa (AN) based on location under reported linkage peaks, previous results in the literature (182 candidate genes), brain expression, biological plausibility, and estrogen responsivity. We employed a case-control design that tested each SNP individually as well as haplotypes derived from these SNPs in 1,085 case individuals with AN diagnoses and 677 control individuals. We also performed separate association analyses using three increasingly restrictive case definitions for AN: all individuals with any subtype of AN (All AN: n = 1,085); individuals with AN with no binge eating behavior (AN with No Binge Eating: n = 687); and individuals with the restricting subtype of AN (Restricting AN: n = 421). After accounting for multiple comparisons, there were no statistically significant associations for any individual SNP or haplotype block with any definition of illness. These results underscore the importance of large samples to yield appropriate power to detect genotypic differences in individuals with AN and also motivate complementary approaches involving Genome-Wide Association (GWA) studies, Copy Number Variation (CNV) analyses, sequencing-based rare variant discovery assays, and pathway-based analysis in order to make up for deficiencies in traditional candidate gene approaches to AN.

Sex chromosomes and male functions: where do new genes go?

The position of a gene in the genome may have important consequences for its function. Therefore, when a new duplicate gene arises, its location may be critical in determining its fate. Our recent work in humans, mouse, and Drosophila provided a test by studying the patterns of duplication in sex chromosome evolution. We revealed a bias in the generation and recruitment of new gene copies involving the X chromosome that has been shaped largely by selection for male germline functions. The gene movement patterns we observed reflect an ongoing process as some of the new genes are very young while others were present before the divergence of humans and mouse. This suggests a continuing redistribution of male-related genes to achieve a more efficient allocation of male functions. This notion should be further tested in organisms employing other sex determination systems or in organisms differing in germline sex chromosome inactivation. It is likely that the selective forces that were detected in these studies are also acting on other types of duplicate genes. As a result, future work elucidating sex chromosome differentiation by other mutational mechanisms will shed light on this important process.

GLUTAMATE RECEPTOR-LIKE genes mediate leaf-to-leaf wound signalling.

Wounded leaves communicate their damage status to one another through a poorly understood process of long-distance signalling. This stimulates the distal production of jasmonates, potent regulators of defence responses. Using non-invasive electrodes we mapped surface potential changes in Arabidopsis thaliana after wounding leaf eight and found that membrane depolarizations correlated with jasmonate signalling domains in undamaged leaves. Furthermore, current injection elicited jasmonoyl-isoleucine accumulation, resulting in a transcriptome enriched in RNAs encoding key jasmonate signalling regulators. From among 34 screened membrane protein mutant lines, mutations in several clade 3 GLUTAMATE RECEPTOR-LIKE genes (GLRs 3.2, 3.3 and 3.6) attenuated wound-induced surface potential changes. Jasmonate-response gene expression in leaves distal to wounds was reduced in a glr3.3 glr3.6 double mutant. This work provides a genetic basis for investigating mechanisms of long-distance wound signalling in plants and indicates that plant genes related to those important for synaptic activity in animals function in organ-to-organ wound signalling.

Reverse transcription quantitative real-time polymerase chain reaction reference genes in the spared nerve injury model of neuropathic pain: validation and literature search.

BACKGROUND: The reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) is a widely used, highly sensitive laboratory technique to rapidly and easily detect, identify and quantify gene expression. Reliable RT-qPCR data necessitates accurate normalization with validated control genes (reference genes) whose expression is constant in all studied conditions. This stability has to be demonstrated.We performed a literature search for studies using quantitative or semi-quantitative PCR in the rat spared nerve injury (SNI) model of neuropathic pain to verify whether any reference genes had previously been validated. We then analyzed the stability over time of 7 commonly used reference genes in the nervous system - specifically in the spinal cord dorsal horn and the dorsal root ganglion (DRG). These were: Actin beta (Actb), Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal proteins 18S (18S), L13a (RPL13a) and L29 (RPL29), hypoxanthine phosphoribosyltransferase 1 (HPRT1) and hydroxymethylbilane synthase (HMBS). We compared the candidate genes and established a stability ranking using the geNorm algorithm. Finally, we assessed the number of reference genes necessary for accurate normalization in this neuropathic pain model. RESULTS: We found GAPDH, HMBS, Actb, HPRT1 and 18S cited as reference genes in literature on studies using the SNI model. Only HPRT1 and 18S had been once previously demonstrated as stable in RT-qPCR arrays. All the genes tested in this study, using the geNorm algorithm, presented gene stability values (M-value) acceptable enough for them to qualify as potential reference genes in both DRG and spinal cord. Using the coefficient of variation, 18S failed the 50% cut-off with a value of 61% in the DRG. The two most stable genes in the dorsal horn were RPL29 and RPL13a; in the DRG they were HPRT1 and Actb. Using a 0.15 cut-off for pairwise variations we found that any pair of stable reference gene was sufficient for the normalization process. CONCLUSIONS: In the rat SNI model, we validated and ranked Actb, RPL29, RPL13a, HMBS, GAPDH, HPRT1 and 18S as good reference genes in the spinal cord. In the DRG, 18S did not fulfill stability criteria. The combination of any two stable reference genes was sufficient to provide an accurate normalization.

Genes involved in the astrocyte-neuron lactate shuttle (ANLS) are specifically regulated in cortical astrocytes following sleep deprivation in mice.

STUDY OBJECTIVES: There is growing evidence indicating that in order to meet the neuronal energy demands, astrocytes provide lactate as an energy substrate for neurons through a mechanism called "astrocyte-neuron lactate shuttle" (ANLS). Since neuronal activity changes dramatically during vigilance states, we hypothesized that the ANLS may be regulated during the sleep-wake cycle. To test this hypothesis we investigated the expression of genes associated with the ANLS specifically in astrocytes following sleep deprivation. Astrocytes were purified by fluorescence-activated cell sorting from transgenic mice expressing the green fluorescent protein (GFP) under the control of the human astrocytic GFAP-promoter. DESIGN: 6-hour instrumental sleep deprivation (TSD). SETTING: Animal sleep research laboratory. PARTICIPANTS: Young (P23-P27) FVB/N-Tg (GFAP-GFP) 14Mes/J (Tg) mice of both sexes and 7-8 week male Tg and FVB/Nj mice. INTERVENTIONS: Basal sleep recordings and sleep deprivation achieved using a modified cage where animals were gently forced to move. MEASUREMENTS AND RESULTS: Since Tg and FVB/Nj mice displayed a similar sleep-wake pattern, we performed a TSD in young Tg mice. Total RNA was extracted from the GFP-positive and GFP-negative cells sorted from cerebral cortex. Quantitative RT-PCR analysis showed that levels of Glut1, α-2-Na/K pump, Glt1, and Ldha mRNAs were significantly increased following TSD in GFP-positive cells. In GFP-negative cells, a tendency to increase, although not significant, was observed for Ldha, Mct2, and α-3-Na/K pump mRNAs. CONCLUSIONS: This study shows that TSD induces the expression of genes associated with ANLS specifically in astrocytes, underlying the important role of astrocytes in the maintenance of the neuro-metabolic coupling across the sleep-wake cycle.

A novel bioinformatics pipeline to discover genes related to arbuscular mycorrhizal symbiosis based on their evolutionary conservation pattern among higher plants.

BACKGROUND: Genes involved in arbuscular mycorrhizal (AM) symbiosis have been identified primarily by mutant screens, followed by identification of the mutated genes (forward genetics). In addition, a number of AM-related genes has been identified by their AM-related expression patterns, and their function has subsequently been elucidated by knock-down or knock-out approaches (reverse genetics). However, genes that are members of functionally redundant gene families, or genes that have a vital function and therefore result in lethal mutant phenotypes, are difficult to identify. If such genes are constitutively expressed and therefore escape differential expression analyses, they remain elusive. The goal of this study was to systematically search for AM-related genes with a bioinformatics strategy that is insensitive to these problems. The central element of our approach is based on the fact that many AM-related genes are conserved only among AM-competent species. RESULTS: Our approach involves genome-wide comparisons at the proteome level of AM-competent host species with non-mycorrhizal species. Using a clustering method we first established orthologous/paralogous relationships and subsequently identified protein clusters that contain members only of the AM-competent species. Proteins of these clusters were then analyzed in an extended set of 16 plant species and ranked based on their relatedness among AM-competent monocot and dicot species, relative to non-mycorrhizal species. In addition, we combined the information on the protein-coding sequence with gene expression data and with promoter analysis. As a result we present a list of yet uncharacterized proteins that show a strongly AM-related pattern of sequence conservation, indicating that the respective genes may have been under selection for a function in AM. Among the top candidates are three genes that encode a small family of similar receptor-like kinases that are related to the S-locus receptor kinases involved in sporophytic self-incompatibility. CONCLUSIONS: We present a new systematic strategy of gene discovery based on conservation of the protein-coding sequence that complements classical forward and reverse genetics. This strategy can be applied to diverse other biological phenomena if species with established genome sequences fall into distinguished groups that differ in a defined functional trait of interest.

GacA-controlled activation of promoters for small RNA genes in Pseudomonas fluorescens.

The Gac/Rsm signal transduction pathway positively regulates secondary metabolism, production of extracellular enzymes, and biocontrol properties of Pseudomonas fluorescens CHA0 via the expression of three noncoding small RNAs, termed RsmX, RsmY, and RsmZ. The architecture and function of the rsmY and rsmZ promoters were studied in vivo. A conserved palindromic upstream activating sequence (UAS) was found to be necessary but not sufficient for rsmY and rsmZ expression and for activation by the response regulator GacA. A poorly conserved linker region located between the UAS and the -10 promoter sequence was also essential for GacA-dependent rsmY and rsmZ expression, suggesting a need for auxiliary transcription factors. One such factor involved in the activation of the rsmZ promoter was identified as the PsrA protein, previously recognized as an activator of the rpoS gene and a repressor of fatty acid degradation. Furthermore, the integration host factor (IHF) protein was found to bind with high affinity to the rsmZ promoter region in vitro, suggesting that DNA bending contributes to the regulated expression of rsmZ. In an rsmXYZ triple mutant, the expression of rsmY and rsmZ was elevated above that found in the wild type. This negative feedback loop appears to involve the translational regulators RsmA and RsmE, whose activity is antagonized by RsmXYZ, and several hypothetical DNA-binding proteins. This highly complex network controls the expression of the three small RNAs in response to cell physiology and cell population densities.