Comparative analysis of 10X Chromium vs. BD Rhapsody whole transcriptome single-cell sequencing technologies in complex human tissues.

The development of single-cell omics tools has enabled scientists to study the tumor microenvironment (TME) in unprecedented detail. However, each of the different techniques may have its unique strengths and limitations. Here we directly compared two commercially available high-throughput single-cell RNA sequencing (scRNA-seq) technologies - droplet-based 10X Chromium vs. microwell-based BD Rhapsody - using paired samples from patients with localized prostate cancer (PCa) undergoing a radical prostatectomy. Although high technical consistency was observed in unraveling the whole transcriptome, the relative abundance of cell populations differed. Cells with low mRNA content such as T cells were underrepresented in the droplet-based system, at least partly due to lower RNA capture rates. In contrast, microwell-based scRNA-seq recovered less cells of epithelial origin. Moreover, we discovered platform-dependent variabilities in mRNA quantification and cell-type marker annotation. Overall, our study provides important information for selection of the appropriate scRNA-seq platform and for the interpretation of published results.

Functional and spatial proteomics profiling reveals intra- and intercellular signaling crosstalk in colorectal cancer.

Precision oncology approaches for patients with colorectal cancer (CRC) continue to lag behind other solid cancers. Functional precision oncology-a strategy that is based on perturbing primary tumor cells from cancer patients-could provide a road forward to personalize treatment. We extend this paradigm to measuring proteome activity landscapes by acquiring quantitative phosphoproteomic data from patient-derived organoids (PDOs). We show that kinase inhibitors induce inhibitor- and patient-specific off-target effects and pathway crosstalk. Reconstruction of the kinase networks revealed that the signaling rewiring is modestly affected by mutations. We show non-genetic heterogeneity of the PDOs and upregulation of stemness and differentiation genes by kinase inhibitors. Using imaging mass-cytometry-based profiling of the primary tumors, we characterize the tumor microenvironment (TME) and determine spatial heterocellular crosstalk and tumor-immune cell interactions. Collectively, we provide a framework for inferring tumor cell intrinsic signaling and external signaling from the TME to inform precision (immuno-) oncology in CRC.

13C tracer analysis reveals the landscape of metabolic checkpoints in human CD8+ T cell differentiation and exhaustion.

Introduction: Naïve T cells remain in an actively maintained state of quiescence until activation by antigenic signals, upon which they start to proliferate and generate effector cells to initiate a functional immune response. Metabolic reprogramming is essential to meet the biosynthetic demands of the differentiation process, and failure to do so can promote the development of hypofunctional exhausted T cells. Methods: Here we used 13C metabolomics and transcriptomics to study the metabolism of CD8+ T cells in their complete course of differentiation from naïve over stem-like memory to effector cells and in exhaustion-inducing conditions. Results: The quiescence of naïve T cells was evident in a profound suppression of glucose oxidation and a decreased expression of ENO1, downstream of which no glycolytic flux was detectable. Moreover, TCA cycle activity was low in naïve T cells and associated with a downregulation of SDH subunits. Upon stimulation and exit from quiescence, the initiation of cell growth and proliferation was accompanied by differential expression of metabolic enzymes and metabolic reprogramming towards aerobic glycolysis with high rates of nutrient uptake, respiration and lactate production. High flux in anabolic pathways imposed a strain on NADH homeostasis, which coincided with engagement of the proline cycle for mitochondrial redox shuttling. With acquisition of effector functions, cells increasingly relied on glycolysis as opposed to oxidative phosphorylation, which was, however, not linked to changes in mitochondrial abundance. In exhaustion, decreased effector function concurred with a reduction in mitochondrial metabolism, glycolysis and amino acid import, and an upregulation of quiescence-associated genes, TXNIP and KLF2, and the T cell suppressive metabolites succinate and itaconate. Discussion: Overall, these results identify multiple metabolic features that regulate quiescence, proliferation and effector function, but also exhaustion of CD8+ T cells during differentiation. Thus, targeting these metabolic checkpoints may be a promising therapeutic strategy for both prevention of exhaustion and promotion of stemness of anti-tumor T cells.

SARS-CoV-2 activates the TLR4/MyD88 pathway in human macrophages: A possible correlation with strong pro-inflammatory responses in severe COVID-19.

Background: Toll-like receptors (TLRs) play a pivotal role in the immunologic response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Exaggerated inflammatory response of innate immune cells, however, may drive morbidity and death in Coronavirus disease 19 (COVID-19). Objective: We investigated the engagement of SARS-CoV-2 with TLR4 in order to better understand how to tackle hyperinflammation in COVID-19. Methods: We combined RNA-sequencing data of human lung tissue and of bronchoalveolar lavage fluid cells derived from COVID-19 patients with functional studies in human macrophages using SARS-CoV-2 spike proteins and viable SARS-CoV-2. Pharmacological inhibitors as well as gene editing with CRISPR/Cas9 were used to delineate the signalling pathways involved. Results: We found TLR4 to be the most abundantly upregulated TLR in human lung tissue irrespective of the underlying pathology. Accordingly, bronchoalveolar lavage fluid cells from patients with severe COVID-19 showed an NF-κB-pathway dominated immune response, whereas they were mostly defined by type I interferon signalling in moderate COVID-19. Mechanistically, we found the Spike ectodomain, but not receptor binding domain monomer to induce TLR4-dependent inflammation in human macrophages. By using pharmacological inhibitors as well as CRISPR/Cas9 deleted macrophages, we identify SARS-CoV-2 to engage canonical TLR4-MyD88 signalling. Importantly, we demonstrate that TLR4 blockage prevents exaggerated inflammatory responses in human macrophages infected with different SARS-CoV-2 variants, including immune escape variants B.1.1.7.-E484K and B.1.1.529 (omicron). Conclusion: Our study critically extends the current knowledge on TLR-mediated hyperinflammatory responses to SARS-CoV-2 in human macrophages, paving the way for novel approaches to tackle severe COVID-19. Take-home message: Our study combining human lung transcriptomics with functional studies in human macrophages clearly supports the design and development of TLR4 - directed therapeutics to mitigate hyperinflammation in severe COVID-19.

Skeletal muscle transcriptomics dissects the pathogenesis of Friedreich's ataxia.

OBJECTIVE: In Friedreich's ataxia (FRDA), the most affected tissues are not accessible to sampling and available transcriptomic findings originate from blood-derived cells and animal models. Herein, we aimed at dissecting for the first time the pathophysiology of FRDA by means of RNA-sequencing in an affected tissue sampled in vivo. METHODS: Skeletal muscle biopsies were collected from seven FRDA patients before and after treatment with recombinant human Erythropoietin (rhuEPO) within a clinical trial. Total RNA extraction, 3'-mRNA library preparation and sequencing were performed according to standard procedures. We tested for differential gene expression with DESeq2 and performed gene set enrichment analysis with respect to control subjects. RESULTS: FRDA transcriptomes showed 1873 genes differentially expressed from controls. Two main signatures emerged: (1) a global downregulation of the mitochondrial transcriptome as well as of ribosome/translational machinery and (2) an upregulation of genes related to transcription and chromatin regulation, especially of repressor terms. Downregulation of the mitochondrial transcriptome was more profound than previously shown in other cellular systems. Furthermore, we observed in FRDA patients a marked upregulation of leptin, the master regulator of energy homeostasis. RhuEPO treatment further enhanced leptin expression. INTERPRETATION: Our findings reflect a double hit in the pathophysiology of FRDA: a transcriptional/translational issue and a profound mitochondrial failure downstream. Leptin upregulation in the skeletal muscle in FRDA may represent a compensatory mechanism of mitochondrial dysfunction, which is amenable to pharmacological boosting. Skeletal muscle transcriptomics is a valuable biomarker to monitor therapeutic interventions in FRDA.

Exosomal mitochondrial tRNAs and miRNAs as potential predictors of inflammation in renal proximal tubular epithelial cells.

Exosomes have emerged as a valuable repository of novel biomarkers for human diseases such as chronic kidney disease (CKD). From a healthy control group, we performed microRNA (miRNA) profiling of urinary exosomes and compared it with a cell culture model of renal proximal tubular epithelial cells (RPTECs). Thereby, a large fraction of abundant urinary exosomal miRNAs could also be detected in exosomes derived from RPTECs, indicating them as a suitable model system for investigation of CKD. We subsequently analyzed exosomes from RPTECs in pro-inflammatory and pro-fibrotic states, mimicking some aspects of CKD. Following cytokine treatment, we observed a significant increase in exosome release and identified 30 dysregulated exosomal miRNAs, predominantly associated with the regulation of pro-inflammatory and pro-fibrotic-related pathways. In addition to miRNAs, we also identified 16 dysregulated exosomal mitochondrial RNAs, highlighting a pivotal role of mitochondria in sensing renal inflammation. Inhibitors of exosome biogenesis and release significantly altered the abundance of selected candidate miRNAs and mitochondrial RNAs, thus suggesting distinct sorting mechanisms of different non-coding RNA (ncRNA) species into exosomes. Hence, these two exosomal ncRNA species might be employed as potential indicators for predicting the pathogenesis of CKD and also might enable effective monitoring of the efficacy of CKD treatment.

Comprehensive characterization of the prostate tumor microenvironment identifies CXCR4/CXCL12 crosstalk as a novel antiangiogenic therapeutic target in prostate cancer.

BACKGROUND: Crosstalk between neoplastic and stromal cells fosters prostate cancer (PCa) progression and dissemination. Insight in cell-to-cell communication networks provides new therapeutic avenues to mold processes that contribute to PCa tumor microenvironment (TME) alterations. Here we performed a detailed characterization of PCa tumor endothelial cells (TEC) to delineate intercellular crosstalk between TEC and the PCa TME. METHODS: TEC isolated from 67 fresh radical prostatectomy (RP) specimens underwent multi-omic ex vivo characterization as well as orthogonal validation of both TEC functions and key markers by immunohistochemistry (IHC) and immunofluorescence (IF). To identify cell-cell interaction targets in TEC, we performed single-cell RNA sequencing (scRNA-seq) in four PCa patients who underwent a RP to catalogue cellular TME composition. Targets were cross-validated using IHC, publicly available datasets, cell culture expriments as well as a PCa xenograft mouse model. RESULTS: Compared to adjacent normal endothelial cells (NEC) bulk RNA-seq analysis revealed upregulation of genes associated with tumor vasculature, collagen modification and extracellular matrix remodeling in TEC. PTGIR, PLAC9, CXCL12 and VDR were identified as TEC markers and confirmed by IF and IHC in an independent patient cohort. By scRNA-seq we identified 27 cell (sub)types, including endothelial cells (EC) with arterial, venous and immature signatures, as well as angiogenic tip EC. A focused molecular analysis revealed that arterial TEC displayed highest CXCL12 mRNA expression levels when compared to all other TME cell (sub)populations and showed a negative prognostic role. Receptor-ligand interaction analysis predicted interactions between arterial TEC derived CXCL12 and its cognate receptor CXCR4 on angiogenic tip EC. CXCL12 was in vitro and in vivo validated as actionable TEC target by highlighting the vessel number- and density- reducing activity of the CXCR4-inhibitor AMD3100 in murine PCa as well as by inhibition of TEC proliferation and migration in vitro. CONCLUSIONS: Overall, our comprehensive analysis identified novel PCa TEC targets and highlights CXCR4/CXCL12 interaction as a potential novel target to interfere with tumor angiogenesis in PCa.

Smooth Muscle Specific Ablation of CXCL12 in Mice Downregulates CXCR7 Associated with Defective Coronary Arteries and Cardiac Hypertrophy.

The chemokine CXCL12 plays a fundamental role in cardiovascular development, cell trafficking, and myocardial repair. Human genome-wide association studies even have identified novel loci downstream of the CXCL12 gene locus associated with coronary artery disease and myocardial infarction. Nevertheless, cell and tissue specific effects of CXCL12 are barely understood. Since we detected high expression of CXCL12 in smooth muscle (SM) cells, we generated a SM22-alpha-Cre driven mouse model to ablate CXCL12 (SM-CXCL12-/-). SM-CXCL12-/- mice revealed high embryonic lethality (50%) with developmental defects, including aberrant topology of coronary arteries. Postnatally, SM-CXCL12-/- mice developed severe cardiac hypertrophy associated with fibrosis, apoptotic cell death, impaired heart function, and severe coronary vascular defects characterized by thinned and dilated arteries. Transcriptome analyses showed specific upregulation of pathways associated with hypertrophic cardiomyopathy, collagen protein network, heart-related proteoglycans, and downregulation of the M2 macrophage modulators. CXCL12 mutants showed endothelial downregulation of the CXCL12 co-receptor CXCR7. Treatment of SM-CXCL12-/- mice with the CXCR7 agonist TC14012 attenuated cardiac hypertrophy associated with increased pERK signaling. Our data suggest a critical role of smooth muscle-specific CXCL12 in arterial development, vessel maturation, and cardiac hypertrophy. Pharmacological stimulation of CXCR7 might be a promising target to attenuate adverse hypertrophic remodeling.

Correction to: Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq data.

It was highlighted that the original article [1] contained a typesetting mistake in the name of Noel Filipe da Cunha Carvalho de Miranda. This was incorrectly captured as Noel Filipe da Cunha Carvahlo de Miranda. It was also highlighted that in Fig. 3C the left panels Y-axis were cropped and in Fig. 5C, CD8 bar was cropped. This Correction article shows the correct Figs. 3 and 5. The original article has been updated.

Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq data.

We introduce quanTIseq, a method to quantify the fractions of ten immune cell types from bulk RNA-sequencing data. quanTIseq was extensively validated in blood and tumor samples using simulated, flow cytometry, and immunohistochemistry data.quanTIseq analysis of 8000 tumor samples revealed that cytotoxic T cell infiltration is more strongly associated with the activation of the CXCR3/CXCL9 axis than with mutational load and that deconvolution-based cell scores have prognostic value in several solid cancers. Finally, we used quanTIseq to show how kinase inhibitors modulate the immune contexture and to reveal immune-cell types that underlie differential patients' responses to checkpoint blockers.Availability: quanTIseq is available at http://icbi.at/quantiseq .

Nuclear receptor NR2F6 inhibition potentiates responses to PD-L1/PD-1 cancer immune checkpoint blockade.

Analyzing mouse tumor models in vivo, human T cells ex vivo, and human lung cancer samples, we provide direct evidence that NR2F6 acts as an immune checkpoint. Genetic ablation of Nr2f6, particularly in combination with established cancer immune checkpoint blockade, efficiently delays tumor progression and improves survival in experimental mouse models. The target genes deregulated in intratumoral T lymphocytes upon genetic ablation of Nr2f6 alone or together with PD-L1 blockade reveal multiple advantageous transcriptional alterations. Acute Nr2f6 silencing in both mouse and human T cells induces hyper-responsiveness that establishes a non-redundant T-cell-inhibitory function of NR2F6. NR2F6 protein expression in T-cell-infiltrating human NSCLC is upregulated in 54% of the cases (n = 303) and significantly correlates with PD-1 and CTLA-4 expression. Our data define NR2F6 as an intracellular immune checkpoint that suppresses adaptive anti-cancer immune responses and set the stage for clinical validation of targeting NR2F6 for next-generation immuno-oncological regimens.

Targeting immune checkpoints potentiates immunoediting and changes the dynamics of tumor evolution.

The cancer immunoediting hypothesis postulates a dual role of the immune system: protecting the host by eliminating tumor cells, and shaping the tumor by editing its genome. Here, we elucidate the impact of evolutionary and immune-related forces on editing the tumor in a mouse model for hypermutated and microsatellite-instable colorectal cancer. Analyses of wild-type and immunodeficient RAG1 knockout mice transplanted with MC38 cells reveal that upregulation of checkpoint molecules and infiltration by Tregs are the major tumor escape mechanisms. Our results show that the effects of immunoediting are weak and that neutral accumulation of mutations dominates. Targeting the PD-1/PD-L1 pathway using immune checkpoint blocker effectively potentiates immunoediting. The immunoediting effects are less pronounced in the CT26 cell line, a non-hypermutated/microsatellite-instable model. Our study demonstrates that neutral evolution is another force that contributes to sculpting the tumor and that checkpoint blockade effectively enforces T-cell-dependent immunoselective pressure.

NR4A3 Suppresses Lymphomagenesis through Induction of Proapoptotic Genes.

Nuclear orphan receptor NR4A1 exerts an essential tumor suppressor function in aggressive lymphomas. In this study, we investigated the hypothesized contribution of the related NR4A family member NR4A3 to lymphomagenesis. In aggressive lymphoma patients, low expression of NR4A3 was associated with poor survival. Ectopic expression or pharmacological activation of NR4A3 in lymphoma cell lines led to a significantly higher proportion of apoptotic cells. In a mouse NSG xenograft model of lymphoma (stably transduced SuDHL4 cells), NR4A3 expression abrogated tumor growth, compared with vector control and uninduced cells that formed massive tumors. Transcript analysis of four different aggressive lymphoma cell lines overexpressing either NR4A3 or NR4A1 revealed that apoptosis was driven similarly by induction of BAK, Puma, BIK, BIM, BID, and Trail. Overall, our results showed that NR4A3 possesses robust tumor suppressor functions of similar impact to NR4A1 in aggressive lymphomas. Cancer Res; 77(9); 2375-86. ©2017 AACR.

Impact of the Chromatin Remodeling Factor CHD1 on Gut Microbiome Composition of Drosophila melanogaster.

The composition of the intestinal microbiota of Drosophila has been studied in some detail in recent years. Environmental, developmental and host-specific genetic factors influence microbiome composition in the fly. Our previous work has indicated that intestinal bacterial load can be affected by chromatin-targeted regulatory mechanisms. Here we studied a potential role of the conserved chromatin assembly and remodeling factor CHD1 in the shaping of the gut microbiome in Drosophila melanogaster. Using high-throughput sequencing of 16S rRNA gene amplicons, we found that Chd1 deletion mutant flies exhibit significantly reduced microbial diversity compared to rescued control strains. Specifically, although Acetobacteraceae dominated the microbiota of both Chd1 wild-type and mutant guts, Chd1 mutants were virtually monoassociated with this bacterial family, whereas in control flies other bacterial taxa constituted ~20% of the microbiome. We further show age-linked differences in microbial load and microbiota composition between Chd1 mutant and control flies. Finally, diet supplementation experiments with Lactobacillus plantarum revealed that, in contrast to wild-type flies, Chd1 mutant flies were unable to maintain higher L. plantarum titres over time. Collectively, these data provide evidence that loss of the chromatin remodeler CHD1 has a major impact on the gut microbiome of Drosophila melanogaster.

Changes in the miRNA-mRNA Regulatory Network Precede Motor Symptoms in a Mouse Model of Multiple System Atrophy: Clinical Implications.

Multiple system atrophy (MSA) is a fatal rapidly progressive α-synucleinopathy, characterized by α-synuclein accumulation in oligodendrocytes. It is accepted that the pathological α-synuclein accumulation in the brain of MSA patients plays a leading role in the disease process, but little is known about the events in the early stages of the disease. In this study we aimed to define potential roles of the miRNA-mRNA regulatory network in the early pre-motor stages of the disease, i.e., downstream of α-synuclein accumulation in oligodendroglia, as assessed in a transgenic mouse model of MSA. We investigated the expression patterns of miRNAs and their mRNA targets in substantia nigra (SN) and striatum, two brain regions that undergo neurodegeneration at a later stage in the MSA model, by microarray and RNA-seq analysis, respectively. Analysis was performed at a time point when α-synuclein accumulation was already present in oligodendrocytes at neuropathological examination, but no neuronal loss nor deficits of motor function had yet occurred. Our data provide a first evidence for the leading role of gene dysregulation associated with deficits in immune and inflammatory responses in the very early, non-symptomatic disease stages of MSA. While dysfunctional homeostasis and oxidative stress were prominent in SN in the early stages of MSA, in striatum differential gene expression in the non-symptomatic phase was linked to oligodendroglial dysfunction, disturbed protein handling, lipid metabolism, transmembrane transport and altered cell death control, respectively. A large number of putative miRNA-mRNAs interaction partners were identified in relation to the control of these processes in the MSA model. Our results support the role of early changes in the miRNA-mRNA regulatory network in the pathogenesis of MSA preceding the clinical onset of the disease. The findings thus contribute to understanding the disease process and are likely to pave the way towards identifying disease biomarkers for early diagnosis of MSA.

Characterization of the immunophenotypes and antigenomes of colorectal cancers reveals distinct tumor escape mechanisms and novel targets for immunotherapy.

BACKGROUND: While large-scale cancer genomic projects are comprehensively characterizing the mutational spectrum of various cancers, so far little attention has been devoted to either define the antigenicity of these mutations or to characterize the immune responses they elicit. Here we present a strategy to characterize the immunophenotypes and the antigen-ome of human colorectal cancer. RESULTS: We apply our strategy to a large colorectal cancer cohort (n = 598) and show that subpopulations of tumor-infiltrating lymphocytes are associated with distinct molecular phenotypes. The characterization of the antigenome shows that a large number of cancer-germline antigens are expressed in all patients. In contrast, neo-antigens are rarely shared between patients, indicating that cancer vaccination requires individualized strategy. Analysis of the genetic basis of the tumors reveals distinct tumor escape mechanisms for the patient subgroups. Hypermutated tumors are depleted of immunosuppressive cells and show upregulation of immunoinhibitory molecules. Non-hypermutated tumors are enriched with immunosuppressive cells, and the expression of immunoinhibitors and MHC molecules is downregulated. Reconstruction of the interaction network of tumor-infiltrating lymphocytes and immunomodulatory molecules followed by a validation with 11 independent cohorts (n = 1,945) identifies BCMA as a novel druggable target. Finally, linear regression modeling identifies major determinants of tumor immunogenicity, which include well-characterized modulators as well as a novel candidate, CCR8, which is then tested in an orthologous immunodeficient mouse model. CONCLUSIONS: The immunophenotypes of the tumors and the cancer antigenome remain widely unexplored, and our findings represent a step toward the development of personalized cancer immunotherapies.

NR4A1-mediated apoptosis suppresses lymphomagenesis and is associated with a favorable cancer-specific survival in patients with aggressive B-cell lymphomas.

NR4A1 (Nur77) and NR4A3 (Nor-1) function as tumor suppressor genes as demonstrated by the rapid development of acute myeloid leukemia in the NR4A1 and NR4A3 knockout mouse. The aim of our study was to investigate NR4A1 and NR4A3 expression and function in lymphoid malignancies. We found a vastly reduced expression of NR4A1 and NR4A3 in chronic lymphocytic B-cell leukemia (71%), in follicular lymphoma (FL, 70%), and in diffuse large B-cell lymphoma (DLBCL, 74%). In aggressive lymphomas (DLBCL and FL grade 3), low NR4A1 expression was significantly associated with a non-germinal center B-cell subtype and with poor overall survival. To investigate the function of NR4A1 in lymphomas, we overexpressed NR4A1 in several lymphoma cell lines. Overexpression of NR4A1 led to a higher proportion of lymphoma cells undergoing apoptosis. To test the tumor suppressor function of NR4A1 in vivo, the stable lentiviral-transduced SuDHL4 lymphoma cell line harboring an inducible NR4A1 construct was further investigated in xenografts. Induction of NR4A1 abrogated tumor growth in the NSG mice, in contrast to vector controls, which formed massive tumors. Our data suggest that NR4A1 has proapoptotic functions in aggressive lymphoma cells and define NR4A1 as a novel gene with tumor suppressor properties involved in lymphomagenesis.

Co-expressed genes prepositioned in spatial neighborhoods stochastically associate with SC35 speckles and RNA polymerase II factories.

Chromosomally separated, co-expressed genes can be in spatial proximity, but there is still debate about how this nuclear organization is achieved. Proposed mechanisms include global genome organization, preferential positioning of chromosome territories, or gene-gene sharing of various nuclear bodies. To investigate this question, we selected a set of genes that were co-expressed upon differentiation of human multipotent stem cells. We applied a novel multi-dimensional analysis procedure which revealed that prior to gene expression, the relative position of these genes was conserved in nuclei. Upon stem cell differentiation and concomitant gene expression, we found that co-expressed genes were closer together. In addition, we found that genes in the same 1-µm-diameter neighborhood associated with either the same splicing speckle or to a lesser extent with the same transcription factory. Dispersal of speckles by overexpression of the serine-arginine (SR) protein kinase cdc2-like kinase Clk2 led to a significant drop in the number of genes in shared neighborhoods. We demonstrate quantitatively that the frequencies of speckle and factory sharing can be explained by assuming stochastic selection of a nuclear body within a restricted sub-volume defined by the original global gene positioning present prior to gene expression. We conclude that the spatial organization of these genes is a two-step process in which transcription-induced association with nuclear bodies enhances and refines a pre-existing global organization.