Christopher R. Bolen

Dynamic expression profiling of type I and type III interferon-stimulated hepatocytes reveals a stable hierarchy of gene expression

Despite activating similar signaling cascades, the type I and type III interferons (IFNs) differ in their ability to antagonize virus replication. However, it is not clear whether these cytokines induce unique antiviral states, particularly in the liver, where the clinically important hepatitis B and C viruses cause persistent infection. Here, clustering and promoter analyses of microarray‐based gene expression profiling were combined with mechanistic studies of signaling pathways to dynamically characterize the transcriptional responses induced by these cytokines in Huh7 hepatoma cells and primary human hepatocytes. Type I and III IFNs differed greatly in their level of interferon‐stimulated gene (ISG) induction with a clearly detectable hierarchy (IFN‐β > IFN‐α > IFN‐λ3 > IFN‐λ1 > IFN‐λ2). Notably, although the hierarchy identified varying numbers of differentially expressed genes when quantified using common statistical thresholds, further analysis of gene expression over multiple timepoints indicated that the individual IFNs do not in fact regulate unique sets of genes. The kinetic profiles of IFN‐induced gene expression were also qualitatively similar with the important exception of IFN‐α. While stimulation with either IFN‐β or IFN‐λs resulted in a similar long‐lasting ISG induction, IFN‐α signaling peaked early after stimulation then declined due to a negative feedback mechanism. The quantitative expression hierarchy and unique kinetics of IFN‐α reveal potential specific roles for individual IFNs in the immune response, and elucidate the mechanism behind previously observed differences in IFN antiviral activity. While current clinical trials are focused on IFN‐λ1 as a potential antiviral therapy, the finding that IFN‐λ3 invariably possesses the highest activity among type III IFNs suggests that this cytokine may have superior clinical activity. (Hepatology 2014;59:1262‐1272)

Non-equivalence of Wnt and R-spondin ligands during Lgr5 + intestinal stem-cell self-renewal

The canonical Wnt/β-catenin signalling pathway governs diverse developmental, homeostatic and pathological processes. Palmitoylated Wnt ligands engage cell-surface frizzled (FZD) receptors and LRP5 and LRP6 co-receptors, enabling β-catenin nuclear translocation and TCF/LEF-dependent gene transactivation. Mutations in Wnt downstream signalling components have revealed diverse functions thought to be carried out by Wnt ligands themselves. However, redundancy between the 19 mammalian Wnt proteins and 10 FZD receptors and Wnt hydrophobicity have made it difficult to attribute these functions directly to Wnt ligands. For example, individual mutations in Wnt ligands have not revealed homeostatic phenotypes in the intestinal epithelium—an archetypal canonical, Wnt pathway-dependent, rapidly self-renewing tissue, the regeneration of which is fueled by proliferative crypt Lgr5+ intestinal stem cells (ISCs). R-spondin ligands (RSPO1–RSPO4) engage distinct LGR4–LGR6, RNF43 and ZNRF3 receptor classes, markedly potentiate canonical Wnt/β-catenin signalling, and induce intestinal organoid growth in vitro and Lgr5+ ISCs in vivo. However, the interchangeability, functional cooperation and relative contributions of Wnt versus RSPO ligands to in vivo canonical Wnt signalling and ISC biology remain unknown. Here we identify the functional roles of Wnt and RSPO ligands in the intestinal crypt stem-cell niche. We show that the default fate of Lgr5+ ISCs is to differentiate, unless both RSPO and Wnt ligands are present. However, gain-of-function studies using RSPO ligands and a new non-lipidated Wnt analogue reveal that these ligands have qualitatively distinct, non-interchangeable roles in ISCs. Wnt proteins are unable to induce Lgr5+ ISC self-renewal, but instead confer a basal competency by maintaining RSPO receptor expression that enables RSPO ligands to actively drive and specify the extent of stem-cell expansion. This functionally non-equivalent yet cooperative interaction between Wnt and RSPO ligands establishes a molecular precedent for regulation of mammalian stem cells by distinct priming and self-renewal factors, with broad implications for precise control of tissue regeneration.

LAMBDA AND ALPHA INTERFERONS INHIBIT HEPATITIS B VIRUS REPLICATION THROUGH A COMMON MOLECULAR MECHANISM BUT WITH DIFFERENT IN VIVO ACTIVITIES

The type III interferons (IFN-lambda1, 2, and 3) induce an antiviral response similar to IFN-alpha/beta, but mediate their activity through a unique receptor. We found that like IFN-alpha/beta, IFN-lambda prevents the assembly of HBV capsids, demonstrating convergence of the two signaling pathways through a single antiviral mechanism. In contrast to IFN-lambda, the structurally related cytokine interleukin (IL)-22 only minimally reduced HBV replication. The transcriptional program activated by IL-22 displayed little similarity to that induced by IFN-lambda, but instead resembled the response elicited by IL-6. We also found that murine IFN-lambda2 had only weak antiviral activity against HBV in the liver of transgenic mice, and that human IFN-lambda2 activity in serum correlated with the sensitivity of the cytokine to proteases. These results demonstrate that the IFN-alpha/beta and IFN-lambda anti-HBV responses operate through a single molecular mechanism, and support the notion that IFN-lambda plays a local, rather than systemic, role in antiviral immunity.

Expression of specific inflammasome gene modules stratifies older individuals into two extreme clinical and immunological states

Low-grade, chronic inflammation has been associated with many diseases of aging, but the mechanisms responsible for producing this inflammation remain unclear. Inflammasomes can drive chronic inflammation in the context of an infectious disease or cellular stress, and they trigger the maturation of interleukin-1β (IL-1β). Here we find that the expression of specific inflammasome gene modules stratifies older individuals into two extremes: those with constitutive expression of IL-1β, nucleotide metabolism dysfunction, elevated oxidative stress, high rates of hypertension and arterial stiffness; and those without constitutive expression of IL-1β, who lack these characteristics. Adenine and N4-acetylcytidine, nucleotide-derived metabolites that are detectable in the blood of the former group, prime and activate the NLRC4 inflammasome, induce the production of IL-1β, activate platelets and neutrophils and elevate blood pressure in mice. In individuals over 85 years of age, the elevated expression of inflammasome gene modules was associated with all-cause mortality. Thus, targeting inflammasome components may ameliorate chronic inflammation and various other age-associated conditions.

Quantitative set analysis for gene expression: a method to quantify gene set differential expression including gene-gene correlations

Enrichment analysis of gene sets is a popular approach that provides a functional interpretation of genome-wide expression data. Existing tests are affected by inter-gene correlations, resulting in a high Type I error. The most widely used test, Gene Set Enrichment Analysis, relies on computationally intensive permutations of sample labels to generate a null distribution that preserves gene–gene correlations. A more recent approach, CAMERA, attempts to correct for these correlations by estimating a variance inflation factor directly from the data. Although these methods generate P-values for detecting gene set activity, they are unable to produce confidence intervals or allow for post hoc comparisons. We have developed a new computational framework for Quantitative Set Analysis of Gene Expression (QuSAGE). QuSAGE accounts for inter-gene correlations, improves the estimation of the variance inflation factor and, rather than evaluating the deviation from a null hypothesis with a P-value, it quantifies gene-set activity with a complete probability density function. From this probability density function, P-values and confidence intervals can be extracted and post hoc analysis can be carried out while maintaining statistical traceability. Compared with Gene Set Enrichment Analysis and CAMERA, QuSAGE exhibits better sensitivity and specificity on real data profiling the response to interferon therapy (in chronic Hepatitis C virus patients) and Influenza A virus infection. QuSAGE is available as an R package, which includes the core functions for the method as well as functions to plot and visualize the results.

Cell subset prediction for blood genomic studies

BackgroundGenome-wide transcriptional profiling of patient blood samples offers a powerful tool to investigate underlying disease mechanisms and personalized treatment decisions. Most studies are based on analysis of total peripheral blood mononuclear cells (PBMCs), a mixed population. In this case, accuracy is inherently limited since cell subset-specific differential expression of gene signatures will be diluted by RNA from other cells. While using specific PBMC subsets for transcriptional profiling would improve our ability to extract knowledge from these data, it is rarely obvious which cell subset(s) will be the most informative.ResultsWe have developed a computational method (Subset Prediction from Enrichment Correlation, SPEC) to predict the cellular source for a pre-defined list of genes (i.e. a gene signature) using only data from total PBMCs. SPEC does not rely on the occurrence of cell subset-specific genes in the signature, but rather takes advantage of correlations with subset-specific genes across a set of samples. Validation using multiple experimental datasets demonstrates that SPEC can accurately identify the source of a gene signature as myeloid or lymphoid, as well as differentiate between B cells, T cells, NK cells and monocytes. Using SPEC, we predict that myeloid cells are the source of the interferon-therapy response gene signature associated with HCV patients who are non-responsive to standard therapy.ConclusionsSPEC is a powerful technique for blood genomic studies. It can help identify specific cell subsets that are important for understanding disease and therapy response. SPEC is widely applicable since only gene expression profiles from total PBMCs are required, and thus it can easily be used to mine the massive amount of existing microarray or RNA-seq data.

Individual heritable differences result in unique cell lymphocyte receptor repertoires of naive and antigen-experienced cells

The adaptive immune systems capability to protect the body requires a highly diverse lymphocyte antigen receptor repertoire. However, the influence of individual genetic and epigenetic differences on these repertoires is not typically measured. By leveraging the unique characteristics of B, CD4+ T and CD8+ T-lymphocyte subsets from monozygotic twins, we quantify the impact of heritable factors on both the V(D)J recombination process and on thymic selection. We show that the resulting biases in both V(D)J usage and N/P addition lengths, which are found in naïve and antigen experienced cells, contribute to significant variation in the CDR3 region. Moreover, we show that the relative usage of V and J gene segments is chromosomally biased, with ∼1.5 times as many rearrangements originating from a single chromosome. These data refine our understanding of the heritable mechanisms affecting the repertoire, and show that biases are evident on a chromosome-wide level.

The blood transcriptional signature of chronic hepatitis C virus is consistent with an ongoing interferon-mediated antiviral response.

Blood transcriptional profiling is a powerful tool for understanding global changes after infection, and may be useful for prognosis and prediction of drug treatment responses. This study characterizes the effects of chronic hepatitis C virus (HCV) infection on gene expression by analyzing blood samples from 10 treatment-naïve HCV patients and 6 healthy volunteers. Differential expression analysis of microarray data from peripheral blood mononuclear cells (PBMCs) identified a 136-gene signature, including 66 genes elevated in infected individuals. Most of the upregulated genes were associated with interferon (IFN) activity (including members of the OAS and MX families, ISG15, and IRF7), suggesting an ongoing immune response. This HCV signature was also found to be consistently enriched in many other viral infection and vaccination datasets. These genes were validated using a second cohort composed of 5 HCV patients and 5 healthy volunteers, confirming the upregulation of the IFN signature. In summary, this is the first study to directly compare blood transcriptional profiles from HCV patients with healthy controls. The results show that chronic HCV infection has a pronounced effect on gene expression in PBMCs of infected individuals, and significantly elevates the expression of a subset of IFN-stimulated genes.

Systems immunology reveals markers of susceptibility to West Nile virus infection.

ABSTRACT West Nile virus (WNV) infection is usually asymptomatic but can cause severe neurological disease and death, particularly in older patients, and how individual variations in immunity contribute to disease severity is not yet defined. Animal studies identified a role for several immunity-related genes that determine the severity of infection. We have integrated systems-level transcriptional and functional data sets from stratified cohorts of subjects with a history of WNV infection to define whether these markers can distinguish susceptibility in a human population. Transcriptional profiles combined with immunophenotyping of primary cells identified a predictive signature of susceptibility that was detectable years after acute infection (67% accuracy), with the most prominent alteration being decreased IL1B induction following ex vivo infection of macrophages with WNV. Deconvolution analysis also determined a significant role for CXCL10 expression in myeloid dendritic cells. This systems analysis identified markers of pathogenic mechanisms and offers insights into potential therapeutic strategies.

Dysregulation of B Cell Repertoire Formation in Myasthenia Gravis Patients Revealed through Deep Sequencing

Myasthenia gravis (MG) is a prototypical B cell-mediated autoimmune disease affecting 20–50 people per 100,000. The majority of patients fall into two clinically distinguishable types based on whether they produce autoantibodies targeting the acetylcholine receptor (AChR-MG) or muscle specific kinase (MuSK-MG). The autoantibodies are pathogenic, but whether their generation is associated with broader defects in the B cell repertoire is unknown. To address this question, we performed deep sequencing of the BCR repertoire of AChR-MG, MuSK-MG, and healthy subjects to generate ∼518,000 unique VH and VL sequences from sorted naive and memory B cell populations. AChR-MG and MuSK-MG subjects displayed distinct gene segment usage biases in both VH and VL sequences within the naive and memory compartments. The memory compartment of AChR-MG was further characterized by reduced positive selection of somatic mutations in the VH CDR and altered VH CDR3 physicochemical properties. The VL repertoire of MuSK-MG was specifically characterized by reduced V-J segment distance in recombined sequences, suggesting diminished VL receptor editing during B cell development. Our results identify large-scale abnormalities in both the naive and memory B cell repertoires. Particular abnormalities were unique to either AChR-MG or MuSK-MG, indicating that the repertoires reflect the distinct properties of the subtypes. These repertoire abnormalities are consistent with previously observed defects in B cell tolerance checkpoints in MG, thereby offering additional insight regarding the impact of tolerance defects on peripheral autoimmune repertoires. These collective findings point toward a deformed B cell repertoire as a fundamental component of MG.