Cardiometabolic disease risk markers are increased following burn injury in children

IntroductionBurn injury in children causes prolonged systemic effects on physiology and metabolism leading to increased morbidity and mortality, yet much remains undefined regarding the metabolic trajectory towards specific health outcomes.MethodsA multi-platform strategy was implemented to evaluate the long-term immuno-metabolic consequences of burn injury combining metabolite, lipoprotein, and cytokine panels. Plasma samples from 36 children aged 4–8 years were collected 3 years after a burn injury together with 21 samples from non-injured age and sex matched controls. Three different 1H Nuclear Magnetic Resonance spectroscopic experiments were applied to capture information on plasma low molecular weight metabolites, lipoproteins, and α-1-acid glycoprotein.ResultsBurn injury was characterized by underlying signatures of hyperglycaemia, hypermetabolism and inflammation, suggesting disruption of multiple pathways relating to glycolysis, tricarboxylic acid cycle, amino acid metabolism and the urea cycle. In addition, very low-density lipoprotein sub-components were significantly reduced in participants with burn injury whereas small-dense low density lipoprotein particles were significantly elevated in the burn injured patient plasma compared to uninjured controls, potentially indicative of modified cardiometabolic risk after a burn. Weighted-node Metabolite Correlation Network Analysis was restricted to the significantly differential features (q <0.05) between the children with and without burn injury and demonstrated a striking disparity in the number of statistical correlations between cytokines, lipoproteins, and small molecular metabolites in the injured groups, with increased correlations between these groups.DiscussionThese findings suggest a ‘metabolic memory’ of burn defined by a signature of interlinked and perturbed immune and metabolic function. Burn injury is associated with a series of adverse metabolic changes that persist chronically and are independent of burn severity and this study demonstrates increased risk of cardiovascular disease in the long-term. These findings highlight a crucial need for improved longer term monitoring of cardiometabolic health in a vulnerable population of children that have undergone burn injury

Outcome after video-assisted thoracoscopic surgery and open pulmonary lobectomy in patients with low VO_{2 max}: a case-matched analysis from the ESTS database

n/

Integrative Modeling of Quantitative Plasma Lipoprotein, Metabolic, and Amino Acid Data Reveals a Multiorgan Pathological Signature of SARS-CoV-2 Infection.

The metabolic effects of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection on human blood plasma were characterized using multiplatform metabolic phenotyping with nuclear magnetic resonance (NMR) spectroscopy and liquid chromatography-mass spectrometry (LC-MS). Quantitative measurements of lipoprotein subfractions, α-1-acid glycoprotein, glucose, and biogenic amines were made on samples from symptomatic coronavirus disease 19 (COVID-19) patients who had tested positive for the SARS-CoV-2 virus (n = 17) and from age- and gender-matched controls (n = 25). Data were analyzed using an orthogonal-projections to latent structures (OPLS) method and used to construct an exceptionally strong (AUROC = 1) hybrid NMR-MS model that enabled detailed metabolic discrimination between the groups and their biochemical relationships. Key discriminant metabolites included markers of inflammation including elevated α-1-acid glycoprotein and an increased kynurenine/tryptophan ratio. There was also an abnormal lipoprotein, glucose, and amino acid signature consistent with diabetes and coronary artery disease (low total and HDL Apolipoprotein A1, low HDL triglycerides, high LDL and VLDL triglycerides), plus multiple highly significant amino acid markers of liver dysfunction (including the elevated glutamine/glutamate and Fischer's ratios) that present themselves as part of a distinct SARS-CoV-2 infection pattern. A multivariate training-test set model was validated using independent samples from additional SARS-CoV-2 positive patients and controls. The predictive model showed a sensitivity of 100% for SARS-CoV-2 positivity. The breadth of the disturbed pathways indicates a systemic signature of SARS-CoV-2 positivity that includes elements of liver dysfunction, dyslipidemia, diabetes, and coronary heart disease risk that are consistent with recent reports that COVID-19 is a systemic disease affecting multiple organs and systems. Metabolights study reference: MTBLS2014

Bacterial and metabolic phenotypes associated with inadequate response to ursodeoxycholic acid treatment in primary biliary cholangitis

Primary biliary cholangitis (PBC) is a chronic cholestatic liver disease with ursodeoxycholic acid (UDCA) as first-line treatment. Poor response to UDCA is associated with a higher risk of progressing to cirrhosis, but the underlying mechanisms are unclear. UDCA modulates the composition of primary and bacterial-derived bile acids (BAs). We characterized the phenotypic response to UDCA based on BA and bacterial profiles of PBC patients treated with UDCA. Patients from the UK-PBC cohort (n = 419) treated with UDCA for a minimum of 12-months were assessed using the Barcelona dynamic response criteria. BAs from serum, urine, and feces were analyzed using Ultra-High-Performance Liquid Chromatography-Mass Spectrometry and fecal bacterial composition measured using 16S rRNA gene sequencing. We identified 191 non-responders, 212 responders, and a subgroup of responders with persistently elevated liver biomarkers (n = 16). Responders had higher fecal secondary and tertiary BAs than non-responders and lower urinary bile acid abundances, with the exception of 12-dehydrocholic acid, which was higher in responders. The sub-group of responders with poor liver function showed lower alpha-diversity evenness, lower abundance of fecal secondary and tertiary BAs than the other groups and lower levels of phyla with BA-deconjugation capacity (Actinobacteriota/Actinomycetota, Desulfobacterota, Verrucomicrobiota) compared to responders. UDCA dynamic response was associated with an increased capacity to generate oxo-/epimerized secondary BAs. 12-dehydrocholic acid is a potential biomarker of treatment response. Lower alpha-diversity and lower abundance of bacteria with BA deconjugation capacity might be associated with an incomplete response to treatment in some patients

A meta-analysis of immune-cell fractions at high resolution reveals novel associations with common phenotypes and health outcomes

Background: Changes in cell-type composition of tissues are associated with a wide range of diseases and environmental risk factors and may be causally implicated in disease development and progression. However, these shifts in cell-type fractions are often of a low magnitude, or involve similar cell subtypes, making their reliable identification challenging. DNA methylation profiling in a tissue like blood is a promising approach to discover shifts in cell-type abundance, yet studies have only been performed at a relatively low cellular resolution and in isolation, limiting their power to detect shifts in tissue composition. Methods: Here we derive a DNA methylation reference matrix for 12 immune-cell types in human blood and extensively validate it with flow-cytometric count data and in whole-genome bisulfite sequencing data of sorted cells. Using this reference matrix, we perform a directional Stouffer and fixed effects meta-analysis comprising 23,053 blood samples from 22 different cohorts, to comprehensively map associations between the 12 immune-cell fractions and common phenotypes. In a separate cohort of 4386 blood samples, we assess associations between immune-cell fractions and health outcomes. Results: Our meta-analysis reveals many associations of cell-type fractions with age, sex, smoking and obesity, many of which we validate with single-cell RNA sequencing. We discover that naïve and regulatory T-cell subsets are higher in women compared to men, while the reverse is true for monocyte, natural killer, basophil, and eosinophil fractions. Decreased natural killer counts associated with smoking, obesity, and stress levels, while an increased count correlates with exercise and sleep. Analysis of health outcomes revealed that increased naïve CD4 + T-cell and N-cell fractions associated with a reduced risk of all-cause mortality independently of all major epidemiological risk factors and baseline co-morbidity. A machine learning predictor built only with immune-cell fractions achieved a C-index value for all-cause mortality of 0.69 (95%CI 0.67–0.72), which increased to 0.83 (0.80–0.86) upon inclusion of epidemiological risk factors and baseline co-morbidity. Conclusions: This work contributes an extensively validated high-resolution DNAm reference matrix for blood, which is made freely available, and uses it to generate a comprehensive map of associations between immune-cell fractions and common phenotypes, including health outcomes

OMICmAge : an integrative multi-omics approach to quantify biological age with electronic medical records

Biological aging is a multifactorial process involving complex interactions of cellular and biochemical processes that is reflected in omic profiles. Using common clinical laboratory measures in ~30,000 individuals from the MGB-Biobank, we developed a robust, predictive biological aging phenotype, EMRAge, that balances clinical biomarkers with overall mortality risk and can be broadly recapitulated across EMRs. We then applied elastic-net regression to model EMRAge with DNA-methylation (DNAm) and multiple omics, generating DNAmEMRAge and OMICmAge, respectively. Both biomarkers demonstrated strong associations with chronic diseases and mortality that outperform current biomarkers across our discovery (MGB-ABC, n=3,451) and validation TruDiagnostic, n=12,666) cohorts. Through the use of epigenetic biomarker proxies, OMICmAge has the unique advantage of expanding the predictive search space to include epigenomic, proteomic, metabolomic, and clinical data while distilling this in a measure with DNAm alone, providing opportunities to identify clinically-relevant interconnections central to the aging process

The evolution of lung cancer and impact of subclonal selection in TRACERx

Lung cancer is the leading cause of cancer-associated mortality worldwide. Here we analysed 1,644 tumour regions sampled at surgery or during follow-up from the first 421 patients with non-small cell lung cancer prospectively enrolled into the TRACERx study. This project aims to decipher lung cancer evolution and address the primary study endpoint: determining the relationship between intratumour heterogeneity and clinical outcome. In lung adenocarcinoma, mutations in 22 out of 40 common cancer genes were under significant subclonal selection, including classical tumour initiators such as TP53 and KRAS. We defined evolutionary dependencies between drivers, mutational processes and whole genome doubling (WGD) events. Despite patients having a history of smoking, 8% of lung adenocarcinomas lacked evidence of tobacco-induced mutagenesis. These tumours also had similar detection rates for EGFR mutations and for RET, ROS1, ALK and MET oncogenic isoforms compared with tumours in never-smokers, which suggests that they have a similar aetiology and pathogenesis. Large subclonal expansions were associated with positive subclonal selection. Patients with tumours harbouring recent subclonal expansions, on the terminus of a phylogenetic branch, had significantly shorter disease-free survival. Subclonal WGD was detected in 19% of tumours, and 10% of tumours harboured multiple subclonal WGDs in parallel. Subclonal, but not truncal, WGD was associated with shorter disease-free survival. Copy number heterogeneity was associated with extrathoracic relapse within 1 year after surgery. These data demonstrate the importance of clonal expansion, WGD and copy number instability in determining the timing and patterns of relapse in non-small cell lung cancer and provide a comprehensive clinical cancer evolutionary data resource

The evolution of non-small cell lung cancer metastases in TRACERx

Metastatic disease is responsible for the majority of cancer-related deaths. We report the longitudinal evolutionary analysis of 126 non-small cell lung cancer (NSCLC) tumours from 421 prospectively recruited patients in TRACERx who developed metastatic disease, compared with a control cohort of 144 non-metastatic tumours. In 25% of cases, metastases diverged early, before the last clonal sweep in the primary tumour, and early divergence was enriched for patients who were smokers at the time of initial diagnosis. Simulations suggested that early metastatic divergence more frequently occurred at smaller tumour diameters (less than 8 mm). Single-region primary tumour sampling resulted in 83% of late divergence cases being misclassified as early, highlighting the importance of extensive primary tumour sampling. Polyclonal dissemination, which was associated with extrathoracic disease recurrence, was found in 32% of cases. Primary lymph node disease contributed to metastatic relapse in less than 20% of cases, representing a hallmark of metastatic potential rather than a route to subsequent recurrences/disease progression. Metastasis-seeding subclones exhibited subclonal expansions within primary tumours, probably reflecting positive selection. Our findings highlight the importance of selection in metastatic clone evolution within untreated primary tumours, the distinction between monoclonal versus polyclonal seeding in dictating site of recurrence, the limitations of current radiological screening approaches for early diverging tumours and the need to develop strategies to target metastasis-seeding subclones before relapse

Genomic–transcriptomic evolution in lung cancer and metastasis

Intratumour heterogeneity (ITH) fuels lung cancer evolution, which leads to immune evasion and resistance to therapy. Here, using paired whole-exome and RNA sequencing data, we investigate intratumour transcriptomic diversity in 354 non-small cell lung cancer tumours from 347 out of the first 421 patients prospectively recruited into the TRACERx study. Analyses of 947 tumour regions, representing both primary and metastatic disease, alongside 96 tumour-adjacent normal tissue samples implicate the transcriptome as a major source of phenotypic variation. Gene expression levels and ITH relate to patterns of positive and negative selection during tumour evolution. We observe frequent copy number-independent allele-specific expression that is linked to epigenomic dysfunction. Allele-specific expression can also result in genomic–transcriptomic parallel evolution, which converges on cancer gene disruption. We extract signatures of RNA single-base substitutions and link their aetiology to the activity of the RNA-editing enzymes ADAR and APOBEC3A, thereby revealing otherwise undetected ongoing APOBEC activity in tumours. Characterizing the transcriptomes of primary–metastatic tumour pairs, we combine multiple machine-learning approaches that leverage genomic and transcriptomic variables to link metastasis-seeding potential to the evolutionary context of mutations and increased proliferation within primary tumour regions. These results highlight the interplay between the genome and transcriptome in influencing ITH, lung cancer evolution and metastasis

Antibodies against endogenous retroviruses promote lung cancer immunotherapy

B cells are frequently found in the margins of solid tumours as organized follicles in ectopic lymphoid organs called tertiary lymphoid structures (TLS). Although TLS have been found to correlate with improved patient survival and response to immune checkpoint blockade (ICB), the underlying mechanisms of this association remain elusive. Here we investigate lung-resident B cell responses in patients from the TRACERx 421 (Tracking Non-Small-Cell Lung Cancer Evolution Through Therapy) and other lung cancer cohorts, and in a recently established immunogenic mouse model for lung adenocarcinoma. We find that both human and mouse lung adenocarcinomas elicit local germinal centre responses and tumour-binding antibodies, and further identify endogenous retrovirus (ERV) envelope glycoproteins as a dominant anti-tumour antibody target. ERV-targeting B cell responses are amplified by ICB in both humans and mice, and by targeted inhibition of KRAS(G12C) in the mouse model. ERV-reactive antibodies exert anti-tumour activity that extends survival in the mouse model, and ERV expression predicts the outcome of ICB in human lung adenocarcinoma. Finally, we find that effective immunotherapy in the mouse model requires CXCL13-dependent TLS formation. Conversely, therapeutic CXCL13 treatment potentiates anti-tumour immunity and synergizes with ICB. Our findings provide a possible mechanistic basis for the association of TLS with immunotherapy response