Mitochondrial DNA Research Articles

Exploring the Mechanisms and Therapeutic Approaches of Mitochondrial Dysfunction in Alzheimer's Disease: An Educational Literature Review.

Alzheimer's disease (AD) presents a significant challenge to global health. It is characterized by progressive cognitive deterioration and increased rates of morbidity and mortality among older adults. Among the various pathophysiologies of AD, mitochondrial dysfunction, encompassing conditions such as increased reactive oxygen production, dysregulated calcium homeostasis, and impaired mitochondrial dynamics, plays a pivotal role. This review comprehensively investigates the mechanisms of mitochondrial dysfunction in AD, focusing on aspects such as glucose metabolism impairment, mitochondrial bioenergetics, calcium signaling, protein tau and amyloid-beta-associated synapse dysfunction, mitophagy, aging, inflammation, mitochondrial DNA, mitochondria-localized microRNAs, genetics, hormones, and the electron transport chain and Krebs cycle. While lecanemab is the only FDA-approved medication to treat AD, we explore various therapeutic modalities for mitigating mitochondrial dysfunction in AD, including antioxidant drugs, antidiabetic agents, acetylcholinesterase inhibitors (FDA-approved to manage symptoms), nutritional supplements, natural products, phenylpropanoids, vaccines, exercise, and other potential treatments.

Unravelling the Molecular Identity of Bulgarian Jumping Plant Lice of the Family Aphalaridae (Hemiptera: Psylloidea)

Psyllids (Hemiptera: Psylloidea) are plant sap-sucking insects whose identification is often difficult for non-experts. Despite the rapid development of DNA barcoding techniques and their widespread use, only a limited number of sequences of psyllids are available in the public databases, and those that are available are often misidentified. Here, we provide 80 sequences of two mitochondrial genes, cytochrome c oxidase I (COI) and cytochrome b (Cytb), for 25 species of Aphalaridae, mainly from Bulgaria. The DNA barcodes for 15 of these species are published for the first time. In cases where standard primers failed to amplify the target gene fragment, we designed new primers that can be used in future studies. The distance-based thresholds for the analysed species were between 0.0015 and 0.3415 for COI and 0.0771 and 0.4721 for Cytb, indicating that the Cytb gene has a higher interspecific divergence, compared to COI, and therefore allows for more accurate species identification. The species delimitation based on DNA barcodes is largely consistent with the differences resulting from morphological and host plant data, demonstrating that the use of DNA barcodes is suitable for successful identification of most aphalarid species studied. The phylogenetic reconstruction based on maximum likelihood and Bayesian inference analyses, while showing similar results at high taxonomic levels to previously published phylogenies, provides additional information on the placement of aphalarids at the species level. The following five species represent new records for Bulgaria: Agonoscena targionii, Aphalara affinis, Colposcenia aliena, Co. bidentata, and Craspedolepta malachitica. Craspedolepta conspersa is reported for the first time from the Czech Republic, while Agonoscena cisti is reported for the first time from Albania.

A new miniaturized species of leaf chameleon, genus Brookesia, from a littoral forest fragment in eastern Madagascar

The number of species of miniaturized chameleons of the genus Brookesia, subgenus Evoluticauda, known to science has witnessed a dramatic increase over the past 15 years, due to the discovery of multiple microendemic species of very strong genetic divergence. So far, no described Evoluticauda species are known from the littoral forest of Madagascar’s east coast, one of the most threatened habitat types of the island. Here, we report on the discovery of a new species of Evoluticauda occurring in the littoral forest at Ankanin’ny Nofy, a touristic site at about sea level and close to Vohibola forest where probably the same species occurs. The new species, Brookesia nofy sp. nov., is sister to B. ramanantsoai which occurs in mid-altitude and highland forests (800–1300 m a.s.l.) at a similar latitude of eastern Madagascar, but differs from this species by an enormous genetic divergence in both mitochondrial and nuclear genes (8.1–9.4% uncorrected pairwise distance in the 16S rRNA gene, 1.8–2.5% in c-mos), a relatively shorter tail, a different arrangement of cephalic ridges, and probably smaller body size and relatively larger head. The discovery of B. nofy sp. nov. emphasizes the importance of conserving Madagascar’s last remaining fragments of littoral forest.

Targeting PRMT3 impairs methylation and oligomerization of HSP60 to boost anti-tumor immunity by activating cGAS/STING signaling

Immune checkpoint blockade (ICB) has emerged as a promising therapeutic option for hepatocellular carcinoma (HCC), but resistance to ICB occurs and patient responses vary. Here, we uncover protein arginine methyltransferase 3 (PRMT3) as a driver for immunotherapy resistance in HCC. We show that PRMT3 expression is induced by ICB-activated T cells via an interferon-gamma (IFNγ)-STAT1 signaling pathway, and higher PRMT3 expression levels correlate with reduced numbers of tumor-infiltrating CD8+ T cells and poorer response to ICB. Genetic depletion or pharmacological inhibition of PRMT3 elicits an influx of T cells into tumors and reduces tumor size in HCC mouse models. Mechanistically, PRMT3 methylates HSP60 at R446 to induce HSP60 oligomerization and maintain mitochondrial homeostasis. Targeting PRMT3-dependent HSP60 methylation disrupts mitochondrial integrity and increases mitochondrial DNA (mtDNA) leakage, which results in cGAS/STING-mediated anti-tumor immunity. Lastly, blocking PRMT3 functions synergize with PD-1 blockade in HCC mouse models. Our study thus identifies PRMT3 as a potential biomarker and therapeutic target to overcome immunotherapy resistance in HCC.

Favorable Outcome after Liver Transplantation in an Infant with Liver Failure Due to Deoxyguanosine Kinase Deficiency

Introduction: Deoxyguanosine Kinase (DGUOK) deficiency is a very rare disorder characterized by liver dysfunction, neurological manifestations, and metabolic disorders secondary to severely reduced mitochondrial DNA content. These patients develop early-onset liver failure, and their liver transplantation (LT) indication remains debatable due to the possibility of neurological involvement. Case Report: We present the case of a 6-month-old female diagnosed with DGUOK deficiency who developed liver failure. At 9 months, she underwent a living-related LT with an initial favorable evolution under immunosuppression therapy with tacrolimus. Four months after LT, she presented two prolonged bacterial and Rotavirus enteritis episodes. She developed classical post-transplant complications (severe renal tubular acidosis type IV, secondary to the high tacrolimus level, and post-transplant lymphoproliferative disease) during these episodes. Her condition deteriorated progressively, with reversible hypotonia and significant weight loss. However, the neurological evaluation did not reveal any signs suggestive of the progression of the underlying disease. A few months later, her clinical features and laboratory parameters improved considerably. Conclusions: This case highlights the unpredictable evolution of children with LT for liver failure due to DGUOK deficiency.

Machine learning algorithm-based biomarker exploration and validation of mitochondria-related diagnostic genes in osteoarthritis.

The role of mitochondria in the pathogenesis of osteoarthritis (OA) is significant. In this study, we aimed to identify diagnostic signature genes associated with OA from a set of mitochondria-related genes (MRGs). First, the gene expression profiles of OA cartilage GSE114007 and GSE57218 were obtained from the Gene Expression Omnibus. And the limma method was used to detect differentially expressed genes (DEGs). Second, the biological functions of the DEGs in OA were investigated using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Wayne plots were employed to visualize the differentially expressed mitochondrial genes (MDEGs) in OA. Subsequently, the LASSO and SVM-RFE algorithms were employed to elucidate potential OA signature genes within the set of MDEGs. As a result, GRPEL and MTFP1 were identified as signature genes. Notably, GRPEL1 exhibited low expression levels in OA samples from both experimental and test group datasets, demonstrating high diagnostic efficacy. Furthermore, RT-qPCR analysis confirmed the reduced expression of Grpel1 in an in vitro OA model. Lastly, ssGSEA analysis revealed alterations in the infiltration abundance of several immune cells in OA cartilage tissue, which exhibited correlation with GRPEL1 expression. Altogether, this study has revealed that GRPEL1 functions as a novel and significant diagnostic indicator for OA by employing two machine learning methodologies. Furthermore, these findings provide fresh perspectives on potential targeted therapeutic interventions in the future.

Cryptic cytoplasmic male sterility-causing gene in the mitochondrial genome of common japonica rice.

Cytoplasmic male sterility (CMS) is an agronomically significant trait that causes dysfunction in pollen and anther development. It is often observed during successive backcrossing between distantly related species. Here, we show that Asian japonica cultivars (Oryza sativa) exhibit CMS when the nucleus is replaced with that of the African rice Oryza glaberrima. The CMS line produced stunted anthers and did not set any seeds. Mitochondrial orf288 RNA was detected in the anthers of CMS lines but not in fertility restorer lines. The mitochondrial genome-edited japonica rice that was depleted of orf288 did not exhibit male sterility when backcrossed with O. glaberrima. These results demonstrate that orf288 is a CMS-causing gene. As orf288 commonly occurs in the mitochondrial genomes of japonica rice, these results indicate that common japonica rice cultivars possess a cryptic CMS-causing gene hidden in their mitochondrial genomes.

Variant load of mitochondrial DNA in single human mesenchymal stem cells

Heteroplasmic mitochondrial DNA (mtDNA) variants accumulate as humans age, particularly in the stem-cell compartments, and are an important contributor to age-related disease. Mitochondrial dysfunction has been observed in osteoporosis and somatic mtDNA pathogenic variants have been observed in animal models of osteoporosis. However, this has never been assessed in the relevant human tissue. Mesenchymal stem cells (MSCs) are the progenitors to many cells of the musculoskeletal system and are critical to skeletal tissues and bone vitality. Investigating mtDNA in MSCs could provide novel insights into the role of mitochondrial dysfunction in osteoporosis. To determine if this is possible, we investigated the landscape of somatic mtDNA variation in MSCs through a combination of fluorescence-activated cell sorting and single-cell next-generation sequencing. Our data show that somatic heteroplasmic variants are present in individual patient-derived MSCs, can reach high heteroplasmic fractions and have the potential to be pathogenic. The identification of somatic heteroplasmic variants in MSCs of patients highlights the potential for mitochondrial dysfunction to contribute to the pathogenesis of osteoporosis.

Photoinduction of Ferroptosis and cGAS-STING Activation by a H2S-Responsive Iridium(III) Complex for Cancer-Specific Therapy.

Triggering ferroptosis represents a promising anticancer therapeutic strategy, but the development of a selective ferroptosis inducer for cancer-specific therapy remains a great challenge. Herein, a H2S-responsive iridium(III) complex NA-Ir has been well-designed as a ferroptosis inducer. NA-Ir could selectively light up H2S-rich cancer cells, primarily localize in mitochondria, intercalate into mitochondrial DNA (mtDNA), and induce mtDNA damage, exhibiting higher anticancer activity under light irradiation. Mechanistic studies showed that NA-Ir-mediated PDT triggered lipid peroxidation and glutathione peroxidase 4 downregulation through ROS production and GSH depletion, resulting in ferroptosis through multiple pathways. Moreover, the intense mtDNA damage can activate the cyclic GMP-AMP synthase-stimulator of the interferon gene (cGAS-STING) pathway, leading to ferritinophagy and further ferroptosis. RNA-sequencing analysis showed that NA-Ir-mediated PDT mainly affects the expression of genes related to ferroptosis, autophagy, and cancer immunity. This study demonstrates the first cancer-specific example with ferroptosis and cGAS-STING activation, which provides a new strategy for multimodal synergistic therapy.

Dominant strain shift in the invasive fall armyworm (Lepidoptera: Noctuidae) populations in Thailand as inferred from mitochondrial COI and nuclear Tpi genes

Abstract The fall armyworm, Spodoptera frugiperda (J. E. Smith, 1797) (Lepidoptera: Noctuidae), is a significant global pest, that exhibits 2 discernible strains, corn strain (CS) and rice strain (RS). After initial detection in the eastern hemisphere in 2016, the dominant strain was identified as RS based only on cytochrome C oxidase subunit I (COI) mitochondrial gene from limited samples from various countries, including Thailand. This study aimed to assess strain and haplotype variation in the S. frugiperda populations in Thailand using both mitochondrial COI and nuclear triosephosphate isomerase (Tpi) genes. Analyses of COI sequences (n = 105) revealed 2 predominant haplotypes, COICSh4 (82.86%) and COIRSh1 (17.14%), and the analyses of Tpi sequences (n = 99) revealed 6 haplotypes, with TpiCa1a (53.53%) being the most prevalent. Of the 98 caterpillar samples, the majority exhibited true CS (83.67%) for both genes. Meanwhile, interstrain hybrids, indicated by gene discordance, accounted for the minority (16.33%). Interestingly, despite the initial dominance of RS during the 2018 outbreak, the current study identified CS as the prevalent strain across all localities in Thailand. These findings suggested a shift in S. frugiperda dynamics in Thailand that was possibly influenced by factors, such as competitive exclusion principle, pesticide usage in rice cultivation, and preferences for corn over rice. Our study suggests a need to reexamine the previous reports of rice-strain dominance in various countries in the eastern hemisphere after the initial invasion.

Comparative Analysis of the Mitochondrial Genome of Eggplant (Solanum melongena L.) to Identify Cytoplasmic Male Sterility Candidate Genes.

Cytoplasmic male sterility (CMS) is important for commercial hybrid seed production. However, it is still not used in eggplant (Solanum melongena L.), and corresponding regulatory genes and mechanisms of action have not been reported. We report CMS line 327A, which was derived from the hybridization between cultivated and wild eggplants. By looking at different stages of anther development under a microscope, we saw that the 327A anther's tapetum layer vacuolized during meiosis, which caused abortion. To investigate the 327A CMS regulatory genes, the mitochondrial genomes of 327A and its maintainer line 327B were assembled de novo. It was found that 15 unique ORFs (Open Reading Frame) were identified in 327A. RT-PCR and RT-QPCAR tests confirmed that orf312a and orf172a, 327A-specific ORFs with a transmembrane domain, were strongly expressed in sterile anthers of 327A. In addition, orf312a has a chimeric structure with the ribosomal protein subunit rpl16. Therefore, orf312a and orf172a can be considered strong candidate genes for CMS. Concurrently, we analyzed the characteristics of CMS to develop a functional molecular marker, CMS312, targeting a future theoretical basis for eggplant CMS three-line molecular breeding.

The complete mitochondrial genome of the banana pathogen Fusarium oxysporum f. sp. cubense M5.

We report the complete mitochondrial genome of a causal agent of banana fusarium wilt isolated in Mexico. The whole set of genes encoding proteins related to respiration and ATP synthesis, rRNAs, tRNAs are enlisted. Two open reading frames of unknown function conserved in Fusarium oxysporum were also identified.

Complete mitochondrial genomes of Dactylogyrus crucifer and Dactylogyrus zandti reveal distinct patterns of codon usage within Dactylogyrus

Monogeneans of the genus Dactylogyrus Diesing, 1850, the largest genus in the family Dactylogyridae, mostly parasitize the gills of cyprinoid hosts; however, only 3 Dactylogyrus’ mitochondrial genomes (mitogenomes) are studied so far. The aim of this research is to extend our understanding of the mitogenomes of Dactylogyrus. We sequenced the mitogenomes of D. crucifer and D. zandti isolated from Rutilus rutilus and Abramis brama orientalis in northwest China, and then we compared these mitogenomes with other monogeneans. We used Illumina NovaSeq to sequence the entire mitochondrial genomes of D. crucifer and D. zandti and characterized the mitogenomes to understand the gene structure, gene identity, the secondary structures of the 22 tRNA genes, and relative synonymous codon usage. We used the analytic Bayesian Information and Maximum Likelihood methods to determine their associated phylogenetic trees. The mitogenomes of D. crucifer and D. zandti were 14,403 and 18,584 bp, respectively. Organization and positioning of these genes were in accordance with Dactylogyrus lamellatus and Dactylogyrus tuba. The nucleotide composition of Dactylogyridae was different from other families of Monogenea, and the A+T count of genus Dactylogyrus (54 – 58.4 %) was lower than other genus species of the family Dactylogyridea (63.9 – 78.4 %) in protein-coding genes. Dactylogyrus members displayed a codon usage bias. The relative synonymous codon used by Dactylogyrus was not conserved and was lower than other monogeneans. The codon use patterns of closely-related species isolated from closely-related hosts were identical. Phylogenetic analyses using mitogenomic dataset produced Dactylogyrus isolated from host subfamily Leuciscinae formed a sister-group. Our results contributed significantly to an increased database of mitogenomes, more than 50 %, for Dactylogyrus that may help future studies of mitochondrial genes and codon uses for the analysis of monogenean phylogenetics.

First Record of the Complete Mitochondrial Genome for the Genus Borbo (Lepidoptera, Hesperiidae): Characterization and Comparative Genomic Analysis

Butterflies of the genus Borbo are mainly distributed in the Oriental and Australian regions and are considered pests of important crops. However, no mitochondrial genomes have been reported for this genus until now, leaving the evolutionary characteristics and differentiation patterns of their mitogenomes unclear. In this study, we present the first complete mitochondrial genome sequence of the rice swift, Borbo cinnara. The circular double-stranded mitogenome was 15,508 bp in length, comprising 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 1 non-coding control region (CR). Among the mitogenomes of Hesperiinae, the ND3 gene was found to be the most variable PCG, while COX1 was the most conserved. Selection pressure analysis revealed that ND3 was under relaxed purifying selection, whereas COX1 was subjected to strong purifying selection. Phylogenetic trees reconstructed using both the Bayesian inference (BI) and maximum likelihood (ML) methods yielded robust and identical topologies, confirming the sister relationship between B. cinnara and Pelopidas mathias at the mitogenome level. Methodologically, this research enriches novel molecular markers for the species identification of butterflies and enhances our understanding of mitogenomic evolution in Lepidoptera.

An independent base composition of each rate class for improved likelihood-based phylogeny estimation; the 5rf model.

The combination of a t ime r eversible Markov process with a "hidden" mixture of g amma distributed relative site rates plus i nvariant sites have become the most favoured options for likelihood and other probabilistic models of nucleotide evolution (e.g., tr4gi which approximates a gamma with four rate classes). However, these models assume a homogeneous and stationary distribution of nucleotide (character or base) frequencies. Here, we explore the potential benefits and pitfalls of allowing each rate category (rate class) of a 4gi mixture model to have its own base frequencies. This is achieved by starting each of the five rate classes, at the tree's r oot, with its own free choice of nucleotide f requencies to create a 4gi5rf model or a 5rf model in shorthand. We assess the practical identifiability of this approach with a BEAST 2 implementation, aiming to determine if it can accurately estimate credibility intervals and expected values for a wide range of plausible parameter values. Practical identifiability, as distinguished from mathematical identifiability, gauges the model's ability to identify parameters in real-world scenarios, as opposed to theoretically with infinite data. One of the most common types of phylogenetic data is mitochondrial DNA (mtDNA) protein coding sequence. It is often assumed current models analyse robustly such data and that higher likelihood/posterior probability models do better. However, this abstract shows that vertebrate mtDNA remains a very difficult type of data to fully model, and that dramatically higher likelihoods do not mean a model is measurably more accurate with respect to recovering key parameters of biological interest (e.g., monophyletic groups, their support and their ages). The 4gi5rf model considerably improves marginal likelihoods and seems to reverse some apparent errors exacerbated by the 4gi model, while introducing others. Problems appear to be linked to non-stationary DNA repair processes that alter the mutation/substitution spectra across lineages and time. We also show such problems are not unique to mtDNA and are encountered in analysing nuclear sequences. Non-stationarity of DNA repair processes mutation/substitution spectra thus pose an active challenge to obtaining reliable inferences of relationships and divergence times near the root of placental mammals, for example. An open source implementation is available under the LGPL 3.0 license in the beastbooster package for BEAST 2, available from https://github.com/rbouckaert/beastbooster .

Association analysis of mitochondrial DNA heteroplasmic variants: Methods and application

We rigorously assessed a comprehensive association testing framework for heteroplasmy, employing both simulated and real-world data. This framework employed a variant allele fraction (VAF) threshold and harnessed multiple gene-based tests for robust identification and association testing of heteroplasmy. Our simulation studies demonstrated that gene-based tests maintained an appropriate type I error rate at α = 0.001. Notably, when 5 % or more heteroplasmic variants within a target region were linked to an outcome, burden-extension tests (including the adaptive burden test, variable threshold burden test, and z-score weighting burden test) outperformed the sequence kernel association test (SKAT) and the original burden test. Applying this framework, we conducted association analyses on whole-blood derived heteroplasmy in 17,507 individuals of African and European ancestries (31 % of African Ancestry, mean age of 62, with 58 % women) with whole genome sequencing data. We performed both cohort- and ancestry-specific association analyses, followed by meta-analysis on both pooled samples and within each ancestry group. Our results suggest that mtDNA-encoded genes/regions are likely to exhibit varying rates in somatic aging, with the notably strong associations observed between heteroplasmy in the RNR1 and RNR2 genes (p < 0.001) and advance aging by the Original Burden test. In contrast, SKAT identified significant associations (p < 0.001) between diabetes and the aggregated effects of heteroplasmy in several protein-coding genes. Further research is warranted to validate these findings. In summary, our proposed statistical framework represents a valuable tool for facilitating association testing of heteroplasmy with disease traits in large human populations.

Molecular species identification from seized dried genital organs: a case study of suspected Bengal monitor lizard (Varanus bengalensis)

The bifurcated hemipenes of any of the four species of Monitor lizards found in India is illegally traded under the name ‘Hatha Jodi’ (a plant root from Martynia annua). A rare plant root is misrepresented as a powerful charm capable of bringing property and good fortune to its possessor. As a result, there is a suspected large-scale poaching of monitors in India to fuel the trade of ‘Hatha Jodi’. With this regard, here we report a case where the Tamil Nadu Forest Department officers seized dried genital organs. We extracted DNA from the seized dried genital organs and amplified a partial fragment of the mitochondrial cytochrome-b gene (cytb). Then, generated DNA sequences were analyzed and confirmed that the dried genital organs were consistent with the Bengal monitor lizard with 99% similarity. In forensic case analysis, dried and processed samples are not preferred for DNA-based analysis. This is due to the low yield of DNA and increased PCR inhibitors, which lead to low success in such experiments. However, the present study demonstrated species identification through DNA analysis from seized dried genital organs and supported law enforcement in successfully prosecuting the case.

Mitochondrial Dysfunction in Aristolochic Acid I-Induced Kidney Diseases: What We Know and What We Do Not Know

Aristolochic acids, compounds derived from Aristolochiaceae plant species, are associated with significant renal nephrotoxicity and carcinogenicity. Aristolochic acid I (AAI), the most predominant and potent of these compounds, is a primary etiological agent in acute and chronic kidney diseases such as Aristolochic Acid Nephropathy (AAN) and Balkan Endemic Nephropathy (BEN). Due to the kidneys’ critical role in xenobiotic excretion, they are the primary organs affected by AAI toxicity. Recent in vitro and in vivo studies have highlighted mitochondrial dysfunction as a crucial factor in the pathogenesis of these kidney diseases. This review provides an update on the recent advances in understanding the causes of acquired mitochondrial dysfunction within the context of AAN and BEN. Key findings include the identification of mitochondrial DNA depletion, loss of mitochondrial membrane potential, and decreased ATP production as significant contributors to kidney damage. Additionally, oxidative stress markers and inflammatory mediators have been implicated in disease progression. Potential therapeutic approaches, such as the use of antioxidants like vitamin C and catalpol, have shown promise in mitigating AAI-induced cytotoxicity. Furthermore, future predictive approaches like pharmacogenomics could pave the way for novel mitochondria-targeted treatments. A comprehensive characterization of mitochondrial function, its underlying molecular mechanisms, and specific biomarkers could offer valuable insights and potential therapeutic options, significantly impacting the current management of AAN and BEN.

The dysregulation of mitochondrial homeostasis-related genes could be involved in the decrease of mtDNA copy number in systemic lupus erythematosus patients.

Systemic lupus erythematosus (SLE) is a chronic multifactorial autoimmune disease. It is now widely demonstrated that oxidative stress (OS) plays an important role in the modulation of the pathogenesis of this disease. Mitochondrial DNA (mtDNA) is highly vulnerable to OS and it is known a decrease of mtDNA copy number in SLE patients. However, to date, it has not been investigated if this decrease is associated with a dysregulation of mitochondrial homeostasis genes. Our aim is to evaluate the amount of mtDNA copy number and the expression of the genes more involved in the mitochondrial homeostasis pathways, in peripheral blood mononuclear cells (PBMCs) of SLE patients and healthy controls. We analysed the amount of mtDNA in PBMCs of 72 SLE patients and 61 healthy controls by qPCR. Then, we investigated the expression variability of TFAM and SIRT1 (biogenesis), MFN1 and MFF (fusion/fission) and PRKN2 (mitophagy) genes in a subgroup of SLE patients and healthy controls. Interestingly, we have observed a highly significant decrease in mtDNA copies in SLE patients compared to healthy controls (P < 0.0001). In addition, we have shown that the expression levels of SIRT1, MFN1 and PRKN2 genes were significantly decreased in SLE patients with respect to healthy controls (P = 0.00001 for SIRT1, P = 0.0150 for MFN1 and P = 0.0009 for PRKN2). Lastly, we have reported a positive correlation between PRKN2 expression level and mtDNA copy number (P = 0.019, r = 0.475). In conclusion, our data confirm the impairment of mtDNA copy number in the disease and show for the first time a dysregulation of the mitochondrial homeostasis genes. These results could provide additional support to the important role of mitochondria in SLE development.

Determining Gene Order Patterns in the Suillus and Boletales through Comparative Analysis of Their Mitogenomes.

Suillus is one of the most important genera of ectomycorrhizal fungi. As a model for studying host specificity, its molecular fragments and nuclear genome have been analyzed. However, its mitochondrial genome has not yet been reported. In this study, we assembled five mitogenomes of Suillus and analyzed and compared their basic characteristics. Owing to the large number of introns as well as intergenic regions, the mitogenomic lengths of species of Suillus were greater than those of other species of Boletales. We identified two main patterns of gene order arrangement in the members of the order Boletales. The Ka/Ks values of 15 protein-coding genes were <1 for the mitochondrial genes of 39 Boletales species, indicating their conserved evolution. Phylogenetic trees, reconstructed using the mitogenomes, indicated that the genus Suillus was monophyletic. Phylogenetic results based on the internal transcribed spacer region and mitogenome were used to confirm the distribution of Suillus placidus in China. The results showed that the mitogenome was superior in distinguishing species compared with a single molecular fragment. This is the first study to investigate the mitogenome of Suillus, enriching the mitogenome information and providing basic data for the phylogeny, resource conservation, and genetic diversity of this genus.