Research article

Bioinformatic analysis of the coding region of the melatonin receptor 1b gene as a reliable DNA marker to resolve interspecific mammal phylogenetic relationships

  • Received: 04 November 2022 Revised: 19 December 2022 Accepted: 22 December 2022 Published: 13 January 2023
  • This research looks into the main DNA markers and the limits of their application in molecular phylogenetic analysis. Melatonin 1B (MTNR1B) receptor genes were analyzed from various biological sources. Based on the coding sequences of this gene, using the class Mammalia as example, phylogenetic reconstructions were made to study the potential of mtnr1b as a DNA marker for phylogenetic relationships investigating. The phylogenetic trees were constructed using NJ, ME and ML methods that establish the evolutionary relationships between different groups of mammals. The resulting topologies were generally in good agreement with topologies established on the basis of morphological and archaeological data as well as with other molecular markers. The present divergences provided a unique opportunity for evolutionary analysis. These results suggest that the coding sequence of the MTNR1B gene can be used as a marker to study the relationships of lower evolutionary levels (order, species) as well as to resolve deeper branches of the phylogenetic tree at the infraclass level.

    Citation: Ekaterina Y. Kasap, Оlga K. Parfenova, Roman V. Kurkin, Dmitry V. Grishin. Bioinformatic analysis of the coding region of the melatonin receptor 1b gene as a reliable DNA marker to resolve interspecific mammal phylogenetic relationships[J]. Mathematical Biosciences and Engineering, 2023, 20(3): 5430-5447. doi: 10.3934/mbe.2023251

    Related Papers:

    [1] Arianna Polverino, Pierpaolo Sorrentino, Matteo Pesoli, Laura Mandolesi . Nutrition and cognition across the lifetime: an overview on epigenetic mechanisms. AIMS Neuroscience, 2021, 8(4): 448-476. doi: 10.3934/Neuroscience.2021024
    [2] Rosa Giacca, Miriana Conte, Alessandro d'Ambrosio, Alvino Bisecco, Renato Docimo, Mario Risi, Manuela Altieri, Riccardo Borgo, Rosario Domenico Melisi, Vittorio Gentile, Antonio Gallo . Use of Transglutaminase 2 mRNA expression in peripheral blood mononuclear cells in patients with Radiologically Isolated Syndrome as a neuroinflammation biomarker: A preliminary study. AIMS Neuroscience, 2025, 12(2): 284-290. doi: 10.3934/Neuroscience.2025015
    [3] Martina Valencia, Odra Santander, Eloísa Torres, Natali Zamora, Fernanda Muñoz, Rodrigo Pascual . Environmental enrichment reverses cerebellar impairments caused by prenatal exposure to a synthetic glucocorticoid. AIMS Neuroscience, 2022, 9(3): 320-344. doi: 10.3934/Neuroscience.2022018
    [4] Mary B. Dratman, Joseph V. Martin . The many faces of thyroxine. AIMS Neuroscience, 2020, 7(1): 17-29. doi: 10.3934/Neuroscience.2020002
    [5] Sara Arciniegas Ruiz, Eliav Tikochinsky, Vardit Rubovitch, Chaim G Pick, Bernard Attali . Contextual fear response is modulated by M-type K+ channels and is associated with subtle structural changes of the axon initial segment in hippocampal GABAergic neurons. AIMS Neuroscience, 2023, 10(1): 33-51. doi: 10.3934/Neuroscience.2023003
    [6] Mohammad Mofatteh . mRNA localization and local translation in neurons. AIMS Neuroscience, 2020, 7(3): 299-310. doi: 10.3934/Neuroscience.2020016
    [7] Mi Zhang, Junjie Ren, Ni Li, Yongyi Li, Linxi Yang, Wenzhuo Wei, Juan Qiu, Xiaochu Zhang, Xiaoming Li . Hypnosis efficacy on nicotine addiction: An analysis of EEG microstates and brain oscillation entropy. AIMS Neuroscience, 2025, 12(1): 15-31. doi: 10.3934/Neuroscience.2025002
    [8] Chacchu Bhattarai, Phanindra Prasad Poudel, Arnab Ghosh, Sneha Guruprasad Kalthur . The RET gene encodes RET protein, which triggers intracellular signaling pathways for enteric neurogenesis, and RET mutation results in Hirschsprung's disease. AIMS Neuroscience, 2022, 9(1): 128-149. doi: 10.3934/Neuroscience.2022008
    [9] Khue Vu Nguyen . Encephalomalacia/gliosis, deep venous thrombosis, and cancer in Arg393His antithrombin Hanoi and the potential impact of the β-amyloid precursor protein (APP) on thrombosis and cancer. AIMS Neuroscience, 2022, 9(2): 175-215. doi: 10.3934/Neuroscience.2022010
    [10] Chris Cadonic, Benedict C. Albensi . Oscillations and NMDA Receptors: Their Interplay Create Memories. AIMS Neuroscience, 2014, 1(1): 52-64. doi: 10.3934/Neuroscience.2014.1.52
  • This research looks into the main DNA markers and the limits of their application in molecular phylogenetic analysis. Melatonin 1B (MTNR1B) receptor genes were analyzed from various biological sources. Based on the coding sequences of this gene, using the class Mammalia as example, phylogenetic reconstructions were made to study the potential of mtnr1b as a DNA marker for phylogenetic relationships investigating. The phylogenetic trees were constructed using NJ, ME and ML methods that establish the evolutionary relationships between different groups of mammals. The resulting topologies were generally in good agreement with topologies established on the basis of morphological and archaeological data as well as with other molecular markers. The present divergences provided a unique opportunity for evolutionary analysis. These results suggest that the coding sequence of the MTNR1B gene can be used as a marker to study the relationships of lower evolutionary levels (order, species) as well as to resolve deeper branches of the phylogenetic tree at the infraclass level.



    The apelin receptor (APLNR) is a type of G-protein-coupled receptor (GPCR), which has garnered significant attention since its sequence similarity to the angiotensin (AT1) receptor was discovered in 1993 [1]. Apelin receptor, coupled with inhibitory G proteins, inhibits adenylate cyclase activity, counteracting the vasoconstrictive effects of angiotensin II and, consequently, modulating blood pressure [2]. Apelin, an APLNR agonist, was discovered in 1998 in bovine stomach extracts [3]. Since then, the intracellular signaling of the apelin/APLNR system has been well studied. This system stimulates signaling pathways, such as phosphoinositide 3-kinase (PI3K)/protein kinase B (Akt) cascades and mitogen-activated protein kinase (MAPK) cascades, which are involved in various cellular signaling pathways, including proliferation and cell death. Akt activation occurs through phosphorylation at serine 473 (S473) by the mammalian target of rapamycin complex 2 (mTORC2), followed by phosphorylation at threonine 308 (T308) by phosphoinositide-dependent protein kinase 1 (PDK1) [4]. Extracellular signal-regulated kinases (ERK)1/2 (p44/42) are members of MAPKs and are activated by upper kinases called MAP kinase kinases via phosphorylation at T202 and Tyrosine 204 (Y204) on ERK1 or T185 and Y187 on ERK2 [5]. In addition to its role in blood pressure regulation [6], the apelin/APLNR system is involved in various mechanisms, including angiogenesis [7], developmental normalization [8],[9] and neuroprotection [10]. Among these, we focused on the neuroprotective effects of the apelin/APLNR system. Cortical neurons transfected with APLNR exhibit neuroprotective properties against NMDA excitotoxicity by activating PI3K/Akt and Raf/ERK cascades [11]. Furthermore, in vivo studies involving tissues such as the mouse hippocampus [11] and retina [12] have demonstrated the neuroprotective effects of APLNR activation in vivo. Pretreatment with apelin reduces neuronal death in the ischemic mouse brain [13]. Evidence suggests that this system may be a potential therapeutic target for neurodegenerative diseases. In our previous study, Apelin knockout mice exhibited accelerated disease progression in amyotrophic lateral sclerosis (ALS) model, whereas increased survival of APLNR-expressing motor neurons was observed in ALS mice [13], suggesting that the APLNR is involved in the maintenance of motor neuron survival. Although the function of the APLNR has been elucidated, its wide range of effective concentrations remain a pharmacological enigma [14],[15]. The factors contributing to this require clarification prior to the use of APLNR as a therapeutic target.

    Post-translational modifications (PTMs) are responsible for the activation or termination of protein function [16]. More than 500 types of PTMs have been reported and registered in the unimod database (http://www.unimod.org/). Among them, phosphorylation is a well-known PTM that is essential for the activation or termination of enzymatic reactions and signal transduction. Glycosylation profoundly affects protein conformation, including folding, localization and stability. Ubiquitination of lysine residues promotes proteasomal degradation, whereas the addition of a small ubiquitin-related modifier (SUMO, i.e., SUMOylation) to lysine residues protects against degradation [17]. The function, localization and stability of GPCRs are also affected by PTMs, including phosphorylation, glycosylation, nitrosylation, tyrosine sulfation, palmitoylation, ubiquitination and SUMOylation [19]. Phosphorylation and glycosylation of the APLNR have been reported and predicted PTMs on the APLNR include SUMOylation and palmitoylation [20]. Although PTMs in most receptors can be inferred from databases and in silico approaches, the precise mechanisms underlying these modifications remain to be elucidated. The addition of functional groups, such as phosphorylation and palmitoylation does not significantly change the molecular weight of proteins, whereas isopeptidation, including ubiquitination and SUMOylation, or glycosylation with a polysaccharide chain, inevitably causes distinct changes in molecular weight. In this study, we investigated the precise signaling mechanism of the APLNR by overexpressing it and explored the potential PTMs involved. We seek to gain a comprehensive understanding of the pharmacological properties of these receptors. We specifically addressed the possibility that certain PTMs of the APLNR, which were not detected in vitro, may stabilize this protein and contribute to its function in vivo.

    (Pyr1)-Apelin-13 was purchased from Peptide Institute (Osaka, Japan). Dulbecco's modified Eagle's medium (DMEM), Blocking One Histo and FluoroKEEPER antifade reagent were purchased from Nacalai Tesque (Kyoto, Japan). ISOGEN II RNA extraction reagent (Nippon Gene, Tokyo, Japan) and ReverTra Ace qPCR RT master mix were provided by Toyobo (Osaka, Japan). Forskolin was supplied by Merk-Sigma (Burlington, MA); polyvinylidene difluoride (PVDF) membrane was obtained from Merk-Millipore (Burlington, MA); Western Lightning Plus-ECL was purchased from PerkinElmer (Waltham, MA); protein A-sepharose 4B resin and EZ-Link NHS-LC-Biotin and Pierce NeutrAvidin UltraLink Resin were purchased from Thermo Scientific (Waltham, MA). PNGase F was provided by New England BioLabs (Ipswich, MA); USP2, Den1 and SENP1 were purchased from LifeSensors (Malvern, PA). All other chemicals used were of the highest commercially available purity.

    Rabbit IgG anti-phospho-PKA substrate (RRXS*/T*) (100G7E), 1:1000 (#9621); anti-Akt (pan), 1:2000 (#3673); anti-phospho-Akt (T308), 1:1000 (#2965); phospho-Akt (S473), 1:2000 (#4060); anti-ERK1/2, 1:2000 (#4695); phospho-ERK1/2 (Thr202/Tyr204), 1:2000 (#4370); Alexa Fluor® 555-conjugated goat IgG anti-mouse IgG (H+L) and F(ab')2 Fragment were purchased from Cell Signaling Technology (Danvers, MA). Mouse IgG anti-human APJ (APLNR) was obtained from R&D systems (Minneapolis, MN), 1:1000 (#MAB8561). Mouse IgG anti-GAPDH, 1:10000 (#014-25524) and anti-FLAG were from FUJIFILM Wako Pure Chemical Corporation (Osaka, Japan), I.F. 1:200, W.B. 1:1000 I.P. 1:500 (#014-22383). Rabbit IgG anti-myelin protein zero was from Abcam (Cambridge, UK), 1:1000 (#ab31851).

    All procedures were performed according to the guidelines of the Japanese Society for Pharmacology and approved by the Committee for Ethical Use of Experimental Animals at Setsunan University. Adult (5–6 weeks old) male C57BL/6 mice were housed at our facilities under a light–dark cycle of 12–12 hr and a humidity of 55% at 23 °C and provided with ad libitum access to food and water.

    HEK293TN cells (System Bioscience, CA, USA) were cultured in DMEM containing 10% fetal bovine serum and 1% penicillin-streptomycin at 37 °C with 5% CO2. Cells reached about 80% confluency and were removed from the plate by Trypsin/EDTA and re-seeded at a density of 1.0 × 105 cells/ml. Under these conditions, cells needed three days to reach about 80% confluency. During the culture period, drugs were added directly to the culture medium. Cells were exposed to (Pyr1)-apelin-13 for 30 min. In some cases, 15 min after the initiation of (Pyr1)-apelin-13 exposure, FSK exposure was performed for 15 min. Cells were exposed to tunicamycin at 1 day in vitro (DIV) for 48 h.

    pcDNA3.1+-DYK-human APLNR (OHu11900) expressing the human APLNR protein with an N-terminal DYKDDDDK (FLAG) tag in mammalian cells (pcDNA-APLNR) was purchased from GenScript Biotech Corporation (Piscataway, NJ). An empty vector, pCDH-CMV-MCS-EF1-GreenPuro (CD513B-1; pCDH-GFP), was purchased from System Bioscience (CA, USA), and the FLAG-tagged APLNR gene from the abovementioned plasmid was subcloned (pCDH-GFP-APLNR). Mammalian cells transfected with this plasmid express APLNR and GFP under the control of the CMV and EF1 promoters, respectively.

    The calcium phosphate transfection method was performed by mixing 1.0 µg of a plasmid containing HEPES-buffered saline and CaCl2 at a final concentration of 86 mM. The DNA mixture was allowed to stand at room temperature for 20 min and then added to the cell culture in a 35 mm dish with 2.0 mL of culture media. Subsequently, the cells were placed in a CO2 incubator for 2 days. For mock transfections, an equal volume of sterile purified water was used instead of the plasmid.

    Total RNA extracted using the ISOGEN II was treated at 65 °C for 5 min and cooled on ice. 1 µg of RNA as template was reacted with reverse transcription using ReverTra Ace qPCR RT master mix to prepare cDNA through at 37 °C for 15 min and 98 °C for 5 min. Gene transcripts were amplified using specific primer sets, as listed below.

    APLNR (207 bp) Forward: 5′-TTCTGCAAGCTCAGCAGCT-3′
    Reverse: 5′-GGTGCGTAACACCATGACAG-3′
    GAPDH (193 bp) Forward: 5′-TCATTGACCTCAACTACATGGTCTA-3′
    Reverse: 5′-ACACCAGTAGACTCCACGACATACT-3′

    Briefly, RT-PCR was performed in three steps. The first step consisted of 1 cycle at 94 °C for 120 s. The second step included 30 cycles at 94 °C for 45 s, 55 °C for 45 s and 72 °C for 30 s. The third step consisted of 1 cycle at 72 °C for 60 s. The PCR products were separated by electrophoresis using 0.8% of TAE agarose gel and visualized using the ImageQuant 400 (GE electronics, Boston, MA).

    Immunoblotting was performed as previously described [21]. Cells were lysed using a buffer composed of 10 mM Tris-HCl (pH 7.5), 0.32 M sucrose, 1 mM EDTA, 1 mM EGTA, 5 mM dithiothreitol, phosphatase inhibitors [10 mM sodium b-glycerophosphate and 10 mM sodium pyrophosphate (PPA), 50 mM sodium fluoride (NaF)] and 1 µg/mL each of protease inhibitors ([p-amidinophenyl] methanesulfonyl fluoride, benzamidine, leupeptin and antipain). In some cases, protease inhibitors were omitted. Proteins in the lysate were solubilized at 100 °C for 10 min, 37 °C for 10, 30 or 60 min or 4 °C for 10 min in a buffer containing 2% sodium dodecyl sulfate (SDS), 5% 2-mercaptoethanol, 10% glycerol and 0.01% bromophenol blue. Unless otherwise specified in the method, the solubilization was performed at 100 °C. Subsequently, 10 µg of proteins were subjected to 10% SDS–polyacrylamide gel electrophoresis and transferred onto PVDF membranes. To block non-specific binding, the membranes were immersed in Tris-buffered saline (TBS) containing 5% skim milk and 0.05% Tween-20 for 1 hr, followed by incubation with a specific primary antibody for 2 hr and with horseradish peroxidase-conjugated secondary antibody for 1 hr at room temperature. Membranes were subjected to detection using Western Lightning Plus-ECL and exposed to X-ray films (Fujifilm, Tokyo, Japan) for visualization. After detecting total protein, the membrane for detecting phosphorylation was soaked in WB stripping solution (Nacalai # 05364-55), gently shaken at room temperature for 30 minutes and used again for antibody reaction.

    Cells on a coverslip were fixed with 4% paraformaldehyde in TBS for 10 min. To enhance membrane permeability, cells were treated with TBS containing 0.1% Triton X-100 for 10 min. After blocking with Blocking One Histo for 10 min, a specific primary antibody was applied in TBS with 0.1% Blocking One and 0.03% Tween-20 and incubated overnight at 4 °C. Cells were then incubated for 2 hr with appropriate combinations of secondary antibodies. After an additional wash with TBS, a coverslip containing immunostained cells was mounted using FluoroKEEPER antifade reagent. Immunostaining images were observed under a confocal laser-scanning microscope (LSM 510; Carl Zeiss, Jena, Germany).

    PNGase F, USP2, Den1 and SENP1 were used to remove N-linked glycosylation, ubiquitination, NEDDylation or SUMOylation, respectively, according to the manufacturer's instructions. Briefly, when the protein required denaturing before reacting with the enzyme, 10 µg of protein was mixed with denaturing buffer containing 0.5% SDS and 40 mM DTT and incubated at 100 °C for 10 min. Subsequently, 10 µg of proteins were mixed with 50 mM PBS containing 1% NP-40 and incubated with 500 units/10 nmol PNGase F, 35 mmol USP2, 1 µmol Den1 or 10 nmol SENP1 at 37 °C for 1 h. The proteins from HEK293TN cells or the mouse spinal cord were solubilized at 37 °C or 100 °C for 10 min respectively and subjected to immunoblotting. In some cases, 100 µg/mL each of protease inhibitors ([p-amidinophenyl]methanesulfonyl fluoride, benzamidine, leupeptin and antipain) were added to the reaction mixture.

    Cells were lysed in RIPA buffer composed of 50 mM Tris-HCl (pH 8.0), 5 mM EGTA, 5 mM EDTA, 50 mM NaF, 10 mM PPA, 1% nonidet P-40 (NP-40), 0.5% deoxycholate and 0.1% SDS. The lysate containing 100 µg of proteins was incubated with 1 µg of anti-FLAG antibody and protein A-sepharose 4B resin by rotation mixer at 4 °C for 120 min. The bound substances were precipitated by centrifugation, and the precipitant was solubilized at 37 °C for 10 min and subjected to immunoblotting.

    Cells were incubated with EZ-Link NHS-LC-Biotin at 4 °C for 60 min and lysed in RIPA buffer. Subsequently, 250 µg of the biotinylated membrane proteins were incubated with Pierce NeutrAvidin UltraLink Resin by rotation mixer at 4 °C for 120 min. The bound substances were precipitated by centrifugation, and the precipitant was solubilized at 37 °C for 10 min and subjected to immunoblotting.

    Densitometric analysis of immunoreactive bands was performed using ImageJ (National Institutes of Health). All data were expressed as means ± SE, and statistical significance was determined using Bonferroni's test.

    The APLNR gene expression was observed when the expression vector was transfected into HEK293TN cells, while GAPDH gene expression was observed under all conditions (Fig. 1a, b). Forskolin (FSK) significantly increased the phosphorylation levels of several protein kinase A (PKA) substrates in both mock and APLNR-expressing cells, whereas (Pyr1)-apelin-13 did not affect phosphorylation levels (Fig. 2a). Furthermore, to investigate the expression of APLNR whether there was a change in its expression with treatment (Pyr1)-apelin-13 or FSK, the results showed that (Pyr1)-apelin-13 and FSK had no effect on APLNR expression, suggesting no membrane expression (Fig. 2b). GPCR-mediated activation of ERK has been shown to be facilitated by G protein-coupled receptor kinase (GRK) and β-arrestin, which are involved in the GPCR biased pathway [14]. In cortical neurons, APLNR activation stimulates PI3K/Akt and Raf/ERK cascades [11]. Relative phosphorylation levels of ERK1/2 were significantly higher in APLNR gene-transfected cells than in mock cells. However, (Pyr1)-apelin-13 did not alter the protein expression and phosphorylation levels of Akt and ERK1/2 (Fig. 3). These results suggest that APLNR protein expression or translocation to the cell membrane may not occur.

    Figure 1.  Overexpression of human APLNR in HEK293TN cells. Cells were cultured without any treatment (none) or transfected with a control (mock) or APLNR expression plasmid, pcDNA3.1-APLNR (APLNR), immediately after plating. After 3 days of transfection, cDNA from each condition was prepared, and a) APLNR or b) GAPDH gene expression was examined by PCR. Black and white arrowheads indicate the desired base pair of APLNR or GAPDH.
    Figure 2.  No significant alteration in phosphorylation levels of PKA substrates following apelin treatment. Cells transfected with a control (mock) or APLNR expression plasmid, pcDNA3.1-APLNR (APLNR), were collected after treatment with or without 1 µM (Pyr1)-apelin-13 for 30 min and/or 25 µM FSK for the last 15 min. a) Cell lysates were subjected to immunoblotting to detect the phosphorylated PKA substrates, GAPDH and APLNR. The black arrowhead indicates the desired molecular weight of the APLNR (42.6 kDa). b) In each lane of a), all band intensities for the phosphorylated PKA substrates or each band intensities for APLNR and GAPDH were quantified, and the ratio of densitometric units (phosphorylated PKA substrates or APLNR/GAPDH) was calculated. The “untreated” group is a sample that has not been treated with any drug, and the “negative (−)” group is a negative control that was added a solvent of drugs. The experiment was repeated three to four times. N.D. non-detectable. N.S. not significant. **P < 0.01 indicates significant differences between groups without FSK and with FSK (Bonferroni's test).
    Figure 3.  No significant alteration in total expression and phosphorylation levels of Akt or ERK following apelin treatment. Cells transfected with a control (mock) or APLNR expression plasmid, pcDNA3.1-APLNR (APLNR), were collected after treatment with or without 1 µM (Pyr1)-apelin-13 for 30 min. Cell lysates were subjected to immunoblotting to detect the total expression and phosphorylation levels of a) Akt or b) ERK1/2. GAPDH was used as an internal control. The experiment was repeated three or four times. The relative phosphorylation level was calculated by dividing the band intensity of each phosphorylation level by the total protein expression. *P < 0.05 and **P < 0.01 indicate significant differences between mock and APLNR transfections for each corresponding condition (Bonferroni's test).

    Using the standard sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) method, we were unable to detect the APLNR protein after heating at 100 °C for 10 min (Fig. 4; third lane from the left). However, upon solubilization at a lower temperature of 37 °C or 4 °C, two bands of approximately 42.6 kDa appeared, corresponding to the desired molecular weight of the APLNR. After incubation at 37 °C for 60 min and 100 °C for 10 min, no bands were observed (Fig. 4). Receptor proteins with a high proportion of hydrophobic amino acid residues may form insoluble aggregates when heated during the SDS solubilization process, and lower temperatures can be more effective in such cases [22]. Apelin receptor may exhibit similar behavior. Treatment with tunicamycin or deglycosidase attenuated the signal intensity of the upper band, which exhibited slower mobility in electrophoresis. However, protease inhibitors did not affect the two bands (Fig. 5). To confirm its subcellular localization, the APLNR gene was subcloned into a plasmid co-expressing GFP. Similarly, in the subcloned plasmid, APLNR bands appeared after solubilization at 37 °C, and the upper band exhibited slightly slower migrating, which intensified and converted to a faster migrating band after tunicamycin treatment. Immunoprecipitation enhanced the intensity of APLNR bands (Fig. 6a). Furthermore, immunostaining revealed the predominant intracellular distribution of APLNR (Fig. 6c). Concentration of the upper band with slow mobility, detected in the biotinylated membrane proteins, using avidin beads resulted in decreased band intensity (Fig. 6b). These findings suggest that although the APLNR protein is overexpressed in cells, only a small amount of protein translocates to the plasma membrane.

    Figure 4.  Solubilization of APLNR in HEK293TN cells. Cells were either cultured without any treatment (None) or transfected with a control (Mock) or APLNR expression plasmid, pcDNA3.1-APLNR (APLNR), for 3 days immediately after plating. Upon protein solubilization by SDS-PAGE, the cell lysates were incubated at 37 °C, 100 °C or 4 °C for 10, 30 or 60 min. GAPDH was used as an internal control. The black arrowhead indicates the desired molecular weight of the APLNR.
    Figure 5.  Inhibition of glycosylation and deglycosylation of APLNR. Left) After exposure to tunicamycin for 48 hr at indicated concentrations, proteins in HEK293TN cells expressing APLNR were solubilized at 37 °C for 10 min. Right) HEK293TN cells transfected with pcDNA3.1-APLNR were lysed and then reacted in the presence of proteolytic enzyme inhibitors or a de-N-glycosylating enzyme, PNGase F. Thereafter, the reaction solution was solubilized at 37 °C for 10 min. GAPDH was detected as an internal control. The arrowhead indicates the desired molecular weight of APLNR.
    Figure 6.  Cellular membrane localization of APLNR in HEK293TN cells. HEK293TN cells were transfected with the pCDH-GFP-APLNR vector. a) The cells were cultured with or without 5 µg/ml tunicamycin for 48 h. Lysates containing 100 µg of each protein were immunoprecipitated with an anti-FLAG antibody, solubilized at 37 °C for 10 min (IP: FLAG) and subjected to immunoblotting using an anti-APLNR antibody. b) Cells were incubated with cellular membrane-impermeant biotin to biotinylate membrane proteins, and a lysate containing 250 µg of biotinylated proteins was precipitated using Neutravidin-beads and solubilized at 37 °C for 10 min (Mem). a, b) The precipitates were subjected to immunoblotting using an anti-FLAG antibody. PC: positive control, cell lysate transfected with pcDNA3.1-APLNR. For the total lysate (input) or PC, 10 µg protein aliquots were used. b) Input and Mem of non-transfected HEK293TN cells were prepared as negative controls (none). c) Cells were fixed and subjected to immunostaining using a FLAG antibody to detect APLNR expression. Bar: 20 µm. A negative control experiment with an empty vector, pCDH-GFP, was also prepared (empty).

    Immunoblot analysis showed the presence of bands corresponding to the molecular weight of APLNR and Protein0, markers of Schwann cells, in all mice spinal cord fractions, particularly in the sacral spinal cord fraction. The bands of the antibody against APLNR were approximately 10 kDa higher than the expected molecular weight (Fig. 7a). Solubilization of proteins from the lumbar spinal cord at a lower temperature of 37 °C resulted in the observation of multiple bands, but none of them corresponded to the desired molecular weight (Fig. 7b). In vitro reactions that cleave PTMs often require prior denaturation of the target protein. However, regardless of whether the protein was denatured or not, deglycosylation of APLNR occurred, causing a slight shift in the band migration towards the lower side with increased mobility. Notably, the deglycosylated band still appeared at a position higher than the expected molecular weight. However, regardless of the degeneration process, the de-ubiquitinating, de-NEDDylation or de-SUMOylation enzymes did not reveal any changes in the band patterns (Fig. 8).

    Figure 7.  Property for APLNR in the mouse spinal cord. Representative immunoblot against APLNR and Protein0 in the spinal cord from three mice. a) Mouse cervical (C), thoracic (T), lumbar (L) and sacral (S) spinal cords were dissected, and the resulting proteins were solubilized at 100 °C and subjected to immunoblotting to detect APLNR, Protein0 and GAPDH. b) Mouse lumbar spinal cords were dissected, and proteins were solubilized at lower (37 °C) or normal (100 °C) temperatures. Lysate of HEK293TN transfected with pcDNA3.1-APLNR solubilized at 37 °C was used as a positive control (PC). White and black arrowheads indicate the main band in the spinal cord and the desired molecular weight of APLNR, respectively.
    Figure 8.  Possible isopeptidation of APLNR in the mouse spinal cord. Mouse spinal cords were dissected and lysed. Lysates were incubated either with or without denaturation and/or enzymes to introduce a) de-glycosylation by PNGaseF, b) de-ubiquitination by USP2, c) de-neddylation by Den1 or d) de-SUMOylation by SENP1 under appropriate enzyme conditions, as described in the Materials and Methods. After the reaction, APLNR and GAPDH were detected using immunoblotting. n = 3 biologically independent mice.

    The purpose of our study was to investigate and elucidate the intracellular signaling mechanisms of the APLNR. Previous studies have suggested that elucidating PTMs of the APLNR is necessary for this purpose. Although some receptors can be functionally characterized by expressing them in HEK293 cells [23], this was not the case for APLNR. Multiple studies have mentioned that the thermosensitive degradation of the APLNR needs to be addressed, but our study is the first to suggest that PMTs may be involved in resolving this issue. Previous studies have reported the instability of APLNR expressed using artificial methods such as gene transfection [24]. Wild-type APLNR has low protein yields and requires mutations for structural studies using X-ray crystallography [25]. This study showed that APLNR protein remains stable at or below 37 °C but denatured at higher temperatures, rendering it undetectable by antibodies (Fig. 4). Hydrophobic proteins such as membrane receptors, aggregate and hold when heated, but aggregation is suppressed at lower temperatures [22]. These differences may be attributed to PTMs. Furthermore, the expressed APLNR appeared to be inadequately localized to the cell membrane, as indicated by its higher concentration in the immunoprecipitate compared to the cell membrane fraction (Fig. 6). As proteins mature via pathways mediated by the endoplasmic reticulum and Golgi apparatus, insufficient PTMs may result in their accumulation in intracellular organelles. In the mouse spinal cord, APLNR was detected at a molecular weight 10 kDa higher than expected and was thermally stable, suggesting its conjugation with other proteins to gain non-denaturing properties (Fig. 7).

    No APLNR isoforms have been reported in the GPCR database (GPCRdb), where GPCR isoforms can be searched (https://gpcrdb.org/protein/apj_human/). The additional 10 kDa molecular weight of APLNR could potentially result from a PTM involving covalent bonds; however, we did not identify the exact modifiers responsible. These modifications are important for improving thermostability and may also affect the functional expression of APLNR. Since ubiquitin is an 8.6 kDa protein, it is a likely candidate for causing the observed molecular weight shift of APLNR. However, since ubiquitinated proteins undergo further polyubiquitination and are degraded by the proteasome system, it is unlikely that polyubiquitination is the PTM of APLNR. SUMO and NEDD are modified at lysine residues by an E1–E3-mediated mechanism similar to ubiquitin; however, unlike ubiquitin, polymerization does not occur [26]. Unfortunately, our prediction did hold true, and neither SUMO nor NEDD was found to be associated with APLNR. Other proteins such as ATG8, ATG12, URM1, UFM1, FAT10 and ISG15 [27] have approximate molecular weights of 10 kDa and undergo isopeptidation. Although these modifications are uncommon, a comprehensive analysis of PTMs for APLNR is required. The glycosylphosphatidylinositol (GPI) anchor is a modification of approximately 10 kDa and is primarily associated with anchoring glycolipids to soluble proteins to lipid bilayer membranes; however, transmembrane proteins do not typically require the GPI anchor [28]. While acetylcholinesterase, alkaline phosphatase and CD90 are known to be modified with GPI anchors, the prediction by NEtGPI (https://services.healthtech.dtu.dk/services/NetGPI-1.1/) suggests that APLNR is not GPI-modified. Nevertheless, of the aforementioned three genes, only acetylcholinesterase is predicted to be GPI-modified, indicating the possibility of GPI modification on APLNR (data not shown). Although the specific PTM that confers stability to APLNR remains unknown, it is evident that APLNR was expressed in an unstable manner under the experimental conditions employed in this study.

    In a previous study, we found that APLNR was expressed in motor neuron cell bodies in the mouse spinal cord [13]. In this study, immunoblotting revealed a substantial number of APLNR-positive bands in the sacral fraction, ranging from the cervical to sacral regions of the mouse spinal cord (Fig. 7). Compared with other spinal cord fractions, the sacral fraction contained peripheral nerve axons (cauda equina). Similarly, the lumbar spinal cord fraction exhibited strong APLNR-positive band intensity (data not shown). Schwann cells, responsible for myelinating peripheral nerve axons, express Gpr126, Gpr44 and LPA1, which are involved in their maturation, myelination and interaction with neuronal axons [29]. In particular, Gpr126 may promote the differentiation of Schwann cells by activating adenylate cyclase and increasing cAMP levels. In addition, Gpr126 is coupled with Gs and Gi proteins, which decreases cAMP concentration and activates Ras and PI3K cascades [29]. The expression of APLNR in Schwann cells suppresses their differentiation by decreasing cAMP levels. Gpr126 is also expressed in motor neuron axons, suggesting that Gpr126 signaling regulates cAMP in both Schwann cells and motor neurons. Therefore, APLNR may be expressed in Schwann cells and/or motor neurons and potentially at the nodes of Ranvier. The GABAB receptor [30] and protease-activated receptor 1 (PAR-1) [31] are present on the nodes of Ranvier and are involved in myelination and conduction blockage during brain hemorrhage.

    Apelin receptor exhibits diverse functions in various cell types, including those in the cardiovascular and nervous systems. Therefore, APLNR in HEK cells used in this study may not have a specific mechanism of action. We believe that activation of ERK may be due to activation of the autophagy pathway to degrade APLNR that cannot reach the cell membrane and accumulates within the cell. In fact, according to Nuria Martinez-Lopez et.al., activation of ERK by autophagosome formation has been reported [32]. (1) Since this is a reaction caused by the agonist expressed in the membrane, we believe that activation via β-arrestin does not contribute as the result that it did not reach the cell membrane this time, but it will reach the membrane in the future. We believe that it is necessary to proceed with analysis of the relationship with PTM. (2) The reagent used for this data was only (Pyr1)-apelin-13, but similar studies were conducted using ML221 as an APLNR agonist, and the results were the same, so the difference between isoforms is I concluded that it would not be possible, but I am thinking of considering isoforms next time. Other studies investigating APLNR in HEK cells have demonstrated the activation of signaling pathways [14]. Transfection experiments using vectors have shown a wide range of effective apelin concentrations. The expression pattern in HEK cells may change in vitro depending on the experimental conditions, such as the plasmid vector used. Therefore, it would be worthwhile to consider transfecting COS7 or CHO-K1 cells, which are commonly used for protein expression after gene transfer, as well as Neuro2A or NSC-34 cells as potential neuronal cell line candidates.

    In vivo experiments have reported variations in the molecular weight of APLNR among different tissues. For example, 43 kDa in the mouse heart [33] and 50 kDa was observed in the spinal cord in this study, suggesting that PTMs differ between central and peripheral tissues. To elucidate the details of APLNR signal transduction, it is essential to establish an assay system that stably expresses APLNR on the cell membrane. To achieve this, we plan to reassess various experimental conditions and continue analyzing the physiological and pathophysiological functions mediated by APLNR.



    [1] V. V. Grechko, Molecular DNA markers in phylogeny and systematics, Russ. J. Genet., 38 (2002), 851–868. https://doi.org/10.1023/A:1016890509443 doi: 10.1023/A:1016890509443
    [2] A. Patwardhan, S. Ray, A. Roy, Molecular markers in phylogenetic studies—A review, J. Phylogen. Evolution Biol., 2 (2014), 1–9. https://doi.org/10.4172/2329-9002.1000131 doi: 10.4172/2329-9002.1000131
    [3] L. Wei, Selection On synonymous Mutations Revealed by 1135 Genomes of Arabidopsis thaliana, Evol. Bioinform. Online, 16 (2020) 1176934320916794. https://doi.org/10.1177/1176934320916794 doi: 10.1177/1176934320916794
    [4] W. M. Brown, M. George, A. C. Wilson, Rapid evolution of animal mitochondrial DNA, PNAS, 76 (1979), 1967–1971. https://doi.org/10.1073/pnas.76.4.1967 doi: 10.1073/pnas.76.4.1967
    [5] A. A. Bannikova, Molecular markers and modern phylogenetics of mammals, Zh. Obshch. Biol., 65 (2004), 278–305. (in Russian)
    [6] A. C. Wilson, R. L. Cann, S. M. Carr, M. George, U. B. Gyllensten, K. M. Helm-Bychowski, et al., Mitochondrial DNA and two perspectives on evolutionary genetics, Biol. J. Linnean Soc., 26 (1985), 375–400. doi.org/10.1111/j.1095-8312.1985.tb02048.x
    [7] G. Pesole, E. Sbisa, G. Preparata, C. Saccone, The evolution of the mitochondrial D-loop region and the origin of modern man, Mol. Biol. Evol., 9 (1992), 587–598. https://doi.org/10.1093/oxfordjournals.molbev.a040747 doi: 10.1093/oxfordjournals.molbev.a040747
    [8] W. M. Brown, E. H. Prager, A. Wang, A. C. Wilson, Mitochondrial DNA sequences of primates: Tempo and mode of evolution, J. Mol. Evol., 18 (1982), 225–239.
    [9] H. Cernohorska, S. Kubickova, O. Kopecna, A. I. Kulemzina, P. L. Perelman, F. F. Elder, et al., Molecular cytogenetic insights to the phylogenetic affinities of the giraffe (Giraffa camelopardalis) and pronghorn (Antilocapra americana), Chromosome Res., 2 (2013), 447–460. https://doi.org/10.1007/s10577-013-9361-0 doi: 10.1007/s10577-013-9361-0
    [10] D. Huchon, O. Madsen, M. Sibbald, K. Ament, M. J. Stanhope, F. Catzeflis, et al., Rodent phylogeny and a timescale for the evolution of glires: Evidence from an extensive taxon sampling using three nuclear genes, Mol. Biol. Evol., 19 (2002), 1053–1065. https://doi.org/10.1093/oxfordjournals.molbev.a004164 doi: 10.1093/oxfordjournals.molbev.a004164
    [11] M. Weksler, Phylogeny of Neotropical oryzomyine rodents (Muridae: Sigmodontinae) based on the nuclear IRBP exon, Mol. Phylogenet. Evol., 29 (2003), 331–349. doi.org/10.1016/S1055-7903(03)00132-5
    [12] E. Zietkiewicz, A. Rafalski, D. Labuda, Genome fingerprinting by simple sequence repeat (SSR)-anchored polymerase chain reaction amplification, Genomics, 20 (1994), 176–183. https://doi.org/10.1006/geno.1994.1151 doi: 10.1006/geno.1994.1151
    [13] H. Nybom, Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants, Mol Ecol., 13 (2004), 1143–1155. https://doi.org/10.1111/j.1365-294X.2004.02141.x doi: 10.1111/j.1365-294X.2004.02141.x
    [14] S. R. Pandi-Perumal, V. Srinivasan, G. J. Maestroni, D. P. Cardinali, B. Poeggeler, R. Hardeland, Melatonin: Nature's most versatile biological signal? FEBS, 273 (2006), 2813–2838. https://doi.org/10.1111/j.1742-4658.2006.05322.x doi: 10.1111/j.1742-4658.2006.05322.x
    [15] A. F. Alamdari, S. Rahnemayan, H. Rajabi, N. Vahed, H. R. K. Kashani, A. Rezabakhsh, et al., Melatonin as a promising modulator of aging related neurodegenerative disorders: Role of microRNAs, Pharmacol. Res., 173 (2021), 105839. https://doi.org/10.1016/j.phrs.2021.105839 doi: 10.1016/j.phrs.2021.105839
    [16] S. Arias-Santiago, J. Aneiros-Fernández, B. Arias-Santiago, M. S. Girón-Prieto, M. Caba-Molina, A. López-Valverde, et al., MTNR1A receptor expression in normal and pathological human salivary glands, Anticancer Res., 32 (2012), 4765–4771.
    [17] Y. Y. Li, H. Wang, Y. Y. Zhang, Melatonin receptor 1B gene rs10830963 C/G polymorphism associated with type 2 diabetes mellitus: An updated meta-analysis of 13,752 participants, Heliyon, 8 (2022), e11786. https://doi.org/10.1016/j.heliyon.2022.e11786 doi: 10.1016/j.heliyon.2022.e11786
    [18] N. R. Sundaresan, M. D. Marcus Leo, J. Subramani, D. Anish, M. Sudhagar, K. A. Ahmed, et al., Expression analysis of melatonin receptor subtypes in the ovary of domestic chicken, Vet. Res. Commun., 33 (2009), 49–56. https://doi.org/10.1007/s11259-008-9071-9 doi: 10.1007/s11259-008-9071-9
    [19] C. Jones, G. Helfer, R. Brandstätter, Melatonin receptor expression in the zebra finch brain and peripheral tissues, Chronobiol. Int., 29 (2012), 189–202. https://doi.org/10.3109/07420528.2011.642912 doi: 10.3109/07420528.2011.642912
    [20] D. Y. Li, D. G. Smith, R. Hardeland, M. Y. Yang, H. L. Xu, L. Zhang, et al., Melatonin receptor genes in vertebrates, Int. J. Mol. Sci., 14 (2013), 11208–11223. https://doi.org/10.3390/ijms140611208 doi: 10.3390/ijms140611208
    [21] C. von Gall, J. H. Stehle, D. R. Weaver, Mammalian melatonin receptors: molecular biology and signal transduction, Cell Tissue Res., 309 (2002), 151–162. https://doi.org/10.1007/s00441-002-0581-4 doi: 10.1007/s00441-002-0581-4
    [22] M. Migaud, A. Daveau, B. Malpaux, MTNR1A melatonin receptors in the ovine premammillary hypothalamus: day-night variation in the expression of the transcripts, Biol. Reprod., 72 (2005), 393–398. https://doi.org/10.1095/biolreprod.104.030064 doi: 10.1095/biolreprod.104.030064
    [23] J. Drew, P. Barrett, J. Mercer, K. Moar, E. Canet, P. Delagrange, et al., Localization of the melatonin-related receptor in the rodent brain and peripheral tissues, J. Neuroendocrinol., 13 (2001), 453–458. https://doi.org/10.1046/j.1365-2826.2001.00651.x doi: 10.1046/j.1365-2826.2001.00651.x
    [24] L. Pinato, D. Ramos, A. Hataka, P. S. Rossignoli, M. D. J. Granado, M. C. Mazzetto, et al., Day/night expression of MT1 and MT2 receptors in hypothalamic nuclei of the primate Sapajus paella, J. Chem. Neuroanat., 81 (2017), 10–17. https://doi.org/10.1016/j.jchemneu.2017.01.005 doi: 10.1016/j.jchemneu.2017.01.005
    [25] H. Ma, J. Kang, W. Fan, H. He, F. Huang, ROR: Nuclear receptor for melatonin or not? Molecules, 26 (2021), 2693. https://doi.org/10.3390/molecules26092693 doi: 10.3390/molecules26092693
    [26] S. D. Huo, R. J. Long, Melatonin receptor (MTNR1A and MTNR2B) expression during the breeding season in the yak (Bos grunniens), Czech J. Anim. Sci., 59 (2014), 140–145. https://doi.org/10.17221/7294-CJAS doi: 10.17221/7294-CJAS
    [27] J. Baker, K. Kimpinski, Role of melatonin in blood pressure regulation: An adjunct anti-hypertensive agent, CEPP, 45 (2018), 755–766. https://doi.org/10.1111/1440-1681.12942 doi: 10.1111/1440-1681.12942
    [28] M. A. Quera-Salva, U. Kilic-Huck, M. F. Vecchierini, Members of the MEL consensus group of the SFRMS. Melatonin (MEL) and its use in circadian rhythm sleep-wake disorders: Recommendations of the French Medical and Research Sleep Society (SFRMS), Rev. Neurol. (Paris)., 177 (2021), 235–244. https://doi.org/10.1016/j.neurol.2020.07.021 doi: 10.1016/j.neurol.2020.07.021
    [29] L. A. Ostrin, Ocular and systemic melatonin and the influence of light exposure, Clin. Exp. Optom., 102 (2019), 99–108. https://doi.org/10.1111/cxo.12824 doi: 10.1111/cxo.12824
    [30] Y. Tao, B. Hu, Z. Ma, H. Li, E. Du, G. Wang, et al., Intravitreous delivery of melatonin affects the retinal neuron survival and visual signal transmission: in vivo and ex vivo study, Drug Deliv., 27 (2020), 1386–1396. https://doi.org/10.1080/10717544.2020.1818882 doi: 10.1080/10717544.2020.1818882
    [31] S. Gurunathan, M. Qasim, M. H. Kang, J. H. Kim, Role and therapeutic potential of melatonin in various type of cancers, Onco. Targets Ther., 14 (2021), 2019–2052. https://doi.org/10.2147/OTT.S298512 doi: 10.2147/OTT.S298512
    [32] Y. J. Guh, T. K. Tamai, T. Yoshimura, The underlying mechanisms of vertebrate seasonal reproduction, Proc. Jpn. Acad. Ser. B Phys. Biol. Sci., 95 (2019), 343–357. https://doi.org/10.2183/pjab.95.025 doi: 10.2183/pjab.95.025
    [33] M. V. Danilova, E. N. Usoltseva, Significance of the pineal gland hormone melatonin in maintaining the health of women of reproductive age (a review), Obstetr. Gynecol. Reproduct., 13 (2019), 337–344. https://doi.org/10.17749/2313-7347.2019.13.4.337-344 doi: 10.17749/2313-7347.2019.13.4.337-344
    [34] C. C. Maganhin, L. F. Fuchs, R. S. Simões, R. M. Oliveira-Filho, M. de Jesus Simões, E. C. Baracat, et al., Effects of melatonin on ovarian follicles, Eur. J. Obstet.Gynecol. Reprod. Biol., 166 (2013), 178–184. https://doi.org/10.1016/j.ejogrb.2012.10.006 doi: 10.1016/j.ejogrb.2012.10.006
    [35] J. S. Heo, S. Pyo, J. Y. Lim, D. W. Yoon, B. Y. Kim, J. H. Kim, et al., Biological effects of melatonin on human adipose‑derived mesenchymal stem cells, Int. J. Mol. Med., 44 (2019), 2234–2244. https://doi.org/10.3892/ijmm.2019.4356 doi: 10.3892/ijmm.2019.4356
    [36] A. Sotthibundhu, P. Phansuwan-Pujito, P. Govitrapong, Melatonin increases proliferation of cultured neural stem cells obtained from adult mouse subventricular zone, J. Pineal. Res., 49 (2010), 291–300. https://doi.org/10.1111/j.1600-079X.2010.00794.x doi: 10.1111/j.1600-079X.2010.00794.x
    [37] G. N. Georgiev, E. Marinova, R. Konakchieva, P. Todorov, Melatonin selectively influences the transcription of pluripotency and differentiation markers in human non-cancer cells, Biotechnol. Biotechnol. Equipment, 33 (2019), 286–293. https://doi.org/10.1080/13102818.2019.1571440 doi: 10.1080/13102818.2019.1571440
    [38] J. Fu, S. D. Zhao, H. J. Liu, Q. H. Yuan, S. M. Liu, Y. M. Zhang, et al., Melatonin promotes proliferation and differentiation of neural stem cells subjected to hypoxia in vitro, J. Pineal. Res., 51 (2011), 104–112. https://doi.org/10.1111/j.1600-079X.2011.00867.x doi: 10.1111/j.1600-079X.2011.00867.x
    [39] C. Tchio, S. K. Musani, A. Quarshie, G. Tosini, Association between MTNR1B polymorphisms and obesity in African American: Findings from the Jackson Heart Study, BMC Med. Genom., 14 (2021), 136. https://doi.org/10.1186/s12920-021-00983-2 doi: 10.1186/s12920-021-00983-2
    [40] R. M. Slominski, R. J. Reiter, N. Schlabritz–Loutsevitch, R. S. Ostrom, A. T. Slominski, Melatonin membrane receptors in peripheral tissues: Distribution and functions, Mol. Cell Endocrinol., 351 (2012), 152–166.
    [41] Q. Xia, Z. X. Chen, Y. C. Wang, Y. S. Ma, F. Zhang, W. Che, et al., Association between the melatonin receptor 1B gene polymorphism on the risk of type 2 diabetes, impaired glucose regulation: A meta-analysis, PLoS One, 7 (2012), e50107. https://doi.org/10.1371/journal.pone.0050107 doi: 10.1371/journal.pone.0050107
    [42] D. Buonfiglio, C. Tchio, I. Furigo, J. Jr. Donato, K. Baba, J. Cipolla-Neto, et al., Removing melatonin receptor type 1 signaling leads to selective leptin resistance in the arcuate nucleus, J. Pineal. Res., 67 (2019), e12580. https://doi.org/10.1111/jpi.12580 doi: 10.1111/jpi.12580
    [43] A. Karamitri, R. Jockers, Melatonin in type 2 diabetes mellitus and obesity, Nat. Rev.Endocrinol., 1 (2019), 105–125. https://doi.org/10.1038/s41574-018-0130-1 doi: 10.1038/s41574-018-0130-1
    [44] H. S. Dashti, C. Vetter, J. M. Lane, M. C. Smith, A. R. Wood, M. N. Weedon, et al., Assessment of MTNR1B Type 2 diabetes genetic risk modification by shift work and morningness-eveningness preference in the UK Biobank, Diabetes, 69 (2020), 259–266. https://doi.org/10.2337/db19-0606 doi: 10.2337/db19-0606
    [45] P. Chaste, N. Clement, O. Mercati, J. L. Guillaume, R. Delorme, H. G. Botros, et al., Identification of pathway-biased and deleterious melatonin receptor mutants in autism spectrum disorders and in the general population, PLoS One, 5 (2010), e11495. https://doi.org/10.1371/journal.pone.0011495 doi: 10.1371/journal.pone.0011495
    [46] D. M. Irwin, T. D. Kocher, A. C. Wilson, Evolution of the cytochrome b gene of mammals, J. Mol. Evol., 32 (1991), 128–144. https://doi.org/10.1007/BF02515385 doi: 10.1007/BF02515385
    [47] S. S. Tobe, A. C. Kitchener, A. M. Linacre, Reconstructing mammalian phylogenies: A detailed comparison of the cytochrome b and cytochrome oxidase subunit I mitochondrial genes, PLoS One, 5 (2010), e14156. https://doi.org/10.1371/journal.pone.0014156 doi: 10.1371/journal.pone.0014156
    [48] J. D. Thompson, D. G. Higgins, T. J. Gibson, CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice, Nucleic Acids Res., 22 (1994), 4673–4680. https://doi.org/10.1093/nar/22.22.4673 doi: 10.1093/nar/22.22.4673
    [49] M. A. Larkin, G. Blackshields, N. P. Brown, R. Chenna, P. A. McGettigan, H. McWilliam, et al., Clustal W and Clustal X version 2.0., Bioinformatics, 23 (2007), 2947–2948. https://doi.org/10.1093/bioinformatics/btm404 doi: 10.1093/bioinformatics/btm404
    [50] K. Tamura, G. Stecher, D. Peterson, A. Filipski, S. Kumar, MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0., Mol. Biol. Evol., 30 (2013), 2725–2729. https://doi.org/10.1093/molbev/mst197 doi: 10.1093/molbev/mst197
    [51] F. Tajima, N. Takezaki, Estimation of evolutionary distance for reconstructing molecular phylogenetic trees, Mol. Biol. Evol., 11 (1994), 278–286. https://doi.org/10.1093/oxfordjournals.molbev.a040109 doi: 10.1093/oxfordjournals.molbev.a040109
    [52] M. Hasegawa, H. Kishino, T. Yano, Dating of the human-ape splitting by a molecular clock of mitochondrial DNA, J. Mol. Evol., 22 (1985), 160–174. https://doi.org/10.1007/BF02101694 doi: 10.1007/BF02101694
    [53] D. Posada, K. A. Crandall, Selecting the best-fit model of nucleotide substitution, Syst. Biol., 50 (2001), 580–601.
    [54] J. Felsenstein, Confidence limits on phylogenies: An approach using the bootstrap, Evolution, 39 (1985), 783–791. https://doi.org/10.1111/j.1558-5646.1985.tb00420.x doi: 10.1111/j.1558-5646.1985.tb00420.x
    [55] W. J. Murphy, P. A. Pevzner, S. J. O'Brien, Mammalian phylogenomics comes of age, Trends Genet., 20 (2004), 631–639. https://doi.org/10.1016/j.tig.2004.09.005 doi: 10.1016/j.tig.2004.09.005
    [56] M. Buckley, Ancient collagen reveals evolutionary history of the endemic South American «ungulates», Proc. R. Soc. B Biol. Sci., 282 (2015), 20142671. https://doi.org/10.1098/rspb.2014.2671 doi: 10.1098/rspb.2014.2671
    [57] A. Sánchez-Gracia, J. Rozas, Divergent evolution and molecular adaptation in the Drosophila odorant-binding protein family: inferences from sequence variation at the OS-E and OS-F genes, BMC Evol. Biol., 8 (2008), 323. https://doi.org/10.1186/1471-2148-8-323 doi: 10.1186/1471-2148-8-323
    [58] C. A. Driscoll, D. W. Macdonald, S. J. O'Brien, From wild animals to domestic pets, an evolutionary view of domestication, PNAS, 106 (2009), 9971–9978. https://doi.org/10.1073/pnas.0901586106 doi: 10.1073/pnas.0901586106
    [59] C. J. Hoskin, M. Higgie, K. R. McDonald, C. Moritz, Reinforcement drives rapid allopatric speciation, Nature, 437 (2005), 1353–1356. https://doi.org/10.1038/nature04004 doi: 10.1038/nature04004
    [60] B. M. Fitzpatrick, J. A. Fordyce, S. Gavrilets, What, if anything, is sympatric speciation? J. Evol. Biol., 21 (2008), 1452–1459. https://doi.org/10.1111/j.1420-9101.2008.01611.x doi: 10.1111/j.1420-9101.2008.01611.x
    [61] W. F. Bottke, D. Vokrouhlický, D. Nesvorný, An asteroid breakup 160 Myr ago as the probable source of the K/T impactor, Nature, 449 (2007), 48—53. https://doi.org/10.1038/nature06070 doi: 10.1038/nature06070
    [62] G. Keller, T. Adatte, S. Gardin, A. Bartolini, S. Bajpai, Main Deccan volcanism phase ends near the K–T boundary: Evidence from the Krishna–Godavari Basin, SE India, Earth Planet. Sci. Lett., 268 (2008), 293–311. doi.org/10.1016/j.epsl.2008.01.015
    [63] R. Nielsen, Statistical tests of selective neutrality in the age of genomics, Heredity, 86 (2001), 641–647. doi.org/10.1046/j.1365-2540.2001.00895.x
    [64] F. Tajima, Statistical method for testing the neutral mutation hypothesis by DNA polymorphism, Genetics, 123 (1989), 585–595.
    [65] K. Tamura, D. Peterson, N. Peterson, G. Stecher, M. Nei, S. Kumar, MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods, Mol. Biol. Evol., 28 (2011), 2731–2739. doi.org/10.1093/molbev/msr121
    [66] M. Nei, S. Kumar, Molecular Evolution and Phylogenetics, NY.: Oxford University Press, (2000), 333.
    [67] R. S. Holmes, L. A. Cox, Comparative structures and evolution of vertebrate lipase H (LIPH) genes and proteins: A relative of the phospholipase A1 gene families, 3 Biotech, 2 (2012), 263–275. doi.org/10.1007/s13205-012-0087-z
    [68] C. Lavialle, G. Cornelis, A. Dupressoir, C. Esnault, O. Heidmann, C. Vernochet, et al., Paleovirology of 'syncytins', retroviral env genes exapted for a role in placentation, Philos. Trans. R. Soc. Lond. B. Biol. Sci., 368 (2013), 20120507. https://doi.org/10.1098/rstb.2012.0507 doi: 10.1098/rstb.2012.0507
    [69] A. Cooper, R. Fortey, Evolutionary explosions and the phylogenetic fuse, Trends Ecol. Evol., 13 (1998), 151–156. https://doi.org/10.1016/s0169-5347(97)01277-9 doi: 10.1016/s0169-5347(97)01277-9
    [70] M. Spaulding, M. A. O'Leary, J. Gatesy, Relationships of Cetacea (Artiodactyla) among mammals: Increased taxon sampling alters interpretations of key fossils and character evolution, PLoS One, 4 (2009), e7062. https://doi.org/10.1371/journal.pone.0007062 doi: 10.1371/journal.pone.0007062
    [71] Z. Luo, In search of the whales' sisters, Nature, 404 (2000), 235–237. https://doi.org/10.1038/35005194 doi: 10.1038/35005194
    [72] R. M. D. Beck, C. Baillie, Improvements in the fossil record may largely resolve current conflicts between morphological and molecular estimates of mammal phylogeny, Proc. R. Soc. B: Biol. Sci., 285 (2018), 20181632. https://doi.org/10.1098/rspb.2018.1632 doi: 10.1098/rspb.2018.1632
  • This article has been cited by:

    1. Gerry T.M. Wagenaar, Gert N. Moll, Advances in the therapeutic potentials of ligands of the apelin receptor APJ, 2025, 00142999, 177302, 10.1016/j.ejphar.2025.177302
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2803) PDF downloads(87) Cited by(0)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog