Search

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.

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.

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.

As miRNAs are associated with normal cellular processes, deregulation of miRNAs is thought to play a causative role in many complex diseases. Nevertheless, the precise contribution of miRNAs in fibrotic lung diseases, especially the idiopathic form (IPF), remains poorly understood. Given the poor response rate of IPF patients to current therapy, new insights into the pathogenic mechanisms controlling lung fibroblasts activation, the key cell type driving the fibrogenic process, are essential to develop new therapeutic strategies for this devastating disease. To identify miRNAs with potential roles in lung fibrogenesis, we performed a genome-wide assessment of miRNA expression in lungs from two different mouse strains known for their distinct susceptibility to develop lung fibrosis after bleomycin exposure. This led to the identification of miR-199a-5p as the best miRNA candidate associated with bleomycin response. Importantly, miR-199a-5p pulmonary expression was also significantly increased in IPF patients (94 IPF versus 83 controls). In particular, levels of miR-199a-5p were selectively increased in myofibroblasts from injured mouse lungs and fibroblastic foci, a histologic feature associated with IPF. Therefore, miR-199a-5p profibrotic effects were further investigated in cultured lung fibroblasts: miR-199a-5p expression was induced upon TGFβ exposure, and ectopic expression of miR-199a-5p was sufficient to promote the pathogenic activation of pulmonary fibroblasts including proliferation, migration, invasion, and differentiation into myofibroblasts. In addition, we demonstrated that miR-199a-5p is a key effector of TGFβ signaling in lung fibroblasts by regulating CAV1, a critical mediator of pulmonary fibrosis. Remarkably, aberrant expression of miR-199a-5p was also found in unilateral ureteral obstruction mouse model of kidney fibrosis, as well as in both bile duct ligation and CCl4-induced mouse models of liver fibrosis, suggesting that dysregulation of miR-199a-5p represents a general mechanism contributing to the fibrotic process. MiR-199a-5p thus behaves as a major regulator of tissue fibrosis with therapeutic potency to treat fibroproliferative diseases.

Idiopathic pulmonary fibrosis (IPF) is a progressive and life threatening disease with median survival of 2.5-3 years. The IPF lung is characterized by abnormal lung remodeling, epithelial cell hyperplasia, myofibroblast foci formation, and extracellular matrix deposition. Analysis of gene expression microarray data revealed that cartilage oligomeric matrix protein (COMP), a non-collagenous extracellular matrix protein is among the most significantly up-regulated genes (Fold change 13, p-value

BACKGROUND:Increased platelet activation in sickle cell disease (SCD) contributes to a state of hypercoagulability and confers a risk of thromboembolic complications. The role for post-transcriptional regulation of the platelet transcriptome by microRNAs (miRNAs) in SCD has not been previously explored. This is the first study to determine whether platelets from SCD exhibit an altered miRNA expression profile. METHODS AND FINDINGS:We analyzed the expression of miRNAs isolated from platelets from a primary cohort (SCD = 19, controls = 10) and a validation cohort (SCD = 7, controls = 7) by hybridizing to the Agilent miRNA microarrays. A dramatic difference in miRNA expression profiles between patients and controls was noted in both cohorts separately. A total of 40 differentially expressed platelet miRNAs were identified as common in both cohorts (p-value 0.05, fold change>2) with 24 miRNAs downregulated. Interestingly, 14 of the 24 downregulated miRNAs were members of three families - miR-329, miR-376 and miR-154 - which localized to the epigenetically regulated, maternally imprinted chromosome 14q32 region. We validated the downregulated miRNAs, miR-376a and miR-409-3p, and an upregulated miR-1225-3p using qRT-PCR. Over-expression of the miR-1225-3p in the Meg01 cells was followed by mRNA expression profiling to identify mRNA targets. This resulted in significant transcriptional repression of 1605 transcripts. A combinatorial approach using Meg01 mRNA expression profiles following miR-1225-3p overexpression, a computational prediction analysis of miRNA target sequences and a previously published set of differentially expressed platelet transcripts from SCD patients, identified three novel platelet mRNA targets: PBXIP1, PLAGL2 and PHF20L1. CONCLUSIONS:We have identified significant differences in functionally active platelet miRNAs in patients with SCD as compared to controls. These data provide an important inventory of differentially expressed miRNAs in SCD patients and an experimental framework for future studies of miRNAs as regulators of biological pathways in platelets.

MicroRNAs (miRNAs) are post-transcriptional regulators that bind to their target mRNAs through base complementarity. Predicting miRNA targets is a challenging task and various studies showed that existing algorithms suffer from high number of false predictions and low to moderate overlap in their predictions. Until recently, very few algorithms considered the dynamic nature of the interactions, including the effect of less specific interactions, the miRNA expression level, and the effect of combinatorial miRNA binding. Addressing these issues can result in a more accurate miRNA:mRNA modeling with many applications, including efficient miRNA-related SNP evaluation. We present a novel thermodynamic model based on the Fermi-Dirac equation that incorporates miRNA expression in the prediction of target occupancy and we show that it improves the performance of two popular single miRNA target finders. Modeling combinatorial miRNA targeting is a natural extension of this model. Two other algorithms show improved prediction efficiency when combinatorial binding models were considered. ComiR (Combinatorial miRNA targeting), a novel algorithm we developed, incorporates the improved predictions of the four target finders into a single probabilistic score using ensemble learning. Combining target scores of multiple miRNAs using ComiR improves predictions over the naïve method for target combination. ComiR scoring scheme can be used for identification of SNPs affecting miRNA binding. As proof of principle, ComiR identified rs17737058 as disruptive to the miR-488-5p:NCOA1 interaction, which we confirmed in vitro. We also found rs17737058 to be significantly associated with decreased bone mineral density (BMD) in two independent cohorts indicating that the miR-488-5p/NCOA1 regulatory axis is likely critical in maintaining BMD in women. With increasing availability of comprehensive high-throughput datasets from patients ComiR is expected to become an essential tool for miRNA-related studies.

BACKGROUND:Idiopathic Pulmonary Fibrosis (IPF) is characterized by profound changes in the lung phenotype including excessive extracellular matrix deposition, myofibroblast foci, alveolar epithelial cell hyperplasia and extensive remodeling. The role of epigenetic changes in determining the lung phenotype in IPF is unknown. In this study we determine whether IPF lungs exhibit an altered global methylation profile. METHODOLOGY/PRINCIPAL FINDINGS:Immunoprecipitated methylated DNA from 12 IPF lungs, 10 lung adenocarcinomas and 10 normal histology lungs was hybridized to Agilent human CpG Islands Microarrays and data analysis was performed using BRB-Array Tools and DAVID Bioinformatics Resources software packages. Array results were validated using the EpiTYPER MassARRAY platform for 3 CpG islands. 625 CpG islands were differentially methylated between IPF and control lungs with an estimated False Discovery Rate less than 5%. The genes associated with the differentially methylated CpG islands are involved in regulation of apoptosis, morphogenesis and cellular biosynthetic processes. The expression of three genes (STK17B, STK3 and HIST1H2AH) with hypomethylated promoters was increased in IPF lungs. Comparison of IPF methylation patterns to lung cancer or control samples, revealed that IPF lungs display an intermediate methylation profile, partly similar to lung cancer and partly similar to control with 402 differentially methylated CpG islands overlapping between IPF and cancer. Despite their similarity to cancer, IPF lungs did not exhibit hypomethylation of long interspersed nuclear element 1 (LINE-1) retrotransposon while lung cancer samples did, suggesting that the global hypomethylation observed in cancer was not typical of IPF. CONCLUSIONS/SIGNIFICANCE:Our results provide evidence that epigenetic changes in IPF are widespread and potentially important. The partial similarity to cancer may signify similar pathogenetic mechanisms while the differences constitute IPF or cancer specific changes. Elucidating the role of these specific changes will potentially allow better understanding of the pathogenesis of IPF.

Idiopathic pulmonary fibrosis (IPF) is a progressive and medically refractory lung disease with a grim prognosis. Although the etiology of IPF remains perplexing, abnormal adaptive immune responses are evident in many afflicted patients. We hypothesized that perturbations of human leukocyte antigen (HLA) allele frequencies, which are often seen among patients with immunologic diseases, may also be present in IPF patients.HLA alleles were determined in subpopulations of IPF and normal subjects using molecular typing methods. HLA-DRB1*15 was over-represented in a discovery cohort of 79 Caucasian IPF subjects who had lung transplantations at the University of Pittsburgh (36.7%) compared to normal reference populations. These findings were prospectively replicated in a validation cohort of 196 additional IPF subjects from four other U.S. medical centers that included both ambulatory patients and lung transplantation recipients. High-resolution typing was used to further define specific HLA-DRB1*15 alleles. DRB1*1501 prevalence in IPF subjects was similar among the 143 ambulatory patients and 132 transplant recipients (31.5% and 34.8%, respectively, p = 0.55). The aggregate prevalence of DRB1*1501 in IPF patients was significantly greater than among 285 healthy controls (33.1% vs. 20.0%, respectively, OR 2.0; 95%CI 1.3-2.9, p = 0.0004). IPF patients with DRB1*1501 (n = 91) tended to have decreased diffusing capacities for carbon monoxide (DL(CO)) compared to the 184 disease subjects who lacked this allele (37.8±1.7% vs. 42.8±1.4%, p = 0.036).DRB1*1501 is more prevalent among IPF patients than normal subjects, and may be associated with greater impairment of gas exchange. These data are novel evidence that immunogenetic processes can play a role in the susceptibility to and/or manifestations of IPF. Findings here of a disease association at the HLA-DR locus have broad pathogenic implications, illustrate a specific chromosomal area for incremental, targeted genomic study, and may identify a distinct clinical phenotype among patients with this enigmatic, morbid lung disease.

BackgroundTransforming growth factor beta 1 (TGFβ1) plays a major role in many lung diseases including lung cancer, pulmonary hypertension, and pulmonary fibrosis. TGFβ1 activates a signal transduction cascade that results in the transcriptional regulation of genes in the nucleus, primarily through the DNA-binding transcription factor SMAD3. The objective of this study is to identify genome-wide scale map of SMAD3 binding targets and the molecular pathways and networks affected by the TGFβ1/SMAD3 signaling in lung epithelial cells.MethodologyWe combined chromatin immunoprecipitation with human promoter region microarrays (ChIP-on-chip) along with gene expression microarrays to study global transcriptional regulation of the TGFβ1/SMAD3 pathway in human A549 alveolar epithelial cells. The molecular pathways and networks associated with TGFβ1/SMAD3 signaling were identified using computational approaches. Validation of selected target gene expression and direct binding of SMAD3 to promoters were performed by quantitative real time RT-PCR and electrophoretic mobility shift assay on A549 and human primary lung epithelial cells.Results and conclusionsKnown TGFβ1 target genes such as SERPINE1, SMAD6, SMAD7, TGFB1 and LTBP3, were found in both ChIP-on-chip and gene expression analyses as well as some previously unrecognized targets such as FOXA2. SMAD3 binding of FOXA2 promoter and changed expression were confirmed. Computational approaches combining ChIP-on-chip and gene expression microarray revealed multiple target molecular pathways affected by the TGFβ1/SMAD3 signaling. Identification of global targets and molecular pathways and networks associated with TGFβ1/SMAD3 signaling allow for a better understanding of the mechanisms that determine epithelial cell phenotypes in fibrogenesis and carcinogenesis as does the discovery of the direct effect of TGFβ1 on FOXA2.

Although the etiology of idiopathic pulmonary fibrosis (IPF) remains perplexing, adaptive immune activation is evident among many afflicted patients. Repeated cycles of antigen-induced proliferation cause T-cells to lose surface expression of CD28, and we hypothesized this process might also occur in IPF.Peripheral blood CD4 T-cells from 89 IPF patients were analyzed by flow cytometry and cytokine multiplex assays, and correlated with clinical events. In comparison to autologous CD4(+)CD28(+)cells, the unusual CD4(+)CD28(null) lymphocytes seen in many IPF patients had discordant expressions of activation markers, more frequently produced cytotoxic mediators perforin (2.4+/-0.8% vs. 60.0+/-7.4%, por=82%, compared to 9+/-9% among those with more extensive CD28 down-regulation (CD4(+)CD28(+)/CD4(total)

BACKGROUND:Chronic obstructive pulmonary disease (COPD) is characterized by both airway remodeling and parenchymal destruction. The identification of unique biomarker patterns associated with airway dominant versus parenchymal dominant patterns would support the existence of unique phenotypes representing independent biologic processes. A cross-sectional study was performed to examine the association of serum biomarkers with radiographic airway and parenchymal phenotypes of COPD. METHODOLOGY/PRINCIPAL FINDINGS:Serum from 234 subjects enrolled in a CT screening cohort was analyzed for 33 cytokines and growth factors using a multiplex protein array. The association of serum markers with forced expiratory volume in one second percent predicted (FEV1%) and quantitative CT measurements of airway thickening and emphysema was assessed with and without stratification for current smoking status. Significant associations were found with several serum inflammatory proteins and measurements of FEV1%, airway thickening, and parenchymal emphysema independent of smoking status. The association of select analytes with airway thickening and emphysema was independent of FEV1%. Furthermore, the relationship between other inflammatory markers and measurements of physiologic obstruction or airway thickening was dependent on current smoking status. CONCLUSIONS/SIGNIFICANCE:Airway and parenchymal phenotypes of COPD are associated with unique systemic serum biomarker profiles. Serum biomarker patterns may provide a more precise classification of the COPD syndrome, provide insights into disease pathogenesis and identify targets for novel patient-specific biological therapies.

Silicosis is a complex lung disease for which no successful treatment is available and therefore lung transplantation is a potential alternative. Tumor necrosis factor alpha (TNFalpha) plays a central role in the pathogenesis of silicosis. TNFalpha signaling is mediated by the transcription factor, Nuclear Factor (NF)-kappaB, which regulates genes controlling several physiological processes including the innate immune responses, cell death, and inflammation. Therefore, inhibition of NF-kappaB activation represents a potential therapeutic strategy for silicosis.In the present work we evaluated the lung transplant database (May 1986-July 2007) at the University of Pittsburgh to study the efficacy of lung transplantation in patients with silicosis (n = 11). We contrasted the overall survival and rate of graft rejection in these patients to that of patients with idiopathic pulmonary fibrosis (IPF, n = 79) that was selected as a control group because survival benefit of lung transplantation has been identified for these patients. At the time of lung transplantation, we found the lungs of silica-exposed subjects to contain multiple foci of inflammatory cells and silicotic nodules with proximal TNFalpha expressing macrophage and NF-kappaB activation in epithelial cells. Patients with silicosis had poor survival (median survival 2.4 yr; confidence interval (CI): 0.16-7.88 yr) compared to IPF patients (5.3 yr; CI: 2.8-15 yr; p = 0.07), and experienced early rejection of their lung grafts (0.9 yr; CI: 0.22-0.9 yr) following lung transplantation (2.4 yr; CI:1.5-3.6 yr; p

MicroRNAs (miRNAs) are short, non-coding RNA regulators of protein coding genes. miRNAs play a very important role in diverse biological processes and various diseases. Many algorithms are able to predict miRNA genes and their targets, but their transcription regulation is still under investigation. It is generally believed that intragenic miRNAs (located in introns or exons of protein coding genes) are co-transcribed with their host genes and most intergenic miRNAs transcribed from their own RNA polymerase II (Pol II) promoter. However, the length of the primary transcripts and promoter organization is currently unknown.We performed Pol II chromatin immunoprecipitation (ChIP)-chip using a custom array surrounding regions of known miRNA genes. To identify the true core transcription start sites of the miRNA genes we developed a new tool (CPPP). We showed that miRNA genes can be transcribed from promoters located several kilobases away and that their promoters share the same general features as those of protein coding genes. Finally, we found evidence that as many as 26% of the intragenic miRNAs may be transcribed from their own unique promoters.miRNA promoters have similar features to those of protein coding genes, but miRNA transcript organization is more complex.

Idiopathic pulmonary fibrosis (IPF) is a chronic progressive fibrotic lung disease associated with substantial morbidity and mortality. The objective of this study was to determine whether there is a peripheral blood protein signature in IPF and whether components of this signature may serve as biomarkers for disease presence and progression.We analyzed the concentrations of 49 proteins in the plasma of 74 patients with IPF and in the plasma of 53 control individuals. We identified a combinatorial signature of five proteins-MMP7, MMP1, MMP8, IGFBP1, and TNFRSF1A-that was sufficient to distinguish patients from controls with a sensitivity of 98.6% (95% confidence interval [CI] 92.7%-100%) and specificity of 98.1% (95% CI 89.9%-100%). Increases in MMP1 and MMP7 were also observed in lung tissue and bronchoalveolar lavage fluid obtained from IPF patients. MMP7 and MMP1 plasma concentrations were not increased in patients with chronic obstructive pulmonary disease or sarcoidosis and distinguished IPF compared to subacute/chronic hypersensitivity pneumonitis, a disease that may mimic IPF, with a sensitivity of 96.3% (95% CI 81.0%-100%) and specificity of 87.2% (95% CI 72.6%-95.7%). We verified our results in an independent validation cohort composed of patients with IPF, familial pulmonary fibrosis, subclinical interstitial lung disease (ILD), as well as with control individuals. MMP7 and MMP1 concentrations were significantly higher in IPF patients compared to controls in this cohort. Furthermore, MMP7 concentrations were elevated in patients with subclinical ILD and negatively correlated with percent predicted forced vital capacity (FVC%) and percent predicted carbon monoxide diffusing capacity (DLCO%).Our experiments provide the first evidence for a peripheral blood protein signature in IPF to our knowledge. The two main components of this signature, MMP7 and MMP1, are overexpressed in the lung microenvironment and distinguish IPF from other chronic lung diseases. Additionally, increased MMP7 concentration may be indicative of asymptomatic ILD and reflect disease progression.

Mechanical ventilation causes ventilator-induced lung injury in animals and humans. Mitogen-activated protein kinases have been implicated in ventilator-induced lung injury though their functional significance remains incomplete. We characterize the role of p38 mitogen-activated protein kinase/mitogen activated protein kinase kinase-3 and c-Jun-NH(2)-terminal kinase-1 in ventilator-induced lung injury and investigate novel independent mechanisms contributing to lung injury during mechanical ventilation.C57/BL6 wild-type mice and mice genetically deleted for mitogen-activated protein kinase kinase-3 (mkk-3(-/-)) or c-Jun-NH(2)-terminal kinase-1 (jnk1(-/-)) were ventilated, and lung injury parameters were assessed. We demonstrate that mkk3(-/-) or jnk1(-/-) mice displayed significantly reduced inflammatory lung injury and apoptosis relative to wild-type mice. Since jnk1(-/-) mice were highly resistant to ventilator-induced lung injury, we performed comprehensive gene expression profiling of ventilated wild-type or jnk1(-/-) mice to identify novel candidate genes which may play critical roles in the pathogenesis of ventilator-induced lung injury. Microarray analysis revealed many novel genes differentially expressed by ventilation including matrix metalloproteinase-8 (MMP8) and GADD45alpha. Functional characterization of MMP8 revealed that mmp8(-/-) mice were sensitized to ventilator-induced lung injury with increased lung vascular permeability.We demonstrate that mitogen-activated protein kinase pathways mediate inflammatory lung injury during ventilator-induced lung injury. C-Jun-NH(2)-terminal kinase was also involved in alveolo-capillary leakage and edema formation, whereas MMP8 inhibited alveolo-capillary protein leakage.

Idiopathic pulmonary fibrosis (IPF) is characterized by the insidious onset of dyspnea or cough. However, a subset of patients has a short duration of symptoms with rapid progression to end-stage disease. In this study, we evaluated clinical and molecular features of "rapid" and "slow" progressors with IPF.26 patients with 24 months of symptoms [slow progressors] were studied. Survival was analyzed by the Kaplan-Meyer method and proportional hazard's model. Lung microarrays and tissue proteins were measured in a subset of patients. No differences were found in age, physiologic impairment and bronchoalveolar lavage (BAL) cellular profile. There were more males (OR = 6.5; CI:1.4-29.5; p = 0.006) and smokers (OR = 3.04; CI:1.1-8.3; p = 0.04) in the rapid progressors group. Survival from the beginning of symptoms was significantly reduced in rapid progressors (HR = 9.0; CI:4.48-18.3; p2-fold increase of active matrix metalloproteinase-9, and induced a higher fibroblast migration compared with slow progressors and controls [238+/-98% versus 123+/-29% (p

Idiopathic pulmonary fibrosis (IPF) is a progressive and lethal disorder characterized by fibroproliferation and excessive accumulation of extracellular matrix in the lung.Using oligonucleotide arrays, we identified osteopontin as one of the genes that significantly distinguishes IPF from normal lungs. Osteopontin was localized to alveolar epithelial cells in IPF lungs and was also significantly elevated in bronchoalveolar lavage from IPF patients. To study the fibrosis-relevant effects of osteopontin we stimulated primary human lung fibroblasts and alveolar epithelial cells (A549) with recombinant osteopontin. Osteopontin induced a significant increase of migration and proliferation in both fibroblasts and epithelial cells. Epithelial growth was inhibited by the pentapeptide Gly-Arg-Gly-Asp-Ser (GRGDS) and antibody to CD44, while fibroproliferation was inhibited by GRGDS and antibody to alphavbeta3 integrin. Fibroblast and epithelial cell migration were inhibited by GRGDS, anti-CD44, and anti-alphavbeta3. In fibroblasts, osteopontin up-regulated tissue inhibitor of metalloprotease-1 and type I collagen, and down-regulated matrix metalloprotease-1 (MMP-1) expression, while in A549 cells it caused up-regulation of MMP-7. In human IPF lungs, osteopontin colocalized with MMP-7 in alveolar epithelial cells, and application of weakest link statistical models to microarray data suggested a significant interaction between osteopontin and MMP-7.Our results provide a potential mechanism by which osteopontin secreted from the alveolar epithelium may exert a profibrotic effect in IPF lungs and highlight osteopontin as a potential target for therapeutic intervention in this incurable disease.