3D-cardiomics: A spatial transcriptional atlas of the mammalian heart.

Understanding the spatial gene expression and regulation in the heart is key to uncovering its developmental and physiological processes, during homeostasis and disease. Numerous techniques exist to gain gene expression and regulation information in organs such as the heart, but few utilize intuitive true-to-life three-dimensional representations to analyze and visualise results. Here we combined transcriptomics with 3D-modelling to interrogate spatial gene expression in the mammalian heart. For this, we microdissected and sequenced transcriptome-wide 18 anatomical sections of the adult mouse heart. Our study has unveiled known and novel genes that display complex spatial expression in the heart sub-compartments. We have also created 3D-cardiomics, an interface for spatial transcriptome analysis and visualization that allows the easy exploration of these data in a 3D model of the heart. 3D-cardiomics is accessible from http://3d-cardiomics.erc.monash.edu/.

MonaGO: a novel gene ontology enrichment analysis visualisation system.

BACKGROUND: Gene ontology (GO) enrichment analysis is frequently undertaken during exploration of various -omics data sets. Despite the wide array of tools available to biologists to perform this analysis, meaningful visualisation of the overrepresented GO in a manner which is easy to interpret is still lacking. RESULTS: Monash Gene Ontology (MonaGO) is a novel web-based visualisation system that provides an intuitive, interactive and responsive interface for performing GO enrichment analysis and visualising the results. MonaGO supports gene lists as well as GO terms as inputs. Visualisation results can be exported as high-resolution images or restored in new sessions, allowing reproducibility of the analysis. An extensive comparison between MonaGO and 11 state-of-the-art GO enrichment visualisation tools based on 9 features revealed that MonaGO is a unique platform that simultaneously allows interactive visualisation within one single output page, directly accessible through a web browser with customisable display options. CONCLUSION: MonaGO combines dynamic clustering and interactive visualisation as well as customisation options to assist biologists in obtaining meaningful representation of overrepresented GO terms, producing simplified outputs in an unbiased manner. MonaGO will facilitate the interpretation of GO analysis and will assist the biologists into the representation of the results.

Enduring epigenetic landmarks define the cancer microenvironment.

The growth and progression of solid tumors involves dynamic cross-talk between cancer epithelium and the surrounding microenvironment. To date, molecular profiling has largely been restricted to the epithelial component of tumors; therefore, features underpinning the persistent protumorigenic phenotype of the tumor microenvironment are unknown. Using whole-genome bisulfite sequencing, we show for the first time that cancer-associated fibroblasts (CAFs) from localized prostate cancer display remarkably distinct and enduring genome-wide changes in DNA methylation, significantly at enhancers and promoters, compared to nonmalignant prostate fibroblasts (NPFs). Differentially methylated regions associated with changes in gene expression have cancer-related functions and accurately distinguish CAFs from NPFs. Remarkably, a subset of changes is shared with prostate cancer epithelial cells, revealing the new concept of tumor-specific epigenome modifications in the tumor and its microenvironment. The distinct methylome of CAFs provides a novel epigenetic hallmark of the cancer microenvironment and promises new biomarkers to improve interpretation of diagnostic samples.

Algorithm for calculating high disease activity in SLE.

BACKGROUND: The ability to identify lupus patients in high disease activity status (HDAS) without knowledge of the SLEDAI could have application in selection of patients for treatment escalation or enrolment in trials. We sought to generate an algorithm that could calculate via model fitting the presence of HDAS using simple demographic and laboratory values. METHODS: We examined the association of high disease activity (HDA) with demographic and laboratory parameters using prospectively collected data. An HDA visit is recorded when SLEDAI-2K ≥10. We utilized the use of combinatorial search to find algorithms to build a mathematical model predictive of HDA. Performance of each algorithm was evaluated using multi-class area under the receiver operating characteristic curve and the final model was compared with the naïve Bayes classifier, and analysed using the confusion matrix for accuracy and misclassification rate. RESULTS: Data on 286 patients, followed for a median of 5.1 years were studied for a total of 5680 visits. Sixteen laboratory parameters were found to be significantly associated with HDA. A total of 216 algorithms were evaluated and the final algorithm chosen was based on seven pathology measures and three demographic variables. It has an accuracy of 88.6% and misclassification rate of 11.4%. When compared with the naïve Bayes classifier [area under the curve (AUC) = 0.663], our algorithm has a better accuracy with AUC = 0.829. CONCLUSION: This study shows that building an accurate model to calculate HDA using routinely available clinical parameters is feasible. Future studies to independently validate the algorithm will be needed to confirm its predictive performance.

View from the heart: cardiac fibroblasts in development, scarring and regeneration.

In the adult, tissue repair after injury is generally compromised by fibrosis, which maintains tissue integrity with scar formation but does not restore normal architecture and function. The process of regeneration is necessary to replace the scar and rebuild normal functioning tissue. Here, we address this problem in the context of heart disease, and discuss the origins and characteristics of cardiac fibroblasts, as well as the crucial role that they play in cardiac development and disease. We discuss the dual nature of cardiac fibroblasts, which can lead to scarring, pathological remodelling and functional deficit, but can also promote heart function in some contexts. Finally, we review current and proposed approaches whereby regeneration could be fostered by interventions that limit scar formation.

Effect of storage duration on cytokine stability in human serum and plasma.

Quantification of analytes such as cytokines in serum samples is intrinsic to translational research in immune diseases. Optimising pre-analytical conditions is critical for ensuring study quality, including evaluation of cytokine stability. We aimed to evaluate the effect on cytokine stability of storage duration prior to freezing of serum, and compare to plasma samples obtained from patients with systemic lupus erythematosus (SLE). Protein stability was analysed by simultaneously quantifying 18 analytes using a custom multi-analyte profile in SLE patient serum and plasma samples that had been prospectively stored at 4 °C for pre-determined periods between 0 and 30 days, prior to freezing. Six analytes were excluded from analysis, because most tested samples were above or below the limit of detection. Amongst the 12 analysed proteins, 11 did not show significant signal degradation. Significant signal degradation was observed from the fourth day of storage for a single analyte, CCL19. Proteins levels were more stable in unseparated serum compared to plasma for most analytes, with the exception of IL-37 which appears slightly more stable in plasma. Based on this, a maximum 3 days of storage at 4 °C for unseparated serum samples is recommended for biobanked samples intended for cytokine analysis in studies of human immune disease.

TrawlerWeb: an online de novo motif discovery tool for next-generation sequencing datasets.

BACKGROUND: A strong focus of the post-genomic era is mining of the non-coding regulatory genome in order to unravel the function of regulatory elements that coordinate gene expression (Nat 489:57-74, 2012; Nat 507:462-70, 2014; Nat 507:455-61, 2014; Nat 518:317-30, 2015). Whole-genome approaches based on next-generation sequencing (NGS) have provided insight into the genomic location of regulatory elements throughout different cell types, organs and organisms. These technologies are now widespread and commonly used in laboratories from various fields of research. This highlights the need for fast and user-friendly software tools dedicated to extracting cis-regulatory information contained in these regulatory regions; for instance transcription factor binding site (TFBS) composition. Ideally, such tools should not require prior programming knowledge to ensure they are accessible for all users. RESULTS: We present TrawlerWeb, a web-based version of the Trawler_standalone tool (Nat Methods 4:563-5, 2007; Nat Protoc 5:323-34, 2010), to allow for the identification of enriched motifs in DNA sequences obtained from next-generation sequencing experiments in order to predict their TFBS composition. TrawlerWeb is designed for online queries with standard options common to web-based motif discovery tools. In addition, TrawlerWeb provides three unique new features: 1) TrawlerWeb allows the input of BED files directly generated from NGS experiments, 2) it automatically generates an input-matched biologically relevant background, and 3) it displays resulting conservation scores for each instance of the motif found in the input sequences, which assists the researcher in prioritising the motifs to validate experimentally. Finally, to date, this web-based version of Trawler_standalone remains the fastest online de novo motif discovery tool compared to other popular web-based software, while generating predictions with high accuracy. CONCLUSIONS: TrawlerWeb provides users with a fast, simple and easy-to-use web interface for de novo motif discovery. This will assist in rapidly analysing NGS datasets that are now being routinely generated. TrawlerWeb is freely available and accessible at: http://trawler.erc.monash.edu.au .

Adaptation to reef habitats through selection on the coral animal and its associated microbiome.

Spatially adjacent habitats on coral reefs can represent highly distinct environments, often harbouring different coral communities. Yet, certain coral species thrive across divergent environments. It is unknown whether the forces of selection are sufficiently strong to overcome the counteracting effects of the typically high gene flow over short distances, and for local adaptation to occur. We screened the coral genome (using restriction site-associated sequencing) and characterized both the dinoflagellate photosymbiont- and tissue-associated prokaryote microbiomes (using metabarcoding) of a reef flat and slope population of the reef-building coral, Pocillopora damicornis, at two locations on Heron Island in the southern Great Barrier Reef. Reef flat and slope populations were separated by <100 m horizontally and ~5 m vertically, and the two study locations were separated by ~1 km. For the coral host, genetic divergence between habitats was much greater than between locations, suggesting limited gene flow between the flat and slope populations. Consistent with environmental selection, outlier loci primarily belonged to the conserved, minimal cellular stress response, likely reflecting adaptation to the different temperature and irradiance regimes on the reef flat and slope. The prokaryote community differed across both habitat and, to a lesser extent, location, whereas the dinoflagellate photosymbionts differed by habitat but not location. The observed intraspecific diversity associated with divergent habitats supports that environmental adaptation involves multiple members of the coral holobiont. Adaptive alleles or microbial associations present in coral populations from the environmentally variable reef flat may provide a source of adaptive variation for assisted evolution approaches, through assisted gene flow, artificial cross-breeding or probiotic inoculations, with the aim to increase climate resilience in the slope populations.

VISIONET: intuitive visualisation of overlapping transcription factor networks, with applications in cardiogenic gene discovery.

BACKGROUND: Existing de novo software platforms have largely overlooked a valuable resource, the expertise of the intended biologist users. Typical data representations such as long gene lists, or highly dense and overlapping transcription factor networks often hinder biologists from relating these results to their expertise. RESULTS: VISIONET, a streamlined visualisation tool built from experimental needs, enables biologists to transform large and dense overlapping transcription factor networks into sparse human-readable graphs via numerically filtering. The VISIONET interface allows users without a computing background to interactively explore and filter their data, and empowers them to apply their specialist knowledge on far more complex and substantial data sets than is currently possible. Applying VISIONET to the Tbx20-Gata4 transcription factor network led to the discovery and validation of Aldh1a2, an essential developmental gene associated with various important cardiac disorders, as a healthy adult cardiac fibroblast gene co-regulated by cardiogenic transcription factors Gata4 and Tbx20. CONCLUSIONS: We demonstrate with experimental validations the utility of VISIONET for expertise-driven gene discovery that opens new experimental directions that would not otherwise have been identified.

Long-term engrafting multilineage hematopoietic cells differentiated from human induced pluripotent stem cells.

Hematopoietic stem cells (HSCs) derived from human induced pluripotent stem cells (iPS cells) have important biomedical applications. We identified differentiation conditions that generate HSCs defined by robust long-term multilineage engraftment in immune-deficient NOD,B6.Prkdcscid Il2rgtm1Wjl/SzJ KitW41/W41 mice. We guided differentiating iPS cells, as embryoid bodies in a defined culture medium supplemented with retinyl acetate, through HOXA-patterned mesoderm to hemogenic endothelium specified by bone morphogenetic protein 4 and vascular endothelial growth factor (VEGF). Removal of VEGF facilitated an efficient endothelial-to-hematopoietic transition, evidenced by release into the culture medium of CD34+ blood cells, which were cryopreserved. Intravenous transplantation of two million thawed CD34+ cells differentiated from four independent iPS cell lines produced multilineage bone marrow engraftment in 25-50% of immune-deficient recipient mice. These functionally defined, multipotent CD34+ hematopoietic cells, designated iPS cell-derived HSCs (iHSCs), produced levels of engraftment similar to those achieved following umbilical cord blood transplantation. Our study provides a step toward the goal of generating HSCs for clinical translation.

Unraveling the evolutionary origin of the complex Nuclear Receptor Element (cNRE), a cis-regulatory module required for preferential expression in the atrial chamber.

Cardiac function requires appropriate proteins in each chamber. Atria requires slow myosin to act as reservoirs, while ventricles demand fast myosin for swift pumping. Myosins are thus under chamber-biased cis-regulation, with myosin gene expression imbalances leading to congenital heart dysfunction. To identify regulatory inputs leading to cardiac chamber-biased expression, we computationally and molecularly dissected the quail Slow Myosin Heavy Chain III (SMyHC III) promoter that drives preferential expression to the atria. We show that SMyHC III gene states are orchestrated by a complex Nuclear Receptor Element (cNRE) of 32 base pairs. Using transgenesis in zebrafish and mice, we demonstrate that preferential atrial expression is achieved by a combinatorial regulatory input composed of atrial activation motifs and ventricular repression motifs. Using comparative genomics, we show that the cNRE might have emerged from an endogenous viral element through infection of an ancestral host germline, revealing an evolutionary pathway to cardiac chamber-specific expression.

Parallel use of human stem cell lung and heart models provide insights for SARS-CoV-2 treatment.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) primarily infects the respiratory tract, but pulmonary and cardiac complications occur in severe coronavirus disease 2019 (COVID-19). To elucidate molecular mechanisms in the lung and heart, we conducted paired experiments in human stem cell-derived lung alveolar type II (AT2) epithelial cell and cardiac cultures infected with SARS-CoV-2. With CRISPR-Cas9-mediated knockout of ACE2, we demonstrated that angiotensin-converting enzyme 2 (ACE2) was essential for SARS-CoV-2 infection of both cell types but that further processing in lung cells required TMPRSS2, while cardiac cells required the endosomal pathway. Host responses were significantly different; transcriptome profiling and phosphoproteomics responses depended strongly on the cell type. We identified several antiviral compounds with distinct antiviral and toxicity profiles in lung AT2 and cardiac cells, highlighting the importance of using several relevant cell types for evaluation of antiviral drugs. Our data provide new insights into rational drug combinations for effective treatment of a virus that affects multiple organ systems.

Spatially resolved transcriptomics in immersive environments.

Spatially resolved transcriptomics is an emerging class of high-throughput technologies that enable biologists to systematically investigate the expression of genes along with spatial information. Upon data acquisition, one major hurdle is the subsequent interpretation and visualization of the datasets acquired. To address this challenge, VR-Cardiomics is presented, which is a novel data visualization system with interactive functionalities designed to help biologists interpret spatially resolved transcriptomic datasets. By implementing the system in two separate immersive environments, fish tank virtual reality (FTVR) and head-mounted display virtual reality (HMD-VR), biologists can interact with the data in novel ways not previously possible, such as visually exploring the gene expression patterns of an organ, and comparing genes based on their 3D expression profiles. Further, a biologist-driven use-case is presented, in which immersive environments facilitate biologists to explore and compare the heart expression profiles of different genes.

Human pluripotent stem cell-derived macrophages host Mycobacterium abscessus infection.

Human macrophages are a natural host of many mycobacterium species, including Mycobacterium abscessus (M. abscessus), an emerging pathogen affecting immunocompromised and cystic fibrosis patients with few available treatments. The search for an effective treatment is hindered by the lack of a tractable in vitro intracellular infection model. Here, we established a reliable model for M. abscessus infection using human pluripotent stem cell-derived macrophages (hPSC-macrophages). hPSC differentiation permitted reproducible generation of functional macrophages that were highly susceptible to M. abscessus infection. Electron microscopy demonstrated that M. abscessus was present in the hPSC-macrophage vacuoles. RNA sequencing analysis revealed a time-dependent host cell response, with differing gene and protein expression patterns post-infection. Engineered tdTOMATO-expressing hPSC-macrophages with GFP-expressing mycobacteria enabled rapid image-based high-throughput analysis of intracellular infection and quantitative assessment of antibiotic efficacy. Our study describes the first to our knowledge hPSC-based model for M. abscessus infection, representing a novel and accessible system for studying pathogen-host interaction and drug discovery.

CD90 Marks a Mesenchymal Program in Human Thymic Epithelial Cells In Vitro and In Vivo.

Thymic epithelium is critical for the structural integrity of the thymus and for T cell development. Within the fully formed thymus, large numbers of hematopoietic cells shape the thymic epithelium into a scaffold-like structure which bears little similarity to classical epithelial layers, such as those observed in the skin, intestine or pancreas. Here, we show that human thymic epithelial cells (TECs) possess an epithelial identity that also incorporates the expression of mesenchymal cell associated genes, whose expression levels vary between medullary and cortical TECs (m/cTECs). Using pluripotent stem cell (PSC) differentiation systems, we identified a unique population of cells that co-expressed the master TEC transcription factor FOXN1, as well as the epithelial associated marker EPCAM and the mesenchymal associated gene CD90. Using the same serum free culture conditions, we also observed co-expression of EPCAM and CD90 on cultured TECs derived from neonatal human thymus in vitro. Single cell RNA-sequencing revealed these cultured TECs possessed an immature mTEC phenotype and expressed epithelial and mesenchymal associated genes, such as EPCAM, CLDN4, CD90 and COL1A1. Importantly, flow cytometry and single cell RNA-sequencing analysis further confirmed the presence of an EPCAM+CD90+ population in the CD45- fraction of neonatal human thymic stromal cells in vivo. Using the human thymus cell atlas, we found that cTECs displayed more pronounced mesenchymal characteristics than mTECs during embryonic development. Collectively, these results suggest human TECs possess a hybrid gene expression program comprising both epithelial and mesenchymal elements, and provide a basis for the further exploration of thymus development from primary tissues and from the in vitro differentiation of PSCs.

Adult mouse fibroblasts retain organ-specific transcriptomic identity.

Organ fibroblasts are essential components of homeostatic and diseased tissues. They participate in sculpting the extracellular matrix, sensing the microenvironment, and communicating with other resident cells. Recent studies have revealed transcriptomic heterogeneity among fibroblasts within and between organs. To dissect the basis of interorgan heterogeneity, we compare the gene expression of murine fibroblasts from different tissues (tail, skin, lung, liver, heart, kidney, and gonads) and show that they display distinct positional and organ-specific transcriptome signatures that reflect their embryonic origins. We demonstrate that expression of genes typically attributed to the surrounding parenchyma by fibroblasts is established in embryonic development and largely maintained in culture, bioengineered tissues and ectopic transplants. Targeted knockdown of key organ-specific transcription factors affects fibroblast functions, in particular genes involved in the modulation of fibrosis and inflammation. In conclusion, our data reveal that adult fibroblasts maintain an embryonic gene expression signature inherited from their organ of origin, thereby increasing our understanding of adult fibroblast heterogeneity. The knowledge of this tissue-specific gene signature may assist in targeting fibrotic diseases in a more precise, organ-specific manner.

Parallel use of pluripotent human stem cell lung and heart models provide new insights for treatment of SARS-CoV-2.

SARS-CoV-2 primarily infects the respiratory tract, but pulmonary and cardiac complications occur in severe COVID-19. To elucidate molecular mechanisms in the lung and heart, we conducted paired experiments in human stem cell-derived lung alveolar type II (AT2) epithelial cell and cardiac cultures infected with SARS-CoV-2. With CRISPR- Cas9 mediated knock-out of ACE2, we demonstrated that angiotensin converting enzyme 2 (ACE2) was essential for SARS-CoV-2 infection of both cell types but further processing in lung cells required TMPRSS2 while cardiac cells required the endosomal pathway. Host responses were significantly different; transcriptome profiling and phosphoproteomics responses depended strongly on the cell type. We identified several antiviral compounds with distinct antiviral and toxicity profiles in lung AT2 and cardiac cells, highlighting the importance of using several relevant cell types for evaluation of antiviral drugs. Our data provide new insights into rational drug combinations for effective treatment of a virus that affects multiple organ systems. One-sentence summary: Rational treatment strategies for SARS-CoV-2 derived from human PSC models.

Towards spatio-temporally resolved developmental cardiac gene regulatory networks in zebrafish.

Heart formation in the zebrafish involves a rapid, complex series of morphogenetic events in three-dimensional space that spans cardiac lineage specification through to chamber formation and maturation. This process is tightly orchestrated by a cardiac gene regulatory network (GRN), which ensures the precise spatio-temporal deployment of genes critical for heart formation. Alterations of the timing or spatial localisation of gene expression can have a significant impact in cardiac ontogeny and may lead to heart malformations. Hence, a better understanding of the cellular and molecular basis of congenital heart disease relies on understanding the behaviour of cardiac GRNs with precise spatiotemporal resolution. Here, we review the recent technical advances that have expanded our capacity to interrogate the cardiac GRN in zebrafish. In particular, we focus on studies utilising high-throughput technologies to systematically dissect gene expression patterns, both temporally and spatially during heart development.

Disease course following High Disease Activity Status revealed patterns in SLE.

BACKGROUND: We sought to examine the disease course of High Disease Activity Status (HDAS) patients and their different disease patterns in a real-world longitudinal cohort. Disease resolution till Lupus Low Disease Activity State (LLDAS) has been a general treatment goal, but there is limited information on this subset of patients who achieve this. METHODS: All consenting patients of the Monash Lupus Cohort who had at least 12 months of observation were included. HDAS was defined as SLEDAI-2K ≥ 10 ever, and HDAS episode as the period from the first HDAS clinic visit until attainment of LLDAS. We examined the associations of different HDAS patterns with the likelihood of damage accrual. RESULTS: Of 342 SLE patients, 151 experienced HDAS at least once, accounting for 298 HDAS episodes. The majority of HDAS patients (76.2%) experienced Recurrent HDAS (> 1 HDAS visit), and a smaller subset (47.7%) had Persistent HDAS (consecutive HDAS visits for longer than 2 months). Recurrent or Persistent HDAS patients were younger at diagnosis and more likely to experience renal or serositis manifestations; persistent HDAS patients were also more likely to experience neurological manifestations. Baseline SLEDAI greater than 10 was associated with longer HDAS episodes. Recurrent and Persistent HDAS were both associated with an increased likelihood of damage accrual. The total duration of HDAS episode greater than 2 years and experiencing multiple HDAS episodes (≥4) was also associated with an increased likelihood of damage accrual (OR 1.80, 95% CI 1.08-2.97, p = 0.02, and OR 3.31, 95% CI 1.66-13.26, p = 0.01, respectively). CONCLUSION: HDAS episodes have a highly variable course. Recurrent and Persistent HDAS, and longer duration of HDAS episodes, increased the risk of damage accrual. In addition to a major signifier of severity in SLE, its resolution to LLDAS can determine the subsequent outcome in SLE patients.

Glucocorticoid gene signatures in systemic lupus erythematosus and the effects of type I interferon: a cross-sectional and in-vitro study.

BACKGROUND: Glucocorticoids, used as a therapy in systemic lupus erythematosus (SLE), interact with the cytoplasmic glucocorticoid receptor to modulate gene transcription. Minimising the use of glucocorticoids is a goal in SLE; however, pharmacological measures to support clinical guidelines are scarce. We evaluated glucocorticoid-regulated genes for their potential use as biomarkers of glucocorticoid exposure in SLE. We examined interactions between changes in gene expression that are induced by glucocorticoids and type I interferon. METHODS: Genes regulated by glucocorticoids and type I interferon were analysed in relation to glucocorticoid exposure in adult patients meeting the American College of Rheumatology criteria for SLE from three cross-sectional cohorts: a local cohort from a tertiary hospital in Melbourne, VIC, Australia, and two public datasets (GSE49454, Hospital de la Conception, Marseille, France, and GSE88884, patients enrolled in a large, multicentre clinical trial). RNA sequencing was done using RNA from healthy donor leucocytes treated with the glucocorticoid dexamethasone, or type I interferon, or both. FINDINGS: Glucocorticoid-regulated genes were analysed in a local SLE cohort (n=18) and public dataset GSE49454 (n=62). Five genes correlated with glucocorticoid dose in both cohorts and were combined to make a glucocorticoid gene signature. Validity of the glucocorticoid gene signature was tested in the public dataset GSE88884 (n=1756). A dose-dependent association was observed with glucocorticoid dose (p<0·0001), and the glucocorticoid gene signature had moderate ability to identify patients taking high-dose glucocorticoid (area under the curve [AUC]=0·77) although was less discriminatory when including all doses (AUC=0·69). We saw no effect of glucocorticoid dose on type I interferon -regulated gene expression. Patients with a high type I interferon gene signature had reduced glucocorticoid gene signature expression compared with patients with a low type I interferon gene signature matched for glucocorticoid dose, suggesting type I interferon inhibits glucocorticoid-stimulated gene expression. In RNA sequencing experiments, type I interferon impaired the expression of glucocorticoid-induced genes, whereas dexamethasone had minimal effect on the expression of type I interferon-stimulated genes. We identified genes regulated by dexamethasone but not affected by type I interferon; combined signatures using these genes also showed moderate ability to distinguish patients taking glucocorticoids. INTERPRETATION: A gene signature for glucocorticoid exposure was identified, but the substantial effect of type I interferon on glucocorticoid-induced genes might limit its application in SLE. These data confirm the insensitivity of type I interferon-regulated genes to glucocorticoids, and together support the concept that type I interferon has a role in glucocorticoid resistance in SLE. FUNDING: Lupus Research Alliance and Australian National Health and Medical Research Council.