Imputation of 3D genome structure by genetic-epigenetic interaction modeling in mice.

Gene expression is known to be affected by interactions between local genetic variation and DNA accessibility, with the latter organized into three-dimensional chromatin structures. Analyses of these interactions have previously been limited, obscuring their regulatory context, and the extent to which they occur throughout the genome. Here, we undertake a genome-scale analysis of these interactions in a genetically diverse population to systematically identify global genetic-epigenetic interaction, and reveal constraints imposed by chromatin structure. We establish the extent and structure of genotype-by-epigenotype interaction using embryonic stem cells derived from Diversity Outbred mice. This mouse population segregates millions of variants from eight inbred founders, enabling precision genetic mapping with extensive genotypic and phenotypic diversity. With 176 samples profiled for genotype, gene expression, and open chromatin, we used regression modeling to infer genetic-epigenetic interactions on a genome-wide scale. Our results demonstrate that statistical interactions between genetic variants and chromatin accessibility are common throughout the genome. We found that these interactions occur within the local area of the affected gene, and that this locality corresponds to topologically associated domains (TADs). The likelihood of interaction was most strongly defined by the three-dimensional (3D) domain structure rather than linear DNA sequence. We show that stable 3D genome structure is an effective tool to guide searches for regulatory elements and, conversely, that regulatory elements in genetically diverse populations provide a means to infer 3D genome structure. We confirmed this finding with CTCF ChIP-seq that revealed strain-specific binding in the inbred founder mice. In stem cells, open chromatin participating in the most significant regression models demonstrated an enrichment for developmental genes and the TAD-forming CTCF-binding complex, providing an opportunity for statistical inference of shifting TAD boundaries operating during early development. These findings provide evidence that genetic and epigenetic factors operate within the context of 3D chromatin structure.

Constitutive expression of interleukin-27 diminishes proinflammatory cytokine production without impairing effector function of engineered T cells.

Immunomodulatory cytokines can alter the tumor microenvironment and promote tumor eradication. Interleukin (IL)-27 is a pleiotropic cytokine that has potential to augment anti-tumor immunity while also facilitating anti-myeloma activity. We engineered human T cells to express a recombinant single-chain (sc)IL-27 and a synthetic antigen receptor targeting the myeloma antigen, B-cell maturation antigen, and evaluated the anti-tumor function of T cells bearing scIL-27 in vitro and in vivo. We discovered that T cells bearing scIL-27 sustained anti-tumor immunity and cytotoxicity yet manifested a profound reduction in pro-inflammatory cytokines granulocyte-macrophage colony-stimulating factor and tumor necrosis factor alpha. IL-27-expressing T cells therefore present a potential avenue to avert treatment-related toxicities commonly associated with engineered T-cell therapy due to the reduced pro-inflammatory cytokine profile.

Variation in histone configurations correlates with gene expression across nine inbred strains of mice.

The Diversity Outbred (DO) mice and their inbred founders are widely used models of human disease. However, although the genetic diversity of these mice has been well documented, their epigenetic diversity has not. Epigenetic modifications, such as histone modifications and DNA methylation, are important regulators of gene expression and, as such, are a critical mechanistic link between genotype and phenotype. Therefore, creating a map of epigenetic modifications in the DO mice and their founders is an important step toward understanding mechanisms of gene regulation and the link to disease in this widely used resource. To this end, we performed a strain survey of epigenetic modifications in hepatocytes of the DO founders. We surveyed four histone modifications (H3K4me1, H3K4me3, H3K27me3, and H3K27ac), as well as DNA methylation. We used ChromHMM to identify 14 chromatin states, each of which represents a distinct combination of the four histone modifications. We found that the epigenetic landscape is highly variable across the DO founders and is associated with variation in gene expression across strains. We found that epigenetic state imputed into a population of DO mice recapitulated the association with gene expression seen in the founders, suggesting that both histone modifications and DNA methylation are highly heritable mechanisms of gene expression regulation. We illustrate how DO gene expression can be aligned with inbred epigenetic states to identify putative cis-regulatory regions. Finally, we provide a data resource that documents strain-specific variation in the chromatin state and DNA methylation in hepatocytes across nine widely used strains of laboratory mice.

LCL161 enhances expansion and survival of engineered anti-tumor T cells but is restricted by death signaling.

Background: The genesis of SMAC mimetic drugs is founded on the observation that many cancers amplify IAP proteins to facilitate their survival, and therefore removal of these pathways would re-sensitize the cells towards apoptosis. It has become increasingly clear that SMAC mimetics also interface with the immune system in a modulatory manner. Suppression of IAP function by SMAC mimetics activates the non-canonical NF-κB pathway which can augment T cell function, opening the possibility of using SMAC mimetics to enhance immunotherapeutics. Methods: We have investigated the SMAC mimetic LCL161, which promotes degradation of cIAP-1 and cIAP-2, as an agent for delivering transient costimulation to engineered BMCA-specific human TAC T cells. In doing so we also sought to understand the cellular and molecular effects of LCL161 on T cell biology. Results: LCL161 activated the non-canonical NF-κB pathway and enhanced antigen-driven TAC T cell proliferation and survival. Transcriptional profiling from TAC T cells treated with LCL161 revealed differential expression of costimulatory and apoptosis-related proteins, namely CD30 and FAIM3. We hypothesized that regulation of these genes by LCL161 may influence the drug's effects on T cells. We reversed the differential expression through genetic engineering and observed impaired costimulation by LCL161, particularly when CD30 was deleted. While LCL161 can provide a costimulatory signal to TAC T cells following exposure to isolated antigen, we did not observe a similar pattern when TAC T cells were stimulated with myeloma cells expressing the target antigen. We questioned whether FasL expression by myeloma cells may antagonize the costimulatory effects of LCL161. Fas-KO TAC T cells displayed superior expansion following antigen stimulation in the presence of LCL161, suggesting a role for Fas-related T cell death in limiting the magnitude of the T cell response to antigen in the presence of LCL161. Conclusions: Our results demonstrate that LCL161 provides costimulation to TAC T cells exposed to antigen alone, however LCL161 did not enhance TAC T cell anti-tumor function when challenged with myeloma cells and may be limited due to sensitization of T cells towards Fas-mediated apoptosis.

Enhlink infers distal and context-specific enhancer-promoter linkages.

Enhancers play a crucial role in regulating gene expression and their functional status can be queried with cell type precision using using single-cell (sc)ATAC-seq. To facilitate analysis of such data, we developed Enhlink, a novel computational approach that leverages single-cell signals to infer linkages between regulatory DNA sequences, such as enhancers and promoters. Enhlink uses an ensemble strategy that integrates cell-level technical covariates to control for batch effects and biological covariates to infer robust condition-specific links and their associated p-values. It can integrate simultaneous gene expression and chromatin accessibility measurements of individual cells profiled by multi-omic experiments for increased specificity. We evaluated Enhlink using simulated and real scATAC-seq data, including those paired with physical enhancer-promoter links enumerated by promoter capture Hi-C and with multi-omic scATAC-/RNA-seq data we generated from the mouse striatum. These examples demonstrated that our method outperforms popular alternative strategies. In conjunction with eQTL analysis, Enhlink revealed a putative super-enhancer regulating key cell type-specific markers of striatal neurons. Taken together, our analyses demonstrate that Enhlink is accurate, powerful, and provides features that can lead to novel biological insights.

Genetic dissection of the pluripotent proteome through multi-omics data integration.

Genetic background drives phenotypic variability in pluripotent stem cells (PSCs). Most studies to date have used transcript abundance as the primary molecular readout of cell state in PSCs. We performed a comprehensive proteogenomics analysis of 190 genetically diverse mouse embryonic stem cell (mESC) lines. The quantitative proteome is highly variable across lines, and we identified pluripotency-associated pathways that were differentially activated in the proteomics data that were not evident in transcriptome data from the same lines. Integration of protein abundance to transcript levels and chromatin accessibility revealed broad co-variation across molecular layers as well as shared and unique drivers of quantitative variation in pluripotency-associated pathways. Quantitative trait locus (QTL) mapping localized the drivers of these multi-omic signatures to genomic hotspots. This study reveals post-transcriptional mechanisms and genetic interactions that underlie quantitative variability in the pluripotent proteome and provides a regulatory map for mESCs that can provide a basis for future mechanistic studies.

T-cell engineered with a fully humanized B-cell maturation antigen-specific T-cell antigen coupler receptor effectively target multiple myeloma.

B-cell maturation antigen (BCMA) is a clinically validated target for multiple myeloma. T-cell engineered with chimeric antigen receptors (CARs) directed against BCMA have demonstrated robust therapeutic activity in clinical trials, but toxicities remain a significant concern for a subset of patients, supporting continued investigation of other engineered T-cell platforms that may offer equal efficacy with an improved toxicity profile. The authors recently described a BCMA-specific, T-cell-centric synthetic antigen receptor, the T-cell antigen coupler (TAC) receptor, that can be used to engineer T-cell with robust anti-myeloma activity. Here the authors describe the creation of a fully humanized BCMA-specific TAC receptor. Single-chain variable fragments (scFvs) were developed from BCMA-specific F(ab)s that were identified in a fully human phage display library. Twenty-four configurations of the F(ab)s were evaluated in a medium-throughput screening using primary T-cell, and a single F(ab), TRAC 3625, emerged as the most robust following in vitro and in vivo evaluation. An optimized BCMA-specific TAC receptor was developed through iterations of the BCMA-TAC design that evaluated a next-generation TAC scaffold sequence, different domains connecting the TAC to the 3625 scFv and different orientations of the TRAC 3625 heavy and light variable regions.

Genetic control of the pluripotency epigenome determines differentiation bias in mouse embryonic stem cells.

Genetically diverse pluripotent stem cells display varied, heritable responses to differentiation cues. Here, we harnessed these disparities through derivation of mouse embryonic stem cells from the BXD genetic reference panel, along with C57BL/6J (B6) and DBA/2J (D2) parental strains, to identify loci regulating cell state transitions. Upon transition to formative pluripotency, B6 stem cells quickly dissolved naïve networks adopting gene expression modules indicative of neuroectoderm lineages, whereas D2 retained aspects of naïve pluripotency. Spontaneous formation of embryoid bodies identified divergent differentiation where B6 showed a propensity toward neuroectoderm and D2 toward definitive endoderm. Genetic mapping identified major trans-acting loci co-regulating chromatin accessibility and gene expression in both naïve and formative pluripotency. These loci distally modulated occupancy of pluripotency factors at hundreds of regulatory elements. One trans-acting locus on Chr 12 primarily impacted chromatin accessibility in embryonic stem cells, while in epiblast-like cells, the same locus subsequently influenced expression of genes enriched for neurogenesis, suggesting early chromatin priming. These results demonstrate genetically determined biases in lineage commitment and identify major regulators of the pluripotency epigenome.

Development of a B-cell maturation antigen-specific T-cell antigen coupler receptor for multiple myeloma.

BACKGROUND AIMS: T cells engineered with synthetic receptors have delivered powerful therapeutic results for patients with relapsed/refractory hematologic malignancies. The authors have recently described the T-cell antigen coupler (TAC) receptor, which co-opts the endogenous T-cell receptor (TCR) and activates engineered T cells in an HLA-independent manner. Here the authors describe the evolution of a next-generation TAC receptor with a focus on developing a TAC-engineered T cell for multiple myeloma. METHODS: To optimize the TAC scaffold, the authors employed a bona fide antigen-binding domain derived from the B-cell maturation antigen-specific monoclonal antibody C11D5.3, which has been used successfully in the clinic. The authors first tested humanized versions of the UCHT1 domain, which is used by the TAC to co-opt the TCR. The authors further discovered that the signal peptide affected surface expression of the TAC receptor. Higher density of the TAC receptor enhanced target binding in vitro, which translated into higher levels of Lck at the immunological synapse and stronger proliferation when only receptor-ligand interactions were present. RESULTS: The authors observed that the humanized UCHT1 improved surface expression and in vivo efficacy. Using TAC T cells derived from both healthy donors and multiple myeloma patients, the authors determined that despite the influence of receptor density on early activation events and effector function, receptor density did not impact late effector functions in vitro, nor did the receptor density affect in vivo efficacy. CONCLUSIONS: The modifications to the TAC scaffold described herein represent an important step in the evolution of this technology, which tolerates a range of expression levels without impacting therapeutic efficacy.

PRD-2 directly regulates casein kinase I and counteracts nonsense-mediated decay in the Neurospora circadian clock.

Circadian clocks in fungi and animals are driven by a functionally conserved transcription-translation feedback loop. In Neurospora crassa, negative feedback is executed by a complex of Frequency (FRQ), FRQ-interacting RNA helicase (FRH), and casein kinase I (CKI), which inhibits the activity of the clock's positive arm, the White Collar Complex (WCC). Here, we show that the prd-2 (period-2) gene, whose mutation is characterized by recessive inheritance of a long 26 hr period phenotype, encodes an RNA-binding protein that stabilizes the ck-1a transcript, resulting in CKI protein levels sufficient for normal rhythmicity. Moreover, by examining the molecular basis for the short circadian period of upf-1prd-6 mutants, we uncovered a strong influence of the Nonsense Mediated Decay pathway on CKI levels. The finding that circadian period defects in two classically derived Neurospora clock mutants each arise from disruption of ck-1a regulation is consistent with circadian period being exquisitely sensitive to levels of casein kinase I.

Naive Pluripotent Stem Cells Exhibit Phenotypic Variability that Is Driven by Genetic Variation.

Variability among pluripotent stem cell (PSC) lines is a prevailing issue that hampers not only experimental reproducibility but also large-scale applications and personalized cell-based therapy. This variability could result from epigenetic and genetic factors that influence stem cell behavior. Naive culture conditions minimize epigenetic fluctuation, potentially overcoming differences in PSC line differentiation potential. Here we derived PSCs from distinct mouse strains under naive conditions and show that lines from distinct genetic backgrounds have divergent differentiation capacity, confirming a major role for genetics in PSC phenotypic variability. This is explained in part through inconsistent activity of extra-cellular signaling, including the Wnt pathway, which is modulated by specific genetic variants. Overall, this study shows that genetic background plays a dominant role in driving phenotypic variability of PSCs.

Mapping the Effects of Genetic Variation on Chromatin State and Gene Expression Reveals Loci That Control Ground State Pluripotency.

Mouse embryonic stem cells (mESCs) cultured in the presence of LIF occupy a ground state with highly active pluripotency-associated transcriptional and epigenetic circuitry. However, ground state pluripotency in some inbred strain backgrounds is unstable in the absence of ERK1/2 and GSK3 inhibition. Using an unbiased genetic approach, we dissect the basis of this divergent response to extracellular cues by profiling gene expression and chromatin accessibility in 170 genetically heterogeneous mESCs. We map thousands of loci affecting chromatin accessibility and/or transcript abundance, including 10 QTL hotspots where genetic variation at a single locus coordinates the regulation of genes throughout the genome. For one hotspot, we identify a single enhancer variant ∼10 kb upstream of Lifr associated with chromatin accessibility and mediating a cascade of molecular events affecting pluripotency. We validate causation through reciprocal allele swaps, demonstrating the functional consequences of noncoding variation in gene regulatory networks that stabilize pluripotent states in vitro.

Unique properties of a subset of human pluripotent stem cells with high capacity for self-renewal.

Archetypal human pluripotent stem cells (hPSC) are widely considered to be equivalent in developmental status to mouse epiblast stem cells, which correspond to pluripotent cells at a late post-implantation stage of embryogenesis. Heterogeneity within hPSC cultures complicates this interspecies comparison. Here we show that a subpopulation of archetypal hPSC enriched for high self-renewal capacity (ESR) has distinct properties relative to the bulk of the population, including a cell cycle with a very low G1 fraction and a metabolomic profile that reflects a combination of oxidative phosphorylation and glycolysis. ESR cells are pluripotent and capable of differentiation into primordial germ cell-like cells. Global DNA methylation levels in the ESR subpopulation are lower than those in mouse epiblast stem cells. Chromatin accessibility analysis revealed a unique set of open chromatin sites in ESR cells. RNA-seq at the subpopulation and single cell levels shows that, unlike mouse epiblast stem cells, the ESR subset of hPSC displays no lineage priming, and that it can be clearly distinguished from gastrulating and extraembryonic cell populations in the primate embryo. ESR hPSC correspond to an earlier stage of post-implantation development than mouse epiblast stem cells.

HELLS and PRDM9 form a pioneer complex to open chromatin at meiotic recombination hot spots.

Chromatin barriers prevent spurious interactions between regulatory elements and DNA-binding proteins. One such barrier, whose mechanism for overcoming is poorly understood, is access to recombination hot spots during meiosis. Here we show that the chromatin remodeler HELLS and DNA-binding protein PRDM9 function together to open chromatin at hot spots and provide access for the DNA double-strand break (DSB) machinery. Recombination hot spots are decorated by a unique combination of histone modifications not found at other regulatory elements. HELLS is recruited to hot spots by PRDM9 and is necessary for both histone modifications and DNA accessibility at hot spots. In male mice lacking HELLS, DSBs are retargeted to other sites of open chromatin, leading to germ cell death and sterility. Together, these data provide a model for hot spot activation in which HELLS and PRDM9 form a pioneer complex to create a unique epigenomic environment of open chromatin, permitting correct placement and repair of DSBs.

Histone methyltransferase PRDM9 is not essential for meiosis in male mice.

A hallmark of meiosis is the rearrangement of parental alleles to ensure genetic diversity in the gametes. These chromosome rearrangements are mediated by the repair of programmed DNA double-strand breaks (DSBs) as genetic crossovers between parental homologs. In mice, humans, and many other mammals, meiotic DSBs occur primarily at hotspots, determined by sequence-specific binding of the PRDM9 protein. Without PRDM9, meiotic DSBs occur near gene promoters and other functional sites. Studies in a limited number of mouse strains showed that functional PRDM9 is required to complete meiosis, but despite its apparent importance, Prdm9 has been repeatedly lost across many animal lineages. Both the reason for mouse sterility in the absence of PRDM9 and the mechanism by which Prdm9 can be lost remain unclear. Here, we explore whether mice can tolerate the loss of Prdm9 By generating Prdm9 functional knockouts in an array of genetic backgrounds, we observe a wide range of fertility phenotypes and ultimately demonstrate that PRDM9 is not required for completion of male meiosis. Although DSBs still form at a common subset of functional sites in all mice lacking PRDM9, meiotic outcomes differ substantially. We speculate that DSBs at functional sites are difficult to repair as a crossover and that by increasing the efficiency of crossover formation at these sites, genetic modifiers of recombination rates can allow for meiotic progression. This model implies that species with a sufficiently high recombination rate may lose Prdm9 yet remain fertile.

Tissue-Specific Trans Regulation of the Mouse Epigenome.

The epigenetic landscape varies greatly among cell types. Although a variety of writers, readers, and erasers of epigenetic features are known, we have little information about the underlying regulatory systems controlling the establishment and maintenance of these features. Here, we have explored how natural genetic variation affects the epigenome in mice. Studying levels of H3K4me3, a histone modification at sites such as promoters, enhancers, and recombination hotspots, we found tissue-specific trans-regulation of H3K4me3 levels in four highly diverse cell types: male germ cells, embryonic stem cells, hepatocytes, and cardiomyocytes. To identify the genetic loci involved, we measured H3K4me3 levels in male germ cells in a mapping population of 59 BXD recombinant inbred lines. We found extensive trans-regulation of H3K4me3 peaks, including six major histone quantitative trait loci (QTL). These chromatin regulatory loci act dominantly to suppress H3K4me3, which at hotspots reduces the likelihood of subsequent DNA double-strand breaks. QTL locations do not correspond with genes encoding enzymes known to metabolize chromatin features. Instead their locations match clusters of zinc finger genes, making these possible candidates that explain the dominant suppression of H3K4me3. Collectively, these data describe an extensive, set of chromatin regulatory loci that control the epigenetic landscape.

Capturing Totipotent Stem Cells.

Minority subpopulations within embryonic stem cell cultures display an expanded developmental potential similar to that of early embryo blastomeres or the early inner cell mass. The ability to isolate and culture totipotent cells capable of giving rise to the entire conceptus would enhance our capacity to study early embryo development, and might enable more efficient generation of chimeric animals for research and organ production for transplantation. Here we review the biological and molecular characterization of cultured cells with developmental potential similar to totipotent blastomeres, and assess recent progress toward the capture and stabilization of the totipotent state in vitro.

The Meiotic Recombination Activator PRDM9 Trimethylates Both H3K36 and H3K4 at Recombination Hotspots In Vivo.

In many mammals, including humans and mice, the zinc finger histone methyltransferase PRDM9 performs the first step in meiotic recombination by specifying the locations of hotspots, the sites of genetic recombination. PRDM9 binds to DNA at hotspots through its zinc finger domain and activates recombination by trimethylating histone H3K4 on adjacent nucleosomes through its PR/SET domain. Recently, the isolated PR/SET domain of PRDM9 was shown capable of also trimethylating H3K36 in vitro, raising the question of whether this reaction occurs in vivo during meiosis, and if so, what its function might be. Here, we show that full-length PRDM9 does trimethylate H3K36 in vivo in mouse spermatocytes. Levels of H3K4me3 and H3K36me3 are highly correlated at hotspots, but mutually exclusive elsewhere. In vitro, we find that although PRDM9 trimethylates H3K36 much more slowly than it does H3K4, PRDM9 is capable of placing both marks on the same histone molecules. In accord with these results, we also show that PRDM9 can trimethylate both K4 and K36 on the same nucleosomes in vivo, but the ratio of K4me3/K36me3 is much higher for the pair of nucleosomes adjacent to the PRDM9 binding site compared to the next pair further away. Importantly, H3K4me3/H3K36me3-double-positive nucleosomes occur only in regions of recombination: hotspots and the pseudoautosomal (PAR) region of the sex chromosomes. These double-positive nucleosomes are dramatically reduced when PRDM9 is absent, showing that this signature is PRDM9-dependent at hotspots; the residual double-positive nucleosomes most likely come from the PRDM9-independent PAR. These results, together with the fact that PRDM9 is the only known mammalian histone methyltransferase with both H3K4 and H3K36 trimethylation activity, suggest that trimethylation of H3K36 plays an important role in the recombination process. Given the known requirement of H3K36me3 for double strand break repair by homologous recombination in somatic cells, we suggest that it may play the same role in meiosis.

Health and population effects of rare gene knockouts in adult humans with related parents.

Examining complete gene knockouts within a viable organism can inform on gene function. We sequenced the exomes of 3222 British adults of Pakistani heritage with high parental relatedness, discovering 1111 rare-variant homozygous genotypes with predicted loss of function (knockouts) in 781 genes. We observed 13.7% fewer homozygous knockout genotypes than we expected, implying an average load of 1.6 recessive-lethal-equivalent loss-of-function (LOF) variants per adult. When genetic data were linked to the individuals' lifelong health records, we observed no significant relationship between gene knockouts and clinical consultation or prescription rate. In this data set, we identified a healthy PRDM9-knockout mother and performed phased genome sequencing on her, her child, and control individuals. Our results show that meiotic recombination sites are localized away from PRDM9-dependent hotspots. Thus, natural LOF variants inform on essential genetic loci and demonstrate PRDM9 redundancy in humans.

Multimer Formation Explains Allelic Suppression of PRDM9 Recombination Hotspots.

Genetic recombination during meiosis functions to increase genetic diversity, promotes elimination of deleterious alleles, and helps assure proper segregation of chromatids. Mammalian recombination events are concentrated at specialized sites, termed hotspots, whose locations are determined by PRDM9, a zinc finger DNA-binding histone methyltransferase. Prdm9 is highly polymorphic with most alleles activating their own set of hotspots. In populations exhibiting high frequencies of heterozygosity, questions remain about the influences different alleles have in heterozygous individuals where the two variant forms of PRDM9 typically do not activate equivalent populations of hotspots. We now find that, in addition to activating its own hotspots, the presence of one Prdm9 allele can modify the activity of hotspots activated by the other allele. PRDM9 function is also dosage sensitive; Prdm9+/- heterozygous null mice have reduced numbers and less active hotspots and increased numbers of aberrant germ cells. In mice carrying two Prdm9 alleles, there is allelic competition; the stronger Prdm9 allele can partially or entirely suppress chromatin modification and recombination at hotspots of the weaker allele. In cell cultures, PRDM9 protein variants form functional heteromeric complexes which can bind hotspots sequences. When a heteromeric complex binds at a hotspot of one PRDM9 variant, the other PRDM9 variant, which would otherwise not bind, can still methylate hotspot nucleosomes. We propose that in heterozygous individuals the underlying molecular mechanism of allelic suppression results from formation of PRDM9 heteromers, where the DNA binding activity of one protein variant dominantly directs recombination initiation towards its own hotspots, effectively titrating down recombination by the other protein variant. In natural populations with many heterozygous individuals, allelic competition will influence the recombination landscape.