Development and validation of a mortality risk prediction model for chronic obstructive pulmonary disease: a cross-sectional study using probabilistic graphical modelling.

Background: Chronic Obstructive Pulmonary Disease (COPD) is a leading cause of mortality. Predicting mortality risk in patients with COPD can be important for disease management strategies. Although all-cause mortality predictors have been developed previously, limited research exists on factors directly affecting COPD-specific mortality. Methods: In a retrospective study, we used probabilistic graphs to analyse clinical cross-sectional data (COPDGene cohort), including demographics, spirometry, quantitative chest imaging, and symptom features, as well as gene expression data. COPDGene recruited current and former smokers, aged 45-80 years with >10 pack-years smoking history, from across the USA (Phase 1, 11/2007-4/2011) and invited them for a follow-up visit (Phase 2, 7/2013-7/2017). ECLIPSE cohort recruited current and former smokers (COPD patients and controls from USA and Europe), aged 45-80 with smoking history >10 pack-years (12/2005-11/2007). We applied graphical models on multi-modal data COPDGene Phase 1 participants to identify factors directly affecting all-cause and COPD-specific mortality (primary outcomes); and on Phase 2 follow-up cohort to identify additional molecular and social factors affecting mortality. We used penalized Cox regression with features selected by the causal graph to build VAPORED, a mortality risk prediction model. VAPORED was compared to existing scores (BODE: BMI, airflow obstruction, dyspnoea, exercise capacity; ADO: age, dyspnoea, airflow obstruction) on the ability to rank individuals by mortality risk, using four evaluation metrics (concordance, concordance probability estimate (CPE), cumulative/dynamic (C/D) area under the receiver operating characteristic curve (AUC), and integrated C/D AUC). The results were validated in ECLIPSE. Findings: Graphical models, applied on the COPDGene Phase 1 samples (n = 8610), identified 11 and 7 variables directly linked to all-cause and COPD-specific mortality, respectively. Although many appear in both models, non-lung comorbidities appear only in the all-cause model, while forced vital capacity (FVC %predicted) appears in COPD-specific mortality model only. Additionally, the graph model of Phase 2 data (n = 3182) identified internet access, CD4 T cells and platelets to be linked to lower mortality risk. Furthermore, using the 7 variables linked to COPD-specific mortality (forced expiratory volume in 1 s/forced vital capacity (FEV1/FVC) ration, FVC %predicted, age, history of pneumonia, oxygen saturation, 6-min walk distance, dyspnoea) we developed VAPORED mortality risk score, which we validated on the ECLIPSE cohort (3-yr all-cause mortality data, n = 2312). VAPORED performed significantly better than ADO, BODE, and updated BODE indices in predicting all-cause mortality in ECLIPSE in terms of concordance (VAPORED [0.719] vs ADO [0.693; FDR p-value 0.014], BODE [0.695; FDR p-value 0.020], and updated BODE [0.694; FDR p-value 0.021]); CPE (VAPORED [0.714] vs ADO [0.673; FDR p-value <0.0001], BODE [0.662; FDR p-value <0.0001], and updated BODE [0.646; FDR p-value <0.0001]); 3-year C/D AUC (VAPORED [0.728] vs ADO [0.702; FDR p-value 0.017], BODE [0.704; FDR p-value 0.021], and updated BODE [0.703; FDR p-value 0.024]); integrated C/D AUC (VAPORED [0.723] vs ADO [0.698; FDR p-value 0.047], BODE [0.695; FDR p-value 0.024], and updated BODE [0.690; FDR p-value 0.021]). Finally, we developed a web tool to help clinicians calculate VAPORED mortality risk and compare it to ADO and BODE predictions. Interpretation: Our work is an important step towards improving our identification of high-risk patients and generating hypotheses of potential biological mechanisms and social factors driving mortality in patients with COPD at the population level. The main limitation of our study is the fact that the analysed datasets consist of older people with extensive smoking history and limited racial diversity. Thus, the results are relevant to high-risk individuals or those diagnosed with COPD and the VAPORED score is validated for them. Funding: This research was supported by NIH [NHLBI, NLM]. The COPDGene study is supported by the COPD Foundation, through grants from AstraZeneca, Bayer Pharmaceuticals, Boehringer Ingelheim, Genentech, GlaxoSmithKline, Novartis, Pfizer and Sunovion.

A hybrid constrained continuous optimization approach for optimal causal discovery from biological data.

MOTIVATION: Understanding causal effects is a fundamental goal of science and underpins our ability to make accurate predictions in unseen settings and conditions. While direct experimentation is the gold standard for measuring and validating causal effects, the field of causal graph theory offers a tantalizing alternative: extracting causal insights from observational data. Theoretical analysis has shown that this is indeed possible, given a large dataset and if certain conditions are met. However, biological datasets, frequently, do not meet such requirements but evaluation of causal discovery algorithms is typically performed on synthetic datasets, which they meet all requirements. Thus, real-life datasets are needed, in which the causal truth is reasonably known. In this work we first construct such a large-scale real-life dataset and then we perform on it a comprehensive benchmarking of various causal discovery methods. RESULTS: We find that the PC algorithm is particularly accurate at estimating causal structure, including the causal direction which is critical for biological applicability. However, PC does only produces cause-effect directionality, but not estimates of causal effects. We propose PC-NOTEARS (PCnt), a hybrid solution, which includes the PC output as an additional constraint inside the NOTEARS optimization. This approach combines PC algorithm's strengths in graph structure prediction with the NOTEARS continuous optimization to estimate causal effects accurately. PCnt achieved best aggregate performance across all structural and effect size metrics. AVAILABILITY AND IMPLEMENTATION: https://github.com/zhu-yh1/PC-NOTEARS.

Identification of factors directly linked to incident chronic obstructive pulmonary disease: A causal graph modeling study.

BACKGROUND: Beyond exposure to cigarette smoking and aging, the factors that influence lung function decline to incident chronic obstructive pulmonary disease (COPD) remain unclear. Advancements have been made in categorizing COPD into emphysema and airway predominant disease subtypes; however, predicting which healthy individuals will progress to COPD is difficult because they can exhibit profoundly different disease trajectories despite similar initial risk factors. This study aimed to identify clinical, genetic, and radiological features that are directly linked-and subsequently predict-abnormal lung function. METHODS AND FINDINGS: We employed graph modeling on 2,643 COPDGene participants (aged 45 to 80 years, 51.25% female, 35.1% African Americans; enrollment 11/2007-4/2011) with smoking history but normal spirometry at study enrollment to identify variables that are directly linked to future lung function abnormalities. We developed logistic regression and random forest predictive models for distinguishing individuals who maintain lung function from those who decline. Of the 131 variables analyzed, 6 were identified as informative to future lung function abnormalities, namely forced expiratory flow in the middle range (FEF25-75%), average lung wall thickness in a 10 mm radius (Pi10), severe emphysema, age, sex, and height. We investigated whether these features predict individuals leaving GOLD 0 status (normal spirometry according to Global Initiative for Obstructive Lung Disease (GOLD) criteria). Linear models, trained with these features, were quite predictive (area under receiver operator characteristic curve or AUROC = 0.75). Random forest predictors performed similarly to logistic regression (AUROC = 0.7), indicating that no significant nonlinear effects were present. The results were externally validated on 150 participants from Specialized Center for Clinically Oriented Research (SCCOR) cohort (aged 45 to 80 years, 52.7% female, 4.7% African Americans; enrollment: 7/2007-12/2012) (AUROC = 0.89). The main limitation of longitudinal studies with 5- and 10-year follow-up is the introduction of mortality bias that disproportionately affects the more severe cases. However, our study focused on spirometrically normal individuals, who have a lower mortality rate. Another limitation is the use of strict criteria to define spirometrically normal individuals, which was unavoidable when studying factors associated with changes in normalized forced expiratory volume in 1 s (FEV1%predicted) or the ratio of FEV1/FVC (forced vital capacity). CONCLUSIONS: This study took an agnostic approach to identify which baseline measurements differentiate and predict the early stages of lung function decline in individuals with previous smoking history. Our analysis suggests that emphysema affects obstruction onset, while airway predominant pathology may play a more important role in future FEV1 (%predicted) decline without obstruction, and FEF25-75% may affect both.

The transcriptome of CD14 + CD163 - HLA-DR low monocytes predicts mortality in Idiopathic Pulmonary Fibrosis.

Rationale: The association between immune-cell-specific transcriptomic profiles and Idiopathic Pulmonary Fibrosis (IPF) mortality is unknown. Objectives: To determine immune-cell-specific transcriptomic profiles associated with IPF mortality. Methods: We profiled peripheral blood mononuclear cells (PBMC) in 18 participants [University of South Florida: IPF, COVID-19, post-COVID-19 Interstitial Lung Disease (Post-COVID-19 ILD), controls] by single-cell RNA sequencing (scRNA-seq) and identified 16 immune-cell-specific transcriptomic profiles. The Scoring Algorithm of Molecular Subphenotypes (SAMS) was used to calculate Up-scores based on these 16 gene profiles. Their association with outcomes was investigated in peripheral blood, Bronchoalveolar Lavage (BAL) and lung tissue of N=416 IPF patients from six cohorts. Findings were validated in an independent IPF, PBMC scRNA-seq dataset (N=38). Measurements and main results: Cox-regression models demonstrated that 230 genes from CD14 + CD163 - HLA-DR low circulating monocytes predicted IPF mortality [Pittsburgh (p=0.02), Chicago (p=0.003)]. PBMC proportions of CD14 + CD163 - HLA-DR low monocytes were higher in progressive versus stable IPF (Yale, 0.13±0.05 versus 0.09±0.05, p=0.034). Receiving operating characteristic identified a 230 gene, Up-score >41.84 (Pittsburgh) predictive of mortality in Chicago (HR: 6.58, 95%CI: 2.15-20.13, p=0.001) and in pooled analysis of BAL cohorts (HR: 2.20, 95%CI: 1.44-3.37, p=0.0003). High-risk patients had decreased expression of the T-cell co-stimulatory genes CD28 , ICOS , ITK and LCK (Pittsburgh and Chicago, p<0.01). 230 gene-up-scores negatively correlated with Forced Vital Capacity (FVC) in IPF lung tissues (LGRC, rho=-0.2, p=0.02). Results were replicated using a subset of 13 genes from the 230-gene signature (pooled PBMC cohorts - HR: 5.34, 95%CI: 2.83-10.06, p<0.0001). Conclusions: The transcriptome of CD14 + CD163 - HLA-DR low monocytes is associated with increased IPF mortality.

A plasma peptidomic signature reveals extracellular matrix remodeling and predicts prognosis in alcohol-associated hepatitis.

BACKGROUND: Alcohol-associated hepatitis (AH) is plagued with high mortality and difficulty in identifying at-risk patients. The extracellular matrix undergoes significant remodeling during inflammatory liver injury and could potentially be used for mortality prediction. METHODS: EDTA plasma samples were collected from patients with AH (n = 62); Model for End-Stage Liver Disease score defined AH severity as moderate (12-20; n = 28) and severe (>20; n = 34). The peptidome data were collected by high resolution, high mass accuracy UPLC-MS. Univariate and multivariate analyses identified differentially abundant peptides, which were used for Gene Ontology, parent protein matrisomal composition, and protease involvement. Machine-learning methods were used to develop mortality predictors. RESULTS: Analysis of plasma peptides from patients with AH and healthy controls identified over 1600 significant peptide features corresponding to 130 proteins. These were enriched for extracellular matrix fragments in AH samples, likely related to the turnover of hepatic-derived proteins. Analysis of moderate versus severe AH peptidomes was dominated by changes in peptides from collagen 1A1 and fibrinogen A proteins. The dominant proteases for the AH peptidome spectrum appear to be CAPN1 and MMP12. Causal graphical modeling identified 3 peptides directly linked to 90-day mortality in >90% of the learned graphs. These peptides improved the accuracy of mortality prediction over the Model for End-Stage Liver Disease score and were used to create a clinically applicable mortality prediction assay. CONCLUSIONS: A signature based on plasma peptidome is a novel, noninvasive method for prognosis stratification in patients with AH. Our results could also lead to new mechanistic and/or surrogate biomarkers to identify new AH mechanisms.

Cellular and transcriptional profiles of peripheral blood mononuclear cells pre-vaccination predict immune response to preventative MUC1 vaccine.

A single arm trial (NCT007773097) and a double-blind, placebo controlled randomized trial ( NCT02134925 ) were conducted in individuals with a history of advanced colonic adenoma to test the safety and immunogenicity of the MUC1 tumor antigen vaccine and its potential to prevent new adenomas. These were the first two trials of a non-viral cancer vaccine administered in the absence of cancer. The vaccine was safe and strongly immunogenic in 43% (NCT007773097) and 25% ( NCT02134925 ) of participants. The lack of response in a significant number of participants suggested, for the first time, that even in a premalignant setting, the immune system may have already been exposed to some level of suppression previously reported only in cancer. Single-cell RNA-sequencing (scRNA-seq) on banked pre-vaccination peripheral blood mononuclear cells (PBMCs) from 16 immune responders and 16 non-responders identified specific cell types, genes, and pathways of a productive vaccine response. Responders had a significantly higher percentage of CD4+ naive T cells pre-vaccination, but a significantly lower percentage of CD8+ T effector memory (TEM) cells and CD16+ monocytes. Differential gene expression (DGE) and transcription factor inference analysis showed a higher level of expression of T cell activation genes, such as Fos and Jun, in CD4+ naive T cells, and pathway analysis showed enriched signaling activity in responders. Furthermore, Bayesian network analysis suggested that these genes were mechanistically connected to response. Our analyses identified several immune mechanisms and candidate biomarkers to be further validated as predictors of immune responses to a preventative cancer vaccine that could facilitate selection of individuals likely to benefit from a vaccine or be used to improve vaccine responses.

Longitudinal multicompartment characterization of host-microbiota interactions in patients with acute respiratory failure.

Critical illness can significantly alter the composition and function of the human microbiome, but few studies have examined these changes over time. Here, we conduct a comprehensive analysis of the oral, lung, and gut microbiota in 479 mechanically ventilated patients (223 females, 256 males) with acute respiratory failure. We use advanced DNA sequencing technologies, including Illumina amplicon sequencing (utilizing 16S and ITS rRNA genes for bacteria and fungi, respectively, in all sample types) and Nanopore metagenomics for lung microbiota. Our results reveal a progressive dysbiosis in all three body compartments, characterized by a reduction in microbial diversity, a decrease in beneficial anaerobes, and an increase in pathogens. We find that clinical factors, such as chronic obstructive pulmonary disease, immunosuppression, and antibiotic exposure, are associated with specific patterns of dysbiosis. Interestingly, unsupervised clustering of lung microbiota diversity and composition by 16S independently predicted survival and performed better than traditional clinical and host-response predictors. These observations are validated in two separate cohorts of COVID-19 patients, highlighting the potential of lung microbiota as valuable prognostic biomarkers in critical care. Understanding these microbiome changes during critical illness points to new opportunities for microbiota-targeted precision medicine interventions.

Pre-vaccination transcriptomic profiles of immune responders to the MUC1 peptide vaccine for colon cancer prevention.

Self-antigens abnormally expressed on tumors, such as MUC1, have been targeted by therapeutic cancer vaccines. We recently assessed in two clinical trials in a preventative setting whether immunity induced with a MUC1 peptide vaccine could reduce high colon cancer risk in individuals with a history of premalignant colon adenomas. In both trials, there were immune responders and non-responders to the vaccine. Here we used PBMC pre-vaccination and 2 weeks after the first vaccine of responders and non-responders selected from both trials to identify early biomarkers of immune response involved in long-term memory generation and prevention of adenoma recurrence. We performed flow cytometry, phosflow, and differential gene expression analyses on PBMCs collected from MUC1 vaccine responders and non-responders pre-vaccination and two weeks after the first of three vaccine doses. MUC1 vaccine responders had higher frequencies of CD4 cells pre-vaccination, increased expression of CD40L on CD8 and CD4 T-cells, and a greater increase in ICOS expression on CD8 T-cells. Differential gene expression analysis revealed that iCOSL, PI3K AKT MTOR, and B-cell signaling pathways are activated early in response to the MUC1 vaccine. We identified six specific transcripts involved in elevated antigen presentation, B-cell activation, and NF-kB1 activation that were directly linked to finding antibody response at week 12. Finally, a model using these transcripts was able to predict non-responders with accuracy. These findings suggest that individuals who can be predicted to respond to the MUC1 vaccine, and potentially other vaccines, have greater readiness in all immune compartments to present and respond to antigens. Predictive biomarkers of MUC1 vaccine response may lead to more effective vaccines tailored to individuals with high risk for cancer but with varying immune fitness.

The Florida Scoring System for stratifying children with suspected Sjögren's disease: a cross-sectional machine learning study.

BACKGROUND: Childhood Sjögren's disease is a rare, underdiagnosed, and poorly-understood condition. By integrating machine learning models on a paediatric cohort in the USA, we aimed to develop a novel system (the Florida Scoring System) for stratifying symptomatic paediatric patients with suspected Sjögren's disease. METHODS: This cross-sectional study was done in symptomatic patients who visited the Department of Pediatric Rheumatology at the University of Florida, FL, USA. Eligible patients were younger than 18 years or had symptom onset before 18 years of age. Patients with confirmed diagnosis of another autoimmune condition or infection with a clear aetiological microorganism were excluded. Eligible patients underwent comprehensive examinations to rule out or diagnose childhood Sjögren's disease. We used latent class analysis with clinical and laboratory variables to detect heterogeneous patient classes. Machine learning models, including random forest, gradient-boosted decision tree, partial least square discriminatory analysis, least absolute shrinkage and selection operator-penalised ordinal regression, artificial neural network, and super learner were used to predict patient classes and rank the importance of variables. Causal graph learning selected key features to build the final Florida Scoring System. The predictors for all models were the clinical and laboratory variables and the outcome was the definition of patient classes. FINDINGS: Between Jan 16, 2018, and April 28, 2022, we screened 448 patients for inclusion. After excluding 205 patients due to symptom onset later than 18 years of age, we recruited 243 patients into our cohort. 26 patients were excluded because of confirmed diagnosis of a disorder other than Sjögren's disease, and 217 patients were included in the final analysis. Median age at diagnosis was 15 years (IQR 11-17). 155 (72%) of 216 patients were female and 61 (28%) were male, 167 (79%) of 212 were White, and 20 (9%) of 213 were Hispanic, Latino, or Spanish. The latent class analysis identified three distinct patient classes: class I (dryness dominant with positive tests, n=27), class II (high symptoms with negative tests, n=98), and class III (low symptoms with negative tests, n=92). Machine learning models accurately predicted patient class and ranked variable importance consistently. The causal graphical model discovered key features for constructing the Florida Scoring System. INTERPRETATION: The Florida Scoring System is a paediatrician-friendly tool that can be used to assist classification and long-term monitoring of suspected childhood Sjögren's disease. The resulting stratification has important implications for clinical management, trial design, and pathobiological research. We found a highly symptomatic patient group with negative serology and diagnostic profiles, which warrants clinical attention. We further revealed that salivary gland ultrasonography can be a non-invasive alternative to minor salivary gland biopsy in children. The Florida Scoring System requires validation in larger prospective paediatric cohorts. FUNDING: National Institute of Dental and Craniofacial Research, National Institute of Arthritis, Musculoskeletal and Skin Diseases, National Heart, Lung, and Blood Institute, and Sjögren's Foundation.

Disentangling Predictors of COPD Mortality with Probabilistic Graphical Models.

Background-Research question: Chronic Obstructive Pulmonary Disease (COPD) is a leading cause of mortality. Predicting mortality risk in COPD patients can be important for disease management strategies. Although scores for all-cause mortality have been developed previously, there is limited research on factors that may directly affect COPD-specific mortality. Study design-Methods: used probabilistic (causal) graphs to analyze clinical baseline COPDGene data, including demographics, spirometry, quantitative chest imaging, and symptom features, as well as gene expression data (from year-5). Results: We identified factors linked to all-cause and COPD-specific mortality. Although many were similar, there were differences in certain comorbidities (all-cause mortality model only) and forced vital capacity (COPD-specific mortality model only). Using our results, we developed VAPORED , a 7-variable COPD-specific mortality risk score, which we validated using the ECLIPSE 3-yr mortality data. We showed that the new model is more accurate than the existing ADO, BODE, and updated BODE indices. Additionally, we identified biological signatures linked to all-cause mortality, including a plasma cell mediated component. Finally, we developed a web page to help clinicians calculate mortality risk using VAPORED, ADO, and BODE indices. Interpretation: Given the importance of predicting COPD-specific and all-cause mortality risk in COPD patients, we showed that probabilistic graphs can identify the features most directly affecting them, and be used to build new, more accurate models of mortality risk. Novel biological features affecting mortality were also identified. This is an important step towards improving our identification of high-risk patients and potential biological mechanisms that drive COPD mortality.

Lung Epithelium Releases Growth Differentiation Factor 15 in Response to Pathogen-mediated Injury.

GDF15 (growth differentiation factor 15) is a stress cytokine with several proposed roles, including support of stress erythropoiesis. Higher circulating GDF15 levels are prognostic of mortality during acute respiratory distress syndrome, but the cellular sources and downstream effects of GDF15 during pathogen-mediated lung injury are unclear. We quantified GDF15 in lower respiratory tract biospecimens and plasma from patients with acute respiratory failure. Publicly available data from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection were reanalyzed. We used mouse models of hemorrhagic acute lung injury mediated by Pseudomonas aeruginosa exoproducts in wild-type mice and mice genetically deficient for Gdf15 or its putative receptor, Gfral. In critically ill humans, plasma levels of GDF15 correlated with lower respiratory tract levels and were higher in nonsurvivors. SARS-CoV-2 infection induced GDF15 expression in human lung epithelium, and lower respiratory tract GDF15 levels were higher in coronavirus disease (COVID-19) nonsurvivors. In mice, intratracheal P. aeruginosa type II secretion system exoproducts were sufficient to induce airspace and plasma release of GDF15, which was attenuated with epithelial-specific deletion of Gdf15. Mice with global Gdf15 deficiency had decreased airspace hemorrhage, an attenuated cytokine profile, and an altered lung transcriptional profile during injury induced by P. aeruginosa type II secretion system exoproducts, which was not recapitulated in mice deficient for Gfral. Airspace GDF15 reconstitution did not significantly modulate key lung cytokine levels but increased circulating erythrocyte counts. Lung epithelium releases GDF15 during pathogen injury, which is associated with plasma levels in humans and mice and can increase erythrocyte counts in mice, suggesting a novel lung-blood communication pathway.

Enabling the clinical application of artificial intelligence in genomics: a perspective of the AMIA Genomics and Translational Bioinformatics Workgroup.

OBJECTIVE: Given the importance AI in genomics and its potential impact on human health, the American Medical Informatics Association-Genomics and Translational Biomedical Informatics (GenTBI) Workgroup developed this assessment of factors that can further enable the clinical application of AI in this space. PROCESS: A list of relevant factors was developed through GenTBI workgroup discussions in multiple in-person and online meetings, along with review of pertinent publications. This list was then summarized and reviewed to achieve consensus among the group members. CONCLUSIONS: Substantial informatics research and development are needed to fully realize the clinical potential of such technologies. The development of larger datasets is crucial to emulating the success AI is achieving in other domains. It is important that AI methods do not exacerbate existing socio-economic, racial, and ethnic disparities. Genomic data standards are critical to effectively scale such technologies across institutions. With so much uncertainty, complexity and novelty in genomics and medicine, and with an evolving regulatory environment, the current focus should be on using these technologies in an interface with clinicians that emphasizes the value each brings to clinical decision-making.

Early Evidence of Chronic Obstructive Pulmonary Disease Obscured by Race-Specific Prediction Equations.

Rationale: The identification of early chronic obstructive pulmonary disease (COPD) is essential to appropriately counsel patients regarding smoking cessation, provide symptomatic treatment, and eventually develop disease-modifying treatments. Disease severity in COPD is defined using race-specific spirometry equations. These may disadvantage non-White individuals in diagnosis and care. Objectives: Determine the impact of race-specific equations on African American (AA) versus non-Hispanic White individuals. Methods: Cross-sectional analyses of the COPDGene (Genetic Epidemiology of Chronic Obstructive Pulmonary Disease) cohort were conducted, comparing non-Hispanic White (n = 6,766) and AA (n = 3,366) participants for COPD manifestations. Measurements and Main Results: Spirometric classifications using race-specific, multiethnic, and "race-reversed" prediction equations (NHANES [National Health and Nutrition Examination Survey] and Global Lung Function Initiative "Other" and "Global") were compared, as were respiratory symptoms, 6-minute-walk distance, computed tomography imaging, respiratory exacerbations, and St. George's Respiratory Questionnaire. Application of different prediction equations to the cohort resulted in different classifications by stage, with NHANES and Global Lung Function Initiative race-specific equations being minimally different, but race-reversed equations moving AA participants to more severe stages and especially between the Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage 0 and preserved ratio impaired spirometry groups. Classification using the established NHANES race-specific equations demonstrated that for each of GOLD stages 1-4, AA participants were younger, had fewer pack-years and more current smoking, but had more exacerbations, shorter 6-minute-walk distance, greater dyspnea, and worse BODE (body mass index, airway obstruction, dyspnea, and exercise capacity) scores and St. George's Respiratory Questionnaire scores. Differences were greatest in GOLD stages 1 and 2. Race-reversed equations reclassified 774 AA participants (43%) from GOLD stage 0 to preserved ratio impaired spirometry. Conclusions: Race-specific equations underestimated disease severity among AA participants. These effects were particularly evident in early disease and may result in late detection of COPD.

CellularPotts.jl: simulating multiscale cellular models in Julia.

SUMMARY: CellularPotts.jl is a software package written in Julia to simulate biological cellular processes such as division, adhesion, and signaling. Accurately modeling and predicting these simple processes is crucial because they facilitate more complex biological phenomena related to important disease states like tumor growth, wound healing, and infection. Here we take advantage of Cellular Potts Modeling to simulate cellular interactions and combine them with differential equations to model dynamic cell signaling patterns. These models are advantageous over other approaches because they retain spatial information about each cell while remaining computationally efficient at larger scales. Users of this package define three key inputs to create valid model definitions: a 2- or 3-dimensional space, a table describing the cells to be positioned in that space, and a list of model penalties that dictate cell behaviors. Models can then be evolved over time to collect statistics, simulated repeatedly to investigate how changing a specific property impacts cellular behavior, and visualized using any of the available plotting libraries in Julia. AVAILABILITY AND IMPLEMENTATION: The CellularPotts.jl package is released under the MIT license and is available at https://github.com/RobertGregg/CellularPotts.jl. An archived version of the code (v0.3.2) at time of submission can also be found at https://doi.org/10.5281/zenodo.10407783.

LEF1 isoforms regulate cellular senescence and aging.

The study of aging and its mechanisms, such as cellular senescence, has provided valuable insights into age-related pathologies, thus contributing to their prevention and treatment. The current abundance of high-throughput data combined with the surge of robust analysis algorithms has facilitated novel ways of identifying underlying pathways that may drive these pathologies. For the purpose of identifying key regulators of lung aging, we performed comparative analyses of transcriptional profiles of aged versus young human subjects and mice, focusing on the common age-related changes in the transcriptional regulation in lung macrophages, T cells, and B immune cells. Importantly, we validated our findings in cell culture assays and human lung samples. Our analysis identified lymphoid enhancer binding factor 1 (LEF1) as an important age-associated regulator of gene expression in all three cell types across different tissues and species. Follow-up experiments showed that the differential expression of long and short LEF1 isoforms is a key regulatory mechanism of cellular senescence. Further examination of lung tissue from patients with idiopathic pulmonary fibrosis, an age-related disease with strong ties to cellular senescence, revealed a stark dysregulation of LEF1. Collectively, our results suggest that LEF1 is a key factor of aging, and its differential regulation is associated with human and murine cellular senescence.

Prognostic Insights from Longitudinal Multicompartment Study of Host-Microbiota Interactions in Critically Ill Patients.

Critical illness can disrupt the composition and function of the microbiome, yet comprehensive longitudinal studies are lacking. We conducted a longitudinal analysis of oral, lung, and gut microbiota in a large cohort of 479 mechanically ventilated patients with acute respiratory failure. Progressive dysbiosis emerged in all three body compartments, characterized by reduced alpha diversity, depletion of obligate anaerobe bacteria, and pathogen enrichment. Clinical variables, including chronic obstructive pulmonary disease, immunosuppression, and antibiotic exposure, shaped dysbiosis. Notably, of the three body compartments, unsupervised clusters of lung microbiota diversity and composition independently predicted survival, transcending clinical predictors, organ dysfunction severity, and host-response sub-phenotypes. These independent associations of lung microbiota may serve as valuable biomarkers for prognostication and treatment decisions in critically ill patients. Insights into the dynamics of the microbiome during critical illness highlight the potential for microbiota-targeted interventions in precision medicine.

Defining trophoblast injury patterns in the transcriptomes of dysfunctional placentas.

Trophoblast injury is central to clinically relevant placenta dysfunction. We hypothesized that the mRNA of primary human trophoblasts, exposed to distinct injuries in vitro, capture transcriptome patterns of placental biopsies obtained from common obstetrical syndromes. We deployed a CIBERSORTx deconvolution method to correlate trophoblastic RNAseq-based expression matrices with the transcriptome of omics-defined placental dysfunction patterns in vivo. We found distinct trophoblast injury patterns in placental biopsies from women with fetal growth restriction and a hypertensive disorder, or in biopsies clustered by their omics analysis. Our RNAseq data are useful for defining the contribution of trophoblast injuries to placental dysfunction syndromes.

Prognostic Insights from Longitudinal Multicompartment Study of Host-Microbiota Interactions in Critically Ill Patients.

Critical illness can disrupt the composition and function of the microbiome, yet comprehensive longitudinal studies are lacking. We conducted a longitudinal analysis of oral, lung, and gut microbiota in a large cohort of 479 mechanically ventilated patients with acute respiratory failure. Progressive dysbiosis emerged in all three body compartments, characterized by reduced alpha diversity, depletion of obligate anaerobe bacteria, and pathogen enrichment. Clinical variables, including chronic obstructive pulmonary disease, immunosuppression, and antibiotic exposure, shaped dysbiosis. Notably, of the three body compartments, unsupervised clusters of lung microbiota diversity and composition independently predicted survival, transcending clinical predictors, organ dysfunction severity, and host-response sub-phenotypes. These independent associations of lung microbiota may serve as valuable biomarkers for prognostication and treatment decisions in critically ill patients. Insights into the dynamics of the microbiome during critical illness highlight the potential for microbiota-targeted interventions in precision medicine.

Integrated unbiased multiomics defines disease-independent placental clusters in common obstetrical syndromes.

BACKGROUND: Placental dysfunction, a root cause of common syndromes affecting human pregnancy, such as preeclampsia (PE), fetal growth restriction (FGR), and spontaneous preterm delivery (sPTD), remains poorly defined. These common, yet clinically disparate obstetrical syndromes share similar placental histopathologic patterns, while individuals within each syndrome present distinct molecular changes, challenging our understanding and hindering our ability to prevent and treat these syndromes. METHODS: Using our extensive biobank, we identified women with severe PE (n = 75), FGR (n = 40), FGR with a hypertensive disorder (FGR + HDP; n = 33), sPTD (n = 72), and two uncomplicated control groups, term (n = 113), and preterm without PE, FGR, or sPTD (n = 16). We used placental biopsies for transcriptomics, proteomics, metabolomics data, and histological evaluation. After conventional pairwise comparison, we deployed an unbiased, AI-based similarity network fusion (SNF) to integrate the datatypes and identify omics-defined placental clusters. We used Bayesian model selection to compare the association between the histopathological features and disease conditions vs SNF clusters. RESULTS: Pairwise, disease-based comparisons exhibited relatively few differences, likely reflecting the heterogeneity of the clinical syndromes. Therefore, we deployed the unbiased, omics-based SNF method. Our analysis resulted in four distinct clusters, which were mostly dominated by a specific syndrome. Notably, the cluster dominated by early-onset PE exhibited strong placental dysfunction patterns, with weaker injury patterns in the cluster dominated by sPTD. The SNF-defined clusters exhibited better correlation with the histopathology than the predefined disease groups. CONCLUSIONS: Our results demonstrate that integrated omics-based SNF distinctively reclassifies placental dysfunction patterns underlying the common obstetrical syndromes, improves our understanding of the pathological processes, and could promote a search for more personalized interventions.

Integrated BATF transcriptional network regulates suppressive intratumoral regulatory T cells.

Human regulatory T cells (Tregs) are crucial regulators of tissue repair, autoimmune diseases, and cancer. However, it is challenging to inhibit the suppressive function of Tregs for cancer therapy without affecting immune homeostasis. Identifying pathways that may distinguish tumor-restricted Tregs is important, yet the transcriptional programs that control intratumoral Treg gene expression, and that are distinct from Tregs in healthy tissues, remain largely unknown. We profiled single-cell transcriptomes of CD4+ T cells in tumors and peripheral blood from patients with head and neck squamous cell carcinomas (HNSCC) and those in nontumor tonsil tissues and peripheral blood from healthy donors. We identified a subpopulation of activated Tregs expressing multiple tumor necrosis factor receptor (TNFR) genes (TNFR+ Tregs) that is highly enriched in the tumor microenvironment (TME) compared with nontumor tissue and the periphery. TNFR+ Tregs are associated with worse prognosis in HNSCC and across multiple solid tumor types. Mechanistically, the transcription factor BATF is a central component of a gene regulatory network that governs key aspects of TNFR+ Tregs. CRISPR-Cas9-mediated BATF knockout in human activated Tregs in conjunction with bulk RNA sequencing, immunophenotyping, and in vitro functional assays corroborated the central role of BATF in limiting excessive activation and promoting the survival of human activated Tregs. Last, we identified a suite of surface molecules reflective of the BATF-driven transcriptional network on intratumoral Tregs in patients with HNSCC. These findings uncover a primary transcriptional regulator of highly suppressive intratumoral Tregs, highlighting potential opportunities for therapeutic intervention in cancer without affecting immune homeostasis.