Systemic lupus erythematosus in American Indian/Alaska natives: Incorporating our new understanding of the biology of trauma

OBJECTIVE

To review the literature regarding systemic lupus erythematosus (SLE) in American Indian/Alaska Native (AI/AN) people and relate prevalence and/or disease severity to our emerging understanding of the biology of trauma and toxic stress.

METHODS

We conducted a search and review of the literature using search terms "lupus and American Indians" "ACEs and disease outcome" "Biology of Adversity" "lupus and ACE scores," " lupus and childhood abuse." These search criteria were entered into Google Scholar and articles retrieved from PubMed, NBCI. This approach yielded a small numbers of papers used throughout this review. We excluded articles that were not published in a peer reviewed journals, as well as editorial commentaries.

RESULTS

In the AI/AN population, SLE shows high prevalence rates and severe disease manifestations, comparable to the African American population. AI/AN populations also have high rates of childhood trauma. Toxic stress and trauma such as those catalogued in the Adverse Childhood Experiences (ACE) study have broad-reaching immunologic and epigenetic effects that are likely to be relevant to our understanding of SLE in AI/AN people.

CONCLUSIONS

AI/AN people have high rates of SLE. These high rates are likely to be driven by many complex factors, not all of which are genetic. Future research is needed to establish (or refute) a causal connection between the biology of adversity and SLE in socially marginalized and historically traumatized populations.

Association of anti-TPM4 autoantibodies with vasculopathic cutaneous manifestations in juvenile dermatomyositis

OBJECTIVES

AECAs are detected in multiple forms of vasculitis or vasculopathy, including JDM. High levels of tropomyosin alpha-4 chain (TPM4) gene expression in cutaneous lesions and TPM4 protein expression in some endothelial cells (ECs) have been proven. Furthermore, the presence of autoantibodies to tropomyosin proteins have been discovered in DM. We therefore investigated whether anti-TPM4 autoantibodies are an AECA in JDM and are correlated with clinical features of JDM.

METHODS

The expression of TPM4 protein in cultured normal human dermal microvascular ECs was investigated by Western blotting. Plasma samples from 63 children with JDM, 50 children with polyarticular JIA (pJIA) and 40 healthy children (HC) were tested for the presence of anti-TPM4 autoantibodies using an ELISA. Clinical features were compared between JDM patients with and without anti-TPM4 autoantibodies.

RESULTS

Autoantibodies to TPM4 were detected in the plasma of 30% of JDM, 2% of pJIA (P < 0.0001) and 0% of HC (P < 0.0001). In JDM, anti-TPM4 autoantibodies were associated with the presence of cutaneous ulcers (53%; P = 0.02), shawl sign rash (47%; P = 0.03), mucous membrane lesions (84%; P = 0.04) and subcutaneous edema (42%; P < 0.05). Anti-TPM4 autoantibodies significantly correlated with the use of intravenous steroids and IVIG therapy in JDM (both P = 0.01). The total number of medications received was higher in patients with anti-TPM4 autoantibodies (P = 0.02).

CONCLUSION

Anti-TPM4 autoantibodies are detected frequently in children with JDM and are novel myositis-associated autoantibodies. Their presence correlates with vasculopathic and other cutaneous manifestations of JDM that may be indicative of more refractory disease.

Comparison of the three-dimensional chromatin structures of adolescent and adult peripheral blood B cells: implications for the study of pediatric autoimmune diseases

Background/Purpose

Knowledge of the 3D genome is essential to elucidate genetic mechanisms driving autoimmune diseases. The 3D genome is distinct for each cell type, and it is uncertain whether cell lines faithfully recapitulate the 3D architecture of primary human cells or whether developmental aspects of the pediatric immune system require use of pediatric samples. We undertook a systematic analysis of B cells and B cell lines to compare 3D genomic features encompassing risk loci for juvenile idiopathic arthritis (JIA), systemic lupus (SLE), and type 1 diabetes (T1D).

Methods

We isolated B cells from healthy individuals, ages 9-17. HiChIP was performed using CTCF antibody, and CTCF peaks were identified. CTCF loops within the pediatric were compared to three datasets: 1) self-called CTCF consensus peaks called within the pediatric samples, 2) ENCODE's publicly available GM12878 CTCF ChIP-seq peaks, and 3) ENCODE's primary B cell CTCF ChIPseq peaks from two adult females. Differential looping was assessed within the pediatric samples and each of the three peak datasets.

Results

The number of consensus peaks called in the pediatric samples was similar to that identified in ENCODE's GM12878 and primary B cell datasets. We observed <1% of loops that demonstrated significantly differential looping between peaks called within the pediatric samples themselves and when called using ENCODE GM12878 peaks . Significant looping differences were even less when comparing loops of the pediatric called peaks to those of the ENCODE primary B cell peaks. When querying loops found in juvenile idiopathic arthritis, type 1 diabetes, or systemic lupus erythematosus risk haplotypes, we observed significant differences in only 2.2%, 1.0%, and 1.3% loops, respectively, when comparing peaks called within the pediatric samples and ENCODE GM12878 dataset. The differences were even less apparent when comparing loops called with the pediatric vs ENCODE adult primary B cell peak datasets.The 3D chromatin architecture in B cells is similar across pediatric, adult, and EBVtransformed cell lines. This conservation of 3D structure includes regions encompassing autoimmune risk haplotypes.

Conclusion

Thus, even for pediatric autoimmune diseases, publicly available adult B cell and cell line datasets may be sufficient for assessing effects exerted in the 3D genomic space.

Analysis of chromatin data supports a role for CD14+ monocytes/macrophages in mediating genetic risk for juvenile idiopathic arthritis

Introduction

Genome wide association studies (GWAS) have identified multiple regions that confer genetic risk for the polyarticular/oligoarticular forms of juvenile idiopathic arthritis (JIA). However, genome-wide scans do not identify the cells impacted by genetic polymorphisms on the risk haplotypes or the genes impacted by those variants. We have shown that genetic variants driving JIA risk are likely to affect both innate and adaptive immune functions. We provide additional evidence that JIA risk variants impact innate immunity.

Materials and methods

We queried publicly available H3K4me1/H3K27ac ChIP-seq data in CD14+ monocytes to determine whether the linkage disequilibrium (LD) blocks incorporating the SNPs that tag JIA risk loci showed enrichment for these epigenetic marks. We also queried monocyte/macrophage GROseq data, a functional readout of active enhancers. We defined the topologically associated domains (TADs) encompassing enhancers on the risk haplotypes and identified genes within those TADs expressed in monocytes. We performed ontology analyses of these genes to identify cellular processes that may be impacted by these variants. We also used whole blood RNAseq data from the Genotype-Tissue Expression (GTEx) data base to determine whether SNPs lying within monocyte GROseq peaks influence plausible target genes within the TADs encompassing the JIA risk haplotypes.

Results

The LD blocks encompassing the JIA genetic risk regions were enriched for H3K4me1/H3K27ac ChIPseq peaks (p=0.00021 and p=0.022) when compared to genome background. Eleven and sixteen JIA were enriched for resting and activated macrophage GROseq peaks, respectively risk regions (p=0.04385 and p=0.00004). We identified 321 expressed genes within the TADs encompassing the JIA haplotypes in human monocytes. Ontological analysis of these genes showed enrichment for multiple immune functions. Finally, we found that SNPs lying within the GROseq peaks are strongly associated with expression levels of plausible target genes in human whole blood.

Conclusions

These findings support the idea that both innate and adaptive immunity are impacted by JIA genetic risk variants.

Health inequities in the rheumatic diseases of childhood

PURPOSE OF REVIEW

To describe differences in disease manifestations and outcomes in pediatric rheumatic diseases as they occur in non-European-descended populations in North America.

RECENT FINDINGS

Differences in disease prevalence, clinical phenotypes, disease course, and outcomes have been described across the spectrum of pediatric-onset rheumatic diseases. Although these differences are commonly explained by differences in genetic risk or access to tertiary healthcare facilities, our emerging understanding of the immunobiology of historical/ongoing trauma suggest a more complex explanation for these observed differences.

SUMMARY

Health inequities as observed in pediatric rheumatic diseases are likely to emerge from a complex interplay between social and biological factors. The important contribution of historical and repetitive trauma deserves further exploration.

Association of anti-HSC70 autoantibodies with cutaneous ulceration and severe disease in juvenile dermatomyositis

OBJECTIVES

JDM is an inflammatory myopathy characterized by prominent vasculopathy. AECAs are frequently detected in inflammatory and autoimmune diseases. We sought to determine whether AECAs correlate with clinical features of JDM, and thus serve as biomarkers to guide therapy or predict outcome.

METHODS

Plasma samples from 63 patients with JDM, 49 patients with polyarticular JIA and 40 juvenile healthy controls were used to detect anti-heat shock cognate 71 kDa protein (HSC70) autoantibodies, a newly identified AECA, in ELISA assays. Clinical features were compared between JDM patients with and without anti-HSC70 autoantibodies.

RESULTS

Anti-HSC70 autoantibodies were detected in 35% of patients with JDM, in 0% of patients with JIA (P < 0.0001) and in 0% of healthy donors (P < 0.0001). Both the presence of cutaneous ulcers (59% vs 17%, P < 0.002) and the use of wheelchairs and/or assistive devices (64% vs 27%, P < 0.007) were strongly associated with anti-HSC70 autoantibodies in JDM. High scores on the severity of myositis damage measures at the time of measurement of anti-HSC70 autoantibodies and an increased number of hospitalizations were also associated with anti-HSC70 autoantibodies. Intravenous immunoglobulin therapy was used more often in anti-HSC70 autoantibody-positive patients.

CONCLUSION

Anti-HCS70 autoantibodies are detected frequently in children with JDM and are novel myositis-associated autoantibodies correlating with disease severity.

Chromatin architecture around stroke haplotypes provides evidence that genetic risk is conferred through vascular cells

Introduction: Genome-wide association studies (GWAS) have identified numerous stroke-associated SNPs. To understand how SNPs affect gene expression related to increased stroke risk, we studied epigenetic landscapes surrounding 26 common, validated stroke-associated loci. Methods: We mapped the SNPs to linkage disequilibrium (LD) blocks and examined H3K27ac, H3K4me1, H3K9ac, and H3K4me3 histone marks and transcription-factor binding-sites in pathologically relevant cell types (hematopoietic and vascular cells). Hi-C data were used to identify topologically associated domains (TADs) encompassing the LD blocks and overlapping genes. Results: Fibroblasts, smooth muscle, and endothelial cells showed significant enrichment for enhancer-associated marks within stroke-associated LD blocks. Genes within encompassing TADs reflected vessel homeostasis, cellular turnover, and enzymatic activity. Conclusions: Stroke-associated genetic variants confer risk predominantly through vascular cells rather than hematopoietic cell types.

Epigenetic landscapes of intracranial aneurysm risk haplotypes implicate enhancer function of endothelial cells and fibroblasts in dysregulated gene expression

BACKGROUND

Genome-wide association studies have identified many single nucleotide polymorphisms (SNPs) associated with increased risk for intracranial aneurysm (IA). However, how such variants affect gene expression within IA is poorly understood. We used publicly-available ChIP-Seq data to study chromatin landscapes surrounding risk loci to determine whether IA-associated SNPs affect functional elements that regulate gene expression in cell types comprising IA tissue.

METHODS

We mapped 16 significant IA-associated SNPs to linkage disequilibrium (LD) blocks within human genome. Using ChIP-Seq data, we examined these regions for presence of H3K4me1, H3K27ac, and H3K9ac histone marks (typically associated with latent/active enhancers). This analysis was conducted in several cell types that are present in IA tissue (endothelial cells, smooth muscle cells, fibroblasts, macrophages, monocytes, neutrophils, T cells, B cells, NK cells). In cell types with significant histone enrichment, we used HiC data to investigate topologically associated domains (TADs) encompassing the LD blocks to identify genes that may be affected by IA-associated variants. Bioinformatics were performed to determine the biological significance of these genes. Genes within HiC-defined TADs were also compared to differentially expressed genes from RNA-seq/microarray studies of IA tissues.

RESULTS

We found that endothelial cells and fibroblasts, rather than smooth muscle or immune cells, have significant enrichment for enhancer marks on IA risk haplotypes (p < 0.05). Bioinformatics demonstrated that genes within TADs subsuming these regions are associated with structural extracellular matrix components and enzymatic activity. The majority of histone marked TADs (83% fibroblasts [IMR90], 77% HUVEC) encompassed at least one differentially expressed gene from IA tissue studies.

CONCLUSIONS

These findings provide evidence that genetic variants associated with IA risk act on endothelial cells and fibroblasts. There is strong circumstantial evidence that this may be mediated through altered enhancer function, as genes in TADs encompassing enhancer marks have also been shown to be differentially expressed in IA tissue. These genes are largely related to organization and regulation of the extracellular matrix. This study builds upon our previous (Poppenberg et al., BMC Med Genomics, 2019) by including a more diverse set of data from additional cell types and by identifying potential affected genes (i.e. those in TADs).

Identification of Circulating Gene Expression Signatures of Intracranial Aneurysm in Peripheral Blood Mononuclear Cells

Peripheral blood mononuclear cells (PBMCs) play an important role in the inflammation that accompanies intracranial aneurysm (IA) pathophysiology. We hypothesized that PBMCs have different transcriptional profiles in patients harboring IAs as compared to IA-free controls, which could be the basis for potential blood-based biomarkers for the disease. To test this, we isolated PBMC RNA from whole blood of 52 subjects (24 with IA, 28 without) and performed next-generation RNA sequencing to obtain their transcriptomes. In a randomly assigned discovery cohort of n = 39 patients, we performed differential expression analysis to define an IA-associated signature of 54 genes (q < 0.05 and an absolute fold-change ≥ 1.3). In the withheld validation dataset, these genes could delineate patients with IAs from controls, as the majority of them still had the same direction of expression difference. Bioinformatics analyses by gene ontology enrichment analysis and Ingenuity Pathway Analysis (IPA) demonstrated enrichment of structural regulation processes, intracellular signaling function, regulation of ion transport, and cell adhesion. IPA analysis showed that these processes were likely coordinated through NF-kB, cytokine signaling, growth factors, and TNF activity. Correlation analysis with aneurysm size and risk assessment metrics showed that 4/54 genes were associated with rupture risk. These findings highlight the potential to develop predictive biomarkers from PBMCs to identify patients harboring IAs.

Broadening our understanding of genetic risk for scleroderma/systemic sclerosis by querying the chromatin architecture surrounding the risk haplotypes

BACKGROUND

Genetic variants in the human leukocyte antigen (HLA) locus contribute to the risk for developing scleroderma/systemic sclerosis (SSc). However, there are other replicated loci that also contribute to genetic risk for SSc, and it is unknown whether genetic risk in these non-HLA loci acts primarily on the vasculature, immune system, fibroblasts, or other relevant cell types. We used the Cistrome database to investigate the epigenetic landscapes surrounding 11 replicated SSc associated loci to determine whether SNPs in these loci may affect regulatory elements and whether they are likely to impact a specific cell type.

METHODS

We mapped 11 replicated SNPs to haplotypes and sought to determine whether there was significant enrichment for H3K27ac and H3K4me1 marks, epigenetic signatures of enhancer function, on these haplotypes. We queried pathologically relevant cell types: B cells, endothelial cells, fibroblasts, monocytes, and T cells. We then identified the topologically associated domains (TADs) that encompass the SSc risk haplotypes in primary T cells to identify the full range of genes that may be influenced by SSc causal SNPs. We used gene ontology analyses of the genes within the TADs to gain insight into immunologic functions that might be affected by SSc causal SNPs.

RESULTS

The SSc-associated haplotypes were enriched (p value < 0.01) for H3K4me1/H3K27ac marks in monocytes. Enrichment of one of the two histone marks was found in B cells, fibroblasts, and T cells. No enrichment was identified in endothelial cells. Ontological analyses of genes within the TADs encompassing the risk haplotypes showed enrichment for regulation of transcription, protein binding, activation of T lymphocytes, and proliferation of immune cells.

CONCLUSIONS

The 11 non-HLA SSc risk haplotypes queried are highly enriched for H3K4me1/H3K27ac-marked regulatory elements in a broad range of immune cells and fibroblasts. Furthermore, in immune cells, the risk haplotypes belong to larger chromatin structures encompassing genes that regulate a wide array of immune processes associated with SSc pathogenesis. Though importance of the vasculature in the pathobiology of SSc is widely accepted, we were unable to find evidence for genetic influences on endothelial cell function in these regions.

Whole blood transcriptome biomarkers of unruptured intracranial aneurysm

BACKGROUND

The rupture of an intracranial aneurysm (IA) causes devastating subarachnoid hemorrhages, yet most IAs remain undiscovered until they rupture. Recently, we found an IA RNA expression signature of circulating neutrophils, and used transcriptome data to build predictive models for unruptured IAs. In this study, we evaluate the feasibility of using whole blood transcriptomes to predict the presence of unruptured IAs.

METHODS

We subjected RNA from peripheral whole blood of 67 patients (34 with unruptured IA, 33 without IA) to next-generation RNA sequencing. Model genes were identified using the least absolute shrinkage and selection operator (LASSO) in a random training cohort (n = 47). These genes were used to train a Gaussian Support Vector Machine (gSVM) model to distinguish patients with IA. The model was applied to an independent testing cohort (n = 20) to evaluate performance by receiver operating characteristic (ROC) curve. Gene ontology and pathway analyses investigated the underlying biology of the model genes.

RESULTS

We identified 18 genes that could distinguish IA patients in a training cohort with 85% accuracy. This SVM model also had 85% accuracy in the testing cohort, with an area under the ROC curve of 0.91. Bioinformatics reflected activation and recruitment of leukocytes, activation of macrophages, and inflammatory response, suggesting that the biomarker captures important processes in IA pathogenesis.

CONCLUSIONS

Circulating whole blood transcriptomes can detect the presence of unruptured IAs. Pending additional testing in larger cohorts, this could serve as a foundation to develop a simple blood-based test to facilitate screening and early detection of IAs.

Classification models using circulating neutrophil transcripts can detect unruptured intracranial aneurysm

Background

Intracranial aneurysms (IAs) are dangerous because of their potential to rupture. We previously found significant RNA expression differences in circulating neutrophils between patients with and without unruptured IAs and trained machine learning models to predict presence of IA using 40 neutrophil transcriptomes. Here, we aim to develop a predictive model for unruptured IA using neutrophil transcriptomes from a larger population and more robust machine learning methods.

Methods

Neutrophil RNA extracted from the blood of 134 patients (55 with IA, 79 IA-free controls) was subjected to next-generation RNA sequencing. In a randomly-selected training cohort (n = 94), the Least Absolute Shrinkage and Selection Operator (LASSO) selected transcripts, from which we constructed prediction models via 4 well-established supervised machine-learning algorithms (K-Nearest Neighbors, Random Forest, and Support Vector Machines with Gaussian and cubic kernels). We tested the models in the remaining samples (n = 40) and assessed model performance by receiver-operating-characteristic (ROC) curves. Real-time quantitative polymerase chain reaction (RT-qPCR) of 9 IA-associated genes was used to verify gene expression in a subset of 49 neutrophil RNA samples. We also examined the potential influence of demographics and comorbidities on model prediction.

Results

Feature selection using LASSO in the training cohort identified 37 IA-associated transcripts. Models trained using these transcripts had a maximum accuracy of 90% in the testing cohort. The testing performance across all methods had an average area under ROC curve (AUC) = 0.97, an improvement over our previous models. The Random Forest model performed best across both training and testing cohorts. RT-qPCR confirmed expression differences in 7 of 9 genes tested. Gene ontology and IPA network analyses performed on the 37 model genes reflected dysregulated inflammation, cell signaling, and apoptosis processes. In our data, demographics and comorbidities did not affect model performance.

Conclusions

We improved upon our previous IA prediction models based on circulating neutrophil transcriptomes by increasing sample size and by implementing LASSO and more robust machine learning methods. Future studies are needed to validate these models in larger cohorts and further investigate effect of covariates.

Utility of CD4 count measurement in the era of universal antiretroviral therapy: an analysis of routine laboratory data in Botswana

OBJECTIVES

National guidelines in Botswana recommend baseline CD4 count measurement and both CD4 and HIV viral load (VL) monitoring post-antiretroviral therapy (ART) initiation. We evaluated the utility of CD4 count measurement in Botswana in the era of universal ART.

METHODS

CD4 and VL data were analysed for HIV-infected adults undergoing CD4 count measurement in 2015-2017 at the Botswana Harvard HIV-Reference Laboratory. We determined (1) the proportion of individuals with advanced HIV disease (CD4 count < 200 cells/µL) at initial CD4 assessment, (2) the proportion with an initial CD4 count ≥ 200 cells/µL experiencing a subsequent decline in CD4 count to < 200 cells/µL, and (3) the proportion of these immunologically failing individuals who had virological failure. Logistic regression modelling examined factors associated with advanced HIV disease. CD4 count trajectories were assessed using locally weighted scatterplot smoothing (LOWESS) regression.

RESULTS

Twenty-five per cent (3571/14 423) of individuals with an initial CD4 assessment during the study period had advanced HIV disease at baseline. Older age [≥ 35 years; adjusted odds ratio (aOR) 1.9; 95% confidence interval (CI) 1.8-2.1] and male sex were associated with advanced HIV disease. Fifty per cent (7163/14 423) of individuals had at least two CD4 counts during the study period. Of those with an initial CD4 count ≥ 200 cells/µL, 4% (180/5061) experienced a decline in CD4 count to < 200 cells/µL; the majority of CD4 count declines were in virologically suppressed individuals and transient.

CONCLUSIONS

One-quarter of HIV-positive individuals in Botswana still present with advanced HIV disease, highlighting the importance of baseline CD4 count measurement to identify this at-risk population. Few with a baseline CD4 count ≥ 200 cells/µL experienced a drop below 200 cells/µL, suggesting limited utility for ongoing CD4 monitoring.

Broadening our understanding of the genetics of Juvenile Idiopathic Arthritis (JIA): Interrogation of three dimensional chromatin structures and genetic regulatory elements within JIA-associated risk loci

Objective

The risk loci for juvenile idiopathic arthritis (JIA) consist of extended haplotypes that include functional elements in addition to canonical coding genes. As with most autoimmune diseases, the risk haplotypes for JIA are highly enriched for H3K4me1/H3K27ac histone marks, epigenetic signatures that typically identify poised or active enhancers. In this study, we test the hypothesis that genetic risk for JIA is exerted through altered enhancer-mediated gene regulation.

Methods

We mined publically available HiC and other chromatin conformation data to determine whether H3K27ac-marked regions in 25 JIA risk loci showed physical evidence of contact with gene promoters. We also used in vitro reporter assays to establish as proof-of-concept the idea that genetic variants in linkage disequilibrium with GWAS-identified tag SNPs alter enhancer function.

Results

All 25 loci examined showed multiple contact sites in the 4 different cell lines that we queried. These regions were characterized by HiC-defined loop structures that included 237 immune-related genes. Using in vitro assays, we found that a 657 bp, H3K4me1/H3K27-marked region within the first intron of IL2RA shows enhancer activity in reporter assays, and this activity is attenuated by SNPs on the IL2RA haplotype that we identified using whole genome sequencing of children with JIA. Similarly, we identified a 1,669 bp sequence in an intergenic region of the IL6R locus where SNPs identified in children with JIA increase enhancer function in reporter assays.

Conclusions

These studies provide evidence that altered enhancer function contributes to genetic risk in JIA. Further studies to identify the specific target genes of genetically altered enhancers are warranted.

Double negative T cells, a potential biomarker for systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is an autoimmune disease that is a challenge to diagnose and treat. There is an urgent need for biomarkers to help define organ involvement, and more effective therapies. A unique population of T cells, the CD3⁺CD4⁻CD8⁻ (DNeg) cells, is significantly increased in lupus patients. Twenty-seven cases (53%) of pediatric SLE patients had elevated DNeg cells in their peripheral blood, which correlated with kidney function (R² = 0.54). Significant infiltration of DNeg cells was observed in both adult and pediatric lupus kidneys by immunofluorescence. For the first time, this study provides direct evidence that DNeg cells facilitate kidney injury in preclinical 8-week-old MRL/lpr lupus mice. In lupus mice, the increase in DNeg cells tracked with worsening disease and correlated with kidney function (R² = 0.85). Our results show that DNeg cells per se can cause kidney dysfunction, increase in number with increase in disease pathology, and could serve as a potential biomarker.

The feasibility of developing biomarkers from peripheral blood mononuclear cell RNAseq data in children with juvenile idiopathic arthritis using machine learning approaches

Background

The response to treatment for juvenile idiopathic arthritis (JIA) can be staged using clinical features. However, objective laboratory biomarkers of remission are still lacking. In this study, we used machine learning to predict JIA activity from transcriptomes from peripheral blood mononuclear cells (PBMCs). We included samples from children with Native American ancestry to determine whether the model maintained validity in an ethnically heterogeneous population.

Methods

Our dataset consisted of 50 samples, 23 from children in remission and 27 from children with an active disease on therapy. Nine of these samples were from children with mixed European/Native American ancestry. We used 4 different machine learning methods to create predictive models in 2 populations: the whole dataset and then the samples from children with exclusively European ancestry.

Results

In both populations, models were able to predict JIA status well, with training accuracies > 74% and testing accuracies > 78%. Performance was better in the whole dataset model. We note a high degree of overlap between genes identified in both populations. Using ingenuity pathway analysis, genes from the whole dataset associated with cell-to-cell signaling and interactions, cell morphology, organismal injury and abnormalities, and protein synthesis.

Conclusions

This study demonstrates it is feasible to use machine learning in conjunction with RNA sequencing of PBMCs to predict JIA stage. Thus, developing objective biomarkers from easy to obtain clinical samples remains an achievable goal.

Epigenetic landscapes suggest that genetic risk for intracranial aneurysm operates on the endothelium

BACKGROUND

Genetics play an important role in intracranial aneurysm (IA) pathophysiology. Genome-wide association studies have identified several single nucleotide polymorphisms (SNPs) that are linked to IA but how they affect disease pathobiology remains poorly understood. We used Encyclopedia of DNA Elements (ENCODE) data to investigate the epigenetic landscapes surrounding genetic risk loci to determine if IA-associated SNPs affect functional elements that regulate gene expression and if those SNPs are most likely to impact a specific type of cells.

METHODS

We mapped 16 highly significant IA-associated SNPs to linkage disequilibrium (LD) blocks within the human genome. Within these regions, we examined the presence of H3K4me1 and H3K27ac histone marks and CCCTC-binding factor (CTCF) and transcription-factor binding sites using chromatin immunoprecipitation-sequencing (ChIP-Seq) data. This analysis was conducted in several cell types relevant to endothelial (human umbilical vein endothelial cells [HUVECs]) and inflammatory (monocytes, neutrophils, and peripheral blood mononuclear cells [PBMCs]) biology. Gene ontology analysis was performed on genes within extended IA-risk regions to understand which biological processes could be affected by IA-risk SNPs. We also evaluated recently published data that showed differential methylation and differential ribonucleic acid (RNA) expression in IA to investigate the correlation between differentially regulated elements and the IA-risk LD blocks.

RESULTS

The IA-associated LD blocks were statistically significantly enriched for H3K4me1 and/or H3K27ac marks (markers of enhancer function) in endothelial cells but not in immune cells. The IA-associated LD blocks also contained more binding sites for CTCF in endothelial cells than monocytes, although not statistically significant. Differentially methylated regions of DNA identified in IA tissue were also present in several IA-risk LD blocks, suggesting SNPs could affect this epigenetic machinery. Gene ontology analysis supports that genes affected by IA-risk SNPs are associated with extracellular matrix reorganization and endopeptidase activity.

CONCLUSION

These findings suggest that known genetic alterations linked to IA risk act on endothelial cell function. These alterations do not correlate with IA-associated gene expression signatures of circulating blood cells, which suggests that such signatures are a secondary response reflecting the presence of IA rather than indicating risk for IA.

Plasma exosomes from children with juvenile dermatomyositis are taken up by human aortic endothelial cells and are associated with altered gene expression in those cells

BACKGROUND

The pathology of juvenile dermatomyositis (JDM) is characterized by prominent vessel wall and perivascular inflammation. This feature of the disease has remained unexplained and under-investigated. We have hypothesized that plasma exosomes, which play an important role in inter-cellular communication, may play a role in the vascular injury associated with JDM.

OBJECTIVE

To characterize the circulating exosomes of children with JDM and determine whether the small RNA cargoes within those exosomes are capable of altering transcriptional programs within endothelial cells.

DESIGN/METHODS

We purified exosomes from plasma samples of children with active, untreated JDM (n = 6) and healthy controls (n = 9). We characterized the small RNA cargoes in JDM and control exosomes by RNA sequencing using the Illumina HiSeq 2500 platform. We then incubated isolated exosomes from healthy controls and children with JDM with cultured human aortic endothelial cells (HAEC) for 24 h. Fluorescence microscopy was used to confirm that both control and JDM exosomes were taken up by HAEC. RNA was then purified from HAEC that had been incubated with either control or JDM exosomes and sequenced on the Illumina platform. Differential expression of mRNAs from HAEC incubated with control or JDM exosomes was ascertained using standard computational methods. Finally, we assessed the degree to which differential gene expression in HAEC could be attributed to the different small RNA cargoes in JDM vs control exosomes using conventional and novel analytic methods.

RESULTS

We identified 10 small RNA molecules that showed differential abundance when we compared JDM and healthy control exosomes. Fluorescence microscopy of labeled exosomes confirmed that both JDM and control exosomes were taken up by HAEC. Differential gene expression analysis revealed 59 genes that showed differential expression between HAEC incubated with JDM exosomes vs HAEC incubated with exosomes from controls. Statistical analysis of gene expression data demonstrated that multiple miRNAs exerted transcriptional control on multiple genes with HAEC.

CONCLUSIONS

Plasma exosomes from children with active, untreated JDM are taken up by HAEC and are associated with alterations in gene expression in those cells. These findings provide new insight into potential mechanisms leading to the targeting of vascular tissue by the immune system in JDM.

Biomarkers from circulating neutrophil transcriptomes have potential to detect unruptured intracranial aneurysms

BACKGROUND

Intracranial aneurysms (IAs) are dangerous because of their potential to rupture and cause deadly subarachnoid hemorrhages. Previously, we found significant RNA expression differences in circulating neutrophils between patients with unruptured IAs and aneurysm-free controls. Searching for circulating biomarkers for unruptured IAs, we tested the feasibility of developing classification algorithms that use neutrophil RNA expression levels from blood samples to predict the presence of an IA.

METHODS

Neutrophil RNA extracted from blood samples from 40 patients (20 with angiography-confirmed unruptured IA, 20 angiography-confirmed IA-free controls) was subjected to next-generation RNA sequencing to obtain neutrophil transcriptomes. In a randomly-selected training cohort of 30 of the 40 samples (15 with IA, 15 controls), we performed differential expression analysis. Significantly differentially expressed transcripts (false discovery rate < 0.05, fold change ≥ 1.5) were used to construct prediction models for IA using four well-known supervised machine-learning approaches (diagonal linear discriminant analysis, cosine nearest neighbors, nearest shrunken centroids, and support vector machines). These models were tested in a testing cohort of the remaining 10 neutrophil samples from the 40 patients (5 with IA, 5 controls), and model performance was assessed by receiver-operating-characteristic (ROC) curves. Real-time quantitative polymerase chain reaction (PCR) was used to corroborate expression differences of a subset of model transcripts in neutrophil samples from a new, separate validation cohort of 10 patients (5 with IA, 5 controls).

RESULTS

The training cohort yielded 26 highly significantly differentially expressed neutrophil transcripts. Models using these transcripts identified IA patients in the testing cohort with accuracy ranging from 0.60 to 0.90. The best performing model was the diagonal linear discriminant analysis classifier (area under the ROC curve = 0.80 and accuracy = 0.90). Six of seven differentially expressed genes we tested were confirmed by quantitative PCR using isolated neutrophils from the separate validation cohort.

CONCLUSIONS

Our findings demonstrate the potential of machine-learning methods to classify IA cases and create predictive models for unruptured IAs using circulating neutrophil transcriptome data. Future studies are needed to replicate these findings in larger cohorts.

Using Chromatin Architecture to Understand the Genetics and Transcriptomics of Juvenile Idiopathic Arthritis

The presence of abnormal gene expression signatures is a well-described feature of the oligoarticular and polyarticular forms of juvenile idiopathic arthritis. In this review, we discuss how new insights into genetic risk for JIA and the three dimensional architecture of the genome may be used to develop a better understanding of the mechanisms driving these gene expression patterns.

Multiple target autoantigens on endothelial cells identified in juvenile dermatomyositis using proteomics

Objective

Although generally classified within the group of inflammatory myopathies, JDM displays many pathological features of vasculitis. Previous work has shown that AECA are abundant in other forms of vasculitis. We therefore investigated whether such antibodies might also be detected in JDM.

Methods

We screened plasma from children with JDM for the presence of AECA by western blotting and 2D gel electrophoresis (2DE) using proteins extracted from human aortic endothelial cells as the substrate. We performed mass spectrometry to identify candidate antigens from 2DE gels and used ELISA to confirm the presence of specific antibodies.

Results

We identified 22 candidate target autoantigens for AECA probed with JDM plasma. Interestingly, 17 of these 22 target antigens were proteins associated with antigen processing and protein trafficking. ELISA confirmed the presence of antibodies to heat shock cognate 71 kDa protein in JDM plasma, particularly in children with active, untreated disease.

Conclusion

Children with JDM express antibodies to autoantigens in endothelial cells. The clinical and pathological significance of such autoantibodies require further investigation.

Soluble inflammatory mediators induce transcriptional re-organization that is independent of dna methylation changes in cultured human chorionic villous trophoblasts

The studies proposed here were undertaken to test the hypothesis that, under specific circumstances (e.g., a strong enough inflammatory stimulus), genes that are repressed at the maternal-fetal interface via DNA methylation might be de-methylated, allowing either a maternal immune response to the semi-allogenic fetus or the onset of early labor. Chorionic trophoblasts (CT) were isolated from fetal membranes, followed by incubation with medium from LPS-activated PBMC or resting PBMC medium for 2 h. RNA and DNA were isolated from the cells for RNA-seq and DNA methylation studies. Two hrs after being exposed to conditioned medium from LPS-activated PBMC, CT showed differential expression of 114 genes, all but 2 of which showed higher expression in the stimulated cells than is the unstimulated cells. We also identified 318 differentially methylated regions (DMRs) that associated with 306 genes (155 protein coding genes) in the two groups, but the observed methylation changes had negligible impact on the observed transcriptional changes in CT.

CT display complex patterns of transcription in response to inflammation. DNA methylation does not appear to be an important regulator of the observed transcriptional changes.

Unique clinical characteristics, autoantibodies and medication use in Native American patients with systemic lupus erythematosus

Objective

Systemic lupus erythematosus (SLE) is a systemic autoimmune disease with varied morbidity and mortality. We assessed clinical presentations, autoantibody specificities and therapeutic interventions in Native American (NA) patients with SLE.

Methods

Patients with SLE meeting 1997 American College of Rheumatology classification criteria (n=3148) were enrolled between 1992 and 2010 in the multiethnic, cross-sectional Lupus Family Registry and Repository. Clinical, demographic and therapeutic information were extracted from medical records using a standardised form and formalised training. Autoantibodies were assessed by indirect immunofluorescence (antinuclear antibodies (ANA) and antidouble-stranded DNA), precipitin (ENA) and ELISA (IgG and IgM anticardiolipins).

Results

NA patients met SLE classification at a younger age (29.89±12.3 years) than European Americans (EA; 32.02±12.87, P=0.0157) and a similar age to African-Americans (AAs) and Hispanics (HIS). More NA patients had concurrent rheumatic diseases or symptoms, such as Raynaud’s phenomenon, interstitial lung disease, Sjӧgren’s syndrome and systemic sclerosis. Compared with EAs, NAs were more likely to have high-titre ANA (≥1:3240; P<0.0001) and had more SLE-associated autoantibodies. Autoantibodies with unknown specificities were more common in NAs (41%) compared with other racial/ethnic groups in this collection (AA: 24%, P=0.0006; EA: 17%, P<0.0001; HIS: 23%, P=0.0050). Fewer NA patients used hydroxychloroquine (68%) compared with others (AA: 74%, P=0.0308; EA: 79%, P=0.0001, HIS: 77%, P=0.0173); this was influenced by lower hydroxychloroquine use in NA patients from Latin America (32%). NA patients had higher rates of methotrexate use (28%) compared with AA (18%, P=0.0006) and HIS patients (14%, P=0.0003), higher azathioprine use (38%) compared with EA patients (30%, P=0.0105) and higher mycophenolate mofetil use (26%) compared with EA (17%, P=0.0012) and HIS patients (11%, P<0.0001).

Conclusions

NA patients are diagnosed with SLE earlier in life and present worse concurrent rheumatic disease symptoms than EA patients. NA patients also are more likely to have expanded autoantibody profiles and precipitins of unknown specificities.

Circulating neutrophil transcriptome may reveal intracranial aneurysm signature

Background

Unruptured intracranial aneurysms (IAs) are typically asymptomatic and undetected except for incidental discovery on imaging. Blood-based diagnostic biomarkers could lead to improvements in IA management. This exploratory study examined circulating neutrophils to determine whether they carry RNA expression signatures of IAs.

Methods

Blood samples were collected from patients receiving cerebral angiography. Eleven samples were collected from patients with IAs and 11 from patients without IAs as controls. Samples from the two groups were paired based on demographics and comorbidities. RNA was extracted from isolated neutrophils and subjected to next-generation RNA sequencing to obtain differential expressions for identification of an IA-associated signature. Bioinformatics analyses, including gene set enrichment analysis and Ingenuity Pathway Analysis, were used to investigate the biological function of all differentially expressed transcripts.

Results

Transcriptome profiling identified 258 differentially expressed transcripts in patients with and without IAs. Expression differences were consistent with peripheral neutrophil activation. An IA-associated RNA expression signature was identified in 82 transcripts (p<0.05, fold-change ≥2). This signature was able to separate patients with and without IAs on hierarchical clustering. Furthermore, in an independent, unpaired, replication cohort of patients with IAs (n = 5) and controls (n = 5), the 82 transcripts separated 9 of 10 patients into their respective groups.

Conclusion

Preliminary findings show that RNA expression from circulating neutrophils carries an IA-associated signature. These findings highlight a potential to use predictive biomarkers from peripheral blood samples to identify patients with IAs.

Chromatin landscapes and genetic risk for juvenile idiopathic arthritis

Background

The transcriptomes of peripheral blood cells in children with juvenile idiopathic arthritis (JIA) have distinct transcriptional aberrations that suggest impairment of transcriptional regulation. To gain a better understanding of this phenomenon, we studied known JIA genetic risk loci, the majority of which are located in non-coding regions, where transcription is regulated and coordinated on a genome-wide basis. We examined human neutrophils and CD4 primary T cells to identify genes and functional elements located within those risk loci.

Methods

We analyzed RNA sequencing (RNA-Seq) data, H3K27ac and H3K4me1 chromatin immunoprecipitation-sequencing (ChIP-Seq) data, and previously published chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) data to characterize the chromatin landscapes within the known JIA-associated risk loci.

Results

In both neutrophils and primary CD4+ T cells, the majority of the JIA-associated linkage disequilibrium (LD) blocks contained H3K27ac and/or H3K4me1 marks. These LD blocks were also binding sites for a small group of transcription factors, particularly in neutrophils. Furthermore, these regions showed abundant intronic and intergenic transcription in neutrophils. In neutrophils, none of the genes that were differentially expressed between untreated patients with JIA and healthy children were located within the JIA-risk LD blocks. In CD4+ T cells, multiple genes, including HLA-DQA1, HLA-DQB2, TRAF1, and IRF1 were associated with the long-distance interacting regions within the LD regions as determined from ChIA-PET data.

Conclusions

These findings suggest that genetic risk contributes to the aberrant transcriptional control observed in JIA. Furthermore, these findings demonstrate the challenges of identifying the actual causal variants within complex genomic/chromatin landscapes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-017-1260-x) contains supplementary material, which is available to authorized users.

Chromatin landscapes and genetic risk in systemic lupus

BACKGROUND

Systemic lupus erythematosus (SLE) is a multi-system, complex disease in which the environment interacts with inherited genes to produce broad phenotypes with inter-individual variability. Of 46 single nucleotide polymorphisms (SNPs) shown to confer genetic risk for SLE in recent genome-wide association studies, 30 lie within noncoding regions of the human genome. We therefore sought to identify and describe the functional elements (aside from genes) located within these regions of interest.

METHODS

We used chromatin immunoprecipitation followed by sequencing to identify epigenetic marks associated with enhancer function in adult neutrophils to determine whether enhancer-associated histone marks were enriched within the linkage disequilibrium (LD) blocks encompassing the 46 SNPs of interest. We also interrogated available data in Roadmap Epigenomics for CD4+ T cells and CD19+ B cells to identify these same elements in lymphoid cells.

RESULTS

All three cell types demonstrated enrichment of enhancer-associated histone marks compared with genomic background within LD blocks encoded by SLE-associated SNPs. In addition, within the promoter regions of these LD blocks, all three cell types demonstrated enrichment for transcription factor binding sites above genomic background. In CD19+ B cells, all but one of the LD blocks of interest were also enriched for enhancer-associated histone marks.

CONCLUSIONS

Much of the genetic risk for SLE lies within or near genomic regions of disease-relevant cells that are enriched for epigenetic marks associated with enhancer function. Elucidating the specific roles of these noncoding elements within these cell-type-specific genomes will be crucial to our understanding of SLE pathogenesis.

Using proteomic and genomic methods to understand JDM

Juvenile dermatomyositis is a complex illness characterized by vascular/perivascular inflammation, primarily in the skin and muscles. In this review, we discuss how proteomic and genomic technologies have expanded our understanding of the immune pathogenesis of this disease. We will also discuss further directions that the field may take to use existing and developing technologies to further our understanding of this often-perplexing disease.

Proceedings of the 2016 Childhood Arthritis and Rheumatology Research Alliance (CARRA) Scientific Meeting : Toronto, Canada. 14-17 April 2016

P1 Serologic evidence of gut-driven systemic inflammation in juvenile idiopathic arthritis Lampros Fotis, Nur Shaikh, Kevin Baszis, Anthony French, Phillip Tarr P2 Oral health and anti-citrullinated peptide antibodies (ACPA) in juvenile idiopathic arthritis Sriharsha Grevich, Peggy Lee, Sarah Ringold, Brian Leroux, Hannah Leahey, Megan Yuasa, Jessica Foster, Jeremy Sokolove, Lauren Lahey, William Robinson, Joshua Newsom, Anne Stevens P3 Novel autoantigens for endothelial cell antibodies in pediatric rheumatic diseases identified by proteomics Rie Karasawa, Mayumi Tamaki, Megumi Tanaka, Toshiko Sato, Kazuo Yudoh, James N. Jarvis P4 Transcriptional profiling reveals monocyte signature associated with JIA patient poor response to methotrexate Halima Moncrieffe, Mark F. Bennett, Monica Tsoras, Lorie Luyrink, Huan Xu, Sampath Prahalad, Paula Morris, Jason Dare, Peter A. Nigrovic, Margalit Rosenkranz, Mara Becker, Kathleen M. O’Neil, Thomas Griffin, Daniel J. Lovell, Alexei A. Grom, Mario Medvedovic, Susan D. Thompson P5 A multi-dimensional genomic map for polyarticular juvenile idiopathic arthritis Lisha Zhu, Kaiyu Jiang, Laiping Wong, Michael J Buck, Yanmin Chen, Halima Moncrieffe, Laura Brungs, Tao Liu, Ting Wang, James N Jarvis P6 Tocilizumab for treatment of children with refractory JIA Khaled Alsaeid, Jasim Alfailakawi, Hamid Alenezi, Hazim Alsaeed P7 Clinical characteristics of the initial patients enrolled in the Childhood Arthritis and Rheumatology Research Alliance (CARRA) Registry Tim Beukelman, Marc Natter, Norm Ilowite, Kelly Mieszkalski, Grendel Burrell, Brian Best, Helen Bristow, Shannon Carr, Anne Dennos, Rachel Kaufmann, Yukiko Kimura, Laura Schanberg P8 Comparative performance of small and large clinical centers in a comprehensive pediatric rheumatology disease registry Peter R Blier P9 Clinical characteristics of children with membranous lupus nephritis: The Childhood Arthritis and Rheumatology Research Alliance Legacy Registry Alexis Boneparth, Scott E. Wenderfer, L. Nandini Moorthy, Suhas M. Radhakrishna, Anna Carmela P. Sagcal-Gironella, Emily von Scheven P10 Rituximab use in pediatric lupus anticoagulant hypoprothrombinemia syndrome - a two center experience Kader Cetin Gedik, Salma Siddique, Cassyanne L. Aguiar, Doruk Erkan P11 Predictors of complementary and alternative medicine use and response in children with musculoskeletal conditions Ezra Cohen, Yvonne Lee, Michelle Dossett, Darshan Mehta, Roger Davis P12 Comparison of pediatric rheumatology and nephrology survey results for the treatment of refractory proliferative lupus nephritis and renal flare in juvenile SLE Mileka Gilbert, Beatrice Goilav, Esra Meidan, Joyce Hsu, Alexis Boneparth, Anabelle Chua, Stacy Ardoin, Scott E. Wenderfer, Emily Von Scheven, Natasha M. Ruth P13 Transitioning lupus patients from pediatric to adult rheumatology Joyce Hui-Yuen, Kader Cetin Gedik, Liza Bermudez, Ashlea Cook, Lisa Imundo, Amy Starr, Andrew Eichenfield, Anca Askanase P14 The systemic juvenile idiopathic arthritis cohort of the Childhood Arthritis & Rheumatology Research Alliance Registry Ginger Janow, Laura E. Schanberg, Soko Setoguchi, Victor Hasselblad, Elizabeth D. Mellins, Rayfel Schneider, Yukiko Kimura, The CARRA Legacy Registry Investigators P15 Results of the pilot study of the Childhood Arthritis and Rheumatology Research Alliance (CARRA) consensus treatment plans for new-onset systemic juvenile idiopathic arthritis Yukiko Kimura, Sriharsha Grevich, Timothy Beukelman, Esi Morgan, T Brent Graham, Maria Ibarra, Yonit Sterba Ruas, Marisa Klein-Gitelman, Karen Onel, Sampath Prahalad, Marilynn Punaro, Sarah Ringold, Dana Toib, Heather Van Mater, Jennifer E. Weiss, Pamela F. Weiss, Kelly Mieszkalski, Laura E. Schanberg P16 A systemic review of pain relief modalities in juvenile idiopathic arthritis: First step in developing a novel decision support intervention Timothy S. H. Kwok, Jacinthe Bisaillon, Christine Smith, Lucie Brosseau, Jennifer Stinson, Adam M. Huber, Ciaran M. Duffy, Karine Toupin April P17 Barriers and facilitators to care retention for pediatric systemic lupus erythematous patients in South Africa: A qualitative study Laura B Lewandowski, Christiaan Scott P18 Evaluating the feasibility of conducting comparative effectiveness studies in juvenile Localized Scleroderma (jLS) Suzanne C. Li, Kathryn S. Torok, C. Egla Rabinovich, Sandy D. Hong, Mara L Becker, Fatma Dedeoglu, Maria F. Ibarra, Polly J Ferguson, Rob C. Fuhbrigge, Katie G. Stewart, Elena Pope, Ronald M. Laxer, Thomas G. Mason, Gloria C. Higgins, Xiaohu Li, Marilynn G. Punaro, George Tomlinson, Eleanor Pullenayegum, John Matelski, Laura Schanberg, Brian M. Feldman P19 Tonsillar histology in patients with periodic fever, aphthous stomatitis, pharyngitis, adenitis (PFAPA) syndrome Kalpana Manthiram, Hernan Correa, Kathryn Edwards P20 Clinical course of juvenile dermatomyositis presenting as skin predominant disease Edward J. Oberle, Michelle Bayer, Dominic O. Co, Hatice Ezgi Baris, Yvonne Chiu, Adam Huber, Susan Kim P21 A Survey of musculoskeletal ultrasound practices of pediatric rheumatologists in North America Edward J Oberle, Timothy Beukelman P22 Assessment, classification and treatment of calcinosis as a complication of juvenile dermatomyositis: A survey of pediatric rheumatologists by the Childhood Arthritis and Rheumatology Research Alliance Amir B. Orandi, Kevin W. Baszis, Vikas Dharnidharka, Mark F. Hoeltzel, for the CARRA JDM Committee P23 CARRA dermatomyositis CTP pilot study Ann Reed, Adam Huber, George Tomlinson, Eleanor Pullenayegum, John Matelski, Y. Ingrid Goh, Laura Schanberg, Brian M. Feldman P24 Unexpectedly high incidences and prolonged disease activity in children with chronic non-bacterial osteomyelitis (CNO) as compared to bacterial osteomyelitis Anja Schnabel, Ursula Range, Gabriele Hahn, Timo Siepmann, Reinhard Berner, Christian Michael Hedrich P25 Juvenile systemic sclerosis cohort within the Childhood Arthritis and Rheumatology Research Alliance (CARRA) Legacy Registry: Follow up characteristics Brandi Stevens, Kathryn S. Torok, Suzanne Li, Nicole Hershey, Megan Curran, Gloria Higgins, Katharine Moore, Egla Rabinovich, Anne M. Stevens, for the CARRA Registry Investigators P26 Development and usability testing of an iPad and desktop psycho-educational game for children with Juvenile Idiopathic Arthritis and their parents Jennifer Stinson, Mark Connelly, Adam Huber, Nadia Luca, Lynn Spiegel, Argerie Tsimicalis, Stephanie Luca, Naweed Tajuddin, Roberta Berard, Julia Barsalou, Sarah Campillo, Paul Dancey, Ciaran Duffy, Brian Feldman, Nicole Johnson, Patrick McGrath, Natalie Shiff, Shirley Tse, Lori Tucker, Charles Victor P27 iCanCope^TM: User-centred design and development of a smartphone app to support self-management for youth with arthritis pain Jennifer Stinson, Chitra Lalloo, Lauren Harris, Joseph Cafazzo, Lynn Spiegel, Brian Feldman, Nadia Luca, Ronald Laxer P28 Accessing pediatric rheumatology care: Despite barriers, few parents prefer telemedicine Danielle R. Bullock, Richard K. Vehe, Lei Zhang, Colleen K. Correll¹ P29 Exploration of factors contributing to time to achieve clinically inactive disease (CID) in juvenile idiopathic arthritis (JIA): A preliminary report Suhas Ganguli, Max Shenberger, Ritesh Korumilli, Beth Gottlieb P30 Pediatric rheumatology referral patterns: Presenting complaints of new patients at a large, urban academic center Martha Rodriguez, Deirdre de Ranieri, Karen Onel, Linda Wagner-Weiner, Melissa Tesher P31 Quality improvement (QI) initiatives in childhood systemic lupus erythematosus (cSLE) Elizabeth Roth Wojcicki, Kristyn L. Maletta, Dominic O. Co, Marsha Malloy, Sarah Thomson, Judyann C. Olson P32 Proliferative lupus nephritis in juvenile SLE: Support from the pediatric nephrology community for the definitions of responsiveness and flare in the 2012 consensus treatment plans Scott E. Wenderfer, Mileka Gilbert, Joyce Hsu, Sangeeta Sule, Tamar B. Rubinstein, Beatrice Goilav, Daryl M. Okamura, Annabelle Chua, Laurence A. Greenbaum, Jerome C. Lane, Emily von Scheven, Stacy P. Ardoin, Natasha M. Ruth P33 The steroid taper app: Making of a mobile app Jennifer M. P. Woo, Marsha M. Malloy, James A. Jegers, Dustin J. Hahn, Mary K. Hintermeyer, Stacey M. Martinetti, Gretchen R. Heckel, Elizabeth L. Roth-Wojcicki, Dominic O. Co

C5a induces caspase-dependent apoptosis in brain vascular endothelial cells in experimental lupus

Blood-brain barrier (BBB) dysfunction complicates central nervous system lupus, an important aspect of systemic lupus erythematosus. To gain insight into the underlying mechanism, vascular corrosion casts of brain were generated from the lupus mouse model, MRL/lpr mice and the MRL/MpJ congenic controls. Scanning electron microscopy of the casts showed loss of vascular endothelial cells in lupus mice compared with controls. Immunostaining revealed a significant increase in caspase 3 expression in the brain vascular endothelial cells, which suggests that apoptosis could be an important mechanism causing cell loss, and thereby loss of BBB integrity. Complement activation occurs in lupus resulting in increased generation of circulating C5a, which caused the endothelial layer to become 'leaky'. In this study, we show that C5a and lupus serum induced apoptosis in cultured human brain microvascular endothelial cells (HBMVECs), whereas selective C5a receptor 1 (C5aR1) antagonist reduced apoptosis in these cells, demonstrating C5a/C5aR1-dependence. Gene expression of initiator caspases, caspase 1 and caspase 8, and pro-apoptotic proteins death-associated protein kinase 1, Fas-associated protein (FADD), cell death-inducing DNA fragmentation factor 45 000 MW subunit A-like effector B (CIDEB) and BCL2-associated X protein were increased in HBMVECs treated with lupus serum or C5a, indicating that both the intrinsic and extrinsic apoptotic pathways could be critical mediators of brain endothelial cell apoptosis in this setting. Overall, our findings suggest that C5a/C5aR1 signalling induces apoptosis through activation of FADD, caspase 8/3 and CIDEB in brain endothelial cells in lupus. Further elucidation of the underlying apoptotic mechanisms mediating the reduced endothelial cell number is important in establishing the potential therapeutic effectiveness of C5aR1 inhibition that could prevent and/or reduce BBB alterations and preserve the physiological function of BBB in central nervous system lupus.

Whole blood expression profiling from the TREAT trial: insights for the pathogenesis of polyarticular juvenile idiopathic arthritis

Background

The Trial of Early Aggressive Therapy in Juvenile Idiopathic Arthritis (TREAT trial) was accompanied by a once-in-a-generation sample collection for translational research. In this paper, we report the results of whole blood gene expression analyses and genomic data-mining designed to cast light on the immunopathogenesis of polyarticular juvenile idiopathic arthritis (JIA).

Methods

TREAT samples and samples from an independent cohort were analyzed on Affymetrix microarrays and compared to healthy controls. Data from the independent cohort were used to validate the TREAT data. Pathways analysis was used to characterize gene expression profiles. Furthermore, we correlated differential gene expression with new information about functional regulatory elements within the genome to develop models of aberrant gene expression in JIA.

Results

There was a strong concordance in gene expression between TREAT samples and the independent cohort. In addition, rheumatoid factor (RF)-positive and RF-negative patients showed only small differences on whole blood expression profiles. Analysis of the combined samples showed 158 genes represented by 176 probes that showed differential expression between TREAT subjects at baseline and healthy controls. None of the differentially expressed genes were encoded within linkage disequilibrium blocks containing single nucleotide polymorphisms known to be associated with risk for JIA. Functional analysis of these genes showed functional associations with multiple processes associated with innate and adaptive immunity, and appeared to reflect overall suppression of STAT1–3/interferon response factor-mediated pathways.

Conclusions

Despite their limitations, whole blood expression profiles clearly distinguish children with polyarticular JIA from healthy controls. Whole blood expression profiles identify several immunologic pathways of biologic relevance that will need to be pursued in homogeneous cell populations in order to clarify mechanisms of pathogenesis.

Trial registration

ClinicalTrials.gov registry #NCT00443430, originally registered 2 March 2007 and last updated 30 May 2013.

Complexity and Specificity of the Neutrophil Transcriptomes in Juvenile Idiopathic Arthritis

NIH projects such as ENCODE and Roadmap Epigenomics have revealed surprising complexity in the transcriptomes of mammalian cells. In this study, we explored transcriptional complexity in human neutrophils, cells generally regarded as nonspecific in their functions and responses. We studied distinct human disease phenotypes and found that, at the gene, gene isoform, and miRNA level, neutrophils exhibit considerable specificity in their transcriptomes. Thus, even cells whose responses are considered non-specific show tailoring of their transcriptional repertoire toward specific physiologic or pathologic contexts. We also found that miRNAs had a global impact on neutrophil transcriptome and are associated with innate immunity in juvenile idiopathic arthritis (JIA). These findings have important implications for our understanding of the link between genes, non-coding transcripts and disease phenotypes.

Dynamic tracking of functional gene modules in treated juvenile idiopathic arthritis

BACKGROUND

We have previously shown that childhood-onset rheumatic diseases show aberrant patterns of gene expression that reflect pathology-associated co-expression networks. In this study, we used novel computational approaches to examine how disease-associated networks are altered in one of the most common rheumatic diseases of childhood, juvenile idiopathic arthritis (JIA).

METHODS

Using whole blood gene expression profiles derived from children in a pediatric rheumatology clinical trial, we used a network approach to understanding the impact of therapy and the underlying biology of response/non-response to therapy.

RESULTS

We demonstrate that therapy for JIA is associated with extensive re-ordering of gene expression networks, even in children who respond inadequately to therapy. Furthermore, we observe distinct differences in the evolution of specific network properties when we compare children who have been treated successfully with those who have inadequate treatment response.

CONCLUSIONS

Despite the inherent noisiness of whole blood gene expression data, our findings demonstrate how therapeutic response might be mapped and understood in pathologically informative cells in a broad range of human inflammatory diseases.

Disease-Associated Single-Nucleotide Polymorphisms From Noncoding Regions in Juvenile Idiopathic Arthritis Are Located Within or Adjacent to Functional Genomic Elements of Human Neutrophils and CD4+ T Cells

OBJECTIVE

Juvenile idiopathic arthritis (JIA) is considered a complex disease in which the environment interacts with inherited genes to produce a phenotype that shows broad interindividual variability. Twenty-four regions of genetic risk for JIA were identified in a recent genome-wide association study (GWAS); however, as is typical of the results of GWAS, most of the regions of genetic risk (22 of 24) were in noncoding regions of the genome. This study was undertaken to identify functional elements (other than genes) that might be located within the regions of genetic risk.

METHODS

We used paired-end RNA sequencing to identify noncoding RNAs (ncRNAs) located within 5 kb of disease-associated single-nucleotide polymorphisms (SNPs). In addition, we used chromatin immunoprecipitation (ChIP) followed by sequencing to identify epigenetic marks associated with enhancer function (H3K4me1 and H3K27ac) in human neutrophils to determine whether enhancer-associated histone marks were enriched in the linkage disequilibrium blocks that encompassed the 22 SNPs identified in the GWAS.

RESULTS

In human neutrophils, we identified H3K4me1 and/or H3K27ac marks in 15 of the 22 regions previously identified as risk loci for JIA. In CD4+ T cells, 18 regions had H3K4me1 and/or H3K27ac marks. In addition, we identified ncRNA transcripts at the rs4705862 and rs6894249 loci in human neutrophils.

CONCLUSION

Much of the genetic risk for JIA lies within or adjacent to regions of neutrophil and CD4+ T cell genomes that carry epigenetic marks associated with enhancer function and/or ncRNA transcripts. These findings are consistent with the hypothesis that JIA is fundamentally a disorder of gene regulation that includes both the innate and the adaptive immune system. Elucidating the specific roles of these noncoding elements within leukocyte genomes will be critical to our understanding of JIA pathogenesis.

RNA sequencing from human neutrophils reveals distinct transcriptional differences associated with chronic inflammatory states

BACKGROUND

The transcriptional complexity of mammalian cells suggests that they have broad abilities to respond to specific environmental stimuli and physiologic contexts. These abilities were not apparent a priori from the structure of mammalian genomes, but have been identified through detailed transcriptome analyses. In this study, we examined the transcriptomes of cells of the innate immune system, human neutrophils, using RNA sequencing (RNAseq).

METHODS

We sequenced poly-A RNA from nine individual samples corresponding to specific phenotypes: three children with active, untreated juvenile idiopathic arthritis (JIA)(AD), three children with the same disease whose disease was inactive on medication (CRM), and three children with cystic fibrosis (CF).

RESULTS

We demonstrate that transcriptomes of neutrophils, typically considered non-specific in their responses and functions, display considerable specificity in their transcriptional repertoires dependent on the pathologic context, and included genes, gene isoforms, and long non-coding RNA transcripts. Furthermore, despite the small sample numbers, these findings demonstrate the potential of RNAseq approaches to biomarker development in rheumatic diseases.

CONCLUSIONS

These data demonstrate the capacity of cells previously considered non-specific in function to adapt their transcriptomes to specific biologic contexts. These data also provide insight into previously unrecognized pathological pathways and show considerable promise for elucidating disease and disease-state specific regulatory networks.

C5a alters blood-brain barrier integrity in a human in vitro model of systemic lupus erythematosus

The blood-brain barrier (BBB) plays a crucial role in brain homeostasis, thereby maintaining the brain environment precise for optimal neuronal function. Its dysfunction is an intriguing complication of systemic lupus erythematosus (SLE). SLE is a systemic autoimmune disorder where neurological complications occur in 5-50% of cases and is associated with impaired BBB integrity. Complement activation occurs in SLE and is an important part of the clinical profile. Our earlier studies demonstrated that C5a generated by complement activation caused the loss of brain endothelial layer integrity in rodents. The goal of the current study was to determine the translational potential of these studies to a human system. To assess this, we used a two dimensional in vitro BBB model constructed using primary human brain microvascular endothelial cells and astroglial cells, which closely emulates the in vivo BBB allowing the assessment of BBB integrity. Increased permeability monitored by changes in transendothelial electrical resistance and cytoskeletal remodelling caused by actin fiber rearrangement were observed when the cells were exposed to lupus serum and C5a, similar to the observations in mice. In addition, our data show that C5a/C5aR1 signalling alters nuclear factor-κB translocation into nucleus and regulates the expression of the tight junction proteins, claudin-5 and zonula occludens 1 in this setting. Our results demonstrate for the first time that C5a regulates BBB integrity in a neuroinflammatory setting where it affects both endothelial and astroglial cells. In addition, we also demonstrate that our previous findings in a mouse model, were emulated in human cells in vitro, bringing the studies one step closer to understanding the translational potential of C5a/C5aR1 blockade as a promising therapeutic strategy in SLE and other neurodegenerative diseases.

Examination of parent-child adjustment in juvenile rheumatic diseases using depression-specific indices of parent and youth functioning

Studies demonstrate a link between parental distress, youth illness appraisals, and depression symptoms in youth with juvenile rheumatic diseases. However, the exclusive use of broadband (i.e. general) measures of parental distress in these studies has resulted in conceptual and clinical imprecision regarding the parent-child adjustment process. Our aim was to reanalyze previously published data (i.e. Wagner et al., 2003) using a depression-specific scale derived from the general adult distress measure in the original study. Parents completed the Brief Symptom Inventory (BSI), youth completed the Child Depression Inventory (CDI), and the Illness Intrusiveness Scale (IIS-C). Thirteen Diagnostic and Statistical Manual of Mental Disorders (Fourth Edition) depression-specific items from the BSI comprised the parent measure of Depressive Symptoms Scale (DS). Consistent with Wagner et al. (2003), adult DS scores were associated with youth CDI scores. However, youth illness appraisals had unique effects on the parent-child depression relation. Elevated youth perceptions of illness intrusiveness amplified the negative effect of parent depressive symptoms on youth depressive symptoms; decreased illness intrusiveness buffered the negative effect of parent depression. The empirical and clinical implications of assessing parent and youth adjustment in a domain-specific manner are discussed.

Whole blood gene expression profiling predicts therapeutic response at six months in patients with polyarticular juvenile idiopathic arthritis

Objective

To determine whether gene expression profiles identified in peripheral whole blood samples could be used to determine therapeutic outcome in a cohort of children with newly diagnosed polyarticular juvenile idiopathic arthritis (JIA).

Methods

Whole blood samples from the Trial of Early Aggressive Therapy (TREAT) in JIA patients were analyzed on Illumina microarrays, and differential gene expression was compared to expression in healthy controls. Microarray results were validated by real-time quantitative polymerase chain reaction in an independent cohort of samples. Pathway analysis software was used to characterize gene expression profiles. Support vector machines were used to develop predictive models for different patient classes.

Results

Differential gene expression profiles for rheumatoid factor (RF)–positive and RF-negative patients were remarkably similar. Pathway analysis revealed a broad range of affected pathways, consistent with current mechanistic theories. Modeling showed that the prognosis at 6 months was strongly linked to gene expression at presentation, irrespective of treatment.

Conclusion

Gene expression is linked to therapeutic outcome, and gene expression in the peripheral blood may be a suitable target for a prognostic test.

Caregiver demand and parent distress in juvenile rheumatic disease: the mediating effect of parent attitude toward illness

Parents of youth with juvenile rheumatic diseases (JRD) often take on illness management responsibilities that can become burdensome, potentially resulting in poor parent adjustment outcomes. However, not all caregivers will experience increased distress as a result of variability in stress appraisals. The current study examined the role of parent illness attitudes in the relation between perceived caregiver demand and parental distress. Youth (N = 70) ages 7-18 years diagnosed with a JRD and their parents were recruited from a pediatric rheumatology clinic. Parents completed measures of caregiver demand, parental distress, and illness attitudes. Hierarchical regression revealed a relationship between caregiver demand and parental distress. A significant relationship was also found between caregiver demand and parent illness attitudes, as well as parent illness attitudes and parental distress. Thus, parent illness attitudes mediated the relationship between caregiver demand and parental distress. Techniques aimed at altering negative illness attitudes may help parents cope with their caregiving responsibilities.

Diagnostic and prognostic potential of gene microarrays in rheumatoid arthritis

Rheumatoid arthritis and juvenile rheumatoid arthritis are histopathologically similar diseases characterized by chronic inflammation of the synovium. The pathogenesis of these diseases is unknown, but the emergence of gene expression profiling provides considerable promise that some of the complex, interconnected immunopathologic events underlying these diseases will soon be better understood. This review will summarize the potential use of gene expression profiling as a diagnostic or prognostic modality, and the potential benefits or limits of such uses. It will conclude with a short discussion of the potential for using gene expression profiling to identify novel targets of therapy in rheumatoid arthritis and related diseases.

Starting ART following cryptococcal meningitis:The optimal time has yet to be defined

Ever since the public sector rollout of antiretroviral therapy (ART) in 2004, the question of the optimal time to start ART following diagnosis of an opportunistic infection has aroused controversy among South African HIV clinicians and researchers.

Genomic characterization of remission in juvenile idiopathic arthritis

Introduction

The attainment of remission has become an important end point for clinical trials in juvenile idiopathic arthritis (JIA), although we do not yet have a full understanding of what remission is at the cell and molecular level.

Methods

Two independent cohorts of patients with JIA and healthy child controls were studied. RNA was prepared separately from peripheral blood mononuclear cells (PBMC) and granulocytes to identify differentially expressed genes using whole genome microarrays. Expression profiling results for selected genes were confirmed by quantitative, real-time polymerase chain reaction (RT-PCR).

Results

We found that remission in JIA induced by either methotrexate (MTX) or MTX plus a TNF inhibitor (etanercept, Et) (MTX + Et) is characterized by numerous differences in gene expression in peripheral blood mononuclear cells and in granulocytes compared with healthy control children; that is, remission is not a restoration of immunologic normalcy. Network analysis of the differentially expressed genes demonstrated that the steroid hormone receptor superfamily member hepatocyte nuclear factor 4 alpha (HNF4α) is a hub in several of the gene networks that distinguished children with arthritis from controls. Confocal microscopy revealed that HNF4a is present in both T lymphocytes and granulocytes, suggesting a previously unsuspected role for this transcription factor in regulating leukocyte function and therapeutic response in JIA.

Conclusions

These findings provide a framework from which to understand therapeutic response in JIA and, furthermore, may be used to develop strategies to increase the frequency with which remission is achieved in adult forms of rheumatoid arthritis.

Utility of the Illness Intrusiveness Scale in parents of children diagnosed with juvenile rheumatic diseases

OBJECTIVE

To examine the factor structure and convergent validity of the Illness Intrusiveness Scale--Parent Version in mother and fathers of children and adolescents ages 7 to 18 (M = 13.56 years, SD = 2.67) diagnosed with a juvenile rheumatic disease.

DESIGN

Parents of 122 children and adolescents (82 girls, 40 boys) completed the Illness Intrusiveness Scale-Parent Version, and both parents and children and adolescents completed measures of functional disability, general distress, and illness uncertainty. An exploratory factor analysis was conducted on the Illness Intrusiveness Scale--Parent Version to identify the factor structure. The factors were then compared with parent- and child-report measures of functional disability, general distress, and uncertainty. Finally, analyses were conducted to determine whether the magnitude of the correlations was significantly different between factors for parents and children and adolescents.

RESULTS

The Illness Intrusiveness Scale--Parent Version was found to have a two-factor structure. The Relationships/Personal Development factor contained items related to self-fulfillment and interactions with others, and the Instrumental factor contained items related to health and work. These factors were found to have good internal consistency and were significantly correlated with measures of parent-reported functional disability and parent- and youth-reported distress and uncertainty. The magnitude of these correlations was also found to differ depending on informant and outcome measure.

CONCLUSION

The Illness Intrusiveness Scale--Parent Version appears to be a valid measure for use in parents of children with juvenile rheumatic disease.

Functional genomics and rheumatoid arthritis: where have we been and where should we go?

Studies in model organisms and humans have begun to reveal the complexity of the transcriptome. In addition to serving as passive templates from which genes are translated, RNA molecules are active, functional elements of the cell whose products can detect, interact with, and modify other transcripts. Gene expression profiling is the method most commonly used thus far to enrich our understanding of the molecular basis of rheumatoid arthritis in adults and juvenile idiopathic arthritis in children. The feasibility of this approach for patient classification (for example, active versus inactive disease, disease subsets) and improving prognosis (for example, response to therapy) has been demonstrated over the past 7 years. Mechanistic understanding of disease-related differences in gene expression must be interpreted in the context of interactions with transcriptional regulatory molecules and epigenetic alterations of the genome. Ongoing work regarding such functional complexities in the human genome will likely bring both insight and surprise to our understanding of rheumatoid arthritis.

Ablation of Stat3 by siRNA alters gene expression profiles in JEG-3 cells: a systems biology approach

Control of inflammation at the maternal-fetal interface is a critical element in mammalian pregnancy. Previous work from our laboratory has shown that Stat3 may be a placental mediator involved in maintaining immunologic homeostasis at the maternal-fetal interface. The aim of the current study is to further elucidate the role of Stat3 in response to inflammation. As ablation of Stat3 in mice results in embryonic lethality, we evaluated the role of Stat3 in vitro using an siRNA approach. Trophoblast-like JEG-3 cells were transfected with an siRNA construct specific to Stat3. Experimental and control cells were exposed to conditioned medium from PHA-activated peripheral blood mononuclear cells and incubated for 45 min. Cells were then collected and RNA isolated for transcriptional profiling using human Affymetrix U133 plus 2.0 GeneChips. Differences in gene expression between control and Stat3-ablated cells were evaluated using conventional statistical methods. Fifty-two genes were detected as up-regulated in conditioned medium in both mock transfected and in Stat3 siRNA transfected JEG-3 cells. Two genes (EPAS1 and RASGEF1B) were up-regulated only in cells transfected with negative control siRNA, while 36 genes were up-regulated only in cells transfected with Stat3 siRNA. Sixty genes were differentially expressed between Stat3 siRNA transfected cells relative to mock transfected cells both in basal and conditioned medium. These included 31 genes up-regulated with Stat3 siRNA transfected cells and 29 genes down-regulated with Stat3 siRNA. Eleven genes were differentially expressed only in basal medium. Seven of these were up-regulated in the presence of Stat3 siRNA and four were down-regulated. Nine genes were differentially expressed only in conditioned medium. Six of these were up-regulated and three down-regulated in the presence of Stat3 siRNA. Off-target effects were excluded in a second set of experiments in which Stat3 mRNA was targeted at a different site and quantitative real-time PCR performed on selected genes derived from the microarray analysis. While some of the genes that showed differential expression between Stat3-ablated cells and mock transfected controls were genes typically associated with immune response (e.g., CCR7 and IRAK1), in silico modeling of the microarray data also revealed complex networks of signaling molecules and molecules associated with cellular metabolism previously seen in transcription factor ablation in model organisms. We conclude thus: Stat3 controls a specific gene set in trophoblast-like JEG-3 cells. While some differentially expressed genes and in silico models of their functions are consistent with the hypothesis that Stat3 plays a role in regulating inflammation, Stat3-mediated response to inflammation appears to also involve complex homeostatic adaptations of a non-immunologic nature.

Gene expression profiling in neutrophils from children with polyarticular juvenile idiopathic arthritis

OBJECTIVE

We have previously reported a defect in neutrophil activation in children with polyarticular juvenile idiopathic arthritis (JIA). The current study was undertaken to determine whether gene expression abnormalities persist in JIA in remission and to use systems biology analysis to elucidate pathologic pathways in polyarticular JIA.

METHODS

We performed gene expression profiling on neutrophils from children with polyarticular JIA. Children were grouped according to disease status. We studied 14 children with active disease who were taking medication, 8 children with clinical remission of disease who were taking medication (CRM status), and 6 children with clinical remission of disease who were not taking medication (CR status). We also studied 13 healthy children whose age ranges overlapped those of the patients.

RESULTS

Neutrophil abnormalities persisted in children with polyarticular JIA even after disease remission was achieved. Children with active disease and those with CRM status showed no differences in expression of specific genes, although they could be separated on cluster analysis. A comparison of children with CR status and healthy control children revealed networks of pro- and antiinflammatory genes that suggested that remission is a state of homeostasis and balance rather than a return to normal immune function. Furthermore, gene overexpression in patients with CR status supports the hypothesis that neutrophils play a role in regulating adaptive immunity in this disease.

CONCLUSION

Neutrophil gene profiling in polyarticular JIA suggests important roles for neutrophils in disease pathogenesis. These findings suggest the presence of complex interactions between innate and adaptive immunity, that are not easily modeled in conventional, linear, reductionist systems.