ALYREF, a novel factor involved in breast carcinogenesis, acts through transcriptional and post-transcriptional mechanisms selectively regulating the short NEAT1 isoform.

The RNA-binding protein ALYREF (THOC4) is involved in transcriptional regulation and nuclear mRNA export, though its role and molecular mode of action in breast carcinogenesis are completely unknown. Here, we identified high ALYREF expression as a factor for poor survival in breast cancer patients. ALYREF significantly influenced cellular growth, apoptosis and mitochondrial energy metabolism in breast cancer cells as well as breast tumorigenesis in orthotopic mouse models. Transcriptional profiling, phenocopy and rescue experiments identified the short isoform of the lncRNA NEAT1 as a molecular trigger for ALYREF effects in breast cancer. Mechanistically, we found that ALYREF binds to the NEAT1 promoter region to enhance the global NEAT1 transcriptional activity. Importantly, by stabilizing CPSF6, a protein that selectively activates the post-transcriptional generation of the short isoform of NEAT1, as well as by direct binding and stabilization of the short isoform of NEAT1, ALYREF selectively fine-tunes the expression of the short NEAT1 isoform. Overall, our study describes ALYREF as a novel factor contributing to breast carcinogenesis and identifies novel molecular mechanisms of regulation the two isoforms of NEAT1.

Complementary omics strategies to dissect p53 signaling networks under nutrient stress.

Signaling trough p53is a major cellular stress response mechanism and increases upon nutrient stresses such as starvation. Here, we show in a human hepatoma cell line that starvation leads to robust nuclear p53 stabilization. Using BioID, we determine the cytoplasmic p53 interaction network within the immediate-early starvation response and show that p53 is dissociated from several metabolic enzymes and the kinase PAK2 for which direct binding with the p53 DNA-binding domain was confirmed with NMR studies. Furthermore, proteomics after p53 immunoprecipitation (RIME) uncovered the nuclear interactome under prolonged starvation, where we confirmed the novel p53 interactors SORBS1 (insulin receptor signaling) and UGP2 (glycogen synthesis). Finally, transcriptomics after p53 re-expression revealed a distinct starvation-specific transcriptome response and suggested previously unknown nutrient-dependent p53 target genes. Together, our complementary approaches delineate several nodes of the p53 signaling cascade upon starvation, shedding new light on the mechanisms of p53 as nutrient stress sensor. Given the central role of p53 in cancer biology and the beneficial effects of fasting in cancer treatment, the identified interaction partners and networks could pinpoint novel pharmacologic targets to fine-tune p53 activity.

ß-catenin regulates FOXP2 transcriptional activity via multiple binding sites.

The transcription factor forkhead box protein P2 (FOXP2) is a highly conserved key regulator of embryonal development. The molecular mechanisms of how FOXP2 regulates embryonal development, however, remain elusive. Using RNA sequencing, we identified the Wnt signaling pathway as key target of FOXP2-dependent transcriptional regulation. Using cell-based assays, we show that FOXP2 transcriptional activity is regulated by the Wnt coregulator ß-catenin and that ß-catenin contacts multiple regions within FOXP2. Using nuclear magnetic resonance spectroscopy, we uncovered the molecular details of these interactions. ß-catenin contacts a disordered FOXP2 region with α-helical propensity via its folded armadillo domain, whereas the intrinsically disordered ß-catenin N terminus and C terminus bind to the conserved FOXP2 DNA-binding domain. Using RNA sequencing, we confirmed that ß-catenin indeed regulates transcriptional activity of FOXP2 and that the FOXP2 α-helical motif acts as a key regulatory element of FOXP2 transcriptional activity. Taken together, our findings provide first insight into novel regulatory interactions and help to understand the intricate mechanisms of FOXP2 function and (mis)-regulation in embryonal development and human diseases. DATABASE: Expression data are available in the GEO database under the accession number GSE138938.

Publisher Correction: Inference of transcription factor binding from cell-free DNA enables tumor subtype prediction and early detection.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Cell-free DNA analysis reveals POLR1D-mediated resistance to bevacizumab in colorectal cancer.

BACKGROUND: Bevacizumab, a monoclonal antibody against soluble VEGFA, is an approved and commonly administered anti-angiogenic drug in patients with metastasized colorectal cancer (mCRC). The survival benefit of anti-VEGF therapy in mCRC patients is limited to a few months, and acquired resistance mechanisms are largely unknown. Here, we employed whole-genome sequencing of plasma DNA to evaluate the tumor genome of patients undergoing treatment with bevacizumab to determine novel aberrations associated with resistance. METHODS: Using longitudinal plasma analyses, we studied the evolution of tumor genomes in a mCRC cohort (n = 150) and conducted analyses of CRC cases from The Cancer Genome Atlas (TCGA) database (n = 619) to identify associations between genomic aberrations and clinical features. We employed whole-genome sequencing to identify the most frequently occurring focal somatic copy number alterations (SCNAs). Using the TCGA data as a comparative and supporting dataset, we defined the minimally amplified overlapping region and studied the mechanistic consequences of copy number gain of the involved genes in this segment. In addition, we established an in vitro cell model and conducted downstream gene expression and cell viability assays to confirm our findings from the patient dataset. RESULTS: We observed a recurrent focal amplification (8.7% of cases) on chromosome 13q12.2. Analysis of CRC cases from the TCGA database suggested that this amplicon is associated with more advanced stages. We confirmed that this 13q12.2 amplicon frequently emerges later during the clinical course of disease. After defining the minimally amplified region, we observed that the amplification and expression of one gene, POLR1D, impacted cell proliferation and resulted in upregulation of VEGFA, an important regulator of angiogenesis which has been implicated in the resistance to bevacizumab treatment. In fact, in several patients, we observed the emergence of this 13q12.2 amplicon under bevacizumab treatment, which was invariably associated with therapy resistance. CONCLUSIONS: Non-invasive analyses of cell-free DNA from patients undergoing treatment with bevacizumab enabled the tracking of evolving tumor genomes and helped identify a recurrent focal SCNA of clinical relevance. Here, we describe a novel resistance mechanism against a widely applied treatment in patients with mCRC which will impact the clinical management of patients.

Inference of transcription factor binding from cell-free DNA enables tumor subtype prediction and early detection.

Deregulation of transcription factors (TFs) is an important driver of tumorigenesis, but non-invasive assays for assessing transcription factor activity are lacking. Here we develop and validate a minimally invasive method for assessing TF activity based on cell-free DNA sequencing and nucleosome footprint analysis. We analyze whole genome sequencing data for >1,000 cell-free DNA samples from cancer patients and healthy controls using a bioinformatics pipeline developed by us that infers accessibility of TF binding sites from cell-free DNA fragmentation patterns. We observe patient-specific as well as tumor-specific patterns, including accurate prediction of tumor subtypes in prostate cancer, with important clinical implications for the management of patients. Furthermore, we show that cell-free DNA TF profiling is capable of detection of early-stage colorectal carcinomas. Our approach for mapping tumor-specific transcription factor binding in vivo based on blood samples makes a key part of the noncoding genome amenable to clinical analysis.

Genome-Wide Analysis of the Nucleosome Landscape in Individuals with Coffin-Siris Syndrome.

The switch/sucrose non-fermenting (SWI/SNF) complex is an ATP-dependent chromatin remodeller that regulates the spacing of nucleosomes and thereby controls gene expression. Heterozygous mutations in genes encoding subunits of the SWI/SNF complex have been reported in individuals with Coffin-Siris syndrome (CSS), with the majority of the mutations in ARID1B. CSS is a rare congenital disorder characterized by facial dysmorphisms, digital anomalies, and variable intellectual disability. We hypothesized that mutations in genes encoding subunits of the ubiquitously expressed SWI/SNF complex may lead to alterations of the nucleosome profiles in different cell types. We performed the first study on CSS-patient samples and investigated the nucleosome landscapes of cell-free DNA (cfDNA) isolated from blood plasma by whole-genome sequencing. In addition, we studied the nucleosome landscapes of CD14+ monocytes from CSS-affected individuals by nucleosome occupancy and methylome-sequencing (NOMe-seq) as well as their expression profiles. In cfDNA of CSS-affected individuals with heterozygous ARID1B mutations, we did not observe major changes in the nucleosome profile around transcription start sites. In CD14+ monocytes, we found few genomic regions with different nucleosome occupancy when compared to controls. RNA-seq analysis of CD14+ monocytes of these individuals detected only few differentially expressed genes, which were not in proximity to any of the identified differential nucleosome-depleted regions. In conclusion, we show that heterozygous mutations in the human SWI/SNF subunit ARID1B do not have a major impact on the nucleosome landscape or gene expression in blood cells. This might be due to functional redundancy, cell-type specificity, or alternative functions of ARID1B.

Untargeted Assessment of Tumor Fractions in Plasma for Monitoring and Prognostication from Metastatic Breast Cancer Patients Undergoing Systemic Treatment.

The aim of this study was to assess the prognostic and predictive value of an untargeted assessment of tumor fractions in the plasma of metastatic breast cancer patients and to compare circulating tumor DNA (ctDNA) with circulating tumor cells (CTC) and conventional tumor markers. In metastatic breast cancer patients (n = 29), tumor fractions in plasma were assessed using the untargeted mFAST-SeqS method from 127 serial blood samples. Resulting z-scores for the ctDNA were compared to tumor fractions established with the recently published ichorCNA algorithm and associated with the clinical outcome. We observed a close correlation between mFAST-SeqS z-scores and ichorCNA ctDNA quantifications. Patients with mFAST-SeqS z-scores above three (34.5%) showed significantly worse overall survival (p = 0.014) and progression-free survival (p = 0.018) compared to patients with lower values. Elevated z-score values were clearly associated with radiologically proven progression. The baseline CTC count, carcinoembryonic antigen (CEA), and cancer antigen (CA)15-5 had no prognostic impact on the outcome of patients in the analyzed cohort. This proof of principle study demonstrates the prognostic impact of ctDNA levels detected with mFAST-SeqS as a very fast and cost-effective means to assess the ctDNA fraction without prior knowledge of the genetic landscape of the tumor. Furthermore, mFAST-SeqS-based ctDNA levels provided an early means of measuring treatment response.

Genomic alterations in plasma DNA from patients with metastasized prostate cancer receiving abiraterone or enzalutamide.

In patients with metastatic castrate-resistant prostate cancer (mCRPC), circulating tumor DNA (ctDNA) analysis offers novel opportunities for the development of non-invasive biomarkers informative of treatment response with novel agents targeting the androgen-receptor (AR) pathway, such as abiraterone or enzalutamide. However, the relationship between ctDNA abundance, detectable somatic genomic alterations and clinical progression of mCRPC remains unexplored. Our study aimed to investigate changes in plasma DNA during disease progression and their associations with clinical variables in mCRPC patients. We analyzed ctDNA in two cohorts including 94 plasma samples from 25 treatment courses (23 patients) and 334 plasma samples from 125 patients, respectively. We conducted whole-genome sequencing (plasma-Seq) for genome-wide profiling of somatic copy number alterations and targeted sequencing of 31 prostate cancer-associated genes. The combination of plasma-Seq with targeted AR analyses identified prostate cancer-related genomic alterations in 16 of 25 (64%) treatment courses in the first cohort, in which we demonstrated that AR amplification does not always correlate with poor abiraterone and enzalutamide therapy outcome. As we observed a wide variability of ctDNA levels, we evaluated ctDNA levels and their association with clinical parameters and included the second, larger cohort for these analyses. Employing altogether 428 longitudinal plasma samples from 148 patients, we identified the presence of bone metastases, increased lactate dehydrogenase and prostate-specific antigen (PSA) as having the strongest association with high ctDNA levels. In summary, ctDNA alterations are observable in the majority of patients with mCRPC and may eventually be useful to guide clinical decision-making in this setting.

Expanded molecular profiling of myxofibrosarcoma reveals potentially actionable targets.

Myxofibrosarcomas are morphologically heterogeneous soft tissue sarcomas lacking a specific immunohistochemical expression profile and recurrent genetic changes. The study was designed to gain further insights into the molecular landscape of myxofibrosarcomas by targeted re-sequencing of known cancer driver hotspot mutations and the analysis of genomewide somatic copy number alterations. A well-defined group of myxofibrosarcomas, including myxofibrosarcomas G1 (n=6), myxofibrosarcomas G3 (n=7), myxofibrosarcomas with morphologically heterogeneous and independently selectable G1 and G3 areas within a tumor (n=8), and myxofibrosarcomas G3 with subsequent tumor recurrence (n=1) or metastatic disease (n=3) were evaluated. Mutational analysis demonstrated mutations in TP53, PTEN, FGFR3, CDKN2A, and RB1. TP53 mutations were seen in 11 (44%) of patients and detected in myxofibrosarcomas G1, G3, with heterogeneous morphology and G3 with subsequent metastases in 1 patient (16%), 3 patients (42%), 2 patients (62.5%), and 3 patients (75%), respectively. Additional mutations were detected in 2 patients, intratumoral mutational heterogeneity in 1 patient. We observed a variety of copy number alterations typical for myxofibrosarcomas, with higher numbers in G3 compared with G1 myxofibrosarcomas. Cluster analysis revealed distinctive features especially in metastatic and recurrent disease. Focal alterations affected CDKN2A, CCND1, CCNE1, EGFR, EPHA3, EPHB1, FGFR1, JUN, NF1, RB1, RET, TP53, and additional novel amplifications in CCNE1, KIT, EGFR, RET, BRAF, NTRK2 were seen in G3 compared with the G1 tumor areas. The total number of focal events in G1 versus G3 tumors differed significantly (P=0.0014). TRIO and RICTOR co-amplification was seen in 8 (44%) G3 and 1 (10%) G1 myxofibrosarcomas and RICTOR amplification alone in 4 (40%) G1 myxofibrosarcomas. TRIO amplification was significantly (P=0.0218) higher in G3 myxofibrosarcomas indicating a late genetic event. These findings support the use of expanded molecular profiling in myxofibrosarcomas to detect drug-able targets to allow patients to participate in basket trials.

Patient monitoring through liquid biopsies using circulating tumor DNA.

Tumors release components such as circulating tumor cells (CTCs), circulating tumor DNA (ctDNA), and tumor-derived extracellular vesicles into the circulation. Multiple studies have demonstrated that molecular information about tumors and metastases can be extracted from these factors, which are therefore frequently referred to as "liquid biopsies." Liquid biopsies allow the longitudinal monitoring of tumor genomes non-invasively and may hence ensure that patients receive appropriate treatments that target the molecular features of their disease. Accordingly, the number of studies employing liquid biopsy based assays has been skyrocketing in the last few years. Here, we focus on three important issues, which are of high relevance for monitoring tumor genomes. First, we analyze the relation between the allele frequency of somatic tumor-specific mutations and the tumor fraction within plasma DNA. Second, we ask how well current tumor evolution models correlate with findings in longitudinal liquid biopsy studies. And, finally, as sensitivity is one of the key challenges of mutation detection, we address the challenge of detecting mutations occurring at very low allele frequencies in plasma DNA.

Autosomal Recessive Keratoderma-Ichthyosis-Deafness (ARKID) Syndrome Is Caused by VPS33B Mutations Affecting Rab Protein Interaction and Collagen Modification.

In this paper, we report three patients with severe palmoplantar keratoderma associated with ichthyosis and sensorineural deafness. Biallelic mutations were found in VPS33B, encoding VPS33B, a Sec1/Munc18 family protein that interacts with Rab11a and Rab25 proteins and is involved in trafficking of the collagen-modifying enzyme LH3. Two patients were homozygous for the missense variant p.Gly131Glu, whereas one patient was compound heterozygous for p.Gly131Glu and the splice site mutation c.240-1G>C, previously reported in patients with arthrogryposis renal dysfunction and cholestasis syndrome. We demonstrated the pathogenicity of variant p.Gly131Glu by assessing the interactions of the mutant VPS33B construct and its ability to traffic LH3. Compared with wild-type VPS33B, the p.Gly131Glu mutant VPS33B had reduced coimmunoprecipitation and colocalization with Rab11a and Rab25 and did not rescue LH3 trafficking. Confirming the cell-based experiments, we found deficient LH3-specific collagen lysine modifications in patients' urine and skin fibroblasts. Additionally, the epidermal ultrastructure of the p.Gly131Glu patients mirrored defects in tamoxifen-inducible VPS33B-deficient Vps33bfl/fl-ERT2 mice. Both patients and murine models revealed an impaired epidermal structure, ascribed to aberrant secretion of lamellar bodies, which are essential for epidermal barrier formation. Our results demonstrate that p.Gly131Glu mutant VPS33B causes an autosomal recessive keratoderma-ichthyosis-deafness syndrome.

mFast-SeqS as a Monitoring and Pre-screening Tool for Tumor-Specific Aneuploidy in Plasma DNA.

Recent progress in the analysis of cell-free DNA fragments (cell-free circulating tumor DNA, ctDNA) now allows monitoring of tumor genomes by non-invasive means. However, previous studies with plasma DNA from patients with cancer demonstrated highly variable allele frequencies of ctDNA. Comprehensive genome-wide analysis of tumor genomes is greatly facilitated when plasma DNA has increased amounts of ctDNA. In order to develop a fast and cost-effective pre-screening method for the identification of plasma samples suitable for further extensive qualitative analysis, we adapted the recently described FAST-SeqS method. We show that our modified FAST-SeqS method (mFAST-SeqS) can be used as a pre-screening tool for an estimation of the ctDNA percentage. Moreover, since the genome-wide mFAST-SeqS z-scores correlate with the actual tumor content in plasma samples, changes in ctDNA levels associated with response to treatment can be easily monitored without prior knowledge of the genetic composition of tumor samples.

Inferring expressed genes by whole-genome sequencing of plasma DNA.

The analysis of cell-free DNA (cfDNA) in plasma represents a rapidly advancing field in medicine. cfDNA consists predominantly of nucleosome-protected DNA shed into the bloodstream by cells undergoing apoptosis. We performed whole-genome sequencing of plasma DNA and identified two discrete regions at transcription start sites (TSSs) where nucleosome occupancy results in different read depth coverage patterns for expressed and silent genes. By employing machine learning for gene classification, we found that the plasma DNA read depth patterns from healthy donors reflected the expression signature of hematopoietic cells. In patients with cancer having metastatic disease, we were able to classify expressed cancer driver genes in regions with somatic copy number gains with high accuracy. We were able to determine the expressed isoform of genes with several TSSs, as confirmed by RNA-seq analysis of the matching primary tumor. Our analyses provide functional information about cells releasing their DNA into the circulation.

Exploring chromosomal abnormalities and genetic changes in uterine smooth muscle tumors.

Smooth muscle tumors of the uterus are a diagnostically challenging group of tumors. Molecular surrogate markers reliably distinguishing between benign and malignant tumors are not available. Therefore, the diagnosis is based on morphologic criteria. The aim was to investigate a well-characterized group of challenging uterine smooth muscle tumors consisting of 20 leiomyomas, 13 leiomyomas with bizarre nuclei, and 14 leiomyosarcomas for copy number alterations, MED12 mutations and FH deletions to search for potential diagnostically useful surrogate markers. MED12 mutations were detected in 47, 15, and 25% of leiomyomas, leiomyomas with bizarre nuclei and leiomyosarcomas, respectively. MED12 mutations in leiomyomas with bizarre nuclei were detected outside the hotspot region. FH-deletions were seen in 27, 30.8, and 25% of leiomyomas, leiomyomas with bizarre nuclei and leiomyosarcomas, respectively. By using copy number alteration profiling a clear separation of leiomyomas, leiomyomas with bizarre nuclei and leiomyosarcomas could not be observed. Copy number alterations revealed clear genetic similarities between leiomyomas with bizarre nuclei and leiomyosarcomas. Leiomyosarcomas showed a similar pattern of gains and losses as leiomyomas with bizarre nuclei, with additional copy number alterations and more homozygous losses and high-level amplifications compared to leiomyomas with bizarre nuclei. In conclusion, this study demonstrates that known FH-deletions, a recurrent molecular change in leiomyomas, occur in morphologically challenging variants of leiomyomas, leiomyomas with bizarre nuclei and leiomyosarcomas. Although MED12 mutations are common in leiomyomas, they infrequently occur in leiomyomas with bizarre nuclei and leiomyosarcomas. The genetic similarities between leiomyomas with bizarre nuclei and leiomyosarcomas raise the intriguing possibility that uterine leiomyomas with bizarre nuclei and leiomyosarcomas are closely related and challenge the traditional concept that leiomyoma with bizarre nuclei is a tumor with just marked 'degenerative' cellular changes. These findings support the hypothesis that tumor progression within uterine smooth muscle tumors might occur.

Whole-genome plasma sequencing reveals focal amplifications as a driving force in metastatic prostate cancer.

Genomic alterations in metastatic prostate cancer remain incompletely characterized. Here we analyse 493 prostate cancer cases from the TCGA database and perform whole-genome plasma sequencing on 95 plasma samples derived from 43 patients with metastatic prostate cancer. From these samples, we identify established driver aberrations in a cancer-related gene in nearly all cases (97.7%), including driver gene fusions (TMPRSS2:ERG), driver focal deletions (PTEN, RYBP and SHQ1) and driver amplifications (AR and MYC). In serial plasma analyses, we observe changes in focal amplifications in 40% of cases. The mean time interval between new amplifications was 26.4 weeks (range: 5-52 weeks), suggesting that they represent rapid adaptations to selection pressure. An increase in neuron-specific enolase is accompanied by clonal pattern changes in the tumour genome, most consistent with subclonal diversification of the tumour. Our findings suggest a high plasticity of prostate cancer genomes with newly occurring focal amplifications as a driving force in progression.

Non-invasive detection of genome-wide somatic copy number alterations by liquid biopsies.

Liquid biopsies, i.e. the analysis of circulating tumor cells (CTCs) or circulating tumor DNA (ctDNA), are evolving into promising clinical tools. Indeed, a plethora of liquid biopsy technologies to deduce non-invasively characteristics of the tumor genome from the peripheral blood have been developed over the last few years. For example, liquid biopsies have been used to assess the tumor burden, to monitor the evolution of tumor genomes, to unravel mechanisms of resistance, to establish the tumor heterogeneity, and for the identification of prognostic and predictive markers. In this review we focus on methods to establish genome-wide profiles of somatic copy number alterations (SCNAs) from plasma DNA and show how they provide novel insights into the biology of cancer and their impact on the management of patients.