High-resolution transcriptomic and epigenetic profiling identifies novel regulators of COPD.

Patients with chronic obstructive pulmonary disease (COPD) are still waiting for curative treatments. Considering its environmental cause, we hypothesized that COPD will be associated with altered epigenetic signaling in lung cells. We generated genome-wide DNA methylation maps at single CpG resolution of primary human lung fibroblasts (HLFs) across COPD stages. We show that the epigenetic landscape is changed early in COPD, with DNA methylation changes occurring predominantly in regulatory regions. RNA sequencing of matched fibroblasts demonstrated dysregulation of genes involved in proliferation, DNA repair, and extracellular matrix organization. Data integration identified 110 candidate regulators of disease phenotypes that were linked to fibroblast repair processes using phenotypic screens. Our study provides high-resolution multi-omic maps of HLFs across COPD stages. We reveal novel transcriptomic and epigenetic signatures associated with COPD onset and progression and identify new candidate regulators involved in the pathogenesis of chronic lung diseases. The presence of various epigenetic factors among the candidates demonstrates that epigenetic regulation in COPD is an exciting research field that holds promise for novel therapeutic avenues for patients.

Transcriptomic and proteomic profiling of young and old mice in the bleomycin model reveals high similarity.

The most common preclinical, in vivo model to study lung fibrosis is the bleomycin-induced lung fibrosis model in 2- to 3-mo-old mice. Although this model resembles key aspects of idiopathic pulmonary fibrosis (IPF), there are limitations in its predictability for the human disease. One of the main differences is the juvenile age of animals that are commonly used in experiments, resembling humans of around 20 yr. Because IPF patients are usually older than 60 yr, aging appears to play an important role in the pathogenesis of lung fibrosis. Therefore, we compared young (3 months) and old mice (21 months) 21 days after intratracheal bleomycin instillation. Analyzing lung transcriptomics (mRNAs and miRNAs) and proteomics, we found most pathways to be similarly regulated in young and old mice. However, old mice show imbalanced protein homeostasis as well as an increased inflammatory state in the fibrotic phase compared to young mice. Comparisons with published human transcriptomic data sets (GSE47460, GSE32537, and GSE24206) revealed that the gene signature of old animals correlates significantly better with IPF patients, and it also turned human healthy individuals better into "IPF patients" using an approach based on predictive disease modeling. Both young and old animals show similar molecular hallmarks of IPF in the bleomycin-induced lung fibrosis model, although old mice more closely resemble several features associated with IPF in comparison to young animals.

An EZH2-dependent transcriptional complex promotes aberrant epithelial remodelling after injury.

Unveiling the molecular mechanisms of tissue remodelling following injury is imperative to elucidate its regenerative capacity and aberrant repair in disease. Using different omics approaches, we identified enhancer of zester homolog 2 (EZH2) as a key regulator of fibrosis in injured lung epithelium. Epithelial injury drives an enrichment of nuclear transforming growth factor-ß-activated kinase 1 (TAK1) that mediates EZH2 phosphorylation to facilitate its liberation from polycomb repressive complex 2 (PRC2). This process results in the establishment of a transcriptional complex of EZH2, RNA-polymerase II (POL2) and nuclear actin, which orchestrates aberrant epithelial repair programmes. The liberation of EZH2 from PRC2 is accompanied by an EZH2-EZH1 switch to preserve H3K27me3 deposition at non-target genes. Loss of epithelial TAK1, EZH2 or blocking nuclear actin influx attenuates the fibrotic cascade and restores respiratory homeostasis. Accordingly, EZH2 inhibition significantly improves outcomes in a pulmonary fibrosis mouse model. Our results reveal an important non-canonical function of EZH2, paving the way for new therapeutic interventions in fibrotic lung diseases.

RORγt inhibitors block both IL-17 and IL-22 conferring a potential advantage over anti-IL-17 alone to treat severe asthma.

BACKGROUND: RORγt is a transcription factor that enables elaboration of Th17-associated cytokines (including IL-17 and IL-22) and is proposed as a pharmacological target for severe asthma. METHODS: IL-17 immunohistochemistry was performed in severe asthma bronchial biopsies (specificity confirmed with in situ hybridization). Primary human small airway epithelial cells in air liquid interface and primary bronchial smooth muscle cells were stimulated with recombinant human IL-17 and/or IL-22 and pro-inflammatory cytokines measured. Balb/c mice were challenged intratracheally with IL-17 and/or IL-22 and airway hyperreactivity, pro-inflammatory cytokines and airway neutrophilia measured. Balb/c mice were sensitized intraperitoneally and challenged intratracheally with house dust mite extract and the effect of either a RORγt inhibitor (BIX119) or an anti-IL-11 antibody assessed on airway hyperreactivity, pro-inflammatory cytokines and airway neutrophilia measured. RESULTS: We confirmed in severe asthma bronchial biopsies both the presence of IL-17-positive lymphocytes and that an IL-17 transcriptome profile in a severe asthma patient sub-population. Both IL-17 and IL-22 stimulated the release of pro-inflammatory cytokine and chemokine release from primary human lung cells and in mice. Furthermore, IL-22 in combination with IL-17, but neither alone, elicits airway hyperresponsiveness (AHR) in naïve mice. A RORγt inhibitor specifically blocked both IL-17 and IL-22, AHR and neutrophilia in a mouse house dust mite model unlike other registered or advanced pipeline modes of action. Full efficacy versus these parameters was associated with 90% inhibition of IL-17 and 50% inhibition of IL-22. In contrast, anti-IL-17 also blocked IL-17, but not IL-22, AHR or neutrophilia. Moreover, the deregulated genes in the lungs from these mice correlated well with deregulated genes from severe asthma biopsies suggesting that this model recapitulates significant severe asthma-relevant biology. Furthermore, these genes were reversed upon RORγt inhibition in the HDM model. Cell deconvolution suggested that the responsible cells were corticosteroid insensitive γδ-T-cells. CONCLUSION: These data strongly suggest that both IL-17 and IL-22 are required for Th2-low endotype associated biology and that a RORγt inhibitor may provide improved clinical benefit in a severe asthma sub-population of patients by blocking both IL-17 and IL-22 biology compared with blocking IL-17 alone.

Recapitulating idiopathic pulmonary fibrosis related alveolar epithelial dysfunction in a human iPSC-derived air-liquid interface model.

Idiopathic pulmonary fibrosis (IPF) is a fatal disease of unknown cause that is characterized by progressive fibrotic lung remodeling. An abnormal emergence of airway epithelial-like cells within the alveolar compartments of the lung, herein termed bronchiolization, is often observed in IPF. However, the origin of this dysfunctional distal lung epithelium remains unknown due to a lack of suitable human model systems. In this study, we established a human induced pluripotent stem cell (iPSC)-derived air-liquid interface (ALI) model of alveolar epithelial type II (ATII)-like cell differentiation that allows us to investigate alveolar epithelial progenitor cell differentiation in vitro. We treated this system with an IPF-relevant cocktail (IPF-RC) to mimic the pro-fibrotic cytokine milieu present in IPF lungs. Stimulation with IPF-RC during differentiation increases secretion of IPF biomarkers and RNA sequencing (RNA-seq) of these cultures reveals significant overlap with human IPF patient data. IPF-RC treatment further impairs ATII differentiation by driving a shift toward an airway epithelial-like expression signature, providing evidence that a pro-fibrotic cytokine environment can influence the proximo-distal differentiation pattern of human lung epithelial cells. In conclusion, we show for the first time, the establishment of a human model system that recapitulates aspects of IPF-associated bronchiolization of the lung epithelium in vitro.

Cell-specific expression of lung disease risk-related genes in the human small airway epithelium.

BACKGROUND: The human small airway epithelium (SAE) plays a central role in the early events in the pathogenesis of most inherited and acquired lung disorders. Little is known about the molecular phenotypes of the specific cell populations comprising the SAE in humans, and the contribution of SAE specific cell populations to the risk for lung diseases. METHODS: Drop-seq single-cell RNA-sequencing was used to characterize the transcriptome of single cells from human SAE of nonsmokers and smokers by bronchoscopic brushing. RESULTS: Eleven distinct cell populations were identified, including major and rare epithelial cells, and immune/inflammatory cells. There was cell type-specific expression of genes relevant to the risk of the inherited pulmonary disorders, genes associated with risk of chronic obstructive pulmonary disease and idiopathic pulmonary fibrosis and (non-mutated) driver genes for lung cancers. Cigarette smoking significantly altered the cell type-specific transcriptomes and disease risk-related genes. CONCLUSIONS: This data provides new insights into the possible contribution of specific lung cells to the pathogenesis of lung disorders.

BAL Cell Gene Expression Is Indicative of Outcome and Airway Basal Cell Involvement in Idiopathic Pulmonary Fibrosis.

RATIONALE: Idiopathic pulmonary fibrosis (IPF) is a fatal disease with a variable and unpredictable course. OBJECTIVES: To determine whether BAL cell gene expression is predictive of survival in IPF. METHODS: This retrospective study analyzed the BAL transcriptome of three independent IPF cohorts: Freiburg (Germany), Siena (Italy), and Leuven (Belgium) including 212 patients. BAL cells from 20 healthy volunteers, 26 patients with sarcoidosis stage III and IV, and 29 patients with chronic obstructive pulmonary disease were used as control subjects. Survival analysis was performed by Cox models and component-wise boosting. Presence of airway basal cells was tested by immunohistochemistry and flow cytometry. MEASUREMENTS AND MAIN RESULTS: A total of 1,582 genes were predictive of mortality in the IPF derivation cohort in univariate analyses adjusted for age and sex at false discovery rate less than 0.05. A nine-gene signature, derived from the discovery cohort (Freiburg), performed well in both replication cohorts, Siena (P < 0.0032) and Leuven (P = 0.0033). nCounter expression analysis confirmed the array results (P < 0.0001). The genes associated with mortality in BAL cells were significantly enriched for genes expressed in airway basal cells. Further analyses by gene expression, flow cytometry, and immunohistochemistry showed an increase in airway basal cells in BAL and tissues of IPF compared with control subjects, but not in chronic obstructive pulmonary disease or sarcoidosis. CONCLUSIONS: Our results identify and validate a BAL signature that predicts mortality in IPF and improves the accuracy of outcome prediction based on clinical parameters. The BAL signature associated with mortality unmasks a potential role for airway basal cells in IPF.

Characterization of an immortalized human small airway basal stem/progenitor cell line with airway region-specific differentiation capacity.

BACKGROUND: The pathology of chronic obstructive pulmonary disease (COPD), idiopathic pulmonary fibrosis (IPF) and most lung cancers involves the small airway epithelium (SAE), the single continuous layer of cells lining the airways ≥ 6th generations. The basal cells (BC) are the stem/progenitor cells of the SAE, responsible for the differentiation into intermediate cells and ciliated, club and mucous cells. To facilitate the study of the biology of the human SAE in health and disease, we immortalized and characterized a normal human SAE basal cell line. METHODS: Small airway basal cells were purified from brushed SAE of a healthy nonsmoker donor with a characteristic normal SAE transcriptome. The BC were immortalized by retrovirus-mediated telomerase reverse transcriptase (TERT) transduction and single cell drug selection. The resulting cell line (hSABCi-NS1.1) was characterized by RNAseq, TaqMan PCR, protein immunofluorescence, differentiation capacity on an air-liquid interface (ALI) culture, transepithelial electrical resistance (TEER), airway region-associated features and response to genetic modification with SPDEF. RESULTS: The hSABCi-NS1.1 single-clone-derived cell line continued to proliferate for > 200 doubling levels and > 70 passages, continuing to maintain basal cell features (TP63+, KRT5+). When cultured on ALI, hSABCi-NS1.1 cells consistently formed tight junctions and differentiated into ciliated, club (SCGB1A1+), mucous (MUC5AC+, MUC5B+), neuroendocrine (CHGA+), ionocyte (FOXI1+) and surfactant protein positive cells (SFTPA+, SFTPB+, SFTPD+), observations confirmed by RNAseq and TaqMan PCR. Annotation enrichment analysis showed that "cilium" and "immunity" were enriched in functions of the top-1500 up-regulated genes. RNAseq reads alignment corroborated expression of CD4, CD74 and MHC-II. Compared to the large airway cell line BCi-NS1.1, differentiated of hSABCi-NS1.1 cells on ALI were enriched with small airway epithelial genes, including surfactant protein genes, LTF and small airway development relevant transcription factors NKX2-1, GATA6, SOX9, HOPX, ID2 and ETV5. Lentivirus-mediated expression of SPDEF in hSABCi-NS1.1 cells induced secretory cell metaplasia, accompanied with characteristic COPD-associated SAE secretory cell changes, including up-regulation of MSMB, CEACAM5 and down-regulation of LTF. CONCLUSIONS: The immortalized hSABCi-NS1.1 cell line has diverse differentiation capacities and retains SAE features, which will be useful for understanding the biology of SAE, the pathogenesis of SAE-related diseases, and testing new pharmacologic agents.

Modeling Pulmonary Disease Pathways Using Recombinant Adeno-Associated Virus 6.2.

Viral vectors have been applied successfully to generate disease-related animal models and to functionally characterize target genes in vivo. However, broader application is still limited by complex vector production, biosafety requirements, and vector-mediated immunogenic responses, possibly interfering with disease-relevant pathways. Here, we describe adeno-associated virus (AAV) variant 6.2 as an ideal vector for lung delivery in mice, overcoming most of the aforementioned limitations. In a proof-of-concept study using AAV6.2 vectors expressing IL-13 and transforming growth factor-ß1 (TGF-ß1), we were able to induce hallmarks of severe asthma and pulmonary fibrosis, respectively. Phenotypic characterization and deep sequencing analysis of the AAV-IL-13 asthma model revealed a characteristic disease signature. Furthermore, suitability of the model for compound testing was also demonstrated by pharmacological intervention studies using an anti-IL-13 antibody and dexamethasone. Similarly, the AAV-TGF-ß1 fibrosis model showed several disease-like pathophenotypes monitored by micro-computed tomography imaging and lung function measurement. Most importantly, analyses using stuffer control vectors demonstrated that in contrast to a common adenovirus-5 vector, AAV6.2 vectors did not induce any measurable inflammation and therefore carry a lower risk of altering relevant readouts. In conclusion, we propose AAV6.2 as an ideal vector system for the functional characterization of target genes in the context of pulmonary diseases in mice.

The FAIR Cookbook - the essential resource for and by FAIR doers.

The notion that data should be Findable, Accessible, Interoperable and Reusable, according to the FAIR Principles, has become a global norm for good data stewardship and a prerequisite for reproducibility. Nowadays, FAIR guides data policy actions and professional practices in the public and private sectors. Despite such global endorsements, however, the FAIR Principles are aspirational, remaining elusive at best, and intimidating at worst. To address the lack of practical guidance, and help with capability gaps, we developed the FAIR Cookbook, an open, online resource of hands-on recipes for "FAIR doers" in the Life Sciences. Created by researchers and data managers professionals in academia, (bio)pharmaceutical companies and information service industries, the FAIR Cookbook covers the key steps in a FAIRification journey, the levels and indicators of FAIRness, the maturity model, the technologies, the tools and the standards available, as well as the skills required, and the challenges to achieve and improve data FAIRness. Part of the ELIXIR ecosystem, and recommended by funders, the FAIR Cookbook is open to contributions of new recipes.

Comparison of normalization methods for Illumina BeadChip HumanHT-12 v3.

BACKGROUND: Normalization of microarrays is a standard practice to account for and minimize effects which are not due to the controlled factors in an experiment. There is an overwhelming number of different methods that can be applied, none of which is ideally suited for all experimental designs. Thus, it is important to identify a normalization method appropriate for the experimental setup under consideration that is neither too negligent nor too stringent. Major aim is to derive optimal results from the underlying experiment. Comparisons of different normalization methods have already been conducted, none of which, to our knowledge, comparing more than a handful of methods. RESULTS: In the present study, 25 different ways of pre-processing Illumina Sentrix BeadChip array data are compared. Among others, methods provided by the BeadStudio software are taken into account. Looking at different statistical measures, we point out the ideal versus the actual observations. Additionally, we compare qRT-PCR measurements of transcripts from different ranges of expression intensities to the respective normalized values of the microarray data. Taking together all different kinds of measures, the ideal method for our dataset is identified. CONCLUSIONS: Pre-processing of microarray gene expression experiments has been shown to influence further downstream analysis to a great extent and thus has to be carefully chosen based on the design of the experiment. This study provides a recommendation for deciding which normalization method is best suited for a particular experimental setup.

A Novel Fibroblast Reporter Cell Line for in vitro Studies of Pulmonary Fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a fatal disease of the lower respiratory tract with restricted therapeutic options. Repetitive injury of the bronchoalveolar epithelium leads to activation of pulmonary fibroblasts, differentiation into myofibroblasts and excessive extracellular matrix (ECM) deposition resulting in aberrant wound repair. However, detailed molecular and cellular mechanisms underlying initiation and progression of fibrotic changes are still elusive. Here, we report the generation of a representative fibroblast reporter cell line (10-4A BFP ) to study pathophysiological mechanisms of IPF in high throughput or high resolution in vitro live cell assays. To this end, we immortalized primary fibroblasts isolated from the distal lung of Sprague-Dawley rats. Molecular and transcriptomic characterization identified clone 10-4A as a matrix fibroblast subpopulation. Mechanical or chemical stimulation induced a reversible fibrotic state comparable to effects observed in primary isolated fibroblasts. Finally, we generated a reporter cell line (10-4A BFP ) to express nuclear blue fluorescent protein (BFP) under the promotor of the myofibroblast marker alpha smooth muscle actin (Acta2) using CRISPR/Cas9 technology. We evaluated the suitability of 10-4A BFP as reporter tool in plate reader assays. In summary, the 10-4A BFP cell line provides a novel tool to study fibrotic processes in vitro to gain new insights into the cellular and molecular processes involved in fibrosis formation and propagation.

Dysregulation of club cell biology in idiopathic pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a progressive, chronic fibrotic lung disease with an irreversible decline of lung function. "Bronchiolization", characterized by ectopic appearance of airway epithelial cells in the alveolar regions, is one of the characteristic features in the IPF lung. Based on the knowledge that club cells are the major epithelial secretory cells in human small airways, and their major secretory product uteroglobin (SCGB1A1) is significantly increased in both serum and epithelial lining fluid of IPF lung, we hypothesize that human airway club cells contribute to the pathogenesis of IPF. By assessing the transcriptomes of the single cells from human lung of control donors and IPF patients, we identified two SCGB1A1+ club cell subpopulations, highly expressing MUC5B, a significant genetic risk factor strongly associated with IPF, and SCGB3A2, a marker heterogeneously expressed in the club cells, respectively. Interestingly, the cellular proportion of SCGB1A1+MUC5B+ club cells was significantly increased in IPF patients, and this club cell subpopulation highly expressed genes related to mucous production and immune cell chemotaxis. In contrast, though the cellular proportion did not change, the molecular phenotype of the SCGB1A1+SCGB3A2high club cell subpopulation was significantly altered in IPF lung, with increased expression of mucins, cytokine and extracellular matrix genes. The single cell transcriptomic analysis reveals the cellular and molecular heterogeneity of club cells, and provide novel insights into the biological functions of club cells in the pathogenesis of IPF.

Intermittent exposure to whole cigarette smoke alters the differentiation of primary small airway epithelial cells in the air-liquid interface culture.

Cigarette smoke (CS) is the leading risk factor to develop COPD. Therefore, the pathologic effects of whole CS on the differentiation of primary small airway epithelial cells (SAEC) were investigated, using cells from three healthy donors and three COPD patients, cultured under ALI (air-liquid interface) conditions. The analysis of the epithelial physiology demonstrated that CS impaired barrier formation and reduced cilia beat activity. Although, COPD-derived ALI cultures preserved some features known from COPD patients, CS-induced effects were similarly pronounced in ALI cultures from patients compared to healthy controls. RNA sequencing analyses revealed the deregulation of marker genes for basal and secretory cells upon CS exposure. The comparison between gene signatures obtained from the in vitro model (CS vs. air) with a published data set from human epithelial brushes (smoker vs. non-smoker) revealed a high degree of similarity between deregulated genes and pathways induced by CS. Taken together, whole cigarette smoke alters the differentiation of small airway basal cells in vitro. The established model showed a good translatability to the situation in vivo. Thus, the model can help to identify and test novel therapeutic approaches to restore the impaired epithelial repair mechanisms in COPD, which is still a high medical need.

Effect of RNA quality on transcript intensity levels in microarray analysis of human post-mortem brain tissues.

BACKGROUND: Large-scale gene expression analysis of post-mortem brain tissue offers unique opportunities for investigating genetic mechanisms of psychiatric and neurodegenerative disorders. On the other hand microarray data analysis associated with these studies is a challenging task. In this publication we address the issue of low RNA quality data and corresponding data analysis strategies. RESULTS: A detailed analysis of effects of post chip RNA quality on the measured abundance of transcripts is presented. Overall Affymetrix GeneChip data (HG-U133_AB and HG-U133_Plus_2.0) derived from ten different brain regions was investigated. Post chip RNA quality being assessed by 5'/3' ratio of housekeeping genes was found to introduce a well pronounced systematic noise into the measured transcript expression levels. According to this study RNA quality effects have: 1) a "random" component which is introduced by the technology and 2) a systematic component which depends on the features of the transcripts and probes. Random components mainly account for numerous negative correlations of low-abundant transcripts. These negative correlations are not reproducible and are mainly introduced by an increased relative level of noise. Three major contributors to the systematic noise component were identified: the first is the probe set distribution, the second is the length of mRNA species, and the third is the stability of mRNA species. Positive correlations reflect the 5'-end to 3'-end direction of mRNA degradation whereas negative correlations result from the compensatory increase in stable and 3'-end probed transcripts. Systematic components affect the expressed transcripts by introducing irrelevant gene correlations and can strongly influence the results of the main experiment. A linear model correcting the effect of RNA quality on measured intensities was introduced. In addition the contribution of a number of pre-mortem and post-mortem attributes to the overall detected RNA quality effect was investigated. Brain pH, duration of agonal stage, post-mortem interval before sampling and donor's age of death within considered limits were found to have no significant contribution. CONCLUSION: Basic conclusions for data analysis in expression profiling study are as follows: 1) testing for RNA quality dependency should be included in the preprocessing of the data; 2) investigating inter-gene correlation without regard to RNA quality effects could be misleading; 3) data normalization procedures relying on housekeeping genes either do not influence the correlation structure (if 3'-end intensities are used) or increase it for negatively correlated transcripts (if 5'-end or median intensities are included in normalization procedure); 4) sample sets should be matched with regard to RNA quality; 5) RMA preprocessing is more sensitive to RNA quality effect, than MAS 5.0.

Key factors for successful data integration in biomarker research.

Integrating a wide range of biomedical data such as that rapidly emerging from the use of next-generation sequencing is expected to have a key role in identifying and qualifying new biomarkers to support precision medicine. Here, we highlight some of the challenges for biomedical data integration and approaches to address them.

Phenocopy--a strategy to qualify chemical compounds during hit-to-lead and/or lead optimization.

A phenocopy is defined as an environmentally induced phenotype of one individual which is identical to the genotype-determined phenotype of another individual. The phenocopy phenomenon has been translated to the drug discovery process as phenotypes produced by the treatment of biological systems with new chemical entities (NCE) may resemble environmentally induced phenotypic modifications. Various new chemical entities exerting inhibition of the kinase activity of Transforming Growth Factor ß Receptor I (TGF-ßR1) were qualified by high-throughput RNA expression profiling. This chemical genomics approach resulted in a precise time-dependent insight to the TGF-ß biology and allowed furthermore a comprehensive analysis of each NCE's off-target effects. The evaluation of off-target effects by the phenocopy approach allows a more accurate and integrated view on optimized compounds, supplementing classical biological evaluation parameters such as potency and selectivity. It has therefore the potential to become a novel method for ranking compounds during various drug discovery phases.

A renewed approach to the nonparametric analysis of replicated microarray experiments.

DNA-microarrays find broad employment in biochemical research. This technology allows the monitoring of the expression levels of thousands of genes at the same time. Often, the goal of a microarray study is to find differentially expressed genes in two different types of tissue, for example normal and cancerous. Multiple hypothesis testing is a useful statistical tool for such studies. One approach using multiple hypothesis testing is nonparametric analysis for replicated microarray experiments. In this paper we present an improved version of this method. We also show how p-values are calculated for all significant genes detected with this testing procedure. All algorithms were implemented in an R-package, and instructions on it's use are included. The package can be downloaded at http://www.statistik.unidortmund.de/de/content/einrichtungen/lehrstuehle/personen/jung.html