Foxm1 regulates cardiomyocyte proliferation in adult zebrafish after cardiac injury.

The regenerative capacity of the mammalian heart is poor, with one potential reason being that adult cardiomyocytes cannot proliferate at sufficient levels to replace lost tissue. During development and neonatal stages, cardiomyocytes can successfully divide under injury conditions; however, as these cells mature their ability to proliferate is lost. Therefore, understanding the regulatory programs that can induce post-mitotic cardiomyocytes into a proliferative state is essential to enhance cardiac regeneration. Here, we report that the forkhead transcription factor Foxm1 is required for cardiomyocyte proliferation after injury through transcriptional regulation of cell cycle genes. Transcriptomic analysis of injured zebrafish hearts revealed that foxm1 expression is increased in border zone cardiomyocytes. Decreased cardiomyocyte proliferation and expression of cell cycle genes in foxm1 mutant hearts was observed, suggesting it is required for cell cycle checkpoints. Subsequent analysis of a candidate Foxm1 target gene, cenpf, revealed that this microtubule and kinetochore binding protein is also required for cardiac regeneration. Moreover, cenpf mutants show increased cardiomyocyte binucleation. Thus, foxm1 and cenpf are required for cardiomyocytes to complete mitosis during zebrafish cardiac regeneration.

An Insulin-Responsive Sensor in the SIRT1 Disordered Region Binds DBC1 and PACS-2 to Control Enzyme Activity.

Current models of SIRT1 enzymatic regulation primarily consider the effects of fluctuating levels of its co-substrate NAD+, which binds to the stably folded catalytic domain. By contrast, the roles of the sizeable disordered N- and C-terminal regions of SIRT1 are largely unexplored. Here we identify an insulin-responsive sensor in the SIRT1 N-terminal region (NTR), comprising an acidic cluster (AC) and a 3-helix bundle (3HB), controlling deacetylase activity. The allosteric assistor DBC1 removes a distal N-terminal shield from the 3-helix bundle, permitting PACS-2 to engage the acidic cluster and the transiently exposed helix 3 of the 3-helix bundle, disrupting its structure and inhibiting catalysis. The SIRT1 activator (STAC) SRT1720 binds and stabilizes the 3-helix bundle, protecting SIRT1 from inhibition by PACS-2. Identification of the SIRT1 insulin-responsive sensor and its engagement by the DBC1 and PACS-2 regulatory hub provides important insight into the roles of disordered regions in enzyme regulation and the mode by which STACs promote metabolic fitness.

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.

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.

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.

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.

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.

An intrinsically interpretable neural network architecture for sequence-to-function learning.

MOTIVATION: Sequence-based deep learning approaches have been shown to predict a multitude of functional genomic readouts, including regions of open chromatin and RNA expression of genes. However, a major limitation of current methods is that model interpretation relies on computationally demanding post hoc analyses, and even then, one can often not explain the internal mechanics of highly parameterized models. Here, we introduce a deep learning architecture called totally interpretable sequence-to-function model (tiSFM). tiSFM improves upon the performance of standard multilayer convolutional models while using fewer parameters. Additionally, while tiSFM is itself technically a multilayer neural network, internal model parameters are intrinsically interpretable in terms of relevant sequence motifs. RESULTS: We analyze published open chromatin measurements across hematopoietic lineage cell-types and demonstrate that tiSFM outperforms a state-of-the-art convolutional neural network model custom-tailored to this dataset. We also show that it correctly identifies context-specific activities of transcription factors with known roles in hematopoietic differentiation, including Pax5 and Ebf1 for B-cells, and Rorc for innate lymphoid cells. tiSFM's model parameters have biologically meaningful interpretations, and we show the utility of our approach on a complex task of predicting the change in epigenetic state as a function of developmental transition. AVAILABILITY AND IMPLEMENTATION: The source code, including scripts for the analysis of key findings, can be found at https://github.com/boooooogey/ATAConv, implemented in Python.

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.

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.

Distinct COPD subtypes in former smokers revealed by gene network perturbation analysis.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) varies significantly in symptomatic and physiologic presentation. Identifying disease subtypes from molecular data, collected from easily accessible blood samples, can help stratify patients and guide disease management and treatment. METHODS: Blood gene expression measured by RNA-sequencing in the COPDGene Study was analyzed using a network perturbation analysis method. Each COPD sample was compared against a learned reference gene network to determine the part that is deregulated. Gene deregulation values were used to cluster the disease samples. RESULTS: The discovery set included 617 former smokers from COPDGene. Four distinct gene network subtypes are identified with significant differences in symptoms, exercise capacity and mortality. These clusters do not necessarily correspond with the levels of lung function impairment and are independently validated in two external cohorts: 769 former smokers from COPDGene and 431 former smokers in the Multi-Ethnic Study of Atherosclerosis (MESA). Additionally, we identify several genes that are significantly deregulated across these subtypes, including DSP and GSTM1, which have been previously associated with COPD through genome-wide association study (GWAS). CONCLUSIONS: The identified subtypes differ in mortality and in their clinical and functional characteristics, underlining the need for multi-dimensional assessment potentially supplemented by selected markers of gene expression. The subtypes were consistent across cohorts and could be used for new patient stratification and disease prognosis.

CausalMGM: an interactive web-based causal discovery tool.

High-throughput sequencing and the availability of large online data repositories (e.g. The Cancer Genome Atlas and Trans-Omics for Precision Medicine) have the potential to revolutionize systems biology by enabling researchers to study interactions between data from different modalities (i.e. genetic, genomic, clinical, behavioral, etc.). Currently, data mining and statistical approaches are confined to identifying correlates in these datasets, but researchers are often interested in identifying cause-and-effect relationships. Causal discovery methods were developed to infer such cause-and-effect relationships from observational data. Though these algorithms have had demonstrated successes in several biomedical applications, they are difficult to use for non-experts. So, there is a need for web-based tools to make causal discovery methods accessible. Here, we present CausalMGM (http://causalmgm.org/), the first web-based causal discovery tool that enables researchers to find cause-and-effect relationships from observational data. Web-based CausalMGM consists of three data analysis tools: (i) feature selection and clustering; (ii) automated identification of cause-and-effect relationships via a graphical model; and (iii) interactive visualization of the learned causal (directed) graph. We demonstrate how CausalMGM enables an end-to-end exploratory analysis of biomedical datasets, giving researchers a clearer picture of its capabilities.

Topographic heterogeneity of lung microbiota in end-stage idiopathic pulmonary fibrosis: the Microbiome in Lung Explants-2 (MiLEs-2) study.

BACKGROUND: Lung microbiota profiles in patients with early idiopathic pulmonary fibrosis (IPF) have been associated with disease progression; however, the topographic heterogeneity of lung microbiota and their roles in advanced IPF are unknown. METHODS: We performed a retrospective, case-control study of explanted lung tissue obtained at the time of lung transplantation or rapid autopsy from patients with IPF and other chronic lung diseases (connective tissue disease-associated interstitial lung disease (CTD-ILD), cystic fibrosis (CF), COPD and donor lungs unsuitable for transplant from Center for Organ Recovery and Education (CORE)). We sampled subpleural tissue and airway-based specimens (bronchial washings and airway tissue) and quantified bacterial load and profiled communities by amplification and sequencing of the 16S rRNA gene. FINDINGS: Explants from 62 patients with IPF, 15 patients with CTD-ILD, 20 patients with CF, 20 patients with COPD and 20 CORE patients were included. Airway-based samples had higher bacterial load compared with distal parenchymal tissue. IPF basilar tissue had much lower bacterial load compared with CF and CORE lungs (p<0.001). No microbial community differences were found between parenchymal tissue samples from different IPF lobes. Dirichlet multinomial models revealed an IPF cluster (29%) with distinct composition, high bacterial load and low alpha diversity, exhibiting higher odds for acute exacerbation or death. INTERPRETATION: IPF explants had low biomass in the distal parenchyma of all three lobes with higher bacterial load in the airways. The discovery of a distinct subgroup of patients with IPF with higher bacterial load and worse clinical outcomes supports investigation of personalised medicine approaches for microbiome-targeted interventions.

Transcriptomics of bronchoalveolar lavage cells identifies new molecular endotypes of sarcoidosis.

BACKGROUND: Sarcoidosis is a multisystem granulomatous disease of unknown origin with a variable and often unpredictable course and pattern of organ involvement. In this study we sought to identify specific bronchoalveolar lavage (BAL) cell gene expression patterns indicative of distinct disease phenotypic traits. METHODS: RNA sequencing by Ion Torrent Proton was performed on BAL cells obtained from 215 well-characterised patients with pulmonary sarcoidosis enrolled in the multicentre Genomic Research in Alpha-1 Antitrypsin Deficiency and Sarcoidosis (GRADS) study. Weighted gene co-expression network analysis and nonparametric statistics were used to analyse genome-wide BAL transcriptome. Validation of results was performed using a microarray expression dataset of an independent sarcoidosis cohort (Freiburg, Germany; n=50). RESULTS: Our supervised analysis found associations between distinct transcriptional programmes and major pulmonary phenotypic manifestations of sarcoidosis including T-helper type 1 (Th1) and Th17 pathways associated with hilar lymphadenopathy, transforming growth factor-ß1 (TGFB1) and mechanistic target of rapamycin (MTOR) signalling with parenchymal involvement, and interleukin (IL)-7 and IL-2 with airway involvement. Our unsupervised analysis revealed gene modules that uncovered four potential sarcoidosis endotypes including hilar lymphadenopathy with increased acute T-cell immune response; extraocular organ involvement with PI3K activation pathways; chronic and multiorgan disease with increased immune response pathways; and multiorgan involvement, with increased IL-1 and IL-18 immune and inflammatory responses. We validated the occurrence of these endotypes using gene expression, pulmonary function tests and cell differentials from Freiburg. CONCLUSION: Taken together, our results identify BAL gene expression programmes that characterise major pulmonary sarcoidosis phenotypes and suggest the presence of distinct disease molecular endotypes.

Causal network perturbations for instance-specific analysis of single cell and disease samples.

MOTIVATION: Complex diseases involve perturbation in multiple pathways and a major challenge in clinical genomics is characterizing pathway perturbations in individual samples. This can lead to patient-specific identification of the underlying mechanism of disease thereby improving diagnosis and personalizing treatment. Existing methods rely on external databases to quantify pathway activity scores. This ignores the data dependencies and that pathways are incomplete or condition-specific. RESULTS: ssNPA is a new approach for subtyping samples based on deregulation of their gene networks. ssNPA learns a causal graph directly from control data. Sample-specific network neighborhood deregulation is quantified via the error incurred in predicting the expression of each gene from its Markov blanket. We evaluate the performance of ssNPA on liver development single-cell RNA-seq data, where the correct cell timing is recovered; and two TCGA datasets, where ssNPA patient clusters have significant survival differences. In all analyses ssNPA consistently outperforms alternative methods, highlighting the advantage of network-based approaches. AVAILABILITY AND IMPLEMENTATION: http://www.benoslab.pitt.edu/Software/ssnpa/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Neurological Complications Acquired During Pediatric Critical Illness: Exploratory "Mixed Graphical Modeling" Analysis Using Serum Biomarker Levels.

OBJECTIVES: Neurologic complications, consisting of the acute development of a neurologic disorder, that is, not present at admission but develops during the course of illness, can be difficult to detect in the PICU due to sedation, neuromuscular blockade, and young age. We evaluated the direct relationships of serum biomarkers and clinical variables to the development of neurologic complications. Analysis was performed using mixed graphical models, a machine learning approach that allows inference of cause-effect associations from continuous and discrete data. DESIGN: Secondary analysis of a previous prospective observational study. SETTING: PICU, single quaternary-care center. PATIENTS: Individuals admitted to the PICU, younger than18 years old, with intravascular access via an indwelling catheter. INTERVENTIONS: None. MEASUREMENTS: About 101 patients were included in this analysis. Serum (days 1-7) was analyzed for glial fibrillary acidic protein, ubiquitin C-terminal hydrolase-L1, and alpha-II spectrin breakdown product 150 utilizing enzyme-linked immunosorbent assays. Serum levels of neuron-specific enolase, myelin basic protein, and S100 calcium binding protein B used in these models were reported previously. Demographic data, use of selected clinical therapies, lengths of stay, and ancillary neurologic testing (head CT, brain MRI, and electroencephalogram) results were recorded. The Mixed Graphical Model-Fast-Causal Inference-Maximum algorithm was applied to the dataset. MAIN RESULTS: About 13 of 101 patients developed a neurologic complication during their critical illness. The mixed graphical model identified peak levels of the neuronal biomarker neuron-specific enolase and ubiquitin C-terminal hydrolase-L1, and the astrocyte biomarker glial fibrillary acidic protein to be the direct causal determinants for the development of a neurologic complication; in contrast, clinical variables including age, sex, length of stay, and primary neurologic diagnosis were not direct causal determinants. CONCLUSIONS: Graphical models that include biomarkers in addition to clinical data are promising methods to evaluate direct relationships in the development of neurologic complications in critically ill children. Future work is required to validate and refine these models further, to determine if they can be used to predict which patients are at risk for/or with early neurologic complications.

Respiratory Tract Dysbiosis Is Associated with Worse Outcomes in Mechanically Ventilated Patients.

Rationale: Host inflammatory responses have been strongly associated with adverse outcomes in critically ill patients, but the biologic underpinnings of such heterogeneous responses have not been defined.Objectives: We examined whether respiratory tract microbiome profiles are associated with host inflammation and clinical outcomes of acute respiratory failure.Methods: We collected oral swabs, endotracheal aspirates (ETAs), and plasma samples from mechanically ventilated patients. We performed 16S ribosomal RNA gene sequencing to characterize upper and lower respiratory tract microbiota and classified patients into host-response subphenotypes on the basis of clinical variables and plasma biomarkers of innate immunity and inflammation. We derived diversity metrics and composition clusters with Dirichlet multinomial models and examined our data for associations with subphenotypes and clinical outcomes.Measurements and Main Results: Oral and ETA microbial communities from 301 mechanically ventilated subjects had substantial heterogeneity in α and ß diversity. Dirichlet multinomial models revealed a cluster with low α diversity and enrichment for pathogens (e.g., high Staphylococcus or Pseudomonadaceae relative abundance) in 35% of ETA samples, associated with a hyperinflammatory subphenotype, worse 30-day survival, and longer time to liberation from mechanical ventilation (adjusted P < 0.05), compared with patients with higher α diversity and relative abundance of typical oral microbiota. Patients with evidence of dysbiosis (low α diversity and low relative abundance of "protective" oral-origin commensal bacteria) in both oral and ETA samples (17%, combined dysbiosis) had significantly worse 30-day survival and longer time to liberation from mechanical ventilation than patients without dysbiosis (55%; adjusted P < 0.05).Conclusions: Respiratory tract dysbiosis may represent an important, modifiable contributor to patient-level heterogeneity in systemic inflammatory responses and clinical outcomes.

An improvement of ComiR algorithm for microRNA target prediction by exploiting coding region sequences of mRNAs.

MicroRNA are small non-coding RNAs that post-transcriptionally regulate the expression levels of messenger RNAs. MicroRNA regulation activity depends on the recognition of binding sites located on mRNA molecules. ComiR is a web tool realized to predict the targets of a set of microRNAs, starting from their expression profile. ComiR was trained with the information regarding binding sites in the 3'utr region, by using a reliable dataset containing the targets of endogenously expressed microRNA in D. melanogaster S2 cells. This dataset was obtained by comparing the results from two different experimental approaches, i.e., inhibition, and immunoprecipitation of the AGO1 protein--a component of the microRNA induced silencing complex.In this work, we tested whether including coding region binding sites in ComiR algorithm improves the performance of the tool in predicting microRNA targets. We focused the analysis on the D. melanogaster species and updated the ComiR underlying database with the currently available releases of mRNA and microRNA sequences. As a result, we find that ComiR algorithm trained with the information related to the coding regions is more efficient in predicting the microRNA targets, with respect to the algorithm trained with 3'utr information. On the other hand, we show that 3'utr based predictions can be seen as complementary to the coding region based predictions, which suggests that both predictions, from 3'utr and coding regions, should be considered in comprehensive analysis.Furthermore, we observed that the lists of targets obtained by analyzing data from one experimental approach only, that is, inhibition or immunoprecipitation of AGO1, are not reliable enough to test the performance of our microRNA target prediction algorithm. Further analysis will be conducted to investigate the effectiveness of the tool with data from other species, provided that validated datasets, as obtained from the comparison of RISC proteins inhibition and immunoprecipitation experiments, will be available for the same samples. Finally, we propose to upgrade the existing ComiR web-tool by including the coding region based trained model, available together with the 3'utr based one.

Mixed graphical models for integrative causal analysis with application to chronic lung disease diagnosis and prognosis.

MOTIVATION: Integration of data from different modalities is a necessary step for multi-scale data analysis in many fields, including biomedical research and systems biology. Directed graphical models offer an attractive tool for this problem because they can represent both the complex, multivariate probability distributions and the causal pathways influencing the system. Graphical models learned from biomedical data can be used for classification, biomarker selection and functional analysis, while revealing the underlying network structure and thus allowing for arbitrary likelihood queries over the data. RESULTS: In this paper, we present and test new methods for finding directed graphs over mixed data types (continuous and discrete variables). We used this new algorithm, CausalMGM, to identify variables directly linked to disease diagnosis and progression in various multi-modal datasets, including clinical datasets from chronic obstructive pulmonary disease (COPD). COPD is the third leading cause of death and a major cause of disability and thus determining the factors that cause longitudinal lung function decline is very important. Applied on a COPD dataset, mixed graphical models were able to confirm and extend previously described causal effects and provide new insights on the factors that potentially affect the longitudinal lung function decline of COPD patients. AVAILABILITY AND IMPLEMENTATION: The CausalMGM package is available on http://www.causalmgm.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Feasibility of lung cancer prediction from low-dose CT scan and smoking factors using causal models.

INTRODUCTION: Low-dose CT (LDCT) is currently used in lung cancer screening of high-risk populations for early lung cancer diagnosis. However, 96% of individuals with detected nodules are false positives. METHODS: In order to develop an efficient early lung cancer predictor from clinical, demographic and LDCT features, we studied a total of 218 subjects with lung cancer or benign nodules. Probabilistic graphical models (PGMs) were used to integrate demographics, clinical data and LDCT features from 92 subjects (training cohort) from the Pittsburgh Lung Screening Study cohort. RESULTS: Learnt PGMs identified three variables directly (causally) linked to malignant nodules and the largest benign nodule and used them to build the Lung Cancer Causal Model (LCCM), which was validated in a separate cohort of 126 subjects. Nodule and vessel numbers and years since the subject quit smoking were sufficient to discriminate malignant from benign nodules. Comparison with existing predictors in the training and validation cohorts showed that (1) incorporating LDCT scan features greatly enhances predictive accuracy; and (2) LCCM improves cancer detection over existing methods, including the Brock parsimonious model (p<0.001). Notably, the number of surrounding vessels, a feature not previously used in predictive models, significantly improves predictive efficiency. Based on the validation cohort results, LCCM is able to identify 30% of the benign nodules without risk of misclassifying cancer nodules. DISCUSSION: LCCM shows promise as a lung cancer predictor as it is significantly improved over existing models. Validated in a larger, prospective study, it may help reduce unnecessary follow-up visits and procedures.