High‐resolution analysis of the conformational transition of pro‐apoptotic Bak at the lipid membrane

Permeabilization of the outer mitochondrial membrane by pore‐forming Bcl2 proteins is a crucial step for the induction of apoptosis. Despite a large set of data suggesting global conformational changes within pro‐apoptotic Bak during pore formation, high‐resolution structural details in a membrane environment remain sparse. Here, we used NMR and HDX‐MS (Hydrogen deuterium exchange mass spectrometry) in lipid nanodiscs to gain important high‐resolution structural insights into the conformational changes of Bak at the membrane that are dependent on a direct activation by BH3‐only proteins. Furthermore, we determined the first high‐resolution structure of the Bak transmembrane helix. Upon activation, α‐helix 1 in the soluble domain of Bak dissociates from the protein and adopts an unfolded and dynamic potentially membrane‐bound state. In line with this finding, comparative protein folding experiments with Bak and anti‐apoptotic BclxL suggest that α‐helix 1 in Bak is a metastable structural element contributing to its pro‐apoptotic features. Consequently, mutagenesis experiments aimed at stabilizing α‐helix 1 yielded Bak variants with delayed pore‐forming activity. These insights will contribute to a better mechanistic understanding of Bak‐mediated membrane permeabilization.     

The effects of intramolecular and intermolecular electrostatic repulsions on the stability and aggregation of NISTmAb revealed by HDX‐MS,DSC, and nanoDSF

The stability and aggregation of NIST monoclonal antibody (NISTmAb) were investigated by hydrogen/deuterium exchange mass spectrometry (HDX‐MS), differential scanning calorimetry (DSC), and nano‐differential scanning fluorimetry (nanoDSF). NISTmAb was prepared in eight formulations at four different pHs (pH 5, 6, 7, and 8) in the presence and absence of 150 mM NaCl and analyzed by the three methods. The HDX‐MS results showed that NISTmAb is more conformationally stable at a pH near its isoelectric point (pI) in the presence of NaCl than a pH far from its pI in the absence of NaCl. The stabilization effects were global and not localized. The midpoint temperature of protein thermal unfolding transition results also showed the C_H2 domain of the protein is more conformationally stable at a pH near its pI. On the other hand, the onset of aggregation temperature results showed that NISTmAb is less prone to aggregate at a pH far from its pI, particularly in the absence of NaCl. These seemingly contradicting results, higher conformational stability yet higher aggregation propensity near the pI than far away from the pI, can be explained by intramolecular and intermolecular electrostatic repulsion using Lumry‐Eyring model, which separates folding/unfolding equilibrium and aggregation event. The further a pH from the pI, the higher the net charge of the protein. The higher net charge leads to greater intramolecular and intermolecular electrostatic repulsions. The greater intramolecular electrostatic repulsion destabilizes the protein and the greater intermolecular electrostatic repulsion prevents aggregation of the protein molecules at pH far from the pI.     

On‐site genetic analysis for species identification using lab‐on‐a‐chip

This paper presents a microfluidic device capable of performing genetic analysis on dung samples to identify White Rhinoceros (Ceratotherium simum). The development of a microfluidic device, which can be used in the field, offers a portable and cost‐effective solution for DNA analysis and species identification to aid conservation efforts. Optimization of the DNA extraction processes produced equivalent yields compared to conventional kit‐based methods within just 5 minutes. The use of a color‐changing loop‐mediated isothermal amplification reaction for simultaneous detection of the cytochrome B sequence of C. simum enabled positive results to be obtained within as little as 30 minutes. Field testing was performed at Knowsley Safari to demonstrate real‐world applicability of the microfluidic device for testing of biological samples.     

Genome‐scale metabolic modeling reveals SARS‐CoV‐2‐induced metabolic changes and antiviral targets

Tremendous progress has been made to control the COVID‐19 pandemic caused by the SARS‐CoV‐2 virus. However, effective therapeutic options are still rare. Drug repurposing and combination represent practical strategies to address this urgent unmet medical need. Viruses, including coronaviruses, are known to hijack host metabolism to facilitate viral proliferation, making targeting host metabolism a promising antiviral approach. Here, we describe an integrated analysis of 12 published in vitro and human patient gene expression datasets on SARS‐CoV‐2 infection using genome‐scale metabolic modeling (GEM), revealing complicated host metabolism reprogramming during SARS‐CoV‐2 infection. We next applied the GEM‐based metabolic transformation algorithm to predict anti‐SARS‐CoV‐2 targets that counteract the virus‐induced metabolic changes. We successfully validated these targets using published drug and genetic screen data and by performing an siRNA assay in Caco‐2 cells. Further generating and analyzing RNA‐sequencing data of remdesivir‐treated Vero E6 cell samples, we predicted metabolic targets acting in combination with remdesivir, an approved anti‐SARS‐CoV‐2 drug. Our study provides clinical data‐supported candidate anti‐SARS‐CoV‐2 targets for future evaluation, demonstrating host metabolism targeting as a promising antiviral strategy.     

Thermodynamic coupling between neighboring binding sites in homo‐oligomeric ligand sensing proteins from mass resolved ligand‐dependent population distributions

Homo‐oligomeric ligand‐activated proteins are ubiquitous in biology. The functions of such molecules are commonly regulated by allosteric coupling between ligand‐binding sites. Understanding the basis for this regulation requires both quantifying the free energy ΔG transduced between sites, and the structural basis by which it is transduced. We consider allostery in three variants of the model ring‐shaped homo‐oligomeric trp RNA‐binding attenuation protein (TRAP). First, we developed a nearest‐neighbor statistical thermodynamic binding model comprising microscopic free energies for ligand binding to isolated sites ΔG ₀, and for coupling between adjacent sites, ΔG _α . Using the resulting partition function (PF) we explored the effects of these parameters on simulated population distributions for the 2 ^N possible liganded states. We then experimentally monitored ligand‐dependent population shifts using conventional spectroscopic and calorimetric methods and using native mass spectrometry (MS). By resolving species with differing numbers of bound ligands by their mass, native MS revealed striking differences in their ligand‐dependent population shifts. Fitting the populations to a binding polynomial derived from the PF yielded coupling free energy terms corresponding to orders of magnitude differences in cooperativity. Uniquely, this approach predicts which of the possible 2 ^N liganded states are populated at different ligand concentrations, providing necessary insights into regulation. The combination of statistical thermodynamic modeling with native MS may provide the thermodynamic foundation for a meaningful understanding of the structure–thermodynamic linkage that drives cooperativity.     

Cryo‐EM structures of pentameric autoinducer‐2 exporter from Escherichia coli reveal its transport mechanism

Bacteria utilize small extracellular molecules to communicate in order to collectively coordinate their behaviors in response to the population density. Autoinducer‐2 (AI‐2), a universal molecule for both intra‐ and inter‐species communication, is involved in the regulation of biofilm formation, virulence, motility, chemotaxis, and antibiotic resistance. While many studies have been devoted to understanding the biosynthesis and sensing of AI‐2, very little information is available on its export. The protein TqsA from Escherichia coli, which belongs to the AI‐2 exporter superfamily, has been shown to export AI‐2. Here, we report the cryogenic electron microscopic structures of two AI‐2 exporters (TqsA and YdiK) from E. coli at 3.35 Å and 2.80 Å resolutions, respectively. Our structures suggest that the AI‐2 exporter exists as a homo‐pentameric complex. In silico molecular docking and native mass spectrometry experiments were employed to demonstrate the interaction between AI‐2 and TqsA, and the results highlight the functional importance of two helical hairpins in substrate binding. We propose that each monomer works as an independent functional unit utilizing an elevator‐type transport mechanism.     

Selective cross‐linking of coinciding protein assemblies by in‐gel cross‐linking mass spectrometry

Cross‐linking mass spectrometry has developed into an important method to study protein structures and interactions. The in‐solution cross‐linking workflows involve time and sample consuming steps and do not provide sensible solutions for differentiating cross‐links obtained from co‐occurring protein oligomers, complexes, or conformers. Here we developed a cross‐linking workflow combining blue native PAGE with in‐gel cross‐linking mass spectrometry (IGX‐MS). This workflow circumvents steps, such as buffer exchange and cross‐linker concentration optimization. Additionally, IGX‐MS enables the parallel analysis of co‐occurring protein complexes using only small amounts of sample. Another benefit of IGX‐MS, demonstrated by experiments on GroEL and purified bovine heart mitochondria, is the substantial reduction of undesired over‐length cross‐links compared to in‐solution cross‐linking. We next used IGX‐MS to investigate the complement components C5, C6, and their hetero‐dimeric C5b6 complex. The obtained cross‐links were used to generate a refined structural model of the complement component C6, resembling C6 in its inactivated state. This finding shows that IGX‐MS can provide new insights into the initial stages of the terminal complement pathway.     

Biochemically validated structural model of the 15‐subunit intraflagellar transport complex IFT‐B

Cilia are ubiquitous eukaryotic organelles impotant for cellular motility, signaling, and sensory reception. Cilium formation requires intraflagellar transport of structural and signaling components and involves 22 different proteins organized into intraflagellar transport (IFT) complexes IFT‐A and IFT‐B that are transported by molecular motors. The IFT‐B complex constitutes the backbone of polymeric IFT trains carrying cargo between the cilium and the cell body. Currently, high‐resolution structures are only available for smaller IFT‐B subcomplexes leaving > 50% structurally uncharacterized. Here, we used Alphafold to structurally model the 15‐subunit IFT‐B complex. The model was validated using cross‐linking/mass‐spectrometry data on reconstituted IFT‐B complexes, X‐ray scattering in solution, diffraction from crystals as well as site‐directed mutagenesis and protein‐binding assays. The IFT‐B structure reveals an elongated and highly flexible complex consistent with cryo‐electron tomographic reconstructions of IFT trains. The IFT‐B complex organizes into IFT‐B1 and IFT‐B2 parts with binding sites for ciliary cargo and the inactive IFT dynein motor, respectively. Interestingly, our results are consistent with two different binding sites for IFT81/74 on IFT88/70/52/46 suggesting the possibility of different structural architectures for the IFT‐B1 complex. Our data present a structural framework to understand IFT‐B complex assembly, function, and ciliopathy variants.     

MicroRNA‐495‐3p diminishes doxorubicin‐induced cardiotoxicity through activating AKT

Doxorubicin (Dox) is a broad‐spectrum antitumour agent; however, its clinical application is impeded due to the cumulative cardiotoxicity. The present study aims to investigate the role and underlying mechanisms of microRNA‐495‐3p (miR‐495‐3p) in Dox‐induced cardiotoxicity. Herein, we found that cardiac miR‐495‐3p expression was significantly decreased in Dox‐treated hearts, and that the miR‐495‐3p agomir could prevent oxidative stress, cell apoptosis, cardiac mass loss, fibrosis and cardiac dysfunction upon Dox stimulation. In contrast, the miR‐495‐3p antagomir dramatically aggravated Dox‐induced cardiotoxicity in mice. Besides, we found that the miR‐495‐3p agomir attenuated, while the miR‐495‐3p antagomir exacerbated Dox‐induced oxidative stress and cellular injury in vitro. Mechanistically, we demonstrated that miR‐495‐3p directly bound to the 3′‐untranslational region of phosphate and tension homology deleted on chromosome ten (PTEN), downregulated PTEN expression and subsequently activated protein kinase B (PKB/AKT) pathway, and that PTEN overexpression or AKT inhibition completely abolished the cardioprotective effects of the miR‐495‐3p agomir. Our study for the first time identify miR‐495‐3p as an endogenous protectant against Dox‐induced cardiotoxicity through activating AKT pathway in vivo and in vitro.     

Tacrolimus alleviates LPS‐induced AKI by inhibiting TLR4/MyD88/NF‐κB signalling in mice

Lipopolysaccharide (LPS)‐induced sepsis‐associated acute kidney injury (SA‐AKI) is a model of clinical serious care syndrome, with high morbidity and mortality. Tacrolimus (TAC), a novel immunosuppressant that inhibits inflammatory response, plays a pivotal role in kidney diseases. In this study, LPS treated mice and cultured podocytes were used as the models of SA‐AKI in vivo and in vitro, respectively. Medium‐ and high‐dose TAC administration significantly attenuated renal function and renal pathological manifestations at 12, 24 and 48 h after LPS treatment in mice. Moreover, the Toll‐like receptor 4 (TLR4)/myeloid differential protein‐88 (MyD88)/nuclear factor‐kappa (NF‐κB) signalling pathway was also dramatically inhibited by medium‐ and high‐dose TAC administration at 12, 24 and 48 h of LPS treatment mice. In addition, TAC reversed LPS‐induced podocyte cytoskeletal injury and podocyte migratory capability. Our findings indicate that TAC has protective effects against LPS‐induced AKI by inhibiting TLR4/MyD88/NF‐κB signalling pathway and podocyte dysfunction, providing another potential therapeutic effects for the LPS‐induced SA‐AKI.     

3,4‐Dihydroxyphenylethanol and 3,4‐dihydroxyphenylacetic acid affect the aggregation process of E46K variant of α‐synuclein at different extent: Insights into the interplay between protein dynamics and catechol effect

Parkinson''s disease (PD) is a chronic multifactorial disease, whose etiology is not completely understood. The amyloid aggregation of α‐synuclein (Syn) is considered a major cause in the development of the disease. The presence of genetic mutations can boost the aggregation of the protein and the likelihood to develop PD. These mutations can lead to early onset (A30P, E46K, and A53T) or late‐onset (H50Q) forms of PD. The disease is also linked to an increase in oxidative stress and altered levels of dopamine metabolites. The molecular interaction of these molecules with Syn has been previously studied, while their effect on the pathological mutant structure and function is not completely clarified. By using biochemical and biophysical approaches, here we have studied the interaction of the familial variant E46K with two dopamine‐derived catechols, 3,4‐dihydroxyphenylacetic acid and 3,4‐dihydroxyphenylethanol. We show that the presence of these catechols causes a decrease in the formation of amyloid fibrils in a dose‐dependent manner. Native‐ and Hydrogen/deuterium exchange‐mass spectrometry (HDX‐MS) provide evidence that this effect is strongly conformation dependent. Indeed, these molecules interact differently with the interconverting conformers of Syn and its familial variant E46K in solution, selecting the most prone‐to‐aggregation one, confining it into an off‐pathway oligomer. These findings suggest that catechols could be a molecular scaffold for the design of compounds potentially useful in the treatment of Parkinson''s disease and related conditions.     

Anti‐inflammatory properties of lemon‐derived extracellular vesicles are achieved through the inhibition of ERK/NF‐κB signalling pathways

Chronic inflammation is associated with the occurrence of several diseases. However, the side effects of anti‐inflammatory drugs prompt the identification of new therapeutic strategies. Plant‐derived extracellular vesicles (PDEVs) are gaining increasing interest in the scientific community for their biological properties. We isolated PDEVs from the juice of Citrus limon L. (LEVs) and characterized their flavonoid, limonoid and lipid contents through reversed‐phase high‐performance liquid chromatography coupled to electrospray ionization quadrupole time‐of‐flight mass spectrometry (RP‐HPLC–ESI‐Q‐TOF‐MS). To investigate whether LEVs have a protective role on the inflammatory process, murine and primary human macrophages were pre‐treated with LEVs for 24 h and then were stimulated with lipopolysaccharide (LPS). We found that pre‐treatment with LEVs decreased gene and protein expression of pro‐inflammatory cytokines, such as IL‐6, IL1‐β and TNF‐α, and reduced the nuclear translocation and phosphorylation of NF‐κB in LPS‐stimulated murine macrophages. The inhibition of NF‐κB activation was associated with the reduction in ERK1‐2 phosphorylation. Furthermore, the ability of LEVs to decrease pro‐inflammatory cytokines and increase anti‐inflammatory molecules was confirmed ex vivo in human primary T lymphocytes. In conclusion, we demonstrated that LEVs exert anti‐inflammatory effects both in vitro and ex vivo by inhibiting the ERK1‐2/NF‐κB signalling pathway.     

Asprin‐loaded strontium‐containing α‐calcium sulphate hemihydrate/nano‐hydroxyapatite composite promotes regeneration of critical bone defects

Our laboratory originally synthesized strontium(Sr)‐containing α‐calcium sulphate hemihydrate/nano‐hydroxyapatite composite (Sr‐α‐CSH/n‐HA) and demonstrated its ability to repair critical bone defects. This study attempted to incorporate aspirin into it to produce a better bone graft material for critical bone defects. After 5% Sr‐α‐CSH was prepared by coprecipitation and hydrothermal methods, it was mixed with aspirin solution of different concentrations (50 μg/ml, 200 μg/ml, 800 μg/ml and 3200 μg/ml) at a fixed liquid‐solid ratio (0.54 v/w) to obtain aspirin‐loaded Sr‐α‐CSH/n‐HA composite. In vitro experiments were performed on the composite extracts. The tibial defects (3 mm*5 mm) in SD rat model were filled with the composite for 4 weeks and 12 weeks to evaluate its osteogenic capacity in vivo. Our results showed its capability of proliferation, migration and osteogenesis of BMSCs in vitro got improved. In vivo treatment with 800 μg/ml aspirin–loaded Sr‐α‐CSH/n‐HA composite led to significantly more new bone formation in the defects compared with Sr‐α‐CSH/n‐HA composite and significantly promoted the expression of osteogenic‐related genes and inhibited osteoclast activity. In general, our research suggests that aspirin‐loaded Sr‐α‐CSH/n‐HA composite may have a greater capacity of repairing tibial defects in SD rats than simple Sr‐α‐CSH/n‐HA composite.     

5‐Fluorouracil efficacy requires anti‐tumor immunity triggered by cancer‐cell‐intrinsic STING

5‐Fluorouracil (5‐FU) is a widely used chemotherapeutic drug, but the mechanisms underlying 5‐FU efficacy in immunocompetent hosts in vivo remain largely elusive. Through modeling 5‐FU response of murine colon and melanoma tumors, we report that effective reduction of tumor burden by 5‐FU is dependent on anti‐tumor immunity triggered by the activation of cancer‐cell‐intrinsic STING. While the loss of STING does not induce 5‐FU resistance in vitro, effective 5‐FU responsiveness in vivo requires cancer‐cell‐intrinsic cGAS, STING, and subsequent type I interferon (IFN) production, as well as IFN‐sensing by bone‐marrow‐derived cells. In the absence of cancer‐cell‐intrinsic STING, a much higher dose of 5‐FU is needed to reduce tumor burden. 5‐FU treatment leads to increased intratumoral T cells, and T‐cell depletion significantly reduces the efficacy of 5‐FU in vivo. In human colorectal specimens, higher STING expression is associated with better survival and responsiveness to chemotherapy. Our results support a model in which 5‐FU triggers cancer‐cell‐initiated anti‐tumor immunity to reduce tumor burden, and our findings could be harnessed to improve therapeutic effectiveness and toxicity for colon and other cancers.     

A NR2E1‐interacting peptide of LSD1 inhibits the proliferation of brain tumour initiating cells

ObjectivesElimination of brain tumour initiating cells (BTICs) is important for the good prognosis of malignant brain tumour treatment. To develop a novel strategy targeting BTICs, we studied NR2E1(TLX) involved self‐renewal mechanism of BTICs and explored the intervention means.Materials and MethodsNR2E1 and its interacting protein‐LSD1 in BTICs were studied by gene interference combined with cell growth, tumour sphere formation, co‐immunoprecipitation and chromatin immunoprecipitation assays. NR2E1 interacting peptide of LSD1 was identified by Amide Hydrogen/Deuterium Exchange and Mass Spectrometry (HDX‐MS) and analysed by in vitro functional assays. The in vivo function of the peptide was examined with intracranial mouse model by transplanting patient‐derived BTICs.ResultsWe found NR2E1 recruits LSD1, a lysine demethylase, to demethylate mono‐ and di‐methylated histone 3 Lys4 (H3K4me/me2) at the Pten promoter and repress its expression, thereby promoting BTIC proliferation. Using Amide Hydrogen/Deuterium Exchange and Mass Spectrometry (HDX‐MS) method, we identified four LSD1 peptides that may interact with NR2E1. One of the peptides, LSD1‐197‐211 that locates at the LSD1 SWIRM domain, strongly inhibited BTIC proliferation by promoting Pten expression through interfering NR2E1 and LSD1 function. Furthermore, overexpression of this peptide in human BTICs can inhibit intracranial tumour formation.ConclusionPeptide LSD1‐197‐211 can repress BTICs by interfering the synergistic function of NR2E1 and LSD1 and may be a promising lead peptide for brain tumour therapy in future.

• NR2E1 and LSD1 complex directly represses Pten in BTICs by demethylating H3K4me and H3K4me2 at its promoter, and thus promotes the proliferation of BTICs.
• LSD1‐197‐211 peptide interferes the interaction between NR2E1 and LSD1 and restores Pten expression, which in turn inhibits the proliferation of BTICs.

     

Flavonol‐mediated stabilization of PIN efflux complexes regulates polar auxin transport

The transport of auxin controls the rate, direction and localization of plant growth and development. The course of auxin transport is defined by the polar subcellular localization of the PIN proteins, a family of auxin efflux transporters. However, little is known about the composition and regulation of the PIN protein complex. Here, using blue‐native PAGE and quantitative mass spectrometry, we identify native PIN core transport units as homo‐ and heteromers assembled from PIN1, PIN2, PIN3, PIN4 and PIN7 subunits only. Furthermore, we show that endogenous flavonols stabilize PIN dimers to regulate auxin efflux in the same way as does the auxin transport inhibitor 1‐naphthylphthalamic acid (NPA). This inhibitory mechanism is counteracted both by the natural auxin indole‐3‐acetic acid and by phosphomimetic amino acids introduced into the PIN1 cytoplasmic domain. Our results lend mechanistic insights into an endogenous control mechanism which regulates PIN function and opens the way for a deeper understanding of the protein environment and regulation of the polar auxin transport complex.     

SARS‐CoV‐2–host proteome interactions for antiviral drug discovery

Treatment options for COVID‐19, caused by SARS‐CoV‐2, remain limited. Understanding viral pathogenesis at the molecular level is critical to develop effective therapy. Some recent studies have explored SARS‐CoV‐2–host interactomes and provided great resources for understanding viral replication. However, host proteins that functionally associate with SARS‐CoV‐2 are localized in the corresponding subnetwork within the comprehensive human interactome. Therefore, constructing a downstream network including all potential viral receptors, host cell proteases, and cofactors is necessary and should be used as an additional criterion for the validation of critical host machineries used for viral processing. This study applied both affinity purification mass spectrometry (AP‐MS) and the complementary proximity‐based labeling MS method (BioID‐MS) on 29 viral ORFs and 18 host proteins with potential roles in viral replication to map the interactions relevant to viral processing. The analysis yields a list of 693 hub proteins sharing interactions with both viral baits and host baits and revealed their biological significance for SARS‐CoV‐2. Those hub proteins then served as a rational resource for drug repurposing via a virtual screening approach. The overall process resulted in the suggested repurposing of 59 compounds for 15 protein targets. Furthermore, antiviral effects of some candidate drugs were observed in vitro validation using image‐based drug screen with infectious SARS‐CoV‐2. In addition, our results suggest that the antiviral activity of methotrexate could be associated with its inhibitory effect on specific protein–protein interactions.     

A new family of CRISPR‐type V nucleases with C‐rich PAM recognition

Most CRISPR‐type V nucleases are stimulated to cleave double‐stranded (ds) DNA targets by a T‐rich PAM, which restricts their targeting range. Here, we identify and characterize a new family of type V RNA‐guided nuclease, Cas12l, that exclusively recognizes a C‐rich (5''‐CCY‐3′) PAM. The organization of genes within its CRISPR locus is similar to type II‐B CRISPR‐Cas9 systems, but both sequence analysis and functional studies establish it as a new family of type V effector. Biochemical experiments show that Cas12l nucleases function optimally between 37 and 52°C, depending on the ortholog, and preferentially cut supercoiled DNA. Like other type V nucleases, it exhibits collateral nonspecific ssDNA and ssRNA cleavage activity that is triggered by ssDNA or dsDNA target recognition. Finally, we show that one family member, Asp2Cas12l, functions in a heterologous cellular environment, altogether, suggesting that this new group of CRISPR‐associated nucleases may be harnessed as genome editing reagents.     

High‐altitude adaptation in vertebrates as revealed by mitochondrial genome analyses

The high‐altitude environment may drive vertebrate evolution in a certain way, and vertebrates living in different altitude environments might have different energy requirements. We hypothesized that the high‐altitude environment might impose different influences on vertebrate mitochondrial genomes (mtDNA). We used selection pressure analyses and PIC (phylogenetic independent contrasts) analysis to detect the evolutionary rate of vertebrate mtDNA protein‐coding genes (PCGs) from different altitudes. The results showed that the ratio of nonsynonymous/synonymous substitutions (dN/dS) in the mtDNA PCGs was significantly higher in high‐altitude vertebrates than in low‐altitude vertebrates. The seven rapidly evolving genes were shared by the high‐altitude vertebrates, and only one positive selection gene (ND5 gene) was detected in the high‐altitude vertebrates. Our results suggest the mtDNA evolutionary rate in high‐altitude vertebrates was higher than in low‐altitude vertebrates as their evolution requires more energy in a high‐altitude environment. Our study demonstrates the high‐altitude environment (low atmospheric O₂ levels) drives vertebrate evolution in mtDNA PCGs.