Unravelling the collateral damage of antibiotics on gut bacteria.

Antibiotics are used to fight pathogens but also target commensal bacteria, disturbing the composition of gut microbiota and causing dysbiosis and disease1. Despite this well-known collateral damage, the activity spectrum of different antibiotic classes on gut bacteria remains poorly characterized. Here we characterize further 144 antibiotics from a previous screen of more than 1,000 drugs on 38 representative human gut microbiome species2. Antibiotic classes exhibited distinct inhibition spectra, including generation dependence for quinolones and phylogeny independence for ß-lactams. Macrolides and tetracyclines, both prototypic bacteriostatic protein synthesis inhibitors, inhibited nearly all commensals tested but also killed several species. Killed bacteria were more readily eliminated from in vitro communities than those inhibited. This species-specific killing activity challenges the long-standing distinction between bactericidal and bacteriostatic antibiotic classes and provides a possible explanation for the strong effect of macrolides on animal3-5 and human6,7 gut microbiomes. To mitigate this collateral damage of macrolides and tetracyclines, we screened for drugs that specifically antagonized the antibiotic activity against abundant Bacteroides species but not against relevant pathogens. Such antidotes selectively protected Bacteroides species from erythromycin treatment in human-stool-derived communities and gnotobiotic mice. These findings illluminate the activity spectra of antibiotics in commensal bacteria and suggest strategies to circumvent their adverse effects on the gut microbiota.

Bioaccumulation of therapeutic drugs by human gut bacteria.

Bacteria in the gut can modulate the availability and efficacy of therapeutic drugs. However, the systematic mapping of the interactions between drugs and bacteria has only started recently1 and the main underlying mechanism proposed is the chemical transformation of drugs by microorganisms (biotransformation). Here we investigated the depletion of 15 structurally diverse drugs by 25 representative strains of gut bacteria. This revealed 70 bacteria-drug interactions, 29 of which had not to our knowledge been reported before. Over half of the new interactions can be ascribed to bioaccumulation; that is, bacteria storing the drug intracellularly without chemically modifying it, and in most cases without the growth of the bacteria being affected. As a case in point, we studied the molecular basis of bioaccumulation of the widely used antidepressant duloxetine by using click chemistry, thermal proteome profiling and metabolomics. We find that duloxetine binds to several metabolic enzymes and changes the metabolite secretion of the respective bacteria. When tested in a defined microbial community of accumulators and non-accumulators, duloxetine markedly altered the composition of the community through metabolic cross-feeding. We further validated our findings in an animal model, showing that bioaccumulating bacteria attenuate the behavioural response of Caenorhabditis elegans to duloxetine. Together, our results show that bioaccumulation by gut bacteria may be a common mechanism that alters drug availability and bacterial metabolism, with implications for microbiota composition, pharmacokinetics, side effects and drug responses, probably in an individual manner.

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.

Consistency across multi-omics layers in a drug-perturbed gut microbial community.

Multi-omics analyses are used in microbiome studies to understand molecular changes in microbial communities exposed to different conditions. However, it is not always clear how much each omics data type contributes to our understanding and whether they are concordant with each other. Here, we map the molecular response of a synthetic community of 32 human gut bacteria to three non-antibiotic drugs by using five omics layers (16S rRNA gene profiling, metagenomics, metatranscriptomics, metaproteomics and metabolomics). We find that all the omics methods with species resolution are highly consistent in estimating relative species abundances. Furthermore, different omics methods complement each other for capturing functional changes. For example, while nearly all the omics data types captured that the antipsychotic drug chlorpromazine selectively inhibits Bacteroidota representatives in the community, the metatranscriptome and metaproteome suggested that the drug induces stress responses related to protein quality control. Metabolomics revealed a decrease in oligosaccharide uptake, likely caused by Bacteroidota depletion. Our study highlights how multi-omics datasets can be utilized to reveal complex molecular responses to external perturbations in microbial communities.

Developmental drugs for sarcoidosis.

Sarcoidosis is a multi-organ granulomatous inflammatory disease of unknown etiology. Over 50% of patients will require treatment at some point in their disease and 10%-30% will develop a chronic progressive disease with pulmonary fibrosis leading to significant morbidity and mortality. Recently published guidelines recommend immunosuppressive therapy for sarcoidosis patients at risk of increased disease-related morbidity and mortality, and in whom disease has negatively impacted quality of life. Prednisone the currently recommended first line therapy is associated with significant toxicity however none of the other guideline recommended steroid sparing therapy is approved by regulatory agencies for use in sarcoidosis, and data in support of their use is weak. For patients with severe refractory disease requiring prolonged therapy, treatment options are limited. The need for expanding treatment options in sarcoidosis has been emphasized. Well conducted large, randomized trials evaluating currently available therapeutic options as well as novel pathways for targeting disease are necessary to better guide treatment decisions. These trials will not be without significant challenges. Sarcoidosis is a rare disease with heterogenous presentation and variable progression and clinical outcome. There are no universally agreed upon biomarkers of disease activity and measurement of outcomes is confounded by the need to balance patient centric measures and objective measures of disease activity. Our paper provides an update on developmental drugs in sarcoidosis and outlines several novel pathways that may be targeted for future drug development. Currently available trials are highlighted and ongoing challenges to drug development and clinical trial design are briefly discussed.

Extensive impact of non-antibiotic drugs on human gut bacteria.

A few commonly used non-antibiotic drugs have recently been associated with changes in gut microbiome composition, but the extent of this phenomenon is unknown. Here, we screened more than 1,000 marketed drugs against 40 representative gut bacterial strains, and found that 24% of the drugs with human targets, including members of all therapeutic classes, inhibited the growth of at least one strain in vitro. Particular classes, such as the chemically diverse antipsychotics, were overrepresented in this group. The effects of human-targeted drugs on gut bacteria are reflected on their antibiotic-like side effects in humans and are concordant with existing human cohort studies. Susceptibility to antibiotics and human-targeted drugs correlates across bacterial species, suggesting common resistance mechanisms, which we verified for some drugs. The potential risk of non-antibiotics promoting antibiotic resistance warrants further exploration. Our results provide a resource for future research on drug-microbiome interactions, opening new paths for side effect control and drug repurposing, and broadening our view of antibiotic resistance.

Alpha-1 antitrypsin inhibits fractalkine-mediated monocyte-lung endothelial cell interactions.

Chronic obstructive pulmonary disease (COPD) is characterized by nonresolving inflammation fueled by breach in the endothelial barrier and leukocyte recruitment into the airspaces. Among the ligand-receptor axes that control leukocyte recruitment, the full-length fractalkine ligand (CX3CL1)-receptor (CX3CR1) ensures homeostatic endothelial-leukocyte interactions. Cigarette smoke (CS) exposure and respiratory pathogens increase expression of endothelial sheddases, such as a-disintegrin-and-metalloproteinase-domain 17 (ADAM17, TACE), inhibited by the anti-protease α-1 antitrypsin (AAT). In the systemic endothelium, TACE cleaves CX3CL1 to release soluble CX3CL1 (sCX3CL1). During CS exposure, it is not known whether AAT inhibits sCX3CL1 shedding and CX3CR1+ leukocyte transendothelial migration across lung microvasculature. We investigated the mechanism of sCX3CL1 shedding, its role in endothelial-monocyte interactions, and AAT effect on these interactions during acute inflammation. We used two, CS and lipopolysaccharide (LPS) models of acute inflammation in transgenic Cx3cr1gfp/gfp mice and primary human endothelial cells and monocytes to study sCX3CL1-mediated CX3CR1+ monocyte adhesion and migration. We measured sCX3CL1 levels in plasma and bronchoalveolar lavage (BALF) of individuals with COPD. Both sCX3CL1 shedding and CX3CR1+ monocytes transendothelial migration were triggered by LPS and CS exposure in mice, and were significantly attenuated by AAT. The inhibition of monocyte-endothelial adhesion and migration by AAT was TACE-dependent. Compared with healthy controls, sCX3CL1 levels were increased in plasma and BALF of individuals with COPD, and were associated with clinical parameters of emphysema. Our results indicate that inhibition of sCX3CL1 as well as AAT augmentation may be effective approaches to decrease excessive monocyte lung recruitment during acute and chronic inflammatory states.NEW & NOTEWORTHY Our novel findings that AAT and other inhibitors of TACE, the sheddase that controls full-length fractalkine (CX3CL1) endothelial expression, may provide fine-tuning of the CX3CL1-CX3CR1 axis specifically involved in endothelial-monocyte cross talk and leukocyte recruitment to the alveolar space, suggests that AAT and inhibitors of sCX3CL1 signaling may be harnessed to reduce lung inflammation.

The Cyanobacterial "Nutraceutical" Phycocyanobilin Inhibits Cysteine Protease Legumain.

The blue biliprotein phycocyanin, produced by photo-autotrophic cyanobacteria including spirulina (Arthrospira) and marketed as a natural food supplement or "nutraceutical," is reported to have anti-inflammatory, antioxidant, immunomodulatory, and anticancer activity. These diverse biological activities have been specifically attributed to the phycocyanin chromophore, phycocyanobilin (PCB). However, the mechanism of action of PCB and the molecular targets responsible for the beneficial properties of PCB are not well understood. We have developed a procedure to rapidly cleave the PCB pigment from phycocyanin by ethanolysis and then characterized it as an electrophilic natural product that interacts covalently with thiol nucleophiles but lacks any appreciable cytotoxicity or antibacterial activity against common pathogens and gut microbes. We then designed alkyne-bearing PCB probes for use in chemical proteomics target deconvolution studies. Target identification and validation revealed the cysteine protease legumain (also known as asparaginyl endopeptidase, AEP) to be a target of PCB. Inhibition of this target may account for PCB's diverse reported biological activities.

Genetic predisposition to sarcoidosis.

Sarcoidosis is a complex systemic disease with clinical heterogeneity based on varying phenotypes and natural history. The detailed etiology of sarcoidosis remains unknown, but genetic predisposition as well as environmental exposures play a significant role in disease pathogenesis. We performed a comprehensive review of germline genetic (DNA) and transcriptomic (RNA) studies of sarcoidosis, including both previous studies and more recent findings. In this review, we provide an assessment of the following: genetic variants in sarcoidosis susceptibility and phenotypes, ancestry- and sex-specific genetic variants in sarcoidosis, shared genetic architecture between sarcoidosis and other diseases, and gene-environment interactions in sarcoidosis. We also highlight the unmet needs in sarcoidosis genetic studies, including the pressing requirement to include diverse populations and have consistent definitions of phenotypes in the sarcoidosis research community to help advance the application of genetic predisposition to sarcoidosis disease risk and manifestations.

Salt-responsive gut commensal modulates TH17 axis and disease.

A Western lifestyle with high salt consumption can lead to hypertension and cardiovascular disease. High salt may additionally drive autoimmunity by inducing T helper 17 (TH17) cells, which can also contribute to hypertension. Induction of TH17 cells depends on gut microbiota; however, the effect of salt on the gut microbiome is unknown. Here we show that high salt intake affects the gut microbiome in mice, particularly by depleting Lactobacillus murinus. Consequently, treatment of mice with L. murinus prevented salt-induced aggravation of actively induced experimental autoimmune encephalomyelitis and salt-sensitive hypertension by modulating TH17 cells. In line with these findings, a moderate high-salt challenge in a pilot study in humans reduced intestinal survival of Lactobacillus spp., increased TH17 cells and increased blood pressure. Our results connect high salt intake to the gut-immune axis and highlight the gut microbiome as a potential therapeutic target to counteract salt-sensitive conditions.

Imaging of traumatic mandibular fractures in young adults using CT-like MRI: a feasibility study.

OBJECTIVES: To assess and compare the diagnostic performance of CT-like images based on a three- dimensional (3D) T1-weighted spoiled gradient-echo sequence (3D T1 GRE) with CT in patients with acute traumatic fractures of the mandible. MATERIALS AND METHODS: Subjects with acute mandibular fractures diagnosed on conventional CT were prospectively recruited and received an additional 3 T MRI with a CT-like 3D T1 GRE sequence. The images were assessed by two radiologists with regard to fracture localization, degree of dislocation, and number of fragments. Bone to soft tissue contrast, diagnostic confidence, artifacts, and overall image quality were rated using a five-point Likert-scale. Agreement of measurements was assessed using an independent t-test. RESULTS: Fourteen subjects and 22 fracture sites were included (26 ± 3.9 years; 4 females, 10 males). All traumatic fractures were accurately detected on CT-like MRI (n = 22, κ 1.00 (95% CI 1.00-1.00)). There was no statistically significant difference in the assessment of the fracture dislocation (axial mean difference (MD) 0.06 mm, p = 0.93, coronal MD, 0.08 mm, p = 0.89 and sagittal MD, 0.04 mm, p = 0.96). The agreement for the fracture classification as well as the inter- and intra-rater agreement was excellent (range κ 0.92-0.98 (95% CI 0.96-0.99)). CONCLUSION: Assessment of mandibular fractures was feasible and accurate using CT-like MRI based on a 3D T1 GRE sequence and is comparable to conventional CT. CLINICAL RELEVANCE: For the assessment of acute mandibular fractures, CT-like MRI might become a useful alternative to CT in order to reduce radiation exposure particularly in young patients.

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.

Prognostic accuracy of a peripheral blood transcriptome signature in chronic hypersensitivity pneumonitis.

The prognostic value of peripheral blood mononuclear cell (PBMC) expression profiles, when used in patients with chronic hypersensitivity pneumonitis (CHP), as an adjunct to traditional clinical assessment is unknown. RNA-seq analysis on PBMC from 37 patients with CHP at initial presentation determined that (1) 74 differentially expressed transcripts at a 10% false discovery rate distinguished those with (n=10) and without (n=27) disease progression, defined as absolute FVC and/or diffusing capacity of the lungs for carbon monoxide (DLCO) decline of ≥10% and increased fibrosis on chest CT images within 24 months, and (2) classification models based on gene expression and clinical factors strongly outperform models based solely on clinical factors (baseline FVC%, DLCO% and chest CT fibrosis).

Towards a mechanistic understanding of reciprocal drug-microbiome interactions.

Broad-spectrum antibiotics target multiple gram-positive and gram-negative bacteria, and can collaterally damage the gut microbiota. Yet, our knowledge of the extent of damage, the antibiotic activity spectra, and the resistance mechanisms of gut microbes is sparse. This limits our ability to mitigate microbiome-facilitated spread of antibiotic resistance. In addition to antibiotics, non-antibiotic drugs affect the human microbiome, as shown by metagenomics as well as in vitro studies. Microbiome-drug interactions are bidirectional, as microbes can also modulate drugs. Chemical modifications of antibiotics mostly function as antimicrobial resistance mechanisms, while metabolism of non-antibiotics can also change the drugs' pharmacodynamic, pharmacokinetic, and toxic properties. Recent studies have started to unravel the extensive capacity of gut microbes to metabolize drugs, the mechanisms, and the relevance of such events for drug treatment. These findings raise the question whether and to which degree these reciprocal drug-microbiome interactions will differ across individuals, and how to take them into account in drug discovery and precision medicine. This review describes recent developments in the field and discusses future study areas that will benefit from systems biology approaches to better understand the mechanistic role of the human gut microbiota in drug actions.

Clinical phenotyping in sarcoidosis using cluster analysis.

BACKGROUND: Most phenotyping paradigms in sarcoidosis are based on expert opinion; however, no paradigm has been widely adopted because of the subjectivity in classification. We hypothesized that cluster analysis could be performed on common clinical variables to define more objective sarcoidosis phenotypes. METHODS: We performed a retrospective cohort study of 554 sarcoidosis cases to identify distinct phenotypes of sarcoidosis based on 29 clinical features. Model-based clustering was performed using the VarSelLCM R package and the Integrated Completed Likelihood (ICL) criteria were used to estimate number of clusters. To identify features associated with cluster membership, features were ranked based on variable importance scores from the VarSelLCM model, and additional univariate tests (Fisher's exact test and one-way ANOVA) were performed using q-values correcting for multiple testing. The Wasfi severity score was also compared between clusters. RESULTS: Cluster analysis resulted in 6 sarcoidosis phenotypes. Salient characteristics for each cluster are as follows: Phenotype (1) supranormal lung function and majority Scadding stage 2/3; phenotype (2) supranormal lung function and majority Scadding stage 0/1; phenotype (3) normal lung function and split Scadding stages between 0/1 and 2/3; phenotype (4) obstructive lung function and majority Scadding stage 2/3; phenotype (5) restrictive lung function and majority Scadding stage 2/3; phenotype (6) mixed obstructive and restrictive lung function and mostly Scadding stage 4. Although there were differences in the percentages, all Scadding stages were encompassed by all of the phenotypes, except for phenotype 1, in which none were Scadding stage 4. Clusters 4, 5, 6 were significantly more likely to have ever been on immunosuppressive treatment and had higher Wasfi disease severity scores. CONCLUSIONS: Cluster analysis produced 6 sarcoidosis phenotypes that demonstrated less severe and severe phenotypes. Phenotypes 1, 2, 3 have less lung function abnormalities, a lower percentage on immunosuppressive treatment and lower Wasfi severity scores. Phenotypes 4, 5, 6 were characterized by lung function abnormalities, more parenchymal abnormalities, an increased percentage on immunosuppressive treatment and higher Wasfi severity scores. These data support using cluster analysis as an objective and clinically useful way to phenotype sarcoidosis subjects and to empower clinicians to identify those with more severe disease versus those who have less severe disease, independent of Scadding stage.

Occupational exposures and sarcoidosis: current understanding and knowledge gaps.

PURPOSE OF REVIEW: Sarcoidosis is an idiopathic granulomatous disease that primarily affects the lungs. Several lines of evidence suggest that occupational exposures are associated with disease risk. This review critically evaluates studies using the Bradford Hill criteria for causation to determine if a causal relationship can be established between occupational exposure and sarcoidosis. RECENT FINDINGS: Large epidemiological studies have proposed multiple occupational exposures associated with sarcoidosis but lack consistency of results. Many convincing studies demonstrate an association between World Trade Center (WTC) dust and sarcoidosis, which illustrates a causal relationship based on the fulfillment of the Bradford Hill criteria. Studies describing an association between silica/metals and sarcoidosis are intriguing but fulfill a limited number of the Bradford Hill criteria and warrant further investigation before a causal relationship can be determined. Finally, we also discuss preliminary studies associating sarcoidosis phenotypes with specific occupational exposures. SUMMARY: Using the Bradford Hill criteria for causation, we demonstrate that WTC dust has a causative relationship with sarcoidosis, which reinforces the theory that sarcoidosis is an exposure-related disease. More research is needed to determine other specific occupational exposures causing disease.

HLA-DPB1 E69 genotype and exposure in beryllium sensitisation and disease.

OBJECTIVES: Human leukocyte antigen-DP beta 1 (HLA-DPB1) with a glutamic acid at the 69th position of the ß chain (E69) genotype and inhalational beryllium exposure individually contribute to risk of chronic beryllium disease (CBD) and beryllium sensitisation (BeS) in exposed individuals. This retrospective nested case-control study assessed the contribution of genetics and exposure in the development of BeS and CBD. METHODS: Workers with BeS (n=444), CBD (n=449) and beryllium-exposed controls (n=890) were enrolled from studies conducted at nuclear weapons and primary beryllium manufacturing facilities. Lifetime-average beryllium exposure estimates were based on workers' job questionnaires and historical and industrial hygienist exposure estimates, blinded to genotype and case status. Genotyping was performed using sequence-specific primer-PCR. Logistic regression models were developed allowing for over-dispersion, adjusting for workforce, race, sex and ethnicity. RESULTS: Having no E69 alleles was associated with lower odds of both CBD and BeS; every additional E69 allele increased odds for CBD and BeS. Increasing exposure was associated with lower odds of BeS. CBD was not associated with exposure as compared to controls, yet the per cent of individuals with CBD versus BeS increased with increasing exposure. No evidence of a gene-by-exposure interaction was found for CBD or BeS. CONCLUSIONS: Risk of CBD increases with E69 allele frequency and increasing exposure, although no gene by environment interaction was found. A decreased risk of BeS with increasing exposure and lack of exposure response in CBD cases may be due to the limitations of reconstructed exposure estimates. Although reducing exposure may not prevent BeS, it may reduce CBD and the associated health effects, especially in those carrying E69 alleles.

Realistic biomarkers from plasma extracellular vesicles for detection of beryllium exposure.

PURPOSE: Exposures related to beryllium (Be) are an enduring concern among workers in the nuclear weapons and other high-tech industries, calling for regular and rigorous biological monitoring. Conventional biomonitoring of Be in urine is not informative of cumulative exposure nor health outcomes. Biomarkers of exposure to Be based on non-invasive biomonitoring could help refine disease risk assessment. In a cohort of workers with Be exposure, we employed blood plasma extracellular vesicles (EVs) to discover novel biomarkers of exposure to Be. METHODS: EVs were isolated from plasma using size-exclusion chromatography and subjected to mass spectrometry-based proteomics. A protein-based classifier was developed using LASSO regression and validated by ELISA. RESULTS: We discovered a dual biomarker signature comprising zymogen granule protein 16B and putative protein FAM10A4 that differentiated between Be-exposed and -unexposed subjects. ELISA-based quantification of the biomarkers in an independent cohort of samples confirmed higher expression of the signature in the Be-exposed group, displaying high predictive accuracy (AUROC = 0.919). Furthermore, the biomarkers efficiently discriminated high- and low-exposure groups (AUROC = 0.749). CONCLUSIONS: This is the first report of EV biomarkers associated with Be exposure and exposure levels. The biomarkers could be implemented in resource-limited settings for Be exposure assessment.

Transcriptomics of bronchoalveolar lavage cells identifies new molecular endotypes of sarcoidosis.

BACKGROUND: Sarcoidosis is a multisystem granulomatous disease of unknown origin with a variable and often unpredictable course and pattern of organ involvement. In this study we sought to identify specific bronchoalveolar lavage (BAL) cell gene expression patterns indicative of distinct disease phenotypic traits. METHODS: RNA sequencing by Ion Torrent Proton was performed on BAL cells obtained from 215 well-characterised patients with pulmonary sarcoidosis enrolled in the multicentre Genomic Research in Alpha-1 Antitrypsin Deficiency and Sarcoidosis (GRADS) study. Weighted gene co-expression network analysis and nonparametric statistics were used to analyse genome-wide BAL transcriptome. Validation of results was performed using a microarray expression dataset of an independent sarcoidosis cohort (Freiburg, Germany; n=50). RESULTS: Our supervised analysis found associations between distinct transcriptional programmes and major pulmonary phenotypic manifestations of sarcoidosis including T-helper type 1 (Th1) and Th17 pathways associated with hilar lymphadenopathy, transforming growth factor-ß1 (TGFB1) and mechanistic target of rapamycin (MTOR) signalling with parenchymal involvement, and interleukin (IL)-7 and IL-2 with airway involvement. Our unsupervised analysis revealed gene modules that uncovered four potential sarcoidosis endotypes including hilar lymphadenopathy with increased acute T-cell immune response; extraocular organ involvement with PI3K activation pathways; chronic and multiorgan disease with increased immune response pathways; and multiorgan involvement, with increased IL-1 and IL-18 immune and inflammatory responses. We validated the occurrence of these endotypes using gene expression, pulmonary function tests and cell differentials from Freiburg. CONCLUSION: Taken together, our results identify BAL gene expression programmes that characterise major pulmonary sarcoidosis phenotypes and suggest the presence of distinct disease molecular endotypes.

Genetics in sarcoidosis.

PURPOSE OF REVIEW: Epidemiological and clinical observations as well as familial clustering support the existence of a genetic predisposition to sarcoidosis. In this article, we review the most recent findings in genetics of sarcoidosis and discuss how the identification of risk alleles may help advancing our understanding of disease etiology and development. RECENT FINDINGS: Genetic studies of sarcoidosis phenotypes have identified novel and ancestry-specific associations. Gene-environment interaction studies highlighted the importance of integrating genetic information when assessing the relationship between sarcoidosis and environmental exposures. A case-control-family study revealed that the heritability of sarcoidosis is only 49%, suggesting the existence of additional important contributors to disease risk. The application of whole-exome sequencing has identified associations with disease activity and prognosis. Finally, gene expression studies of circulating immune cells have identified shared and unique pathways between sarcoidosis and other granulomatous diseases. SUMMARY: Sarcoidosis genetic research has led to the identification of a number of associations with both sarcoidoses per se and disease phenotypes. Newer sequencing technologies are likely to increase the number of genetic variants associated with sarcoidosis. However, studying phenotypically and ethnically homogeneous patient subsets remains critically important regardless of the genetic approach used.