Obesity Surgery Mortality Risk Score as a Predictor for Intensive Care Unit Admission in Patients Undergoing Laparoscopic Bariatric Surgery.

Background: Laparoscopic bariatric surgery provides many benefits including lower postoperative pain scores, reduced opioid consumption, shorter hospital stays, and improved quality of recovery. However, the anaesthetic management of obese patients requires caution in determining postoperative risk and in planning adequate postoperative pathways. Currently, there are no specific indications for intensive care unit (ICU) admission in this surgical population and most decisions are made on a case-by-case basis. The aim of this study is to investigate whether Obesity Surgery Mortality Risk Score (OS-MRS) is able to predict ICU admission in patients undergoing laparoscopic bariatric surgery (LBS). Methods: We retrospectively reviewed data of patients who underwent LBS during a 2-year period (2017-2019). The collected data included demographics, comorbidities and surgery-related variables. Postoperative ICU admission was decided via bariatric anaesthesiologists' evaluations, based on the high risk of postoperative cardiac or respiratory complications. Anaesthesia protocol was standardized. Logistic regression was used for statistical analysis. Results: ICU admission was required in 2% (n = 15) of the 763 patients. The intermediate risk group of the OS-MRS was detected in 84% of patients, while the American Society of Anaesthesiologists class III was reported in 80% of patients. A greater OS-MRS (p = 0.01), advanced age (p = 0.04), male gender (p = 0.001), longer duration of surgery (p = 0.0001), increased number of patient comorbidities (p = 0.002), and previous abdominal surgeries (p = 0.003) were predictive factors for ICU admission. Conclusions: ICU admission in obese patients undergoing LBS is predicted by OS-MRS together with age, male gender, number of comorbidities, previous abdominal surgeries, and duration of surgery.

Heteroplasmic mitochondrial DNA variants in cardiovascular diseases.

Mitochondria are implicated in the pathogenesis of cardiovascular diseases (CVDs) but the reasons for this are not well understood. Maternally-inherited population variants of mitochondrial DNA (mtDNA) which affect all mtDNA molecules (homoplasmic) are associated with cardiometabolic traits and the risk of developing cardiovascular disease. However, it is not known whether mtDNA mutations only affecting a proportion of mtDNA molecules (heteroplasmic) also play a role. To address this question, we performed a high-depth (~1000-fold) mtDNA sequencing of blood DNA in 1,399 individuals with hypertension (HTN), 1,946 with ischemic heart disease (IHD), 2,146 with ischemic stroke (IS), and 723 healthy controls. We show that the per individual burden of heteroplasmic single nucleotide variants (mtSNVs) increases with age. The age-effect was stronger for low-level heteroplasmies (heteroplasmic fraction, HF, 5-10%), likely reflecting acquired somatic events based on trinucleotide mutational signatures. After correcting for age and other confounders, intermediate heteroplasmies (HF 10-95%) were more common in hypertension, particularly involving non-synonymous variants altering the amino acid sequence of essential respiratory chain proteins. These findings raise the possibility that heteroplasmic mtSNVs play a role in the pathophysiology of hypertension.

Single-molecule mitochondrial DNA sequencing shows no evidence of CpG methylation in human cells and tissues.

Methylation on CpG residues is one of the most important epigenetic modifications of nuclear DNA, regulating gene expression. Methylation of mitochondrial DNA (mtDNA) has been studied using whole genome bisulfite sequencing (WGBS), but recent evidence has uncovered technical issues which introduce a potential bias during methylation quantification. Here, we validate the technical concerns of WGBS, and develop and assess the accuracy of a new protocol for mtDNA nucleotide variant-specific methylation using single-molecule Oxford Nanopore Sequencing (ONS). Our approach circumvents confounders by enriching for full-length molecules over nuclear DNA. Variant calling analysis against showed that 99.5% of homoplasmic mtDNA variants can be reliably identified providing there is adequate sequencing depth. We show that some of the mtDNA methylation signal detected by ONS is due to sequence-specific false positives introduced by the technique. The residual signal was observed across several human primary and cancer cell lines and multiple human tissues, but was always below the error threshold modelled using negative controls. We conclude that there is no evidence for CpG methylation in human mtDNA, thus resolving previous controversies. Additionally, we developed a reliable protocol to study epigenetic modifications of mtDNA at single-molecule and single-base resolution, with potential applications beyond CpG methylation.

Mitochondrial DNA heteroplasmy is modulated during oocyte development propagating mutation transmission.

Heteroplasmic mitochondrial DNA (mtDNA) mutations are a common cause of inherited disease, but a few recurrent mutations account for the vast majority of new families. The reasons for this are not known. We studied heteroplasmic mice transmitting m.5024C>T corresponding to a human pathogenic mutation. Analyzing 1167 mother-pup pairs, we show that m.5024C>T is preferentially transmitted from low to higher levels but does not reach homoplasmy. Single-cell analysis of the developing mouse oocytes showed the preferential increase in mutant over wild-type mtDNA in the absence of cell division. A similar inheritance pattern is seen in human pedigrees transmitting several pathogenic mtDNA mutations. In m.5024C>T mice, this can be explained by the preferential propagation of mtDNA during oocyte maturation, counterbalanced by purifying selection against high heteroplasmy levels. This could explain how a disadvantageous mutation in a carrier increases to levels that cause disease but fails to fixate, causing multigenerational heteroplasmic mtDNA disorders.

Use of whole genome sequencing to determine genetic basis of suspected mitochondrial disorders: cohort study.

OBJECTIVE: To determine whether whole genome sequencing can be used to define the molecular basis of suspected mitochondrial disease. DESIGN: Cohort study. SETTING: National Health Service, England, including secondary and tertiary care. PARTICIPANTS: 345 patients with suspected mitochondrial disorders recruited to the 100 000 Genomes Project in England between 2015 and 2018. INTERVENTION: Short read whole genome sequencing was performed. Nuclear variants were prioritised on the basis of gene panels chosen according to phenotypes, ClinVar pathogenic/likely pathogenic variants, and the top 10 prioritised variants from Exomiser. Mitochondrial DNA variants were called using an in-house pipeline and compared with a list of pathogenic variants. Copy number variants and short tandem repeats for 13 neurological disorders were also analysed. American College of Medical Genetics guidelines were followed for classification of variants. MAIN OUTCOME MEASURE: Definite or probable genetic diagnosis. RESULTS: A definite or probable genetic diagnosis was identified in 98/319 (31%) families, with an additional 6 (2%) possible diagnoses. Fourteen of the diagnoses (4% of the 319 families) explained only part of the clinical features. A total of 95 different genes were implicated. Of 104 families given a diagnosis, 39 (38%) had a mitochondrial diagnosis and 65 (63%) had a non-mitochondrial diagnosis. CONCLUSION: Whole genome sequencing is a useful diagnostic test in patients with suspected mitochondrial disorders, yielding a diagnosis in a further 31% after exclusion of common causes. Most diagnoses were non-mitochondrial disorders and included developmental disorders with intellectual disability, epileptic encephalopathies, other metabolic disorders, cardiomyopathies, and leukodystrophies. These would have been missed if a targeted approach was taken, and some have specific treatments.

Mitochondrial DNA variants modulate N-formylmethionine, proteostasis and risk of late-onset human diseases.

Mitochondrial DNA (mtDNA) variants influence the risk of late-onset human diseases, but the reasons for this are poorly understood. Undertaking a hypothesis-free analysis of 5,689 blood-derived biomarkers with mtDNA variants in 16,220 healthy donors, here we show that variants defining mtDNA haplogroups Uk and H4 modulate the level of circulating N-formylmethionine (fMet), which initiates mitochondrial protein translation. In human cytoplasmic hybrid (cybrid) lines, fMet modulated both mitochondrial and cytosolic proteins on multiple levels, through transcription, post-translational modification and proteolysis by an N-degron pathway, abolishing known differences between mtDNA haplogroups. In a further 11,966 individuals, fMet levels contributed to all-cause mortality and the disease risk of several common cardiovascular disorders. Together, these findings indicate that fMet plays a key role in common age-related disease through pleiotropic effects on cell proteostasis.

An atlas of mitochondrial DNA genotype-phenotype associations in the UK Biobank.

Mitochondrial DNA (mtDNA) variation in common diseases has been underexplored, partly due to a lack of genotype calling and quality-control procedures. Developing an at-scale workflow for mtDNA variant analyses, we show correlations between nuclear and mitochondrial genomic structures within subpopulations of Great Britain and establish a UK Biobank reference atlas of mtDNA-phenotype associations. A total of 260 mtDNA-phenotype associations were new (P < 1 × 10-5), including rs2853822 /m.8655 C>T (MT-ATP6) with type 2 diabetes, rs878966690 /m.13117 A>G (MT-ND5) with multiple sclerosis, 6 mtDNA associations with adult height, 24 mtDNA associations with 2 liver biomarkers and 16 mtDNA associations with parameters of renal function. Rare-variant gene-based tests implicated complex I genes modulating mean corpuscular volume and mean corpuscular hemoglobin. Seven traits had both rare and common mtDNA associations, where rare variants tended to have larger effects than common variants. Our work illustrates the value of studying mtDNA variants in common complex diseases and lays foundations for future large-scale mtDNA association studies.

Metabolic shift underlies recovery in reversible infantile respiratory chain deficiency.

Reversible infantile respiratory chain deficiency (RIRCD) is a rare mitochondrial myopathy leading to severe metabolic disturbances in infants, which recover spontaneously after 6-months of age. RIRCD is associated with the homoplasmic m.14674T>C mitochondrial DNA mutation; however, only ~ 1/100 carriers develop the disease. We studied 27 affected and 15 unaffected individuals from 19 families and found additional heterozygous mutations in nuclear genes interacting with mt-tRNAGlu including EARS2 and TRMU in the majority of affected individuals, but not in healthy carriers of m.14674T>C, supporting a digenic inheritance. Our transcriptomic and proteomic analysis of patient muscle suggests a stepwise mechanism where first, the integrated stress response associated with increased FGF21 and GDF15 expression enhances the metabolism modulated by serine biosynthesis, one carbon metabolism, TCA lipid oxidation and amino acid availability, while in the second step mTOR activation leads to increased mitochondrial biogenesis. Our data suggest that the spontaneous recovery in infants with digenic mutations may be modulated by the above described changes. Similar mechanisms may explain the variable penetrance and tissue specificity of other mtDNA mutations and highlight the potential role of amino acids in improving mitochondrial disease.

Genomic basis for RNA alterations in cancer.

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

A Pan-cancer Transcriptome Analysis Reveals Pervasive Regulation through Alternative Promoters.

Most human protein-coding genes are regulated by multiple, distinct promoters, suggesting that the choice of promoter is as important as its level of transcriptional activity. However, while a global change in transcription is recognized as a defining feature of cancer, the contribution of alternative promoters still remains largely unexplored. Here, we infer active promoters using RNA-seq data from 18,468 cancer and normal samples, demonstrating that alternative promoters are a major contributor to context-specific regulation of transcription. We find that promoters are deregulated across tissues, cancer types, and patients, affecting known cancer genes and novel candidates. For genes with independently regulated promoters, we demonstrate that promoter activity provides a more accurate predictor of patient survival than gene expression. Our study suggests that a dynamic landscape of active promoters shapes the cancer transcriptome, opening new diagnostic avenues and opportunities to further explore the interplay of regulatory mechanisms with transcriptional aberrations in cancer.

HmtDB 2016: data update, a better performing query system and human mitochondrial DNA haplogroup predictor.

The HmtDB resource hosts a database of human mitochondrial genome sequences from individuals with healthy and disease phenotypes. The database is intended to support both population geneticists as well as clinicians undertaking the task to assess the pathogenicity of specific mtDNA mutations. The wide application of next-generation sequencing (NGS) has provided an enormous volume of high-resolution data at a low price, increasing the availability of human mitochondrial sequencing data, which called for a cogent and significant expansion of HmtDB data content that has more than tripled in the current release. We here describe additional novel features, including: (i) a complete, user-friendly restyling of the web interface, (ii) links to the command-line stand-alone and web versions of the MToolBox package, an up-to-date tool to reconstruct and analyze human mitochondrial DNA from NGS data and (iii) the implementation of the Reconstructed Sapiens Reference Sequence (RSRS) as mitochondrial reference sequence. The overall update renders HmtDB an even more handy and useful resource as it enables a more rapid data access, processing and analysis. HmtDB is accessible at http://www.hmtdb.uniba.it/.

A multi-parametric workflow for the prioritization of mitochondrial DNA variants of clinical interest.

Assigning a pathogenic role to mitochondrial DNA (mtDNA) variants and unveiling the potential involvement of the mitochondrial genome in diseases are challenging tasks in human medicine. Assuming that rare variants are more likely to be damaging, we designed a phylogeny-based prioritization workflow to obtain a reliable pool of candidate variants for further investigations. The prioritization workflow relies on an exhaustive functional annotation through the mtDNA extraction pipeline MToolBox and includes Macro Haplogroup Consensus Sequences to filter out fixed evolutionary variants and report rare or private variants, the nucleotide variability as reported in HmtDB and the disease score based on several predictors of pathogenicity for non-synonymous variants. Cutoffs for both the disease score as well as for the nucleotide variability index were established with the aim to discriminate sequence variants contributing to defective phenotypes. The workflow was validated on mitochondrial sequences from Leber's Hereditary Optic Neuropathy affected individuals, successfully identifying 23 variants including the majority of the known causative ones. The application of the prioritization workflow to cancer datasets allowed to trim down the number of candidate for subsequent functional analyses, unveiling among these a high percentage of somatic variants. Prioritization criteria were implemented in both standalone ( http://sourceforge.net/projects/mtoolbox/ ) and web version ( https://mseqdr.org/mtoolbox.php ) of MToolBox.

Defining needs and goals of post-ICU care for trauma patients: preliminary study.

BACKGROUND: The aim of this study was to assess the long-term physical and psychological disabilities and their economic impact in severe trauma survivors. METHODS: Adult patients with Injury Severity Score >15 and Abbreviated Injury Scale ≤3 admitted to the ICU of a Level 1 trauma centre in the Lazio Region and discharged alive from hospital underwent a structured interview 12-24 months after the event. Self-reported somatic symptoms, autonomy, anxiety and depression were evaluated using a Likert-type Scale, Barthel Index and Hospital Anxiety and Depression Score (HADS), respectively. Patients' working and economic status were also investigated. RESULTS: A total of 32/58 patients matching the inclusion criteria were included in the final analysis. Eighteen patients (56%) reported at least a partial restriction in daily activities. Most common symptoms included muscle or joint pain, fatigue, and headache. All patients were receiving rehabilitation 1-2 years after the event. Fifty-eight percent of the patients spent more than €3600/year from their family budget for rehabilitation and medical care, however only 25% were receiving financial support from regional social services and 44% were unemployed at the time of the interview. Thirty patients (94%) had HADS Depression Score≥11. CONCLUSION: Survivors of severe trauma in our cohort had limited autonomy and need long-term rehabilitation. Most of them rely on private healthcare services with a significant financial impact on their family budget. Almost all patients had moderate to severe depression. Future post-ICU counseling services should facilitate access to rehabilitation and psychological support for these patients.

High-resolution genomic profiling of thyroid lesions uncovers preferential copy number gains affecting mitochondrial biogenesis loci in the oncocytic variants.

Oncocytic change is the result of aberrant mitochondrial hyperplasia, which may occur in both neoplastic and non-neoplastic cells and is not infrequent in the thyroid. Despite being a well-characterized histologic phenotype, the molecular causes underlying such a distinctive cellular change are poorly understood. To identify potential genetic causes for the oncocytic phenotype in thyroid, we analyzed copy number alterations in a set of oncocytic (n=21) and non-oncocytic (n=20) thyroid lesions by high-resolution microarray-based comparative genomic hybridization (aCGH). Each group comprised lesions of diverse histologic types, including hyperplastic nodules, adenomas and carcinomas. Unsupervised hierarchical clustering of categorical aCGH data resulted in two distinct branches, one of which was significantly enriched for samples with the oncocytic phenotype, regardless of histologic type. Analysis of aCGH events showed that the oncocytic group harbored a significantly higher number of genes involved in copy number gains, when compared to that of conventional thyroid lesions. Functional annotation demonstrated an enrichment for copy number gains that affect genes encoding activators of mitochondrial biogenesis in oncocytic cases but not in their non-oncocytic counterparts. Taken together, our data suggest that genomic alterations may represent additional/alternative mechanisms underlying the development of the oncocytic phenotype in the thyroid.

A comprehensive characterization of mitochondrial DNA mutations in glioblastoma multiforme.

Glioblastoma multiforme (GBM) is the most malignant brain cancer in adults, with a poor prognosis, whose molecular stratification still represents a challenge in pathology and clinics. On the other hand, mitochondrial DNA (mtDNA) mutations have been found in most tumors as modifiers of the bioenergetics state, albeit in GBM a characterization of the mtDNA status is lacking to date. Here, a characterization of the burden of mtDNA mutations in GBM samples was performed. First, investigation of tumor-specific vs. non tumor-specific mutations was carried out with the MToolBox bioinformatics pipeline by analyzing 45 matched tumor/blood samples, from whole genome or whole exome sequencing datasets obtained from The Cancer Genome Atlas (TCGA) consortium. Additionally, the entire mtDNA sequence was obtained in a dataset of 104 fresh-frozen GBM samples. Mitochondrial mutations with potential pathogenic interest were prioritized based on heteroplasmic fraction, nucleotide variability, and in silico prediction of pathogenicity. A preliminary biochemical analysis of the activity of mitochondrial respiratory complexes was also performed on fresh-frozen GBM samples. Although a high number of mutations was detected, we report that the large majority of them does not pass the prioritization filters. Therefore, a relatively limited burden of pathogenic mutations is indeed carried by GBM, which did not appear to determine a general impairment of the respiratory chain. This article is part of a Directed Issue entitled: Energy Metabolism Disorders and Therapies.

Extraction and annotation of human mitochondrial genomes from 1000 Genomes Whole Exome Sequencing data.

BACKGROUND: Whole Exome Sequencing (WES) is one of the most used and cost-effective next generation technologies that allows sequencing of all nuclear exons. Off-target regions may be captured if they present high sequence similarity with baits. Bioinformatics tools have been optimized to retrieve a large amount of WES off-target mitochondrial DNA (mtDNA), by exploiting the aspecificity of probes, partially overlapping to Nuclear mitochondrial Sequences (NumtS). The 1000 Genomes project represents one of the widest resources to extract mtDNA sequences from WES data, considering the large effort the scientific community is undertaking to reconstruct human population history using mtDNA as marker, and the involvement of mtDNA in pathology. RESULTS: A previously published pipeline aimed at assembling mitochondrial genomes from off-target WES reads and further improved to detect insertions and deletions (indels) and heteroplasmy in a dataset of 1242 samples from the 1000 Genomes project, enabled to obtain a nearly complete mitochondrial genome from 943 samples (76% analyzed exomes). The robustness of our computational strategy was highlighted by the reduction of reads amount recognized as mitochondrial in the original annotation produced by the Consortium, due to NumtS filtering. CONCLUSIONS: To the best of our knowledge, this is likely the most extended population-scale mitochondrial genotyping in humans enriched with the estimation of heteroplasmies.

MToolBox: a highly automated pipeline for heteroplasmy annotation and prioritization analysis of human mitochondrial variants in high-throughput sequencing.

MOTIVATION: The increasing availability of mitochondria-targeted and off-target sequencing data in whole-exome and whole-genome sequencing studies (WXS and WGS) has risen the demand of effective pipelines to accurately measure heteroplasmy and to easily recognize the most functionally important mitochondrial variants among a huge number of candidates. To this purpose, we developed MToolBox, a highly automated pipeline to reconstruct and analyze human mitochondrial DNA from high-throughput sequencing data. RESULTS: MToolBox implements an effective computational strategy for mitochondrial genomes assembling and haplogroup assignment also including a prioritization analysis of detected variants. MToolBox provides a Variant Call Format file featuring, for the first time, allele-specific heteroplasmy and annotation files with prioritized variants. MToolBox was tested on simulated samples and applied on 1000 Genomes WXS datasets. AVAILABILITY AND IMPLEMENTATION: MToolBox package is available at https://sourceforge.net/projects/mtoolbox/.

Evidence of association of human papillomavirus with prognosis worsening in glioblastoma multiforme.

BACKGROUND: Glioblastoma multiforme (GBM) is the most malignant brain tumor in adults, but its etiology still remains unknown. Recently, a role of viruses such as cytomegalovirus and JC virus in gliomagenesis has been suggested. Since human papillomavirus (HPV) is considered the most common oncogenic virus in humans, we evaluated its occurrence in GBM samples. MATERIAL AND METHODS: Fifty-two formalin-fixed paraffin-embedded primary glioblastoma specimens were retrospectively analyzed. The presence of HPV genome on tumor DNA was assessed by MY/GP nested PCR. Confirmation of HPV detection was obtained by chromogenic in situ hybridization (CISH) and immunohistochemistry (IHC) with an antibody directed against the L1 capsidic protein. Finally, univariate and multivariate proportional-hazards models were used to compare the risk of death among HPV-positive and HPV-negative patients. RESULTS: Strikingly, viral DNA was detected after PCR in 12 cases (23%). HPV16 genome was present in 25% infected samples, whereas the remaining samples tested positive for HPV6. CISH confirmed positivity in all infected samples for which enough material was available. Moreover, IHC positivity suggested that production of viral proteins from HPV genome is an ongoing process in GBM cancer cells. Finally an association between HPV infection and a worse prognosis was found in patients upon age stratification with a univariate analysis (HR, 2.10; 95% CI, 1.00-4.44; log-rank P = .045). CONCLUSIONS: HPV infection status may be considered an independent prognostic factor in GBM patients and suggests that prevention may be considered, should HPV be recognized as a causative agent in gliomagenesis.