Predicting functional gene links from phylogenetic-statistical analyses of whole genomes

An important element of the developing field of proteomics is to understand protein-protein interactions and other functional links amongst genes. Across-species correlation methods for detecting functional links work on the premise that functionally linked proteins will tend to show a common pattern of presence and absence across a range of genomes. We describe a maximum likelihood statistical model for predicting functional gene linkages. The method detects independent instances of the correlated gain or loss of pairs of proteins on phylogenetic trees, reducing the high rates of false positives observed in conventional across-species methods that do not explicitly incorporate a phylogeny. We show, in a dataset of 10,551 protein pairs, that the phylogenetic method improves by up to 35% on across-species analyses at identifying known functionally linked proteins. The method shows that protein pairs with at least two to three correlated events of gain or loss are almost certainly functionally linked. Contingent evolution, in which one gene's presence or absence depends upon the presence of another, can also be detected phylogenetically, and may identify genes whose functional significance depends upon its interaction with other genes. Incorporating phylogenetic information improves the prediction of functional linkages. The improvement derives from having a lower rate of false positives and from detecting trends that across-species analyses miss. Phylogenetic methods can easily be incorporated into the screening of large-scale bioinformatics datasets to identify sets of protein links and to characterise gene networks.

Constrained models of evolution lead to improved prediction of functional linkage from correlated gain and loss of genes

Motivation: We compare phylogenetic approaches for inferring functional gene links. The approaches detect independent instances of the correlated gain and loss of pairs of genes from species' genomes. We investigate the effect on results of basing evidence of correlations on two phylogenetic approaches, Dollo parsminony and maximum likelihood (ML). We further examine the effect of constraining the ML model by fixing the rate of gene gain at a low value, rather than estimating it from the data. Results: We detect correlated evolution among a test set of pairs of yeast (Saccharomyces cerevisiae) genes, with a case study of 21 eukaryotic genomes and test data derived from known yeast protein complexes. If the rate at which genes are gained is constrained to be low, ML achieves by far the best results at detecting known functional links. The model then has fewer parameters but it is more realistic by preventing genes from being gained more than once. Availability: BayesTraits by M. Pagel and A. Meade, and a script to configure and repeatedly launch it by D. Barker and M. Pagel, are available at http://www.evolution.reading.ac.uk .

Characterization of para-Nitrophenol-degrading bacterial communities in river water by using functional markers and stable isotope probing

Microbial degradation is a major determinant of the fate of pollutants in the environment. para-Nitrophenol (PNP) is an EPA listed priority pollutant with a wide environmental distribution, but little is known about the microorganisms that degrade it in the environment. We studied the diversity of active PNP-degrading bacterial populations in river water using a novel functional marker approach coupled with [13C6]PNP stable isotope probing (SIP). Culturing together with culture-independent terminal restriction fragment length polymorphism analysis of 16S rRNA gene amplicons identified Pseudomonas syringae to be the major driver of PNP degradation in river water microcosms. This was confirmed by SIP-pyrosequencing of amplified 16S rRNA. Similarly, functional gene analysis showed that degradation followed the Gram-negative bacterial pathway and involved pnpA from Pseudomonas spp. However, analysis of maleylacetate reductase (encoded by mar), an enzyme common to late stages of both Gram-negative and Gram-positive bacterial PNP degradation pathways, identified a diverse assemblage of bacteria associated with PNP degradation, suggesting that mar has limited use as a specific marker of PNP biodegradation. Both the pnpA and mar genes were detected in a PNP-degrading isolate, P. syringae AKHD2, which was isolated from river water. Our results suggest that PNP-degrading cultures of Pseudomonas spp. are representative of environmental PNP-degrading populations.

Variations in the human cannabinoid receptor (CNR1) gene modulate striatal responses to happy faces.

Happy facial expressions are innate social rewards and evoke a response in the striatum, a region known for its role in reward processing in rats, primates and humans. The cannabinoid receptor 1 (CNR1) is the best-characterized molecule of the endocannabinoid system, involved in processing rewards. We hypothesized that genetic variation in human CNR1 gene would predict differences in the striatal response to happy faces. In a 3T functional magnetic resonance imaging (fMRI) scanning study on 19 Caucasian volunteers, we report that four single nucleotide polymorphisms (SNPs) in the CNR1 locus modulate differential striatal response to happy but not to disgust faces. This suggests a role for the variations of the CNR1 gene in underlying social reward responsivity. Future studies should aim to replicate this finding with a balanced design in a larger sample, but these preliminary results suggest neural responsivity to emotional and socially rewarding stimuli varies as a function of CNR1 genotype. This has implications for medical conditions involving hypo-responsivity to emotional and social stimuli, such as autism.

Use of global gene expression patterns in mechanistic studies of oestrogen action in MCF7 human breast cancer cells

Over the years, the MCF7 human breast cancer cell line has provided a model system for the study of cellular and molecular mechanisms in oestrogen regulation of cell proliferation and in progression to oestrogen and antioestrogen independent growth. Global gene expression profiling has shown that oestrogen action in MCF7 cells involves the coordinated regulation of hundreds of genes across a wide range of functional groupings and that more genes are down regulated than upregulated. Adaptation to long-term oestrogen deprivation, which results in loss of oestrogen-responsive growth, involves alterations to gene patterns not only at early time points (0-4 weeks) but continuing through to later times (20-55 weeks), and even involves alterations to patterns of oestrogen-regulated gene expression. Only 48% of the genes which were regulated >= 2-fold by oestradiol in oestrogen-responsive cells retained this responsiveness after long-term oestrogen deprivation but other genes developed de novo oestrogen regulation. Long-term exposure to fulvestrant, which resulted in loss of growth inhibition by the antioestrogen, resulted in some very large fold changes in gene expression up to 10,000-fold. Comparison of gene profiles produced by environmental chemicals with oestrogenic properties showed that each ligand gave its own unique expression profile which suggests that environmental oestrogens entering the human breast may give rise to a more complex web of interference in cell function than simply mimicking oestrogen action at inappropriate times. (C) 2009 Elsevier Ltd. All rights reserved.

EfeUOB (YcdNOB) is a tripartite, acid-induced and CpxAR-regulated, low-pH Fe2+ transporter that is cryptic in Escherichia coli K-12 but functional in E-coli O157 : H7

Escherichia coli possesses iron transporters specific for either Fe2+ or Fe3+. Although Fe2+ is far more soluble than Fe3+, it rapidly oxidizes aerobically at pH >= 7. Thus, FeoAB, the major Fe2+ transporter of E. coli, operates anaerobically. However, Fe2+ remains stable aerobically under acidic conditions, although a low-pH Fe2+ importer has not been previously identified. Here we show that ycdNOB (efeUOB) specifies the first such transporter. efeUOB is repressed at high pH by CpxAR, and is Fe2+-Fur repressed. EfeU is homologous to the high-affinity iron permease, Ftr1p, of Saccharomyces cerevisiae and other fungi. EfeO is periplasmic with a cupredoxin N-terminal domain; EfeB is also periplasmic and is haem peroxidase-like. All three Efe proteins are required for Efe function. The efeU gene of E. coli K-12 is cryptic due to a frameshift mutation - repair of the single-base-pair deletion generates a functional EfeUOB system. In contrast, the efeUOB operon of the enterohaemorrhagic strain, O157:1147, lacks any frameshift and is functional. A 'wild-type' K-12 strain bearing a functional EfeUOB displays a major growth advantage under aerobic, low-pH, low-iron conditions when a competing metal is provided. Fe-55 transport assays confirm the ferrous iron specificity of EfeUOB.

Rapid parallel expression in E-coli and insect cells: Analysis of five lef gene products of the Autographa californica multiple nuclear polyhedrosis virus (AcMNPV)

A number of strategies are emerging for the high throughput (HTP) expression of recombinant proteins to enable structural and functional study. Here we describe a workable HTP strategy based on parallel protein expression in E. coli and insect cells. Using this system we provide comparative expression data for five proteins derived from the Autographa californica polyhedrosis virus genome that vary in amino acid composition and in molecular weight. Although the proteins are part of a set of factors known to be required for viral late gene expression, the precise function of three of the five, late expression factors (lefs) 6, 7 and 10, is unknown. Rapid expression and characterisation has allowed the determination of their ability to bind DNA and shown a cellular location consistent with their properties. Our data point to the utility of a parallel expression strategy to rapidly obtain workable protein expression levels from many open reading frames (ORFs).

ms1, a novel stress-responsive, muscle-specific gene that is up-regulated in the early stages of pressure overload-induced left ventricular hypertrophy

We have identified and characterised a cDNA encoding a novel gene, designated myocyte stress 1 (ms1), that is up-regulated within 1 h in the left ventricle following the application of pressure overload by aortic banding in the rat. The deduced ms1 protein of 317 amino acids contains several putative functional motifs, including a region that is evolutionarily conserved. Distribution analysis indicates that rat ms1 mRNA expression is predominantly expressed in striated muscle and progressively increases in the left ventricle from embryo to adulthood. These findings suggest that rust may be important in striated muscle biology and the development of pressure-induced left ventricular hypertrophy. (C) 2002 Published by Elsevier Science B.V. on behalf of the Federation of European Biochemical Societies.

Influence of the Photorhabdus luminescens phosphomannose isomerase gene, manA, on mannose utilization, exopolysaccharide structure, and biofilm formation

Extracellular polysaccharide (EPS) is produced by diverse bacterial pathogens and fulfills assorted roles, including providing a structural matrix for biofilm formation and more specific functions in virulence, such as protection against immune defenses. We report here the first investigation of some of the genes important for biofilm formation in Photorhabdus luminescens and demonstrate the key role of the phosphomannose isomerase gene, manA, in the structure of functional EPS. Phenotypic analyses of a manA-deficient mutant showed the importance of EPS in motility, insect virulence, and biofilm formation on abiotic surfaces as well as the requirement of this gene for the use of mannose as the sole carbon source. Conversely, this defect had no apparent impact on symbiosis with the heterorhabditid nematode vector. A more detailed analysis of biofilm formation revealed that the manA mutant was able to attach to surfaces with the same efficiency as that of the wild-type strain but could not develop the more extended biofilm matrix structures. A compositional analysis of P. luminescens EPS reveals how the manA mutation has a major effect on the formation of a complete, branched EPS.

Assessment of motor symptoms and functional impact in prodromal and early Huntington disease.

The Functional Rating Scale Taskforce for pre-Huntington Disease (FuRST-pHD) is a multinational, multidisciplinary initiative with the goal of developing a data-driven, comprehensive, psychometrically sound, rating scale for assessing symptoms and functional ability in prodromal and early Huntington disease (HD) gene expansion carriers. The process involves input from numerous sources to identify relevant symptom domains, including HD individuals, caregivers, and experts from a variety of fields, as well as knowledge gained from the analysis of data from ongoing large-scale studies in HD using existing clinical scales. This is an iterative process in which an ongoing series of field tests in prodromal (prHD) and early HD individuals provides the team with data on which to make decisions regarding which questions should undergo further development or testing and which should be excluded. We report here the development and assessment of the first iteration of interview questions aimed to assess functional impact of motor manifestations in prHD and early HD individuals.

Prediction and functional analysis of native disorder in proteins from the Three Kingdoms of Life

An automatic method for recognizing natively disordered regions from amino acid sequence is described and benchmarked against predictors that were assessed at the latest critical assessment of techniques for protein structure prediction (CASP) experiment. The method attains a Wilcoxon score of 90.0, which represents a statistically significant improvement on the methods evaluated on the same targets at CASP. The classifier, DISOPRED2, was used to estimate the frequency of native disorder in several representative genomes from the three kingdoms of life. Putative, long (>30 residue) disordered segments are found to occur in 2.0% of archaean, 4.2% of eubacterial and 33.0% of eukaryotic proteins. The function of proteins with long predicted regions of disorder was investigated using the gene ontology annotations supplied with the Saccharomyces genome database. The analysis of the yeast proteome suggests that proteins containing disorder are often located in the cell nucleus and are involved in the regulation of transcription and cell signalling. The results also indicate that native disorder is associated with the molecular functions of kinase activity and nucleic acid binding.

Profiling dehydrin gene sequence and physiological parameters in drought tolerant and susceptible spring wheat cultivars.

Physiological and yield traits such as stomatal conductance (mmol m-2s-1), Leaf relative water content (RWC %) and grain yield per plant were studied in a separate experiment. Results revealed that five out of sixteen cultivars viz. Anmol, Moomal, Sarsabz, Bhitai and Pavan, appeared to be relatively more drought tolerant. Based on morphophysiological results, studies were continued to look at these cultivars for drought tolerance at molecular level. Initially, four well recognized primers for dehydrin genes (DHNs) responsible for drought induction in T. durum L., T. aestivum L. and O. sativa L. were used for profiling gene sequence of sixteen wheat cultivars. The primers amplified the DHN genes variably like Primer WDHN13 (T. aestivum L.) amplified the DHN gene in only seven cultivars whereas primer TdDHN15 (T. durum L.) amplified all the sixteen cultivars with even different DNA banding patterns some showing second weaker DNA bands. Third primer TdDHN16 (T. durum L.) has shown entirely different PCR amplification prototype, specially showing two strong DNA bands while fourth primer RAB16C (O. sativa L.) failed to amplify DHN gene in any of the cultivars. Examination of DNA sequences revealed several interesting features. First, it identified the two exon/one intron structure of this gene (complete sequences were not shown), a feature not previously described in the two database cDNA sequences available from T. aestivum L. (gi|21850). Secondly, the analysis identified several single nucleotide polymorphisms (SNPs), positions in gene sequence. Although complete gene sequence was not obtained for all the cultivars, yet there were a total of 38 variable positions in exonic (coding region) sequence, from a total gene length of 453 nucleotides. Matrix of SNP shows these 37 positions with individual sequence at positions given for each of the 14 cultivars (sequence of two cultivars was not obtained) included in this analysis. It demonstrated a considerable diversity for this gene with only three cultivars i.e. TJ-83, Marvi and TD-1 being similar to the consensus sequence. All other cultivars showed a unique combination of SNPs. In order to prove a functional link between these polymorphisms and drought tolerance in wheat, it would be necessary to conduct a more detailed study involving directed mutation of this gene and DHN gene expression.

Endosomal proteolysis regulates calcitonin gene-related peptide responses in mesenteric arteries

Background and Purpose: Calcitonin gene‐related peptide (CGRP) is a potent vasodilator, implicated in the pathogenesis of migraine. CGRP activates a receptor complex comprising, calcitonin receptor‐like receptor (CLR) and receptor activity‐modifying protein 1 (RAMP1). In vitro studies indicate recycling of CLR•RAMP1 is regulated by degradation of CGRP in early endosomes by endothelin‐converting enzyme‐1 (ECE‐1). However, it is not known if ECE‐1 regulates the resensitization of CGRP‐induced responses in functional arterial tissue. Experimental Approach: CLR, ECE‐1a‐d and RAMP1 expression in rat mesenteric artery smooth muscle cells (RMA‐SMCs) and mesenteric arteries was analyzed by RT‐PCR and by immunofluorescence and confocal microscopy. CGRP‐induced signaling in cells was examined by measuring cAMP production and ERK activation. CGRP‐induced relaxation of arteries was measured by isometric wire myography. ECE‐1 was inhibited using the specific inhibitor, SM‐19712. Key Results: RMA‐SMCs and arteries contained mRNA for CLR, ECE‐1a‐d and RAMP1. ECE‐1 was present in early endosomes of RMA‐SMCs and in the smooth muscle layer of arteries. CGRP induced endothelium‐independent relaxation of arteries. ECE‐1 inhibition had no effect on initial CGRP‐induced responses but reduced cAMP generation in RMA‐SMCs and vasodilation in mesenteric arteries responses to subsequent CGRP challenges. Conclusions and Implications: ECE‐1 regulates the resensitization of responses to CGRP in RMA‐SMCs and mesenteric arteries. CGRP‐induced relaxation does not involve endothelium‐derived pathways. This is the first report of ECE‐1 regulating CGRP responses in SMCs and arteries. ECE‐1 inhibitors may attenuate an important vasodilatory pathway, implicated in primary headaches and may represent a new therapeutic approach for the treatment of migraine.

Adeno-associated virus-8-Mediated intravenous transfer of myostatin propeptide leads to systemic functional improvements of slow but not fast muscle

Myostatin is a member of the transformating growth factor-_ (TGF-_) superfamily of proteins and is produced almost exclusively in skeletal muscle tissue, where it is secreted and circulates as a serum protein. Myostatin acts as a negative regulator of muscle mass through the canonical SMAD2/3/4 signaling pathway. Naturally occurring myostatin mutants exhibit a ‘double muscling’ phenotype in which muscle mass is dramatically increased as a result of both hypertrophy and hyperplasia. Myostatin is naturally inhibited by its own propeptide; therefore, we assessed the impact of adeno associated virus-8 (AAV8) myostatin propeptide vectors when systemically introduced in MF-1 mice. We noted a significant systemic increase in muscle mass in both slow and fast muscle phenotypes, with no evidence of hyperplasia; however, the nuclei-to- cytoplasm ratio in all myofiber types was significantly reduced. An increase in muscle mass in slow (soleus) muscle led to an increase in force output; however, an increase in fast (extensor digitorum longus [EDL]) muscle mass did not increase force output. These results suggest that the use of gene therapeutic regimens of myostatin inhibition for age-related or disease-related muscle loss may have muscle-specific effects.

A functional link between bone morphogenetic proteins and insulin-like peptide 3 signaling in modulating ovarian androgen production

Bone morphogenetic proteins (BMP) are firmly implicated as intra-ovarian regulators of follicle development and steroidogenesis. Here we report a microarray analysis showing that treatment of cultured bovine theca cells (TC) with BMP6 significantly (>2-fold; P<0.01) up- or down-regulated expression of 445 genes. Insulin-like peptide 3 (INSL3) was the most heavily down-regulated gene (-43-fold) with CYP17A1 and other key transcripts involved in TC steroidogenesis including LHCGR, INHA, STAR, CYP11A1 and HSD3B1 also down-regulated. BMP6 also reduced expression of NR5A1 encoding steroidogenic factor-1 known to target the promoter regions of the aforementioned genes. Real-time PCR confirmed these findings and also revealed a marked reduction in expression of INSL3 receptor (RXFP2). Secretion of INSL3 protein and androstenedione were also suppressed suggesting a functional link between BMP and INSL3 pathways in controlling androgen synthesis. RNAi-mediated knockdown of INSL3 reduced INSL3 mRNA and secreted protein level (75 and 94%, respectively) and elicited a 77% reduction in CYP17A1 mRNA level and 83% reduction in androstenedione secretion. Knockdown of RXFP2 also reduced CYP17A1 mRNA level (81%) and androstenedione secretion (88%). Conversely, treatment with exogenous (human) INSL3 increased androstenedione secretion ~2-fold. The CYP17 inhibitor abiraterone abolished androgen secretion and reduced expression of both INSL3 and RXFP2. Collectively, these findings indicate a positive autoregulatory role for INSL3 signaling in maintaining thecal androgen production, and visa versa. Moreover, BMP6-induced suppression of thecal androgen synthesis may be mediated, at least in part, by reduced INSL3-RXFP2 signaling.