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.

Molecular networks in atopic mothers impact the risk of infant atopy.

BACKGROUND: The prevalence of atopic diseases has increased with atopic dermatitis (AD) as the earliest manifestation. We assessed if molecular risk factors in atopic mothers influence their infants' susceptibility to an atopic disease. METHODS: Pregnant women and their infants with (n = 174, high-risk) or without (n = 126, low-risk) parental atopy were enrolled in a prospective birth cohort. Global differentially methylated regions (DMRs) were determined in atopic (n = 92) and non-atopic (n = 82) mothers. Principal component analysis was used to predict atopy risk in children dependent on maternal atopy. Genome-wide transcriptomic analyses were performed in paired atopic (n = 20) and non-atopic (n = 15) mothers and cord blood. Integrative genomic analyses were conducted to define methylation-gene expression relationships. RESULTS: Atopic dermatitis was more prevalent in high-risk compared to low-risk children by age 2. Differential methylation analyses identified 165 DMRs distinguishing atopic from non-atopic mothers. Inclusion of DMRs in addition to maternal atopy significantly increased the odds ratio to develop AD in children from 2.56 to 4.26. In atopic compared to non-atopic mothers, 139 differentially expressed genes (DEGs) were identified significantly enriched of genes within the interferon signaling pathway. Expression quantitative trait methylation analyses dependent on maternal atopy identified 29 DEGs controlled by 136 trans-acting methylation marks, some located near transcription factors. Differential expression for the same nine genes, including MX1 and IFI6 within the interferon pathway, was identified in atopic and non-atopic mothers and high-risk and low-risk children. CONCLUSION: These data suggest that in utero epigenetic and transcriptomic mechanisms predominantly involving the interferon pathway may impact and predict the development of infant atopy.

Creating an ignorance-base: Exploring known unknowns in the scientific literature.

BACKGROUND: Scientific discovery progresses by exploring new and uncharted territory. More specifically, it advances by a process of transforming unknown unknowns first into known unknowns, and then into knowns. Over the last few decades, researchers have developed many knowledge bases to capture and connect the knowns, which has enabled topic exploration and contextualization of experimental results. But recognizing the unknowns is also critical for finding the most pertinent questions and their answers. Prior work on known unknowns has sought to understand them, annotate them, and automate their identification. However, no knowledge-bases yet exist to capture these unknowns, and little work has focused on how scientists might use them to trace a given topic or experimental result in search of open questions and new avenues for exploration. We show here that a knowledge base of unknowns can be connected to ontologically grounded biomedical knowledge to accelerate research in the field of prenatal nutrition. RESULTS: We present the first ignorance-base, a knowledge-base created by combining classifiers to recognize ignorance statements (statements of missing or incomplete knowledge that imply a goal for knowledge) and biomedical concepts over the prenatal nutrition literature. This knowledge-base places biomedical concepts mentioned in the literature in context with the ignorance statements authors have made about them. Using our system, researchers interested in the topic of vitamin D and prenatal health were able to uncover three new avenues for exploration (immune system, respiratory system, and brain development) by searching for concepts enriched in ignorance statements. These were buried among the many standard enriched concepts. Additionally, we used the ignorance-base to enrich concepts connected to a gene list associated with vitamin D and spontaneous preterm birth and found an emerging topic of study (brain development) in an implied field (neuroscience). The researchers could look to the field of neuroscience for potential answers to the ignorance statements. CONCLUSION: Our goal is to help students, researchers, funders, and publishers better understand the state of our collective scientific ignorance (known unknowns) in order to help accelerate research through the continued illumination of and focus on the known unknowns and their respective goals for scientific knowledge.

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).

CHD4 is essential for transcriptional repression and lineage progression in B lymphopoiesis.

Cell lineage specification is a tightly regulated process that is dependent on appropriate expression of lineage and developmental stage-specific transcriptional programs. Here, we show that Chromodomain Helicase DNA-binding protein 4 (CHD4), a major ATPase/helicase subunit of Nucleosome Remodeling and Deacetylase Complexes (NuRD) in lymphocytes, is essential for specification of the early B cell lineage transcriptional program. In the absence of CHD4 in B cell progenitors in vivo, development of these cells is arrested at an early pro-B-like stage that is unresponsive to IL-7 receptor signaling and unable to efficiently complete V(D)J rearrangements at Igh loci. Our studies confirm that chromatin accessibility and transcription of thousands of gene loci are controlled dynamically by CHD4 during early B cell development. Strikingly, CHD4-deficient pro-B cells express transcripts of many non-B cell lineage genes, including genes that are characteristic of other hematopoietic lineages, neuronal cells, and the CNS, lung, pancreas, and other cell types. We conclude that CHD4 inhibits inappropriate transcription in pro-B cells. Together, our data demonstrate the importance of CHD4 in establishing and maintaining an appropriate transcriptome in early B lymphopoiesis via chromatin accessibility.

The Kinome of Human Alveolar Type II and Basal Cells, and Its Reprogramming in Lung Cancer.

The discovery of mutant tyrosine kinases as oncogenic drivers of lung adenocarcinomas has changed the basic understanding of lung cancer development and therapy. Yet, expressed kinases (kinome) in lung cancer progenitor cells, as well as whether kinase expression and the overall kinome changes or is reprogrammed upon transformation, is incompletely understood. We hypothesized that the kinome differs between lung cancer progenitor cells, alveolar type II cells (ATII), and basal cells (BC) and that their respective kinomes undergo distinct lineage-specific reprogramming to adenocarcinomas and squamous cell carcinomas upon transformation. We performed RNA sequencing on freshly isolated human ATII, BC, and lung cancer cell lines to define the kinome in nontransformed cells and transformed cells. Our studies identified a unique kinome for ATII and BC and changes in their kinome upon transformation to their respective carcinomas.

Identification of shared TCR sequences from T cells in human breast cancer using emulsion RT-PCR.

Infiltration of T cells in breast tumors correlates with improved survival of patients with breast cancer, despite relatively few mutations in these tumors. To determine if T-cell specificity can be harnessed to augment immunotherapies of breast cancer, we sought to identify the alpha-beta paired T-cell receptors (TCRs) of tumor-infiltrating lymphocytes shared between multiple patients. Because TCRs function as heterodimeric proteins, we used an emulsion-based RT-PCR assay to link and amplify TCR pairs. Using this assay on engineered T-cell hybridomas, we observed ∼85% accurate pairing fidelity, although TCR recovery frequency varied. When we applied this technique to patient samples, we found that for any given TCR pair, the dominant alpha- or beta-binding partner comprised ∼90% of the total binding partners. Analysis of TCR sequences from primary tumors showed about fourfold more overlap in tumor-involved relative to tumor-free sentinel lymph nodes. Additionally, comparison of sequences from both tumors of a patient with bilateral breast cancer showed 10% overlap. Finally, we identified a panel of unique TCRs shared between patients' tumors and peripheral blood that were not found in the peripheral blood of controls. These TCRs encoded a range of V, J, and complementarity determining region 3 (CDR3) sequences on the alpha-chain, and displayed restricted V-beta use. The nucleotides encoding these shared TCR CDR3s varied, suggesting immune selection of this response. Harnessing these T cells may provide practical strategies to improve the shared antigen-specific response to breast cancer.

The transcription factors GATA2 and microphthalmia-associated transcription factor regulate Hdc gene expression in mast cells and are required for IgE/mast cell-mediated anaphylaxis.

BACKGROUND: Histamine is a critical mediator of IgE/mast cell-mediated anaphylaxis. Histamine is synthesized by decarboxylating the amino acid histidine, a reaction catalyzed by the histidine decarboxylase (Hdc) gene-encoded enzyme HDC. However, regulation of the Hdc gene in mast cells is poorly understood. OBJECTIVE: We sought to investigate the in vivo regulation of IgE/mast cell-mediated anaphylaxis by the transcription factors GATA2 and microphthalmia-associated transcription factor (MITF) and the mechanisms by which GATA2 and MITF regulate Hdc gene expression in mouse and human mast cells. METHODS: Mice deficient in the transcription factors Gata2, aryl hydrocarbon receptor (Ahr), aryl hydrocarbon receptor repressor (Ahrr), or basic helix-loop-helix family member E40 (Bhlhe40) were assessed for anaphylactic reactions. Chromatin immunoprecipitation sequencing analysis identified putative Hdc enhancers. Luciferase reporter transcription assay confirmed enhancer activities of putative enhancers in the Hdc gene. The short hairpin RNA knockdown approach was used to determine the role of MITF in regulating mouse and human HDC gene expression. RESULTS: Connective tissue mast cell-specific Gata2-deficient mice did not have IgE/mast cell-mediated anaphylaxis. GATA2 induced the expression of Mitf, Ahr, Ahrr, and Bhlhe40 in mast cells. MITF, but not AHR, AHRR, or BHLHE40, was required for anaphylaxis. MITF bound to an enhancer located 8.8 kb upstream of the transcription start site of the Hdc gene and directed enhancer activity. MITF overexpression largely restored Hdc gene expression in the Gata2-deficient mast cells. In the human mast cell line LAD2, MITF was required for the HDC gene expression and histamine synthesis. CONCLUSION: The transcription factors GATA2 and MITF regulate Hdc gene expression in mast cells and are required for IgE/mast cell-mediated anaphylaxis.

Deletion of c-FLIP from CD11bhi Macrophages Prevents Development of Bleomycin-induced Lung Fibrosis.

Idiopathic pulmonary fibrosis is a progressive lung disease with complex pathophysiology and fatal prognosis. Macrophages (MΦ) contribute to the development of lung fibrosis; however, the underlying mechanisms and specific MΦ subsets involved remain unclear. During lung injury, two subsets of lung MΦ coexist: Siglec-Fhi resident alveolar MΦ and a mixed population of CD11bhi MΦ that primarily mature from immigrating monocytes. Using a novel inducible transgenic system driven by a fragment of the human CD68 promoter, we targeted deletion of the antiapoptotic protein cellular FADD-like IL-1ß-converting enzyme-inhibitory protein (c-FLIP) to CD11bhi MΦ. Upon loss of c-FLIP, CD11bhi MΦ became susceptible to cell death. Using this system, we were able to show that eliminating CD11bhi MΦ present 7-14 days after bleomycin injury was sufficient to protect mice from fibrosis. RNA-seq analysis of lung MΦ present during this time showed that CD11bhi MΦ, but not Siglec-Fhi MΦ, expressed high levels of profibrotic chemokines and growth factors. Human MΦ from patients with idiopathic pulmonary fibrosis expressed many of the same profibrotic chemokines identified in murine CD11bhi MΦ. Elimination of monocyte-derived MΦ may help in the treatment of fibrosis. We identify c-FLIP and the associated extrinsic cell death program as a potential pathway through which these profibrotic MΦ may be pharmacologically targeted.

Systems biology of facial development: contributions of ectoderm and mesenchyme.

The rapid increase in gene-centric biological knowledge coupled with analytic approaches for genomewide data integration provides an opportunity to develop systems-level understanding of facial development. Experimental analyses have demonstrated the importance of signaling between the surface ectoderm and the underlying mesenchyme are coordinating facial patterning. However, current transcriptome data from the developing vertebrate face is dominated by the mesenchymal component, and the contributions of the ectoderm are not easily identified. We have generated transcriptome datasets from critical periods of mouse face formation that enable gene expression to be analyzed with respect to time, prominence, and tissue layer. Notably, by separating the ectoderm and mesenchyme we considerably improved the sensitivity compared to data obtained from whole prominences, with more genes detected over a wider dynamic range. From these data we generated a detailed description of ectoderm-specific developmental programs, including pan-ectodermal programs, prominence- specific programs and their temporal dynamics. The genes and pathways represented in these programs provide mechanistic insights into several aspects of ectodermal development. We also used these data to identify co-expression modules specific to facial development. We then used 14 co-expression modules enriched for genes involved in orofacial clefts to make specific mechanistic predictions about genes involved in tongue specification, in nasal process patterning and in jaw development. Our multidimensional gene expression dataset is a unique resource for systems analysis of the developing face; our co-expression modules are a resource for predicting functions of poorly annotated genes, or for predicting roles for genes that have yet to be studied in the context of facial development; and our analytic approaches provide a paradigm for analysis of other complex developmental programs.

Molecular Profile of Tumor-Specific CD8+ T Cell Hypofunction in a Transplantable Murine Cancer Model.

Mechanisms of self-tolerance often result in CD8(+) tumor-infiltrating lymphocytes (TIL) with a hypofunctional phenotype incapable of tumor clearance. Using a transplantable colon carcinoma model, we found that CD8(+) T cells became tolerized in <24 h in an established tumor environment. To define the collective impact of pathways suppressing TIL function, we compared genome-wide mRNA expression of tumor-specific CD8(+) T cells from the tumor and periphery. Notably, gene expression induced during TIL hypofunction more closely resembled self-tolerance than viral exhaustion. Differential gene expression was refined to identify a core set of genes that defined hypofunctional TIL; these data comprise the first molecular profile of tumor-specific TIL that are naturally responding and represent a polyclonal repertoire. The molecular profile of TIL was further dissected to determine the extent of overlap and distinction between pathways that collectively restrict T cell functions. As suggested by the molecular profile of TIL, protein expression of inhibitory receptor LAG-3 was differentially regulated throughout prolonged late-G1/early-S phase of the cell cycle. Our data may accelerate efficient identification of combination therapies to boost anti-tumor function of TIL specifically against tumor cells.

Regulation of toll-like receptor signaling by the SF3a mRNA splicing complex.

The innate immune response plays a key role in fighting infection by activating inflammation and stimulating the adaptive immune response. However, chronic activation of innate immunity can contribute to the pathogenesis of many diseases with an inflammatory component. Thus, various negatively acting factors turn off innate immunity subsequent to its activation to ensure that inflammation is self-limiting and to prevent inflammatory disease. These negatively acting pathways include the production of inhibitory acting alternate proteins encoded by alternative mRNA splice forms of genes in Toll-like receptor (TLR) signaling pathways. We previously found that the SF3a mRNA splicing complex was required for a robust innate immune response; SF3a acts to promote inflammation in part by inhibiting the production of a negatively acting splice form of the TLR signaling adaptor MyD88. Here we inhibit SF3a1 using RNAi and subsequently perform an RNAseq study to identify the full complement of genes and splicing events regulated by SF3a in murine macrophages. Surprisingly, in macrophages, SF3a has significant preference for mRNA splicing events within innate immune signaling pathways compared with other biological pathways, thereby affecting the splicing of specific genes in the TLR signaling pathway to modulate the innate immune response.

Three Unique Interstitial Macrophages in the Murine Lung at Steady State.

The current paradigm in macrophage biology is that some tissues mainly contain macrophages from embryonic origin, such as microglia in the brain, whereas other tissues contain postnatal-derived macrophages, such as the gut. However, in the lung and in other organs, such as the skin, there are both embryonic and postnatal-derived macrophages. In this study, we demonstrate in the steady-state lung that the mononuclear phagocyte system is comprised of three newly identified interstitial macrophages (IMs), alveolar macrophages, dendritic cells, and few extravascular monocytes. We focused on similarities and differences between the three IM subtypes, specifically, their phenotype, location, transcriptional signature, phagocytic capacity, turnover, and lack of survival dependency on fractalkine receptor, CX3CR1. Pulmonary IMs were located in the bronchial interstitium but not the alveolar interstitium. At the transcriptional level, all three IMs displayed a macrophage signature and phenotype. All IMs expressed MER proto-oncogene, tyrosine kinase, CD64, CD11b, and CX3CR1, and were further distinguished by differences in cell surface protein expression of CD206, Lyve-1, CD11c, CCR2, and MHC class II, along with the absence of Ly6C, Ly6G, and Siglec F. Most intriguingly, in addition to the lung, similar phenotypic populations of IMs were observed in other nonlymphoid organs, perhaps highlighting conserved functions throughout the body. These findings promote future research to track four distinct pulmonary macrophages and decipher the division of labor that exists between them.

Cell Origin Dictates Programming of Resident versus Recruited Macrophages during Acute Lung Injury.

Two populations of alveolar macrophages (AMs) coexist in the inflamed lung: resident AMs that arise during embryogenesis, and recruited AMs that originate postnatally from circulating monocytes. The objective of this study was to determine whether origin or environment dictates the transcriptional, metabolic, and functional programming of these two ontologically distinct populations over the time course of acute inflammation. RNA sequencing demonstrated marked transcriptional differences between resident and recruited AMs affecting three main areas: proliferation, inflammatory signaling, and metabolism. Functional assays and metabolomic studies confirmed these differences and demonstrated that resident AMs proliferate locally and are governed by increased tricarboxylic acid cycle and amino acid metabolism. Conversely, recruited AMs produce inflammatory cytokines in association with increased glycolytic and arginine metabolism. Collectively, the data show that even though they coexist in the same environment, inflammatory macrophage subsets have distinct immunometabolic programs and perform specialized functions during inflammation that are associated with their cellular origin.

Transcriptome analysis highlights the conserved difference between embryonic and postnatal-derived alveolar macrophages.

Alveolar macrophages (AMs) reside on the luminal surfaces of the airways and alveoli where they maintain host defense and promote alveolar homeostasis by ingesting inhaled particulates and regulating inflammatory responses. Recent studies have demonstrated that AMs populate the lungs during embryogenesis and self-renew throughout life with minimal replacement by circulating monocytes, except under extreme conditions of depletion or radiation injury. Here we demonstrate that on a global scale, environment appears to dictate AM development and function. Indeed, transcriptome analysis of embryonic host-derived and postnatal donor-derived AMs coexisting within the same mouse demonstrated >98% correlation and overall functional analyses were similar. However, we also identified several genes whose expression was dictated by origin rather than environment. The most differentially expressed gene not altered by environment was Marco, a gene recently demonstrated to have enhancer activity in embryonic-derived but not postnatal-derived tissue macrophages. Overall, we show that under homeostatic conditions, the environment largely dictates the programming and function of AMs, whereas the expression of a small number of genes remains linked to the origin of the cell.

Unexpected effects of different genetic backgrounds on identification of genomic rearrangements via whole-genome next generation sequencing.

BACKGROUND: Whole genome next generation sequencing (NGS) is increasingly employed to detect genomic rearrangements in cancer genomes, especially in lymphoid malignancies. We recently established a unique mouse model by specifically deleting a key non-homologous end-joining DNA repair gene, Xrcc4, and a cell cycle checkpoint gene, Trp53, in germinal center B cells. This mouse model spontaneously develops mature B cell lymphomas (termed G1XP lymphomas). RESULTS: Here, we attempt to employ whole genome NGS to identify novel structural rearrangements, in particular inter-chromosomal translocations (CTXs), in these G1XP lymphomas. We sequenced six lymphoma samples, aligned our NGS data with mouse reference genome (in C57BL/6J (B6) background) and identified CTXs using CREST algorithm. Surprisingly, we detected widespread CTXs in both lymphomas and wildtype control samples, majority of which were false positive and attributable to different genetic backgrounds. In addition, we validated our NGS pipeline by sequencing multiple control samples from distinct tissues of different genetic backgrounds of mouse (B6 vs non-B6). Lastly, our studies showed that widespread false positive CTXs can be generated by simply aligning sequences from different genetic backgrounds of mouse. CONCLUSIONS: We conclude that mapping and alignment with reference genome might not be a preferred method for analyzing whole-genome NGS data obtained from a genetic background different from reference genome. Given the complex genetic background of different mouse strains or the heterogeneity of cancer genomes in human patients, in order to minimize such systematic artifacts and uncover novel CTXs, a preferred method might be de novo assembly of personalized normal control genome and cancer cell genome, instead of mapping and aligning NGS data to mouse or human reference genome. Thus, our studies have critical impact on the manner of data analysis for cancer genomics.

Bioinformatics-driven discovery of rational combination for overcoming EGFR-mutant lung cancer resistance to EGFR therapy.

MOTIVATION: Non-small-cell lung cancer (NSCLC) is the leading cause of cancer death in the United States. Targeted tyrosine kinase inhibitors (TKIs) directed against the epidermal growth factor receptor (EGFR) have been widely and successfully used in treating NSCLC patients with activating EGFR mutations. Unfortunately, the duration of response is short-lived, and all patients eventually relapse by acquiring resistance mechanisms. RESULT: We performed an integrative systems biology approach to determine essential kinases that drive EGFR-TKI resistance in cancer cell lines. We used a series of bioinformatics methods to analyze and integrate the functional genetics screen and RNA-seq data to identify a set of kinases that are critical in survival and proliferation in these TKI-resistant lines. By connecting the essential kinases to compounds using a novel kinase connectivity map (K-Map), we identified and validated bosutinib as an effective compound that could inhibit proliferation and induce apoptosis in TKI-resistant lines. A rational combination of bosutinib and gefitinib showed additive and synergistic effects in cancer cell lines resistant to EGFR TKI alone. CONCLUSIONS: We have demonstrated a bioinformatics-driven discovery roadmap for drug repurposing and development in overcoming resistance in EGFR-mutant NSCLC, which could be generalized to other cancer types in the era of personalized medicine. AVAILABILITY AND IMPLEMENTATION: K-Map can be accessible at: http://tanlab.ucdenver.edu/kMap. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Rationale and design of the prognostic transcriptomic signature in fibrotic hypersensitivity pneumonitis (PREDICT) study.

Hypersensitivity pneumonitis is an immunologically mediated form of lung disease, resulting from inhalational exposure to a large variety of antigens. A subgroup of patients with fibrotic hypersensitivity pneumonitis (FHP) develop symptomatic, functional and radiographic disease progression. Mortality occurs primarily from respiratory failure as a result of progressive and self-sustaining lung injury that often occurs despite immunosuppression and removal of the inciting antigen. The development and validation of a prognostic transcriptomic signature for FHP (PREDICT-HP) is an observational multicentre cohort study designed to explore a transcriptomic signature from peripheral blood mononuclear cells in patients with FHP that is predictive of disease progression. This article describes the design and rationale of the PREDICT-HP study. This study will enrol ∼135 patients with FHP at approximately seven academic medical sites. Participants with a confirmed diagnosis of FHP are followed over 24 months and undergo physical examinations, self-administered questionnaires, chest computed tomography, pulmonary function tests, a 6-min walk test and blood testing for transcriptomic analyses. At each 6-month follow-up visit the study will assess the participants' clinical course and clinical events including hospitalisations and respiratory exacerbations. The PREDICT study has the potential to enhance our ability to predict disease progression and fundamentally advance our understanding of the pathobiology of FHP disease progression.

Expression of cilium-associated genes defines novel molecular subtypes of idiopathic pulmonary fibrosis.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is an untreatable lung disease with a median survival of only 3-5 years that is diagnosed using a combination of clinical, radiographic and pathologic criteria. Histologically, IPF is characterised by usual interstitial pneumonia (UIP), a fibrosing interstitial pneumonia with a pattern of heterogeneous, subpleural regions of fibrotic and remodelled lung. We hypothesised that gene expression profiles of lung tissue may identify molecular subtypes of disease that could classify subtypes of IPF/UIP that have clinical implications. METHODS AND FINDINGS: We collected transcriptional profiles on lung tissue from 119 patients with IPF/UIP and 50 non-diseased controls. Differential expression of individual transcripts was identified using an analysis of covariance (ANCOVA) model incorporating the clinical diagnosis of each patient as well as age, gender and smoking status. Validation was performed in an independent cohort of 111 IPF/UIP and 39 non-diseased controls. Our analysis identified two subtypes of IPF/UIP based on a strong molecular signature associated with expression of genes previously associated with fibrosis (matrix metalloproteinases, osteopontin, keratins), cilium genes and genes with unknown function. We demonstrate that elevated expression of cilium genes is associated with more extensive microscopic honeycombing and higher expression of both the airway mucin gene MUC5B and the metalloproteinase MMP7, a gene recently implicated in attenuating ciliated cell differentiation during wound repair. CONCLUSIONS: Expression of cilium genes appears to identify two unique molecular phenotypes of IPF/UIP. The different molecular profiles may be relevant to therapeutic responsiveness in patients with IPF/UIP.

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.