On the Use of Gene Ontology Annotations to Assess Functional Similarity among Orthologs and Paralogs: A Short Report

A recent paper (Nehrt et al., PLoS Comput. Biol. 7:e1002073, 2011) has proposed a metric for the “functional similarity” between two genes that uses only the Gene Ontology (GO) annotations directly derived from published experimental results. Applying this metric, the authors concluded that paralogous genes within the mouse genome or the human genome are more functionally similar on average than orthologous genes between these genomes, an unexpected result with broad implications if true. We suggest, based on both theoretical and empirical considerations, that this proposed metric should not be interpreted as a functional similarity, and therefore cannot be used to support any conclusions about the “ortholog conjecture” (or, more properly, the “ortholog functional conservation hypothesis”). First, we reexamine the case studies presented by Nehrt et al. as examples of orthologs with divergent functions, and come to a very different conclusion: they actually exemplify how GO annotations for orthologous genes provide complementary information about conserved biological functions. We then show that there is a global ascertainment bias in the experiment-based GO annotations for human and mouse genes: particular types of experiments tend to be performed in different model organisms. We conclude that the reported statistical differences in annotations between pairs of orthologous genes do not reflect differences in biological function, but rather complementarity in experimental approaches. Our results underscore two general considerations for researchers proposing novel types of analysis based on the GO: 1) that GO annotations are often incomplete, potentially in a biased manner, and subject to an “open world assumption” (absence of an annotation does not imply absence of a function), and 2) that conclusions drawn from a novel, large-scale GO analysis should whenever possible be supported by careful, in-depth examination of examples, to help ensure the conclusions have a justifiable biological basis

The nuclear receptors of Biomphalaria glabrata and Lottia gigantea: Implications for developing new model organisms

© 2015 Kaur et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are creditedNuclear receptors (NRs) are transcription regulators involved in an array of diverse physiological functions including key roles in endocrine and metabolic function. The aim of this study was to identify nuclear receptors in the fully sequenced genome of the gastropod snail, Biomphalaria glabrata, intermediate host for Schistosoma mansoni and compare these to known vertebrate NRs, with a view to assessing the snail's potential as a invertebrate model organism for endocrine function, both as a prospective new test organism and to elucidate the fundamental genetic and mechanistic causes of disease. For comparative purposes, the genome of a second gastropod, the owl limpet, Lottia gigantea was also investigated for nuclear receptors. Thirty-nine and thirty-three putative NRs were identified from the B. glabrata and L. gigantea genomes respectively, based on the presence of a conserved DNA-binding domain and/or ligand-binding domain. Nuclear receptor transcript expression was confirmed and sequences were subjected to a comparative phylogenetic analysis, which demonstrated that these molluscs have representatives of all the major NR subfamilies (1-6). Many of the identified NRs are conserved between vertebrates and invertebrates, however differences exist, most notably, the absence of receptors of Group 3C, which includes some of the vertebrate endocrine hormone targets. The mollusc genomes also contain NR homologues that are present in insects and nematodes but not in vertebrates, such as Group 1J (HR48/DAF12/HR96). The identification of many shared receptors between humans and molluscs indicates the potential for molluscs as model organisms; however the absence of several steroid hormone receptors indicates snail endocrine systems are fundamentally different.The National Centre for the Replacement, Refinement and Reduction of Animals in Research, Grant Ref:G0900802 to CSJ, LRN, SJ & EJR [www.nc3rs.org.uk]

Mechanisms for the Evolution of a Derived Function in the Ancestral Glucocorticoid Receptor

Understanding the genetic, structural, and biophysical mechanisms that caused protein functions to evolve is a central goal of molecular evolutionary studies. Ancestral sequence reconstruction (ASR) offers an experimental approach to these questions. Here we use ASR to shed light on the earliest functions and evolution of the glucocorticoid receptor (GR), a steroid-activated transcription factor that plays a key role in the regulation of vertebrate physiology. Prior work showed that GR and its paralog, the mineralocorticoid receptor (MR), duplicated from a common ancestor roughly 450 million years ago; the ancestral functions were largely conserved in the MR lineage, but the functions of GRs—reduced sensitivity to all hormones and increased selectivity for glucocorticoids—are derived. Although the mechanisms for the evolution of glucocorticoid specificity have been identified, how reduced sensitivity evolved has not yet been studied. Here we report on the reconstruction of the deepest ancestor in the GR lineage (AncGR1) and demonstrate that GR's reduced sensitivity evolved before the acquisition of restricted hormone specificity, shortly after the GR–MR split. Using site-directed mutagenesis, X-ray crystallography, and computational analyses of protein stability to recapitulate and determine the effects of historical mutations, we show that AncGR1's reduced ligand sensitivity evolved primarily due to three key substitutions. Two large-effect mutations weakened hydrogen bonds and van der Waals interactions within the ancestral protein, reducing its stability. The degenerative effect of these two mutations is extremely strong, but a third permissive substitution, which has no apparent effect on function in the ancestral background and is likely to have occurred first, buffered the effects of the destabilizing mutations. Taken together, our results highlight the potentially creative role of substitutions that partially degrade protein structure and function and reinforce the importance of permissive mutations in protein evolution

Protein Evolution by Molecular Tinkering: Diversification of the Nuclear Receptor Superfamily from a Ligand-Dependent Ancestor

Phylogenetic reconstruction of the structure and function of the ancestor of the nuclear receptor protein family reveals how functional diversity evolves by subtle tinkering with an ancestral template

Evolution of a Signaling Nexus Constrained by Protein Interfaces and Conformational States

Heterotrimeric G proteins act as the physical nexus between numerous receptors that respond to extracellular signals and proteins that drive the cytoplasmic response. The Gα subunit of the G protein, in particular, is highly constrained due to its many interactions with proteins that control or react to its conformational state. Various organisms contain differing sets of Gα-interacting proteins, clearly indicating that shifts in sequence and associated Gα functionality were acquired over time. These numerous interactions constrained much of Gα evolution; yet Gα has diversified, through poorly understood processes, into several functionally specialized classes, each with a unique set of interacting proteins. Applying a synthetic sequence-based approach to mammalian Gα subunits, we established a set of seventy-five evolutionarily important class-distinctive residues, sites where a single Gα class is differentiated from the three other classes. We tested the hypothesis that shifts at these sites are important for class-specific functionality. Importantly, we mapped known and well-studied class-specific functionalities from all four mammalian classes to sixteen of our class-distinctive sites, validating the hypothesis. Our results show how unique functionality can evolve through the recruitment of residues that were ancestrally functional. We also studied acquisition of functionalities by following these evolutionarily important sites in non-mammalian organisms. Our results suggest that many class-distinctive sites were established early on in eukaryotic diversification and were critical for the establishment of new Gα classes, whereas others arose in punctuated bursts throughout metazoan evolution. These Gα class-distinctive residues are rational targets for future structural and functional studies

MPI-PHYLIP: Parallelizing Computationally Intensive Phylogenetic Analysis Routines for the Analysis of Large Protein Families

Background: Phylogenetic study of protein sequences provides unique and valuable insights into the molecular and genetic basis of important medical and epidemiological problems as well as insights about the origins and development of physiological features in present day organisms. Consensus phylogenies based on the bootstrap and other resampling methods play a crucial part in analyzing the robustness of the trees produced for these analyses. Methodology: Our focus was to increase the number of bootstrap replications that can be performed on large protein datasets using the maximum parsimony, distance matrix, and maximum likelihood methods. We have modified the PHYLIP package using MPI to enable large-scale phylogenetic study of protein sequences, using a statistically robust number of bootstrapped datasets, to be performed in a moderate amount of time. This paper discusses the methodology used to parallelize the PHYLIP programs and reports the performance of the parallel PHYLIP programs that are relevant to the study of protein evolution on several protein datasets. Conclusions: Calculations that currently take a few days on a state of the art desktop workstation are reduced to calculations that can be performed over lunchtime on a modern parallel computer. Of the three protein methods tested, the maximum likelihood method scales the best, followed by the distance method, and then the maximum parsimony method. However, the maximum likelihood method requires significant memory resources, which limits its application to mor

Peatland Initiation, Carbon Accumulation, and 2 ka Depth in the James Bay Lowland and Adjacent Regions

Copyright © 2014 University of Colorado at Boulder, Institute of Arctic and Alpine ResearchPeatlands surrounding Hudson and James Bays form the second largest peatland complex in the world and contain major stores of soil carbon (C). This study utilized a transect of eight ombrotrophic peat cores from remote regions of central and northern Ontario to quantify the magnitude and rate of C accumulation since peatland initiation and for the past 2000 calendar years before present (2 ka). These new data were supplemented by 17 millennially resolved chronologies from a literature review covering the Boreal Shield, Hudson Plains, and Taiga Shield bordering Hudson and James Bays. Peatlands initiated in central and northern Ontario by 7.8 ka following deglaciation and isostatic emergence of northern areas to above sea level. Total C accumulated since inception averaged 109.7 ± (std. dev.) 36.2 kg C m–2. Approximately 40% of total soil C has accumulated since 2 ka at an average apparent rate of 20.2 ± 6.9 g C m–2 yr–1. The 2 ka depths correlate significantly and positively with modern gridded climate estimates for mean annual precipitation, mean annual air temperature, growing degree-days > 0 °C, and photosynthetically active radiation integrated over days > 0 °C. There are significantly shallower depths in permafrost peatlands. Vertical peat accumulation was likely constrained by temperature, growing season length, and photosynthetically active radiation over the last 2 ka in the Hudson Bay Lowlands and surrounding regions.US National Science Foundatio

Methane Fluxes in Tidal Marshes of the Conterminous United States

Methane (CH4) is a potent greenhouse gas (GHG) with atmospheric concentrations that have nearly tripled since pre-industrial times. Wetlands account for a large share of global CH4 emissions, yet the magnitude and factors controlling CH4 fluxes in tidal wetlands remain uncertain. We synthesized CH4 flux data from 100 chamber and 9 eddy covariance (EC) sites across tidal marshes in the conterminous United States to assess controlling factors and improve predictions of CH4 emissions. This effort included creating an open-source database of chamber-based GHG fluxes (https://doi.org/10.25573/serc.14227085). Annual fluxes across chamber and EC sites averaged 26 ± 53 g CH4 m−2 year−1, with a median of 3.9 g CH4 m−2 year−1, and only 25% of sites exceeding 18 g CH4 m−2 year−1. The highest fluxes were observed at fresh-oligohaline sites with daily maximum temperature normals (MATmax) above 25.6°C. These were followed by frequently inundated low and mid-fresh-oligohaline marshes with MATmax ≤25.6°C, and mesohaline sites with MATmax \u3e19°C. Quantile regressions of paired chamber CH4 flux and porewater biogeochemistry revealed that the 90th percentile of fluxes fell below 5 ± 3 nmol m−2 s−1 at sulfate concentrations \u3e4.7 ± 0.6 mM, porewater salinity \u3e21 ± 2 psu, or surface water salinity \u3e15 ± 3 psu. Across sites, salinity was the dominant predictor of annual CH4 fluxes, while within sites, temperature, gross primary productivity (GPP), and tidal height controlled variability at diel and seasonal scales. At the diel scale, GPP preceded temperature in importance for predicting CH4 flux changes, while the opposite was observed at the seasonal scale. Water levels influenced the timing and pathway of diel CH4 fluxes, with pulsed releases of stored CH4 at low to rising tide. This study provides data and methods to improve tidal marsh CH4 emission estimates, support blue carbon assessments, and refine national and global GHG inventories

Evolution of ligand specificity in vertebrate corticosteroid receptors

<p>Abstract</p> <p>Background</p> <p>Corticosteroid receptors include mineralocorticoid (MR) and glucocorticoid (GR) receptors. Teleost fishes have a single MR and duplicate GRs that show variable sensitivities to mineralocorticoids and glucocorticoids. How these receptors compare functionally to tetrapod MR and GR, and the evolutionary significance of maintaining two GRs, remains unclear.</p> <p>Results</p> <p>We used up to seven steroids (including aldosterone, cortisol and 11-deoxycorticosterone [DOC]) to compare the ligand specificity of the ligand binding domains of corticosteroid receptors between a mammal (<it>Mus musculus</it>) and the midshipman fish (<it>Porichthys notatus</it>), a teleost model for steroid regulation of neural and behavioral plasticity. Variation in mineralocorticoid sensitivity was considered in a broader phylogenetic context by examining the aldosterone sensitivity of MR and GRs from the distantly related daffodil cichlid (<it>Neolamprologus pulcher</it>), another teleost model for neurobehavioral plasticity. Both teleost species had a single MR and duplicate GRs. All MRs were sensitive to DOC, consistent with the hypothesis that DOC was the initial ligand of the ancestral MR. Variation in GR steroid-specificity corresponds to nine identified amino acid residue substitutions rather than phylogenetic relationships based on receptor sequences.</p> <p>Conclusion</p> <p>The mineralocorticoid sensitivity of duplicate GRs in teleosts is highly labile in the context of their evolutionary phylogeny, a property that likely led to neo-functionalization and maintenance of two GRs.</p

The Origins of Novel Protein Interactions during Animal Opsin Evolution

Background. Biologists are gaining an increased understanding of the genetic bases of phenotypic change during evolution. Nevertheless, the origins of phenotypes mediated by novel protein-protein interactions remain largely undocumented. Methodology/Principle Findings. Here we analyze the evolution of opsin visual pigment proteins from the genomes of early branching animals, including a new class of opsins from Cnidaria. We combine these data with existing knowledge of the molecular basis of opsin function in a rigorous phylogenetic framework. We identify adaptive amino acid substitutions in duplicated opsin genes that correlate with a diversification of physiological pathways mediated by different protein-protein interactions. Conclusions/Significance. This study documents how gene duplication events early in the history of animals followed by adaptive structural mutations increased organismal complexity by adding novel protein-protein interactions that underlie different physiological pathways. These pathways are central to vision and other photo-reactive phenotypes in most extant animals. Similar evolutionary processes may have been a work in generating other metazoan sensory systems and other physiological processes mediated by signal transduction