Glucose dissociates DDX21 dimers to regulate mRNA splicing and tissue differentiation.

Glucose is a universal bioenergy source; however, its role in controlling protein interactions is unappreciated, as are its actions during differentiation-associated intracellular glucose elevation. Azido-glucose click chemistry identified glucose binding to a variety of RNA binding proteins (RBPs), including the DDX21 RNA helicase, which was found to be essential for epidermal differentiation. Glucose bound the ATP-binding domain of DDX21, altering protein conformation, inhibiting helicase activity, and dissociating DDX21 dimers. Glucose elevation during differentiation was associated with DDX21 re-localization from the nucleolus to the nucleoplasm where DDX21 assembled into larger protein complexes containing RNA splicing factors. DDX21 localized to specific SCUGSDGC motif in mRNA introns in a glucose-dependent manner and promoted the splicing of key pro-differentiation genes, including GRHL3, KLF4, OVOL1, and RBPJ. These findings uncover a biochemical mechanism of action for glucose in modulating the dimerization and function of an RNA helicase essential for tissue differentiation.

Multimodal Analysis of Composition and Spatial Architecture in Human Squamous Cell Carcinoma.

To define the cellular composition and architecture of cutaneous squamous cell carcinoma (cSCC), we combined single-cell RNA sequencing with spatial transcriptomics and multiplexed ion beam imaging from a series of human cSCCs and matched normal skin. cSCC exhibited four tumor subpopulations, three recapitulating normal epidermal states, and a tumor-specific keratinocyte (TSK) population unique to cancer, which localized to a fibrovascular niche. Integration of single-cell and spatial data mapped ligand-receptor networks to specific cell types, revealing TSK cells as a hub for intercellular communication. Multiple features of potential immunosuppression were observed, including T regulatory cell (Treg) co-localization with CD8 T cells in compartmentalized tumor stroma. Finally, single-cell characterization of human tumor xenografts and in vivo CRISPR screens identified essential roles for specific tumor subpopulation-enriched gene networks in tumorigenesis. These data define cSCC tumor and stromal cell subpopulations, the spatial niches where they interact, and the communicating gene networks that they engage in cancer.

Defining a Cancer Dependency Map.

Most human epithelial tumors harbor numerous alterations, making it difficult to predict which genes are required for tumor survival. To systematically identify cancer dependencies, we analyzed 501 genome-scale loss-of-function screens performed in diverse human cancer cell lines. We developed DEMETER, an analytical framework that segregates on- from off-target effects of RNAi. 769 genes were differentially required in subsets of these cell lines at a threshold of six SDs from the mean. We found predictive models for 426 dependencies (55%) by nonlinear regression modeling considering 66,646 molecular features. Many dependencies fall into a limited number of classes, and unexpectedly, in 82% of models, the top biomarkers were expression based. We demonstrated the basis behind one such predictive model linking hypermethylation of the UBB ubiquitin gene to a dependency on UBC. Together, these observations provide a foundation for a cancer dependency map that facilitates the prioritization of therapeutic targets.

Long Neural Genes Harbor Recurrent DNA Break Clusters in Neural Stem/Progenitor Cells.

Repair of DNA double-strand breaks (DSBs) by non-homologous end joining is critical for neural development, and brain cells frequently contain somatic genomic variations that might involve DSB intermediates. We now use an unbiased, high-throughput approach to identify genomic regions harboring recurrent DSBs in primary neural stem/progenitor cells (NSPCs). We identify 27 recurrent DSB clusters (RDCs), and remarkably, all occur within gene bodies. Most of these NSPC RDCs were detected only upon mild, aphidicolin-induced replication stress, providing a nucleotide-resolution view of replication-associated genomic fragile sites. The vast majority of RDCs occur in long, transcribed, and late-replicating genes. Moreover, almost 90% of identified RDC-containing genes are involved in synapse function and/or neural cell adhesion, with a substantial fraction also implicated in tumor suppression and/or mental disorders. Our characterization of NSPC RDCs reveals a basis of gene fragility and suggests potential impacts of DNA breaks on neurodevelopment and neural functions.

Chromosomal Loop Domains Direct the Recombination of Antigen Receptor Genes.

RAG initiates antibody V(D)J recombination in developing lymphocytes by generating "on-target" DNA breaks at matched pairs of bona fide recombination signal sequences (RSSs). We employ bait RAG-generated breaks in endogenous or ectopically inserted RSS pairs to identify huge numbers of RAG "off-target" breaks. Such breaks occur at the simple CAC motif that defines the RSS cleavage site and are largely confined within convergent CTCF-binding element (CBE)-flanked loop domains containing bait RSS pairs. Marked orientation dependence of RAG off-target activity within loops spanning up to 2 megabases implies involvement of linear tracking. In this regard, major RAG off-targets in chromosomal translocations occur as convergent RSS pairs at enhancers within a loop. Finally, deletion of a CBE-based IgH locus element disrupts V(D)J recombination domains and, correspondingly, alters RAG on- and off-target distributions within IgH. Our findings reveal how RAG activity is developmentally focused and implicate mechanisms by which chromatin domains harness biological processes within them.

Sequence-Intrinsic Mechanisms that Target AID Mutational Outcomes on Antibody Genes.

In activated B lymphocytes, AID initiates antibody variable (V) exon somatic hypermutation (SHM) for affinity maturation in germinal centers (GCs) and IgH switch (S) region DNA breaks (DSBs) for class-switch recombination (CSR). To resolve long-standing questions, we have developed an in vivo assay to study AID targeting of passenger sequences replacing a V exon. First, we find AID targets SHM hotspots within V exon and S region passengers at similar frequencies and that the normal SHM process frequently generates deletions, indicating that SHM and CSR employ the same mechanism. Second, AID mutates targets in diverse non-Ig passengers in GC B cells at levels similar to those of V exons, definitively establishing the V exon location as "privileged" for SHM. Finally, Peyer's patch GC B cells generate a reservoir of V exons that are highly mutated before selection for affinity maturation. We discuss the implications of these findings for harnessing antibody diversification mechanisms.

Convergent transcription at intragenic super-enhancers targets AID-initiated genomic instability.

Activation-induced cytidine deaminase (AID) initiates both somatic hypermutation (SHM) for antibody affinity maturation and DNA breakage for antibody class switch recombination (CSR) via transcription-dependent cytidine deamination of single-stranded DNA targets. Though largely specific for immunoglobulin genes, AID also acts on a limited set of off-targets, generating oncogenic translocations and mutations that contribute to B cell lymphoma. How AID is recruited to off-targets has been a long-standing mystery. Based on deep GRO-seq studies of mouse and human B lineage cells activated for CSR or SHM, we report that most robust AID off-target translocations occur within highly focal regions of target genes in which sense and antisense transcription converge. Moreover, we found that such AID-targeting "convergent" transcription arises from antisense transcription that emanates from super-enhancers within sense transcribed gene bodies. Our findings provide an explanation for AID off-targeting to a small subset of mostly lineage-specific genes in activated B cells.

PROBER identifies proteins associated with programmable sequence-specific DNA in living cells.

DNA-protein interactions mediate physiologic gene regulation and may be altered by DNA variants linked to polygenic disease. To enhance the speed and signal-to-noise ratio (SNR) in the identification and quantification of proteins associated with specific DNA sequences in living cells, we developed proximal biotinylation by episomal recruitment (PROBER). PROBER uses high-copy episomes to amplify SNR, and proximity proteomics (BioID) to identify the transcription factors and additional gene regulators associated with short DNA sequences of interest. PROBER quantified both constitutive and inducible association of transcription factors and corresponding chromatin regulators to target DNA sequences and binding quantitative trait loci due to single-nucleotide variants. PROBER identified alterations in regulator associations due to cancer hotspot mutations in the hTERT promoter, indicating that these mutations increase promoter association with specific gene activators. PROBER provides an approach to rapidly identify proteins associated with specific DNA sequences and their variants in living cells.

Orientation-specific joining of AID-initiated DNA breaks promotes antibody class switching.

During B-cell development, RAG endonuclease cleaves immunoglobulin heavy chain (IgH) V, D, and J gene segments and orchestrates their fusion as deletional events that assemble a V(D)J exon in the same transcriptional orientation as adjacent Cµ constant region exons. In mice, six additional sets of constant region exons (CHs) lie 100-200 kilobases downstream in the same transcriptional orientation as V(D)J and Cµ exons. Long repetitive switch (S) regions precede Cµ and downstream CHs. In mature B cells, class switch recombination (CSR) generates different antibody classes by replacing Cµ with a downstream CH (ref. 2). Activation-induced cytidine deaminase (AID) initiates CSR by promoting deamination lesions within Sµ and a downstream acceptor S region; these lesions are converted into DNA double-strand breaks (DSBs) by general DNA repair factors. Productive CSR must occur in a deletional orientation by joining the upstream end of an Sµ DSB to the downstream end of an acceptor S-region DSB. However, the relative frequency of deletional to inversional CSR junctions has not been measured. Thus, whether orientation-specific joining is a programmed mechanistic feature of CSR as it is for V(D)J recombination and, if so, how this is achieved is unknown. To address this question, we adapt high-throughput genome-wide translocation sequencing into a highly sensitive DSB end-joining assay and apply it to endogenous AID-initiated S-region DSBs in mouse B cells. We show that CSR is programmed to occur in a productive deletional orientation and does so via an unprecedented mechanism that involves in cis Igh organizational features in combination with frequent S-region DSBs initiated by AID. We further implicate ATM-dependent DSB-response factors in enforcing this mechanism and provide an explanation of why CSR is so reliant on the 53BP1 DSB-response factor.

Sequence intrinsic somatic mutation mechanisms contribute to affinity maturation of VRC01-class HIV-1 broadly neutralizing antibodies.

Variable regions of Ig chains provide the antigen recognition portion of B-cell receptors and derivative antibodies. Ig heavy-chain variable region exons are assembled developmentally from V, D, J gene segments. Each variable region contains three antigen-contacting complementarity-determining regions (CDRs), with CDR1 and CDR2 encoded by the V segment and CDR3 encoded by the V(D)J junction region. Antigen-stimulated germinal center (GC) B cells undergo somatic hypermutation (SHM) of V(D)J exons followed by selection for SHMs that increase antigen-binding affinity. Some HIV-1-infected human subjects develop broadly neutralizing antibodies (bnAbs), such as the potent VRC01-class bnAbs, that neutralize diverse HIV-1 strains. Mature VRC01-class bnAbs, including VRC-PG04, accumulate very high SHM levels, a property that hinders development of vaccine strategies to elicit them. Because many VRC01-class bnAb SHMs are not required for broad neutralization, high overall SHM may be required to achieve certain functional SHMs. To elucidate such requirements, we used a V(D)J passenger allele system to assay, in mouse GC B cells, sequence-intrinsic SHM-targeting rates of nucleotides across substrates representing maturation stages of human VRC-PG04. We identify rate-limiting SHM positions for VRC-PG04 maturation, as well as SHM hotspots and intrinsically frequent deletions associated with SHM. We find that mature VRC-PG04 has low SHM capability due to hotspot saturation but also demonstrate that generation of new SHM hotspots and saturation of existing hotspot regions (e.g., CDR3) does not majorly influence intrinsic SHM in unmutated portions of VRC-PG04 progenitor sequences. We discuss implications of our findings for bnAb affinity maturation mechanisms.

Microbial colonization influences early B-lineage development in the gut lamina propria.

The RAG1/RAG2 endonuclease (RAG) initiates the V(D)J recombination reaction that assembles immunoglobulin heavy (IgH) and light (IgL) chain variable region exons from germline gene segments to generate primary antibody repertoires. IgH V(D)J assembly occurs in progenitor (pro-) B cells followed by that of IgL in precursor (pre-) B cells. Expression of IgH µ and IgL (Igκ or Igλ) chains generates IgM, which is expressed on immature B cells as the B-cell antigen-binding receptor (BCR). Rag expression can continue in immature B cells, allowing continued Igκ V(D)J recombination that replaces the initial VκJκ exon with one that generates a new specificity. This 'receptor editing' process, which can also lead to Igλ V(D)J recombination and expression, provides a mechanism whereby antigen encounter at the Rag-expressing immature B-cell stage helps shape pre-immune BCR repertoires. As the major site of postnatal B-cell development, the bone marrow is the principal location of primary immunoglobulin repertoire diversification in mice. Here we report that early B-cell development also occurs within the mouse intestinal lamina propria (LP), where the associated V(D)J recombination/receptor editing processes modulate primary LP immunoglobulin repertoires. At weanling age in normally housed mice, the LP contains a population of Rag-expressing B-lineage cells that harbour intermediates indicative of ongoing V(D)J recombination and which contain cells with pro-B, pre-B and editing phenotypes. Consistent with LP-specific receptor editing, Rag-expressing LP B-lineage cells have similar VH repertoires, but significantly different Vκ repertoires, compared to those of Rag2-expressing bone marrow counterparts. Moreover, colonization of germ-free mice leads to an increased ratio of Igλ-expressing versus Igκ-expressing B cells specifically in the LP. We conclude that B-cell development occurs in the intestinal mucosa, where it is regulated by extracellular signals from commensal microbes that influence gut immunoglobulin repertoires.

Transcription-associated processes cause DNA double-strand breaks and translocations in neural stem/progenitor cells.

High-throughput, genome-wide translocation sequencing (HTGTS) studies of activated B cells have revealed that DNA double-strand breaks (DSBs) capable of translocating to defined bait DSBs are enriched around the transcription start sites (TSSs) of active genes. We used the HTGTS approach to investigate whether a similar phenomenon occurs in primary neural stem/progenitor cells (NSPCs). We report that breakpoint junctions indeed are enriched around TSSs that were determined to be active by global run-on sequencing analyses of NSPCs. Comparative analyses of transcription profiles in NSPCs and B cells revealed that the great majority of TSS-proximal junctions occurred in genes commonly expressed in both cell types, possibly because this common set has higher transcription levels on average than genes transcribed in only one or the other cell type. In the latter context, among all actively transcribed genes containing translocation junctions in NSPCs, those with junctions located within 2 kb of the TSS show a significantly higher transcription rate on average than genes with junctions in the gene body located at distances greater than 2 kb from the TSS. Finally, analysis of repair junction signatures of TSS-associated translocations in wild-type versus classical nonhomologous end-joining (C-NHEJ)-deficient NSPCs reveals that both C-NHEJ and alternative end-joining pathways can generate translocations by joining TSS-proximal DSBs to DSBs on other chromosomes. Our studies show that the generation of transcription-associated DSBs is conserved across divergent cell types.

Developmental propagation of V(D)J recombination-associated DNA breaks and translocations in mature B cells via dicentric chromosomes.

Mature IgM(+) B-cell lymphomas that arise in certain ataxia telangiectasia-mutated (ATM)-deficient compound mutant mice harbor translocations that fuse V(D)J recombination-initiated IgH double-strand breaks (DSBs) on chromosome 12 to sequences downstream of c-myc on chromosome 15, generating dicentric chromosomes and c-myc amplification via a breakage-fusion-bridge mechanism. As V(D)J recombination DSBs occur in developing progenitor B cells in the bone marrow, we sought to elucidate a mechanism by which such DSBs contribute to oncogenic translocations/amplifications in mature B cells. For this purpose, we applied high-throughput genome-wide translocation sequencing to study the fate of introduced c-myc DSBs in splenic IgM(+) B cells stimulated for activation-induced cytidine deaminase (AID)-dependent IgH class switch recombination (CSR). We found frequent translocations of c-myc DSBs to AID-initiated DSBs in IgH switch regions in wild-type and ATM-deficient B cells. However, c-myc also translocated frequently to newly generated DSBs within a 35-Mb region downstream of IgH in ATM-deficient, but not wild-type, CSR-activated B cells. Moreover, we found such DSBs and translocations in activated B cells that did not express AID or undergo CSR. Our findings indicate that ATM deficiency leads to formation of chromosome 12 dicentrics via recombination-activating gene-initiated IgH DSBs in progenitor B cells and that these dicentrics can be propagated developmentally into mature B cells where they generate new DSBs downstream of IgH via breakage-fusion-bridge cycles. We propose that dicentrics formed by joining V(D)J recombination-associated IgH DSBs to DSBs downstream of c-myc in ATM-deficient B lineage cells similarly contribute to c-myc amplification and mature B-cell lymphomas.

IgH class switching exploits a general property of two DNA breaks to be joined in cis over long chromosomal distances.

Antibody class switch recombination (CSR) in B lymphocytes joins two DNA double-strand breaks (DSBs) lying 100-200 kb apart within switch (S) regions in the immunoglobulin heavy-chain locus (IgH). CSR-activated B lymphocytes generate multiple S-region DSBs in the donor Sµ and in a downstream acceptor S region, with a DSB in Sµ being joined to a DSB in the acceptor S region at sufficient frequency to drive CSR in a large fraction of activated B cells. Such frequent joining of widely separated CSR DSBs could be promoted by IgH-specific or B-cell-specific processes or by general aspects of chromosome architecture and DSB repair. Previously, we found that B cells with two yeast I-SceI endonuclease targets in place of Sγ1 undergo I-SceI-dependent class switching from IgM to IgG1 at 5-10% of normal levels. Now, we report that B cells in which Sγ1 is replaced with a 28 I-SceI target array, designed to increase I-SceI DSB frequency, undergo I-SceI-dependent class switching at almost normal levels. High-throughput genome-wide translocation sequencing revealed that I-SceI-generated DSBs introduced in cis at Sµ and Sγ1 sites are joined together in T cells at levels similar to those of B cells. Such high joining levels also occurred between I-SceI-generated DSBs within c-myc and I-SceI- or CRISPR/Cas9-generated DSBs 100 kb downstream within Pvt1 in B cells or fibroblasts, respectively. We suggest that CSR exploits a general propensity of intrachromosomal DSBs separated by several hundred kilobases to be frequently joined together and discuss the relevance of this finding for recurrent interstitial deletions in cancer.

Analyzing RNA-Protein Interactions by Cross-Link Rates and CLIP-seq Libraries.

UV cross-linking-based methods are the most common tool to explore in vivo RNA-protein interactions. UV cross-linking enables the freezing of direct interactions in the cell, which can then be mapped by high-throughput sequencing through a family of methods termed CLIP-seq. CLIP-seq measures the distribution of cross-link events by purifying a protein of interest and sequencing the covalently bound RNA fragments. However, there are disagreements and ambiguities as to which proteins are RNA-binding proteins and what interactions are significant as all proteins contact all RNAs at some frequency. Here we describe a protocol for both determining RNA-protein interactions through a combination of RNA library preparation and the measurement of absolute cross-link rates, which helps determine what proteins are RNA-binding proteins and what interactions are significant. This protocol, comprising an updated form of the easyCLIP protocol, describes guidelines for RNA library preparation, oligo and protein standard construction, and the measurement of cross-link rates. These methods are easily visualizable through their fluorescent labels and can be adapted to study RNA-binding properties of both functional, high affinity RNA-binding proteins, and the accidental RNA interactions of non-RNA-binding proteins. © 2023 Wiley Periodicals LLC. Basic Protocol 1: RNA library construction Basic Protocol 2: Determining UV cross-link rates Support Protocol 1: Cross-linking and lysing cells Support Protocol 2: Adapter preparation Support Protocol 3: Preparation of cross-linked RBP standard.

Genome-scale functional genomics identify genes preferentially essential for multiple myeloma cells compared to other neoplasias.

Clinical progress in multiple myeloma (MM), an incurable plasma cell (PC) neoplasia, has been driven by therapies that have limited applications beyond MM/PC neoplasias and do not target specific oncogenic mutations in MM. Instead, these agents target pathways critical for PC biology yet largely dispensable for malignant or normal cells of most other lineages. Here we systematically characterized the lineage-preferential molecular dependencies of MM through genome-scale clustered regularly interspaced short palindromic repeats (CRISPR) studies in 19 MM versus hundreds of non-MM lines and identified 116 genes whose disruption more significantly affects MM cell fitness compared with other malignancies. These genes, some known, others not previously linked to MM, encode transcription factors, chromatin modifiers, endoplasmic reticulum components, metabolic regulators or signaling molecules. Most of these genes are not among the top amplified, overexpressed or mutated in MM. Functional genomics approaches thus define new therapeutic targets in MM not readily identifiable by standard genomic, transcriptional or epigenetic profiling analyses.

A cis-regulatory lexicon of DNA motif combinations mediating cell-type-specific gene regulation.

Gene expression is controlled by transcription factors (TFs) that bind cognate DNA motif sequences in cis-regulatory elements (CREs). The combinations of DNA motifs acting within homeostasis and disease, however, are unclear. Gene expression, chromatin accessibility, TF footprinting, and H3K27ac-dependent DNA looping data were generated and a random-forest-based model was applied to identify 7,531 cell-type-specific cis-regulatory modules (CRMs) across 15 diploid human cell types. A co-enrichment framework within CRMs nominated 838 cell-type-specific, recurrent heterotypic DNA motif combinations (DMCs), which were functionally validated using massively parallel reporter assays. Cancer cells engaged DMCs linked to neoplasia-enabling processes operative in normal cells while also activating new DMCs only seen in the neoplastic state. This integrative approach identifies cell-type-specific cis-regulatory combinatorial DNA motifs in diverse normal and diseased human cells and represents a general framework for deciphering cis-regulatory sequence logic in gene regulation.

A genome-wide atlas of co-essential modules assigns function to uncharacterized genes.

A central question in the post-genomic era is how genes interact to form biological pathways. Measurements of gene dependency across hundreds of cell lines have been used to cluster genes into 'co-essential' pathways, but this approach has been limited by ubiquitous false positives. In the present study, we develop a statistical method that enables robust identification of gene co-essentiality and yields a genome-wide set of functional modules. This atlas recapitulates diverse pathways and protein complexes, and predicts the functions of 108 uncharacterized genes. Validating top predictions, we show that TMEM189 encodes plasmanylethanolamine desaturase, a key enzyme for plasmalogen synthesis. We also show that C15orf57 encodes a protein that binds the AP2 complex, localizes to clathrin-coated pits and enables efficient transferrin uptake. Finally, we provide an interactive webtool for the community to explore our results, which establish co-essentiality profiling as a powerful resource for biological pathway identification and discovery of new gene functions.

Small-Molecule and CRISPR Screening Converge to Reveal Receptor Tyrosine Kinase Dependencies in Pediatric Rhabdoid Tumors.

Cancer is often seen as a disease of mutations and chromosomal abnormalities. However, some cancers, including pediatric rhabdoid tumors (RTs), lack recurrent alterations targetable by current drugs and need alternative, informed therapeutic options. To nominate potential targets, we performed a high-throughput small-molecule screen complemented by a genome-scale CRISPR-Cas9 gene-knockout screen in a large number of RT and control cell lines. These approaches converged to reveal several receptor tyrosine kinases (RTKs) as therapeutic targets, with RTK inhibition effective in suppressing RT cell growth in vitro and against a xenograft model in vivo. RT cell lines highly express and activate (phosphorylate) different RTKs, creating dependency without mutation or amplification. Downstream of RTK signaling, we identified PTPN11, encoding the pro-growth signaling protein SHP2, as a shared dependency across all RT cell lines. This study demonstrates that large-scale perturbational screening can uncover vulnerabilities in cancers with "quiet" genomes.