Methylome and transcriptome signature of bronchoalveolar cells from multiple sclerosis patients in relation to smoking

BACKGROUND

Despite compelling evidence that cigarette smoking impacts the risk of developing multiple sclerosis (MS), little is known about smoking-associated changes in the primary exposed lung cells of patients.

OBJECTIVES

We aimed to examine molecular changes occurring in bronchoalveolar lavage (BAL) cells from MS patients in relation to smoking and in comparison to healthy controls (HCs).

METHODS

We profiled DNA methylation in BAL cells from female MS (n = 17) and HC (n = 22) individuals, using Illumina Infinium EPIC and performed RNA-sequencing in non-smokers.

RESULTS

The most prominent changes were found in relation to smoking, with 1376 CpG sites (adjusted P < 0.05) differing between MS smokers and non-smokers. Approximately 30% of the affected genes overlapped with smoking-associated changes in HC, leading to a strong common smoking signature in both MS and HC after gene ontology analysis. Smoking in MS patients resulted in additional discrete changes related to neuronal processes. Methylome and transcriptome analyses in non-smokers suggest that BAL cells from MS patients display very subtle (not reaching adjusted P < 0.05) but concordant changes in genes connected to reduced transcriptional/translational processes and enhanced cellular motility.

CONCLUSIONS

Our study provides insights into the impact of smoking on lung inflammation and immunopathogenesis of MS.

Hospital diagnosed pneumonia before age 20 years and multiple sclerosis risk

Introduction

Respiratory inflammation has been proposed as a risk factor for MS. This study aims to determine if hospital-diagnosed pneumonia in adolescence (before age 20 years) is associated with subsequent multiple sclerosis (MS).

Methods

This case-control study included incident MS cases after age 20 years identified using the Swedish national registers. Cases were matched with 10 general population controls by age, sex and region. Pneumonia diagnoses were identified between 0–5, 6–10, 11–15 and 16–20 years of age. Conditional logistic regression models adjusted for infectious mononucleosis (IM) and education calculated ORs with 95% CIs. Urinary tract infections (UTIs), a common complication of MS, before age 20 years were included as a control diagnosis for reverse causation.

Results

There were 6109 cases and 49 479 controls included. Pneumonia diagnosed between age 11–15 years was associated with subsequent MS (adj OR 2.00, 95% CI 1.22 to 3.27). Although not statistically significant, sensitivity analyses showed similar magnitude associations of pneumonia between age 11–15 years and MS. No statistically significant associations with MS for pneumonia at other age groups were observed. Adjustment for IM had no notable effect on associations, but was statistically significantly associated with MS. UTIs were not associated with MS.

Conclusion

Pneumonia at 11–15 years of age was associated with MS, suggesting a possible role for inflammation of the respiratory system in the aetiology of MS during a period of susceptibility in adolescence. Further research on respiratory infections prior to MS onset should be conducted to replicate this finding and determine explanatory causal mechanisms.

Tobacco smoking induces changes in true DNA methylation, hydroxymethylation and gene expression in bronchoalveolar lavage cells

Background

While smoking is known to associate with development of multiple diseases, the underlying mechanisms are still poorly understood. Tobacco smoking can modify the chemical integrity of DNA leading to changes in transcriptional activity, partly through an altered epigenetic state. We aimed to investigate the impact of smoking on lung cells collected from bronchoalveolar lavage (BAL).

Methods

We profiled changes in DNA methylation (5mC) and its oxidised form hydroxymethylation (5hmC) using conventional bisulphite (BS) treatment and oxidative bisulphite treatment with Illumina Infinium MethylationEPIC BeadChip, and examined gene expression by RNA-seq in healthy smokers.

Findings

We identified 1667 total 5mC + 5hmC, 1756 5mC and 67 5hmC differentially methylated positions (DMPs) between smokers and non-smokers (FDR-adjusted P <.05, absolute Δβ >0.15). Both 5mC DMPs and to a lesser extent 5mC + 5hmC were predominantly hypomethylated. In contrast, almost all 5hmC DMPs were hypermethylated, supporting the hypothesis that smoking-associated oxidative stress can lead to DNA demethylation, via the established sequential oxidation of which 5hmC is the first step. While we confirmed differential methylation of previously reported smoking-associated 5mC + 5hmC CpGs using former generations of BeadChips in alveolar macrophages, the large majority of identified DMPs, 5mC + 5hmC (1639/1667), 5mC (1738/1756), and 5hmC (67/67), have not been previously reported. Most of these novel smoking-associating sites are specific to the EPIC BeadChip and, interestingly, many of them are associated to FANTOM5 enhancers. Transcriptional changes affecting 633 transcripts were consistent with DNA methylation profiles and converged to alteration of genes involved in migration, signalling and inflammatory response of immune cells.

Interpretation

Collectively, these findings suggest that tobacco smoke exposure epigenetically modifies BAL cells, possibly involving a continuous active demethylation and subsequent increased activity of inflammatory processes in the lungs.

Fund

The study was supported by the Swedish Research Council, the Swedish Heart-Lung Foundation, the Stockholm County Council (ALF), the King Gustav's and Queen Victoria's Freemasons' Foundation, Knut and Alice Wallenberg Foundation, Neuro Sweden, and the Swedish MS foundation.

Expression of MATE1, P-gp, OCTN1 and OCTN2, in epithelial and immune cells in the lung of COPD and healthy individuals

Background

Several inhaled drugs are dependent on organic cation transporters to cross cell membranes. To further evaluate their potential to impact on inhaled drug disposition, the localization of MATE1, P-gp, OCTN1 and OCTN2 were investigated in human lung.

Methods

Transporter proteins were analysed by immunohistochemistry in lung tissue from healthy subjects and COPD patients. Transporter mRNA was analysed by qPCR in lung tissue and in bronchoalveolar lavage (BAL) cells from smokers and non-smokers.

Results

We demonstrate for the first time MATE1 protein expression in the lung with localization to the apical side of bronchial and bronchiolar epithelial cells. Interestingly, MATE1 was strongly expressed in alveolar macrophages as demonstrated both in lung tissue and in BAL cells, and in inflammatory cells including CD3 positive T cells. P-gp, OCTN1 and OCTN2 were also expressed in the alveolar epithelial cells and in inflammatory cells including alveolar macrophages. In BAL cells from smokers, MATE1 and P-gp mRNA expression was significantly lower compared to cells from non-smokers whereas no difference was observed between COPD patients and healthy subjects. THP-1 cells were evaluated as a model for alveolar macrophages but did not reflect the transporter expression observed in BAL cells.

Conclusions

We conclude that MATE1, P-gp, OCTN1 and OCTN2 are expressed in pulmonary lung epithelium, in alveolar macrophages and in other inflammatory cells. This is important to consider in the development of drugs treating pulmonary disease as the transporters may impact drug disposition in the lung and consequently affect pharmacological efficacy and toxicity.

Electronic supplementary material

The online version of this article (10.1186/s12931-018-0760-9) contains supplementary material, which is available to authorized users.

The lung microbiota in early rheumatoid arthritis and autoimmunity

BACKGROUND

Airway abnormalities and lung tissue citrullination are found in both rheumatoid arthritis (RA) patients and individuals at-risk for disease development. This suggests the possibility that the lung could be a site of autoimmunity generation in RA, perhaps in response to microbiota changes. We therefore sought to test whether the RA lung microbiome contains distinct taxonomic features associated with local and/or systemic autoimmunity.

METHODS

16S rRNA gene high-throughput sequencing was utilized to compare the bacterial community composition of bronchoalveolar lavage fluid (BAL) in patients with early, disease-modifying anti-rheumatic drugs (DMARD)-naïve RA, patients with lung sarcoidosis, and healthy control subjects. Samples were further assessed for the presence and levels of anti-citrullinated peptide antibodies (including fine specificities) in both BAL and serum.

RESULTS

The BAL microbiota of RA patients was significantly less diverse and abundant when compared to healthy controls, but similar to sarcoidosis patients. This distal airway dysbiosis was attributed to the reduced presence of several genus (i.e., Actynomyces and Burkhordelia) as well as reported periodontopathic taxa, including Treponema, Prevotella, and Porphyromonas. While multiple clades correlated with local and systemic levels of autoantibodies, the genus Pseudonocardia and various related OTUs were the only taxa overrepresented in RA BAL and correlated with higher disease activity and erosions.

CONCLUSIONS

Distal airway dysbiosis is present in untreated early RA and similar to that detected in sarcoidosis lung inflammation. This community perturbation, which correlates with local and systemic autoimmune/inflammatory changes, may potentially drive initiation of RA in a proportion of cases.

Distinctive Regulatory T Cells and Altered Cytokine Profile Locally in the Airways of Young Smokers with Normal Lung Function

Smoking influences the immune system in different ways and, hypothetically, effects on pulmonary effector and regulatory T cells emerge as potentially detrimental. Therefore, we characterized the frequencies and characteristics of CD4+ and CD8+ T cell subsets in the blood and lungs of young tobacco smokers. Bronchoalveolar lavage (BAL) and peripheral blood were obtained from healthy moderate smokers (n = 18; 2–24 pack-years) and never-smokers (n = 15), all with normal lung function. Cells were stimulated ex vivo and key intracellular cytokines (IFNγ, IL-17, IL-10 and TNFα) and transcription factors (Foxp3, T-bet and Helios) were analyzed using flow cytometry. Our results indicate that smoking is associated with a decline in lung IL-17+ CD4+ T cells, increased IFNγ+ CD8+ T cells and these alterations relate to the history of daily cigarette consumption. There is an increased fraction of Foxp3+ regulatory T cells being Helios- in the lungs of smokers. Cytokine production is mainly confined to the Helios- T cells, both in regulatory and effector subsets. Moreover, we detected a decline of Helios+Foxp3- postulated regulatory CD8+ T cells in smokers. These alterations in the immune system are likely to increase risk for infection and may have implications for autoimmune processes initiated in the lungs among tobacco smokers.

T-cell activation and HLA-regulated response to smoking in the deep airways of patients with multiple sclerosis

Cigarette smoking is a risk factor for multiple sclerosis (MS), and the risk is further multiplied for HLA-DRB1*15(+) smokers. To define the smoke-induced immune responses in the lung we performed bronchoscopy with bronchoalveolar lavage (BAL) on smokers and non-smokers, both MS-patients and healthy volunteers. In the BAL, non-smokers with MS showed an increased preformed CD40L expression in CD4(+) T-cells while smokers displayed an increase in proliferating (Ki-67(+)) T-cells. In addition, our results confirm that smoking induces an increase of alveolar macrophages in BAL, and further defined a significant attenuation of this response in carriers of the HLA-DRB1*15 allele, in both MS patients and healthy controls. This first systematic investigation of the immune response in the lungs of smokers and non-smokers diagnosed with MS, thus suggests an MS-associated lung T-cell phenotype, involvement of a specific T-cell response to smoke, and a genetic regulation of the macrophage response.

Parent-of-origin effects implicate epigenetic regulation of experimental autoimmune encephalomyelitis and identify imprinted Dlk1 as a novel risk gene

Parent-of-origin effects comprise a range of genetic and epigenetic mechanisms of inheritance. Recently, detection of such effects implicated epigenetic mechanisms in the etiology of multiple sclerosis (MS), a chronic inflammatory disease of the central nervous system. We here sought to dissect the magnitude and the type of parent-of-origin effects in the pathogenesis of experimental neuroinflammation under controlled environmental conditions. We investigated inheritance of an MS-like disease in rat, experimental autoimmune encephalomyelitis (EAE), using a backcross strategy designed to identify the parental origin of disease-predisposing alleles. A striking 37–54% of all detected disease-predisposing loci depended on parental transmission. Additionally, the Y chromosome from the susceptible strain contributed to disease susceptibility. Accounting for parent-of-origin enabled more powerful and precise identification of novel risk factors and increased the disease variance explained by the identified factors by 2-4-fold. The majority of loci displayed an imprinting–like pattern whereby a gene expressed only from the maternal or paternal copy exerts an effect. In particular, a locus on chromosome 6 comprises a well-known cluster of imprinted genes including the paternally expressed Dlk1, an atypical Notch ligand. Disease-predisposing alleles at the locus conferred lower Dlk1 expression in rats and, together with data from transgenic overexpressing Dlk1 mice, demonstrate that reduced Dlk1 drives more severe disease and modulates adaptive immune reactions in EAE. Our findings suggest a significant epigenetic contribution to the etiology of EAE. Incorporating these effects enables more powerful and precise identification of novel risk factors with diagnostic and prognostic implications for complex disease.

Even with recent progress in determining the genetic basis of complex diseases, the issue of ‘missing heritability’ remains and its potential sources are frequently speculated about but rarely explained. Parent-of-origin effects might contribute to the ‘missing heritability’ and involve genetic and epigenetic mechanisms of inheritance. Our study is the first that establishes (i) the magnitude and (ii) the type of parent-of-origin effects in the pathogenesis of a multiple sclerosis-like disease, experimental autoimmune encephalomyelitis (EAE) in rat, using a strategy designed to identify genes that confer risk only when inherited from either mother or father. A striking 37-54% of all risk loci depended on parental origin. Accounting for parent-of-origin enabled more powerful and precise identification of novel risk factors for EAE, such as the imprinted Dlk1gene. Disease-predisposing alleles conferred lower Dlk1 expression in rats and transgenic Dlk1 mice demonstrated that lower Dlk1 drives more severe EAE and modulates adaptive immune responses. Because parental-origin effects are epigenetically regulated, our data implicate a contributory role for epigenetic mechanisms in complex diseases. Considering parent-of-origin effects in complex disease has enabled more powerful and precise identification of novel risk factors.

Combining genetic mapping with genome-wide expression in experimental autoimmune encephalomyelitis highlights a gene network enriched for T cell functions and candidate genes regulating autoimmunity

The experimental autoimmune encephalomyelitis (EAE) is an autoimmune disease of the central nervous system commonly used to study multiple sclerosis (MS). We combined clinical EAE phenotypes with genome-wide expression profiling in spleens from 150 backcross rats between susceptible DA and resistant PVG rat strains during the chronic EAE phase. This enabled correlation of transcripts with genotypes, other transcripts and clinical EAE phenotypes and implicated potential genetic causes and pathways in EAE. We detected 2285 expression quantitative trait loci (eQTLs). Sixty out of 599 cis-eQTLs overlapped well-known EAE QTLs and constitute positional candidate genes, including Ifit1 (Eae7), Atg7 (Eae20-22), Klrc3 (eEae22) and Mfsd4 (Eae17). A trans-eQTL that overlaps Eae23a regulated a large number of small RNAs and implicates a master regulator of transcription. We defined several disease-correlated networks enriched for pathways involved in cell-mediated immunity. They include C-type lectins, G protein coupled receptors, mitogen-activated protein kinases, transmembrane proteins, suppressors of transcription (Jundp2 and Nr1d1) and STAT transcription factors (Stat4) involved in interferon signaling. The most significant network was enriched for T cell functions, similar to genetic findings in MS, and revealed both established and novel gene interactions. Transcripts in the network have been associated with T cell proliferation and differentiation, the TCR signaling and regulation of regulatory T cells. A number of network genes and their family members have been associated with MS and/or other autoimmune diseases. Combining disease and genome-wide expression phenotypes provides a link between disease risk genes and distinct molecular pathways that are dysregulated during chronic autoimmune inflammation.

Expression of Ccl11 associates with immune response modulation and protection against neuroinflammation in rats

Multiple sclerosis (MS) is a polygenic disease characterized by inflammation and demyelination in the central nervous system (CNS), which can be modeled in experimental autoimmune encephalomyelitis (EAE). The Eae18b locus on rat chromosome 10 has previously been linked to regulation of beta-chemokine expression and severity of EAE. Moreover, the homologous chemokine cluster in humans showed evidence of association with susceptibility to MS. We here established a congenic rat strain with Eae18b locus containing a chemokine cluster (Ccl2, Ccl7, Ccl11, Ccl12 and Ccl1) from the EAE- resistant PVG rat strain on the susceptible DA background and utilized myelin oligodendrocyte glycoprotein (MOG)-induced EAE to characterize the mechanisms underlying the genetic regulation. Congenic rats developed a milder disease compared to the susceptible DA strain, and this was reflected in decreased demyelination and in reduced recruitment of inflammatory cells to the brain. The congenic strain also showed significantly increased Ccl11 mRNA expression in draining lymph nodes and spinal cord after EAE induction. In the lymph nodes, macrophages were the main producers of CCL11, whereas macrophages and lymphocytes expressed the main CCL11 receptor, namely CCR3. Accordingly, the congenic strain also showed significantly increased Ccr3 mRNA expression in lymph nodes. In the CNS, the main producers of CCL11 were neurons, whereas CCR3 was detected on neurons and CSF producing ependymal cells. This corresponded to increased levels of CCL11 protein in the cerebrospinal fluid of the congenic rats. Increased intrathecal production of CCL11 in congenic rats was accompanied by a tighter blood brain barrier, reflected by more occludin⁺ blood vessels. In addition, the congenic strain showed a reduced antigen specific response and a predominant anti-inflammatory Th2 phenotype. These results indicate novel mechanisms in the genetic regulation of neuroinflammation.

A resource for the simultaneous high-resolution mapping of multiple quantitative trait loci in rats: the NIH heterogeneous stock

The laboratory rat (Rattus norvegicus) is a key tool for the study of medicine and pharmacology for human health. A large database of phenotypes for integrated fields such as cardiovascular, neuroscience, and exercise physiology exists in the literature. However, the molecular characterization of the genetic loci that give rise to variation in these traits has proven to be difficult. Here we show how one obstacle to progress, the fine-mapping of quantitative trait loci (QTL), can be overcome by using an outbred population of rats. By use of a genetically heterogeneous stock of rats, we map a locus contributing to variation in a fear-related measure (two-way active avoidance in the shuttle box) to a region on chromosome 5 containing nine genes. By establishing a protocol measuring multiple phenotypes including immunology, neuroinflammation, and hematology, as well as cardiovascular, metabolic, and behavioral traits, we establish the rat HS as a new resource for the fine-mapping of QTLs contributing to variation in complex traits of biomedical relevance.