Post-duplication charge evolution of phosphoglucose isomerases in teleost fishes through weak selection on many amino acid sites

When a foreign gene meets its native counterpart: computational biophysics analysis of two PgiC loci in the grass Festuca ovina

Duplicative horizontal gene transfer may bring two previously separated homologous genes together, which may raise questions about the interplay between the gene products. One such gene pair is the “native” PgiC1 and “foreign” PgiC2 in the perennial grass Festuca ovina. Both PgiC1 and PgiC2 encode cytosolic phosphoglucose isomerase, a dimeric enzyme whose proper binding is functionally essential. Here, we use biophysical simulations to explore the inter-monomer binding of the two homodimers and the heterodimer that can be produced by PgiC1 and PgiC2 in F. ovina. Using simulated native-state ensembles, we examine the structural properties and binding tightness of the dimers. In addition, we investigate their ability to withstand dissociation when pulled by a force. Our results suggest that the inter-monomer binding is tighter in the PgiC2 than the PgiC1 homodimer, which could explain the more frequent occurrence of the foreign PgiC2 homodimer in dry habitats. We further find that the PgiC1 and PgiC2 monomers are compatible with heterodimer formation; the computed binding tightness is comparable to that of the PgiC1 homodimer. Enhanced homodimer stability and capability of heterodimer formation with PgiC1 are properties of PgiC2 that may contribute to the retaining of the otherwise redundant PgiC2 gene.

An inter-switch between hydrophobic and charged amino acids generated druggable small molecule binding pocket in chemokine paralog CXCL3

Multigene families such as chemokines arose as a result of gene duplication events, followed by mutations and selection. GRO chemokines are three duplicated CXCL genes, comprising of CXCL1, CXCL2 and CXCL3 proteins. Comparative structural analysis of the two closely related paralog chemokines CXCL2 and CXCL3 in the current study indicated a variable electrostatic surface between them, and a specific hydrophobic pocket on the surface of CXCL3 that can bind naphthalene derivatives. Combined fluorescence and NMR analyses revealed that CXCL3 monomer can specifically bind to ANS (8-Anilinonaphthalene-1-sulfonic acid) with a stoichiometry of 1:1 by involving the residues belonging to the structural elements 3₁₀ helix and the α-helix. A close observation of the surfaces of these paralogs suggested that such a hydrophobic pocket is a resultant of inter-switch between a charged and a hydrophobic residue on the primary sequence of the two paralog proteins. Interestingly, the hydrophobic pocket is in the vicinity of GAG binding region of CXCL3, a molecular determinant in leukocyte trafficking. Such unique pockets/patches on specific chemokine surfaces can be exploited to design the naphthalene/small molecule based inhibitors against GAG binding to regulate their molecular interactions during the onset and progression of various types of cancers and inflammatory diseases.

Mechanistic insights into molecular evolution of species-specific differential glycosaminoglycan binding surfaces in growth-related oncogene chemokines

Chemokines are chemotactic cytokines involved in leucocyte trafficking to infected tissue. Growth-related oncogene (GRO) chemokines namely CXCL1, CXCL2 and CXCL3 are neutrophil activating chemokines sharing a conserved three-dimensional structure, but encompassing functional diversity due to gene duplication and evolutionary events. However, the evolutionary mechanisms including selection pressures involved in diversification of GRO genes have not yet been characterized. Here, we performed comprehensive evolutionary analysis of GRO genes among different mammalian species. Phylogenetic analysis illustrated a species-specific evolution pattern. Selection analysis evidenced that these genes have undergone concerted evolution. Seventeen positively selected sites were obtained, although the majority of the protein is under purifying selection. Interestingly, these positively selected sites are more concentrated on the C-terminal/glycosaminoglycan (GAG) binding and dimerization segment compared to receptor binding domain. Substitution rate analysis confirmed the C-terminal domain of GRO genes as the highest substituted segment. Further, structural analysis established that the nucleotide alterations in the GAG binding domain are the source of surface charge modulation, thus generating the differential GAG binding surfaces and multiple binding sites as per evolutionary pressure, although the helical surface is primordial for GAG binding. Indeed, such variable electrostatic surfaces are crucial to regulate chemokine gradient formation during a host's defence against pathogens and also explain the significance of chemokine promiscuity.

Evolution of estrogen receptors in ray-finned fish and their comparative responses to estrogenic substances

In vertebrates, estrogens play fundamental roles in regulating reproductive activities through estrogen receptors (ESRs), and disruption of estrogen signaling is now of global concern for both wildlife and human health. To date, ESRs of only a limited number of species have been characterized. We investigated the functional diversity and molecular basis or ligand sensitivity of ESRs among ray-finned fish species (Actinopterygii), the most variable group within vertebrates. We cloned and characterized ESRs from several key species in the evolution of ray-finned fish including bichir (Polypteriformes, ESR1 and ESR2) at the basal lineage of ray-finned fish, and arowana (Osteoglossiformes, ESR1 and ESR2b) and eel (Anguilliformes, ESR1, ESR2a and ESR2b) both belonging to ancient early-branching lineages of teleosts, and suggest that ESR2a and ESR2b emerged through teleost-specific whole genome duplication, but an ESR1 paralogue has been lost in the early lineage of euteleost fish species. All cloned ESR isoforms showed similar responses to endogenous and synthetic steroidal estrogens, but they responded differently to non-steroidal estrogenic endocrine disrupting chemicals (EDCs) (e.g., ESR2a exhibits a weaker reporter activity compared with ESR2b). We show that variation in ligand sensitivity of ESRs can be attributed to phylogeny among species of different taxonomic groups in ray-finned fish. The molecular information provided contributes both to understanding of the comparative role of ESRs in the reproductive biology of fish and their comparative responses to EDCs.

Chromosomal-Level Assembly of the Asian Seabass Genome Using Long Sequence Reads and Multi-layered Scaffolding

We report here the ~670 Mb genome assembly of the Asian seabass (Lates calcarifer), a tropical marine teleost. We used long-read sequencing augmented by transcriptomics, optical and genetic mapping along with shared synteny from closely related fish species to derive a chromosome-level assembly with a contig N50 size over 1 Mb and scaffold N50 size over 25 Mb that span ~90% of the genome. The population structure of L. calcarifer species complex was analyzed by re-sequencing 61 individuals representing various regions across the species' native range. SNP analyses identified high levels of genetic diversity and confirmed earlier indications of a population stratification comprising three clades with signs of admixture apparent in the South-East Asian population. The quality of the Asian seabass genome assembly far exceeds that of any other fish species, and will serve as a new standard for fish genomics.

Evolutionary significance and diversification of the phosphoglucose isomerase genes in vertebrates

Background

Phosphoglucose isomerase (PGI) genes are important multifunctional proteins whose evolution has, until now, not been well elucidated because of the limited number of completely sequenced genomes. Although the multifunctionality of this gene family has been considered as an original and innate characteristic, PGI genes may have acquired novel functions through changes in coding sequences and exon/intron structure, which are known to lead to functional divergence after gene duplication. A whole-genome comparative approach was used to estimate the rates of molecular evolution of this protein family.

Results

The results confirm the presence of two isoforms in teleost fishes and only one variant in all other vertebrates. Phylogenetic reconstructions grouped the PGI genes into five main groups: lungfishes/coelacanth/cartilaginous fishes, teleost fishes, amphibians, reptiles/birds and mammals, with the teleost group being subdivided into two subclades comprising PGI1 and PGI2. This PGI partitioning into groups is consistent with the synteny and molecular evolution results based on the estimation of the ratios of nonsynonymous to synonymous changes (Ka/Ks) and divergence rates between both PGI paralogs and orthologs. Teleost PGI2 shares more similarity with the variant found in all other vertebrates, suggesting that it has less evolved than PGI1 relative to the PGI of common vertebrate ancestor.

Conclusions

The diversification of PGI genes into PGI1 and PGI2 is consistent with a teleost-specific duplication before the radiation of this lineage, and after its split from the other infraclasses of ray-finned fishes. The low average Ka/Ks ratios within teleost and mammalian lineages suggest that both PGI1 and PGI2 are functionally constrained by purifying selection and may, therefore, have the same functions. By contrast, the high average Ka/Ks ratios and divergence rates within reptiles and birds indicate that PGI may be involved in different functions. The synteny analyses show that the genomic region harbouring PGI genes has independently undergone genomic rearrangements in mammals versus the reptile/bird lineage in particular, which may have contributed to the actual functional diversification of this gene family.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-015-1683-x) contains supplementary material, which is available to authorized users.

The circadian clock of teleost fish: a comparative analysis reveals distinct fates for duplicated genes

The circadian clock is a central oscillator that coordinates endogenous rhythms. Members of six gene families underlie the metabolic machinery of this system. Although this machinery appears to correspond to a highly conserved genetic system in metazoans, it has been recognized that vertebrates possess a more diverse gene inventory than that of non-vertebrates. This difference could have originated in the two successive rounds of whole-genome duplications that took place in the common ancestor of the group. Teleost fish underwent an extra event of whole-genome duplication, which is thought to have provided an abundance of raw genetic material for the biological innovations that facilitated the radiation of the group. In this study, we assessed the relative contributions of whole-genome duplication and small-scale gene duplication to generate the repertoire of genes associated with the circadian clock of teleost fish. To achieve this goal, we annotated genes from six gene families associated with the circadian clock in eight teleost fish species, and we reconstructed their evolutionary history by inferring phylogenetic relationships. Our comparative analysis indicated that teleost species possess a variable repertoire of genes related to the circadian clock gene families and that the actual diversity of these genes has been shaped by a variety of phenomena, such as the complete deletion of ohnologs, the differential retention of genes, and lineage-specific gene duplications. From a functional perspective, the subfunctionalization of two ohnolog genes (PER1a and PER1b) in zebrafish highlights the power of whole-genome duplications to generate biological diversity.

Identification, tissue distribution and functional characterization of the ghrelin receptor in West African lungfish, Protopterus annectens

We identified two ghrelin receptor isoforms, the ghrelin receptor type-1a (GHS-R1a) and its alternative splice form (GHS-R1b) for West African lungfish, Protopterus annectens. Lungfish GHS-R1a and 1b comprised 361 and 281 amino acids, respectively. Lungfish GHS-R1a showed the highest identity to coelacanth GHS-R1a (80.4%). The highest expression of GHS-R1a mRNAs was seen in the brain, liver, ovary, heart, intestine, and gills. GHS-R1b mRNAs were also detected in the same tissues with GHS-R1a, but their expression level was 1/20 that of GHS-R1a. In human embryonic kidney 293 cells transiently expressing lungfish GHS-R1a, rat and bullfrog ghrelin, and two GHS-R1a agonists, GHRP-6 and hexarelin, increased intracellular Ca(2+) concentrations. The intensity of the Ca(2+) increases induced by GHS-R1a agonists was twice when compared to that induced by ghrelin, although the median effective doses (ED50) were similar, suggesting a long-lasting effect of GHS-R1a agonists with similar affinity. We also examined changes in the GHS-R gene expression during an eight-week estivation. Body weight was slightly lowered, but plasma sodium and glucose concentrations decreased; plasma urea concentration increased significantly 4weeks after the start of estivation. Overall, expression of GHS-R1a mRNA decreased, but changes in GHS-R1b mRNA expression were inconsistent with those of GHS-R1a during estivation, suggesting an involvement of GHS-R in energy homeostasis, as seen in mammals. Our results suggest that the ghrelin-GHS-R1a system is present in this lungfish although ghrelin has not yet been found. The structure of GHS-R1a is closer to that of tetrapods than Actinopterygian fish, indicating a process of evolution that follows the Crossopterygii such as coelacanth.

Positive selection in glycolysis among Australasian stick insects

BACKGROUND

The glycolytic pathway is central to cellular energy production. Selection on individual enzymes within glycolysis, particularly phosphoglucose isomerase (Pgi), has been associated with metabolic performance in numerous organisms. Nonetheless, how whole energy-producing pathways evolve to allow organisms to thrive in different environments and adopt new lifestyles remains little explored. The Lanceocercata radiation of Australasian stick insects includes transitions from tropical to temperate climates, lowland to alpine habitats, and winged to wingless forms. This permits a broad investigation to determine which steps within glycolysis and what sites within enzymes are the targets of positive selection. To address these questions we obtained transcript sequences from seven core glycolysis enzymes, including two Pgi paralogues, from 29 Lanceocercata species.

RESULTS

Using maximum likelihood methods a signature of positive selection was inferred in two core glycolysis enzymes. Pgi and Glyceraldehyde 3-phosphate dehydrogenase (Gaphd) genes both encode enzymes linking glycolysis to the pentose phosphate pathway. Positive selection among Pgi paralogues and orthologues predominately targets amino acids with residues exposed to the protein's surface, where changes in physical properties may alter enzyme performance.

CONCLUSION

Our results suggest that, for Lancerocercata stick insects, adaptation to new stressful lifestyles requires a balance between maintaining cellular energy production, efficiently exploiting different energy storage pools and compensating for stress-induced oxidative damage.

Insights into the evolution of the CSP gene family through the integration of evolutionary analysis and comparative protein modeling

Insect chemical communication and chemosensory systems rely on proteins coded by several gene families. Here, we have combined protein modeling with evolutionary analysis in order to study the evolution and structure of chemosensory proteins (CSPs) within arthropods and, more specifically, in ants by using the data available from sequenced genomes. Ants and other social insects are especially interesting model systems for the study of chemosensation, as they communicate in a highly complex social context and much of their communication relies on chemicals. Our ant protein models show how this complexity has shaped CSP evolution; the proteins are highly modifiable by their size, surface charge and binding pocket. Based on these findings, we divide ant CSPs into three groups: typical insect CSPs, an ancient 5-helical CSP and hymenopteran CSPs with a small binding pocket, and suggest that these groups likely serve different functions. The hymenopteran CSPs have duplicated repeatedly in individual ant lineages. In these CSPs, positive selection has driven surface charge changes, an observation which has possible implications for the interaction between CSPs and ligands or odorant receptors. Our phylogenetic analysis shows that within the Arthropoda the only highly conserved gene is the ancient 5-helical CSP, which is likely involved in an essential ubiquitous function rather than chemosensation. During insect evolution, the 6-helical CSPs have diverged and perform chemosensory functions among others. Our results contribute to the general knowledge of the structural differences between proteins underlying chemosensation and highlight those protein properties which have been affected by adaptive evolution.

Whole-genome duplication and the functional diversification of teleost fish hemoglobins

Subsequent to the two rounds of whole-genome duplication that occurred in the common ancestor of vertebrates, a third genome duplication occurred in the stem lineage of teleost fishes. This teleost-specific genome duplication (TGD) is thought to have provided genetic raw materials for the physiological, morphological, and behavioral diversification of this highly speciose group. The extreme physiological versatility of teleost fish is manifest in their diversity of blood–gas transport traits, which reflects the myriad solutions that have evolved to maintain tissue O₂ delivery in the face of changing metabolic demands and environmental O₂ availability during different ontogenetic stages. During the course of development, regulatory changes in blood–O₂ transport are mediated by the expression of multiple, functionally distinct hemoglobin (Hb) isoforms that meet the particular O₂-transport challenges encountered by the developing embryo or fetus (in viviparous or oviparous species) and in free-swimming larvae and adults. The main objective of the present study was to assess the relative contributions of whole-genome duplication, large-scale segmental duplication, and small-scale gene duplication in producing the extraordinary functional diversity of teleost Hbs. To accomplish this, we integrated phylogenetic reconstructions with analyses of conserved synteny to characterize the genomic organization and evolutionary history of the globin gene clusters of teleosts. These results were then integrated with available experimental data on functional properties and developmental patterns of stage-specific gene expression. Our results indicate that multiple α- and β-globin genes were present in the common ancestor of gars (order Lepisoteiformes) and teleosts. The comparative genomic analysis revealed that teleosts possess a dual set of TGD-derived globin gene clusters, each of which has undergone lineage-specific changes in gene content via repeated duplication and deletion events. Phylogenetic reconstructions revealed that paralogous genes convergently evolved similar functional properties in different teleost lineages. Consistent with other recent studies of globin gene family evolution in vertebrates, our results revealed evidence for repeated evolutionary transitions in the developmental regulation of Hb synthesis.

Testis-specific glyceraldehyde-3-phosphate dehydrogenase: origin and evolution

Background

Glyceraldehyde-3-phosphate dehydrogenase (GAPD) catalyses one of the glycolytic reactions and is also involved in a number of non-glycolytic processes, such as endocytosis, DNA excision repair, and induction of apoptosis. Mammals are known to possess two homologous GAPD isoenzymes: GAPD-1, a well-studied protein found in all somatic cells, and GAPD-2, which is expressed solely in testis. GAPD-2 supplies energy required for the movement of spermatozoa and is tightly bound to the sperm tail cytoskeleton by the additional N-terminal proline-rich domain absent in GAPD-1. In this study we investigate the evolutionary history of GAPD and gain some insights into specialization of GAPD-2 as a testis-specific protein.

Results

A dataset of GAPD sequences was assembled from public databases and used for phylogeny reconstruction by means of the Bayesian method. Since resolution in some clades of the obtained tree was too low, syntenic analysis was carried out to define the evolutionary history of GAPD more precisely. The performed selection tests showed that selective pressure varies across lineages and isoenzymes, as well as across different regions of the same sequences.

Conclusions

The obtained results suggest that GAPD-1 and GAPD-2 emerged after duplication during the early evolution of chordates. GAPD-2 was subsequently lost by most lineages except lizards, mammals, as well as cartilaginous and bony fishes. In reptilians and mammals, GAPD-2 specialized to a testis-specific protein and acquired the novel N-terminal proline-rich domain anchoring the protein in the sperm tail cytoskeleton. This domain is likely to have originated by exonization of a microsatellite genomic region. Recognition of the proline-rich domain by cytoskeletal proteins seems to be unspecific. Besides testis, GAPD-2 of lizards was also found in some regenerating tissues, but it lacks the proline-rich domain due to tissue-specific alternative splicing.

cDNA identification, comparison and phylogenetic aspects of lombricine kinase from two oligochaete species

Creatine kinase and arginine kinase are the typical representatives of an eight-member phosphagen kinase family, which play important roles in the cellular energy metabolism of animals. The phylum Annelida underwent a series of evolutionary processes that resulted in rapid divergence and radiation of these enzymes, producing the greatest diversity of the phosphagen kinases within this phylum. Lombricine kinase (EC 2.7.3.5) is one of such enzymes and sequence information is rather limited compared to other phosphagen kinases. This study presents data on the cDNA sequences of lombricine kinase from two oligochaete species, the California blackworm (Lumbriculus variegatus) and the sludge worm (Tubifex tubifex). The deduced amino acid sequences are analyzed and compared with other selected phosphagen kinases, including two additional lombricine kinase sequences extracted from DNA databases and provide further insights in the evolution and position of these enzymes within the phosphagen kinase family. The data confirms the presence of a deleted region within the flexible loop (the GS region) of all six examined lombricine kinases. A phylogenetic analysis of these six lombricine kinases clearly positions the enzymes together in a small subcluster within the larger creatine kinase (EC 2.7.3.2) clade.

Temporal pattern of loss/persistence of duplicate genes involved in signal transduction and metabolic pathways after teleost-specific genome duplication

Background

Recent genomic studies have revealed a teleost-specific third-round whole genome duplication (3R-WGD) event occurred in a common ancestor of teleost fishes. However, it is unclear how the genes duplicated in this event were lost or persisted during the diversification of teleosts, and therefore, how many of the duplicated genes contribute to the genetic differences among teleosts. This subject is also important for understanding the process of vertebrate evolution through WGD events. We applied a comparative evolutionary approach to this question by focusing on the genes involved in long-term potentiation, taste and olfactory transduction, and the tricarboxylic acid cycle, based on the whole genome sequences of four teleosts; zebrafish, medaka, stickleback, and green spotted puffer fish.

Results

We applied a state-of-the-art method of maximum-likelihood phylogenetic inference and conserved synteny analyses to each of 130 genes involved in the above biological systems of human. These analyses identified 116 orthologous gene groups between teleosts and tetrapods, and 45 pairs of 3R-WGD-derived duplicate genes among them. This suggests that more than half [(45×2)/(116+45)] = 56.5%) of the loci, probably more than ten thousand genes, present in a common ancestor of the four teleosts were still duplicated after the 3R-WGD. The estimated temporal pattern of gene loss suggested that, after the 3R-WGD, many (71/116) of the duplicated genes were rapidly lost during the initial 75 million years (MY), whereas on average more than half (27.3/45) of the duplicated genes remaining in the ancestor of the four teleosts (45/116) have persisted for about 275 MY. The 3R-WGD-derived duplicates that have persisted for a long evolutionary periods of time had significantly larger number of interacting partners and longer length of protein coding sequence, implying that they tend to be more multifunctional than the singletons after the 3R-WGD.

Conclusion

We have shown firstly the temporal pattern of gene loss process after 3R-WGD on the basis of teleost phylogeny and divergence time frameworks. The 3R-WGD-derived duplicates have not undergone constant exponential decay, suggesting that selection favoured the long-term persistence of a subset of duplicates that tend to be multi-functional. On the basis of these results obtained from the analysis of 116 orthologous gene groups, we propose that more than ten thousand of 3R-WGD-derived duplicates have experienced lineage-specific evolution, that is, the differential sub-/neo-functionalization or secondary loss between lineages, and contributed to teleost diversity.

Evolution of multiple phosphodiesterase isoforms in stickleback involved in cAMP signal transduction pathway

Background

Duplicate genes are considered to have evolved through the partitioning of ancestral functions among duplicates (subfunctionalization) and/or the acquisition of novel functions from a beneficial mutation (neofunctionalization). Additionally, an increase in gene dosage resulting from duplication may also confer an advantageous effect, as has been suggested for histone, tRNA, and rRNA genes. Currently, there is little understanding of the effect of increased gene dosage on subcellular networks like signal transduction pathways. Addressing this issue may provide further insights into the evolution by gene duplication.

Results

We analyzed the evolution of multiple stickleback phosphodiesterase (PDE, EC: 3.1.4.17) 1C genes involved in the cyclic nucleotide signaling pathway. Stickleback has 8–9 copies of this gene, whereas only one or two loci exist in other model vertebrates. Our phylogenetic and synteny analyses suggested that the multiple PDE1C genes in stickleback were generated by repeated duplications of >100-kbp chromosome segments. Sequence evolution analysis did not provide strong evidence for neofunctionalization in the coding sequences of stickleback PDE1C isoforms. On the other hand, gene expression analysis suggested that the derived isoforms acquired expression in new organs, implying their neofunctionalization in terms of expression patterns. In addition, at least seven isoforms of the stickleback PDE1C were co-expressed with olfactory-type G-proteins in the nose, suggesting that PDE1C dosage is increased in the stickleback olfactory transduction (OT) pathway. In silico simulations of OT implied that the increased PDE1C dosage extends the longevity of the depolarization signals of the olfactory receptor neuron.

Conclusion

The predicted effect of the increase in PDE1C products on the OT pathway may play an important role in stickleback behavior and ecology. However, this possibility should be empirically examined. Our analyses imply that an increase in gene product sometimes has a significant, yet unexpected, effect on the functions of subcellular networks.

Electric charge divergence in proteins: insights into the evolution of their three-dimensional properties

Previous studies of protein evolution have identified important mutations in various proteins that affect a small number of residues, but dramatically alter protein function. However, the evolutionary process underlying the three-dimensional protein properties, which are determined by a much larger number of residues, remains unclear. Based on a comparative evolutionary analysis of teleost phosphoglucose isomerases (PGIs; EC 5.3.1.9), we previously demonstrated that the relatively weak selection on many amino acid sites has played an important role in the evolution of protein electric charge as a model of three-dimensional protein properties. To ascertain the generality of this finding, we sought further evidence of this type of protein evolution. For this purpose, we analyzed the vertebrate isoforms of fructose-1,6-bisphosphate aldolase (ALD; EC 4.1.2.13), for which electric charges are known to have diverged after gene duplication. The results showed that the divergence in electric charge between the ALD isoforms was also driven by weak selection on many amino acid sites, as in PGI, confirming the generality of earlier findings. To obtain further insights, ALD and PGI were compared to the proteins pancreatic ribonuclease (EC 3.1.27.5) and triose-phosphate isomerase (EC 5.3.1.1), for which electric charges likely evolved through a well-defined mode of molecular evolution; namely, strong selection on specific amino acid sites. Comparison of the number and composition of amino acids on the protein surface suggested that the absolute number of evolutionarily changeable amino acids in a protein affects the strength of selection pressure acting on individual amino acid sites.