Next Article in Journal
The Photoperiod Regulates Granulosa Cell Apoptosis through the FSH-Nodal/ALK7 Signaling Pathway in Phodopus sungorus
Next Article in Special Issue
Comparative Analyses Reveal the Genetic Mechanism of Ambergris Production in the Sperm Whale Based on the Chromosome-Level Genome
Previous Article in Journal
Informing Implementation: Perspectives from the Australian University Community Regarding an Animal Assisted Intervention
Previous Article in Special Issue
Complete Mitochondrial Genome, Genetic Diversity and Phylogenetic Analysis of Pingpu Yellow Chicken (Gallus gallus)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 12
Issue 24
10.3390/ani12243571
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rubbing Salt in the Wound: Molecular Evolutionary Analysis of Pain-Related Genes Reveals the Pain Adaptation of Cetaceans in Seawater

by
Xiaoyue Ding
†,
Fangfang Yu
†,
Xiaofang He
,
Shixia Xu
,
Guang Yang
and
Wenhua Ren
*
Jiangsu Key Laboratory for Biodiversity and Biotechnology, College of Life Sciences, Nanjing Normal University, Nanjing 210000, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Animals 2022, 12(24), 3571; https://doi.org/10.3390/ani12243571
Submission received: 25 October 2022 / Revised: 26 November 2022 / Accepted: 14 December 2022 / Published: 16 December 2022
(This article belongs to the Special Issue Adaptive Evolution and Trait Formation of Animals)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Cetaceans are aquatic mammals that evolved a series of specializations for life in an aquatic habitat, such as lack of distal hindlimbs, loss of hair, and derivation of echolocation. Notably, in the face of high salinity of seawater, the molecular mechanism of adaptation to pain in cetacean is still unclear. In this study, we performed molecular evolutionary analyses of genes related to pain perception (pain) and pain relief (analgesia) in selected representatives of mammals to explore the molecular mechanisms of pain adaptation in cetaceans to ‘rubbing salt in the wound’. Relaxed selection, positive selection, and convergent and specific amino acid substitutions were identified within cetacean’s pain-related genes, showing that the adaptation of mammals to a seawater environment might also include molecular evolution towards greater sensitivity to pain and more effective analgesia. Our study could have implications for diagnosis and treatment of human pain.

Abstract

Pain, usually caused by a strong or disruptive stimulus, is an unpleasant sensation that serves as a warning to organisms. To adapt to extreme environments, some terrestrial animals have evolved to be inherently insensitive to pain. Cetaceans are known as supposedly indifferent to pain from soft tissue injury representatives of marine mammals. However, the molecular mechanisms that explain how cetaceans are adapted to pain in response to seawater environment remain unclear. Here, we performed a molecular evolutionary analysis of pain-related genes in selected representatives of cetaceans. ASIC4 gene was identified to be pseudogenized in all odontocetes (toothed whales) except from Physeter macrocephalus (sperm whales), and relaxed selection of this gene was detected in toothed whales with pseudogenized ASIC4. In addition, positive selection was detected in pain perception (i.e., ASIC3, ANO1, CCK, and SCN9A) and analgesia (i.e., ASIC3, ANO1, CCK, and SCN9A) genes among the examined cetaceans. In this study, potential convergent amino acid substitutions within predicted proteins were found among the examined cetaceans and other terrestrial mammals, inhabiting extreme environments (e.g., V441I of TRPV1 in cetaceans and naked mole rats). Moreover, specific amino acid substitutions within predicted sequences of several proteins were found in the studied representatives of cetaceans (e.g., F56L and D163A of ASIC3, E88G of GRK2, and F159L of OPRD1). Most of the substitutions were located within important functional domains of proteins, affecting their protein functions. The above evidence suggests that cetaceans might have undergone adaptive molecular evolution in pain-related genes through different evolutionary patterns to adapt to pain, resulting in greater sensitivity to pain and more effective analgesia. This study could have implications for diagnosis and treatment of human pain.

1. Introduction

Pain is our body’s alarm system, warning us about dangers in the environment, wounds, and the presence of diseases. It is necessarily unpleasant, which enables protection by driving immediate attention, action, and adaptive learning [1]. Thus, the ability to detect harmful stimuli is critical for the survival of organisms. Pain is not only the result of a rigid hardware system, but also the interactions between highly plastic molecules and circuits [2]. Nociceptive stimuli that cause pain are sensed by special afferent nerve fibers called primary sensory neurons or nociceptors [3]. At the end of primary afferent neurons, four main nociceptors, i.e., acid-sensitive ion channels (ASICs), G protein-coupled receptors (GPCRs), P2X receptors (P2RXs), and transient receptor potential channels (TRP) are present [4,5]. In response to mechanical, thermal, and chemical stimuli, nociceptors cause membrane depolarization, which is called as generator potential. The initiated generator potentials are then converted into action potentials (AP) by ion channels, especially Nav1.7 channels, in adjacent areas. Next, the action potentials are transmitted further along the afferent sensory fibers wrapped by the epineurium. These fibers (whose body is located in the dorsal root ganglion) terminate in the spinal cord, where they communicate with the secondary neurons of the spinothalamic tract through synapses. Finally, the nociceptive signals are transmitted to the high-level brain center via the secondary neurons [6].
The pain response can be divided into physiological, psychological, and behavioral aspects, of which the latter one is the most easily to recognize, which is manifested as avoidance, crying, shouting, and so on [7,8,9]. However, during the process of adaptation to extreme environments some animals have lost their sense of pain or become insensitive to it in certain situations. For example, the African naked mole-rat (Heterocephalus glaber) acquired insensitivity to pain caused by acid stimuli. It was found that cutaneous nerves of the species are lacking in substance P (SP) and calcitonin gene-related peptide (CGRP) that are associated with nociceptive primary afferent fibers [10,11]. Another study indicated a mutation in TrkA gene as putative mechanism of pain insensitivity in the species [12]. Other example is pallid bat (Antrozous pallidus) that feeds on various species of scorpions. The species can endure multiple scorpion’s stings during hunting without behavioral pain response (e.g., avoidance and vocalizing) [13]. It was suggested that resistance to pain caused by scorpion venom might be acquired via amino acid substitutions in some proteins of Nav1.7 channel (e.g., E284G) [14].
So far, studies have investigated the mechanisms behind pain adaptation in mammals that are adapted to extreme terrestrial environments, but there are few studies on the mechanisms behind pain adaptation in marine mammals. Salinity is an important feature that distinguishes the marine from the terrestrial environment. The present in the saltwater sodium chloride, magnesium chloride and magnesium sulfate can directly stimulate the pain sensation in sea mammals when injury occurs [15]. In addition, salt has a strong water absorption function, which causes the edema around the wound to be more serious, further increasing the compression of pain nerves and causing acute pain [16]. Seals and sea lions that were heat-marked showed wound licking behavior and they stay in the water for a time about 12% shorter than before wounds acquisition [17,18]. It indicates that marine mammals can sense pain, and increased pain sensation of wounds caused by its contact with salt water. Moreover, bottlenose dolphins (Tursiops 2runcates) respond to painful stimuli with protest/withdrawal movements following a needle stick, suggesting that dolphins can perceive pain. However, bottlenose dolphins show apparent indifference to pain after soft tissue injury [19]. This also suggests that cetaceans may have adapted to pain in response to the marine environment. However, the molecular mechanism of this adaptation in cetaceans is still unclear.
At present, the most commonly used analgesia (the treatment of pain) methods are acupuncture and opioid therapy, which are mainly based on raising the pain threshold and attenuating/blocking pain signaling. In this study, we defined the genes reported to be involved in pain perception/sensation as pain genes or pain perception/sensation gene, and to be involved in analgesia as analgesia genes. The above pain perception/sensation and analgesia genes are collectively referred to as pain-related genes. We performed molecular evolutionary analyses of genes related to pain sensation and analgesia in representative mammals to explore the molecular mechanism of pain adaptation in cetaceans. This study could have implications for diagnosis and treatment of human pain.

2. Materials and Methods

2.1. Selected Species, Candidate Genes and Sequence Acquisition

Our research included the sequences of 31 pain-related genes, encompassing 18 pain perception genes: SCN9A, SCN10A, PRDM12, ANO1, CCK, MARK1, MARK2, RUNX1, NGF, NTRK1, CACNA1G, TRPV1, TRPV4, P2RX3, P2RX4, ASIC1, ASIC3 and ASIC4 and 13 analgesia genes: CA8, GRK2, ARRB2, KCNA1, KCNA2, TRAAK, KCNIP3, KCNJ10, ANO3, OPRM1, OPRD1, OPRK1, and OPRL1. All the selected genes have been proven to be associated with the pain (Supplementary Table S1 [20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47]). In total, 43 representatives of mammalian species were examined, including: 17 Cetacea, 4 Artiodactyla, 8 Carnivora, 2 Perissodactyla, 3 Chiroptera, 3 Primates, 1 Scandentia, 3 Rodentia, 1 Sirenia, and 1 Proboscidea. Most of the gene sequences were obtained from NCBI (http://www.ncbi.nlm.nih.gov/ accessed on 15 March 2022) and OrthoMaM [48] (https://orthomam.mbb.cnrs.fr/ accessed on 15 March 2022) (Supplementary Table S2). The gene sequences that are not available from the online databases, were extracted directly from well-annotated genomes stored in NCBI by our custom script. In addition, a set of genes of comparable size, 30 housekeeping genes (ACTB, ACTG1, ARF4, ATP6V1F, BECN1, CALM2, EIF4A2, ERH, HBOA, HINT1, HPCAL1, ID3, LASP1, LDHA, NFKBIA, PABPC1, PGK1, RPL11, RPL32, RPLP1, RPS18, RPS19, RPS25, RPS27A, SEC61B, SEPTIN7, UBB, YWHAB, YWHAZ, and ZFP36L1) of the above 43 species were used as a control set of genes to prove that the analysis results of pain-related genes are not randomly generated and are unusual. Housekeeping genes are being expressed stably in all conditions and cells, essentially, belong to cellular maintenance pathways and be conserved. Control genes sequences of human were downloaded from NCBI, and gene sequences of other species were blasted from genomes by using gene sequences of human as queries. The information of genomes was shown in Supplementary Table S3. The nucleotide and amino acid sequences (translated by MEGA) of each gene were aligned by MUSCLE in MEGA [49] and trimmed manually.

2.2. Screening for and Validating Pseudogenes

If premature stop codons, frame-shifting indels, and/or splicing site mutations were found in genes by alignment, then they were identified as potential pseudogenes and validated by PCR. Genomic DNA was extracted from cetacean muscle collected from baiji dolphin (Lipotes vexillifer) and Yangtze Finless Porpoise (Neophocaena asiaeorientalis) using the standard phenol/chloroform extraction method followed by ethanol precipitation [50]. All cetacean samples used in this study were collected from dead individuals in the wild with accordance to the ethical guidelines and legal requirements in China. DNA concentration and quality were measured by a UV spectrophotometer (Thermo Scientific of the U.S., NANODROP 2000). Primers (Supplementary Table S4) were designed in the conserved regions at both ends of the mutational locus in cetaceans by the online Primer-BLAST tool of NCBI. PCR amplification was performed according to the protocol of TIANGEN of the China (KT121221 2 × Taq Plus). PCR products were detected by agarose gel electrophoresis, and the product with the most distinct band was sent to Sangon Biotech for sequencing.

2.3. Selective Pressure Analysis

Comparison of ω (dN/dS), the ratio of nonsynonymous (dN) to synonymous (dS) substitution rates, among branches is an indication of the form and intensity of natural gene selection, where ω > 1, ω = 1 and ω < 1 indicate positive, neutral and purifying selection, respectively. The codon-based maximum likelihood models implemented in the CODEML program in PAML 4.7a [51] were employed to estimate the rates of synonymous (dS) and nonsynonymous substitutions (dN), as well as the dN/dS ratio (ω).
For pseudogenes, five models were used to determine the selection pressure in different branches. Model A assumes that all branches have the same ω values, model B is similar to model A, but ω is fixed to 1 (ω = 1). Model A versus model B was used to detect an overall selection pressure for pseudogenes. Model C assumes that particular pseudogenized branches have a common ω2, model D is similar to model C, but ω2 is fixed to 1 (ω2 = 1). Model C versus model D was used to detect the selection pressure for each pseudegenized branches. Finally, model E, which allowed different branches to have their own ω independently, was compared with model C to assess the divergence in the selective pressure across the phylogenetic branches [52].
In order to explore the evolutionary pattern of pain-related genes in cetaceans, we used the Branch Model and Branch-site Model in the PAML software package to perform selection pressure analysis on the selected genes. All the phylogenetic trees used in this study were obtained from the Timetree [53] website (http://www.timetree.org/ accessed on 20 March 2022). Foreground branches (tested as the pedigree being selected) and background branches (other pedigrees) are required to be defined a priori in the branch model and branch site model. The two-ratio model was treated as an alternative hypothesis for the branch model, which allows for the foreground and background branches have different ω values. The one-ratio model was used as the corresponding null hypothesis, which assumed that all evolutionary branches have the same ω value. In the branch-site model, the alternative hypothesis allowed each codon of the foreground branch to have its own ω value and allowed its ω to be greater than 1 (positive selection model: 0 < ω0 < 1, ω1 = 1 and ω2 ≥ 1); the null hypothesis (neutral model: 0 < ω0 < 1, ω1 = 1 and ω2 = 1) did not allow positive selection to appear. Then, we used the likelihood ratio test (LRT) with a χ2 distribution to determine which models were statistically different from the null model (threshold of p < 0.05). Moreover, the p-values of all genes were multiple-calibrated by false discovery rate (FDR) using the method of Benjamini–Hochberg [54], namely FDR (padjust) < 0.05. Finally, we used Bayes empirical Bayes (BEB) analysis to determine sites under positive selection with posterior probabilities ≥ 0.8 [55].

2.4. Labeling Positive Selection Sites on the Three-Dimensional Structure of Proteins

To emphasize the significance of positively selected sites in function of protein, we mapped these sites onto the predicted 3D structure of proteins. First, we used protein sequences of the bottlenose dolphin to predict the 3D structure by I-TASSER (https://zhanglab.ccmb.med.umich.edu/services/ accessed on 5 April 2022). Then, we used the Pfam (http://pfam.xfam.org/ accessed on 8 April 2022) website to search for domains of each protein. Finally, the EzMol (http://www.sbg.bio.ic.ac.uk/ezmol/ accessed on 17 April 2022) was used to annotate the positively selective sites and functional domains in the obtained 3D structure.

2.5. Identification of Convergent Amino Acids and Cetaceans-Specific Amino Acid Substitutions

To explore whether there are convergent molecular mechanisms among species that adapt to pain and whether cetaceans have their own unique molecular evolutionary mechanisms, we used MEGA 6 software and Alignment of Fasparser 2.10.0 software [56] to compare the amino acid sequences of pain-related genes, where set target group = 1, another group > 0.85, and ignore gaps. We defined that, for an amino acid site, if multiple species adapted to pain or cetaceans shared an amino acid that was unique to most other species, then that amino acid was convergent and related to pain insensitivity or cetaceans-specific amino acid substitutions.

2.6. Predicting the Effect of Cetacean-Specific Amino Acid Substitutions on Protein Function

We used Polyphen2 [57] (Phenotyping Version2, http://genetics.bwh.harvard.edu/pph2/ accessed on 5 May 2022) and PROVEAN [58] (Protein variation effect analyzer, http://provean.jcvi.org/seq_submit.php accessed on 5 May 2022) to predict the effect of cetacean-specific amino acid substitutions on protein function. The human protein sequences and cetacean-specific amino acid substitutions were inputted into the corresponding to protein sequences boxes and substitutions/mutation boxes of Polyphen2 and PROVEAN. The Polyphen2 score over the threshold is considered to affect protein function; a PROVEAN score less than or equal to the threshold −2.5 is considered deleterious (Deleterious), and a score higher than −2.5 is considered neutral, meaning the mutation may not affect the protein’s function.

3. Results

3.1. Pseudogenization of ASCI4 in Toothed Whales

Of the 31 pain-related genes examined in this study, ASCI4 was identified to be pseudogenized in all toothed whales except the sperm whale due to premature stop codons detected in the first exon of the sequences (Figure 1, Supplementary Figure S1). The other 30 genes were intact in all the cetaceans studied. In a control set of 30 housekeeping genes, we did not detect pseudogenization.

3.2. Relaxed Selection in ASIC4 Detected in Toothed Whales

Then, we calculated a series of ω values to determine whether the ASIC4 gene has experienced relaxed selection pressure (Table 1). The average ω across the tree was estimated to be 0.087 in model A, which was significantly lower than in model B (ω = 1.000), indicating strong purifying selection on ASIC4 (p = 0; Table 1). Model C allowed toothed whales with pseudogenized ASIC4 having ω2 (ω2 = 0.200), fitting the data significantly better than model A (ω = 0.087; p = 0.004 in model C vs. A; Table 1). This indicates the relaxed selection pressure on ASIC4 gene in toothed whales. Model C (ω2 = 0.200) of the pseudogene branch was significantly different from model D (ω2 =1.000; p = 6.31221 × 10−8 in model C vs. D; Table 1), suggesting that the selective pressure on ASIC4 was not completely relaxed. Finally, model E significantly fitted the data compared with model C (p = 6.15578 × 10−12 in model C vs. E; Table 1), indicating a divergence in the selective pressure across the tree.

3.3. Cetacean Lineages Displayed Stronger Positive Selection than Unreported Pain-Adaption Branches and Reported Pain-Adaption Non-Cetacean Branches

Three positive selection instances in three genes (two pain perception genes, one analgesia gene) were detected in the mammalian dataset using the branch model and all were found in cetaceans, especially in baleen whales (Supplementary Table S5). In turns, ten positive selection instances in six pain perception genes and four positive selection instances in four analgesia genes were identified in the mammalian dataset using the branch-site model, when compared with a control set of housekeeping genes (no positive selection was detected). Among these 14 positive selection instances, cetacean lineages (7/14) displayed stronger positive selection than unreported pain-adaption branches (4/14) and reported pain-adaption non-cetacean branches (3/14) (Figure 2, Supplementary Table S6). In addition, positive selection with pain sensation genes were detected in selectively pain-insensitive species such as bats, naked mole rats, cows, and mice (for which no pain-related reports were reported), but they were mainly concentrated in cetacean branches (4/10), especially baleen whales and sperm whales (the only toothed whale with no pseudogene detected on ASIC4). Among the analgesia genes, we only detected positive selection signals in the cetacean lineages of three (ARRB2, KCNK4, and OPRL1), and positive selection signals were detected in OPRM1 in pigs. In other words, analgesic genes (3/4) showed more signals in cetacean lineages than pain perception genes (4/10) (Figure 2, Supplementary Table S6). Furthermore, most positive selection sites were located at or near the protein’s functional domain. For example, the 118th position of the ARRB2 protein is located in the Arrestin_N domain; the 103rd and 104th positions of the KCNK4 protein are located in the Ion_trans_2 domain; the 225th position of the OPRL1 protein is located in the 7tm_1 domain (Figure 3).

3.4. Potential Molecular Convergence among Species Adapted to Pain

We divided cetaceans, other marine mammals, bats, and naked mole rats into four groups of mammals that may be adapted to pain. Seven potential convergent amino acids in four genes among them were identified, when compared with a control set of housekeeping genes (no such molecular convergence). ASIC3 has two convergent amino acids, R226K of bats and naked mole rats, and I384M of cetaceans and bats; L70P of bats and naked mole rats, and V441I of cetaceans and naked mole rats were found in TRPV1; SCN9A has two convergent amino acids, V169I and A825S of bats and naked mole rats; C295S of cetaceans and bats was found in ANO3. Among them, ASIC3, TRPV1, and SCN9A are pain sensation genes and ANO3 is an analgesia gene (Figure 4).

3.5. Cetaceans-Specific Amino Acid Substitutions

We also found cetacean-specific amino acid substitutions, when compared with a control set of housekeeping genes (no such cetaceans-specific amino acid substitutions). Nine cetacean-specific amino acid substitutions were identified, of which ASIC3 F56L, ASIC3 D163A, TRPV1 D471G, and NGF T202A were detected in pain perception genes, and GRK2 E88G, OPRD1 F159L, KCNIP3 A127S, KCNJ10 G350R, and ANO3 V635I were detected in analgesic genes (Figure 5A). Furthermore, ASIC3 D163A, GRK2 E88G, and KCNJ10 G350R were predicted by Polyphen2 (scores > threshold of 0.85) to potentially affect protein function (Figure 5B, Supplementary Table S7); TRPV1 D471G, NGF T202A, and OPRD1 F159L were predicted by PROVEAN (scores < threshold of −2.5) to potentially affect protein function (Figure 5C, Supplementary Table S7).

4. Discussion

Pain is currently one of the most common disease symptoms in mammals. In 2018, the International Association for the Study of Pain (IASP) defined pain as “an unpleasant sensory and emotional experience associated or similarly associated with tissue damage or potential tissue damage.” From a neurobiological point of view, pain can be divided into three categories: nociceptive pain caused by noxious stimuli, inflammatory pain associated with tissue damage and immune cell infiltration, and pathological pain due to damage to the nervous system or its dysfunction [59]. Lack of nociceptive pain is a problem. There are reports on its molecular genetic cause—e.g., congenital insensitivity or indifference to pain due to in-frame deletions or rare missense of SCN9A, which encodes the Nav1.7, a voltage-gated sodium channel [60]. Another example is that mutations in the loss of function of nerve growth factor and the neurotrophic tyrosine kinase receptor TrkA typically lead to self-mutilation, bone fractures, multiple scars, joint deformities, amputations, and early death [61].
In this study, we discovered a pseudogenization of the acid-sensitive ion channel 4 gene (ASIC4) in all toothed whales except sperm whales. Furthermore, we showed that ASIC4 experiences relaxed selection pressure on the branch on which we detected pseudogenization via selection pressure analysis. Notably, ASIC4 is the newest member of the ASIC gene family of neuronal proton-gated cation channels. ASIC4 itself is inactive and cannot form proton-gated ion channel complexes [62]. To some extent, it can promote the degradation of ASIC1 or inhibits its expression [63]. ASIC1 is the most abundant proton-gated ion channel gene in mammals. It is involved in the transmission of pain signals, and its degradation or reduced expression can lead to insensitivity to pain [64]. Therefore, the relaxed selection of ASIC4 in toothed whales may have deprived it of its inhibitory effect on ASIC1, thereby making toothed whales more sensitive to pain. While examining the sperm whale and baleen whale lineages, we detected positive selection sites in important functional domains of the pain protein (i.e., ASIC3, ANO1, CCK, and SCN9A), also suggesting that sperm and baleen whales may also be more sensitive to pain. That is to say, different cetaceans may have evolved increased pain sensitivity through different paths, which is consistent with reports that bottlenose dolphins respond to painful stimuli through protest/withdrawal movements after acupuncture [19]. We also detected positive selection in pain sensation genes of species such as bats, naked mole rats, cows, and mice, which may be related to selective pain insensitivity. For example, primary afferent nociceptors in naked mole-rats are insensitive to acid stimuli, but nociceptors of them respond vigorously to capsaicin, and capsaicin-sensitive acid transduction by sensory neurons (TRPV1) is expressed in sensory neurons [11].
In addition, we specifically detected positive selection in three analgesic genes (ARRB2, KCNK4 and OPRL1) in cetaceans. ARRB2 is a member of the arrestin family. Its N-terminal domain can bind to phosphorylated GPCRs to form trimers and block pain signaling, exerting analgesic effects [65]. Thus, the 118th positive selection site (located in the N-terminal domain) of ARRB2 in cetaceans may strengthen the negative regulation of pain signaling by ARRB2, thereby enhancing analgesia after suffering pain. It seems to explain that bottlenose dolphins show apparent indifference to pain after soft tissue injury [19].
Convergent evolution has always been a hot spot in the field of evolutionary biology. We identified seven potential convergent amino acids in four genes of species adapted to pain; most of these convergent amino acids were located in important functional domains of proteins, suggesting that different species that are adapted to pain may have undergone convergent evolution at the molecular level. Among the four genes, SCN9A-mapped mutations in pain disorders have been studied fully—e.g., p.Arg 222 His, evaluated from patients from a large family with early-onset pain symptoms, a familial mutation (I136V) and two sporadic mutations (I848T, V1316A) identified in patients with primary erythromelalgia (PE) [66,67,68]. The convergent amino acid substitutions (V169I and A825S) identified in this study, similar to these mutations, are located in the DI/S1 and DII important domains, respectively. Thus, these two convergent sites likely affect the function of SCN9A and contribute to pain insensitivity in bats and naked mole rats. Finally, we also found nine cetacean-specific amino acid substitutions and six of them were identified to affect the function of proteins. This means that cetaceans may have experienced a unique evolutionary trajectory different from other species adapted to pain.

5. Conclusions

Overall, the above evidence reveals that in response to seawater environmental conditions, cetaceans might have undergone adaptive molecular evolution in pain and analgesic traits through different evolutionary patterns. On the one hand, cetaceans may be more sensitive to pain in certain situations, which is conducive to perceiving danger and avoiding risks, and on the other hand, they may be more effective in analgesia, which is beneficial for them to relieve pain after suffering “pain”.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ani12243571/s1, Table S1: 31 pain-related genes; Table S2. Information of genes sequences in this study; Table S3. Information of genome with Blast; Table S4. Primer sequence information; Table S5. Branch detected positive selection for pain-related genes by the branch model; Table S6. Branch detected positive selection for pain-related genes by the branch site model; Table S7. Predicted deleterious score by PROVEAN and PolyPhen2; Figure S1. Identification of ASIC4 gene of L.vexillifer and N. asiaeorientalis. The purplish red box indicates the stop codon.

Author Contributions

Conceptualization, X.D. and W.R.; formal analysis, X.D. and F.Y.; methodology, X.D., F.Y., and X.H.; writing—original draft preparation, X.D.; writing—review and editing, W.R. and G.Y.; project administration, W.R.; resources, W.R.; funding acquisition, W.R., G.Y., and S.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the National Natural Science Foundation of China (grant nos. 31872219, 31370401, 32030011, 31630071, 31772448), National Key Programme of Research and Development, Ministry of Science and Technology of China (grant no. 2022YFF1301600), and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Institutional Review Board Statement

All animal study protocol was approved by the Animal Ethical and Welfare Committee (Approval No: IACUC-20200512) of Nanjing Normal University (NNU).

Informed Consent Statement

Not applicable.

Data Availability Statement

Most of the data supporting the findings of this study can be found within the supplementary files, some other data are available from the corresponding author upon reasonable request.

Acknowledgments

We thank Simin Chai, Xu Zhou, Tianzhen Wu, and Xing Liu for their suggestion and technical support during this study.

Conflicts of Interest

All authors declare no conflicts of interest.

References

  1. Tracey, I.; Woolf, C.J.; Andrews, N.A. Composite Pain Biomarker Signatures for Objective Assessment and Effective Treatment. Neuron 2019, 101, 783–800. [Google Scholar] [CrossRef] [Green Version]
  2. Apkarian, A.V.; Bushnell, M.C.; Treede, R.-D.; Zubieta, J.-K. Human brain mechanisms of pain perception and regulation in health and disease. Eur. J. Pain 2005, 9, 463. [Google Scholar] [CrossRef]
  3. Dubin, A.E.; Patapoutian, A. Nociceptors: The sensors of the pain pathway. J. Clin. Investig. 2010, 120, 3760–3772. [Google Scholar] [CrossRef] [Green Version]
  4. Ji, R.-R.; Xu, Z.-Z.; Gao, Y.-J. Emerging targets in neuroinflammation-driven chronic pain. Nat. Rev. Drug Discov. 2014, 13, 533–548. [Google Scholar] [CrossRef] [Green Version]
  5. Moore, C.; Gupta, R.; Jordt, S.-E.; Chen, Y.; Liedtke, W.B. Regulation of Pain and Itch by TRP Channels. Neurosci. Bull. 2018, 34, 120–142. [Google Scholar] [CrossRef] [Green Version]
  6. King, G.F.; Vetter, I. No Gain, No Pain: NaV1.7 as an Analgesic Target. ACS Chem. Neurosci. 2014, 5, 749–751. [Google Scholar] [CrossRef]
  7. Philips, H.C.; Jahanshahi, M. The components of pain behaviour report. Behav. Res. Ther. 1986, 24, 117–125. [Google Scholar] [CrossRef] [PubMed]
  8. Philips, H.C. Avoidance behaviour and its role in sustaining chronic pain. Behav. Res. Ther. 1987, 25, 273–279. [Google Scholar] [CrossRef]
  9. Gigliuto, C.; De Gregori, M.; Malafoglia, V.; Raffaeli, W.; Compagnone, C.; Visai, L.; Petrini, P.; Avanzini, M.A.; Muscoli, C.; Viganò, J.; et al. Pain assessment in animal models: Do we need further studies? J. Pain Res. 2014, 7, 227–236. [Google Scholar] [CrossRef] [Green Version]
  10. Cao, Y.Q.; Mantyh, P.W.; Carlson, E.J.; Gillespie, A.-M.; Epstein, C.J.; Basbaum, A.I. Primary afferent tachykinins are required to experience moderate to intense pain. Nature 1998, 392, 390–394. [Google Scholar] [CrossRef]
  11. Park, T.J.; Lu, Y.; Jüttner, R.; Smith, E.S.J.; Hu, J.; Brand, A.; Wetzel, C.; Milenkovic, N.; Erdmann, B.; Heppenstall, P.A.; et al. Selective inflammatory pain insensitivity in the African naked mole-rat (Heterocephalus glaber). PLoS Biol. 2008, 6, e13. [Google Scholar] [CrossRef] [PubMed]
  12. Omerbašić, D.; Smith, E.S.J.; Moroni, M.; Homfeld, J.; Eigenbrod, O.; Bennett, N.C.; Reznick, J.; Faulkes, C.G.; Selbach, M.; Lewin, G.R. Hypofunctional TrkA Accounts for the Absence of Pain Sensitization in the African Naked Mole-Rat. Cell Rep. 2016, 17, 748–758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Gurevitz, M. Mapping of scorpion toxin receptor sites at voltage-gated sodium channels. Toxicon 2012, 60, 502–511. [Google Scholar] [CrossRef]
  14. Hopp, B.H.; Arvidson, R.S.; Adams, M.E.; Razak, K.A. Arizona bark scorpion venom resistance in the pallid bat, Antrozous pallidus. PLoS ONE 2017, 12, e0183215. [Google Scholar] [CrossRef] [Green Version]
  15. Ahern, G.P.; Brooks, I.M.; Miyares, R.L.; Wang, X.-b. Extracellular Cations Sensitize and Gate Capsaicin Receptor TRPV1 Modulating Pain Signaling. J. Neurosci. 2005, 25, 5109. [Google Scholar] [CrossRef] [Green Version]
  16. Longhurst, H.; Cicardi, M. Hereditary angio-oedema. Lancet 2012, 379, 474–481. [Google Scholar] [CrossRef] [Green Version]
  17. Walker, K.A. Behavioural and Physiological Responses of Steller SEA Lions to Invasive Marking Techniques: Evidence of Post-operative Pain. Ph.D. Thesis, University of British Columbia, Vancouver, BC, Canada, 2010. [Google Scholar]
  18. Walker, K.A.; Mellish, J.E.; Weary, D.M. Effects of hot-iron branding on heart rate, breathing rate and behaviour of anaesthetised Steller sea lions. Vet. Rec. 2011, 169, 363. [Google Scholar] [CrossRef]
  19. Zasloff, M. Observations on the Remarkable (and Mysterious) Wound-Healing Process of the Bottlenose Dolphin. J. Investig. Dermatol. 2011, 131, 2503–2505. [Google Scholar] [CrossRef] [Green Version]
  20. Xue, Y.; Chidiac, C.; Herault, Y.; Gaveriaux-Ruff, C. Pain behavior in SCN9A (Nav1.7) and SCN10A (Nav1.8) mutant rodent models. Neurosci. Lett. 2021, 753, 135844. [Google Scholar] [CrossRef]
  21. Landy, M.A.; Goyal, M.; Casey, K.M.; Liu, C.; Lai, H.C. Loss of Prdm12 during development, but not in mature nociceptors, causes defects in pain sensation. Cell Rep. 2021, 34, 108913. [Google Scholar] [CrossRef]
  22. Liu, Y.; Liu, Z.; Wang, K. The Ca2+-activated chloride channel ANO1/TMEM16A: An emerging therapeutic target for epithelium-originated diseases? Acta Pharm. Sin. B 2021, 11, 1412–1433. [Google Scholar] [CrossRef] [PubMed]
  23. McRoberts, J.W. Cholecystokinin and pain: A review. Anesthesia Prog. 1986, 33, 87–90. [Google Scholar]
  24. Yu, L.-N.; Sun, L.-H.; Wang, M.; Yan, M. Research progress of the role and mechanism of extracellular signal-regulated protein kinase 5 (ERK5) pathway in pathological pain. J. Zhejiang Univ. Sci. B 2016, 17, 733–741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Obata, K.; Noguchi, K. MAPK activation in nociceptive neurons and pain hypersensitivity. Life Sci. 2004, 74, 2643–2653. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, J.W.; Stifani, S. Roles of Runx Genes in Nervous System Development. In RUNX Proteins in Development and Cancer; Groner, Y., Ito, Y., Liu, P., Neil, J.C., Speck, N.A., van Wijnen, A., Eds.; Springer: Singapore, 2017; pp. 103–116. [Google Scholar]
  27. Testa, G.; Cattaneo, A.; Capsoni, S. Understanding pain perception through genetic painlessness diseases: The role of NGF and proNGF. Pharmacol. Res. 2021, 169, 105662. [Google Scholar] [CrossRef]
  28. Indo, Y. Genetics of congenital insensitivity to pain with anhidrosis (CIPA) or hereditary sensory and autonomic neuropathy type IV. Clin. Auton. Res. 2002, 12, I20–I32. [Google Scholar] [CrossRef]
  29. Coutelier, M.; Blesneac, I.; Monteil, A.; Monin, M.-L.; Ando, K.; Mundwiller, E.; Brusco, A.; Le Ber, I.; Anheim, M.; Castrioto, A.; et al. A Recurrent Mutation in CACNA1G Alters Cav3.1 T-Type Calcium-Channel Conduction and Causes Autosomal-Dominant Cerebellar Ataxia. Am. J. Hum. Genet. 2015, 97, 726–737. [Google Scholar] [CrossRef] [Green Version]
  30. Davis, J.B.; Gray, J.; Gunthorpe, M.J.; Hatcher, J.P.; Davey, P.T.; Overend, P.; Harries, M.H.; Latcham, J.; Clapham, C.; Atkinson, K.; et al. Vanilloid receptor-1 is essential for inflammatory thermal hyperalgesia. Nature 2000, 405, 183–187. [Google Scholar] [CrossRef]
  31. Wirkner, K.; Sperlagh, B.; Illes, P. P2X3 Receptor Involvement in Pain States. Mol. Neurobiol. 2007, 36, 165–183. [Google Scholar] [CrossRef]
  32. Zhang, W.-J.; Luo, H.-L.; Zhu, Z.-M. The role of P2X4 receptors in chronic pain: A potential pharmacological target. Biomed. Pharmacother. 2020, 129, 110447. [Google Scholar] [CrossRef]
  33. Chauhan, S.A.; Sahoo, C.G.; Dikhit, R.M.; Das, P. Acid-Sensing Ion Channels Structural Aspects, Pathophysiological Importance and Experimental Mutational Data Available Across Various Species to Target Human ASIC1. Curr. Drug Targets 2019, 20, 111–121. [Google Scholar] [CrossRef] [PubMed]
  34. Wemmie, J.A.; Taugher, R.J.; Kreple, C.J. Acid-sensing ion channels in pain and disease. Nat. Rev. Neurosci. 2013, 14, 461–471. [Google Scholar] [CrossRef] [PubMed]
  35. Fu, E.S.; Erasso, D.M.; Zhuang, G.Z.; Upadhyay, U.; Ozdemir, M.; Wiltshire, T.; Sarantopoulos, K.D.; Smith, S.B.; Maixner, W.; Martin, E.R.; et al. Impact of human CA8 on thermal antinociception in relation to morphine equivalence in mice. Neuroreport 2017, 28, 1215–1220. [Google Scholar] [CrossRef] [PubMed]
  36. Wang, H.; Heijnen, C.J.; van Velthoven, C.T.J.; Willemen, H.L.D.M.; Ishikawa, Y.; Zhang, X.; Sood, A.K.; Vroon, A.; Eijkelkamp, N.; Kavelaars, A. Balancing GRK2 and EPAC1 levels prevents and relieves chronic pain. J. Clin. Investig. 2013, 123, 5023–5034. [Google Scholar] [CrossRef] [Green Version]
  37. Ozberk, D.; Haywood, A.; Sutherland, H.G.; Yu, C.; Albury, C.L.; Zunk, M.; George, R.; Good, P.; Griffiths, L.R.; Hardy, J.; et al. Association of polymorphisms in ARRB2 and clinical response to methadone for pain in advanced cancer. Pharmacogenomics 2022, 23, 281–289. [Google Scholar] [CrossRef]
  38. Imbrici, P.; Altamura, C.; Gualandi, F.; Mangiatordi, G.F.; Neri, M.; De Maria, G.; Ferlini, A.; Padovani, A.; D’Adamo, M.C.; Nicolotti, O.; et al. A novel KCNA1 mutation in a patient with paroxysmal ataxia, myokymia, painful contractures and metabolic dysfunctions. Mol. Cell. Neurosci. 2017, 83, 6–12. [Google Scholar] [CrossRef]
  39. Peng, Y.; Zhang, Q.; Cheng, H.; Yan, G.; Xing, C. Upregulation of ubiquitin conjugating enzyme E2B (Ube2b) ameliorates neuropathic pain by regulating Kcna2 (potassium voltage-gated channel subfamily A member 2) in primary afferent neurons. Bioengineered 2021, 12, 7470–7480. [Google Scholar] [CrossRef]
  40. Noël, J.; Zimmermann, K.; Busserolles, J.; Deval, E.; Alloui, A.; Diochot, S.; Guy, N.; Borsotto, M.; Reeh, P.; Eschalier, A.; et al. The mechano-activated K+ channels TRAAK and TREK-1 control both warm and cold perception. EMBO J. 2009, 28, 1308–1318. [Google Scholar] [CrossRef] [Green Version]
  41. Guo, Y.-P.; Zhi, Y.-R.; Liu, T.-T.; Wang, Y.; Zhang, Y. Global Gene Knockout of Kcnip3 Enhances Pain Sensitivity and Exacerbates Negative Emotions in Rats. Front. Mol. Neurosci. 2019, 12, 5. [Google Scholar] [CrossRef] [Green Version]
  42. Mannerak, M.A.; Lashkarivand, A.; Eide, P.K. Trigeminal neuralgia and genetics: A systematic review. Mol. Pain 2021, 17, 17448069211016139. [Google Scholar] [CrossRef]
  43. Ślęczkowska, M.; Almomani, R.; Marchi, M.; de Greef, B.T.A.; Sopacua, M.; Hoeijmakers, J.G.J.; Lindsey, P.; Salvi, E.; Bönhof, G.J.; Ziegler, D.; et al. Peripheral Ion Channel Gene Screening in Painful- and Painless-Diabetic Neuropathy. Int. J. Mol. Sci. 2022, 23, 7190. [Google Scholar] [CrossRef] [PubMed]
  44. Campa, D.; Gioia, A.; Tomei, A.; Poli, P.; Barale, R. Association of ABCB1/MDR1 and OPRM1 Gene Polymorphisms With Morphine Pain Relief. Clin. Pharmacol. Ther. 2008, 83, 559–566. [Google Scholar] [CrossRef] [PubMed]
  45. Seo, E.-J.; Efferth, T.; Panossian, A. Curcumin downregulates expression of opioid-related nociceptin receptor gene (OPRL1) in isolated neuroglia cells. Phytomedicine 2018, 50, 285–299. [Google Scholar] [CrossRef] [PubMed]
  46. Ho, K.W.D.; Wallace, M.R.; Staud, R.; Fillingim, R.B. OPRM1, OPRK1, and COMT genetic polymorphisms associated with opioid effects on experimental pain: A randomized, double-blind, placebo-controlled study. Pharm. J. 2020, 20, 471–481. [Google Scholar] [CrossRef]
  47. Nielsen, L.M.; Christrup, L.L.; Sato, H.; Drewes, A.M.; Olesen, A.E. Genetic Influences of OPRM1, OPRD1 and COMT on Morphine Analgesia in a Multi-Modal, Multi-Tissue Human Experimental Pain Model. Basic Clin. Pharmacol. Toxicol. 2017, 121, 6–12. [Google Scholar] [CrossRef] [Green Version]
  48. Scornavacca, C.; Belkhir, K.; Lopez, J.; Dernat, R.; Delsuc, F.; Douzery, E.J.P.; Ranwez, V. OrthoMaM v10: Scaling-Up Orthologous Coding Sequence and Exon Alignments with More than One Hundred Mammalian Genomes. Mol. Biol. Evol. 2019, 36, 861–862. [Google Scholar] [CrossRef] [Green Version]
  49. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [Green Version]
  50. Green, M.R.; Sambrook, J. Molecular Cloning: A Laboratory Manual, 4th ed.; Cold Spring Harbor Laboratory Pr: New York, NY, USA, 2012; Volume 3. [Google Scholar]
  51. Yang, Z. PAML 4: Phylogenetic Analysis by Maximum Likelihood. Mol. Biol. Evol. 2007, 24, 1586–1591. [Google Scholar] [CrossRef] [Green Version]
  52. Zhao, H.; Zhou, Y.; Pinto, C.M.; Charles-Dominique, P.; Galindo-González, J.; Zhang, S.; Zhang, J. Evolution of the sweet taste receptor gene Tas1r2 in bats. Mol. Biol. Evol. 2010, 27, 2642–2650. [Google Scholar] [CrossRef] [Green Version]
  53. Kumar, S.; Stecher, G.; Suleski, M.; Hedges, S.B. TimeTree: A Resource for Timelines, Timetrees, and Divergence Times. Mol. Biol. Evol. 2017, 34, 1812–1819. [Google Scholar] [CrossRef]
  54. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Society. Ser. B Methodol. 1995, 57, 289–300. [Google Scholar] [CrossRef]
  55. Yang, Z.; Wong, W.S.W.; Nielsen, R. Bayes Empirical Bayes Inference of Amino Acid Sites Under Positive Selection. Mol. Biol. Evol. 2005, 22, 1107–1118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Sun, Y.B. FasParser: A package for manipulating sequence data. Zool. Res. 2017, 38, 110–112. [Google Scholar] [PubMed] [Green Version]
  57. Adzhubei, I.A.; Schmidt, S.; Peshkin, L.; Ramensky, V.E.; Gerasimova, A.; Bork, P.; Kondrashov, A.S.; Sunyaev, S.R. A method and server for predicting damaging missense mutations. Nat. Methods 2010, 7, 248–249. [Google Scholar] [CrossRef] [Green Version]
  58. Choi, Y.; Chan, A.P. PROVEAN web server: A tool to predict the functional effect of amino acid substitutions and indels. Bioinformatics 2015, 31, 2745–2747. [Google Scholar] [CrossRef] [Green Version]
  59. Woolf, C.J. What is this thing called pain? J. Clin. Investig. 2010, 120, 3742–3744. [Google Scholar] [CrossRef]
  60. Raouf, R.; Quick, K.; Wood, J.N. Pain as a channelopathy. J. Clin. Investig. 2010, 120, 3745–3752. [Google Scholar] [CrossRef]
  61. Reilly, M.M.; Shy, M.E. Diagnosis and new treatments in genetic neuropathies. J. Neurol. Neurosurg. Psychiatry 2009, 80, 1304. [Google Scholar] [CrossRef]
  62. Gründer, S.; Geisler, H.-S.; Rainier, S.; Fink, J.K. Acid-sensing ion channel (ASIC) 4 gene: Physical mapping, genomic organisation, and evaluation as a candidate for paroxysmal dystonia. Eur. J. Hum. Genet. 2001, 9, 672–676. [Google Scholar] [CrossRef]
  63. Donier, E.; Rugiero, F.; Jacob, C.; Wood, J.N. Regulation of ASIC activity by ASIC4—New insights into ASIC channel function revealed by a yeast two-hybrid assay. Eur. J. Neurosci. 2008, 28, 74–86. [Google Scholar] [CrossRef]
  64. Lingueglia, E. Acid-Sensing Ion Channels (ASICs) in pain. Biol. Aujourdhui 2014, 208, 13. [Google Scholar] [CrossRef] [PubMed]
  65. Chen, G.; Xie, R.-G.; Gao, Y.-J.; Xu, Z.-Z.; Zhao, L.-X.; Bang, S.; Berta, T.; Park, C.-K.; Lay, M.; Chen, W.; et al. β-arrestin-2 regulates NMDA receptor function in spinal lamina II neurons and duration of persistent pain. Nat. Commun. 2016, 7, 12531. [Google Scholar] [CrossRef] [Green Version]
  66. Goldberg, Y.P.; MacFarlane, J.; MacDonald, M.L.; Thompson, J.; Dube, M.P.; Mattice, M.; Fraser, R.; Young, C.; Hossain, S.; Pape, T.; et al. Loss-of-function mutations in the Nav1.7 gene underlie congenital indifference to pain in multiple human populations. Clin. Genet. 2007, 71, 311–319. [Google Scholar] [CrossRef]
  67. Wu, M.-T.; Huang, P.-Y.; Yen, C.-T.; Chen, C.-C.; Lee, M.-J. A novel SCN9A mutation responsible for primary erythromelalgia and is resistant to the treatment of sodium channel blockers. PLoS ONE 2013, 8, e55212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Han, C.; Yang, Y.; te Morsche, R.H.; Drenth, J.P.H.; Politei, J.M.; Waxman, S.G.; Dib-Hajj, S.D. Familial gain-of-function Nav1.9 mutation in a painful channelopathy. J. Neurol. Neurosurg. Psychiatry 2017, 88, 233. [Google Scholar] [CrossRef] [PubMed]
Animals 12 03571 g001 550
Figure 1. Pseudogenizations were detected in the ASIC4 of almost all the toothed whales. A schematic diagram of the gene structure of ASIC4 in the bottlenose dolphin is shown (on the top). Orange rectangles indicate exons and the black vertical bar indicates the location of premature stop codons. All 10 toothed whales in this study, except for sperm whales, had premature stop codons in the first exon of ASIC4. Blue stars represent toothed whale species; blue triangles represent baleen whale species. Purplish red font represents species with detected pseudogenes, and purplish red boxes are stop codons.
Figure 1. Pseudogenizations were detected in the ASIC4 of almost all the toothed whales. A schematic diagram of the gene structure of ASIC4 in the bottlenose dolphin is shown (on the top). Orange rectangles indicate exons and the black vertical bar indicates the location of premature stop codons. All 10 toothed whales in this study, except for sperm whales, had premature stop codons in the first exon of ASIC4. Blue stars represent toothed whale species; blue triangles represent baleen whale species. Purplish red font represents species with detected pseudogenes, and purplish red boxes are stop codons.
Animals 12 03571 g001
Animals 12 03571 g002 550
Figure 2. Positive selection detected in the mammalian dataset using the branch-site model. Blue circles represent cetacean species. Reported pain-adapted species are shown in purplish red font; unreported pain-adapted species are shown in black font. The lineages for which positive selection was detected on pain perception genes are shown with vertical bars of different colors: ASIC3 (red), ANO1 (orange), CCK (purple), SCN9A (blue), SCN10A (green), TRPV1 (pink). The branches for which positive selection was detected on analgesia genes are indicated by diamonds of different colors: ARRB2 (red), KCNK4 (orange), OPRL1 (purple), OPRM1 (blue).
Figure 2. Positive selection detected in the mammalian dataset using the branch-site model. Blue circles represent cetacean species. Reported pain-adapted species are shown in purplish red font; unreported pain-adapted species are shown in black font. The lineages for which positive selection was detected on pain perception genes are shown with vertical bars of different colors: ASIC3 (red), ANO1 (orange), CCK (purple), SCN9A (blue), SCN10A (green), TRPV1 (pink). The branches for which positive selection was detected on analgesia genes are indicated by diamonds of different colors: ARRB2 (red), KCNK4 (orange), OPRL1 (purple), OPRM1 (blue).
Animals 12 03571 g002
Animals 12 03571 g003 550
Figure 3. The distribution of positive selection sites on the three-dimensional structure. (A–C) represent the distribution of positive selection sites in ARRB2, KCNK4, and OPRL1, respectively. These genes were specifically detected positive selection signals in cetacean lineages. The yellow represents the domain located by positive selection sites, the purplish red represents the positive selection sites, and the corresponding numbers represent the amino acid positions.
Figure 3. The distribution of positive selection sites on the three-dimensional structure. (A–C) represent the distribution of positive selection sites in ARRB2, KCNK4, and OPRL1, respectively. These genes were specifically detected positive selection signals in cetacean lineages. The yellow represents the domain located by positive selection sites, the purplish red represents the positive selection sites, and the corresponding numbers represent the amino acid positions.
Animals 12 03571 g003
Animals 12 03571 g004 550
Figure 4. Potential convergent amino acids in species adapted to pain. Four groups of mammals adapted to pain are shown with circles of different colors: cetaceans (blue), other marine mammals (orange), bats (green), naked mole rats (pink). Amino acids in pink and purple are convergent amino acids associated with pain sensation and analgesia genes, respectively. These convergent amino acids are annotated on the protein domains (on the top). Amino acid positions are at the human scale.
Figure 4. Potential convergent amino acids in species adapted to pain. Four groups of mammals adapted to pain are shown with circles of different colors: cetaceans (blue), other marine mammals (orange), bats (green), naked mole rats (pink). Amino acids in pink and purple are convergent amino acids associated with pain sensation and analgesia genes, respectively. These convergent amino acids are annotated on the protein domains (on the top). Amino acid positions are at the human scale.
Animals 12 03571 g004
Animals 12 03571 g005 550
Figure 5. Cetaceans-specific amino acid substitutions. (A) Cetacean-specific amino acid substitutions identified in pain-related genes. Blue circles represent cetacean species. Amino acids in pink and purple are cetacean-specific amino acid substitutions identified in pain perception and analgesic genes, respectively. Amino acid positions are at the human scale. (B) Bar graph of the score of Polyphen2 predicting the effect of cetacean-specific amino acids on protein function. The purplish red bar graph shows that the scores are more than the threshold of 0.85, meaning that the protein’s function may be affected. (C) Bar graph of the PROVEAN score predicting the effect of cetacean-specific amino acids on protein function. The purplish red bar graph shows that the scores are less than the threshold of −2.5, meaning that the protein’s function may be affected.
Figure 5. Cetaceans-specific amino acid substitutions. (A) Cetacean-specific amino acid substitutions identified in pain-related genes. Blue circles represent cetacean species. Amino acids in pink and purple are cetacean-specific amino acid substitutions identified in pain perception and analgesic genes, respectively. Amino acid positions are at the human scale. (B) Bar graph of the score of Polyphen2 predicting the effect of cetacean-specific amino acids on protein function. The purplish red bar graph shows that the scores are more than the threshold of 0.85, meaning that the protein’s function may be affected. (C) Bar graph of the PROVEAN score predicting the effect of cetacean-specific amino acids on protein function. The purplish red bar graph shows that the scores are less than the threshold of −2.5, meaning that the protein’s function may be affected.
Animals 12 03571 g005
Table
Table 1. Analysis of selective pressure relaxation in the ASIC4 gene.
Table 1. Analysis of selective pressure relaxation in the ASIC4 gene.
Modelsω−lnLModels Compared2Δ(lnL)p Value
One ratio ω(A)0.087−11,012.355
One ratio ω = 1(B)1.000−12,064.363B vs. A2104.0170
The branches with pseudogenized ASIC4 have ω2, others ω1€ω1 = 0.084
ω2 = 0.200		−11,008.339A vs. C8.0310.004
The branches with pseudogenized ASIC4 have ω2 = 1, others ω1(D)ω1 = 0.084
ω2 = 1.000		−11,022.971D vs. C29.2656.31221 × 10−8
Each branch has its own €(E) −10,906.493C vs. E203.6926.15578 × 10−12
0: The p-value approaches 0, but is displayed as 0 because the value is too small.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ding, X.; Yu, F.; He, X.; Xu, S.; Yang, G.; Ren, W. Rubbing Salt in the Wound: Molecular Evolutionary Analysis of Pain-Related Genes Reveals the Pain Adaptation of Cetaceans in Seawater. Animals 2022, 12, 3571. https://doi.org/10.3390/ani12243571

AMA Style

Ding X, Yu F, He X, Xu S, Yang G, Ren W. Rubbing Salt in the Wound: Molecular Evolutionary Analysis of Pain-Related Genes Reveals the Pain Adaptation of Cetaceans in Seawater. Animals. 2022; 12(24):3571. https://doi.org/10.3390/ani12243571

Chicago/Turabian Style

Ding, Xiaoyue, Fangfang Yu, Xiaofang He, Shixia Xu, Guang Yang, and Wenhua Ren. 2022. "Rubbing Salt in the Wound: Molecular Evolutionary Analysis of Pain-Related Genes Reveals the Pain Adaptation of Cetaceans in Seawater" Animals 12, no. 24: 3571. https://doi.org/10.3390/ani12243571

APA Style

Ding, X., Yu, F., He, X., Xu, S., Yang, G., & Ren, W. (2022). Rubbing Salt in the Wound: Molecular Evolutionary Analysis of Pain-Related Genes Reveals the Pain Adaptation of Cetaceans in Seawater. Animals, 12(24), 3571. https://doi.org/10.3390/ani12243571

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop