Insect mitochondrial genomes (mitogenomes) are usually represented by a conserved gene order. Whiteflies exhibit gene rearrangement in their mitogenomes; however, understanding how nucleotide substitution rates shape gene rearrangement in whiteflies is unclear due to the limited number of mitogenomes. Additionally, the mechanisms by which selection pressure drives adaptations in mitochondrial genes in the two subfamilies of whiteflies are not yet known. Here, we analyzed 18 whitefly mitogenomes, including one newly generated mitogenome, to compare nucleotide substitution rates, selection pressure, and gene arrangements. The newly generated mitogenome is reported along with reannotation of Pealius mori and comparisons to other whitefly mitogenomes. Comparative studies on nucleotide composition of 18 whiteflies revealed the positive GC skewness, confirming the reversal of strand asymmetry. We found 11 rearranged gene orders within two subfamilies of whiteflies with 8-18 breakpoints of gene rearrangements. Members of the subfamily Aleyrodinae exhibit more complex pathways in the evolution of gene order as compared to the subfamily Aleurodicinae. Our findings also revealed that the increase or reduction of nucleotide substitution rates does not have an impact on any of the gene rearrangement scenarios depicting neutral correlation. Selection pressure analysis revealed that the mitogenomes from members of both the subfamilies Aleurodicinae and Aleyrodinae are characterized by intense purifying selection pressure.

BACKGROUND: Mitochondrial (MT) dysfunction is a hallmark of liver diseases. However, the effects of functional variants such as protein truncating variants (PTVs) in MT-related genes on the risk of liver diseases have not been extensively explored.
METHODS: We extracted 60,928 PTVs across 2466 MT-related nucleus genes using whole-exome sequencing data obtained from 442,603 participants in the UK Biobank. We examined their associations with liver dysfunction that represented by the liver-related biomarkers and the risks of chronic liver diseases and liver-related mortality.
RESULTS: 96.10% of the total participants carried at least one PTV. We identified 866 PTVs that were positively associated with liver dysfunction at the threshold of P value < 8.21e - 07. The coding genes of these PTVs were mainly enriched in pathways related to lipid, fatty acid, amino acid, and carbohydrate metabolisms. The 866 PTVs were presented in 1.07% (4721) of participants. Compared with participants who did not carry any of the PTVs, the carriers had a 5.33-fold (95% CI 4.15-6.85), 2.82-fold (1.69-4.72), and 4.41-fold (3.04-6.41) increased risk for fibrosis and cirrhosis of liver, liver cancer, and liver disease-related mortality, respectively. These adverse effects were consistent across subgroups based on age, sex, body mass index, smoking status, and presence of hypertension, diabetes, dyslipidemia, and metabolic syndrome.
CONCLUSIONS: Our findings revealed a significant impact of PTVs in MT-related genes on liver disease risk, highlighting the importance of these variants in identifying populations at risk of liver diseases and facilitating early clinical interventions.

Physella acuta is a freshwater snail native to North America. Understanding the phylogeography and genetic structure of P. acuta will help elucidate its evolution. In this study, we used mitochondrial (COI and 16S rDNA) and nuclear (ITS1) markers to identify the species and examine its genetic diversity, population structure, and demographic history of P. acuta in Thailand. Phylogenetic and network analyses of P. acuta in Thailand pertained to clade A, which exhibits a global distribution. Analysis of the genetic structure of the population revealed that the majority of pairwise comparisons showed no genetic dissimilarity. An isolation-by-distance test indicates no significant correlation between genetic and geographical distances among P. acuta populations, suggesting that gene flow is not restricted by distance. Demographic history and haplotype network analyses suggest a population expansion of P. acuta, as evidenced by the star-like structure detected in the median-joining network. Based on these results, we concluded that P. acuta in Thailand showed gene flow and recent population expansion. Our findings provide fundamental insights into the genetic variation of P. acuta in Thailand.

Locomotor preferences and habitat types may drive animal evolution. In this study, we speculated that locomotor preference and habitat type may have diverse influences on Bovidae mitochondrial genes. We used selection pressure and statistical analysis to explore the evolution of mitochondrial DNA (mtDNA) protein-coding genes (PCGs) from diverse locomotor preferences and habitat types. Our study demonstrates that locomotor preference (energy demand) drives the evolution of Bovidae in mtDNA PCGs. The habitat types had no significant effect on the rate of evolution in Bovidae mitochondrial genes. Our study provides deep insight into the adaptation of Bovidae.

Background: Colon adenocarcinoma (COAD) has increasing incidence and is one of the most common malignant tumors. The mitochondria involved in cell energy metabolism, oxygen free radical generation, and cell apoptosis play important roles in tumorigenesis and progression. The relationship between mitochondrial genes and COAD remains largely unknown. Methods: COAD data including 512 samples were set out from the UCSC Xena database. The nuclear mitochondrial-related genes (NMRGs)-related risk prognostic model and prognostic nomogram were constructed, and NMRGs-related gene mutation and the immune environment were analyzed using bioinformatics methods. Then, a liver metastasis model of colorectal cancer was constructed and protein expression was detected using Western blot assay. Results: A prognostic model for COAD was constructed. Comparing the prognostic model dataset and the validation dataset showed considerable correlation in both risk grouping and prognosis. Based on the risk score (RS) model, the samples of the prognostic dataset were divided into high risk group and low risk group. Moreover, pathologic N and T stage and tumor recurrence in the two risk groups were significantly different. The four prognostic factors, including age and pathologic T stage in the nomogram survival model also showed excellent predictive performance. An optimal combination of nine differentially expressed NMRGs was finally obtained, including LARS2, PARS2, ETHE1, LRPPRC, TMEM70, AARS2, ACAD9, VARS2, and ATP8A2. The high-RS group had more inflamed immune features, including T and CD4+ memory cell activation. Besides, mitochondria-associated LRPPRC and LARS2 expression levels were increased in vivo xenograft construction and liver metastases assays. Conclusion: This study established a comprehensive prognostic model for COAD, incorporating nine genes associated with nuclear-mitochondrial functions. This model demonstrates superior predictive performance across four prognostic factors: age, pathological T stage, tumor recurrence, and overall prognosis. It is anticipated to be an effective model for enhancing the prognosis and treatment of COAD.

Numerous studies have highlighted the pivotal role of mitochondria-related genes (MRGs) in the initiation and progression of glioblastoma (GBM). However, the specific contributions of MRGs coding proteins to GBM pathology remain incompletely elucidated. The identification of prognostic MRGs in GBM holds promise for the development of personalized targeted therapies and the enhancement of patient prognosis. We combined differential expression with univariate Cox regression analysis to screen prognosis-associated MRGs in GBM. Based on the nine MRGs, the hazard ratio model was conducted using a multivariate Cox regression algorithm. SHC-related survival, pathway, and immune analyses in GBM cohorts were obtained from the Biomarker Exploration of the Solid Tumor database. The proliferation and migration of U87 cells were measured by CCK-8 and transwell assay. Apoptosis in U87 cells was evaluated using flow cytometry. Confocal microscopy was employed to measure mitochondrial reactive oxygen species (ROS) levels and morphology. The expression levels of SHC1 and other relevant proteins were examined via western blotting. We screened 15 prognosis-associated MRGs and constructed a 9 MRGs-based model. Validation of the model\'s risk score confirmed its efficacy in predicting the prognosis of patients with GBM. Furthermore, analysis revealed that SHC1, a constituent MRG of the prognostic model, was upregulated and implicated in the progression, migration, and immune infiltration of GBM. In vitro experiments elucidated that p66Shc, the longest isoform of SHC1, modulates mitochondrial ROS production and morphology, consequently promoting the proliferation and migration of U87 cells. The 9 MRGs-based prognostic model could predict the prognosis of GBM. SHC1 was upregulated and correlated with the prognosis of patients by involvement in immune infiltration. Furthermore, in vitro experiments demonstrated that p66Shc promotes U87 cell proliferation and migration by mediating mitochondrial ROS production. Thus, p66Shc may serve as a promising biomarker and therapeutic target for GBM.

OBJECTIVE: The purpose of this study was to identify and analyze the mitochondrial genes associated with sepsis patients in order to elucidate the underlying mechanism of sepsis immunity and provide new ideas for the clinical treatment of sepsis.
METHODS: The hospitalized cases of sepsis (n = 20) and systemic inflammatory response syndrome (SIRS) (n = 12) admitted to the Emergency Intensive Care Unit (EICU) of the Affiliated Hospital of Southwest Medical University from January 2019 to December 2019 were collected consecutively. RNA-seq was used to sequence the RNA (mRNA) of peripheral blood cells. Bioinformatics techniques were used to screen and identify differentially expressed RNAs, with an absolute value of fold change (FC) greater than or equal to 1.2 and a false discovery rate (FDR) less than 0.05. At the same time, mitochondrial genes were obtained from the MitoCarta 3.0 database. Differential genes were then intersected with mitochondrial genes. The resulting crossover genes were subjected to GO, KEGG, and PPI analysis. Subsequently, the GSE65682 dataset was downloaded from the GEO database for survival analysis to assess the prognostic value of core genes, and GSE67652 was downloaded for ROC curve analysis to validate the diagnostic value of core genes. Finally, the localization of core genes was clarified through 10X single-cell sequencing.
RESULTS: The crossing of 314 sepsis differential genes and 1136 mitochondrial genes yielded 28 genes. GO and KEGG analysis showed that the crossover genes were mainly involved in the mitochondrion, mitochondrial matrix, and mitochondrial inner membrane. Survival analysis screened four genes that were significantly negatively associated with the prognosis of sepsis, namely FIS1, FKBP8, GLRX5, and GUK1. A comparison of peripheral blood RNA-seq results between the sepsis group and the SIRS group showed that the expression levels of these four genes were significantly decreased in the sepsis group compared to the SIRS group. ROC curve analysis based on GSE67652 indicates these four genes\' high sensitivity and specificity for sepsis detection. Additionally, single-cell RNA sequencing found that the core genes were mainly expressed in macrophages, T cells, and B cells.
CONCLUSIONS: Mitochondria-related genes (FIS1, FKBP8, GLRX5, GUK1) were underexpressed in the sepsis group, negatively correlated with survival, and mainly distributed in immune cells. This finding may guide studying the immune-related mechanisms of sepsis. This study protocol was reviewed by the Ethics Committee of the Affiliated Hospital of Southwest Medical University (ethics number: KY2018029), the clinical trial registration number is ChiCTR1900021261, and the registration date is February 4, 2019.

BACKGROUND: Maternally inherited hearing loss has been associated with mitochondrial genes, including MT-RNR1, MT-TL1, MT-TS1, MT-TK and MT-TE. Among these genes, MT-RNR1 is known to be a hotspot for pathogenic variants related to aminoglycoside ototoxicity and nonsyndromic hearing loss. However, the frequency and spectrum of variants in these genes, particularly in multi-ethnic hearing loss patients from Southwestern China, are still not fully understood.
METHODS: In this study, we enrolled 460 hearing loss patients from various ethnic backgrounds (Han, Yi, Dai, Hani, etc.) in Southwestern China. Next-generation sequencing was used to analyze the mitochondrial MT-RNR1, MT-TL1, MT-TS1, MT-TK and MT-TE genes. Subsequently, bioinformatical methods were employed to evaluate the identified variants.
RESULTS: Among the patients with hearing loss, we identified 70 variants in MT-RNR1 (78.6 %, 55/70), MT-TL1 (4.3 %, 3/70), MT-TS1 (4.3 %, 3/70), MT-TK (7.1 %, 5/70) and MT-TE (5.7 %, 4/70) genes. We found that 15 variants were associated with hearing loss, including m.1555 A > G and m.1095 T > C. Additionally, we discovered three reported mitochondrial variants (m.676 G > A, m.7465 insC, and m.7474 A > G) newly correlated with hearing loss. Notably, certain pathogenic variants, such as m.1555 A > G, displayed non-consistent distributions among the multi-ethnic patients with hearing loss. Furthermore, the number of variants associated with hearing loss was higher in the Sinitic group (n = 181) and Tibeto-Burman group (n = 215) compared to the Kra-Dai group (n = 38) and Hmong-Mien group (n = 26).
CONCLUSIONS: This present study revealed the distribution of mitochondrial variants linked to hearing loss across various ethnic groups in Southwestern China. These data suggest a potential correlation between the distribution of mitochondrial variants associated with hearing loss and ethnic genetic backgrounds.

BACKGROUND: Evidence suggests that hepatocyte mitochondrial dysfunction leads to abnormal lipid metabolism, redox imbalance, and programmed cell death, driving the onset and progression of non-alcoholic steatohepatitis (NASH). Identifying hub mitochondrial genes linked to NASH may unveil potential therapeutic targets.
METHODS: Mitochondrial hub genes implicated in NASH were identified via analysis using 134 algorithms.
RESULTS: The Random Forest algorithm (RF), the most effective among the 134 algorithms, identified three genes: Aldo-keto reductase family 1 member B10 (AKR1B10), thymidylate synthase (TYMS), and triggering receptor expressed in myeloid cell 2 (TREM2). They were upregulated and positively associated with genes promoting inflammation, genes involved in lipid synthesis, fibrosis, and nonalcoholic steatohepatitis activity scores in patients with NASH. Moreover, using these three genes, patients with NASH were accurately categorized into cluster 1, exhibiting heightened disease severity, and cluster 2, distinguished by milder disease activity.
CONCLUSIONS: These three genes are pivotal mitochondrial genes implicated in NASH progression.

Although genome-scale data generation is becoming more tractable for phylogenetics, there are large quantities of single gene fragment data in public repositories and such data are still being generated. We therefore investigated whether single mitochondrial genes are suitable proxies for phylogenetic reconstruction as compared to the application of full mitogenomes. With near complete taxon sampling for the southern African dwarf chameleons (Bradypodion), we estimated and compared phylogenies for the complete mitogenome with topologies generated from individual mitochondrial genes and various combinations of these genes. Our results show that the topologies produced by single genes (ND2, ND4, ND5, COI, and COIII) were analogous to the complete mitogenome, suggesting that these genes may be reliable markers for generating mitochondrial phylogenies in lieu of generating entire mitogenomes. In contrast, the short fragment of 16S commonly used in herpetological systematics, produced a topology quite dissimilar to the complete mitogenome and its concatenation with ND2 weakened the resolution of ND2. We therefore recommend the avoidance of this 16S fragment in future phylogenetic work.