Search

Objectives: COVID-19 and systemic sclerosis (SSc) share multiple similarities in their clinical manifestations, alterations in immune response, and therapeutic options. These resemblances have also been identified in other immune-mediated inflammatory diseases where a common genetic component has been found. Thus, we decided to evaluate for the first time this shared genetic architecture with SSc., Methods: For this study, we retrieved genomic data from two European-ancestry cohorts: 2,597 856 individuals from The COVID-19 Host Genetics Initiative consortium, and 26 679 individuals from the largest genomic scan in SSc. We performed a cross-trait meta-analyses including >9.3 million SNPs. Finally, we conducted functional annotation to prioritize potential causal genes and performed drug repurposing analysis., Results: Our results revealed a total of 19 non-HLA pleiotropic loci, including 2 novel associations for both conditions (BMP1 and PPARG), and 12 emerging as new shared loci. Functional annotation of these regions underscored their potential regulatory role and identified potential causal genes, many of which are implicated in fibrotic and inflammatory pathways. Remarkably, we observed an antagonistic pleiotropy model of the IFN signalling between COVID-19 and SSc, including the well-known TYK2 P1104A missense variant, showing a protective effect for SSc while being a risk factor for COVID-19, along with two additional novel pleiotropic associations (IRF8 and SENP7). Finally, our findings provide new therapeutic options that could potentially benefit both conditions., Conclusion: Our study confirms the genetic resemblance between susceptibility to and severity of COVID-19 and SSc, revealing a novel common genetic contribution affecting fibrotic and immune pathways., (© The Author(s) 2025. Published by Oxford University Press on behalf of the British Society for Rheumatology.)

Objective: Nongenetic factors influence systemic sclerosis (SSc) pathogenesis, underscoring epigenetics as a relevant contributor to the disease. We aimed to unravel DNA methylation abnormalities associated with SSc through an epigenome-wide association study., Methods: We analyzed DNA methylation data from whole-blood samples in 179 patients with SSc and 241 unaffected individuals to identify differentially methylated positions (DMPs) with a false discovery rate (FDR) <0.05. These results were further integrated with RNA sequencing data from the same patients to assess their functional consequence. Additionally, we examined the impact of DNA methylation changes on transcription factors and analyzed the relationship between alterations of the methylation and gene expression profile and serum proteins levels., Results: This analysis yielded 525 DMPs enriched in immune-related pathways, with leukocyte cell-cell adhesion being the most significant (FDR = 4.91 × 10 -9 ), prioritizing integrins as they were exposed by integrating methylome and transcriptome data. Furthermore, through this integrative approach, we observed an enrichment of neutrophil-related pathways, highlighting this myeloid cell type as a relevant contributor in SSc pathogenesis. In addition, we uncovered novel profibrotic and proinflammatory mechanisms involved in the disease. Finally, the altered epigenetic and transcriptomic signature revealed an increased activity of CCAAT/enhancer-binding protein transcription factor family in SSc, which is crucial in the myeloid lineage development., Conclusion: Our findings uncover the impaired epigenetic regulation of the disease and its impact on gene expression, identifying new molecules for potential clinical applications and improving our understanding of SSc pathogenesis., (© 2024 The Author(s). Arthritis & Rheumatology published by Wiley Periodicals LLC on behalf of American College of Rheumatology.)

Objectives: Systemic sclerosis (SSc) is a heterogeneous disease, complicating its management. Its complexity and the insufficiency of clinical manifestations alone to delineate homogeneous patient groups further challenge this task. However, autoantibodies could serve as relevant markers for the pathophysiological mechanisms driving the disease. Identifying specific immunological mechanisms based on patients' serological statuses might facilitate a deeper understanding of the diversity of the disease., Methods: A cohort of 206 patients with SSc enrolled in the PRECISESADS cross-sectional study was examined. Patients were stratified based on their anti-centromere (ACA) and anti-SCL70 (SCL70) antibody statuses. Comprehensive omics analyses including transcriptomic, flow cytometric, cytokine and metabolomic data were analysed to characterise the differences between these patient groups., Results: Patients with SCL70 antibodies showed severe clinical features such as diffuse cutaneous sclerosis and pulmonary fibrosis and were biologically distinguished by unique transcriptomic profiles. They exhibit a pro-inflammatory and fibrotic signature associated with impaired tissue remodelling and increased carnitine metabolism. Conversely, ACA-positive patients exhibited an immunomodulation and tissue homeostasis signature and increased phospholipid metabolism., Conclusions: Patients with SSc display varying biological profiles based on their serological status. The findings highlight the potential utility of serological status as a discriminating factor in disease severity and suggest its relevance in tailoring treatment strategies and future research directions., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2024. Re-use permitted under CC BY. Published by BMJ on behalf of EULAR.)

The heterogeneity of systemic lupus erythematosus (SLE) can be explained by epigenetic alterations that disrupt transcriptional programs mediating environmental and genetic risk. This study evaluated the epigenetic contribution to SLE heterogeneity considering molecular and serological subtypes, genetics and transcriptional status, followed by drug target discovery. We performed a stratified epigenome-wide association studies of whole blood DNA methylation from 213 SLE patients and 221 controls. Methylation quantitative trait loci analyses, cytokine and transcription factor activity - epigenetic associations and methylation-expression correlations were conducted. New drug targets were searched for based on differentially methylated genes. In a stratified approach, a total of 974 differential methylation CpG sites with dependency on molecular subtypes and autoantibody profiles were found. Mediation analyses suggested that SLE-associated SNPs in the HLA region exert their risk through DNA methylation changes. Novel genetic variants regulating DNAm in disease or in specific molecular contexts were identified. The epigenetic landscapes showed strong association with transcription factor activity and cytokine levels, conditioned by the molecular context. Epigenetic signals were enriched in known and novel drug targets for SLE. This study reveals possible genetic drivers and consequences of epigenetic variability on SLE heterogeneity and disentangles the DNAm mediation role on SLE genetic risk and novel disease-specific meQTLs. Finally, novel targets for drug development were discovered., (© 2024. The Author(s).)

Background: Giant cell arteritis is an age-related vasculitis that mainly affects the aorta and its branches in individuals aged 50 years and older. Current options for diagnosis and treatment are scarce, highlighting the need to better understand its underlying pathogenesis. Genome-wide association studies (GWAS) have emerged as a powerful tool for unravelling the pathogenic mechanisms involved in complex diseases. We aimed to characterise the genetic basis of giant cell arteritis by performing the largest GWAS of this vasculitis to date and to assess the functional consequences and clinical implications of identified risk loci., Methods: We collected and meta-analysed genomic data from patients with giant cell arteritis and healthy controls of European ancestry from ten cohorts across Europe and North America. Eligible patients required confirmation of giant cell arteritis diagnosis by positive temporal artery biopsy, positive temporal artery doppler ultrasonography, or imaging techniques confirming large-vessel vasculitis. We assessed the functional consequences of loci associated with giant cell arteritis using cell enrichment analysis, fine-mapping, and causal gene prioritisation. We also performed a drug repurposing analysis and developed a polygenic risk score to explore the clinical implications of our findings., Findings: We included a total of 3498 patients with giant cell arteritis and 15 550 controls. We identified three novel loci associated with risk of giant cell arteritis. Two loci, MFGE8 (rs8029053; p=4·96 × 10 -8 ; OR 1·19 [95% CI 1·12-1·26]) and VTN (rs704; p=2·75 × 10 -9 ; OR 0·84 [0·79-0·89]), were related to angiogenesis pathways and the third locus, CCDC25 (rs11782624; p=1·28 × 10 -8 ; OR 1·18 [1·12-1·25]), was related to neutrophil extracellular traps (NETs). We also found an association between this vasculitis and HLA region and PLG. Variants associated with giant cell arteritis seemed to fulfil a specific regulatory role in crucial immune cell types. Furthermore, we identified several drugs that could represent promising candidates for treatment of this disease. The polygenic risk score model was able to identify individuals at increased risk of developing giant cell arteritis (90th percentile OR 2·87 [95% CI 2·15-3·82]; p=1·73 × 10 -13 )., Interpretation: We have found several additional loci associated with giant cell arteritis, highlighting the crucial role of angiogenesis in disease susceptibility. Our study represents a step forward in the translation of genomic findings to clinical practice in giant cell arteritis, proposing new treatments and a method to measure genetic predisposition to this vasculitis., Funding: Institute of Health Carlos III, Spanish Ministry of Science and Innovation, UK Medical Research Council, and National Institute for Health and Care Research., Competing Interests: Declaration of interests MCC reports support from the Spanish Ministry of Science and Innovation (PID2020-114909RB-I00), Vasculitis Foundation, Agency for the Management of University and Research Grants (2021 SGR 01561), and Kiniksa Pharmaceuticals; consulting fees or honoraria from GSK, CSL Vifor, AbbVie, and AstraZeneca; support for attending meetings from Kiniksa Pharmaceuticals; and participation on a data safety monitoring board or advisory board for GSK, CSL Vifor, and AstraZeneca. GE has acted as a consultant for GSK, AstraZeneca, Sobi, Novartis, Boehringer, and CSL Vifor. AWM reports support from the UK Medical Research Council (MRC), National Institute for Health and Care Research (NIHR), Leeds Care, and Roche Products; and consulting fees or honoraria from CSL Vifor and AstraZeneca. PL reports grants or contracts from the Federal Ministry of Education and Research, German Research Society, German Society for Rheumatology, John Grube Foundation, and CSL Vifor; consulting fees or honoraria from GSK, CSL Vifor, AstraZeneca, Bristol Myers Squibb, Boehringer Ingelheim, Forum für medizinische Fortbildung, Janssen, Rheumaakademie, and UCB; support for attending meetings from CSL Vifor; and participation on a data safety monitoring board or advisory board for GSK, CSL Vifor, AbbVie, and Novartis. TW reports consulting fees or honoraria from AbbVie, AstraZeneca, Lilly, UCB, and Novartis; and participation on a data safety monitoring board or advisory board for AbbVie, AstraZeneca, Lilly, UCB, Novartis, and Fresenius. MAG-G reports honoraria from GSK. NK reports grants or contracts from Bristol Myers Squibb, AbbVie, and Sanofi; and consulting fees or honoraria from Roche, Otsuka, GSK, and Mallinckrodt. CAL reports grants or contracts from Bristol Myers Squibb and support from the National Institutes of Health. PAM reports grants, contracts, or consulting fees from AbbVie, Amgen, AstraZeneca, ArGenx, Boehringer Ingelheim, Bristol Myers Squibb, Cabaletta, CSL Behring, Eicos, Electra, Forbius, Genentech–Roche, GSK, HiBio, InflaRx, Janssen, Jubilant, Kyverna, MiroBio, Neutrolis, Novartis, NS Pharma, Q32, Regeneron, Sanofi, Sparrow, Takeda, and Vistera; royalties or licenses from UpToDate; and stock or stock options from Kyverna, Q32, and Sparrow. SLM reports grants or contracts from MRC, NIHR, and CSL Vifor; consulting fees from Roche, Sanofi, AbbVie, AstraZeneca, and Pfizer; payment or honoraria for lectures or educational events from Roche, Pfizer, UCB, CSL Vifor, Fresenius Kabi, and Novartis; support for attending meetings from Pfizer; participation on a data safety monitoring board or advisory board for Collaboration for Leadership in Applied Health Research and Care, Haywood Foundation, and GC-SheaLD; a leadership or fiduciary role in the British Society for Rheumatology Clinical Affairs Committee; participation as an investigator on industry-sponsored clinical trials for Sanofi; and infrastructure support from MRC. LB reports payment or honoraria for lectures, presentations, speakers bureaus, manuscript writing, or educational events from Instrumentation Laboratory SPA and AbbVie; and support for attending meetings from AbbVie and Novartis. EB reports payments or honoraria from EULAR and received grants from the Dutch Arthritis Society DAS and the EU/EFPIA/Innovative Medicines Initiative 2 Joint Undertaking Immune-Image grant no 831514. EB is member of the board of the non-profit organisation, Auto-immune Research Hub, in the Netherlands. All other authors declare no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Ltd. This is an Open Access article under the CC BY 4.0 license. Published by Elsevier Ltd.. All rights reserved.)

Objective: Erdheim-Chester disease (ECD) is rare histiocytosis with a wide range of clinical manifestations. Somatic mutations are key to the pathogenesis of the disease; however, the relationship between germline genetic variants and ECD has not been examined so far. The present study aims to explore the inherited genetic component of ECD by performing the first genome-wide association study., Methods: After quality controls, a cohort of 255 patients with ECD and 7,471 healthy donors was included in this study. Afterward, a logistic regression followed by in silico functional annotation was performed., Results: A signal at the 18q12.3 genomic region was identified as a new susceptibility locus for ECD (P = 2.75 × 10 -11 ; Odds Ratio = 2.09). This association was annotated to the SETBP1 gene, which is involved in clonal haematopoiesis. Functional annotation of this region and of the identified suggestive signals revealed additional genes that could be potentially involved in the pathogenesis of the disease., Conclusion: Overall, this work demonstrates that germline genetic variants can impact on the development of ECD and suggests new pathways with a potential pathogenic role., (© 2023 The Authors. Arthritis & Rheumatology published by Wiley Periodicals LLC on behalf of American College of Rheumatology.)

Immune-mediated inflammatory diseases (IMIDs) comprise a complex group of pathologies with diverse etiologies and clinical manifestations. In particular, omics technologies have remodeled our understanding of a set of IMIDs such as systemic autoimmune rheumatic diseases (SARDs), generating vast amounts of data on the genome, epigenome, transcriptome, proteome and metabolome of immune cells and SARDs patients. However, the integration of omics data to advance our knowledge of these diseases is challenging, requiring advanced bioinformatic tools. This review explores different multi-omic integrative tools for refining previous research, exploring the biological relevance of datasets within different contexts, or translating omics results into clinical advances. We also discuss relevant multi-omic studies in SARDs research and the potential of omics data from available repositories to complement ongoing investigation in this field., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2023 The Authors.)

Systemic sclerosis (SSc) is a complex disease that affects the connective tissue, causing fibrosis. SSc patients show altered immune cell composition and activation in the peripheral blood (PB). PB monocytes (Mos) are recruited into tissues where they differentiate into macrophages, which are directly involved in fibrosis. To understand the role of CD14 + PB Mos in SSc, a single-cell transcriptome analysis (scRNA-seq) was conducted on 8 SSc patients and 8 controls. Using unsupervised clustering methods, CD14 + cells were assigned to 11 clusters, which added granularity to the known monocyte subsets: classical (cMos), intermediate (iMos) and non-classical Mos (ncMos) or type 2 dendritic cells. NcMos were significantly overrepresented in SSc patients and showed an active IFN-signature and increased expression levels of PTGES, in addition to monocyte motility and adhesion markers. We identified a SSc-related cluster of IRF7+ STAT1+ iMos with an aberrant IFN-response. Finally, a depletion of M2 polarised cMos in SSc was observed. Our results highlighted the potential of PB Mos as biomarkers for SSc and provided new possibilities for putative drug targets for modulating the innate immune response in SSc., (Copyright © 2023 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Objectives: The number of susceptibility loci currently associated with vasculitis is lower than in other immune-mediated diseases due in part to small cohort sizes, a consequence of the low prevalence of vasculitides. This study aimed to identify new genetic risk loci for the main systemic vasculitides through a comprehensive analysis of their genetic overlap., Methods: Genome-wide data from 8467 patients with any of the main forms of vasculitis and 29 795 healthy controls were meta-analysed using ASSET. Pleiotropic variants were functionally annotated and linked to their target genes. Prioritised genes were queried in DrugBank to identify potentially repositionable drugs for the treatment of vasculitis., Results: Sixteen variants were independently associated with two or more vasculitides, 15 of them representing new shared risk loci. Two of these pleiotropic signals, located close to CTLA4 and CPLX1 , emerged as novel genetic risk loci in vasculitis. Most of these polymorphisms appeared to affect vasculitis by regulating gene expression. In this regard, for some of these common signals, potential causal genes were prioritised based on functional annotation, including CTLA4 , RNF145 , IL12B , IL5 , IRF1 , IFNGR1 , PTK2B , TRIM35 , EGR2 and ETS2 , each of which has key roles in inflammation. In addition, drug repositioning analysis showed that several drugs, including abatacept and ustekinumab, could be potentially repurposed in the management of the analysed vasculitides., Conclusions: We identified new shared risk loci with functional impact in vasculitis and pinpointed potential causal genes, some of which could represent promising targets for the treatment of vasculitis., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2023. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Objective: Systemic sclerosis (SSc) is a complex autoimmune disease with a strong genetic component. However, most of the genes associated with the disease are still unknown because associated variants affect mostly noncoding intergenic elements of the genome. We used functional genomics to translate the genetic findings into a better understanding of the disease., Methods: Promoter capture Hi-C and RNA-sequencing experiments were performed in CD4+ T cells and CD14+ monocytes from 10 SSc patients and 5 healthy controls to link SSc-associated variants with their target genes, followed by differential expression and differential interaction analyses between cell types., Results: We linked SSc-associated loci to 39 new potential target genes and confirmed 7 previously known SSc-associated genes. We highlight novel causal genes, such as CXCR5, as the most probable candidate gene for the DDX6 locus. Some previously known SSc-associated genes, such as IRF8, STAT4, and CD247, showed cell type-specific interactions. We also identified 15 potential drug targets already in use in other similar immune-mediated diseases that could be repurposed for SSc treatment. Furthermore, we observed that interactions were directly correlated with the expression of important genes implicated in cell type-specific pathways and found evidence that chromatin conformation is associated with genotype., Conclusion: Our study revealed potential causal genes for SSc-associated loci, some of them acting in a cell type-specific manner, suggesting novel biologic mechanisms that might mediate SSc pathogenesis., (© 2022 The Authors. Arthritis & Rheumatology published by Wiley Periodicals LLC on behalf of American College of Rheumatology.)

Objectives: rs76428106-C, a low frequency polymorphism that affects the splicing of the FLT3 gene, has recently been associated with several seropositive autoimmune diseases. Here, we aimed to evaluate the potential implication of rs76428106-C in the susceptibility to systemic sclerosis (SSc)., Methods: We analysed a total of 26 598 European ancestry individuals, 9063 SSc and 17 535 healthy controls, to test the association between FLT3 rs76428106-C and SSc and its different subphenotypes. Genotype data of rs76428106 were obtained by imputation of already available genome-wide association study data and analysed by logistic regression analysis., Results: In accordance with that observed in other autoimmune disorders, the FLT3 rs76428106-C allele was significantly increased [P-value = 2.03 × 10-3, odds ratio (OR) = 1.34] in SSc patients compared with healthy controls. A similar risk effect was found when the main SSc clinical and serological subgroups were compared with controls. When comparing SSc patients with and without digital ulcers (DU), the rs76428106-C frequency was significantly increased in DU-positive SSc patients in comparison with DU-negative patients (P-value = 0.036, OR = 2.16)., Conclusion: This study is the first to report an association between rs76428176-C and SSc. Our results support the role of FLT3 as a relevant gene in seropositive immune-mediated diseases and a potential biomarker for SSc microangiopathy., (© The Author(s) 2022. Published by Oxford University Press on behalf of the British Society for Rheumatology.)

Rheumatoid arthritis (RA) is a highly heritable complex disease with unknown etiology. Multi-ancestry genetic research of RA promises to improve power to detect genetic signals, fine-mapping resolution and performances of polygenic risk scores (PRS). Here, we present a large-scale genome-wide association study (GWAS) of RA, which includes 276,020 samples from five ancestral groups. We conducted a multi-ancestry meta-analysis and identified 124 loci (P < 5 × 10 -8 ), of which 34 are novel. Candidate genes at the novel loci suggest essential roles of the immune system (for example, TNIP2 and TNFRSF11A) and joint tissues (for example, WISP1) in RA etiology. Multi-ancestry fine-mapping identified putatively causal variants with biological insights (for example, LEF1). Moreover, PRS based on multi-ancestry GWAS outperformed PRS based on single-ancestry GWAS and had comparable performance between populations of European and East Asian ancestries. Our study provides several insights into the etiology of RA and improves the genetic predictability of RA., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Obesity contributes to a chronic proinflammatory state, which is a known risk factor to develop immune-mediated diseases. However, its role in systemic sclerosis (SSc) remains to be elucidated. Therefore, we conducted a two-sample mendelian randomization (2SMR) study to analyze the effect of three body fat distribution parameters in SSc. As instrumental variables, we used the allele effects described for single nucleotide polymorphisms (SNPs) in different genome-wide association studies (GWAS) for SSc, body mass index (BMI), waist-to-hip ratio (WHR) and WHR adjusted for BMI (WHRadjBMI). We performed local (pHESS) and genome-wide (LDSC) genetic correlation analyses between each of the traits and SSc and we applied several Mendelian randomization (MR) methods (i.e., random effects inverse-variance weight, MR-Egger regression, MR pleiotropy residual sum and outlier method and a multivariable model). Our results show no genetic correlation or causal relationship between any of these traits and SSc. Nevertheless, we observed a negative causal association between WHRadjBMI and SSc, which might be due to the effect of gastrointestinal complications suffered by the majority of SSc patients. In conclusion, reverse causality might be an especially difficult confounding factor to define the effect of obesity in the onset of SSc.

Chagas disease is an infection caused by the parasite Trypanosoma cruzi, endemic in Latino America. Leveraging the three-way admixture between Native American (AMR), European (EUR) and African (AFR) populations in Latin Americans, we aimed to better understand the genetic basis of Chagas disease by performing an admixture mapping study in a Colombian population. A two-stage study was conducted, and subjects were classified as seropositive and seronegative for T. cruzi. In stage 1, global and local ancestries were estimated using reference data from the 1000 Genomes Project (1KGP), and local ancestry associations were performed by logistic regression models. The AMR ancestry showed a protective association with Chagas disease within the major histocompatibility complex region [Odds ratio (OR) = 0.74, 95% confidence interval (CI) = 0.66-0.83, lowest P-value = 4.53 × 10-8]. The fine mapping assessment on imputed genotypes combining data from stage 1 and 2 from an independent Colombian cohort, revealed nominally associated variants in high linkage disequilibrium with the top signal (rs2032134, OR = 0.93, 95% CI = 0.90-0.97, P-value = 3.54 × 10-4) in the previously associated locus. To assess ancestry-specific adaptive signals, a selective sweep scan in an AMR reference population from 1KGP together with an in silico functional analysis highlighted the Tripartite Motif family and the human leukocyte antigen genes, with crucial role in the immune response against pathogens. Furthermore, these analyses emphasized the macrophages, neutrophils and eosinophils, as key players in the defense against T. cruzi. This first admixture mapping study in Chagas disease provided novel insights underlying the host immune response in the pathogenesis of this neglected disease., (© The Author(s) 2021. Published by Oxford University Press.)

Eukaryotic cells can survive the loss of their mitochondrial genome, but consequently suffer from severe growth defects. 'Petite yeasts', characterized by mitochondrial genome loss, are instrumental for studying mitochondrial function and physiology. However, the molecular cause of their reduced growth rate remains an open question. Here we show that petite cells suffer from an insufficient capacity to synthesize glutamate, glutamine, leucine and arginine, negatively impacting their growth. Using a combination of molecular genetics and omics approaches, we demonstrate the evolution of fast growth overcomes these amino acid deficiencies, by alleviating a perturbation in mitochondrial iron metabolism and by restoring a defect in the mitochondrial tricarboxylic acid cycle, caused by aconitase inhibition. Our results hence explain the slow growth of mitochondrial genome-deficient cells with a partial auxotrophy in four amino acids that results from distorted iron metabolism and an inhibited tricarboxylic acid cycle., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

A recent genome-wide association study (GWAS) identified a locus in chromosome 11 associated with the chronic cardiac form of Chagas disease. Here we aimed to elucidate the potential functional mechanism underlying this genetic association by analyzing the correlation among single nucleotide polymorphisms (SNPs) and DNA methylation (DNAm) levels as cis methylation quantitative trait loci (cis-mQTL) within this region. A total of 2,611 SNPs were tested against 2,647 DNAm sites, in a subset of 37 chronic Chagas cardiomyopathy patients and 20 asymptomatic individuals from the GWAS. We identified 6,958 significant cis-mQTLs (False Discovery Rate [FDR]<0.05) at 1 Mb each side of the GWAS leading variant, where six of them potentially modulate the expression of the SAC3D1 gene, the reported gene in the previous GWAS. In addition, a total of 268 cis-mQTLs showed differential methylation between chronic Chagas cardiomyopathy patients and asymptomatic individuals. The most significant cis-mQTLs mapped in the gene bodies of POLA2 (FDR = 1.04x10-11), PLAAT3 (FDR = 7.22x10-03), and CCDC88B (FDR = 1.89x10-02) that have been associated with cardiovascular and hematological traits in previous studies. One of the most relevant interactions correlated with hypermethylation of CCDC88B. This gene is involved in the inflammatory response, and its methylation and expression levels have been previously reported in Chagas cardiomyopathy. Our findings support the functional relevance of the previously associated genomic region, highlighting the regulation of novel genes that could play a role in the chronic cardiac form of the disease., Competing Interests: The authors have declared that no competing interests exist.

Characterizing genetic influences on DNA methylation (DNAm) provides an opportunity to understand mechanisms underpinning gene regulation and disease. In the present study, we describe results of DNAm quantitative trait locus (mQTL) analyses on 32,851 participants, identifying genetic variants associated with DNAm at 420,509 DNAm sites in blood. We present a database of >270,000 independent mQTLs, of which 8.5% comprise long-range (trans) associations. Identified mQTL associations explain 15-17% of the additive genetic variance of DNAm. We show that the genetic architecture of DNAm levels is highly polygenic. Using shared genetic control between distal DNAm sites, we constructed networks, identifying 405 discrete genomic communities enriched for genomic annotations and complex traits. Shared genetic variants are associated with both DNAm levels and complex diseases, but only in a minority of cases do these associations reflect causal relationships from DNAm to trait or vice versa, indicating a more complex genotype-phenotype map than previously anticipated., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Objective: The greatest genetic effect reported for systemic sclerosis (SSc) lies in the major histocompatibility complex (MHC) locus. Leveraging the largest SSc genome-wide association study, we aimed to fine-map this region to identify novel human leucocyte antigen (HLA) genetic variants associated with SSc susceptibility and its main clinical and serological subtypes., Methods: 9095 patients with SSc and 17 584 controls genome-wide genotyped were used to impute and test single-nucleotide polymorphisms (SNPs) across the MHC, classical HLA alleles and their composite amino acid residues. Additionally, patients were stratified according to their clinical and serological status, namely, limited cutaneous systemic sclerosis (lcSSc), diffuse cutaneous systemic sclerosis (dcSSc), anticentromere (ACA), antitopoisomerase (ATA) and anti-RNApolIII autoantibodies (ARA)., Results: Sequential conditional analyses showed nine SNPs, nine classical alleles and seven amino acids that modelled the observed associations with SSc. This confirmed previously reported associations with HLA-DRB1*11:04 and HLA-DPB1*13:01 , and revealed a novel association of HLA-B*08:01 . Stratified analyses showed specific associations of HLA-DQA1*02:01 with lcSSc, and an exclusive association of HLA-DQA1*05:01 with dcSSc. Similarly, private associations were detected in HLA-DRB1*08:01 and confirmed the previously reported association of HLA-DRB1*07:01 with ACA-positive patients, as opposed to the HLA-DPA1*02:01 and HLA-DQB1*03:01 alleles associated with ATA presentation., Conclusions: This study confirms the contribution of HLA class II and reveals a novel association of HLA class I with SSc, suggesting novel pathways of disease pathogenesis. Furthermore, we describe specific HLA associations with SSc clinical and serological subtypes that could serve as biomarkers of disease severity and progression., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2021. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Objective: To identify the genetic variants that affect gene expression (expression quantitative trait loci [eQTLs]) in systemic sclerosis (SSc) and to investigate their role in the pathogenesis of the disease., Methods: We performed an eQTL analysis using whole-blood sequencing data from 333 SSc patients and 524 controls and integrated them with SSc genome-wide association study (GWAS) data. We integrated our findings from expression modeling, differential expression analysis, and transcription factor binding site enrichment with key clinical features of SSc., Results: We detected 49,123 validated cis-eQTLs from 4,539 SSc-associated single-nucleotide polymorphisms (SNPs) (P GWAS < 10 -5 ). A total of 1,436 genes were within 1 Mb of the 4,539 SSc-associated SNPs. Of those 1,436 genes, 565 were detected as having ≥1 eQTL with an SSc-associated SNP. We developed a strategy to prioritize disease-associated genes based on their expression variance explained by SSc eQTLs (r 2 > 0.05). As a result, 233 candidates were identified, 134 (58%) of them associated with hallmarks of SSc and 105 (45%) of them differentially expressed in the blood cells, skin, or lung tissue of SSc patients. Transcription factor binding site analysis revealed enriched motifs of 24 transcription factors (5%) among SSc eQTLs, 5 of which were found to be differentially regulated in the blood cells (ELF1 and MGA), skin (KLF4 and ID4), and lungs (TBX4) of SSc patients. Ten candidate genes (4%) can be targeted by approved medications for immune-mediated diseases, of which only 3 have been tested in clinical trials in patients with SSc., Conclusion: The findings of the present study indicate a new layer to the molecular complexity of SSc, contributing to a better understanding of the pathogenesis of the disease., (© 2021, American College of Rheumatology.)

Objective: Clinical heterogeneity, a hallmark of systemic autoimmune diseases, impedes early diagnosis and effective treatment, issues that may be addressed if patients could be classified into groups defined by molecular pattern. This study was undertaken to identify molecular clusters for reclassifying systemic autoimmune diseases independently of clinical diagnosis., Methods: Unsupervised clustering of integrated whole blood transcriptome and methylome cross-sectional data on 955 patients with 7 systemic autoimmune diseases and 267 healthy controls was undertaken. In addition, an inception cohort was prospectively followed up for 6 or 14 months to validate the results and analyze whether or not cluster assignment changed over time., Results: Four clusters were identified and validated. Three were pathologic, representing "inflammatory," "lymphoid," and "interferon" patterns. Each included all diagnoses and was defined by genetic, clinical, serologic, and cellular features. A fourth cluster with no specific molecular pattern was associated with low disease activity and included healthy controls. A longitudinal and independent inception cohort showed a relapse-remission pattern, where patients remained in their pathologic cluster, moving only to the healthy one, thus showing that the molecular clusters remained stable over time and that single pathogenic molecular signatures characterized each individual patient., Conclusion: Patients with systemic autoimmune diseases can be jointly stratified into 3 stable disease clusters with specific molecular patterns differentiating different molecular disease mechanisms. These results have important implications for future clinical trials and the study of nonresponse to therapy, marking a paradigm shift in our view of systemic autoimmune diseases., (© 2020, American College of Rheumatology.)

Objectives: Genomic Risk Scores (GRS) successfully demonstrated the ability of genetics to identify those individuals at high risk for complex traits including immune-mediated inflammatory diseases (IMIDs). We aimed to test the performance of GRS in the prediction of risk for systemic sclerosis (SSc) for the first time., Methods: Allelic effects were obtained from the largest SSc Genome-Wide Association Study (GWAS) to date (9 095 SSc and 17 584 healthy controls with European ancestry). The best-fitting GRS was identified under the additive model in an independent cohort that comprised 400 patients with SSc and 571 controls. Additionally, GRS for clinical subtypes (limited cutaneous SSc and diffuse cutaneous SSc) and serological subtypes (anti-topoisomerase positive (ATA+) and anti-centromere positive (ACA+)) were generated. We combined the estimated GRS with demographic and immunological parameters in a multivariate generalised linear model., Results: The best-fitting SSc GRS included 33 single nucleotide polymorphisms (SNPs) and discriminated between patients with SSc and controls (area under the receiver operating characteristic (ROC) curve (AUC)=0.673). Moreover, the GRS differentiated between SSc and other IMIDs, such as rheumatoid arthritis and Sjögren's syndrome. Finally, the combination of GRS with age and immune cell counts significantly increased the performance of the model (AUC=0.787). While the SSc GRS was not able to discriminate between ATA+ and ACA+ patients (AUC<0.5), the serological subtype GRS, which was based on the allelic effects observed for the comparison between ACA+ and ATA+ patients, reached an AUC=0.693., Conclusions: GRS was successfully implemented in SSc. The model discriminated between patients with SSc and controls or other IMIDs, confirming the potential of GRS to support early and differential diagnosis for SSc., Competing Interests: Competing interests: LB-C: none; GV-M: none; MK: none; MA-H: none; ELI: none; International SSc Group: none; PRECISESADS Consortium: none; MEAl-R: none; LB: none; JM: none., (© Author(s) (or their employer(s)) 2021. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Objectives: The analysis of annotated transcripts from genome-wide expression studies may help to understand the pathogenesis of complex diseases, such as systemic sclerosis (SSc). We performed a whole blood (WB) transcriptome analysis on RNA collected in the context of the European PRECISESADS project, aiming at characterising the pathways that differentiate SSc from controls and that are reproducible in geographically diverse populations., Methods: Samples from 162 patients and 252 controls were collected in RNA stabilisers. Cases and controls were divided into a discovery (n=79+163; Southern Europe) and validation cohort (n=83+89; Central-Western Europe). RNA sequencing was performed by an Illumina assay. Functional annotations of Reactome pathways were performed with the Functional Analysis of Individual Microarray Expression (FAIME) algorithm. In parallel, immunophenotyping of 28 circulating cell populations was performed. We tested the presence of differentially expressed genes/pathways and the correlation between absolute cell counts and RNA transcripts/FAIME scores in regression models. Results significant in both populations were considered as replicated., Results: Overall, 15 224 genes and 1277 functional pathways were available; of these, 99 and 225 were significant in both sets. Among replicated pathways, we found a deregulation in type-I interferon, Toll-like receptor cascade, tumour suppressor p53 protein function, platelet degranulation and activation. RNA transcripts or FAIME scores were jointly correlated with cell subtypes with strong geographical differences; neutrophils were the major determinant of gene expression in SSc-WB samples., Conclusions: We discovered a set of differentially expressed genes/pathways validated in two independent sets of patients with SSc, highlighting a number of deregulated processes that have relevance for the pathogenesis of autoimmunity and SSc., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2020. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Mixed Connective Tissue Disease (MCTD) is a rare complex systemic autoimmune disease (SAD) characterized by the presence of increased levels of anti-U1 ribonucleoprotein autoantibodies and signs and symptoms that resemble other SADs such as systemic sclerosis (SSc), rheumatoid arthritis (RA), and systemic lupus erythematosus (SLE). Due to its low prevalence, this disease has been very poorly studied at the molecular level. We performed for the first time an epigenome-wide association study interrogating DNA methylation data obtained with the Infinium MethylationEPIC array from whole blood samples in 31 patients diagnosed with MCTD and 255 healthy subjects. We observed a pervasive hypomethylation involving 170 genes enriched for immune-related function such as those involved in type I interferon signaling pathways or in negative regulation of viral genome replication. We mostly identified epigenetic signals at genes previously implicated in other SADs, for example MX1, PARP9, DDX60 , or IFI44L , for which we also observed that MCTD patients exhibit higher DNA methylation variability compared with controls, suggesting that these sites might be involved in plastic immune responses that are relevant to the disease. Through methylation quantitative trait locus (meQTL) analysis we identified widespread local genetic effects influencing DNA methylation variability at MCTD-associated sites. Interestingly, for IRF7, IFI44 genes, and the HLA region we have evidence that they could be exerting a genetic risk on MCTD mediated through DNA methylation changes. Comparison of MCTD-associated epigenome with patients diagnosed with SLE, or Sjögren's Syndrome, reveals a common interferon-related epigenetic signature, however we find substantial epigenetic differences when compared with patients diagnosed with rheumatoid arthritis and systemic sclerosis. Furthermore, we show that MCTD-associated CpGs are potential epigenetic biomarkers with high diagnostic value. Our study serves to reveal new genes and pathways involved in MCTD, to illustrate the important role of epigenetic modifications in MCTD pathology, in mediating the interaction between different genetic and environmental MCTD risk factors, and as potential biomarkers of SADs.

Objective: Immune-mediated inflammatory diseases (IMIDs) are heterogeneous and complex conditions with overlapping clinical symptoms and elevated familial aggregation, which suggests the existence of a shared genetic component. In order to identify this genetic background in a systematic fashion, we performed the first cross-disease genome-wide meta-analysis in systemic seropositive rheumatic diseases, namely, systemic sclerosis, systemic lupus erythematosus, rheumatoid arthritis and idiopathic inflammatory myopathies., Methods: We meta-analysed ~6.5 million single nucleotide polymorphisms in 11 678 cases and 19 704 non-affected controls of European descent populations. The functional roles of the associated variants were interrogated using publicly available databases., Results: Our analysis revealed five shared genome-wide significant independent loci that had not been previously associated with these diseases: NAB1 , KPNA4-ARL14 , DGQK , LIMK1 and PRR12 . All of these loci are related with immune processes such as interferon and epidermal growth factor signalling, response to methotrexate, cytoskeleton dynamics and coagulation cascade. Remarkably, several of the associated loci are known key players in autoimmunity, which supports the validity of our results. All the associated variants showed significant functional enrichment in DNase hypersensitivity sites, chromatin states and histone marks in relevant immune cells, including shared expression quantitative trait loci . Additionally, our results were significantly enriched in drugs that are being tested for the treatment of the diseases under study., Conclusions: We have identified shared new risk loci with functional value across diseases and pinpoint new potential candidate loci that could be further investigated. Our results highlight the potential of drug repositioning among related systemic seropositive rheumatic IMIDs., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2019. No commercial re-use. See rights and permissions. Published by BMJ.)

Background: In recent years, research has consistently proven the occurrence of genetic overlap across autoimmune diseases, which supports the existence of common pathogenic mechanisms in autoimmunity. The objective of this study was to further investigate this shared genetic component., Methods: For this purpose, we performed a cross-disease meta-analysis of Immunochip data from 37,159 patients diagnosed with a seropositive autoimmune disease (11,489 celiac disease (CeD), 15,523 rheumatoid arthritis (RA), 3477 systemic sclerosis (SSc), and 6670 type 1 diabetes (T1D)) and 22,308 healthy controls of European origin using the R package ASSET., Results: We identified 38 risk variants shared by at least two of the conditions analyzed, five of which represent new pleiotropic loci in autoimmunity. We also identified six novel genome-wide associations for the diseases studied. Cell-specific functional annotations and biological pathway enrichment analyses suggested that pleiotropic variants may act by deregulating gene expression in different subsets of T cells, especially Th17 and regulatory T cells. Finally, drug repositioning analysis evidenced several drugs that could represent promising candidates for CeD, RA, SSc, and T1D treatment., Conclusions: In this study, we have been able to advance in the knowledge of the genetic overlap existing in autoimmunity, thus shedding light on common molecular mechanisms of disease and suggesting novel drug targets that could be explored for the treatment of the autoimmune diseases studied.

The Revised International Prognostic Score and some somatic mutations in myelodysplastic syndrome (MDS) are independently associated with transformation to acute myeloid leukemia (AML). Immunity has also been implicated in the pathogenesis of MDS, although the underlying mechanism remains unclear. We performed a SNP array on chromosome 6 in CD34 + purified blasts from 19 patients diagnosed with advanced MDS and 8 patients with other myeloid malignancies to evaluate the presence of loss of heterozygosity (LOH) in HLA and its impact on disease progression. Three patients had acquired copy-neutral LOH (CN-LOH) on 6p arms, which may disrupt antigen presentation and act as a mechanism for immune system evasion. Interestingly, these patients had previously been classified at low risk of AML progression, and the poor outcome cannot be explained by the acquisition of adverse mutations. LOH HLA was not detected in the remaining 24 patients, who all had adverse risk factors. In summary, the clinical outcome of patients with advanced MDS might be influenced by HLA allelic loss, wich allows subclonal expansions to evade cytotoxic-T and NK cell attack. CN-LOH HLA may therefore be a factor favoring MDS progression to AML independently of the somatic tumor mutation load., Competing Interests: CONFLICTS OF INTEREST No potential conflicts of interest were disclosed.

Background: Non-small cell lung cancer (NSCLC) is the most common cause of cancer-related deaths worldwide and is primarily treated with radiation, surgery, and platinum-based drugs like cisplatin and carboplatin. The major challenge in the treatment of NSCLC patients is intrinsic or acquired resistance to chemotherapy. Molecular markers predicting the outcome of the patients are urgently needed., Methods: Here, we employed patient-derived xenografts (PDXs) to detect predictive methylation biomarkers for platin-based therapies. We used MeDIP-Seq to generate genome-wide DNA methylation profiles of 22 PDXs, their parental primary NSCLC, and their corresponding normal tissues and complemented the data with gene expression analyses of the same tissues. Candidate biomarkers were validated with quantitative methylation-specific PCRs (qMSP) in an independent cohort., Results: Comprehensive analyses revealed that differential methylation patterns are highly similar, enriched in PDXs and lung tumor-specific when comparing differences in methylation between PDXs versus primary NSCLC. We identified a set of 40 candidate regions with methylation correlated to carboplatin response and corresponding inverse gene expression pattern even before therapy. This analysis led to the identification of a promoter CpG island methylation of LDL receptor-related protein 12 (LRP12) associated with increased resistance to carboplatin. Validation in an independent patient cohort (n = 35) confirmed that LRP12 methylation status is predictive for therapeutic response of NSCLC patients to platin therapy with a sensitivity of 80% and a specificity of 84% (p < 0.01). Similarly, we find a shorter survival time for patients with LRP12 hypermethylation in the TCGA data set for NSCLC (lung adenocarcinoma)., Conclusions: Using an epigenome-wide sequencing approach, we find differential methylation patterns from primary lung cancer and PDX-derived cancers to be very similar, albeit with a lower degree of differential methylation in primary tumors. We identify LRP12 DNA methylation as a powerful predictive marker for carboplatin resistance. These findings outline a platform for the identification of epigenetic therapy resistance biomarkers based on PDX NSCLC models.

Immune-checkpoint inhibitors show encouraging results in cancer treatment, but the clinical benefit is limited exclusively to a subset of patients. We analyzed the density and composition of tumor T-cell infiltration in non-small-cell lung carcinoma (NSCLC) in relation to PD-L1 and HLA class I (HLA-I) expression. We found that positive HLA-I expression, independently on PD-L1 status, is the key factor determining the increased density of the immune infiltrate. When both markers were analyzed simultaneously, we identified four phenotypes of HLA-I and PD-L1 co-expression. They demonstrated different patterns of tumor infiltration and clinicopathologic characteristics, including the tumor size and lymphatic spread. All HLA-I+/PD-L1+ tumors had a high degree of intratumoral infiltration with CD8+T-lymphocytes, whereas HLA-I loss was associated with a significantly reduced number of tumor infiltrating T-lymphocytes mostly restrained in the stroma surrounding the tumor nest. HLA-I-negative/PD-L1-positive tumors had bigger size (T) and lower grade of infiltration with CD8+T-cells. It represents a cancer immune escape phenotype that combines two independent mechanisms of immune evasion: loss of HLA-I and upregulation of PD-L1. Using GCH-array analysis of human lung cancer cell lines we found that the loss of heterozygosity (LOH) with complete or partial deletion of HLA-I genes is the principal mechanism of HLA-I alterations. This irreversible defect, which could potentially decrease the clinical efficacy of lung cancer immunotherapy, appears to be underestimated. In conclusion, our results suggest that the analysis of HLA-I is very important for the selection of potential responders to cancer immunotherapy., Competing Interests: CONFLICTS OF INTEREST The authors have no conflicts of interest to disclose.

Background: Watery diarrhea is the cardinal symptom of lymphocytic colitis (LC). We have previously shown that colonic Na malabsorption is one of the major pathologic alterations of LC and found evidence for an epithelial barrier defect. On these grounds, this study aimed to identify the inherent mechanisms of this epithelial barrier dysfunction and its regulatory features., Methods: Epithelial resistance (R epi ) was determined by one-path impedance spectroscopy and 3 H-mannitol fluxes were performed on biopsies from sigmoid colon in miniaturized Ussing chambers. Tight junction proteins were analyzed by Western blot and confocal microscopy. Inflammatory signaling was characterized in HT-29/B6 cells. Apoptosis and mucosal surface parameters were quantified morphologically., Results: R epi was reduced to 53% and 3 H-mannitol fluxes increased 1.7-fold in LC due to lower expression of claudin-4, -5, and -8 and altered subcellular claudin-5 and -8 distributions off the tight junction. TNFα and IFNγ could mimic subcellular redistribution in HT-29/B6 cells, a process which was independent on MLCK activation. Epithelial apoptosis did not contribute to barrier dysfunction in LC and mucosal surface area was unchanged., Conclusions: Epithelial barrier dysfunction in LC occurs through downregulation of claudin-4, -5, and -8, and redistribution of claudin-5 and -8 off the tight junction, which contributes to diarrhea by a leak-flux mechanism. The key effector cytokines TNFα and IFNγ turned out to be the trigger for redistribution of claudin-5 and -8. Thus, alongside sodium malabsorption, leak-flux is yet another important diarrheal mechanism in LC.

Colorectal cancer (CRC) arising in Lynch syndrome (LS) comprises tumours with constitutional mutations in DNA mismatch repair genes. There is still a lack of whole-genome and transcriptome studies of LS-CRC to address questions about similarities and differences in mutation and gene expression characteristics between LS-CRC and sporadic CRC, about the molecular heterogeneity of LS-CRC, and about specific mechanisms of LS-CRC genesis linked to dysfunctional mismatch repair in LS colonic mucosa and the possible role of immune editing. Here, we provide a first molecular characterization of LS tumours and of matched tumour-distant reference colonic mucosa based on whole-genome DNA-sequencing and RNA-sequencing analyses. Our data support two subgroups of LS-CRCs, G1 and G2, whereby G1 tumours show a higher number of somatic mutations, a higher amount of microsatellite slippage, and a different mutation spectrum. The gene expression phenotypes support this difference. Reference mucosa of G1 shows a strong immune response associated with the expression of HLA and immune checkpoint genes and the invasion of CD4+ T cells. Such an immune response is not observed in LS tumours, G2 reference and normal (non-Lynch) mucosa, and sporadic CRC. We hypothesize that G1 tumours are edited for escape from a highly immunogenic microenvironment via loss of HLA presentation and T-cell exhaustion. In contrast, G2 tumours seem to develop in a less immunogenic microenvironment where tumour-promoting inflammation parallels tumourigenesis. Larger studies on non-neoplastic mucosa tissue of mutation carriers are required to better understand the early phases of emerging tumours. Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd., (Copyright © 2017 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.)

The plasticity of myeloid cells is illustrated by a diversity of functions including their role as effectors of innate immunity as macrophages (MACs) and bone remodelling as osteoclasts (OCs). TET2, a methylcytosine dioxygenase highly expressed in these cells and frequently mutated in myeloid leukemias, may be a key contributor to this plasticity. Through transcriptomic and epigenomic analyses, we investigated 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC) and gene expression changes in two divergent terminal myeloid differentiation processes, namely MAC and OC differentiation. MACs and OCs undergo highly similar 5hmC and 5mC changes, despite their wide differences in gene expression. Many TET2- and thymine-DNA glycosylase (TDG)-dependent 5mC and 5hmC changes directly activate the common terminal myeloid differentiation programme. However, the acquisition of differential features between MACs and OCs also depends on TET2/TDG. In fact, 5mC oxidation precedes differential histone modification changes between MACs and OCs. TET2 and TDG downregulation impairs the acquisition of such differential histone modification and expression patterns at MAC-/OC-specific genes. We prove that the histone H3K4 methyltransferase SETD1A is differentially recruited between MACs and OCs in a TET2-dependent manner. We demonstrate a novel role of these enzymes in the establishment of specific elements of identity and function in terminal myeloid differentiation., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The epigenetic writer lysine-specific demethylase 1 (LSD1) is aberrantly upregulated in many cancer types and its overexpression correlates with poor survival and tumor progression. In this study, we analysed LSD1 function in non-small cell lung cancer adenocarcinomas. Expression profiling of 182 cases of lung adenocarcinoma proved a significant correlation of LSD1 overexpression with lung adenocarcinoma progression and metastasis. KRAS-mutated lung cancer cell clones were stably silenced for LSD1 expression. RNA-seq and comprehensive pathway analysis revealed, that genes related to a recently described non-canonical integrin β3 pathway, were significantly downregulated by LSD1 silencing. Hence, invasion and self-renewal capabilities were strongly decreased. Notably, this novel defined LSD1/integrin β3 axis, was also detected in human lung adenocarcinoma specimens. Furthermore, the linkage of LSD1 to an altered expression pattern of lung-lineage specific transcription factors and genes, which are involved in alveolar epithelial differentiation, was demonstrated. Thus, our findings point to a LSD1-integrin β3 axis, conferring attributes of invasiveness and tumor progression to lung adenocarcinoma.

The cellular response to heat stress is an ancient and evolutionarily highly conserved defence mechanism characterised by the transcriptional up-regulation of cyto-protective genes and a partial inhibition of splicing. These features closely resemble the proteotoxic stress response during tumor development. The bromodomain protein BRD4 has been identified as an integral member of the oxidative stress as well as of the inflammatory response, mainly due to its role in the transcriptional regulation process. In addition, there are also several lines of evidence implicating BRD4 in the splicing process. Using RNA-sequencing we found a significant increase in splicing inhibition, in particular intron retentions (IR), following heat treatment in BRD4-depleted cells. This leads to a decrease of mRNA abundancy of the affected transcripts, most likely due to premature termination codons. Subsequent experiments revealed that BRD4 interacts with the heat shock factor 1 (HSF1) such that under heat stress BRD4 is recruited to nuclear stress bodies and non-coding SatIII RNA transcripts are up-regulated. These findings implicate BRD4 as an important regulator of splicing during heat stress. Our data which links BRD4 to the stress induced splicing process may provide novel mechanisms of BRD4 inhibitors in regard to anti-cancer therapies., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Unlabelled: Deletions at the C-terminus of the proto-oncogene protein c-Src kinase are found in the viral oncogene protein v-Src as well as in some advanced human colon cancers. They are associated with increased kinase activity and cellular invasiveness. Here, we analyzed the mRNA expression signature of a constitutively active C-terminal mutant of c-Src, c-Src(mt), in comparison with its wild-type protein, c-Src(wt), in the human non-transformed breast epithelial cell line MCF-10A. We demonstrated previously that the mutant altered migratory and metastatic properties. Genome-wide transcriptome analysis revealed that c-Src(mt) de-regulated the expression levels of approximately 430 mRNAs whose gene products are mainly involved in the cellular processes of migration and adhesion, apoptosis and protein synthesis. 82.9% of these genes have previously been linked to cellular migration, while the others play roles in RNA transport and splicing processes, for instance. Consistent with the transcriptome data, cells expressing c-Src(mt), but not those expressing c-Src(wt), showed the capacity to metastasize into the lungs of mice in vivo. The mRNA expression profile of c-Src(mt)-expressing cells shows significant overlap with that of various primary human tumor samples, possibly reflecting elevated Src activity in some cancerous cells. Expression of c-Src(mt) led to elevated migratory potential. We used this model system to analyze the transcriptional changes associated with an invasive cellular phenotype. These genes and pathways de-regulated by c-Src(mt) may provide suitable biomarkers or targets of therapeutic approaches for metastatic cells., Database: This project was submitted to the National Center for Biotechnology Information BioProject under ID PRJNA288540. The Illumina RNA-Seq reads are available in the National Center for Biotechnology Information Sequence Read Archive under study ID SRP060008 with accession numbers SRS977414 for MCF-10A cells, SRS977717 for mock cells, SRS978053 for c-Src(wt) cells and SRS978046 for c-Src(mt) cells., (© 2016 Federation of European Biochemical Societies.)

In up to 30% of patients with colorectal adenomatous polyposis, no germline mutation in the known genes APC, causing familial adenomatous polyposis, MUTYH, causing MUTYH-associated polyposis, and POLE or POLD1, causing Polymerase-Proofreading-associated polyposis can be identified, although a hereditary etiology is likely. To uncover new causative genes, exome sequencing was performed using DNA from leukocytes and a total of 12 colorectal adenomas from seven unrelated patients with unexplained sporadic adenomatous polyposis. For data analysis and variant filtering, an established bioinformatics pipeline including in-house tools was applied. Variants were filtered for rare truncating point mutations and copy-number variants assuming a dominant, recessive, or tumor suppressor model of inheritance. Subsequently, targeted sequence analysis of the most promising candidate genes was performed in a validation cohort of 191 unrelated patients. All relevant variants were validated by Sanger sequencing. The analysis of exome sequencing data resulted in the identification of rare loss-of-function germline mutations in three promising candidate genes (DSC2, PIEZO1, ZSWIM7). In the validation cohort, further variants predicted to be pathogenic were identified in DSC2 and PIEZO1. According to the somatic mutation spectra, the adenomas in this patient cohort follow the classical pathways of colorectal tumorigenesis. The present study identified three candidate genes which might represent rare causes for a predisposition to colorectal adenoma formation. Especially PIEZO1 (FAM38A) and ZSWIM7 (SWS1) warrant further exploration. To evaluate the clinical relevance of these genes, investigation of larger patient cohorts and functional studies are required.

Background: Lymphocytic colitis (LC) causes watery diarrhea. We aimed to identify mechanisms of altered Na absorption and regulatory inputs in patients with LC by examining the epithelial Na channel (ENaC) function as the predominant Na transport system in human distal colon., Methods: Epithelial Na channel function and regulation was analyzed in biopsies from sigmoid colon of patients with LC and in rat distal colon in Ussing chambers. ENaC-subunit expression was measured by real-time PCR and RNA sequencing. Correction factors for subepithelial resistance contributions were determined by impedance spectroscopy. Upstream regulators in LC were determined by RNA sequencing., Results: Epithelial Na channel-mediated electrogenic Na transport was inhibited despite aldosterone stimulation in human sigmoid colon of patients with LC. The increase in γ-ENaC mRNA expression in response to aldosterone was MEK1/2-dependently reduced in LC, since it could be restored toward normal by MEK1/2 inhibition through U0126. Parallel experiments for identification of signaling in rat distal colon established MEK1/2 to be activated by a cytokine cocktail of TNFα, IFNγ, and IL-15, which were identified as the most important regulators in the upstream regulator analysis in LC., Conclusions: In the sigmoid colon of patients with LC, the key effector cytokines TNFα, IFNγ, and IL-15 inhibited γ-ENaC upregulation in response to aldosterone through a MEK1/2-mediated pathway. This prevents ENaC to reach its maximum transport capacity and results in Na malabsorption which contributes to diarrhea.

Background: In 30-50% of patients with colorectal adenomatous polyposis, no germline mutation in the known genes APC, causing familial adenomatous polyposis, MUTYH, causing MUTYH-associated polyposis, or POLE or POLD1, causing polymerase-proofreading-associated polyposis can be identified, although a hereditary aetiology is likely. This study aimed to explore the impact of APC mutational mosaicism in unexplained polyposis., Methods: To comprehensively screen for somatic low-level APC mosaicism, high-coverage next-generation sequencing of the APC gene was performed using DNA from leucocytes and a total of 53 colorectal tumours from 20 unrelated patients with unexplained sporadic adenomatous polyposis. APC mosaicism was assumed if the same loss-of-function APC mutation was present in ≥ 2 anatomically separated colorectal adenomas/carcinomas per patient. All mutations were validated using diverse methods., Results: In 25% (5/20) of patients, somatic mosaicism of a pathogenic APC mutation was identified as underlying cause of the disease. In 2/5 cases, the mosaic level in leucocyte DNA was slightly below the sensitivity threshold of Sanger sequencing; while in 3/5 cases, the allelic fraction was either very low (0.1-1%) or no mutations were detectable. The majority of mosaic mutations were located outside the somatic mutation cluster region of the gene., Conclusions: The present data indicate a high prevalence of pathogenic mosaic APC mutations below the detection thresholds of routine diagnostics in adenomatous polyposis, even if high-coverage sequencing of leucocyte DNA alone is taken into account. This has important implications for both routine work-up and strategies to identify new causative genes in this patient group., (Published by the BMJ Publishing Group Limited. For permission to use (where not already granted under a licence) please go to http://www.bmj.com/company/products-services/rights-and-licensing/)

Background: The detection of BRAFV600 mutations in patients with metastatic melanoma is important because of the availability of BRAF inhibitor therapy. However, the clinical relevance of the frequency of BRAFV600 mutant alleles is unclear., Patients and Methods: Allele frequencies of BRAFV600 mutations were analyzed by ultra-deep next-generation sequencing in formalin-fixed, paraffin-embedded melanoma tissue (75 primary melanomas and 88 matched metastases). In a second study, pretreatment specimens from 76 patients who received BRAF inhibitors were retrospectively analyzed, and BRAFV600 allele frequencies were correlated with therapeutic results., Results: Thirty-five patients had concordantly BRAF-positive and 36 (48%) patients had concordantly BRAF-negative primary melanomas and matched metastases, and four patients had discordant samples with low allele frequencies (3.4-5.2%). Twenty-six of 35 patients with concordant samples had BRAFV600E mutations, three of whom had additional mutations (V600K in two patients and V600R in one) and nine patients had exclusively non-V600E mutations (V600K in eight patients and V600E -c.1799&#95;1800TG > AA- in one patient). The frequency of mutated BRAFV600 alleles was similar in the primary melanoma and matched metastasis in 27/35 patients, but differed by >3-fold in 8/35 of samples. BRAFV600E allele frequencies in pretreatment tumor specimens were not significantly correlated with treatment outcomes in 76 patients with metastatic melanoma who were treated with BRAF inhibitors., Conclusions: BRAFV600 mutation status and allele frequency is consistent in the majority of primary melanomas and matched metastases. A small subgroup of patients has double mutations. BRAFV600 allele frequencies are not correlated with the response to BRAF inhibitors.

Background: While the song of all songbirds is controlled by the same neural circuit, the hormone dependence of singing behavior varies greatly between species. For this reason, songbirds are ideal organisms to study ultimate and proximate mechanisms of hormone-dependent behavior and neuronal plasticity., Results: We present the high quality assembly and annotation of a female 1.2-Gbp canary genome. Whole genome alignments between the canary and 13 genomes throughout the bird taxa show a much-conserved synteny, whereas at the single-base resolution there are considerable species differences. These differences impact small sequence motifs like transcription factor binding sites such as estrogen response elements and androgen response elements. To relate these species-specific response elements to the hormone-sensitivity of the canary singing behavior, we identify seasonal testosterone-sensitive transcriptomes of major song-related brain regions, HVC and RA, and find the seasonal gene networks related to neuronal differentiation only in the HVC. Testosterone-sensitive up-regulated gene networks of HVC of singing males concerned neuronal differentiation. Among the testosterone-regulated genes of canary HVC, 20% lack estrogen response elements and 4 to 8% lack androgen response elements in orthologous promoters in the zebra finch., Conclusions: The canary genome sequence and complementary expression analysis reveal intra-regional evolutionary changes in a multi-regional neural circuit controlling seasonal singing behavior and identify gene evolution related to the hormone-sensitivity of this seasonal singing behavior. Such genes that are testosterone- and estrogen-sensitive specifically in the canary and that are involved in rewiring of neurons might be crucial for seasonal re-differentiation of HVC underlying seasonal song patterning.

Background: Classic bladder exstrophy (CBE) is the most common form of the bladder exstrophy and epispadias complex. Previously, we and others have identified four patients with a duplication of 22q11.21 among a total of 96 unrelated CBE patients., Methods: Here, we investigated whether this chromosomal aberration was commonly associated with CBE/bladder exstrophy and epispadias complex in an extended case-control sample. Multiplex ligation-dependent probe amplification and microarray-based analysis were used to identify 22q11.21 duplications in 244 unrelated bladder exstrophy and epispadias complex patients (including 217 CBE patients) and 665 healthy controls., Results: New duplications of variable size were identified in four CBE patients and one control. Pooling of our previous and present data (eight duplications in 313 CBE patients) yielded a combined odds ratio of 31.86 (95% confidence interval, 4.24-1407.97). Array-based sequence capture and high-throughput targeted re-sequencing established that all breakpoints resided within the low-copy repeats 22A to 22D. Comparison of the eight duplications revealed a 414 kb phenocritical region harboring 12 validated RefSeq genes. Characterization of these 12 candidate genes through whole-mount in situ hybridization of mouse embryos at embryonic day 9.5 suggested that CRKL, THAP7, and LZTR1 are CBE candidate genes., Conclusion: Our data suggest that duplication of 22q11.21 increases CBE risk and implicate a phenocritical region in disease formation., (© 2014 Wiley Periodicals, Inc.)

MiRNAs are discussed as diagnostic and therapeutic molecules. However, effective miRNA drug treatments with miRNAs are, so far, hampered by the complexity of the miRNA networks. To identify potential miRNA drugs in colorectal cancer, we profiled miRNA and mRNA expression in matching normal, tumor and metastasis tissues of eight patients by Illumina sequencing. We validated six miRNAs in a large tissue screen containing 16 additional tumor entities and identified miRNA-1, miRNA-129, miRNA-497 and miRNA-215 as constantly de-regulated within the majority of cancers. Of these, we investigated miRNA-1 as representative in a systems-biology simulation of cellular cancer models implemented in PyBioS and assessed the effects of depletion as well as overexpression in terms of miRNA-1 as a potential treatment option. In this system, miRNA-1 treatment reverted the disease phenotype with different effectiveness among the patients. Scoring the gene expression changes obtained through mRNA-Seq from the same patients we show that the combination of deep sequencing and systems biological modeling can help to identify patient-specific responses to miRNA treatments. We present this data as guideline for future pre-clinical assessments of new and personalized therapeutic options.

Androgens and oestrogens have been implicated in prostatic carcinogenesis and tumour progression. Although the actions of androgens have been studied extensively, the mechanisms underlying oestrogen signalling in prostate cancer are not fully understood. In the present study, we analyzed the effect of androgens and oestrogens on the expression of anterior gradient 2 (AGR2) and anterior gradient 3 (AGR3), comprising two highly-related genes encoding secretory proteins that are expressed in prostate cancer and one of which (AGR2) has been associated with tumour metastasis. Quantitative reverse-transcriptase PCR and western blot analysis showed androgen induction of AGR2 and AGR3 in three androgen receptor positive cell lines, starting at concentrations of 0.1 nm. Both AGR genes were also transcriptionally activated by ≥ 5 nM oestradiol but not by isotype selective or nonselective oestrogen receptor agonists in DUCaP cells that harbour a high-level of wild-type androgen receptor. A functional androgen receptor but not oestrogen receptor turned out to be required for both androgen and oestrogen regulation. This pattern of androgen and oestrogen regulation was confirmed in VCaP cells and was also observed for FKBP5, a well-characterized androgen-regulated gene. Genome-wide chromatin-immunoprecipitation studies coupled with deep sequencing identified androgen receptor binding sites localized in the distal promoter and intron regions of the AGR2 and AGR3 genes, respectively. The androgen responsiveness of these enhancers was verified by luciferase reporter gene assays and site-directed mutagenesis analysis. Androgen treatment also induced p300 and RNA Pol II recruitment to androgen receptor enhancers of AGR2 and initiated local chromatin remodelling and the formation of RNA Pol II-containing androgen receptor transcription complexes., (© 2013 The Authors Journal compilation © 2013 FEBS.)

Unlabelled: Prostate cancer is the second most common cancer among men worldwide. Alterations in the DNA methylation pattern can be one of the leading causes for prostate cancer formation. This study is the first high-throughput sequencing study investigating genome-wide DNA methylation patterns in a large cohort of 51 tumor and 53 benign prostate samples using methylated DNA immunoprecipitation sequencing. Comparative analyses identified more than 147,000 cancer-associated epigenetic alterations. In addition, global methylation patterns show significant differences based on the TMPRSS2-ERG rearrangement status. We propose the hypermethylation of miR-26a as an alternative pathway of ERG rearrangement-independent EZH2 activation. The observed increase in differential methylation events in fusion-negative tumors can explain the tumorigenic process in the absence of genomic rearrangements., Significance: In contrast to TMPRSS2-ERG -rearranged tumors, the pathomechanism for gene fusion-negative tumors is completely unclear. Using a sequencing-based approach, our work uncovers significant global epigenetic alterations in TMPRSS2-ERG gene fusion-negative tumors and provides a mechanistic explanation for the tumor formation process., (©2012 AACR.)

Background: The heat shock protein 90 (Hsp90) is required for the stability of many signalling kinases. As a target for cancer therapy it allows the simultaneous inhibition of several signalling pathways. However, its inhibition in healthy cells could also lead to severe side effects. This is the first comprehensive analysis of the response to Hsp90 inhibition at the kinome level., Methods: We quantitatively profiled the effects of Hsp90 inhibition by geldanamycin on the kinome of one primary (Hs68) and three tumour cell lines (SW480, U2OS, A549) by affinity proteomics based on immobilized broad spectrum kinase inhibitors ("kinobeads"). To identify affected pathways we used the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway classification. We combined Hsp90 and proteasome inhibition to identify Hsp90 substrates in Hs68 and SW480 cells. The mutational status of kinases from the used cell lines was determined using next-generation sequencing. A mutation of Hsp90 candidate client RIPK2 was mapped onto its structure., Results: We measured relative abundances of > 140 protein kinases from the four cell lines in response to geldanamycin treatment and identified many new potential Hsp90 substrates. These kinases represent diverse families and cellular functions, with a strong representation of pathways involved in tumour progression like the BMP, MAPK and TGF-beta signalling cascades. Co-treatment with the proteasome inhibitor MG132 enabled us to classify 64 kinases as true Hsp90 clients. Finally, mutations in 7 kinases correlate with an altered response to Hsp90 inhibition. Structural modelling of the candidate client RIPK2 suggests an impact of the mutation on a proposed Hsp90 binding domain., Conclusions: We propose a high confidence list of Hsp90 kinase clients, which provides new opportunities for targeted and combinatorial cancer treatment and diagnostic applications.

Background: Massively parallel sequencing technologies have brought an enormous increase in sequencing throughput. However, these technologies need to be further improved with regard to reproducibility and applicability to clinical samples and settings., Methods: Using identification of genetic variations in prostate cancer as an example we address three crucial challenges in the field of targeted re-sequencing: Small nucleotide variation (SNV) detection in samples of formalin-fixed paraffin embedded (FFPE) tissue material, minimal amount of input sample and sampling in view of tissue heterogeneity., Results: We show that FFPE tissue material can supplement for fresh frozen tissues for the detection of SNVs and that solution-based enrichment experiments can be accomplished with small amounts of DNA with only minimal effects on enrichment uniformity and data variance.Finally, we address the question whether the heterogeneity of a tumor is reflected by different genetic alterations, e.g. different foci of a tumor display different genomic patterns. We show that the tumor heterogeneity plays an important role for the detection of copy number variations., Conclusions: The application of high throughput sequencing technologies in cancer genomics opens up a new dimension for the identification of disease mechanisms. In particular the ability to use small amounts of FFPE samples available from surgical tumor resections and histopathological examinations facilitates the collection of precious tissue materials. However, care needs to be taken in regard to the locations of the biopsies, which can have an influence on the prediction of copy number variations. Bearing these technological challenges in mind will significantly improve many large-scale sequencing studies and will - in the long term - result in a more reliable prediction of individual cancer therapies.

Aims: A shift in primary carbon metabolism is the fastest response to oxidative stress. Induced within seconds, it precedes transcriptional regulation, and produces reducing equivalents in form of NADPH within the pentose phosphate pathway (PPP)., Results: Here, we provide evidence for a regulatory signaling function of this metabolic transition in yeast. Several PPP-deficiencies caused abnormal accumulation of intermediate metabolites during the stress response. These PPP-deficient strains had strong growth deficits on media containing oxidants, but we observed that part of their oxidant-phenotypes were not attributable to the production of NADPH equivalents. This pointed to a second, yet unknown role of the PPP in the antioxidant response. Comparing transcriptome profiles obtained by RNA sequencing, we found gene expression profiles that resembled oxidative conditions when PPP activity was increased. Vice versa, deletion of PPP enzymes disturbed and delayed mRNA and protein expression during the antioxidant response., Innovation: Thus, the transient activation of the PPP is a metabolic signal required for balancing and timing gene expression upon an oxidative burst., Conclusion: Consequently, dynamic rearrangements in central carbon metabolism seem to be of major importance for eukaryotic redox sensing, and represent a novel class of dynamic gene expression regulators.

Alternative pre-mRNA splicing is regarded as a pivotal mechanism for generating proteome diversity and complexity from a limited inventory of mammalian genes. Aberrant splicing has been described as a predisposing factor for a number of diseases, but very little is known about its role in chronic inflammation. In this study, we systematically screened 149 splicing factors and 145 potential intron retention events for occurrence and differential expression in inflammatory bowel diseases (IBD). As a result, we identified 47 splicing factors and 33 intron retention events that were differentially regulated in mucosal tissue of IBD patients at transcript level. Despite the fact that Crohn's disease and ulcerative colitis, two subtypes of IBD, share the expression patterns of splicing factors and intron retention events in the majority of cases, we observed significant differences. To investigate these subtype-specific changes in detail we determined the expression levels of seven splicing factors (DUSP11, HNRPAB, HNRPH3, SLU7, SFR2IP, SFPQ, SF3B14) and three intron retention events (PARC, IER3, FGD2) in a cohort of 165 patients with inflammatory diseases of the colon (120 with IBD) and 30 healthy controls by real time PCR (TaqMan). This study demonstrates the potential impact of regulated splicing factors on subsequent regulated intron retention in the pathogenesis of chronic inflammation, exemplified by IBD., (Copyright © 2010 Elsevier GmbH. All rights reserved.)

The development of cancer is characterized by the joined occurrence of alterations on different levels--from single nucleotide changes via structural and copy number variations to epigenetic alterations. With the advent of advanced technologies such as next generation sequencing, we have now the tools in hands to put some light on complex processes and recognize systematic patterns that develop throughout cancer progression. The combination of single hypothesis-driven experiments with a system-wide genetic view enables us to prove so far not addressable questions such as the influence of DNA methylation on gene expression or the disruption of genome homeostasis by structural variations and miRNA expression patterns. Out of this enormous amount of information, specific biomarkers for cancer progression have been discovered, which pave the way for the development of new therapeutic strategies. Here, we will review the status quo of integrative cancer genomic approaches, give an overview over the power of next generation sequencing technologies in oncology, and outline future perspective. Both sides--clinical as well as basic research aspects--will be considered., (© Springer Science+Business Media, LLC 2011)