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.

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.

Transcriptional changes of the aging lung.

Aging is a natural process associated with declined organ function and higher susceptibility to developing chronic diseases. A systemic single-cell type-based study provides a unique opportunity to understand the mechanisms behind age-related pathologies. Here, we use single-cell gene expression analysis comparing healthy young and aged human lungs from nonsmoker donors to investigate age-related transcriptional changes. Our data suggest that aging has a heterogenous effect on lung cells, as some populations are more transcriptionally dynamic while others remain stable in aged individuals. We found that monocytes and alveolar macrophages were the most transcriptionally affected populations. These changes were related to inflammation and regulation of the immune response. Additionally, we calculated the LungAge score, which reveals the diversity of lung cell types during aging. Changes in DNA damage repair, fatty acid metabolism, and inflammation are essential for age prediction. Finally, we quantified the senescence score in aged lungs and found that the more biased cells toward senescence are immune and progenitor cells. Our study provides a comprehensive and systemic analysis of the molecular signatures of lung aging. Our LungAge signature can be used to predict molecular signatures of physiological aging and to detect common signatures of age-related lung diseases.

LEF1 isoforms regulate cellular senescence and aging.

Background: 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. Methods: With the focus on 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. Results: We 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 (IPF), an age-related disease with strong ties to cellular senescence, we demonstrated a stark dysregulation of LEF1. Conclusions: Collectively, our results suggest that the LEF1 is a key factor of aging, and its differential regulation is associated with human and murine cellular senescence.

Early events marking lung fibroblast transition to profibrotic state in idiopathic pulmonary fibrosis.

BACKGROUND: Idiopathic Pulmonary Fibrosis (IPF) is an age-associated progressive lung disease with accumulation of scar tissue impairing gas exchange. Previous high-throughput studies elucidated the role of cellular heterogeneity and molecular pathways in advanced disease. However, critical pathogenic pathways occurring in the transition of fibroblasts from normal to profibrotic have been largely overlooked. METHODS: We used single cell transcriptomics (scRNA-seq) from lungs of healthy controls and IPF patients (lower and upper lobes). We identified fibroblast subclusters, genes and pathways associated with early disease. Immunofluorescence assays validated the role of MOXD1 early in fibrosis. RESULTS: We identified four distinct fibroblast subgroups, including one marking the normal-to-profibrotic state transition. Our results show for the first time that global downregulation of ribosomal proteins and significant upregulation of the majority of copper-binding proteins, including MOXD1, mark the IPF transition. We find no significant differences in gene expression in IPF upper and lower lobe samples, which were selected to have low and high degree of fibrosis, respectively. CONCLUSIONS: Early events during IPF onset in fibroblasts include dysregulation of ribosomal and copper-binding proteins. Fibroblasts in early stage IPF may have already acquired a profibrotic phenotype while hallmarks of advanced disease, including fibroblast foci and honeycomb formation, are still not evident. The new transitional fibroblasts we discover could prove very important for studying the role of fibroblast plasticity in disease progression and help develop early diagnosis tools and therapeutic interventions targeting earlier disease states.

Causal Discovery in High-dimensional, Multicollinear Datasets.

As the cost of high-throughput genomic sequencing technology declines, its application in clinical research becomes increasingly popular. The collected datasets often contain tens or hundreds of thousands of biological features that need to be mined to extract meaningful information. One area of particular interest is discovering underlying causal mechanisms of disease outcomes. Over the past few decades, causal discovery algorithms have been developed and expanded to infer such relationships. However, these algorithms suffer from the curse of dimensionality and multicollinearity. A recently introduced, non-orthogonal, general empirical Bayes approach to matrix factorization has been demonstrated to successfully infer latent factors with interpretable structures from observed variables. We hypothesize that applying this strategy to causal discovery algorithms can solve both the high dimensionality and collinearity problems, inherent to most biomedical datasets. We evaluate this strategy on simulated data and apply it to two real-world datasets. In a breast cancer dataset, we identified important survival-associated latent factors and biologically meaningful enriched pathways within factors related to important clinical features. In a SARS-CoV-2 dataset, we were able to predict whether a patient (1) had Covid-19 and (2) would enter the ICU. Furthermore, we were able to associate factors with known Covid-19 related biological pathways.

Myofibroblast transcriptome indicates SFRP2hi fibroblast progenitors in systemic sclerosis skin.

Skin and lung fibrosis in systemic sclerosis (SSc) is driven by myofibroblasts, alpha-smooth muscle actin expressing cells. The number of myofibroblasts in SSc skin correlates with the modified Rodnan skin score, the most widely used clinical measure of skin disease severity. Murine fibrosis models indicate that myofibroblasts can arise from a variety of different cell types, but their origin in SSc skin has remained uncertain. Utilizing single cell RNA-sequencing, we define different dermal fibroblast populations and transcriptome changes, comparing SSc to healthy dermal fibroblasts. Here, we show that SSc dermal myofibroblasts arise in two steps from an SFRP2hi/DPP4-expressing progenitor fibroblast population. In the first step, SSc fibroblasts show globally upregulated expression of transcriptome markers, such as PRSS23 and THBS1. A subset of these cells shows markers indicating that they are proliferating. Only a fraction of SFRP2hi SSc fibroblasts differentiate into myofibroblasts, as shown by expression of additional markers, SFRP4 and FNDC1. Bioinformatics analysis of the SSc fibroblast transcriptomes implicated upstream transcription factors, including FOSL2, RUNX1, STAT1, FOXP1, IRF7 and CREB3L1, as well as SMAD3, driving SSc myofibroblast differentiation.