Search

Understanding factors that impact CAR T cell expansion in the clinic is crucial to improving its therapeutic success. Here the authors document heterogeneity in the clonal dynamics of CAR-T cells by tracking individual clones using the endogenous TCR and integration sites and provide further insights into the role of transcriptional states in clonal kinetics.

Abstract Previous advances have identified immune pathways associated with inflammatory skin diseases, leading to the development of targeted therapies. However, there is a lack of molecular approaches that delineate these pathways at the individual patient level for personalized diagnostic and therapeutic guidance. Here, we conduct a cross-comparison of expression profiles from multiple inflammatory skin diseases to identify gene modules defining relevant immune pathways. Seven modules are identified, representing key immune pathways: Th17, Th2, Th1, Type I IFNs, neutrophilic, macrophagic, and eosinophilic. These modules allow the development of a molecular map with high diagnostic efficacy for inflammatory skin diseases and clinico-pathologically undetermined cases. Aligning dominant modules with treatment targets offers a rational framework for treatment selection, improving response rates in both treatment-naïve patients and non-responders to targeted therapies. Overall, our approach offers precision medicine for inflammatory skin diseases, utilizing transcriptional modules to support diagnosis and guide personalized treatment selection.

Abstract As robust cellular responses are important for protection against dengue, this phase 2 study evaluated the kinetics and phenotype of T cell responses induced by TAK-003, a live-attenuated tetravalent dengue vaccine, in 4–16-year-old living in dengue-endemic countries (NCT02948829). Two hundred participants received TAK-003 on Days 1 and 90. Interferon-gamma (IFN-γ) enzyme-linked immunospot assay [ELISPOT] and intracellular cytokine staining were used to analyze T cell response and functionality, using peptide pools representing non-structural (NS) proteins NS3 and NS5 matching DENV-1, -2, -3, and -4 and DENV-2 NS1. One month after the second TAK-003 dose (Day 120), IFN-γ ELISPOT T cell response rates against any peptide pool were 97.1% (95% CI: 93.4% to 99.1%), and similar for baseline dengue seropositive (96.0%) and seronegative (98.6%) participants. IFN-γ ELISPOT T cell response rates at Day 120 were 79.8%, 90.2%, 77.3%, and 74.0%, against DENV-1, -2, -3, and -4, respectively, and remained elevated through 3 years post-vaccination. Multifunctional CD4 and CD8 T cell responses against DENV-2 NS peptides were observed, independent of baseline serostatus: CD8 T cells typically secreted IFN-γ and TNF-α whereas CD4 T cells secreted ≥ 2 of IFN-γ, IL-2 and TNF-α cytokines. NAb titers and seropositivity rates remained substantially elevated through 3 years post-vaccination. Overall, TAK-003 was well tolerated and elicited durable T cell responses against all four DENV serotypes irrespective of baseline serostatus in children and adolescents aged 4–16 years living in dengue-endemic countries. TAK-003-elicited CD4 and CD8 T cells were multifunctional and persisted up to 3 years post-vaccination.

IntroductionSystemic lupus erythematosus (SLE) is characterized by dysregulated humoral immunity, leading to the generation of autoreactive B cells that can differentiate both within and outside of lymph node (LN) follicles.MethodsHere, we employed spatial transcriptomics and multiplex imaging to investigate the follicular immune landscaping and the in situ transcriptomic profile in LNs from SLE individuals.ResultsOur spatial transcriptomic analysis revealed robust type I IFN and plasma cell signatures in SLE compared to reactive, control follicles. Cell deconvolution revealed that follicular T cell subsets are mainly affected by the type I IFN fingerprint of SLE follicles. Dysregulation of TFH differentiation was documented by i) the significant reduction of Bcl6hi TFH cells, ii) the reduced cell density of potential IL-4 producing TFH cell subsets associated with the impaired transcriptomic signature of follicular IL-4 signaling and iii) the loss of their correlation with GC-B cells. This profile was accompanied by a marked reduction of Bcl6hi B cells and an enrichment of extrafollicular CD19hiCD11chiTbethi, age-associated B cells (ABCs), known for their autoreactive potential. The increased prevalence of follicular IL-21hi cells further reveals a hyperactive microenvironment in SLE compared to control.DiscussionTaken together, our findings highlight the altered immunological landscape of SLE follicles, likely fueled by potent inflammatory signals such as sustained type I IFN and/or IL-21 signaling. Our work provides novel insights into the spatial molecular and cellular signatures of SLE follicular B and TFH cell dynamics, and points to druggable targets to restore immune tolerance and enhance vaccine responses in SLE patients.

Abstract Broadly neutralizing antibodies (bNAbs) are promising candidates for the treatment and prevention of HIV-1 infections. Despite their critical importance, automatic detection of HIV-1 bNAbs from immune repertoires is still lacking. Here, we develop a straightforward computational method for the Rapid Automatic Identification of bNAbs (RAIN) based on machine learning methods. In contrast to other approaches, which use one-hot encoding amino acid sequences or structural alignment for prediction, RAIN uses a combination of selected sequence-based features for the accurate prediction of HIV-1 bNAbs. We demonstrate the performance of our approach on non-biased, experimentally obtained and sequenced BCR repertoires from HIV-1 immune donors. RAIN processing leads to the successful identification of distinct HIV-1 bNAbs targeting the CD4-binding site of the envelope glycoprotein. In addition, we validate the identified bNAbs using an in vitro neutralization assay and we solve the structure of one of them in complex with the soluble native-like heterotrimeric envelope glycoprotein by single-particle cryo-electron microscopy (cryo-EM). Overall, we propose a method to facilitate and accelerate HIV-1 bNAbs discovery from non-selected immune repertoires.

Abstract Interleukin (IL)-26 is a TH17 cytokine with known antimicrobial and pro-inflammatory functions. However, the precise role of IL-26 in the context of pathogenic TH17 responses is unknown. Here we identify a population of blood TH17 intermediates that produce high levels of IL-26 and differentiate into IL-17A-producing TH17 cells upon TGF-β1 exposure. By combining single cell RNA sequencing, TCR sequencing and spatial transcriptomics we show that this process occurs in psoriatic skin. In fact, IL-26+ TH17 intermediates infiltrating psoriatic skin induce TGF-β1 expression in basal keratinocytes and thereby promote their own differentiation into IL-17A-producing cells. Thus, our study identifies IL-26-producing cells as an early differentiation stage of TH17 cells that infiltrates psoriatic skin and controls its own maturation into IL17A-producing TH17 cells, via epithelial crosstalk involving paracrine production of TGF-β1.

Approximately 50% of Merkel cell carcinoma (MCC) patients facing this highly aggressive skin cancer initially respond positively to PD-1-based immunotherapy. Nevertheless, the recurrence of MCC post-immunotherapy emphasizes the pressing need for more effective treatments. Recent research has highlighted Cyclin-dependent kinases 4 and 6 (CDK4/6) as pivotal cell cycle regulators gaining prominence in cancer studies. This study reveals that the CDK4/6 inhibitor, palbociclib can enhance PD-L1 gene transcription and surface expression in MCC cells by activating HIF2α. Inhibiting HIF2α with TC-S7009 effectively counteracts palbociclib-induced PD-L1 transcription and significantly intensifies cell death in MCC. Simultaneously, co-targeting CDK4/6 and HIF2α boosts ROS levels while suppressing SLC7A11, a key regulator of cellular redox balance, promoting ferroptosis- a form of immunogenic cell death linked to iron. Considering the rising importance of immunogenic cell death in immunotherapy, this strategy holds promise for improving future MCC treatments, markedly increasing immunogenic cell death across various MCC cell lines, thus advancing cancer immunotherapy.

Measurement(s) Transcriptomics • Hemagglutination Inhibition Assay • IgG IgM IgA Total Measurement • Virus-neutralizing Antibody • ELISA Technology Type(s) Microarray • RNA sequencing • Hemagglutination Inhibition Assay • ELISA • Microneutralization Assay • serum neutralization of viral infectivity assay Sample Characteristic - Organism Homo sapiens

Abstract We consider an increasingly popular study design where single-cell RNA-seq data are collected from multiple individuals and the question of interest is to find genes that are differentially expressed between two groups of individuals. Towards this end, we propose a statistical method named IDEAS (individual level differential expression analysis for scRNA-seq). For each gene, IDEAS summarizes its expression in each individual by a distribution and then assesses whether these individual-specific distributions are different between two groups of individuals. We apply IDEAS to assess gene expression differences of autism patients versus controls and COVID-19 patients with mild versus severe symptoms.

Advanced prostate malignancies are a leading cause of cancer-related deaths in men, in large part due to our incomplete understanding of cellular drivers of disease progression. We investigate prostate cancer cell dynamics at single-cell resolution from disease onset to the development of androgen independence in an in vivo murine model. We observe an expansion of a castration-resistant intermediate luminal cell type that correlates with treatment resistance and poor prognosis in human patients. Moreover, transformed epithelial cells and associated fibroblasts create a microenvironment conducive to pro-tumorigenic immune infiltration, which is partially androgen responsive. Androgen-independent prostate cancer leads to significant diversification of intermediate luminal cell populations characterized by a range of androgen signaling activity, which is inversely correlated with proliferation and mRNA translation. Accordingly, distinct epithelial populations are exquisitely sensitive to translation inhibition, which leads to epithelial cell death, loss of pro-tumorigenic signaling, and decreased tumor heterogeneity. Our findings reveal a complex tumor environment largely dominated by castration-resistant luminal cells and immunosuppressive infiltrates.

Background Merkel cell carcinoma (MCC) often responds to PD-1 pathway blockade, regardless of tumor-viral status (~80% of cases driven by the Merkel cell polyomavirus (MCPyV)). Prior studies have characterized tumor-specific T cell responses to MCPyV, which have typically been CD8, but little is known about the T cell response to UV-induced neoantigens.Methods A patient in her mid-50s with virus-negative (VN) MCC developed large liver metastases after a brief initial response to chemotherapy. She received anti-PD-L1 (avelumab) and had a partial response within 4 weeks. Whole exome sequencing (WES) was performed to determine potential neoantigen peptides. Characterization of peripheral blood neoantigen T cell responses was evaluated via interferon-gamma (IFNγ) ELISpot, flow cytometry and single-cell RNA sequencing. Tumor-resident T cells were characterized by multiplexed immunohistochemistry.Results WES identified 1027 tumor-specific somatic mutations, similar to the published average of 1121 for VN-MCCs. Peptide prediction with a binding cut-off of ≤100 nM resulted in 77 peptides that were synthesized for T cell assays. Although peptides were predicted based on class I HLAs, we identified circulating CD4 T cells targeting 5 of 77 neoantigens. In contrast, no neoantigen-specific CD8 T cell responses were detected. Neoantigen-specific CD4 T cells were undetectable in blood before anti-PD-L1 therapy but became readily detectible shortly after starting therapy. T cells produced robust IFNγ when stimulated by neoantigen (mutant) peptides but not by the normal (wild-type) peptides. Single cell RNAseq showed neoantigen-reactive T cells expressed the Th1-associated transcription factor (T-bet) and associated cytokines. These CD4 T cells did not significantly exhibit cytotoxicity or non-Th1 markers. Within the pretreatment tumor, resident CD4 T cells were also Th1-skewed and expressed T-bet.Conclusions We identified and characterized tumor-specific Th1-skewed CD4 T cells targeting multiple neoantigens in a patient who experienced a profound and durable partial response to anti-PD-L1 therapy. To our knowledge, this is the first report of neoantigen-specific T cell responses in MCC. Although CD4 and CD8 T cells recognizing viral tumor antigens are often detectible in virus-positive MCC, only CD4 T cells recognizing neoantigens were detected in this patient. These findings suggest that CD4 T cells can play an important role in the response to anti-PD-(L)1 therapy.

Optimal T follicular helper (Tfh) cells function is important to promote the development of germinal centers and maturation of high affinity antigen-specific B cells. We have found that the expression of CXCR3 defines distinct Tfh subsets: CXCR3+ Th1-like Tfh cells mainly producing single IFN-γ and dual IL-21/IFN-γ and CXCR3- Th2-like Tfh cells mainly producing single IL-4 and dual IL-21/IL-4 cytokines. CXCR3- Th2-like Tfhs are significantly reduced during ongoing HIV replication. While the percentage of Th2-like Tfh cells correlates with that of total and cycling HIV-specific B cells, the percentage of CXCR3+ Th1-like Tfhs correlates with HIV-specific B cells expressing T-bet and CXCR3. Of note, only IL-4 and IL-21 cytokines boosted efficient maturation of HIV-specific B cells while IFN-γ induced expression of T-bet and CXCR3 in B cells. Interestingly, total and HIV-specific CXCR3+ B cells showed lower rate of somatic hypermutation, as compared to CXCR3- B cells. Therefore, the imbalance in Th2/Th1-like Tfhs affects B cell responses in viremic HIV infection.

To control HIV infection there is a need for vaccines to induce broad, potent and long-term B and T cell immune responses. With the objective to accelerate and maintain the induction of substantial levels of HIV-1 Env-specific antibodies and, at the same time, to enhance balanced CD4 and CD8 T cell responses, we evaluated the effect of concurrent administration of MF59-adjuvanted Env protein together with DNA or NYVAC vectors at priming to establish if early administration of Env leads to early induction of antibody responses. The primary goal was to assess the immunogenicity endpoint at week 26. Secondary endpoints were (i) to determine the quality of responses with regard to RV144 correlates of protection and (ii) to explore a potential impact of two late boosts. In this study, five different prime/boost vaccination regimens were tested in rhesus macaques. Animals received priming immunizations with either NYVAC or DNA alone or in combination with Env protein, followed by NYVAC + protein or DNA + protein boosts. All regimens induced broad, polyfunctional and well-balanced CD4 and CD8 T cell responses, with DNA-primed regimens eliciting higher response rates and magnitudes than NYVAC-primed regimens. Very high plasma binding IgG titers including V1/V2 specific antibodies, modest antibody-dependent cellular cytotoxicity (ADCC) and moderate neutralization activity were observed. Of note, early administration of the MF59-adjuvanted Env protein in parallel with DNA priming leads to more rapid elicitation of humoral responses, without negatively affecting the cellular responses, while responses were rapidly boosted after repeated immunizations, indicating the induction of a robust memory response. In conclusion, our findings support the use of the Env protein component during priming in the context of an heterologous immunization regimen with a DNA and/or NYVAC vector as an optimized immunization protocol against HIV infection.

ABSTRACT Heavy exposure to Mycobacterium tuberculosis, the etiologic agent of tuberculosis (TB) and among the top infectious killers worldwide, results in infection that is cleared, contained, or progresses to disease. Some heavily exposed tuberculosis contacts show no evidence of infection using the tuberculin skin test (TST) and interferon gamma release assay (IGRA); yet the mechanisms underlying this “resister” (RSTR) phenotype are unclear. To identify transcriptional responses that distinguish RSTR monocytes, we performed transcriptome sequencing (RNA-seq) on monocytes isolated from heavily exposed household contacts in Uganda and gold miners in South Africa after ex vivo M. tuberculosis infection. Gene set enrichment analysis (GSEA) revealed several gene pathways that were consistently enriched in response to M. tuberculosis among RSTR subjects compared to controls with positive TST/IGRA testing (latent TB infection [LTBI]) across Uganda and South Africa. The most significantly enriched gene set in which expression was increased in RSTR relative to LTBI M. tuberculosis-infected monocytes was the tumor necrosis factor alpha (TNF-α) signaling pathway whose core enrichment (leading edge) substantially overlapped across RSTR populations. These leading-edge genes included candidate resistance genes (ABCA1 and DUSP2) with significantly increased expression among Uganda RSTRs (false-discovery rate [FDR],

Abstract Background New strategies are needed to reduce the incidence of malaria, and promising approaches include the development of vaccines and monoclonal antibodies (mAbs) that target the circumsporozoite protein (CSP). To select the best candidates and speed development, it is essential to standardize preclinical assays to measure the potency of such interventions in animal models. Methods Two assay configurations were studied using transgenic Plasmodium berghei expressing Plasmodium falciparum full-length circumsporozoite protein. The assays measured (1) reduction in parasite infection of the liver (liver burden) following an intravenous (i.v) administration of sporozoites and (2) protection from parasitaemia following mosquito bite challenge. Two human CSP mAbs, AB311 and AB317, were compared for their ability to inhibit infection. Multiple independent experiments were conducted to define assay variability and resultant impact on the ability to discriminate differences in mAb functional activity. Results Overall, the assays produced highly consistent results in that all individual experiments showed greater functional activity for AB317 compared to AB311 as calculated by the dose required for 50% inhibition (ID50) as well as the serum concentration required for 50% inhibition (IC50). The data were then used to model experimental designs with adequate statistical power to rigorously screen, compare, and rank order novel anti-CSP mAbs. Conclusion The results indicate that in vivo assays described here can provide reliable information for comparing the functional activity of mAbs. The results also provide guidance regarding selection of the appropriate experimental design, dose selection, and group sizes.

Broadly neutralizing antibodies (bNAbs) are promising agents to prevent HIV infection and achieve HIV remission without antiretroviral therapy (ART). As with ART, bNAb combinations are likely needed to cover HIV's extensive diversity. Not all bNAbs are identical in terms of their breadth, potency, and in vivo longevity (half-life). Given these differences, it is important to optimally select the composition, or dose ratio, of combination bNAb therapies for future clinical studies. We developed a model that synthesizes 1) pharmacokinetics, 2) potency against a wide HIV diversity, 3) interaction models for how drugs work together, and 4) correlates that translate in vitro potency to clinical protection. We found optimization requires drug-specific balances between potency, longevity, and interaction type. As an example, tradeoffs between longevity and potency are shown by comparing a combination therapy to a bi-specific antibody (a single protein merging both bNAbs) that takes the better potency but the worse longevity of the two components. Then, we illustrate a realistic dose ratio optimization of a triple combination of VRC07, 3BNC117, and 10-1074 bNAbs. We apply protection estimates derived from both a non-human primate (NHP) challenge study meta-analysis and the human antibody mediated prevention (AMP) trials. In both cases, we find a 2:1:1 dose emphasizing VRC07 is nearly optimal. Our approach can be immediately applied to optimize the next generation of combination antibody prevention and cure studies.

The Fbw7 ubiquitin ligase targets many proteins for proteasomal degradation, which include oncogenic transcription factors (TFs) (e.g., c-Myc, c-Jun, and Notch). Fbw7 is a tumor suppressor and tumors often contain mutations in FBXW7, the gene that encodes Fbw7. The complexity of its substrate network has obscured the mechanisms of Fbw7-associated tumorigenesis, yet this understanding is needed for developing therapies. We used an integrated approach employing RNA-Seq and high-resolution mapping (cleavage under target and release using nuclease) of histone modifications and TF occupancy (c-Jun and c-Myc) to examine the combinatorial effects of misregulated Fbw7 substrates in colorectal cancer (CRC) cells with engineered tumor-associated FBXW7 null or missense mutations. Both Fbw7 mutations caused widespread transcriptional changes associated with active chromatin and altered TF occupancy: some were common to both Fbw7 mutant cell lines, whereas others were mutation specific. We identified loci where both Jun and Myc were coregulated by Fbw7, suggesting that substrates may have synergistic effects. One coregulated gene was CIITA, the master regulator of MHC Class II gene expression. Fbw7 loss increased MHC Class II expression and Fbw7 mutations were correlated with increased CIITA expression in TCGA colorectal tumors and cell lines, which may have immunotherapeutic implications for Fbw7-associated cancers. Analogous studies in neural stem cells in which FBXW7 had been acutely deleted closely mirrored the results in CRC cells. Gene set enrichment analyses revealed Fbw7-associated pathways that were conserved across both cell types that may reflect fundamental Fbw7 functions. These analyses provide a framework for understanding normal and neoplastic context-specific Fbw7 functions.

Background: In a phase 3 trial in African infants and children, the RTS,S/AS01 vaccine (GSK) showed moderate efficacy against clinical malaria. We sought to further understand RTS,S/AS01-induced immune responses associated with vaccine protection. Methods: Applying the blood transcriptional module (BTM) framework, we characterized the transcriptomic response to RTS,S/AS01 vaccination in antigen-stimulated (and vehicle control) peripheral blood mononuclear cells sampled from a subset of trial participants at baseline and month 3 (1-month post-third dose). Using a matched case–control study design, we evaluated which of these ‘RTS,S/AS01 signature BTMs’ associated with malaria case status in RTS,S/AS01 vaccinees. Antigen-specific T-cell responses were analyzed by flow cytometry. We also performed a cross-study correlates analysis where we assessed the generalizability of our findings across three controlled human malaria infection studies of healthy, malaria-naive adult RTS,S/AS01 recipients. Results: RTS,S/AS01 vaccination was associated with downregulation of B-cell and monocyte-related BTMs and upregulation of T-cell-related BTMs, as well as higher month 3 (vs. baseline) circumsporozoite protein-specific CD4+ T-cell responses. There were few RTS,S/AS01-associated BTMs whose month 3 levels correlated with malaria risk. In contrast, baseline levels of BTMs associated with dendritic cells and with monocytes (among others) correlated with malaria risk. The baseline dendritic cell- and monocyte-related BTM correlations with malaria risk appeared to generalize to healthy, malaria-naive adults. Conclusions: A prevaccination transcriptomic signature associates with malaria in RTS,S/AS01-vaccinated African children, and elements of this signature may be broadly generalizable. The consistent presence of monocyte-related modules suggests that certain monocyte subsets may inhibit protective RTS,S/AS01-induced responses. Funding: Funding was obtained from the NIH-NIAID (R01AI095789), NIH-NIAID (U19AI128914), PATH Malaria Vaccine Initiative (MVI), and Ministerio de Economía y Competitividad (Instituto de Salud Carlos III, PI11/00423 and PI14/01422). The RNA-seq project has been funded in whole or in part with Federal funds from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services, under grant number U19AI110818 to the Broad Institute. This study was also supported by the Vaccine Statistical Support (Bill and Melinda Gates Foundation award INV-008576/OPP1154739 to R.G.). C.D. was the recipient of a Ramon y Cajal Contract from the Ministerio de Economía y Competitividad (RYC-2008-02631). G.M. was the recipient of a Sara Borrell–ISCIII fellowship (CD010/00156) and work was performed with the support of Department of Health, Catalan Government grant (SLT006/17/00109). This research is part of the ISGlobal’s Program on the Molecular Mechanisms of Malaria which is partially supported by the Fundación Ramón Areces and we acknowledge support from the Spanish Ministry of Science and Innovation through the ‘Centro de Excelencia Severo Ochoa 2019–2023’ Program (CEX2018-000806-S), and support from the Generalitat de Catalunya through the CERCA Program.

RTS,S/AS01 (GSK) is the world’s first malaria vaccine. However, despite initial efficacy of almost 70% over the first 6 months of follow-up, efficacy waned over time. A deeper understanding of the immune features that contribute to RTS,S/AS01-mediated protection could be beneficial for further vaccine development. In two recent controlled human malaria infection (CHMI) trials of the RTS,S/AS01 vaccine in malaria-naïve adults, MAL068 and MAL071, vaccine efficacy against patent parasitemia ranged from 44% to 87% across studies and arms (each study included a standard RTS,S/AS01 arm with three vaccine doses delivered in four-week-intervals, as well as an alternative arm with a modified version of this regimen). In each trial, RTS,S/AS01 immunogenicity was interrogated using a broad range of immunological assays, assessing cellular and humoral immune parameters as well as gene expression. Here, we used a predictive modeling framework to identify immune biomarkers measured at day-of-challenge that could predict sterile protection against malaria infection. Using cross-validation on MAL068 data (either the standard RTS,S/AS01 arm alone, or across both the standard RTS,S/AS01 arm and the alternative arm), top-performing univariate models identified variables related to Fc effector functions and titer of antibodies that bind to the central repeat region (NANP6) of CSP as the most predictive variables; all NANP6-related variables consistently associated with protection. In cross-study prediction analyses of MAL071 outcomes (the standard RTS,S/AS01 arm), top-performing univariate models again identified variables related to Fc effector functions of NANP6-targeting antibodies as highly predictive. We found little benefit–with this dataset–in terms of improved prediction accuracy in bivariate models vs. univariate models. These findings await validation in children living in malaria-endemic regions, and in vaccinees administered a fourth RTS,S/AS01 dose. Our findings support a “quality as well as quantity” hypothesis for RTS,S/AS01-elicited antibodies against NANP6, implying that malaria vaccine clinical trials should assess both titer and Fc effector functions of anti-NANP6 antibodies.

Summary: We recently reported that the risk of sexually acquired HIV-1 infection is increased significantly by variants in the gene encoding CD101, a protein thought to modify inflammatory responses. Using blood samples from individuals with and without these variants, we demonstrate that CD101 variants modify the prevalence of circulating inflammatory cell types and show that CD101 variants are associated with increased proinflammatory cytokine production by circulating T cells. One category of CD101 variants is associated with a reduced capacity of regulatory T cells to suppress T cell cytokine production, resulting in a reduction in the baseline level of immune quiescence. These data are supported by transcriptomics data revealing alterations in the intrinsic regulation of antiviral pathways and HIV resistance genes in individuals with CD101 variants. Our data support the hypothesis that CD101 contributes to homeostatic regulation of bystander inflammation, with CD101 variants altering heterosexual HIV-1 acquisition by facilitating increased prevalence and altered function of T cell subsets.

Background Clinical benefit from programmed cell death 1 receptor (PD-1) inhibitors relies on reinvigoration of endogenous antitumor immunity. Nonetheless, robust immunological markers, based on circulating immune cell subsets associated with therapeutic efficacy are yet to be validated.Methods We isolated peripheral blood mononuclear cell from three independent cohorts of melanoma and Merkel cell carcinoma patients treated with PD-1 inhibitor, at baseline and longitudinally after therapy. Using multiparameter flow cytometry and cell sorting, we isolated four subsets of CD8+ T cells, based on PD-1 and TIGIT expression profiles. We performed phenotypic characterization, T cell receptor sequencing, targeted transcriptomic analysis and antitumor reactivity assays to thoroughly characterize each of these subsets.Results We documented that the frequency of circulating PD-1+TIGIT+ (DPOS) CD8+ T-cells after 1 month of anti-PD-1 therapy was associated with clinical response and overall survival. This DPOS T-cell population was enriched in highly activated T-cells, tumor-specific and emerging T-cell clonotypes and T lymphocytes overexpressing CXCR5, a key marker of the CD8 cytotoxic follicular T cell population. Additionally, transcriptomic profiling defined a specific gene signature for this population as well as the overexpression of specific pathways associated with the therapeutic response.Conclusions Our results provide a convincing rationale for monitoring this PD-1+TIGIT+ circulating population as an early cellular-based marker of therapeutic response to anti-PD-1 therapy.

DNA vectors have been widely used as a priming of poxvirus vaccine in prime/boost regimens. Whether the number of DNA impacts qualitatively or quantitatively the immune response is not fully explored. With the aim to reinforce T-cell responses by optimizing the prime-boost regimen, the multicentric EV03/ANRS VAC20 phase I/II trial, randomized 147 HIV-negative volunteers to either 3xDNA plus 1xNYVAC (weeks 0, 4, 8 plus 24; n = 74) or to 2xDNA plus 2xNYVAC (weeks 0, 4 plus 20, 24; n = 73) groups. T-cell responses (IFN-γ ELISPOT) to at least one peptide pool were higher in the 3xDNA than the 2xDNA groups (91% and 80% of vaccinees) (P = 0.049). In the 3xDNA arm, 26 (37%) recipients developed a broader T-cell response (Env plus at least to one of the Gag, Pol, Nef pools) than in the 2xDNA (15; 22%) arms (primary endpoint; P = 0.047) with a higher magnitude against Env (at week 26) (P

The RV144 vaccine trial showed reduced risk of HIV-1 acquisition, but mechanisms underlying protection are poorly understood. Here, Fourati et al. assess the transcriptomic profile of blood collected from 223 vaccinees and 40 placebo recipients and identify IRF7 as a mediator of protection.

Summary: High-throughput single-cell RNA sequencing (scRNA-seq) has become a frequently used tool to assess immune cell heterogeneity. Recently, the combined measurement of RNA and protein expression was developed, commonly known as cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq). Acquisition of protein expression data along with transcriptome data resolves some of the limitations inherent to only assessing transcripts but also nearly doubles the sequencing read depth required per single cell. Furthermore, there is still a paucity of analysis tools to visualize combined transcript-protein datasets. Here, we describe a targeted transcriptomics approach that combines an analysis of over 400 genes with simultaneous measurement of over 40 proteins on 2 × 104 cells in a single experiment. This targeted approach requires only about one-tenth of the read depth compared to a whole-transcriptome approach while retaining high sensitivity for low abundance transcripts. To analyze these multi-omic datasets, we adapted one-dimensional soli expression by nonlinear stochastic embedding (One-SENSE) for intuitive visualization of protein-transcript relationships on a single-cell level. : Mair et al. describe a targeted transcriptomics approach combined with surface protein measurement to capture immune cell heterogeneity at a low sequencing depth. One-SENSE is used as a visualization tool to intuitively explore the relationship of protein and transcript expression on the single-cell level. Keywords: single-cell RNA sequencing, multi-omic, AbSeq, high-dimensional cytometry, human immunology, One-SENSE, targeted transcriptomics, barcoded antibody, Rhapsody

A central tenet of reproducible research is that scientific results are published along with the underlying data and software code necessary to reproduce and verify the findings. A host of tools and software have been released that facilitate such work-flows and scientific journals have increasingly demanded that code and primary data be made available with publications. There has been little practical advice on implementing reproducible research work-flows for large ’omics’ or systems biology data sets used by teams of analysts working in collaboration. In such instances it is important to ensure all analysts use the same version of a data set for their analyses. Yet, instantiating relational databases and standard operating procedures can be unwieldy, with high "startup" costs and poor adherence to procedures when they deviate substantially from an analyst’s usual work-flow. Ideally a reproducible research work-flow should fit naturally into an individual’s existing work-flow, with minimal disruption. Here, we provide an overview of how we have leveraged popular open source tools, including Bioconductor, Rmarkdown, git version control, R, and specifically R’s package system combined with a new tool DataPackageR, to implement a lightweight reproducible research work-flow for preprocessing large data sets, suitable for sharing among small-to-medium sized teams of computational scientists. Our primary contribution is the DataPackageR tool, which decouples time-consuming data processing from data analysis while leaving a traceable record of how raw data is processed into analysis-ready data sets. The software ensures packaged data objects are properly documented and performs checksum verification of these along with basic package version management, and importantly, leaves a record of data processing code in the form of package vignettes. Our group has implemented this work-flow to manage, analyze and report on pre-clinical immunological trial data from multi-center, multi-assay studies for the past three years.

Malaria continues to be one of mankind's most devastating diseases despite the many and varied efforts to combat it. Indispensable for malaria elimination and eventual eradication is the development of effective vaccines. Controlled human malaria infection (CHMI) is an invaluable tool for vaccine efficacy assessment and investigation of early immunological and molecular responses against Plasmodium falciparum infection. Here, we investigated gene expression changes following CHMI using RNA-Seq. Peripheral blood samples were collected in Bagamoyo, Tanzania, from ten adults who were injected intradermally (ID) with 2.5x104 aseptic, purified, cryopreserved P. falciparum sporozoites (Sanaria® PfSPZ Challenge). A total of 2,758 genes were identified as differentially expressed following CHMI. Transcriptional changes were most pronounced on day 5 after inoculation, during the clinically silent liver phase. A secondary analysis, grouping the volunteers according to their prepatent period duration, identified 265 genes whose expression levels were linked to time of blood stage parasitemia detection. Gene modules associated with these 265 genes were linked to regulation of transcription, cell cycle, phosphatidylinositol signaling and erythrocyte development. Our study showed that in malaria pre-exposed volunteers, parasite prepatent period in each individual is linked to magnitude and timing of early gene expression changes after ID CHMI.

The RV144 Phase III clinical trial with ALVAC-HIV prime and AIDSVAX B/E subtypes CRF01_AE (A244) and B (MN) gp120 boost vaccine regime in Thailand provided a foundation for the future development of improved vaccine strategies that may afford protection against the human immunodeficiency virus type 1 (HIV-1). Results from this trial showed that immune responses directed against specific regions V1V2 of the viral envelope (Env) glycoprotein gp120 of HIV-1, were inversely correlated to the risk of HIV-1 infection. Due to the low production of gp120 proteins in CHO cells (2-20 mg/L), cleavage sites in V1V2 loops (A244) and V3 loop (MN) causing heterogeneous antigen products, it was an urgent need to generate CHO cells harboring A244 gp120 with high production yields and an additional, homogenous and uncleaved subtype B gp120 protein to replace MN used in RV144 for the future clinical trials. Here we describe the generation of Chinese Hamster Ovary (CHO) cell lines stably expressing vaccine HIV-1 Env antigens for these purposes: one expressing an HIV-1 subtype CRF01_AE A244 Env gp120 protein (A244.AE) and one expressing an HIV-1 subtype B 6240 Env gp120 protein (6240.B) suitable for possible future manufacturing of Phase I clinical trial materials with cell culture expression levels of over 100 mg/L. The antigenic profiles of the molecules were elucidated by comprehensive approaches including analysis with a panel of well-characterized monoclonal antibodies recognizing critical epitopes using Biacore and ELISA, and glycosylation analysis by mass spectrometry, which confirmed previously identified glycosylation sites and revealed unknown sites of O-linked and N-linked glycosylations at non-consensus motifs. Overall, the vaccines given with MF59 adjuvant induced higher and more rapid antibody (Ab) responses as well as higher Ab avidity than groups given with aluminum hydroxide. Also, bivalent proteins (A244.AE and 6240.B) formulated with MF59 elicited distinct V2-specific Abs to the epitope previously shown to correlate with decreased risk of HIV-1 infection in the RV144 trial. All together, these results provide critical information allowing the consideration of these candidate gp120 proteins for future clinical evaluations in combination with a potent adjuvant.

CD4 T cells harboring HIV-1/SIV represent a formidable hurdle to eradicating infection, and yet their detailed phenotype remains unknown. Here we integrate two single-cell technologies, flow cytometry and highly multiplexed quantitative RT-PCR, to characterize SIV-infected CD4 T cells directly ex vivo. Within individual cells, we correlate the cellular phenotype, in terms of host protein and RNA expression, with stages of the viral life cycle defined by combinatorial expression of viral RNAs. Spliced RNA+ infected cells display multiple memory and activation phenotypes, indicating virus production by diverse CD4 T cell subsets. In most (but not all) cells, progressive infection accompanies post-transcriptional downregulation of CD4 protein, while surface MHC class I is largely retained. Interferon-stimulated genes were also commonly upregulated. Thus, we demonstrate that combined quantitation of transcriptional and post-transcriptional regulation at the single-cell level informs in vivo mechanisms of viral replication and immune evasion.

To evaluate the role of V3-specific IgG antibodies (Abs) in the RV144 clinical HIV vaccine trial, which reduced HIV-1 infection by 31.2%, the anti-V3 Ab response was assessed. Vaccinees' V3 Abs were highly cross-reactive with cyclic V3 peptides (cV3s) from diverse virus subtypes. Sieve analysis of CRF01_AE breakthrough viruses from 43 vaccine- and 66 placebo-recipients demonstrated an estimated vaccine efficacy of 85% against viruses with amino acids mismatching the vaccine at V3 site 317 (p = 0.004) and 52% against viruses matching the vaccine at V3 site 307 (p = 0.004). This analysis was supported by data showing that vaccinees' plasma Abs were less reactive with I307 when replaced with residues found more often in vaccinees' breakthrough viruses. Simultaneously, viruses with mutations at F317 were less infectious, possibly due to the contribution of F317 to optimal formation of the V3 hydrophobic core. These data suggest that RV144-induced V3-specific Abs imposed immune pressure on infecting viruses and inform efforts to design an HIV vaccine.

Although the HVTN 505 DNA/recombinant adenovirus type 5 vector HIV-1 vaccine trial showed no overall efficacy, analysis of breakthrough HIV-1 sequences in participants can help determine whether vaccine-induced immune responses impacted viruses that caused infection. We analyzed 480 HIV-1 genomes sampled from 27 vaccine and 20 placebo recipients and found that intra-host HIV-1 diversity was significantly lower in vaccine recipients (P ≤ 0.04, Q-values ≤ 0.09) in Gag, Pol, Vif and envelope glycoprotein gp120 (Env-gp120). Furthermore, Env-gp120 sequences from vaccine recipients were significantly more distant from the subtype B vaccine insert than sequences from placebo recipients (P = 0.01, Q-value = 0.12). These vaccine effects were associated with signatures mapping to CD4 binding site and CD4-induced monoclonal antibody footprints. These results suggest either (i) no vaccine efficacy to block acquisition of any viral genotype but vaccine-accelerated Env evolution post-acquisition; or (ii) vaccine efficacy against HIV-1s with Env sequences closest to the vaccine insert combined with increased acquisition due to other factors, potentially including the vaccine vector.

Improved antigenicity against HIV-1 envelope (Env) protein is needed to elicit vaccine-induced protective immunity in humans. Here we describe the first tests in non-human primates (NHPs) of Env gp140 protein fused to a humanized anti-LOX-1 recombinant antibody for delivering Env directly to LOX-1-bearing antigen presenting cells, especially dendritic cells (DC). LOX-1, or 1ectin-like oxidized low-density lipoprotein (LDL) receptor-1, is expressed on various antigen presenting cells and endothelial cells, and is involved in promoting humoral immune responses. The anti-LOX-1 Env gp140 fusion protein was tested for priming immune responses and boosting responses in animals primed with replication competent NYVAC-KC Env gp140 vaccinia virus. Anti-LOX-1 Env gp140 vaccination elicited robust cellular and humoral responses when used for either priming or boosting immunity. Co-administration with Poly ICLC, a TLR3 agonist, was superior to GLA, a TLR4 agonist. Both CD4+ and CD8+ Env-specific T cell responses were elicited by anti-LOX-1 Env gp140, but in particular the CD4+ T cells were multifunctional and directed to multiple epitopes. Serum IgG and IgA antibody responses induced by anti-LOX-1 Env gp140 against various gp140 domains were cross-reactive across HIV-1 clades; however, the sera neutralized only HIV-1 bearing sequences most similar to the clade C 96ZM651 Env gp140 carried by the anti-LOX-1 vehicle. These data, as well as the safety of this protein vaccine, justify further exploration of this DC-targeting vaccine approach for protective immunity against HIV-1.

The RV144 clinical trial showed the partial efficacy of a vaccine regimen with an estimated vaccine efficacy (VE) of 31% for protecting low-risk Thai volunteers against acquisition of HIV-1. The impact of vaccine-induced immune responses can be investigated through sieve analysis of HIV-1 breakthrough infections (infected vaccine and placebo recipients). A V1/V2-targeted comparison of the genomes of HIV-1 breakthrough viruses identified two V2 amino acid sites that differed between the vaccine and placebo groups. Here we extended the V1/V2 analysis to the entire HIV-1 genome using an array of methods based on individual sites, k-mers and genes/proteins. We identified 56 amino acid sites or "signatures" and 119 k-mers that differed between the vaccine and placebo groups. Of those, 19 sites and 38 k-mers were located in the regions comprising the RV144 vaccine (Env-gp120, Gag, and Pro). The nine signature sites in Env-gp120 were significantly enriched for known antibody-associated sites (p = 0.0021). In particular, site 317 in the third variable loop (V3) overlapped with a hotspot of antibody recognition, and sites 369 and 424 were linked to CD4 binding site neutralization. The identified signature sites significantly covaried with other sites across the genome (mean = 32.1) more than did non-signature sites (mean = 0.9) (p < 0.0001), suggesting functional and/or structural relevance of the signature sites. Since signature sites were not preferentially restricted to the vaccine immunogens and because most of the associations were insignificant following correction for multiple testing, we predict that few of the genetic differences are strongly linked to the RV144 vaccine-induced immune pressure. In addition to presenting results of the first complete-genome analysis of the breakthrough infections in the RV144 trial, this work describes a set of statistical methods and tools applicable to analysis of breakthrough infection genomes in general vaccine efficacy trials for diverse pathogens.

Tenofovir gel is being evaluated for vaginal and rectal pre-exposure prophylaxis against HIV transmission. Because this is a new prevention strategy, we broadly assessed its effects on the mucosa. In MTN-007, a phase-1, randomized, double-blinded rectal microbicide trial, we used systems genomics/proteomics to determine the effect of tenofovir 1% gel, nonoxynol-9 2% gel, placebo gel or no treatment on rectal biopsies (15 subjects/arm). We also treated primary vaginal epithelial cells from four healthy women with tenofovir in vitro. After seven days of administration, tenofovir 1% gel had broad-ranging effects on the rectal mucosa, which were more pronounced than, but different from, those of the detergent nonoxynol-9. Tenofovir suppressed anti-inflammatory mediators, increased T cell densities, caused mitochondrial dysfunction, altered regulatory pathways of cell differentiation and survival, and stimulated epithelial cell proliferation. The breadth of mucosal changes induced by tenofovir indicates that its safety over longer-term topical use should be carefully monitored. Clinical trial registration: NCT01232803.

Flow cytometry is used increasingly in clinical research for cancer, immunology and vaccines. Technological advances in cytometry instrumentation are increasing the size and dimensionality of data sets, posing a challenge for traditional data management and analysis. Automated analysis methods, despite a general consensus of their importance to the future of the field, have been slow to gain widespread adoption. Here we present OpenCyto, a new BioConductor infrastructure and data analysis framework designed to lower the barrier of entry to automated flow data analysis algorithms by addressing key areas that we believe have held back wider adoption of automated approaches. OpenCyto supports end-to-end data analysis that is robust and reproducible while generating results that are easy to interpret. We have improved the existing, widely used core BioConductor flow cytometry infrastructure by allowing analysis to scale in a memory efficient manner to the large flow data sets that arise in clinical trials, and integrating domain-specific knowledge as part of the pipeline through the hierarchical relationships among cell populations. Pipelines are defined through a text-based csv file, limiting the need to write data-specific code, and are data agnostic to simplify repetitive analysis for core facilities. We demonstrate how to analyze two large cytometry data sets: an intracellular cytokine staining (ICS) data set from a published HIV vaccine trial focused on detecting rare, antigen-specific T-cell populations, where we identify a new subset of CD8 T-cells with a vaccine-regimen specific response that could not be identified through manual analysis, and a CyTOF T-cell phenotyping data set where a large staining panel and many cell populations are a challenge for traditional analysis. The substantial improvements to the core BioConductor flow cytometry packages give OpenCyto the potential for wide adoption. It can rapidly leverage new developments in computational cytometry and facilitate reproducible analysis in a unified environment.

Advances in high-throughput, single cell gene expression are allowing interrogation of cell heterogeneity. However, there is concern that the cell cycle phase of a cell might bias characterizations of gene expression at the single-cell level. We assess the effect of cell cycle phase on gene expression in single cells by measuring 333 genes in 930 cells across three phases and three cell lines. We determine each cell's phase non-invasively without chemical arrest and use it as a covariate in tests of differential expression. We observe bi-modal gene expression, a previously-described phenomenon, wherein the expression of otherwise abundant genes is either strongly positive, or undetectable within individual cells. This bi-modality is likely both biologically and technically driven. Irrespective of its source, we show that it should be modeled to draw accurate inferences from single cell expression experiments. To this end, we propose a semi-continuous modeling framework based on the generalized linear model, and use it to characterize genes with consistent cell cycle effects across three cell lines. Our new computational framework improves the detection of previously characterized cell-cycle genes compared to approaches that do not account for the bi-modality of single-cell data. We use our semi-continuous modelling framework to estimate single cell gene co-expression networks. These networks suggest that in addition to having phase-dependent shifts in expression (when averaged over many cells), some, but not all, canonical cell cycle genes tend to be co-expressed in groups in single cells. We estimate the amount of single cell expression variability attributable to the cell cycle. We find that the cell cycle explains only 5%-17% of expression variability, suggesting that the cell cycle will not tend to be a large nuisance factor in analysis of the single cell transcriptome.

Male circumcision provides partial protection against multiple sexually transmitted infections (STIs), including HIV, but the mechanisms are not fully understood. To examine potential vulnerabilities in foreskin epithelial structure, we used Wilcoxon paired tests adjusted using the false discovery rate method to compare inner and outer foreskin samples from 20 healthy, sexually active Peruvian males who have sex with males or transgender females, ages 21-29, at elevated risk of HIV infection. No evidence of epithelial microtrauma was identified, as assessed by keratinocyte activation, fibronectin deposition, or parakeratosis. However, multiple suprabasal tight junction differences were identified: 1) inner foreskin stratum corneum was thinner than outer (p = 0.035); 2) claudin 1 had extended membrane-bound localization throughout inner epidermis stratum spinosum (p = 0.035); 3) membrane-bound claudin 4 was absent from inner foreskin stratum granulosum (p = 0.035); and 4) occludin had increased membrane deposition in inner foreskin stratum granulosum (p = 0.042) versus outer. Together, this suggests subclinical inflammation and paracellular transport modifications to the inner foreskin. A setting of inflammation was further supported by inner foreskin epithelial explant cultures secreting higher levels of GM-CSF (p = 0.029), IP-10 (p = 0.035) and RANTES (p = 0.022) than outer foreskin, and also containing an increased density of CCR5+ and CD4+ CCR5+ cells (p = 0.022). Inner foreskin dermis also secreted more RANTES than outer (p = 0.036), and had increased density of CCR5+ cells (p = 0.022). In conclusion, subclinical changes to the inner foreskin of sexually active males may support an inflammatory state, with availability of target cells for HIV infection and modifications to epidermal barriers, potentially explaining the benefits of circumcision for STI prevention.

Neutralizing and non-neutralizing antibodies to linear epitopes on HIV-1 envelope glycoproteins have potential to mediate antiviral effector functions that could be beneficial to vaccine-induced protection. Here, plasma IgG responses were assessed in three HIV-1 gp120 vaccine efficacy trials (RV144, Vax003, Vax004) and in HIV-1-infected individuals by using arrays of overlapping peptides spanning the entire consensus gp160 of all major genetic subtypes and circulating recombinant forms (CRFs) of the virus. In RV144, where 31.2% efficacy against HIV-1 infection was seen, dominant responses targeted the C1, V2, V3 and C5 regions of gp120. An analysis of RV144 case-control samples showed that IgG to V2 CRF01_AE significantly inversely correlated with infection risk (OR= 0.54, p=0.0042), as did the response to other V2 subtypes (OR=0.60-0.63, p=0.016-0.025). The response to V3 CRF01_AE also inversely correlated with infection risk but only in vaccine recipients who had lower levels of other antibodies, especially Env-specific plasma IgA (OR=0.49, p=0.007) and neutralizing antibodies (OR=0.5, p=0.008). Responses to C1 and C5 showed no significant correlation with infection risk. In Vax003 and Vax004, where no significant protection was seen, serum IgG responses targeted the same epitopes as in RV144 with the exception of an additional C1 reactivity in Vax003 and infrequent V2 reactivity in Vax004. In HIV-1 infected subjects, dominant responses targeted the V3 and C5 regions of gp120, as well as the immunodominant domain, heptad repeat 1 (HR-1) and membrane proximal external region (MPER) of gp41. These results highlight the presence of several dominant linear B cell epitopes on the HIV-1 envelope glycoproteins. They also generate the hypothesis that IgG to linear epitopes in the V2 and V3 regions of gp120 are part of a complex interplay of immune responses that contributed to protection in RV144.

The RV144 clinical trial of a prime/boost immunizing regimen using recombinant canary pox (ALVAC-HIV) and two gp120 proteins (AIDSVAX B and E) was previously shown to have a 31.2% efficacy rate. Plasma specimens from vaccine and placebo recipients were used in an extensive set of assays to identify correlates of HIV-1 infection risk. Of six primary variables that were studied, only one displayed a significant inverse correlation with risk of infection: the antibody (Ab) response to a fusion protein containing the V1 and V2 regions of gp120 (gp70-V1V2). This finding prompted a thorough examination of the results generated with the complete panel of 13 assays measuring various V2 Abs in the stored plasma used in the initial pilot studies and those used in the subsequent case-control study. The studies revealed that the ALVAC-HIV/AIDSVAX vaccine induced V2-specific Abs that cross-react with multiple HIV-1 subgroups and recognize both conformational and linear epitopes. The conformational epitope was present on gp70-V1V2, while the predominant linear V2 epitope mapped to residues 165-178, immediately N-terminal to the putative α4β7 binding motif in the mid-loop region of V2. Odds ratios (ORs) were calculated to compare the risk of infection with data from 12 V2 assays, and in 11 of these, the ORs were ≤1, reaching statistical significance for two of the variables: Ab responses to gp70-V1V2 and to overlapping V2 linear peptides. It remains to be determined whether anti-V2 Ab responses were directly responsible for the reduced infection rate in RV144 and whether anti-V2 Abs will prove to be important with other candidate HIV vaccines that show efficacy, however, the results support continued dissection of Ab responses to the V2 region which may illuminate mechanisms of protection from HIV-1 infection and may facilitate the development of an effective HIV-1 vaccine.

We describe a model-based method, PING, for predicting nucleosome positions in MNase-Seq and MNase- or sonicated-ChIP-Seq data. PING compares favorably to NPS and TemplateFilter in scalability, accuracy and robustness to low read density. To demonstrate that PING predictions from widely available sonicated data can have sufficient spatial resolution to be to be useful for biological inference, we use Illumina H3K4me1 ChIP-seq data to detect changes in nucleosome positioning around transcription factor binding sites due to tamoxifen stimulation, to discriminate functional and non-functional transcription factor binding sites more effectively than with enrichment profiles, and to confirm that the pioneer transcription factor Foxa2 associates with the accessible major groove of nucleosomal DNA.

ChIP-Seq has become the standard method for genome-wide profiling DNA association of transcription factors. To simplify analyzing and interpreting ChIP-Seq data, which typically involves using multiple applications, we describe an integrated, open source, R-based analysis pipeline. The pipeline addresses data input, peak detection, sequence and motif analysis, visualization, and data export, and can readily be extended via other R and Bioconductor packages. Using a standard multicore computer, it can be used with datasets consisting of tens of thousands of enriched regions. We demonstrate its effectiveness on published human ChIP-Seq datasets for FOXA1, ER, CTCF and STAT1, where it detected co-occurring motifs that were consistent with the literature but not detected by other methods. Our pipeline provides the first complete set of Bioconductor tools for sequence and motif analysis of ChIP-Seq and ChIP-chip data.