Downregulation of exhausted cytotoxic T cells in gene expression networks of multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) presents with fever, inflammation and pathology of multiple organs in individuals under 21 years of age in the weeks following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Although an autoimmune pathogenesis has been proposed, the genes, pathways and cell types causal to this new disease remain unknown. Here we perform RNA sequencing of blood from patients with MIS-C and controls to find disease-associated genes clustered in a co-expression module annotated to CD56dimCD57+ natural killer (NK) cells and exhausted CD8+ T cells. A similar transcriptome signature is replicated in an independent cohort of Kawasaki disease (KD), the related condition after which MIS-C was initially named. Probing a probabilistic causal network previously constructed from over 1,000 blood transcriptomes both validates the structure of this module and reveals nine key regulators, including TBX21, a central coordinator of exhausted CD8+ T cell differentiation. Together, this unbiased, transcriptome-wide survey implicates downregulation of NK cells and cytotoxic T cell exhaustion in the pathogenesis of MIS-C.

A Proinflammatory Cytokine Network Profile in Th1/Type 1 Effector B Cells Delineates a Common Group of Patients in Four Systemic Autoimmune Diseases

OBJECTIVE

The effector T cell and B cell cytokine networks have been implicated in the pathogenesis of systemic autoimmune diseases, but the association of these cytokine networks with the heterogeneity of clinical manifestations and immune profiles has not been carefully examined. This study was undertaken to examine whether cytokine profiles can delineate distinct groups of patients in 4 systemic autoimmune diseases (systemic lupus erythematosus, Sjögren's syndrome, rheumatoid arthritis, and systemic sclerosis).

METHODS

A total of 179 patients and 48 healthy volunteers were enrolled in the multicenter cross-sectional PRECISE Systemic Autoimmune Diseases (PRECISESADS) study. Multi-low-dimensional omics data (cytokines, autoantibodies, circulating immune cells) were examined. Coculture experiments were performed to test the impact of the cytokine microenvironment on T cell/B cell cross-talk.

RESULTS

A proinflammatory cytokine profile defined by high levels of CXCL10, interleukin-6 (IL-6), IL-2, and tumor necrosis factor characterized a distinct group of patients in the 4 systemic autoimmune diseases. In each disease, this proinflammatory cluster was associated with a specific circulating immune cell signature, more severe disease, and higher levels of autoantibodies, suggesting an uncontrolled proinflammatory Th1 immune response. We observed in vitro that B cells reinforce Th1 differentiation and naive T cell proliferation, leading to the induction of type 1 effector B cells and IgG production. This process was associated with an increase in CXCL10, IL-6, IL-2, and interferon-γ production.

CONCLUSION

This composite analysis brings new insights into human B cell functional heterogeneity based on T cell/B cell cross-talk, and proposes a better stratification of patients with systemic autoimmune diseases, suggesting that combined biomarkers would be of great value for the design of personalized treatments.

Continental-scale genomic analysis suggests shared post-admixture adaptation in the Americas

American populations are one of the most interesting examples of recently admixed groups, where ancestral components from three major continental human groups (Africans, Eurasians and Native Americans) have admixed within the last 15 generations. Recently, several genetic surveys focusing on thousands of individuals shed light on the geography, chronology and relevance of these events. However, even though gene flow could drive adaptive evolution, it is unclear whether and how natural selection acted on the resulting genetic variation in the Americas.

In this study, we analysed the patterns of local ancestry of genomic fragments in genome-wide data for ~ 6000 admixed individuals from 10 American countries. In doing so, we identified regions characterized by a divergent ancestry profile (DAP), in which a significant over or under ancestral representation is evident.

Our results highlighted a series of genomic regions with DAPs associated with immune system response and relevant medical traits, with the longest DAP region encompassing the human leukocyte antigen locus. Furthermore, we found that DAP regions are enriched in genes linked to cancer-related traits and autoimmune diseases. Then, analysing the biological impact of these regions, we showed that natural selection could have acted preferentially towards variants located in coding and non-coding transcripts and characterized by a high deleteriousness score.

Taken together, our analyses suggest that shared patterns of post admixture adaptation occurred at a continental scale in the Americas, affecting more often functional and impactful genomic variants.

A comprehensive database for integrated analysis of omics data in autoimmune diseases

BACKGROUND

Autoimmune diseases are heterogeneous pathologies with difficult diagnosis and few therapeutic options. In the last decade, several omics studies have provided significant insights into the molecular mechanisms of these diseases. Nevertheless, data from different cohorts and pathologies are stored independently in public repositories and a unified resource is imperative to assist researchers in this field.

RESULTS

Here, we present Autoimmune Diseases Explorer ( https://adex.genyo.es ), a database that integrates 82 curated transcriptomics and methylation studies covering 5609 samples for some of the most common autoimmune diseases. The database provides, in an easy-to-use environment, advanced data analysis and statistical methods for exploring omics datasets, including meta-analysis, differential expression or pathway analysis.

CONCLUSIONS

This is the first omics database focused on autoimmune diseases. This resource incorporates homogeneously processed data to facilitate integrative analyses among studies.

A new molecular classification to drive precision treatment strategies in primary Sjögren's syndrome

There is currently no approved treatment for primary Sjögren's syndrome, a disease that primarily affects adult women. The difficulty in developing effective therapies is -in part- because of the heterogeneity in the clinical manifestation and pathophysiology of the disease. Finding common molecular signatures among patient subgroups could improve our understanding of disease etiology, and facilitate the development of targeted therapeutics. Here, we report, in a cross-sectional cohort, a molecular classification scheme for Sjögren's syndrome patients based on the multi-omic profiling of whole blood samples from a European cohort of over 300 patients, and a similar number of age and gender-matched healthy volunteers. Using transcriptomic, genomic, epigenetic, cytokine expression and flow cytometry data, combined with clinical parameters, we identify four groups of patients with distinct patterns of immune dysregulation. The biomarkers we identify can be used by machine learning classifiers to sort future patients into subgroups, allowing the re-evaluation of response to treatments in clinical trials.

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Data processing workflow for large-scale immune monitoring studies by mass cytometry

Mass cytometry is a powerful tool for deep immune monitoring studies. To ensure maximal data quality, a careful experimental and analytical design is required. However even in well-controlled experiments variability caused by either operator or instrument can introduce artifacts that need to be corrected or removed from the data. Here we present a data processing pipeline which ensures the minimization of experimental artifacts and batch effects, while improving data quality. Data preprocessing and quality controls are carried out using an R pipeline and packages like CATALYST for bead-normalization and debarcoding, flowAI and flowCut for signal anomaly cleaning, AOF for files quality control, flowClean and flowDensity for gating, CytoNorm for batch normalization and FlowSOM and UMAP for data exploration. As proper experimental design is key in obtaining good quality events, we also include the sample processing protocol used to generate the data. Both, analysis and experimental pipelines are easy to scale-up, thus the workflow presented here is particularly suitable for large-scale, multicenter, multibatch and retrospective studies.

"Precision Medicine in Autoimmune Diseases: Fact or Fiction"

Much is said about precision medicine, but its real significance and the possibility of making it a real possibility is far from certain. Several studies in each of the autoimmune diseases have provided important insight into molecular pathways but the use of molecular studies, particularly those looking into transcriptome pathways, have seldom approached the possibility of using the data for disease stratification and then for prediction, or diagnosis. Only the type I interferon signature has been considered in the use of this signature for therapeutic purposes, particularly in the case of systemic lupus erythematosus. Here, the authors provide an update on precision medicine, what can be translated into clinical practice, and what do single-cell molecular studies provide to our knowledge in autoimmune diseases, focusing on a few examples. The main message being that we should try to move from precision medicine of established disease to preventive medicine in order to predict the development of disease.

The Role of BANK1 in B Cell Signaling and Disease

The B cell scaffold protein with ankyrin repeats (BANK1) is expressed primarily in B cells and with multiple but discrete roles in B cell signaling, including B cell receptor signaling, CD40-related signaling, and Toll-like receptor (TLR) signaling. The gene for BANK1, located in chromosome 4, has been found to contain genetic variants that are associated with several autoimmune diseases and also other complex phenotypes, in particular, with systemic lupus erythematosus. Common genetic variants are associated with changes in BANK1 expression in B cells, while rare variants modify their capacity to bind efferent effectors during signaling. A BANK1-deficient model has shown the importance of BANK1 during TLR7 and TLR9 signaling and has confirmed its role in the disease. Still, much needs to be done to fully understand the function of BANK1, but the main conclusion is that it may be the link between different signaling functions within the B cells and they may act to synergize the various pathways within a cell. With this review, we hope to enhance the interest in this molecule.

Integrative Analysis Reveals a Molecular Stratification of Systemic Autoimmune Diseases

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.

Age-dependent impact of the major common genetic risk factor for COVID-19 on severity and mortality

BACKGROUND

There is considerable variability in COVID-19 outcomes amongst younger adults-and some of this variation may be due to genetic predisposition. We characterized the clinical implications of the major genetic risk factor for COVID-19 severity, and its age-dependent effect, using individual-level data in a large international multi-centre consortium.

METHOD

The major common COVID-19 genetic risk factor is a chromosome 3 locus, tagged by the marker rs10490770. We combined individual level data for 13,424 COVID-19 positive patients (N=6,689 hospitalized) from 17 cohorts in nine countries to assess the association of this genetic marker with mortality, COVID-19-related complications and laboratory values. We next examined if the magnitude of these associations varied by age and were independent from known clinical COVID-19 risk factors.

FINDINGS

We found that rs10490770 risk allele carriers experienced an increased risk of all-cause mortality (hazard ratio [HR] 1·4, 95% confidence interval [CI] 1·2-1·6) and COVID-19 related mortality (HR 1·5, 95%CI 1·3-1·8). Risk allele carriers had increased odds of several COVID-19 complications: severe respiratory failure (odds ratio [OR] 2·0, 95%CI 1·6-2·6), venous thromboembolism (OR 1·7, 95%CI 1·2-2·4), and hepatic injury (OR 1·6, 95%CI 1·2-2·0). Risk allele carriers ≤ 60 years had higher odds of death or severe respiratory failure (OR 2·6, 95%CI 1·8-3·9) compared to those > 60 years OR 1·5 (95%CI 1·3-1·9, interaction p-value=0·04). Amongst individuals ≤ 60 years who died or experienced severe respiratory COVID-19 outcome, we found that 31·8% (95%CI 27·6-36·2) were risk variant carriers, compared to 13·9% (95%CI 12·6-15·2%) of those not experiencing these outcomes. Prediction of death or severe respiratory failure among those ≤ 60 years improved when including the risk allele (AUC 0·82 vs 0·84, p=0·016) and the prediction ability of rs10490770 risk allele was similar to, or better than, most established clinical risk factors.

INTERPRETATION

The major common COVID-19 risk locus on chromosome 3 is associated with increased risks of morbidity and mortality-and these are more pronounced amongst individuals ≤ 60 years. The effect on COVID-19 severity was similar to, or larger than most established risk factors, suggesting potential implications for clinical risk management.

FUNDING

Funding was obtained by each of the participating cohorts individually.

Stabilization of Human Whole Blood Samples for Multicenter and Retrospective Immunophenotyping Studies

Whole blood is often collected for large-scale immune monitoring studies to track changes in cell frequencies and responses using flow (FC) or mass cytometry (MC). In order to preserve sample composition and phenotype, blood samples should be analyzed within 24 h after bleeding, restricting the recruitment, analysis protocols, as well as biobanking. Herein, we have evaluated two whole blood preservation protocols that allow rapid sample processing and long-term stability. Two fixation buffers were used, Phosphoflow Fix and Lyse (BD) and Proteomic Stabilizer (PROT) to fix and freeze whole blood samples for up to 6 months. After analysis by an 8-plex panel by FC and a 26-plex panel by MC, manual gating of circulating leukocyte populations and cytokines was performed. Additionally, we tested the stability of a single sample over a 13-months period using 45 consecutive aliquots and a 34-plex panel by MC. We observed high correlation and low bias toward any cell population when comparing fresh and 6 months frozen blood with FC and MC. This correlation was confirmed by hierarchical clustering. Low coefficients of variation (CV) across studied time points indicate good sample preservation for up to 6 months. Cytokine detection stability was confirmed by low CVs, with some differences between fresh and fixed conditions. Thirteen months regular follow-up of PROT samples showed remarkable sample stability. Whole blood can be preserved for phenotyping and cytokine-response studies provided the careful selection of a compatible antibody panel. However, possible changes in cell morphology, differences in antibody affinity, and changes in cytokine-positive cell frequencies when compared to fresh blood should be considered. Our setting constitutes a valuable tool for multicentric and retrospective studies. © 2020 International Society for Advancement of Cytometry.

Cytotoxic lymphocytes are dysregulated in multisystem inflammatory syndrome in children

Multisystem inflammatory syndrome in children (MIS-C) presents with fever, inflammation and multiple organ involvement in individuals under 21 years following severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. To identify genes, pathways and cell types driving MIS-C, we sequenced the blood transcriptomes of MIS-C cases, pediatric cases of coronavirus disease 2019, and healthy controls. We define a MIS-C transcriptional signature partially shared with the transcriptional response to SARS-CoV-2 infection and with the signature of Kawasaki disease, a clinically similar condition. By projecting the MIS-C signature onto a co-expression network, we identified disease gene modules and found genes downregulated in MIS-C clustered in a module enriched for the transcriptional signatures of exhausted CD8 + T-cells and CD56 dim CD57 + NK cells. Bayesian network analyses revealed nine key regulators of this module, including TBX21 , a central coordinator of exhausted CD8 + T-cell differentiation. Together, these findings suggest dysregulated cytotoxic lymphocyte response to SARS-Cov-2 infection in MIS-C.

Proteomic analysis in lupus mice identifies Coronin-1A as a potential biomarker for lupus nephritis

BACKGROUND

Approximately 50% of systemic lupus erythematosus (SLE) patients develop nephritis, which is among the most severe and frequent complications of the disease and a leading cause of morbidity and mortality. Despite intensive research, there are still no reliable lupus nephritis (LN) markers in clinical use that can assess renal damage and activity with a high sensitivity and specificity. To this end, the aim of this study was to identify new clinically relevant tissue-specific protein biomarkers and possible underlying molecular mechanisms associated with renal involvement in SLE, using mass spectrometry (MS)-based proteomics.

METHODS

Kidneys were harvested from female triple congenic B6.NZMsle1/sle2/sle3 lupus mice model, and the respective sex- and age-matched C57BL/6 control mice at 12, 24 and 36 weeks of age, representing pre-symptomatic, established and end-stage LN, respectively. Proteins were extracted from kidneys, purified, reduced, alkylated and digested by trypsin. Purified peptides were separated by liquid chromatography and analysed by high-resolution MS. Data were processed by the Progenesis QIp software, and functional annotation analysis was performed using DAVID bioinformatics resources. Immunofluorescence and multiple reaction monitoring (MRM) MS methods were used to confirm prospective biomarkers in SLE mouse strains as well as human serum samples.

RESULTS

Proteomic profiling of kidney tissues from SLE and control mice resulted in the identification of more than 3800 unique proteins. Pathway analysis revealed a number of dysregulated molecular pathways that may be mechanistically involved in renal pathology, including phagosome and proximal tubule bicarbonate reclamation pathways. Proteomic analysis supported by human transcriptomic data and pathway analysis revealed Coronin-1A, Ubiquitin-like protein ISG15, and Rho GDP-dissociation inhibitor 2, as potential LN biomarkers. These results were further validated in other SLE mouse strains using MRM-MS. Most importantly, experiments in humans showed that measurement of Coronin-1A in human sera using MRM-MS can segregate LN patients from SLE patients without nephritis with a high sensitivity (100%) and specificity (100%).

CONCLUSIONS

These preliminary findings suggest that serum Coronin-1A may serve as a promising non-invasive biomarker for LN and, upon validation in larger cohorts, may be employed in the future as a screening test for renal disease in SLE patients.

Key steps and methods in the experimental design and data analysis of highly multi-parametric flow and mass cytometry

High-dimensional, single-cell cell technologies revolutionized the way to study biological systems, and polychromatic flow cytometry (FC) and mass cytometry (MC) are two of the drivers of this revolution. As up to 30-50 dimensions respectively can be measured per single-cell, they allow deep phenotyping combined with cellular functions studies, like cytokine production or protein phosphorylation. In parallel, the bioinformatics field develops algorithms that are able to process incoming data and extract the most useful and meaningful biological information. However, the success of automated analysis tools depends on the generation of high-quality data. In this review we present the most recent FC and MC computational approaches that are used to prepare, process and interpret high-content cytometry data. We also underscore proper experimental design as a key step for obtaining good quality data.

A survey of gene expression meta-analysis: methods and applications

The increasing use of high-throughput gene expression quantification technologies over the last two decades and the fact that most of the published studies are stored in public databases has triggered an explosion of studies available through public repositories. All this information offers an invaluable resource for reuse to generate new knowledge and scientific findings. In this context, great interest has been focused on meta-analysis methods to integrate and jointly analyze different gene expression datasets. In this work, we describe the main steps in the gene expression meta-analysis, from data preparation to the state-of-the art statistical methods. We also analyze the main types of applications and problems that can be approached in gene expression meta-analysis studies and provide a comparative overview of the available software and bioinformatics tools. Moreover, a practical guide for choosing the most appropriate method in each case is also provided.

ImaGEO: integrative gene expression meta-analysis from GEO database

SUMMARY

The Gene Expression Omnibus (GEO) database provides an invaluable resource of publicly available gene expression data that can be integrated and analyzed to derive new hypothesis and knowledge. In this context, gene expression meta-analysis (geMAs) is increasingly used in several fields to improve study reproducibility and discovering robust biomarkers. Nevertheless, integrating data is not straightforward without bioinformatics expertise. Here, we present ImaGEO, a web tool for geMAs that implements a complete and comprehensive meta-analysis workflow starting from GEO dataset identifiers. The application integrates GEO datasets, applies different meta-analysis techniques and provides functional analysis results in an easy-to-use environment. ImaGEO is a powerful and useful resource that allows researchers to integrate and perform meta-analysis of GEO datasets to lead robust findings for biomarker discovery studies.

AVAILABILITY AND IMPLEMENTATION

ImaGEO is accessible at http://bioinfo.genyo.es/imageo/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Discovering new metabolite alterations in primary sjögren's syndrome in urinary and plasma samples using an HPLC-ESI-QTOF-MS methodology

Sjögren's Syndrome (SjS) is a complex autoimmune disease characterized by the affection of the exocrine glands and the involvement of multiple organs. Although a greater number of biomarker studies have been carried out in recent years, the origin and pathogenesis are not yet well known and therefore there is a need to continue studying this pathology. This work aims to find metabolic changes in biological samples (plasma and urine), which could help identify the metabolic pathways affected by the SjS pathogenesis. The samples collected from SjS patients and healthy volunteers were analyzed by a fingerprinting metabolomic approach based on HPLC-ESI-QTOF-MS methodology. After feature pre-selection by univariate statistical tests, an integrated PLS-DA model using data from urine and plasma was constructed obtaining a good classification between cases and controls (AUROC = 0.839 ± 0.021). 31 and 38 metabolites in plasma and urine, respectively, showed significant differences between healthy volunteers and SjS patients and were proposed for their identification. From them, 12 plasma and 24 urinary metabolites could be annotated. In general, the main metabolic pathways altered in SjS patients were related to the metabolism of phospholipids, fatty acids, and amino acids, specially tryptophan, proline and phenylalanine.

The Genomic Impact of European Colonization of the Americas

The human genetic diversity of the Americas has been affected by several events of gene flow that have continued since the colonial era and the Atlantic slave trade. Moreover, multiple waves of migration followed by local admixture occurred in the last two centuries, the impact of which has been largely unexplored. Here, we compiled a genome-wide dataset of ∼12,000 individuals from twelve American countries and ∼6,000 individuals from worldwide populations and applied haplotype-based methods to investigate how historical movements from outside the New World affected (1) the genetic structure, (2) the admixture profile, (3) the demographic history, and (4) sex-biased gene-flow dynamics of the Americas. We revealed a high degree of complexity underlying the genetic contribution of European and African populations in North and South America, from both geographic and temporal perspectives, identifying previously unreported sources related to Italy, the Middle East, and to specific regions of Africa.

Differential Treatments Based on Drug-induced Gene Expression Signatures and Longitudinal Systemic Lupus Erythematosus Stratification

Systemic lupus erythematosus (SLE) is a heterogeneous disease with unpredictable patterns of activity. Patients with similar activity levels may have different prognosis and molecular abnormalities. In this study, we aimed to measure the main differences in drug-induced gene expression signatures across SLE patients and to evaluate the potential for clinical data to build a machine learning classifier able to predict the SLE subset for individual patients. SLE transcriptomic data from two cohorts were compared with drug-induced gene signatures from the CLUE database to compute a connectivity score that reflects the capability of a drug to revert the patient signatures. Patient stratification based on drug connectivity scores revealed robust clusters of SLE patients identical to the clusters previously obtained through longitudinal gene expression data, implying that differential treatment depends on the cluster to which patients belongs. The best drug candidates found, mTOR inhibitors or those reducing oxidative stress, showed stronger cluster specificity. We report that drug patterns for reverting disease gene expression follow the cell-specificity of the disease clusters. We used 2 cohorts to train and test a logistic regression model that we employed to classify patients from 3 independent cohorts into the SLE subsets and provide a clinically useful model to predict subset assignment and drug efficacy.

Epigenome-Wide Comparative Study Reveals Key Differences Between Mixed Connective Tissue Disease and Related Systemic Autoimmune Diseases

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.

The TREX1 Dinosaur Bites the Brain through the LINE

In this issue of Cell Stem Cell, Thomas et al. (2017) define the nature of accumulated ssDNA present in the neuron and astrocyte cytoplasm of TREX1 mutated stem cell-derived organoids. Accumulated ssDNAs are derived from LINE-1 endogenous retroelements, providing new clues as to the development of Aicardi-Goutières syndrome in the neural system.

Exploring Impact of Rare Variation in Systemic Lupus Erythematosus by a Genome Wide Imputation Approach

The importance of low frequency and rare variation in complex disease genetics is difficult to estimate in patient populations. Genome-wide association studies are therefore, underpowered to detect rare variation. We have used a combined approach of genome-wide-based imputation with a highly stringent sequence kernel association (SKAT) test and a case-control burden test. We identified 98 candidate genes containing rare variation that in aggregate show association with SLE many of which have recognized immunological function, but also function and expression related to relevant tissues such as the joints, skin, blood or central nervous system. In addition we also find that there is a significant enrichment of genes annotated for disease-causing mutations in the OMIM database, suggesting that in complex diseases such as SLE, such mutations may be involved in subtle or combined phenotypes or could accelerate specific organ abnormalities found in the disease. We here provide an important resource of candidate genes for SLE.