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

Evolution of ion channels in cetaceans: a natural experiment in the tree of life

Cetaceans represent a natural experiment within the tree of life in which a lineage changed from terrestrial to aquatic habitats. This shift involved phenotypic modifications, representing an opportunity to explore the genetic bases of phenotypic diversity. Among the different molecular systems that maintain cellular homeostasis, ion channels are crucial for the proper physiological functioning of all living species. This study aims to explore the evolution of ion channels during the evolutionary history of cetaceans. To do so, we created a bioinformatic pipeline to annotate the repertoire of ion channels in the genome of the species included in our sampling. Our main results show that cetaceans have, on average, fewer protein-coding genes and a higher percentage of annotated ion channels than non-cetacean mammals. Signals of positive selection were detected in ion channels related to the heart, locomotion, visual and neurological phenotypes. Interestingly, we predict that the NaV1.5 ion channel of most toothed whales (odontocetes) is sensitive to tetrodotoxin, similar to NaV1.7, given the presence of tyrosine instead of cysteine, in a specific position of the ion channel. Finally, the gene turnover rate of the cetacean crown group is more than three times faster than that of non-cetacean mammals.

Contraction of the ROS Scavenging Enzyme Glutathione S-Transferase Gene Family in Cetaceans

Cetaceans are a group of marine mammals whose ancestors were adaptated for life on land. Life in an aquatic environment poses many challenges for air-breathing mammals. Diving marine mammals have adapted to rapid reoxygenation and reactive oxygen species (ROS)-mediated reperfusion injury. Here, we considered the evolution of the glutathione transferase (GST) gene family which has important roles in the detoxification of endogenously-derived ROS and environmental pollutants. We characterized the cytosolic GST gene family in 21 mammalian species; cetaceans, sirenians, pinnipeds, and their terrestrial relatives. All seven GST classes were identified, showing that GSTs are ubiquitous in mammals. Some GST genes are the product of lineage-specific duplications and losses, in line with a birth-and-death evolutionary model. We detected sites with signatures of positive selection that possibly influence GST structure and function, suggesting that adaptive evolution of GST genes is important for defending mammals from various types of noxious environmental compounds. We also found evidence for loss of alpha and mu GST subclass genes in cetacean lineages. Notably, cetaceans have retained a homolog of at least one of the genes GSTA1, GSTA4, and GSTM1; GSTs that are present in both the cytosol and mitochondria. The observed variation in number and selection pressure on GST genes suggest that the gene family structure is dynamic within cetaceans.

Beluga whale pVHL enhances HIF-2α activity via inducing HIF-2α proteasomal degradation under hypoxia

Aquatic mammals, such as cetaceans experience various depths, with accordingly diverse oxygenation, thus, cetaceans have developed adaptations for hypoxia, but mechanisms underlying this tolerance to low oxygen are unclear. Here we analyzed VHL and HIF-2α, in the hypoxia signaling pathway. Variations in VHL are greater than HIF-2α between cetaceans and terrestrial mammals, and beluga whale VHL (BW-VHL) promotes HIF-2α degradation under hypoxia. BW-VHL catalyzes BW-HIF-2α to form K48-linked poly-ubiquitin chains mainly at the lysine 429 of BW-HIF-2α (K429) and induces BW-HIF-2α for proteasomal degradation. W100 within BW-VHL is a key site for BW-VHL functionally and BW-VHL enhances transcriptional activity of BW-HIF-2α under hypoxia. Our data therefore reveal that BW-VHL has a unique function that may contribute to hypoxic adaptation.

Selection on different genes with equivalent functions: the convergence story told by Hox genes along the evolution of aquatic mammalian lineages

Background

Convergent evolution has been a challenging topic for decades, being cetaceans, pinnipeds and sirenians textbook examples of three independent origins of equivalent phenotypes. These mammalian lineages acquired similar anatomical features correlated to an aquatic life, and remarkably differ from their terrestrial counterparts. Whether their molecular evolutionary history also involved similar genetic mechanisms underlying such morphological convergence nevertheless remained unknown. To test for the existence of convergent molecular signatures, we studied the molecular evolution of Hox genes in these three aquatic mammalian lineages, comparing their patterns to terrestrial mammals. Hox genes are transcription factors that play a pivotal role in specifying embryonic regional identity of nearly any bilateral animal, and are recognized major agents for diversification of body plans.

Results

We detected few signatures of positive selection on Hox genes across the three aquatic mammalian lineages and verified that purifying selection prevails in these sequences, as expected for pleiotropic genes. Genes found as being positively selected differ across the aquatic mammalian lineages, but we identified a substantial overlap of their developmental functions. Such pattern likely resides on the duplication history of Hox genes, which probably provided different possible evolutionary routes for achieving the same phenotypic solution.

Conclusions

Our results indicate that convergence occurred at a functional level of Hox genes along three independent origins of aquatic mammals. This conclusion reinforces the idea that different changes in developmental genes may lead to similar phenotypes, probably due to the redundancy provided by the participation of Hox paralogous genes in several developmental functions.

Electronic supplementary material

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

Evolutionary Genetics of Hypoxia Tolerance in Cetaceans during Diving

Hypoxia was a major challenge faced by cetaceans during the course of secondary aquatic adaptation. Although physiological traits of hypoxia tolerance in cetaceans have been well characterized, the underlying molecular mechanisms remain unknown. We investigated the sequences of 17 hypoxia-tolerance-related genes in representative cetaceans to provide a comprehensive insight into the genetic basis of hypoxia tolerance in these animals. Genes involved in carrying and transporting oxygen in the blood and muscle (hemoglobin-α and β, myoglobin), and genes involved in the regulation of vasoconstriction (endothelin-1, -2, and -3; endothelin receptor type A and B; adrenergic receptor α-1D; and arginine vasopressin) appear to have undergone adaptive evolution, evidence for positive selection on their particular sites, and radical physiochemical property changes of selected condons. Interestingly, “long-diving” cetaceans had relatively higher ω (d_N/d_S) values than “short-diving” cetaceans for the hemoglobin β gene, indicating divergent selective pressure presented in cetacean lineages with different diving abilities. Additionally, parallel positive selection or amino acid changes (ADRA1D: P50A, A53G, AVPR1B: I/V270T) among animals exposed to different hypoxia habitats reflect functional convergence or similar genetic mechanisms of hypoxia tolerance. In summary, positive selection, divergent selective pressures, and parallel evolution at the molecular level provided some new insights into the genetic adaptation of hypoxia tolerance.

Beta-globin gene evolution in the ruminants: evidence for an ancient origin of sheep haplotype B

Domestic sheep (Ovis aries) can be divided into two groups with significantly different responses to hypoxic environments, determined by two allelic beta-globin haplotypes. Haplotype A is very similar to the goat beta-globin locus, whereas haplotype B has a deletion spanning four globin genes, including beta-C globin, which encodes a globin with high oxygen affinity. We surveyed the beta-globin locus using resequencing data from 70 domestic sheep from 42 worldwide breeds and three Ovis canadensis and two Ovis dalli individuals. Haplotype B has an allele frequency of 71.4% in O. aries and was homozygous (BB) in all five wild sheep. This shared ancestry indicates haplotype B is at least 2-3 million years old. Approximately 40 kb of the sequence flanking the ~37-kb haplotype B deletion had unexpectedly low identity between haplotypes A and B. Phylogenetic analysis showed that the divergent region of sheep haplotype B is remarkably distinct from the beta-globin loci in goat and cattle but still groups with the Ruminantia. We hypothesize that this divergent ~40-kb region in haplotype B may be from an unknown ancestral ruminant and was maintained in the lineage to O. aries, but not other Bovidae, evolving independently of haplotype A. Alternatively, the ~40-kb sequence in haplotype B was more recently acquired by an ancestor of sheep from an unknown non-Bovidae ruminant, replacing part of haplotype A. Haplotype B has a lower nucleotide diversity than does haplotype A, suggesting a recent bottleneck, whereas the higher frequency of haplotype B suggests a subsequent spread through the global population of O. aries.

Myoglobin oxygen affinity in aquatic and terrestrial birds and mammals

Myoglobin (Mb) is an oxygen binding protein found in vertebrate skeletal muscle, where it facilitates intracellular transport and storage of oxygen. This protein has evolved to suit unique physiological needs in the muscle of diving vertebrates that express Mb at much greater concentrations than their terrestrial counterparts. In this study, we characterized Mb oxygen affinity (P50) from 25 species of aquatic and terrestrial birds and mammals. Among diving species, we tested for correlations between Mb P50 and routine dive duration. Across all species examined, Mb P50 ranged from 2.40 to 4.85 mmHg. The mean P50 of Mb from terrestrial ungulates was 3.72±0.15 mmHg (range 3.70-3.74 mmHg). The P50 of cetaceans was similar to terrestrial ungulates ranging from 3.54 to 3.82 mmHg, with the exception of the melon-headed whale, which had a significantly higher P50 of 4.85 mmHg. Among pinnipeds, the P50 ranged from 3.23 to 3.81 mmHg and showed a trend for higher oxygen affinity in species with longer dive durations. Among diving birds, the P50 ranged from 2.40 to 3.36 mmHg and also showed a trend of higher affinities in species with longer dive durations. In pinnipeds and birds, low Mb P50 was associated with species whose muscles are metabolically active under hypoxic conditions associated with aerobic dives. Given the broad range of potential globin oxygen affinities, Mb P50 from diverse vertebrate species appears constrained within a relatively narrow range. High Mb oxygen affinity within this range may be adaptive for some vertebrates that make prolonged dives.