Pervasive structural heterogeneity rewires glioblastoma chromosomes to sustain patient-specific transcriptional programs.

Glioblastoma multiforme (GBM) encompasses brain malignancies marked by phenotypic and transcriptional heterogeneity thought to render these tumors aggressive, resistant to therapy, and inevitably recurrent. However, little is known about how the spatial organization of GBM genomes underlies this heterogeneity and its effects. Here, we compile a cohort of 28 patient-derived glioblastoma stem cell-like lines (GSCs) known to reflect the properties of their tumor-of-origin; six of these were primary-relapse tumor pairs from the same patient. We generate and analyze 5 kbp-resolution chromosome conformation capture (Hi-C) data from all GSCs to systematically map thousands of standalone and complex structural variants (SVs) and the multitude of neoloops arising as a result. By combining Hi-C, histone modification, and gene expression data with chromatin folding simulations, we explain how the pervasive, uneven, and idiosyncratic occurrence of neoloops sustains tumor-specific transcriptional programs via the formation of new enhancer-promoter contacts. We also show how even moderately recurrent neoloops can relate to patient-specific vulnerabilities. Together, our data provide a resource for dissecting GBM biology and heterogeneity, as well as for informing therapeutic approaches.

Heterogeneity and molecular landscape of melanoma: implications for targeted therapy.

Uveal cancer (UM) offers a complex molecular landscape characterized by substantial heterogeneity, both on the genetic and epigenetic levels. This heterogeneity plays a critical position in shaping the behavior and response to therapy for this uncommon ocular malignancy. Targeted treatments with gene-specific therapeutic molecules may prove useful in overcoming radiation resistance, however, the diverse molecular makeups of UM call for a patient-specific approach in therapy procedures. We need to understand the intricate molecular landscape of UM to develop targeted treatments customized to each patient's specific genetic mutations. One of the promising approaches is using liquid biopsies, such as circulating tumor cells (CTCs) and circulating tumor DNA (ctDNA), for detecting and monitoring the disease at the early stages. These non-invasive methods can help us identify the most effective treatment strategies for each patient. Single-cellular is a brand-new analysis platform that gives treasured insights into diagnosis, prognosis, and remedy. The incorporation of this data with known clinical and genomics information will give a better understanding of the complicated molecular mechanisms that UM diseases exploit. In this review, we focused on the heterogeneity and molecular panorama of UM, and to achieve this goal, the authors conducted an exhaustive literature evaluation spanning 1998 to 2023, using keywords like "uveal melanoma, "heterogeneity". "Targeted therapies"," "CTCs," and "single-cellular analysis".

Interpretation of 10 years of Alzheimer's disease genetic findings in the perspective of statistical heterogeneity.

Common genetic variants and susceptibility loci associated with Alzheimer's disease (AD) have been discovered through large-scale genome-wide association studies (GWAS), GWAS by proxy (GWAX) and meta-analysis of GWAS and GWAX (GWAS+GWAX). However, due to the very low repeatability of AD susceptibility loci and the low heritability of AD, these AD genetic findings have been questioned. We summarize AD genetic findings from the past 10 years and provide a new interpretation of these findings in the context of statistical heterogeneity. We discovered that only 17% of AD risk loci demonstrated reproducibility with a genome-wide significance of P < 5.00E-08 across all AD GWAS and GWAS+GWAX datasets. We highlighted that the AD GWAS+GWAX with the largest sample size failed to identify the most significant signals, the maximum number of genome-wide significant genetic variants or maximum heritability. Additionally, we identified widespread statistical heterogeneity in AD GWAS+GWAX datasets, but not in AD GWAS datasets. We consider that statistical heterogeneity may have attenuated the statistical power in AD GWAS+GWAX and may contribute to explaining the low repeatability (17%) of genome-wide significant AD susceptibility loci and the decreased AD heritability (40-2%) as the sample size increased. Importantly, evidence supports the idea that a decrease in statistical heterogeneity facilitates the identification of genome-wide significant genetic loci and contributes to an increase in AD heritability. Collectively, current AD GWAX and GWAS+GWAX findings should be meticulously assessed and warrant additional investigation, and AD GWAS+GWAX should employ multiple meta-analysis methods, such as random-effects inverse variance-weighted meta-analysis, which is designed specifically for statistical heterogeneity.

Integrating spatial and single-cell transcriptomics reveals tumor heterogeneity and intercellular networks in colorectal cancer.

Single cell RNA sequencing (scRNA-seq), a powerful tool for studying the tumor microenvironment (TME), does not preserve/provide spatial information on tissue morphology and cellular interactions. To understand the crosstalk between diverse cellular components in proximity in the TME, we performed scRNA-seq coupled with spatial transcriptomic (ST) assay to profile 41,700 cells from three colorectal cancer (CRC) tumor-normal-blood pairs. Standalone scRNA-seq analyses revealed eight major cell populations, including B cells, T cells, Monocytes, NK cells, Epithelial cells, Fibroblasts, Mast cells, Endothelial cells. After the identification of malignant cells from epithelial cells, we observed seven subtypes of malignant cells that reflect heterogeneous status in tumor, including tumor_CAV1, tumor_ATF3_JUN | FOS, tumor_ZEB2, tumor_VIM, tumor_WSB1, tumor_LXN, and tumor_PGM1. By transferring the cellular annotations obtained by scRNA-seq to ST spots, we annotated four regions in a cryosection from CRC patients, including tumor, stroma, immune infiltration, and colon epithelium regions. Furthermore, we observed intensive intercellular interactions between stroma and tumor regions which were extremely proximal in the cryosection. In particular, one pair of ligands and receptors (C5AR1 and RPS19) was inferred to play key roles in the crosstalk of stroma and tumor regions. For the tumor region, a typical feature of TMSB4X-high expression was identified, which could be a potential marker of CRC. The stroma region was found to be characterized by VIM-high expression, suggesting it fostered a stromal niche in the TME. Collectively, single cell and spatial analysis in our study reveal the tumor heterogeneity and molecular interactions in CRC TME, which provides insights into the mechanisms underlying CRC progression and may contribute to the development of anticancer therapies targeting on non-tumor components, such as the extracellular matrix (ECM) in CRC. The typical genes we identified may facilitate to new molecular subtypes of CRC.

Polygenic Risk Score Reveals Genetic Heterogeneity of Alzheimer's Disease between the Chinese and European Populations.

BACKGROUND: The polygenic risk score (PRS) aggregates the effects of numerous genetic variants associated with a condition across the human genome and may help to predict late-onset Alzheimer's disease (LOAD). Most of the current PRS studies on Alzheimer's disease (AD) have been conducted in Caucasian ancestry populations, while it is less studied in Chinese. OBJECTIVE: To establish and examine the validity of Chinese PRS, and explore its racial heterogeneity. DESIGN: We constructed a PRS using both discovery (N = 2012) and independent validation samples (N = 1008) from Chinese population. The associations between PRS and age at onset of LOAD or cerebrospinal fluid (CSF) biomarkers were assessed. We also replicated the PRS in an independent replication cohort with CSF data and constructed an alternative PRS using European weights. SETTING: Multi-center genetics study. PARTICIPANTS: A total of 3020 subjects were included in the study. MEASUREMENTS: PRS was calculated using genome-wide association studies data and evaluated the performance alone (PRSnoAPOE) and with other predictors (full model: LOAD ~ PRSnoAPOE + APOE+ sex + age) by measuring the area under the receiver operating curve (AUC). RESULTS: PRS of the full model achieved the highest AUC of 84.0% (95% CI = 81.4-86.5) with pT< 0.5, compared with the model containing APOE alone (61.0%). The AUC of PRS with pT<5e-8 was 77.8% in the PRSnoAPOE model, 81.5% in the full model, and only ranged from 67.5% to 75.1% in the PRS with the European weights model. A higher PRS was significantly associated with an earlier age at onset (P <0.001). The PRS also performed well in the replication cohort of the full model (AUC=83.1%, 95% CI = 74.3-92.0). The CSF biomarkers of Aß42 and the ratio of Aß42/Aß40 were significantly inversely associated with the PRS, while p-Tau181 showed a positive association. CONCLUSIONS: This finding suggests that PRS reveal genetic heterogeneity and higher prediction accuracy of the PRS for AD can be achieved using a base dataset and validation within the same ethnicity. The effective PRS model has the clinical potential to predict individuals at risk of developing LOAD at a given age and with abnormal levels of CSF biomarkers in the Chinese population.

Understanding the genomic heterogeneity of North African Imazighen: from broad to microgeographical perspectives.

The strategic location of North Africa has led to cultural and demographic shifts, shaping its genetic structure. Historical migrations brought different genetic components that are evident in present-day North African genomes, along with autochthonous components. The Imazighen (plural of Amazigh) are believed to be the descendants of autochthonous North Africans and speak various Amazigh languages, which belong to the Afro-Asiatic language family. However, the arrival of different human groups, especially during the Arab conquest, caused cultural and linguistic changes in local populations, increasing their heterogeneity. We aim to characterize the genetic structure of the region, using the largest Amazigh dataset to date and other reference samples. Our findings indicate microgeographical genetic heterogeneity among Amazigh populations, modeled by various admixture waves and different effective population sizes. A first admixture wave is detected group-wide around the twelfth century, whereas a second wave appears in some Amazigh groups around the nineteenth century. These events involved populations with higher genetic ancestry from south of the Sahara compared to the current North Africans. A plausible explanation would be the historical trans-Saharan slave trade, which lasted from the Roman times to the nineteenth century. Furthermore, our investigation shows that assortative mating in North Africa has been rare.

Unveiling the Dynamics behind Glioblastoma Multiforme Single-Cell Data Heterogeneity.

Glioblastoma Multiforme is a brain tumor distinguished by its aggressiveness. We suggested that this aggressiveness leads single-cell RNA-sequence data (scRNA-seq) to span a representative portion of the cancer attractors domain. This conjecture allowed us to interpret the scRNA-seq heterogeneity as reflecting a representative trajectory within the attractor's domain. We considered factors such as genomic instability to characterize the cancer dynamics through stochastic fixed points. The fixed points were derived from centroids obtained through various clustering methods to verify our method sensitivity. This methodological foundation is based upon sample and time average equivalence, assigning an interpretative value to the data cluster centroids and supporting parameters estimation. We used stochastic simulations to reproduce the dynamics, and our results showed an alignment between experimental and simulated dataset centroids. We also computed the Waddington landscape, which provided a visual framework for validating the centroids and standard deviations as characterizations of cancer attractors. Additionally, we examined the stability and transitions between attractors and revealed a potential interplay between subtypes. These transitions might be related to cancer recurrence and progression, connecting the molecular mechanisms of cancer heterogeneity with statistical properties of gene expression dynamics. Our work advances the modeling of gene expression dynamics and paves the way for personalized therapeutic interventions.

Mutation rate heterogeneity at the sub-gene scale due to local DNA hypomethylation.

Local mutation rates in human are highly heterogeneous, with known variability at the scale of megabase-sized chromosomal domains, and, on the other extreme, at the scale of oligonucleotides. The intermediate, kilobase-scale heterogeneity in mutation risk is less well characterized. Here, by analyzing thousands of somatic genomes, we studied mutation risk gradients along gene bodies, representing a genomic scale spanning roughly 1-10 kb, hypothesizing that different mutational mechanisms are differently distributed across gene segments. The main heterogeneity concerns several kilobases at the transcription start site and further downstream into 5' ends of gene bodies; these are commonly hypomutated with several mutational signatures, most prominently the ubiquitous C > T changes at CpG dinucleotides. The width and shape of this mutational coldspot at 5' gene ends is variable across genes, and corresponds to variable interval of lowered DNA methylation depending on gene activity level and regulation. Such hypomutated loci, at 5' gene ends or elsewhere, correspond to DNA hypomethylation that can associate with various landmarks, including intragenic enhancers, Polycomb-marked regions, or chromatin loop anchor points. Tissue-specific DNA hypomethylation begets tissue-specific local hypomutation. Of note, direction of mutation risk is inverted for AID/APOBEC3 cytosine deaminase activity, whose signatures are enriched in hypomethylated regions.

Patterns of Chromosomal Instability and Clonal Heterogeneity in Luminal B Breast Cancer: A Pilot Study.

Chromosomal instability (CIN), defined by variations in the number or structure of chromosomes from cell to cell, is recognized as a distinctive characteristic of cancer associated with the ability of tumors to adapt to challenging environments. CIN has been recognized as a source of genetic variation that leads to clonal heterogeneity (CH). Recent findings suggest a potential association between CIN and CH with the prognosis of BC patients, particularly in tumors expressing the epidermal growth factor receptor 2 (HER2+). In fact, information on the role of CIN in other BC subtypes, including luminal B BC, is limited. Additionally, it remains unknown whether CIN in luminal B BC tumors, above a specific threshold, could have a detrimental effect on the growth of human tumors or whether low or intermediate CIN levels could be linked to a more favorable BC patient prognosis when contrasted with elevated levels. Clarifying these relationships could have a substantial impact on risk stratification and the development of future therapeutic strategies aimed at targeting CIN in BC. This study aimed to assess CIN and CH in tumor tissue samples from ten patients with luminal B BC and compare them with established clinicopathological parameters. The results of this study reveal that luminal B BC patients exhibit intermediate CIN and stable aneuploidy, both of which correlate with lymphovascular invasion. Our results also provide valuable preliminary data that could contribute to the understanding of the implications of CIN and CH in risk stratification and the development of future therapeutic strategies in BC.

A Multiplexed Approach to Assess Small Cell Lung Cancer Subtype Heterogeneity in Primary and Patient-Derived Tumor Samples.

Unlocking the heterogeneity of cancers is crucial for developing therapeutic approaches that effectively eradicate disease. As our understanding of markers specific to cancer subclones or subtypes expands, there is a growing demand for advanced technologies that enable the simultaneous investigation of multiple targets within an individual tumor sample. Indeed, multiplex approaches offer distinct benefits, particularly when tumor specimens are small and scarce. Here we describe the utility of two fluorescence-based multiplex approaches; fluorescent Western blots, and multiplex immunohistochemistry (Opal™) staining to interrogate heterogeneity, using small cell lung cancer as an example. Critically, the coupling of Opal™ staining with advanced image quantitation, permits the dissection of cancer cell phenotypes at a single cell level. These approaches can be applied to patient biopsies and/or patient-derived xenograft (PDX) models and serve as powerful methodologies for assessing tumor cell heterogeneity in response to therapy or between metastatic lesions across diverse tissue sites.

Quantifying intra-tumoral genetic heterogeneity of glioblastoma toward precision medicine using MRI and a data-inclusive machine learning algorithm.

BACKGROUND AND OBJECTIVE: Glioblastoma (GBM) is one of the most aggressive and lethal human cancers. Intra-tumoral genetic heterogeneity poses a significant challenge for treatment. Biopsy is invasive, which motivates the development of non-invasive, MRI-based machine learning (ML) models to quantify intra-tumoral genetic heterogeneity for each patient. This capability holds great promise for enabling better therapeutic selection to improve patient outcome. METHODS: We proposed a novel Weakly Supervised Ordinal Support Vector Machine (WSO-SVM) to predict regional genetic alteration status within each GBM tumor using MRI. WSO-SVM was applied to a unique dataset of 318 image-localized biopsies with spatially matched multiparametric MRI from 74 GBM patients. The model was trained to predict the regional genetic alteration of three GBM driver genes (EGFR, PDGFRA and PTEN) based on features extracted from the corresponding region of five MRI contrast images. For comparison, a variety of existing ML algorithms were also applied. Classification accuracy of each gene were compared between the different algorithms. The SHapley Additive exPlanations (SHAP) method was further applied to compute contribution scores of different contrast images. Finally, the trained WSO-SVM was used to generate prediction maps within the tumoral area of each patient to help visualize the intra-tumoral genetic heterogeneity. RESULTS: WSO-SVM achieved 0.80 accuracy, 0.79 sensitivity, and 0.81 specificity for classifying EGFR; 0.71 accuracy, 0.70 sensitivity, and 0.72 specificity for classifying PDGFRA; 0.80 accuracy, 0.78 sensitivity, and 0.83 specificity for classifying PTEN; these results significantly outperformed the existing ML algorithms. Using SHAP, we found that the relative contributions of the five contrast images differ between genes, which are consistent with findings in the literature. The prediction maps revealed extensive intra-tumoral region-to-region heterogeneity within each individual tumor in terms of the alteration status of the three genes. CONCLUSIONS: This study demonstrated the feasibility of using MRI and WSO-SVM to enable non-invasive prediction of intra-tumoral regional genetic alteration for each GBM patient, which can inform future adaptive therapies for individualized oncology.

Better together against genetic heterogeneity: A sex-combined joint main and interaction analysis of 290 quantitative traits in the UK Biobank.

Genetic effects can be sex-specific, particularly for traits such as testosterone, a sex hormone. While sex-stratified analysis provides easily interpretable sex-specific effect size estimates, the presence of sex-differences in SNP effect implies a SNP×sex interaction. This suggests the usage of the often overlooked joint test, testing for an SNP's main and SNP×sex interaction effects simultaneously. Notably, even without individual-level data, the joint test statistic can be derived from sex-stratified summary statistics through an omnibus meta-analysis. Utilizing the available sex-stratified summary statistics of the UK Biobank, we performed such omnibus meta-analyses for 290 quantitative traits. Results revealed that this approach is robust to genetic effect heterogeneity and can outperform the traditional sex-stratified or sex-combined main effect-only tests. Therefore, we advocate using the omnibus meta-analysis that captures both the main and interaction effects. Subsequent sex-stratified analysis should be conducted for sex-specific effect size estimation and interpretation.

Molecular Heterogeneity in Leiomyosarcoma and Implications for Personalised Medicine.

OPINION STATEMENT: Leiomyosarcoma (LMS) is one of the more common subtypes of soft tissue sarcomas (STS), accounting for about 20% of cases. Differences in anatomical location, risk of recurrence and histomorphological variants contribute to the substantial clinical heterogeneity in survival outcomes and therapy responses observed in patients. There is therefore a need to move away from the current one-size-fits-all treatment approach towards a personalised strategy tailored for individual patients. Over the past decade, tissue profiling studies have revealed key genomic features and an additional layer of molecular heterogeneity among patients, with potential utility for optimal risk stratification and biomarker-matched therapies. Furthermore, recent studies investigating intratumour heterogeneity and tumour evolution patterns in LMS suggest some key features that may need to be taken into consideration when designing treatment strategies and clinical trials. Moving forward, national and international collaborative efforts to aggregate expertise, data, resources and tools are needed to achieve a step change in improving patient survival outcomes in this disease of unmet need.

Bacterial phenotypic heterogeneity through the lens of single-cell RNA sequencing.

Bacterial transcription is not monolithic. Microbes exist in a wide variety of cell states that help them adapt to their environment, acquire and produce essential nutrients, and engage in both competition and cooperation with their neighbors. While we typically think of bacterial adaptation as a group behavior, where all cells respond in unison, there is often a mixture of phenotypic responses within a bacterial population, where distinct cell types arise. A primary phenomenon driving these distinct cell states is transcriptional heterogeneity. Given that bacterial mRNA transcripts are extremely short-lived compared to eukaryotes, their transcriptional state is closely associated with their physiology, and thus the transcriptome of a bacterial cell acts as a snapshot of the behavior of that bacterium. Therefore, the application of single-cell transcriptomics to microbial populations will provide novel insight into cellular differentiation and bacterial ecology. In this review, we provide an overview of transcriptional heterogeneity in microbial systems, discuss the findings already provided by single-cell approaches, and plot new avenues of inquiry in transcriptional regulation, cellular biology, and mechanisms of heterogeneity that are made possible when microbial communities are analyzed at single-cell resolution.

Treatments, prognostic factors, and genetic heterogeneity in advanced cholangiocarcinoma: A multicenter real-world study.

BACKGROUND AND AIMS: Cholangiocarcinoma (CCA), a rare and aggressive hepatobiliary malignancy, presents significant clinical management challenges. Despite rising incidence and evolving treatment options, prognosis remains poor, motivating the exploration of real-world data for enhanced understanding and patient care. METHODS: This multicenter study analyzed data from 120 metastatic CCA patients at three institutions from 2016 to 2023. Kaplan-Meier curves assessed overall survival (OS), while univariate and multivariate analyses evaluated links between clinical variables (age, gender, tumor site, metastatic burden, ECOG performance status, response to first-line chemotherapy) and OS. Genetic profiling was conducted selectively. RESULTS: Enrolled patients had a median age of 68.5 years, with intrahepatic tumors predominant in 79 cases (65.8%). Among 85 patients treated with first-line chemotherapy, cisplatin and gemcitabine (41.1%) was the most common regimen. Notably, one-third received no systemic treatment. After a median 14-month follow-up, 81 CCA-related deaths occurred, with a median survival of 13.1 months. Two clinical variables independently predicted survival: response to first-line chemotherapy (disease control vs. no disease control; HR: 0.27; 95% CI: 0.14-0.50; p < 0.0001) and metastatic involvement (>1 site vs. 1 site; HR: 1.99; 95% CI: 1.04-3.80; p = 0.0366). The three most common genetic alterations involved the ARID1A, tp53, and CDKN2A genes. CONCLUSIONS: Advanced CCA displays aggressive clinical behavior, emphasizing the need for treatments beyond chemotherapy. Genetic diversity supports potential personalized therapies. Collaborative research and deeper CCA biology understanding are crucial to enhance patient outcomes in this challenging malignancy.

Statistical methods for assessing the effects of de novo variants on birth defects.

With the development of next-generation sequencing technology, de novo variants (DNVs) with deleterious effects can be identified and investigated for their effects on birth defects such as congenital heart disease (CHD). However, statistical power is still limited for such studies because of the small sample size due to the high cost of recruiting and sequencing samples and the low occurrence of DNVs. DNV analysis is further complicated by genetic heterogeneity across diseased individuals. Therefore, it is critical to jointly analyze DNVs with other types of genomic/biological information to improve statistical power to identify genes associated with birth defects. In this review, we discuss the general workflow, recent developments in statistical methods, and future directions for DNV analysis.

Patient genetic heterogeneities acting as indicators of post-operative pain and opioid requirement in orthopedic surgery: A systematic review.

INTRODUCTION: Orthopedic surgical procedures are expected to increase annually, making it imperative to understand the correlations between patient genetic makeup and post-operative pain levels. METHODS: We performed a systematic literature review using PubMed and Cochrane databases in accordance with Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines. A total of 299 articles were initially selected, 20 articles remained after title and abstract review, and nine articles were selected for inclusion upon full text review. RESULTS: Genetic risk factors identified included the A allele of the 5HT2A gene single nucleotide polymorphism, the AA genotype of the ADRB2 gene, the CG genotype of the IL6 gene, the genotypes CT and TT of the NTRK1 gene, genotypes AA and GA of the OPRM gene, and the AA and GA genotypes of the COMT gene. Additional studies in the review discuss statistical significance of other variants of the COMT gene. CONCLUSION: There have been genetic association studies performed on the patient heterogeneity and its relationship on patient pain levels, but more data need to be collected to understand the clinical utility of stratifying patients based on genomic sequence.

Genetic heterogeneity and respiratory chain enzyme analysis in pediatric Indian patients with mitochondrial disorder: Report of novel variants in POLG1 gene and their functional implication using molecular dynamic simulation.

Mitochondrial disorders are a heterogeneous group of disorders caused by mutations in the mitochondrial DNA or in nuclear genes encoding the mitochondrial proteins and subunits. Polymerase Gamma (POLG) is a nuclear gene and mutation in the POLG gene are one of the major causes of inherited mitochondrial disorders. In this study, 15 pediatric patients, with a wide spectrum of clinical phenotypes were screened using blood samples (n = 15) and muscle samples (n = 4). Respiratory chain enzyme analysis in the muscle samples revealed multi-complex deficiencies with Complex I deficiency present in (1/4) patients, Complex II (2/4), Complex III (3/4) and Complex IV (2/4) patients. Multiple large deletions were observed in 4/15 patients using LR-PCR. Whole exome sequencing (WES) revealed a compound heterozygous mutation consisting of a POLG1 novel variant (NP_002684.1:p.Trp261X) and a missense variant (NP_002684.1:p. Leu304Arg) in one patient and another patient harboring a novel homozygous POLG1 variant (NP_002684.1:p. Phe750Val). These variants (NP_002684.1:p. Leu304Arg) and (NP_002684.1:p. Phe750Val) and their interactions with DNA were modelled using molecular docking and molecular dynamics (MD) simulation studies. The protein conformation was analyzed as root mean square deviation (RMSD), root mean square fluctuation (RMSF) which showed local fluctuations in the mutants compared to the wildtype. However, Solvent Accessible Surface Area (SASA) significantly increased for NP_002684.1:p.Leu304Arg and decreased in NP_002684.1:p.Phe750Val mutants. Further, Contact Order analysis indicated that the Aromatic-sulfur interactions were destabilizing in the mutants. Overall, these in-silico analysis has revealed a destabilizing mutations suggesting pathogenic variants in POLG1 gene.