Clonal selection of stable aneuploidies in progenitor cells drives high-prevalence tumorigenesis.

Chromosome gains and losses are a frequent feature of human cancers. However, how these aberrations can outweigh the detrimental effects of aneuploidy remains unclear. An initial comparison of existing chromosomal instability (CIN) mouse models suggests that aneuploidy accumulates to low levels in these animals. We therefore developed a novel mouse model that enables unprecedented levels of chromosome missegregation in the adult animal. At the earliest stages of T-cell development, cells with random chromosome gains and/or losses are selected against, but CIN eventually results in the expansion of progenitors with clonal chromosomal imbalances. Clonal selection leads to the development of T-cell lymphomas with stereotypic karyotypes in which chromosome 15, containing the Myc oncogene, is gained with high prevalence. Expressing human MYC from chromosome 6 (MYCChr6) is sufficient to change the karyotype of these lymphomas to include universal chromosome 6 gains. Interestingly, while chromosome 15 is still gained in MYCChr6 tumors after genetic ablation of the endogenous Myc locus, this chromosome is not efficiently gained after deletion of one copy of Rad21, suggesting a synergistic effect of both MYC and RAD21 in driving chromosome 15 gains. Our results show that the initial detrimental effects of random missegregation are outbalanced by clonal selection, which is dictated by the chromosomal location and nature of certain genes and is sufficient to drive cancer with high prevalence.

Transformation of Accessible Chromatin and 3D Nucleome Underlies Lineage Commitment of Early T Cells.

How chromatin reorganization coordinates differentiation and lineage commitment from hematopoietic stem and progenitor cells (HSPCs) to mature immune cells has not been well understood. Here, we carried out an integrative analysis of chromatin accessibility, topologically associating domains, AB compartments, and gene expression from HSPCs to CD4+CD8+ T cells. We found that abrupt genome-wide changes at all three levels of chromatin organization occur during the transition from double-negative stage 2 (DN2) to DN3, accompanying the T lineage commitment. The transcription factor BCL11B, a critical regulator of T cell commitment, is associated with increased chromatin interaction, and Bcl11b deletion compromised chromatin interaction at its target genes. We propose that these large-scale and concerted changes in chromatin organization present an energy barrier to prevent the cell from reversing its fate to earlier stages or redirecting to alternatives and thus lock the cell fate into the T lineages.

A Muscle-Specific Enhancer RNA Mediates Cohesin Recruitment and Regulates Transcription In trans.

The enhancer regions of the myogenic master regulator MyoD give rise to at least two enhancer RNAs. Core enhancer eRNA (CEeRNA) regulates transcription of the adjacent MyoD gene, whereas DRReRNA affects expression of Myogenin in trans. We found that DRReRNA is recruited at the Myogenin locus, where it colocalizes with Myogenin nascent transcripts. DRReRNA associates with the cohesin complex, and this association correlates with its transactivating properties. Despite being expressed in undifferentiated cells, cohesin is not loaded on Myogenin until the cells start expressing DRReRNA, which is then required for cohesin chromatin recruitment and maintenance. Functionally, depletion of either cohesin or DRReRNA reduces chromatin accessibility, prevents Myogenin activation, and hinders muscle cell differentiation. Thus, DRReRNA ensures spatially appropriate cohesin loading in trans to regulate gene expression.

Characterization of genetically defined sporadic and hereditary type 1 papillary renal cell carcinoma cell lines.

Renal cell carcinoma (RCC) is not a single disease but is made up of several different histologically defined subtypes that are associated with distinct genetic alterations which require subtype specific management and treatment. Papillary renal cell carcinoma (pRCC) is the second most common subtype after conventional/clear cell RCC (ccRCC), representing ~20% of cases, and is subcategorized into type 1 and type 2 pRCC. It is important for preclinical studies to have cell lines that accurately represent each specific RCC subtype. This study characterizes seven cell lines derived from both primary and metastatic sites of type 1 pRCC, including the first cell line derived from a hereditary papillary renal carcinoma (HPRC)-associated tumor. Complete or partial gain of chromosome 7 was observed in all cell lines and other common gains of chromosomes 16, 17, or 20 were seen in several cell lines. Activating mutations of MET were present in three cell lines that all demonstrated increased MET phosphorylation in response to HGF and abrogation of MET phosphorylation in response to MET inhibitors. CDKN2A loss due to mutation or gene deletion, associated with poor outcomes in type 1 pRCC patients, was observed in all cell line models. Six cell lines formed tumor xenografts in athymic nude mice and thus provide in vivo models of type 1 pRCC. These type 1 pRCC cell lines provide a comprehensive representation of the genetic alterations associated with pRCC that will give insight into the biology of this disease and be ideal preclinical models for therapeutic studies.

Trac-looping measures genome structure and chromatin accessibility.

Long-range chromatin interactions play critical roles in genome organization and regulation of transcription. We now report transposase-mediated analysis of chromatin looping (Trac-looping) for simultaneous detection of multiscale genome-wide chromatin interactions among regulatory elements and chromatin accessibility. With this technique, a bivalent oligonucleotide linker is inserted between two interacting regions such that the chromatin interactions are captured without prior chromatin fragmentation and proximity-based ligation. Application of Trac-looping to human CD4+ T cells revealed substantial reorganization of enhancer-promoter interactions associated with changes in gene expression after T cell receptor stimulation.

Novel renal medullary carcinoma cell lines, UOK353 and UOK360, provide preclinical tools to identify new therapeutic treatments.

Renal medullary carcinoma (RMC) is a rare, aggressive disease that predominantly afflicts individuals of African or Mediterranean descent with sickle cell trait. RMC comprises 1% of all renal cell carcinoma diagnoses with a median overall survival of 13 months. Patients are typically young (median age-22) and male (male:female ratio of 2:1) and tumors are characterized by complete loss of expression of the SMARCB1 tumor suppressor protein. Due to the low incidence of RMC and the disease's aggressiveness, treatment decisions are often based on case reports. Thus, it is critical to develop preclinical models of RMC to better understand the pathogenesis of this disease and to identify effective forms of therapy. Two novel cell line models, UOK353 and UOK360, were derived from primary RMCs that both demonstrated the characteristic SMARCB1 loss. Both cell lines overexpressed EZH2 and other members of the polycomb repressive complex and EZH2 inhibition in RMC tumor spheroids resulted in decreased viability. High throughput drug screening of both cell lines revealed several additional candidate compounds, including bortezomib that had both in vitro and in vivo antitumor activity. The activity of bortezomib was shown to be partially dependent on increased oxidative stress as addition of the N-acetyl cysteine antioxidant reduced the effect on cell proliferation. Combining bortezomib and cisplatin further decreased cell viability both in vitro and in vivo that single agent bortezomib treatment. The UOK353 and UOK360 cell lines represent novel preclinical models for the development of effective forms of therapy for RMC patients.

Single-cell genetic analysis of clonal dynamics in colorectal adenomas indicates CDX2 gain as a predictor of recurrence.

Colorectal adenomas are common precancerous lesions with the potential for malignant transformation to colorectal adenocarcinoma. Endoscopic polypectomy provides an opportunity for cancer prevention; however, recurrence rates are high. We collected formalin-fixed paraffin-embedded tissue of 15 primary adenomas with recurrence, 15 adenomas without recurrence, and 14 matched pair samples (primary adenoma and the corresponding recurrent adenoma). The samples were analysed by array-comparative genomic hybridisation (aCGH) and single-cell multiplex interphase fluorescence in situ hybridisation (miFISH) to understand clonal evolution, to examine the dynamics of copy number alterations (CNAs) and to identify molecular markers for recurrence prediction. The miFISH probe panel consisted of 14 colorectal carcinogenesis-relevant genes (COX2, PIK3CA, APC, CLIC1, EGFR, MYC, CCND1, CDX2, CDH1, TP53, HER2, SMAD7, SMAD4 and ZNF217), and a centromere probe (CEP10). The aCGH analysis confirmed the genetic landscape typical for colorectal tumorigenesis, that is, CNAs of chromosomes 7, 13q, 18 and 20q. Focal aberrations (≤10 Mbp) were mapped to chromosome bands 6p22.1-p21.33 (33.3%), 7q22.1 (31.4%) and 16q21 (29.4%). MiFISH detected gains of EGFR (23.6%), CDX2 (21.8%) and ZNF217 (18.2%). Most adenomas exhibited a major clone population which was accompanied by multiple smaller clone populations. Gains of CDX2 were exclusively seen in primary adenomas with recurrence (25%) compared to primary adenomas without recurrence (0%). Generation of phylogenetic trees for matched pair samples revealed four distinct patterns of clonal dynamics. In conclusion, adenoma development and recurrence are complex genetic processes driven by multiple CNAs whose evaluations by miFISH, with emphasis on CDX2, might serve as a predictor of recurrence.

Near-tetraploid cancer cells show chromosome instability triggered by replication stress and exhibit enhanced invasiveness.

A considerable proportion of tumors exhibit aneuploid karyotypes, likely resulting from the progressive loss of chromosomes after whole-genome duplication. Here, by using isogenic diploid and near-tetraploid (4N) single-cell-derived clones from the same parental cell lines, we aimed at exploring how polyploidization affects cellular functions and how tetraploidy generates chromosome instability. Gene expression profiling in 4N clones revealed a significant enrichment of transcripts involved in cell cycle and DNA replication. Increased levels of replication stress in 4N cells resulted in DNA damage, impaired proliferation caused by a cell cycle delay during S phase, and higher sensitivity to S phase checkpoint inhibitors. In fact, increased levels of replication stress were also observed in nontransformed, proliferative posttetraploid RPE1 cells. Additionally, replication stress promoted higher levels of intercellular genomic heterogeneity and ongoing genomic instability, which could be explained by high rates of mitotic defects, and was alleviated by the supplementation of exogenous nucleosides. Finally, our data found that 4N cancer cells displayed increased migratory and invasive capacity, both in vitro and in primary colorectal tumors, indicating that tetraploidy can promote aggressive cancer cell behavior.-Wangsa, D., Quintanilla, I., Torabi, K., Vila-Casadesús, M., Ercilla, A., Klus, G., Yuce, Z., Galofré, C., Cuatrecasas, M., Lozano, J. J., Agell, N., Cimini, D., Castells, A., Ried, T., Camps, J. Near-tetraploid cancer cells show chromosome instability triggered by replication stress and exhibit enhanced invasiveness.

HiCTMap: Detection and analysis of chromosome territory structure and position by high-throughput imaging.

The spatial organization of chromosomes in the nuclear space is an extensively studied field that relies on measurements of structural features and 3D positions of chromosomes with high precision and robustness. However, no tools are currently available to image and analyze chromosome territories in a high-throughput format. Here, we have developed High-throughput Chromosome Territory Mapping (HiCTMap), a method for the robust and rapid analysis of 2D and 3D chromosome territory positioning in mammalian cells. HiCTMap is a high-throughput imaging-based chromosome detection method which enables routine analysis of chromosome structure and nuclear position. Using an optimized FISH staining protocol in a 384-well plate format in conjunction with a bespoke automated image analysis workflow, HiCTMap faithfully detects chromosome territories and their position in 2D and 3D in a large population of cells per experimental condition. We apply this novel technique to visualize chromosomes 18, X, and Y in male and female primary human skin fibroblasts, and show accurate detection of the correct number of chromosomes in the respective genotypes. Given the ability to visualize and quantitatively analyze large numbers of nuclei, we use HiCTMap to measure chromosome territory area and volume with high precision and determine the radial position of chromosome territories using either centroid or equidistant-shell analysis. The HiCTMap protocol is also compatible with RNA FISH as demonstrated by simultaneous labeling of X chromosomes and Xist RNA in female cells. We suggest HiCTMap will be a useful tool for routine precision mapping of chromosome territories in a wide range of cell types and tissues.

The evolution of single cell-derived colorectal cancer cell lines is dominated by the continued selection of tumor-specific genomic imbalances, despite random chromosomal instability.

Intratumor heterogeneity is a major challenge in cancer treatment. To decipher patterns of chromosomal heterogeneity, we analyzed six colorectal cancer cell lines by multiplex interphase FISH (miFISH). The mismatch-repair-deficient cell lines DLD-1 and HCT116 had the most stable copy numbers, whereas aneuploid cell lines (HT-29, SW480, SW620 and H508) displayed a higher degree of instability. We subsequently assessed the clonal evolution of single cells in two colorectal carcinoma cell lines, SW480 and HT-29, which both have aneuploid karyotypes but different degrees of chromosomal instability. The clonal compositions of the single cell-derived daughter lines, as assessed by miFISH, differed for HT-29 and SW480. Daughters of HT-29 were stable, clonal, with little heterogeneity. Daughters of SW480 were more heterogeneous, with the single cell-derived daughter lines separating into two distinct populations with different ploidy (hyper-diploid and near-triploid), morphology, gene expression and tumorigenicity. To better understand the evolutionary trajectory for the two SW480 populations, we constructed phylogenetic trees which showed ongoing instability in the daughter lines. When analyzing the evolutionary development over time, most single cell-derived daughter lines maintained their major clonal pattern, with the exception of one daughter line that showed a switch involving a loss of APC. Our meticulous analysis of the clonal evolution and composition of these colorectal cancer models shows that all chromosomes are subject to segregation errors, however, specific net genomic imbalances are maintained. Karyotype evolution is driven by the necessity to arrive at and maintain a specific plateau of chromosomal copy numbers as the drivers of carcinogenesis.

Transcription-dependent radial distribution of TCF7L2 regulated genes in chromosome territories.

Human chromosomes occupy distinct territories in the interphase nucleus. Such chromosome territories (CTs) are positioned according to gene density. Gene-rich CTs are generally located in the center of the nucleus, while gene-poor CTs are positioned more towards the nuclear periphery. However, the association between gene expression levels and the radial positioning of genes within the CT is still under debate. In the present study, we performed three-dimensional fluorescence in situ hybridization experiments in the colorectal cancer cell lines DLD-1 and LoVo using whole chromosome painting probes for chromosomes 8 and 11 and BAC clones targeting four genes with different expression levels assessed by gene expression arrays and RT-PCR. Our results confirmed that the two over-expressed genes, MYC on chromosome 8 and CCND1 on chromosome 11, are located significantly further away from the center of the CT compared to under-expressed genes on the same chromosomes, i.e., DLC1 and SCN3B. When CCND1 expression was reduced after silencing the major transcription factor of the WNT/ß-catenin signaling pathway, TCF7L2, the gene was repositioned and mostly detected in the interior of the CT. Thus, we suggest a non-random distribution in which over-expressed genes are located more towards the periphery of the respective CTs.

Phylogenetic analysis of multiple FISH markers in oral tongue squamous cell carcinoma suggests that a diverse distribution of copy number changes is associated with poor prognosis.

Oral tongue squamous cell carcinoma (OTSCC) is associated with poor prognosis. To improve prognostication, we analyzed four gene probes (TERC, CCND1, EGFR and TP53) and the centromere probe CEP4 as a marker of chromosomal instability, using fluorescence in situ hybridization (FISH) in single cells from the tumors of sixty-five OTSCC patients (Stage I, n = 15; Stage II, n = 30; Stage III, n = 7; Stage IV, n = 13). Unsupervised hierarchical clustering of the FISH data distinguished three clusters related to smoking status. Copy number increases of all five markers were found to be correlated to non-smoking habits, while smokers in this cohort had low-level copy number gains. Using the phylogenetic modeling software FISHtrees, we constructed models of tumor progression for each patient based on the four gene probes. Then, we derived test statistics on the models that are significant predictors of disease-free and overall survival, independent of tumor stage and smoking status in multivariate analysis. The patients whose tumors were modeled as progressing by a more diverse distribution of copy number changes across the four genes have poorer prognosis. This is consistent with the view that multiple genetic pathways need to become deregulated in order for cancer to progress.

Inferring models of multiscale copy number evolution for single-tumor phylogenetics.

MOTIVATION: Phylogenetic algorithms have begun to see widespread use in cancer research to reconstruct processes of evolution in tumor progression. Developing reliable phylogenies for tumor data requires quantitative models of cancer evolution that include the unusual genetic mechanisms by which tumors evolve, such as chromosome abnormalities, and allow for heterogeneity between tumor types and individual patients. Previous work on inferring phylogenies of single tumors by copy number evolution assumed models of uniform rates of genomic gain and loss across different genomic sites and scales, a substantial oversimplification necessitated by a lack of algorithms and quantitative parameters for fitting to more realistic tumor evolution models. RESULTS: We propose a framework for inferring models of tumor progression from single-cell gene copy number data, including variable rates for different gain and loss events. We propose a new algorithm for identification of most parsimonious combinations of single gene and single chromosome events. We extend it via dynamic programming to include genome duplications. We implement an expectation maximization (EM)-like method to estimate mutation-specific and tumor-specific event rates concurrently with tree reconstruction. Application of our algorithms to real cervical cancer data identifies key genomic events in disease progression consistent with prior literature. Classification experiments on cervical and tongue cancer datasets lead to improved prediction accuracy for the metastasis of primary cervical cancers and for tongue cancer survival. AVAILABILITY AND IMPLEMENTATION: Our software (FISHtrees) and two datasets are available at ftp://ftp.ncbi.nlm.nih.gov/pub/FISHtrees.

Inflammation-mediated genetic and epigenetic alterations drive cancer development in the neighboring epithelium upon stromal abrogation of TGF-ß signaling.

Deletion of tumor suppressor genes in stromal fibroblasts induces epithelial cancer development, suggesting an important role of stroma in epithelial homoeostasis. However, the underlying mechanisms remain to be elucidated. Here we report that deletion of the gene encoding TGFß receptor 2 (Tgfbr2) in the stromal fibroblasts (Tgfbr2(fspKO)) induces inflammation and significant DNA damage in the neighboring epithelia of the forestomach. This results in loss or down-regulation of cyclin-dependent kinase inhibitors p15, p16, and p21, which contribute to the development of invasive squamous cell carcinoma (SCC). Anti-inflammation treatment restored p21 expression, delayed tumorigenesis, and increased survival of Tgfbr2(fspKO) mice. Our data demonstrate for the first time that inflammation is a critical player in the epigenetic silencing of p21 in tumor progression. Examination of human esophageal SCC showed a down-regulation of TGFß receptor 2 (TßRII) in the stromal fibroblasts, as well as increased inflammation, DNA damage, and loss or decreased p15/p16 expression. Our study suggests anti-inflammation may be a new therapeutic option in treating human SCCs with down-regulation of TßRII in the stroma.

Loss of lamin B1 results in prolongation of S phase and decondensation of chromosome territories.

Nuclear lamin B1 (LMNB1) constitutes one of the major structural proteins in the lamina mesh. We silenced the expression of LMNB1 by RNA interference in the colon cancer cell line DLD-1 and showed a dramatic redistribution of H3K27me3 from the periphery to a more homogeneous nuclear dispersion. In addition, we observed telomere attrition and an increased frequency of micronuclei and nuclear blebs. By 3D-FISH analyses, we demonstrated that the volume and surface of chromosome territories were significantly larger in LMNB1-depleted cells, suggesting that LMNB1 is required to maintain chromatin condensation in interphase nuclei. These changes led to a prolonged S phase due to activation of Chk1. Finally, silencing of LMNB1 resulted in extensive changes in alternative splicing of multiple genes and in a higher number of enlarged nuclear speckles. Taken together, our results suggest a mechanistic role of the nuclear lamina in the organization of chromosome territories, maintenance of genome integrity and proper gene splicing.

Concurrent V(D)J recombination and DNA end instability increase interchromosomal trans-rearrangements in ATM-deficient thymocytes.

During the CD4(-)CD8(-) (DN) stage of T-cell development, RAG-dependent DNA breaks and V(D)J recombination occur at three T-cell receptor (TCR) loci: TCRß, TCRγ and TCRδ. During this stage, abnormal trans-rearrangements also take place between TCR loci, occurring at increased frequency in absence of the DNA damage response mediator ataxia telangiectasia mutated (ATM). Here, we use this model of physiologic trans-rearrangement to study factors that predispose to rearrangement and the role of ATM in preventing chromosomal translocations. The frequency of DN thymocytes with DNA damage foci at multiple TCR loci simultaneously is increased 2- to 3-fold in the absence of ATM. However, trans-rearrangement is increased 10 000- to 100 000-fold, indicating that ATM function extends beyond timely resolution of DNA breaks. RAG-mediated synaptic complex formation occurs between recombination signal sequences with unequal 12 and 23 base spacer sequences (12/23 rule). TCR trans-rearrangements violate this rule, as we observed similar frequencies of 12/23 and aberrant 12/12 or 23/23 recombination products. This suggests that trans-rearrangements are not the result of trans-synaptic complex formation, but they are instead because of unstable cis synaptic complexes that form simultaneously at distinct TCR loci. Thus, ATM suppresses trans-rearrangement primarily through stabilization of DNA breaks at TCR loci.

Haplotype-resolved karyotype construction from Hi-C data using refLinker.

Chromosomal aberrations are prevalent in cancer genomes, yet it remains challenging to resolve the long-range structure of rearranged chromosomes. A key problem is to determine the chromosomal origin of rearranged genomic segments, which requires chromosome-length haplotype information. Here we describe refLinker, a new computational method for whole-chromosome haplotype inference using external reference panels and Hi-C. We show that refLinker ensures consistent long-range phasing accuracy in both diploid human genomes and aneuploid cancers, including regions with loss-of-heterozygosity and high-level focal amplification. We further demonstrate the feasibility of complex genome reconstruction using haplotype-specific Hi-C contacts, revealing new karyotype features in two widely studied cancer cell lines. Together, these findings provide a new framework for the complete resolution of long-range chromosome structure in complex genomes and highlight the unique advantages of Hi-C data for reconstructing cancer genomes with chromosome-scale continuity.

Extrachromosomal Amplification of Human Papillomavirus Episomes Is a Mechanism of Cervical Carcinogenesis.

SIGNIFICANCE: Multimers of the HPV genome are generated in cervical tumors replicating as extrachromosomal episomes, which is associated with deletion and rearrangement of the HPV genome and provides a mechanism for oncogenesis without integration.

Intratumoral Heterogeneity and Clonal Evolution Induced by HPV Integration.

The human papillomavirus (HPV) genome is integrated into host DNA in most HPV-positive cancers, but the consequences for chromosomal integrity are unknown. Continuous long-read sequencing of oropharyngeal cancers and cancer cell lines identified a previously undescribed form of structural variation, "heterocateny," characterized by diverse, interrelated, and repetitive patterns of concatemerized virus and host DNA segments within a cancer. Unique breakpoints shared across structural variants facilitated stepwise reconstruction of their evolution from a common molecular ancestor. This analysis revealed that virus and virus-host concatemers are unstable and, upon insertion into and excision from chromosomes, facilitate capture, amplification, and recombination of host DNA and chromosomal rearrangements. Evidence of heterocateny was detected in extrachromosomal and intrachromosomal DNA. These findings indicate that heterocateny is driven by the dynamic, aberrant replication and recombination of an oncogenic DNA virus, thereby extending known consequences of HPV integration to include promotion of intratumoral heterogeneity and clonal evolution. SIGNIFICANCE: Long-read sequencing of HPV-positive cancers revealed "heterocateny," a previously unreported form of genomic structural variation characterized by heterogeneous, interrelated, and repetitive genomic rearrangements within a tumor. Heterocateny is driven by unstable concatemerized HPV genomes, which facilitate capture, rearrangement, and amplification of host DNA, and promotes intratumoral heterogeneity and clonal evolution. See related commentary by McBride and White, p. 814. This article is highlighted in the In This Issue feature, p. 799.

ATR inhibition augments the efficacy of lurbinectedin in small-cell lung cancer.

Small-cell lung cancer (SCLC) is the most lethal type of lung cancer. Specifically, MYC-driven non-neuroendocrine SCLC is particularly resistant to standard therapies. Lurbinectedin was recently approved for the treatment of relapsed SCLC, but combinatorial approaches are needed to increase the depth and duration of responses to lurbinectedin. Using high-throughput screens, we found inhibitors of ataxia telangiectasia mutated and rad3 related (ATR) as the most effective agents for augmenting lurbinectedin efficacy. First-in-class ATR inhibitor berzosertib synergized with lurbinectedin in multiple SCLC cell lines, organoid, and in vivo models. Mechanistically, ATR inhibition abrogated S-phase arrest induced by lurbinectedin and forced cell cycle progression causing mitotic catastrophe and cell death. High CDKN1A/p21 expression was associated with decreased synergy due to G1 arrest, while increased levels of ERCC5/XPG were predictive of increased combination efficacy. Importantly, MYC-driven non-neuroendocrine tumors which are resistant to first-line therapies show reduced CDKN1A/p21 expression and increased ERCC5/XPG indicating they are primed for response to lurbinectedin-berzosertib combination. The combination is being assessed in a clinical trial NCT04802174.