Compartment-specific protein interactions in beryllium lung disease.

The study provides insights into proteins that may be relevant in BeS and CBD. It provides a framework to investigate the global changes in lung compartment-specific inflammatory cells to better understand the potential interplay of proteins in CBD. https://bit.ly/3PLNTXC.

Genome-wide association study identifies multiple HLA loci for sarcoidosis susceptibility.

Sarcoidosis is a complex systemic disease. Our study aimed to (1) identify novel alleles associated with sarcoidosis susceptibility; (2) provide an in-depth evaluation of HLA alleles and sarcoidosis susceptibility and (3) integrate genetic and transcription data to identify risk loci that may more directly impact disease pathogenesis. We report a genome-wide association study of 1335 sarcoidosis cases and 1264 controls of European descent (EA) and investigate associated alleles in a study of African Americans (AA: 1487 cases and 1504 controls). The EA and AA cohort was recruited from multiple United States sites. HLA alleles were imputed and tested for association with sarcoidosis susceptibility. Expression quantitative locus and colocalization analysis were performed using a subset of subjects with transcriptome data. Forty-nine SNPs in the HLA region in HLA-DRA, -DRB9, -DRB5, -DQA1 and BRD2 genes were significantly associated with sarcoidosis susceptibility in EA, rs3129888 was also a risk variant for sarcoidosis in AA. Classical HLA alleles DRB1*0101, DQA1*0101 and DQB1*0501, which are highly correlated, were also associated with sarcoidosis. rs3135287 near HLA-DRA was associated with HLA-DRA expression in peripheral blood mononuclear cells and bronchoalveolar lavage from subjects and lung tissue and whole blood from GTEx. We identified six novel SNPs (out of the seven SNPs representing the 49 significant SNPs) and nine HLA alleles associated with sarcoidosis susceptibility in the largest EA population. We also replicated our findings in an AA population. Our study reiterates the potential role of antigen recognition and/or presentation HLA class II genes in sarcoidosis pathogenesis.

Standardization of flow cytometry and cell sorting to enable a transcriptomic analysis in a multi-site sarcoidosis study.

The contribution and regulation of various CD4+ T cell lineages that occur with remitting vs progressive courses in sarcoidosis are poorly understood. We developed a multiparameter flow cytometry panel to sort these CD4+ T cell lineages followed by measurement of their functional potential using RNA-sequencing analysis at six-month intervals across multiple study sites. To obtain good quality RNA for sequencing, we relied on chemokine receptor expression to identify and sort lineages. To minimize gene expression changes induced by perturbations of T cells and avoid protein denaturation caused by freeze/thaw cycles, we optimized our protocols using freshly isolated samples at each study site. To accomplish this study, we had to overcome significant standardization challenges across multiple sites. Here, we detail standardization considerations for cell processing, flow staining, data acquisition, sorting parameters, and RNA quality control analysis that were performed as part of the NIH-sponsored, multi-center study, BRonchoscopy at Initial sarcoidosis diagnosis Targeting longitudinal Endpoints (BRITE). After several rounds of iterative optimization, we identified the following aspects as critical for successful standardization: 1) alignment of PMT voltages across sites using CS&T/rainbow bead technology; 2) a single template created in the cytometer program that was used by all sites to gate cell populations during data acquisition and cell sorting; 3) use of standardized lyophilized flow cytometry staining cocktails to reduce technical error during processing; 4) development and implementation of a standardized Manual of Procedures. After standardization of cell sorting, we were able to determine the minimum number of sorted cells necessary for next generation sequencing through analysis of RNA quality and quantity from sorted T cell populations. Overall, we found that implementing a multi-parameter cell sorting with RNA-seq analysis clinical study across multiple study sites requires iteratively tested standardized procedures to ensure comparable and high-quality results.

Multi-Omic Signatures of Sarcoidosis and Progression in Bronchoalveolar Lavage Cells.

Introduction: Sarcoidosis is a heterogeneous, granulomatous disease that can prove difficult to diagnose, with no accurate biomarkers of disease progression. Therefore, we profiled and integrated the DNA methylome, mRNAs, and microRNAs to identify molecular changes associated with sarcoidosis and disease progression that might illuminate underlying mechanisms of disease and potential genomic biomarkers. Methods: Bronchoalveolar lavage cells from 64 sarcoidosis subjects and 16 healthy controls were used. DNA methylation was profiled on Illumina HumanMethylationEPIC arrays, mRNA by RNA-sequencing, and miRNAs by small RNA-sequencing. Linear models were fit to test for effect of diagnosis and phenotype, adjusting for age, sex, and smoking. We built a supervised multi-omics model using a subset of features from each dataset. Results: We identified 46,812 CpGs, 1,842 mRNAs, and 5 miRNAs associated with sarcoidosis versus controls and 1 mRNA, SEPP1 - a protein that supplies selenium to cells, associated with disease progression. Our integrated model emphasized the prominence of the PI3K/AKT1 pathway in sarcoidosis, which is important in T cell and mTOR function. Novel immune related genes and miRNAs including LYST, RGS14, SLFN12L, and hsa-miR-199b-5p, distinguished sarcoidosis from controls. Our integrated model also demonstrated differential expression/methylation of IL20RB, ABCC11, SFSWAP, AGBL4, miR-146a-3p, and miR-378b between non-progressive and progressive sarcoidosis. Conclusions: Leveraging the DNA methylome, transcriptome, and miRNA-sequencing in sarcoidosis BAL cells, we detected widespread molecular changes associated with disease, many which are involved in immune response. These molecules may serve as diagnostic/prognostic biomarkers and/or drug targets, although future testing will be required for confirmation.

Multiomic Signatures of Chronic Beryllium Disease Bronchoalveolar Lavage Cells Relate to T-Cell Function and Innate Immunity.

Chronic beryllium disease (CBD) is a Th1 granulomatous lung disease preceded by sensitization to beryllium (BeS). We profiled the methylome, transcriptome, and selected proteins in the lung to identify molecular signatures and networks associated with BeS and CBD. BAL cell DNA and RNA were profiled using microarrays from CBD (n = 30), BeS (n = 30), and control subjects (n = 12). BAL fluid proteins were measured using Olink Immune Response Panel proteins from CBD (n = 22) and BeS (n = 22) subjects. Linear models identified features associated with CBD, adjusting for covariation and batch effects. Multiomic integration methods identified correlated features between datasets. We identified 1,546 differentially expressed genes in CBD versus control subjects and 204 in CBD versus BeS. Of the 101 shared transcripts, 24 have significant cis relationships between gene expression and DNA methylation, assessed using expression quantitative trait methylation analysis, including genes not previously identified in CBD. A multiomic model of top DNA methylation and gene expression features demonstrated that the first component separated CBD from other samples and the second component separated control subjects from remaining samples. The top features on component one were enriched for T-lymphocyte function, and the top features on component two were enriched for innate immune signaling. We identified six differentially abundant proteins in CBD versus BeS, with two (SIT1 and SH2D1A) selected as important RNA features in the multiomic model. Our integrated analysis of DNA methylation, gene expression, and proteins in the lung identified multiomic signatures of CBD that differentiated it from BeS and control subjects.

Genomic biomarkers in chronic beryllium disease and sarcoidosis.

Background Previous gene expression studies have identified genes IFNγ, TNFα, RNase 3, CXCL9, and CD55 as potential biomarkers for sarcoidosis and/or chronic beryllium disease (CBD). We hypothesized that differential expression of these genes could function as diagnostic biomarkers for sarcoidosis and CBD, and prognostic biomarkers for sarcoidosis. Study Design/Methods We performed RT-qPCR on whole blood samples from CBD (n = 132), beryllium sensitized (BeS) (n = 109), and sarcoidosis (n = 99) cases and non-diseased controls (n = 97) to determine differential expression of target genes. We then performed logistic regression modeling and generated ROC curves to determine which genes could most accurately differentiate: 1) CBD versus sarcoidosis 2) CBD versus BeS 3) sarcoidosis versus controls 4) non-progressive versus progressive sarcoidosis. Results CD55 and TNFα were significantly upregulated, while CXCL9 was significantly downregulated in CBD compared to sarcoidosis (p < 0.05). The ROC curve from the logistic regression model demonstrated high discriminatory ability of the combination of CD55, TNFα, and CXCL9 to distinguish between CBD and sarcoidosis with an AUC of 0.98. CD55 and TNFα were significantly downregulated in sarcoidosis compared to controls (p < 0.05). The ROC curve from the model showed a reasonable discriminatory ability of CD55 and TNFα to distinguish between sarcoidosis and controls with an AUC of 0.86. There was no combination of genes that could accurately differentiate between CBD and BeS or sarcoidosis phenotypes. Interpretation CD55, TNFα and CXCL9 expression levels can accurately differentiate between CBD and sarcoidosis, while CD55 and TNFα expression levels can accurately differentiate sarcoidosis and controls.

DNA Methylation Changes in Lung Immune Cells Are Associated with Granulomatous Lung Disease.

Epigenetic marks are likely to explain variability of response to antigen in granulomatous lung disease. The objective of this study was to identify DNA methylation and gene expression changes associated with chronic beryllium disease (CBD) and sarcoidosis in lung cells obtained by BAL. BAL cells from CBD (n = 8), beryllium-sensitized (n = 8), sarcoidosis (n = 8), and additional progressive sarcoidosis (n = 9) and remitting (n = 15) sarcoidosis were profiled on the Illumina 450k methylation and Affymetrix/Agilent gene expression microarrays. Statistical analyses were performed to identify DNA methylation and gene expression changes associated with CBD, sarcoidosis, and disease progression in sarcoidosis. DNA methylation array findings were validated by pyrosequencing. We identified 52,860 significant (P < 0.005 and q < 0.05) CpGs associated with CBD; 2,726 CpGs near 1,944 unique genes have greater than 25% methylation change. A total of 69% of differentially methylated genes are significantly (q < 0.05) differentially expressed in CBD, with many canonical inverse relationships of methylation and expression in genes critical to T-helper cell type 1 differentiation, chemokines and their receptors, and other genes involved in immunity. Testing of these CBD-associated CpGs in sarcoidosis reveals that methylation changes only approach significance, but are methylated in the same direction, suggesting similarities between the two diseases with more heterogeneity in sarcoidosis that limits power with the current sample size. Analysis of progressive versus remitting sarcoidosis identified 15,215 CpGs (P < 0.005 and q < 0.05), but only 801 of them have greater than 5% methylation change, demonstrating that DNA methylation marks of disease progression changes are more subtle. Our study highlights the significance of epigenetic marks in lung immune response in granulomatous lung disease.

Polymorphism of FCGR3A gene in chronic beryllium disease.

Previously we showed that alveolar macrophages (AMs) from patients with chronic beryllium disease (CBD) and beryllium sensitization (BeS) demonstrated significantly greater cell surface CD16 (encoded by the FCGR3A gene) than controls. We hypothesized that these differences were related to polymorphisms in the FCGR3A gene. This study was to determine the association between FCGR3A polymorphisms in CBD, BeS versus controls as well as clinical data, providing potential information about disease pathogenesis, risk, and activity. A total of 189 CBD/154 BeS/150 controls (92 Be-exposed non-diseased and 58 healthy controls) were included in this study. Sequence-specific primers polymerase chain reaction (PCR-SSP) was used to determine FCGR3A 158V/F polymorphisms. We found significantly higher frequencies of the 158V allele (OR: 1.60 (CI: 1.17-2.19), p = 0.004) and 158VV homozygotes (OR: 2.97 (CI: 1.48-5.97) p = 0.007) in CBD versus controls. No differences were found in the frequencies of FCGR3A alleles or genotypes between BeS versus controls and CBD versus BeS. Average changes in exercise testing maximum workload (Wlm), maximum oxygen consumption (VO2m), and diffusion capacity of carbon monoxide (DLCO) demonstrated greater decline over time in those CBD cases with the 158VV gene, modeled between 10 and 40 years from first beryllium exposure. The FCGR3A V158F polymorphism is associated with CBD compared to BeS and controls and may impact lung function in CBD.

Gene-environment interactions influence airways function in laboratory animal workers.

BACKGROUND: Most diseases, including asthma, result from the interaction between environmental exposures and genetic variants. Functional variants of CD14 negatively affect lung function in farm workers and children exposed to animal allergens and endotoxin. OBJECTIVE: We hypothesized that CD14 polymorphisms interact with inhaled endotoxin, mouse allergen, or both to decrease airways function in laboratory animal workers. METHODS: Three hundred sixty-nine Caucasian workers completed a symptom and work exposure questionnaire, skin prick testing, and spirometry. Individual exposure estimates for endotoxin and murine allergen were calculated by weighting task-based breathing zone concentrations by time reported for each task and length of time in the current job. Real-time PCR was used to assess CD14/-1619, -550, and -159 alleles. Multiple linear regression predicting airways function included an interaction term between genotype and exposure. RESULTS: Workers at the highest quartile of the natural log-transformed cumulative endotoxin exposure and with the endotoxin-responsive CD14/-1619 G allele had significantly lower FEV(1) and forced expiratory flow, midexpiratory phase (FEF(25-75)) percent predicted compared with workers with an AA genotype, with no significant differences noted at lower endotoxin levels for either genotype. The gene-environment effect was marked for atopic workers. Laboratory animal allergy, murine allergen exposure, CD14/-159 or -550 genotype, and a gene-exposure interaction term for these genotypes and exposures did not predict changes in lung function. CONCLUSIONS: A significant gene-environment interaction affects airways function in laboratory animal workers. More highly endotoxin-exposed workers with CD14/-1619G alleles have significantly lower FEV(1) and FEF(25-75) percent predicted than those with CD14/-1619AA alleles. Atopic workers are particularly affected by cumulative endotoxin exposures.