Peripheral blood cellular dynamics of rheumatoid arthritis treatment informs about efficacy of response to disease modifying drugs.

Rheumatoid arthritis (RA) is an autoimmune disease characterized by systemic inflammation and is mediated by multiple immune cell types. In this work, we aimed to determine the relevance of changes in cell proportions in peripheral blood mononuclear cells (PBMCs) during the development of disease and following treatment. Samples from healthy blood donors, newly diagnosed RA patients, and established RA patients that had an inadequate response to MTX and were about to start tumor necrosis factor inhibitors (TNFi) treatment were collected before and after 3 months of treatment. We used in parallel a computational deconvolution approach based on RNA expression and flow cytometry to determine the relative cell-type frequencies. Cell-type frequencies from deconvolution of gene expression indicate that monocytes (both classical and non-classical) and CD4+ cells (Th1 and Th2) were increased in RA patients compared to controls, while NK cells and B cells (naïve and mature) were significantly decreased in RA patients. Treatment with MTX caused a decrease in B cells (memory and plasma cell), and a decrease in CD4 Th cells (Th1 and Th17), while treatment with TNFi resulted in a significant increase in the population of B cells. Characterization of the RNA expression patterns found that most of the differentially expressed genes in RA subjects after treatment can be explained by changes in cell frequencies (98% and 74% respectively for MTX and TNFi).

Molecular signature of methotrexate response among rheumatoid arthritis patients.

Background: Methotrexate (MTX) is the first line treatment for rheumatoid arthritis (RA), but failure of satisfying treatment response occurs in a significant proportion of patients. Here we present a longitudinal multi-omics study aimed at detecting molecular and cellular processes in peripheral blood associated with a successful methotrexate treatment of rheumatoid arthritis. Methods: Eighty newly diagnosed patients with RA underwent clinical assessment and donated blood before initiation of MTX, and 3 months into treatment. Flow cytometry was used to describe cell types and presence of activation markers in peripheral blood, the expression of 51 proteins was measured in serum or plasma, and RNA sequencing was performed in peripheral blood mononuclear cells (PBMC). Response to treatment after 3 months was determined using the EULAR response criteria. We assessed the changes in biological phenotypes during treatment, and whether these changes differed between responders and non-responders with regression analysis. By using measurements from baseline, we also tried to find biomarkers of future MTX response or, alternatively, to predict MTX response. Results: Among the MTX responders, (Good or Moderate according to EULAR treatment response classification, n = 60, 75%), we observed changes in 29 partly overlapping cell types proportions, levels of 13 proteins and expression of 38 genes during treatment. These changes were in most cases suppressions that were stronger among responders compared to non-responders. Within responders to treatment, we observed a suppression of FOXP3 gene expression, reduction of immunoglobulin gene expression and suppression of genes involved in cell proliferation. The proportion of many HLA-DR expressing T-cell populations were suppressed in all patients irrespective of clinical response, and the proportion of many IL21R+ T-cells were reduced exclusively in non-responders. Using only the baseline measurements we could not detect any biomarkers or prediction models that could predict response to MTX. Conclusion: We conclude that a deep molecular and cellular phenotyping of peripheral blood cells in RA patients treated with methotrexate can reveal previously not recognized differences between responders and non-responders during 3 months of treatment with MTX. This may contribute to the understanding of MTX mode of action and explain non-responsiveness to MTX therapy.

Molecular and Functional Characterization of Human Intestinal Organoids and Monolayers for Modeling Epithelial Barrier.

BACKGROUND: Patient-derived organoid (PDO) models offer potential to transform drug discovery for inflammatory bowel disease (IBD) but are limited by inconsistencies with differentiation and functional characterization. We profiled molecular and cellular features across a range of intestinal organoid models and examined differentiation and establishment of a functional epithelial barrier. METHODS: Patient-derived organoids or monolayers were generated from control or IBD patient-derived colon or ileum and were molecularly or functionally profiled. Biological or technical replicates were examined for transcriptional responses under conditions of expansion or differentiation. Cell-type composition was determined by deconvolution of cell-associated gene signatures and histological features. Differentiated control or IBD-derived monolayers were examined for establishment of transepithelial electrical resistance (TEER), loss of barrier integrity in response to a cocktail of interferon (IFN)-γ and tumor necrosis factor (TNF)-α, and prevention of cytokine-induced barrier disruption by the JAK inhibitor, tofacitinib. RESULTS: In response to differentiation media, intestinal organoids and monolayers displayed gene expression patterns consistent with maturation of epithelial cell types found in the human gut. Upon differentiation, both colon- and ileum-derived monolayers formed functional barriers, with sustained TEER. Barrier integrity was compromised by inflammatory cytokines IFN-γ and TNF-α, and damage was inhibited in a dose-dependent manner by tofacitinib. CONCLUSIONS: We describe the generation and characterization of human colonic or ileal organoid models capable of functional differentiation to mature epithelial cell types. In monolayer culture, these cells formed a robust epithelial barrier with sustained TEER and responses to pharmacological modulation. Our findings demonstrate that control and IBD patient-derived organoids possess consistent transcriptional and functional profiles that can enable development of epithelial-targeted therapies.

Ageing and interferon gamma response drive the phenotype of neutrophils in the inflamed joint.

OBJECTIVE: Neutrophils are typically the most abundant leucocyte in arthritic synovial fluid. We sought to understand changes that occur in neutrophils as they migrate from blood to joint. METHODS: We performed RNA sequencing of neutrophils from healthy human blood, arthritic blood and arthritic synovial fluid, comparing transcriptional signatures with those from murine K/BxN serum transfer arthritis. We employed mass cytometry to quantify protein expression and sought to reproduce the synovial fluid phenotype ex vivo in cultured healthy blood neutrophils. RESULTS: Blood neutrophils from healthy donors and patients with active arthritis showed largely similar transcriptional signatures. By contrast, synovial fluid neutrophils exhibited more than 1600 differentially expressed genes. Gene signatures identified a prominent response to interferon gamma (IFN-γ), as well as to tumour necrosis factor, interleukin-6 and hypoxia, in both humans and mice. Mass cytometry confirmed that healthy and arthritic donor blood neutrophils are largely indistinguishable but revealed a range of neutrophil phenotypes in synovial fluid defined by downregulation of CXCR1 and upregulation of FcγRI, HLA-DR, PD-L1, ICAM-1 and CXCR4. Reproduction of key elements of this signature in cultured blood neutrophils required both IFN-γ and prolonged culture. CONCLUSIONS: Circulating neutrophils from patients with arthritis resemble those from healthy controls, but joint fluid cells exhibit a network of changes, conserved across species, that implicate IFN-γ response and ageing as complementary drivers of the synovial fluid neutrophil phenotype.

Early prediction of clinical response to anti-TNF treatment using multi-omics and machine learning in rheumatoid arthritis.

OBJECTIVES: Advances in immunotherapy by blocking TNF have remarkably improved treatment outcomes for Rheumatoid arthritis (RA) patients. Although treatment specifically targets TNF, the downstream mechanisms of immune suppression are not completely understood. The aim of this study was to detect biomarkers and expression signatures of treatment response to TNF inhibition. METHODS: Peripheral blood mononuclear cells (PBMCs) from 39 female patients were collected before anti-TNF treatment initiation (day 0) and after 3 months. The study cohort included patients previously treated with MTX who failed to respond adequately. Response to treatment was defined based on the EULAR criteria and classified 23 patients as responders and 16 as non-responders. We investigated differences in gene expression in PBMCs, the proportion of cell types and cell phenotypes in peripheral blood using flow cytometry and the level of proteins in plasma. Finally, we used machine learning models to predict non-response to anti-TNF treatment. RESULTS: The gene expression analysis in baseline samples revealed notably higher expression of the gene EPPK1 in future responders. We detected the suppression of genes and proteins following treatment, including suppressed expression of the T cell inhibitor gene CHI3L1 and its protein YKL-40. The gene expression results were replicated in an independent cohort. Finally, machine learning models mainly based on transcriptomic data showed high predictive utility in classifying non-response to anti-TNF treatment in RA. CONCLUSIONS: Our integrative multi-omics analyses identified new biomarkers for the prediction of response, found pathways influenced by treatment and suggested new predictive models of anti-TNF treatment in RA patients.

What Is the Persistence to Methotrexate in Rheumatoid Arthritis, and Does Machine Learning Outperform Hypothesis-Based Approaches to Its Prediction?

OBJECTIVE: The objectives of this study were to assess the 1-year persistence to methotrexate (MTX) initiated as the first ever conventional synthetic disease-modifying antirheumatic drug in new-onset rheumatoid arthritis (RA) and to investigate the marginal gains and robustness of the results by increasing the number and nature of covariates and by using data-driven, instead of hypothesis-based, methods to predict this persistence. METHODS: Through the Swedish Rheumatology Quality Register, linked to other data sources, we identified a cohort of 5475 patients with new-onset RA in 2006-2016 who were starting MTX monotherapy as their first disease-modifying antirheumatic drug. Data on phenotype at diagnosis and demographics were combined with increasingly detailed data on medical disease history and medication use in four increasingly complex data sets (48-4162 covariates). We performed manual model building using logistic regression. We also performed five different machine learning (ML) methods and combined the ML results into an ensemble model. We calculated the area under the receiver operating characteristic curve (AUROC) and made calibration plots. We trained on 90% of the data, and tested the models on a holdout data set. RESULTS: Of the 5475 patients, 3834 (70%) remained on MTX monotherapy 1 year after treatment start. Clinical RA disease activity and baseline characteristics were most strongly associated with the outcome. The best manual model had an AUROC of 0.66 (95% confidence interval [CI] 0.60-0.71). For the ML methods, Lasso regression performed best (AUROC = 0.67; 95% CI 0.62-0.71). CONCLUSION: Approximately two thirds of patients with early RA who start MTX remain on this therapy 1 year later. Predicting this persistence remains a challenge, whether using hypothesis-based or ML models, and may yet require additional types of data.

Distinct clinical phenotypes for Crohn's disease derived from patient surveys.

BACKGROUND: Defining clinical phenotypes provides opportunities for new diagnostics and may provide insights into early intervention and disease prevention. There is increasing evidence that patient-derived health data may contain information that complements traditional methods of clinical phenotyping. The utility of these data for defining meaningful phenotypic groups is of great interest because social media and online resources make it possible to query large cohorts of patients with health conditions. METHODS: We evaluated the degree to which patient-reported categorical data is useful for discovering subclinical phenotypes and evaluated its utility for discovering new measures of disease severity, treatment response and genetic architecture. Specifically, we examined the responses of 1961 patients with inflammatory bowel disease to questionnaires in search of sub-phenotypes. We applied machine learning methods to identify novel subtypes of Crohn's disease and studied their associations with drug responses. RESULTS: Using the patients' self-reported information, we identified two subpopulations of Crohn's disease; these subpopulations differ in disease severity, associations with smoking, and genetic transmission patterns. We also identified distinct features of drug response for the two Crohn's disease subtypes. These subtypes show a trend towards differential genotype signatures. CONCLUSION: Our findings suggest that patient-defined data can have unplanned utility for defining disease subtypes and may be useful for guiding treatment approaches.

Myod1 and GR coordinate myofiber-specific transcriptional enhancers.

Skeletal muscle is a dynamic tissue the size of which can be remodeled through the concerted actions of various cues. Here, we investigated the skeletal muscle transcriptional program and identified key tissue-specific regulatory genetic elements. Our results show that Myod1 is bound to numerous skeletal muscle enhancers in collaboration with the glucocorticoid receptor (GR) to control gene expression. Remarkably, transcriptional activation controlled by these factors occurs through direct contacts with the promoter region of target genes, via the CpG-bound transcription factor Nrf1, and the formation of Ctcf-anchored chromatin loops, in a myofiber-specific manner. Moreover, we demonstrate that GR negatively controls muscle mass and strength in mice by down-regulating anabolic pathways. Taken together, our data establish Myod1, GR and Nrf1 as key players of muscle-specific enhancer-promoter communication that orchestrate myofiber size regulation.

Single-cell analyses of Crohn's disease tissues reveal intestinal intraepithelial T cells heterogeneity and altered subset distributions.

Crohn's disease (CD) is a chronic transmural inflammation of intestinal segments caused by dysregulated interaction between microbiome and gut immune system. Here, we profile, via multiple single-cell technologies, T cells purified from the intestinal epithelium and lamina propria (LP) from terminal ileum resections of adult severe CD cases. We find that intraepithelial lymphocytes (IEL) contain several unique T cell subsets, including NKp30+γδT cells expressing RORγt and producing IL-26 upon NKp30 engagement. Further analyses comparing tissues from non-inflamed and inflamed regions of patients with CD versus healthy controls show increased activated TH17 but decreased CD8+T, γδT, TFH and Treg cells in inflamed tissues. Similar analyses of LP find increased CD8+, as well as reduced CD4+T cells with an elevated TH17 over Treg/TFH ratio. Our analyses of CD tissues thus suggest a potential link, pending additional validations, between transmural inflammation, reduced IEL γδT cells and altered spatial distribution of IEL and LP T cell subsets.

Publisher Correction: Subsets of ILC3-ILC1-like cells generate a diversity spectrum of innate lymphoid cells in human mucosal tissues.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Subsets of ILC3-ILC1-like cells generate a diversity spectrum of innate lymphoid cells in human mucosal tissues.

Innate lymphoid cells (ILCs) are tissue-resident lymphocytes categorized on the basis of their core regulatory programs and the expression of signature cytokines. Human ILC3s that produce the cytokine interleukin-22 convert into ILC1-like cells that produce interferon-γ in vitro, but whether this conversion occurs in vivo remains unclear. In the present study we found that ILC3s and ILC1s in human tonsils represented the ends of a spectrum that included additional discrete subsets. RNA velocity analysis identified an intermediate ILC3-ILC1 cluster, which had strong directionality toward ILC1s. In humanized mice, the acquisition of ILC1 features by ILC3s showed tissue dependency. Chromatin studies indicated that the transcription factors Aiolos and T-bet cooperated to repress regulatory elements active in ILC3s. A transitional ILC3-ILC1 population was also detected in the human intestine. We conclude that ILC3s undergo conversion into ILC1-like cells in human tissues in vivo, and that tissue factors and Aiolos were required for this process.

A BMP/activin A chimera is superior to native BMPs and induces bone repair in nonhuman primates when delivered in a composite matrix.

Bone morphogenetic protein (BMP)/carriers approved for orthopedic procedures achieve efficacy superior or equivalent to autograft bone. However, required supraphysiological BMP concentrations have been associated with potential local and systemic adverse events. Suboptimal BMP/receptor binding and rapid BMP release from approved carriers may contribute to these outcomes. To address these issues and improve efficacy, we engineered chimeras with increased receptor binding by substituting BMP-6 and activin A receptor binding domains into BMP-2 and optimized a carrier for chimera retention and tissue ingrowth. BV-265, a BMP-2/BMP-6/activin A chimera, demonstrated increased binding affinity to BMP receptors, including activin-like kinase-2 (ALK2) critical for bone formation in people. BV-265 increased BMP intracellular signaling, osteogenic activity, and expression of bone-related genes in murine and human cells to a greater extent than BMP-2 and was not inhibited by BMP antagonist noggin or gremlin. BV-265 induced larger ectopic bone nodules in rats compared to BMP-2 and was superior to BMP-2, BMP-2/6, and other chimeras in nonhuman primate bone repair models. A composite matrix (CM) containing calcium-deficient hydroxyapatite granules suspended in a macroporous, fenestrated, polymer mesh-reinforced recombinant human type I collagen matrix demonstrated improved BV-265 retention, minimal inflammation, and enhanced handling. BV-265/CM was efficacious in nonhuman primate bone repair models at concentrations ranging from 1/10 to 1/30 of the BMP-2/absorbable collagen sponge (ACS) concentration approved for clinical use. Initial toxicology studies were negative. These results support evaluations of BV-265/CM as an alternative to BMP-2/ACS in clinical trials for orthopedic conditions requiring augmented healing.

High-dimensional immune phenotyping and transcriptional analyses reveal robust recovery of viable human immune and epithelial cells from frozen gastrointestinal tissue.

Simultaneous analyses of peripheral and mucosal immune compartments can yield insight into the pathogenesis of mucosal-associated diseases. Although methods to preserve peripheral immune cells are well established, studies involving mucosal immune cells have been hampered by lack of simple storage techniques. We provide a cryopreservation protocol allowing for storage of gastrointestinal (GI) tissue with preservation of viability and functionality of both immune and epithelial cells. These methods will facilitate translational studies allowing for batch analysis of mucosal tissue to investigate disease pathogenesis, biomarker discovery and treatment responsiveness.

Human CCL3L1 copy number variation, gene expression, and the role of the CCL3L1-CCR5 axis in lung function.

Background: The CCL3L1-CCR5 signaling axis is important in a number of inflammatory responses, including macrophage function, and T-cell-dependent immune responses. Small molecule CCR5 antagonists exist, including the approved antiretroviral drug maraviroc, and therapeutic monoclonal antibodies are in development. Repositioning of drugs and targets into new disease areas can accelerate the availability of new therapies and substantially reduce costs. As it has been shown that drug targets with genetic evidence supporting their involvement in the disease are more likely to be successful in clinical development, using genetic association studies to identify new target repurposing opportunities could be fruitful. Here we investigate the potential of perturbation of the CCL3L1-CCR5 axis as treatment for respiratory disease. Europeans typically carry between 0 and 5 copies of CCL3L1 and this multi-allelic variation is not detected by widely used genome-wide single nucleotide polymorphism studies.  Methods: We directly measured the complex structural variation of CCL3L1 using the Paralogue Ratio Test and imputed (with validation) CCR5del32 genotypes in 5,000 individuals from UK Biobank, selected from the extremes of the lung function distribution, and analysed DNA and RNAseq data for CCL3L1 from the 1000 Genomes Project. Results: We confirmed the gene dosage effect of CCL3L1 copy number on CCL3L1 mRNA expression levels.  We found no evidence for association of CCL3L1 copy number or CCR5del32 genotype with lung function. Conclusions: These results suggest that repositioning CCR5 antagonists is unlikely to be successful for the treatment of airflow obstruction.

IRAK4 kinase activity controls Toll-like receptor-induced inflammation through the transcription factor IRF5 in primary human monocytes.

Interleukin-1 receptor-associated kinase 4 (IRAK4) plays a critical role in innate immune signaling by Toll-like receptors (TLRs), and loss of IRAK4 activity in mice and humans increases susceptibility to bacterial infections and causes defects in TLR and IL1 ligand sensing. However, the mechanism by which IRAK4 activity regulates the production of downstream inflammatory cytokines is unclear. Using transcriptomic and biochemical analyses of human monocytes treated with a highly potent and selective inhibitor of IRAK4, we show that IRAK4 kinase activity controls the activation of interferon regulatory factor 5 (IRF5), a transcription factor implicated in the pathogenesis of multiple autoimmune diseases. Following TLR7/8 stimulation by its agonist R848, chemical inhibition of IRAK4 abolished IRF5 translocation to the nucleus and thus prevented IRF5 binding to and activation of the promoters of inflammatory cytokines in human monocytes. We also found that IKKß, an upstream IRF5 activator, is phosphorylated in response to the agonist-induced TLR signaling. Of note, IRAK4 inhibition blocked IKKß phosphorylation but did not block the nuclear translocation of NFκB, which was surprising, given the canonical role of IKKß in phosphorylating IκB to allow NFκB activation. Moreover, pharmacological inhibition of either IKKß or the serine/threonine protein kinase TAK1 in monocytes blocked TLR-induced cytokine production and IRF5 translocation to the nucleus, but not nuclear translocation of NFκB. Taken together, our data suggest a mechanism by which IRAK4 activity regulates TAK1 and IKKß activation, leading to the nuclear translocation of IRF5 and induction of inflammatory cytokines in human monocytes.

Defining the inflammatory signature of human lung explant tissue in the presence and absence of glucocorticoid.

BACKGROUND:  Airway inflammation is a feature of many respiratory diseases and there is a need for newer, more effective anti-inflammatory compounds. The aim of this study was to develop an ex vivo human lung explant model which can be used to help study the mechanisms underlying inflammatory responses and which can provide a tool to aid drug discovery for inflammatory respiratory diseases such as asthma and COPD. METHOD:  Parenchymal lung tissue from 6 individual donors was dissected and cultured with two pro-inflammatory stimuli, lipopolysaccharide (LPS) (1 µg/ml) and interleukin-1 beta (IL-1ß) (10 ng/ml) in the presence or absence of dexamethasone (1 µM).  Inflammatory responses were assessed using Luminex analysis of tissue culture supernatants to measure levels of 21 chemokines, growth factors and cytokines. RESULTS:  A robust and reproducible inflammatory signal was detected across all donors for 12 of the analytes measured following LPS stimulation with a modest fold increase (<2-fold) in levels of CCL22, IL-4, and IL-2; increases of 2-4-fold in levels of CXCL8, VEGF and IL-6 and increases >4-fold in CCL3, CCL4, GM-CSF, IL-10, TNF-α and IL-1ß.  The inflammatory signal induced by IL-1ß stimulation was less than that observed with LPS but resulted in elevated levels of 7 analytes (CXCL8, CCL3, CCL4, GM-CSF, IL-6, IL-10 and TNF-α).  The inflammatory responses induced by both stimulations was supressed by dexamethasone for the majority of analytes. CONCLUSIONS:  These data provide proof of concept that this ex vivo human lung explant model is responsive to inflammatory signals and could be used to investigate the anti-inflammatory effects of existing and novel compounds.  In addition this model could be used to help define the mechanisms and pathways involved in development of inflammatory airway disease. ABBREVIATIONS:  COPD: Chronic Obstructive Pulmonary Disease; ICS: inhaled corticosteroids; LPS: lipopolysaccharide; IL-1ß: interleukin-1 beta; PSF: penicillin, streptomycin and fungizone.

Stepwise Distributed Open Innovation Contests for Software Development: Acceleration of Genome-Wide Association Analysis.

Background: The association of differing genotypes with disease-related phenotypic traits offers great potential to both help identify new therapeutic targets and support stratification of patients who would gain the greatest benefit from specific drug classes. Development of low-cost genotyping and sequencing has made collecting large-scale genotyping data routine in population and therapeutic intervention studies. In addition, a range of new technologies is being used to capture numerous new and complex phenotypic descriptors. As a result, genotype and phenotype datasets have grown exponentially. Genome-wide association studies associate genotypes and phenotypes using methods such as logistic regression. As existing tools for association analysis limit the efficiency by which value can be extracted from increasing volumes of data, there is a pressing need for new software tools that can accelerate association analyses on large genotype-phenotype datasets. Results: Using open innovation (OI) and contest-based crowdsourcing, the logistic regression analysis in a leading, community-standard genetics software package (PLINK 1.07) was substantially accelerated. OI allowed us to do this in <6 months by providing rapid access to highly skilled programmers with specialized, difficult-to-find skill sets. Through a crowd-based contest a combination of computational, numeric, and algorithmic approaches was identified that accelerated the logistic regression in PLINK 1.07 by 18- to 45-fold. Combining contest-derived logistic regression code with coarse-grained parallelization, multithreading, and associated changes to data initialization code further developed through distributed innovation, we achieved an end-to-end speedup of 591-fold for a data set size of 6678 subjects by 645 863 variants, compared to PLINK 1.07's logistic regression. This represents a reduction in run time from 4.8 hours to 29 seconds. Accelerated logistic regression code developed in this project has been incorporated into the PLINK2 project. Conclusions: Using iterative competition-based OI, we have developed a new, faster implementation of logistic regression for genome-wide association studies analysis. We present lessons learned and recommendations on running a successful OI process for bioinformatics.

Genetic variants affecting cross-sectional lung function in adults show little or no effect on longitudinal lung function decline.

BACKGROUND: Genome-wide association studies have identified numerous genetic regions that influence cross-sectional lung function. Longitudinal decline in lung function also includes a heritable component but the genetic determinants have yet to be defined. OBJECTIVES: We aimed to determine whether regions associated with cross-sectional lung function were also associated with longitudinal decline and to seek novel variants which influence decline. METHODS: We analysed genome-wide data from 4167 individuals from the Busselton Health Study cohort, who had undergone spirometry (12 695 observations across eight time points). A mixed model was fitted and weighted risk scores were calculated for the joint effect of 26 known regions on baseline and longitudinal changes in FEV1 and FEV1/FVC. Potential additional regions of interest were identified and followed up in two independent cohorts. RESULTS: The 26 regions previously associated with cross-sectional lung function jointly showed a strong effect on baseline lung function (p=4.44×10-16 for FEV1/FVC) but no effect on longitudinal decline (p=0.160 for FEV1/FVC). This was replicated in an independent cohort. 39 additional regions of interest (48 variants) were identified; these associations were not replicated in two further cohorts. CONCLUSIONS: Previously identified genetic variants jointly have a strong effect on cross-sectional lung function in adults but little or no effect on the rate of decline of lung function. It is possible that they influence COPD risk through lung development. Although no genetic variants have yet been associated with lung function decline at stringent genome-wide significance, longitudinal change in lung function is heritable suggesting that there is scope for future discoveries.

Lysine Methyltransferase SETD7 (SET7/9) Regulates ROS Signaling through mitochondria and NFE2L2/ARE pathway.

Reactive oxygen species (ROS) homeostasis requires stringent regulation. ROS imbalance, especially ROS accumulation, has profound implications in various disease pathogenesis. Lysine methylation of histone and non-histone proteins has been implicated in various cellular responses. The main objective of this study is to investigate the role of SET domain containing lysine methyltransferase SETD7 (SET7/9) in the regulation of ROS-mediated signaling. Here we report that inhibition of SETD7 with siRNA or a SETD7 small molecule inhibitor in both macrophages and a human bronchial epithelial cell line (Beas-2B) were able to counter NF-ĸB-induced oxidative stress and pro-inflammatory cytokine production. Meanwhile, inhibition of SETD7 elevates mitochondria antioxidant functions via negative regulation of PPARGC1A and NFE2L2. Using a co-expression system and purified proteins, we detected direct interaction between SETD7 and NFE2L2. These results indicate that lysine methylation by SETD7 is important for the fine-tuning of ROS signaling through its regulation on pro-inflammatory responses, mitochondrial function and the NFE2L2/ARE pathway. Up-regulation of multiple antioxidant genes and improved ROS clearance by inhibition of SETD7 suggests the potential benefit of targeting SETD7 in treating ROS-associated diseases.

Whole exome re-sequencing implicates CCDC38 and cilia structure and function in resistance to smoking related airflow obstruction.

Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV1) despite many years of heavy smoking. We hypothesised that these "resistant smokers" may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re-sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the "resistant smokers" and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.34 × 10(-4)) in CCDC38, a gene which has previously been reported to show association with FEV1/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke.