Insights from Drosophila on mitochondrial complex I

Behavioral and Transcriptomic Analyses in the Indoxacarb Response of a Non-Target Damselfly Species

Ischnura senegalensis, which widely spreads in paddy fields, has the potential to be used as a natural predator of insect pests. However, the application of insecticides in the field could pose a threat to the survival of I. senegalensis. Among these pesticides, indoxacarb, an oxadiazine insecticide, is renowned for its broad-spectrum efficacy against numerous insect pests. In this study, we examined the toxicity of indoxacarb towards the larvae of I. senegalensis. Behavioral experiments and transcriptome analyses were conducted under indoxacarb treatments. Results revealed that indoxacarb induced abnormal body gestures and significant locomotory impairments, which could ultimately reduce the survival rate of the larvae in their natural habitat. Moreover, transcriptome analyses indicated that genes related to muscle function were significantly affected. Interestingly, at lower concentrations of indoxacarb (0.004 mg/L), the larvae seem to detoxify the indoxacarb with the aid of the cytochrome P450 gene. However, under higher concentrations (0.4 mg/L), the sensory abilities of the larvae were significantly diminished, and they were unable to degrade the toxicity of indoxacarb. Our study underscores the importance of carefully evaluating the impact of insecticides on non-target predatory insects before their widespread application.

Generation of Mutants from the 57B Region of Drosophila melanogaster

The 57B region of Drosophila melanogaster includes a cluster of the three homeobox genes orthopedia (otp), Drosophila Retinal homeobox (DRx), and homeobrain (hbn). In an attempt to isolate mutants for these genes, we performed an EMS mutagenesis and isolated lethal mutants from the 57B region, among them mutants for otp, DRx, and hbn. With the help of two newly generated deletions from the 57B region, we mapped additional mutants to specific chromosomal intervals and identified several of these mutants from the 57B region molecularly. In addition, we generated mutants for CG15651 and RIC-3 by gene targeting and mutants for the genes CG9344, CG15649, CG15650, and ND-B14.7 using the CRISPR/Cas9 system. We determined the lethality period during development for most isolated mutants. In total, we analysed alleles from nine different genes from the 57B region of Drosophila, which could now be used to further explore the functions of the corresponding genes in the future.

The effects of maternal care on the developmental transcriptome and metatranscriptome of a wild bee

Maternal care acts as a strong environmental stimulus that can induce phenotypic plasticity in animals and may also alter their microbial communities through development. Here, we characterize the developmental metatranscriptome of the small carpenter bee, Ceratina calcarata, across developmental stages and in the presence or absence of mothers. Maternal care had the most influence during early development, with the greatest number and magnitude of differentially expressed genes between maternal care treatments, and enrichment for transcription factors regulating immune response in motherless early larvae. Metatranscriptomic data revealed fungi to be the most abundant group in the microbiome, with Aspergillus the most abundant in early larvae raised without mothers. Finally, integrative analysis between host transcriptome and metatranscriptome highlights several fungi correlating with developmental and immunity genes. Our results provide characterizations of the influence of maternal care on gene expression and the microbiome through development in a wild bee.

Phospholipids can regulate complex I assembly independent of their role in maintaining mitochondrial membrane integrity

Several phospholipid (PL) molecules are intertwined with some mitochondrial complex I (CI) subunits in the membrane domain of CI, but their function is unclear. We report that when the Drosophila melanogaster ortholog of the intramitochondrial PL transporter, STARD7, is severely disrupted, assembly of the oxidative phosphorylation (OXPHOS) system is impaired, and the biogenesis of several CI subcomplexes is hampered. However, intriguingly, a restrained knockdown of STARD7 impairs the incorporation of NDUFS5 and NDUFA1 into the proximal part of the CI membrane domain without directly affecting the incorporation of subunits in the distal part of the membrane domain, OXPHOS complexes already assembled, or mitochondrial cristae integrity. Importantly, the restrained knockdown of STARD7 appears to induce a modest amount of cardiolipin remodeling, indicating that there could be some alteration in the composition of the mitochondrial phospholipidome. We conclude that PLs can regulate CI biogenesis independent of their role in maintaining mitochondrial membrane integrity.

A new metabolic model of Drosophila melanogaster and the integrative analysis of Parkinson's disease

High conservation of the disease-associated genes between flies and humans facilitates the common use of Drosophila melanogaster to study metabolic disorders under controlled laboratory conditions. However, metabolic modeling studies are highly limited for this organism. We here report a comprehensively curated genome-scale metabolic network model of Drosophila using an orthology-based approach. The gene coverage and metabolic information of the draft model derived from a reference human model were expanded via Drosophila-specific KEGG and MetaCyc databases, with several curation steps to avoid metabolic redundancy and stoichiometric inconsistency. Furthermore, we performed literature-based curations to improve gene-reaction associations, subcellular metabolite locations, and various metabolic pathways. The performance of the resulting Drosophila model (8,230 reactions, 6,990 metabolites, and 2,388 genes), iDrosophila1 (https://github.com/SysBioGTU/iDrosophila), was assessed using flux balance analysis in comparison with the other currently available fly models leading to superior or comparable results. We also evaluated the transcriptome-based prediction capacity of iDrosophila1, where differential metabolic pathways during Parkinson's disease could be successfully elucidated. Overall, iDrosophila1 is promising to investigate system-level metabolic alterations in response to genetic and environmental perturbations.

Single-cell RNA-sequencing reveals the transcriptional landscape of ND-42 mediated spermatid elongation via mitochondrial derivative maintenance in Drosophila testes

During spermatogenesis, mitochondria extend along the whole length of spermatid tail and offer a structural platform for microtubule reorganization and synchronized spermatid individualization, that eventually helps to generate mature sperm in Drosophila. However, the regulatory mechanism of spermatid mitochondria during elongation remains largely unknown. Herein, we demonstrated that NADH dehydrogenase (ubiquinone) 42 kDa subunit (ND-42) was essential for male fertility and spermatid elongation in Drosophila. Moreover, ND-42 depletion led to mitochondrial disorders in Drosophila testes. Based on single-cell RNA-sequencing (scRNA-seq), we identified 15 distinct cell clusters, including several unanticipated transitional subpopulations or differentiative stages for testicular germ cell complexity in Drosophila testes. Enrichments of the transcriptional regulatory network in the late-stage cell populations revealed key roles of ND-42 in mitochondria and its related biological processes during spermatid elongation. Notably, we demonstrated that ND-42 depletion led to maintenance defects of the major mitochondrial derivative and the minor mitochondrial derivative by affecting mitochondrial membrane potential and mitochondrial-encoded genes. Our study proposes a novel regulatory mechanism of ND-42 for spermatid mitochondrial derivative maintenance, contributing to a better understanding of spermatid elongation.

Mitochondrial Neurodegeneration: Lessons from Drosophila melanogaster Models

The fruit fly-i.e., Drosophila melanogaster-has proven to be a very useful model for the understanding of basic physiological processes, such as development or ageing. The availability of straightforward genetic tools that can be used to produce engineered individuals makes this model extremely interesting for the understanding of the mechanisms underlying genetic diseases in physiological models. Mitochondrial diseases are a group of yet-incurable genetic disorders characterized by the malfunction of the oxidative phosphorylation system (OXPHOS), which is the highly conserved energy transformation system present in mitochondria. The generation of D. melanogaster models of mitochondrial disease started relatively recently but has already provided relevant information about the molecular mechanisms and pathological consequences of mitochondrial dysfunction. Here, we provide an overview of such models and highlight the relevance of D. melanogaster as a model to study mitochondrial disorders.

Drosophila as a diet discovery tool for treating amino acid disorders

Amino acid disorders (AADs) are a large group of rare inherited conditions that collectively impact one in 6500 live births, often resulting in rapid neurological decline and death during infancy. For several AADs, including phenylketonuria, dietary modification prevents physiological deterioration and ameliorates symptoms. Despite this remarkable potential for treatment success, dietary therapy for most AADs remains largely unexplored. Although animal models have provided novel insights into AAD mechanisms, few have been used for therapeutic diet discovery. Here, we find that of all the animal models, Drosophila is particularly well suited for nutrigenomic disease modelling, having amino acid pathways conserved with humans, exceptional genetic tractability, and the unique availability of a synthetic customisable diet.

Mutations in Complex I of the Mitochondrial Electron-Transport Chain Sensitize the Fruit Fly (Drosophila melanogaster) to Ether and Non-Ether Volatile Anesthetics

The mitochondrial electron transport chain (mETC) contains molecular targets of volatile general anesthetics (VGAs), which places carriers of mutations at risk for anesthetic complications. The ND-23⁶⁰¹¹⁴ and mt:ND2^del1 lines of fruit flies (Drosophila melanogaster) that carry mutations in core subunits of Complex I of the mETC replicate numerous characteristics of Leigh syndrome (LS) caused by orthologous mutations in mammals and serve as models of LS. ND-23⁶⁰¹¹⁴ flies are behaviorally hypersensitive to volatile anesthetic ethers and develop an age- and oxygen-dependent anesthetic-induced neurotoxicity (AiN) phenotype after exposure to isoflurane but not to the related anesthetic sevoflurane. The goal of this paper was to investigate whether the alkane volatile anesthetic halothane and other mutations in Complex I and in Complexes II-V of the mETC cause AiN. We found that (i) ND-23⁶⁰¹¹⁴ and mt:ND2^del1 were susceptible to toxicity from halothane; (ii) in wild-type flies, halothane was toxic under anoxic conditions; (iii) alleles of accessory subunits of Complex I predisposed to AiN; and (iv) mutations in Complexes II-V did not result in an AiN phenotype. We conclude that AiN is neither limited to ether anesthetics nor exclusive to mutations in core subunits of Complex I.

Cryo-EM structures of mitochondrial respiratory complex I from Drosophila melanogaster

Respiratory complex I powers ATP synthesis by oxidative phosphorylation, exploiting the energy from NADH oxidation by ubiquinone to drive protons across an energy-transducing membrane. Drosophila melanogaster is a candidate model organism for complex I due to its high evolutionary conservation with the mammalian enzyme, well-developed genetic toolkit, and complex physiology for studies in specific cell types and tissues. Here, we isolate complex I from Drosophila and determine its structure, revealing a 43-subunit assembly with high structural homology to its 45-subunit mammalian counterpart, including a hitherto unknown homologue to subunit NDUFA3. The major conformational state of the Drosophila enzyme is the mammalian-type 'ready-to-go' active resting state, with a fully ordered and enclosed ubiquinone-binding site, but a subtly altered global conformation related to changes in subunit ND6. The mammalian-type 'deactive' pronounced resting state is not observed: in two minor states, the ubiquinone-binding site is unchanged, but a deactive-type π-bulge is present in ND6-TMH3. Our detailed structural knowledge of Drosophila complex I provides a foundation for new approaches to disentangle mechanisms of complex I catalysis and regulation in bioenergetics and physiology.

The membrane domain of respiratory complex I accumulates during muscle aging in Drosophila melanogaster

The boot-shaped respiratory complex I (CI) consists of a mitochondrial matrix and membrane domain organized into N-, Q- and P-modules. The N-module is the most distal part of the matrix domain, whereas the Q-module is situated between the N-module and the membrane domain. The proton-pumping P-module is situated in the membrane domain. We explored the effect of aging on the disintegration of CI and its constituent subcomplexes and modules in Drosophila flight muscles. We find that the fully-assembled complex remains largely intact in aged flies. And while the effect of aging on the stability of many Q- and N-module subunits in subcomplexes was stochastic, NDUFS3 was consistently down-regulated in subcomplexes with age. This was associated with an accumulation of many P-module subunits in subcomplexes. The potential significance of these studies is that genetic manipulations aimed at boosting, perhaps, a few CI subunits may suffice to restore the whole CI biosynthesis pathway during muscle aging.

Proteomic and metabolomic analyses reveal the novel targets of spermine for alleviating diabetic cardiomyopathy in type II diabetic mice

Diabetic cardiomyopathy (DCM) is one of the most serious complications of diabetes. Recent cardiology studies suggest that spermine has a cardioprotective effect. Here, we used proteomic and metabolomic analyses to reveal the underlying research targets in a type II diabetic (T2D) mouse model treated with spermine. Left ventricular tissues from nine mice (Control group, three; T2D group, three; T2D+SP group, three) were excised and analyzed. Quantitative analysis of the global proteome and metabolome was performed using the 4D label-free technique and untargeted metabolomics, respectively, and differentially expressed proteins (DEPs) and metabolites were used to perform bioinformatic analyses. A total of 169 DEPs were identified in T2D/Control group, including 115 upregulated and 54 downregulated proteins. Furthermore, 16 DEPs were identified in T2D+SP/T2D group, where these DEPs were found highly enriched in the cellular, metabolic processes, biological regulation, response to stimulus, and immune system process. The results of association analysis between proteomics and metabolomics showed that SP could affect the production of 51 metabolites by regulating the expression of 16 DEPs in the T2D+SP/T2D group. We also found that PRKG1 was closely related to the expressions of 10 overlapping metabolites between db/db and SP-treated mice. Our findings provide insights into the underlying mechanisms for DCM and suggest the potential applicability of utilizing spermine on protecting against DCM-associated cardiac function deterioration.

Mitochondrial electron transport chain defects modify Parkinson's disease phenotypes in a Drosophila model

INTRODUCTION

Mitochondrial defects have been implicated in Parkinson's disease (PD) since complex I poisons were found to cause accelerated parkinsonism in young people in the early 1980s. More evidence of mitochondrial involvement arose when many of the genes whose mutations caused inherited PD were discovered to be subcellularly localized to mitochondria or have mitochondrial functions. However, the details of how mitochondrial dysfunction might impact or cause PD remain unclear. The aim of our study was to better understand mitochondrial dysfunction in PD by evaluating mitochondrial respiratory complex mutations in a Drosophila melanogaster (fruit fly) model of PD.

METHODS

We have conducted a targeted heterozygous enhancer/suppressor screen using Drosophila mutations within mitochondrial electron transport chain (ETC) genes against a null PD mutation in parkin. The interactions were assessed by climbing assays at 2-5 days as an indicator of motor function. A strong enhancer mutation in COX5A was examined further for L-dopa rescue, oxygen consumption, mitochondrial content, and reactive oxygen species. A later timepoint of 16-20 days was also investigated for both COX5A and a suppressor mutation in cyclope. Generalized Linear Models and similar statistical tests were used to verify significance of the findings.

RESULTS

We have discovered that mutations in individual genes for subunits within the mitochondrial respiratory complexes have interactions with parkin, while others do not, irrespective of complex. One intriguing mutation in a complex IV subunit (cyclope) shows a suppressor rescue effect at early time points, improving the gross motor defects caused by the PD mutation, providing a strong candidate for drug discovery. Most mutations, however, show varying degrees of enhancement or slight suppression of the PD phenotypes. Thus, individual mitochondrial mutations within different oxidative phosphorylation complexes have different interactions with PD with regard to degree and direction. Upon further investigation of the strongest enhancer (COX5A), the mechanism by which these interactions occur initially does not appear to be based on defects in ATP production, but rather may be related to increased levels of reactive oxygen species.

CONCLUSIONS

Our work highlights some key subunits potentially involved in mechanisms underlying PD pathogenesis, implicating ETC complexes other than complex I in PD.

IDH2-mediated regulation of the biogenesis of the oxidative phosphorylation system

Several subunits in the matrix domain of mitochondrial complex I (CI) have been posited to be redox sensors for CI, but how elevated levels of reactive oxygen species (ROS) impinge on CI assembly is unknown. We report that genetic disruption of the mitochondrial NADPH-generating enzyme, isocitrate dehydrogenase 2 (IDH2), in Drosophila flight muscles results in elevated ROS levels and impairment of assembly of the oxidative phosphorylation system (OXPHOS). Mechanistically, this begins with an inhibition of biosynthesis of the matrix domain of CI and progresses to involve multiple OXPHOS complexes. Despite activation of multiple compensatory mechanisms, including enhanced coenzyme Q biosynthesis and the mitochondrial unfolded protein response, ferroptotic cell death ensues. Disruption of enzymes that eliminate hydrogen peroxide, but not those that eliminate the superoxide radical, recapitulates the phenotype, thereby implicating hydrogen peroxide as the signaling molecule involved. Thus, IDH2 modulates the assembly of the matrix domain of CI and ultimately that of the entire OXPHOS.

Physiological Signaling Functions of Reactive Oxygen Species in Stem Cells: From Flies to Man

Reactive oxygen species (ROS), superoxide anion and hydrogen peroxide, are generated as byproducts of oxidative phosphorylation in the mitochondria or via cell signaling-induced NADPH oxidases in the cytosol. In the recent two decades, a plethora of studies established that elevated ROS levels generated by oxidative eustress are crucial physiological mediators of many cellular and developmental processes. In this review, we discuss the mechanisms of ROS generation and regulation, current understanding of ROS functions in the maintenance of adult and embryonic stem cells, as well as in the process of cell reprogramming to a pluripotent state. Recently discovered cell-non-autonomous ROS functions mediated by growth factors are crucial for controlling cell differentiation and cellular immune response in Drosophila. Importantly, many physiological functions of ROS discovered in Drosophila may allow for deciphering and understanding analogous processes in human, which could potentially lead to the development of novel therapeutic approaches in ROS-associated diseases treatment.

Genome of the pincer wasp Gonatopus flavifemur reveals unique venom evolution and a dual adaptation to parasitism and predation

Background

Hymenoptera comprise extremely diverse insect species with extensive variation in their life histories. The Dryinidae, a family of solitary wasps of Hymenoptera, have evolved innovations that allow them to hunt using venom and a pair of chelae developed from the fore legs that can grasp prey. Dryinidae larvae are also parasitoids of Auchenorrhyncha, a group including common pests such as planthoppers and leafhoppers. Both of these traits make them effective and valuable for pest control, but little is yet known about the genetic basis of its dual adaptation to parasitism and predation.

Results

We sequenced and assembled a high-quality genome of the dryinid wasp Gonatopus flavifemur, which at 636.5 Mb is larger than most hymenopterans. The expansion of transposable elements, especially DNA transposons, is a major contributor to the genome size enlargement. Our genome-wide screens reveal a number of positively selected genes and rapidly evolving proteins involved in energy production and motor activity, which may contribute to the predatory adaptation of dryinid wasp. We further show that three female-biased, reproductive-associated yellow genes, in response to the prey feeding behavior, are significantly elevated in adult females, which may facilitate the egg production. Venom is a powerful weapon for dryinid wasp during parasitism and predation. We therefore analyze the transcriptomes of venom glands and describe specific expansions in venom Idgf-like genes and neprilysin-like genes. Furthermore, we find the LWS2-opsin gene is exclusively expressed in male G. flavifemur, which may contribute to partner searching and mating.

Conclusions

Our results provide new insights into the genome evolution, predatory adaptation, venom evolution, and sex-biased genes in G. flavifemur, and present genomic resources for future in-depth comparative analyses of hymenopterans that may benefit pest control.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01081-6.

Dissecting the concordant and disparate roles of NDUFAF3 and NDUFAF4 in mitochondrial complex I biogenesis

Distinct sub-assemblies (modules) of mitochondrial complex I (CI) are assembled with the assistance of CI Assembly Factors (CIAFs) through mechanisms that are incompletely defined. Here, using genetic analyses in Drosophila, we report that when either of the CIAFs – NDUFAF3 or NDUFAF4 – is disrupted, biogenesis of the Q-, N-, and P_P-b-modules of CI is impaired. This is due, at least in part, to the compromised integration of NDUFS3 and NDUFS5 into the Q-, and P_P-b-modules, respectively, coupled with a destabilization of another CIAF, TIMMDC1, in assembly intermediates. Notably, forced expression of NDUFAF4 rescues the biogenesis defects in the Q-module and some aspects of the defects in the P_P-b-module of CI when NDUFAF3 is disrupted. Altogether, our studies furnish new fundamental insights into the mechanism by which NDUFAF3 and NDUFAF4 regulate CI assembly and raises the possibility that certain point mutations in NDUFAF3 may be rescued by overexpression of NDUFAF4.

•

Disruption of NDUFAF3 and NDUFAF4 in Drosophila muscles destabilizes TIMMDC1

•

NDUFAF3 and NDUFAF4 regulate biogenesis of the N, Q, and Pp modules

•

NDUFAF4 ameliorates some of the CI biogenesis defects in NDUFAF3 mutants

Proteomics and transcriptomics jointly identify the key role of oxidative phosphorylation in fluoride-induced myocardial mitochondrial dysfunction in rats

The regulation of mitochondrial function, which is dominated by oxidative phosphorylation (OXPHOs), is important in fluoride induced cardiovascular disease. Based on the previous study of fluoride-induced mitochondrial structure and membrane potential abnormalities, this study integrated ITRAQ protein quantification and RNA-Seq methods to analyze the sequencing data of rat myocardial tissue under fluoride exposure (0, 30, 60 and 90 mg/L). A total of 22 differentially expressed genes associated with the OXPHOs pathway were screened by Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) co-enrichment analysis, and were localizated by Interaction Network and calculated inter-genes and inter-omics correlations by Pearson correlation. In general, fluoride exposure can down-regulate genes related OXPHOs, particularly affecting the assembly of the complex I including Ndufa10, resulting in abnormal mitochondrial ATP synthesis and reduced myocardial energy supply. Most importantly, this study shows that the enriched information from the proteomics can explain the change process of energy production, but the specific molecules involved in energy supply cannot be obtained via transcriptomics information alone. Based on the results of transcriptional and protein analysis, our findings contribute to an innovative understanding of the pathways and molecular changes of myocardial injury induced by fluorosis.

An antibody toolbox to track complex I assembly defines AIF's mitochondrial function

An ability to comprehensively track the assembly intermediates (AIs) of complex I (CI) biogenesis in Drosophila will enable the characterization of the precise mechanism(s) by which various CI regulators modulate CI assembly. Accordingly, we generated 21 novel antibodies to various mitochondrial proteins and used this resource to characterize the mechanism by which apoptosis-inducing factor (AIF) regulates CI biogenesis by tracking the AI profile observed when AIF expression is impaired. We find that when the AIF–Mia40 translocation complex is disrupted, the part of CI that transfers electrons to ubiquinone is synthesized but fails to progress in the CI biosynthetic pathway. This is associated with a reduction in intramitochondrial accumulation of the Mia40 substrate, MIC19. Importantly, knockdown of either MIC19 or MIC60, components of the mitochondrial contact site and cristae organizing system (MICOS), fully recapitulates the AI profile observed when AIF is inhibited. Thus, AIF’s effect on CI assembly is principally due to compromised intramitochondrial transport of the MICOS complex.

Downregulation of the tyrosine degradation pathway extends Drosophila lifespan

Aging is characterized by extensive metabolic reprogramming. To identify metabolic pathways associated with aging, we analyzed age-dependent changes in the metabolomes of long-lived Drosophila melanogaster. Among the metabolites that changed, levels of tyrosine were increased with age in long-lived flies. We demonstrate that the levels of enzymes in the tyrosine degradation pathway increase with age in wild-type flies. Whole-body and neuronal-specific downregulation of enzymes in the tyrosine degradation pathway significantly extends Drosophila lifespan, causes alterations of metabolites associated with increased lifespan, and upregulates the levels of tyrosine-derived neuromediators. Moreover, feeding wild-type flies with tyrosine increased their lifespan. Mechanistically, we show that suppression of ETC complex I drives the upregulation of enzymes in the tyrosine degradation pathway, an effect that can be rescued by tigecycline, an FDA-approved drug that specifically suppresses mitochondrial translation. In addition, tyrosine supplementation partially rescued lifespan of flies with ETC complex I suppression. Altogether, our study highlights the tyrosine degradation pathway as a regulator of longevity.

Evolution of Oxidative Phosphorylation (OXPHOS) Genes Reflecting the Evolutionary and Life Histories of Fig Wasps (Hymenoptera, Chalcidoidea)

Fig wasps are a peculiar group of insects which, for millions of years, have inhabited the enclosed syconia of fig trees. Considering the relatively closed and dark environment of fig syconia, we hypothesize that the fig wasps’ oxidative phosphorylation (OXPHOS) pathway, which is the main oxygen consumption and adenosine triphosphate (ATP) production system, may have adaptively evolved. In this study, we manually annotated the OXPHOS genes of 11 species of fig wasps, and compared the evolutionary patterns of OXPHOS genes for six pollinators and five non-pollinators. Thirteen mitochondrial protein-coding genes and 30 nuclear-coding single-copy orthologous genes were used to analyze the amino acid substitution rate and natural selection. The results showed high amino acid substitution rates of both mitochondrial and nuclear OXPHOS genes in fig wasps, implying the co-evolution of mitochondrial and nuclear genes. Our results further revealed that the OXPHOS-related genes evolved significantly faster in pollinators than in non-pollinators, and five genes had significant positive selection signals in the pollinator lineage, indicating that OXPHOS genes play an important role in the adaptation of pollinators. This study can help us understand the relationship between gene evolution and environmental adaptation.