Introduction Over the last decade, the field of systems vaccinology has emerged, in which high throughput transcriptomics and other omics assays are used to probe changes of the innate and adaptive immune system in response to vaccination. The goal of this study was to benchmark key technical and analytical parameters of RNA sequencing (RNA-seq) in the context of a multi-site, double-blind randomized vaccine clinical trial. Methods We collected longitudinal peripheral blood mononuclear cell (PBMC) samples from 10 subjects before and after vaccination with a live attenuated Francisella tularensis vaccine and performed RNA-Seq at two different sites using aliquots from the same sample to generate two replicate datasets (5 time points for 50 samples each). We evaluated the impact of (i) filtering lowly-expressed genes, (ii) using external RNA controls, (iii) fold change and false discovery rate (FDR) filtering, (iv) read length, and (v) sequencing depth on differential expressed genes (DEGs) concordance between replicate datasets. Using synthetic mRNA spike-ins, we developed a method for empirically establishing minimal read-count thresholds for maintaining fold change accuracy on a per-experiment basis. We defined a reference PBMC transcriptome by pooling sequence data and established the impact of sequencing depth and gene filtering on transcriptome representation. Lastly, we modeled statistical power to detect DEGs for a range of sample sizes, effect sizes, and sequencing depths. Results and Discussion Our results showed that (i) filtering lowly-expressed genes is recommended to improve fold-change accuracy and inter-site agreement, if possible guided by mRNA spike-ins (ii) read length did not have a major impact on DEG detection, (iii) applying fold-change cutoffs for DEG detection reduced inter-set agreement and should be used with caution, if at all, (iv) reduction in sequencing depth had a minimal impact on statistical power but reduced the identifiable fraction of the PBMC transcriptome, (v) after sample size, effect size (i.e. the magnitude of fold change) was the most important driver of statistical power to detect DEG. The results from this study provide RNA sequencing benchmarks and guidelines for planning future similar vaccine studies.

West Nile virus (WNV) is a neurotropic mosquito-borne flavivirus of global importance. Neuroinvasive WNV infection results in encephalitis and can lead to prolonged neurological impairment or death. Type I interferon (IFN-I) is crucial for promoting antiviral defenses through the induction of antiviral effectors, which function to restrict viral replication and spread. However, our understanding of the antiviral response to WNV infection is mostly derived from analysis of bulk cell populations. It is becoming increasingly apparent that substantial heterogeneity in cellular processes exists among individual cells, even within a seemingly homogenous cell population. Here, we present WNV-inclusive single-cell RNA sequencing (scRNA-seq), an approach to examine the transcriptional variation and viral RNA burden across single cells. We observed that only a few cells within the bulk population displayed robust transcription of IFN-β mRNA, and this did not appear to depend on viral RNA abundance within the same cell. Furthermore, we observed considerable transcriptional heterogeneity in the IFN-I response, with genes displaying high unimodal and bimodal expression patterns. Broadly, IFN-stimulated genes negatively correlated with viral RNA abundance, corresponding with a precipitous decline in expression in cells with high viral RNA levels. Altogether, we demonstrated the feasibility and utility of WNV-inclusive scRNA-seq as a high-throughput technique for single-cell transcriptomics and WNV RNA detection. This approach can be implemented in other models to provide insights into the cellular features of protective immunity and identify novel therapeutic targets.IMPORTANCE West Nile virus (WNV) is a clinically relevant pathogen responsible for recurrent epidemics of neuroinvasive disease. Type I interferon is essential for promoting an antiviral response against WNV infection; however, it is unclear how heterogeneity in the antiviral response at the single-cell level impacts viral control. Specifically, conventional approaches lack the ability to distinguish differences across cells with varying viral abundance. The significance of our research is to demonstrate a new technique for studying WNV infection at the single-cell level. We discovered extensive variation in antiviral gene expression and viral abundance across cells. This protocol can be applied to primary cells or in vivo models to better understand the underlying cellular heterogeneity following WNV infection for the development of targeted therapeutic strategies.

Conventional immunization strategies will likely be insufficient for the development of a broadly neutralizing antibody (bnAb) vaccine for HIV or other difficult pathogens because of the immunological hurdles posed, including B cell immunodominance and germinal center (GC) quantity and quality. We found that two independent methods of slow delivery immunization of rhesus monkeys (RMs) resulted in more robust T follicular helper (TFH) cell responses and GC B cells with improved Env-binding, tracked by longitudinal fine needle aspirates. Improved GCs correlated with the development of >20-fold higher titers of autologous nAbs. Using a new RM genomic immunoglobulin locus reference, we identified differential IgV gene use between immunization modalities. Ab mapping demonstrated targeting of immunodominant non-neutralizing epitopes by conventional bolus-immunized animals, whereas slow delivery-immunized animals targeted a more diverse set of epitopes. Thus, alternative immunization strategies can enhance nAb development by altering GCs and modulating the immunodominance of non-neutralizing epitopes.

Breast cancer is the second leading cause of cancer-related deaths in women in the United States. Triple-negative breast cancer constitutes a subset of breast cancer that is associated with higher rates of relapse, decreased survival, and limited therapeutic options for patients afflicted with this type of breast cancer. Mammalian orthoreovirus (reovirus) selectively infects and kills transformed cells, and a serotype 3 reovirus is in clinical trials to assess its efficacy as an oncolytic agent against several cancers. It is unclear if reovirus serotypes differentially infect and kill triple-negative breast cancer cells and if reovirus-induced cytotoxicity of breast cancer cells can be enhanced by modulating the activity of host molecules and pathways. Here, we generated reassortant reoviruses by forward genetics with enhanced infective and cytotoxic properties in triple-negative breast cancer cells. From a high-throughput screen of small-molecule inhibitors, we identified topoisomerase inhibitors as a class of drugs that enhance reovirus infectivity and cytotoxicity of triple-negative breast cancer cells. Treatment of triple-negative breast cancer cells with topoisomerase inhibitors activates DNA damage response pathways, and reovirus infection induces robust production of type III, but not type I, interferon (IFN). Although type I and type III IFNs can activate STAT1 and STAT2, triple-negative breast cancer cellular proliferation is only negatively affected by type I IFN. Together, these data show that reassortant viruses with a novel genetic composition generated by forward genetics in combination with topoisomerase inhibitors more efficiently infect and kill triple-negative breast cancer cells. IMPORTANCE Patients afflicted by triple-negative breast cancer have decreased survival and limited therapeutic options. Reovirus infection results in cell death of a variety of cancers, but it is unknown if different reovirus types lead to triple-negative breast cancer cell death. In this study, we generated two novel reoviruses that more efficiently infect and kill triple-negative breast cancer cells. We show that infection in the presence of DNA-damaging agents enhances infection and triple-negative breast cancer cell killing by reovirus. These data suggest that a combination of a genetically engineered oncolytic reovirus and topoisomerase inhibitors may provide a potent therapeutic option for patients afflicted with triple-negative breast cancer.

Whole-genome sequencing and comparative genomic analysis of immune-related genes of Cercocebus atys and Macaca mulatta identify candidate genes, such as ICAM2 and TLR4, that may explain the AIDS resistance of C. atys. Sooty mangabeys are a primate model for non-pathogenic simian immunodeficiency virus (SIV) infection. Guido Silvestri and colleagues sequenced the genome of a captive sooty mangabey from the Yerkes National Primate Research Center in Atlanta, Georgia. They compared this with the genomic sequences of other primates and suggest several candidate genes that may influence differences in SIV pathogenicity across nonhuman primates. The team also analysed genetic variation in the genomes of 10 sooty mangabeys from the same centre and the transcriptomes of 14 tissues. This work provides a resource for comparative genomics and studies on SIV pathogenicity. In contrast to infections with human immunodeficiency virus (HIV) in humans and simian immunodeficiency virus (SIV) in macaques, SIV infection of a natural host, sooty mangabeys (Cercocebus atys), is non-pathogenic despite high viraemia1. Here we sequenced and assembled the genome of a captive sooty mangabey. We conducted genome-wide comparative analyses of transcript assemblies from C. atys and AIDS-susceptible species, such as humans and macaques, to identify candidates for host genetic factors that influence susceptibility. We identified several immune-related genes in the genome of C. atys that show substantial sequence divergence from macaques or humans. One of these sequence divergences, a C-terminal frameshift in the toll-like receptor-4 (TLR4) gene of C. atys, is associated with a blunted in vitro response to TLR-4 ligands. In addition, we found a major structural change in exons 3–4 of the immune-regulatory protein intercellular adhesion molecule 2 (ICAM-2); expression of this variant leads to reduced cell surface expression of ICAM-2. These data provide a resource for comparative genomic studies of HIV and/or SIV pathogenesis and may help to elucidate the mechanisms by which SIV-infected sooty mangabeys avoid AIDS.

PROBLEM:In women, the use of progestin-based contraception may increase the risk of vaginal HIV acquisition. We previously showed in macaques that there is a significantly higher simian-human immunodeficiency virus (SHIV) acquisition rate in the luteal phase of the menstrual cycle, which presents a naturally high-progesterone state, and this may be attributable to altered expression of innate immune factors. We hypothesized that progestin-based contraception, especially depot medroxyprogesterone acetate (DMPA), would, in a similar way, affect mucosal immune factors that influence HIV acquisition risk.METHOD OF STUDY:We used a pig-tailed macaque model to evaluate the effects of two progestin-based contraceptives, DMPA, and levonorgestrel (LNG)/ethinyl estradiol (EE)-based combined oral contraceptives (COCs), on innate mucosal factors. We compared the vaginal epithelial thickness data from previous studies and used cytokine profiling and microarray analysis to evaluate contraception-induced molecular changes in the vagina.RESULTS:The administration of DMPA caused a reduction in the thickness of the vaginal epithelium relative to that of the follicular or luteal phase. DMPA also induced a significant increase in vaginal levels of the anti-inflammatory cytokine IL-10. Both DMPA- and LNG-based contraception induced a signature of gene expression similar to that of the luteal phase, only more exacerbated, including widespread downregulation of antiviral genes.CONCLUSION:The use of progestin-based contraception might engender a milieu that poses an increased risk of HIV acquisition as compared to both the luteal and follicular phases of the menstrual cycle.

B cells play a critical role in the immune response by producing antibodies, which display remarkable diversity. Here we describe a bioinformatic pipeline, BALDR (BCR Assignment of Lineage using De novo Reconstruction) that accurately reconstructs the paired heavy and light chain immunoglobulin gene sequences from Illumina single-cell RNA-seq data. BALDR was accurate for clonotype identification in human and rhesus macaque influenza vaccine and simian immunodeficiency virus vaccine induced vaccine-induced plasmablasts and naïve and antigen-specific memory B cells. BALDR enables matching of clonotype identity with single-cell transcriptional information in B cell lineages and will have broad application in the fields of vaccines, human immunodeficiency virus broadly neutralizing antibody development, and cancer.BALDR is available at https://github.com/BosingerLab/BALDR .

Nat. Genet. 49, 238–248 (2017); published online 9 January 2017; corrected after print 20 March 2017 In the version of this article initially published, the legend to Figure 4c stated that only one proband without SMCHD1 mutation was tested for D4Z4 methylation pattern. However, three probands and one affected family member without SMCHD1 mutation were tested, as shown in the figure.

The dynamics and molecular mechanisms underlying vaccine immunity in early childhood remain poorly understood. Here we applied systems approaches to investigate the innate and adaptive responses to trivalent inactivated influenza vaccine (TIV) and MF59-adjuvanted TIV (ATIV) in 90 14- to 24-mo-old healthy children. MF59 enhanced the magnitude and kinetics of serum antibody titers following vaccination, and induced a greater frequency of vaccine specific, multicytokine-producing CD4(+) T cells. Compared with transcriptional responses to TIV vaccination previously reported in adults, responses to TIV in infants were markedly attenuated, limited to genes regulating antiviral and antigen presentation pathways, and observed only in a subset of vaccinees. In contrast, transcriptional responses to ATIV boost were more homogenous and robust. Interestingly, a day 1 gene signature characteristic of the innate response (antiviral IFN genes, dendritic cell, and monocyte responses) correlated with hemagglutination at day 28. These findings demonstrate that MF59 enhances the magnitude, kinetics, and consistency of the innate and adaptive response to vaccination with the seasonal influenza vaccine during early childhood, and identify potential molecular correlates of antibody responses.

Respiratory syncytial viruses encode a nonstructural protein (NS1) that interferes with type I and III interferon and other antiviral responses. Proteomic studies were conducted on human A549 type II alveolar epithelial cells and type I interferon-deficient Vero cells (African green monkey kidney cells) infected with wild-type and NS1-deficient clones of human respiratory syncytial virus to identify other potential pathway and molecular targets of NS1 interference. These analyses included two-dimensional differential gel electrophoresis and quantitative Western blotting. Surprisingly, NS1 was found to suppress the induction of manganese superoxide dismutase (SOD2) expression in A549 cells and to a much lesser degree Vero cells in response to infection. Because SOD2 is not directly inducible by type I interferons, it served as a marker to probe the impact of NS1 on signaling of other cytokines known to induce SOD2 expression and/or indirect effects of type I interferon signaling. Deductive analysis of results obtained from cell infection and cytokine stimulation studies indicated that interferon-γ signaling was a potential target of NS1, possibly as a result of modulation of STAT1 levels. However, this was not sufficient to explain the magnitude of the impact of NS1 on SOD2 induction in A549 cells. Vero cell infection experiments indicated that NS1 targeted a component of the type I interferon response that does not directly induce SOD2 expression but is required to induce another initiator of SOD2 expression. STAT2 was ruled out as a target of NS1 interference using quantitative Western blot analysis of infected A549 cells, but data were obtained to indicate that STAT1 was one of a number of potential targets of NS1. A label-free mass spectrometry-based quantitative approach is proposed as a means of more definitive identification of NS1 targets.

The folate pathway plays a crucial role in the regeneration and repair of the adult CNS after injury. Here, we have shown in rodents that such repair occurs at least in part through DNA methylation. In animals with combined spinal cord and sciatic nerve injury, folate-mediated CNS axon regeneration was found to depend on injury-related induction of the high-affinity folate receptor 1 (Folr1). The activity of folate was dependent on its activation by the enzyme dihydrofolate reductase (Dhfr) and a functional methylation cycle. The effect of folate on the regeneration of afferent spinal neurons was biphasic and dose dependent and correlated closely over its dose range with global and gene-specific DNA methylation and with expression of both the folate receptor Folr1 and the de novo DNA methyltransferases. These data implicate an epigenetic mechanism in CNS repair. Folic acid and possibly other nontoxic dietary methyl donors may therefore be useful in clinical interventions to promote brain and spinal cord healing. If indeed the benefit of folate is mediated by epigenetic mechanisms that promote endogenous axonal regeneration, this provides possible avenues for new pharmacologic approaches to treating CNS injuries.