Genomic Identification of Low-Density Granulocytes and Analysis of Their Role in the Pathogenesis of Systemic Lupus Erythematosus.

Systemic lupus erythematosus (SLE) is an autoimmune disease characterized by the presence of low-density granulocytes (LDGs) with a heightened capacity for spontaneous NETosis, but the contribution of LDGs to SLE pathogenesis remains unclear. To characterize LDGs in human SLE, gene expression profiles derived from isolated LDGs were characterized by weighted gene coexpression network analysis, and a 92-gene module was identified. The LDG gene signature was enriched in genes related to neutrophil degranulation and cell cycle regulation. This signature was assessed in gene expression datasets from two large-scale SLE clinical trials to study associations between LDG enrichment, SLE manifestations, and treatment regimens. LDG enrichment in the blood was associated with corticosteroid treatment as well as anti-dsDNA, low serum complement, renal manifestations, and vasculitis, but the latter two of these associations were dependent on concomitant corticosteroid treatment. In addition, LDG enrichment was associated with enrichment of gene signatures induced by type I IFN and TNF irrespective of corticosteroid treatment. Notably, LDG enrichment was not found in numerous tissues affected by SLE. Comparison with relevant reference datasets indicated that LDG enrichment is likely reflective of increased granulopoiesis in the bone marrow and not peripheral neutrophil activation. The results have uncovered important determinants of the appearance of LDGs in SLE and have emphasized the likely role of LDGs in specific aspects of lupus pathogenesis.

Machine learning approaches to predict lupus disease activity from gene expression data.

The integration of gene expression data to predict systemic lupus erythematosus (SLE) disease activity is a significant challenge because of the high degree of heterogeneity among patients and study cohorts, especially those collected on different microarray platforms. Here we deployed machine learning approaches to integrate gene expression data from three SLE data sets and used it to classify patients as having active or inactive disease as characterized by standard clinical composite outcome measures. Both raw whole blood gene expression data and informative gene modules generated by Weighted Gene Co-expression Network Analysis from purified leukocyte populations were employed with various classification algorithms. Classifiers were evaluated by 10-fold cross-validation across three combined data sets or by training and testing in independent data sets, the latter of which amplified the effects of technical variation. A random forest classifier achieved a peak classification accuracy of 83 percent under 10-fold cross-validation, but its performance could be severely affected by technical variation among data sets. The use of gene modules rather than raw gene expression was more robust, achieving classification accuracies of approximately 70 percent regardless of how the training and testing sets were formed. Fine-tuning the algorithms and parameter sets may generate sufficient accuracy to be informative as a standalone estimate of disease activity.

Potential strain-dependent mechanisms defining matrix alignment in healing tendons.

Tendon mechanical function after injury and healing is largely determined by its underlying collagen structure, which in turn is dependent on the degree of mechanical loading experienced during healing. Experimental studies have shown seemingly conflicting outcomes: although collagen content steadily increases with increasing loads, collagen alignment peaks at an intermediate load. Herein, we explored potential collagen remodeling mechanisms that could give rise to this structural divergence in response to strain. We adapted an established agent-based model of collagen remodeling in order to simulate various strain-dependent cell and collagen interactions that govern long-term collagen content and fiber alignment. Our simulation results show two collagen remodeling mechanisms that give rise to divergent collagen content and alignment in healing tendons: (1) strain-induced collagen fiber damage in concert with increased rates of deposition at higher strains, or (2) strain-dependent rates of enzymatic degradation. These model predictions identify critical future experiments needed to isolate each mechanism's specific contribution to the structure of healing tendons.

Identification of alterations in macrophage activation associated with disease activity in systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is characterized by abnormalities in B cell and T cell function, but the role of disturbances in the activation status of macrophages (Mϕ) has not been well described in human patients. To address this, gene expression profiles from isolated lymphoid and myeloid populations were analyzed to identify differentially expressed (DE) genes between healthy controls and patients with either inactive or active SLE. While hundreds of DE genes were identified in B and T cells of active SLE patients, there were no DE genes found in B or T cells from patients with inactive SLE compared to healthy controls. In contrast, large numbers of DE genes were found in myeloid cells (MC) from both active and inactive SLE patients. Among the DE genes were several known to play roles in Mϕ activation and polarization, including the M1 genes STAT1 and SOCS3 and the M2 genes STAT3, STAT6, and CD163. M1-associated genes were far more frequent in data sets from active versus inactive SLE patients. To characterize the relationship between Mϕ activation and disease activity in greater detail, weighted gene co-expression network analysis (WGCNA) was used to identify modules of genes associated with clinical activity in SLE patients. Among these were disease activity-correlated modules containing activation signatures of predominantly M1-associated genes. No disease activity-correlated modules were enriched in M2-associated genes. Pathway and upstream regulator analysis of DE genes from both active and inactive SLE MC were cross-referenced with high-scoring hits from the drug discovery Library of Integrated Network-based Cellular Signatures (LINCS) to identify new strategies to treat both stages of SLE. A machine learning approach employing MC gene modules and a generalized linear model was able to predict the disease activity status in unrelated gene expression data sets. In summary, altered MC gene expression is characteristic of both active and inactive SLE. However, disease activity is associated with an alteration in the activation of MC, with a bias toward the M1 proinflammatory phenotype. These data suggest that while hyperactivity of B cells and T cells is associated with active SLE, MC potentially direct flare-ups and remission by altering their activation status toward the M1 state.