The globin gene repertoire of lampreys: convergent evolution of hemoglobin and myoglobin in jawed and jawless vertebrates

The knockout of cytoglobin 1 in zebrafish (Danio rerio) alters lipid metabolism, iron homeostasis and oxidative stress response

Cytoglobin (Cygb) is an evolutionary ancient heme protein with yet unclear physiological function(s). Mammalian Cygb is ubiquitously expressed in all tissues and is proposed to be involved in reactive oxygen species (ROS) detoxification, nitric oxide (NO) metabolism and lipid-based signaling processes. Loss-of-function studies in mouse associate Cygb with apoptosis, inflammation, fibrosis, cardiovascular dysfunction or oncogenesis. In zebrafish (Danio rerio), two cygb genes exist, cytoglobin 1 (cygb1) and cytoglobin 2 (cygb2). Both have different coordination states and distinct expression sites within zebrafish tissues. The biological roles of the cygb paralogs are largely uncharacterized. We used a CRISPR/Cas9 genome editing approach and generated a knockout of the penta-coordinated cygb1 for in vivo analysis. Adult male cygb1 knockouts develop phenotypic abnormalities, including weight loss. To identify the molecular mechanisms underlying the occurrence of these phenotypes and differentiate between function and effect of the knockout we compared the transcriptomes of cygb1 knockout at different ages to age-matched wild-type zebrafish. We found that immune regulatory and cell cycle regulatory transcripts (e.g. tp53) were up-regulated in the cygb1 knockout liver. Additionally, the expression of transcripts involved in lipid metabolism and transport, the antioxidative defense and iron homeostasis was affected in the cygb1 knockout. Cygb1 may function as an anti-inflammatory and cytoprotective factor in zebrafish liver, and may be involved in lipid-, iron-, and ROS-dependent signaling.

Ancient genomic linkage couples metabolism with erythroid development

Generation of mature cells from progenitors requires tight coupling of differentiation and metabolism. During erythropoiesis, erythroblasts are required to massively upregulate globin synthesis then clear extraneous material and enucleate to produce erythrocytes1-3. Nprl3 has remained in synteny with the α-globin genes for >500 million years4, and harbours the majority of the α-globin enhancers5. Nprl3 is a highly conserved inhibitor of mTORC1, which controls cellular metabolism. However, whether Nprl3 itself serves an erythroid role is unknown. Here, we show that Nprl3 is a key regulator of erythroid metabolism. Using Nprl3-deficient fetal liver and adult competitive bone marrow - fetal liver chimeras, we show that NprI3 is required for sufficient erythropoiesis. Loss of Nprl3 elevates mTORC1 signalling, suppresses autophagy and disrupts erythroblast glycolysis and redox control. Human CD34+ progenitors lacking NPRL3 produce fewer enucleated cells and demonstrate dysregulated mTORC1 signalling in response to nutrient availability and erythropoietin. Finally, we show that the α-globin enhancers upregulate NprI3 expression, and that this activity is necessary for optimal erythropoiesis. Therefore, the anciently conserved linkage of NprI3, α-globin and their associated enhancers has enabled coupling of metabolic and developmental control in erythroid cells. This may enable erythropoiesis to adapt to fluctuating nutritional and environmental conditions.

Molecular evolution of the hemoglobin gene family across vertebrates

Adaptation to various altitudes and oxygen levels is a major aspect of vertebrate evolution. Hemoglobin is an erythrocyte protein belonging to the globin superfamily, and the α-, β-globin genes of jawed vertebrates encode tetrameric ((α₂β₂) hemoglobin, which contributes to aerobic metabolism by delivering oxygen from the respiratory exchange surfaces into cells. However, there are various gaps in knowledge regarding hemoglobin gene evolution, including patterns in cartilaginous fish and the roles of gene conversion in various taxa. Hence, we evaluated the evolutionary history of the vertebrate hemoglobin gene family by analyses of 97 species representing all classes of vertebrates. By genome-wide analyses, we extracted 879 hemoglobin sequences. Members of the hemoglobin gene family were conserved in birds and reptiles but variable in mammals, amphibians, and teleosts. Gene motifs, structures, and synteny were relatively well-conserved among vertebrates. Our results revealed that purifying selection contributed substantially to the evolution of all vertebrate hemoglobin genes, with mean d_N/d_S (ω) values ranging from 0.057 in teleosts to 0.359 in reptiles. In general, after the fish-specific genome duplication, the teleost hemoglobin genes showed variation in rates of evolution, and the β-globin genes showed relatively high ω values after a gene transposition event in amniotes. We also observed that the frequency of gene conversion was high in amniotes, with fewer hemoglobin genes and higher rates of evolution. Collectively, our findings provide detail insight into complex evolutionary processes shaping the vertebrate hemoglobin gene family, involving gene duplication, gene loss, purifying selection, and gene conversion.

Structural evolution of an amphibian-specific globin: A computational evolutionary biochemistry approach

Studies on the globin family are continuously revealing insights into the mechanisms of gene and protein evolution. The rise of a new globin gene type in Pelobatoidea and Neobatrachia (Amphibia:Anura) from an α-globin precursor provides the opportunity to investigate the genetic and physical mechanisms underlying the origin of new protein structural and functional properties. This amphibian-specific globin (globin A/GbA) discovered in the heart of Rana catesbeiana is a monomer. As the ancestral oligomeric state of α-globins is a homodimer, we inferred that the ancestral state was lost somewhere in the GbA lineage. Here, we combined computational molecular evolution with structural bioinformatics to determine the extent to which the loss of the homodimeric state is pervasive in the GbA clade. We also characterized the loci of GbA genes in Bufo bufo. We found two GbA clades in Neobatrachia. One was deleted in Ranidae, but retained and expanded to yield a new globin cluster in Bufonidae species. Loss of the ancestral oligomeric state seems to be pervasive in the GbA clade. However, a taxonomic sampling that includes more Pelobatoidea, as well as early Neobatrachia, lineages would be necessary to determine the oligomeric state of the last common ancestor of all GbA. The evidence presented here points out a possible loss of oligomerization in Pelobatoidea GbA as a result of amino acid substitutions that weaken the homodimeric state. In contrast, the loss of oligomerization in both Neobatrachia GbA clades was linked to independent deletions that disrupted many packing contacts at the homodimer interface.

Whole genome sequencing of the fast-swimming Southern bluefin tuna (Thunnus maccoyii)

The economically important Southern bluefin tuna (Thunnus maccoyii) is a world-famous fast-swimming fish, but its genomic information is limited. Here, we performed whole genome sequencing and assembled a draft genome for Southern bluefin tuna, aiming to generate useful genetic data for comparative functional prediction. The final genome assembly is 806.54 Mb, with scaffold and contig N50 values of 3.31 Mb and 67.38 kb, respectively. Genome completeness was evaluated to be 95.8%. The assembled genome contained 23,403 protein-coding genes and 236.1 Mb of repeat sequences (accounting for 29.27% of the entire assembly). Comparative genomics analyses of this fast-swimming tuna revealed that it had more than twice as many hemoglobin genes (18) as other relatively slow-moving fishes (such as seahorse, sunfish, and tongue sole). These hemoglobin genes are mainly localized in two big clusters (termed as “MNˮ and “LAˮ respectively), which is consistent with other reported fishes. However, Thr39 of beta-hemoglobin in the MN cluster, conserved in other fishes, was mutated as cysteine in tunas including the Southern bluefin tuna. Since hemoglobins are reported to transport oxygen efficiently for aerobic respiration, our genomic data suggest that both high copy numbers of hemoglobin genes and an adjusted function of the beta-hemoglobin may support the fast-swimming activity of tunas. In summary, we produced a primary genome assembly and predicted hemoglobin-related roles for the fast-swimming Southern bluefin tuna.

The More, the Merrier? Multiple Myoglobin Genes in Fish Species, Especially in Gray Bichir (Polypterus senegalus) and Reedfish (Erpetoichthys calabaricus)

The members of the globin superfamily are a classical model system to investigate gene evolution and their fates as well as the diversity of protein function. One of the best-known globins is myoglobin (Mb), which is mainly expressed in heart muscle and transports oxygen from the sarcolemma to the mitochondria. Most vertebrates harbor a single copy of the myoglobin gene, but some fish species have multiple myoglobin genes. Phylogenetic analyses indicate an independent emergence of multiple myoglobin genes, whereby the origin is mostly the last common ancestor of each order. By analyzing different transcriptome data sets, we found at least 15 multiple myoglobin genes in the polypterid gray bichir (Polypterus senegalus) and reedfish (Erpetoichthys calabaricus). In reedfish, the myoglobin genes are expressed in a broad range of tissues but show very different expression values. In contrast, the Mb genes of the gray bichir show a rather scattered expression pattern; only a few Mb genes were found expressed in the analyzed tissues. Both, gray bichir and reedfish possess lungs which enable them to inhabit shallow and swampy waters throughout tropical Africa with frequently fluctuating and low oxygen concentrations. The myoglobin repertoire probably reflects the molecular adaptation to these conditions. The sequence divergence, the substitution rate, and the different expression pattern of multiple myoglobin genes in gray bichir and reedfish imply different functions, probably through sub- and neofunctionalization during evolution.

Myoglobin primary structure reveals multiple convergent transitions to semi-aquatic life in the world's smallest mammalian divers

The speciose mammalian order Eulipotyphla (moles, shrews, hedgehogs, solenodons) combines an unusual diversity of semi-aquatic, semi-fossorial, and fossorial forms that arose from terrestrial forbearers. However, our understanding of the ecomorphological pathways leading to these lifestyles has been confounded by a fragmentary fossil record, unresolved phylogenetic relationships, and potential morphological convergence, calling for novel approaches. The net surface charge of the oxygen-storing muscle protein myoglobin (ZMb), which can be readily determined from its primary structure, provides an objective target to address this question due to mechanistic linkages with myoglobin concentration. Here, we generate a comprehensive 71 species molecular phylogeny that resolves previously intractable intra-family relationships and then ancestrally reconstruct ZMb evolution to identify ancient lifestyle transitions based on protein sequence alone. Our phylogenetically informed analyses confidently resolve fossorial habits having evolved twice in talpid moles and reveal five independent secondary aquatic transitions in the order housing the world's smallest endothermic divers.

Globins in the marine annelid Platynereis dumerilii shed new light on hemoglobin evolution in bilaterians

BACKGROUND

How vascular systems and their respiratory pigments evolved is still debated. While many animals present a vascular system, hemoglobin exists as a blood pigment only in a few groups (vertebrates, annelids, a few arthropod and mollusk species). Hemoglobins are formed of globin sub-units, belonging to multigene families, in various multimeric assemblages. It was so far unclear whether hemoglobin families from different bilaterian groups had a common origin.

RESULTS

To unravel globin evolution in bilaterians, we studied the marine annelid Platynereis dumerilii, a species with a slow evolving genome. Platynereis exhibits a closed vascular system filled with extracellular hemoglobin. Platynereis genome and transcriptomes reveal a family of 19 globins, nine of which are predicted to be extracellular. Extracellular globins are produced by specialized cells lining the vessels of the segmental appendages of the worm, serving as gills, and thus likely participate in the assembly of a previously characterized annelid-specific giant hemoglobin. Extracellular globin mRNAs are absent in smaller juveniles, accumulate considerably in growing and more active worms and peak in swarming adults, as the need for O2 culminates. Next, we conducted a metazoan-wide phylogenetic analysis of globins using data from complete genomes. We establish that five globin genes (stem globins) were present in the last common ancestor of bilaterians. Based on these results, we propose a new nomenclature of globins, with five clades. All five ancestral stem-globin clades are retained in some spiralians, while some clades disappeared early in deuterostome and ecdysozoan evolution. All known bilaterian blood globin families are grouped in a single clade (clade I) together with intracellular globins of bilaterians devoid of red blood.

CONCLUSIONS

We uncover a complex "pre-blood" evolution of globins, with an early gene radiation in ancestral bilaterians. Circulating hemoglobins in various bilaterian groups evolved convergently, presumably in correlation with animal size and activity. However, all hemoglobins derive from a clade I globin, or cytoglobin, probably involved in intracellular O2 transit and regulation. The annelid Platynereis is remarkable in having a large family of extracellular blood globins, while retaining all clades of ancestral bilaterian globins.

The rise and fall of globins in the amphibia

The globin gene repertoire of gnathostome vertebrates is dictated by differential retention and loss of nine paralogous genes: androglobin, neuroglobin, globin X, cytoglobin, globin Y, myoglobin, globin E, and the α- and β-globins. We report the globin gene repertoire of three orders of modern amphibians: Anura, Caudata, and Gymnophiona. Combining phylogenetic and conserved synteny analysis, we show that myoglobin and globin E were lost only in the Batrachia clade, but retained in Gymnophiona. The major amphibian groups also retained different paralogous copies of globin X. None of the amphibian presented α^D-globin gene. Nevertheless, two clades of β-globins are present in all amphibians, indicating that the amphibian ancestor possessed two paralogous proto β-globins. We also show that orthologs of the gene coding for the monomeric hemoglobin found in the heart of Rana catesbeiana are present in Neobatrachia and Pelobatoidea species we analyzed. We suggest that these genes might perform myoglobin- and globin E-related functions. We conclude that the repertoire of globin genes in amphibians is dictated by both retention and loss of the paralogous genes cited above and the rise of a new globin gene through co-option of an α-globin, possibly facilitated by a prior event of transposition.

Lessons from the post-genomic era: Globin diversity beyond oxygen binding and transport

Vertebrate hemoglobin (Hb) and myoglobin (Mb) were among the first proteins whose structures and sequences were determined over 50 years ago. In the subsequent pregenomic period, numerous related proteins came to light in plants, invertebrates and bacteria, that shared the myoglobin fold, a signature sequence motif characteristic of a 3-on-3 α-helical sandwich. Concomitantly, eukaryote and bacterial globins with a truncated 2-on-2 α-helical fold were discovered. Genomic information over the last 20 years has dramatically expanded the list of known globins, demonstrating their existence in a limited number of archaeal genomes, a majority of bacterial genomes and an overwhelming majority of eukaryote genomes. In vertebrates, 6 additional globin types were identified, namely neuroglobin (Ngb), cytoglobin (Cygb), globin E (GbE), globin X (GbX), globin Y (GbY) and androglobin (Adgb). Furthermore, functions beyond the familiar oxygen transport and storage have been discovered within the vertebrate globin family, including NO metabolism, peroxidase activity, scavenging of free radicals, and signaling functions. The extension of the knowledge on globin functions suggests that the original roles of bacterial globins must have been enzymatic, involved in defense against NO toxicity, and perhaps also as sensors of O₂, regulating taxis away or towards high O₂ concentrations. In this review, we aimed to discuss the evolution and remarkable functional diversity of vertebrate globins with particular focus on the variety of non-canonical expression sites of mammalian globins and their according impressive variability of atypical functions.

An evolutionarily ancient mechanism for regulation of hemoglobin expression in vertebrate red cells

The oxygen transport function of hemoglobin (HB) is thought to have arisen ∼500 million years ago, roughly coinciding with the divergence between jawless (Agnatha) and jawed (Gnathostomata) vertebrates. Intriguingly, extant HBs of jawless and jawed vertebrates were shown to have evolved twice, and independently, from different ancestral globin proteins. This raises the question of whether erythroid-specific expression of HB also evolved twice independently. In all jawed vertebrates studied to date, one of the HB gene clusters is linked to the widely expressed NPRL3 gene. Here we show that the nprl3-linked hb locus of a jawless vertebrate, the river lamprey (Lampetra fluviatilis), shares a range of structural and functional properties with the equivalent jawed vertebrate HB locus. Functional analysis demonstrates that an erythroid-specific enhancer is located in intron 7 of lamprey nprl3, which corresponds to the NPRL3 intron 7 MCS-R1 enhancer of jawed vertebrates. Collectively, our findings signify the presence of an nprl3-linked multiglobin gene locus, which contains a remote enhancer that drives globin expression in erythroid cells, before the divergence of jawless and jawed vertebrates. Different globin genes from this ancestral cluster evolved in the current NPRL3-linked HB genes in jawless and jawed vertebrates. This provides an explanation of the enigma of how, in different species, globin genes linked to the same adjacent gene could undergo convergent evolution.

Adaptations to environmental change: Globin superfamily evolution in Antarctic fishes

The ancient origins and functional versatility of globins make them ideal subjects for studying physiological adaptation to environmental change. Our goals in this review are to describe the evolution of the vertebrate globin gene superfamily and to explore the structure/function relationships of hemoglobin, myoglobin, neuroglobin and cytoglobin in teleost fishes. We focus on the globins of Antarctic notothenioids, emphasizing their adaptive features as inferred from comparisons with human proteins. We dedicate this review to Guido di Prisco, our co-author, colleague, friend, and husband of C.V. Ever thoughtful, creative, and enthusiastic, Guido spearheaded study of the structure, function, and evolution of the hemoglobins of polar fishes - this review is testimony to his wide-ranging contributions. Throughout his career, Guido inspired younger scientists to embrace polar biological research, and he challenged researchers of all ages to explore evolutionary adaptation in the context of global climate change. Beyond his scientific contributions, we will miss his warmth, his culture, and his great intellect. Guido has left an outstanding legacy, one that will continue to inspire us and our research.

Evolution of nodal and nodal-related genes and the putative composition of the heterodimers that trigger the nodal pathway in vertebrates

Nodal is a signaling molecule that belongs to the transforming growth factor-β superfamily that plays key roles during the early stages of development of animals. In vertebrates Nodal forms an heterodimer with a GDF1/3 protein to activate the Nodal pathway. Vertebrates have a paralog of nodal in their genomes labeled Nodal-related, but the evolutionary history of these genes is a matter of debate, mainly because of the presence of a variable numbers of genes in the vertebrate genomes sequenced so far. Thus, the goal of this study was to investigate the evolutionary history of the Nodal and Nodal-related genes with an emphasis in tracking changes in the number of genes among vertebrates. Our results show the presence of two gene lineages (Nodal and Nodal-related) that can be traced back to the ancestor of jawed vertebrates. These lineages have undergone processes of differential retention and lineage-specific expansions. Our results imply that Nodal and Nodal-related duplicated at the latest in the ancestor of gnathostomes, and they still retain a significant level of functional redundancy. By comparing the evolution of the Nodal/Nodal-related with GDF1/3 gene family, it is possible to infer that there are several types of heterodimers that can trigger the Nodal pathway among vertebrates.

Genomic Insights into the Adaptive Convergent Evolution

Adaptive convergent evolution, which refers to the same or similar phenotypes produced by species from independent lineages under similar selective pressures, has been widely examined for a long time. Accumulating studies on the adaptive convergent evolution have been reported from many different perspectives (cellular, anatomical, morphological, physiological, biochemical, and behavioral). Recent advances in the genomic technologies have demonstrated that adaptive convergence can arise from specific genetic mechanisms in different hierarchies, ranging from the same nucleotide or amino acid substitutions to the biological functions or pathways. Among these genetic mechanisms, the same amino acid changes in protein-coding genes play an important role in adaptive phenotypic convergence. Methods for detecting adaptive convergence at the protein sequence level have been constantly debated and developed. Here, we review recent progress on using genomic approaches to evaluate the genetic mechanisms of adaptive convergent evolution, summarize the research methods for identifying adaptive amino acid convergence, and discuss the future perspectives for researching adaptive convergent evolu-tion.

Globin E is a myoglobin-related, respiratory protein highly expressed in lungfish oocytes

Globins are a classical model system for the studies of protein evolution and function. Recent studies have shown that – besides the well-known haemoglobin and myoglobin – additional globin-types occur in vertebrates that serve different functions. Globin E (GbE) was originally identified as an eye-specific protein of birds that is distantly related to myoglobin. GbE is also present in turtles and the coelacanth but appeared to have been lost in other vertebrates. Here, we show that GbE additionally occurs in lungfish, the closest living relatives of the tetrapods. Each lungfish species harbours multiple (≥5) GbE gene copies. Surprisingly, GbE is exclusively and highly expressed in oocytes, with mRNA levels that exceed that of myoglobin in the heart. Thus, GbE is the first known oocyte-specific globin in vertebrates. No GbE transcripts were found in the ovary or egg transcriptomes of other vertebrates, suggesting a lungfish-specific function. Spectroscopic analysis and kinetic studies of recombinant GbE1 of the South American lungfish Lepidosiren paradoxa revealed a typical pentacoordinate globin with myoglobin-like O₂-binding kinetics, indicating similar functions. Our findings suggest that the multiple copies of GbE evolved to enhance O₂-supply in the developing embryo of lungfish, analogous to the embryonic and fetal haemoglobins of other vertebrates. In evolution, GbE must have changed its expression site from oocytes to eyes, or vice versa.

Genome sequence of walking catfish (Clarias batrachus) provides insights into terrestrial adaptation

Background

Walking catfish (Clarias batrachus) is a freshwater fish capable of air-breathing and locomotion on land. It usually inhabits various low-oxygen habitats, burrows inside the mudflat, and sometimes “walks” to search for suitable environments during summer. It has evolved accessory air-breathing organs for respiring air and corresponding mechanisms to survive in such challenging environments. Thereby, it serves as a great model for understanding adaptations to terrestrial life.

Results

Comparative genomics with channel catfish (Ictalurus punctatus) revealed specific adaptations of C. batrachus in DNA repair, enzyme activator activity, and small GTPase regulator activity. Comparative analysis with 11 non-air-breathing fish species suggested adaptive evolution in gene expression and nitrogenous waste metabolic processes. Further, myoglobin, olfactory receptor related to class A G protein-coupled receptor 1, and sulfotransferase 6b1 genes were found to be expanded in the air-breathing walking catfish genome, with 15, 15, and 12 copies, respectively, compared to non-air-breathing fishes that possess only 1–2 copies of these genes. Additionally, we sequenced and compared the transcriptomes of the gill and the air-breathing organ to characterize the mechanism of aerial respiration involved in elastic fiber formation, oxygen binding and transport, angiogenesis, ion homeostasis and acid-base balance. The hemoglobin genes were expressed dramatically higher in the air-breathing organ than in the gill of walking catfish.

Conclusions

This study provides an important genomic resource for understanding the adaptive mechanisms of walking catfish to terrestrial environments. It is possible that the coupling of enhanced abilities for oxygen storage and oxygen transport through genomic expansion of myoglobin genes and transcriptomic up-regulation of hemoglobin and angiogenesis-related genes are important components of the molecular basis for adaptation of this aquatic species to terrestrial life.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5355-9) contains supplementary material, which is available to authorized users.

Gene Turnover and Diversification of the α- and β-Globin Gene Families in Sauropsid Vertebrates

The genes that encode the α- and β-chain subunits of vertebrate hemoglobin have served as a model system for elucidating general principles of gene family evolution, but little is known about patterns of evolution in amniotes other than mammals and birds. Here, we report a comparative genomic analysis of the α- and β-globin gene clusters in sauropsids (archosaurs and nonavian reptiles). The objectives were to characterize changes in the size and membership composition of the α- and β-globin gene families within and among the major sauropsid lineages, to reconstruct the evolutionary history of the sauropsid α- and β-globin genes, to resolve orthologous relationships, and to reconstruct evolutionary changes in the developmental regulation of gene expression. Our comparisons revealed contrasting patterns of evolution in the unlinked α- and β-globin gene clusters. In the α-globin gene cluster, which has remained in the ancestral chromosomal location, evolutionary changes in gene content are attributable to the differential retention of paralogous gene copies that were present in the common ancestor of tetrapods. In the β-globin gene cluster, which was translocated to a new chromosomal location, evolutionary changes in gene content are attributable to differential gene gains (via lineage-specific duplication events) and gene losses (via lineage-specific deletions and inactivations). Consequently, all major groups of amniotes possess unique repertoires of embryonic and postnatally expressed β-type globin genes that diversified independently in each lineage. These independently derived β-type globins descend from a pair of tandemly linked paralogs in the most recent common ancestor of sauropsids.

Evolution of hemoglobin loci and their regulatory elements

Across the expanse of vertebrate evolution, each species produces multiple forms of hemoglobin in erythroid cells at appropriate times and in the proper amounts. The multiple hemoglobins are encoded in two globin gene clusters in almost all species. One globin gene cluster, linked to the gene NPRL3, is preserved in all vertebrates, including a gene cluster encoding the highly divergent globins from jawless vertebrates. This preservation of synteny may reflect the presence of a powerful enhancer of globin gene expression in the NPRL3 gene. Despite substantial divergence in noncoding DNA sequences among mammals, several epigenetic features of the globin gene regulatory regions are preserved across vertebrates. The preserved features include multiple DNase hypersensitive sites, at least one of which is an enhancer, and binding by key lineage-restricted transcription factors such as GATA1 and TAL1, which in turn recruit coactivators such as P300 that catalyze acetylation of histones. The maps of epigenetic features are strongly correlated with activity in gene regulation, and resources for accessing and visualizing such maps are readily available to the community of researchers and students.

Convergent evolution of caffeine in plants by co-option of exapted ancestral enzymes

Convergent evolution is a process that has occurred throughout the tree of life, but the historical genetic and biochemical context promoting the repeated independent origins of a trait is rarely understood. The well-known stimulant caffeine, and its xanthine alkaloid precursors, has evolved multiple times in flowering plant history for various roles in plant defense and pollination. We have shown that convergent caffeine production, surprisingly, has evolved by two previously unknown biochemical pathways in chocolate, citrus, and guaraná plants using either caffeine synthase- or xanthine methyltransferase-like enzymes. However, the pathway and enzyme lineage used by any given plant species is not predictable from phylogenetic relatedness alone. Ancestral sequence resurrection reveals that this convergence was facilitated by co-option of genes maintained over 100 million y for alternative biochemical roles. The ancient enzymes of the Citrus lineage were exapted for reactions currently used for various steps of caffeine biosynthesis and required very few mutations to acquire modern-day enzymatic characteristics, allowing for the evolution of a complete pathway. Future studies aimed at manipulating caffeine content of plants will require the use of different approaches given the metabolic and genetic diversity revealed by this study.

Functional diversification of sea lamprey globins in evolution and development

Agnathans have a globin repertoire that markedly differs from that of jawed (gnathostome) vertebrates. The sea lamprey (Petromyzon marinus) harbors at least 18 hemoglobin, two myoglobin, two globin X, and one cytoglobin genes. However, agnathan hemoglobins and myoglobins are not orthologous to their cognates in jawed vertebrates. Thus, blood-based O₂ transport and muscle-based O₂ storage proteins emerged twice in vertebrates from a tissue-globin ancestor. Notably, the sea lamprey displays three switches in hemoglobin expression in its life cycle, analogous to hemoglobin switching in vertebrates. To study the functional changes associated with the evolution and ontogenesis of distinct globin types, we determined O₂ binding equilibria, type of quaternary assembly, and nitrite reductase enzymatic activities of one adult (aHb5a) and one embryonic/larval hemoglobin (aHb6), myoglobin (aMb1) and cytoglobin (Cygb) of the sea lamprey. We found clear functional differentiation among globin types expressed at different developmental stages and in different tissues. Cygb and aMb1 have high O₂ affinity and nitrite reductase activity, while the two hemoglobins display low O₂ affinity and nitrite reductase activity. Cygb and aHb6 but not aHb5a show cooperative O₂ binding, correlating with increased stability of dimers, as shown by gel filtration and molecular modeling. The high O₂-affinity and the lack of cooperativity confirm the identity of the sea lamprey aMb1 as O₂ storage protein of the muscle. The dimeric structure and O₂-binding properties of sea lamprey and mammalian Cygb were very similar, suggesting a conservation of function since their divergence around 500million years ago.

Gene Duplication and Evolutionary Innovations in Hemoglobin-Oxygen Transport

During vertebrate evolution, duplicated hemoglobin (Hb) genes diverged with respect to functional properties as well as the developmental timing of expression. For example, the subfamilies of genes that encode the different subunit chains of Hb are ontogenetically regulated such that functionally distinct Hb isoforms are expressed during different developmental stages. In some vertebrate taxa, functional differentiation between co-expressed Hb isoforms may also contribute to physiologically important divisions of labor.

CRH peptide evolution occurred in three phases: Evidence from characterizing sea lamprey CRH system members

The corticotropin releasing hormone (CRH) system, which includes the CRH family of peptides, their receptors (CRHRs) and a binding protein (CRHBP), has been strongly conserved throughout vertebrate evolution. The identification of invertebrate homologues suggests this system evolved over 500 million years ago. However, the early vertebrate evolution of the CRH system is not understood. Current theory indicates that agnathans (hagfishes and lampreys) are monophyletic with a conservative evolution over the past 500million years and occupy a position at the root of vertebrate phylogeny. We isolated the cDNAs for three CRH family members, two CRHRs and a CRHBP from the sea lamprey, Petromyzon marinus. Two of the CRH peptides are related to the CRH/urotensin-1 (UI) lineage, whereas the other is a urocortin (Ucn) 3 orthologue. The predicted amino acid identity of CRH and UI is 61% but they possess distinct motifs indicative of each peptide, suggesting an early divergence of the two genes. Based on our findings we propose the CRH peptides evolved in at least 3 distinct phases. The first occurring prior to the agnathans gave rise to the CRH/UI-like and Ucn2/3-like paralogous lineages. The second was a partial sub-genomic duplication of the ancestral CRH/UI-like gene, but not the Ucn2/3-like gene, giving rise to the CRH and UI (Ucn) lineages. The third event which resulted in the appearance of Ucn2 and Ucn3 must have occurred after the evolution of the cartilaginous fishes. Interestingly, unlike other vertebrate CRHRs, we were unable to classify our two P. marinus receptors (designated CRHRα and CRHRβ) as either type 1 or type 2, indicating that this split evolved later in vertebrate evolution. A single CRHBP gene was found suggesting that either this gene has not been affected by the vertebrate genome duplications or there have been a series of paralogous gene deletions. This study suggests that P. marinus possess a functional CRH system that differs from that of the gnathostomes and may represent a model for the earliest functioning CRH system in vertebrates.

Insight into the phylogenetic relationship and structural features of vertebrate myoglobin family

Myoglobin (Mb) is studied to clarify the structure-function relationships in protein science. In this work, we report the results of a comparative analysis of amino acid sequences from 298 vertebrate Mbs. Forty-one high conserved residues were identified and seven of them were invariants [E18, G25, F43, V68, L72, H93 (proximal histidine) and H97]. E18 is the only invariant amino acid residue located out of the heme-pocket and Xe-cavities playing a role in interaction between the A and E-helices. A comparative analysis of several parameters related to amino acid composition shows an increase of average mass, accessible surface area and volume per residue from Actinopterygii to Mammalia and Aves. This may be due to an increased number of bulky residues reducing the non-specific cavities volume and thus improving the oxygen flow between the heme site and the outside of the protein. Finally, the phylogenetic analyses of Mb in vertebrates are consistent with an evolution that runs with the diversification of the species, but in which several episodes of gene duplication and lost have occurred, less frequently in the ancestors of great taxons, cartilaginous fishes and non-avian reptiles, most frequently in ray-finned fishes and mammals, and very frequently in birds.

Distribution of Cytoglobin in the Mouse Brain

Cytoglobin (Cygb) is a vertebrate globin with so far poorly defined function. It is expressed in the fibroblast cell-lineage but has also been found in neurons. Here we provide, using immunohistochemistry, a detailed study on the distribution of Cygb in the mouse brain. While Cygb is a cytoplasmic protein in active cells of the supportive tissue, in neurons it is located in the cytoplasm and the nucleus. We found the expression of Cygb in all brain regions, although only a fraction of the neurons was Cygb-positive. Signals were of different intensity ranging from faint to very intense. Telencephalic neurons in all laminae of the cerebral cortex (CCo), in the olfactory bulb (in particular periglomerular cells), in the hippocampal formation (strongly stained pyramidal cells with long processes), basal ganglia (scattered multipolar neurons in the dorsal striatum, dorsal and ventral pallidum (VP)), and in the amygdala (neurons with unlabeled processes) were labeled by the antibody. In the diencephalon, we observed Cygb-positive neurons of moderate intensity in various nuclei of the dorsal thalamus, in the hypothalamus, metathalamus (geniculate nuclei), epithalamus with strong labeling of habenular nucleus neurons and no labeling of pineal cells, and in the ventral thalamus. Tegmental neurons stood out by strongly stained somata with long processes in, e.g., the laterodorsal nucleus. In the tectum, faintly labeled neurons and fibers were detected in the superior colliculus (SC). The cerebellum exhibited unlabeled Purkinje-neurons but signs of strong afferent cortical innervation. Neurons in the gray matter of the spinal cord showed moderate immunofluorescence. Peripheral ganglia were not labeled by the antibody. The Meynert-fascicle and the olfactory and optic nerves/tracts were the only Cygb-immunoreactive (Cygb-IR) fiber systems. Notably, we found a remarkable level of colocalization of Cygb and neuronal nitric oxide (NO)-synthase in neurons, which supports a functional association.

Convergent evolution of hemoglobin switching in jawed and jawless vertebrates

Background

During development, humans and other jawed vertebrates (Gnathostomata) express distinct hemoglobin genes, resulting in different hemoglobin tetramers. Embryonic and fetal hemoglobin have higher oxygen affinities than the adult hemoglobin, sustaining the oxygen demand of the developing organism. Little is known about the expression of hemoglobins during development of jawless vertebrates (Agnatha).

Results

We identified three hemoglobin switches in the life cycle of the sea lamprey. Three hemoglobin genes are specifically expressed in the embryo, four genes in the filter feeding larva (ammocoete), and nine genes correspond to the adult hemoglobin chains. During the development from the parasitic to the reproductive adult, the composition of hemoglobin changes again, with a massive increase of chain aHb1. A single hemoglobin chain is expressed constitutively in all stages. We further showed the differential expression of other globin genes: Myoglobin 1 is most highly expressed in the reproductive adult, myoglobin 2 expression peaks in the larva. Globin X1 is restricted to the embryo; globin X2 was only found in the reproductive adult. Cytoglobin is expressed at low levels throughout the life cycle.

Conclusion

Because the hemoglobins of jawed and jawless vertebrates evolved independently from a common globin ancestor, hemoglobin switching must also have evolved convergently in these taxa. Notably, the ontogeny of sea lamprey hemoglobins essentially recapitulates their phylogeny, with the embryonic hemoglobins emerging first, followed by the evolution of larval and adult hemoglobins.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0597-0) contains supplementary material, which is available to authorized users.

Early Evolution of Vertebrate Mybs: An Integrative Perspective Combining Synteny, Phylogenetic, and Gene Expression Analyses

The genes in the Myb superfamily encode for three related transcription factors in most vertebrates, A-, B-, and c-Myb, with functionally distinct roles, whereas most invertebrates have a single Myb. B-Myb plays an essential role in cell division and cell cycle progression, c-Myb is involved in hematopoiesis, and A-Myb is involved in spermatogenesis and regulating expression of pachytene PIWI interacting RNAs, a class of small RNAs involved in posttranscriptional gene regulation and the maintenance of reproductive tissues. Comparisons between teleost fish and tetrapods suggest that the emergence and functional divergence of the Myb genes were linked to the two rounds of whole-genome duplication early in vertebrate evolution. We combined phylogenetic, synteny, structural, and gene expression analyses of the Myb paralogs from elephant shark and lampreys with data from 12 bony vertebrates to reconstruct the early evolution of vertebrate Mybs. Phylogenetic and synteny analyses suggest that the elephant shark and Japanese lamprey have copies of the A-, B-, and c-Myb genes, implying their origin could be traced back to the common ancestor of lampreys and gnathostomes. However, structural and gene expression analyses suggest that their functional roles diverged between gnathostomes and cyclostomes. In particular, we did not detect A-Myb expression in testis suggesting that the involvement of A-Myb in the pachytene PIWI interacting RNA pathway is probably a gnathostome-specific innovation. We speculate that the secondary loss of a central domain in lamprey A-Myb underlies the functional differences between the cyclostome and gnathostome A-Myb proteins.

The Full Globin Repertoire of Turtles Provides Insights into Vertebrate Globin Evolution and Functions

Globins are small heme proteins that play an important role in oxygen supply, but may also have other functions. Globins offer a unique opportunity to study the functional evolution of genes and proteins. We have characterized the globin repertoire of two different turtle species: the Chinese softshell turtle (Pelodiscus sinensis) and the western painted turtle (Chrysemys picta bellii). In the genomes of both species, we have identified eight distinct globin types: hemoglobin (Hb), myoglobin, neuroglobin, cytoglobin, globin E, globin X, globin Y, and androglobin. Therefore, along with the coelacanth, turtles are so far the only known vertebrates with a full globin repertoire. This fact allows for the first time a comparative analysis of the expression of all eight globins in a single species. Phylogenetic analysis showed an early divergence of neuroglobin and globin X before the radiation of vertebrates. Among the other globins, cytoglobin diverged first, and there is a close relationship between myoglobin and globin E; the position of globin Y is not resolved. The globin E gene was selectively lost in the green anole, and the genes coding for globin X and globin Y were deleted in chicken. Quantitative real-time reverse transcription polymerase chain reaction experiments revealed that myoglobin, neuroglobin, and globin E are highly expressed with tissue-specific patterns, which are in line with their roles in the oxidative metabolism of the striated muscles, the brain, and the retina, respectively. Histochemical analyses showed high levels of globin E in the pigment epithelium of the eye. Globin E probably has a myoglobin-like role in transporting O₂ across the pigment epithelium to supply in the metabolically highly active retina.

Gene turnover in the avian globin gene families and evolutionary changes in hemoglobin isoform expression

The apparent stasis in the evolution of avian chromosomes suggests that birds may have experienced relatively low rates of gene gain and loss in multigene families. To investigate this possibility and to explore the phenotypic consequences of variation in gene copy number, we examined evolutionary changes in the families of genes that encode the α- and β-type subunits of hemoglobin (Hb), the tetrameric α2β2 protein responsible for blood-O2 transport. A comparative genomic analysis of 52 bird species revealed that the size and membership composition of the α- and β-globin gene families have remained remarkably constant during approximately 100 My of avian evolution. Most interspecific variation in gene content is attributable to multiple independent inactivations of the α(D)-globin gene, which encodes the α-chain subunit of a functionally distinct Hb isoform (HbD) that is expressed in both embryonic and definitive erythrocytes. Due to consistent differences in O2-binding properties between HbD and the major adult-expressed Hb isoform, HbA (which incorporates products of the α(A)-globin gene), recurrent losses of α(D)-globin contribute to among-species variation in blood-O2 affinity. Analysis of HbA/HbD expression levels in the red blood cells of 122 bird species revealed high variability among lineages and strong phylogenetic signal. In comparison with the homologous gene clusters in mammals, the low retention rate for lineage-specific gene duplicates in the avian globin gene clusters suggests that the developmental regulation of Hb synthesis in birds may be more highly conserved, with orthologous genes having similar stage-specific expression profiles and similar functional properties in disparate taxa.

Ancient Duplications and Expression Divergence in the Globin Gene Superfamily of Vertebrates: Insights from the Elephant Shark Genome and Transcriptome

Comparative analyses of vertebrate genomes continue to uncover a surprising diversity of genes in the globin gene superfamily, some of which have very restricted phyletic distributions despite their antiquity. Genomic analysis of the globin gene repertoire of cartilaginous fish (Chondrichthyes) should be especially informative about the duplicative origins and ancestral functions of vertebrate globins, as divergence between Chondrichthyes and bony vertebrates represents the most basal split within the jawed vertebrates. Here, we report a comparative genomic analysis of the vertebrate globin gene family that includes the complete globin gene repertoire of the elephant shark (Callorhinchus milii). Using genomic sequence data from representatives of all major vertebrate classes, integrated analyses of conserved synteny and phylogenetic relationships revealed that the last common ancestor of vertebrates possessed a repertoire of at least seven globin genes: single copies of androglobin and neuroglobin, four paralogous copies of globin X, and the single-copy progenitor of the entire set of vertebrate-specific globins. Combined with expression data, the genomic inventory of elephant shark globins yielded four especially surprising findings: 1) there is no trace of the neuroglobin gene (a highly conserved gene that is present in all other jawed vertebrates that have been examined to date), 2) myoglobin is highly expressed in heart, but not in skeletal muscle (reflecting a possible ancestral condition in vertebrates with single-circuit circulatory systems), 3) elephant shark possesses two highly divergent globin X paralogs, one of which is preferentially expressed in gonads, and 4) elephant shark possesses two structurally distinct α-globin paralogs, one of which is preferentially expressed in the brain. Expression profiles of elephant shark globin genes reveal distinct specializations of function relative to orthologs in bony vertebrates and suggest hypotheses about ancestral functions of vertebrate globins.