Tumor-associated macrophages trigger MAIT cell dysfunction at the HCC invasive margin.

Mucosal-associated invariant T (MAIT) cells represent an abundant innate-like T cell subtype in the human liver. MAIT cells are assigned crucial roles in regulating immunity and inflammation, yet their role in liver cancer remains elusive. Here, we present a MAIT cell-centered profiling of hepatocellular carcinoma (HCC) using scRNA-seq, flow cytometry, and co-detection by indexing (CODEX) imaging of paired patient samples. These analyses highlight the heterogeneity and dysfunctionality of MAIT cells in HCC and their defective capacity to infiltrate liver tumors. Machine-learning tools were used to dissect the spatial cellular interaction network within the MAIT cell neighborhood. Co-localization in the adjacent liver and interaction between niche-occupying CSF1R+PD-L1+ tumor-associated macrophages (TAMs) and MAIT cells was identified as a key regulatory element of MAIT cell dysfunction. Perturbation of this cell-cell interaction in ex vivo co-culture studies using patient samples and murine models reinvigorated MAIT cell cytotoxicity. These studies suggest that aPD-1/aPD-L1 therapies target MAIT cells in HCC patients.

Transcription factor Tox2 is required for metabolic adaptation and tissue residency of ILC3 in the gut.

Group 3 innate lymphoid cells (ILC3) are the major subset of gut-resident ILC with essential roles in infections and tissue repair, but how they adapt to the gut environment to maintain tissue residency is unclear. We report that Tox2 is critical for gut ILC3 maintenance and function. Gut ILC3 highly expressed Tox2, and depletion of Tox2 markedly decreased ILC3 in gut but not at central sites, resulting in defective control of Citrobacter rodentium infection. Single-cell transcriptional profiling revealed decreased expression of Hexokinase-2 in Tox2-deficient gut ILC3. Consistent with the requirement for hexokinases in glycolysis, Tox2-/- ILC3 displayed decreased ability to utilize glycolysis for protein translation. Ectopic expression of Hexokinase-2 rescued Tox2-/- gut ILC3 defects. Hypoxia and interleukin (IL)-17A each induced Tox2 expression in ILC3, suggesting a mechanism by which ILC3 adjusts to fluctuating environments by programming glycolytic metabolism. Our results reveal the requirement for Tox2 to support the metabolic adaptation of ILC3 within the gastrointestinal tract.

Genome-scale exon perturbation screens uncover exons critical for cell fitness.

CRISPR-Cas technology has transformed functional genomics, yet understanding of how individual exons differentially shape cellular phenotypes remains limited. Here, we optimized and conducted massively parallel exon deletion and splice-site mutation screens in human cell lines to identify exons that regulate cellular fitness. Fitness-promoting exons are prevalent in essential and highly expressed genes and commonly overlap with protein domains and interaction interfaces. Conversely, fitness-suppressing exons are enriched in nonessential genes, exhibiting lower inclusion levels, and overlap with intrinsically disordered regions and disease-associated mutations. In-depth mechanistic investigation of the screen-hit TAF5 alternative exon-8 revealed that its inclusion is required for assembly of the TFIID general transcription initiation complex, thereby regulating global gene expression output. Collectively, our orthogonal exon perturbation screens established a comprehensive repository of phenotypically important exons and uncovered regulatory mechanisms governing cellular fitness and gene expression.

hRpn13 shapes the proteome and transcriptome through epigenetic factors HDAC8, PADI4, and transcription factor NF-κB p50.

The anti-cancer target hRpn13 is a proteasome substrate receptor. However, hRpn13-targeting molecules do not impair its interaction with proteasomes or ubiquitin, suggesting other critical cellular activities. We find that hRpn13 depletion causes correlated proteomic and transcriptomic changes, with pronounced effects in myeloma cells for cytoskeletal and immune response proteins and bone-marrow-specific arginine deiminase PADI4. Moreover, a PROTAC against hRpn13 co-depletes PADI4, histone deacetylase HDAC8, and DNA methyltransferase MGMT. PADI4 binds and citrullinates hRpn13 and proteasomes, and proteasomes from PADI4-inhibited myeloma cells exhibit reduced peptidase activity. When off proteasomes, hRpn13 can bind HDAC8, and this interaction inhibits HDAC8 activity. Further linking hRpn13 to transcription, its loss reduces nuclear factor κB (NF-κB) transcription factor p50, which proteasomes generate by cleaving its precursor protein. NF-κB inhibition depletes hRpn13 interactors PADI4 and HDAC8. Altogether, we find that hRpn13 acts dually in protein degradation and expression and that proteasome constituency and, in turn, regulation varies by cell type.

Comprehensive mapping of cell fates in microsatellite unstable cancer cells supports dual targeting of WRN and ATR.

Addiction to the WRN helicase is a unique vulnerability of human cancers with high levels of microsatellite instability (MSI-H). However, while prolonged loss of WRN ultimately leads to cell death, little is known about how MSI-H cancers initially respond to acute loss of WRN-knowledge that would be helpful for informing clinical development of WRN targeting therapy, predicting possible resistance mechanisms, and identifying useful biomarkers of successful WRN inhibition. Here, we report the construction of an inducible ligand-mediated degradation system in which the stability of endogenous WRN protein can be rapidly and specifically tuned, enabling us to track the complete sequence of cellular events elicited by acute loss of WRN function. We found that WRN degradation leads to immediate accrual of DNA damage in a replication-dependent manner that curiously did not robustly engage checkpoint mechanisms to halt DNA synthesis. As a result, WRN-degraded MSI-H cancer cells accumulate DNA damage across multiple replicative cycles and undergo successive rounds of increasingly aberrant mitoses, ultimately triggering cell death. Of potential therapeutic importance, we found no evidence of any generalized mechanism by which MSI-H cancers could adapt to near-complete loss of WRN. However, under conditions of partial WRN degradation, addition of low-dose ATR inhibitor significantly increased their combined efficacy to levels approaching full inactivation of WRN. Overall, our results provide the first comprehensive view of molecular events linking upstream inhibition of WRN to subsequent cell death and suggest that dual targeting of WRN and ATR might be a useful strategy for treating MSI-H cancers.

The dystonia gene THAP1 controls DNA double-strand break repair choice.

The Shieldin complex shields double-strand DNA breaks (DSBs) from nucleolytic resection. Curiously, the penultimate Shieldin component, SHLD1, is one of the least abundant mammalian proteins. Here, we report that the transcription factors THAP1, YY1, and HCF1 bind directly to the SHLD1 promoter, where they cooperatively maintain the low basal expression of SHLD1, thereby ensuring a proper balance between end protection and resection during DSB repair. The loss of THAP1-dependent SHLD1 expression confers cross-resistance to poly (ADP-ribose) polymerase (PARP) inhibitor and cisplatin in BRCA1-deficient cells and shorter progression-free survival in ovarian cancer patients. Moreover, the embryonic lethality and PARPi sensitivity of BRCA1-deficient mice is rescued by ablation of SHLD1. Our study uncovers a transcriptional network that directly controls DSB repair choice and suggests a potential link between DNA damage and pathogenic THAP1 mutations, found in patients with the neurodevelopmental movement disorder adult-onset torsion dystonia type 6.

DNAJC9 prevents CENP-A mislocalization and chromosomal instability by maintaining the fidelity of histone supply chains.

The centromeric histone H3 variant CENP-A is overexpressed in many cancers. The mislocalization of CENP-A to noncentromeric regions contributes to chromosomal instability (CIN), a hallmark of cancer. However, pathways that promote or prevent CENP-A mislocalization remain poorly defined. Here, we performed a genome-wide RNAi screen for regulators of CENP-A localization which identified DNAJC9, a J-domain protein implicated in histone H3-H4 protein folding, as a factor restricting CENP-A mislocalization. Cells lacking DNAJC9 exhibit mislocalization of CENP-A throughout the genome, and CIN phenotypes. Global interactome analysis showed that DNAJC9 depletion promotes the interaction of CENP-A with the DNA-replication-associated histone chaperone MCM2. CENP-A mislocalization upon DNAJC9 depletion was dependent on MCM2, defining MCM2 as a driver of CENP-A deposition at ectopic sites when H3-H4 supply chains are disrupted. Cells depleted for histone H3.3, also exhibit CENP-A mislocalization. In summary, we have defined novel factors that prevent mislocalization of CENP-A, and demonstrated that the integrity of H3-H4 supply chains regulated by histone chaperones such as DNAJC9 restrict CENP-A mislocalization and CIN.

Targeted long-read sequencing of the Ewing sarcoma 6p25.1 susceptibility locus identifies germline-somatic interactions with EWSR1-FLI1 binding.

Ewing sarcoma (EwS) is a rare bone and soft tissue malignancy driven by chromosomal translocations encoding chimeric transcription factors, such as EWSR1-FLI1, that bind GGAA motifs forming novel enhancers that alter nearby expression. We propose that germline microsatellite variation at the 6p25.1 EwS susceptibility locus could impact downstream gene expression and EwS biology. We performed targeted long-read sequencing of EwS blood DNA to characterize variation and genomic features important for EWSR1-FLI1 binding. We identified 50 microsatellite alleles at 6p25.1 and observed that EwS-affected individuals had longer alleles (>135 bp) with more GGAA repeats. The 6p25.1 GGAA microsatellite showed chromatin features of an EWSR1-FLI1 enhancer and regulated expression of RREB1, a transcription factor associated with RAS/MAPK signaling. RREB1 knockdown reduced proliferation and clonogenic potential and reduced expression of cell cycle and DNA replication genes. Our integrative analysis at 6p25.1 details increased binding of longer GGAA microsatellite alleles with acquired EWSR-FLI1 to promote Ewing sarcomagenesis by RREB1-mediated proliferation.

Saturation genome editing of 11 codons and exon 13 of BRCA2 coupled with chemotherapeutic drug response accurately determines pathogenicity of variants.

The unknown pathogenicity of a significant number of variants found in cancer-related genes is attributed to limited epidemiological data, resulting in their classification as variant of uncertain significance (VUS). To date, Breast Cancer gene-2 (BRCA2) has the highest number of VUSs, which has necessitated the development of several robust functional assays to determine their functional significance. Here we report the use of a humanized-mouse embryonic stem cell (mESC) line expressing a single copy of the human BRCA2 for a CRISPR-Cas9-based high-throughput functional assay. As a proof-of-principle, we have saturated 11 codons encoded by BRCA2 exons 3, 18, 19 and all possible single-nucleotide variants in exon 13 and multiplexed these variants for their functional categorization. Specifically, we used a pool of 180-mer single-stranded donor DNA to generate all possible combination of variants. Using a high throughput sequencing-based approach, we show a significant drop in the frequency of non-functional variants, whereas functional variants are enriched in the pool of the cells. We further demonstrate the response of these variants to the DNA-damaging agents, cisplatin and olaparib, allowing us to use cellular survival and drug response as parameters for variant classification. Using this approach, we have categorized 599 BRCA2 variants including 93-single nucleotide variants (SNVs) across the 11 codons, of which 28 are reported in ClinVar. We also functionally categorized 252 SNVs from exon 13 into 188 functional and 60 non-functional variants, demonstrating that saturation genome editing (SGE) coupled with drug sensitivity assays can enhance functional annotation of BRCA2 VUS.

Massively parallel reporter assays and variant scoring identified functional variants and target genes for melanoma loci and highlighted cell-type specificity.

The most recent genome-wide association study (GWAS) of cutaneous melanoma identified 54 risk-associated loci, but functional variants and their target genes for most have not been established. Here, we performed massively parallel reporter assays (MPRAs) by using malignant melanoma and normal melanocyte cells and further integrated multi-layer annotation to systematically prioritize functional variants and susceptibility genes from these GWAS loci. Of 1,992 risk-associated variants tested in MPRAs, we identified 285 from 42 loci (78% of the known loci) displaying significant allelic transcriptional activities in either cell type (FDR < 1%). We further characterized MPRA-significant variants by motif prediction, epigenomic annotation, and statistical/functional fine-mapping to create integrative variant scores, which prioritized one to six plausible candidate variants per locus for the 42 loci and nominated a single variant for 43% of these loci. Overlaying the MPRA-significant variants with genome-wide significant expression or methylation quantitative trait loci (eQTLs or meQTLs, respectively) from melanocytes or melanomas identified candidate susceptibility genes for 60% of variants (172 of 285 variants). CRISPRi of top-scoring variants validated their cis-regulatory effect on the eQTL target genes, MAFF (22q13.1) and GPRC5A (12p13.1). Finally, we identified 36 melanoma-specific and 45 melanocyte-specific MPRA-significant variants, a subset of which are linked to cell-type-specific target genes. Analyses of transcription factor availability in MPRA datasets and variant-transcription-factor interaction in eQTL datasets highlighted the roles of transcription factors in cell-type-specific variant functionality. In conclusion, MPRAs along with variant scoring effectively prioritized plausible candidates for most melanoma GWAS loci and highlighted cellular contexts where the susceptibility variants are functional.

The histone H3/H4 chaperone CHAF1B prevents the mislocalization of CENP-A for chromosomal stability.

Restricting the localization of the evolutionarily conserved centromeric histone H3 variant CENP-A to centromeres prevents chromosomal instability (CIN). The mislocalization of CENP-A to non-centromeric regions contributes to CIN in yeasts, flies and human cells. Even though overexpression and mislocalization of CENP-A have been reported in cancers, the mechanisms responsible for its mislocalization remain poorly understood. Here, we used an imaging-based high-throughput RNAi screen to identify factors that prevent mislocalization of overexpressed YFP-tagged CENP-A (YFP-CENP-A) in HeLa cells. Among the top five candidates in the screen - the depletion of which showed increased nuclear YFP-CENP-A fluorescence - were the histone chaperones CHAF1B (or p60) and CHAF1A (or p150). Follow-up validation and characterization experiments showed that CHAF1B-depleted cells exhibited CENP-A mislocalization, CIN phenotypes and increased enrichment of CENP-A in chromatin fractions. The depletion of DAXX, a histone H3.3 chaperone, suppressed CENP-A mislocalization and CIN in CHAF1B-depleted cells. We propose that in CHAF1B-depleted cells, DAXX promotes mislocalization of the overexpressed CENP-A to non-centromeric regions, resulting in CIN. In summary, we identified regulators of CENP-A localization and defined a role for CHAF1B in preventing DAXX-dependent CENP-A mislocalization and CIN.

An engineered cell line with a hRpn1-attached handle to isolate proteasomes.

Regulated protein degradation in eukaryotes is performed by the 26S proteasome, which contains a 19-subunit regulatory particle (RP) that binds, processes, and translocates substrates to a 28-subunit hollow core particle (CP) where proteolysis occurs. In addition to its intrinsic subunits, myriad proteins interact with the proteasome transiently, including factors that assist and/or regulate its degradative activities. Efforts to identify proteasome-interacting components and/or to solve its structure have relied on over-expression of a tagged plasmid, establishing stable cell lines, or laborious purification protocols to isolate native proteasomes from cells. Here, we describe an engineered human cell line, derived from colon cancer HCT116 cells, with a biotin handle on the RP subunit hRpn1/PSMD2 (proteasome 26S subunit, non-ATPase 2) for purification of 26S proteasomes. A 75-residue sequence from Propionibacterium shermanii that is biotinylated in mammalian cells was added following a tobacco etch virus protease cut site at the C terminus of hRpn1. We tested and found that 26S proteasomes can be isolated from this modified HCT116 cell line by using a simple purification protocol. More specifically, biotinylated proteasomes were purified from the cell lysates by using neutravidin agarose resin and released from the resin following incubation with tobacco etch virus protease. The purified proteasomes had equivalent activity in degrading a model ubiquitinated substrate, namely ubiquitinated p53, compared to commercially available bovine proteasomes that were purified by fractionation. In conclusion, advantages of this approach to obtain 26S proteasomes over others is the simple purification protocol and that all cellular proteins, including the tagged hRpn1 subunit, remain at endogenous stoichiometry.

A UVB-responsive common variant at chromosome band 7p21.1 confers tanning response and melanoma risk via regulation of the aryl hydrocarbon receptor, AHR.

Genome-wide association studies (GWASs) have identified a melanoma-associated locus on chromosome band 7p21.1 with rs117132860 as the lead SNP and a secondary independent signal marked by rs73069846. rs117132860 is also associated with tanning ability and cutaneous squamous cell carcinoma (cSCC). Because ultraviolet radiation (UVR) is a key environmental exposure for all three traits, we investigated the mechanisms by which this locus contributes to melanoma risk, focusing on cellular response to UVR. Fine-mapping of melanoma GWASs identified four independent sets of candidate causal variants. A GWAS region-focused Capture-C study of primary melanocytes identified physical interactions between two causal sets and the promoter of the aryl hydrocarbon receptor (AHR). Subsequent chromatin state annotation, eQTL, and luciferase assays identified rs117132860 as a functional variant and reinforced AHR as a likely causal gene. Because AHR plays critical roles in cellular response to dioxin and UVR, we explored links between this SNP and AHR expression after both 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and ultraviolet B (UVB) exposure. Allele-specific AHR binding to rs117132860-G was enhanced following both, consistent with predicted weakened AHR binding to the risk/poor-tanning rs117132860-A allele, and allele-preferential AHR expression driven from the protective rs117132860-G allele was observed following UVB exposure. Small deletions surrounding rs117132860 introduced via CRISPR abrogates AHR binding, reduces melanocyte cell growth, and prolongs growth arrest following UVB exposure. These data suggest AHR is a melanoma susceptibility gene at the 7p21.1 risk locus and rs117132860 is a functional variant within a UVB-responsive element, leading to allelic AHR expression and altering melanocyte growth phenotypes upon exposure.

A structurally conserved site in AUP1 binds the E2 enzyme UBE2G2 and is essential for ER-associated degradation.

Endoplasmic reticulum-associated degradation (ERAD) is a protein quality control pathway of fundamental importance to cellular homeostasis. Although multiple ERAD pathways exist for targeting topologically distinct substrates, all pathways require substrate ubiquitination. Here, we characterize a key role for the UBE2G2 Binding Region (G2BR) of the ERAD accessory protein ancient ubiquitous protein 1 (AUP1) in ERAD pathways. This 27-amino acid (aa) region of AUP1 binds with high specificity and low nanomolar affinity to the backside of the ERAD ubiquitin-conjugating enzyme (E2) UBE2G2. The structure of the AUP1 G2BR (G2BRAUP1) in complex with UBE2G2 reveals an interface that includes a network of salt bridges, hydrogen bonds, and hydrophobic interactions essential for AUP1 function in cells. The G2BRAUP1 shares significant structural conservation with the G2BR found in the E3 ubiquitin ligase gp78 and in vitro can similarly allosterically activate ubiquitination in conjunction with ERAD E3s. In cells, AUP1 is uniquely required to maintain normal levels of UBE2G2; this is due to G2BRAUP1 binding to the E2 and preventing its rapid degradation. In addition, the G2BRAUP1 is required for both ER membrane recruitment of UBE2G2 and for its activation at the ER membrane. Thus, by binding to the backside of a critical ERAD E2, G2BRAUP1 plays multiple critical roles in ERAD.

Beyond editing to writing large genomes.

Recent exponential advances in genome sequencing and engineering technologies have enabled an unprecedented level of interrogation into the impact of DNA variation (genotype) on cellular function (phenotype). Furthermore, these advances have also prompted realistic discussion of writing and radically re-writing complex genomes. In this Perspective, we detail the motivation for large-scale engineering, discuss the progress made from such projects in bacteria and yeast and describe how various genome-engineering technologies will contribute to this effort. Finally, we describe the features of an ideal platform and provide a roadmap to facilitate the efficient writing of large genomes.

dbGuide: a database of functionally validated guide RNAs for genome editing in human and mouse cells.

With the technology's accessibility and ease of use, CRISPR has been employed widely in many different organisms and experimental settings. As a result, thousands of publications have used CRISPR to make specific genetic perturbations, establishing in itself a resource of validated guide RNA sequences. While numerous computational tools to assist in the design and identification of candidate guide RNAs exist, these are still just at best predictions and generally, researchers inevitably will test multiple sequences for functional activity. Here, we present dbGuide (https://sgrnascorer.cancer.gov/dbguide), a database of functionally validated guide RNA sequences for CRISPR/Cas9-based knockout in human and mouse. Our database not only contains computationally determined candidate guide RNA sequences, but of even greater value, over 4000 sequences which have been functionally validated either through direct amplicon sequencing or manual curation of literature from over 1000 publications. Finally, our established framework will allow for continual addition of newly published and experimentally validated guide RNA sequences for CRISPR/Cas9-based knockout as well as incorporation of sequences from different gene editing systems, additional species and other types of site-specific functionalities such as base editing, gene activation, repression and epigenetic modification.

Protonation-Dependent Sequencing of 5-Formylcytidine in RNA.

Chemical modification of cytidine in noncoding RNAs plays a key role in regulating translation and disease. However, the distribution and dynamics of many of these modifications remain unknown due to a lack of sensitive site-specific sequencing technologies. Here, we report a protonation-dependent sequencing reaction for the detection of 5-formylcytidine (5fC) and 5-carboxycytidine (5caC) in RNA. First, we evaluate how protonation combined with electron-withdrawing substituents alters the molecular orbital energies and reduction of modified cytidine nucleosides, highlighting 5fC and 5caC as reactive species. Next, we apply this reaction to detect these modifications in synthetic oligonucleotides as well as endogenous human transfer RNA (tRNA). Finally, we demonstrate the utility of our method to characterize a patient-derived model of 5fC deficiency, where it enables facile monitoring of both pathogenic loss and exogenous rescue of NSUN3-dependent 5fC within the wobble base of human mitochondrial tRNAMet. These studies showcase the ability of protonation to enhance the reactivity and sensitive detection of 5fC in RNA and more broadly provide a molecular foundation for using optimized sequencing reactions to better understand the role of oxidized RNA cytidine residues in diseases.

Comparison of Cas9 activators in multiple species.

Several programmable transcription factors exist based on the versatile Cas9 protein, yet their relative potency and effectiveness across various cell types and species remain unexplored. Here, we compare Cas9 activator systems and examine their ability to induce robust gene expression in several human, mouse, and fly cell lines. We also explore the potential for improved activation through the combination of the most potent activator systems, and we assess the role of cooperativity in maximizing gene expression.

Unraveling CRISPR-Cas9 genome engineering parameters via a library-on-library approach.

We developed an in vivo library-on-library methodology to simultaneously assess single guide RNA (sgRNA) activity across ∼1,400 genomic loci. Assaying across multiple human cell types and end-processing enzymes as well as two Cas9 orthologs, we unraveled underlying nucleotide sequence and epigenetic parameters. Our results and software (http://crispr.med.harvard.edu/sgRNAScorer) enable improved design of reagents, shed light on mechanisms of genome targeting, and provide a generalizable framework to study nucleic acid-nucleic acid interactions and biochemistry in high throughput.

Cas9 gRNA engineering for genome editing, activation and repression.

We demonstrate that by altering the length of Cas9-associated guide RNA (gRNA) we were able to control Cas9 nuclease activity and simultaneously perform genome editing and transcriptional regulation with a single Cas9 protein. We exploited these principles to engineer mammalian synthetic circuits with combined transcriptional regulation and kill functions governed by a single multifunctional Cas9 protein.