Gene duplication, genome duplication, and the functional diversification of vertebrate globins

Haemoglobin Gene Repertoire in Teleost and Cichlid Fishes Shaped by Gene Duplications and Genome Rearrangements

Haemoglobin is a key molecule for oxygen transport in vertebrates. It exhibits remarkable gene diversity in teleost fishes, reflecting adaptation to various aquatic environments. In this study, we present the dynamic evolution of haemoglobin subunit genes based on a comparison of high-quality genome assemblies of 24 vertebrate species, including 17 teleosts (of which six are cichlids). Our findings indicate that teleost genomes contain a range of haemoglobin genes, from as few as five in fugu to as many as 43 in salmon, with the latter being the largest repertoire found in vertebrates. We find evidence that the teleost ancestor had at least four Hbα and three or four Hbβ subunit genes, and that the current gene diversity emerged during teleost radiation, driven primarily by (tandem) gene duplications, genome compaction, and rearrangement dynamics. We provide insights into the genomic organisation of haemoglobin clusters in different teleost species. We further show that the evolution of paralogous rhbdf1 genes flanking both teleost clusters (LA and MN) supports the hypothesis for the origin of the LA cluster by rearrangement within teleosts, rather than by the teleost specific whole-genome duplication. We specifically focus on cichlid fishes, where adaptation to low oxygen environment plays role in species diversification. Our analysis of six cichlid genomes, including Pungu maclareni from the Barombi Mbo crater lake, for which we sequenced a representative genome, reveals 18-32 copies of the Hb genes, and elevated rates of non-synonymous substitutions compared to other teleosts. Overall, this work facilitates a deeper understanding of how haemoglobin genes contribute to the adaptive potential of teleosts.

Selection and Gene Duplication Associated With High-Elevation Diversification in Pristimantis, the Largest Terrestrial Vertebrate Genus

The genus Pristimantis diversified in the tropical Andes mountains and is the most speciose genus of terrestrial vertebrates. Pristimantis are notable among frogs in that they thrive at high elevations (>2,000 m) and are direct developers without a tadpole stage. Despite their ecological significance, little is known about the genetic and physiological traits enabling their success. We conducted transcriptomic analysis on seven Pristimantis species sampled across elevations in the Ecuadorean Andes to explore three hypotheses for their success: (i) unique genes are under selection relative to all other frogs, (ii) common selection occurs across all direct developers, or (iii) common selection occurs across all high-elevation frog clades. Comparative analysis with 34 frog species revealed unique positive selection in Pristimantis genes related to aerobic respiration, hemostasis, signaling, cellular transportation of proteins and ions, and immunity. Additionally, we detected positive selection across all direct developers for genes associated with oxygenase activity and metal ion binding. While many genes under selection in Pristimantis were not positively selected in other high-elevation frog species, we identified some shared genes and pathways linked to lipid metabolism, innate immunity, and cellular redox processes. We observed more positive selection in duplicated- versus single-copy genes, while relaxed purifying selection was prevalent in single-copy genes. Notably, copy number of an innate immunity complement gene was positively correlated with Pristimantis species elevation. Our findings contribute novel insights into the genetic basis of adaptation in Pristimantis and provide a foundation for future studies on the evolutionary mechanisms leading to direct development and coping with high elevations.

The unique allosteric property of crocodilian haemoglobin elucidated by cryo-EM

The principal effect controlling the oxygen affinity of vertebrate haemoglobins (Hbs) is the allosteric switch between R and T forms with relatively high and low oxygen affinity respectively. Uniquely among jawed vertebrates, crocodilians possess Hb that shows a profound drop in oxygen affinity in the presence of bicarbonate ions. This allows them to stay underwater for extended periods by consuming almost all the oxygen present in the blood-stream, as metabolism releases carbon dioxide, whose conversion to bicarbonate and hydrogen ions is catalysed by carbonic anhydrase. Despite the apparent universal utility of bicarbonate as an allosteric regulator of Hb, this property evolved only in crocodilians. We report here the molecular structures of both human and a crocodilian Hb in the deoxy and liganded states, solved by cryo-electron microscopy. We reveal the precise interactions between two bicarbonate ions and the crocodilian protein at symmetry-related sites found only in the T state. No other known effector of vertebrate Hbs binds anywhere near these sites.

Globin phylogeny, evolution and function, the newest update

Our globin census update allows us to refine our vision of globin origin, evolution, and structure to function relationship in the context of the currently accepted tree of life. The modern globin domain originates as a single domain, three-over-three α-helical folded structure before the diversification of the kingdoms of life (Bacteria, Archaea, Eukarya). Together with the diversification of prokaryotes, three monophyletic globin families (M, S, and T) emerged, most likely in Proteobacteria and Actinobacteria, displaying specific sequence and structural features, and spread by vertical and horizontal gene transfer, most probably already present in the last universal common ancestor (LUCA). Non-globin domains were added, and eventually lost again, creating multi-domain structures in key branches of M- (FHb and Adgb) and the vast majority of S globins, which with their coevolved multi-domain architectures, have predominantly "sensor" functions. Single domain T-family globins diverged into four major groups and most likely display functions related to reactive nitrogen and oxygen species (RNOS) chemistry, as well as oxygen storage/transport which drives the evolution of its major branches with their characteristic key distal residues (B10, E11, E7, and G8). M-family evolution also lead to distinctive major types (FHb and Fgb, Ngb, Adgb, GbX vertebrate Gbs), and shows the shift from high oxygen affinity controlled by TyrB10-Gln/AsnE11 likely related to RNOS chemistry in microorganisms, to a moderate oxygen affinity storage/transport function controlled by hydrophobic B10/E11-HisE7 in multicellular animals.

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.

Different structural variant prediction tools yield considerably different results in Caenorhabditis elegans

The accurate characterization of structural variation is crucial for our understanding of how large chromosomal alterations affect phenotypic differences and contribute to genome evolution. Whole-genome sequencing is a popular approach for identifying structural variants, but the accuracy of popular tools remains unclear due to the limitations of existing benchmarks. Moreover, the performance of these tools for predicting variants in non-human genomes is less certain, as most tools were developed and benchmarked using data from the human genome. To evaluate the use of long-read data for the validation of short-read structural variant calls, the agreement between predictions from a short-read ensemble learning method and long-read tools were compared using real and simulated data from Caenorhabditis elegans. The results obtained from simulated data indicate that the best performing tool is contingent on the type and size of the variant, as well as the sequencing depth of coverage. These results also highlight the need for reference datasets generated from real data that can be used as 'ground truth' in benchmarks.

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.

A Rag GTPase dimer code defines the regulation of mTORC1 by amino acids

Amino acid availability controls mTORC1 activity via a heterodimeric Rag GTPase complex that functions as a scaffold at the lysosomal surface, bringing together mTORC1 with its activators and effectors. Mammalian cells express four Rag proteins (RagA–D) that form dimers composed of RagA/B bound to RagC/D. Traditionally, the Rag paralogue pairs (RagA/B and RagC/D) are referred to as functionally redundant, with the four dimer combinations used interchangeably in most studies. Here, by using genetically modified cell lines that express single Rag heterodimers, we uncover a Rag dimer code that determines how amino acids regulate mTORC1. First, RagC/D differentially define the substrate specificity downstream of mTORC1, with RagD promoting phosphorylation of its lysosomal substrates TFEB/TFE3, while both Rags are involved in the phosphorylation of non-lysosomal substrates such as S6K. Mechanistically, RagD recruits mTORC1 more potently to lysosomes through increased affinity to the anchoring LAMTOR complex. Furthermore, RagA/B specify the signalling response to amino acid removal, with RagB-expressing cells maintaining lysosomal and active mTORC1 even upon starvation. Overall, our findings reveal key qualitative differences between Rag paralogues in the regulation of mTORC1, and underscore Rag gene duplication and diversification as a potentially impactful event in mammalian evolution.

Gollwitzer, Grützmacher et al. and Figlia et al. establish that the various Rag GTPase genes and isoforms differentially regulate mTORC1 activity and distinctly modulate the responsiveness of mammalian cells to amino acid availability.

Segregation of α- and β-Globin Gene Cluster in Vertebrate Evolution: Chance or Necessity?

The review is devoted to the patterns of evolution of α- and β-globin gene domains. A hypothesis is presented according to which segregation of the ancestral cluster of α/β-globin genes in Amniota occurred due to the performance by α-globins and β-globins of non-canonical functions not related to oxygen transport.

Comparative genome analysis of the monogonont marine rotifer Brachionus manjavacas Australian strain: Potential application for ecotoxicology and environmental genomics

This is the first study to analyze the whole-genome sequence of B. manjavacas Australian (Aus.) strain through combination of Oxford Nanopore long-read seq, resulting in a total length of 108.1 Mb and 75 contigs. Genome-wide detoxification related gene families in B. manjavacas Aus. strain were comparatively analyzed with those previously identified in other Brachionus spp., including B. manjavacas German (Ger.) strain. Most of the subfamilies in detoxification related families (CYPs, GSTs, and ABCs) were highly conserved and confirmed orthologous relationship with Brachionus spp., and with accumulation of genome data, clear differences between genomic repertoires were demonstrated the marine and the freshwater species. Furthermore, strain-specific genetic variations were present between the Aus. and Ger. strains of B. manjavacas. This whole-genome analysis provides in-depth review on the genomic structural differences for detoxification-related gene families and further provides useful information for comparative ecotoxicological studies and evolution of detoxification mechanisms in Brachionus spp.

Evolutionary Adaptations of Parasitic Flatworms to Different Oxygen Tensions

During the evolution of the Earth, the increase in the atmospheric concentration of oxygen gave rise to the development of organisms with aerobic metabolism, which utilized this molecule as the ultimate electron acceptor, whereas other organisms maintained an anaerobic metabolism. Platyhelminthes exhibit both aerobic and anaerobic metabolism depending on the availability of oxygen in their environment and/or due to differential oxygen tensions during certain stages of their life cycle. As these organisms do not have a circulatory system, gas exchange occurs by the passive diffusion through their body wall. Consequently, the flatworms developed several adaptations related to the oxygen gradient that is established between the aerobic tegument and the cellular parenchyma that is mostly anaerobic. Because of the aerobic metabolism, hydrogen peroxide (H₂O₂) is produced in abundance. Catalase usually scavenges H₂O₂ in mammals; however, this enzyme is absent in parasitic platyhelminths. Thus, the architecture of the antioxidant systems is different, depending primarily on the superoxide dismutase, glutathione peroxidase, and peroxiredoxin enzymes represented mainly in the tegument. Here, we discuss the adaptations that parasitic flatworms have developed to be able to transit from the different metabolic conditions to those they are exposed to during their life cycle.

Changes in hemoglobin function and isoform expression during embryonic development in the American alligator, Alligator mississippiensis

In the developing embryos of egg-laying vertebrates, O₂ flux takes place across a fixed surface area of the eggshell and the chorioallantoic membrane. In the case of crocodilians, the developing embryo may experience a decrease in O₂ flux when the nest becomes hypoxic, which may cause compensatory adjustments in blood O₂ transport. However, whether the switch from embryonic to adult hemoglobin isoforms (isoHbs) plays some role in these adjustments is unknown. Here, we provide a detailed characterization of the developmental switch of isoHb synthesis in the American alligator, Alligator mississippiensis. We examined the in vitro functional properties and subunit composition of purified alligator isoHbs expressed during embryonic developmental stages in normoxia and hypoxia (10% O₂). We found distinct patterns of isoHb expression in alligator embryos at different stages of development, but these patterns were not affected by hypoxia. Specifically, alligator embryos expressed two main isoHbs: HbI, prevalent at early developmental stages, with a high O₂ affinity and high ATP sensitivity, and HbII, prevalent at later stages and identical to the adult protein, with a low O₂ affinity and high CO₂ sensitivity. These results indicate that whole blood O₂ affinity is mainly regulated by ATP in the early embryo and by CO₂ and bicarbonate from the late embryo until adult life, but the developmental regulation of isoHb expression is not affected by hypoxia exposure.

Lineage-Specific Genes and Family Expansions in Dictyostelid Genomes Display Expression Bias and Evolutionary Diversification during Development

Gene duplications generate new genes that can contribute to expression changes and the evolution of new functions. Genomes often consist of gene families that undergo expansions, some of which occur in specific lineages that reflect recent adaptive diversification. In this study, lineage-specific genes and gene family expansions were studied across five dictyostelid species to determine when and how they are expressed during multicellular development. Lineage-specific genes were found to be enriched among genes with biased expression (predominant expression in one developmental stage) in each species and at most developmental time points, suggesting independent functional innovations of new genes throughout the phylogeny. Biased duplicate genes had greater expression divergence than their orthologs and paralogs, consistent with subfunctionalization or neofunctionalization. Lineage-specific expansions in particular had biased genes with both molecular signals of positive selection and high expression, suggesting adaptive genetic and transcriptional diversification following duplication. Our results present insights into the potential contributions of lineage-specific genes and families in generating species-specific phenotypes during multicellular development in dictyostelids.

Whole-Genome Duplications and the Diversification of the Globin-X Genes of Vertebrates

Globin-X (GbX) is an enigmatic member of the vertebrate globin gene family with a wide phyletic distribution that spans protostomes and deuterostomes. Unlike canonical globins such as hemoglobins and myoglobins, functional data suggest that GbX does not have a primary respiratory function. Instead, evidence suggests that the monomeric, membrane-bound GbX may play a role in cellular signaling or protection against the oxidation of membrane lipids. Recently released genomes from key vertebrates provide an excellent opportunity to address questions about the early stages of the evolution of GbX in vertebrates. We integrate bioinformatics, synteny, and phylogenetic analyses to characterize the diversity of GbX genes in nonteleost ray-finned fishes, resolve relationships between the GbX genes of cartilaginous fish and bony vertebrates, and demonstrate that the GbX genes of cyclostomes and gnathostomes derive from independent duplications. Our study highlights the role that whole-genome duplications (WGDs) have played in expanding the repertoire of genes in vertebrate genomes. Our results indicate that GbX paralogs have a remarkably high rate of retention following WGDs relative to other globin genes and provide an evolutionary framework for interpreting results of experiments that examine functional properties of GbX and patterns of tissue-specific expression. By identifying GbX paralogs that are products of different WGDs, our results can guide the design of experimental work to explore whether gene duplicates that originate via WGDs have evolved novel functional properties or expression profiles relative to singleton or tandemly duplicated copies of GbX.

Hypoxic and nitrosative stress conditions modulate expression of myoglobin genes in a carcinogenic hepatobiliary trematode, Clonorchis sinensis

Despite recent evidence suggesting that adult trematodes require oxygen for the generation of bioenergy and eggshells, information on the molecular mechanism by which the parasites acquire oxygen remains largely elusive. In this study, the structural and expressional features of globin genes identified in Clonorchis sinensis, a carcinogenic trematode parasite that invades the hypoxic biliary tracts of mammalian hosts, were investigated to gain insight into the molecules that enable oxygen metabolism. The number of globin paralogs substantially differed among parasitic platyhelminths, ranging from one to five genes, and the C. sinensis genome encoded at least five globin genes. The expression of these Clonorchis genes, named CsMb (CsMb1—CsMb3), CsNgb, and CsGbX, according to their preferential similarity patterns toward respective globin subfamilies, exponentially increased in the worms coinciding with their sexual maturation, after being downregulated in early juveniles compared to those in metacercariae. The CsMb1 protein was detected throughout the parenchymal region of adult worms as well as in excretory-secretory products, whereas the other proteins were localized exclusively in the sexual organs and intrauterine eggs. Stimuli generated by exogenous oxygen, nitric oxide (NO), and nitrite as well as co-incubation with human cholangiocytes variously affected globin gene expression in live C. sinensis adults. Together with the specific histological distributions, these hypoxia-induced patterns may suggest that oxygen molecules transported by CsMb1 from host environments are provided to cells in the parenchyma and intrauterine eggs/sex organs of the worms for energy metabolism and/or, more importantly, eggshell formation by CsMb1 and CsMb3, respectively. Other globin homologs are likely to perform non-respiratory functions. Based on the responsive expression profile against nitrosative stress, an oxygenated form of secreted CsMb1 is suggested to play a pivotal role in parasite survival by scavenging NO generated by host immune cells via its NO dioxygenase activity.

Trematode parasites that invade mammalian tissues have long been believed to produce bioenergy via anaerobic respiration in their definitive hosts. However, recent studies have revealed that these parasites require considerable amounts of oxygen for the generation of hard eggshells during sexual reproduction as well as energy metabolism. Despite these findings, information on the biological mechanisms and relevant molecules responsible for oxygen uptake in the host environment remains largely elusive. Clonorchis sinensis is a carcinogenic trematode parasite that causes clonorchiasis in humans by infecting the bile ducts. Here, we investigated globin genes/proteins in the liver fluke. The genome of C. sinensis encoded at least five globin paralogs (CsMb1, CsMb2, CsMb3, CsNgb, and CsGbX). Temporal expression of these globin genes coincided with the sexual maturation of C. sinensis. Based on the histological localities and induction profiles upon hypoxia, it could be postulated that the oxygen molecules transported by CsMb1 from host environments are provided to cells in the parenchyma and intrauterine eggs/sex organs of the worms by CsMb1 and CsMb3, respectively, for energy metabolism and eggshell formation. Other globin homologs were likely to perform non-respiratory functions. In addition, the oxygenated form of secreted CsMb1 seemed to participate in the scavenging of nitric oxide generated by host immune cells via its nitric oxide dioxygenase activity to increase the survival of the parasite.

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.

Co-Regulated Genes and Gene Clusters

There are many co-regulated genes in eukaryotic cells. The coordinated activation or repression of such genes occurs at specific stages of differentiation, or under the influence of external stimuli. As a rule, co-regulated genes are dispersed in the genome. However, there are also gene clusters, which contain paralogous genes that encode proteins with similar functions. In this aspect, they differ significantly from bacterial operons containing functionally linked genes that are not paralogs. In this review, we discuss the reasons for the existence of gene clusters in vertebrate cells and propose that clustering is necessary to ensure the possibility of selective activation of one of several similar genes.

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.

Recent genome duplications facilitate the phenotypic diversity of Hb repertoire in the Cyprinidae

Whole-genome duplications (WGDs) are an important contributor to phenotypic innovations in evolutionary history. The diversity of blood oxygen transport traits is the perfect reflection of physiological versatility for evolutionary success among vertebrates. In this study, the evolutionary changes of hemoglobin (Hb) repertoire driven by the recent genome duplications were detected in representative Cyprinidae fish, including eight diploid and four tetraploid species. Comparative genomic analysis revealed a substantial variation in both membership composition and intragenomic organization of Hb genes in these species. Phylogenetic reconstruction analyses were conducted to characterize the evolutionary history of these genes. Data were integrated with the expression profiles of the genes during ontogeny. Our results indicated that genome duplications facilitated the phenotypic diversity of the Hb gene family; each was associated with species-specific changes in gene content via gene loss and fusion after genome duplications. This led to repeated evolutionary transitions in the ontogenic regulation of Hb gene expression. Our results revealed that genome duplications helped to generate phenotypic changes in Cyprinidae Hb systems.

Evolutionary History of the Globin Gene Family in Annelids

Animals depend on the sequential oxidation of organic molecules to survive; thus, oxygen-carrying/transporting proteins play a fundamental role in aerobic metabolism. Globins are the most common and widespread group of respiratory proteins. They can be divided into three types: circulating intracellular, noncirculating intracellular, and extracellular, all of which have been reported in annelids. The diversity of oxygen transport proteins has been underestimated across metazoans. We probed 250 annelid transcriptomes in search of globin diversity in order to elucidate the evolutionary history of this gene family within this phylum. We report two new globin types in annelids, namely androglobins and cytoglobins. Although cytoglobins and myoglobins from vertebrates and from invertebrates are referred to by the same name, our data show they are not genuine orthologs. Our phylogenetic analyses show that extracellular globins from annelids are more closely related to extracellular globins from other metazoans than to the intracellular globins of annelids. Broadly, our findings indicate that multiple gene duplication and neo-functionalization events shaped the evolutionary history of the globin family.

Phylogenic Determinants of Cardiovascular Frailty, Focus on Hemodynamics and Arterial Smooth Muscle Cells

The evolution of the circulatory system from invertebrates to mammals has involved the passage from an open system to a closed in-parallel system via a closed in-series system, accompanying the increasing complexity and efficiency of life’s biological functions. The archaic heart enables pulsatile motion waves of hemolymph in invertebrates, and the in-series circulation in fish occurs with only an endothelium, whereas mural smooth muscle cells appear later. The present review focuses on evolution of the circulatory system. In particular, we address how and why this evolution took place from a closed, flowing, longitudinal conductance at low pressure to a flowing, highly pressurized and bifurcating arterial compartment. However, although arterial pressure was the latest acquired hemodynamic variable, the general teleonomy of the evolution of species is the differentiation of individual organ function, supported by specific fueling allowing and favoring partial metabolic autonomy. This was achieved via the establishment of an active contractile tone in resistance arteries, which permitted the regulation of blood supply to specific organ activities via its localized function-dependent inhibition (active vasodilation). The global resistance to viscous blood flow is the peripheral increase in frictional forces caused by the tonic change in arterial and arteriolar radius, which backscatter as systemic arterial blood pressure. Consequently, the arterial pressure gradient from circulating blood to the adventitial interstitium generates the unidirectional outward radial advective conductance of plasma solutes across the wall of conductance arteries. This hemodynamic evolution was accompanied by important changes in arterial wall structure, supported by smooth muscle cell functional plasticity, including contractility, matrix synthesis and proliferation, endocytosis and phagocytosis, etc. These adaptive phenotypic shifts are due to epigenetic regulation, mainly related to mechanotransduction. These paradigms actively participate in cardio-arterial pathologies such as atheroma, valve disease, heart failure, aneurysms, hypertension, and physiological aging.

Structural Phylogenetics with Confidence

For evaluating the deepest evolutionary relationships among proteins, sequence similarity is too low for application of sequence-based homology search or phylogenetic methods. In such cases, comparison of protein structures, which are often better conserved than sequences, may provide an alternative means of uncovering deep evolutionary signal. Although major protein structure databases such as SCOP and CATH hierarchically group protein structures, they do not describe the specific evolutionary relationships within a hierarchical level. Structural phylogenies have the potential to fill this gap. However, it is difficult to assess evolutionary relationships derived from structural phylogenies without some means of assessing confidence in such trees. We therefore address two shortcomings in the application of structural data to deep phylogeny. First, we examine whether phylogenies derived from pairwise structural comparisons are sensitive to differences in protein length and shape. We find that structural phylogenetics is best employed where structures have very similar lengths, and that shape fluctuations generated during molecular dynamics simulations impact pairwise comparisons, but not so drastically as to eliminate evolutionary signal. Second, we address the absence of statistical support for structural phylogeny. We present a method for assessing confidence in a structural phylogeny using shape fluctuations generated via molecular dynamics or Monte Carlo simulations of proteins. Our approach will aid the evolutionary reconstruction of relationships across structurally defined protein superfamilies. With the Protein Data Bank now containing in excess of 158,000 entries (December 2019), we predict that structural phylogenetics will become a useful tool for ordering the protein universe.

Unravelling the Complex Duplication History of Deuterostome Glycerol Transporters

Transmembrane glycerol transport is an ancient biophysical property that evolved in selected subfamilies of water channel (aquaporin) proteins. Here, we conducted broad level genome (>550) and transcriptome (>300) analyses to unravel the duplication history of the glycerol-transporting channels (glps) in Deuterostomia. We found that tandem duplication (TD) was the major mechanism of gene expansion in echinoderms and hemichordates, which, together with whole genome duplications (WGD) in the chordate lineage, continued to shape the genomic repertoires in craniates. Molecular phylogenies indicated that aqp3-like and aqp13-like channels were the probable stem subfamilies in craniates, with WGD generating aqp9 and aqp10 in gnathostomes but aqp7 arising through TD in Osteichthyes. We uncovered separate examples of gene translocations, gene conversion, and concerted evolution in humans, teleosts, and starfishes, with DNA transposons the likely drivers of gene rearrangements in paleotetraploid salmonids. Currently, gene copy numbers and BLAST are poor predictors of orthologous relationships due to asymmetric glp gene evolution in the different lineages. Such asymmetries can impact estimations of divergence times by millions of years. Experimental investigations of the salmonid channels demonstrated that approximately half of the 20 ancestral paralogs are functional, with neofunctionalization occurring at the transcriptional level rather than the protein transport properties. The combined findings resolve the origins and diversification of glps over >800 million years old and thus form the novel basis for proposing a pandeuterostome glp gene nomenclature.

Multiplicity and Polymorphism of Fish Hemoglobins

The diversity of fish hemoglobins and the association with oxygen availability and physiological requirements during the life cycle has attracted scientists since the first report on multiple hemoglobin in fishes (Buhler and Shanks 1959). The functional heterogeneity of the fish hemoglobins enables many species to tolerate hypoxic conditions and exhausting swimming, but also to maintain the gas pressure in the swim bladder at large depths. The hemoglobin repertoire has further increased in various species displaying polymorphic hemoglobin variants differing in oxygen binding properties. The multiplicity of fish hemoglobins as particularly found in the tetraploid salmonids strongly contrasts with the complete loss of hemoglobins in Antarctic icefishes and illustrates the adaptive radiation in the oxygen transport of this successful vertebrate group.

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.

Prevalence of hemoglobinopathies in the Brazilian adult population: National Health Survey 2014-2015

OBJECTIVE

To describe the prevalence of hemoglobinopathies in the Brazilian adult population, according to laboratory tests from the National Health Survey.

METHODS

A descriptive study was carried out with National Health Survey laboratory data collected between 2014 and 2015. The hemoglobinopathies test was performed using the High Performance Liquid Chromatography method. The results of the individual tests were interpreted as providing normal, homozygous or heterozygous results for S, C and D hemoglobin, in addition to other possible hemoglobinopathies. Prevalence of hemoglobinopathies according to gender, skin color, region, age and schooling was estimated.

RESULTS

Hemoglobinopathies were present in 3.7% of the population. The main ones were the sickle cell trait (2.49%), thalassemia minor (0.30%) and suspected thalassemia major (0.80%). In relation to the sickle cell trait and suspected thalassemia major, there was a statistically significant difference for the skin color variable (p<0.05). The prevalences found for sickle cell trait according to skin color was: 4.1% among dark-skinned blacks, 3.6% among light-skinned blacks, 1.2% among whites, and 1.7% among others.

CONCLUSION

The most prevalent hemoglobinopathies were the sickle cell trait and minor thalassemia, and were predominate among light- and dark-skinned black people. The study helps in identifying hemoglobinopathies and in genetic counseling in pre-conception.

Ancient Hybridization and Adaptive Introgression of an Invadolysin Gene in Schistosome Parasites

Introgression among parasite species has the potential to transfer traits of biomedical importance across species boundaries. The parasitic blood fluke Schistosoma haematobium causes urogenital schistosomiasis in humans across sub-Saharan Africa. Hybridization with other schistosome species is assumed to occur commonly, because genetic crosses between S. haematobium and livestock schistosomes, including S. bovis, can be staged in the laboratory, and sequencing of mtDNA and rDNA amplified from microscopic miracidia larvae frequently reveals markers from different species. However, the frequency, direction, age, and genomic consequences of hybridization are unknown. We hatched miracidia from eggs and sequenced the exomes from 96 individual S. haematobium miracidia from infected patients from Niger and the Zanzibar archipelago. These data revealed no evidence for contemporary hybridization between S. bovis and S. haematobium in our samples. However, all Nigerien S. haematobium genomes sampled show hybrid ancestry, with 3.3-8.2% of their nuclear genomes derived from S. bovis, providing evidence of an ancient introgression event that occurred at least 108-613 generations ago. Some S. bovis-derived alleles have spread to high frequency or reached fixation and show strong signatures of directional selection; the strongest signal spans a single gene in the invadolysin gene family (Chr. 4). Our results suggest that S. bovis/S. haematobium hybridization occurs rarely but demonstrate profound consequences of ancient introgression from a livestock parasite into the genome of S. haematobium, the most prevalent schistosome species infecting humans.

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.

Cluster expansion of apolipoprotein D (ApoD) genes in teleost fishes

BACKGROUND

Gene and genome duplication play important roles in the evolution of gene function. Compared to individual duplicated genes, gene clusters attract particular attention considering their frequent associations with innovation and adaptation. Here, we report for the first time the expansion of the apolipoprotein D (ApoD) ligand-transporter genes in a cluster manner specific to teleost fishes.

RESULTS

Based on comparative genomic and transcriptomic analyses, protein 3D structure comparison, positive selection detection and breakpoints detection, the single ApoD gene in the ancestor expanded into two clusters following a dynamic evolutionary pattern in teleost fishes. Orthologous genes show conserved expression patterns, whereas lineage-specific duplicated genes show tissue-specific expression patterns and even evolve new gene expression profiles. Positive selection occurred in branches before and after gene duplication, especially for lineage-specific duplicated genes. Cluster analyses based on protein 3D structure comparisons, especially comparisons of the four loops at the opening side, show gene duplication-segregating patterns. Duplicated ApoD genes are predicted to be associated with forkhead transcription factors and MAPK genes. ApoD clusters are located next to the breakpoints of genome rearrangements.

CONCLUSIONS

Here, we report the expansion of ApoD genes specific to teleost fishes in a cluster manner for the first time. Neofunctionalization and subfunctionalization were observed at both the protein and expression levels after duplication. Evidence from different aspects-i.e., abnormal expression-induced disease in humans, fish-specific expansion, predicted associations with forkhead transcription factors and MAPK genes, specific expression patterns in tissues related to sexual selection and adaptation, duplicated genes under positive selection and their location next to the breakpoints of genome rearrangements-suggests the potentially advantageous roles of ApoD genes in teleost fishes. The cluster expansion of ApoD genes specific to teleost fishes provides thus an ideal evo-devo model for studying gene duplication, cluster maintenance and new gene function emergence.

Dynamics based clustering of globin family members

A methodology to cluster proteins based on their dynamics’ similarity is presented. For each pair of proteins from a dataset, the structures are superimposed, and the Anisotropic Network Model modes of motions are calculated. The twelve slowest modes from each protein are matched using a local mode alignment algorithm based on the local sequence alignment algorithm of Smith–Waterman. The dynamical similarity distance matrix is calculated based on the top scoring matches of each pair and the proteins are clustered using a hierarchical clustering algorithm. The utility of this method is exemplified on a dataset of protein chains from the globin family and a dataset of tetrameric hemoglobins. The results demonstrate the effect of the quaternary structure of globin members on their intrinsic dynamics and show good ability to distinguish between different states of hemoglobin, revealing the dynamical relations between them.

Dynamical comparison between myoglobin and hemoglobin

Myoglobin and hemoglobin are globular hemeproteins, when the former is a monomer and the latter a heterotetramer. Despite the structural similarity of myoglobin to α and β subunits of hemoglobin, there is a functional difference between the two proteins, owing to the quaternary structure of hemoglobin. The effect of the quaternary structure of hemoglobin on the intrinsic dynamics of its subunits is explored by dynamical comparison of the two proteins. Anisotropic Network Model modes of motion were calculated for hemoglobin and myoglobin. Dynamical comparison between the proteins was performed using global and local Anisotropic Network Model mode alignment algorithms based on the algorithms of Smith-Waterman and Needleman-Wunsch for sequence comparison. The results indicate that the quaternary structure of hemoglobin substantially alters the intrinsic dynamics of its subunits, an effect that may contribute to the functional difference between the two proteins. Local dynamics similarity between the proteins is still observed at the major exit route of the ligand.

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.

The Diversity of Mammalian Hemoproteins and Microbial Heme Scavengers Is Shaped by an Arms Race for Iron Piracy

Iron is an essential micronutrient for most living species. In mammals, hemoglobin (Hb) stores more than two thirds of the body's iron content. In the bloodstream, haptoglobin (Hp) and hemopexin (Hpx) sequester free Hb or heme. Pathogenic microorganisms usually acquire iron from their hosts and have evolved complex systems of iron piracy to circumvent nutritional immunity. Herein, we performed an evolutionary analysis of genes coding for mammalian heme-binding proteins and heme-scavengers in pathogen species. The underlying hypothesis is that these molecules are engaged in a molecular arms race. We show that positive selection drove the evolution of mammalian Hb and Hpx. Positively selected sites in Hb are located at the interaction surface with Neisseria meningitidis heme scavenger HpuA and with Staphylococcus aureus iron-regulated surface determinant B (IsdB). In turn, positively selected sites in HpuA and IsdB are located in the flexible protein regions that contact Hb. A residue in Hb (S45H) was also selected on the Caprinae branch. This site stabilizes the interaction with Trypanosoma brucei hemoglobin-haptoglobin (HbHp) receptor (TbHpHbR), a molecule that also mediates trypanosome lytic factor (TLF) entry. In TbHpHbR, positive selection drove the evolution of a variant (L210S) which allows evasion from TLF but reduces affinity for HbHp. Finally, selected sites in Hpx are located at the interaction surface with the Haemophilus influenzae hemophore HxuA, which in turn displays fast evolving sites at the Hpx-binding interface. These results shed light into host-pathogens conflicts and establish the importance of nutritional immunity as an evolutionary force.

The Reprimo Gene Family: A Novel Gene Lineage in Gastric Cancer with Tumor Suppressive Properties

The reprimo (RPRM) gene family is a group of single exon genes present exclusively within the vertebrate lineage. Two out of three members of this family are present in humans: RPRM and RPRM-Like (RPRML). RPRM induces cell cycle arrest at G2/M in response to p53 expression. Loss-of-expression of RPRM is related to increased cell proliferation and growth in gastric cancer. This evidence suggests that RPRM has tumor suppressive properties. However, the molecular mechanisms and signaling partners by which RPRM exerts its functions remain unknown. Moreover, scarce studies have attempted to characterize RPRML, and its functionality is unclear. Herein, we highlight the role of the RPRM gene family in gastric carcinogenesis, as well as its potential applications in clinical settings. In addition, we summarize the current knowledge on the phylogeny and expression patterns of this family of genes in embryonic zebrafish and adult humans. Strikingly, in both species, RPRM is expressed primarily in the digestive tract, blood vessels and central nervous system, supporting the use of zebrafish for further functional characterization of RPRM. Finally, drawing on embryonic and adult expression patterns, we address the potential relevance of RPRM and RPRML in cancer. Active investigation or analytical research in the coming years should contribute to novel translational applications of this poorly understood gene family as potential biomarkers and development of novel cancer therapies.

Evolution and Expression of Tissue Globins in Ray-Finned Fishes

The globin gene family encodes oxygen-binding hemeproteins conserved across the major branches of multicellular life. The origins and evolutionary histories of complete globin repertoires have been established for many vertebrates, but there remain major knowledge gaps for ray-finned fish. Therefore, we used phylogenetic, comparative genomic and gene expression analyses to discover and characterize canonical “non-blood” globin family members (i.e., myoglobin, cytoglobin, neuroglobin, globin-X, and globin-Y) across multiple ray-finned fish lineages, revealing novel gene duplicates (paralogs) conserved from whole genome duplication (WGD) and small-scale duplication events. Our key findings were that: (1) globin-X paralogs in teleosts have been retained from the teleost-specific WGD, (2) functional paralogs of cytoglobin, neuroglobin, and globin-X, but not myoglobin, have been conserved from the salmonid-specific WGD, (3) triplicate lineage-specific myoglobin paralogs are conserved in arowanas (Osteoglossiformes), which arose by tandem duplication and diverged under positive selection, (4) globin-Y is retained in multiple early branching fish lineages that diverged before teleosts, and (5) marked variation in tissue-specific expression of globin gene repertoires exists across ray-finned fish evolution, including several previously uncharacterized sites of expression. In this respect, our data provide an interesting link between myoglobin expression and the evolution of air breathing in teleosts. Together, our findings demonstrate great-unrecognized diversity in the repertoire and expression of nonblood globins that has arisen during ray-finned fish evolution.

Functional explanation and the problem of functional equivalence

The legitimacy of functional explanations in biology is threatened by a problem first identified by Hempel: the problem of functional equivalence. In order for the prevalence of a trait to be explained by its function, the function would have to explain why that very trait is prevalent and not some other functionally equivalent trait. But functions alone cannot meet this explanatory demand. I argue that this is a problem not only for Nagelian deductive-nomological models but also for etiological models of functional explanation. I contrast these models with a dual model of adaptive explanation and design explanation. This dual model largely circumvents the problem of functional equivalence, but divests functions of much explanatory power.

Evolution of the α2-adrenoreceptors in vertebrates: ADRA2D is absent in mammals and crocodiles

Evolutionary studies of genes that have been functionally characterized and whose variation has been associated with pathological conditions represent an opportunity to understand the genetic basis of pathologies. α₂-Adrenoreceptors (ADRA2) are a class of G protein-coupled receptors that regulate several physiological processes including blood pressure, platelet aggregation, insulin secretion, lipolysis, and neurotransmitter release. This gene family has been extensively studied from a molecular/physiological perspective, yet much less is known about its evolutionary history. Accordingly, the goal of this study was to investigate the evolutionary history of α₂-adrenoreceptors (ADRA2) in vertebrates. Our results show that in addition to the three well-recognized α₂-adrenoreceptor genes (ADRA2A, ADRA2B and ADRA2C), we recovered a clade that corresponds to the fourth member of the α₂-adrenoreceptor gene family (ADRA2D). We also recovered a clade that possesses two ADRA2 sequences found in two lamprey species. Furthermore, our results show that mammals and crocodiles are characterized by possessing three α₂-adrenoreceptor genes, whereas all other vertebrate groups possess the full repertoire of α₂-adrenoreceptor genes. Among vertebrates ADRA2D seems to be a dispensable gene, as it was lost two independent times during the evolutionary history of the group. Additionally, we found that most examined species possess the most common alleles described for humans; however, there are cases in which non-human mammals possess the alternative variant. Finally, transcript abundance profiles revealed that during the early evolutionary history of gnathostomes, the expression of ADRA2D in different taxonomic groups became specialized to different tissues, but in the ancestor of sarcopterygians this specialization would have been lost.

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.

The Phenotypic Plasticity of Duplicated Genes in Saccharomyces cerevisiae and the Origin of Adaptations

Gene and genome duplication are the major sources of biological innovations in plants and animals. Functional and transcriptional divergence between the copies after gene duplication has been considered the main driver of innovations . However, here we show that increased phenotypic plasticity after duplication plays a more major role than thought before in the origin of adaptations. We perform an exhaustive analysis of the transcriptional alterations of duplicated genes in the unicellular eukaryote Saccharomyces cerevisiae when challenged with five different environmental stresses. Analysis of the transcriptomes of yeast shows that gene duplication increases the transcriptional response to environmental changes, with duplicated genes exhibiting signatures of adaptive transcriptional patterns in response to stress. The mechanism of duplication matters, with whole-genome duplicates being more transcriptionally altered than small-scale duplicates. The predominant transcriptional pattern follows the classic theory of evolution by gene duplication; with one gene copy remaining unaltered under stress, while its sister copy presents large transcriptional plasticity and a prominent role in adaptation. Moreover, we find additional transcriptional profiles that are suggestive of neo- and subfunctionalization of duplicate gene copies. These patterns are strongly correlated with the functional dependencies and sequence divergence profiles of gene copies. We show that, unlike singletons, duplicates respond more specifically to stress, supporting the role of natural selection in the transcriptional plasticity of duplicates. Our results reveal the underlying transcriptional complexity of duplicated genes and its role in the origin of adaptations.

Copy Number Variants in Alzheimer's Disease

Alzheimer's disease (AD) is a devastating disease mainly afflicting elderly people, characterized by decreased cognition, loss of memory, and eventually death. Although risk and deterministic genes are known, major genetics research programs are underway to gain further insights into the inheritance of AD. In the last years, in particular, new developments in genome-wide scanning methodologies have enabled the association of a number of previously uncharacterized copy number variants (CNVs, gain or loss of DNA) in AD. Because of the exceedingly large number of studies performed, it has become difficult for geneticists as well as clinicians to systematically follow, evaluate, and interpret the growing number of (sometime conflicting) CNVs implicated in AD. In this review, after a brief introduction of this type of structural variation, and a description of available databases, computational analyses, and technologies involved, we provide a systematic review of all published data showing statistical and scientific significance of pathogenic CNVs and discuss the role they might play in AD.

The mammalian skeletal muscle DHPR has larger Ca2+ conductance and is phylogenetically ancient to the early ray-finned fish sterlet (Acipenser ruthenus)

The L-type Ca²⁺ channel or dihydropyridine receptor (DHPR) in vertebrate skeletal muscle is responsible for sensing sarcolemmal depolarizations and transducing this signal to the sarcoplasmic Ca²⁺ release channel RyR1 via conformational coupling to initiate muscle contraction. During this excitation-contraction (EC) coupling process there is a slow Ca²⁺ current through the mammalian DHPR which is fully missing in euteleost fishes. In contrast to ancestral evolutionary stages where skeletal muscle EC coupling is still depended on Ca²⁺-induced Ca²⁺-release (CICR), it is possible that the DHPR Ca²⁺ conductivity during mammalian (conformational) EC coupling was retained as an evolutionary remnant (vestigiality). Here, we wanted to test the hypothesis that due to the lack of evolutionary pressure in post-CICR species skeletal muscle DHPR Ca²⁺ conductivity gradually reduced as evolution progressed. Interestingly, we identified that the DHPR of the early ray-finned fish sterlet (Acipenser ruthenus) is phylogenetically positioned above the mammalian rabbit DHPR which retained robust Ca²⁺ conductivity, but below the euteleost zebrafish DHPR which completely lost Ca²⁺ conductivity. Remarkably, our results revealed that sterlet DHPR still retained the Ca²⁺ conductivity but currents are significantly reduced compared to rabbit. This decrease is due to lower DHPR membrane expression similar to zebrafish, as well as due to reduced channel open probability (P_o). In both these fish species the lower DHPR expression density is partially compensated by higher efficacy of DHPR-RyR1 coupling. The complete loss of P_o in zebrafish and other euteleost species was presumably based on the teleost specific 3rd round of genome duplication (Ts3R). Ts3R headed into the appearance of two skeletal muscle DHPR isoforms which finally, together with the radiation of the euteleost clade, fully lost the P_o.

Hemoglobin-oxygen affinity in high-altitude vertebrates: is there evidence for an adaptive trend?

In air-breathing vertebrates at high altitude, fine-tuned adjustments in hemoglobin (Hb)-O2 affinity provide an energetically efficient means of mitigating the effects of arterial hypoxemia. However, it is not always clear whether an increased or decreased Hb-O2 affinity should be expected to improve tissue O2 delivery under different degrees of hypoxia, due to the inherent trade-off between arterial O2 loading and peripheral O2 unloading. Theoretical results indicate that the optimal Hb-O2 affinity varies as a non-linear function of environmental O2 availability, and the threshold elevation at which an increased Hb-O2 affinity becomes advantageous depends on the magnitude of diffusion limitation (the extent to which O2 equilibration at the blood-gas interface is limited by the kinetics of O2 exchange). This body of theory provides a framework for interpreting the possible adaptive significance of evolved changes in Hb-O2 affinity in vertebrates that have colonized high-altitude environments. To evaluate the evidence for an empirical generalization and to test theoretical predictions, I synthesized comparative data in a phylogenetic framework to assess the strength of the relationship between Hb-O2 affinity and native elevation in mammals and birds. Evidence for a general trend in mammals is equivocal, but there is a remarkably strong positive relationship between Hb-O2 affinity and native elevation in birds. Evolved changes in Hb function in high-altitude birds provide one of the most compelling examples of convergent biochemical adaptation in vertebrates.

The Primary Structure of β(I)-Chain of Hemoglobin from Snake Sindhi Krait (Bungarus sindanus sindanus)

The amino acid sequence of β(I)-globin chain from Sindhi Krait (Bungarus sindanus sindanus) was determined to study the molecular evolution among snakes. The hemoglobin was isolated from the red blood cells and was analyzed by ion-exchange chromatography (IEX). The crude globin was subjected to reversed phased-high performance liquid chromatography (RP-HPLC) using C4 column. The N-terminal sequences of intact globin chains and tryptic peptides were determined by Edman degradation in a pulsed liquid gas phase sequencer using an online Phenylthiohydantoin analyzer. Sindhi Krait is expected to express three hemoglobin components that are composed of β(II), β(I), α(D) and α(A)-globin chains, as apparent by IEX, RP-HPLC and N-terminal sequence analyses. Sequence alignment and phylogenetic analyses of β(I) globin chain from Sindhi Krait showed closest relationship with β(I) globin chain from Rattlesnake, Water snake and Indigo snake. Interestingly, comparison of primary sequence of β(I) globin chain of Sindhi Krait with human β chain revealed 63 % similarity along with the retention of all heme contact points. Variations among the two sequences were prominent at αβ contact points and in regions directly not important for function.

Appetite regulation in Schizothorax prenanti by three CART genes

In recent years, cocaine- and amphetamine-regulated transcript (CART) has received much attention as mediators of appetite regulation in mammals. However, the involvement of CART in the feeding behavior of teleosts has not been well understood. In this study, three distinct CARTs were cloned from the Schizothorax prenanti (S. prenanti). Real-time quantitative PCR were applied to characterize the tissue distribution and appetite regulatory effects of CARTs in S. prenanti. The S. prenanti CART-1, CART-2 and CART-3 full-length cDNA sequences were 597 bp, 694 bp and 749 bp in length, encoding the peptides of 125, 120 and 104 amino acid residues, respectively. All the S. prenanti CARTs consisted of three exons and two introns. Tissue distribution analysis showed that the high mRNA levels of S. prenanti CART-1 were observed in the telencephalon and eye, followed by the hypothalamus, myelencephalon, and mesencephalon. The S. prenanti CART-2 mRNA was mainly found in the mesencephalon, hypothalamus, telencephalon and myelencephalon. The S. prenanti CART-3 mRNA was widely distributed among the tissues, with the high levels in the hypothalamus and foregut. In the periprandial experiment, all three CARTs mRNA expressions in the hypothalamus were highly elevated after a meal, suggesting that CARTs are postprandial satiety signals. In the fasting experiment, all three CARTs mRNA expressions decreased after fasting and increased after refeeding, suggesting that CARTs might be involved in regulation of appetite in the S. prenanti.

Contribution of a mutational hot spot to hemoglobin adaptation in high-altitude Andean house wrens

A key question in evolutionary genetics is why certain mutations or certain types of mutation make disproportionate contributions to adaptive phenotypic evolution. In principle, the preferential fixation of particular mutations could stem directly from variation in the underlying rate of mutation to function-altering alleles. However, the influence of mutation bias on the genetic architecture of phenotypic evolution is difficult to evaluate because data on rates of mutation to function-altering alleles are seldom available. Here, we report the discovery that a single point mutation at a highly mutable site in the β(A)-globin gene has contributed to an evolutionary change in hemoglobin (Hb) function in high-altitude Andean house wrens (Troglodytes aedon). Results of experiments on native Hb variants and engineered, recombinant Hb mutants demonstrate that a nonsynonymous mutation at a CpG dinucleotide in the β(A)-globin gene is responsible for an evolved difference in Hb-O2 affinity between high- and low-altitude house wren populations. Moreover, patterns of genomic differentiation between high- and low-altitude populations suggest that altitudinal differentiation in allele frequencies at the causal amino acid polymorphism reflects a history of spatially varying selection. The experimental results highlight the influence of mutation rate on the genetic basis of phenotypic evolution by demonstrating that a large-effect allele at a highly mutable CpG site has promoted physiological differentiation in blood O2 transport capacity between house wren populations that are native to different elevations.