Unraveling the functional divergence of membrane-bound pyrophosphatases

Inorganic pyrophosphatases (PPases) are enzymes that hydrolyze pyrophosphate (PPi)which is produced as a byproduct in many important growth related processes e.g. in the biosynthesis of DNA, proteins and lipids. PPases can be either soluble or membranebound. Membrane-bound PPases (mPPases) are ion transporters that couple the energy released during PPi hydrolysis to Na+ or H+ transport. When I started the project, only three Na+-transporting mPPases were known to exist. In this study, I aimed to confirm if Na+-transport is a common function of mPPases. Furthermore, the amino acid residues responsible for determining the transporter specificity were unknown. I constructed a phylogenetic tree for mPPases and selected the representative bacterial and archaeal mPPases to be investigated. I expressed different prokaryotic mPPases in Escherichia coli, isolated these as inverted membrane vesicles and characterized their functions. In the first project I identified four new Na+-PPases, two K+-dependent H+-PPases and one K+-independent mPPase. The residues determining the transporter specificity were identified by site-directed mutagenesis. I showed that the conserved glutamate residues are important for specificity, though are not the only residues that influence it. This research clarified the ion transport specificities throughout the mPPase phylogenetic tree, and revealed that Na+ transport is a widespread function of mPPases. In addition, it became clear that the transporter specificity can be predicted from the amino acid sequence in combination with a phylogenetic analysis. In the second project, I identified a novel class of mPPases, which is capable of transporting both Na+ and H+ ions and is mainly found in bacteria of the human gastrointestinal tract. The physiological role of these novel enzymes may be to help the bacteria survive in the demanding conditions of the host. In the third project, I characterized the Chlorobium limicola Na+-PPase and found that this and related mPPases are able to transport H+ ions at subphysiological Na+ concentrations. In addition, the H+-transport activity was shown to be a common function of all studied Na+-PPases at low Na+ concentrations. I observed that mutating gate-lysine to asparagine eliminated the H+ but not the Na+ ion transport function, indicating the important role of the residue in the transport of H+. In the fourth project, I characterized the unknown and evolutionary divergent mPPase clade of the phylogenetic tree. The enzymes belonging to this clade are able to transport H+ ions and, based on their sequence, were expected to be K+- and Na+-independent. The sequences of membrane-bound PPase are usually highly conserved, but the enzymes belonging to this clade are more divergent and usually contain 100−150 extra amino acid residues compared to other known mPPases. Despite the vast sequence differences, these mPPases have the full set of important residues and, surprisingly, are regulated by Na+ and K+ ions. These enzymes are mainly of bacterial origin.

Duplication and functional divergence in the chalcone synthase gene family of Asteraceae: evolution with substrate change and catalytic simplification.

Plant-specific polyketide synthase genes constitute a gene superfamily, including universal chalcone synthase [CHS; malonyl-CoA:4-coumaroyl-CoA malonyltransferase (cyclizing) (EC 2.3.1.74)] genes, sporadically distributed stilbene synthase (SS) genes, and atypical, as-yet-uncharacterized CHS-like genes. We have recently isolated from Gerbera hybrida (Asteraceae) an unusual CHS-like gene, GCHS2, which codes for an enzyme with structural and enzymatic properties as well as ontogenetic distribution distinct from both CHS and SS. Here, we show that the GCHS2-like function is encoded in the Gerbera genome by a family of at least three transcriptionally active genes. Conservation within the GCHS2 family was exploited with selective PCR to study the occurrence of GCHS2-like genes in other Asteraceae. Parsimony analysis of the amplified sequences together with CHS-like genes isolated from other taxa of angiosperm subclass Asteridae suggests that GCHS2 has evolved from CHS via a gene duplication event that occurred before the diversification of the Asteraceae. Enzyme activity analysis of proteins produced in vitro indicates that the GCHS2 reaction is a non-SS variant of the CHS reaction, with both different substrate specificity (to benzoyl-CoA) and a truncated catalytic profile. Together with the recent results of Durbin et al. [Durbin, M. L., Learn, G. H., Jr., Huttley, G. A. & Clegg, M. T. (1995) Proc. Natl. Acad. Sci. USA 92, 3338-3342], our study confirms a gene duplication-based model that explains how various related functions have arisen from CHS during plant evolution.

Developmental and Environmental Regulation of Aquaporin Gene Expression across Populus Species: Divergence or Redundancy?

Aquaporins (AQPs) are membrane channels belonging to the major intrinsic proteins family and are known for their ability to facilitate water movement. While in Populus trichocarpa, AQP proteins form a large family encompassing fifty-five genes, most of the experimental work focused on a few genes or subfamilies. The current work was undertaken to develop a comprehensive picture of the whole AQP gene family in Populus species by delineating gene expression domain and distinguishing responsiveness to developmental and environmental cues. Since duplication events amplified the poplar AQP family, we addressed the question of expression redundancy between gene duplicates. On these purposes, we carried a meta-analysis of all publicly available Affymetrix experiments. Our in-silico strategy controlled for previously identified biases in cross-species transcriptomics, a necessary step for any comparative transcriptomics based on multispecies design chips. Three poplar AQPs were not supported by any expression data, even in a large collection of situations (abiotic and biotic constraints, temporal oscillations and mutants). The expression of 11 AQPs was never or poorly regulated whatever the wideness of their expression domain and their expression level. Our work highlighted that PtTIP1;4 was the most responsive gene of the AQP family. A high functional divergence between gene duplicates was detected across species and in response to tested cues, except for the root-expressed PtTIP2;3/PtTIP2;4 pair exhibiting 80% convergent responses. Our meta-analysis assessed key features of aquaporin expression which had remained hidden in single experiments, such as expression wideness, response specificity and genotype and environment interactions. By consolidating expression profiles using independent experimental series, we showed that the large expansion of AQP family in poplar was accompanied with a strong divergence of gene expression, even if some cases of functional redundancy could be suspected.

The UlaG protein familly defines novel structural and functional motifs grafted on an ancient RNase fold

Background: Bacterial populations are highly successful at colonizing new habitats and adapting to changing environmental conditions, partly due to their capacity to evolve novel virulence and metabolic pathways in response to stress conditions and to shuffle them by horizontal gene transfer (HGT). A common theme in the evolution of new functions consists of gene duplication followed by functional divergence. UlaG, a unique manganese-dependent metallo-b-lactamase (MBL) enzyme involved in L-ascorbate metabolism by commensal and symbiotic enterobacteria, provides a model for the study of the emergence of new catalytic activities from the modification of an ancient fold. Furthermore, UlaG is the founding member of the so-called UlaG-like (UlaGL) protein family, a recently established and poorly characterized family comprising divalent (and perhaps trivalent)metal-binding MBLs that catalyze transformations on phosphorylated sugars and nucleotides. Results: Here we combined protein structure-guided and sequence-only molecular phylogenetic analyses to dissect the molecular evolution of UlaG and to study its phylogenomic distribution, its relatedness with present-day UlaGL protein sequences and functional conservation. Phylogenetic analyses indicate that UlaGL sequences are present in Bacteria and Archaea, with bona fide orthologs found mainly in mammalian and plant-associated Gramnegative and Gram-positive bacteria. The incongruence between the UlaGL tree and known species trees indicates exchange by HGT and suggests that the UlaGL-encoding genes provided a growth advantage under changing conditions. Our search for more distantly related protein sequences aided by structural homology has uncovered that UlaGL sequences have a common evolutionary origin with present-day RNA processing and metabolizing MBL enzymes widespread in Bacteria, Archaea, and Eukarya. This observation suggests an ancient origin for the UlaGL family within the broader trunk of the MBL superfamily by duplication, neofunctionalization and fixation. Conclusions: Our results suggest that the forerunner of UlaG was present as an RNA metabolizing enzyme in the last common ancestor, and that the modern descendants of that ancestral gene have a wide phylogenetic distribution and functional roles. We propose that the UlaGL family evolved new metabolic roles among bacterial and possibly archeal phyla in the setting of a close association with metazoans, such as in the mammalian gastrointestinal tract or in animal and plant pathogens, as well as in environmental settings. Accordingly, the major evolutionary forces shaping the UlaGL family include vertical inheritance and lineage-specific duplication and acquisition of novel metabolic functions, followed by HGT and numerous lineage-specific gene loss events.

Functional and evolutionary analysis of DXL1, a non-essential gene encoding a 1-deoxy-D-xylulose 5-phosphate synthase like protein in Arabidopsis thaliana

The synthesis of 1-deoxy-D-xylulose 5-phosphate (DXP), catalyzed by the enzyme DXP synthase (DXS), represents a key regulatory step of the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway for isoprenoid biosynthesis. In plants DXS is encoded by small multigene families that can be classified into, at least, three specialized subfamilies. Arabidopsis thaliana contains three genes encoding proteins with similarity to DXS, including the well-known DXS1/CLA1 gene, which clusters within subfamily I. The remaining proteins, initially named DXS2 and DXS3, have not yet been characterized. Here we report the expression and functional analysis of A. thaliana DXS2. Unexpectedly, the expression of DXS2 failed to rescue Escherichia coli and A. thaliana mutants defective in DXS activity. Coherently, we found that DXS activity was negligible in vitro, being renamed as DXL1 following recent nomenclature recommendation. DXL1 is targeted to plastids as DXS1, but shows a distinct expression pattern. The phenotypic analysis of a DXL1 defective mutant revealed that the function of the encoded protein is not essential for growth and development. Evolutionary analyses indicated that DXL1 emerged from DXS1 through a recent duplication apparently specific of the Brassicaceae lineage. Divergent selective constraints would have affected a significant fraction of sites after diversification of the paralogues. Furthermore, amino acids subjected to divergent selection and likely critical for functional divergence through the acquisition of a novel, although not yet known, biochemical function, were identified. Our results provide with the first evidences of functional specialization at both the regulatory and biochemical level within the plant DXS family.

Responses of aquatic insect functional diversity to landscape changes in Atlantic forest

The Atlantic Forest domain, one of the 25 world's hotspots for biodiversity, has experienced dramatic changes in its landscape. While the loss of species diversity is well documented, functional diversity has not received the same amount of attention. In this study, we evaluated functional diversity of insects in streams utilizing three indices: functional diversity (FD), functional dispersion (FDis), and functional divergence (FDiv), seeking to understand the roles of three predictor sets in explaining functional patterns: (1) bioclimatic and landscape variables; (2) spatial variables; and (3) local environmental variables. We determined the amount of variation in different measures of functional diversity that was explained by each predictor set and their interplays using variation partitioning. Our study showed that variation in functional diversity is better explained by a set of variables linked to different scales dependent on spatial structures, indicating the importance of landscape and mainly environmental variables in the functional organization of aquatic insect communities, and that the relative importance of predictor sets depends on the indices considered. Variation in FD was better explained by the interplay among the three predictor sets and by local environmental variables, whereas variation in FDis was better explained by spatial variables and by the interplay between environmental and spatial variables. Variation in FDiv was not significantly explained by any predictors. Our study adds more evidence on the harmful effects caused by landscape changes on biodiversity in the Atlantic Forest, suggesting that these effects also influence the functional organization of stream insect communities. © 2013 The Author(s) Journal compilation © 2013 by The Association for Tropical Biology and Conservation.

Cladoceran functional assemblages and functional diversity indices from Mallusjärvi and Tiiläänjärvi (southern Finland) sediment cores

will be submitted by the author

Evidence for an episodic model of protein sequence evolution.

The evolution of protein function appears to involve alternating periods of conservative evolution and of relatively rapid change. Evidence for such episodic evolution, consistent with some theoretical expectations, comes from the application of increasingly sophisticated models of evolution to large sequence datasets. We present here some of the recent methods to detect functional shifts, using amino acid or codon models. Both provide evidence for punctual shifts in patterns of amino acid conservation, including the fixation of key changes by positive selection. Although a link to gene duplication, a presumed source of functional changes, has been difficult to establish, this episodic model appears to apply to a wide variety of proteins and organisms.

Large-scale analysis of orthologs and paralogs under covarion-like and constant-but-different models of amino acid evolution.

Functional divergence between homologous proteins is expected to affect amino acid sequences in two main ways, which can be considered as proxies of biochemical divergence: a "covarion-like" pattern of correlated changes in evolutionary rates, and switches in conserved residues ("conserved but different"). Although these patterns have been used in case studies, a large-scale analysis is needed to estimate their frequency and distribution. We use a phylogenomic framework of animal genes to answer three questions: 1) What is the prevalence of such patterns? 2) Can we link such patterns at the amino acid level with selection inferred at the codon level? 3) Are patterns different between paralogs and orthologs? We find that covarion-like patterns are more frequently detected than "constant but different," but that only the latter are correlated with signal for positive selection. Finally, there is no obvious difference in patterns between orthologs and paralogs.

Comparative genomic analysis of the odorant-binding protein family in 12 Drosophila genomes: purifying selection and birth-and-death evolution

Background: Chemoreception is a widespread mechanism that is involved in critical biologic processes, including individual and social behavior. The insect peripheral olfactory system comprises three major multigene families: the olfactory receptor (Or), the gustatory receptor (Gr), and the odorant-binding protein (OBP) families. Members of the latter family establish the first contact with the odorants, and thus constitute the first step in the chemosensory transduction pathway.Results: Comparative analysis of the OBP family in 12 Drosophila genomes allowed the identification of 595 genes that encode putative functional and nonfunctional members in extant species, with 43 gene gains and 28 gene losses (15 deletions and 13 pseudogenization events). The evolution of this family shows tandem gene duplication events, progressive divergence in DNA and amino acid sequence, and prevalence of pseudogenization events in external branches of the phylogenetic tree. We observed that the OBP arrangement in clusters is maintained across the Drosophila species and that purifying selection governs the evolution of the family; nevertheless, OBP genes differ in their functional constraints levels. Finally, we detect that the OBP repertoire evolves more rapidly in the specialist lineages of the Drosophila melanogaster group (D. sechellia and D. erecta) than in their closest generalists.Conclusion: Overall, the evolution of the OBP multigene family is consistent with the birth-and-death model. We also found that members of this family exhibit different functional constraints, which is indicative of some functional divergence, and that they might be involved in some of the specialization processes that occurred through the diversification of the Drosophila genus.

Comparative genomics of Hedgehog loci in chordates and the origin of Shh regulatory novelties

The origin and evolution of the complex regulatory landscapes of some vertebrate developmental genes, often spanning hundreds of Kbp and including neighboring genes, remain poorly understood. The Sonic Hedgehog (Shh) genomic regulatory block (GRB) is one of the best functionally characterized examples, with several discrete enhancers reported within its introns, vast upstream gene-free region and neighboring genes (Lmbr1 and Rnf32). To investigate the origin and evolution of this GRB, we sequenced and characterized the Hedgehog (Hh) loci from three invertebrate chordate amphioxus species, which share several early expression domains with Shh. Using phylogenetic footprinting within and between chordate lineages, and reporter assays in zebrafish probing >30 Kbp of amphioxus Hh, we report large sequence and functional divergence between both groups. In addition, we show that the linkage of Shh to Lmbr1 and Rnf32, necessary for the unique gnatostomate-specific Shh limb expression, is a vertebrate novelty occurred between the two whole-genome duplications.

Molecular Evolution of Metazoan Hypoxia-inducing Factors

Most metazoans rely on aerobic energy production, which is dependent on adequate oxygen supply. In the case of reduced oxygen supply (hypoxia), the most profound changes in gene expression are mediated by transcription factors named hypoxia-inducible factors (HIF alpha). These proteins are post-translationally regulated by prolyl-4-hydroxylase (PHD) enzymes that are direct “sensors” of cellular oxygen levels. This thesis examines the molecular evolution of metazoan HIF systems. In early metazoans the HIF system emerged from pre-existing PHD oxygen sensors and early bHLH-PAS transcription factors. In invertebrates our analysis revealed an unexpected diversity of PHD genes and HIF alpha sequence characteristics. An early branching vertebrate, the epaulette shark (Hemiscyllium ocellatum) was chosen for sequencing and hypoxia preconditioning studies of HIF alpha and PHD genes. As no quantitative PCR reference genes were available, this thesis includes the first study of reference genes in cartilaginous fish species. Applying multiple statistical analysis we also discoveredthat commonly used reference gene software may perform poorly with some data sets. Novel reference genes allowed accurate measurements of the mRNAlevels of the studied target genes. Cartilaginous fishes have three genomic duplicates of both HIF alpha and PHD genes like mammals and teleost fishes. Combining functional divergence and selection analyses it was possible to describe how sequence changes in both HIF alpha and PHD duplicates may have contributed to the differential oxygen sensitivityof HIF alphas. Additionally, novel teleost HIF-1 alpha sequences were produced and used to reveal the molecular evolution of HIF-1 alpha in this lineage rich with hypoxia tolerant species.

Structural Aspects of the Distinct Biochemical Properties of Glutaredoxin 1 and Glutaredoxin 2 from Saccharomyces cerevisiae

Glutaredoxins (Grxs) are small (9-12 kDa) heat-stable proteins that are ubiquitously distributed. In Saccharomyces cerevisiae, seven Grx enzymes have been identified. Two of them (yGrx1 and yGrx2) are dithiolic, possessing a conserved Cys-Pro-Tyr-Cys motif. Here, we show that yGrx2 has a specific activity 15 times higher than that of yGrx1, although these two oxidoreductases share 64% identity and 85% similarity with respect to their amino acid sequences. Further characterization of the enzymatic activities through two-substrate kinetics analysis revealed that yGrx2 possesses a lower Km for glutathione and a higher turnover than yGrx1. To better comprehend these biochemical differences, the pK(a) of the N-terminal active-site cysteines (Cys27) of these two proteins and of the yGrx2-C30S mutant were determined. Since the pK(a) values of the yGrx1 and yGix2 Cys27 residues are very similar, these parameters cannot account for the difference observed between their specific activities. Therefore, crystal structures of yGrx2 in the oxidized form and with a glutathionyl mixed disulfide were determined at resolutions of 2.05 and 1.91 angstrom, respectively. Comparisons of yGrx2 structures with the recently determined structures of yGrx1 provided insights into their remarkable functional divergence. We hypothesize that the substitutions of Ser23 and Gln52 in yGrx1 by Ala23 and Glu52 in yGrx2 modify the capability of the active-site C-terminal cysteine to attack the mixed disulfide between the N-terminal active-site cysteine and the glutathione molecule. Mutagenesis studies supported this hypothesis. The observed structural and functional differences between yGrx1 and yGrx2 may reflect variations in substrate specificity. (C) 2008 Elsevier Ltd. All rights reserved.