Breathing Air and Living Underwater: Molecular Evolution of Genes Related to Antioxidant Response in Cetaceans and Pinnipeds

Cetaceans and pinnipeds are lineages of mammals that have independently returned to the aquatic environment, acquiring varying degrees of dependence on it while sharing adaptations for underwater living. Here, we focused on one critical adaptation from both groups, their ability to withstand the ischemia and reperfusion experienced during apnea diving, which can lead to the production of reactive oxygen species (ROS) and subsequent oxidative damage. Previous studies have shown that cetaceans and pinnipeds possess efficient antioxidant enzymes that protect against ROS. In this study, we investigated the molecular evolution of key antioxidant enzyme genes (CAT, GPX3, GSR, PRDX1, PRDX3, and SOD1) and the ROS-producing gene XDH, in cetaceans and pinnipeds lineages. We used the ratio of non-synonymous (dN) to synonymous (dS) substitutions as a measure to identify signatures of adaptive molecular evolution in these genes within and between the two lineages. Additionally, we performed protein modeling and variant impact analyzes to assess the functional consequences of observed mutations. Our findings revealed distinct selective regimes between aquatic and terrestrial mammals in five of the examined genes, including divergences within cetacean and pinniped lineages, between ancestral and recent lineages and between crowns groups. We identified specific sites under positive selection unique to Cetacea and Pinnipedia, with one site showing evidence of convergent evolution in species known for their long and deep-diving capacities. Notably, many sites under adaptive selection exhibited radical changes in amino acid properties, with some being damaging mutations in human variations, but with no apparent detrimental impacts on aquatic mammals. In conclusion, our study provides insights into the adaptive changes that have occurred in the antioxidant systems of aquatic mammals throughout their evolutionary history. We observed both distinctive features within each group of Cetacea and Pinnipedia and instances of convergence. These findings highlight the dynamic nature of the antioxidant system in response to challenges of the aquatic environment and provide a foundation for further investigations into the molecular mechanisms underlying these adaptations.

Molecular Footprints on Osmoregulation-Related Genes Associated with Freshwater Colonization by Cetaceans and Sirenians

The genetic basis underlying adaptive physiological mechanisms has been extensively explored in mammals after colonizing the seas. However, independent lineages of aquatic mammals exhibit complex patterns of secondary colonization in freshwater environments. This change in habitat represents new osmotic challenges, and additional changes in key systems, such as the osmoregulatory system, are expected. Here, we studied the selective regime on coding and regulatory regions of 20 genes related to the osmoregulation system in strict aquatic mammals from independent evolutionary lineages, cetaceans, and sirenians, with representatives in marine and freshwater aquatic environments. We identified positive selection signals in genes encoding the protein vasopressin (AVP) in mammalian lineages with secondary colonization in the fluvial environment and in aquaporins for lineages inhabiting the marine and fluvial environments. A greater number of sites with positive selection signals were found for the dolphin species compared to the Amazonian manatee. Only the AQP5 and AVP genes showed selection signals in more than one independent lineage of these mammals. Furthermore, the vasopressin gene tree indicates greater similarity in river dolphin sequences despite the independence of their lineages based on the species tree. Patterns of distribution and enrichment of Transcription Factors in the promoter regions of target genes were analyzed and appear to be phylogenetically conserved among sister species. We found accelerated evolution signs in genes ACE, AQP1, AQP5, AQP7, AVP, NPP4, and NPR1 for the fluvial mammals. Together, these results allow a greater understanding of the molecular bases of the evolution of genes responsible for osmotic control in aquatic mammals.

Patterns of enrichment and acceleration in evolutionary rates of promoters suggest a role of regulatory regions in cetacean gigantism

BACKGROUND

Cetaceans (whales, porpoises, and dolphins) are a lineage of aquatic mammals from which some species became giants. Only recently, gigantism has been investigated from the molecular point of view. Studies focused mainly on coding regions, and no data on the influence of regulatory regions on gigantism in this group was available. Accordingly, we investigated the molecular evolution of non-coding regulatory regions of genes already described in the literature for association with size in mammals, focusing mainly on the promoter regions. For this, we used Ciiider and phyloP tools. Ciiider identifies significantly enriched transcription factor binding sites, and phyloP estimates the molecular evolution rate of the promoter.

RESULTS

We found evidence of enrichment of transcription binding factors related to large body size, with distinct patterns between giant and non-giant cetaceans in the IGFBP7 and NCAPG promoters, in which repressive agents are present in small cetaceans and those that stimulate transcription, in giant cetaceans. In addition, we found evidence of acceleration in the IGF2, IGFBP2, IGFBP7, and ZFAT promoters.

CONCLUSION

Our results indicate that regulatory regions may also influence cetaceans' body size, providing candidate genes for future research to understand the molecular basis of the largest living animals.

The molecular evolution of genes previously associated with large sizes reveals possible pathways to cetacean gigantism

Cetaceans are a group of aquatic mammals with the largest body sizes among living animals, including giant representatives such as blue and fin whales. To understand the genetic bases of gigantism in cetaceans, we performed molecular evolutionary analyses on five genes (GHSR, IGF2, IGFBP2, IGFBP7, and EGF) from the growth hormone/insulin-like growth factor axis, and four genes (ZFAT, EGF, LCORL, and PLAG1) previously described as related to the size of species evolutionarily close to cetaceans, such as pigs, cows, and sheep. Our dataset comprised 19 species of cetaceans, seven of which are classified as giants because they exceed 10 m in length. Our results revealed signs of positive selection in genes from the growth hormone/insulin-like growth factor axis and also in those related to body increase in cetacean-related species. In addition, pseudogenization of the EGF gene was detected in the lineage of toothless cetaceans, Mysticeti. Our results suggest the action of positive selection on gigantism in genes that act both in body augmentation and in mitigating its consequences, such as cancer suppression when involved in processes such as division, migration, and cell development control.

Positive selection in multiple salivary gland proteins of Anophelinae reveals potential targets for vector control

Anopheles is a genus belonging to the Culicidae family, which has great medical importance due to its role as a vector of Plasmodium, the causative agent of malaria. Great focus has been given to the salivary gland proteins (SGPs) group from Anopheles' functional genomics. This class of proteins is essential to blood-feeding behavior as they have attributes such as vasodilators and anti-clotting properties. Recently, a comprehensive review on Anopheles SGPs was performed; however, the authors did not deeply explore the adaptive molecular evolution of these genes. In this context, this work aimed to perform a more detailed analysis of the adaptive molecular evolution of SGPs in Anopheles, carrying out positive selection and gene family evolution analysis on 824 SGPs. Our results show that most SGPs have positively selected codon sites that can be used as targets in developing new strategies for vector control and that younger SGPs evolve at a faster rate than older SGPs. Notably, we could not find any evidence of an accelerated shift in SGPs' rates of gene gain and loss compared with other proteins, as suggested in previous works.

Expansions and contractions in gene families of independently-evolved blood-feeding insects

BACKGROUND

The blood-feeding behavior evolved multiple times in Insecta lineages and it represents an excellent opportunity to study patterns of convergent molecular evolution regarding this habit. In insects the expansion of some gene families is linked with blood-feeding behavior, but a wide study comparing the evolution of these gene families among different lineages is still missing. Here we gathered genomic data from six independently-evolved hematophagous lineages, aiming to identify convergent expansions and/or contractions of gene families in hematophagous lineages of insects.

RESULTS

We found four rapidly evolving gene families shared by at least two hematophagous independently-evolved lineages, including a heat-shock and a chemosensory protein. On the expression of these four rapidly evolving gene families we found more genes expressed in mated individuals compared with virgin individuals in rapidly-expanded families and more genes expressed in non-blood-feeding individuals compared with blood-feeding individuals in rapidly-contracted families.

CONCLUSION

Our results reveal a new set of candidate genes to be explored in further analysis to help the development of new strategies to deal with blood-feeding vectors and also presents a new perspective to study the evolution of hematophagy identifying convergent molecular patterns.

Gene fusion of heterophyletic gamma-globin genes in platyrrhine primates

We performed phylogenetic analyses of HBG genes to assess its origin and interspecific variation among primates. Our analyses showed variation in HBG genes copy number ranging from one to three, some of them pseudogenes. For platyrrhines HBG genes, phylogenetic reconstructions of flanking regions recovered orthologous clades with distinct topologies for 5' and 3' flanking regions. The 5' region originated in the common ancestor of platyrrhines but the 3' region had an anthropoid origin. We hypothesize that the platyrrhine HBG genes of 5' and 3' heterophyletic origins arose from subsequent fusions of the (earlier) platyrrhine 5' portion and the (later) anthropoid 3' portion.

Molecular evolution of HoxA13 and the multiple origins of limbless morphologies in amphibians and reptiles

Developmental processes and their results, morphological characters, are inherited through transmission of genes regulating development. While there is ample evidence that cis-regulatory elements tend to be modular, with sequence segments dedicated to different roles, the situation for proteins is less clear, being particularly complex for transcription factors with multiple functions. Some motifs mediating protein-protein interactions may be exclusive to particular developmental roles, but it is also possible that motifs are mostly shared among different processes. Here we focus on HoxA13, a protein essential for limb development. We asked whether the HoxA13 amino acid sequence evolved similarly in three limbless clades: Gymnophiona, Amphisbaenia and Serpentes. We explored variation in ω (dN/dS) using a maximum-likelihood framework and HoxA13sequences from 47 species. Comparisons of evolutionary models provided low ω global values and no evidence that HoxA13 experienced relaxed selection in limbless clades. Branch-site models failed to detect evidence for positive selection acting on any site along branches of Amphisbaena and Gymnophiona, while three sites were identified in Serpentes. Examination of alignments did not reveal consistent sequence differences between limbed and limbless species. We conclude that HoxA13 has no modules exclusive to limb development, which may be explained by its involvement in multiple developmental processes.

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.

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.

Site fidelity of Sotalia guianensis (Cetacea: Delphinidae) in Sepetiba Bay, Rio de Janeiro, Brazil

Estuarine dolphins, Sotalia guianensis (Van Bénéden, 1864), were identified in Sepetiba Bay (22º35'S, 44º03'W), state of Rio de Janeiro, Brazil, between February 1994 and July 2006 using the photo-identification method. A total of 72 surveys resulted in the identification and cataloguing of 217 individual dolphins. Seventy-nine (36%) were resighted at least once during the study period. The average interval between resightings was 40,6 ± 31,1 months. Large numbers of mother-and-calf pairs were recorded during the 12 years and new individuals were frequently encountered in the photographic records. These results suggest that part of the dolphin population is more consistently found in Sepetiba Bay, probably for feeding and reproduction, while other individuals use the Bay only sporadically.