Engineering TALE-linked deaminases to facilitate precision adenine base editing in mitochondrial DNA.

DddA-derived cytosine base editors (DdCBEs) and transcription activator-like effector (TALE)-linked deaminases (TALEDs) catalyze targeted base editing of mitochondrial DNA (mtDNA) in eukaryotic cells, a method useful for modeling of mitochondrial genetic disorders and developing novel therapeutic modalities. Here, we report that A-to-G-editing TALEDs but not C-to-T-editing DdCBEs induce tens of thousands of transcriptome-wide off-target edits in human cells. To avoid these unwanted RNA edits, we engineered the substrate-binding site in TadA8e, the deoxy-adenine deaminase in TALEDs, and created TALED variants with fine-tuned deaminase activity. Our engineered TALED variants not only reduced RNA off-target edits by >99% but also minimized off-target mtDNA mutations and bystander edits at a target site. Unlike wild-type versions, our TALED variants were not cytotoxic and did not cause developmental arrest of mouse embryos. As a result, we obtained mice with pathogenic mtDNA mutations, associated with Leigh syndrome, which showed reduced heart rates.

Mapping information-rich genotype-phenotype landscapes with genome-scale Perturb-seq.

A central goal of genetics is to define the relationships between genotypes and phenotypes. High-content phenotypic screens such as Perturb-seq (CRISPR-based screens with single-cell RNA-sequencing readouts) enable massively parallel functional genomic mapping but, to date, have been used at limited scales. Here, we perform genome-scale Perturb-seq targeting all expressed genes with CRISPR interference (CRISPRi) across >2.5 million human cells. We use transcriptional phenotypes to predict the function of poorly characterized genes, uncovering new regulators of ribosome biogenesis (including CCDC86, ZNF236, and SPATA5L1), transcription (C7orf26), and mitochondrial respiration (TMEM242). In addition to assigning gene function, single-cell transcriptional phenotypes allow for in-depth dissection of complex cellular phenomena-from RNA processing to differentiation. We leverage this ability to systematically identify genetic drivers and consequences of aneuploidy and to discover an unanticipated layer of stress-specific regulation of the mitochondrial genome. Our information-rich genotype-phenotype map reveals a multidimensional portrait of gene and cellular function.

DNA segregation in mitochondria and beyond: insights from the trypanosomal tripartite attachment complex.

The tripartite attachment complex (TAC) of the single mitochondrion of trypanosomes allows precise segregation of its single nucleoid mitochondrial genome during cytokinesis. It couples the segregation of the duplicated mitochondrial genome to the segregation of the basal bodies of the flagella. Here, we provide a model of the molecular architecture of the TAC that explains how its eight essential subunits connect the basal body, across the mitochondrial membranes, with the mitochondrial genome. We also discuss how the TAC subunits are imported into the mitochondrion and how they assemble to form a new TAC. Finally, we present a comparative analysis of the trypanosomal TAC with open and closed mitotic spindles, which reveals conserved concepts between these diverse DNA segregation systems.

mitoSplitter: A mitochondrial variants-based method for efficient demultiplexing of pooled single-cell RNA-seq.

Single-cell RNA-seq (scRNA-seq) analysis of multiple samples separately can be costly and lead to batch effects. Exogenous barcodes or genome-wide RNA mutations can be used to demultiplex pooled scRNA-seq data, but they are experimentally or computationally challenging and limited in scope. Mitochondrial genomes are small but diverse, providing concise genotype information. We developed "mitoSplitter," an algorithm that demultiplexes samples using mitochondrial RNA (mtRNA) variants, and demonstrated that mtRNA variants can be used to demultiplex large-scale scRNA-seq data. Using affordable computational resources, mitoSplitter can accurately analyze 10 samples and 60,000 cells in 6 h. To avoid the batch effects from separated experiments, we applied mitoSplitter to analyze the responses of five non-small cell lung cancer cell lines to BET (Bromodomain and extraterminal) chemical degradation in a multiplexed fashion. We found the synthetic lethality of TOP2A inhibition and BET chemical degradation in BET inhibitor-resistant cells. The result indicates that mitoSplitter can accelerate the application of scRNA-seq assays in biomedical research.

AMPK-FOXO-IP3R signaling pathway mediates neurological and developmental defects caused by mitochondrial DNA mutations.

Pathological mutations in human mitochondrial genomes (mtDNA) can cause a series of neurological, behavioral, and developmental defects, but the underlying molecular mechanisms are poorly understood. We show here that the energy-sensing adenosine monophosphate (AMP)-activated protein kinase (AMPK) signaling pathway plays a key role in mediating similar defects caused by different mtDNA mutations in Caenorhabditis elegans, including loss or reduction of osmotic, chemical and olfactory sensing, locomotion, and associative learning and memory, as well as increased embryonic lethality. mtDNA mutations cause reduced ATP (adenosine triphosphate) levels, activation of C. elegans AMPK AAK-2, and nuclear translocation of the FOXO transcription factor DAF-16. Activated DAF-16 up-regulates the expression of inositol triphosphate receptor ITR-1, an endoplasmic reticulum calcium channel, leading to increased basal cytosolic Ca2+ levels, decreased neuronal responsiveness, compromised synapses, and increased embryonic death. Treatment of mtDNA mutants with vitamin MK-4 restores cellular ATP and cytosolic Ca2+ levels, improves synaptic development, and suppresses sensory and behavioral defects and embryonic death. Our study provides crucial mechanistic insights into neuronal and developmental defects caused by mtDNA mutations and will improve understanding and treatment of related mitochondrial diseases.

Self-assembly of multi-component mitochondrial nucleoids via phase separation.

Mitochondria contain an autonomous and spatially segregated genome. The organizational unit of their genome is the nucleoid, which consists of mitochondrial DNA (mtDNA) and associated architectural proteins. Here, we show that phase separation is the primary physical mechanism for assembly and size control of the mitochondrial nucleoid (mt-nucleoid). The major mtDNA-binding protein TFAM spontaneously phase separates in vitro via weak, multivalent interactions into droplets with slow internal dynamics. TFAM and mtDNA form heterogenous, viscoelastic structures in vitro, which recapitulate the dynamics and behavior of mt-nucleoids in vivo. Mt-nucleoids coalesce into larger droplets in response to various forms of cellular stress, as evidenced by the enlarged and transcriptionally active nucleoids in mitochondria from patients with the premature aging disorder Hutchinson-Gilford Progeria Syndrome (HGPS). Our results point to phase separation as an evolutionarily conserved mechanism of genome organization.

Neither cardiac mitochondrial DNA variation nor copy number contribute to congenital heart disease risk.

The well-established manifestation of mitochondrial mutations in functional cardiac disease (e.g., mitochondrial cardiomyopathy) prompted the hypothesis that mitochondrial DNA (mtDNA) sequence and/or copy number (mtDNAcn) variation contribute to cardiac defects in congenital heart disease (CHD). MtDNAcns were calculated and rare, non-synonymous mtDNA mutations were identified in 1,837 CHD-affected proband-parent trios, 116 CHD-affected singletons, and 114 paired cardiovascular tissue/blood samples. The variant allele fraction (VAF) of heteroplasmic variants in mitochondrial RNA from 257 CHD cardiovascular tissue samples was also calculated. On average, mtDNA from blood had 0.14 rare variants and 52.9 mtDNA copies per nuclear genome per proband. No variation with parental age at proband birth or CHD-affected proband age was seen. mtDNAcns in valve/vessel tissue (320 ± 70) were lower than in atrial tissue (1,080 ± 320, p = 6.8E-21), which were lower than in ventricle tissue (1,340 ± 280, p = 1.4E-4). The frequency of rare variants in CHD-affected individual DNA was indistinguishable from the frequency in an unaffected cohort, and proband mtDNAcns did not vary from those of CHD cohort parents. In both the CHD and the comparison cohorts, mtDNAcns were significantly correlated between mother-child, father-child, and mother-father. mtDNAcns among people with European (mean = 52.0), African (53.0), and Asian haplogroups (53.5) were calculated and were significantly different for European and Asian haplogroups (p = 2.6E-3). Variant heteroplasmic fraction (HF) in blood correlated well with paired cardiovascular tissue HF (r = 0.975) and RNA VAF (r = 0.953), which suggests blood HF is a reasonable proxy for HF in heart tissue. We conclude that mtDNA mutations and mtDNAcns are unlikely to contribute significantly to CHD risk.

Master graph: an essential integrated assembly model for the plant mitogenome based on a graph-based framework.

Unlike the typical single circular structure of most animal mitochondrial genomes (mitogenome), the drastic structural variation of plant mitogenomes is a result of a mixture of molecules of various sizes and structures. Obtaining the full panoramic plant mitogenome is still considered a roadblock in evolutionary biology. In this study, we developed a graph-based sequence assembly toolkit (GSAT) to construct the pan-structural landscape of plant mitogenome with high-quality mitochondrial master graphs (MMGs) for model species including rice (Oryza sativa) and thale cress (Arabidopsis thaliana). The rice and thale cress MMGs have total lengths of 346 562 and 358 041 bp, including 9 and 6 contigs and 12 and 8 links, respectively, and could be further divided into 6 and 3 minimum master circles and 4 and 2 minimum secondary circles separately. The nuclear mitochondrial DNA segments (NUMTs) in thale cress strongly affected the frequency evaluation of the homologous structures in the mitogenome, while the effects of NUMTs in rice were relatively weak. The mitochondrial plastid DNA segments (MTPTs) in both species had no effects on the assessment of the MMGs. All potential recombinant structures were evaluated, and the findings revealed that all, except for nuclear-homologous structures, MMG structures are present at a much higher frequency than non-MMG structures are. Investigations of potential circular and linear molecules further supported multiple dominant structures in the mitogenomes and could be completely summarized in the MMG. Our study provided an efficient and accurate model for assembling and applying graph-based plant mitogenomes to assess their pan-structural variations.

Modeling mitochondrial DNA diseases: from base editing to pluripotent stem-cell-derived organoids.

Mitochondrial DNA (mtDNA) diseases are multi-systemic disorders caused by mutations affecting a fraction or the entirety of mtDNA copies. Currently, there are no approved therapies for the majority of mtDNA diseases. Challenges associated with engineering mtDNA have in fact hindered the study of mtDNA defects. Despite these difficulties, it has been possible to develop valuable cellular and animal models of mtDNA diseases. Here, we describe recent advances in base editing of mtDNA and the generation of three-dimensional organoids from patient-derived human-induced pluripotent stem cells (iPSCs). Together with already available modeling tools, the combination of these novel technologies could allow determining the impact of specific mtDNA mutations in distinct human cell types and might help uncover how mtDNA mutation load segregates during tissue organization. iPSC-derived organoids could also represent a platform for the identification of treatment strategies and for probing the in vitro effectiveness of mtDNA gene therapies. These studies have the potential to increase our mechanistic understanding of mtDNA diseases and may open the way to highly needed and personalized therapeutic interventions.

Ancient endosymbiont-mediated transmission of a selfish gene provides a model for overcoming barriers to gene transfer into animal mitochondrial genomes.

In contrast to bilaterian animals, non-bilaterian mitochondrial genomes contain atypical genes, often attributed to horizontal gene transfer (HGT) as an ad hoc explanation. Although prevalent in plants, HGT into animal mitochondrial genomes is rare, lacking suitable explanatory models for their occurrence. HGT of the mismatch DNA repair gene (mtMutS) from giant viruses to octocoral (soft corals and their kin) mitochondrial genomes provides a model for how barriers to HGT to animal mitochondria may be overcome. A review of the available literature suggests that this HGT was mediated by an alveolate endosymbiont infected with a lysogenic phycodnavirus that enabled insertion of the homing endonuclease containing mtMutS into octocoral mitochondrial genomes. We posit that homing endonuclease domains and similar selfish elements play a crucial role in such inter-domain gene transfers. Understanding the role of selfish genetic elements in HGT has the potential to aid development of tools for manipulating animal mitochondrial DNA.

Mesoscale structure-function relationships in mitochondrial transcriptional condensates.

In live cells, phase separation is thought to organize macromolecules into membraneless structures known as biomolecular condensates. Here, we reconstituted transcription in condensates from purified mitochondrial components using optimized in vitro reaction conditions to probe the structure-function relationships of biomolecular condensates. We find that the core components of the mt-transcription machinery form multiphasic, viscoelastic condensates in vitro. Strikingly, the rates of condensate-mediated transcription are substantially lower than in solution. The condensate-mediated decrease in transcriptional rates is associated with the formation of vesicle-like structures that are driven by the production and accumulation of RNA during transcription. The generation of RNA alters the global phase behavior and organization of transcription components within condensates. Coarse-grained simulations of mesoscale structures at equilibrium show that the components stably assemble into multiphasic condensates and that the vesicles formed in vitro are the result of dynamical arrest. Overall, our findings illustrate the complex phase behavior of transcribing, multicomponent condensates, and they highlight the intimate, bidirectional interplay of structure and function in transcriptional condensates.

Single p197 molecules of the mitochondrial genome segregation system of Trypanosoma brucei determine the distance between basal body and outer membrane.

The tripartite attachment complex (TAC) couples the segregation of the single unit mitochondrial DNA of trypanosomes with the basal body (BB) of the flagellum. Here, we studied the architecture of the exclusion zone filament (EZF) of the TAC, the only known component of which is p197, that connects the BB with the mitochondrial outer membrane (OM). We show that p197 has three domains that are all essential for mitochondrial DNA inheritance. The C terminus of p197 interacts with the mature and probasal body (pro-BB), whereas its N terminus binds to the peripheral OM protein TAC65. The large central region of p197 has a high α-helical content and likely acts as a flexible spacer. Ultrastructure expansion microscopy (U-ExM) of cell lines exclusively expressing p197 versions of different lengths that contain both N- and C-terminal epitope tags demonstrates that full-length p197 alone can bridge the ∼270-nm distance between the BB and the cytosolic face of the OM. Thus U-ExM allows the localization of distinct domains within the same molecules and suggests that p197 is the TAC subunit most proximal to the BB. In addition, U-ExM revealed that p197 acts as a spacer molecule, as two shorter versions of p197, with the repeat domain either removed or replaced by the central domain of the Trypanosoma cruzi p197 ortholog reduced the distance between the BB and the OM in proportion to their predicted molecular weight.

Targeted base editing in the mitochondrial genome of Arabidopsis thaliana.

Beyond their well-known role in respiration, mitochondria of land plants contain biologically essential and/or agriculturally important genes whose function and regulation are not fully understood. Until recently, it has been difficult to analyze these genes or, in the case of crops, to improve their functions, due to a lack of methods for stably modifying plant mitochondrial genomes. In rice, rapeseed, and Arabidopsis thaliana, mitochondria-targeting transcription activator-like effector nucleases (mitoTALENs) have recently been used to disrupt targeted genes in an inheritable and stable manner. However, this technique can also induce large deletions around the targeted sites, as well as cause ectopic homologous recombinations, which can change the sequences and gene order of mitochondrial genomes. Here, we used mitochondria-targeting TALEN-based cytidine deaminase to successfully substitute targeted C:G pairs with T:A pairs in the mitochondrial genomes of plantlets of A. thaliana without causing deletions or changes in genome structure. Expression vectors of the base editor genes were stably introduced into the nuclear genome by the easy-to-use floral dipping method. Some T1 plants had apparent homoplasmic substitutions that were stably inherited by seed progenies, independently of the inheritance of nuclear-introduced genes. As a demonstration of the method, we used it to restore the growth of an organelle transcript processing 87 (otp87) mutant that is defective in the editing of RNA transcripts of the mitochondrial atp1 gene and to identify bases in atp1 that affect the efficiency of RNA editing by OTP87.

Assembly and comparative analysis of the complete mitochondrial genome of white towel gourd (Luffa cylindrica).

Mitochondria play an important role in the energy production of plant cells through independent genetic systems. This study has aimed to assemble and annotate the functions of the mitochondrial (mt) genome of Luffa cylindrica. The mt genome of L. cylindrica contained two chromosomes with lengths of 380,879 bp and 67,982 bp, respectively. Seventy-seven genes including 39 protein-coding genes, 34 tRNA genes, 3 rRNA genes, and 1 pseudogene, were identified. About 90.63% of the codons ended with A or U bases, and 98.63% of monomers contained A/T, which contributed to the high A/T content (55.91%) of the complete mt genome. Six genes (ATP8, CCMFC, NAD4, RPL10, RPL5 and RPS4) showed positive selection. Phylogenetic analysis indicates that L. cylindrica is closely related to L. acutangula. The present results provide the mt genome of L. cylindrica, which may facilitate possible genetic variation, evolutionary, and molecular breeding studies of L. cylindrica.

The role of lysine acetylation in the function of mitochondrial ribosomal protein L12.

Mitochondria play a central role in energy production and cellular metabolism. Mitochondria contain their own small genome (mitochondrial DNA, mtDNA) that carries the genetic instructions for proteins required for ATP synthesis. The mitochondrial proteome, including the mitochondrial transcriptional machinery, is subject to post-translational modifications (PTMs), including acetylation and phosphorylation. We set out to determine whether PTMs of proteins associated with mtDNA may provide a potential mechanism for the regulation of mitochondrial gene expression. Here, we focus on mitochondrial ribosomal protein L12 (MRPL12), which is thought to stabilize mitochondrial RNA polymerase (POLRMT) and promote transcription. Numerous acetylation sites of MRPL12 were identified by mass spectrometry. We employed amino acid mimics of the acetylated (lysine to glutamine mutants) and deacetylated (lysine to arginine mutants) versions of MRPL12 to interrogate the role of lysine acetylation in transcription initiation in vitro and mitochondrial gene expression in HeLa cells. MRPL12 acetyl and deacetyl protein mimics were purified and assessed for their ability to impact mtDNA promoter binding of POLRMT. We analyzed mtDNA content and mitochondrial transcript levels in HeLa cells upon overexpression of acetyl and deacetyl mimics of MRPL12. Our results suggest that MRPL12 single-site acetyl mimics do not change the mtDNA promoter binding ability of POLRMT or mtDNA content in HeLa cells. Individual acetyl mimics may have modest effects on mitochondrial transcript levels. We found that the mitochondrial deacetylase, Sirtuin 3, is capable of deacetylating MRPL12 in vitro, suggesting a potential role for dynamic acetylation controlling MRPL12 function in a role outside of the regulation of gene expression.

Assembly and comparative analysis of the initial complete mitochondrial genome of Primulina hunanensis (Gesneriaceae): a cave-dwelling endangered plant.

BACKGROUND: Primulina hunanensis, a troglobitic plant within the Primulina genus of Gesneriaceae family, exhibits robust resilience to arid conditions and holds great horticultural potential as an ornamental plant. The work of chloroplast genome (cpDNA) has been recently accomplished, however, the mitochondrial genome (mtDNA) that is crucial for plant evolution has not been reported. RESULTS: In this study, we sequenced and assembled the P. hunanensis complete mtDNA, and elucidated its evolutionary and phylogenetic relationships. The assembled mtDNA spans 575,242 bp with 43.54% GC content, encompassing 60 genes, including 37 protein-coding genes (PCGs), 20 tRNA genes, and 3 rRNA genes. Notably, high number of repetitive sequences in the mtDNA and substantial sequence translocation from chloroplasts to mitochondria were observed. To determine the evolutionary and taxonomic positioning of P. hunanensis, a phylogenetic tree was constructed using mitochondrial PCGs from P. hunanensis and 32 other taxa. Furthermore, an exploration of PCGs relative synonymous codon usage, identification of RNA editing events, and an investigation of collinearity with closely related species were conducted. CONCLUSIONS: This study reports the initial assembly and annotation of P. hunanensis mtDNA, contributing to the limited mtDNA repository for Gesneriaceae plants and advancing our understanding of their evolution for improved utilization and conservation.

Mitochondrial genome features and systematic evolution of diospyros kaki thunb 'Taishuu'.

BACKGROUND: 'Taishuu' has a crisp texture, abundant juice, and sweet flavor with hints of cantaloupe. The availability of mitochondrial genome data of Diospyros species is far from the known number of species. RESULTS: The sequencing data were assembled into a closed circular mitochondrial chromosome with a 421,308 bp length and a 45.79% GC content. The mitochondrial genome comprised 40 protein-coding, 24 tRNA, and three rRNA genes. The most common codons for arginine (Arg), proline (Pro), glycine (Gly), tryptophan (Trp), valine (Val), alanine (Ala), and leucine (Leu) were AGA, CCA, GGA, UGG, GUA, GCA, and CUA, respectively. The start codon for cox1 and nad4L protein-coding genes was ACG (ATG), whereas the remaining protein-coding genes started with ATG. There are four types of stop codons: CGA, TAA, TAG, and TGA, with TAA being the most frequently used stop codon (45.24%). In the D. kaki Thunb. 'Taishuu' mitochondrial genome, a total of 645 repeat sequences were identified, including 125 SSRs, 7 tandem repeats, and 513 dispersed repeats. Collinearity analysis revealed a close relationship between D. kaki Thunb. 'Taishuu' and Diospyros oleifera, with conserved homologous gene fragments shared among these species in large regions of the mitochondrial genome. The protein-coding genes ccmB and nad4L were observed to undergo positive selection. Analysis of homologous sequences between chloroplasts and mitochondria identified 28 homologous segments, with a total length of 24,075 bp, accounting for 5.71% of the mitochondrial genome. These homologous segments contain 8 annotated genes, including 6 tRNA genes and 2 protein-coding genes (rrn18 and ccmC). There are 23 homologous genes between chloroplasts and nuclei. Mitochondria, chloroplasts, and nuclei share two homologous genes, which are trnV-GAC and trnW-CCA. CONCLUSION: In conclusion, a high-quality chromosome-level draft genome for D. kaki was generated in this study, which will contribute to further studies of major economic traits in the genus Diospyros.

Insights into the multi-chromosomal mitochondrial genome structure of the xero-halophytic plant Haloxylon Ammodendron (C.A.Mey.) Bunge ex Fenzl.

BACKGROUND: Haloxylon ammodendron holds significance as an ecological plant, showcasing remarkable adaptability to desert conditions, halophytic environments, and sand fixation. With its potential for carbon sequestration, it emerges as a promising candidate for environmental sustainability. Furthermore, it serves as a valuable C4 plant model, offering insights into the genetic foundations of extreme drought tolerance. Despite the availability of plastid and nuclear genomes, the absence of a mitochondrial genome (mitogenome or mtDNA) hinders a comprehensive understanding of its its mtDNA structure, organization, and phylogenetic implications. RESULTS: In the present study, the mitochondrial genome of H. ammodendron was assembled and annotated, resulting in a multi-chromosomal configuration with two circular chromosomes. The mtDNA measured 210,149 bp in length and contained 31 protein-coding genes, 18 tRNA and three rRNA. Our analysis identified a total of 66 simple sequence repeats along with 27 tandem repeats, 312 forward repeats, and 303 palindromic repeats were found. Notably, 17 sequence fragments displayed homology between the mtDNA and chloroplast genome (cpDNA), spanning 5233 bp, accounting for 2.49% of the total mitogenome size. Additionally, we predicted 337 RNA editing sites, all of the C-to-U conversion type. Phylogenetic inference confidently placed H. ammodendron in the Amaranthacea family and its close relative, Suaeda glacum. CONCLUSIONS: H. ammodendron mtDNA showed a multi-chromosomal structure with two fully circularized molecules. This newly characterized mtDNA represents a valuable resource for gaining insights into the basis of mtDNA structure variation within Caryophyllales and the evolution of land plants, contributing to their identification, and classification.

Intraspecific and interspecific variations in the synonymous codon usage in mitochondrial genomes of 8 pleurotus strains.

In this study, we investigated the codon bias of twelve mitochondrial core protein coding genes (PCGs) in eight Pleurotus strains, two of which are from the same species. The results revealed that the codons of all Pleurotus strains had a preference for ending in A/T. Furthermore, the correlation between codon base compositions and codon adaptation index (CAI), codon bias index (CBI) and frequency of optimal codons (FOP) indices was also detected, implying the influence of base composition on codon bias. The two P. ostreatus species were found to have differences in various base bias indicators. The average effective number of codons (ENC) of mitochondrial core PCGs of Pleurotus was found to be less than 35, indicating strong codon preference of mitochondrial core PCGs of Pleurotus. The neutrality plot analysis and PR2-Bias plot analysis further suggested that natural selection plays an important role in Pleurotus codon bias. Additionally, six to ten optimal codons (ΔRSCU > 0.08 and RSCU > 1) were identified in eight Pleurotus strains, with UGU and ACU being the most widely used optimal codons in Pleurotus. Finally, based on the combined mitochondrial sequence and RSCU value, the genetic relationship between different Pleurotus strains was deduced, showing large variations between them. This research has improved our understanding of synonymous codon usage characteristics and evolution of this important fungal group.

Novel insight into lepidopteran phylogenetics from the mitochondrial genome of the apple fruit moth of the family Argyresthiidae.

BACKGROUND: The order Lepidoptera has an abundance of species, including both agriculturally beneficial and detrimental insects. Molecular data has been used to investigate the phylogenetic relationships of major subdivisions in Lepidoptera, which has enhanced our understanding of the evolutionary relationships at the family and superfamily levels. However, the phylogenetic placement of many superfamilies and/or families in this order is still unknown. In this study, we determine the systematic status of the family Argyresthiidae within Lepidoptera and explore its phylogenetic affinities and implications for the evolution of the order. We describe the first mitochondrial (mt) genome from a member of Argyresthiidae, the apple fruit moth Argyresthia conjugella. The insect is an important pest on apples in Fennoscandia, as it switches hosts when the main host fails to produce crops. RESULTS: The mt genome of A. conjugella contains 16,044 bp and encodes all 37 genes commonly found in insect mt genomes, including 13 protein-coding genes (PCGs), two ribosomal RNAs, 22 transfer RNAs, and a large control region (1101 bp). The nucleotide composition was extremely AT-rich (82%). All detected PCGs (13) began with an ATN codon and terminated with a TAA stop codon, except the start codon in cox1 is ATT. All 22 tRNAs had cloverleaf secondary structures, except trnS1, where one of the dihydrouridine (DHU) arms is missing, reflecting potential differences in gene expression. When compared to the mt genomes of 507 other Lepidoptera representing 18 superfamilies and 42 families, phylogenomic analyses found that A. conjugella had the closest relationship with the Plutellidae family (Yponomeutoidea-super family). We also detected a sister relationship between Yponomeutoidea and the superfamily Tineidae. CONCLUSIONS: Our results underline the potential importance of mt genomes in comparative genomic analyses of Lepidoptera species and provide valuable evolutionary insight across the tree of Lepidoptera species.