Temporal chromatin accessibility changes define transcriptional states essential for osteosarcoma metastasis.

The metastasis-invasion cascade describes the series of steps required for a cancer cell to successfully spread from its primary tumor and ultimately grow within a secondary organ. Despite metastasis being a dynamic, multistep process, most omics studies to date have focused on comparing primary tumors to the metastatic deposits that define end-stage disease. This static approach means we lack information about the genomic and epigenomic changes that occur during the majority of tumor progression. One particularly understudied phase of tumor progression is metastatic colonization, during which cells must adapt to the new microenvironment of the secondary organ. Through temporal profiling of chromatin accessibility and gene expression in vivo, we identify dynamic changes in the epigenome that occur as osteosarcoma tumors form and grow within the lung microenvironment. Furthermore, we show through paired in vivo and in vitro CRISPR drop-out screens and pharmacological validation that the upstream transcription factors represent a class of metastasis-specific dependency genes. While current models depict lung colonization as a discrete step within the metastatic cascade, our study shows it is a defined trajectory through multiple epigenetic states, revealing new therapeutic opportunities undetectable with standard approaches.

Ex vivo screen identifies CDK12 as a metastatic vulnerability in osteosarcoma.

Despite progress in intensification of therapy, outcomes for patients with metastatic osteosarcoma (OS) have not improved in thirty years. We developed a system that enabled preclinical screening of compounds against metastatic OS cells in the context of the native lung microenvironment. Using this strategy to screen a library of epigenetically targeted compounds, we identified inhibitors of CDK12 to be most effective, reducing OS cell outgrowth in the lung by more than 90% at submicromolar doses. We found that knockout of CDK12 in an in vivo model of lung metastasis significantly decreased the ability of OS to colonize the lung. CDK12 inhibition led to defects in transcription elongation in a gene length- and expression-dependent manner. These effects were accompanied by defects in RNA processing and altered the expression of genes involved in transcription regulation and the DNA damage response. We further identified OS models that differ in their sensitivity to CDK12 inhibition in the lung and provided evidence that upregulated MYC levels may mediate these differences. Our studies provided a framework for rapid preclinical testing of compounds with antimetastatic activity and highlighted CDK12 as a potential therapeutic target in OS.

Mismatch repair-signature mutations activate gene enhancers across human colorectal cancer epigenomes.

Commonly-mutated genes have been found for many cancers, but less is known about mutations in cis-regulatory elements. We leverage gains in tumor-specific enhancer activity, coupled with allele-biased mutation detection from H3K27ac ChIP-seq data, to pinpoint potential enhancer-activating mutations in colorectal cancer (CRC). Analysis of a genetically-diverse cohort of CRC specimens revealed that microsatellite instable (MSI) samples have a high indel rate within active enhancers. Enhancers with indels show evidence of positive selection, increased target gene expression, and a subset is highly recurrent. The indels affect short homopolymer tracts of A/T and increase affinity for FOX transcription factors. We further demonstrate that signature mismatch-repair (MMR) mutations activate enhancers using a xenograft tumor metastasis model, where mutations are induced naturally via CRISPR/Cas9 inactivation of MLH1 prior to tumor cell injection. Our results suggest that MMR signature mutations activate enhancers in CRC tumor epigenomes to provide a selective advantage.

Corrigendum: Positively selected enhancer elements endow osteosarcoma cells with metastatic competence.

This corrects the article DOI: 10.1038/nm.4475.

Positively selected enhancer elements endow osteosarcoma cells with metastatic competence.

Metastasis results from a complex set of traits acquired by tumor cells, distinct from those necessary for tumorigenesis. Here, we investigate the contribution of enhancer elements to the metastatic phenotype of osteosarcoma. Through epigenomic profiling, we identify substantial differences in enhancer activity between primary and metastatic human tumors and between near isogenic pairs of highly lung metastatic and nonmetastatic osteosarcoma cell lines. We term these regions metastatic variant enhancer loci (Met-VELs). Met-VELs drive coordinated waves of gene expression during metastatic colonization of the lung. Met-VELs cluster nonrandomly in the genome, indicating that activity of these enhancers and expression of their associated gene targets are positively selected. As evidence of this causal association, osteosarcoma lung metastasis is inhibited by global interruptions of Met-VEL-associated gene expression via pharmacologic BET inhibition, by knockdown of AP-1 transcription factors that occupy Met-VELs, and by knockdown or functional inhibition of individual genes activated by Met-VELs, such as that encoding coagulation factor III/tissue factor (F3). We further show that genetic deletion of a single Met-VEL at the F3 locus blocks metastatic cell outgrowth in the lung. These findings indicate that Met-VELs and the genes they regulate play a functional role in metastasis and may be suitable targets for antimetastatic therapies.

The thiirane-based selective MT1-MMP/MMP2 inhibitor ND-322 reduces melanoma tumor growth and delays metastatic dissemination.

MT1-MMP and MMP2 have been implicated as pro-tumorigenic and pro-metastatic factors in a wide variety of cancers including melanoma. We have previously demonstrated that MT1-MMP is highly expressed in melanoma where it promotes melanoma cell invasion and metastasis in part through the activation of its target MMP2. Given the accessibility of MMPs, as they are either secreted (e.g. MMP2) or membrane-tethered (e.g. MT1-MMP), they represent ideal targets for specific inhibition via small molecules. Here we show that the novel small-molecule inhibitor ND-322 with high selectivity for MT1-MMP and MMP2, effectively inhibits MT1-MMP and MMP2 activity resulting in reduced in vitro melanoma cell growth, migration and invasion. Importantly, these inhibitory effects lead to significant reduction of melanoma tumor growth and metastasis. We further show that while cell migration and invasion could be similarly hampered by specific inhibition of either MT1-MMP or MMP2 via shRNAs, the growth inhibitory activity of ND-322 could only be mirrored by specific inhibition of MT1-MMP. These data support ND-322 as a novel effective inhibitor capable of counteracting both MT1-MMP and MMP2, two key proteases involved in melanoma growth and metastasis. ND-322 may therefore represent a new inhibitor in the repertoire of treatments against melanoma.

Modeling disease risk through analysis of physical interactions between genetic variants within chromatin regulatory circuitry.

SNPs associated with disease susceptibility often reside in enhancer clusters, or super-enhancers. Constituents of these enhancer clusters cooperate to regulate target genes and often extend beyond the linkage disequilibrium (LD) blocks containing risk SNPs identified in genome-wide association studies (GWAS). We identified 'outside variants', defined as SNPs in weak LD with GWAS risk SNPs that physically interact with risk SNPs as part of a target gene's regulatory circuitry. These outside variants further explain variation in target gene expression beyond that explained by GWAS-associated SNPs. Additionally, the clinical risk associated with GWAS SNPs is considerably modified by the genotype of outside variants. Collectively, these findings suggest a potential model in which outside variants and GWAS SNPs that physically interact in 3D chromatin collude to influence target transcript levels as well as clinical risk. This model offers an additional hypothesis for the source of missing heritability for complex traits.

Transcription elongation factor ELL2 drives Ig secretory-specific mRNA production and the unfolded protein response.

Differentiation of B cells into Ab-secreting cells induces changes in gene transcription, IgH RNA processing, the unfolded protein response (UPR), and cell architecture. The transcription elongation factor eleven nineteen lysine-rich leukemia gene (ELL2) stimulates the processing of the secreted form of the IgH mRNA from the H chain gene. Mice (mus musculus) with the ELL2 gene floxed in either exon 1 or exon 3 were constructed and crossed to CD19-driven cre/CD19(+). The B cell-specific ELL2 conditional knockouts (cKOs; ell2(loxp/loxp) CD19(cre/+)) exhibit curtailed humoral responses both in 4-hydroxy-3-nitrophenyl acetyl-Ficoll and in 4-hydroxy-3-nitrophenyl acetyl-keyhole limpet hemocyanin immunized animals; recall responses were also diminished. The number of immature and recirculating B cells in the bone marrow is increased in the cKOs, whereas plasma cells in spleen are reduced relative to control animals. There are fewer IgG1 Ab-producing cells in the bone marrow of cKOs. LPS ex vivo-stimulated B220(lo)CD138(+) cells from ELL2-deficient mouse spleens are 4-fold less abundant than from control splenic B cells; have a paucity of secreted IgH; and have distended, abnormal-appearing endoplasmic reticulum. IRE1α is efficiently phosphorylated, but the amounts of Ig κ, ATF6, BiP, Cyclin B2, OcaB (BOB1, Pou2af1), and XBP1 mRNAs, unspliced and spliced, are severely reduced in ELL2-deficient cells. ELL2 enhances the expression of BCMA (also known as Tnfrsf17), which is important for long-term survival. Transcription yields from the cyclin B2 and the canonical UPR promoter elements are upregulated by ELL2 cDNA. Thus, ELL2 is important for many aspects of Ab secretion, XBP1 expression, and the UPR.

Plasma cell formation, secretion, and persistence: the short and the long of it.

B cells can be activated by cognate antigen, anti-B-cell receptor antibody, complement receptors, or polyclonal stimulators like lipopolysaccharide; the overall result is a large shift in RNA processing to the secretory-specific form of immunoglobulin (Ig) heavy chain mRNA and an upregulation of Igh mRNA amounts. Associated with this shift is the large-scale induction of Ig protein synthesis and the unfolded protein response to accommodate the massive quantity of secretory Ig that results. Stimulation to secretion also produces major structural accommodations and stress, with extensive generation of endoplasmic reticulum and Golgi as part of the cellular architecture. Reactive oxygen species can lead to either activation or apoptosis based on context and the high or low oxygen tension surrounding the cells. Transcription elongation factor ELL2 plays an important role in the induction of Ig secretory mRNA production, the unfolded protein response, and gene expression during hypoxia. After antigen stimulation, activated B cells from either the marginal zones or follicles can produce short-lived antibody secreting cells; it is not clear whether cells from both locations can become long-lived plasma cells. Autophagy is necessary for plasma cell long-term survival through the elimination of some of the accumulated damage to the ER from producing so much protein. Survival signals from the bone marrow stromal cells also contribute to plasma cell longevity, with BCMA serving a potentially unique survival role. Integrating the various information pathways converging on the plasma cell is crucial to the development of their long-lived, productive immune response.

Expanding the diversity of mycobacteriophages: insights into genome architecture and evolution.

Mycobacteriophages are viruses that infect mycobacterial hosts such as Mycobacterium smegmatis and Mycobacterium tuberculosis. All mycobacteriophages characterized to date are dsDNA tailed phages, and have either siphoviral or myoviral morphotypes. However, their genetic diversity is considerable, and although sixty-two genomes have been sequenced and comparatively analyzed, these likely represent only a small portion of the diversity of the mycobacteriophage population at large. Here we report the isolation, sequencing and comparative genomic analysis of 18 new mycobacteriophages isolated from geographically distinct locations within the United States. Although no clear correlation between location and genome type can be discerned, these genomes expand our knowledge of mycobacteriophage diversity and enhance our understanding of the roles of mobile elements in viral evolution. Expansion of the number of mycobacteriophages grouped within Cluster A provides insights into the basis of immune specificity in these temperate phages, and we also describe a novel example of apparent immunity theft. The isolation and genomic analysis of bacteriophages by freshman college students provides an example of an authentic research experience for novice scientists.