Distinct baseline immune characteristics associated with responses to conjugated and unconjugated pneumococcal polysaccharide vaccines in older adults.

Pneumococcal infections cause serious illness and death among older adults. The capsular polysaccharide vaccine PPSV23 and conjugated alternative PCV13 can prevent these infections; yet, underlying immunological responses and baseline predictors remain unknown. We vaccinated 39 older adults (>60 years) with PPSV23 or PCV13 and observed comparable antibody responses (day 28) and plasmablast transcriptional responses (day 10); however, the baseline predictors were distinct. Analyses of baseline flow cytometry and bulk and single-cell RNA-sequencing data revealed a baseline phenotype specifically associated with weaker PCV13 responses, which was characterized by increased expression of cytotoxicity-associated genes, increased frequencies of CD16+ natural killer cells and interleukin-17-producing helper T cells and a decreased frequency of type 1 helper T cells. Men displayed this phenotype more robustly and mounted weaker PCV13 responses than women. Baseline expression levels of a distinct gene set predicted PPSV23 responses. This pneumococcal precision vaccinology study in older adults uncovered distinct baseline predictors that might transform vaccination strategies and initiate novel interventions.

Sestrins induce natural killer function in senescent-like CD8+ T cells.

Aging is associated with remodeling of the immune system to enable the maintenance of life-long immunity. In the CD8+ T cell compartment, aging results in the expansion of highly differentiated cells that exhibit characteristics of cellular senescence. Here we found that CD27-CD28-CD8+ T cells lost the signaling activity of the T cell antigen receptor (TCR) and expressed a protein complex containing the agonistic natural killer (NK) receptor NKG2D and the NK adaptor molecule DAP12, which promoted cytotoxicity against cells that expressed NKG2D ligands. Immunoprecipitation and imaging cytometry indicated that the NKG2D-DAP12 complex was associated with sestrin 2. The genetic inhibition of sestrin 2 resulted in decreased expression of NKG2D and DAP12 and restored TCR signaling in senescent-like CD27-CD28-CD8+ T cells. Therefore, during aging, sestrins induce the reprogramming of non-proliferative senescent-like CD27-CD28-CD8+ T cells to acquire a broad-spectrum, innate-like killing activity.

H3K4me3 breadth is linked to cell identity and transcriptional consistency.

Trimethylation of histone H3 at lysine 4 (H3K4me3) is a chromatin modification known to mark the transcription start sites of active genes. Here, we show that H3K4me3 domains that spread more broadly over genes in a given cell type preferentially mark genes that are essential for the identity and function of that cell type. Using the broadest H3K4me3 domains as a discovery tool in neural progenitor cells, we identify novel regulators of these cells. Machine learning models reveal that the broadest H3K4me3 domains represent a distinct entity, characterized by increased marks of elongation. The broadest H3K4me3 domains also have more paused polymerase at their promoters, suggesting a unique transcriptional output. Indeed, genes marked by the broadest H3K4me3 domains exhibit enhanced transcriptional consistency and [corrected] increased transcriptional levels, and perturbation of H3K4me3 breadth leads to changes in transcriptional consistency. Thus, H3K4me3 breadth contains information that could ensure transcriptional precision at key cell identity/function genes.

Baseline immune states (BIS) associated with vaccine responsiveness and factors that shape the BIS.

Vaccines are among the greatest inventions in medicine, leading to the elimination or control of numerous diseases, including smallpox, polio, measles, rubella, and, most recently, COVID-19. Yet, the effectiveness of vaccines varies among individuals. In fact, while some recipients mount a robust response to vaccination that protects them from the disease, others fail to respond. Multiple clinical and epidemiological factors contribute to this heterogeneity in responsiveness. Systems immunology studies fueled by advances in single-cell biology have been instrumental in uncovering pre-vaccination immune cell types and genomic features (i.e., the baseline immune state, BIS) that have been associated with vaccine responsiveness. Here, we review clinical factors that shape the BIS, and the characteristics of the BIS associated with responsiveness to frequently studied vaccines (i.e., influenza, COVID-19, bacterial pneumonia, malaria). Finally, we discuss potential strategies to enhance vaccine responsiveness in high-risk groups, focusing specifically on older adults.

CoRE-ATAC: A deep learning model for the functional classification of regulatory elements from single cell and bulk ATAC-seq data.

Cis-Regulatory elements (cis-REs) include promoters, enhancers, and insulators that regulate gene expression programs via binding of transcription factors. ATAC-seq technology effectively identifies active cis-REs in a given cell type (including from single cells) by mapping accessible chromatin at base-pair resolution. However, these maps are not immediately useful for inferring specific functions of cis-REs. For this purpose, we developed a deep learning framework (CoRE-ATAC) with novel data encoders that integrate DNA sequence (reference or personal genotypes) with ATAC-seq cut sites and read pileups. CoRE-ATAC was trained on 4 cell types (n = 6 samples/replicates) and accurately predicted known cis-RE functions from 7 cell types (n = 40 samples) that were not used in model training (mean average precision = 0.80, mean F1 score = 0.70). CoRE-ATAC enhancer predictions from 19 human islet samples coincided with genetically modulated gain/loss of enhancer activity, which was confirmed by massively parallel reporter assays (MPRAs). Finally, CoRE-ATAC effectively inferred cis-RE function from aggregate single nucleus ATAC-seq (snATAC) data from human blood-derived immune cells that overlapped with known functional annotations in sorted immune cells, which established the efficacy of these models to study cis-RE functions of rare cells without the need for cell sorting. ATAC-seq maps from primary human cells reveal individual- and cell-specific variation in cis-RE activity. CoRE-ATAC increases the functional resolution of these maps, a critical step for studying regulatory disruptions behind diseases.

A Common Type 2 Diabetes Risk Variant Potentiates Activity of an Evolutionarily Conserved Islet Stretch Enhancer and Increases C2CD4A and C2CD4B Expression.

Genome-wide association studies (GWASs) and functional genomics approaches implicate enhancer disruption in islet dysfunction and type 2 diabetes (T2D) risk. We applied genetic fine-mapping and functional (epi)genomic approaches to a T2D- and proinsulin-associated 15q22.2 locus to identify a most likely causal variant, determine its direction of effect, and elucidate plausible target genes. Fine-mapping and conditional analyses of proinsulin levels of 8,635 non-diabetic individuals from the METSIM study support a single association signal represented by a cluster of 16 strongly associated (p < 10-17) variants in high linkage disequilibrium (r2 > 0.8) with the GWAS index SNP rs7172432. These variants reside in an evolutionarily and functionally conserved islet and ß cell stretch or super enhancer; the most strongly associated variant (rs7163757, p = 3 × 10-19) overlaps a conserved islet open chromatin site. DNA sequence containing the rs7163757 risk allele displayed 2-fold higher enhancer activity than the non-risk allele in reporter assays (p < 0.01) and was differentially bound by ß cell nuclear extract proteins. Transcription factor NFAT specifically potentiated risk-allele enhancer activity and altered patterns of nuclear protein binding to the risk allele in vitro, suggesting that it could be a factor mediating risk-allele effects. Finally, the rs7163757 proinsulin-raising and T2D risk allele (C) was associated with increased expression of C2CD4B, and possibly C2CD4A, both of which were induced by inflammatory cytokines, in human islets. Together, these data suggest that rs7163757 contributes to genetic risk of islet dysfunction and T2D by increasing NFAT-mediated islet enhancer activity and modulating C2CD4B, and possibly C2CD4A, expression in (patho)physiologic states.

V-SVA: an R Shiny application for detecting and annotating hidden sources of variation in single-cell RNA-seq data.

SUMMARY: Single-cell RNA-sequencing (scRNA-seq) technology enables studying gene expression programs from individual cells. However, these data are subject to diverse sources of variation, including 'unwanted' variation that needs to be removed in downstream analyses (e.g. batch effects) and 'wanted' or biological sources of variation (e.g. variation associated with a cell type) that needs to be precisely described. Surrogate variable analysis (SVA)-based algorithms, are commonly used for batch correction and more recently for studying 'wanted' variation in scRNA-seq data. However, interpreting whether these variables are biologically meaningful or stemming from technical reasons remains a challenge. To facilitate the interpretation of surrogate variables detected by algorithms including IA-SVA, SVA or ZINB-WaVE, we developed an R Shiny application [Visual Surrogate Variable Analysis (V-SVA)] that provides a web-browser interface for the identification and annotation of hidden sources of variation in scRNA-seq data. This interactive framework includes tools for discovery of genes associated with detected sources of variation, gene annotation using publicly available databases and gene sets, and data visualization using dimension reduction methods. AVAILABILITY AND IMPLEMENTATION: The V-SVA Shiny application is publicly hosted at https://vsva.jax.org/ and the source code is freely available at https://github.com/nlawlor/V-SVA. CONTACT: leed13@miamioh.edu or duygu.ucar@jax.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

BiFET: sequencing Bias-free transcription factor Footprint Enrichment Test.

Transcription factor (TF) footprinting uncovers putative protein-DNA binding via combined analyses of chromatin accessibility patterns and their underlying TF sequence motifs. TF footprints are frequently used to identify TFs that regulate activities of cell/condition-specific genomic regions (target loci) in comparison to control regions (background loci) using standard enrichment tests. However, there is a strong association between the chromatin accessibility level and the GC content of a locus and the number and types of TF footprints that can be detected at this site. Traditional enrichment tests (e.g. hypergeometric) do not account for this bias and inflate false positive associations. Therefore, we developed a novel post-processing method, Bias-free Footprint Enrichment Test (BiFET), that corrects for the biases arising from the differences in chromatin accessibility levels and GC contents between target and background loci in footprint enrichment analyses. We applied BiFET on TF footprint calls obtained from EndoC-ßH1 ATAC-seq samples using three different algorithms (CENTIPEDE, HINT-BC and PIQ) and showed BiFET's ability to increase power and reduce false positive rate when compared to hypergeometric test. Furthermore, we used BiFET to study TF footprints from human PBMC and pancreatic islet ATAC-seq samples to show its utility to identify putative TFs associated with cell-type-specific loci.

Genomics of Islet (Dys)function and Type 2 Diabetes.

Pancreatic islet dysfunction and beta cell failure are hallmarks of type 2 diabetes mellitus (T2DM) pathogenesis. In this review, we discuss how genome-wide association studies (GWASs) and recent developments in islet (epi)genome and transcriptome profiling (particularly single cell analyses) are providing novel insights into the genetic, environmental, and cellular contributions to islet (dys)function and T2DM pathogenesis. Moving forward, study designs that interrogate and model genetic variation [e.g., allelic profiling and (epi)genome editing] will be critical to dissect the molecular genetics of T2DM pathogenesis, to build next-generation cellular and animal models, and to develop precision medicine approaches to detect, treat, and prevent islet (dys)function and T2DM.

The lethal sex gap: COVID-19.

While Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is disrupting lives across the globe for everyone, it has a more devastating impact on the health of older adults, especially that of older men. This pandemic has highlighted the crucial importance of considering an individual's age and biological sex in the clinic in addition to other confounding diseases (Kuchel, G.A, J Am Geriatr Soc, 67, 203, 2019, Tannenbaum, C., Nature, 575 451-458, 2009) As an interdisciplinary team of scientists in immunology, hematology, genomics, bioinformatics, and geriatrics, we have been studying how age and sex shape the human immune system. Herein we reflect on how our recent findings on the alterations of the immune system in aging might contribute to our current understanding of COVID-19 infection rate and disease risk.

FoxO6 regulates memory consolidation and synaptic function.

The FoxO family of transcription factors is known to slow aging downstream from the insulin/IGF (insulin-like growth factor) signaling pathway. The most recently discovered FoxO isoform in mammals, FoxO6, is highly enriched in the adult hippocampus. However, the importance of FoxO factors in cognition is largely unknown. Here we generated mice lacking FoxO6 and found that these mice display normal learning but impaired memory consolidation in contextual fear conditioning and novel object recognition. Using stereotactic injection of viruses into the hippocampus of adult wild-type mice, we found that FoxO6 activity in the adult hippocampus is required for memory consolidation. Genome-wide approaches revealed that FoxO6 regulates a program of genes involved in synaptic function upon learning in the hippocampus. Consistently, FoxO6 deficiency results in decreased dendritic spine density in hippocampal neurons in vitro and in vivo. Thus, FoxO6 may promote memory consolidation by regulating a program coordinating neuronal connectivity in the hippocampus, which could have important implications for physiological and pathological age-dependent decline in memory.

The tandem duplicator phenotype as a distinct genomic configuration in cancer.

Next-generation sequencing studies have revealed genome-wide structural variation patterns in cancer, such as chromothripsis and chromoplexy, that do not engage a single discernable driver mutation, and whose clinical relevance is unclear. We devised a robust genomic metric able to identify cancers with a chromotype called tandem duplicator phenotype (TDP) characterized by frequent and distributed tandem duplications (TDs). Enriched only in triple-negative breast cancer (TNBC) and in ovarian, endometrial, and liver cancers, TDP tumors conjointly exhibit tumor protein p53 (TP53) mutations, disruption of breast cancer 1 (BRCA1), and increased expression of DNA replication genes pointing at rereplication in a defective checkpoint environment as a plausible causal mechanism. The resultant TDs in TDP augment global oncogene expression and disrupt tumor suppressor genes. Importantly, the TDP strongly correlates with cisplatin sensitivity in both TNBC cell lines and primary patient-derived xenografts. We conclude that the TDP is a common cancer chromotype that coordinately alters oncogene/tumor suppressor expression with potential as a marker for chemotherapeutic response.

QuIN: A Web Server for Querying and Visualizing Chromatin Interaction Networks.

UNLABELLED: Recent studies of the human genome have indicated that regulatory elements (e.g. promoters and enhancers) at distal genomic locations can interact with each other via chromatin folding and affect gene expression levels. Genomic technologies for mapping interactions between DNA regions, e.g., ChIA-PET and HiC, can generate genome-wide maps of interactions between regulatory elements. These interaction datasets are important resources to infer distal gene targets of non-coding regulatory elements and to facilitate prioritization of critical loci for important cellular functions. With the increasing diversity and complexity of genomic information and public ontologies, making sense of these datasets demands integrative and easy-to-use software tools. Moreover, network representation of chromatin interaction maps enables effective data visualization, integration, and mining. Currently, there is no software that can take full advantage of network theory approaches for the analysis of chromatin interaction datasets. To fill this gap, we developed a web-based application, QuIN, which enables: 1) building and visualizing chromatin interaction networks, 2) annotating networks with user-provided private and publicly available functional genomics and interaction datasets, 3) querying network components based on gene name or chromosome location, and 4) utilizing network based measures to identify and prioritize critical regulatory targets and their direct and indirect interactions. AVAILABILITY: QuIN's web server is available at http://quin.jax.org QuIN is developed in Java and JavaScript, utilizing an Apache Tomcat web server and MySQL database and the source code is available under the GPLV3 license available on GitHub: https://github.com/UcarLab/QuIN/.

Transgenerational epigenetic inheritance of longevity in Caenorhabditis elegans.

Chromatin modifiers regulate lifespan in several organisms, raising the question of whether changes in chromatin states in the parental generation could be incompletely reprogrammed in the next generation and thereby affect the lifespan of descendants. The histone H3 lysine 4 trimethylation (H3K4me3) complex, composed of ASH-2, WDR-5 and the histone methyltransferase SET-2, regulates Caenorhabditis elegans lifespan. Here we show that deficiencies in the H3K4me3 chromatin modifiers ASH-2, WDR-5 or SET-2 in the parental generation extend the lifespan of descendants up until the third generation. The transgenerational inheritance of lifespan extension by members of the ASH-2 complex is dependent on the H3K4me3 demethylase RBR-2, and requires the presence of a functioning germline in the descendants. Transgenerational inheritance of lifespan is specific for the H3K4me3 methylation complex and is associated with epigenetic changes in gene expression. Thus, manipulation of specific chromatin modifiers only in parents can induce an epigenetic memory of longevity in descendants.

Developmental enhancers are marked independently of zygotic Nodal signals in Xenopus.

To determine the hierarchy of transcriptional regulation within the in vivo vertebrate embryo, we examined whether developmental enhancers were influenced by Nodal signaling during early embryogenesis in Xenopus tropicalis. We find that developmental enhancers, defined by the active enhancer chromatin marks H3K4me1 and H3K27ac, are established as early as blastula stage and that Smad2/3 only strongly associates with these regions at gastrula stages. Significantly, when we perturb Nodal signaling using the drug SB431542, most enhancers remain marked, including at genes known to be sensitive to Nodal signaling. Overall, as enhancers are in an active conformation prior to Nodal signaling and are established independently of Nodal signaling, we suggest that many developmental enhancers are marked maternally, prior to exposure to extrinsic signals.

Transcriptional Regulation of the Pancreatic Islet: Implications for Islet Function.

Islets of Langerhans contain multiple hormone-producing endocrine cells controlling glucose homeostasis. Transcription establishes and maintains islet cellular fates and identities. Genetic and environmental disruption of islet transcription triggers cellular dysfunction and disease. Early transcriptional regulation studies of specific islet genes, including insulin (INS) and the transcription factor PDX1, identified the first cis-regulatory DNA sequences and trans-acting factors governing islet function. Here, we review how human islet "omics" studies are reshaping our understanding of transcriptional regulation in islet (dys)function and diabetes. First, we highlight the expansion of islet transcript number, form, and function and of DNA transcriptional regulatory elements controlling their production. Next, we cover islet transcriptional effects of genetic and environmental perturbation. Finally, we discuss how these studies' emerging insights should empower our diabetes research community to build mechanistic understanding of diabetes pathophysiology and to equip clinicians with tailored, precision medicine options to prevent and treat islet dysfunction and diabetes.

Multi-omic human pancreatic islet endoplasmic reticulum and cytokine stress response mapping provides type 2 diabetes genetic insights.

Endoplasmic reticulum (ER) and inflammatory stress responses contribute to islet dysfunction in type 2 diabetes (T2D). Comprehensive genomic understanding of these human islet stress responses and whether T2D-associated genetic variants modulate them is lacking. Here, comparative transcriptome and epigenome analyses of human islets exposed ex vivo to these stressors revealed 30% of expressed genes and 14% of islet cis-regulatory elements (CREs) as stress responsive, modulated largely in an ER- or cytokine-specific fashion. T2D variants overlapped 86 stress-responsive CREs, including 21 induced by ER stress. We linked the rs6917676-T T2D risk allele to increased islet ER-stress-responsive CRE accessibility and allele-specific ß cell nuclear factor binding. MAP3K5, the ER-stress-responsive putative rs6917676 T2D effector gene, promoted stress-induced ß cell apoptosis. Supporting its pro-diabetogenic role, MAP3K5 expression correlated inversely with human islet ß cell abundance and was elevated in T2D ß cells. This study provides genome-wide insights into human islet stress responses and context-specific T2D variant effects.

Combinatorial chromatin modification patterns in the human genome revealed by subspace clustering.

Chromatin modifications, such as post-translational modification of histone proteins and incorporation of histone variants, play an important role in regulating gene expression. Joint analyses of multiple histone modification maps are starting to reveal combinatorial patterns of modifications that are associated with functional DNA elements, providing support to the 'histone code' hypothesis. However, due to the lack of analytical methods, only a small number of chromatin modification patterns have been discovered so far. Here, we introduce a scalable subspace clustering algorithm, coherent and shifted bicluster identification (CoSBI), to exhaustively identify the set of combinatorial modification patterns across a given epigenome. Performance comparisons demonstrate that CoSBI can generate biclusters with higher intra-cluster coherency and biological relevance. We apply our algorithm to a compendium of 39 genome-wide chromatin modification maps in human CD4(+) T cells. We identify 843 combinatorial patterns that recur at >0.1% of the genome. A total of 19 chromatin modifications are observed in the combinatorial patterns, 10 of which occur in more than half of the patterns. We also identify combinatorial modification signatures for eight classes of functional DNA elements. Application of CoSBI to epigenome maps of different cells and developmental stages will aid in understanding how chromatin structure helps regulate gene expression.

Transcriptional activation of Jun and Fos members of the AP-1 complex is a conserved signature of immune aging that contributes to inflammaging.

Diverse mouse strains have different health and life spans, mimicking the diversity among humans. To capture conserved aging signatures, we studied long-lived C57BL/6J and short-lived NZO/HILtJ mouse strains by profiling transcriptomes and epigenomes of immune cells from peripheral blood and the spleen from young and old mice. Transcriptional activation of the AP-1 transcription factor complex, particularly Fos, Junb, and Jun genes, was the most significant and conserved aging signature across tissues and strains. ATAC-seq data analyses showed that the chromatin around these genes was more accessible with age and there were significantly more binding sites for these TFs with age across all studied tissues, targeting pro-inflammatory molecules including Il6. Age-related increases in binding sites of JUN and FOS factors were also conserved in human peripheral blood ATAC-seq data. Single-cell RNA-seq data from the mouse aging cell atlas Tabula Muris Senis showed that the expression of these genes increased with age in B, T, NK cells, and macrophages, with macrophages from old mice expressing these molecules more abundantly than other cells. Functional data showed that upon myeloid cell activation via poly(I:C), the levels of JUN protein and its binding activity increased more significantly in spleen cells from old compared to young mice. In addition, upon activation, old cells produced more IL6 compared to young cells. In sum, we showed that the aging-related transcriptional activation of Jun and Fos family members in AP-1 complex is conserved across immune tissues and long- and short-living mouse strains, possibly contributing to increased inflammation with age.