A single-array preprocessing method for estimating full-resolution raw copy numbers from all Affymetrix genotyping arrays including GenomeWideSNP 5 & 6

Multi-omics profiling reveals key factors involved in Ewing sarcoma metastasis

Ewing sarcoma (EWS) is the second most common bone tumor affecting children and young adults, with dismal outcomes for patients with metastasis at diagnosis. Mechanisms leading to metastasis remain poorly understood. To deepen our knowledge on EWS progression, we have profiled tumors and metastases from a spontaneous metastasis mouse model using a multi-omics approach. Combining transcriptomics, proteomics, and methylomics analyses, we identified signaling cascades and candidate genes enriched in metastases that could be modulating aggressiveness in EWS. Phenotypical validation of two of these candidates, cyclic AMP-responsive element-binding protein 1 (CREB1) and lipoxygenase homology domain-containing protein 1 (LOXHD1), showed an association with migration and clonogenic abilities. Moreover, previously described CREB1 downstream targets were present amongst the metastatic-enriched results. Combining the different omics datasets, we identified FYVE, RhoGEF, and PH domain-containing protein 4 (FGD4) as a CREB1 target interconnecting the different EWS biological layers (RNA, protein and methylation status) and whose high expression is associated with worse clinical outcome. Further studies will provide insight into EWS metastasis mechanisms and ultimately improve survival rates for EWS patients.

Frequent copy number variants in a cohort of Mexican-Mestizo individuals

BACKGROUND

The human genome presents variation at distinct levels, copy number variants (CNVs) are DNA segments of variable lengths that range from several base pairs to megabases and are present at a variable number of copies in human genomes. Common CNVs have no apparent influence on the phenotype; however, some rare CNVs have been associated with phenotypic traits, depending on their size and gene content. CNVs are detected by microarrays of different densities and are generally visualized, and their frequencies analysed using the HapMap as default reference population. Nevertheless, this default reference is inadequate when the samples analysed are from people from Mexico, since population with a Hispanic genetic background are minimally represented. In this work, we describe the variation in the frequencies of four common CNVs in Mexican-Mestizo individuals.

RESULTS

In a cohort of 147 unrelated Mexican-Mestizo individuals, we found that the common CNVs 2p11.2 (99.6%), 8p11.22 (54.5%), 14q32.33 (100%), and 15q11.2 (71.1%) appeared with unexpectedly high frequencies when contrasted with the HapMap reference (ChAS). Yet, while when comparing to an ethnically related reference population, these differences were significantly reduced or even disappeared.

CONCLUSION

The findings in this work contribute to (1) a better description of the CNVs characteristics of the Mexican Mestizo population and enhance the knowledge of genome variation in different ethnic groups. (2) emphasize the importance of contrasting CNVs identified in studied individuals against a reference group that-as best as possible-share the same ethnicity while keeping this relevant information in mind when conducting CNV studies at the population or clinical level.

A family study implicates GBE1 in the etiology of autism spectrum disorder

Autism spectrum disorders (ASD) are neurodevelopmental disorders with an estimated heritability of >60%. Family-based genetic studies of ASD have generally focused on multiple small kindreds, searching for de novo variants of major effect. We hypothesized that molecular genetic analysis of large multiplex families would enable the identification of variants of milder effects. We studied a large multigenerational family of European ancestry with multiple family members affected with ASD or the broader autism phenotype (BAP). We identified a rare heterozygous variant in the gene encoding 1,4-ɑ-glucan branching enzyme 1 (GBE1) that was present in seven of seven individuals with ASD, nine of ten individuals with the BAP, and none of four tested unaffected individuals. We genotyped a community-acquired cohort of 389 individuals with ASD and identified three additional probands. Cascade analysis demonstrated that the variant was present in 11 of 13 individuals with familial ASD/BAP and neither of the two tested unaffected individuals in these three families, also of European ancestry. The variant was not enriched in the combined UK10K ASD cohorts of European ancestry but heterozygous GBE1 deletion was overrepresented in large ASD cohorts, collectively suggesting an association between GBE1 and ASD.

Deletions in VANGL1 are a risk factor for antibody-mediated kidney disease

We identify an intronic deletion in VANGL1 that predisposes to renal injury in high risk populations through a kidney-intrinsic process. Half of all SLE patients develop nephritis, yet the predisposing mechanisms to kidney damage remain poorly understood. There is limited evidence of genetic contribution to specific organ involvement in SLE.1,2 We identify a large deletion in intron 7 of Van Gogh Like 1 (VANGL1), which associates with nephritis in SLE patients. The same deletion occurs at increased frequency in an indigenous population (Tiwi Islanders) with 10-fold higher rates of kidney disease compared with non-indigenous populations. Vangl1 hemizygosity in mice results in spontaneous IgA and IgG deposition within the glomerular mesangium in the absence of autoimmune nephritis. Serum transfer into B cell-deficient Vangl1+/- mice results in mesangial IgG deposition indicating that Ig deposits occur in a kidney-intrinsic fashion in the absence of Vangl1. These results suggest that Vangl1 acts in the kidney to prevent Ig deposits and its deficiency may trigger nephritis in individuals with SLE.

A genome‑wide scan to locate regions associated with familial vesicoureteral reflux

Vesicoureteral reflux (VUR) is a congenital malformation carrying a high risk of recurrent urinary tract infections (UTI) and, at worst, chronic renal failure. Familial clustering implies a genetic etiology, but studies during the past few decades have demonstrated a causal gene variant in <10% of patients with VUR. The aim of the present study was to search for fully or partially shared ancestral haplotypes in 14 families from south-western Sweden with at least three affected members. High-density single nucleotide polymorphism microarray was used for genotyping prior to analysis with a compatibility matching method developed in-house, and the analysis of copy number variations (CNV). No single unique haplotype was revealed to be shared by the families, thereby excluding a common ancestry and founder mutations as a probable cause of VUR. After evaluation of haplotypes shared by subsets of families, a haplotype shared by nine families was found to be of particular interest. This haplotype, located at chromosomal region 4q21.21, harbours two tentative candidate genes (bone morphogenetic protein 3 and fibroblast growth factor 5), both expressed in metanephros and with known functions during nephrogenesis. As to CNV, only one family had a specific CNV shared by all affected members. This was a focal deletion at 5q31.1 including follistatin-like 4, a gene without a previous known connection to VUR. These data demonstrated the genetic heterogeneity of VUR and indicated that an interaction of environmental and genetic factors, including non-coding and epigenetic regulators, all contribute to the complexity of VUR.

Dissecting the heterogeneity of the alternative polyadenylation profiles in triple-negative breast cancers

Background: Triple-negative breast cancer (TNBC) is an aggressive malignancy with high heterogeneity. However, the alternative polyadenylation (APA) profiles of TNBC remain unknown. Here, we aimed to define the characteristics of the APA events at post-transcription level among TNBCs. Methods: Using transcriptome microarray data, we analyzed APA profiles of 165 TNBC samples and 33 paired normal tissues. A pooled short hairpin RNA screen targeting 23 core cleavage and polyadenylation (C/P) genes was used to identify key C/P factors. Results: We established an unconventional APA subtyping system composed of four stable subtypes: 1) luminal androgen receptor (LAR), 2) mesenchymal-like immune-activated (MLIA), 3) basal-like (BL), 4) suppressed (S) subtypes. Patients in the S subtype had the worst disease-free survival comparing to other patients (log-rank p = 0.021). Enriched clinically actionable pathways and putative therapeutic APA events were analyzed among each APA subtype. Furthermore, CPSF1 and PABPN1 were identified as the master C/P factors in regulating APA events and TNBC proliferation. The depletion of CPSF1 or PABPN1 weakened cell proliferation, enhanced apoptosis, resulted in cell cycle redistribution and a reversion of APA events of genes associated with tumorigenesis, proliferation, metastasis and chemosensitivity in breast cancer. Conclusions: Our findings advance the understanding of tumor heterogeneity regulation in APA and yield new insights into therapeutic target identification in TNBC.

Combinations of common SNPs of the transporter gene ABCB1 influence apparent bioavailability, but not renal elimination of oral digoxin

Effects of different genotypes on the pharmacokinetics of probe substrates may support their use as phenotyping agents for the activity of the respective enzyme or transporter. Digoxin is recommended as a probe substrate to assess the activity of the transporter P-glycoprotein (P-gp) in humans. Current studies on the individual effects of three commonly investigated single nucleotide polymorphisms (SNPs) of the ABCB1 gene encoding P-gp (C1236T, G2677T/A, and C3435T) on digoxin pharmacokinetics are inconclusive. Since SNPs are in incomplete linkage disequilibrium, considering combinations of these SNPs might be necessary to assess the role of polymorphisms in digoxin pharmacokinetics accurately. In this study, the relationship between SNP combinations and digoxin pharmacokinetics was explored via a population pharmacokinetic approach in 40 volunteers who received oral doses of 0.5 mg digoxin. Concerning the SNPs 1236/2677/3435, the following combinations were evaluated: CGC, CGT, and TTT. Carriers of CGC/CGT and TTT/TTT had 35% higher apparent bioavailability compared to the reference group CGC/CGC, while no difference was seen in CGC/TTT carriers. No significant effect on renal clearance was observed. The population pharmacokinetic model supports the use of oral digoxin as a phenotyping substrate of intestinal P-gp, but not to assess renal P-gp activity.

The Landscape of Somatic Copy Number Alterations in Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is the sixth most common malignancy worldwide. Somatic copy number alterations (CNAs) play a significant role in the development of this lethal cancer. In this study, we present a meta-analysis of CNAs for a total of 1,395 HNSCC samples. Publicly available R packages and in-house scripts were used for genomic array data processing, including normalization, segmentation and CNA calling. We detected 125 regions of significant gains or losses using GISTIC algorithm and found several potential driver genes in these regions. The incidence of chromothripsis in HNSCC was estimated to be 6%, and the chromosome pulverization hotspot regions were detected. We determined 323 genomic locations significantly enriched for breakpoints, which indicate HNSCC-specific genomic instability regions. Unsupervised clustering of genome-wide CNA data revealed a sub-cluster predominantly composed of nasopharynx tumors and presented a large proportion of HPV-positive samples. These results will facilitate the discovery of therapeutic candidates and extend our molecular understanding of HNSCC.

Genome-wide somatic copy number alteration analysis and database construction for cervical cancer

Cervical cancer is a common gynecological malignancy with high incidence and mortality. Somatic copy number alterations (CNAs) play an important role in identifying tumor suppressor genes and oncogenes and are a useful diagnostic indicator for many cancer types. However, the genomic landscape of CNAs in cervical cancer has not yet been comprehensively characterized. In the present study, we collected 974 cervical cancer samples from different data sources. All samples were analyzed by genomic arrays to obtain high-resolution CNAs. Focal genomic regions with CNA events and potential cancer driver genes were identified by GISTIC2.0. Meanwhile, we constructed a comprehensive cervical cancer database by PHP and self-written Perl and R scripts. In total, 54 recurrent regions of amplification and deletion were detected. Frequently altered tumor suppressor genes were found in these regions, including PIK3CA, ERBB2, EP300 and FBXW7. CNA hotspots and related enriched functional categories were also identified. The incidence of chromothripsis in cervical cancer was estimated to be 6.06%, and the chromosome pulverization hotspot regions were detected. Based on the curated data, we developed CNAdbCC (http://cailab.labshare.cn/CNAdbCC/), a comprehensive database for copy number alterations in cervical cancer. We provide a user-friendly Web interface for data mining and visualization. It is the most comprehensive public database devoted exclusively to genomic alterations in cervical cancer. These results extend our molecular understanding of cervical cancer. The database will enable researchers to explore specific CNA patterns in this lethal cancer and facilitate the discovery of therapeutic candidates.

Comprehensive identification and characterization of somatic copy number alterations in triple‑negative breast cancer

Triple‑negative breast cancer (TNBC) accounts for ~15% of all breast cancer diagnoses each year. Patients with TNBC tend to have a higher risk for early relapse and a worse prognosis. TNBC is characterized by extensive somatic copy number alterations (CNAs). However, the DNA CNA profile of TNBC remains to be extensively investigated. The present study assessed the genomic profile of CNAs in 201 TNBC samples, aiming to identify recurrent CNAs that may drive the pathogenesis of TNBC. In total, 123 regions of significant amplification and deletion were detected using the Genomic Identification of Significant Targets in Cancer algorithm, and potential driver genes for TNBC were identified. A total of 31 samples exhibited signs of chromothripsis and revealed chromosome pulverization hotspot regions. The present study further determined 199 genomic locations that were significantly enriched for breakpoints, which indicated TNBC‑specific genomic instability regions. Unsupervised hierarchical clustering of tumors resulted in three main subgroups that exhibited distinct CNA profiles, which may reveal the heterogeneity of molecular mechanisms in TNBC subgroups. These results will extend the molecular understanding of TNBC and will facilitate the discovery of therapeutic and diagnostic target candidates.

A high-resolution 3D epigenomic map reveals insights into the creation of the prostate cancer transcriptome

To better understand the impact of chromatin structure on regulation of the prostate cancer transcriptome, we develop high-resolution chromatin interaction maps in normal and prostate cancer cells using in situ Hi-C. By combining the in situ Hi-C data with active and repressive histone marks, CTCF binding sites, nucleosome-depleted regions, and transcriptome profiling, we identify topologically associating domains (TADs) that change in size and epigenetic states between normal and prostate cancer cells. Moreover, we identify normal and prostate cancer-specific enhancer-promoter loops and involved transcription factors. For example, we show that FOXA1 is enriched in prostate cancer-specific enhancer-promoter loop anchors. We also find that the chromatin structure surrounding the androgen receptor (AR) locus is altered in the prostate cancer cells with many cancer-specific enhancer-promoter loops. This creation of 3D epigenomic maps enables a better understanding of prostate cancer biology and mechanisms of gene regulation.

In prostate cancer, chromatin structure can impact the transcriptome. Here, the authors develop high resolution chromatin interaction maps in prostate cancer cells using in situ Hi-C, revealing prostate cancer-specific TADs and enhancer-promoter loops surrounding the androgen receptor (AR) locus.

A Clinical Drug-Drug Interaction Study Assessing a Novel Drug Transporter Phenotyping Cocktail With Adefovir, Sitagliptin, Metformin, Pitavastatin, and Digoxin

A new probe drug cocktail containing substrates of important drug transporters was tested for mutual interactions in a clinical trial. The cocktail consisted of (predominant transporter; primary phenotyping metric): 10 mg adefovir-dipivoxil (OAT1; renal clearance (CL_R )), 100 mg sitagliptin (OAT3; CL_R ), 500 mg metformin (several renal transporters; CL_R ), 2 mg pitavastatin (OATP1B1; clearance/F), and 0.5 mg digoxin (intestinal P-gp, renal P-gp, and OATP4C1; peak plasma concentration (C_max ) and CL_R ). Using a randomized six-period, open change-over design, single oral doses were administrated either concomitantly or separately to 24 healthy male and female volunteers. Phenotyping metrics were evaluated by noncompartmental analysis and compared between periods by the standard average bioequivalence approach (boundaries for ratios 0.80-1.25). Primary metrics supported the absence of relevant interactions, whereas secondary metrics suggested that mainly adefovir was a victim of minor drug-drug interactions (DDIs). All drugs were well tolerated. This cocktail may be another useful tool to assess transporter-based DDIs in vivo.

Differential Expression of DNA Repair Genes in Prognostically-Favorable versus Unfavorable Uveal Melanoma

Expression of DNA repair genes was studied in uveal melanoma (UM) in order to identify genes that may play a role in metastases formation. We searched for genes that are differentially expressed between tumors with a favorable and unfavorable prognosis. Gene-expression profiling was performed on 64 primary UM from the Leiden University Medical Center (LUMC), Leiden, The Netherlands. The expression of 121 genes encoding proteins involved in DNA repair pathways was analyzed: a total of 44 genes differed between disomy 3 and monosomy 3 tumors. Results were validated in a cohort from Genoa and Paris and the The Cancer Genome Atlas (TCGA) cohort. Expression of the PRKDC, WDR48, XPC, and BAP1 genes was significantly associated with clinical outcome after validation. PRKDC was highly expressed in metastasizing UM (p < 0.001), whereas WDR48, XPC, and BAP1 were lowly expressed (p < 0.001, p = 0.006, p = 0.003, respectively). Low expression of WDR48 and XPC was related to a large tumor diameter (p = 0.01 and p = 0.004, respectively), and a mixed/epithelioid cell type (p = 0.007 and p = 0.03, respectively). We conclude that the expression of WDR48, XPC, and BAP1 is significantly lower in UM with an unfavorable prognosis, while these tumors have a significantly higher expression of PRKDC. Pharmacological inhibition of DNA-PKcs resulted in decreased survival of UM cells. PRKDC may be involved in proliferation, invasion and metastasis of UM cells. Unraveling the role of DNA repair genes may enhance our understanding of UM biology and result in the identification of new therapeutic targets.

Contribution of epigenetic variation to adaptation in Arabidopsis

In plants, transgenerational inheritance of some epialleles has been demonstrated but it remains controversial whether epigenetic variation is subject to selection and contributes to adaptation. Simulating selection in a rapidly changing environment, we compare phenotypic traits and epigenetic variation between Arabidopsis thaliana populations grown for five generations under selection and their genetically nearly identical ancestors. Selected populations of two distinct genotypes show significant differences in flowering time and plant architecture, which are maintained for at least 2–3 generations in the absence of selection. While we cannot detect consistent genetic changes, we observe a reduction of epigenetic diversity and changes in the methylation state of about 50,000 cytosines, some of which are associated with phenotypic changes. Thus, we propose that epigenetic variation is subject to selection and can contribute to rapid adaptive responses, although the extent to which epigenetics plays a role in adaptation is still unclear.

Whether plant epigenetic variation is subject to selection and contributes to adaptation is under debate. Here, the authors compare DNA methylation and phenotypes of Arabidopsis lines subject to simulated selection and their nearly isogenic ancestors and provide evidence that epigenetic variation contributes to adaptive responses.

Fatty acid oxidation is associated with proliferation and prognosis in breast and other cancers

Background

Altered cellular metabolism is a hallmark of cancer but the association between utilisation of particular metabolic pathways in tumours and patient outcome is poorly understood. We sought to investigate the association between fatty acid metabolism and outcome in breast and other cancers.

Methods

Cox regression analysis and Gene Set Enrichment Analysis (GSEA) of a gene expression dataset from primary breast tumours with well annotated clinical and survival information was used to identify genesets associated with outcome. A geneset representing fatty acid oxidation (FAO) was then examined in other datasets. A doxycycline-inducible breast cancer cell line model overexpressing the rate-limiting enzyme in FAO, carnitine palmitoyl transferase 1A (CPT1A) was generated and analysed to confirm the association between FAO and cancer-associated characteristics in vitro.

Results

We identified a gene expression signature composed of 19 genes associated with fatty acid oxidation (FAO) that was significantly associated with patient outcome. We validated this observation in eight independent breast cancer datasets, and also observed the FAO signature to be prognostic in other cancer types. Furthermore, the FAO signature expression was significantly downregulated in tumours, compared to normal tissues from a variety of anatomic origins. In breast cancer, the expression of CPT1A was higher in oestrogen receptor (ER)-positive, compared to ER-negative tumours and cell lines. Importantly, overexpression of CPT1A significantly decreased the proliferation and wound healing migration rates of MDA-MB231 breast cancer cells, compared to basal expression control.

Conclusions

Our findings suggest that FAO is downregulated in multiple tumour types, and activation of this pathway may lower cancer cell proliferation, and is associated with improved outcomes in some cancers.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-4626-9) contains supplementary material, which is available to authorized users.

Integrated genomics and functional validation identifies malignant cell specific dependencies in triple negative breast cancer

Triple negative breast cancers (TNBCs) lack recurrent targetable driver mutations but demonstrate frequent copy number aberrations (CNAs). Here, we describe an integrative genomic and RNAi-based approach that identifies and validates gene addictions in TNBCs. CNAs and gene expression alterations are integrated and genes scored for pre-specified target features revealing 130 candidate genes. We test functional dependence on each of these genes using RNAi in breast cancer and non-malignant cells, validating malignant cell selective dependence upon 37 of 130 genes. Further analysis reveals a cluster of 13 TNBC addiction genes frequently co-upregulated that includes genes regulating cell cycle checkpoints, DNA damage response, and malignant cell selective mitotic genes. We validate the mechanism of addiction to a potential drug target: the mitotic kinesin family member C1 (KIFC1/HSET), essential for successful bipolar division of centrosome-amplified malignant cells and develop a potential selection biomarker to identify patients with tumors exhibiting centrosome amplification.

Triple Angiokinase Inhibitor Nintedanib Directly Inhibits Tumor Cell Growth and Induces Tumor Shrinkage via Blocking Oncogenic Receptor Tyrosine Kinases

The triple-angiokinase inhibitor nintedanib is an orally available, potent, and selective inhibitor of tumor angiogenesis by blocking the tyrosine kinase activities of vascular endothelial growth factor receptor (VEGFR) 1-3, platelet-derived growth factor receptor (PDGFR)-α and -β, and fibroblast growth factor receptor (FGFR) 1-3. Nintedanib has received regulatory approval as second-line treatment of adenocarcinoma non-small cell lung cancer (NSCLC), in combination with docetaxel. In addition, nintedanib has been approved for the treatment of idiopathic lung fibrosis. Here we report the results from a broad kinase screen that identified additional kinases as targets for nintedanib in the low nanomolar range. Several of these kinases are known to be mutated or overexpressed and are involved in tumor development (discoidin domain receptor family, member 1 and 2, tropomyosin receptor kinase A (TRKA) and C, rearranged during transfection proto-oncogene [RET proto oncogene]), as well as in fibrotic diseases (e.g., DDRs). In tumor cell lines displaying molecular alterations in potential nintedanib targets, the inhibitor demonstrates direct antiproliferative effects: in the NSCLC cell line NCI-H1703 carrying a PDGFRα amplification (ampl.); the gastric cancer cell line KatoIII and the breast cancer cell line MFM223, both driven by a FGFR2 amplification; AN3CA (endometrial carcinoma) bearing a mutated FGFR2; the acute myeloid leukemia cell lines MOLM-13 and MV-4-11-B with FLT3 mutations; and the NSCLC adenocarcinoma LC-2/ad harboring a CCDC6-RET fusion. Potent kinase inhibition does not, however, strictly translate into antiproliferative activity, as demonstrated in the TRKA-dependent cell lines CUTO-3 and KM-12. Importantly, nintedanib treatment of NCI-H1703 tumor xenografts triggered effective tumor shrinkage, indicating a direct effect on the tumor cells in addition to the antiangiogenic effect on the tumor stroma. These findings will be instructive in guiding future genome-based clinical trials of nintedanib.

ITGB1-dependent upregulation of Caveolin-1 switches TGFβ signalling from tumour-suppressive to oncogenic in prostate cancer

Caveolin-1 (CAV1) is over-expressed in prostate cancer (PCa) and is associated with adverse prognosis, but the molecular mechanisms linking CAV1 expression to disease progression are poorly understood. Extensive gene expression correlation analysis, quantitative multiplex imaging of clinical samples, and analysis of the CAV1-dependent transcriptome, supported that CAV1 re-programmes TGFβ signalling from tumour suppressive to oncogenic (i.e. induction of SLUG, PAI-1 and suppression of CDH1, DSP, CDKN1A). Supporting such a role, CAV1 knockdown led to growth arrest and inhibition of cell invasion in prostate cancer cell lines. Rationalized RNAi screening and high-content microscopy in search for CAV1 upstream regulators revealed integrin beta1 (ITGB1) and integrin associated proteins as CAV1 regulators. Our work suggests TGFβ signalling and beta1 integrins as potential therapeutic targets in PCa over-expressing CAV1, and contributes to better understand the paradoxical dual role of TGFβ in tumour biology.

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.

The cancer immunoediting hypothesis assumes the immune system sculpts the cancer genome. Here the authors show, in a mouse model, that neutral evolution outweighs the effects of immunoselection and that immune checkpoint blockade potentiates the immunoediting, switching the system to non-neutral evolution.

A Network-Based Integrative Workflow to Unravel Mechanisms Underlying Disease Progression

Unraveling mechanisms underlying diseases has motivated the development of systems biology approaches. The key challenges for the development of mathematical models and computational tool are (1) the size of molecular networks, (2) the nonlinear nature of spatio-temporal interactions, and (3) feedback loops in the structure of interaction networks. We here propose an integrative workflow that combines structural analyses of networks, high-throughput data, and mechanistic modeling. As an illustration of the workflow, we use prostate cancer as a case study with the aim of identifying key functional components associated with primary to metastasis transitions. Analysis carried out by the workflow revealed that HOXD10, BCL2, and PGR are the most important factors affected in primary prostate samples, whereas, in the metastatic state, STAT3, JUN, and JUNB are playing a central role. The identified key elements of each network are validated using patient survival analysis. The workflow presented here allows experimentalists to use heterogeneous data sources for the identification of diagnostic and prognostic signatures.

Integrative analysis of genome-wide gene copy number changes and gene expression in non-small cell lung cancer

Non-small cell lung cancer (NSCLC) represents a genomically unstable cancer type with extensive copy number aberrations. The relationship of gene copy number alterations and subsequent mRNA levels has only fragmentarily been described. The aim of this study was to conduct a genome-wide analysis of gene copy number gains and corresponding gene expression levels in a clinically well annotated NSCLC patient cohort (n = 190) and their association with survival. While more than half of all analyzed gene copy number-gene expression pairs showed statistically significant correlations (10,296 of 18,756 genes), high correlations, with a correlation coefficient >0.7, were obtained only in a subset of 301 genes (1.6%), including KRAS, EGFR and MDM2. Higher correlation coefficients were associated with higher copy number and expression levels. Strong correlations were frequently based on few tumors with high copy number gains and correspondingly increased mRNA expression. Among the highly correlating genes, GO groups associated with posttranslational protein modifications were particularly frequent, including ubiquitination and neddylation. In a meta-analysis including 1,779 patients we found that survival associated genes were overrepresented among highly correlating genes (61 of the 301 highly correlating genes, FDR adjusted p<0.05). Among them are the chaperone CCT2, the core complex protein NUP107 and the ubiquitination and neddylation associated protein CAND1. In conclusion, in a comprehensive analysis we described a distinct set of highly correlating genes. These genes were found to be overrepresented among survival-associated genes based on gene expression in a large collection of publicly available datasets.

An indicator cell assay for blood-based diagnostics

We have established proof of principle for the Indicator Cell Assay Platform™ (iCAP™), a broadly applicable tool for blood-based diagnostics that uses specifically-selected, standardized cells as biosensors, relying on their innate ability to integrate and respond to diverse signals present in patients' blood. To develop an assay, indicator cells are exposed in vitro to serum from case or control subjects and their global differential response patterns are used to train reliable, disease classifiers based on a small number of features. In a feasibility study, the iCAP detected pre-symptomatic disease in a murine model of amyotrophic lateral sclerosis (ALS) with 94% accuracy (p-Value = 3.81E-6) and correctly identified samples from a murine Huntington's disease model as non-carriers of ALS. Beyond the mouse model, in a preliminary human disease study, the iCAP detected early stage Alzheimer's disease with 72% cross-validated accuracy (p-Value = 3.10E-3). For both assays, iCAP features were enriched for disease-related genes, supporting the assay's relevance for disease research.

Integrative clustering reveals a novel split in the luminal A subtype of breast cancer with impact on outcome

BACKGROUND

Breast cancer is a heterogeneous disease at the clinical and molecular level. In this study we integrate classifications extracted from five different molecular levels in order to identify integrated subtypes.

METHODS

Tumor tissue from 425 patients with primary breast cancer from the Oslo2 study was cut and blended, and divided into fractions for DNA, RNA and protein isolation and metabolomics, allowing the acquisition of representative and comparable molecular data. Patients were stratified into groups based on their tumor characteristics from five different molecular levels, using various clustering methods. Finally, all previously identified and newly determined subgroups were combined in a multilevel classification using a "cluster-of-clusters" approach with consensus clustering.

RESULTS

Based on DNA copy number data, tumors were categorized into three groups according to the complex arm aberration index. mRNA expression profiles divided tumors into five molecular subgroups according to PAM50 subtyping, and clustering based on microRNA expression revealed four subgroups. Reverse-phase protein array data divided tumors into five subgroups. Hierarchical clustering of tumor metabolic profiles revealed three clusters. Combining DNA copy number and mRNA expression classified tumors into seven clusters based on pathway activity levels, and tumors were classified into ten subtypes using integrative clustering. The final consensus clustering that incorporated all aforementioned subtypes revealed six major groups. Five corresponded well with the mRNA subtypes, while a sixth group resulted from a split of the luminal A subtype; these tumors belonged to distinct microRNA clusters. Gain-of-function studies using MCF-7 cells showed that microRNAs differentially expressed between the luminal A clusters were important for cancer cell survival. These microRNAs were used to validate the split in luminal A tumors in four independent breast cancer cohorts. In two cohorts the microRNAs divided tumors into subgroups with significantly different outcomes, and in another a trend was observed.

CONCLUSIONS

The six integrated subtypes identified confirm the heterogeneity of breast cancer and show that finer subdivisions of subtypes are evident. Increasing knowledge of the heterogeneity of the luminal A subtype may add pivotal information to guide therapeutic choices, evidently bringing us closer to improved treatment for this largest subgroup of breast cancer.

p120-catenin prevents multinucleation through control of MKLP1-dependent RhoA activity during cytokinesis

Spatiotemporal activation of RhoA and actomyosin contraction underpins cellular adhesion and division. Loss of cell-cell adhesion and chromosomal instability are cardinal events that drive tumour progression. Here, we show that p120-catenin (p120) not only controls cell-cell adhesion, but also acts as a critical regulator of cytokinesis. We find that p120 regulates actomyosin contractility through concomitant binding to RhoA and the centralspindlin component MKLP1, independent of cadherin association. In anaphase, p120 is enriched at the cleavage furrow where it binds MKLP1 to spatially control RhoA GTPase cycling. Binding of p120 to MKLP1 during cytokinesis depends on the N-terminal coiled-coil domain of p120 isoform 1A. Importantly, clinical data show that loss of p120 expression is a common event in breast cancer that strongly correlates with multinucleation and adverse patient survival. In summary, our study identifies p120 loss as a driver event of chromosomal instability in cancer.

Upregulation of HLA Expression in Primary Uveal Melanoma by Infiltrating Leukocytes

Introduction

Uveal melanoma (UM) with an inflammatory phenotype, characterized by infiltrating leukocytes and increased human leukocyte antigen (HLA) expression, carry an increased risk of death due to metastases. These tumors should be ideal for T-cell based therapies, yet it is not clear why prognostically-infaust tumors have a high HLA expression. We set out to determine whether the level of HLA molecules in UM is associated with other genetic factors, HLA transcriptional regulators, or microenvironmental factors.

Methods

28 enucleated UM were used to study HLA class I and II expression, and several regulators of HLA by immunohistochemistry, PCR microarray, qPCR and chromosome SNP-array. Fresh tumor samples of eight primary UM and four metastases were compared to their corresponding xenograft in SCID mice, using a PCR microarray and SNP array.

Results

Increased expression levels of HLA class I and II showed no dosage effect of chromosome 6p, but, as expected, were associated with monosomy of chromosome 3. Increased HLA class I and II protein levels were positively associated with their gene expression and with raised levels of the peptide-loading gene TAP1, and HLA transcriptional regulators IRF1, IRF8, CIITA, and NLRC5, revealing a higher transcriptional activity in prognostically-bad tumors. Implantation of fresh human tumor samples into SCID mice led to a loss of infiltrating leukocytes, and to a decreased expression of HLA class I and II genes, and their regulators.

Conclusion

Our data provides evidence for a proper functioning HLA regulatory system in UM, offering a target for T-cell based therapies.

Genomic landscape of liposarcoma

Liposarcoma (LPS) is the most common type of soft tissue sarcoma accounting for 20% of all adult sarcomas. Due to absence of clinically effective treatment options in inoperable situations and resistance to chemotherapeutics, a critical need exists to identify novel therapeutic targets. We analyzed LPS genomic landscape using SNP arrays, whole exome sequencing and targeted exome sequencing to uncover the genomic information for development of specific anti-cancer targets. SNP array analysis indicated known amplified genes (MDM2, CDK4, HMGA2) and important novel genes (UAP1, MIR557, LAMA4, CPM, IGF2, ERBB3, IGF1R). Carboxypeptidase M (CPM), recurrently amplified gene in well-differentiated/de-differentiated LPS was noted as a putative oncogene involved in the EGFR pathway. Notable deletions were found at chromosome 1p (RUNX3, ARID1A), chromosome 11q (ATM, CHEK1) and chromosome 13q14.2 (MIR15A, MIR16-1). Significantly and recurrently mutated genes (false discovery rate < 0.05) included PLEC (27%), MXRA5 (21%), FAT3 (24%), NF1 (20%), MDC1 (10%), TP53 (7%) and CHEK2 (6%). Further, in vitro and in vivo functional studies provided evidence for the tumor suppressor role for Neurofibromin 1 (NF1) gene in different subtypes of LPS. Pathway analysis of recurrent mutations demonstrated signaling through MAPK, JAK-STAT, Wnt, ErbB, axon guidance, apoptosis, DNA damage repair and cell cycle pathways were involved in liposarcomagenesis. Interestingly, we also found mutational and copy number heterogeneity within a primary LPS tumor signifying the importance of multi-region sequencing for cancer-genome guided therapy. In summary, these findings provide insight into the genomic complexity of LPS and highlight potential druggable pathways for targeted therapeutic approach.

CTLPScanner: a web server for chromothripsis-like pattern detection

Chromothripsis is a recently observed phenomenon in cancer cells in which one or several chromosomes shatter into pieces with subsequent inaccurate reassembly and clonal propagation. This type of event generates a potentially vast number of mutations within a relatively short-time period, and has been considered as a new paradigm in cancer development. Despite recent advances, much work is still required to better understand the molecular mechanisms of this phenomenon, and thus an easy-to-use tool is in urgent need for automatically detecting and annotating chromothripsis. Here we present CTLPScanner, a web server for detection of chromothripsis-like pattern (CTLP) in genomic array data. The output interface presents intuitive graphical representations of detected chromosome pulverization region, as well as detailed results in table format. CTLPScanner also provides additional information for associated genes in chromothripsis region to help identify the potential candidates involved in tumorigenesis. To assist in performing meta-data analysis, we integrated over 50 000 pre-processed genomic arrays from The Cancer Genome Atlas and Gene Expression Omnibus into CTLPScanner. The server allows users to explore the presence of chromothripsis signatures from public data resources, without carrying out any local data processing. CTLPScanner is freely available at http://cgma.scu.edu.cn/CTLPScanner/.

Three-dimensional disorganization of the cancer genome occurs coincident with long-range genetic and epigenetic alterations

A three-dimensional chromatin state underpins the structural and functional basis of the genome by bringing regulatory elements and genes into close spatial proximity to ensure proper, cell-type–specific gene expression profiles. Here, we performed Hi-C chromosome conformation capture sequencing to investigate how three-dimensional chromatin organization is disrupted in the context of copy-number variation, long-range epigenetic remodeling, and atypical gene expression programs in prostate cancer. We find that cancer cells retain the ability to segment their genomes into megabase-sized topologically associated domains (TADs); however, these domains are generally smaller due to establishment of additional domain boundaries. Interestingly, a large proportion of the new cancer-specific domain boundaries occur at regions that display copy-number variation. Notably, a common deletion on 17p13.1 in prostate cancer spanning the TP53 tumor suppressor locus results in bifurcation of a single TAD into two distinct smaller TADs. Change in domain structure is also accompanied by novel cancer-specific chromatin interactions within the TADs that are enriched at regulatory elements such as enhancers, promoters, and insulators, and associated with alterations in gene expression. We also show that differential chromatin interactions across regulatory regions occur within long-range epigenetically activated or silenced regions of concordant gene activation or repression in prostate cancer. Finally, we present a novel visualization tool that enables integrated exploration of Hi-C interaction data, the transcriptome, and epigenome. This study provides new insights into the relationship between long-range epigenetic and genomic dysregulation and changes in higher-order chromatin interactions in cancer.

Abundant copy-number loss of CYCLOPS and STOP genes in gastric adenocarcinoma

BACKGROUND

Gastric cancer, a leading cause of cancer death worldwide, has been little studied compared with other cancers that impose similar health burdens. Our goal is to assess genomic copy-number loss and the possible functional consequences and therapeutic implications thereof across a large series of gastric adenocarcinomas.

METHODS

We used high-density single-nucleotide polymorphism microarrays to determine patterns of copy-number loss and allelic imbalance in 74 gastric adenocarcinomas. We investigated whether suppressor of tumorigenesis and/or proliferation (STOP) genes are associated with genomic copy-number loss. We also analyzed the extent to which copy-number loss affects Copy-number alterations Yielding Cancer Liabilities Owing to Partial losS (CYCLOPS) genes-genes that may be attractive targets for therapeutic inhibition when partially deleted.

RESULTS

The proportion of the genome subject to copy-number loss varies considerably from tumor to tumor, with a median of 5.5 %, and a mean of 12 % (range 0-58.5 %). On average, 91 STOP genes were subject to copy-number loss per tumor (median 35, range 0-452), and STOP genes tended to have lower copy-number compared with the rest of the genes. Furthermore, on average, 1.6 CYCLOPS genes per tumor were both subject to copy-number loss and downregulated, and 51.4 % of the tumors had at least one such gene.

CONCLUSIONS

The enrichment of STOP genes in regions of copy-number loss indicates that their deletion may contribute to gastric carcinogenesis. Furthermore, the presence of several deleted and downregulated CYCLOPS genes in some tumors suggests potential therapeutic targets in these tumors.

Genomic landscape of liposarcoma

Liposarcoma (LPS) is the most common type of soft tissue sarcoma accounting for 20% of all adult sarcomas. Due to absence of clinically effective treatment options in inoperable situations and resistance to chemotherapeutics, a critical need exists to identify novel therapeutic targets. We analyzed LPS genomic landscape using SNP arrays, whole exome sequencing and targeted exome sequencing to uncover the genomic information for development of specific anti-cancer targets. SNP array analysis indicated known amplified genes (MDM2, CDK4, HMGA2) and important novel genes (UAP1, MIR557, LAMA4, CPM, IGF2, ERBB3, IGF1R). Carboxypeptidase M (CPM), recurrently amplified gene in well-differentiated/de-differentiated LPS was noted as a putative oncogene involved in the EGFR pathway. Notable deletions were found at chromosome 1p (RUNX3, ARID1A), chromosome 11q (ATM, CHEK1) and chromosome 13q14.2 (MIR15A, MIR16-1). Significantly and recurrently mutated genes (false discovery rate < 0.05) included PLEC (27%), MXRA5 (21%), FAT3 (24%), NF1 (20%), MDC1 (10%), TP53 (7%) and CHEK2 (6%). Further, in vitro and in vivo functional studies provided evidence for the tumor suppressor role for Neurofibromin 1 (NF1) gene in different subtypes of LPS. Pathway analysis of recurrent mutations demonstrated signaling through MAPK, JAK-STAT, Wnt, ErbB, axon guidance, apoptosis, DNA damage repair and cell cycle pathways were involved in liposarcomagenesis. Interestingly, we also found mutational and copy number heterogeneity within a primary LPS tumor signifying the importance of multi-region sequencing for cancer-genome guided therapy. In summary, these findings provide insight into the genomic complexity of LPS and highlight potential druggable pathways for targeted therapeutic approach.

Identifying rhesus macaque gene orthologs using heterospecific human CNV probes

We used the Affymetrix(®) Genome-Wide Human SNP Array 6.0 to identify heterospecific markers and compare copy number and structural genomic variation between humans and rhesus macaques. Over 200,000 human copy number variation (CNV) probes were mapped to a Chinese and an Indian rhesus macaque sample. Observed genomic rearrangements and synteny were in agreement with the results of a previously published genomic comparison between humans and rhesus macaques. Comparisons between each of the two rhesus macaques and humans yielded 206 regions with copy numbers that differed by at least two fold in the Indian rhesus macaque and human, 32 in the Chinese rhesus macaque and human, and 147 in both rhesus macaques. The detailed genomic map and preliminary CNV data are useful for better understanding genetic variation in rhesus macaques, identifying derived changes in human CNVs that may have evolved by selection, and determining the suitability of rhesus macaques as human models for particular biomedical studies.

miR-17-92 fine-tunes MYC expression and function to ensure optimal B cell lymphoma growth

The synergism between c-MYC and miR-17-19b, a truncated version of the miR-17-92 cluster, is well-documented during tumor initiation. However, little is known about miR-17-19b function in established cancers. Here we investigate the role of miR-17-19b in c-MYC-driven lymphomas by integrating SILAC-based quantitative proteomics, transcriptomics and 3′ untranslated region (UTR) analysis upon miR-17-19b overexpression. We identify over one hundred miR-17-19b targets, of which 40% are co-regulated by c-MYC. Downregulation of a new miR-17/20 target, checkpoint kinase 2 (Chek2), increases the recruitment of HuR to c-MYC transcripts, resulting in the inhibition of c-MYC translation and thus interfering with in vivo tumor growth. Hence, in established lymphomas, miR-17-19b fine-tunes c-MYC activity through a tight control of its function and expression, ultimately ensuring cancer cell homeostasis. Our data highlight the plasticity of miRNA function, reflecting changes in the mRNA landscape and 3′ UTR shortening at different stages of tumorigenesis.

The synergism between c-MYC and miR-17-19b plays an important role in lymphoma initiation. In this study, the authors identify a panel of targets co-regulated by miR-17-19b and in MYC-driven lymphoma and unravel the molecular mechanism through which miR-17-19b inhibits MYC translation.

A novel genomic alteration of LSAMP associates with aggressive prostate cancer in African American men

Evaluation of cancer genomes in global context is of great interest in light of changing ethnic distribution of the world population. We focused our study on men of African ancestry because of their disproportionately higher rate of prostate cancer (CaP) incidence and mortality. We present a systematic whole genome analyses, revealing alterations that differentiate African American (AA) and Caucasian American (CA) CaP genomes. We discovered a recurrent deletion on chromosome 3q13.31 centering on the LSAMP locus that was prevalent in tumors from AA men (cumulative analyses of 435 patients: whole genome sequence, 14; FISH evaluations, 101; and SNP array, 320 patients). Notably, carriers of this deletion experienced more rapid disease progression. In contrast, PTEN and ERG common driver alterations in CaP were significantly lower in AA prostate tumors compared to prostate tumors from CA. Moreover, the frequency of inter-chromosomal rearrangements was significantly higher in AA than CA tumors. These findings reveal differentially distributed somatic mutations in CaP across ancestral groups, which have implications for precision medicine strategies.

Patterns of somatic uniparental disomy identify novel tumor suppressor genes in colorectal cancer

Colorectal cancer (CRC) is characterized by specific patterns of copy number alterations (CNAs), which helped with the identification of driver oncogenes and tumor suppressor genes (TSGs). More recently, the usage of single nucleotide polymorphism arrays provided information of copy number neutral loss of heterozygosity, thus suggesting the occurrence of somatic uniparental disomy (UPD) and uniparental polysomy (UPP) events. The aim of this study is to establish an integrative profiling of recurrent UPDs/UPPs and CNAs in sporadic CRC. Our results indicate that regions showing high frequencies of UPD/UPP mostly coincide with regions typically involved in genomic losses. Among them, chromosome arms 3p, 5q, 9q, 10q, 14q, 17p, 17q, 20p, 21q and 22q preferentially showed UPDs/UPPs over genomic losses suggesting that tumor cells must maintain the disomic state of certain genes to favor cellular fitness. A meta-analysis using over 300 samples from The Cancer Genome Atlas confirmed our findings. Several regions affected by recurrent UPDs/UPPs contain well-known TSGs, as well as novel candidates such as ARID1A, DLC1, TCF7L2 and DMBT1. In addition, VCAN, FLT4, SFRP1 and GAS7 were also frequently involved in regions of UPD/UPP and displayed high levels of methylation. Finally, sequencing and fluorescence in situ hybridization analysis of the gene APC underlined that a somatic UPD event might represent the second hit to achieve biallelic inactivation of this TSG in colorectal tumors. In summary, our data define a profile of somatic UPDs/UPPs in sporadic CRC and highlights the importance of these events as a mechanism to achieve the inactivation of TSGs.

Allelic loss of 9p21.3 is a prognostic factor in 1p/19q codeleted anaplastic gliomas

OBJECTIVES

We aimed to study the potential clinical relevance of 9p allelic loss, with or without copy number variation, in 1p/19q codeleted anaplastic oligodendroglial tumors (AOTs).

METHODS

This study enrolled 216 patients with 1p/19q codeleted AOT. The prognostic value of 9p allelic loss was investigated using a French nation-wide prospective registry, POLA (prise en charge des tumeurs oligodendrogliales anaplasiques) and high-density single nucleotide polymorphism arrays. We validated our results using the Repository of Molecular Brain Neoplasia Data (REMBRANDT) dataset.

RESULTS

The minimal common region of allelic loss in chromosome arm 9p was 9p21.3. Allelic loss of 9p21.3, detected in 41.7% of tumors, was associated with shorter progression-free and overall survival rates in univariate (p = 0.008 and p < 0.001, respectively) and multivariate analyses (p = 0.009 and p = 0.009, respectively). This finding was validated in the REMBRANDT dataset in univariate and multivariate analysis (p = 0.01 and p = 0.01, respectively).

CONCLUSION

Our study highlights a novel potential prognostic biomarker in 1p/19q codeleted AOT. Further prospective studies are warranted to investigate our finding.

Genome-Wide Loss of Heterozygosity and DNA Copy Number Aberration in HPV-Negative Oral Squamous Cell Carcinoma and Their Associations with Disease-Specific Survival

Oral squamous cell cancer of the oral cavity and oropharynx (OSCC) is associated with high case-fatality. For reasons that are largely unknown, patients with the same clinical and pathologic staging have heterogeneous response to treatment and different probability of recurrence and survival, with patients with Human Papillomavirus (HPV)-positive oropharyngeal tumors having the most favorable survival. To gain insight into the complexity of OSCC and to identify potential chromosomal changes that may be associated with OSCC mortality, we used Affymtrix 6.0 SNP arrays to examine paired DNA from peripheral blood and tumor cell populations isolated by laser capture microdissection to assess genome-wide loss of heterozygosity (LOH) and DNA copy number aberration (CNA) and their associations with risk factors, tumor characteristics, and oral cancer-specific mortality among 75 patients with HPV-negative OSCC. We found a highly heterogeneous and complex genomic landscape of HPV-negative tumors, and identified regions in 4q, 8p, 9p and 11q that seem to play an important role in oral cancer biology and survival from this disease. If confirmed, these findings could assist in designing personalized treatment or in the creation of models to predict survival in patients with HPV-negative OSCC.

Family-Based Benchmarking of Copy Number Variation Detection Software

The analysis of structural variants, in particular of copy-number variations (CNVs), has proven valuable in unraveling the genetic basis of human diseases. Hence, a large number of algorithms have been developed for the detection of CNVs in SNP array signal intensity data. Using the European and African HapMap trio data, we undertook a comparative evaluation of six commonly used CNV detection software tools, namely Affymetrix Power Tools (APT), QuantiSNP, PennCNV, GLAD, R-gada and VEGA, and assessed their level of pair-wise prediction concordance. The tool-specific CNV prediction accuracy was assessed in silico by way of intra-familial validation. Software tools differed greatly in terms of the number and length of the CNVs predicted as well as the number of markers included in a CNV. All software tools predicted substantially more deletions than duplications. Intra-familial validation revealed consistently low levels of prediction accuracy as measured by the proportion of validated CNVs (34-60%). Moreover, up to 20% of apparent family-based validations were found to be due to chance alone. Software using Hidden Markov models (HMM) showed a trend to predict fewer CNVs than segmentation-based algorithms albeit with greater validity. PennCNV yielded the highest prediction accuracy (60.9%). Finally, the pairwise concordance of CNV prediction was found to vary widely with the software tools involved. We recommend HMM-based software, in particular PennCNV, rather than segmentation-based algorithms when validity is the primary concern of CNV detection. QuantiSNP may be used as an additional tool to detect sets of CNVs not detectable by the other tools. Our study also reemphasizes the need for laboratory-based validation, such as qPCR, of CNVs predicted in silico.

Progesterone receptor modulates ERα action in breast cancer

Progesterone receptor (PR) expression is used as a biomarker of oestrogen receptor-α (ERα) function and breast cancer prognosis. Here we show that PR is not merely an ERα-induced gene target, but is also an ERα-associated protein that modulates its behaviour. In the presence of agonist ligands, PR associates with ERα to direct ERα chromatin binding events within breast cancer cells, resulting in a unique gene expression programme that is associated with good clinical outcome. Progesterone inhibited oestrogen-mediated growth of ERα(+) cell line xenografts and primary ERα(+) breast tumour explants, and had increased anti-proliferative effects when coupled with an ERα antagonist. Copy number loss of PGR, the gene coding for PR, is a common feature in ERα(+) breast cancers, explaining lower PR levels in a subset of cases. Our findings indicate that PR functions as a molecular rheostat to control ERα chromatin binding and transcriptional activity, which has important implications for prognosis and therapeutic interventions.

TAFFYS: An Integrated Tool for Comprehensive Analysis of Genomic Aberrations in Tumor Samples

Background

Tumor single nucleotide polymorphism (SNP) array is a common platform for investigating the cancer genomic aberration and the functionally important altered genes. Original SNP array signals are usually corrupted by noise, and need to be de-convoluted into absolute copy number profile by analytical methods. Unfortunately, in contrast with the popularity of tumor Affymetrix SNP array, the methods that are specifically designed for this platform are still limited. The complicated characteristics of noise in signals is one of the difficulties for dissecting tumor Affymetrix SNP array data, as they inevitably blur the distinction between aberrations and create an obstacle for the copy number aberration (CNA) identification.

Results

We propose a tool named TAFFYS for comprehensive analysis of tumor Affymetrix SNP array data. TAFFYS introduce a wavelet-based de-noising approach and copy number-specific signal variance model for suppressing and modelling the noise in signals. Then a hidden Markov model is employed for copy number inference. Finally, by using the absolute copy number profile, statistical significance of each aberration region is calculated in term of different aberration types, including amplification, deletion and loss of heterozygosity (LOH). The result shows that copy number specific-variance model and wavelet de-noising algorithm fits well with the Affymetrix SNP array signals, leading to more accurate estimation for diluted tumor sample (even with only 30% of cancer cells) than other existed methods. Results of examinations also demonstrate a good compatibility and extensibility for different Affymetrix SNP array platforms. Application on the 35 breast tumor samples shows that TAFFYS can automatically dissect the tumor samples and reveal statistically significant aberration regions where cancer-related genes locate.

Conclusions

TAFFYS provide an efficient and convenient tool for identifying the copy number alteration and allelic imbalance and assessing the recurrent aberrations for the tumor Affymetrix SNP array data.

Pan-cancer analysis of genomic scar signatures associated with homologous recombination deficiency suggests novel indications for existing cancer drugs

BACKGROUND

Ovarian and triple-negative breast cancers with BRCA1 or BRCA2 loss are highly sensitive to treatment with PARP inhibitors and platinum-based cytotoxic agents and show an accumulation of genomic scars in the form of gross DNA copy number aberrations. Cancers without BRCA1 or BRCA2 loss but with accumulation of similar genomic scars also show increased sensitivity to platinum-based chemotherapy. Therefore, reliable biomarkers to identify DNA repair-deficient cancers prior to treatment may be useful for directing patients to platinum chemotherapy and possibly PARP inhibitors. Recently, three SNP array-based signatures of chromosomal instability were published that each quantitate a distinct type of genomic scar considered likely to be caused by improper DNA repair. They measure telomeric allelic imbalance (named NtAI), large scale transition (named LST), and loss of heterozygosity (named HRD-LOH), and it is suggested that these signatures may act as biomarkers for the state of DNA repair deficiency in a given cancer.

RESULTS

We explored the pan-cancer distribution of scores of the three signatures utilizing a panel of 5371 tumors representing 15 cancer types from The Cancer Genome Atlas, and found a good correlation between scores of the three signatures (Spearman's ρ 0.73-0.87). In addition we found that cancer types ordinarily receiving platinum as standard of care have higher median scores of all three signatures. Interestingly, we also found that smaller subpopulations of high-scoring tumors exist in most cancer types, including those for which platinum chemotherapy is not standard therapy.

CONCLUSIONS

Within several cancer types that are not ordinarily treated with platinum chemotherapy, we identified tumors with high levels of the three genomic biomarkers. These tumors represent identifiable subtypes of patients which may be strong candidates for clinical trials with PARP inhibitors or platinum-based chemotherapeutic regimens.

Clonal status of actionable driver events and the timing of mutational processes in cancer evolution

Deciphering whether actionable driver mutations are found in all or a subset of tumor cells will likely be required to improve drug development and precision medicine strategies. We analyzed nine cancer types to determine the subclonal frequencies of driver events, to time mutational processes during cancer evolution, and to identify drivers of subclonal expansions. Although mutations in known driver genes typically occurred early in cancer evolution, we also identified later subclonal "actionable" mutations, including BRAF (V600E), IDH1 (R132H), PIK3CA (E545K), EGFR (L858R), and KRAS (G12D), which may compromise the efficacy of targeted therapy approaches. More than 20% of IDH1 mutations in glioblastomas, and 15% of mutations in genes in the PI3K (phosphatidylinositol 3-kinase)-AKT-mTOR (mammalian target of rapamycin) signaling axis across all tumor types were subclonal. Mutations in the RAS-MEK (mitogen-activated protein kinase kinase) signaling axis were less likely to be subclonal than mutations in genes associated with PI3K-AKT-mTOR signaling. Analysis of late mutations revealed a link between APOBEC-mediated mutagenesis and the acquisition of subclonal driver mutations and uncovered putative cancer genes involved in subclonal expansions, including CTNNA2 and ATXN1. Our results provide a pan-cancer census of driver events within the context of intratumor heterogeneity and reveal patterns of tumor evolution across cancers. The frequent presence of subclonal driver mutations suggests the need to stratify targeted therapy response according to the proportion of tumor cells in which the driver is identified.

Whole-genome and multisector exome sequencing of primary and post-treatment glioblastoma reveals patterns of tumor evolution

Glioblastoma (GBM) is a prototypical heterogeneous brain tumor refractory to conventional therapy. A small residual population of cells escapes surgery and chemoradiation, resulting in a typically fatal tumor recurrence ∼7 mo after diagnosis. Understanding the molecular architecture of this residual population is critical for the development of successful therapies. We used whole-genome sequencing and whole-exome sequencing of multiple sectors from primary and paired recurrent GBM tumors to reconstruct the genomic profile of residual, therapy resistant tumor initiating cells. We found that genetic alteration of the p53 pathway is a primary molecular event predictive of a high number of subclonal mutations in glioblastoma. The genomic road leading to recurrence is highly idiosyncratic but can be broadly classified into linear recurrences that share extensive genetic similarity with the primary tumor and can be directly traced to one of its specific sectors, and divergent recurrences that share few genetic alterations with the primary tumor and originate from cells that branched off early during tumorigenesis. Our study provides mechanistic insights into how genetic alterations in primary tumors impact the ensuing evolution of tumor cells and the emergence of subclonal heterogeneity.

BCL11A is a triple-negative breast cancer gene with critical functions in stem and progenitor cells

Triple-negative breast cancer (TNBC) has poor prognostic outcome compared with other types of breast cancer. The molecular and cellular mechanisms underlying TNBC pathology are not fully understood. Here, we report that the transcription factor BCL11A is overexpressed in TNBC including basal-like breast cancer (BLBC) and that its genomic locus is amplified in up to 38% of BLBC tumours. Exogenous BCL11A overexpression promotes tumour formation, whereas its knockdown in TNBC cell lines suppresses their tumourigenic potential in xenograft models. In the DMBA-induced tumour model, Bcl11a deletion substantially decreases tumour formation, even in p53-null cells and inactivation of Bcl11a in established tumours causes their regression. At the cellular level, Bcl11a deletion causes a reduction in the number of mammary epithelial stem and progenitor cells. Thus, BCL11A has an important role in TNBC and normal mammary epithelial cells. This study highlights the importance of further investigation of BCL11A in TNBC-targeted therapies.

Single cell-derived clonal analysis of human glioblastoma links functional and genomic heterogeneity

Glioblastoma (GBM) is a cancer comprised of morphologically, genetically, and phenotypically diverse cells. However, an understanding of the functional significance of intratumoral heterogeneity is lacking. We devised a method to isolate and functionally profile tumorigenic clones from patient glioblastoma samples. Individual clones demonstrated unique proliferation and differentiation abilities. Importantly, naïve patient tumors included clones that were temozolomide resistant, indicating that resistance to conventional GBM therapy can preexist in untreated tumors at a clonal level. Further, candidate therapies for resistant clones were detected with clone-specific drug screening. Genomic analyses revealed genes and pathways that associate with specific functional behavior of single clones. Our results suggest that functional clonal profiling used to identify tumorigenic and drug-resistant tumor clones will lead to the discovery of new GBM clone-specific treatment strategies.

A tumor DNA complex aberration index is an independent predictor of survival in breast and ovarian cancer

Complex focal chromosomal rearrangements in cancer genomes, also called "firestorms", can be scored from DNA copy number data. The complex arm-wise aberration index (CAAI) is a score that captures DNA copy number alterations that appear as focal complex events in tumors, and has potential prognostic value in breast cancer. This study aimed to validate this DNA-based prognostic index in breast cancer and test for the first time its potential prognostic value in ovarian cancer. Copy number alteration (CNA) data from 1950 breast carcinomas (METABRIC cohort) and 508 high-grade serous ovarian carcinomas (TCGA dataset) were analyzed. Cases were classified as CAAI positive if at least one complex focal event was scored. Complex alterations were frequently localized on chromosome 8p (n = 159), 17q (n = 176) and 11q (n = 251). CAAI events on 11q were most frequent in estrogen receptor positive (ER+) cases and on 17q in estrogen receptor negative (ER-) cases. We found only a modest correlation between CAAI and the overall rate of genomic instability (GII) and number of breakpoints (r = 0.27 and r = 0.42, p < 0.001). Breast cancer specific survival (BCSS), overall survival (OS) and ovarian cancer progression free survival (PFS) were used as clinical end points in Cox proportional hazard model survival analyses. CAAI positive breast cancers (43%) had higher mortality: hazard ratio (HR) of 1.94 (95%CI, 1.62-2.32) for BCSS, and of 1.49 (95%CI, 1.30-1.71) for OS. Representations of the 70-gene and the 21-gene predictors were compared with CAAI in multivariable models and CAAI was independently significant with a Cox adjusted HR of 1.56 (95%CI, 1.23-1.99) for ER+ and 1.55 (95%CI, 1.11-2.18) for ER- disease. None of the expression-based predictors were prognostic in the ER- subset. We found that a model including CAAI and the two expression-based prognostic signatures outperformed a model including the 21-gene and 70-gene signatures but excluding CAAI. Inclusion of CAAI in the clinical prognostication tool PREDICT significantly improved its performance. CAAI positive ovarian cancers (52%) also had worse prognosis: HRs of 1.3 (95%CI, 1.1-1.7) for PFS and 1.3 (95%CI, 1.1-1.6) for OS. This study validates CAAI as an independent predictor of survival in both ER+ and ER- breast cancer and reveals a significant prognostic value for CAAI in high-grade serous ovarian cancer.

Sequenza: allele-specific copy number and mutation profiles from tumor sequencing data

BACKGROUND

Exome or whole-genome deep sequencing of tumor DNA along with paired normal DNA can potentially provide a detailed picture of the somatic mutations that characterize the tumor. However, analysis of such sequence data can be complicated by the presence of normal cells in the tumor specimen, by intratumor heterogeneity, and by the sheer size of the raw data. In particular, determination of copy number variations from exome sequencing data alone has proven difficult; thus, single nucleotide polymorphism (SNP) arrays have often been used for this task. Recently, algorithms to estimate absolute, but not allele-specific, copy number profiles from tumor sequencing data have been described.

MATERIALS AND METHODS

We developed Sequenza, a software package that uses paired tumor-normal DNA sequencing data to estimate tumor cellularity and ploidy, and to calculate allele-specific copy number profiles and mutation profiles. We applied Sequenza, as well as two previously published algorithms, to exome sequence data from 30 tumors from The Cancer Genome Atlas. We assessed the performance of these algorithms by comparing their results with those generated using matched SNP arrays and processed by the allele-specific copy number analysis of tumors (ASCAT) algorithm.

RESULTS

Comparison between Sequenza/exome and SNP/ASCAT revealed strong correlation in cellularity (Pearson's r = 0.90) and ploidy estimates (r = 0.42, or r = 0.94 after manual inspecting alternative solutions). This performance was noticeably superior to previously published algorithms. In addition, in artificial data simulating normal-tumor admixtures, Sequenza detected the correct ploidy in samples with tumor content as low as 30%.

CONCLUSIONS

The agreement between Sequenza and SNP array-based copy number profiles suggests that exome sequencing alone is sufficient not only for identifying small scale mutations but also for estimating cellularity and inferring DNA copy number aberrations.

MPAgenomics: an R package for multi-patient analysis of genomic markers

BACKGROUND

Last generations of Single Nucleotide Polymorphism (SNP) arrays allow to study copy-number variations in addition to genotyping measures.

RESULTS

MPAgenomics, standing for multi-patient analysis (MPA) of genomic markers, is an R-package devoted to: (i) efficient segmentation and (ii) selection of genomic markers from multi-patient copy number and SNP data profiles. It provides wrappers from commonly used packages to streamline their repeated (sometimes difficult) manipulation, offering an easy-to-use pipeline for beginners in R.The segmentation of successive multiple profiles (finding losses and gains) is performed with an automatic choice of parameters involved in the wrapped packages. Considering multiple profiles in the same time, MPAgenomics wraps efficient penalized regression methods to select relevant markers associated with a given outcome.

CONCLUSIONS

MPAgenomics provides an easy tool to analyze data from SNP arrays in R. The R-package MPAgenomics is available on CRAN.

Mayer-Rokitansky-Küster-Hauser syndrome discordance in monozygotic twins: matrix metalloproteinase 14, low-density lipoprotein receptor-related protein 10, extracellular matrix, and neoangiogenesis genes identified as candidate genes in a tissue-specific mosaicism

OBJECTIVE

To find a potential underlying cause for Mayer-Rokitansky-Küster-Hauser syndrome (MRKHS) discordance in monozygotic twins.

DESIGN

Prospective comparative study.

SETTING

University hospital.

PATIENT(S)

Our study genetically analyzed 5 MRKHS-discordant monozygotic twin pairs with the unique opportunity to include saliva and rudimentary uterine tissue.

INTERVENTION(S)

Blood, saliva, or rudimentary uterine tissue from five MRKHS-discordant twins was analyzed and compared between twin pairs as well as within the same individual where applicable. We used copy number variations (CNVs) to identify differences.

MAIN OUTCOME MEASURE(S)

CNVs in blood, rudimentary uterine tissue, and saliva, network analysis, and review of the literature.

RESULT(S)

One duplication found in the affected twin included two genes, matrix metalloproteinase 14 (MMP14) and low-density lipoprotein receptor-related protein 10 (LRP10), which have known functions in the embryonic development of the uterus and endometrium. The duplicated region was detected in rudimentary uterine tissue from the same individual but not in saliva, making a tissue-specific mosaicism a possible explanation for twin discordance. Additional network analysis revealed important connections to differentially expressed genes from previous studies. These genes encode several molecules involved in extracellular matrix (ECM) remodeling and neoangiogenesis.

CONCLUSION(S)

MMP-14, LRP-10, ECM, and neoangiogenesis genes are identified as candidate genes in a tissue-specific mosaicism. The detected clusters provide evidence of deficient vascularization during uterine development and/or disturbed reorganization of ECM components, potentially during müllerian duct elongation signaled by the embryologically relevant phosphatidylinositol 3-kinase/protein kinase B pathway. Therefore, we consider these genes to be new candidates in the manifestation of MRKHS.

arrayMap 2014: an updated cancer genome resource

Somatic copy number aberrations (CNA) represent a mutation type encountered in the majority of cancer genomes. Here, we present the 2014 edition of arrayMap (http://www.arraymap.org), a publicly accessible collection of pre-processed oncogenomic array data sets and CNA profiles, representing a vast range of human malignancies. Since the initial release, we have enhanced this resource both in content and especially with regard to data mining support. The 2014 release of arrayMap contains more than 64 000 genomic array data sets, representing about 250 tumor diagnoses. Data sets included in arrayMap have been assembled from public repositories as well as additional resources, and integrated by applying custom processing pipelines. Online tools have been upgraded for a more flexible array data visualization, including options for processing user provided, non-public data sets. Data integration has been improved by mapping to multiple editions of the human reference genome, with the majority of the data now being available for the UCSC hg18 as well as GRCh37 versions. The large amount of tumor CNA data in arrayMap can be freely downloaded by users to promote data mining projects, and to explore special events such as chromothripsis-like genome patterns.