Long-term retention of antigens in germinal centers is controlled by the spatial organization of the follicular dendritic cell network.

Germinal centers (GCs) require sustained availability of antigens to promote antibody affinity maturation against pathogens and vaccines. A key source of antigens for GC B cells are immune complexes (ICs) displayed on follicular dendritic cells (FDCs). Here we show that FDC spatial organization regulates antigen dynamics in the GC. We identify heterogeneity within the FDC network. While the entire light zone (LZ) FDC network captures ICs initially, only the central cells of the network function as the antigen reservoir, where different antigens arriving from subsequent immunizations colocalize. Mechanistically, central LZ FDCs constitutively express subtly higher CR2 membrane densities than peripheral LZ FDCs, which strongly increases the IC retention half-life. Even though repeated immunizations gradually saturate central FDCs, B cell responses remain efficient because new antigens partially displace old ones. These results reveal the principles shaping antigen display on FDCs during the GC reaction.

Regulation of the RNAPII Pool Is Integral to the DNA Damage Response.

In response to transcription-blocking DNA damage, cells orchestrate a multi-pronged reaction, involving transcription-coupled DNA repair, degradation of RNA polymerase II (RNAPII), and genome-wide transcription shutdown. Here, we provide insight into how these responses are connected by the finding that ubiquitylation of RNAPII itself, at a single lysine (RPB1 K1268), is the focal point for DNA-damage-response coordination. K1268 ubiquitylation affects DNA repair and signals RNAPII degradation, essential for surviving genotoxic insult. RNAPII degradation results in a shutdown of transcriptional initiation, in the absence of which cells display dramatic transcriptome alterations. Additionally, regulation of RNAPII stability is central to transcription recovery-persistent RNAPII depletion underlies the failure of this process in Cockayne syndrome B cells. These data expose regulation of global RNAPII levels as integral to the cellular DNA-damage response and open the intriguing possibility that RNAPII pool size generally affects cell-specific transcription programs in genome instability disorders and even normal cells.

Deterministic Evolutionary Trajectories Influence Primary Tumor Growth: TRACERx Renal.

The evolutionary features of clear-cell renal cell carcinoma (ccRCC) have not been systematically studied to date. We analyzed 1,206 primary tumor regions from 101 patients recruited into the multi-center prospective study, TRACERx Renal. We observe up to 30 driver events per tumor and show that subclonal diversification is associated with known prognostic parameters. By resolving the patterns of driver event ordering, co-occurrence, and mutual exclusivity at clone level, we show the deterministic nature of clonal evolution. ccRCC can be grouped into seven evolutionary subtypes, ranging from tumors characterized by early fixation of multiple mutational and copy number drivers and rapid metastases to highly branched tumors with >10 subclonal drivers and extensive parallel evolution associated with attenuated progression. We identify genetic diversity and chromosomal complexity as determinants of patient outcome. Our insights reconcile the variable clinical behavior of ccRCC and suggest evolutionary potential as a biomarker for both intervention and surveillance.

UV Irradiation Induces a Non-coding RNA that Functionally Opposes the Protein Encoded by the Same Gene.

The transcription-related DNA damage response was analyzed on a genome-wide scale with great spatial and temporal resolution. Upon UV irradiation, a slowdown of transcript elongation and restriction of gene activity to the promoter-proximal ∼25 kb is observed. This is associated with a shift from expression of long mRNAs to shorter isoforms, incorporating alternative last exons (ALEs) that are more proximal to the transcription start site. Notably, this includes a shift from a protein-coding ASCC3 mRNA to a shorter ALE isoform of which the RNA, rather than an encoded protein, is critical for the eventual recovery of transcription. The non-coding ASCC3 isoform counteracts the function of the protein-coding isoform, indicating crosstalk between them. Thus, the ASCC3 gene expresses both coding and non-coding transcript isoforms with opposite effects on transcription recovery after UV-induced DNA damage.

Branching topology of the human embryo transcriptome revealed by Entropy Sort Feature Weighting.

Analysis of single cell transcriptomics (scRNA-seq) data is typically performed after subsetting to highly variable genes (HVGs). Here, we show that Entropy Sorting provides an alternative mathematical framework for feature selection. On synthetic datasets, continuous Entropy Sort Feature Weighting (cESFW) outperforms HVG selection in distinguishing cell-state-specific genes. We apply cESFW to six merged scRNA-seq datasets spanning human early embryo development. Without smoothing or augmenting the raw counts matrices, cESFW generates a high-resolution embedding displaying coherent developmental progression from eight-cell to post-implantation stages and delineating 15 distinct cell states. The embedding highlights sequential lineage decisions during blastocyst development, while unsupervised clustering identifies branch point populations obscured in previous analyses. The first branching region, where morula cells become specified for inner cell mass or trophectoderm, includes cells previously asserted to lack a developmental trajectory. We quantify the relatedness of different pluripotent stem cell cultures to distinct embryo cell types and identify marker genes of naïve and primed pluripotency. Finally, by revealing genes with dynamic lineage-specific expression, we provide markers for staging progression from morula to blastocyst.

Regulation of intestinal immunity and tissue repair by enteric glia.

Tissue maintenance and repair depend on the integrated activity of multiple cell types1. Whereas the contributions of epithelial2,3, immune4,5 and stromal cells6,7 in intestinal tissue integrity are well understood, the role of intrinsic neuroglia networks remains largely unknown. Here we uncover important roles of enteric glial cells (EGCs) in intestinal homeostasis, immunity and tissue repair. We demonstrate that infection of mice with Heligmosomoides polygyrus leads to enteric gliosis and the upregulation of an interferon gamma (IFNγ) gene signature. IFNγ-dependent gene modules were also induced in EGCs from patients with inflammatory bowel disease8. Single-cell transcriptomics analysis of the tunica muscularis showed that glia-specific abrogation of IFNγ signalling leads to tissue-wide activation of pro-inflammatory transcriptional programs. Furthermore, disruption of the IFNγ-EGC signalling axis enhanced the inflammatory and granulomatous response of the tunica muscularis to helminths. Mechanistically, we show that the upregulation of Cxcl10 is an early immediate response of EGCs to IFNγ signalling and provide evidence that this chemokine and the downstream amplification of IFNγ signalling in the tunica muscularis are required for a measured inflammatory response to helminths and resolution of the granulomatous pathology. Our study demonstrates that IFNγ signalling in enteric glia is central to intestinal homeostasis and reveals critical roles of the IFNγ-EGC-CXCL10 axis in immune response and tissue repair after infectious challenge.

Neuronal programming by microbiota regulates intestinal physiology.

Neural control of the function of visceral organs is essential for homeostasis and health. Intestinal peristalsis is critical for digestive physiology and host defence, and is often dysregulated in gastrointestinal disorders1. Luminal factors, such as diet and microbiota, regulate neurogenic programs of gut motility2-5, but the underlying molecular mechanisms remain unclear. Here we show that the transcription factor aryl hydrocarbon receptor (AHR) functions as a biosensor in intestinal neural circuits, linking their functional output to the microbial environment of the gut lumen. Using nuclear RNA sequencing of mouse enteric neurons that represent distinct intestinal segments and microbiota states, we demonstrate that the intrinsic neural networks of the colon exhibit unique transcriptional profiles that are controlled by the combined effects of host genetic programs and microbial colonization. Microbiota-induced expression of AHR in neurons of the distal gastrointestinal tract enables these neurons to respond to the luminal environment and to induce expression of neuron-specific effector mechanisms. Neuron-specific deletion of Ahr, or constitutive overexpression of its negative feedback regulator CYP1A1, results in reduced peristaltic activity of the colon, similar to that observed in microbiota-depleted mice. Finally, expression of Ahr in the enteric neurons of mice treated with antibiotics partially restores intestinal motility. Together, our experiments identify AHR signalling in enteric neurons as a regulatory node that integrates the luminal environment with the physiological output of intestinal neural circuits to maintain gut homeostasis and health.

Phosphoproteomic identification of ULK substrates reveals VPS15-dependent ULK/VPS34 interplay in the regulation of autophagy.

Autophagy is a process through which intracellular cargoes are catabolised inside lysosomes. It involves the formation of autophagosomes initiated by the serine/threonine kinase ULK and class III PI3 kinase VPS34 complexes. Here, unbiased phosphoproteomics screens in mouse embryonic fibroblasts deleted for Ulk1/2 reveal that ULK loss significantly alters the phosphoproteome, with novel high confidence substrates identified including VPS34 complex member VPS15 and AMPK complex subunit PRKAG2. We identify six ULK-dependent phosphorylation sites on VPS15, mutation of which reduces autophagosome formation in cells and VPS34 activity in vitro. Mutation of serine 861, the major VPS15 phosphosite, decreases both autophagy initiation and autophagic flux. Analysis of VPS15 knockout cells reveals two novel ULK-dependent phenotypes downstream of VPS15 removal that can be partially recapitulated by chronic VPS34 inhibition, starvation-independent accumulation of ULK substrates and kinase activity-regulated recruitment of autophagy proteins to ubiquitin-positive structures.

KLF17 promotes human naïve pluripotency but is not required for its establishment.

Current knowledge of the transcriptional regulation of human pluripotency is incomplete, with lack of interspecies conservation observed. Single-cell transcriptomics analysis of human embryos previously enabled us to identify transcription factors, including the zinc-finger protein KLF17, that are enriched in the human epiblast and naïve human embryonic stem cells (hESCs). Here, we show that KLF17 is expressed coincident with the known pluripotency-associated factors NANOG and SOX2 across human blastocyst development. We investigate the function of KLF17 using primed and naïve hESCs for gain- and loss-of-function analyses. We find that ectopic expression of KLF17 in primed hESCs is sufficient to induce a naïve-like transcriptome and that KLF17 can drive transgene-mediated resetting to naïve pluripotency. This implies a role for KLF17 in establishing naïve pluripotency. However, CRISPR-Cas9-mediated knockout studies reveal that KLF17 is not required for naïve pluripotency acquisition in vitro. Transcriptome analysis of naïve hESCs identifies subtle effects on metabolism and signalling pathways following KLF17 loss of function, and possible redundancy with other KLF paralogues. Overall, we show that KLF17 is sufficient, but not necessary, for naïve pluripotency under the given in vitro conditions.

Mutation of cancer driver MLL2 results in transcription stress and genome instability.

Genome instability is a recurring feature of tumorigenesis. Mutation in MLL2, encoding a histone methyltransferase, is a driver in numerous different cancer types, but the mechanism is unclear. Here, we present evidence that MLL2 mutation results in genome instability. Mouse cells in which MLL2 gene deletion can be induced display elevated levels of sister chromatid exchange, gross chromosomal aberrations, 53BP1 foci, and micronuclei. Human MLL2 knockout cells are characterized by genome instability as well. Interestingly, MLL2 interacts with RNA polymerase II (RNAPII) and RECQL5, and, although MLL2 mutated cells have normal overall H3K4me levels in genes, nucleosomes in the immediate vicinity of RNAPII are hypomethylated. Importantly, MLL2 mutated cells display signs of substantial transcription stress, and the most affected genes overlap with early replicating fragile sites, show elevated levels of γH2AX, and suffer frequent mutation. The requirement for MLL2 in the maintenance of genome stability in genes helps explain its widespread role in cancer and points to transcription stress as a strong driver in tumorigenesis.

Control of skeletal morphogenesis by the Hippo-YAP/TAZ pathway.

The Hippo-YAP/TAZ pathway is an important regulator of tissue growth, but can also control cell fate or tissue morphogenesis. Here, we investigate the function of the Hippo pathway during the development of cartilage, which forms the majority of the skeleton. Previously, YAP was proposed to inhibit skeletal size by repressing chondrocyte proliferation and differentiation. We find that, in vitro, Yap/Taz double knockout impairs murine chondrocyte proliferation, whereas constitutively nuclear nls-YAP5SA accelerates proliferation, in line with the canonical role of this pathway in most tissues. However, in vivo, cartilage-specific knockout of Yap/Taz does not prevent chondrocyte proliferation, differentiation or skeletal growth, but rather results in various skeletal deformities including cleft palate. Cartilage-specific expression of nls-YAP5SA or knockout of Lats1/2 do not increase cartilage growth, but instead lead to catastrophic malformations resembling chondrodysplasia or achondrogenesis. Physiological YAP target genes in cartilage include Ctgf, Cyr61 and several matrix remodelling enzymes. Thus, YAP/TAZ activity controls chondrocyte proliferation in vitro, possibly reflecting a regenerative response, but is dispensable for chondrocyte proliferation in vivo, and instead functions to control cartilage morphogenesis via regulation of the extracellular matrix.

Extracellular matrix anisotropy is determined by TFAP2C-dependent regulation of cell collisions.

The isotropic or anisotropic organization of biological extracellular matrices has important consequences for tissue function. We study emergent anisotropy using fibroblasts that generate varying degrees of matrix alignment from uniform starting conditions. This reveals that the early migratory paths of fibroblasts are correlated with subsequent matrix organization. Combined experimentation and adaptation of Vicsek modelling demonstrates that the reorientation of cells relative to each other following collision plays a role in generating matrix anisotropy. We term this behaviour 'cell collision guidance'. The transcription factor TFAP2C regulates cell collision guidance in part by controlling the expression of RND3. RND3 localizes to cell-cell collision zones where it downregulates actomyosin activity. Cell collision guidance fails without this mechanism in place, leading to isotropic matrix generation. The cross-referencing of alignment and TFAP2C gene expression signatures against existing datasets enables the identification and validation of several classes of pharmacological agents that disrupt matrix anisotropy.

Analysis of RNA polymerase II ubiquitylation and proteasomal degradation.

Transcribing RNA polymerase II (RNAPII) is decorated by a plethora of post-translational modifications that mark different stages of transcription. One important modification is RNAPII ubiquitylation, which occurs in response to numerous different stimuli that cause RNAPII stalling, such as DNA damaging agents, RNAPII inhibitors, or depletion of the nucleotide pool. Stalled RNAPII triggers a so-called "last resort pathway", which involves RNAPII poly-ubiquitylation and proteasome-mediated degradation. Different approaches have been described to study RNAPII poly-ubiquitylation and degradation, each method with its own advantages and caveats. Here, we describe optimised strategies for detecting ubiquitylated RNAPII and studying its degradation, but these protocols are suitable for studying other ubiquitylated proteins as well.

Integrin signalling regulates YAP and TAZ to control skin homeostasis.

The skin is a squamous epithelium that is continuously renewed by a population of basal layer stem/progenitor cells and can heal wounds. Here, we show that the transcription regulators YAP and TAZ localise to the nucleus in the basal layer of skin and are elevated upon wound healing. Skin-specific deletion of both YAP and TAZ in adult mice slows proliferation of basal layer cells, leads to hair loss and impairs regeneration after wounding. Contact with the basal extracellular matrix and consequent integrin-Src signalling is a key determinant of the nuclear localisation of YAP/TAZ in basal layer cells and in skin tumours. Contact with the basement membrane is lost in differentiating daughter cells, where YAP and TAZ become mostly cytoplasmic. In other types of squamous epithelia and squamous cell carcinomas, a similar control mechanism is present. By contrast, columnar epithelia differentiate an apical domain that recruits CRB3, Merlin (also known as NF2), KIBRA (also known as WWC1) and SAV1 to induce Hippo signalling and retain YAP/TAZ in the cytoplasm despite contact with the basal layer extracellular matrix. When columnar epithelial tumours lose their apical domain and become invasive, YAP/TAZ becomes nuclear and tumour growth becomes sensitive to the Src inhibitor Dasatinib.

Cockayne syndrome B protein regulates recruitment of the Elongin A ubiquitin ligase to sites of DNA damage.

Elongin A performs dual functions as the transcriptionally active subunit of RNA polymerase II (Pol II) elongation factor Elongin and as the substrate recognition subunit of a Cullin-RING E3 ubiquitin ligase that ubiquitylates Pol II in response to DNA damage. Assembly of the Elongin A ubiquitin ligase and its recruitment to sites of DNA damage is a tightly regulated process induced by DNA-damaging agents and α-amanitin, a drug that induces Pol II stalling. In this study, we demonstrate (i) that Elongin A and the ubiquitin ligase subunit CUL5 associate in cells with the Cockayne syndrome B (CSB) protein and (ii) that this interaction is also induced by DNA-damaging agents and α-amanitin. In addition, we present evidence that the CSB protein promotes stable recruitment of the Elongin A ubiquitin ligase to sites of DNA damage. Our findings are consistent with the model that the Elongin A ubiquitin ligase and the CSB protein function together in a common pathway in response to Pol II stalling and DNA damage.

A ubiquitylation site in Cockayne syndrome B required for repair of oxidative DNA damage, but not for transcription-coupled nucleotide excision repair.

Cockayne syndrome B (CSB), best known for its role in transcription-coupled nucleotide excision repair (TC-NER), contains a ubiquitin-binding domain (UBD), but the functional connection between protein ubiquitylation and this UBD remains unclear. Here, we show that CSB is regulated via site-specific ubiquitylation. Mass spectrometry analysis of CSB identified lysine (K) 991 as a ubiquitylation site. Intriguingly, mutation of this residue (K991R) does not affect CSB's catalytic activity or protein stability, but greatly affects genome stability, even in the absence of induced DNA damage. Moreover, cells expressing CSB K991R are sensitive to oxidative DNA damage, but proficient for TC-NER. K991 becomes ubiquitylated upon oxidative DNA damage, and while CSB K991R is recruited normally to such damage, it fails to dissociate in a timely manner, suggesting a requirement for K991 ubiquitylation in CSB activation. Interestingly, deletion of CSB's UBD gives rise to oxidative damage sensitivity as well, while CSB ΔUBD and CSB K991R affects expression of overlapping groups of genes, further indicating a functional connection. Together, these results shed new light on the regulation of CSB, with K991R representing an important separation-of-function-mutation in this multi-functional protein.

TLR and TNF-R1 activation of the MKK3/MKK6-p38α axis in macrophages is mediated by TPL-2 kinase.

Previous studies suggested that Toll-like receptor (TLR) stimulation of the p38α MAP kinase (MAPK) is mediated by transforming growth factor-ß-activated kinase 1 (TAK1) activation of MAPK kinases, MKK3, MKK4 and MKK6. We used quantitative mass spectrometry to monitor tumour progression locus 2 (TPL-2)-dependent protein phosphorylation following TLR4 stimulation with lipopolysaccharide, comparing macrophages from wild-type mice and Map3k8(D270A/D270A) mice expressing catalytically inactive TPL-2 (MAP3K8). In addition to the established TPL-2 substrates MKK1/2, TPL-2 kinase activity was required to phosphorylate the activation loops of MKK3/6, but not of MKK4. MKK3/6 activation required IκB kinase (IKK) phosphorylation of the TPL-2 binding partner nuclear factor κ-light-chain-enhancer of activated B cells (NF-κB1) p105, similar to MKK1/2 activation. Tumour necrosis factor (TNF) stimulation of MKK3/6 phosphorylation was similarly dependent on TPL-2 catalytic activity and IKK phosphorylation of NF-κB1 p105. Owing to redundancy of MKK3/6 with MKK4, Map3k8(D270A) mutation only fractionally decreased lipopolysaccharide activation of p38α. TNF activation of p38α, which is mediated predominantly via MKK3/6, was substantially reduced. TPL-2 catalytic activity was also required for MKK3/6 and p38α activation following macrophage stimulation with Mycobacterium tuberculosis and Listeria monocytogenes Our experiments demonstrate that the IKK/NF-κB1 p105/TPL-2 signalling pathway, downstream of TAK1, regulates MKK3/6 and p38α activation in macrophages in inflammation.

Systematic diet composition swap in a mouse genome-scale metabolic model reveals determinants of obesogenic diet metabolism in liver cancer.

Dietary nutrient availability and gene expression, together, influence tissue metabolic activity. Here, we explore whether altering dietary nutrient composition in the context of mouse liver cancer suffices to overcome chronic gene expression changes that arise from tumorigenesis and western-style diet (WD). We construct a mouse genome-scale metabolic model and estimate metabolic fluxes in liver tumors and non-tumoral tissue after computationally varying the composition of input diet. This approach, called Systematic Diet Composition Swap (SyDiCoS), revealed that, compared to a control diet, WD increases production of glycerol and succinate irrespective of specific tissue gene expression patterns. Conversely, differences in fatty acid utilization pathways between tumor and non-tumor liver are amplified with WD by both dietary carbohydrates and lipids together. Our data suggest that combined dietary component modifications may be required to normalize the distinctive metabolic patterns that underlie selective targeting of tumor metabolism.

A branching model of lineage differentiation underpinning the neurogenic potential of enteric glia.

Glial cells have been proposed as a source of neural progenitors, but the mechanisms underpinning the neurogenic potential of adult glia are not known. Using single cell transcriptomic profiling, we show that enteric glial cells represent a cell state attained by autonomic neural crest cells as they transition along a linear differentiation trajectory that allows them to retain neurogenic potential while acquiring mature glial functions. Key neurogenic loci in early enteric nervous system progenitors remain in open chromatin configuration in mature enteric glia, thus facilitating neuronal differentiation under appropriate conditions. Molecular profiling and gene targeting of enteric glial cells in a cell culture model of enteric neurogenesis and a gut injury model demonstrate that neuronal differentiation of glia is driven by transcriptional programs employed in vivo by early progenitors. Our work provides mechanistic insight into the regulatory landscape underpinning the development of intestinal neural circuits and generates a platform for advancing glial cells as therapeutic agents for the treatment of neural deficits.