Search

Abstract Spatial barcoding-based transcriptomic (ST) data require deconvolution for cellular-level downstream analysis. Here we present SDePER, a hybrid machine learning and regression method to deconvolve ST data using reference single-cell RNA sequencing (scRNA-seq) data. SDePER tackles platform effects between ST and scRNA-seq data, ensuring a linear relationship between them while addressing sparsity and spatial correlations in cell types across capture spots. SDePER estimates cell-type proportions, enabling enhanced resolution tissue mapping by imputing cell-type compositions and gene expressions at unmeasured locations. Applications to simulated data and four real datasets showed SDePER’s superior accuracy and robustness over existing methods.

Summary: Matrix stiffening by lysyl oxidase-like 2 (LOXL2)-mediated collagen cross-linking is proposed as a core feedforward mechanism that promotes fibrogenesis. Failure in clinical trials of simtuzumab (the humanized version of AB0023, a monoclonal antibody against human LOXL2) suggested that targeting LOXL2 may not have disease relevance; however, target engagement was not directly evaluated. We compare the spatial transcriptome of active human lung fibrogenesis sites with different human cell culture models to identify a disease-relevant model. Within the selected model, we then evaluate AB0023, identifying that it does not inhibit collagen cross-linking or reduce tissue stiffness, nor does it inhibit LOXL2 catalytic activity. In contrast, it does potently inhibit angiogenesis consistent with an alternative, non-enzymatic mechanism of action. Thus, AB0023 is anti-angiogenic but does not inhibit LOXL2 catalytic activity, collagen cross-linking, or tissue stiffening. These findings have implications for the interpretation of the lack of efficacy of simtuzumab in clinical trials of fibrotic diseases.

Objectives Mycophenolate mofetil (MMF) and azathioprine (AZA) are immunomodulatory treatments in interstitial lung disease (ILD). This systematic review aimed to evaluate the efficacy of MMF or AZA on pulmonary function in ILD.Design Population included any ILD diagnosis, intervention included MMF or AZA treatment, outcome was delta change from baseline in per cent predicted forced vital capacity (%FVC) and gas transfer (diffusion lung capacity of carbon monoxide, %DLco). The primary endpoint compared outcomes relative to placebo comparator, the secondary endpoint assessed outcomes in treated groups only.Eligibility criteria Randomised controlled trials (RCTs) and prospective observational studies were included. No language restrictions were applied. Retrospective studies and studies with high-dose concomitant steroids were excluded.Data synthesis The systematic search was performed on 9 May. Meta-analyses according to drug and outcome were specified with random effects, I2 evaluated heterogeneity and Grading of Recommendations, Assessment, Development and Evaluation evaluated certainty of evidence. Primary endpoint analysis was restricted to RCT design, secondary endpoint included subgroup analysis according to prospective observational or RCT design.Results A total of 2831 publications were screened, 12 were suitable for quantitative synthesis. Three MMF RCTs were included with no significant effect on the primary endpoints (%FVC 2.94, 95% CI −4.00 to 9.88, I2=79.3%; %DLco −2.03, 95% CI −4.38 to 0.32, I2=0.0%). An overall 2.03% change from baseline in %FVC (95% CI 0.65 to 3.42, I2=0.0%) was observed in MMF, and RCT subgroup summary estimated a 4.42% change from baseline in %DLCO (95% CI 2.05 to 6.79, I2=0.0%). AZA studies were limited. All estimates were considered very low certainty evidence.Conclusions There were limited RCTs of MMF or AZA and their benefit in ILD was of very low certainty. MMF may support preservation of pulmonary function, yet confidence in the effect was weak. To support high certainty evidence, RCTs should be designed to directly assess MMF efficacy in ILD.PROSPERO registration number CRD42023423223.

IntroductionChronic rejection is a major complication post-transplantation. Within lung transplantation, chronic rejection was considered as airway centred. Chronic Lung Allograft Dysfunction (CLAD), defined to cover all late chronic complications, makes it more difficult to understand chronic rejection from an immunological perspective. This study investigated the true nature, timing and location of chronic rejection as a whole, within mouse lung transplantation.Methods40 mice underwent an orthotopic left lung transplantation, were sacrificed at day 70 and evaluated by histology and in vivo µCT. For timing and location of rejection, extra grafts were sacrificed at day 7, 35, 56 and investigated by ex vivo µCT or single cell RNA (scRNA) profiling.ResultsChronic rejection originated as innate inflammation around small arteries evolving toward adaptive organization with subsequent end-arterial fibrosis and obliterans. Subsequently, venous and pleural infiltration appeared, followed by airway related bronchiolar folding and rarely bronchiolitis obliterans was observed. Ex vivo µCT and scRNA profiling validated the time, location and sequence of events with endothelial destruction and activation as primary onset.ConclusionAgainst the current belief, chronic rejection in lung transplantation may start as an arterial response, followed by responses in venules, pleura, and, only in the late stage, bronchioles, as may be seen in some but not all patients with CLAD.

Alveolarization ensures sufficient lung surface area for gas exchange, and during bulk alveolarization in mice (postnatal day [P] 4.5–14.5), alpha-smooth muscle actin (SMA)+ myofibroblasts accumulate, secrete elastin, and lay down alveolar septum. Herein, we delineate the dynamics of the lineage of early postnatal SMA+ myofibroblasts during and after bulk alveolarization and in response to lung injury. SMA+ lung myofibroblasts first appear at ∼ P2.5 and proliferate robustly. Lineage tracing shows that, at P14.5 and over the next few days, the vast majority of SMA+ myofibroblasts downregulate smooth muscle cell markers and undergo apoptosis. Of note, ∼8% of these dedifferentiated cells and another ∼1% of SMA+ myofibroblasts persist to adulthood. Single cell RNA sequencing analysis of the persistent SMA− cells and SMA+ myofibroblasts in the adult lung reveals distinct gene expression profiles. For instance, dedifferentiated SMA− cells exhibit higher levels of tissue remodeling genes. Most interestingly, these dedifferentiated early postnatal myofibroblasts re-express SMA upon exposure of the adult lung to hypoxia or the pro-fibrotic drug bleomycin. However, unlike during alveolarization, these cells that re-express SMA do not proliferate with hypoxia. In sum, dedifferentiated early postnatal myofibroblasts are a previously undescribed cell type in the adult lung and redifferentiate in response to injury.

Abstract The activation and accumulation of lung fibroblasts resulting in aberrant deposition of extracellular matrix components, is a pathogenic hallmark of Idiopathic Pulmonary Fibrosis, a lethal and incurable disease. In this report, increased expression of TKS5, a scaffold protein essential for the formation of podosomes, was detected in the lung tissue of Idiopathic Pulmonary Fibrosis patients and bleomycin-treated mice. Τhe profibrotic milieu is found to induce TKS5 expression and the formation of prominent podosome rosettes in lung fibroblasts, that are retained ex vivo, culminating in increased extracellular matrix invasion. Tks5 +/- mice are found resistant to bleomycin-induced pulmonary fibrosis, largely attributed to diminished podosome formation in fibroblasts and decreased extracellular matrix invasion. As computationally predicted, inhibition of src kinase is shown to potently attenuate podosome formation in lung fibroblasts and extracellular matrix invasion, and bleomycin-induced pulmonary fibrosis, suggesting pharmacological targeting of podosomes as a very promising therapeutic option in pulmonary fibrosis.

Abstract Over the last century, outbreaks and pandemics have occurred with disturbing regularity, necessitating advance preparation and large-scale, coordinated response. Here, we developed a machine learning predictive model of disease severity and length of hospitalization for COVID-19, which can be utilized as a platform for future unknown viral outbreaks. We combined untargeted metabolomics on plasma data obtained from COVID-19 patients (n = 111) during hospitalization and healthy controls (n = 342), clinical and comorbidity data (n = 508) to build this patient triage platform, which consists of three parts: (i) the clinical decision tree, which amongst other biomarkers showed that patients with increased eosinophils have worse disease prognosis and can serve as a new potential biomarker with high accuracy (AUC = 0.974), (ii) the estimation of patient hospitalization length with ± 5 days error (R2 = 0.9765) and (iii) the prediction of the disease severity and the need of patient transfer to the intensive care unit. We report a significant decrease in serotonin levels in patients who needed positive airway pressure oxygen and/or were intubated. Furthermore, 5-hydroxy tryptophan, allantoin, and glucuronic acid metabolites were increased in COVID-19 patients and collectively they can serve as biomarkers to predict disease progression. The ability to quickly identify which patients will develop life-threatening illness would allow the efficient allocation of medical resources and implementation of the most effective medical interventions. We would advocate that the same approach could be utilized in future viral outbreaks to help hospitals triage patients more effectively and improve patient outcomes while optimizing healthcare resources.

Abstract Background Single-cell RNA sequencing (scRNA-seq) technology has enabled assessment of transcriptome-wide changes at single-cell resolution. Due to the heterogeneity in environmental exposure and genetic background across subjects, subject effect contributes to the major source of variation in scRNA-seq data with multiple subjects, which severely confounds cell type specific differential expression (DE) analysis. Moreover, dropout events are prevalent in scRNA-seq data, leading to excessive number of zeroes in the data, which further aggravates the challenge in DE analysis. Results We developed iDESC to detect cell type specific DE genes between two groups of subjects in scRNA-seq data. iDESC uses a zero-inflated negative binomial mixed model to consider both subject effect and dropouts. The prevalence of dropout events (dropout rate) was demonstrated to be dependent on gene expression level, which is modeled by pooling information across genes. Subject effect is modeled as a random effect in the log-mean of the negative binomial component. We evaluated and compared the performance of iDESC with eleven existing DE analysis methods. Using simulated data, we demonstrated that iDESC had well-controlled type I error and higher power compared to the existing methods. Applications of those methods with well-controlled type I error to three real scRNA-seq datasets from the same tissue and disease showed that the results of iDESC achieved the best consistency between datasets and the best disease relevance. Conclusions iDESC was able to achieve more accurate and robust DE analysis results by separating subject effect from disease effect with consideration of dropouts to identify DE genes, suggesting the importance of considering subject effect and dropouts in the DE analysis of scRNA-seq data with multiple subjects.

Abstract Chronic obstructive pulmonary disease (COPD) and interstitial lung disease (ILD) are clinically and molecularly heterogeneous diseases. We utilized clustering and integrative network analyses to elucidate roles for microRNAs (miRNAs) and miRNA isoforms (isomiRs) in COPD and ILD pathogenesis. Short RNA sequencing was performed on 351 lung tissue samples of COPD (n = 145), ILD (n = 144) and controls (n = 64). Five distinct subclusters of samples were identified including 1 COPD-predominant cluster and 2 ILD-predominant clusters which associated with different clinical measurements of disease severity. Utilizing 262 samples with gene expression and SNP microarrays, we built disease-specific genetic and expression networks to predict key miRNA regulators of gene expression. Members of miR-449/34 family, known to promote airway differentiation by repressing the Notch pathway, were among the top connected miRNAs in both COPD and ILD networks. Genes associated with miR-449/34 members in the disease networks were enriched among genes that increase in expression with airway differentiation at an air–liquid interface. A highly expressed isomiR containing a novel seed sequence was identified at the miR-34c-5p locus. 47% of the anticorrelated predicted targets for this isomiR were distinct from the canonical seed sequence for miR-34c-5p. Overexpression of the canonical miR-34c-5p and the miR-34c-5p isomiR with an alternative seed sequence down-regulated NOTCH1 and NOTCH4. However, only overexpression of the isomiR down-regulated genes involved in Ras signaling such as CRKL and GRB2. Overall, these findings elucidate molecular heterogeneity inherent across COPD and ILD patients and further suggest roles for miR-34c in regulating disease-associated gene-expression.

Abstract Three and a half years after the pandemic outbreak, now that WHO has formally declared that the emergency is over, COVID-19 is still a significant global issue. Here, we focus on recent developments in genetic and genomic research on COVID-19, and we give an outlook on state-of-the-art therapeutical approaches, as the pandemic is gradually transitioning to an endemic situation. The sequencing and characterization of rare alleles in different populations has made it possible to identify numerous genes that affect either susceptibility to COVID-19 or the severity of the disease. These findings provide a beginning to new avenues and pan-ethnic therapeutic approaches, as well as to potential genetic screening protocols. The causative virus, SARS-CoV-2, is still in the spotlight, but novel threatening virus could appear anywhere at any time. Therefore, continued vigilance and further research is warranted. We also note emphatically that to prevent future pandemics and other world-wide health crises, it is imperative to capitalize on what we have learnt from COVID-19: specifically, regarding its origins, the world’s response, and insufficient preparedness. This requires unprecedented international collaboration and timely data sharing for the coordination of effective response and the rapid implementation of containment measures.

Introduction Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial pneumonia marked by progressive lung fibrosis and a poor prognosis. Recent studies have highlighted the potential role of infection in the pathogenesis of IPF, and a prior association of the HLA-DQB1 gene with idiopathic fibrotic interstitial pneumonia (including IPF) has been reported. Owing to the important role that the human leukocyte antigen (HLA) region plays in the immune response, here we evaluated if HLA genetic variation was associated specifically with IPF risk. Methods We performed a meta-analysis of associations of the HLA region with IPF risk in individuals of European ancestry from seven independent case–control studies of IPF (comprising 5159 cases and 27 459 controls, including a prior study of fibrotic interstitial pneumonia). Single nucleotide polymorphisms, classical HLA alleles and amino acids were analysed and signals meeting a region-wide association threshold of p90% were considered significant. We sought to replicate the previously reported HLA-DQB1 association in the subset of studies independent of the original report. Results The meta-analysis of all seven studies identified four significant independent single nucleotide polymorphisms associated with IPF risk. However, none met the posterior probability for replication criterion. The HLA-DQB1 association was not replicated in the independent IPF studies. Conclusion Variation in the HLA region was not consistently associated with risk in studies of IPF. However, this does not preclude the possibility that other genomic regions linked to the immune response may be involved in the aetiology of IPF.

Abstract Engineered whole lungs may one day expand therapeutic options for patients with end-stage lung disease. However, the feasibility of ex vivo lung regeneration remains limited by the inability to recapitulate mature, functional alveolar epithelium. Here, we modulate multimodal components of the alveolar epithelial type 2 cell (AEC2) niche in decellularized lung scaffolds in order to guide AEC2 behavior for epithelial regeneration. First, endothelial cells coordinate with fibroblasts, in the presence of soluble growth and maturation factors, to promote alveolar scaffold population with surfactant-secreting AEC2s. Subsequent withdrawal of Wnt and FGF agonism synergizes with tidal-magnitude mechanical strain to induce the differentiation of AEC2s to squamous type 1 AECs (AEC1s) in cultured alveoli, in situ. These results outline a rational strategy to engineer an epithelium of AEC2s and AEC1s contained within epithelial-mesenchymal-endothelial alveolar-like units, and highlight the critical interplay amongst cellular, biochemical, and mechanical niche cues within the reconstituting alveolus.

RationaleCOPD is characterized by chronic airway inflammation, small airways changes, with disappearance and obstruction, and also distal/alveolar destruction (emphysema). The chronology by which these three features evolve with altered mucosal immunity remains elusive. This study assessed the mucosal immune defense in human control and end-stage COPD lungs, by detailed microCT and RNA transcriptomic analysis of diversely affected zones.MethodsIn 11 control (non-used donors) and 11 COPD (end-stage) explant frozen lungs, 4 cylinders/cores were processed per lung for microCT and tissue transcriptomics. MicroCT was used to quantify tissue percentage and alveolar surface density to classify the COPD cores in mild, moderate and severe alveolar destruction groups, as well as to quantify terminal bronchioles in each group. Transcriptomics of each core assessed fold changes in innate and adaptive cells and pathway enrichment score between control and COPD cores. Immunostainings of immune cells were performed for validation.ResultsIn mildly affected zones, decreased defensins and increased mucus production were observed, along CD8+ T cell accumulation and activation of the IgA pathway. In more severely affected zones, CD68+ myeloid antigen-presenting cells, CD4+ T cells and B cells, as well as MHCII and IgA pathway genes were upregulated. In contrast, terminal bronchioles were decreased in all COPD cores.ConclusionSpatial investigation of end-stage COPD lungs show that mucosal defense dysregulation with decreased defensins and increased mucus and IgA responses, start concomitantly with CD8+ T-cell accumulation in mild emphysema zones, where terminal bronchioles are already decreased. In contrast, adaptive Th and B cell activation is observed in areas with more advanced tissue destruction. This study suggests that in COPD innate immune alterations occur early in the tissue destruction process, which affects both the alveoli and the terminal bronchioles, before the onset of an adaptive immune response.

The functional role of airway basal cells has not been comprehensively studied in idiopathic pulmonary fibrosis (IPF). Here, the authors show that airway basal cells of IPF patients display a distinct phenotype, are profibrotic if transplanted to mice and that fibrosis can be ameliorated by Src iinhibitors.

Finding disease-relevant tissues and cell types can facilitate the identification and investigation of functional genes and variants. In particular, cell type proportions can serve as potential disease predictive biomarkers. In this manuscript, we introduce a novel statistical framework, cell-type Wide Association Study (cWAS), that integrates genetic data with transcriptomics data to identify cell types whose genetically regulated proportions (GRPs) are disease/trait-associated. On simulated and real GWAS data, cWAS showed good statistical power with newly identified significant GRP associations in disease-associated tissues. More specifically, GRPs of endothelial and myofibroblasts in lung tissue were associated with Idiopathic Pulmonary Fibrosis and Chronic Obstructive Pulmonary Disease, respectively. For breast cancer, the GRP of blood CD8+ T cells was negatively associated with breast cancer (BC) risk as well as survival. Overall, cWAS is a powerful tool to reveal cell types associated with complex diseases mediated by GRPs.

Abstract The increased resolution of single-cell RNA-sequencing technologies has led to major breakthroughs and improved our understanding of the normal and pathologic conditions of multiple tissues and organs. In the study of parenchymal lung disease, single-cell RNA-sequencing has better delineated known cell populations and identified novel cells and changes in cellular phenotypes and gene expression patterns associated with disease. In this review, we aim to highlight the advances and insights that have been made possible by applying these technologies to two seemingly very different lung diseases: fibrotic interstitial lung diseases, a group of relentlessly progressive lung diseases leading to pulmonary fibrosis, and COVID-19 pneumonia, an acute viral disease with life-threatening complications, including pulmonary fibrosis. We discuss changes in cell populations and gene expression, highlighting potential common features, such as alveolar cell epithelial injury and aberrant repair and monocyte-derived macrophage populations, as well as relevance and implications to mechanisms of disease and future directions.

Summary: Most cell types in multicellular organisms can perform multiple functions. However, not all functions can be optimally performed simultaneously by the same cells. Functions incompatible at the level of individual cells can be performed at the cell population level, where cells divide labor and specialize in different functions. Division of labor can arise due to instruction by tissue environment or through self-organization. Here, we develop a computational framework to investigate the contribution of these mechanisms to division of labor within a cell-type population. By optimizing collective cellular task performance under trade-offs, we find that distinguishable expression patterns can emerge from cell-cell interactions versus instructive signals. We propose a method to construct ligand-receptor networks between specialist cells and use it to infer division-of-labor mechanisms from single-cell RNA sequencing (RNA-seq) and spatial transcriptomics data of stromal, epithelial, and immune cells. Our framework can be used to characterize the complexity of cell interactions within tissues.

Mitochondrial dysfunction at birth predicts bronchopulmonary dysplasia (BPD) in extremely low–birth weight (ELBW) infants. Recently, nebulized thyroid hormone (TH), given as triiodothyronine (T3), was noted to decrease pulmonary fibrosis in adult animals through improved mitochondrial function. In this study, we tested the hypothesis that TH may have similar effects on hyperoxia-induced neonatal lung injury and mitochondrial dysfunction by testing whether i.n. T3 decreases neonatal hyperoxic lung injury in newborn mice; whether T3 improves mitochondrial function in lung homogenates, neonatal murine lung fibroblasts (NMLFs), and umbilical cord–derived mesenchymal stem cells (UC-MSCs) obtained from ELBW infants; and whether neonatal hypothyroxinemia is associated with BPD in ELBW infants. We found that inhaled T3 (given i.n.) attenuated hyperoxia-induced lung injury and mitochondrial dysfunction in newborn mice. T3 also reduced bioenergetic deficits in UC-MSCs obtained from both infants with no or mild BPD and those with moderate to severe BPD. T3 also increased the content of peroxisome proliferator–activated receptor γ coactivator 1α in lung homogenates of mice exposed to hyperoxia as well as mitochondrial potential in both NMLFs and UC-MSCs. ELBW infants who died or developed moderate to severe BPD had lower total T4 (TT4) compared with survivors with no or mild BPD. In conclusion, TH signaling and function may play a critical role in neonatal lung injury, and inhaled T3 supplementation may be useful as a therapeutic strategy for BPD.

Abstract Single-cell RNA-sequencing data has revolutionized our ability to understand of the patterns of cell–cell and ligand–receptor connectivity that influence the function of tissues and organs. However, the quantification and visualization of these patterns in a way that informs tissue biology are major computational and epistemological challenges. Here, we present Connectome, a software package for R which facilitates rapid calculation and interactive exploration of cell–cell signaling network topologies contained in single-cell RNA-sequencing data. Connectome can be used with any reference set of known ligand–receptor mechanisms. It has built-in functionality to facilitate differential and comparative connectomics, in which signaling networks are compared between tissue systems. Connectome focuses on computational and graphical tools designed to analyze and explore cell–cell connectivity patterns across disparate single-cell datasets and reveal biologic insight. We present approaches to quantify focused network topologies and discuss some of the biologic theory leading to their design.

SARS-CoV-2 infection can lead to progressive pathology in patients with COVID-19, but information for this disease progression is sparse. Here the authors use multi-omics approach to profile the immune responses of patients, assessing immune repertoire and effects of tocilizumab treatments, to find a dyssynchrony between innate and adaptive immunity in progressive COVID-19.

Chronic obstructive pulmonary disease is a leading cause of death worldwide, while our understanding of cell-specific mechanisms underlying its pathobiology remains incomplete. Here the authors perform scRNA-seq of human lung tissue to identify transcriptional changes in alveolar niche cells associated with the disease.

Idiopathic pulmonary fibrosis (IPF) is a progressive and ultimately fatal disease. Recent findings have shown a marked metabolic reprogramming associated with changes in mitochondrial homeostasis and autophagy during pulmonary fibrosis. The microRNA-33 (miR-33) family of microRNAs (miRNAs) encoded within the introns of sterol regulatory element binding protein (SREBP) genes are master regulators of sterol and fatty acid (FA) metabolism. miR-33 controls macrophage immunometabolic response and enhances mitochondrial biogenesis, FA oxidation, and cholesterol efflux. Here, we show that miR-33 levels are increased in bronchoalveolar lavage (BAL) cells isolated from patients with IPF compared with healthy controls. We demonstrate that specific genetic ablation of miR-33 in macrophages protects against bleomycin-induced pulmonary fibrosis. The absence of miR-33 in macrophages improves mitochondrial homeostasis and increases autophagy while decreasing inflammatory response after bleomycin injury. Notably, pharmacological inhibition of miR-33 in macrophages via administration of anti–miR-33 peptide nucleic acids (PNA-33) attenuates fibrosis in different in vivo and ex vivo mice and human models of pulmonary fibrosis. These studies elucidate a major role of miR-33 in macrophages in the regulation of pulmonary fibrosis and uncover a potentially novel therapeutic approach to treat this disease.

The pluripotency factor KLF4 has been described as pro-fibrotic or anti-fibrotic in various diseases. Herein, the authors show that during lung fibrosis, these distinct effects can be attributed to mesenchymal cell-type specific functions of KLF4.

Abstract Background Recent development of single cell sequencing technologies has made it possible to identify genes with different expression (DE) levels at the cell type level between different groups of samples. In this article, we propose to borrow information through known biological networks to increase statistical power to identify differentially expressed genes (DEGs). Results We develop MRFscRNAseq, which is based on a Markov random field (MRF) model to appropriately accommodate gene network information as well as dependencies among cell types to identify cell-type specific DEGs. We implement an Expectation-Maximization (EM) algorithm with mean field-like approximation to estimate model parameters and a Gibbs sampler to infer DE status. Simulation study shows that our method has better power to detect cell-type specific DEGs than conventional methods while appropriately controlling type I error rate. The usefulness of our method is demonstrated through its application to study the pathogenesis and biological processes of idiopathic pulmonary fibrosis (IPF) using a single-cell RNA-sequencing (scRNA-seq) data set, which contains 18,150 protein-coding genes across 38 cell types on lung tissues from 32 IPF patients and 28 normal controls. Conclusions The proposed MRF model is implemented in the R package MRFscRNAseq available on GitHub. By utilizing gene-gene and cell-cell networks, our method increases statistical power to detect differentially expressed genes from scRNA-seq data.

Summary: Background: MicroRNAs are non-coding RNAs that negatively regulate gene networks. Previously, we reported that systemically delivered miR-29 mimic MRG-201 reduced fibrosis in animal models, supporting the consideration of miR-29-based therapies for idiopathic pulmonary fibrosis (IPF). Methods: We generated MRG-229, a next-generation miR-29 mimic based on MRG-201 with improved chemical stability due to additional sugar modifications and conjugation with the internalization moiety BiPPB (PDGFbetaR-specific bicyclic peptide)1. We investigated the anti-fibrotic efficacy of MRG-229 on TGF-β1 treated human lung fibroblasts (NHLFs), human precision cut lung slices (hPCLS), and in vivo bleomycin studies; toxicology was assessed in two animal models, rats, and non-human primates. Finally, we examined miR-29b levels in a cohort of 46 and 213 patients with IPF diagnosis recruited from Yale and Nottingham Universities (Profile Cohort), respectively. Findings: The peptide-conjugated MRG-229 mimic decreased expression of pro-fibrotic genes and reduced collagen production in each model. In bleomycin-treated mice, the peptide-conjugated MRG-229 mimic downregulated profibrotic gene programs at doses more than ten-fold lower than the original compound. In rats and non-human primates, the peptide-conjugated MRG-229 mimic was well tolerated at clinically relevant doses with no adverse findings observed. In human peripheral blood from IPF patients decreased miR-29 concentrations were associated with increased mortality in two cohorts potentially identified as a target population for treatment. Interpretation: Collectively, our results provide support for the development of the peptide-conjugated MRG-229 mimic as a potential therapy in humans with IPF. Funding: This work was supported by NIH NHLBI grants UH3HL123886, R01HL127349, R01HL141852, U01HL145567.

Studies comparing single cell RNA-Seq (scRNA-Seq) data between conditions mainly focus on differences in the proportion of cell types or on differentially expressed genes. In many cases these differences are driven by changes in cell interactions which are challenging to infer without spatial information. To determine cell-cell interactions that differ between conditions we developed the Cell Interaction Network Inference (CINS) pipeline. CINS combines Bayesian network analysis with regression-based modeling to identify differential cell type interactions and the proteins that underlie them. We tested CINS on a disease case control and on an aging mouse dataset. In both cases CINS correctly identifies cell type interactions and the ligands involved in these interactions improving on prior methods suggested for cell interaction predictions. We performed additional mouse aging scRNA-Seq experiments which further support the interactions identified by CINS.

Patients with chronic lung disease (CLD) have an increased risk for severe coronavirus disease-19 and poor outcomes. Here the authors compare the transcriptomes of single cells isolated from healthy and CLD lungs to identify molecular characteristics of lung cells that may account for worse COVID-19 outcomes in these patients.

Background: We have previously shown that SHP2 downregulation may predispose fibroblasts to differentiate into myofibroblasts and proposed a role for SHP2 downregulation in the pathogenesis of idiopathic pulmonary fibrosis (IPF). Recent data have shown that SHP2 localizes to the mitochondrial intercristae, and its overexpression enhances mitochondrial metabolism leading to oxidative stress and senescence. Objective: To determine the effect of SHP2 on fibrotic responses. Methods and Results: Primary mouse lung fibroblasts derived from mice carrying a conditional knock-in mutation (D61G/+), rendering the SHP2 catalytic domain constitutively active, had reduced proliferation (1.6-fold, p < 0.05), migration (2-fold, p < 0.05), as well as reduced responsiveness of TGFB-1 induced fibroblasts-to-myofibroblasts differentiation, compared to wild-type ones. Electron microscope analysis revealed that SHP2 D61G/+ mouse lung fibroblasts were characterized by mitochondrial abnormalities, including swollen mitochondria with disrupted electron-lucent cristae and an increased number of autophagosomes compared to wild-type ones. SHP2 D61G/+ MLFs exhibited increased protein levels of autophagy markers, including LC3B-II and p-62, evidence that was confirmed by immunofluorescence analysis. Mitochondrial function analysis revealed that stable (genotype D61G/+) overexpression of SHP2 led to impaired mitochondrial function, as assessed by decreased mitochondrial membrane potential (1.29-fold, p < 0.05), coupling efficiency (1.82 fold, p < 0.05), oxygen consumption rate (1.9-fold, p < 0.05), and increased reactive oxygen species production both at baseline (1.75-fold, p < 0.05) and following H2O2 stimulation (1.63-fold, p < 0.05) compared to wild-type ones (SHP2+/+). SHP2 D61G/+ mouse lung fibroblasts showed enhanced AMPK activity, as well as decreased activation of the mTORC1 signaling pathway, potentially leading to ineffective mitochondrial metabolism and increased autophagy. Conclusions: SHP2 attenuates fibrotic responses in fibroblast cell lines through negative regulation of mitochondrial metabolism and induction of autophagy. SHP2 activation may represent a promising therapeutic strategy for patients with fibrotic lung diseases.

Abstract The recent coronavirus disease (COVID-19), caused by SARS-CoV-2, is inarguably the most challenging coronavirus outbreak relative to the previous outbreaks involving SARS-CoV and MERS-CoV. With the number of COVID-19 cases now exceeding 2 million worldwide, it is apparent that (i) transmission of SARS-CoV-2 is very high and (ii) there are large variations in disease severity, one component of which may be genetic variability in the response to the virus. Controlling current rates of infection and combating future waves require a better understanding of the routes of exposure to SARS-CoV-2 and the underlying genomic susceptibility to this disease. In this mini-review, we highlight possible genetic determinants of COVID-19 and the contribution of aerosol exposure as a potentially important transmission route of SARS-CoV-2.

Collagen production by lung cells is critical to maintain organ architecture but can also drive pathological scarring. Here the authors perform single cell RNA sequencing of collagen-producing lung cells identifying a subset of pathologic fibroblasts characterized by Cthrc1 expression which are concentrated within fibroblastic foci in fibrotic lungs and show a pro-fibrotic phenotype.

Abstract Background Numerous publications suggest an association between herpes virus infection and idiopathic pulmonary fibrosis (IPF). These reports have employed immunohistochemistry, in situ hybridization and/or PCR, which are susceptible to specificity artifacts. Methods We investigated the possible association between IPF and viral RNA expression using next-generation sequencing, which has the potential to provide a high degree of both sensitivity and specificity. We quantified viral RNA expression for 740 viruses in 28 IPF patient lung biopsy samples and 20 controls. Key RNA-seq results were confirmed using Real-time RT-PCR for select viruses (EBV, HCV, herpesvirus saimiri and HERV-K). Results We identified sporadic low-level evidence of viral infections in our lung tissue specimens, but did not find a statistical difference for expression of any virus, including EBV, herpesvirus saimiri and HERV-K, between IPF and control lungs. Conclusions To the best of our knowledge, this is the first publication that employs RNA-seq to assess whether viral infections are linked to the pathogenesis of IPF. Our results do not address the role of viral infection in acute exacerbations of IPF, however, this analysis patently did not support an association between herpes virus detection and IPF.

Extracellular matrix (ECM) stiffening with downstream activation of mechanosensitive pathways is strongly implicated in fibrosis. We previously reported that altered collagen nanoarchitecture is a key determinant of pathogenetic ECM structure-function in human fibrosis (Jones et al., 2018). Here, through human tissue, bioinformatic and ex vivo studies we provide evidence that hypoxia-inducible factor (HIF) pathway activation is a critical pathway for this process regardless of the oxygen status (pseudohypoxia). Whilst TGFβ increased the rate of fibrillar collagen synthesis, HIF pathway activation was required to dysregulate post-translational modification of fibrillar collagen, promoting pyridinoline cross-linking, altering collagen nanostructure, and increasing tissue stiffness. In vitro, knockdown of Factor Inhibiting HIF (FIH), which modulates HIF activity, or oxidative stress caused pseudohypoxic HIF activation in the normal fibroblasts. By contrast, endogenous FIH activity was reduced in fibroblasts from patients with lung fibrosis in association with significantly increased normoxic HIF pathway activation. In human lung fibrosis tissue, HIF-mediated signalling was increased at sites of active fibrogenesis whilst subpopulations of human lung fibrosis mesenchymal cells had increases in both HIF and oxidative stress scores. Our data demonstrate that oxidative stress can drive pseudohypoxic HIF pathway activation which is a critical regulator of pathogenetic collagen structure-function in fibrosis.

Asthma is a common disease with profoundly variable natural history and patient morbidity. Heterogeneity has long been appreciated, and much work has focused on identifying subgroups of patients with similar pathobiological underpinnings. Previous studies of the Severe Asthma Research Program (SARP) cohort linked gene expression changes to specific clinical and physiologic characteristics. While invaluable for hypothesis generation, these data include extensive candidate gene lists that complicate target identification and validation. In this analysis, we performed unsupervised clustering of the SARP cohort using bronchial epithelial cell gene expression data, identifying a transcriptional signature for participants suffering exacerbation-prone asthma with impaired lung function. Clinically, participants in this asthma cluster exhibited a mixed inflammatory process and bore transcriptional hallmarks of NF-κB and activator protein 1 (AP-1) activation, despite high corticosteroid exposure. Using supervised machine learning, we found a set of 31 genes that classified patients with high accuracy and could reconstitute clinical and transcriptional hallmarks of our patient clustering in an external cohort. Of these genes, IL18R1 (IL-18 Receptor 1) negatively associated with lung function and was highly expressed in the most severe patient cluster. We validated IL18R1 protein expression in lung tissue and identified downstream NF-κB and AP-1 activity, supporting IL-18 signaling in severe asthma pathogenesis and highlighting this approach for gene and pathway discovery.

The development of an in vitro system for the study of lung vascular disease is critical to understanding human pathologies. Conventional culture systems fail to fully recapitulate native microenvironmental conditions and are typically limited in their ability to represent human pathophysiology for the study of disease and drug mechanisms. Whole organ decellularization provides a means to developing a construct that recapitulates structural, mechanical, and biological features of a complete vascular structure. Here, we developed a culture protocol to improve endothelial cell coverage in whole lung scaffolds and used single-cell RNA-sequencing analysis to explore the impact of decellularized whole lung scaffolds on endothelial phenotypes and functions in a biomimetic bioreactor system. Intriguingly, we found that the phenotype and functional signals of primary pulmonary microvascular revert back—at least partially—toward native lung endothelium. Additionally, human induced pluripotent stem cell-derived endothelium cultured in decellularized lung systems start to gain various native human endothelial phenotypes. Vascular barrier function was partially restored, while small capillaries remained patent in endothelial cell-repopulated lungs. To evaluate the ability of the engineered endothelium to modulate permeability in response to exogenous stimuli, lipopolysaccharide (LPS) was introduced into repopulated lungs to simulate acute lung injury. After LPS treatment, proinflammatory signals were significantly increased and the vascular barrier was impaired. Taken together, these results demonstrate a novel platform that recapitulates some pulmonary microvascular functions and phenotypes at a whole organ level. This development may help pave the way for using the whole organ engineering approach to model vascular diseases.

Hypoxia adversely affects the pulmonary circulation of mammals, including vasoconstriction leading to elevated pulmonary arterial pressures. The clinical importance of changes in the structure and function of the large, elastic pulmonary arteries is gaining increased attention, particularly regarding impact in multiple chronic cardiopulmonary conditions. We establish a multi-disciplinary workflow to understand better transcriptional, microstructural, and functional changes of the pulmonary artery in response to sustained hypoxia and how these changes inter-relate. We exposed adult male C57BL/6J mice to normoxic or hypoxic (FiO2 10%) conditions. Excised pulmonary arteries were profiled transcriptionally using single cell RNA sequencing, imaged with multiphoton microscopy to determine microstructural features under in vivo relevant multiaxial loading, and phenotyped biomechanically to quantify associated changes in material stiffness and vasoactive capacity. Pulmonary arteries of hypoxic mice exhibited an increased material stiffness that was likely due to collagen remodeling rather than excessive deposition (fibrosis), a change in smooth muscle cell phenotype reflected by decreased contractility and altered orientation aligning these cells in the same direction as the remodeled collagen fibers, endothelial proliferation likely representing endothelial-to-mesenchymal transitioning, and a network of cell-type specific transcriptomic changes that drove these changes. These many changes resulted in a system-level increase in pulmonary arterial pulse wave velocity, which may drive a positive feedback loop exacerbating all changes. These findings demonstrate the power of a multi-scale genetic-functional assay. They also highlight the need for systems-level analyses to determine which of the many changes are clinically significant and may be potential therapeutic targets.

This data is a curated collection of visual images of gene expression patterns from the pre- and post-natal mouse lung, accompanied by associated mRNA probe sequences and RNA-Seq expression profiles. Mammalian lungs undergo significant growth and cellular differentiation before and after the transition to breathing air. Documenting normal lung development is an important step in understanding abnormal lung development, as well as the challenges faced during a preterm birth. Images in this dataset indicate the spatial distribution of mRNA transcripts for over 500 different genes that are active during lung development, as initially determined via RNA-Seq. Images were systematically acquired using high-throughput in situ hybridization with non-radioactive digoxigenin-labeled mRNA probes across mouse lungs from developmental time points E16.5, E18.5, P7, and P28. The dataset was produced as part of The Molecular Atlas of Lung Development Program (LungMAP) and is hosted at https://lungmap.net. This manuscript describes the nature of the data and the protocols for generating the dataset.

Abstract Long‐noncoding RNAs (lncRNAs) have numerous biological functions controlling cell differentiation and tissue development. The knowledge about the role of lncRNAs in human lungs remains limited. Here we found the regulatory role of the terminal differentiation‐induced lncRNA (TINCR) in bronchial cell differentiation. RNA in situ hybridization revealed that TINCR was mainly expressed in bronchial epithelial cells in normal human lung. We performed RNA sequencing analysis of normal human bronchial epithelial cells (NHBECs) with or without TINCR inhibition and found the differential expression of 603 genes, which were enriched for cell adhesion and migration, wound healing, extracellular matrix organization, tissue development and differentiation. To investigate the role of TINCR in the differentiation of NHBECs, we employed air–liquid interface culture and 3D organoid formation assay. TINCR was upregulated during differentiation, loss of TINCR significantly induced an early basal‐like cell phenotype (TP63) and a ciliated cell differentiation (FOXJ1) in late phase and TINCR overexpression suppressed basal cell phenotype and the differentiation toward to ciliated cells. Critical regulators of differentiation such as SOX2 and NOTCH genes (NOTCH1, HES1, and JAG1) were significantly upregulated by TINCR inhibition and downregulated by TINCR overexpression. RNA immunoprecipitation assay revealed that TINCR was required for the direct bindings of Staufen1 protein to SOX2, HES1, and JAG1 mRNA. Loss of Staufen1 induced TP63, SOX2, NOTCH1, HES1, and JAG1 mRNA expressions, which TINCR overexpression suppressed partially. In conclusion, TINCR is a novel regular of bronchial cell differentiation, affecting downstream regulators such as SOX2 and NOTCH genes, potentially in coordination with Staufen1.

The pathogenesis of chronic obstructive pulmonary disease (COPD) involves aberrant responses to cellular stress caused by chronic cigarette smoke (CS) exposure. However, not all smokers develop COPD and the critical mechanisms that regulate cellular stress responses to increase COPD susceptibility are not understood. Because microRNAs are well-known regulators of cellular stress responses, we evaluated microRNA expression arrays performed on distal parenchymal lung tissue samples from 172 subjects with and without COPD. We identified miR-24-3p as the microRNA that best correlated with radiographic emphysema and validated this finding in multiple cohorts. In a CS exposure mouse model, inhibition of miR-24-3p increased susceptibility to apoptosis, including alveolar type II epithelial cell apoptosis, and emphysema severity. In lung epithelial cells, miR-24-3p suppressed apoptosis through the BH3-only protein BIM and suppressed homology-directed DNA repair and the DNA repair protein BRCA1. Finally, we found BIM and BRCA1 were increased in COPD lung tissue, and BIM and BRCA1 expression inversely correlated with miR-24-3p. We concluded that miR-24-3p, a regulator of the cellular response to DNA damage, is decreased in COPD, and decreased miR-24-3p increases susceptibility to emphysema through increased BIM and apoptosis.

Antibodies are central effectors of the adaptive immune response, widespread used therapeutics, but also potentially disease-causing biomolecules. Antibody folding catalysts in the plasma cell are incompletely defined. Idiopathic pulmonary fibrosis (IPF) is a fatal chronic lung disease with increasingly recognized autoimmune features. We found elevated expression of FK506-binding protein 11 (FKBP11) in IPF lungs where FKBP11 specifically localized to antibody-producing plasma cells. Suggesting a general role in plasma cells, plasma cell-specific FKBP11 expression was equally observed in lymphatic tissues, and in vitro B cell to plasma cell differentiation was accompanied by induction of FKBP11 expression. Recombinant human FKBP11 was able to refold IgG antibody in vitro and inhibited by FK506, strongly supporting a function as antibody peptidyl-prolyl cis-trans isomerase. Induction of ER stress in cell lines demonstrated induction of FKBP11 in the context of the unfolded protein response in an X-box-binding protein 1 (XBP1)-dependent manner. While deficiency of FKBP11 increased susceptibility to ER stress-mediated cell death in an alveolar epithelial cell line, FKBP11 knockdown in an antibody-producing hybridoma cell line neither induced cell death nor decreased expression or secretion of IgG antibody. Similarly, antibody secretion by the same hybridoma cell line was not affected by knockdown of the established antibody peptidyl-prolyl isomerase cyclophilin B. The results are consistent with FKBP11 as a novel XBP1-regulated antibody peptidyl-prolyl cis-trans isomerase and indicate significant redundancy in the ER-resident folding machinery of antibody-producing hybridoma cells.

Summary: Cell-based therapies have shown promise for treating myriad chronic pulmonary diseases through direct application of epithelial progenitors or by way of engineered tissue grafts or whole organs. To elucidate environmental effects on epithelial regenerative outcomes in vitro, here, we isolate and culture a population of pharmacologically expanded basal cells (peBCs) from rat tracheas. At peak basal marker expression, we simultaneously split peBCs into four in vitro platforms: organoid, air-liquid interface (ALI), engineered trachea, and engineered lung. Following differentiation, these samples are evaluated using single-cell RNA sequencing (scRNA-seq) and computational pipelines are developed to compare samples both globally and at the population level. A sample of native rat tracheal epithelium is also evaluated by scRNA-seq as a control for engineered epithelium. Overall, this work identifies platform-specific effects that support the use of engineered models to achieve the most physiologic differential outcomes in pulmonary epithelial regenerative applications. : Greaney et al. compare pulmonary epithelial regeneration across multiple modalities in vitro, finding that decellularized scaffolds achieved the most physiologic differentiation over more artificial platforms. scRNA-seq enables high-resolution comparison between engineered and native cell populations, thereby better gauging progress toward the generation of a tissue that may function on implantation. Keywords: Pulmonary Cell Biology, Regenerative Medicine, Tissue Engineering, Primary Epithelial Cells

Abstract Background Lung Cancer is occasionally observed in patients with Idiopathic Pulmonary Fibrosis (IPF). We sought to describe the epidemiologic and clinical characteristics of lung cancer for patients with IPF and other interstitial lung disease (ILD) using institutional and statewide data registries. Methods We conducted a retrospective analysis of IPF and non-IPF ILD patients from the ILD center registry, to compare with lung cancer registries at the University of Pittsburgh as well as with population data of lung cancer obtained from Pennsylvania Department of Health between 2000 and 2015. Results Among 1108 IPF patients, 31 patients were identified with IPF and lung cancer. The age-adjusted standard incidence ratio of lung cancer was 3.34 (with IPF) and 2.3 (with non-IPF ILD) (between-group Hazard ratio = 1.4, p = 0.3). Lung cancer worsened the mortality of IPF (p

Abstract Background Idiopathic pulmonary fibrosis is a fatal lung disease characterized by a progressive formation of fibroblastic foci in the interstitium. This disease is strongly associated with telomere dysfunction but the extent of telomere shortening and consequent chromosomal damage within IPF lungs and with regional disease severity remains unknown. Methods Explanted IPF lungs (n = 10) were collected from transplant surgeries with six samples per lung analysed to capture the regional heterogeneity ranging from mild to severe disease. Non-used donor lungs (n = 6) were collected as “healthy” controls. Structural changes related to disease severity (microCT surface density), relative telomere length (real-time qPCR), and quantitative histology of chromosomal damage (γ-H2A.X) and extracellular matrix (elastin, total collagen, collagen 1, and collagen 3) were measured. A multivariate linear mixed-effects model controlling for subject was used to identify association of disease severity or fibrotic markers with telomere length and chromosomal damage. Results We observed shorter telomere length (p = 0.001) and increased chromosomal damage (p = 0.018) in IPF lungs compared to controls. In IPF lungs, telomere length was associated with total collagen (p

Abstract Background Telomere shortening has been associated with several lung diseases. However, telomere length is generally measured in peripheral blood leucocytes rather than in lung tissue, where disease occurs. Consequently, telomere dynamics have not been established for the normal human lung nor for diseased lung tissue. We hypothesized an age- and disease-dependent shortening of lung tissue telomeres. Methods At time of (re-)transplantation or autopsy, 70 explant lungs were collected: from unused donors (normal, n = 13) and patients with cystic fibrosis (CF, n = 12), chronic obstructive pulmonary disease (COPD, n = 11), chronic hypersensitivity pneumonitis (cHP, n = 9), bronchiolitis obliterans syndrome (BOS) after prior transplantation (n = 11) and restrictive allograft syndrome (RAS) after prior transplantation (n = 14). Lungs were inflated, frozen and then scanned using CT. Four tissue cores from distinct lung regions were sampled for analysis. Disease severity was evaluated using CT and micro CT imaging. DNA was extracted from the samples and average relative telomere length (RTL) was determined using real-time qPCR. Results The normal lungs showed a decrease in RTL with age (p

Vascular smooth muscle cell (VSMC) dysfunction is a common feature of thoracic aortic aneurysms (TAAs). Here, Lino Cardenas and colleagues show that the formation of a HDAC9-MALAT1-BRG1 complex promotes VSMC dysfunction in TAA by epigenetically altering the expression of key components of the cytoskeleton in VSMCs.

Regenerative capacity in vital organs is limited by fibrosis propensity. Idiopathic pulmonary fibrosis (IPF), a progressive lung disease linked with aging, is a classic example. In this study, we show that in flow cytometry, immunoblots (IB) and in lung sections, FLIP levels can be regulated, in vivo and in vitro, through SIRT1 activity inhibition by CMH (4-(4-Chloro-2-methylphenoxy)-N-hydroxybutanamide), a small molecule that, as we determined here by structural biology calculations, docked into its nonhistone substrate Ku70-binding site. Ku70 immunoprecipitations and immunoblots confirmed our theory that Ku70-deacetylation, Ku70/FLIP complex, myofibroblast resistance to apoptosis, cell survival, and lung fibrosis in bleomycin-treated mice, are reduced and regulated by CMH. Thus, small molecules associated with SIRT1-mediated regulation of Ku70 deacetylation, affecting FLIP stabilization in fibrotic-lung myofibroblasts, may be a useful strategy, enabling tissue regeneration.

Idiopathic pulmonary fibrosis (IPF) is a lethal fibrotic lung disease characterized by aberrant remodeling of the lung parenchyma with extensive changes to the phenotypes of all lung resident cells. The introduction of transcriptomics, genome scale profiling of thousands of RNA transcripts, caused a significant inversion in IPF research. Instead of generating hypotheses based on animal models of disease, or biological plausibility, with limited validation in humans, investigators were able to generate hypotheses based on unbiased molecular analysis of human samples and then use animal models of disease to test their hypotheses. In this review, we describe the insights made from transcriptomic analysis of human IPF samples. We describe how transcriptomic studies led to identification of novel genes and pathways involved in the human IPF lung such as: matrix metalloproteinases, WNT pathway, epithelial genes, role of microRNAs among others, as well as conceptual insights such as the involvement of developmental pathways and deep shifts in epithelial and fibroblast phenotypes. The impact of lung and transcriptomic studies on disease classification, endotype discovery, and reproducible biomarkers is also described in detail. Despite these impressive achievements, the impact of transcriptomic studies has been limited because they analyzed bulk tissue and did not address the cellular and spatial heterogeneity of the IPF lung. We discuss new emerging technologies and applications, such as single-cell RNAseq and microenvironment analysis that may address cellular and spatial heterogeneity. We end by making the point that most current tissue collections and resources are not amenable to analysis using the novel technologies. To take advantage of the new opportunities, we need new efforts of sample collections, this time focused on access to all the microenvironments and cells in the IPF lung.

The Dynamic Regulatory Events Miner (DREM) software reconstructs dynamic regulatory networks by integrating static protein-DNA interaction data with time series gene expression data. In recent years, several additional types of high-throughput time series data have been profiled when studying biological processes including time series miRNA expression, proteomics, epigenomics and single cell RNA-Seq. Combining all available time series and static datasets in a unified model remains an important challenge and goal. To address this challenge we have developed a new version of DREM termed interactive DREM (iDREM). iDREM provides support for all data types mentioned above and combines them with existing interaction data to reconstruct networks that can lead to novel hypotheses on the function and timing of regulators. Users can interactively visualize and query the resulting model. We showcase the functionality of the new tool by applying it to microglia developmental data from multiple labs.

Investigation of the Th1 immune response in sarcoidosis CD4+ T cells has revealed reduced proliferative capacity and cytokine expression upon TCR stimulation. In other disease models, such cellular dysfunction has been associated with a step-wise, progressive loss of T cell function that results from chronic antigenic stimulation. T cell exhaustion is defined by decreased cytokine production upon TCR activation, decreased proliferation, increased expression of inhibitory cell surface receptors, and increased susceptibility to apoptosis. We characterized sarcoidosis CD4+ T cell immune function in systemic and local environments among subjects undergoing disease progression compared to those experiencing disease resolution. Spontaneous and TCR-stimulated Th1 cytokine expression and proliferation assays were performed in 53 sarcoidosis subjects and 30 healthy controls. PD-1 expression and apoptosis were assessed by flow cytometry. Compared to healthy controls, sarcoidosis CD4+ T cells demonstrated reductions in Th1 cytokine expression, proliferative capacity (p

Biological systems are increasingly being studied by high throughput profiling of molecular data over time. Determining the set of time points to sample in studies that profile several different types of molecular data is still challenging. Here we present the Time Point Selection (TPS) method that solves this combinatorial problem in a principled and practical way. TPS utilizes expression data from a small set of genes sampled at a high rate. As we show by applying TPS to study mouse lung development, the points selected by TPS can be used to reconstruct an accurate representation for the expression values of the non selected points. Further, even though the selection is only based on gene expression, these points are also appropriate for representing a much larger set of protein, miRNA and DNA methylation changes over time. TPS can thus serve as a key design strategy for high throughput time series experiments. Supporting Website: www.sb.cs.cmu.edu/TPS