Tangent normalization for somatic copy-number inference in cancer genome analysis.

MOTIVATION: Somatic copy-number alterations (SCNAs) play an important role in cancer development. Systematic noise in sequencing and array data present a significant challenge to the inference of SCNAs for cancer genome analyses. As part of The Cancer Genome Atlas, the Broad Institute Genome Characterization Center developed the Tangent normalization method to generate copy-number profiles using data from single-nucleotide polymorphism (SNP) arrays and whole-exome sequencing (WES) technologies for over 10 000 pairs of tumors and matched normal samples. Here, we describe the Tangent method, which uses a unique linear combination of normal samples as a reference for each tumor sample, to subtract systematic errors that vary across samples. We also describe a modification of Tangent, called Pseudo-Tangent, which enables denoising through comparisons between tumor profiles when few normal samples are available. RESULTS: Tangent normalization substantially increases signal-to-noise ratios (SNRs) compared to conventional normalization methods in both SNP array and WES analyses. Tangent and Pseudo-Tangent normalizations improve the SNR by reducing noise with minimal effect on signal and exceed the contribution of other steps in the analysis such as choice of segmentation algorithm. Tangent and Pseudo-Tangent are broadly applicable and enable more accurate inference of SCNAs from DNA sequencing and array data. AVAILABILITY AND IMPLEMENTATION: Tangent is available at https://github.com/broadinstitute/tangent and as a Docker image (https://hub.docker.com/r/broadinstitute/tangent). Tangent is also the normalization method for the copy-number pipeline in Genome Analysis Toolkit 4 (GATK4). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Higher-Order Inter-chromosomal Hubs Shape 3D Genome Organization in the Nucleus.

Eukaryotic genomes are packaged into a 3-dimensional structure in the nucleus. Current methods for studying genome-wide structure are based on proximity ligation. However, this approach can fail to detect known structures, such as interactions with nuclear bodies, because these DNA regions can be too far apart to directly ligate. Accordingly, our overall understanding of genome organization remains incomplete. Here, we develop split-pool recognition of interactions by tag extension (SPRITE), a method that enables genome-wide detection of higher-order interactions within the nucleus. Using SPRITE, we recapitulate known structures identified by proximity ligation and identify additional interactions occurring across larger distances, including two hubs of inter-chromosomal interactions that are arranged around the nucleolus and nuclear speckles. We show that a substantial fraction of the genome exhibits preferential organization relative to these nuclear bodies. Our results generate a global model whereby nuclear bodies act as inter-chromosomal hubs that shape the overall packaging of DNA in the nucleus.

Comparative Analysis of Immune Cells Reveals a Conserved Regulatory Lexicon.

Most well-characterized enhancers are deeply conserved. In contrast, genome-wide comparative studies of steady-state systems showed that only a small fraction of active enhancers are conserved. To better understand conservation of enhancer activity, we used a comparative genomics approach that integrates temporal expression and epigenetic profiles in an innate immune system. We found that gene expression programs diverge among mildly induced genes, while being highly conserved for strongly induced genes. The fraction of conserved enhancers varies greatly across gene expression programs, with induced genes and early-response genes, in particular, being regulated by a higher fraction of conserved enhancers. Clustering of conserved accessible DNA sequences within enhancers resulted in over 60 sequence motifs including motifs for known factors, as well as many with unknown function. We further show that the number of instances of these motifs is a strong predictor of the responsiveness of a gene to pathogen detection.

High Epstein-Barr Virus Load and Genomic Diversity Are Associated with Generation of gp350-Specific Neutralizing Antibodies following Acute Infectious Mononucleosis.

The Epstein-Barr virus (EBV) gp350 glycoprotein interacts with the cellular receptor to mediate viral entry and is thought to be the major target for neutralizing antibodies. To better understand the role of EBV-specific antibodies in the control of viral replication and the evolution of sequence diversity, we measured EBV gp350-specific antibody responses and sequenced the gp350 gene in samples obtained from individuals experiencing primary EBV infection (acute infectious mononucleosis [AIM]) and again 6 months later (during convalescence [CONV]). EBV gp350-specific IgG was detected in the sera of 17 (71%) of 24 individuals at the time of AIM and all 24 (100%) individuals during CONV; binding antibody titers increased from AIM through CONV, reaching levels equivalent to those in age-matched, chronically infected individuals. Antibody-dependent cell-mediated phagocytosis (ADCP) was rarely detected during AIM (4 of 24 individuals; 17%) but was commonly detected during CONV (19 of 24 individuals; 79%). The majority (83%) of samples taken during AIM neutralized infection of primary B cells; all samples obtained at 6 months postdiagnosis neutralized EBV infection of cultured and primary target cells. Deep sequencing revealed interpatient gp350 sequence variation but conservation of the CR2-binding site. The levels of gp350-specific neutralizing activity directly correlated with higher peripheral blood EBV DNA levels during AIM and a greater evolution of diversity in gp350 nucleotide sequences from AIM to CONV. In summary, we conclude that the viral load and EBV gp350 diversity during early infection are associated with the development of neutralizing antibody responses following AIM. IMPORTANCE: Antibodies against viral surface proteins can blunt the spread of viral infection by coating viral particles, mediating uptake by immune cells, or blocking interaction with host cell receptors, making them a desirable component of a sterilizing vaccine. The EBV surface protein gp350 is a major target for antibodies. We report the detection of EBV gp350-specific antibodies capable of neutralizing EBV infection in vitro The majority of gp350-directed vaccines focus on glycoproteins from lab-adapted strains, which may poorly reflect primary viral envelope diversity. We report some of the first primary gp350 sequences, noting that the gp350 host receptor binding site is remarkably stable across patients and time. However, changes in overall gene diversity were detectable during infection. Patients with higher peripheral blood viral loads in primary infection and greater changes in viral diversity generated more efficient antibodies. Our findings provide insight into the generation of functional antibodies, necessary for vaccine development.

Expression profiles of 151 pediatric low-grade gliomas reveal molecular differences associated with location and histological subtype.

BACKGROUND: Pediatric low-grade gliomas (PLGGs), the most frequent pediatric brain tumor, comprise a heterogeneous group of diseases. Recent genomic analyses suggest that these tumors are mostly driven by mitogene-activated protein kinase (MAPK) pathway alterations. However, little is known about the molecular characteristics inherent to their clinical and histological heterogeneity. METHODS: We performed gene expression profiling on 151 paraffin-embedded PLGGs from different locations, ages, and histologies. Using unsupervised and supervised analyses, we compared molecular features with age, location, histology, and BRAF genomic status. We compared molecular differences with normal pediatric brain expression profiles to observe whether those patterns were mirrored in normal brain. RESULTS: Unsupervised clustering distinguished 3 molecular groups that correlated with location in the brain and histological subtype. "Not otherwise specified" (NOS) tumors did not constitute a unified class. Supratentorial pilocytic astrocytomas (PAs) were significantly enriched with genes involved in pathways related to inflammatory activity compared with infratentorial tumors. Differences based on tumor location were not mirrored in location-dependent differences in expression within normal brain tissue. We identified significant differences between supratentorial PAs and diffuse astrocytomas as well as between supratentorial PAs and dysembryoplastic neuroepithelial tumors but not between supratentorial PAs and gangliogliomas. Similar expression patterns were observed between childhood and adolescent PAs. We identified differences between BRAF-duplicated and V600E-mutated tumors but not between primary and recurrent PLGGs. CONCLUSION: Expression profiling of PLGGs reveals significant differences associated with tumor location, histology, and BRAF genomic status. Supratentorial PAs, in particular, are enriched in inflammatory pathways that appear to be tumor-related.

Functional annotation of native enhancers with a Cas9-histone demethylase fusion.

Understanding of mammalian enhancers is limited by the lack of a technology to rapidly and thoroughly test the cell type-specific function. Here, we use a nuclease-deficient Cas9 (dCas9)-histone demethylase fusion to functionally characterize previously described and new enhancer elements for their roles in the embryonic stem cell state. Further, we distinguish the mechanism of action of dCas9-LSD1 at enhancers from previous dCas9-effectors.

HIV type 1 (HIV-1) proviral reservoirs decay continuously under sustained virologic control in HIV-1-infected children who received early treatment.

BACKGROUND: Early initiation of combination antiretroviral therapy (cART) to human immunodeficiency virus type 1 (HIV-1)-infected infants controls HIV-1 replication and reduces mortality. METHODS: Plasma viremia (lower limit of detection, <2 copies/mL), T-cell activation, HIV-1-specific immune responses, and the persistence of cells carrying replication-competent virus were quantified during long-term effective combination antiretroviral therapy (cART) in 4 perinatally HIV-1-infected youth who received treatment early (the ET group) and 4 who received treatment late (the LT group). Decay in peripheral blood mononuclear cell (PBMC) proviral DNA levels was also measured over time in the ET youth. RESULTS: Plasma viremia was not detected in any ET youth but was detected in all LT youth (median, 8 copies/mL; P = .03). PBMC proviral load was significantly lower in ET youth (median, 7 copies per million PBMCs) than in LT youth (median, 181 copies; P = .03). Replication-competent virus was recovered from all LT youth but only 1 ET youth. Decay in proviral DNA was noted in all 4 ET youth in association with limited T-cell activation and with absent to minimal HIV-1-specific immune responses. CONCLUSIONS: Initiation of early effective cART during infancy significantly limits circulating levels of proviral and replication-competent HIV-1 and promotes continuous decay of viral reservoirs. Continued cART with reduction in HIV-1 reservoirs over time may facilitate HIV-1 eradication strategies.

Epstein-Barr virus latent membrane protein 1 genetic variability in peripheral blood B cells and oropharyngeal fluids.

UNLABELLED: We report the diversity of latent membrane protein 1 (LMP1) gene founder sequences and the level of Epstein-Barr virus (EBV) genome variability over time and across anatomic compartments by using virus genomes amplified directly from oropharyngeal wash specimens and peripheral blood B cells during acute infection and convalescence. The intrahost nucleotide variability of the founder virus was 0.02% across the region sequences, and diversity increased significantly over time in the oropharyngeal compartment (P = 0.004). The LMP1 region showing the greatest level of variability in both compartments, and over time, was concentrated within the functional carboxyl-terminal activating regions 2 and 3 (CTAR2 and CTAR3). Interestingly, a deletion in a proline-rich repeat region (amino acids 274 to 289) of EBV commonly reported in EBV sequenced from cancer specimens was not observed in acute infectious mononucleosis (AIM) patients. Taken together, these data highlight the diversity in circulating EBV genomes and its potential importance in disease pathogenesis and vaccine design. IMPORTANCE: This study is among the first to leverage an improved high-throughput deep-sequencing methodology to investigate directly from patient samples the degree of diversity in Epstein-Barr virus (EBV) populations and the extent to which viral genome diversity develops over time in the infected host. Significant variability of circulating EBV latent membrane protein 1 (LMP1) gene sequences was observed between cellular and oral wash samples, and this variability increased over time in oral wash samples. The significance of EBV genetic diversity in transmission and disease pathogenesis are discussed.

SQSTM1 is a pathogenic target of 5q copy number gains in kidney cancer.

Clear cell renal cell carcinoma (ccRCC) is the most common form of kidney cancer and is often linked to loss of chromosome 3p, which harbors the VHL tumor suppressor gene, loss of chromosome 14q, which includes HIF1A, and gain of chromosome 5q. The relevant target(s) on chromosome 5q is not known. Here, we show that 5q amplification leads to overexpression of the SQSTM1 oncogene in ccRCC lines and tumors. Overexpression of SQSTM1 in ccRCC lines promoted resistance to redox stress and increased soft agar growth, while downregulation of SQSTM1 decreased resistance to redox stress, impaired cellular fitness, and decreased tumor formation. Therefore, the selection pressure to amplify 5q in ccRCC is driven, at least partly, by SQSTM1.

Coordinate activation of Shh and PI3K signaling in PTEN-deficient glioblastoma: new therapeutic opportunities.

In glioblastoma, phosphatidylinositol 3-kinase (PI3K) signaling is frequently activated by loss of the tumor suppressor phosphatase and tensin homolog (PTEN). However, it is not known whether inhibiting PI3K represents a selective and effective approach for treatment. We interrogated large databases and found that sonic hedgehog (SHH) signaling is activated in PTEN-deficient glioblastoma. We demonstrate that the SHH and PI3K pathways synergize to promote tumor growth and viability in human PTEN-deficient glioblastomas. A combination of PI3K and SHH signaling inhibitors not only suppressed the activation of both pathways but also abrogated S6 kinase (S6K) signaling. Accordingly, targeting both pathways simultaneously resulted in mitotic catastrophe and tumor apoptosis and markedly reduced the growth of PTEN-deficient glioblastomas in vitro and in vivo. The drugs tested here appear to be safe in humans; therefore, this combination may provide a new targeted treatment for glioblastoma.

Pan-cancer patterns of somatic copy number alteration.

Determining how somatic copy number alterations (SCNAs) promote cancer is an important goal. We characterized SCNA patterns in 4,934 cancers from The Cancer Genome Atlas Pan-Cancer data set. Whole-genome doubling, observed in 37% of cancers, was associated with higher rates of every other type of SCNA, TP53 mutations, CCNE1 amplifications and alterations of the PPP2R complex. SCNAs that were internal to chromosomes tended to be shorter than telomere-bounded SCNAs, suggesting different mechanisms underlying their generation. Significantly recurrent focal SCNAs were observed in 140 regions, including 102 without known oncogene or tumor suppressor gene targets and 50 with significantly mutated genes. Amplified regions without known oncogenes were enriched for genes involved in epigenetic regulation. When levels of genomic disruption were accounted for, 7% of region pairs were anticorrelated, and these regions tended to encompass genes whose proteins physically interact, suggesting related functions. These results provide insights into mechanisms of generation and functional consequences of cancer-related SCNAs.

Genomic analysis of diffuse pediatric low-grade gliomas identifies recurrent oncogenic truncating rearrangements in the transcription factor MYBL1.

Pediatric low-grade gliomas (PLGGs) are among the most common solid tumors in children but, apart from BRAF kinase mutations or duplications in specific subclasses, few genetic driver events are known. Diffuse PLGGs comprise a set of uncommon subtypes that exhibit invasive growth and are therefore especially challenging clinically. We performed high-resolution copy-number analysis on 44 formalin-fixed, paraffin-embedded diffuse PLGGs to identify recurrent alterations. Diffuse PLGGs exhibited fewer such alterations than adult low-grade gliomas, but we identified several significantly recurrent events. The most significant event, 8q13.1 gain, was observed in 28% of diffuse astrocytoma grade IIs and resulted in partial duplication of the transcription factor MYBL1 with truncation of its C-terminal negative-regulatory domain. A similar recurrent deletion-truncation breakpoint was identified in two angiocentric gliomas in the related gene v-myb avian myeloblastosis viral oncogene homolog (MYB) on 6q23.3. Whole-genome sequencing of a MYBL1-rearranged diffuse astrocytoma grade II demonstrated MYBL1 tandem duplication and few other events. Truncated MYBL1 transcripts identified in this tumor induced anchorage-independent growth in 3T3 cells and tumor formation in nude mice. Truncated transcripts were also expressed in two additional tumors with MYBL1 partial duplication. Our results define clinically relevant molecular subclasses of diffuse PLGGs and highlight a potential role for the MYB family in the biology of low-grade gliomas.

Integrative molecular analysis of intrahepatic cholangiocarcinoma reveals 2 classes that have different outcomes.

BACKGROUND & AIMS: Cholangiocarcinoma, the second most common liver cancer, can be classified as intrahepatic cholangiocarcinoma (ICC) or extrahepatic cholangiocarcinoma. We performed an integrative genomic analysis of ICC samples from a large series of patients. METHODS: We performed a gene expression profile, high-density single-nucleotide polymorphism array, and mutation analyses using formalin-fixed ICC samples from 149 patients. Associations with clinicopathologic traits and patient outcomes were examined for 119 cases. Class discovery was based on a non-negative matrix factorization algorithm and significant copy number variations were identified by Genomic Identification of Significant Targets in Cancer (GISTIC) analysis. Gene set enrichment analysis was used to identify signaling pathways activated in specific molecular classes of tumors, and to analyze their genomic overlap with hepatocellular carcinoma (HCC). RESULTS: We identified 2 main biological classes of ICC. The inflammation class (38% of ICCs) is characterized by activation of inflammatory signaling pathways, overexpression of cytokines, and STAT3 activation. The proliferation class (62%) is characterized by activation of oncogenic signaling pathways (including RAS, mitogen-activated protein kinase, and MET), DNA amplifications at 11q13.2, deletions at 14q22.1, mutations in KRAS and BRAF, and gene expression signatures previously associated with poor outcomes for patients with HCC. Copy number variation-based clustering was able to refine these molecular groups further. We identified high-level amplifications in 5 regions, including 1p13 (9%) and 11q13.2 (4%), and several focal deletions, such as 9p21.3 (18%) and 14q22.1 (12% in coding regions for the SAV1 tumor suppressor). In a complementary approach, we identified a gene expression signature that was associated with reduced survival times of patients with ICC; this signature was enriched in the proliferation class (P < .001). CONCLUSIONS: We used an integrative genomic analysis to identify 2 classes of ICC. The proliferation class has specific copy number alterations, activation of oncogenic pathways, and is associated with worse outcome. Different classes of ICC, based on molecular features, therefore might require different treatment approaches.