Evolutionary genetics of pulmonary anatomical adaptations in deep-diving cetaceans.

BACKGROUND: Cetaceans, having experienced prolonged adaptation to aquatic environments, have undergone evolutionary changes in their respiratory systems. This process of evolution has resulted in the emergence of distinctive phenotypic traits, notably the abundance of elastic fibers and thickened alveolar walls in their lungs, which may facilitate alveolar collapse during diving. This structure helps selective exchange of oxygen and carbon dioxide, while minimizing nitrogen exchange, thereby reducing the risk of DCS. Nevertheless, the scientific inquiry into the mechanisms through which these unique phenotypic characteristics govern the diving behavior of marine mammals, including cetaceans, remains unresolved. RESULTS: This study entails an evolutionary analysis of 42 genes associated with pulmonary fibrosis across 45 mammalian species. Twenty-one genes in cetaceans exhibited accelerated evolution, featuring specific amino acid substitutions in 14 of them. Primarily linked to the development of the respiratory system and lung morphological construction, these genes play a crucial role. Moreover, among marine mammals, we identified eight genes undergoing positive selection, and the evolutionary rates of three genes significantly correlated with diving depth. Specifically, the SFTPC gene exhibited convergent amino acid substitutions. Through in vitro cellular experiments, we illustrated that convergent amino acid site mutations in SFTPC contribute positively to pulmonary fibrosis in marine mammals, and the presence of this phenotype can induce deep alveolar collapse during diving, thereby reducing the risk of DCS during diving. CONCLUSIONS: The study unveils pivotal genetic signals in cetaceans and other marine mammals, arising through evolution. These genetic signals may influence lung characteristics in marine mammals and have been linked to a reduced risk of developing DCS. Moreover, the research serves as a valuable reference for delving deeper into human diving physiology.

Association analysis of dopaminergic degeneration and the neutrophil-to-lymphocyte ratio in Parkinson's disease.

Introduction: Peripheral inflammatory responses are suggested to play a major role in the pathophysiology of Parkinson's disease (PD). The neutrophil-to-lymphocyte ratio (NLR), a new recognized biomarker, can reflect peripheral inflammation in PD. However, the association between the NLR and dopaminergic degeneration in PD remains unclear. Methods: In this retrospective study, 101 enrolled PD patients were categorized into early-stage and advanced-stage PD based on the Hoehn and Yahr (HY) scale. We evaluated the clinical characteristics, peripheral immune profile, and 11C-CFT striatal dopamine transporter (DAT) binding levels. Linear regression analyses were employed to assess the associations between NLR and striatal DAT levels at different stages in PD patients. Results: Covariate-controlled regression analysis revealed that higher NLR was significantly associated with lower DAT levels in the caudate (ß = -0.27, p = 0.003) and the putamen (ß = -0.27, p = 0.011). Moreover, in the early-stage PD subgroup, a similar association was observed (caudate: ß = -0.37, p = 0.013; putamen: ß = -0.45, p = 0.005). The lymphocytes count was correlated positively with the striatal DAT levels in the Spearman correlation analysis whether in total patients (caudate: ρ = 0.25, p = 0.013; putamen: ρ = 0.22, p = 0.026) or in the early-stage subgroup (caudate: ρ = 0.31, p = 0.023, putamen: ρ = 0.34, p = 0.011). Conclusion: Dopaminergic degeneration is associated with peripheral inflammation in PD. The NLR, a widely used inflammatory marker, may have the potential to reflect the degree of dopaminergic degeneration in individuals with early-stage PD.

Neuromelanin-sensitive magnetic resonance imaging: Possibilities and promises as an imaging biomarker for Parkinson's disease.

Parkinson's disease (PD) is an age-related progressive neurodegenerative disorder characterized by both motor and non-motor symptoms resulting from the death of dopaminergic neurons in the substantia nigra pars compacta (SNpc) and noradrenergic neurons in the locus coeruleus (LC). The current diagnosis of PD primarily relies on motor symptoms, often leading to diagnoses in advanced stages, where a significant portion of SNpc dopamine neurons has already succumbed. Therefore, the identification of imaging biomarkers for early-stage PD diagnosis and disease progression monitoring is imperative. Recent studies propose that neuromelanin-sensitive magnetic resonance imaging (NM-MRI) holds promise as an imaging biomarker. In this review, we summarize the latest findings concerning NM-MRI characteristics at various stages in patients with PD and those with atypical parkinsonism. In conclusion, alterations in neuromelanin within the LC are associated with non-motor symptoms and prove to be a reliable imaging biomarker in the prodromal phase of PD. Furthermore, NM-MRI demonstrates efficacy in differentiating progressive supranuclear palsy (PSP) from PD and multiple system atrophy with predominant parkinsonism. The spatial patterns of changes in the SNpc can be indicative of PD progression and aid in distinguishing between PSP and synucleinopathies. We recommend that patients with PD and individuals at risk for PD undergo regular NM-MRI examinations. This technology holds the potential for widespread use in PD diagnosis.

Adaptive changes in BMAL2 with increased locomotion associated with the evolution of unihemispheric slow-wave sleep in mammals.

Marine mammals, especially cetaceans, have evolved a very special form of sleep characterized by unihemispheric slow-wave sleep (USWS) and a negligible amount or complete absence of rapid-eye-movement sleep; however, the underlying genetic mechanisms remain unclear. Here, we detected unique, significant selection signatures in basic helix-loop-helix ARNT like 2 (BMAL2; also called ARNTL2), a key circadian regulator, in marine mammal lineages, and identified two nonsynonymous amino acid substitutions (K204E and K346Q) in the important PER-ARNT-SIM domain of cetacean BMAL2 via sequence comparison with other mammals. In vitro assays revealed that these cetacean-specific mutations specifically enhanced the response to E-box-like enhancer and consequently promoted the transcriptional activation of PER2, which is closely linked to sleep regulation. The increased PER2 expression, which was further confirmed both in vitro and in vivo, is beneficial for allowing cetaceans to maintain continuous movement and alertness during sleep. Concordantly, the locomotor activities of zebrafish overexpressing the cetacean-specific mutant bmal2 were significantly higher than the zebrafish overexpressing the wild-type gene. Subsequently, transcriptome analyses revealed that cetacean-specific mutations caused the upregulation of arousal-related genes and the downregulation of several sleep-promoting genes, which is consistent with the need to maintain hemispheric arousal during USWS. Our findings suggest a potential close relationship between adaptive changes in BMAL2 and the remarkable adaptation of USWS and may provide novel insights into the genetic basis of the evolution of animal sleep.

Evaluation of risk factors for impulse control disorder in Parkinson's disease in northern China.

Introduction: Impulse control disorder (ICD) is a common non-motor symptom of Parkinson's disease (PD), but its risk factors are still controversial. This study aimed to determine the prevalence of ICD in northern China and analyze the risk factors associated with ICD, multiple ICDs, and four subtypes. Methods: A total of 285 PD patients were enrolled in this study. Each patient was screened using the Questionnaire for Impulse and Compulsive Control Disorders (QUIP). Stepwise regression analysis was performed to identify independent risk factors, and a prediction model was developed. Results: The prevalence of ICD in the study population was 11.6%. Stepwise regression analysis showed that ICD was associated with disease duration, motor symptoms, dyskinesia, depression, REM sleep behavior disorder (RBD) and cognitive decline; multiple ICDs were related to coffee history, motor symptoms, dyskinesia, depression, apathy and RBD. The prediction model demonstrated good performance with AUC values of 0.93, 0.88, and 0.66 on the balanced train set, balanced test set, and the original imbalanced data set, respectively. Conclusion: The risk factors for PD-ICD are complex and influenced by regional economic and cultural backgrounds. Clarifying these factors and developing predictive models can help to delay or even prevent the development of ICD through early screening and intervention.

Molecular evolution of wound healing-related genes during cetacean secondary aquatic adaptation.

The marine environment presents challenges for wound healing in cetaceans, despite their remarkable recovery abilities with minimal infections or complications. However, the molecular mechanism underlying this efficient wound healing remains underexplored. To better understand the molecular mechanisms behind wound healing in cetaceans, we investigated the evolutionary patterns of 37 wound healing-related genes in representative mammals. We found wound healing-related genes experience adaptive evolution in cetaceans: (1) Three extrinsic coagulation pathway-related genes-tissue factor (F3), coagulation factor VII (F7), and coagulation factor X (F10)-are subject to positive selection in cetaceans, which might promote efficient hemostasis after injury; positive selection in transforming growth factor-beta 2 (TGF-ß2), transforming growth factor-beta 3 (TGF-ß3), and platelet-derived growth factor D (PDGFD), which play immunological roles in wound healing, may help cetaceans enhance inflammatory response and tissue debridement. (2) Coagulation factor XII (F12) is the initiation factor in the intrinsic coagulation pathway. It had a premature stop codon mutation and was subjected to selective stress relaxation in cetaceans, suggesting that the early termination of F12 may help cetaceans avoid the risk of vascular blockage during diving. (3) Fibrinogen alpha chain (FGA) and FIII, which were detected to contain the specific amino acid substitutions in marine mammals, indicating similar evolutionary mechanisms might exist among marine mammals to maintain strong wound-healing ability. Thus, our research provides further impetus to study the evolution of the wound healing system in cetaceans and other marine mammals, extending knowledge of preventing coagulation disorder and atherosclerosis in humans.

Determination of Vitamin K1, MK-4, MK-7, and D Levels in Human Serum of Postmenopausal Osteoporosis Women Based on High Stability LC-MS/MS: MK-7 May Be a New Marker of Bone Metabolism.

INTRODUCTION: Vitamin K (VK) as well as vitamin D (VD) plays an important role in osteoporosis. Vitamin K1 (VK1), vitamin K2 (VK2, menaquinone-4 (MK-4), and menaquinone-7 (MK-7)) are significant for the metabolism of skeletal muscle. 25-hydroxyvitamin D2 and 25-hydroxyvitamin D3 (25(OH)D2 and 25(OH)D3) reflect circulating VD levels. More sensitive measurements remain to be developed. In the present study, a new high-performance liquid chromatography-tandem mass spectrometry (LC-MS/MS) method was developed for the determination of VK1, VK2 (MK-4 and MK-7), as well as 25(OH)D2 and 25(OH)D3 levels in human serum. METHODS: We developed a simple LC-MS/MS method for the determination of VK1, MK-4, MK-7, 25(OH)D2, and 25(OH)D3 levels in human serum and validated the method in a study cohort of 200 patients divided into the premenopausal women group and postmenopausal osteoporosis patient group. RESULTS: The overall precision (coefficient of variation) ranged from 2.66 to 10.11% in the specified working range (0.05-5 ng/mL) for VK1, MK-4, and MK-7. Serum VK1, MK-4, and MK-7 levels in postmenopausal women with osteoporosis were 1.187 ± 0.094 ng/mL, 0.058 ± 0.009 ng/mL, and 0.885 ± 0.064 ng/mL, respectively (mean ± standard deviation). Serum VK1, MK-4, and MK-7 levels in premenopausal women were 1.143 ± 0.103 ng/mL, 0.028 ± 0.003 ng/mL, and 1.553 ± 0.226 ng/mL, respectively. Serum 25(OH)D2 and 25(OH)D3 levels in postmenopausal women with osteoporosis were 0.757 ± 0.056 ng/mL and 11.72 ± 0.632 ng/mL, respectively. Serum 25(OH)D2 and 25(OH)D3 levels in premenopausal were 1.793 ± 0.575 ng/mL and 12.42 ± 1.069 ng/mL, respectively. CONCLUSION: A new LC-MS/MS method for determination of serum VK and VD levels was evaluated and validated. MK-7 in plasma decreased earlier than VD in postmenopausal osteoporosis patients. MK-7 status is significantly associated with osteoporosis and could be considered a predictable biomarker in the diagnosis of osteoporosis in postmenopausal women.

The evolution of the discrete multirenculate kidney in mammals from ecological and molecular perspectives.

Mammals have developed different kinds of renal structures during evolution, yet the origin of the renal structural phenotypes and the molecular mechanisms underlying their adaptive evolution remains unclear. Here, we reconstructed the ancestral state of the renal structures across mammals and found that the unilobar kidney was the ancestral character in mammals. The subsequent correlation analyses between renal phenotypes and life history traits revealed that species with a larger body or in aquatic habitats tend to have evolved discrete multirenculate kidneys. To explore the molecular convergent mechanisms among mammals with this most distinct renal structure, the discrete multirenculate kidney, we used 45 genes related to duplex/multiplex kidney diseases to compare the evolutions of species with discrete multirenculate kidneys and with other renal phenotypes. Twelve rapidly evolving genes that were functionally enriched in cilium assembly and centrosome were identified in species with discrete multirenculate kidneys, suggesting that these genes played key roles in the evolution of discrete multirenculate kidneys. In addition, positive selection was detected in six crucial genes which are mainly involved in epithelial tube morphogenesis and regulation of neurogenesis. Finally, 12 convergent amino acid substitutions, six of which are in crucial domain of proteins, were shared by two or more lineages with discrete multirenculate kidneys. These findings could provide some novel insights into the origin and evolution of renal structures across mammals and the pathogenesis of renal diseases in humans.

Development of a matrix-based candidate reference material for human glycated albumin using isotope-dilution liquid chromatography/tandem mass spectrometry.

Glycated albumin (GA) in human serum is tested clinically as a short-term indicator for glucose monitoring. Here, we evaluated a candidate serum reference material (RM) at three different GA concentrations to help standardize serum GA measurements. Both albumin and GA were quantitatively determined using isotope-dilution liquid chromatography/tandem mass spectrometry with lysine-4,4,5,5-D4·2HCl (D4-lysine) and Nε-l3C6-(l-deoxy-d-fructose-1-yl)-l-lysine (13C6-DOF-lysine) as internal standards and lysine and synthetic DOF-lysine as calibration standards. The method was evaluated with the RM, JCCRM611-1, from the Reference Material Institute for Clinical Chemistry Standards. The homogeneity and stability of the candidate RMs were examined using a commercial biochemical analyzer. Fifteen units were randomly selected, and statistical analysis showed no inhomogeneity. The candidate RMs were stable for at least 6 months at -80 °C. The coefficients of variation (CVs) for the JCCRM611-1 RM ranged from 3.2% to 2.3%, and the biases ranged from 4.12% to -1.84%. GA was tested at low, medium, and high concentrations, which were quantified as 249.53 ± 13.29, 408.02 ± 11.70, and 637.22 ± 17.03 mmol/mol, respectively. The overall CVs ranged from 0.99% to 2.51%. The candidate RMs can potentially be used to develop a traceability chain to improve the accuracy of GA measurements.

Evolutionary genetics of flipper forelimb and hindlimb loss from limb development-related genes in cetaceans.

BACKGROUND: Cetacean hindlimbs were lost and their forelimb changed into flippers characterized by webbed digits and hyperphalangy, thus allowing them to adapt to a completely aquatic environment. However, the underlying molecular mechanism behind cetacean limb development remains poorly understood. RESULTS: In the present study, we explored the evolution of 16 limb-related genes and their cis-regulatory elements in cetaceans and compared them with that of other mammals. TBX5, a forelimb specific expression gene, was identified to have been under accelerated evolution in the ancestral branches of cetaceans. In addition, 32 cetacean-specific changes were examined in the SHH signaling network (SHH, PTCH1, TBX5, BMPs and SMO), within which mutations could yield webbed digits or an additional phalange. These findings thus suggest that the SHH signaling network regulates cetacean flipper formation. By contrast, the regulatory activity of the SHH gene enhancer-ZRS in cetaceans-was significantly lower than in mice, which is consistent with the cessation of SHH gene expression in the hindlimb bud during cetacean embryonic development. It was suggested that the decreased SHH activity regulated by enhancer ZRS might be one of the reasons for hindlimb degeneration in cetaceans. Interestingly, a parallel / convergent site (D42G) and a rapidly evolving CNE were identified in marine mammals in FGF10 and GREM1, respectively, and shown to be essential to restrict limb bud size; this is molecular evidence explaining the convergence of flipper-forelimb and shortening or degeneration of hindlimbs in marine mammals. CONCLUSIONS: We did evolutionary analyses of 16 limb-related genes and their cis-regulatory elements in cetaceans and compared them with those of other mammals to provide novel insights into the molecular basis of flipper forelimb and hindlimb loss in cetaceans.

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

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.

Analysis of Influencing Factors of Apathy in Patients with Parkinson's Disease.

BACKGROUND: Apathy is a common non-motor symptom of Parkinson's disease (PD). The influencing factors of apathy are currently controversial. This study aimed to describe the clinical characteristics of PD-associated apathy and to analyze the associated risk factors. METHODS: Two hundred patients diagnosed with PD were selected. Included patients were divided into an apathetic group and a non-apathetic group. Demographic and clinical data, motor symptoms, non-motor symptoms and medication use of the two groups were assessed. RESULTS: The incidence of apathy was 69%. Demographic and clinical data, motor symptoms, non-motor symptoms and medications use were statistically significant. CONCLUSIONS: PD patients with more severe motor symptoms, cognitive impairment, depression, anxiety, RBD, excessive daytime sleep, fatigue, low education level, long disease course, poor quality of life and lower DA dosage are more prone to apathy. Cognitive function, quality of life, educational level, DA and LEDD are independent risk factors for apathy.

Protein and miRNA profile of circulating extracellular vesicles in patients with primary sclerosing cholangitis.

Primary sclerosing cholangitis (PSC) is an idiopathic and heterogenous cholestatic liver disease characterized by chronic inflammation and fibrosis of the biliary tree. Currently, no effective therapies are available for this condition, whose incidence is rising. At present, specificity and sensitivity of current serum markers used to diagnose PSC are limited and often unreliable. In this study, we characterize circulating extracellular vesicles and provide supporting data on their potential use as novel surrogate biomarkers for PSC. EVs are membrane surrounded structures, 100-1000 nm in size, released by cells under various conditions and which carry a variety of bioactive molecules, including small non-coding RNAs, lipids and proteins. In recent years, a large body of evidence has pointed to diagnostic implications of EVs and relative cargo in various human diseases. We isolated EVs from serum of well-characterized patients with PSC or control subjects by differential centrifugation and size-exclusion chromatography. A complete characterization identified elevated levels of circulating EVs in PSC patients compared to healthy control subjects (2000 vs. 500 Calcein-FITC + EVs/µL). Tissue and cell specificity of circulating EVs was assessed by identification of liver-specific markers and cholangiocyte marker CK-19. Further molecular characterization identified 282 proteins that were differentially regulated in PSC-derived compared to healthy control-EVs. Among those, IL-13Ra1 was the most significantly and differentially expressed protein in PSC-derived EVs and correlated with the degree of liver fibrosis. In addition to protein profiling, we performed a miRNA-sequencing analysis which identified 11 among established, liver-specific (e.g., miR-122 and miR-192) and novel miRNAs. One of the newly identified miRNAs, miR-4645-3p, was significantly up-regulated fourfold in PSC-derived EVs compared to circulating EVs isolated from healthy controls. This study provides supporting evidence of the potential role of circulating EVs and associated protein and miRNA cargo as surrogate noninvasive and reliable biomarker for PSC.

Insights into the Evolution of Spermatogenesis-Related Ubiquitin-Proteasome System Genes in Abdominal Testicular Laurasiatherians.

During embryonic development in mammals, the testicles generally descend into the scrotum, making the testicular temperature 2-4 °C lower than the core temperature via heat exchange and clearance, and thus more beneficial for normal spermatogenesis. Failure to descend, known as cryptorchidism, carries a series of risks such as infertility and testicular cancer. However, some mammals have evolved abdominal testes while maintaining healthy reproduction. To explore the underlying molecular mechanism, we conducted comparative genomic analyses and functional assays on the spermatogenesis-related ubiquitin-proteasome system (UPS) genes essential to sperm formation in representative laurasiatherians. Here, positive selection and rapid evolution of spermatogenesis-related UPS genes were identified in the abdominal testicular laurasiatherians. Moreover, potential convergent amino acids were found between distantly related species with similar abdominal testicles and functional analyses showed RNF8 (V437I) in abdominal testicular species (437I) has a stronger ubiquitination ability, which suggests that the mammals with abdominal testes might exhibit enhanced sperm cell histone clearance to maintain sperm formation. This evidence implies that, in response to "cryptorchidism injury", spermatogenesis-related UPS genes in the abdominal testicular species might have undergone adaptive evolution to stabilize sperm formation. Thus, our study could provide some novel insights into the reproductive adaptation in abdominal testicular mammals.

Comparative genomics reveals molecular mechanisms underlying health and reproduction in cryptorchid mammals.

BACKGROUND: Mammals have wide variations in testicular position, with scrotal testes in some species and ascrotal testes in others. Although cryptorchidism is hazardous to human health, some mammalian taxa are natural cryptorchids. However, the evolution of testicular position and the molecular mechanisms underlying the maintenance of health, including reproductive health, in ascrotal mammals are not clear. RESULTS: In the present study, comparative genomics and evolutionary analyses revealed that genes associated with the extracellular matrix and muscle, contributing to the development of the gubernaculum, were involved in the evolution of testicular position in mammals. Moreover, genes related to testicular position were significantly associated with spermatogenesis and sperm fertility. These genes showed rapid evolution and the signature of positive selection, with specific substitutions in ascrotal mammals. Genes associated with testicular position were significantly enriched in functions and pathways related to cancer, DNA repair, DNA replication, and autophagy. CONCLUSIONS: Our results revealed that alterations in gubernaculum development contributed to the evolution of testicular position in mammals and provided the first support for two hypotheses for variation in testicular position in mammals, the "cooling hypothesis", which proposes that the scrotum provides a cool environment for acutely heat-sensitive sperm and the "training hypothesis", which proposes that the scrotum develops the sperm by exposing them to an exterior environment. Further, we identified cancer resistance and DNA repair as potential protective mechanisms in natural cryptorchids. These findings provide general insights into cryptorchidism and have implications for health and infertility both in humans and domestic mammals.

Comparative analyses of aging-related genes in long-lived mammals provide insights into natural longevity.

Extreme longevity has evolved multiple times during the evolution of mammals, yet its underlying molecular mechanisms remain largely underexplored. Here, we compared the evolution of 115 aging-related genes in 11 long-lived species and 25 mammals with non-increased lifespan (control group) in the hopes of better understanding the common molecular mechanisms behind longevity. We identified 16 unique positively selected genes and 23 rapidly evolving genes in long-lived species, which included nine genes involved in regulating lifespan through the insulin/IGF-1 signaling (IIS) pathway and 11 genes highly enriched in immune-response-related pathways, suggesting that the IIS pathway and immune response play a particularly important role in exceptional mammalian longevity. Interestingly, 11 genes related to cancer progression, including four positively selected genes and seven genes with convergent amino acid changes, were shared by two or more long-lived lineages, indicating that long-lived mammals might have evolved convergent or similar mechanisms of cancer resistance that extended their lifespan. This suggestion was further corroborated by our identification of 12 robust candidates for longevity-related genes closely related to cancer.

Over-expression of the bottlenose dolphin Hoxd13 gene in zebrafish provides new insights into the cetacean flipper formation.

Cetaceans have evolved elongated soft-tissue flipper with digits made of hyperphalangy. Cetaceans were found to have 2-3 more alanine residues in Hoxd13 than other mammals, which were suggested to be related to their flipper. However, how Hoxd13 regulates other genes and induces hyperphalangy in cetaceans remain poorly understood. Here, we overexpressed the bottlenose dolphin Hoxd13 in zebrafish (Danio rerio). Combined with transcriptome data and evolutionary analyses, our results revealed that the Wingless/Integrated (Wnt) and Hedgehog signaling pathways and multiple genes might regulate hyperphalangy development in cetaceans. Meanwhile, the Notch and mitogen-activated protein kinase (Mapk) signaling pathways and Fibroblast growth factor receptor 1 (Fgfr1) are probably correlated with interdigital tissues retained in the cetacean flipper. In conclusion, this is the first study to use a transgenic zebrafish to explore the molecular evolution of Hoxd13 in cetaceans, and it provides new insights into cetacean flipper formation.