Evolution of lysine-specific demethylase 1 and REST corepressor gene families and their molecular interaction

Lysine-specific demethylase 1A (LSD1) binds to the REST corepressor (RCOR) protein family of corepressors to erase transcriptionally active marks on histones. Functional diversity in these complexes depends on the type of RCOR included, which modulates the catalytic activity of the complex. Here, we studied the duplicative history of the RCOR and LSD gene families and analyzed the evolution of their interaction. We found that RCOR genes are the product of the two rounds of whole-genome duplications that occurred early in vertebrate evolution. In contrast, the origin of the LSD genes traces back before to the divergence of animals and plants. Using bioinformatics tools, we show that the RCOR and LSD1 interaction precedes the RCOR repertoire expansion that occurred in the last common ancestor of jawed vertebrates. Overall, we trace LSD1-RCOR complex evolution and propose that animal non-model species offer advantages in addressing questions about the molecular biology of this epigenetic complex.

Discovery of BbX transcription factor in the patagonian blennie: Exploring expression changes following combined bacterial and thermal stress exposure

High-Mobility Group (HMG) proteins are involved in different processes such as transcription, replication, DNA repair, and immune response. The role of HMG proteins in the immune response of fish has been studied mainly for HMGB1, where its expression can be induced by the stimulation of viral/bacterial PAMPs and can act as a proinflammatory mediator and as a global regulator of transcription in response to temperature. However, for BbX this role remains to be discovered. In this work, we identified the BbX of E. maclovinus and evaluated the temporal expression levels after simultaneous challenge with P. salmonis and thermal stress. Phylogenetic analysis does not significantly deviate from the expected organismal relationships suggesting orthologous relationships and that BbX was present in the common ancestor of the group. BbX mRNA expression levels were very high in the intestinal tissue of E. maclovinus (foregut, midgut, and hindgut). Nevertheless, the protein levels analyzed by WB showed the highest levels of BbX protein in the liver (constitutive expression). On the other hand, the mRNA expression levels of BbX in the liver of E. maclovinus injected with P. salmonis and subjected to thermal stress showed an increase at days 16 and 20 in all treatments applied at 12 °C and 18 °C. Meanwhile, the protein levels quantified by WB showed a statistically significant increase in the HMG-Bbx at all experimental times (4, 8, 12, 16, and 20 dpi). However, at 4 dpi the HMG-Bbx protein levels were much higher than the other days evaluated. The results suggest that BbX protein may be implicated in the response mechanism to temperature and bacterial stimulation in the foregut, midgut, hindgut, and liver, according to our findings at the level of mRNA and protein. Furthermore, our WB analysis suggests an effect of P. salmonis on the expression of this protein that can be observed in condition C+ 12 °C compared to C- 12 °C. Then, there is an effect of temperature that can be evidenced in the condition AM 18 °C and SM 18 °C, compared to AB 18 °C and SB 18 °C at 4, 8, and 12 dpi. We found not differences in the levels of this protein if the thermal stress is achieved through acclimatization or shock. More research is necessary to clarify the importance of this type of HMG in the immune response and thermal tolerance in fish.

Genomic insights into the mystery of mouse mummies on the summits of Atacama volcanoes

Our understanding of the limits of animal life is continually revised by scientific exploration of extreme environments. Here we report the discovery of mummified cadavers of leaf-eared mice, Phyllotis vaccarum, from the summits of three different Andean volcanoes at elevations 6,029-6,233 m above sea level in the Puna de Atacama in Chile and Argentina. Such extreme elevations were previously assumed to be completely uninhabitable by mammals. In combination with a live-captured specimen of the same species from the nearby summit of Volcán Llullaillaco (6,739 m)1, the summit mummies represent the highest altitude physical records of mammals in the world. We also report a chromosome-level genome assembly for P. vaccarum that, in combination with a whole-genome re-sequencing analysis and radiocarbon dating analysis, provides insights into the provenance and antiquity of the summit mice. Radiocarbon data indicate that the most ancient of the mummies are, at most, a few centuries old. Genomic polymorphism data revealed a high degree of continuity between the summit mice and conspecifics from lower elevations in the surrounding Altiplano. Genomic data also revealed equal numbers of males and females among the summit mice and evidence of close kinship between some individuals from the same summits. These findings bolster evidence for resident populations of Phyllotis at elevations >6,000 m and challenge assumptions about the environmental limits of vertebrate life and the physiological tolerances of small mammals.

Extreme high-elevation mammal surveys reveal unexpectedly high upper range limits of Andean mice

In the world's highest mountain ranges, uncertainty about the upper elevational range limits of alpine animals represents a critical knowledge gap regarding the environmental limits of life and presents a problem for detecting range shifts in response to climate change. Here we report results of mountaineering mammal surveys in the Central Andes, which led to the discovery of multiple species of mice living at extreme elevations that far surpass previously assumed range limits for mammals. We live-trapped small mammals from ecologically diverse sites spanning >6700 m of vertical relief, from the desert coast of northern Chile to the summits of the highest volcanoes in the Andes. We used molecular sequence data and whole-genome sequence data to confirm the identities of species that represent new elevational records and to test hypotheses regarding species limits. These discoveries contribute to a new appreciation of the environmental limits of vertebrate life.

How Many Sirtuin Genes Are Out There? Evolution of Sirtuin Genes in Vertebrates With a Description of a New Family Member

Studying the evolutionary history of gene families is a challenging and exciting task with a wide range of implications. In addition to exploring fundamental questions about the origin and evolution of genes, disentangling their evolution is also critical to those who do functional/structural studies to allow a deeper and more precise interpretation of their results in an evolutionary context. The sirtuin gene family is a group of genes that are involved in a variety of biological functions mostly related to aging. Their duplicative history is an open question, as well as the definition of the repertoire of sirtuin genes among vertebrates. Our results show a well-resolved phylogeny that represents an improvement in our understanding of the duplicative history of the sirtuin gene family. We identified a new sirtuin gene family member (SIRT3.2) that was apparently lost in the last common ancestor of amniotes but retained in all other groups of jawed vertebrates. According to our experimental analyses, elephant shark SIRT3.2 protein is located in mitochondria, the overexpression of which leads to an increase in cellular levels of ATP. Moreover, in vitro analysis demonstrated that it has deacetylase activity being modulated in a similar way to mammalian SIRT3. Our results indicate that there are at least eight sirtuin paralogs among vertebrates and that all of them can be traced back to the last common ancestor of the group that existed between 676 and 615 millions of years ago.

Sequence and structural conservation reveal fingerprint residues in TRP channels

Transient receptor potential (TRP) proteins are a large family of cation-selective channels, surpassed in variety only by voltage-gated potassium channels. Detailed molecular mechanisms governing how membrane voltage, ligand binding, or temperature can induce conformational changes promoting the open state in TRP channels are still a matter of debate. Aiming to unveil distinctive structural features common to the transmembrane domains within the TRP family, we performed phylogenetic reconstruction, sequence statistics, and structural analysis over a large set of TRP channel genes. Here, we report an exceptionally conserved set of residues. This fingerprint is composed of twelve residues localized at equivalent three-dimensional positions in TRP channels from the different subtypes. Moreover, these amino acids are arranged in three groups, connected by a set of aromatics located at the core of the transmembrane structure. We hypothesize that differences in the connectivity between these different groups of residues harbor the apparent differences in coupling strategies used by TRP subgroups.

Evolution of the DAN gene family in vertebrates

The DAN gene family (DAN, Differential screening-selected gene Aberrant in Neuroblastoma) is a group of genes that is expressed during development and plays fundamental roles in limb bud formation and digitation, kidney formation and morphogenesis and left-right axis specification. During adulthood the expression of these genes are associated with diseases, including cancer. Although most of the attention to this group of genes has been dedicated to understanding its role in physiology and development, its evolutionary history remains poorly understood. Thus, the goal of this study is to investigate the evolutionary history of the DAN gene family in vertebrates, with the objective of complementing the already abundant physiological information with an evolutionary context. Our results recovered the monophyly of all DAN gene family members and divide them into five main groups. In addition to the well-known DAN genes, our phylogenetic results revealed the presence of two new DAN gene lineages; one is only retained in cephalochordates, whereas the other one (GREM3) was only identified in cartilaginous fish, holostean fish, and coelacanth. According to the phyletic distribution of the genes, the ancestor of gnathostomes possessed a repertoire of eight DAN genes, and during the radiation of the group GREM1, GREM2, SOST, SOSTDC1, and NBL1 were retained in all major groups, whereas, GREM3, CER1, and DAND5 were differentially lost.

Evolutionary history of the vertebrate Piwi gene family

PIWIs are regulatory proteins that belong to the Argonaute family. Piwis are primarily expressed in gonads and protect the germline against the mobilization and propagation of transposable elements (TEs) through transcriptional gene silencing. Vertebrate genomes encode up to four Piwi genes: Piwil1, Piwil2, Piwil3 and Piwil4, but their duplication history is unresolved. We leveraged phylogenetics, synteny and expression analyses to address this void. Our phylogenetic analysis suggests Piwil1 and Piwil2 were retained in all vertebrate members. Piwil4 was the result of Piwil1 duplication in the ancestor of gnathostomes, but was independently lost in ray-finned fishes and birds. Further, Piwil3 was derived from a tandem Piwil1 duplication in the common ancestor of marsupial and placental mammals, but was secondarily lost in Atlantogenata (Xenarthra and Afrotheria) and some rodents. The evolutionary rate of Piwil3 is considerably faster than any Piwi among all lineages, but an explanation is lacking. Our expression analyses suggest Piwi expression has mostly been constrained to gonads throughout vertebrate evolution. Vertebrate evolution is marked by two early rounds of whole genome duplication and many multigene families are linked to these events. However, our analyses suggest Piwi expansion was independent of whole genome duplications.

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.

Positive selection and gene duplications in tumour suppressor genes reveal clues about how cetaceans resist cancer

Cetaceans are the longest-living species of mammals and the largest in the history of the planet. They have developed mechanisms against diseases such cancer, although the underlying molecular bases of these remain unknown. The goal of this study was to investigate the role of natural selection in the evolution of 1077 tumour suppressor genes (TSGs) in cetaceans. We used a comparative genomic approach to analyse two sources of molecular variation in the form of d_N/d_S rates and gene copy number variation. We found a signal of positive selection in the ancestor of cetaceans within the CXCR2 gene, an important regulator of DNA damage, tumour dissemination and immune system. Further, in the ancestor of baleen whales, we found six genes exhibiting positive selection relating to diseases such as breast carcinoma, lung neoplasm (ADAMTS8) and leukaemia (ANXA1). The TSGs turnover rate (gene gain and loss) was almost 2.4-fold higher in cetaceans when compared with other mammals, and notably even faster in baleen whales. The molecular variants in TSGs found in baleen whales, combined with the faster gene turnover rate, could have favoured the evolution of their particular traits of anti-cancer resistance, gigantism and longevity. Additionally, we report 71 genes with duplications, of which 11 genes are linked to longevity (e.g. NOTCH3 and SIK1) and are important regulators of senescence, cell proliferation and metabolism. Overall, these results provide evolutionary evidence that natural selection in TSGs could act on species with large body sizes and extended lifespan, providing novel insights into the genetic basis of disease resistance.

The Globin Gene Family in Arthropods: Evolution and Functional Diversity

Globins are small heme-proteins that reversibly bind oxygen. Their most prominent roles in vertebrates are the transport and storage of O₂ for oxidative energy metabolism, but recent research has suggested alternative, non-respiratory globin functions. In the species-rich and ecologically highly diverse taxon of arthropods, the copper-containing hemocyanin is considered the main respiratory protein. However, recent studies have suggested the presence of globin genes and their proteins in arthropod taxa, including model species like Drosophila. To systematically assess the taxonomic distribution, evolution and diversity of globins in arthropods, we systematically searched transcriptome and genome sequence data and found a conserved, widespread occurrence of three globin classes in arthropods: hemoglobin-like (HbL), globin X (GbX), and globin X-like (GbXL) protein lineages. These globin types were previously identified in protostome and deuterostome animals including vertebrates, suggesting their early ancestry in Metazoa. The HbL genes show multiple, lineage-specific gene duplications in all major arthropod clades. Some HbL genes (e.g., Glob2 and 3 of Drosophila) display particularly fast substitution rates, possibly indicating the evolution of novel functions, e.g., in spermatogenesis. In contrast, arthropod GbX and GbXL globin genes show high evolutionary stability: GbXL is represented by a single-copy gene in all arthropod groups except Brachycera, and representatives of the GbX clade are present in all examined taxa except holometabolan insects. GbX and GbXL both show a brain-specific expression. Most arthropod GbX and GbXL proteins, but also some HbL variants, include sequence motifs indicative of potential N-terminal acylation (i.e., N-myristoylation, 3C-palmitoylation). All arthropods except for the brachyceran Diptera harbor at least one such potentially acylated globin copy, confirming the hypothesis of an essential, conserved globin function associated with the cell membrane. In contrast to other animals, the fourth ancient globin lineage, represented by neuroglobin, appears to be absent in arthropods, and the putative arthropod orthologs of the fifth metazoan globin lineage, androglobin, lack a recognizable globin domain. Thus, the remarkable evolutionary stability of some globin variants is contrasted by occasional dynamic gene multiplication or even loss of otherwise strongly conserved globin lineages in arthropod phylogeny.

Discovery of the world's highest-dwelling mammal

Environmental limits of animal life are invariably revised when the animals themselves are investigated in their natural habitats. Here we report results of a scientific mountaineering expedition to survey the high-altitude rodent fauna of Volcán Llullaillaco in the Puna de Atacama of northern Chile, an effort motivated by video documentation of mice (genus Phyllotis) at a record altitude of 6,205 m. Among numerous trapping records at altitudes of >5,000 m, we captured a specimen of the yellow-rumped leaf-eared mouse (Phyllotis xanthopygus rupestris) on the very summit of Llullaillaco at 6,739 m. This summit specimen represents an altitudinal world record for mammals, far surpassing all specimen-based records from the Himalayas and other mountain ranges. This discovery suggests that we may have generally underestimated the altitudinal range limits and physiological tolerances of small mammals simply because the world's high summits remain relatively unexplored by biologists.

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.

The Reprimo gene family member, reprimo-like (rprml), is required for blood development in embryonic zebrafish

The Reprimo gene family comprises a group of single-exon genes for which their physiological function remains poorly understood. Heretofore, mammalian Reprimo (RPRM) has been described as a putative p53-dependent tumor suppressor gene that functions at the G2/M cell cycle checkpoint. Another family member, Reprimo-like (RPRML), has not yet an established role in physiology or pathology. Importantly, RPRML expression pattern is conserved between zebrafish and human species. Here, using CRISPR-Cas9 and antisense morpholino oligonucleotides, we disrupt the expression of rprml in zebrafish and demonstrate that its loss leads to impaired definitive hematopoiesis. The formation of hemangioblasts and the primitive wave of hematopoiesis occur normally in absence of rprml. Later in development there is a significant reduction in erythroid-myeloid precursors (EMP) at the posterior blood island (PBI) and a significant decline of definitive hematopoietic stem/progenitor cells (HSPCs). Furthermore, loss of rprml also increases the activity of caspase-3 in endothelial cells within the caudal hematopoietic tissue (CHT), the first perivascular niche where HSPCs reside during zebrafish embryonic development. Herein, we report an essential role for rprml during hematovascular development in zebrafish embryos, specifically during the definitive waves of hematopoiesis, indicating for the first time a physiological role for the rprml gene.

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.

Phylogenetic evidence for independent origins of GDF1 and GDF3 genes in anurans and mammals

Growth differentiation factors 1 (GDF1) and 3 (GDF3) are members of the transforming growth factor superfamily (TGF-β) that is involved in fundamental early-developmental processes that are conserved across vertebrates. The evolutionary history of these genes is still under debate due to ambiguous definitions of homologous relationships among vertebrates. Thus, the goal of this study was to unravel the evolution of the GDF1 and GDF3 genes of vertebrates, emphasizing the understanding of homologous relationships and their evolutionary origin. Our results revealed that the GDF1 and GDF3 genes found in anurans and mammals are the products of independent duplication events of an ancestral gene in the ancestor of each of these lineages. The main implication of this result is that the GDF1 and GDF3 genes of anurans and mammals are not 1:1 orthologs. In other words, genes that participate in fundamental processes during early development have been reinvented two independent times during the evolutionary history of tetrapods.

An association between differential expression and genetic divergence in the Patagonian olive mouse (Abrothrix olivacea)

Recent molecular studies have found striking differences between desert-adapted species and model mammals regarding water conservation. In particular, aquaporin 4, a classical gene involved in water regulation of model species, is absent or not expressed in the kidneys of desert-adapted species. To further understand the molecular response to water availability, we studied the Patagonian olive mouse Abrothrix olivacea, a species with an unusually broad ecological tolerance that exhibits a great urine concentration capability. The species is able to occupy both the arid Patagonian steppe and the Valdivian and Magellanic forests. We sampled 95 olive mouse specimens from four localities (two in the steppe and two in the forests) and analysed both phenotypic variables and transcriptomic data to investigate the response of this species to the contrasting environmental conditions. The relative size of the kidney and the ratio of urine to plasma concentrations were, as expected, negatively correlated with annual rainfall. Expression analyses uncovered nearly 3,000 genes that were differentially expressed between steppe and forest samples and indicated that this species resorts to the "classical" gene pathways for water regulation. Differential expression across biomes also involves genes that involved in immune and detoxification functions. Overall, genes that were differentially expressed showed a slight tendency to be more divergent and to display an excess of intermediate allele frequencies, relative to the remaining loci. Our results indicate that both differential expression in pathways involved in water conservation and geographical allelic variation are important in the occupation of contrasting habitats by the Patagonian olive mouse.

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 of dopamine receptors: phylogenetic evidence suggests a later origin of the DRD2l and DRD4rs dopamine receptor gene lineages

Dopamine receptors are integral membrane proteins whose endogenous ligand is dopamine. They play a fundamental role in the central nervous system and dysfunction of dopaminergic neurotransmission is responsible for the generation of a variety of neuropsychiatric disorders. From an evolutionary standpoint, phylogenetic relationships among the DRD₁ class of dopamine receptors are still a matter of debate as in the literature different tree topologies have been proposed. In contrast, phylogenetic relationships among the DRD₂ group of receptors are well understood. Understanding the time of origin of the different dopamine receptors is also an issue that needs further study, especially for the genes that have restricted phyletic distributions (e.g., DRD_2l and DRD_4rs). Thus, the goal of this study was to investigate the evolution of dopamine receptors, with emphasis on shedding light on the phylogenetic relationships among the D₁ class of dopamine receptors and the time of origin of the DRD_2l and DRD_4rs gene lineages. Our results recovered the monophyly of the two groups of dopamine receptors. Within the DRD₁ group the monophyly of each paralog was recovered with strong support, and phylogenetic relationships among them were well resolved. Within the DRD₁ class of dopamine receptors we recovered the sister group relationship between the DRD_1C and DRD_1E, and this clade was recovered sister to a cyclostome sequence. The DRD₁ clade was recovered sister to the aforementioned clade, and the group containing DRD₅ receptors was sister to all other DRD₁ paralogs. In agreement with the literature, among the DRD₂ class of receptors, DRD₂ was recovered sister to DRD₃, whereas DRD₄ was sister to the DRD₂/DRD₃ clade. According to our phylogenetic tree, the DRD_2l and DRD_4rs gene lineages would have originated in the ancestor of gnathostomes between 615 and 473 mya. Conservation of sequences required for dopaminergic neurotransmission and small changes in regulatory regions suggest a functional refinement of the dopaminergic pathways along evolution.

Expression of RPRM/rprm in the Olfactory System of Embryonic Zebrafish (Danio rerio)

The Reprimo (RPRM) family is composed of highly conserved single-exon genes. The expression pattern of this gene family has been recently described during zebrafish (Danio rerio) embryogenesis, and primarily locates in the nervous system. Its most characterized member, RPRM, which duplicated to give rise rprma and rprmb in the fish lineage, is known to act as a tumor-suppressor gene in mammalian models. Here, we describe in detail the spatiotemporal expression of three rprm genes (rprma, rprmb, and rprml) within distinct anatomical structures in the developing peripheral and central nervous system. In the zebrafish, rprma mRNA is expressed in the olfactory placodes (OP) and olfactory epithelium (OE), rprmb is observed in the tectum opticum (TeO) and trigeminal ganglion (Tg), whereas rprml is found primarily in the telencephalon (Tel). At protein level, RPRM is present in a subset of cells in the OP, and neurons in the OE, TeO, hindbrain and sensory peripheral structures. Most importantly, the expression of RPRM has been conserved between teleosts and mammals. Thus, we provide a reference dataset describing the expression patterns of RPRM gene products during zebrafish and mouse development as a first step to approach the physiological role of the RPRM gene family.

Progressive erosion of the Relaxin1 gene in bovids

The relaxin/insulin-like (RLN/INSL) gene family is a group of genes that encode peptide hormones involved in a variety of physiological functions related to reproduction. Previous studies have shown that relaxin plays a key role in widening of the pubic bone during labor and in gamete maturation. Because of these functions, studying the evolution of RLN1, the gene encoding for relaxin, is relevant in livestock species, most of which belong in the group Laurasiatheria, which includes cow, pig, horse, goat, and sheep in addition to bats, cetaceans and carnivores. Experimental evidence suggests that cows do not synthesize relaxin, but respond to it, and sheep apparently have a truncated RLN1 gene. Thus, we made use of genome sequence data to characterize the genomic locus of the RLN1 gene in Laurasiatherian mammals to better understand how cows lost the ability to synthesize this peptide. We found that all ruminants in our study (cow, giraffe, goat, sheep and Tibetan antelope) lack a functional RLN1 gene, and document the progressive loss of RLN1 in the lineage leading to cows. Our analyses indicate that 1 - all ruminants have lost all key regulatory elements upstream of the first exon, 2 - giraffe, goat, sheep and Tibetan antelope have multiple inactivating mutations in the RLN1 pseudogene, and 3 - the cow genome has lost all traces of RLN1. The 5' regulatory sequence plays a key role in activating expression, and the loss of this sequence would impair synthesis of mRNA. Our results suggest that changes in regulatory sequence preceded mutations in coding sequence and highlight the importance of these regions in maintaining proper gene function. In addition, we found that all bovids examined posses copies of the relaxin receptors, which explains why they are able to respond to relaxin despite their inability to produce it.

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.

Reprimo tissue-specific expression pattern is conserved between zebrafish and human

Reprimo (RPRM), a member of the RPRM gene family, is a tumor-suppressor gene involved in the regulation of the p53-mediated cell cycle arrest at G2/M. RPRM has been associated with malignant tumor progression and proposed as a potential biomarker for early cancer detection. However, the expression and role of RPRM, as well as its family, are poorly understood and their physiology is as yet unstudied. In this scenario, a model system like the zebrafish could serve to dissect the role of the RPRM family members in vivo. Phylogenetic analysis reveals that RPRM and RPRML have been differentially retained by most species throughout vertebrate evolution, yet RPRM3 has been retained only in a small group of distantly related species, including zebrafish. Herein, we characterized the spatiotemporal expression of RPRM (present in zebrafish as an infraclass duplication rprma/rprmb), RPRML and RPRM3 in the zebrafish. By whole-mount in situ hybridization (WISH) and fluorescent in situ hybridization (FISH), we demonstrate that rprm (rprma/rprmb) and rprml show a similar spatiotemporal expression profile during zebrafish development. At early developmental stages rprmb is expressed in somites. After one day post-fertilization, rprm (rprma/rprmb) and rprml are expressed in the notochord, brain, blood vessels and digestive tube. On the other hand, rprm3 shows the most unique expression profile, being expressed only in the central nervous system (CNS). We assessed the expression patterns of RPRM gene transcripts in adult zebrafish and human RPRM protein product in tissue samples by RT-qPCR and immunohistochemistry (IHC) staining, respectively. Strikingly, tissue-specific expression patterns of the RPRM transcripts and protein are conserved between zebrafish and humans. We propose the zebrafish as a powerful tool to elucidate the both physiological and pathological roles of the RPRM gene family.

Evolution of gremlin 2 in cetartiodactyl mammals: gene loss coincides with lack of upper jaw incisors in ruminants

Understanding the processes that give rise to genomic variability in extant species is an active area of research within evolutionary biology. With the availability of whole genome sequences, it is possible to quantify different forms of variability such as variation in gene copy number, which has been described as an important source of genetic variability and in consequence of phenotypic variability. Most of the research on this topic has been focused on understanding the biological significance of gene duplication, and less attention has been given to the evolutionary role of gene loss. Gremlin 2 is a member of the DAN gene family and plays a significant role in tooth development by blocking the ligand-signaling pathway of BMP2 and BMP4. The goal of this study was to investigate the evolutionary history of gremlin 2 in cetartiodactyl mammals, a group that possesses highly divergent teeth morphology. Results from our analyses indicate that gremlin 2 has experienced a mixture of gene loss, gene duplication, and rate acceleration. Although the last common ancestor of cetartiodactyls possessed a single gene copy, pigs and camels are the only cetartiodactyl groups that have retained gremlin 2. According to the phyletic distribution of this gene and synteny analyses, we propose that gremlin 2 was lost in the common ancestor of ruminants and cetaceans between 56.3 and 63.5 million years ago as a product of a chromosomal rearrangement. Our analyses also indicate that the rate of evolution of gremlin 2 has been accelerated in the two groups that have retained this gene. Additionally, the lack of this gene could explain the high diversity of teeth among cetartiodactyl mammals; specifically, the presence of this gene could act as a biological constraint. Thus, our results support the notions that gene loss is a way to increase phenotypic diversity and that gremlin 2 is a dispensable gene, at least in cetartiodactyl mammals.

Evolution of the β-adrenoreceptors in vertebrates

The study of the evolutionary history of genes related to human disease lies at the interface of evolution and medicine. These studies provide the evolutionary context on which medical researchers should work, and are also useful in providing information to suggest further genetic experiments, especially in model species where genetic manipulations can be made. Here we studied the evolution of the β-adrenoreceptor gene family in vertebrates with the aim of adding an evolutionary framework to the already abundant physiological information. Our results show that in addition to the three already described vertebrate β-adrenoreceptor genes there is an additional group containing cyclostome sequences. We suggest that β-adrenoreceptors diversified as a product of the two whole genome duplications that occurred in the ancestor of vertebrates. Gene expression patterns are in general consistent across species, suggesting that expression dynamics were established early in the evolutionary history of vertebrates, and have been maintained since then. Finally, amino acid polymorphisms that are associated to pathological conditions in humans appear to be common in non-human mammals, suggesting that the phenotypic effects of these mutations depend on epistatic interaction with other positions. The evolutionary analysis of the β-adrenoreceptors delivers new insights about the diversity of these receptors in vertebrates, the evolution of the expression patterns and a comparative perspective regarding the polymorphisms that in humans are linked to pathological conditions.

Evolutionary history of the reprimo tumor suppressor gene family in vertebrates with a description of a new reprimo gene lineage

Genes related to human diseases should be natural targets for evolutionary studies, since they could provide clues regarding the genetic bases of pathologies and potential treatments. Here we studied the evolution of the reprimo gene family, a group of tumor-suppressor genes that are implicated in p53-mediated cell cycle arrest. These genes, especially the reprimo duplicate located on human chromosome 2, have been associated with epigenetic modifications correlated with transcriptional silencing and cancer progression. We demonstrate the presence of a third reprimo lineage that, together with the reprimo and reprimo-like genes, appears to have been differentially retained during the evolutionary history of vertebrates. We present evidence that these reprimo lineages originated early in vertebrate evolution and expanded as a result of the two rounds of whole genome duplications that occurred in the last common ancestor of vertebrates. The reprimo gene has been lost in birds, and the third reprimo gene lineage has been retained in only a few distantly related species, such as coelacanth and gar. Expression analyses revealed that the reprimo paralogs are mainly expressed in the nervous system. Different vertebrate lineages have retained different reprimo paralogs, and even in species that have retained multiple copies, only one of them is heavily expressed.

Lineage-Specific Expansion of the Chalcone Synthase Gene Family in Rosids

Rosids are a monophyletic group that includes approximately 70,000 species in 140 families, and they are found in a variety of habitats and life forms. Many important crops such as fruit trees and legumes are rosids. The evolutionary success of this group may have been influenced by their ability to produce flavonoids, secondary metabolites that are synthetized through a branch of the phenylpropanoid pathway where chalcone synthase is a key enzyme. In this work, we studied the evolution of the chalcone synthase gene family in 12 species belonging to the rosid clade. Our results show that the last common ancestor of the rosid clade possessed six chalcone synthase gene lineages that were differentially retained during the evolutionary history of the group. In fact, of the six gene lineages that were present in the last common ancestor, 7 species retained 2 of them, whereas the other 5 only retained one gene lineage. We also show that one of the gene lineages was disproportionately expanded in species that belonged to the order Fabales (soybean, barrel medic and Lotus japonicas). Based on the available literature, we suggest that this gene lineage possesses stress-related biological functions (e.g., response to UV light, pathogen defense). We propose that the observed expansion of this clade was a result of a selective pressure to increase the amount of enzymes involved in the production of phenylpropanoid pathway-derived secondary metabolites, which is consistent with the hypothesis that suggested that lineage-specific expansions fuel plant adaptation.

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.

Three crocodilian genomes reveal ancestral patterns of evolution among archosaurs

To provide context for the diversification of archosaurs--the group that includes crocodilians, dinosaurs, and birds--we generated draft genomes of three crocodilians: Alligator mississippiensis (the American alligator), Crocodylus porosus (the saltwater crocodile), and Gavialis gangeticus (the Indian gharial). We observed an exceptionally slow rate of genome evolution within crocodilians at all levels, including nucleotide substitutions, indels, transposable element content and movement, gene family evolution, and chromosomal synteny. When placed within the context of related taxa including birds and turtles, this suggests that the common ancestor of all of these taxa also exhibited slow genome evolution and that the comparatively rapid evolution is derived in birds. The data also provided the opportunity to analyze heterozygosity in crocodilians, which indicates a likely reduction in population size for all three taxa through the Pleistocene. Finally, these data combined with newly published bird genomes allowed us to reconstruct the partial genome of the common ancestor of archosaurs, thereby providing a tool to investigate the genetic starting material of crocodilians, birds, and dinosaurs.

Comparative genomics reveals insights into avian genome evolution and adaptation

Birds are the most species-rich class of tetrapod vertebrates and have wide relevance across many research fields. We explored bird macroevolution using full genomes from 48 avian species representing all major extant clades. The avian genome is principally characterized by its constrained size, which predominantly arose because of lineage-specific erosion of repetitive elements, large segmental deletions, and gene loss. Avian genomes furthermore show a remarkably high degree of evolutionary stasis at the levels of nucleotide sequence, gene synteny, and chromosomal structure. Despite this pattern of conservation, we detected many non-neutral evolutionary changes in protein-coding genes and noncoding regions. These analyses reveal that pan-avian genomic diversity covaries with adaptations to different lifestyles and convergent evolution of traits.

Evolution of the relaxin/insulin-like gene family in anthropoid primates

The relaxin/insulin-like gene family includes signaling molecules that perform a variety of physiological roles mostly related to reproduction and neuroendocrine regulation. Several previous studies have focused on the evolutionary history of relaxin genes in anthropoid primates, with particular attention on resolving the duplication history of RLN1 and RLN2 genes, which are found as duplicates only in apes. These studies have revealed that the RLN1 and RLN2 paralogs in apes have a more complex history than their phyletic distribution would suggest. In this regard, alternative scenarios have been proposed to explain the timing of duplication, and the history of gene gain and loss along the organismal tree. In this article, we revisit the question and specifically reconstruct phylogenies based on coding and noncoding sequence in anthropoid primates to readdress the timing of the duplication event giving rise to RLN1 and RLN2 in apes. Results from our phylogenetic analyses based on noncoding sequence revealed that the duplication event that gave rise to the RLN1 and RLN2 occurred in the last common ancestor of catarrhine primates, between ∼44.2 and 29.6 Ma, and not in the last common ancestor of apes or anthropoids, as previously suggested. Comparative analyses based on coding and noncoding sequence suggests an event of convergent evolution at the sequence level between co-ortholog genes, the single-copy RLN gene found in New World monkeys and the RLN1 gene of apes, where changes in a fraction of the convergent sites appear to be driven by positive selection.

Increased rate of hair keratin gene loss in the cetacean lineage

BACKGROUND

Hair represents an evolutionary innovation that appeared early on mammalian evolutionary history, and presumably contributed significantly to the rapid radiation of the group. An interesting event in hair evolution has been its secondary loss in some mammalian groups, such as cetaceans, whose hairless phenotype appears to be an adaptive response to better meet the environmental conditions. To determine whether different repertoire of keratin genes among mammals can potentially explain the phenotypic hair features of different lineages, we characterized the type I and II clusters of alpha keratins from eight mammalian species, including the hairless dolphin and minke whale representing the order Cetacea.

RESULTS

We combined the available genomic information with phylogenetic analysis to conduct a comprehensive analysis of the evolutionary patterns of keratin gene clusters. We found that both type I and II gene clusters are fairly conserved among the terrestrial mammals included in this study, with lineage specific gene duplication and gene loss. Nevertheless, there is also evidence for an increased rate of pseudogenization in the cetacean lineage when compared to their terrestrial relatives, especially among the hair type keratins.

CONCLUSIONS

Here we present a comprehensive characterization of alpha-keratin genes among mammals and elucidate the mechanisms involved in the evolution of this gene family. We identified lineage-specific gene duplications and gene loss among the Laurasiatherian and Euarchontoglires species included in the study. Interestingly, cetaceans present an increased loss of hair-type keratin genes when compared to other terrestrial mammals. As suggested by the 'less-is-more' hypothesis, we do not rule out the possibility that the gene loss of hair-type keratin genes in these species might be associated to the hairless phenotype and could have been adaptive in response to new selective pressures imposed by the colonization of a new habitat. Our study provides support for the idea that pseudogenes are not simply 'genomic fossils' but instead have adaptive roles during the evolutionary process.

Genomic organization and differential signature of positive selection in the alpha and beta globin gene clusters in two cetacean species

The hemoglobin of jawed vertebrates is a heterotetramer protein that contains two α- and two β-chains, which are encoded by members of α- and β-globin gene families. Given the hemoglobin role in mediating an adaptive response to chronic hypoxia, it is likely that this molecule may have experienced a selective pressure during the evolution of cetaceans, which have to deal with hypoxia tolerance during prolonged diving. This selective pressure could have generated a complex history of gene turnover in these clusters and/or changes in protein structure themselves. Accordingly, we aimed to characterize the genomic organization of α- and β-globin gene clusters in two cetacean species and to detect a possible role of positive selection on them using a phylogenetic framework. Maximum likelihood and Bayesian phylogeny reconstructions revealed that both cetacean species had retained a similar complement of putatively functional genes. For the α-globin gene cluster, the killer whale presents a complement of genes composed of HBZ, HBK, and two functional copies of HBA and HBQ genes, whereas the dolphin possesses HBZ, HBK, HBA and HBQ genes, and one HBA pseudogene. For the β-globin gene cluster, both species retained a complement of four genes, two early expressed genes—HBE and HBH—and two adult expressed genes—HBD and HBB. Our natural selection analysis detected two positively selected sites in the HBB gene (56 and 62) and four in HBA (15, 21, 49, 120). Interestingly, only the genes that are expressed during the adulthood showed the signature of positive selection.

Characterization of the kidney transcriptome of the South American olive mouse Abrothrix olivacea

BACKGROUND

The olive mouse Abrothrix olivacea is a cricetid rodent of the subfamily Sigmodontinae that inhabits a wide range of contrasting environments in southern South America, from aridlands to temperate rainforests. Along its distribution, it presents different geographic forms that make the olive mouse a good focal case for the study of geographical variation in response to environmental variation. We chose to characterize the kidney transcriptome because this organ has been shown to be associated with multiple physiological processes, including water reabsorption.

RESULTS

Transcriptomes of thirteen kidneys from individuals from Argentina and Chile were sequenced using Illumina technology in order to obtain a kidney reference transcriptome. After combining the reads produced for each sample, we explored three assembly strategies to obtain the best reconstruction of transcripts, TrinityNorm and DigiNorm, which include its own normalization algorithms for redundant reads removal, and Multireads, which simply consist on the assembly of the joined reads. We found that Multireads strategy produces a less fragmented assembly than normalization algorithms but recovers fewer number of genes. In general, about 15000 genes were annotated, of which almost half had at least one coding sequence reconstructed at 99% of its length. We also built a list of highly expressed genes, of which several are involved in water conservation under laboratory conditions using mouse models.

CONCLUSION

Based on our assembly results, Trinity's in silico normalization is the best algorithm in terms of cost-benefit returns; however, our results also indicate that normalization should be avoided if complete or nearly complete coding sequences of genes are desired. Given that this work is the first to characterize the transcriptome of any member of Sigmodontinae, a subfamily of cricetid rodents with about 400 living species, it will provide valuable resources for future ecological and evolutionary genomic analyses.

The complete mitochondrial genome of the land snail Cornu aspersum (Helicidae: Mollusca): intra-specific divergence of protein-coding genes and phylogenetic considerations within Euthyneura

The complete sequences of three mitochondrial genomes from the land snail Cornu aspersum were determined. The mitogenome has a length of 14050 bp, and it encodes 13 protein-coding genes, 22 transfer RNA genes and two ribosomal RNA genes. It also includes nine small intergene spacers, and a large AT-rich intergenic spacer. The intra-specific divergence analysis revealed that COX1 has the lower genetic differentiation, while the most divergent genes were NADH1, NADH3 and NADH4. With the exception of Euhadra herklotsi, the structural comparisons showed the same gene order within the family Helicidae, and nearly identical gene organization to that found in order Pulmonata. Phylogenetic reconstruction recovered Basommatophora as polyphyletic group, whereas Eupulmonata and Pulmonata as paraphyletic groups. Bayesian and Maximum Likelihood analyses showed that C. aspersum is a close relative of Cepaea nemoralis, and with the other Helicidae species form a sister group of Albinaria caerulea, supporting the monophyly of the Stylommatophora clade.

How to Make a Dolphin: Molecular Signature of Positive Selection in Cetacean Genome

Cetaceans are unique in being the only mammals completely adapted to an aquatic environment. This adaptation has required complex changes and sometimes a complete restructuring of physiology, behavior and morphology. Identifying genes that have been subjected to selection pressure during cetacean evolution would greatly enhance our knowledge of the ways in which genetic variation in this mammalian order has been shaped by natural selection. Here, we performed a genome-wide scan for positive selection in the dolphin lineage. We employed models of codon substitution that account for variation of selective pressure over branches on the tree and across sites in a sequence. We analyzed 7,859 nuclear-coding ortholog genes and using a series of likelihood ratio tests (LRTs), we identified 376 genes (4.8%) with molecular signatures of positive selection in the dolphin lineage. We used the cow as the sister group and compared estimates of selection in the cetacean genome to this using the same methods. This allowed us to define which genes have been exclusively under positive selection in the dolphin lineage. The enrichment analysis found that the identified positively selected genes are significantly over-represented for three exclusive functional categories only in the dolphin lineage: segment specification, mesoderm development and system development. Of particular interest for cetacean adaptation to an aquatic life are the following GeneOntology targets under positive selection: genes related to kidney, heart, lung, eye, ear and nervous system development.

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

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

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

The functional diversification of the vertebrate globin gene superfamily provides an especially vivid illustration of the role of gene duplication and whole-genome duplication in promoting evolutionary innovation. For example, key globin proteins that evolved specialized functions in various aspects of oxidative metabolism and oxygen signaling pathways (hemoglobin [Hb], myoglobin [Mb], and cytoglobin [Cygb]) trace their origins to two whole-genome duplication events in the stem lineage of vertebrates. The retention of the proto-Hb and Mb genes in the ancestor of jawed vertebrates permitted a physiological division of labor between the oxygen-carrier function of Hb and the oxygen-storage function of Mb. In the Hb gene lineage, a subsequent tandem gene duplication gave rise to the proto α- and β-globin genes, which permitted the formation of multimeric Hbs composed of unlike subunits (α(2)β(2)). The evolution of this heteromeric quaternary structure was central to the emergence of Hb as a specialized oxygen-transport protein because it provided a mechanism for cooperative oxygen-binding and allosteric regulatory control. Subsequent rounds of duplication and divergence have produced diverse repertoires of α- and β-like globin genes that are ontogenetically regulated such that functionally distinct Hb isoforms are expressed during different stages of prenatal development and postnatal life. In the ancestor of jawless fishes, the proto Mb and Hb genes appear to have been secondarily lost, and the Cygb homolog evolved a specialized respiratory function in blood-oxygen transport. Phylogenetic and comparative genomic analyses of the vertebrate globin gene superfamily have revealed numerous instances in which paralogous globins have convergently evolved similar expression patterns and/or similar functional specializations in different organismal lineages.

Evolution of the globin gene family in deuterostomes: lineage-specific patterns of diversification and attrition

In the Metazoa, globin proteins display an underlying unity in tertiary structure that belies an extraordinary diversity in primary structures, biochemical properties, and physiological functions. Phylogenetic reconstructions can reveal which of these functions represent novel, lineage-specific innovations, and which represent ancestral functions that are shared with homologous globin proteins in other eukaryotes and even prokaryotes. To date, our understanding of globin diversity in deuterostomes has been hindered by a dearth of genomic sequence data from the Ambulacraria (echinoderms + hemichordates), the sister group of chordates, and the phylum Xenacoelomorpha, which includes xenoturbellids, acoelomorphs, and nemertodermatids. Here, we report the results of a phylogenetic and comparative genomic analysis of the globin gene repertoire of deuterostomes. We first characterized the globin genes of the acorn worm, Saccoglossus kowalevskii, a representative of the phylum Hemichordata. We then integrated genomic sequence data from the acorn worm into a comprehensive analysis of conserved synteny and phylogenetic relationships among globin genes from representatives of the eight lineages that comprise the superphylum Deuterostomia. The primary aims were 1) to unravel the evolutionary history of the globin gene superfamily in deuterostomes and 2) to use the estimated phylogeny to gain insights into the functional evolution of deuterostome globins. Results of our analyses indicate that the deuterostome common ancestor possessed a repertoire of at least four distinct globin paralogs and that different subsets of these ancestral genes have been retained in each of the descendant organismal lineages. In each major deuterostome group, a different subset of ancestral precursor genes underwent lineage-specific expansions of functional diversity through repeated rounds of gene duplication and divergence. By integrating results of the phylogenetic analysis with available functional data, we discovered that circulating oxygen-transport hemoglobins evolved independently in several deuterostome lineages and that intracellular nerve globins evolved independently in chordates and acoelomorph worms.

Gene turnover and differential retention in the relaxin/insulin-like gene family in primates

The relaxin/insulin-like gene family is related to the insulin gene family, and includes two separate types of peptides: relaxins (RLNs) and insulin-like peptides (INSLs) that perform a variety of physiological roles including testicular descent, growth and differentiation of the mammary glands, trophoblast development, and cell differentiation. In vertebrates, these genes are found on three separate genomic loci, and in mammals, variation in the number and nature of genes in this family is mostly restricted to the Relaxin Family Locus B. For example, this locus contains a single copy of RLN in platypus and opossum, whereas it contains copies of the INSL6, INSL4, RLN2 and RLN1 genes in human and chimp. The main objective of this research is to characterize changes in the size and membership composition of the RLN/INSL gene family in primates, reconstruct the history of the RLN/INSL genes of primates, and test competing evolutionary scenarios regarding the origin of INSL4 and of the duplicated copies of the RLN gene of apes. Our results show that the relaxin/INSL-like gene family of primates has had a more dynamic evolutionary history than previously thought, including several examples of gene duplications and losses which are consistent with the predictions of the birth-and-death model of gene family evolution. In particular, we found that the differential retention of relatively old paralogs played a key role in shaping the gene complement of this family in primates. Two examples of this phenomenon are the origin of the INSL4 gene of catarrhines (the group that includes Old World monkeys and apes), and of the duplicate RLN1 and RLN2 paralogs of apes. In the case of INSL4, comparative genomics and phylogenetic analyses indicate that the origin of this gene, which was thought to represent a catarrhine-specific evolutionary innovation, is as old as the split between carnivores and primates, which took place approximately 97 million years ago. In addition, in the case of the RLN1 and RLN2 genes of apes our phylogenetic trees and topology tests indicate that the duplication that gave rise to these two genes maps to the last common ancestor of anthropoid primates. All these genomic changes in gene complement, which are particularly prevalent among anthropoid primates, might be linked to the many physiological and anatomical changes found in this group. Given the various roles of members of the RLN/INSL-like gene family in reproductive biology, it might be that changes in this gene family are associated to changes in reproductive traits.

Silencing, positive selection and parallel evolution: busy history of primate cytochromes C

Cytochrome c (cyt c) participates in two crucial cellular processes, energy production and apoptosis, and unsurprisingly is a highly conserved protein. However, previous studies have reported for the primate lineage (i) loss of the paralogous testis isoform, (ii) an acceleration and then a deceleration of the amino acid replacement rate of the cyt c somatic isoform, and (iii) atypical biochemical behavior of human cyt c. To gain insight into the cause of these major evolutionary events, we have retraced the history of cyt c loci among primates. For testis cyt c, all primate sequences examined carry the same nonsense mutation, which suggests that silencing occurred before the primates diversified. For somatic cyt c, maximum parsimony, maximum likelihood, and Bayesian phylogenetic analyses yielded the same tree topology. The evolutionary analyses show that a fast accumulation of non-synonymous mutations (suggesting positive selection) occurred specifically on the anthropoid lineage root and then continued in parallel on the early catarrhini and platyrrhini stems. Analysis of evolutionary changes using the 3D structure suggests they are focused on the respiratory chain rather than on apoptosis or other cyt c functions. In agreement with previous biochemical studies, our results suggest that silencing of the cyt c testis isoform could be linked with the decrease of primate reproduction rate. Finally, the evolution of cyt c in the two sister anthropoid groups leads us to propose that somatic cyt c evolution may be related both to COX evolution and to the convergent brain and body mass enlargement in these two anthropoid clades.

Whole-genome duplications spurred the functional diversification of the globin gene superfamily in vertebrates

It has been hypothesized that two successive rounds of whole-genome duplication (WGD) in the stem lineage of vertebrates provided genetic raw materials for the evolutionary innovation of many vertebrate-specific features. However, it has seldom been possible to trace such innovations to specific functional differences between paralogous gene products that derive from a WGD event. Here, we report genomic evidence for a direct link between WGD and key physiological innovations in the vertebrate oxygen transport system. Specifically, we demonstrate that key globin proteins that evolved specialized functions in different aspects of oxidative metabolism (hemoglobin, myoglobin, and cytoglobin) represent paralogous products of two WGD events in the vertebrate common ancestor. Analysis of conserved macrosynteny between the genomes of vertebrates and amphioxus (subphylum Cephalochordata) revealed that homologous chromosomal segments defined by myoglobin + globin-E, cytoglobin, and the α-globin gene cluster each descend from the same linkage group in the reconstructed proto-karyotype of the chordate common ancestor. The physiological division of labor between the oxygen transport function of hemoglobin and the oxygen storage function of myoglobin played a pivotal role in the evolution of aerobic energy metabolism, supporting the hypothesis that WGDs helped fuel key innovations in vertebrate evolution.

Differential loss and retention of cytoglobin, myoglobin, and globin-E during the radiation of vertebrates

If rates of postduplication gene retention are positively correlated with levels of functional constraint, then gene duplicates that have been retained in a restricted number of taxonomic lineages would be expected to exhibit relatively low levels of sequence conservation. Paradoxical patterns are presented by gene duplicates that have been retained in a small number of taxa but which are nonetheless subject to strong purifying selection relative to paralogous members of the same multigene family. This pattern suggests that such genes may have been co-opted for novel, lineage-specific functions. One possible example involves the enigmatic globin-E gene (GbE), which appears to be exclusively restricted to birds. Available data indicate that this gene is expressed exclusively in the avian eye, but its physiological function remains a mystery. In contrast to the highly restricted phyletic distribution of GbE, the overwhelming majority of jawed vertebrates (gnathostomes) possess copies of the related cytoglobin (Cygb) and myoglobin (Mb) genes. The purpose of the present study was 1) to assess the phyletic distribution of the Cygb, Mb, and GbE genes among vertebrates, 2) to elucidate the duplicative origins and evolutionary histories of these three genes, and 3) to evaluate the relative levels of functional constraint of these genes based on comparative sequence analysis. To accomplish these objectives, we conducted a combined phylogenetic and comparative genomic analysis involving taxa that represent each of the major lineages of gnathostome vertebrates. Results of synteny comparisons and phylogenetic topology tests revealed that GbE is clearly not the product of a recent, bird-specific duplication event. Instead, GbE originated via duplication of a proto-Mb gene in the stem lineage of gnathostomes. Unlike the Mb gene, which has been retained in all major gnathostome lineages other than amphibians, the GbE gene has been retained only in the lineage leading to modern birds and has been independently lost in at least four major lineages: teleost fish, amphibians, mammals, and nonavian reptiles. Despite the restricted phyletic distribution of this gene, our results indicate that GbE is one of the most highly conserved globins in the avian genome.