Selective deployment of transcription factor paralogs with submaximal strength facilitates gene regulation in the immune system.

In multicellular organisms, duplicated genes can diverge through tissue-specific gene expression patterns, as exemplified by highly regulated expression of RUNX transcription factor paralogs with apparent functional redundancy. Here we asked what cell-type-specific biologies might be supported by the selective expression of RUNX paralogs during Langerhans cell and inducible regulatory T cell differentiation. We uncovered functional nonequivalence between RUNX paralogs. Selective expression of native paralogs allowed integration of transcription factor activity with extrinsic signals, while non-native paralogs enforced differentiation even in the absence of exogenous inducers. DNA binding affinity was controlled by divergent amino acids within the otherwise highly conserved RUNT domain and evolutionary reconstruction suggested convergence of RUNT domain residues toward submaximal strength. Hence, the selective expression of gene duplicates in specialized cell types can synergize with the acquisition of functional differences to enable appropriate gene expression, lineage choice and differentiation in the mammalian immune system.

Conserved phosphorylation hotspots in eukaryotic protein domain families.

Protein phosphorylation is the best characterized post-translational modification that regulates almost all cellular processes through diverse mechanisms such as changing protein conformations, interactions, and localization. While the inventory for phosphorylation sites across different species has rapidly expanded, their functional role remains poorly investigated. Here, we combine 537,321 phosphosites from 40 eukaryotic species to identify highly conserved phosphorylation hotspot regions within domain families. Mapping these regions onto structural data reveals that they are often found at interfaces, near catalytic residues and tend to harbor functionally important phosphosites. Notably, functional studies of a phospho-deficient mutant in the C-terminal hotspot region within the ribosomal S11 domain in the yeast ribosomal protein uS11 shows impaired growth and defective cytoplasmic 20S pre-rRNA processing at 16 °C and 20 °C. Altogether, our study identifies phosphorylation hotspots for 162 protein domains suggestive of an ancient role for the control of diverse eukaryotic domain families.

iProteinDB: An Integrative Database of Drosophila Post-translational Modifications.

Post-translational modification (PTM) serves as a regulatory mechanism for protein function, influencing their stability, interactions, activity and localization, and is critical in many signaling pathways. The best characterized PTM is phosphorylation, whereby a phosphate is added to an acceptor residue, most commonly serine, threonine and tyrosine in metazoans. As proteins are often phosphorylated at multiple sites, identifying those sites that are important for function is a challenging problem. Considering that any given phosphorylation site might be non-functional, prioritizing evolutionarily conserved phosphosites provides a general strategy to identify the putative functional sites. To facilitate the identification of conserved phosphosites, we generated a large-scale phosphoproteomics dataset from Drosophila embryos collected from six closely-related species. We built iProteinDB (https://www.flyrnai.org/tools/iproteindb/), a resource integrating these data with other high-throughput PTM datasets, including vertebrates, and manually curated information for Drosophila At iProteinDB, scientists can view the PTM landscape for any Drosophila protein and identify predicted functional phosphosites based on a comparative analysis of data from closely-related Drosophila species. Further, iProteinDB enables comparison of PTM data from Drosophila to that of orthologous proteins from other model organisms, including human, mouse, rat, Xenopus tropicalis, Danio rerio, and Caenorhabditis elegans.

Absence of Selection for Quantum Coherence in the Fenna-Matthews-Olson Complex: A Combined Evolutionary and Excitonic Study.

We present a study on the evolution of the Fenna-Matthews-Olson bacterial photosynthetic pigment-protein complex. This protein complex functions as an antenna. It transports absorbed photons-excitons-to a reaction center where photosynthetic reactions initiate. The efficiency of exciton transport is therefore fundamental for the photosynthetic bacterium's survival. We have reconstructed an ancestor of the complex to establish whether coherence in the exciton transport was selected for or optimized over time. We have also investigated the role of optimizing free energy variation upon folding in evolution. We studied whether mutations which connect the ancestor to current day species were stabilizing or destabilizing from a thermodynamic viewpoint. From this study, we established that most of these mutations were thermodynamically neutral. Furthermore, we did not see a large change in exciton transport efficiency or coherence, and thus our results predict that exciton coherence was not specifically selected for.

Evolution of protein phosphorylation across 18 fungal species.

Living organisms have evolved protein phosphorylation, a rapid and versatile mechanism that drives signaling and regulates protein function. We report the phosphoproteomes of 18 fungal species and a phylogenetic-based approach to study phosphosite evolution. We observe rapid divergence, with only a small fraction of phosphosites conserved over hundreds of millions of years. Relative to recently acquired phosphosites, ancient sites are enriched at protein interfaces and are more likely to be functionally important, as we show for sites on H2A1 and eIF4E. We also observe a change in phosphorylation motif frequencies and kinase activities that coincides with the whole-genome duplication event. Our results provide an evolutionary history for phosphosites and suggest that rapid evolution of phosphorylation can contribute strongly to phenotypic diversity.

Evolutionary diversification of retinoic acid receptor ligand-binding pocket structure by molecular tinkering.

Whole genome duplications (WGDs) have been classically associated with the origin of evolutionary novelties and the so-called duplication-degeneration-complementation model describes the possible fates of genes after duplication. However, how sequence divergence effectively allows functional changes between gene duplicates is still unclear. In the vertebrate lineage, two rounds of WGDs took place, giving rise to paralogous gene copies observed for many gene families. For the retinoic acid receptors (RARs), for example, which are members of the nuclear hormone receptor (NR) superfamily, a unique ancestral gene has been duplicated resulting in three vertebrate paralogues: RARα, RARß and RARγ. It has previously been shown that this single ancestral RAR was neofunctionalized to give rise to a larger substrate specificity range in the RARs of extant jawed vertebrates (also called gnathostomes). To understand RAR diversification, the members of the cyclostomes (lamprey and hagfish), jawless vertebrates representing the extant sister group of gnathostomes, provide an intermediate situation and thus allow the characterization of the evolutionary steps that shaped RAR ligand-binding properties following the WGDs. In this study, we assessed the ligand-binding specificity of cyclostome RARs and found that their ligand-binding pockets resemble those of gnathostome RARα and RARß. In contrast, none of the cyclostome receptors studied showed any RARγ-like specificity. Together, our results suggest that cyclostome RARs cover only a portion of the specificity repertoire of the ancestral gnathostome RARs and indicate that the establishment of ligand-binding specificity was a stepwise event. This iterative process thus provides a rare example for the diversification of receptor-ligand interactions of NRs following WGDs.

Correction to 'Evolutionary diversification of retinoic acid receptor ligand-binding pocket structure by molecular tinkering'.

[This corrects the article DOI: 10.1098/rsos.150484.].

CATH: comprehensive structural and functional annotations for genome sequences.

The latest version of the CATH-Gene3D protein structure classification database (4.0, http://www.cathdb.info) provides annotations for over 235,000 protein domain structures and includes 25 million domain predictions. This article provides an update on the major developments in the 2 years since the last publication in this journal including: significant improvements to the predictive power of our functional families (FunFams); the release of our 'current' putative domain assignments (CATH-B); a new, strictly non-redundant data set of CATH domains suitable for homology benchmarking experiments (CATH-40) and a number of improvements to the web pages.

Origin of the response to adrenal and sex steroids: Roles of promiscuity and co-evolution of enzymes and steroid receptors.

Many responses to adrenal and sex steroids are mediated by receptors that belong to the nuclear receptor family of transcription factors. We investigated the co-evolution of these vertebrate steroid receptors and the enzymes that synthesize adrenal and sex steroids through data mining of genomes from cephalochordates [amphioxus], cyclostomes [lampreys, hagfish], chondrichthyes [sharks, rays, skates], actinopterygii [ray-finned fish], sarcopterygii [coelacanths, lungfishes and terrestrial vertebrates]. An ancestor of the estrogen receptor and 3-ketosteroid receptors evolved in amphioxus. A corticoid receptor and a progesterone receptor evolved in cyclostomes, and an androgen receptor evolved in gnathostomes. Amphioxus contains CYP11, CYP17, CYP19, 3ß/Δ5-4-HSD and 17ß-HSD14, which suffice for the synthesis of estradiol and Δ5-androstenediol. Amphioxus also contains CYP27, which catalyzes the synthesis of 27-hydroxy-cholesterol, another estrogen. Lamprey contains, in addition, CYP21, which catalyzes the synthesis of 11-deoxycortisol. Chondrichthyes contain, in addition, CYP11A, CYP11C, CYP17A1, CYP17A2. Coelacanth also contains CYP11C1, the current descendent from a common ancestor with modern land vertebrate CYP11B genes, which catalyze the synthesis of cortisol, corticosterone and aldosterone. Interestingly, CYP11B2, aldosterone synthase, evolved from separate gene duplications in at least old world monkeys and two suborders of rodents. Sciurognathi (including mice and rats) and Hystricomorpha (including guinea pigs). Thus, steroid receptors and steroidogenic enzymes co-evolved at key transitions in the evolution of vertebrates. Together, this suite of receptors and enzymes through their roles in transcriptional regulation of reproduction, development, homeostasis and the response to stress contributed to the evolutionary diversification of vertebrates. This article is part of a Special Issue entitled 'Steroid/Sterol signaling'.

Epithelial sodium transport and its control by aldosterone: the story of our internal environment revisited.

Transcription and translation require a high concentration of potassium across the entire tree of life. The conservation of a high intracellular potassium was an absolute requirement for the evolution of life on Earth. This was achieved by the interplay of P- and V-ATPases that can set up electrochemical gradients across the cell membrane, an energetically costly process requiring the synthesis of ATP by F-ATPases. In animals, the control of an extracellular compartment was achieved by the emergence of multicellular organisms able to produce tight epithelial barriers creating a stable extracellular milieu. Finally, the adaptation to a terrestrian environment was achieved by the evolution of distinct regulatory pathways allowing salt and water conservation. In this review we emphasize the critical and dual role of Na(+)-K(+)-ATPase in the control of the ionic composition of the extracellular fluid and the renin-angiotensin-aldosterone system (RAAS) in salt and water conservation in vertebrates. The action of aldosterone on transepithelial sodium transport by activation of the epithelial sodium channel (ENaC) at the apical membrane and that of Na(+)-K(+)-ATPase at the basolateral membrane may have evolved in lungfish before the emergence of tetrapods. Finally, we discuss the implication of RAAS in the origin of the present pandemia of hypertension and its associated cardiovascular diseases.

Understanding the functional difference between growth arrest-specific protein 6 and protein S: an evolutionary approach.

Although protein S (PROS1) and growth arrest-specific protein 6 (GAS6) proteins are homologous with a high degree of structural similarity, they are functionally different. The objectives of this study were to identify the evolutionary origins from which these functional differences arose. Bioinformatics methods were used to estimate the evolutionary divergence time and to detect the amino acid residues under functional divergence between GAS6 and PROS1. The properties of these residues were analysed in the light of their three-dimensional structures, such as their stability effects, the identification of electrostatic patches and the identification potential protein-protein interaction. The divergence between GAS6 and PROS1 probably occurred during the whole-genome duplications in vertebrates. A total of 78 amino acid sites were identified to be under functional divergence. One of these sites, Asn463, is involved in N-glycosylation in GAS6, but is mutated in PROS1, preventing this post-translational modification. Sites experiencing functional divergence tend to express a greater diversity of stabilizing/destabilizing effects than sites that do not experience such functional divergence. Three electrostatic patches in the LG1/LG2 domains were found to differ between GAS6 and PROS1. Finally, a surface responsible for protein-protein interactions was identified. These results may help researchers to analyse disease-causing mutations in the light of evolutionary and structural constraints, and link genetic pathology to clinical phenotypes.

Positive selection during the evolution of the blood coagulation factors in the context of their disease-causing mutations.

Blood coagulation occurs through a cascade of enzymes and cofactors that produces a fibrin clot, while otherwise maintaining hemostasis. The 11 human coagulation factors (FG, FII-FXIII) have been identified across all vertebrates, suggesting that they emerged with the first vertebrates around 500 Ma. Human FVIII, FIX, and FXI are associated with thousands of disease-causing mutations. Here, we evaluated the strength of selective pressures on the 14 genes coding for the 11 factors during vertebrate evolution, and compared these with human mutations in FVIII, FIX, and FXI. Positive selection was identified for fibrinogen (FG), FIII, FVIII, FIX, and FX in the mammalian Primates and Laurasiatheria and the Sauropsida (reptiles and birds). This showed that the coagulation system in vertebrates was under strong selective pressures, perhaps to adapt against blood-invading pathogens. The comparison of these results with disease-causing mutations reported in FVIII, FIX, and FXI showed that the number of disease-causing mutations, and the probability of positive selection were inversely related to each other. It was concluded that when a site was under positive selection, it was less likely to be associated with disease-causing mutations. In contrast, sites under negative selection were more likely to be associated with disease-causing mutations and be destabilizing. A residue-by-residue comparison of the FVIII, FIX, and FXI sequence alignments confirmed this. This improved understanding of evolutionary changes in FVIII, FIX, and FXI provided greater insight into disease-causing mutations, and better assessments of the codon sites that may be mutated in applications of gene therapy.

Stability-activity tradeoffs constrain the adaptive evolution of RubisCO.

A well-known case of evolutionary adaptation is that of ribulose-1,5-bisphosphate carboxylase (RubisCO), the enzyme responsible for fixation of CO2 during photosynthesis. Although the majority of plants use the ancestral C3 photosynthetic pathway, many flowering plants have evolved a derived pathway named C4 photosynthesis. The latter concentrates CO2, and C4 RubisCOs consequently have lower specificity for, and faster turnover of, CO2. The C4 forms result from convergent evolution in multiple clades, with substitutions at a small number of sites under positive selection. To understand the physical constraints on these evolutionary changes, we reconstructed in silico ancestral sequences and 3D structures of RubisCO from a large group of related C3 and C4 species. We were able to precisely track their past evolutionary trajectories, identify mutations on each branch of the phylogeny, and evaluate their stability effect. We show that RubisCO evolution has been constrained by stability-activity tradeoffs similar in character to those previously identified in laboratory-based experiments. The C4 properties require a subset of several ancestral destabilizing mutations, which from their location in the structure are inferred to mainly be involved in enhancing conformational flexibility of the open-closed transition in the catalytic cycle. These mutations are near, but not in, the active site or at intersubunit interfaces. The C3 to C4 transition is preceded by a sustained period in which stability of the enzyme is increased, creating the capacity to accept the functionally necessary destabilizing mutations, and is immediately followed by compensatory mutations that restore global stability.

Selectome update: quality control and computational improvements to a database of positive selection.

Selectome (http://selectome.unil.ch/) is a database of positive selection, based on a branch-site likelihood test. This model estimates the number of nonsynonymous substitutions (dN) and synonymous substitutions (dS) to evaluate the variation in selective pressure (dN/dS ratio) over branches and over sites. Since the original release of Selectome, we have benchmarked and implemented a thorough quality control procedure on multiple sequence alignments, aiming to provide minimum false-positive results. We have also improved the computational efficiency of the branch-site test implementation, allowing larger data sets and more frequent updates. Release 6 of Selectome includes all gene trees from Ensembl for Primates and Glires, as well as a large set of vertebrate gene trees. A total of 6810 gene trees have some evidence of positive selection. Finally, the web interface has been improved to be more responsive and to facilitate searches and browsing.

Gene3D: Multi-domain annotations for protein sequence and comparative genome analysis.

Gene3D (http://gene3d.biochem.ucl.ac.uk) is a database of protein domain structure annotations for protein sequences. Domains are predicted using a library of profile HMMs from 2738 CATH superfamilies. Gene3D assigns domain annotations to Ensembl and UniProt sequence sets including >6000 cellular genomes and >20 million unique protein sequences. This represents an increase of 45% in the number of protein sequences since our last publication. Thanks to improvements in the underlying data and pipeline, we see large increases in the domain coverage of sequences. We have expanded this coverage by integrating Pfam and SUPERFAMILY domain annotations, and we now resolve domain overlaps to provide highly comprehensive composite multi-domain architectures. To make these data more accessible for comparative genome analyses, we have developed novel search algorithms for searching genomes to identify related multi-domain architectures. In addition to providing domain family annotations, we have now developed a pipeline for 3D homology modelling of domains in Gene3D. This has been applied to the human genome and will be rolled out to other major organisms over the next year.

Vitellogenin underwent subfunctionalization to acquire caste and behavioral specific expression in the harvester ant Pogonomyrmex barbatus.

The reproductive ground plan hypothesis (RGPH) proposes that the physiological pathways regulating reproduction were co-opted to regulate worker division of labor. Support for this hypothesis in honeybees is provided by studies demonstrating that the reproductive potential of workers, assessed by the levels of vitellogenin (Vg), is linked to task performance. Interestingly, contrary to honeybees that have a single Vg ortholog and potentially fertile nurses, the genome of the harvester ant Pogonomyrmex barbatus harbors two Vg genes (Pb_Vg1 and Pb_Vg2) and nurses produce infertile trophic eggs. P. barbatus, thus, provides a unique model to investigate whether Vg duplication in ants was followed by subfunctionalization to acquire reproductive and non-reproductive functions and whether Vg reproductive function was co-opted to regulate behavior in sterile workers. To investigate these questions, we compared the expression patterns of P. barbatus Vg genes and analyzed the phylogenetic relationships and molecular evolution of Vg genes in ants. qRT-PCRs revealed that Pb_Vg1 is more highly expressed in queens compared to workers and in nurses compared to foragers. By contrast, the level of expression of Pb_Vg2 was higher in foragers than in nurses and queens. Phylogenetic analyses show that a first duplication of the ancestral Vg gene occurred after the divergence between the poneroid and formicoid clades and subsequent duplications occurred in the lineages leading to Solenopsis invicta, Linepithema humile and Acromyrmex echinatior. The initial duplication resulted in two Vg gene subfamilies preferentially expressed in queens and nurses (subfamily A) or in foraging workers (subfamily B). Finally, molecular evolution analyses show that the subfamily A experienced positive selection, while the subfamily B showed overall relaxation of purifying selection. Our results suggest that in P. barbatus the Vg gene underwent subfunctionalization after duplication to acquire caste- and behavior- specific expression associated with reproductive and non-reproductive functions, supporting the validity of the RGPH in ants.

Residue mutations and their impact on protein structure and function: detecting beneficial and pathogenic changes.

The present review focuses on the evolution of proteins and the impact of amino acid mutations on function from a structural perspective. Proteins evolve under the law of natural selection and undergo alternating periods of conservative evolution and of relatively rapid change. The likelihood of mutations being fixed in the genome depends on various factors, such as the fitness of the phenotype or the position of the residues in the three-dimensional structure. For example, co-evolution of residues located close together in three-dimensional space can occur to preserve global stability. Whereas point mutations can fine-tune the protein function, residue insertions and deletions ('decorations' at the structural level) can sometimes modify functional sites and protein interactions more dramatically. We discuss recent developments and tools to identify such episodic mutations, and examine their applications in medical research. Such tools have been tested on simulated data and applied to real data such as viruses or animal sequences. Traditionally, there has been little if any cross-talk between the fields of protein biophysics, protein structure-function and molecular evolution. However, the last several years have seen some exciting developments in combining these approaches to obtain an in-depth understanding of how proteins evolve. For example, a better understanding of how structural constraints affect protein evolution will greatly help us to optimize our models of sequence evolution. The present review explores this new synthesis of perspectives.

New functional families (FunFams) in CATH to improve the mapping of conserved functional sites to 3D structures.

CATH version 3.5 (Class, Architecture, Topology, Homology, available at http://www.cathdb.info/) contains 173 536 domains, 2626 homologous superfamilies and 1313 fold groups. When focusing on structural genomics (SG) structures, we observe that the number of new folds for CATH v3.5 is slightly less than for previous releases, and this observation suggests that we may now know the majority of folds that are easily accessible to structure determination. We have improved the accuracy of our functional family (FunFams) sub-classification method and the CATH sequence domain search facility has been extended to provide FunFam annotations for each domain. The CATH website has been redesigned. We have improved the display of functional data and of conserved sequence features associated with FunFams within each CATH superfamily.

Resolving the ortholog conjecture: orthologs tend to be weakly, but significantly, more similar in function than paralogs.

The function of most proteins is not determined experimentally, but is extrapolated from homologs. According to the "ortholog conjecture", or standard model of phylogenomics, protein function changes rapidly after duplication, leading to paralogs with different functions, while orthologs retain the ancestral function. We report here that a comparison of experimentally supported functional annotations among homologs from 13 genomes mostly supports this model. We show that to analyze GO annotation effectively, several confounding factors need to be controlled: authorship bias, variation of GO term frequency among species, variation of background similarity among species pairs, and propagated annotation bias. After controlling for these biases, we observe that orthologs have generally more similar functional annotations than paralogs. This is especially strong for sub-cellular localization. We observe only a weak decrease in functional similarity with increasing sequence divergence. These findings hold over a large diversity of species; notably orthologs from model organisms such as E. coli, yeast or mouse have conserved function with human proteins.

Evolution of the epithelial sodium channel and the sodium pump as limiting factors of aldosterone action on sodium transport.

Despite large changes in salt intake, the mammalian kidney is able to maintain the extracellular sodium concentration and osmolarity within very narrow margins, thereby controlling blood volume and blood pressure. In the aldosterone-sensitive distal nephron (ASDN), aldosterone tightly controls the activities of epithelial sodium channel (ENaC) and Na,K-ATPase, the two limiting factors in establishing transepithelial sodium transport. It has been proposed that the ENaC/degenerin gene family is restricted to Metazoans, whereas the α- and ß-subunits of Na,K-ATPase have homologous genes in prokaryotes. This raises the question of the emergence of osmolarity control. By exploring recent genomic data of diverse organisms, we found that: 1) ENaC/degenerin exists in all of the Metazoans screened, including nonbilaterians and, by extension, was already present in ancestors of Metazoa; 2) ENaC/degenerin is also present in Naegleria gruberi, an eukaryotic microbe, consistent with either a vertical inheritance from the last common ancestor of Eukaryotes or a lateral transfer between Naegleria and Metazoan ancestors; and 3) The Na,K-ATPase ß-subunit is restricted to Holozoa, the taxon that includes animals and their closest single-cell relatives. Since the ß-subunit of Na,K-ATPase plays a key role in targeting the α-subunit to the plasma membrane and has an additional function in the formation of cell junctions, we propose that the emergence of Na,K-ATPase, together with ENaC/degenerin, is linked to the development of multicellularity in the Metazoan kingdom. The establishment of multicellularity and the associated extracellular compartment ("internal milieu") precedes the emergence of other key elements of the aldosterone signaling pathway.