Distinct autoantibody profiles across checkpoint inhibitor types and toxicities

Immune checkpoint inhibitors (ICI) are increasingly used in combination. To understand the effects of different ICI categories, we characterized changes in circulating autoantibodies in patients enrolled in the E4412 trial (NCT01896999) of brentuximab vedotin (BV) plus ipilimumab, BV plus nivolumab, or BV plus ipilimumab-nivolumab for Hodgkin Lymphoma. Cycle 2 Day 1 (C2D1) autoantibody levels were compared to pre-treatment baseline. Across 112 autoantibodies tested, we generally observed increases in ipilimumab-containing regimens, with decreases noted in the nivolumab arm. Among 15 autoantibodies with significant changes at C2D1, all nivolumab cases exhibited decreases, with more than 90% of ipilimumab-exposed cases showing increases. Autoantibody profiles also showed differences according to immune-related adverse event (irAE) type, with rash generally featuring increases and liver toxicity demonstrating decreases. We conclude that dynamic autoantibody profiles may differ according to ICI category and irAE type. These findings may have relevance to clinical monitoring and irAE treatment.

Systemic immune parameters after prior radiation therapy in patients receiving immune checkpoint inhibitors

INTRODUCTION

Preclinical studies have demonstrated the ability of radiation therapy (RT) to augment immune response and tumor control by immune checkpoint inhibitors (ICI). However, numerous clinical trials combining RT and ICI have yielded relatively disappointing results. To improve understanding of optimal use of these therapies, we assessed systemic immune effects of prior RT in patients receiving ICI.

METHODS AND MATERIALS

Pre- and post-ICI blood samples were collected from patients enrolled in a prospective immunotherapy biospecimen protocol. Mutiplex panels of 40 cytokines and 120 autoantibodies (Ab) were analyzed. We identified differences in these parameters according to receipt, timing, and type of prior RT. We calculated P values using the Pearson product-moment correlation coefficient and false discovery rate (FDR) using the Benjamini-Hochberg Procedure.

RESULTS

Among 277 total patients, 69 (25%) received RT in the 6 months prior to ICI initiation. Among RT-treated patients, 23 (33%) received stereotactic RT, and 33 (48%) received curative intent RT. There was no significant difference in demographics or type of immunotherapy between patients according to prior RT exposure. Baseline complement C8 Ab and MIP-1d/CCL15 were significantly higher among patients with prior RT. For MIP-1d/CCL15, only prior stereotactic RT was associated with significant differences.

CONCLUSIONS

Prior RT is associated with few changes in systemic immune parameters in patients receiving ICI. The underlying mechanisms and optimal approach to harnessing the potential synergy of RT and ICI require further prospective clinical investigation.

T-cell tolerant fraction as a predictor of immune-related adverse events

BACKGROUND

Immune checkpoint inhibitor (ICI) therapies may cause unpredictable and potentially severe autoimmune toxicities termed immune-related adverse events (irAEs). Because T cells mediate ICI effects, T cell profiling may provide insight into the risk of irAEs. Here we evaluate a novel metric-the T-cell tolerant fraction-as a predictor of future irAEs.

METHODS

We examined T-cell receptor beta (TRB) locus sequencing from baseline pretreatment samples from an institutional registry and previously published studies. For each patient, we used TRB sequences to calculate the T-cell tolerant fraction, which was then assessed as a predictor of future irAEs (classified as Common Terminology Criteria for Adverse Event grade 0-1 vs grade ≥2). We then compared the tolerant fraction to TRB clonality and diversity. Finally, the tolerant fraction was assessed on (1) T cells enriched against napsin A, a potential autoantigen of irAEs; (2) thymic versus peripheral blood T cells; and (3) TRBs specific for various infections and autoimmune diseases.

RESULTS

A total of 77 patients with cancer (22 from an institutional registry and 55 from published studies) receiving ICI therapy (43 CTLA4, 19 PD1/PDL1, 15 combination CTLA4+PD1/PDL1) were included in the study. The tolerant fraction was significantly lower in cases with clinically significant irAEs (p<0.001) and had an area under the receiver operating curve (AUC) of 0.79. The tolerant fraction was lower for each ICI treatment category, reaching statistical significance for CTLA4 (p<0.001) and demonstrating non-significant trends for PD1/PDL1 (p=0.21) and combination ICI (p=0.18). The tolerant fraction for T cells enriched against napsin A was lower than other samples. The tolerant fraction was also lower in thymic versus peripheral blood samples, and lower in some (multiple sclerosis) but not other (type 1 diabetes) autoimmune diseases. In our study cohort, TRB clonality had an AUC of 0.62, and TRB diversity had an AUC of 0.60 for predicting irAEs.

CONCLUSIONS

Among patients receiving ICI, the baseline T-cell tolerant fraction may serve as a predictor of clinically significant irAEs.

Loss-of-function variants in SAT1 cause X-linked childhood-onset systemic lupus erythematosus

OBJECTIVES

Families that contain multiple siblings affected with childhood onset of systemic lupus erythematosus (SLE) likely have strong genetic predispositions. We performed whole exome sequencing (WES) to identify familial rare risk variants and to assess their effects in lupus.

METHODS

Sanger sequencing validated the two ultra-rare, predicted pathogenic risk variants discovered by WES and identified additional variants in 562 additional patients with SLE. Effects of a splice site variant and a frameshift variant were assessed using a Minigene assay and CRISPR/Cas9-mediated knock-in (KI) mice, respectively.

RESULTS

The two familial ultra-rare, predicted loss-of-function (LOF) SAT1 variants exhibited X-linked recessive Mendelian inheritance in two unrelated African-American families. Each LOF variant was transmitted from the heterozygous unaffected mother to her two sons with childhood-onset SLE. The p.Asp40Tyr variant affected a splice donor site causing deleterious transcripts. The young hemizygous male and homozygous female Sat1 p.Glu92Leufs*6 KI mice spontaneously developed splenomegaly, enlarged glomeruli with leucocyte infiltration, proteinuria and elevated expression of type I interferon-inducible genes. SAT1 is highly expressed in neutrophils and encodes spermidine/spermine-N1-acetyltransferase 1 (SSAT1), a rate-limiting enzyme in polyamine catabolism. Young male KI mice exhibited neutrophil defects and decreased proportions of Foxp3 +CD4+ T-cell subsets. Circulating neutrophil counts and proportions of Foxp3 +CD4+ T cells correlated with decreased plasma levels of spermine in treatment-naive, incipient SLE patients.

CONCLUSIONS

We identified two novel SAT1 LOF variants, showed the ability of the frameshift variant to confer murine lupus, highlighted the pathogenic role of dysregulated polyamine catabolism and identified SAT1 LOF variants as new monogenic causes for SLE.

Association between Antibiotic Exposure and Systemic Immune Parameters in Cancer Patients Receiving Checkpoint Inhibitor Therapy

Antibiotic administration is associated with worse clinical outcomes and changes to the gut microbiome in cancer patients receiving immune checkpoint inhibitors (ICI). However, the effects of antibiotics on systemic immune function are unknown. We, therefore, evaluated antibiotic exposure, therapeutic responses, and multiplex panels of 40 serum cytokines and 124 antibodies at baseline and six weeks after ICI initiation, with p < 0.05 and false discovery rate (FDR) < 0.2 considered significant. A total of 251 patients were included, of whom the 135 (54%) who received antibiotics had lower response rates and shorter survival. Patients who received antibiotics prior to ICI initiation had modestly but significantly lower baseline levels of nucleolin, MDA5, c-reactive protein, and liver cytosol antigen type 1 (LC1) antibodies, as well as higher levels of heparin sulfate and Matrigel antibodies. After ICI initiation, antibiotic-treated patients had significantly lower levels of MDA5, CENP.B, and nucleolin antibodies. Although there were no clear differences in cytokines in the overall cohort, in the lung cancer subset (53% of the study population), we observed differences in IFN-γ, IL-8, and macrophage inflammatory proteins. In ICI-treated patients, antibiotic exposure is associated with changes in certain antibodies and cytokines. Understanding the relationship between these factors may improve the clinical management of patients receiving ICI.

Outcome and Immune Correlates of a Phase II Trial of High-Dose Interleukin-2 and Stereotactic Ablative Radiotherapy for Metastatic Renal Cell Carcinoma

PURPOSE

This phase II clinical trial evaluated whether the addition of stereotactic ablative radiotherapy (SAbR), which may promote tumor antigen presentation, improves the overall response rate (ORR) to high-dose IL2 (HD IL2) in metastatic renal cell carcinoma (mRCC).

PATIENTS AND METHODS

Patients with pathologic evidence of clear cell renal cell carcinoma (RCC) and radiographic evidence of metastasis were enrolled in this single-arm trial and were treated with SAbR, followed by HD IL2. ORR was assessed based on nonirradiated metastases. Secondary endpoints included overall survival (OS), progression-free survival (PFS), toxicity, and treatment-related tumor-specific immune response. Correlative studies involved whole-exome and transcriptome sequencing, T-cell receptor sequencing, cytokine analysis, and mass cytometry on patient samples.

RESULTS

Thirty ethnically diverse mRCC patients were enrolled. A median of two metastases were treated with SAbR. Among 25 patients evaluable by RECIST v1.1, ORR was 16% with 8% complete responses. Median OS was 37 months. Treatment-related adverse events (AE) included 22 grade ≥3 events that were not dissimilar from HD IL2 alone. There were no grade 5 AEs. A correlation was observed between SAbR to lung metastases and improved PFS (P = 0.0165). Clinical benefit correlated with frameshift mutational load, mast cell tumor infiltration, decreased circulating tumor-associated T-cell clones, and T-cell clonal expansion. Higher regulatory/CD8⁺ T-cell ratios at baseline in the tumor and periphery correlated with no clinical benefit.

CONCLUSIONS

Adding SAbR did not improve the response rate to HD IL2 in patients with mRCC in this study. Tissue analyses suggest a possible correlation between frameshift mutation load as well as tumor immune infiltrates and clinical outcomes.

Humoral and cellular correlates of a novel immune-related adverse event and its treatment

Immune-related adverse events (irAE) may affect almost any organ system and occur at any point during treatment with immune checkpoint inhibitors (ICI). We present a patient with advanced lung cancer receiving antiprogrammed death 1 checkpoint inhibitor who developed a delayed-onset visual irAE treated with corticosteroids. Through assessment of longitudinal biospecimens, we analyzed serial autoantibodies, cytokines, and cellular populations. Months after ICI initiation and preceding clinical toxicity, the patient developed broad increases in cytokines (most notably interleukin-6 (IL-6), interferon-γ (IFNγ), C-X-C motif chemokine ligand 2 (CXCL2), and C-C motif chemokine ligand 17 (CCL17)), autoantibodies (including anti-angiotensin receptor, α-actin, and amyloid), CD8 T cells, and plasmablasts. Such changes were not observed in healthy controls and ICI-treated patients without irAE. Administration of corticosteroids resulted in immediate and profound decreases in cytokines, autoantibodies, and inflammatory cells. This case highlights the potential for late-onset changes in humoral and cellular immunity in patients receiving ICI. It also demonstrates the biologic effects of corticosteroids on these parameters. Application of humoral and cellular immune biomarkers across ICI populations may inform toxicity monitoring and management.

Association between body mass index, dosing strategy, and efficacy of immune checkpoint inhibitors

Background

Increased body mass index (BMI) has been associated with improved response to immune checkpoint inhibitors (ICIs) in multiple cancer types. We evaluated associations between BMI, ICI dosing strategy, and clinical outcomes.

Methods

We abstracted clinical data on patients with cancer treated with ICI, including age, sex, cancer type, BMI, ICI type, dosing strategy (weight-based or fixed), radiographic response, overall survival (OS), and progression-free survival (PFS). We compared clinical outcomes between low-BMI and high-BMI populations using Kaplan-Meier curves, Cox regressions, and Pearson product-moment correlation coefficients.

Results

A total of 297 patients were enrolled, of whom 40% were women and 59% were overweight (BMI≥25). Of these, 204 (69%) received fixed and 93 (31%) received weight-based ICI dosing. In the overall cohort, overweight BMI was associated with improved PFS (HR 0.69; 95% CI 0.51 to 0.94; p=0.02) and had a trend toward improved OS (HR 0.77; 95% CI 0.57 to 1.04; p=0.08). For both endpoints, improved outcomes in the overweight population were limited to patients who received weight-based ICI dosing (PFS HR 0.53; p=0.04 for weight-based; vs HR 0.79; p=0.2 for fixed dosing) (OS HR 0.56; p=0.03 for weight-based; vs HR 0.89; p=0.54 for fixed dosing). In multivariable analysis, BMI was not associated with PFS or OS. However, the interaction of BMI≥25 and weight-based dosing had a trend toward association with PFS (HR 0.53; 95% CI 0.26 to 1.10; p=0.09) and was associated with OS (HR 0.50; 95% CI 0.25 to 0.99; p=0.05). Patients with BMI<25 tended to have better outcomes with fixed-dose compared with weight-based ICI, while patients with BMI≥25 tended to have better outcomes with weight-based ICI, although these differences did not achieve statistical significance. There was no association between radiographic response and BMI with fixed-dose ICI (p=0.97), but a near-significant trend with weight-based ICI (p=0.1). In subset analyses, the association between BMI, ICI dosing strategy, and clinical outcomes appeared limited to men.

Conclusions

The clinical benefit of ICI in high-BMI populations appears limited to individuals receiving weight-based ICI dosing. Further research into optimal ICI dosing strategies may be warranted.

IRF1 governs the differential interferon-stimulated gene responses in human monocytes and macrophages by regulating chromatin accessibility

Myeloid lineage cells use TLRs to recognize and respond to diverse microbial ligands. Although unique transcription factors dictate the outcome of specific TLR signaling, whether lineage-specific differences exist to further modulate the quality of TLR-induced inflammation remains unclear. Comprehensive analysis of global gene transcription in human monocytes, monocyte-derived macrophages, and monocyte-derived dendritic cells stimulated with various TLR ligands identifies multiple lineage-specific, TLR-responsive gene programs. Monocytes are hyperresponsive to TLR7/8 stimulation that correlates with the higher expression of the receptors. While macrophages and monocytes express similar levels of TLR4, macrophages, but not monocytes, upregulate interferon-stimulated genes (ISGs) in response to TLR4 stimulation. We find that TLR4 signaling in macrophages uniquely engages transcription factor IRF1, which facilitates the opening of ISG loci for transcription. This study provides a critical mechanistic basis for lineage-specific TLR responses and uncovers IRF1 as a master regulator for the ISG transcriptional program in human macrophages.

Deep sequencing reveals a DAP1 regulatory haplotype that potentiates autoimmunity in systemic lupus erythematosus

BACKGROUND

Systemic lupus erythematosus (SLE) is a clinically heterogeneous autoimmune disease characterized by the development of anti-nuclear antibodies. Susceptibility to SLE is multifactorial, with a combination of genetic and environmental risk factors contributing to disease development. Like other polygenic diseases, a significant proportion of estimated SLE heritability is not accounted for by common disease alleles analyzed by SNP array-based GWASs. Death-associated protein 1 (DAP1) was implicated as a candidate gene in a previous familial linkage study of SLE and rheumatoid arthritis, but the association has not been explored further.

RESULTS

We perform deep sequencing across the DAP1 genomic segment in 2032 SLE patients, and healthy controls, and discover a low-frequency functional haplotype strongly associated with SLE risk in multiple ethnicities. We find multiple cis-eQTLs embedded in a risk haplotype that progressively downregulates DAP1 transcription in immune cells. Decreased DAP1 transcription results in reduced DAP1 protein in peripheral blood mononuclear cells, monocytes, and lymphoblastoid cell lines, leading to enhanced autophagic flux in immune cells expressing the DAP1 risk haplotype. Patients with DAP1 risk allele exhibit significantly higher autoantibody titers and altered expression of the immune system, autophagy, and apoptosis pathway transcripts, indicating that the DAP1 risk allele mediates enhanced autophagy, leading to the survival of autoreactive lymphocytes and increased autoantibody.

CONCLUSIONS

We demonstrate how targeted sequencing captures low-frequency functional risk alleles that are missed by SNP array-based studies. SLE patients with the DAP1 genotype have distinct autoantibody and transcription profiles, supporting the dissection of SLE heterogeneity by genetic analysis.

Transcriptional profiling identifies caspase-1 as a T cell-intrinsic regulator of Th17 differentiation

Dendritic cells (DCs) are critical for the differentiation of pathogen-specific CD4 T cells. However, to what extent innate cues from DCs dictate transcriptional changes in T cells remains elusive. Here, we used DCs stimulated with specific pathogens to prime CD4 T cells in vitro and found that these T cells express unique transcriptional profiles dictated by the nature of the priming pathogen. More specifically, the transcriptome of in vitro C. rodentium–primed Th17 cells resembled that of Th17 cells primed following infection in vivo but was remarkably distinct from cytokine-polarized Th17 cells. We identified caspase-1 as a unique gene up-regulated only in pathogen-primed Th17 cells and discovered a critical role for T cell–intrinsic caspase-1, independent of inflammasome, in optimal priming of Th17 responses. T cells lacking caspase-1 failed to induce colitis or confer protection against C. rodentium infection due to suboptimal Th17 cell differentiation in vivo. This study underlines the importance of DC-mediated priming in identifying novel regulators of T cell differentiation.

sncRNA-1 Is a Small Noncoding RNA Produced by Mycobacterium tuberculosis in Infected Cells That Positively Regulates Genes Coupled to Oleic Acid Biosynthesis

Nearly one third of the world’s population is infected with Mycobacterium tuberculosis (Mtb). While much work has focused on the role of different Mtb encoded proteins in pathogenesis, recent studies have revealed that Mtb also transcribes many noncoding RNAs whose functions remain poorly characterized. We performed RNA sequencing and identified a subset of Mtb H37Rv-encoded small RNAs (<30 nts in length) that were produced in infected macrophages. Designated as smaller noncoding RNAs (sncRNAs), three of these predominated the read counts. Each of the three, sncRNA-1, sncRNA-6, and sncRNA-8 had surrounding sequences with predicted stable secondary RNA stem loops. Site-directed mutagenesis of the precursor sequences suggest the existence of a hairpin loop dependent RNA processing mechanism. A functional assessment of sncRNA-1 suggested that it positively regulated two mycobacterial transcripts involved in oleic acid biosynthesis. Complementary loss- and gain- of-function approaches revealed that sncRNA-1 positively supports Mtb growth and survival in nutrient-depleted cultures as well as in infected macrophages. Overall, the findings reveal that Mtb produces sncRNAs in infected cells, with sncRNA-1 modulating mycobacterial gene expression including genes coupled to oleic acid biogenesis.

Statin Intolerance, Anti-HMGCR Antibodies, and Immune Checkpoint Inhibitor-Associated Myositis: A "Two-Hit" Autoimmune Toxicity or Clinical Predisposition?

Immune-related adverse events induced by immune checkpoint inhibitor (ICI) therapy may affect diverse organ systems, including skeletal and cardiac muscle. ICI-associated myositis may result in substantial morbidity and occasional mortality. We present a case of a patient with advanced non-small cell lung cancer who developed grade 4 myositis with concurrent myocarditis early after initiation of anti-programmed death ligand 1 therapy (durvalumab). Autoantibody analysis revealed marked increases in anti-3-hydroxy-3-methylglutaryl-coenzyme A reductase antibody levels that preceded clinical toxicity, and further increased during toxicity. Notably, the patient had a history of intolerable statin myopathy, which had resolved clinically after statin discontinuation and prior to ICI initiation. This case demonstrates a potential association between statin exposure, autoantibodies, and ICI-associated myositis.

Longitudinal Immune and Genomic Monitoring Reveals Signatures of Immune-related Adverse Events in Cancer Patients Receiving Checkpoint Inhibitor Therapy

Despite the remarkable success of immune checkpoint inhibitor (ICI) therapy, a significant number of patients develop severe and unpredictable immune-related adverse events (irAEs) affecting a wide variety of organs. Concerns over irAE have led to the exclusion of patients with autoimmune disease from ICI clinical trials. Role of host genetic and immune factors in mediating irAEs remain unclear and it is not clear if the manifestations of irAEs is associated with response to therapy. Here, we use multi-faceted approach to characterize changes in host immune system in 200 patients receiving ICI therapy at baseline and post immunotherapy. In addition, we assessed genetic predisposition to autoimmunity using the Illumina SNP array and via targeted resequencing of over 150 immunoregulatory loci including the HLA region. In this meeting, we will present our initial genetic data, serum cytokine and autoantibody profiles and RNA sequencing on peripheral blood mononuclear cells (PBMCs) at baseline and post immunotherapy in patients with and without irAEs. Our preliminary findings suggest that patients who developed irAEs have lower baseline levels and greater post-treatment increases in key interferon gamma inducible cytokine/chemokine levels. Our autoantibody profiling data reveals unique sets of autoantibodies associated with specific irAEs. We will present a comprehensive analysis of immune and genetic correlates of irAEs and response to therapy. We hope our study can help gain insight in to the mechanisms underlying irAE and to identify biomarker signatures predictive of irAEs and/or response. These findings may ultimately help identify high-risk patients, customize therapy, expand use of immunotherapy and prevent toxicities.

Late-Onset Immunotherapy Toxicity and Delayed Autoantibody Changes: Checkpoint Inhibitor-Induced Raynaud's-Like Phenomenon

Immune checkpoint inhibitor (ICI)‐induced immune‐related adverse events (irAEs) may affect almost any organ system and occur at any point during therapy. Autoantibody analysis may provide insight into the mechanism, nature, and timing of these events. We report a case of ICI‐induced late‐onset Raynaud's‐like phenomenon in a patient receiving combination immunotherapy. A 53‐year‐old woman with advanced non‐small lung cancer received combination anti‐cytotoxic T‐lymphocyte antigen 4 and anti‐programmed death 1 ICI therapy. She developed early (hypophysitis at 4 months) and late (Raynaud's at >20 months) irAEs. Longitudinal assessment of 124 autoantibodies was correlated with toxicity. Although autoantibody levels were generally stable for the first 18 months of therapy, shortly before the development of Raynaud's, a marked increase in multiple autoantibodies was observed. This case highlights the potential for delayed autoimmune toxicities and provides potential biologic insights into the dynamic nature of these events.

Key Points

A patient treated with dual anti‐PD1 and anti‐CTLA4 therapy developed Raynaud's‐like signs and symptoms more than 18 months after starting therapy.

In this case, autoantibody changes became apparent shortly before onset of clinical toxicity.

This case highlights the potential for late‐onset immune‐related adverse events checkpoint inhibitors, requiring continuous clinical vigilance.

The optimal duration of checkpoint inhibitor therapy in patients with profound and prolonged responses remains unclear.

Autoantibody analysis may provide insight into the mechanism, nature, and timing of immune‐related adverse events. This case report describes a case of immune checkpoint inhibitor‐induced late‐onset Raynaud's‐like phenomenon in a patient receiving combination immunotherapy.

Toll-Like Receptor 9 Deficiency Breaks Tolerance to RNA-Associated Antigens and Up-Regulates Toll-Like Receptor 7 Protein in Sle1 Mice

Objective

Toll‐like receptors (TLRs) 7 and 9 are important innate signaling molecules with opposing roles in the development and progression of systemic lupus erythematosus (SLE). While multiple studies support the notion of a dependency on TLR‐7 for disease development, genetic ablation of TLR‐9 results in severe disease with glomerulonephritis (GN) by a largely unknown mechanism. This study was undertaken to examine the suppressive role of TLR‐9 in the development of severe lupus in a mouse model.

Methods

We crossed Sle1 lupus‐prone mice with TLR‐9–deficient mice to generate Sle1TLR‐9^−/− mice. Mice ages 4.5–6.5 months were evaluated for severe autoimmunity by assessing splenomegaly, GN, immune cell populations, autoantibody and total Ig profiles, kidney dendritic cell (DC) function, and TLR‐7 protein expression. Mice ages 8–10 weeks were used for functional B cell studies, Ig profiling, and determination of TLR‐7 expression.

Results

Sle1TLR‐9^−/− mice developed severe disease similar to TLR‐9–deficient MRL and Nba2 models. Sle1TLR‐9^−/− mouse B cells produced more class‐switched antibodies, and the autoantibody repertoire was skewed toward RNA‐containing antigens. GN in these mice was associated with DC infiltration, and purified Sle1TLR‐9^−/− mouse renal DCs were more efficient at TLR‐7–dependent antigen presentation and expressed higher levels of TLR‐7 protein. Importantly, this increase in TLR‐7 expression occurred prior to disease development, indicating a role in the initiation stages of tissue destruction.

Conclusion

The increase in TLR‐7–reactive immune complexes, and the concomitant enhanced expression of their receptor, promotes inflammation and disease in Sle1TLR9^−/− mice.

Serum Antibody Profiling Healthy Toddlers Reveals Low, Intermediate and High Response Cohorts

During their first two years of life, infant/toddlers produce their own antibodies in response to diverse environmental exposures, infections, and vaccinations. The specificity of these antibodies towards self-antigens, infections, and vaccine components and how they vary among individuals is poorly characterized. We developed a new antigen array comprising autoantigens, infectious agents, and vaccine antigens to assess the serum antibody specificities among a cohort of &gt;150 health toddlers. An analysis of their responses reveals a stratification of healthy toddlers into 3 groups, low, intermediate, and high responders. Among those 16% that are high responders, their serum IgGs were able to bind a diverse array of self-antigens as well as infectious agents. Longitudinal follow-up suggests this pattern is relatively stable over time. Comparing clinical data reveals a significant correlation with the high IgG responder cohort and a family history of asthma and maternal gestational diabetes. Targeted DNA sequencing in the high responder group revealed a strong genetic association signal at the HLA locus, with polymorphisms at this locus associated with high ANA and IgG titers to many antigens. Findings from this study may provide insights into the natural history of human autoantibody formation.

Amino acid signatures of HLA Class-I and II molecules are strongly associated with SLE susceptibility and autoantibody production in Eastern Asians

Human leukocyte antigen (HLA) is a key genetic factor conferring risk of systemic lupus erythematosus (SLE), but precise independent localization of HLA effects is extremely challenging. As a result, the contribution of specific HLA alleles and amino-acid residues to the overall risk of SLE and to risk of specific autoantibodies are far from completely understood. Here, we dissected (a) overall SLE association signals across HLA, (b) HLA-peptide interaction, and (c) residue-autoantibody association. Classical alleles, SNPs, and amino-acid residues of eight HLA genes were imputed across 4,915 SLE cases and 13,513 controls from Eastern Asia. We performed association followed by conditional analysis across HLA, assessing both overall SLE risk and risk of autoantibody production. DR15 alleles HLA-DRB1*15:01 (P = 1.4x10^-27, odds ratio (OR) = 1.57) and HLA-DQB1*06:02 (P = 7.4x10^-23, OR = 1.55) formed the most significant haplotype (OR = 2.33). Conditioned protein-residue signals were stronger than allele signals and mapped predominantly to HLA-DRB1 residue 13 (P = 2.2x10^-75) and its proxy position 11 (P = 1.1x10^-67), followed by HLA-DRB1-37 (P = 4.5x10^-24). After conditioning on HLA-DRB1, novel associations at HLA-A-70 (P = 1.4x10^-8), HLA-DPB1-35 (P = 9.0x10^-16), HLA-DQB1-37 (P = 2.7x10^-14), and HLA-B-9 (P = 6.5x10^-15) emerged. Together, these seven residues increased the proportion of explained heritability due to HLA to 2.6%. Risk residues for both overall disease and hallmark autoantibodies (i.e., nRNP: DRB1-11, P = 2.0x10^-14; DRB1-13, P = 2.9x10^-13; DRB1-30, P = 3.9x10^-14) localized to the peptide-binding groove of HLA-DRB1. Enrichment for specific amino-acid characteristics in the peptide-binding groove correlated with overall SLE risk and with autoantibody presence. Risk residues were in primarily negatively charged side-chains, in contrast with rheumatoid arthritis. We identified novel SLE signals in HLA Class I loci (HLA-A, HLA-B), and localized primary Class II signals to five residues in HLA-DRB1, HLA-DPB1, and HLA-DQB1. These findings provide insights about the mechanisms by which the risk residues interact with each other to produce autoantibodies and are involved in SLE pathophysiology.

The Human leukocyte antigen (HLA) region is a key genetic factor conferring risk of systemic lupus erythematosus (SLE). In spite of multiple SLE association signals identified in the HLA region, only amino-acid residues within HLA-DRB1 have been specifically described previously. In this study, we performed an imputation-based analysis on individuals with East Asian ancestry, and characterized SLE risk within the HLA region for all involved independent genes (HLA-DRB1, HLA-DPB1, HLA-DQB1, HLA-A, and HLA-B). Furthermore, we identified a characteristic SLE risk residue signature as well as a pattern of specific nRNP and Ro/La autoantibody residues located in the peptide-binding grooves, suggesting their key involvement in autoantibody production.

T cell-intrinsic IL-1R signaling licenses effector cytokine production by memory CD4 T cells

Innate cytokines are critical drivers of priming and differentiation of naive CD4 T cells, but their functions in memory T cell response are largely undefined. Here we show that IL-1 acts as a licensing signal to permit effector cytokine production by pre-committed Th1 (IFN-γ), Th2 (IL-13, IL-4, and IL-5) and Th17 (IL-17A, IL-17F, and IL-22) lineage cells. This licensing function of IL-1 is conserved across effector CD4 T cells generated by diverse immunological insults. IL-1R signaling stabilizes cytokine transcripts to enable productive and rapid effector functions. We also demonstrate that successful lineage commitment does not translate into productive effector functions in the absence of IL-1R signaling. Acute abrogation of IL-1R signaling in vivo results in reduced IL-17A production by intestinal Th17 cells. These results extend the role of innate cytokines beyond CD4 T cell priming and establish IL-1 as a licensing signal for memory CD4 T cell function.

CD4 T cell polarizations and functions are regulated by cytokines from innate cells. Here the authors show that IL-1 deficiency does not impair the differentiation of Th1, Th2 and Th17, but IL-1 signaling is required for maintaining the expressions of their respective key cytokines to ‘license’ the functions of these T cell subsets.

A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.

The impact of Staphylococcus aureus genomic variation on clinical phenotype of children with acute hematogenous osteomyelitis

Background

Children with acute hematogenous osteomyelitis (AHO) have a broad spectrum of illness ranging from mild to severe. The purpose of this study is to evaluate the impact of genomic variation of Staphylococcus aureus on clinical phenotype of affected children and determine which virulence genes correlate with severity of illness.

Methods

De novo whole genome sequencing was conducted for a strain of Community Acquired Methicillin Resistant Staphylococcus aureus (CA-MRSA), using PacBio Hierarchical Genome Assembly Process (HGAP) from 6 Single Molecule Real Time (SMRT) Cells, as a reference for DNA library assembly of 71 Staphylococcus aureus isolates from children with AHO. Virulence gene annotation was based on exhaustive literature review and genomic data in NCBI for Staphylococcus aureus. Clinical phenotype was assessed using a validated severity score. Kruskal-Wallis rank sum test determined association between clinical severity and virulence gene presence using False Discovery Rate (FDR), significance <0.01.

Results

PacBio produced an assembled genome of 2,898,306 bp and 2054 Open Reading Frames (ORFs). Annotation confirmed 201 virulence genes. Statistical analysis of gene presence by clinical severity found 40 genes significantly associated with severity of illness (FDR ≤0.009). MRSA isolates encoded a significantly greater number of virulence genes than did MSSA (p < 0.0001). Phylogenetic analysis by maximum likelihood (PAML) demonstrated the relatedness of genomic distance to clinical phenotype.

Conclusions

The Staphylococcus aureus genome contains virulence genes which are significantly associated with severity of illness in children with osteomyelitis. This study introduces a novel reference strain and detailed annotation of Staphylococcus aureus virulence genes. While this study does not address bacterial gene expression, a platform is created for future transcriptome investigations to elucidate the complex mechanisms involved in childhood osteomyelitis.

Association of Novel ALX4 Gene Polymorphisms with Antidepressant Treatment Response: Findings from the CO-MED Trial

Genome-wide association studies (GWAS) were conducted in participants of the CO-MED (Combining Medications to Enhance Depression Outcomes) trial, a randomized, 3-treatment arm clinical trial of major depressive disorder (MDD) designed to identify markers of differential treatment outcome (response and remission). The QIDS-SR (Quick Inventory of Depressive Symptomatology, Self-Reported version) was used to measure response at week 6 (QIDS-SR ≤5) and remission at week 12 (QIDS-SR ≤6 and ≤8 at the last two study visits). Three treatment groups (escitalopram monotherapy, escitalopram + bupropion, and venlafaxine + mirtazapine) were evaluated. GWAS identified a potentially regulatory SNP rs10769025 in the ALX4 gene on chromosome 11 with a strong association (p value = 9.85925E-08) with response to escitalopram monotherapy in Caucasians. Further, haplotype analysis on 7 ALX4 variants showed that a regulatory haplotype CAAACTG was significantly associated (odds ratio = 3.4, p = 2.00E-04) with response to escitalopram monotherapy at week 6. Ingenuity pathway analyses in the present study suggest that ALX4 has an indirect connection with antidepressant gene pathways in MDD, which may account for the genetic association with treatment outcome. Functional genomics studies to investigate the role of ALX4 in antidepressant treatment outcome will be an interesting future direction.

Profiling Serum Antibody Specificities in Infants Reveals a Significant Number with Autoreactive Antibodies

The antibody repertoire of an infant develops in response to infections, environmental exposures, and vaccinations. By adulthood, many will produce antibodies that react against self-antigens, often causal to autoimmune diseases such as lupus erythematosus. Limited scientific literature exists as to whether such autoreactive antibodies develop early in infancy. For this reason, the antibody specificities in the serum of healthy infants at 1 and 2 years of age was analyzed. Screening in a cohort of 82 infants revealed that 23 (28%) had moderate to high titered antibodies directed to diverse self-antigens. These numbers are consistent with observed ANA positivity in adults. Ongoing custom targeted DNA sequencing analysis will assess the genetic load of known autoimmune risk alleles in the ANA positive infants. Relationships between the antibody specificities and the genetic risk alleles assembled for each infant will be presented. The comparisons will include clinical information; sex, ethnicity, growth records, vaccination status, infectious history, antibiotic and antiviral treatments, disease status, and family history. Genomic analysis of infants with autoantibodies may facilitate implementation of a new wellness screen to identify antibody positive infants that have significant genetic predisposition and therefore are at risk for developing diverse immune system abnormalities including autoimmune disorders.

T cell intrinsic IL-1R signaling licenses effector cytokine production by memory CD4 T cells

Innate cytokines are critical drivers of priming and differentiation of naïve CD4 T cells but their role in memory T cell function is largely undefined. We discovered that while TCR and CD28 ligation are primary requirements for memory T cell reactivation, IL-1 acts as a licensing signal to permit effector cytokine production by pre-committed Th1 (IFNg), Th2 (IL-13, IL-4 and IL-5) and Th17 (IL-17A, IL-17F and IL-22) lineage cells. This licensing function of IL-1 is conserved across effector CD4 T cells generated by diverse pathogenic insults and different routes of infections. We found that Il1r−/− CD4 T cells from spleen as well as small intestine can differentiate in vivo into Th1, Th2 and Th17 lineages at the steady state. However, the same CD4 T cells are unable to secrete effector cytokines following TCR ligation. Our data demonstrates that successful lineage commitment does not translate into productive effector function in the absence of IL-1R signaling. Acute abrogation of IL-1R signaling in vivo also resulted in reduced in situ IL-17A production by primed intestinal Th17 cells. Consistently, reactivation of circulating memory CD4 T cells from healthy human donors with autologous DCs in the presence of hIL-1R antagonist (Anakinra) resulted in significantly diminished IL-17A production. Mechanistically, abrogation of IL-1R signaling during reactivation led to rapid decay of T cell cytokine transcripts. Conversely, the presence of recombinant IL-1b during T cell reactivation rescued destabilized cytokine transcripts in a p38 dependent fashion. Together, these results uncover a novel layer of regulation of memory CD4 T cells and establish a broadly conserved role for IL-1 in memory CD4 T cell effector function.

Transcriptional profiling of pathogen-specific CD4 T cells reveals key events that regulate Th17 cell priming and differentiation

CD4 T cells play a central role in adaptive immunity through secretion of specific cytokines which exert effector functions against invading pathogens. They also drive B cell activation and coordinate CD8 T cell responses. Productive CD4 T cell responses are determined not only by the specificity of the T cell receptor, but also by innate immune cues provided by dendritic cells (DCs) which shape the nature of the effector lineage commitment. However, to what extent pathogen-recognizing DCs prime naïve T cells on a global transcriptional level remains elusive, because of the inability to track and identify pathogen-specific T cells during cognate DC-T cell interaction. Here, we analyzed transcriptional profiles of pathogen-specific mouse Th17 cells utilizing an in vitro approach, where naïve CD4 T cells were differentiated by mature DCs in a highly antigen-specific manner. These pathogen-specific Th17 cells exhibit enormous lineage heterogeneity and resemble Th17 cells found in vivo following infections. We further discovered that DC-mediated priming not only generates a distinct transcriptional landscape in pathogen-specific Th17 cells, but also warrants their acquisition of memory gene signature and metabolic switch from glycolysis to oxidative phosphorylation. These results indicate that pathogen-activated dendritic cells play a critical role in CD4 T cell transcription programming that extends beyond determining lineage commitment. By mining this dataset, we have identified a previously unappreciated gene in promoting Th17 lineage commitment. Conceptually, our study underlines the necessity of investigating physiologically differentiated CD4 T cells in order to reveal novel regulators of T cell biology.

Identification of Btnl2 as a Key Gene in Suppressing Autoimmunity: New Players in SLE pathogenesis

Systemic lupus erythematosus (SLE) is a chronic multisystem autoimmune disease with extensive clinical heterogeneity, affecting mainly females. It is initiated with breach in immune tolerance and as a result one of the main characteristics of SLE is the production of autoantibodies to a wide spectrum of self-antigens. Sle1 is a potent susceptibility locus for spontaneous systemic autoimmunity derived from the NZM2410 mouse strain. The NZW-derived suppressive modifier locus, Sles1, specifically prevents the spontaneous loss in tolerance mediated by the B6.Sle1 congenic. Fine mapping combined with phenotypic characterization of series of recombinant mouse strains have mapped Sles1 to an approximately 638 Kb segment on murine chromosome 17. This revealed extensive epistatic gene interactions within the Sles1 region and suggested Btnl2 and H2 genes as the strongest candidates. We generated BAC transgenic mice by introducing BAC harboring B6 allele of Btnl2 alone or Btnl2 and H2 combined allelic region on to Sle1.Sles1 mice to identify the gene driving the suppressive phenotype. Phenotypic characterization of aging BAC transgenic mice is completed and results strongly indicate Btnl2 as the key gene within the 638 Kb region driving the suppression. Additionally we have performed targeted resequencing on B6 and B6.Sles1 mice and have identified functional polymorphisms within the Btnl2 gene between the two mice strains. We are currently generating and analyzing RNA-seq data on our BAC transgenic mice to gain insight in to the mechanisms on how Btnl2 is suppressing autoimmunity. This study confirms Btnl2 as the key driver gene in suppressing loss of tolerance in SLE pathogenesis.

Chromatin accessibility modulates the induction of innate immunity in cells of myeloid lineage

The Toll-like receptors (TLRs) of the innate immune system play a key role in the recognition of pathogens. TLRs signal through the recruitment of specific adaptor molecules, leading to activation of the transcription factors, such as AP-1, NF-κB and interferon regulatory factors (IRFs). Although TLRs are mainly expressed in immune cells, a systematic comparison of immune cells activation programs induced by specific individual TLRs is not yet available. Here, we used whole transcriptome sequencing to evaluate global changes in gene expression on monocytes, monocyte-derived macrophages (MDM) and monocyte-derived dendritic cells (MDDC) of five human donors before and after 18 hours stimulation through the TLR2, TLR4, TLR7 and TLR7/8 agonists (Pam3CSK4, LPS, R837 and R848). Cluster analysis demonstrated eight distinct clusters of unresponsive or responsive genes by the four TLR agonists in monocytes, MDM and MDDC. TLR7/8 stimulation by R848 led to the most robust cytokine and chemokine transcription in all cell types, especially in monocytes. Interestingly, we found that interferon-stimulated genes were strongly induced only by TLR4 activation in MDM and only by TLR7 and TLR8 activation in monocytes. Unexpectedly, there was no difference of activation of MAPK, NF-κB and IRF by LPS or R848 in these cell types. Surprisingly, we revealed that chromatin accessibility is significantly correlated with genes that are induced by specific TLR agonists on specific cell type(s), using ATAC-seq. Therefore, our results suggest in spite of robust early signal transduction events, factor that regulate chromatin accessibility play a critical role in regulation of inflammatory gene expression by TLR agonists in cells of the innate immune system.

Distinct patterns of innate immune activation by clinical isolates of respiratory syncytial virus

Respiratory syncytial virus (RSV) is a major respiratory pathogen of infants and young children. Multiple strains of both subgroup A and B viruses circulate during each seasonal epidemic. Genetic heterogeneity among RSV genomes, in large part due to the error prone RNA-dependent, RNA polymerase, could mediate variations in pathogenicity. We evaluated clinical strains of RSV for their ability to induce the innate immune response. Subgroup B viruses were used to infect human pulmonary epithelial cells (A549) and primary monocyte-derived human macrophages (MDM) from a variety of donors. Secretions of IL-6 and CCL5 (RANTES) from infected cells were measured following infection. Host and viral transcriptome expression were assessed using RNA-SEQ technology and the genomic sequences of several clinical isolates were determined. There were dramatic differences in the induction of IL-6 and CCL5 in both A549 cells and MDM infected with a variety of clinical isolates of RSV. Transcriptome analyses revealed that the pattern of innate immune activation in MDM was virus-specific and host-specific. Specifically, viruses that induced high levels of secreted IL-6 and CCL5 tended to induce cellular innate immune pathways whereas viruses that induced relatively low level of IL-6 or CCL5 did not induce or suppressed innate immune gene expression. Activation of the host innate immune response mapped to variations in the RSV G gene and the M2-1 gene. Viral transcriptome data indicated that there was a gradient of transcription across the RSV genome though in some strains, RSV G was the expressed in the highest amounts at late times post-infection. Clinical strains of RSV differ in cytokine/chemokine induction and in induction and suppression of host genes expression suggesting that these viruses may have inherent differences in virulence potential. Identification of the genetic elements responsible for these differences may lead to novel approaches to antiviral agents and vaccines.

RAD51 interconnects between DNA replication, DNA repair and immunity

RAD51, a multifunctional protein, plays a central role in DNA replication and homologous recombination repair, and is known to be involved in cancer development. We identified a novel role for RAD51 in innate immune response signaling. Defects in RAD51 lead to the accumulation of self-DNA in the cytoplasm, triggering a STING-mediated innate immune response after replication stress and DNA damage. In the absence of RAD51, the unprotected newly replicated genome is degraded by the exonuclease activity of MRE11, and the fragmented nascent DNA accumulates in the cytosol, initiating an innate immune response. Our data suggest that in addition to playing roles in homologous recombination-mediated DNA double-strand break repair and replication fork processing, RAD51 is also implicated in the suppression of innate immunity. Thus, our study reveals a previously uncharacterized role of RAD51 in initiating immune signaling, placing it at the hub of new interconnections between DNA replication, DNA repair, and immunity.

A missense variant in NCF1 is associated with susceptibility to multiple autoimmune diseases

Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease with a strong genetic component characterized by autoantibody production and a type I interferon signature. Here we report a missense variant (g.74779296G>A; p.Arg90His) in NCF1, encoding the p47^phox subunit of the phagocyte NADPH oxidase (NOX2), as the putative underlying causal variant that drives a strong SLE-associated signal detected by the Immunochip in the GTF2IRD1-GTF2I region at 7q11.23 with a complex genomic structure. We show that the p.Arg90His substitution, which is reported to cause reduced reactive oxygen species (ROS) production, predisposes to SLE (odds ratio (OR) = 3.47 in Asians (P_meta = 3.1 × 10^-104), OR = 2.61 in European Americans, OR = 2.02 in African Americans) and other autoimmune diseases, including primary Sjögren's syndrome (OR = 2.45 in Chinese, OR = 2.35 in European Americans) and rheumatoid arthritis (OR = 1.65 in Koreans). Additionally, decreased and increased copy numbers of NCF1 predispose to and protect against SLE, respectively. Our data highlight the pathogenic role of reduced NOX2-derived ROS levels in autoimmune diseases.

Identification of a Systemic Lupus Erythematosus Risk Locus Spanning ATG16L2, FCHSD2, and P2RY2 in Koreans

OBJECTIVE

Systemic lupus erythematosus (SLE) is a chronic autoimmune disorder whose etiology is incompletely understood, but likely involves environmental triggers in genetically susceptible individuals. Using an unbiased genome-wide association (GWA) scan and replication analysis, we sought to identify the genetic loci associated with SLE in a Korean population.

METHODS

A total of 1,174 SLE cases and 4,246 population controls from Korea were genotyped and analyzed with a GWA scan to identify single-nucleotide polymorphisms (SNPs) significantly associated with SLE, after strict quality control measures were applied. For select variants, replication of SLE risk loci was tested in an independent data set of 1,416 SLE cases and 1,145 population controls from Korea and China.

RESULTS

Eleven regions outside the HLA exceeded the genome-wide significance level (P = 5 × 10(-8) ). A novel SNP-SLE association was identified between FCHSD2 and P2RY2, peaking at rs11235667 (P = 1.03 × 10(-8) , odds ratio [OR] 0.59) on a 33-kb haplotype upstream of ATG16L2. In the independent replication data set, the SNP rs11235667 continued to show a significant association with SLE (replication meta-analysis P = 0.001, overall meta-analysis P = 6.67 × 10(-11) ; OR 0.63). Within the HLA region, the SNP-SLE association peaked in the class II region at rs116727542, with multiple independent effects observed in this region. Classic HLA allele imputation analysis identified HLA-DRB1*1501 and HLA-DQB1*0602, each highly correlated with one another, as most strongly associated with SLE. Ten previously established SLE risk loci were replicated: STAT1-STAT4, TNFSF4, TNFAIP3, IKZF1, HIP1, IRF5, BLK, WDFY4, ETS1, and IRAK1-MECP2. Of these loci, previously unreported, independent second risk effects of SNPs in TNFAIP3 and TNFSF4, as well as differences in the association with a putative causal variant in the WDFY4 region, were identified.

CONCLUSION

Further studies are needed to identify true SLE risk effects in other loci suggestive of a significant association, and to identify the causal variants in the regions of ATG16L2, FCHSD2, and P2RY2.

Mycobacterium tuberculosis produces small noncoding RNAs during infections to regulate mycobacterial and eukaryotic gene expression

Mycobacterium tuberculosis (Mtb) presently infects 1/3 of the World’s population, causing both active and latent pulmonary tuberculosis (TB). The pathogenic mechanisms used by Mtb remain poorly defined. We used a small RNA sequencing strategy on RNA extracted from Mtb-infected THP-1 macrophages to characterize the changes in eukaryotic microRNAs (miRs) and identify any mycobacterially-encoded small RNAs. A large number of mammalian microRNAs were differentially regulated during the infection course. Target prediction programs indicated that the mRNAs regulated by the miRs were involved in immune processes, metabolic pathways, cell communication, and developmental events. In addition, 35 Mtb-encoded small RNAs (18–30 nucleotides, mRs) were discovered. Several of the Mtb-encoded mRs were detected in lung tissue biopsies from TB infected monkeys. The 100–150 nucleotides surrounding the 35 Mtb-encoded small RNAs had extensive secondary RNA folding capabilities, including hairpins and antisense stretches more characteristic of eukaryotic pre-microRNAs. Two such small RNAs, Mtb-mR-1 and Mtb-mR-6 were transcribed in Mycobacterium avium complex and Mtb in a hairpin loop and antisense sequence-dependent manner. These mRs were not expressed in Mycobacterium smegmatis, revealing a restricted expression to pathogenic strains. Mtb-encoded mR-1 and mR-6 modulated both prokaryotic and eukaryotic gene expression. Taken together, our findings reveal potentially novel pathogenic processes utilized by Mtb during infections.

Regulatory polymorphisms modulate the expression of HLA class II molecules and promote autoimmunity

Targeted sequencing of sixteen SLE risk loci among 1349 Caucasian cases and controls produced a comprehensive dataset of the variations causing susceptibility to systemic lupus erythematosus (SLE). Two independent disease association signals in the HLA-D region identified two regulatory regions containing 3562 polymorphisms that modified thirty-seven transcription factor binding sites. These extensive functional variations are a new and potent facet of HLA polymorphism. Variations modifying the consensus binding motifs of IRF4 and CTCF in the XL9 regulatory complex modified the transcription of HLA-DRB1, HLA-DQA1 and HLA-DQB1 in a chromosome-specific manner, resulting in a 2.5-fold increase in the surface expression of HLA-DR and DQ molecules on dendritic cells with SLE risk genotypes, which increases to over 4-fold after stimulation. Similar analyses of fifteen other SLE risk loci identified 1206 functional variants tightly linked with disease-associated SNPs and demonstrated that common disease alleles contain multiple causal variants modulating multiple immune system genes.

Cutting Edge: Inhibiting TBK1 by Compound II Ameliorates Autoimmune Disease in Mice

TANK-binding kinase 1 (TBK1) is a serine/threonine protein kinase that plays a crucial role in innate immunity. Enhanced TBK1 function is associated with autoimmune diseases and cancer, implicating the potential benefit of therapeutically targeting TBK1. In this article, we examined a recently identified TBK1 inhibitor Compound II on treating autoimmune diseases. We found that Compound II is a potent and specific inhibitor of TBK1-mediated IFN response. Compound II inhibited polyinosinic-polycytidylic acid-induced immune activation in vitro and in vivo. Compound II treatment also ameliorated autoimmune disease phenotypes of Trex1(-/-) mice, increased mouse survival, and dampened the IFN gene signature in TREX1 mutant patient lymphoblasts. In addition, we found that TBK1 gene expression is elevated in systemic lupus erythematosus patient cells, and systemic lupus erythematosus cells with high IFN signature responded well to Compound II treatment. Together, our findings provided critical experimental evidence for inhibiting TBK1 with Compound II as an effective treatment for TREX1-associated autoimmune diseases and potentially other interferonopathies.

Genomic Heterogeneity of Methicillin Resistant Staphylococcus aureus Associated with Variation in Severity of Illness among Children with Acute Hematogenous Osteomyelitis

Introduction

The association between severity of illness of children with osteomyelitis caused by Methicillin-resistant Staphylococcus aureus (MRSA) and genomic variation of the causative organism has not been previously investigated. The purpose of this study is to assess genomic heterogeneity among MRSA isolates from children with osteomyelitis who have diverse severity of illness.

Materials and Methods

Children with osteomyelitis were prospectively studied between 2010 and 2011. Severity of illness of the affected children was determined from clinical and laboratory parameters. MRSA isolates were analyzed with next generation sequencing (NGS) and optical mapping. Sequence data was used for multi-locus sequence typing (MLST), phylogenetic analysis by maximum likelihood (PAML), and identification of virulence genes and single nucleotide polymorphisms (SNP) relative to reference strains.

Results

The twelve children studied demonstrated severity of illness scores ranging from 0 (mild) to 9 (severe). All isolates were USA300, ST 8, SCC mec IVa MRSA by MLST. The isolates differed from reference strains by 2 insertions (40 Kb each) and 2 deletions (10 and 25 Kb) but had no rearrangements or copy number variations. There was a higher occurrence of virulence genes among study isolates when compared to the reference strains (p = 0.0124). There were an average of 11 nonsynonymous SNPs per strain. PAML demonstrated heterogeneity of study isolates from each other and from the reference strains.

Discussion

Genomic heterogeneity exists among MRSA isolates causing osteomyelitis among children in a single community. These variations may play a role in the pathogenesis of variation in clinical severity among these children.

B cell-intrinsic CD84 and Ly108 maintain germinal center B cell tolerance

Signaling lymphocyte activation molecules (SLAMs) play an integral role in immune regulation. Polymorphisms in the SLAM family receptors are implicated in human and mouse model of lupus disease. The lupus-associated, somatically mutated, and class-switched pathogenic autoantibodies are generated in spontaneously developed germinal centers (GCs) in secondary lymphoid organs. The role and mechanism of B cell-intrinsic expression of polymorphic SLAM receptors that affect B cell tolerance at the GC checkpoint are not clear. In this study, we generated several bacterial artificial chromosome-transgenic mice that overexpress C57BL/6 (B6) alleles of different SLAM family genes on an autoimmune-prone B6.Sle1b background. B6.Sle1b mice overexpressing B6-derived Ly108 and CD84 exhibit a significant reduction in the spontaneously developed GC response and autoantibody production compared with B6.Sle1b mice. These data suggest a prominent role for Sle1b-derived Ly108 and CD84 in altering the GC checkpoint. We further confirm that expression of lupus-associated CD84 and Ly108 specifically on GC B cells in B6.Sle1b mice is sufficient to break B cell tolerance, leading to an increase in autoantibody production. In addition, we observe that B6.Sle1b B cells have reduced BCR signaling and a lower frequency of B cell-T cell conjugates; the reverse is seen in B6.Sle1b mice overexpressing B6 alleles of CD84 and Ly108. Finally, we find a significant decrease in apoptotic GC B cells in B6.Sle1b mice compared with B6 controls. Our study establishes a central role for GC B cell-specific CD84 and Ly108 expression in maintaining B cell tolerance in GCs and in preventing autoimmunity.

T/B-cell interactions are more transient in response to weak stimuli in SLE-prone mice

Changes in immune function during the course of systemic lupus erythematosus (SLE) are well characterized. Class-switched antinuclear antibodies are the hallmark of SLE, and T/B-cell interactions are thus critical. However, changes in immune function contributing to disease susceptibility are unknown. Here, we have analyzed primary T and B cells from a mouse model of SLE prior to the onset of disease. To allow cognate T-cell activation with low affinity, we have developed a lower potency peptide ligand for the OTII TCR. T- and B-cell couples formed less frequently and retained their polarity less efficiently preferentially in response to low-affinity stimulation in SLE-prone mice. This matched decreased recruitment of actin and Vav1 and an enhanced PKCΘ recruitment to the cellular interface in T cells. The induction of the GC B-cell marker GL7 was increased in T/B cell couples from SLE-prone mice when the T-cell numbers were limited. However, the overall gene expression changes were marginal. Taken together, the enhanced cell-couple transience may allow a more efficient sampling of a large number of T/B cell couples, preferentially in response to limiting stimuli, therefore enhancing the immune reactivity in the development of SLE.

Whole transcriptome RNA-seq analysis: tumorigenesis and metastasis of melanoma

Melanoma is the most malignant cutaneous cancer and causes over 9000 deaths annually. Because fatality rates from malignant melanoma (MM) increase dramatically upon metastasis, we investigated tumorigenesis and metastasis of MM in transcriptome analyses of three distinct cell lines that correspond with the stages of MM pathogenesis: the normal stage (HEMn-LP), the onset of MM (A375), and the metastasis stage (A2058). Using next-generation sequencing (NGS) technology, we detected asymmetrical expression of genes among the three cell lines, notably on chromosomes 9, 11, 12, and 14, suggesting their involvement in tumorigenesis and metastasis of MM. These genes were clustered into 41 categories based on their expression patterns, and their biological functions were analyzed using Ingenuity Pathway Analysis. In the top cancer-associated category, HIF1A, IL8, TERT, ONECUT1, and FOXA1 directly interacted with either transcription factors or cytokines that are known to be involved in the tumorigenesis or metastasis of other malignant tumors. The present data suggest that cytokine regulatory pathways in macrophages predominate over other pathways during the pathogenesis of MM. This study provides new targets for the downstream mechanistic studies of the tumorigenesis and metastasis of MM and demonstrates a new strategy for studies of the progression of other malignant cancers.

Two functional lupus-associated BLK promoter variants control cell-type- and developmental-stage-specific transcription

Efforts to identify lupus-associated causal variants in the FAM167A/BLK locus on 8p21 are hampered by highly associated noncausal variants. In this report, we used a trans-population mapping and sequencing strategy to identify a common variant (rs922483) in the proximal BLK promoter and a tri-allelic variant (rs1382568) in the upstream alternative BLK promoter as putative causal variants for association with systemic lupus erythematosus. The risk allele (T) at rs922483 reduced proximal promoter activity and modulated alternative promoter usage. Allelic differences at rs1382568 resulted in altered promoter activity in B progenitor cell lines. Thus, our results demonstrated that both lupus-associated functional variants contribute to the autoimmune disease association by modulating transcription of BLK in B cells and thus potentially altering immune responses.

Allelic heterogeneity in NCF2 associated with systemic lupus erythematosus (SLE) susceptibility across four ethnic populations

Recent reports have associated NCF2, encoding a core component of the multi-protein NADPH oxidase (NADPHO), with systemic lupus erythematosus (SLE) susceptibility in individuals of European ancestry. To identify ethnicity-specific and -robust variants within NCF2, we assessed 145 SNPs in and around the NCF2 gene in 5325 cases and 21 866 controls of European-American (EA), African-American (AA), Hispanic (HS) and Korean (KR) ancestry. Subsequent imputation, conditional, haplotype and bioinformatic analyses identified seven potentially functional SLE-predisposing variants. Association with non-synonymous rs17849502, previously reported in EA, was detected in EA, HS and AA (P(EA) = 1.01 × 10(-54), PHS = 3.68 × 10(-10), P(AA) = 0.03); synonymous rs17849501 was similarly significant. These SNPs were monomorphic in KR. Novel associations were detected with coding variants at rs35937854 in AA (PAA = 1.49 × 10(-9)), and rs13306575 in HS and KR (P(HS) = 7.04 × 10(-7), P(KR) = 3.30 × 10(-3)). In KR, a 3-SNP haplotype was significantly associated (P = 4.20 × 10(-7)), implying that SLE predisposing variants were tagged. Significant SNP-SNP interaction (P = 0.02) was detected between rs13306575 and rs17849502 in HS, and a dramatically increased risk (OR = 6.55) with a risk allele at each locus. Molecular modeling predicts that these non-synonymous mutations could disrupt NADPHO complex assembly. The risk allele of rs17849501, located in a conserved transcriptional regulatory region, increased reporter gene activity, suggesting in vivo enhancer function. Our results not only establish allelic heterogeneity within NCF2 associated with SLE, but also emphasize the utility of multi-ethnic cohorts to identify predisposing variants explaining additional phenotypic variance ('missing heritability') of complex diseases like SLE.

Transcriptome dynamics during human erythroid differentiation and development

To explore the mechanisms controlling erythroid differentiation and development, we analyzed the genome-wide transcription dynamics occurring during the differentiation of human embryonic stem cells (HESCs) into the erythroid lineage and development of embryonic to adult erythropoiesis using high throughput sequencing technology. HESCs and erythroid cells at three developmental stages: ESER (embryonic), FLER (fetal), and PBER (adult) were analyzed. Our findings revealed that the number of expressed genes decreased during differentiation, whereas the total expression intensity increased. At each of the three transitions (HESCs-ESERs, ESERs-FLERs, and FLERs-PBERs), many differentially expressed genes were observed, which were involved in maintaining pluripotency, early erythroid specification, rapid cell growth, and cell-cell adhesion and interaction. We also discovered dynamic networks and their central nodes in each transition. Our study provides a fundamental basis for further investigation of erythroid differentiation and development, and has implications in using ESERs for transfusion product in clinical settings.

SLE peripheral blood B cell, T cell and myeloid cell transcriptomes display unique profiles and each subset contributes to the interferon signature

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease that is characterized by defective immune tolerance combined with immune cell hyperactivity resulting in the production of pathogenic autoantibodies. Previous gene expression studies employing whole blood or peripheral blood mononuclear cells (PBMC) have demonstrated that a majority of patients with active disease have increased expression of type I interferon (IFN) inducible transcripts known as the IFN signature. The goal of the current study was to assess the gene expression profiles of isolated leukocyte subsets obtained from SLE patients. Subsets including CD19⁺ B lymphocytes, CD3⁺CD4⁺ T lymphocytes and CD33⁺ myeloid cells were simultaneously sorted from PBMC. The SLE transcriptomes were assessed for differentially expressed genes as compared to healthy controls. SLE CD33⁺ myeloid cells exhibited the greatest number of differentially expressed genes at 208 transcripts, SLE B cells expressed 174 transcripts and SLE CD3⁺CD4⁺ T cells expressed 92 transcripts. Only 4.4% (21) of the 474 total transcripts, many associated with the IFN signature, were shared by all three subsets. Transcriptional profiles translated into increased protein expression for CD38, CD63, CD107a and CD169. Moreover, these studies demonstrated that both SLE lymphoid and myeloid subsets expressed elevated transcripts for cytosolic RNA and DNA sensors and downstream effectors mediating IFN and cytokine production. Prolonged upregulation of nucleic acid sensing pathways could modulate immune effector functions and initiate or contribute to the systemic inflammation observed in SLE.

Complete genome analysis of three Acinetobacter baumannii clinical isolates in China for insight into the diversification of drug resistance elements

Background

The emergence and rapid spreading of multidrug-resistant Acinetobacter baumannii strains has become a major health threat worldwide. To better understand the genetic recombination related with the acquisition of drug-resistant elements during bacterial infection, we performed complete genome analysis on three newly isolated multidrug-resistant A. baumannii strains from Beijing using next-generation sequencing technology.

Methodologies/Principal Findings

Whole genome comparison revealed that all 3 strains share some common drug resistant elements including carbapenem-resistant bla_OXA-23 and tetracycline (tet) resistance islands, but the genome structures are diversified among strains. Various genomic islands intersperse on the genome with transposons and insertions, reflecting the recombination flexibility during the acquisition of the resistant elements. The blood-isolated BJAB07104 and ascites-isolated BJAB0868 exhibit high similarity on their genome structure with most of the global clone II strains, suggesting these two strains belong to the dominant outbreak strains prevalent worldwide. A large resistance island (RI) of about 121-kb, carrying a cluster of resistance-related genes, was inserted into the ATPase gene on BJAB07104 and BJAB0868 genomes. A 78-kb insertion element carrying tra-locus and bla_OXA-23 island, can be either inserted into one of the tniB gene in the 121-kb RI on the chromosome, or transformed to conjugative plasmid in the two BJAB strains. The third strains of this study, BJAB0715, which was isolated from spinal fluid, exhibit much more divergence compared with above two strains. It harbors multiple drug-resistance elements including a truncated AbaR-22-like RI on its genome. One of the unique features of this strain is that it carries both bla_OXA-23 and bla_OXA-58 genes on its genome. Besides, an Acinetobacter lwoffii adeABC efflux element was found inserted into the ATPase position in BJAB0715.

Conclusions

Our comparative analysis on currently completed Acinetobacter baumannii genomes revealed extensive and dynamic genome organizations, which may facilitate the bacteria to acquire drug-resistance elements into their genomes.