Gender specific airway gene expression in COPD sub-phenotypes supports a role of mitochondria and of different types of leukocytes.

Chronic obstructive pulmonary disease (COPD) is a destructive inflammatory disease and the genes expressed within the lung are crucial to its pathophysiology. We have determined the RNAseq transcriptome of bronchial brush cells from 312 stringently defined ex-smoker patients. Compared to healthy controls there were for males 40 differentially expressed genes (DEGs) and 73 DEGs for females with only 26 genes shared. The gene ontology (GO) term "response to bacterium" was shared, with several different DEGs contributing in males and females. Strongly upregulated genes TCN1 and CYP1B1 were unique to males and females, respectively. For male emphysema (E)-dominant and airway disease (A)-dominant COPD (defined by computed tomography) the term "response to stress" was found for both sub-phenotypes, but this included distinct up-regulated genes for the E-sub-phenotype (neutrophil-related CSF3R, CXCL1, MNDA) and for the A-sub-phenotype (macrophage-related KLF4, F3, CD36). In E-dominant disease, a cluster of mitochondria-encoded (MT) genes forms a signature, able to identify patients with emphysema features in a confirmation cohort. The MT-CO2 gene is upregulated transcriptionally in bronchial epithelial cells with the copy number essentially unchanged. Both MT-CO2 and the neutrophil chemoattractant CXCL1 are induced by reactive oxygen in bronchial epithelial cells. Of the female DEGs unique for E- and A-dominant COPD, 88% were detected in females only. In E-dominant disease we found a pronounced expression of mast cell-associated DEGs TPSB2, TPSAB1 and CPA3. The differential genes discovered in this study point towards involvement of different types of leukocytes in the E- and A-dominant COPD sub-phenotypes in males and females.

Lung microbiome composition and bronchial epithelial gene expression in patients with COPD versus healthy individuals: a bacterial 16S rRNA gene sequencing and host transcriptomic analysis.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is associated with airway inflammation and bacterial dysbiosis. The relationship between the airway microbiome and bronchial gene expression in COPD is poorly understood. We aimed to identify differences in the airway microbiome from bronchial brushings in patients with COPD and healthy individuals and to investigate whether any distinguishing bacteria are related to bronchial gene expression. METHODS: For this 16S rRNA gene sequencing and host transcriptomic analysis, individuals aged 45-75 years with mild-to-moderate COPD either receiving or not receiving inhaled corticosteroids and healthy individuals in the same age group were recruited as part of the Emphysema versus Airways Disease (EvA) consortium from nine centres in the UK, Germany, Italy, Poland, and Hungary. Individuals underwent clinical characterisation, spirometry, CT scans, and bronchoscopy. From bronchoscopic bronchial brush samples, we obtained the microbial profiles using 16S rRNA gene sequencing and gene expression using the RNA-Seq technique. We analysed bacterial genera relative abundance and the associations between genus abundance and clinical characteristics or between genus abundance and host lung transcriptional signals in patients with COPD versus healthy individuals, and in patients with COPD with versus without inhaled corticosteroids treatment. FINDINGS: Between February, 2009, and March, 2012, we obtained brush samples from 574 individuals. We used 546 of 574 samples for analysis, including 207 from healthy individuals and 339 from patients with COPD (192 with inhaled corticosteroids and 147 without). The bacterial genera that most strongly distinguished patients with COPD from healthy individuals were Prevotella (median relative abundance 33·5%, IQR 14·5-49·4, in patients with COPD vs 47·7%, 31·1-60·7, in healthy individuals; p<0·0001), Streptococcus (8·6%, 3·8-15·8, vs 5·3%, 3·0-10·1; p<0·0001), and Moraxella (0·05%, 0·02-0·14, vs 0·02%, 0-0·07; p<0·0001). Prevotella abundance was inversely related to COPD severity in terms of symptoms and positively related to lung function and exercise capacity. 446 samples had assessable RNA-seq data, 257 from patients with COPD (136 with inhaled corticosteroids and 121 without) and 189 from healthy individuals. No significant associations were observed between lung transcriptional signals from bronchial brushings and abundance of bacterial genera in patients with COPD without inhaled corticosteroids treatment and in healthy individuals. In patients with COPD treated with inhaled corticosteroids, Prevotella abundance was positively associated with expression of epithelial genes involved in tight junction promotion and Moraxella abundance was associated with expression of the IL-17 and TNF inflammatory pathways. INTERPRETATION: With increasing severity of COPD, the airway microbiome is associated with decreased abundance of Prevotella and increased abundance of Moraxella in concert with downregulation of genes promoting epithelial defence and upregulation of pro-inflammatory genes associated with inhaled corticosteroids use. Our work provides further insight in understanding the relationship between microbiome alteration and host inflammatory response, which might lead to novel therapeutic strategies for COPD. FUNDING: EU Seventh Framework Programme, National Institute for Health Research.

Subsets of CD1c+ DCs: Dendritic Cell Versus Monocyte Lineage.

Currently three bona fide dendritic cell (DC) types are distinguished in human blood. Herein we focus on type 2 DCs (DC2s) and compare the three defining markers CD1c, CD172, and CD301. When using CD1c to define DC2s, a CD14+ and a CD14- subset can be detected. The CD14+ subset shares features with monocytes, and this includes substantially higher expression levels for CD64, CD115, CD163, and S100A8/9. We review the current knowledge of these CD1c+CD14+ cells as compared to the CD1c+CD14- cells with respect to phenotype, function, transcriptomics, and ontogeny. Here, we discuss informative mutations, which suggest that two populations have different developmental requirements. In addition, we cover subsets of CD11c+CD8- DC2s in the mouse, where CLEC12A+ESAMlow cells, as compared to the CLEC12A-ESAMhigh subset, also express higher levels of monocyte-associated markers CD14, CD3, and CD115. Finally, we summarize, for both man and mouse, the data on lower antigen presentation and higher cytokine production in the monocyte-marker expressing DC2 subset, which demonstrate that the DC2 subsets are also functionally distinct.

Epigenetics in non-classical monocytes support their pro-inflammatory gene expression.

Non-classical human monocytes are characterized by high-level expression of cytokines like TNF, but the mechanisms involved are elusive. We have identified miRNAs and CpG-methylation sites that are unique to non-classical monocytes, defined via CD14 and CD16 expression levels. For down-regulated miRNAs that are linked to up-regulated mRNAs the dominant gene ontology term was intracellular signal transduction. This included down-regulated miRNA-20a-5p and miRNA-106b-5p, which both are linked to increased mRNA for the TRIM8 signaling molecule. Methylation analysis revealed 16 hypo-methylated CpG sites upstream of 14 differentially increased mRNAs including 2 sites upstream of TRIM8. Consistent with a positive role in signal transduction, high TRIM8 levels went along with high basal TNF mRNA levels in non-classical monocytes. Since cytokine expression levels in monocytes strongly increase after stimulation with toll-like-receptor ligands, we have analyzed non-classical monocytes (defined via slan expression) after stimulation with lipopolysaccharide (LPS). LPS-stimulated cells continued to have low miRNA-20a and miRNA-106b and high TRIM8 mRNA levels and they showed a 10-fold increase in TNF mRNA. These data suggest that decreased miRNAs and CpG hypo-methylation is linked to enhanced expression of TRIM8 and that this can contribute to the increased TNF levels in non-classical human monocytes.

Blood eosinophil count and airway epithelial transcriptome relationships in COPD versus asthma.

BACKGROUND: Whether the clinical or pathophysiologic significance of the "treatable trait" high blood eosinophil count in COPD is the same as for asthma remains controversial. We sought to determine the relationship between the blood eosinophil count, clinical characteristics and gene expression from bronchial brushings in COPD and asthma. METHODS: Subjects were recruited into a COPD (emphysema versus airway disease [EvA]) or asthma cohort (Unbiased BIOmarkers in PREDiction of respiratory disease outcomes, U-BIOPRED). We determined gene expression using RNAseq in EvA (n = 283) and Affymetrix microarrays in U-BIOPRED (n = 85). We ran linear regression analysis of the bronchial brushings transcriptional signal versus blood eosinophil counts as well as differential expression using a blood eosinophil > 200 cells/µL as a cut-off. The false discovery rate was controlled at 1% (with continuous values) and 5% (with dichotomized values). RESULTS: There were no differences in age, gender, lung function, exercise capacity and quantitative computed tomography between eosinophilic versus noneosinophilic COPD cases. Total serum IgE was increased in eosinophilic asthma and COPD. In EvA, there were 12 genes with a statistically significant positive association with the linear blood eosinophil count, whereas in U-BIOPRED, 1197 genes showed significant associations (266 positive and 931 negative). The transcriptome showed little overlap between genes and pathways associated with blood eosinophil counts in asthma versus COPD. Only CST1 was common to eosinophilic asthma and COPD and was replicated in independent cohorts. CONCLUSION: Despite shared "treatable traits" between asthma and COPD, the molecular mechanisms underlying these clinical entities are predominately different.

Resistome analyses of sputum from COPD and healthy subjects reveals bacterial load-related prevalence of target genes.

BACKGROUND: Antibiotic resistance is a major global threat. We hypothesised that the chronic obstructive pulmonary disease (COPD) airway is a reservoir of antimicrobial resistance genes (ARGs) that associate with microbiome-specific COPD subgroups. OBJECTIVE: To determine the resistance gene profiles in respiratory samples from COPD patients and healthy volunteers. METHODS: Quantitative PCR targeting 279 specific ARGs was used to profile the resistomes in sputum from subjects with COPD at stable, exacerbation and recovery visits (n=55; COPD-BEAT study), healthy controls with (n=7) or without (n=22) exposure to antibiotics in the preceding 12 months (EXCEED study) and in bronchial brush samples from COPD (n=8) and healthy controls (n=7) (EvA study). RESULTS: ARG mean (SEM) prevalence was greater in stable COPD samples (35.2 (1.6)) than in healthy controls (27.6 (1.7); p=0.004) and correlated with total bacterial abundance (r2=0.23; p<0.001). Prevalence of ARG positive signals in individuals was not related to COPD symptoms, lung function or their changes at exacerbation. In the COPD subgroups designated High γProteobacteria and High Firmicutes, ARG prevalence was not different at stable state but significantly declined from stable through exacerbation to recovery in the former (p=0.011) without changes in total bacterial abundance. The ARG patterns were similar in COPD versus health, COPD microbiome-subgroups and between sputum and bronchoscopic samples independent of antibiotic exposure in the last 12 months. CONCLUSIONS: ARGs are highly prevalent in sputum, broadly in proportion to bacterial abundance in both healthy and COPD subjects. Thus, COPD appears to be an ARG reservoir due to high levels of bacterial colonisation.

High degree of polyclonality hinders somatic mutation calling in lung brush samples of COPD cases and controls.

Chronic obstructive pulmonary disease (COPD) is induced by cigarette smoking and characterized by inflammation of airway tissue. Since smokers with COPD have a higher risk of developing lung cancer than those without, we hypothesized that they carry more mutations in affected tissue. We called somatic mutations in airway brush samples from medium-coverage whole genome sequencing data from healthy never and ex-smokers (n = 8), as well as from ex-smokers with variable degrees of COPD (n = 4). Owing to the limited concordance of resulting calls between the applied tools we built a consensus, a strategy that was validated with high accuracy for cancer data. However, consensus calls showed little promise of representing true positives due to low mappability of corresponding sequence reads and high overlap with positions harbouring known genetic polymorphisms. A targeted re-sequencing approach suggested that only few mutations would survive stringent verification testing and that our data did not allow the inference of any difference in the mutational load of bronchial brush samples between former smoking COPD cases and controls. High polyclonality in airway brush samples renders medium-depth sequencing insufficient to provide the resolution to detect somatic mutations. Deep sequencing data of airway biopsies are needed to tackle the question.

6-Sulfo LacNAc (Slan) as a Marker for Non-classical Monocytes.

Monocytes are subdivided into three subsets, which have different phenotypic and functional characteristics and different roles in inflammation and malignancy. When in man CD14 and CD16 monoclonal antibodies are used to define these subsets, then the distinction of non-classical CD14low and intermediate CD14high monocytes requires setting a gate in what is a gradually changing level of CD14 expression. In the search for an additional marker to better dissect the two subsets we have explored the marker 6-sulfo LacNAc (slan). Slan is a carbohydrate residue originally described to be expressed on the cell surface of a type of dendritic cell in human blood. We elaborate herein that the features of slan+ cells are congruent with the features of CD16+ non-classical monocytes and that slan is a candidate marker for definition of non-classical monocytes. The use of this marker may help in studying the role of non-classical monocytes in health and in diagnosis and monitoring of disease.

Influence of lung CT changes in chronic obstructive pulmonary disease (COPD) on the human lung microbiome.

BACKGROUND: Changes in microbial community composition in the lung of patients suffering from moderate to severe COPD have been well documented. However, knowledge about specific microbiome structures in the human lung associated with CT defined abnormalities is limited. METHODS: Bacterial community composition derived from brush samples from lungs of 16 patients suffering from different CT defined subtypes of COPD and 9 healthy subjects was analyzed using a cultivation independent barcoding approach applying 454-pyrosequencing of 16S rRNA gene fragment amplicons. RESULTS: We could show that bacterial community composition in patients with changes in CT (either airway or emphysema type changes, designated as severe subtypes) was different from community composition in lungs of patients without visible changes in CT as well as from healthy subjects (designated as mild COPD subtype and control group) (PC1, Padj = 0.002). Higher abundance of Prevotella in samples from patients with mild COPD subtype and from controls and of Streptococcus in the severe subtype cases mainly contributed to the separation of bacterial communities of subjects. No significant effects of treatment with inhaled glucocorticoids on bacterial community composition were detected within COPD cases with and without abnormalities in CT in PCoA. Co-occurrence analysis suggests the presence of networks of co-occurring bacteria. Four communities of positively correlated bacteria were revealed. The microbial communities can clearly be distinguished by their associations with the CT defined disease phenotype. CONCLUSION: Our findings indicate that CT detectable structural changes in the lung of COPD patients, which we termed severe subtypes, are associated with alterations in bacterial communities, which may induce further changes in the interaction between microbes and host cells. This might result in a changed interplay with the host immune system.

Reprint of: MicroRNA profiling of human intermediate monocytes.

Among the three human monocyte subsets, intermediate CD14++CD16+ monocytes have been characterized as particularly proinflammatory cells in experimental studies and as potential biomarkers of cardiovascular risk in clinical cohorts. To further substantiate the distinct role of intermediate monocytes within human monocyte heterogeneity, we assessed subset-specific expression of miRNAs as central epigenetic regulators of gene expression. We hypothesized that intermediate monocytes have a distinct miRNA profile compared to classical and non-classical monocytes. By using small RNA-seq we analyzed 662 miRNAs in the three monocyte subsets. We identified 38 miRNAs that are differentially expressed in intermediate monocytes compared to both classical and non-classical monocytes with a p value of <10-10, of which two miRNAs - miR-6087 (upregulated) and miR-150-5p (downregulated) - differed in their expression more than ten-fold. Pathway analysis of the 38 differentially expressed miRNAs linked intermediate monocytes to distinct biological processes such as gene regulation, cell differentiation, toll-like receptor signaling as well as antigen processing and presentation. Moreover, differentially expressed miRNAs were connected to those genes that we previously identified as markers of intermediate monocytes. In aggregation, we provide first genome-wide miRNA data in the context of monocyte heterogeneity, which substantiate the concept of monocyte trichotomy in human immunity. The identification of miRNAs that are specific for intermediate monocytes may allow to develop strategies, which particularly target this cell population while sparing the other two subsets.

MicroRNA profiling of human intermediate monocytes.

Among the three human monocyte subsets, intermediate CD14++CD16+ monocytes have been characterized as particularly proinflammatory cells in experimental studies and as potential biomarkers of cardiovascular risk in clinical cohorts. To further substantiate the distinct role of intermediate monocytes within human monocyte heterogeneity, we assessed subset-specific expression of miRNAs as central epigenetic regulators of gene expression. We hypothesized that intermediate monocytes have a distinct miRNA profile compared to classical and non-classical monocytes. By using small RNA-seq we analyzed 662 miRNAs in the three monocyte subsets. We identified 38 miRNAs that are differentially expressed in intermediate monocytes compared to both classical and non-classical monocytes with a p value of <10-10, of which two miRNAs - miR-6087 (upregulated) and miR-150-5p (downregulated) - differed in their expression more than ten-fold. Pathway analysis of the 38 differentially expressed miRNAs linked intermediate monocytes to distinct biological processes such as gene regulation, cell differentiation, toll-like receptor signaling as well as antigen processing and presentation. Moreover, differentially expressed miRNAs were connected to those genes that we previously identified as markers of intermediate monocytes. In aggregation, we provide first genome-wide miRNA data in the context of monocyte heterogeneity, which substantiate the concept of monocyte trichotomy in human immunity. The identification of miRNAs that are specific for intermediate monocytes may allow to develop strategies, which particularly target this cell population while sparing the other two subsets.

Emphysema- and airway-dominant COPD phenotypes defined by standardised quantitative computed tomography.

EvA (Emphysema versus Airway disease) is a multicentre project to study mechanisms and identify biomarkers of emphysema and airway disease in chronic obstructive pulmonary disease (COPD). The objective of this study was to delineate objectively imaging-based emphysema-dominant and airway disease-dominant phenotypes using quantitative computed tomography (QCT) indices, standardised with a novel phantom-based approach.441 subjects with COPD (Global Initiative for Chronic Obstructive Lung Disease (GOLD) stages 1-3) were assessed in terms of clinical and physiological measurements, laboratory testing and standardised QCT indices of emphysema and airway wall geometry.QCT indices were influenced by scanner non-conformity, but standardisation significantly reduced variability (p<0.001) and led to more robust phenotypes. Four imaging-derived phenotypes were identified, reflecting "emphysema-dominant", "airway disease-dominant", "mixed" disease and "mild" disease. The emphysema-dominant group had significantly higher lung volumes, lower gas transfer coefficient, lower oxygen (PO2 ) and carbon dioxide (PCO2 ) tensions, higher haemoglobin and higher blood leukocyte numbers than the airway disease-dominant group.The utility of QCT for phenotyping in the setting of an international multicentre study is improved by standardisation. QCT indices of emphysema and airway disease can delineate within a population of patients with COPD, phenotypic groups that have typical clinical features known to be associated with emphysema-dominant and airway-dominant disease.

slan-defined subsets of CD16-positive monocytes: impact of granulomatous inflammation and M-CSF receptor mutation.

Human monocytes are subdivided into classical, intermediate, and nonclassical subsets, but there is no unequivocal strategy to dissect the latter 2 cell types. We show herein that the cell surface marker 6-sulfo LacNAc (slan) can define slan-positive CD14(+)CD16(++) nonclassical monocytes and slan-negative CD14(++)CD16(+) intermediate monocytes. Gene expression profiling confirms that slan-negative intermediate monocytes show highest expression levels of major histocompatibility complex class II genes, whereas a differential ubiquitin signature is a novel feature of the slan approach. In unsupervised hierarchical clustering, the slan-positive nonclassical monocytes cluster with monocytes and are clearly distinct from CD1c(+) dendritic cells. In clinical studies, we show a selective increase of the slan-negative intermediate monocytes to >100 cells per microliter in patients with sarcoidosis and a fivefold depletion of the slan-positive monocytes in patients with hereditary diffuse leukoencephalopathy with axonal spheroids (HDLS), which is caused by macrophage colony-stimulating factor (M-CSF) receptor mutations. These data demonstrate that the slan-based definition of CD16-positive monocyte subsets is informative in molecular studies and in clinical settings.

Blood Monocytes and Their Subsets: Established Features and Open Questions.

In contrast to the past reliance on morphology, the identification and enumeration of blood monocytes are nowadays done with monoclonal antibodies and flow cytometry and this allows for subdivision into classical, intermediate, and non-classical monocytes. Using specific cell surface markers, dendritic cells in blood can be segregated from these monocytes. While in the past, changes in monocyte numbers as determined in standard hematology counters have not had any relevant clinical impact, the subset analysis now has uncovered informative changes that may be used in management of disease.

Comparison of two different strategies for human monocyte subsets gating within the large-scale prospective CARE FOR HOMe Study.

Monocytes are heterogeneous cells consisting of (at least) three subsets: classical, intermediate, and nonclassical monocytes. Correct enumeration of cell counts necessitates well-defined gating strategies, which are essentially based upon CD14 and CD16 expression. For the delineation of intermediate from nonclassical monocytes, a "rectangular gating (RG) strategy" and a "trapezoid gating (TG) strategy" have been proposed. We compared the two gating strategies in a well-defined clinical cohort of patients with chronic kidney disease (CKD). Within the ongoing CARE FOR HOMe study, monocyte subsets were reanalyzed in 416 CKD patients, who were followed 3.6 ± 1.6 years for the occurrence of a cardiovascular event. Gating was performed by either RG or TG. We analyzed the expression of surface markers, and compared the predictive role of cell counts of monocyte subsets, as defined by RG and TG, respectively. With both gating strategies, higher intermediate monocyte counts predicted the cardiovascular endpoint in Kaplan-Meier analyses (P < 0.001 with RG; P < 0.001 with TG). After correction for confounders, intermediate monocyte counts remained independent predictors in Cox-Regression analyses (HR = 1.013 [95% CI: 1.006-1.020; P < 0.001] with RG; HR = 1.015 [95% CI: 1.006-1.024; P = 0.001] with TG). NRI was 3.9% when reclassifying patients from quartiles of intermediate monocyte counts with RG strategy toward quartiles of intermediate monocytes counts with TG strategy. In expression analysis, those monocytes which are defined as intermediate monocytes by the RG strategy and as nonclassical monocytes by the TG strategy share characteristics of both subsets. In conclusion, intermediate monocytes were independent predictors of cardiovascular outcome irrespective of the applied gating strategy. Future studies should aim to identify markers that allow for an unequivocal definition of intermediate monocytes, which may further improve their power to predict cardiovascular events.

Mechanisms of Hypoxic Up-Regulation of Versican Gene Expression in Macrophages.

Hypoxia is a hallmark of many pathological tissues. Macrophages accumulate in hypoxic sites and up-regulate a range of hypoxia-inducible genes. The matrix proteoglycan versican has been identified as one such gene, but the mechanisms responsible for hypoxic induction are not fully characterised. Here we investigate the up-regulation of versican by hypoxia in primary human monocyte-derived macrophages (HMDM), and, intriguingly, show that versican mRNA is up-regulated much more highly (>600 fold) by long term hypoxia (5 days) than by 1 day of hypoxia (48 fold). We report that versican mRNA decay rates are not affected by hypoxia, demonstrating that hypoxic induction of versican mRNA is mediated by increased transcription. Deletion analysis of the promoter identified two regions required for high level promoter activity of luciferase reporter constructs in human macrophages. The hypoxia-inducible transcription factor HIF-1 has previously been implicated as a key potential regulator of versican expression in hypoxia, however our data suggest that HIF-1 up-regulation is unlikely to be principally responsible for the high levels of induction observed in HMDM. Treatment of HMDM with two distinct specific inhibitors of Phosphoinositide 3-kinase (PI3K), LY290042 and wortmannin, significantly reduced induction of versican mRNA by hypoxia and provides evidence of a role for PI3K in hypoxic up-regulation of versican expression.