Concordance between an FDA-approved companion diagnostic and an alternative assay kit for assessing homologous recombination deficiency in ovarian cancer.

OBJECTIVE: Authors evaluated the performance of a commercially available next-generation sequencing assay kit; this was based on genomic content from Illumina's TruSight™ Oncology 500 research assay that identifies BRCA variants and proprietary algorithms licensed from Myriad and, with additional genomic content, measures the homologous recombination deficiency (HRD) genomic instability score (GIS) in tumor tissue (TSO 500 HRD assay). METHODS: Data from the TSO 500 HRD assay were compared with data from the Myriad MyChoice®CDx PLUS assay (Myriad assay). Prevalence rates for overall HRD status and BRCA mutations (a deleterious or suspected deleterious BRCA1 or BRCA2 mutation or both) and assay agreement rates for HRD GIS and BRCA analysis were assessed in ovarian tumor samples. Pearson correlations of the continuous HRD GIS and analytic sensitivity and specificity were evaluated. RESULTS: The prevalence of overall HRD positivity was 51.2% (TSO 500 HRD assay) versus 49.2% (Myriad assay) and the prevalence of BRCA mutations was 27.6% (TSO 500 HRD assay) versus 25.5% (Myriad assay). After post-processing optimization, concordance of the HRD GIS was 0.980 in all samples and 0.976 in the non-BRCA mutation cohort; the area under the receiver operating characteristic curve was 0.995 and 0.992, respectively. CONCLUSIONS: Comparison between the Illumina and Myriad assays showed that overall HRD status, the individual components of BRCA analysis, and HRD GIS detection results were highly concordant (>93%), suggesting the TSO 500 HRD assay will approach the analytical accuracy of the FDA-approved Myriad assay.

The interplay between cancer type, panel size and tumor mutational burden threshold in patient selection for cancer immunotherapy.

The tumor mutational burden (TMB) is increasingly recognized as an emerging biomarker that predicts improved outcomes or response to immune checkpoint inhibitors in cancer. A multitude of technical and biological factors make it difficult to compare TMB values across platforms, histologies, and treatments. Here, we present a mechanistic model that explains the association between panel size, histology, and TMB threshold with panel performance and survival outcome and demonstrate the limitations of existing methods utilized to harmonize TMB across platforms.

The CDX1-microRNA-215 axis regulates colorectal cancer stem cell differentiation.

The transcription factor caudal-type homeobox 1 (CDX1) is a key regulator of differentiation in the normal colon and in colorectal cancer (CRC). CDX1 activates the expression of enterocyte genes, but it is not clear how the concomitant silencing of stem cell genes is achieved. MicroRNAs (miRNAs) are important mediators of gene repression and have been implicated in tumor suppression and carcinogenesis, but the roles of miRNAs in differentiation, particularly in CRC, remain poorly understood. Here, we identified microRNA-215 (miR-215) as a direct transcriptional target of CDX1 by using high-throughput small RNA sequencing to profile miRNA expression in two pairs of CRC cell lines: CDX1-low HCT116 and HCT116 with stable CDX1 overexpression, and CDX1-high LS174T and LS174T with stable CDX1 knockdown. Validation of candidate miRNAs identified by RNA-seq in a larger cell-line panel revealed miR-215 to be most significantly correlated with CDX1 expression. Quantitative ChIP-PCR and promoter luciferase assays confirmed that CDX1 directly activates miR-215 transcription. miR-215 expression is depleted in FACS-enriched cancer stem cells compared with unsorted samples. Overexpression of miR-215 in poorly differentiated cell lines causes a decrease in clonogenicity, whereas miR-215 knockdown increases clonogenicity and impairs differentiation in CDX1-high cell lines. We identified the genome-wide targets of miR-215 and found that miR-215 mediates the repression of cell cycle and stemness genes downstream of CDX1. In particular, the miR-215 target gene BMI1 has been shown to promote stemness and self-renewal and to vary inversely with CDX1. Our work situates miR-215 as a link between CDX1 expression and BMI1 repression that governs differentiation in CRC.

Somatic mutations in murine models of leukemia and lymphoma: Disease specificity and clinical relevance.

Malignant transformation is a multistep process that is dictated by the acquisition of multiple genomic aberrations that provide growth and survival advantage. During the post genomic era, high throughput genomic sequencing has advanced exponentially, leading to identification of countless cancer associated mutations with potential for targeted therapy. Mouse models of cancer serve as excellent tools to examine the functionality of gene mutations and their contribution to the malignant process. However, it remains unclear whether the genetic events that occur during transformation are similar in mice and humans. To address that, we chose several transgenic mouse models of hematopoietic malignancies and identified acquired mutations in these mice by means of targeted re-sequencing of known cancer-associated genes as well as whole exome sequencing. We found that mutations that are typically found in acute myeloid leukemia or T cell acute lymphoblastic leukemia patients are also common in mouse models of the respective disease. Moreover, we found that the most frequent mutations found in a mouse model of lymphoma occur in a set of epigenetic modifier genes, implicating this pathway in the generation of lymphoma. These results demonstrate that genetically engineered mouse models (GEMM) mimic the genetic evolution of human cancer and serve as excellent platforms for target discovery and validation.

EWS-FLI1 employs an E2F switch to drive target gene expression.

Cell cycle progression is orchestrated by E2F factors. We previously reported that in ETS-driven cancers of the bone and prostate, activating E2F3 cooperates with ETS on target promoters. The mechanism of target co-regulation remained unknown. Using RNAi and time-resolved chromatin-immunoprecipitation in Ewing sarcoma we report replacement of E2F3/pRB by constitutively expressed repressive E2F4/p130 complexes on target genes upon EWS-FLI1 modulation. Using mathematical modeling we interrogated four alternative explanatory models for the observed EWS-FLI1/E2F3 cooperation based on longitudinal E2F target and regulating transcription factor expression analysis. Bayesian model selection revealed the formation of a synergistic complex between EWS-FLI1 and E2F3 as the by far most likely mechanism explaining the observed kinetics of E2F target induction. Consequently we propose that aberrant cell cycle activation in Ewing sarcoma is due to the de-repression of E2F targets as a consequence of transcriptional induction and physical recruitment of E2F3 by EWS-FLI1 replacing E2F4 on their target promoters.

Oncogenic ETS fusions deregulate E2F3 target genes in Ewing sarcoma and prostate cancer.

Deregulated E2F transcription factor activity occurs in the vast majority of human tumors and has been solidly implicated in disturbances of cell cycle control, proliferation, and apoptosis. Aberrant E2F regulatory activity is often caused by impairment of control through pRB function, but little is known about the interplay of other oncoproteins with E2F. Here we show that ETS transcription factor fusions resulting from disease driving rearrangements in Ewing sarcoma (ES) and prostate cancer (PC) are one such class of oncoproteins. We performed an integrative study of genome-wide DNA-binding and transcription data in EWSR1/FLI1 expressing ES and TMPRSS2/ERG containing PC cells. Supported by promoter activity and mutation analyses, we demonstrate that a large fraction of E2F3 target genes are synergistically coregulated by these aberrant ETS proteins. We propose that the oncogenic effect of ETS fusion oncoproteins is in part mediated by the disruptive effect of the E2F-ETS interaction on cell cycle control. Additionally, a detailed analysis of the regulatory targets of the characteristic EWSR1/FLI1 fusion in ES identifies two functionally distinct gene sets. While synergistic regulation in concert with E2F in the promoter of target genes has a generally activating effect, EWSR1/FLI1 binding independent of E2F3 is predominantly associated with repressed differentiation genes. Thus, EWSR1/FLI1 appears to promote oncogenesis by simultaneously promoting cell proliferation and perturbing differentiation.

A chromatin structure-based model accurately predicts DNA replication timing in human cells.

The metazoan genome is replicated in precise cell lineage-specific temporal order. However, the mechanism controlling this orchestrated process is poorly understood as no molecular mechanisms have been identified that actively regulate the firing sequence of genome replication. Here, we develop a mechanistic model of genome replication capable of predicting, with accuracy rivaling experimental repeats, observed empirical replication timing program in humans. In our model, replication is initiated in an uncoordinated (time-stochastic) manner at well-defined sites. The model contains, in addition to the choice of the genomic landmark that localizes initiation, only a single adjustable parameter of direct biological relevance: the number of replication forks. We find that DNase-hypersensitive sites are optimal and independent determinants of DNA replication initiation. We demonstrate that the DNA replication timing program in human cells is a robust emergent phenomenon that, by its very nature, does not require a regulatory mechanism determining a proper replication initiation firing sequence.

Genome-wide depletion of replication initiation events in highly transcribed regions.

This report investigates the mechanisms by which mammalian cells coordinate DNA replication with transcription and chromatin assembly. In yeast, DNA replication initiates within nucleosome-free regions, but studies in mammalian cells have not revealed a similar relationship. Here, we have used genome-wide massively parallel sequencing to map replication initiation events, thereby creating a database of all replication initiation sites within nonrepetitive DNA in two human cell lines. Mining this database revealed that genomic regions transcribed at moderate levels were generally associated with high replication initiation frequency. In genomic regions with high rates of transcription, very few replication initiation events were detected. High-resolution mapping of replication initiation sites showed that replication initiation events were absent from transcription start sites but were highly enriched in adjacent, downstream sequences. Methylation of CpG sequences strongly affected the location of replication initiation events, whereas histone modifications had minimal effects. These observations suggest that high levels of transcription interfere with formation of pre-replication protein complexes. Data presented here identify replication initiation sites throughout the genome, providing a foundation for further analyses of DNA-replication dynamics and cell-cycle progression.

Enforced expression of Lin28b leads to impaired T-cell development, release of inflammatory cytokines, and peripheral T-cell lymphoma.

LIN28A and LIN28B, the mammalian homologs of lin-28, are implicated in malignant transformation in part because of their ability to promote degradation of the let-7 family of miRs. In the present study, we show that overexpression of Lin28b in vivo leads to an aggressive peripheral T-cell lymphoma (PTCL) characterized by widespread infiltration of parenchymal organs with malignant CD4(+) cells. Similar to patients with PTCL, Lin28b-transgenic mice show signs of inflammation such as eosinophilia, increased C-reactive protein, release of inflammatory cytokines, and pleural effusion. The PTCLs that develop in Lin28b mice are derived from activated T cells and show decreased let-7 expression, increased Il6 expression, activation of NF-κB, and infiltration of B cells, all resulting in an inflammatory microenvironment. In addition, LIN28B is overexpressed 7.5-fold in PTCL patient samples compared with activated CD4(+) cells. The results of the present study demonstrate for the first time that Lin28b can transform primary cells in vivo, identify a previously unsuspected link between Lin28b and PTCL, and provide a unique animal model for the study of PTCL biology and therapy.

A methyl-deviator epigenotype of estrogen receptor-positive breast carcinoma is associated with malignant biology.

We broadly profiled DNA methylation in breast cancers (n = 351) and benign parenchyma (n = 47) for correspondence with disease phenotype, using FFPE diagnostic surgical pathology specimens. Exploratory analysis revealed a distinctive primary invasive carcinoma subclass featuring extreme global methylation deviation. Subsequently, we tested the correlation between methylation remodeling pervasiveness and malignant biological features. A methyl deviation index (MDI) was calculated for each lesion relative to terminal ductal-lobular unit baseline, and group comparisons revealed that high-grade and short-survival estrogen receptor-positive (ER(+)) cancers manifest a significantly higher MDI than low-grade and long-survival ER(+) cancers. In contrast, ER(-) cancers display a significantly lower MDI, revealing a striking epigenomic distinction between cancer hormone receptor subtypes. Kaplan-Meier survival curves of MDI-based risk classes showed significant divergence between low- and high-risk groups. MDI showed superior prognostic performance to crude methylation levels, and MDI retained prognostic significance (P < 0.01) in Cox multivariate analysis, including clinical stage and pathological grade. Most MDI targets individually are significant markers of ER(+) cancer survival. Lymphoid and mesenchymal indexes were not substantially different between ER(+) and ER(-) groups and do not explain MDI dichotomy. However, the mesenchymal index was associated with ER(+) cancer survival, and a high lymphoid index was associated with medullary carcinoma. Finally, a comparison between metastases and primary tumors suggests methylation patterns are established early and maintained through disease progression for both ER(+) and ER(-) tumors.

Dynamic bookmarking of primary response genes by p300 and RNA polymerase II complexes.

Profiling the dynamic interaction of p300 with proximal promoters of human T cells identified a class of genes that rapidly coassemble p300 and RNA polymerase II (pol II) following mitogen stimulation. Several of these p300 targets are immediate early genes, including FOS, implicating a prominent role for p300 in the control of primary genetic responses. The recruitment of p300 and pol II rapidly transitions to the assembly of several elongation factors, including the positive transcriptional elongation factor (P-TEFb), the bromodomain-containing protein (BRD4), and the elongin-like eleven nineteen lysine-rich leukemia protein (ELL). However, transcription at many of these rapidly induced genes is transient, wherein swift departure of P-TEFb, BRD4, and ELL coincides with termination of transcriptional elongation. Unexpectedly, both p300 and pol II remain accumulated or "bookmarked" at the proximal promoter long after transcription has terminated, demarking a clear mechanistic separation between the recruitment and elongation phases of transcription in vivo. The bookmarked pol II is depleted of both serine-2 and serine-5 phosphorylation of its C-terminal domain and remains proximally positioned at the promoter for hours. Surprisingly, these p300/pol II bookmarked genes can be readily reactivated, and elongation factors can be reassembled by subsequent addition of nonmitogenic agents that, alone, have minimal effects on transcription in the absence of prior preconditioning by mitogen stimulation. These findings suggest that p300 is likely to play an important role in biological processes in which transcriptional bookmarking or preconditioning influences cellular growth and development through the dynamic priming of genes for response to rechallenge by secondary stimuli.

Author Correction: Interaction between the microbiome and TP53 in human lung cancer.

Following publication of the original paper [1], the authors submitted a new Additional file 5 to replace the one containing formatting issues. The updated Additional file 5 is published in this correction.

Whole chromosome alterations predict survival in high-risk neuroblastoma without MYCN amplification.

PURPOSE: Patients with stage IV neuroblastoma over the age of 500 days without MYCN amplification have a survival rate of <30% and there are currently no reliable means of predicting which of these patients will survive or succumb to the disease. The goal of this study is to develop a DNA copy number-based prognostic profile for these patients. EXPERIMENTAL DESIGN: We have used comparative genomic hybridization to identify genome copy number changes that can predict outcome in patients with stage IV neuroblastoma without MYCN amplification. RESULTS: A strong correlation of patient survival with the presence of whole chromosome changes (WCC >or=2) was observed, even in the group of patients older than 500 days at time of diagnosis. This novel prognostic marker showed a significant dependence on the date of diagnosis; patients with WCC >or=2 diagnosed after 1998 had a significantly higher probability of survival compared with those diagnosed earlier. At the same time, no such time dependence was found among the samples with WCC <2, suggesting that medical progress patients in recent years has particularly benefited those patients with a stage IV non-MYCN-amplified disease if WCC >or=2 were present. CONCLUSIONS: In this pilot study, we present a novel prognostic marker for survival of high-risk neuroblastoma patients over the age of 500 days without MYCN amplification and diagnosed after 1998. Further validation study is required to establish this risk stratification for these patients.

An integrated cross-platform prognosis study on neuroblastoma patients.

There have been several reports about the potential for predicting prognosis of neuroblastoma patients using microarray gene expression profiling of the tumors. However these studies have revealed an apparent diversity in the identity of the genes in their predictive signatures. To test the contribution of the platform to this discrepancy we applied the z-scoring method to minimize the impact of platform and combine gene expression profiles of neuroblastoma (NB) tumors from two different platforms, cDNA and Affymetrix. A total of 12442 genes were common to both cDNA and Affymetrix arrays in our data set. Two-way ANOVA analysis was applied to the combined data set for assessing the relative effect of prognosis and platform on gene expression. We found that 26.6% (3307) of the genes had significant impact on survival. There was no significant impact of microarray platform on expression after application of z-scoring standardization procedure. Artificial neural network (ANN) analysis of the combined data set in a leave-one-out prediction strategy correctly predicted the outcome for 90% of the samples. Hierarchical clustering analysis using the top-ranked 160 genes showed the great separation of two clusters, and the majority of matched samples from the different platforms were clustered next to each other. The ANN classifier trained with our combined cross-platform data for these 160 genes could predict the prognosis of 102 independent test samples with 71% accuracy. Furthermore it correctly predicted the outcome for 85/102 (83%) NB patients through the leave-one-out cross-validation approach. Our study showed that gene expression studies performed in different platforms could be integrated for prognosis analysis after removing variation resulting from different platforms.

Credentialing preclinical pediatric xenograft models using gene expression and tissue microarray analysis.

Human tumor xenografts have been used extensively for rapid screening of the efficacy of anticancer drugs for the past 35 years. The selection of appropriate xenograft models for drug testing has been largely empirical and has not incorporated a similarity to the tumor type of origin at the molecular level. This study is the first comprehensive analysis of the transcriptome of a large set of pediatric xenografts, which are currently used for preclinical drug testing. Suitable models representing the tumor type of origin were identified. It was found that the characteristic expression patterns of the primary tumors were maintained in the corresponding xenografts for the majority of samples. Because a prerequisite for developing rationally designed drugs is that the target is expressed at the protein level, we developed tissue arrays from these xenografts and corroborated that high mRNA levels yielded high protein levels for two tested genes. The web database and availability of tissue arrays will allow for the rapid confirmation of the expression of potential targets at both the mRNA and the protein level for molecularly targeted agents. The database will facilitate the identification of tumor markers predictive of response to tested agents as well as the discovery of new molecular targets.

A Genome-Wide Screen Reveals a Role for the HIR Histone Chaperone Complex in Preventing Mislocalization of Budding Yeast CENP-A.

Centromeric localization of the evolutionarily conserved centromere-specific histone H3 variant CENP-A (Cse4 in yeast) is essential for faithful chromosome segregation. Overexpression and mislocalization of CENP-A lead to chromosome segregation defects in yeast, flies, and human cells. Overexpression of CENP-A has been observed in human cancers; however, the molecular mechanisms preventing CENP-A mislocalization are not fully understood. Here, we used a genome-wide synthetic genetic array (SGA) to identify gene deletions that exhibit synthetic dosage lethality (SDL) when Cse4 is overexpressed. Deletion for genes encoding the replication-independent histone chaperone HIR complex (HIR1, HIR2, HIR3, HPC2) and a Cse4-specific E3 ubiquitin ligase, PSH1, showed highest SDL. We defined a role for Hir2 in proteolysis of Cse4 that prevents mislocalization of Cse4 to noncentromeric regions for genome stability. Hir2 interacts with Cse4 in vivo, and hir2∆ strains exhibit defects in Cse4 proteolysis and stabilization of chromatin-bound Cse4 Mislocalization of Cse4 to noncentromeric regions with a preferential enrichment at promoter regions was observed in hir2∆ strains. We determined that Hir2 facilitates the interaction of Cse4 with Psh1, and that defects in Psh1-mediated proteolysis contribute to increased Cse4 stability and mislocalization of Cse4 in the hir2∆ strain. In summary, our genome-wide screen provides insights into pathways that regulate proteolysis of Cse4 and defines a novel role for the HIR complex in preventing mislocalization of Cse4 by facilitating proteolysis of Cse4, thereby promoting genome stability.

Interaction between the microbiome and TP53 in human lung cancer.

BACKGROUND: Lung cancer is the leading cancer diagnosis worldwide and the number one cause of cancer deaths. Exposure to cigarette smoke, the primary risk factor in lung cancer, reduces epithelial barrier integrity and increases susceptibility to infections. Herein, we hypothesize that somatic mutations together with cigarette smoke generate a dysbiotic microbiota that is associated with lung carcinogenesis. Using lung tissue from 33 controls and 143 cancer cases, we conduct 16S ribosomal RNA (rRNA) bacterial gene sequencing, with RNA-sequencing data from lung cancer cases in The Cancer Genome Atlas serving as the validation cohort. RESULTS: Overall, we demonstrate a lower alpha diversity in normal lung as compared to non-tumor adjacent or tumor tissue. In squamous cell carcinoma specifically, a separate group of taxa are identified, in which Acidovorax is enriched in smokers. Acidovorax temporans is identified within tumor sections by fluorescent in situ hybridization and confirmed by two separate 16S rRNA strategies. Further, these taxa, including Acidovorax, exhibit higher abundance among the subset of squamous cell carcinoma cases with TP53 mutations, an association not seen in adenocarcinomas. CONCLUSIONS: The results of this comprehensive study show both microbiome-gene and microbiome-exposure interactions in squamous cell carcinoma lung cancer tissue. Specifically, tumors harboring TP53 mutations, which can impair epithelial function, have a unique bacterial consortium that is higher in relative abundance in smoking-associated tumors of this type. Given the significant need for clinical diagnostic tools in lung cancer, this study may provide novel biomarkers for early detection.

Functional profiling of uncommon VCAM1 promoter polymorphisms prevalent in African American populations.

Multiple variants of the vascular adhesion molecule-1 (VCAM1) promoter show increased nucleotide heterozygosity in the African American population. Using a novel transfection-based transcriptional pathway profiling method, we show that select uncommon variants are functionally hyperactive. Eight candidate VCAM1 promoter haplotypes comprising 13 previously identified SNPs were assessed for response to known mitogens. Activity was correlated with bioinformatic analysis of hyper- and hyporesponsive variants to identify the gain or loss of haplotype-specific transcription factor binding site (TFBS). Using this approach, a low frequency regulatory allele (c.-540A>G; dbSNP rs3783605:A>G), found in a hyperactive VCAM1 promoter haplotype, was shown to create a candidate binding site for ETS2 that was confirmed in vivo by chromatin immunoprecipitation. This report provides the first functional evaluation of VCAM1 promoter polymorphisms and establishes a hypothetical foundation for investigation of their role in the pathogenesis of VCAM1-associated diseases that disproportionately afflict African Americans, including thromboembolic diseases, asthma, and multiple myeloma.

Early B-cell factor, E2A, and Pax-5 cooperate to activate the early B cell-specific mb-1 promoter.

Previous studies have suggested that the early-B-cell-specific mb-1(Igalpha) promoter is regulated by EBF and Pax-5. Here, we used in vivo footprinting assays to detect occupation of binding sites in endogenous mb-1 promoters at various stages of B-cell differentiation. In addition to EBF and Pax-5 binding sites, we detected occupancy of a consensus binding site for E2A proteins (E box) in pre-B cells. EBF and E box sites are crucial for promoter function in transfected pre-B cells, and EBF and E2A proteins synergistically activated the promoter in transfected HeLa cells. Other data suggest that EBF and E box sites are less important for promoter function at later stages of differentiation, whereas binding sites for Pax-5 (and its Ets ternary complex partners) are required for promoter function in all mb-1-expressing cells. Using DNA microarrays, we found that expression of endogenous mb-1 transcripts correlates most closely with EBF expression and negatively with Id1, an inhibitor of E2A protein function, further linking regulation of the mb-1 gene with EBF and E2A. Together, our studies demonstrate the complexity of factors regulating tissue-specific transcription and support the concept that EBF, E2A, and Pax-5 cooperate to activate target genes in early B-cell development.

Analysis of comparative genomic hybridization data on cDNA microarrays.

We present a detailed method to analyze DNA copy number data generated on cDNA microarrays. A web interface is made available for those steps in the workflow that are not typically used in gene expression analysis so that these steps can be carried out online. The end result of the analysis is a list of p-values for the presence of genomic gains or losses for each sample individually or an average p-value, which we show is useful to identify recurrent genomic imbalances.