Nasal epithelium as a proxy for bronchial epithelium for smoking-induced gene expression and expression Quantitative Trait Loci

Improved Annotation of Asthma Gene Variants with Cell Type Deconvolution of Nasal and Lung Expression Quantitative Trait Loci

Asthma is a genetically complex inflammatory airway disease associated with more than 200 SNPs. However, the functional effects of many asthma-associated SNPs in lung and airway epithelial samples are unknown. Here, we aimed to conduct expression quantitative trait loci (eQTL) analysis using a meta-analysis of nasal and lung samples. We hypothesize that incorporating cell type proportions of airway and lung samples enhances eQTL analysis outcomes. Nasal brush (n = 792) and lung tissue (n = 1,087) samples were investigated separately. Initially, a general eQTL analysis identified genetic variants associated with gene expression levels. Estimated cell type proportions were adjusted based on the Human Lung Cell Atlas. In addition, the presence of significant interaction effects between asthma-associated SNPs and each cell type proportion was explored and considered evidence for cell type-associated eQTL. In nasal brush and lung parenchyma samples, 44 and 116 asthma-associated SNPs were identified as eQTL. Adjusting for cell type proportions revealed eQTL for an additional 17 genes (e.g., FCER1G, CD200R1, and GABBR2) and 16 genes (e.g., CYP2C8, SLC9A2, and SGCD) in nose and lung, respectively. Moreover, we identified eQTL for nine SNPs annotated to genes such as VASP, FOXA3, and PCDHB12 displayed significant interactions with cell type proportions of club, goblet, and alveolar macrophages. Our findings demonstrate increased power for identifying eQTL among asthma-associated SNPs by considering cell type proportion of the bulk RNA-sequencing data from nasal and lung tissues. Integration of cell type deconvolution and eQTL analysis enhances our understanding of asthma genetics and cellular mechanisms, uncovering potential therapeutic targets for personalized interventions.

Nasal epithelial gene expression identifies relevant asthma endotypes in the ATLANTIS study

INTRODUCTION

Asthma is an inflammatory airways disease encompassing multiple phenotypes and endotypes. Several studies suggested gene expression in nasal epithelium to serve as a proxy for bronchial epithelium, being a non-invasive approach to investigate lung diseases. We hypothesised that molecular differences in upper airway epithelium reflect asthma-associated differences in the lower airways and are associated with clinical expression of asthma.

METHODS

We analysed nasal epithelial gene expression data from 369 patients with asthma and 58 non-asthmatic controls from the Assessment of Small Airways Involvement in Asthma study. Unsupervised hierarchical clustering was performed on asthma-associated genes. Asthma-associated gene signatures were replicated in independent cohorts with nasal and bronchial brushes data by comparing Gene Set Variation Analysis scores between asthma patients and non-asthmatic controls.

RESULTS

We identified 67 higher expressed and 59 lower expressed genes in nasal epithelium from asthma patients compared with controls (false discovery rate<0.05), including CLCA1, CST1 and POSTN, genes well known to reflect asthma in bronchial airway epithelium. Hierarchical clustering revealed several molecular asthma endotypes with distinct clinical characteristics, including an endotype with higher blood and sputum eosinophils, high fractional exhaled nitric oxide, and more severe small airway dysfunction, as reflected by lower forced expiratory flow at 50%. In an independent cohort, we demonstrated that genes higher expressed in the nasal epithelium reflect asthma-associated changes in the lower airways.

CONCLUSION

Our results show that the nasal epithelial gene expression profile reflects asthma-related processes in the lower airways. We suggest that nasal epithelium may be a useful non-invasive tool to identify asthma endotypes and may advance personalised management of the disease.

A bronchial gene signature specific for severe COPD that is retained in the nose

Introduction

A subset of COPD patients develops advanced disease with severe airflow obstruction, hyperinflation and extensive emphysema. We propose that the pathogenesis in these patients differs from mild-moderate COPD and is reflected by bronchial gene expression. The aim of the present study was to identify a unique bronchial epithelial gene signature for severe COPD patients.

Methods

We obtained RNA sequencing data from bronchial brushes from 123 ex-smokers with severe COPD, 23 with mild-moderate COPD and 23 non-COPD controls. We identified genes specific to severe COPD by comparing severe COPD to non-COPD controls, followed by removing genes that were also differentially expressed between mild-moderate COPD and non-COPD controls. Next, we performed a pathway analysis on these genes and evaluated whether this signature is retained in matched nasal brushings.

Results

We identified 219 genes uniquely differentially expressed in severe COPD. Interaction network analysis identified VEGFA and FN1 as the key genes with the most interactions. Genes were involved in extracellular matrix regulation, collagen binding and the immune response. Of interest were 10 genes (VEGFA, DCN, SPARC, COL6A2, MGP, CYR61, ANXA6, LGALS1, C1QA and C1QB) directly connected to fibronectin 1 (FN1). Most of these genes were lower expressed in severe COPD and showed the same effect in nasal brushings.

Conclusions

We found a unique severe COPD bronchial gene signature with key roles for VEGFA and FN1, which was retained in the upper airways. This supports the hypothesis that severe COPD, at least partly, comprises a different pathology and supports the potential for biomarker development based on nasal brushes in COPD.

The role of DNA methylation in personalized medicine for immune-related diseases

Epigenetics functions as a bridge between host genetic & environmental factors, aiding in human health and diseases. Many immune-related diseases, including infectious and allergic diseases, have been linked to epigenetic mechanisms, particularly DNA methylation. In this review, we summarized an updated overview of DNA methylation and its importance in personalized medicine, and demonstrated that DNA methylation has excellent potential for disease prevention, diagnosis, and treatment in a personalized manner. The future implications and limitations of the DNA methylation study have also been well-discussed.

Smoking induces shifts in cellular composition and transcriptome within the bronchial mucus barrier

BACKGROUND AND OBJECTIVE

Smoking disturbs the bronchial-mucus-barrier. This study assesses the cellular composition and gene expression shifts of the bronchial-mucus-barrier with smoking to understand the mechanism of mucosal damage by cigarette smoke exposure. We explore whether single-cell-RNA-sequencing (scRNA-seq) based cellular deconvolution (CD) can predict cell-type composition in RNA-seq data.

METHODS

RNA-seq data of bronchial biopsies from three cohorts were analysed using CD. The cohorts included 56 participants with chronic obstructive pulmonary disease [COPD] (38 smokers; 18 ex-smokers), 77 participants without COPD (40 never-smokers; 37 smokers) and 16 participants who stopped smoking for 1 year (11 COPD and 5 non-COPD-smokers). Differential gene expression was used to investigate gene expression shifts. The CD-derived goblet cell ratios were validated by correlating with staining-derived goblet cell ratios from the COPD cohort. Statistics were done in the R software (false discovery rate p-value < 0.05).

RESULTS

Both CD methods indicate a shift in bronchial-mucus-barrier cell composition towards goblet cells in COPD and non-COPD-smokers compared to ex- and never-smokers. It shows that the effect was reversible within a year of smoking cessation. A reduction of ciliated and basal cells was observed with current smoking, which resolved following smoking cessation. The expression of mucin and sodium channel (ENaC) genes, but not chloride channel genes, were altered in COPD and current smokers compared to never smokers or ex-smokers. The goblet cell-derived staining scores correlate with CD-derived goblet cell ratios.

CONCLUSION

Smoking alters bronchial-mucus-barrier cell composition, transcriptome and increases mucus production. This effect is partly reversible within a year of smoking cessation. CD methodology can predict goblet-cell percentages from RNA-seq.

Multi-omic association study identifies DNA methylation-mediated genotype and smoking exposure effects on lung function in children living in urban settings

Impaired lung function in early life is associated with the subsequent development of chronic respiratory disease. Most genetic associations with lung function have been identified in adults of European descent and therefore may not represent those most relevant to pediatric populations and populations of different ancestries. In this study, we performed genome-wide association analyses of lung function in a multiethnic cohort of children (n = 1,035) living in low-income urban neighborhoods. We identified one novel locus at the TDRD9 gene in chromosome 14q32.33 associated with percent predicted forced expiratory volume in one second (FEV1) (p = 2.4x10-9; βz = -0.31, 95% CI = -0.41- -0.21). Mendelian randomization and mediation analyses revealed that this genetic effect on FEV1 was partially mediated by DNA methylation levels at this locus in airway epithelial cells, which were also associated with environmental tobacco smoke exposure (p = 0.015). Promoter-enhancer interactions in airway epithelial cells revealed chromatin interaction loops between FEV1-associated variants in TDRD9 and the promoter region of the PPP1R13B gene, a stimulator of p53-mediated apoptosis. Expression of PPP1R13B in airway epithelial cells was significantly associated the FEV1 risk alleles (p = 1.3x10-5; β = 0.12, 95% CI = 0.06-0.17). These combined results highlight a potential novel mechanism for reduced lung function in urban youth resulting from both genetics and smoking exposure.

The establishment of COPD organoids to study host-pathogen interaction reveals enhanced viral fitness of SARS-CoV-2 in bronchi

Chronic obstructive pulmonary disease (COPD) is characterised by airflow limitation and infective exacerbations, however, in-vitro model systems for the study of host-pathogen interaction at the individual level are lacking. Here, we describe the establishment of nasopharyngeal and bronchial organoids from healthy individuals and COPD that recapitulate disease at the individual level. In contrast to healthy organoids, goblet cell hyperplasia and reduced ciliary beat frequency were observed in COPD organoids, hallmark features of the disease. Single-cell transcriptomics uncovered evidence for altered cellular differentiation trajectories in COPD organoids. SARS-CoV-2 infection of COPD organoids revealed more productive replication in bronchi, the key site of infection in severe COVID-19. Viral and bacterial exposure of organoids induced greater pro-inflammatory responses in COPD organoids. In summary, we present an organoid model that recapitulates the in vivo physiological lung microenvironment at the individual level and is amenable to the study of host-pathogen interaction and emerging infectious disease.

Organotypic sinonasal airway culture systems are predictive of the mucociliary phenotype produced by bronchial airway epithelial cells

Differentiated air-liquid interface models are the current standard to assess the mucociliary phenotype using clinically-derived samples in a controlled environment. However, obtaining basal progenitor airway epithelial cells (AEC) from the lungs is invasive and resource-intensive. Hence, we applied a tissue engineering approach to generate organotypic sinonasal AEC (nAEC) epithelia to determine whether they are predictive of bronchial AEC (bAEC) models. Basal progenitor AEC were isolated from healthy participants using a cytological brushing method and differentiated into epithelia on transwells until the mucociliary phenotype was observed. Tissue architecture was assessed using H&E and alcian blue/Verhoeff-Van Gieson staining, immunofluorescence (for cilia via acetylated α-tubulin labelling) and scanning electron microscopy. Differentiation and the formation of tight-junctions were monitored over the culture period (day 1-32) by quantifying trans-epithelial electrical resistance. End point (day 32) tight junction protein expression was assessed using Western blot analysis of ZO-1, Occludin-1 and Claudin-1. Reverse transcription qPCR-array was used to assess immunomodulatory and autophagy-specific transcript profiles. All outcome measures were assessed using R-statistical software. Mucociliary architecture was comparable for nAEC and bAEC-derived cultures, e.g. cell density P = 0.55, epithelial height P = 0.88 and cilia abundance P = 0.41. Trans-epithelial electrical resistance measures were distinct from day 1-14, converged over days 16-32, and were statistically similar over the entire culture period (global P < 0.001). This agreed with end-point (day 32) measures of tight junction protein abundance which were non-significant for each analyte (P > 0.05). Transcript analysis for inflammatory markers demonstrated significant variation between nAEC and bAEC epithelial cultures, and favoured increased abundance in the nAEC model (e.g. TGFβ and IL-1β; P < 0.05). Conversely, the abundance of autophagy-related transcripts were comparable and the range of outcome measures for either model exhibited a considerably more confined uncertainty distribution than those observed for the inflammatory markers. Organotypic air-liquid interface models of nAEC are predictive of outcomes related to barrier function, mucociliary architecture and autophagy gene activity in corresponding bAEC models. However, inflammatory markers exhibited wide variation which may be explained by the sentinel immunological surveillance role of the sinonasal epithelium.

Lipid mediators are detectable in the nasal epithelium and differ by asthma status in female subjects

BACKGROUND

Lipid mediators, bioactive products of polyunsaturated fatty acid metabolism, contribute to inflammation initiation and resolution in allergic diseases; however, their presence in lung-related biosamples has not been fully described.

OBJECTIVE

We aimed to quantify lipid mediators in the nasal airway epithelium and characterize preliminary associations with asthma.

METHODS

Using liquid chromatography-mass spectrometry, we conducted a pilot study to quantify 56 lipid mediators from nasal epithelial samples collected from 11 female participants of an outpatient asthma clinic and community controls (aged 30-55 years). We examined the presence of each compound using descriptive statistics to test whether lipid mediators could distinguish subjects with asthma (n = 8) from control subjects (n = 3) using linear regression and partial least squares discriminant analysis.

RESULTS

Fifteen lipid mediators were detectable in all samples, including resolvin (Rv) D5 (RvD5), with the highest median concentrations (in pg/μg protein) of 13-HODE (126.481), 15-HETE (32.869), and 13-OxoODE (13.251). From linear regression adjusted for age, prostaglandin E2 (PGE2) had a trend (P < .1) for higher concentrations in patients with severe asthma compared to controls (mean difference, 0.95; 95% confidence interval, -0.04 to 1.95). Asthma patients had higher scores on principal component 3 compared to controls (mean difference, 2.42; 95% confidence interval, 0.89 to 3.96), which represented lower levels of proresolving 15-HEPE, 19,20-DiHDPA, RvD5, 14-HDHA, 17-HDHA, and 13-HOTrE. Most of these compounds were best at discriminating asthma cases from controls in partial least squares discriminant analysis.

CONCLUSION

Lipid mediators are detectable in the nasal epithelium, and their levels distinguish asthma cases from controls.

Increased SARS-CoV-2 Infection, Protease, and Inflammatory Responses in Chronic Obstructive Pulmonary Disease Primary Bronchial Epithelial Cells Defined with Single-Cell RNA Sequencing

Rationale: Patients with chronic obstructive pulmonary disease (COPD) develop more severe coronavirus disease (COVID-19); however, it is unclear whether they are more susceptible to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and what mechanisms are responsible for severe disease. Objectives: To determine whether SARS-CoV-2 inoculated primary bronchial epithelial cells (pBECs) from patients with COPD support greater infection and elucidate the effects and mechanisms involved. Methods: We performed single-cell RNA sequencing analysis on differentiated pBECs from healthy subjects and patients with COPD 7 days after SARS-CoV-2 inoculation. We correlated changes with viral titers, proinflammatory responses, and IFN production. Measurements and Main Results: Single-cell RNA sequencing revealed that COPD pBECs had 24-fold greater infection than healthy cells, which was supported by plaque assays. Club/goblet and basal cells were the predominant populations infected and expressed mRNAs involved in viral replication. Proteases involved in SARS-CoV-2 entry/infection (TMPRSS2 and CTSB) were increased, and protease inhibitors (serpins) were downregulated more so in COPD. Inflammatory cytokines linked to COPD exacerbations and severe COVID-19 were increased, whereas IFN responses were blunted. Coexpression analysis revealed a prominent population of club/goblet cells with high type 1/2 IFN responses that were important drivers of immune responses to infection in both healthy and COPD pBECs. Therapeutic inhibition of proteases and inflammatory imbalances reduced viral titers and cytokine responses, particularly in COPD pBECs. Conclusions: COPD pBECs are more susceptible to SARS-CoV-2 infection because of increases in coreceptor expression and protease imbalances and have greater inflammatory responses. A prominent cluster of IFN-responsive club/goblet cells emerges during infection, which may be important drivers of immunity. Therapeutic interventions suppress SARS-CoV-2 replication and consequent inflammation.

Early transcriptional responses of bronchial epithelial cells to whole cigarette smoke mirror those of in-vivo exposed human bronchial mucosa

BACKGROUND

Despite the well-known detrimental effects of cigarette smoke (CS), little is known about the complex gene expression dynamics in the early stages after exposure. This study aims to investigate early transcriptomic responses following CS exposure of airway epithelial cells in culture and compare these to those found in human CS exposure studies.

METHODS

Primary bronchial epithelial cells (PBEC) were differentiated at the air-liquid interface (ALI) and exposed to whole CS. Bulk RNA-sequencing was performed at 1 h, 4 h, and 24 h hereafter, followed by differential gene expression analysis. Results were additionally compared to data retrieved from human CS studies.

RESULTS

ALI-PBEC gene expression in response to CS was most significantly changed at 4 h after exposure. Early transcriptomic changes (1 h, 4 h post CS exposure) were related to oxidative stress, xenobiotic metabolism, higher expression of immediate early genes and pro-inflammatory pathways (i.e., Nrf2, AP-1, AhR). At 24 h, ferroptosis-associated genes were significantly increased, whereas PRKN, involved in removing dysfunctional mitochondria, was downregulated. Importantly, the transcriptome dynamics of the current study mirrored in-vivo human studies of acute CS exposure, chronic smokers, and inversely mirrored smoking cessation.

CONCLUSION

These findings show that early after CS exposure xenobiotic metabolism and pro-inflammatory pathways were activated, followed by activation of the ferroptosis-related cell death pathway. Moreover, significant overlap between these transcriptomic responses in the in-vitro model and human in-vivo studies was found, with an early response of ciliated cells. These results provide validation for the use of ALI-PBEC cultures to study the human lung epithelial response to inhaled toxicants.

Multi-omics colocalization with genome-wide association studies reveals a context-specific genetic mechanism at a childhood onset asthma risk locus

BACKGROUND

Genome-wide association studies (GWASs) have identified thousands of variants associated with asthma and other complex diseases. However, the functional effects of most of these variants are unknown. Moreover, GWASs do not provide context-specific information on cell types or environmental factors that affect specific disease risks and outcomes. To address these limitations, we used an upper airway epithelial cell (AEC) culture model to assess transcriptional and epigenetic responses to rhinovirus (RV), an asthma-promoting pathogen, and provide context-specific functional annotations to variants discovered in GWASs of asthma.

METHODS

Genome-wide genetic, gene expression, and DNA methylation data in vehicle- and RV-treated upper AECs were collected from 104 individuals who had a diagnosis of airway disease (n=66) or were healthy participants (n=38). We mapped cis expression and methylation quantitative trait loci (cis-eQTLs and cis-meQTLs, respectively) in each treatment condition (RV and vehicle) in AECs from these individuals. A Bayesian test for colocalization between AEC molecular QTLs and adult onset asthma and childhood onset asthma GWAS SNPs, and a multi-ethnic GWAS of asthma, was used to assign the function to variants associated with asthma. We used Mendelian randomization to demonstrate DNA methylation effects on gene expression at asthma colocalized loci.

RESULTS

Asthma and allergic disease-associated GWAS SNPs were specifically enriched among molecular QTLs in AECs, but not in GWASs from non-immune diseases, and in AEC eQTLs, but not among eQTLs from other tissues. Colocalization analyses of AEC QTLs with asthma GWAS variants revealed potential molecular mechanisms of asthma, including QTLs at the TSLP locus that were common to both the RV and vehicle treatments and to both childhood onset and adult onset asthma, as well as QTLs at the 17q12-21 asthma locus that were specific to RV exposure and childhood onset asthma, consistent with clinical and epidemiological studies of these loci.

CONCLUSIONS

This study provides evidence of functional effects for asthma risk variants in AECs and insight into RV-mediated transcriptional and epigenetic response mechanisms that modulate genetic effects in the airway and risk for asthma.

Stability in eosinophil categorisation during subsequent severe exacerbations of COPD

Background

The blood eosinophil count has been shown to be a promising biomarker for establishing personalised treatment strategies to reduce corticosteroid use, either inhaled or systemic, in chronic obstructive pulmonary disease (COPD). Eosinophil levels seem relatively stable over time in stable state, but little is known whether this is also true in subsequent severe acute exacerbations of COPD (AECOPD).

Aims and objectives

To determine the stability in eosinophil categorisation between two subsequent severe AECOPDs employing frequently used cut-off levels.

Methods

During two subsequent severe AECOPDs, blood eosinophil counts were determined at admission to the hospital in 237 patients in the Cohort of Mortality and Inflammation in COPD Study. The following four cut-off levels were analysed: absolute counts of eosinophils ≥0.2×10⁹/L (200 cells/µL) and ≥0.3×10⁹/L (300 cells/µL) and relative eosinophil percentage of ≥2% and ≥3% of total leucocyte count. Categorisations were considered stable if during the second AECOPD their blood eosinophil status led to the same classification: eosinophilic or not.

Results

Depending on the used cut-off, the overall stability in eosinophil categorisation varied between 70% and 85% during two subsequent AECOPDs. From patients who were eosinophilic at the first AECOPD, 34%–45% remained eosinophilic at the subsequent AECOPD, while 9%–21% of patients being non-eosinophilic at the first AECOPD became eosinophilic at the subsequent AECOPD.

Conclusions

The eosinophil variability leads to category changes in subsequent AECOPDs, which limits the eosinophil categorisation stability. Therefore, measurement of eosinophils at each new exacerbation seems warranted.

Single-nucleotide polymorphism rs2070600 regulates AGER splicing and the sputum levels of the COPD biomarker soluble receptor for advanced glycation end-products

The COPD susceptibility SNP rs2070600 affects the levels of the COPD biomarker sRAGE in sputum as well as splicing of AGER. Moreover, @PouwelsScience et al. demonstrate large differences in sRAGE levels between serum and sputum. https://bit.ly/3t0pJtK.

Current-Smoking alters Gene Expression and DNA Methylation in the Nasal Epithelium of Asthmatics

Current-smoking contributes to worsened asthma prognosis, more severe symptoms and limits the beneficial effects of corticosteroids. As the nasal epithelium can reflect smoking-induced changes in the lower airways, it is a relevant source to investigate changes in gene expression and DNA methylation. This study explores gene expression and DNA methylation changes in current and ex-smokers with asthma. Matched gene expression and epigenome-wide DNA methylation samples collected from nasal brushings of 55 patients enrolled in a clinical trial investigation of current and ex-smoker asthma patients were analysed. Differential gene expression and DNA methylation analyses were conducted comparing current- vs ex-smokers. Expression quantitative trait methylation (eQTM) analysis was completed to explore smoking relevant genes by CpG sites that differ between current and ex-smokers. To investigate the relevance of the smoking-associated DNA methylation changes for the lower airways, significant CpG sites were explored in bronchial biopsies from patients who had stopped smoking. 809 genes and 18,814 CpG sites were differentially associated with current-smoking in the nose. The cis-eQTM analysis uncovered 171 CpG sites whose methylation status associated with smoking-related gene expression, including AHRR, ALDH3A1, CYP1A1 and CYP1B1. Methylation status of CpG sites altered by current-smoking reversed with one-year smoking cessation. We confirm current-smoking alters epigenetic patterns and affects gene expression in the nasal epithelium of asthma patients, which is partially reversible in bronhcial biopsies after smoking cessation. We demonstrate the ability to discern molecular changes in the nasal epithelium, presenting this as a tool in future investigations into disease-relevant effects of tobacco smoke.

Significant functional differences in differentiated Conditionally Reprogrammed (CRC)- and Feeder-free Dual SMAD inhibited-expanded human nasal epithelial cells

BACKGROUND

Patient-derived airway cells differentiated at Air Liquid Interface (ALI) are valuable models for Cystic fibrosis (CF) precision therapy. Different culture expansion methods have been established to extend expansion capacity of airway basal cells, while retaining functional airway epithelium physiology. Considerable variation in response to CFTR modulators is observed in cultures even within the same CFTR genotype and despite the use of similar ALI culture techniques. We aimed to address culture expansion method impact on differentiation.

METHODS

Nasal epithelial brushings from 14 individuals (CF=9; non-CF=5) were collected, then equally divided and expanded under conditional reprogramming culture (CRC) and feeder-serum-free "dual-SMAD inhibition" (SMADi) methods. Expanded cells from each culture were differentiated with proprietary PneumaCult™-ALI media. Morphology (Immunofluorescence), global proteomics (LC-MS/MS) and function (barrier integrity, cilia motility, and ion transport) were compared in CRC^ALI and SMADi^ALI under basal and CFTR corrector treated (VX-809) conditions.

RESULTS

No significant difference in the structural morphology or baseline global proteomics profile were observed. Barrier integrity and cilia motility were significantly different, despite no difference in cell junction morphology or cilia abundance. Epithelial Sodium Channels and Calcium-activated Chloride Channel activity did not differ but CFTR mediated chloride currents were significantly reduced in SMADi^ALI compare to their CRC^ALI counterparts.

CONCLUSION

Alteration of cellular physiological function in vitro were more prominent than structural and differentiation potential in airway ALI. Since initial expansion culture conditions significantly influence CFTR activity, this could lead to false conclusions if data from different labs are compared against each other without specific reference ranges.

Acute cigarette smoke-induced eQTL affects formyl peptide receptor expression and lung function

Background and objective

Cigarette smoking is one of the most prevalent causes of preventable deaths worldwide, leading to chronic diseases, including chronic obstructive pulmonary disease (COPD). Cigarette smoke is known to induce significant transcriptional modifications throughout the respiratory tract. However, it is largely unknown how genetic profiles influence the smoking‐related transcriptional changes and how changes in gene expression translate into altered alveolar epithelial repair responses.

Methods

We performed a candidate‐based acute cigarette smoke‐induced eQTL study, investigating the association between SNP and differential gene expression of FPR family members in bronchial epithelial cells isolated 24 h after smoking and after 48 h without smoking. The effects FPR1 on lung epithelial integrity and repair upon damage in the presence and absence of cigarette smoke were studied in CRISPR‐Cas9‐generated lung epithelial knockout cells.

Results

One significant (FDR < 0.05) inducible eQTL (rs3212855) was identified that induced a >2‐fold change in gene expression. The minor allele of rs3212855 was associated with significantly higher gene expression of FPR1, FPR2 and FPR3 upon smoking. Importantly, the minor allele of rs3212855 was also associated with lower lung function. Alveolar epithelial FPR1 knockout cells were protected against CSE‐induced reduction in repair capacity upon wounding.

Conclusion

We identified a novel smoking‐related inducible eQTL that is associated with a smoke‐induced increase in the expression of FPR1, FPR2 and FPR3, and with lowered lung function. in vitro FPR1 down‐regulation protects against smoke‐induced reduction in lung epithelial repair.

Strong correlation between air-liquid interface cultures and in vivo transcriptomics of nasal brush biopsy

Air-liquid interface (ALI) cultures are ex vivo models that are used extensively to study the epithelium of patients with chronic respiratory diseases. However, the in vitro conditions impose a milieu different from that encountered in the patient in vivo, and the degree to which this alters gene expression remains unclear. In this study we employed RNA sequencing to compare the transcriptome of fresh brushings of nasal epithelial cells with that of ALI-cultured epithelial cells from the same patients. We observed a strong correlation between cells cultured at the ALI and cells obtained from the brushed nasal epithelia: 96% of expressed genes showed similar expression profiles, although there was greater similarity between the brushed samples. We observed that while the ALI model provides an excellent representation of the in vivo airway epithelial transcriptome for mechanistic studies, several pathways are affected by the change in milieu.

Genetic profiling for disease stratification in chronic obstructive pulmonary disease and asthma

PURPOSE OF REVIEW

In asthma and chronic obstructive pulmonary disease (COPD), the movement towards genetic profiling with a push towards 'personalized medicine' has been hindered by complex environment--gene interactions and lack of tools to identify clear causal genetic traits. In this review, we will discuss the need for genetic profiling in asthma and COPD, what methods are currently used in the clinics and the recent finding using new sequencing methods.

RECENT FINDINGS

Over the past 10-15 years, genome-wide association studies analysis of common variants has provide little in the way of new genetic profiling markers for asthma and COPD. Whole exome/genome sequencing has provided a new method to identify lowly abundant alleles, which might have a much higher impact. Although, low population numbers due to high costs has hindered early studies, recent studies have reached genome wide significance.

SUMMARY

The use of genetic profiling of COPD in the clinic is current limited to the identification of Alpha-1 antitrypsin deficiency, while being absent in asthma. Advances in sequencing technology provide new avenues to identify disease causes or therapy response altering variants that in the short-term will allow for the development of screening procedures for disease to identify patients at risk of developing asthma or COPD.

Human airway epithelial extracellular vesicle miRNA signature is altered upon asthma development

BACKGROUND

miRNAs are master regulators of signaling pathways critically involved in asthma and are transferred between cells in extracellular vesicles (EV). We aimed to investigate whether the miRNA content of EV secreted by primary normal human bronchial epithelial cells (NHBE) is altered upon asthma development.

METHODS

NHBE cells were cultured at air-liquid interface and treated with interleukin (IL)-13 to induce an asthma-like phenotype. EV isolations by precipitation from basal culture medium or apical surface wash were characterized by nanoparticle tracking analysis, transmission electron microscopy, and Western blot, and EV-associated miRNAs were identified by a RT-qPCR-based profiling. Significant candidates were confirmed in EVs isolated by size-exclusion chromatography from nasal lavages of children with mild-to-moderate (n = 8) or severe asthma (n = 9), and healthy controls (n = 9).

RESULTS

NHBE cells secrete EVs to the apical and basal side. 47 miRNAs were expressed in EVs and 16 thereof were significantly altered in basal EV upon IL-13 treatment. Expression of miRNAs could be confirmed in EVs from human nasal lavages. Of note, levels of miR-92b, miR-210, and miR-34a significantly correlated with lung function parameters in children (FEV₁ FVC_%pred and FEF_25-75%pred ), thus lower sEV-miRNA levels in nasal lavages associated with airway obstruction. Subsequent ingenuity pathway analysis predicted the miRNAs to regulate Th2 polarization and dendritic cell maturation.

CONCLUSION

Our data indicate that secretion of miRNAs in EVs from the airway epithelium, in particular miR-34a, miR-92b, and miR-210, might be involved in the early development of a Th2 response in the airways and asthma.

Nasal DNA methylation profiling of asthma and rhinitis

BACKGROUND

Epigenetic signatures in the nasal epithelium, which is a primary interface with the environment and an accessible proxy for the bronchial epithelium, might provide insights into mechanisms of allergic disease.

OBJECTIVE

We aimed to identify and interpret methylation signatures in nasal epithelial brushes associated with rhinitis and asthma.

METHODS

Nasal epithelial brushes were obtained from 455 children at the 16-year follow-up of the Dutch Prevention and Incidence of Asthma and Mite Allergy birth cohort study. Epigenome-wide association studies were performed on children with asthma, rhinitis, and asthma and/or rhinitis (AsRh) by using logistic regression, and the top results were replicated in 2 independent cohorts of African American and Puerto Rican children. Significant CpG sites were related to environmental exposures (pets, active and passive smoking, and molds) during secondary school and were correlated with gene expression by RNA-sequencing (n = 244).

RESULTS

The epigenome-wide association studies identified CpG sites significantly associated with rhinitis (n = 81) and AsRh (n = 75), but not with asthma. We significantly replicated 62 of 81 CpG sites with rhinitis and 60 of 75 with AsRh, as well as 1 CpG site with asthma. Methylation of cg03565274 was negatively associated with AsRh and positively associated with exposure to pets during secondary school. DNA methylation signals associated with AsRh were mainly driven by specific IgE-positive subjects. DNA methylation related to gene transcripts that were enriched for immune pathways and expressed in immune and epithelial cells. Nasal CpG sites performed well in predicting AsRh.

CONCLUSIONS

We identified replicable DNA methylation profiles of asthma and rhinitis in nasal brushes. Exposure to pets may affect nasal epithelial methylation in relation to asthma and rhinitis.

Gene network approach reveals co-expression patterns in nasal and bronchial epithelium

Nasal gene expression profiling is a new approach to investigate the airway epithelium as a biomarker to study the activity and treatment responses of obstructive pulmonary diseases. We investigated to what extent gene expression profiling of nasal brushings is similar to that of bronchial brushings. We performed genome wide gene expression profiling on matched nasal and bronchial epithelial brushes from 77 respiratory healthy individuals. To investigate differences and similarities among regulatory modules, network analysis was performed on correlated, differentially expressed and smoking-related genes using Gaussian Graphical Models. Between nasal and bronchial brushes, 619 genes were correlated and 1692 genes were differentially expressed (false discovery rate <0.05, |Fold-change|>2). Network analysis of correlated genes showed pro-inflammatory pathways to be similar between the two locations. Focusing on smoking-related genes, cytochrome-P450 pathway related genes were found to be similar, supporting the concept of a detoxifying response to tobacco exposure throughout the airways. In contrast, cilia-related pathways were decreased in nasal compared to bronchial brushes when focusing on differentially expressed genes. Collectively, while there are substantial differences in gene expression between nasal and bronchial brushes, we also found similarities, especially in the response to the external factors such as smoking.

Cigarette smoke exposure alters phosphodiesterases in human structural lung cells

Cigarette smoke (CS), a highly complex mixture containing more than 4,000 compounds, causes aberrant cell responses leading to tissue damage around the airways and alveoli, which underlies various lung diseases. Phosphodiesterases (PDEs) are a family of enzymes that hydrolyze cyclic nucleotides. PDE inhibition induces bronchodilation, reduces the activation and recruitment of inflammatory cells, and the release of various cytokines. Currently, the selective PDE4 inhibitor roflumilast is an approved add-on treatment for patients with severe chronic obstructive pulmonary disease with chronic bronchitis and a history of frequent exacerbations. Additional selective PDE inhibitors are being tested in preclinical and clinical studies. However, the effect of chronic CS exposure on the expression of PDEs is unknown. Using mRNA isolated from nasal and bronchial brushes and lung tissues of never smokers and current smokers, we compared the gene expression of 25 PDE coding genes. Additionally, the expression and distribution of PDE3A and PDE4D in human lung tissues was examined. This study reveals that chronic CS exposure modulates the expression of various PDE members. Thus, CS exposure may change the levels of intracellular cyclic nucleotides and thereby impact the efficiency of PDE-targeted therapies.

AGER expression and alternative splicing in bronchial biopsies of smokers and never smokers

Cigarette smoking is one of the major risk factors for the development of chronic obstructive pulmonary disease (COPD). Evidence is accumulating that Receptor for Advanced Glycation-End products (RAGE)-signaling is a key pathway in the pathophysiology of COPD. To date, it is unknown how smoking affects RAGE expression. In the current study, we investigated the effect of smoking on AGER, the gene encoding RAGE, expression and on alternative splicing of AGER. To this end, we conducted RNA-Seq on bronchial biopsies for asymptomatic smokers (n = 36) and never smokers (n = 40). Total AGER gene expression was accessed using DESeq2, while alternative splicing was investigated by measuring the number of specific split reads spanning exon-exon junctions and the total split reads. One of the major isoforms of RAGE is endogenous soluble (es) RAGE, an anti-inflammatory decoy receptor, making up for approximately 10% of the total amount of soluble (s)RAGE. We found that smokers show decreased total gene expression of AGER in bronchial biopsies, while the relative abundance of the esRAGE isoform is increased. Furthermore, no difference in the serum levels of total sRAGE were observed between smokers and non-smokers. Our data indicates that smoking initiates a protective anti-inflammatory mechanism with decreased expression of the pro-inflammatory gene AGER and increased relative abundance of the anti-inflammatory isoform esRAGE.

RNA-sequencing across three matched tissues reveals shared and tissue-specific gene expression and pathway signatures of COPD

BACKGROUND

Multiple gene expression studies have been performed separately in peripheral blood, lung, and airway tissues to study COPD. We performed RNA-sequencing gene expression profiling of large-airway epithelium, alveolar macrophage and peripheral blood samples from the same subset of COPD cases and controls from the COPDGene study who underwent bronchoscopy at a single center. Using statistical and gene set enrichment approaches, we sought to improve the understanding of COPD by studying gene sets and pathways across these tissues, beyond the individual genomic determinants.

METHODS

We performed differential expression analysis using RNA-seq data obtained from 63 samples from 21 COPD cases and controls (includes four non-smokers) via the R package DESeq2. We tested associations between gene expression and variables related to lung function, smoking history, and CT scan measures of emphysema and airway disease. We examined the correlation of differential gene expression across the tissues and phenotypes, hypothesizing that this would reveal preserved and private gene expression signatures. We performed gene set enrichment analyses using curated databases and findings from prior COPD studies to provide biological and disease relevance.

RESULTS

The known smoking-related genes CYP1B1 and AHRR were among the top differential expression results for smoking status in the large-airway epithelium data. We observed a significant overlap of genes primarily across large-airway and macrophage results for smoking and airway disease phenotypes. We did not observe specific genes differentially expressed in all three tissues for any of the phenotypes. However, we did observe hemostasis and immune signaling pathways in the overlaps across all three tissues for emphysema, and amyloid and telomere-related pathways for smoking. In peripheral blood, the emphysema results were enriched for B cell related genes previously identified in lung tissue studies.

CONCLUSIONS

Our integrative analyses across COPD-relevant tissues and prior studies revealed shared and tissue-specific disease biology. These replicated and novel findings in the airway and peripheral blood have highlighted candidate genes and pathways for COPD pathogenesis.