Chromothripsis followed by circular recombination drives oncogene amplification in human cancer

The mechanisms behind the evolution of complex genomic amplifications in cancer have remained largely unclear. Using whole-genome sequencing data of the pediatric tumor neuroblastoma, we here identified a type of amplification, termed 'seismic amplification', that is characterized by multiple rearrangements and discontinuous copy number levels. Overall, seismic amplifications occurred in 9.9% (274 of 2,756) of cases across 38 cancer types, and were associated with massively increased copy numbers and elevated oncogene expression. Reconstruction of the development of seismic amplification showed a stepwise evolution, starting with a chromothripsis event, followed by formation of circular extrachromosomal DNA that subsequently underwent repetitive rounds of circular recombination. The resulting amplicons persisted as extrachromosomal DNA circles or had reintegrated into the genome in overt tumors. Together, our data indicate that the sequential occurrence of chromothripsis and circular recombination drives oncogene amplification and overexpression in a substantial fraction of human malignancies.

Quantification of Translocation-Specific ctDNA Provides an Integrating Parameter for Early Assessment of Treatment Response and Risk Stratification in Ewing Sarcoma

PURPOSE

We evaluated the predictive and prognostic value of circulating tumor DNA (ctDNA) in patients with Ewing sarcoma (EWS) treated in the EWING2008 trial.

EXPERIMENTAL DESIGN

Plasma samples from 102 patients with EWS enrolled in the EWING2008 trial were obtained before and during induction chemotherapy. Genomic EWSR1 fusion sequence spanning primers and probes were used for highly specific and sensitive quantification of the levels of ctDNA by digital droplet PCR. ctDNA levels were correlated to established clinical risk factors and outcome parameters.

RESULTS

Pretreatment ctDNA copy numbers were correlated with event-free and overall survival. The reduction in ctDNA levels below the detection limit was observed in most cases after only two blocks of vincristine, ifosfamide, doxorubicin, and etoposide (VIDE) induction chemotherapy. The persistence of ctDNA after two VIDE blocks was a strong predictor of poor outcomes. ctDNA levels correlated well with most established clinical risk factors; an inverse correlation was found only for the histologic response to induction therapy. ctDNA levels did not provide simple representations of tumor volume, but integrated information from various tumor characteristics represented an independent EWS tumor marker with predictive and prognostic value.

CONCLUSIONS

ctDNA copy number in the plasma of patients with EWS is a quantifiable parameter for early risk stratification and can be used as a dynamic noninvasive biomarker for early prediction of treatment response and outcome of patients.

Multimodal analysis of cell-free DNA whole-genome sequencing for pediatric cancers with low mutational burden

Sequencing of cell-free DNA in the blood of cancer patients (liquid biopsy) provides attractive opportunities for early diagnosis, assessment of treatment response, and minimally invasive disease monitoring. To unlock liquid biopsy analysis for pediatric tumors with few genetic aberrations, we introduce an integrated genetic/epigenetic analysis method and demonstrate its utility on 241 deep whole-genome sequencing profiles of 95 patients with Ewing sarcoma and 31 patients with other pediatric sarcomas. Our method achieves sensitive detection and classification of circulating tumor DNA in peripheral blood independent of any genetic alterations. Moreover, we benchmark different metrics for cell-free DNA fragmentation analysis, and we introduce the LIQUORICE algorithm for detecting circulating tumor DNA based on cancer-specific chromatin signatures. Finally, we combine several fragmentation-based metrics into an integrated machine learning classifier for liquid biopsy analysis that exploits widespread epigenetic deregulation and is tailored to cancers with low mutation rates. Clinical associations highlight the potential value of cfDNA fragmentation patterns as prognostic biomarkers in Ewing sarcoma. In summary, our study provides a comprehensive analysis of circulating tumor DNA beyond recurrent genetic aberrations, and it renders the benefits of liquid biopsy more readily accessible for childhood cancers.

Genome-wide Screens Implicate Loss of Cullin Ring Ligase 3 in Persistent Proliferation and Genome Instability in TP53-Deficient Cells

TP53 deficiency is the most common alteration in cancer; however, this alone is typically insufficient to drive tumorigenesis. To identify genes promoting tumorigenesis in combination with TP53 deficiency, we perform genome-wide CRISPR-Cas9 knockout screens coupled with proliferation and transformation assays in isogenic cell lines. Loss of several known tumor suppressors enhances cellular proliferation and transformation. Loss of neddylation pathway genes promotes uncontrolled proliferation exclusively in TP53-deficient cells. Combined loss of CUL3 and TP53 activates an oncogenic transcriptional program governed by the nuclear factor κB (NF-κB), AP-1, and transforming growth factor β (TGF-β) pathways. This program maintains persistent cellular proliferation, induces partial epithelial to mesenchymal transition, and increases DNA damage, genomic instability, and chromosomal rearrangements. Our findings reveal CUL3 loss as a key event stimulating persistent proliferation in TP53-deficient cells. These findings may be clinically relevant, since TP53-CUL3-deficient cells are highly sensitive to ataxia telangiectasia mutated (ATM) inhibition, exposing a vulnerability that could be exploited for cancer treatment.

•

Mixed-effect models with MEMcrispR applied to CRISPR screen analyses

•

Knockout of neddylation genes increases persistent proliferation in TP53^−/− cells

•

TP53^−/−,CUL3^−/− cells exhibit persistent proliferation and partial EMT phenotype

•

TP53^−/−,CUL3^−/− cells show increased DNA damage and display sensitivity to ATM inhibition

Molecular Evolution of Early-Onset Prostate Cancer Identifies Molecular Risk Markers and Clinical Trajectories

Identifying the earliest somatic changes in prostate cancer can give important insights into tumor evolution and aids in stratifying high- from low-risk disease. We integrated whole genome, transcriptome and methylome analysis of early-onset prostate cancers (diagnosis ≤55 years). Characterization across 292 prostate cancer genomes revealed age-related genomic alterations and a clock-like enzymatic-driven mutational process contributing to the earliest mutations in prostate cancer patients. Our integrative analysis identified four molecular subgroups, including a particularly aggressive subgroup with recurrent duplications associated with increased expression of ESRP1, which we validate in 12,000 tissue microarray tumors. Finally, we combined the patterns of molecular co-occurrence and risk-based subgroup information to deconvolve the molecular and clinical trajectories of prostate cancer from single patient samples.

Identification of ZCCHC8 as fusion partner of ROS1 in a case of congenital glioblastoma multiforme with a t(6;12)(q21;q24.3)

Congenital gliobastoma multiforme (GBM) is rare and little is known about the molecular defects underlying the initiation and progression of this tumor type. We present a case of congenital GBM analyzed by conventional cytogenetics, fluorescence in situ hybridization, array comparative genomic hybridization and next generation sequencing. On cytogenetic analysis we detected a reciprocal translocation t(6;12)(q21;q24.3). By sequencing, the translocation was shown to form a fusion between the 5' region of ZCCHC8 and the 3' region of ROS1. RT-PCR analyses confirmed the existence of an in-frame fusion transcript with ZCCHC8 exons 1-3 joined to ROS1 exons 36-43. In addition to the ZCCHC8-ROS1 fusion, we detected a deletion in the short arm of chromosome 9, including homozygous loss of the CDKN2A/2B locus in 9p21.3 and heterozygous deletion of the HAUS6 gene in 9p22.1. The latter encodes a protein involved in faithful chromosome segregation by regulating microtubule nucleation and its deletion might be associated with the marked subclonal changes of ploidy observed in the tumor. This report adds the ZCCHC8-ROS1 fusion as oncogenic driver in GBM and supports the role of ROS1 activation in the pathogenesis of a subset of GBM. © 2016 Wiley Periodicals, Inc.

MCM3AP and POMP Mutations Cause a DNA-Repair and DNA-Damage-Signaling Defect in an Immunodeficient Child

Immunodeficiency patients with DNA repair defects exhibit radiosensitivity and proneness to leukemia/lymphoma formation. Though progress has been made in identifying the underlying mutations, in most patients the genetic basis is unknown. Two de novo mutated candidate genes, MCM3AP encoding germinal center-associated nuclear protein (GANP) and POMP encoding proteasome maturation protein (POMP), were identified by whole-exome sequencing (WES) and confirmed by Sanger sequencing in a child with complex phenotype displaying immunodeficiency, genomic instability, skin changes, and myelodysplasia. GANP was previously described to promote B-cell maturation by nuclear targeting of activation-induced cytidine deaminase (AID) and to control AID-dependent hyperrecombination. POMP is required for 20S proteasome assembly and, thus, for efficient NF-κB signaling. Patient-derived cells were characterized by impaired homologous recombination, moderate radio- and cross-linker sensitivity associated with accumulation of damage, impaired DNA damage-induced NF-κB signaling, and reduced nuclear AID levels. Complementation by wild-type (WT)-GANP normalized DNA repair and WT-POMP rescued defective NF-κB signaling. In conclusion, we identified for the first time mutations in MCM3AP and POMP in an immunodeficiency patient. These mutations lead to cooperative effects on DNA recombination and damage signaling. Digenic/polygenic mutations may constitute a novel genetic basis in immunodeficiency patients with DNA repair defects.

An integrated map of structural variation in 2,504 human genomes

Structural variants are implicated in numerous diseases and make up the majority of varying nucleotides among human genomes. Here we describe an integrated set of eight structural variant classes comprising both balanced and unbalanced variants, which we constructed using short-read DNA sequencing data and statistically phased onto haplotype blocks in 26 human populations. Analysing this set, we identify numerous gene-intersecting structural variants exhibiting population stratification and describe naturally occurring homozygous gene knockouts that suggest the dispensability of a variety of human genes. We demonstrate that structural variants are enriched on haplotypes identified by genome-wide association studies and exhibit enrichment for expression quantitative trait loci. Additionally, we uncover appreciable levels of structural variant complexity at different scales, including genic loci subject to clusters of repeated rearrangement and complex structural variants with multiple breakpoints likely to have formed through individual mutational events. Our catalogue will enhance future studies into structural variant demography, functional impact and disease association.

A cell-based model system links chromothripsis with hyperploidy

A remarkable observation emerging from recent cancer genome analyses is the identification of chromothripsis as a one-off genomic catastrophe, resulting in massive somatic DNA structural rearrangements (SRs). Largely due to lack of suitable model systems, the mechanistic basis of chromothripsis has remained elusive. We developed an integrative method termed “complex alterations after selection and transformation (CAST),” enabling efficient in vitro generation of complex DNA rearrangements including chromothripsis, using cell perturbations coupled with a strong selection barrier followed by massively parallel sequencing. We employed this methodology to characterize catastrophic SR formation processes, their temporal sequence, and their impact on gene expression and cell division. Our in vitro system uncovered a propensity of chromothripsis to occur in cells with damaged telomeres, and in particular in hyperploid cells. Analysis of primary medulloblastoma cancer genomes verified the link between hyperploidy and chromothripsis in vivo. CAST provides the foundation for mechanistic dissection of complex DNA rearrangement processes.

Pediatric T-cell lymphoblastic leukemia evolves into relapse by clonal selection, acquisition of mutations and promoter hypomethylation

Relapsed precursor T-cell acute lymphoblastic leukemia is characterized by resistance against chemotherapy and is frequently fatal. We aimed at understanding the molecular mechanisms resulting in relapse of T-cell acute lymphoblastic leukemia and analyzed 13 patients at first diagnosis, remission and relapse by whole exome sequencing, targeted ultra-deep sequencing, multiplex ligation dependent probe amplification and DNA methylation array. Compared to primary T-cell acute lymphoblastic leukemia, in relapse the number of single nucleotide variants and small insertions and deletions approximately doubled from 11.5 to 26. Targeted ultra-deep sequencing sensitively detected subclones that were selected for in relapse. The mutational pattern defined two types of relapses. While both are characterized by selection of subclones and acquisition of novel mutations, 'type 1' relapse derives from the primary leukemia whereas 'type 2' relapse originates from a common pre-leukemic ancestor. Relapse-specific changes included activation of the nucleotidase NT5C2 resulting in resistance to chemotherapy and mutations of epigenetic modulators, exemplified by SUZ12, WHSC1 and SMARCA4. While mutations present in primary leukemia and in relapse were enriched for known drivers of leukemia, relapse-specific changes revealed an association with general cancer-promoting mechanisms. This study thus identifies mechanisms that drive progression of pediatric T-cell acute lymphoblastic leukemia to relapse and may explain the characteristic treatment resistance of this condition.

Assembly and diploid architecture of an individual human genome via single-molecule technologies

We present the first comprehensive analysis of a diploid human genome that combines single-molecule sequencing with single-molecule genome maps. Our hybrid assembly markedly improves upon the contiguity observed from traditional shotgun sequencing approaches, with scaffold N50 values approaching 30 Mb, and we identified complex structural variants (SVs) missed by other high-throughput approaches. Furthermore, by combining Illumina short-read data with long reads, we phased both single-nucleotide variants and SVs, generating haplotypes with over 99% consistency with previous trio-based studies. Our work shows that it is now possible to integrate single-molecule and high-throughput sequence data to generate de novo assembled genomes that approach reference quality.

Genomics and drug profiling of fatal TCF3-HLF-positive acute lymphoblastic leukemia identifies recurrent mutation patterns and therapeutic options

TCF3-HLF-fusion positive acute lymphoblastic leukemia (ALL) is currently incurable. Employing an integrated approach, we uncovered distinct mutation, gene expression, and drug response profiles in TCF3-HLF-positive and treatment-responsive TCF3-PBX1-positive ALL. Recurrent intragenic deletions of PAX5 or VPREB1 were identified in constellation with TCF3-HLF. Moreover somatic mutations in the non-translocated allele of TCF3 and a reduction of PAX5 gene dosage in TCF3-HLF ALL suggest cooperation within a restricted genetic context. The enrichment for stem cell and myeloid features in the TCF3-HLF signature may reflect reprogramming by TCF3-HLF of a lymphoid-committed cell of origin towards a hybrid, drug-resistant hematopoietic state. Drug response profiling of matched patient-derived xenografts revealed a distinct profile for TCF3-HLF ALL with resistance to conventional chemotherapeutics, but sensitivity towards glucocorticoids, anthracyclines and agents in clinical development. Striking on-target sensitivity was achieved with the BCL2-specific inhibitor venetoclax (ABT-199). This integrated approach thus provides alternative treatment options for this deadly disease.

Analysis of deletion breakpoints from 1,092 humans reveals details of mutation mechanisms

Investigating genomic structural variants at basepair resolution is crucial for understanding their formation mechanisms. We identify and analyze 8,943 deletion breakpoints in 1,092 samples from the 1000 Genomes Project. We find breakpoints have more nearby SNPs and indels than the genomic average, likely a consequence of relaxed selection. By investigating the correlation of breakpoints with DNA methylation, Hi-C interactions, and histone marks and the substitution patterns of nucleotides near them, we find that breakpoints with the signature of non-allelic homologous recombination (NAHR) are associated with open chromatin. We hypothesize that some NAHR deletions occur without DNA replication and cell division, in embryonic and germline cells. In contrast, breakpoints associated with non-homologous (NH) mechanisms often have sequence micro-insertions, templated from later replicating genomic sites, spaced at two characteristic distances from the breakpoint. These micro-insertions are consistent with template-switching events and suggest a particular spatiotemporal configuration for DNA during the events.

A common microdeletion affecting a hippocampus- and amygdala-specific isoform of tryptophan hydroxylase 2 is not associated with affective disorders

OBJECTIVES

Copy number variants (CNVs) have been shown to affect susceptibility for neuropsychiatric disorders. To date, studies implicating the serotonergic system in complex conditions have just focused on single nucleotide polymorphisms (SNPs). We therefore sought to identify novel common genetic copy number polymorphisms affecting genes of the serotonergic system, and to assess their putative role in bipolar affective disorder (BPAD) and major depressive disorder (MDD).

METHODS

A selection of 41 genes of the serotonergic system encoding receptors, the serotonin transporter, metabolic enzymes and chaperones were investigated using a paired-end mapping (PEM) approach on next-generation sequencing data from the pilot project of the 1000 Genomes Project. For association testing, 593 patients with MDD, 1,145 patients with BPAD, and 1,738 healthy controls were included in the study.

RESULTS

PEM led to the identification of a microdeletion in the gene encoding tryptophan hydroxylase 2 (TPH2), affecting an amygdala- and hippocampus-specific isoform. It was not associated with BPAD or MDD using a case-control association approach.

CONCLUSIONS

We did not find evidence for a role of the TPH2 microdeletion in the pathoetiology of affective disorders. Further studies examining its putative role in behavioral traits regulated by the limbic system are warranted.

Enhancer hijacking activates GFI1 family oncogenes in medulloblastoma

Medulloblastoma is a highly malignant paediatric brain tumour currently treated with a combination of surgery, radiation and chemotherapy, posing a considerable burden of toxicity to the developing child. Genomics has illuminated the extensive intertumoral heterogeneity of medulloblastoma, identifying four distinct molecular subgroups. Group 3 and group 4 subgroup medulloblastomas account for most paediatric cases; yet, oncogenic drivers for these subtypes remain largely unidentified. Here we describe a series of prevalent, highly disparate genomic structural variants, restricted to groups 3 and 4, resulting in specific and mutually exclusive activation of the growth factor independent 1 family proto-oncogenes, GFI1 and GFI1B. Somatic structural variants juxtapose GFI1 or GFI1B coding sequences proximal to active enhancer elements, including super-enhancers, instigating oncogenic activity. Our results, supported by evidence from mouse models, identify GFI1 and GFI1B as prominent medulloblastoma oncogenes and implicate 'enhancer hijacking' as an efficient mechanism driving oncogene activation in a childhood cancer.

SPAG7 is a candidate gene for the periodic fever, aphthous stomatitis, pharyngitis and adenopathy (PFAPA) syndrome

Periodic fever, aphthous stomatitis, pharyngitis and adenopathy (PFAPA) syndrome is an auto-inflammatory disease for which a genetic basis has been postulated. Nevertheless, in contrast to the other periodic fever syndromes, no candidate genes have yet been identified. By cloning, following long insert size paired-end sequencing, of a de novo chromosomal translocation t(10;17)(q11.2;p13) in a patient with typical PFAPA syndrome lacking mutations in genes associated with other periodic fever syndromes we identified SPAG7 as a candidate gene for PFAPA. SPAG7 protein is expressed in tissues affected by PFAPA and has been functionally linked to antiviral and inflammatory responses. Haploinsufficiency of SPAG7 due to a microdeletion at the translocation breakpoint leading to loss of exons 2-7 from one allele was associated with PFAPA in the index. Sequence analyses of SPAG7 in additional patients with PFAPA point to genetic heterogeneity or alternative mechanisms of SPAG7 deregulation, such as somatic or epigenetic changes.

Recurrent somatic alterations of FGFR1 and NTRK2 in pilocytic astrocytoma

Pilocytic astrocytoma, the most common childhood brain tumor, is typically associated with mitogen-activated protein kinase (MAPK) pathway alterations. Surgically inaccessible midline tumors are therapeutically challenging, showing sustained tendency for progression and often becoming a chronic disease with substantial morbidities. Here we describe whole-genome sequencing of 96 pilocytic astrocytomas, with matched RNA sequencing (n = 73), conducted by the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. We identified recurrent activating mutations in FGFR1 and PTPN11 and new NTRK2 fusion genes in non-cerebellar tumors. New BRAF-activating changes were also observed. MAPK pathway alterations affected all tumors analyzed, with no other significant mutations identified, indicating that pilocytic astrocytoma is predominantly a single-pathway disease. Notably, we identified the same FGFR1 mutations in a subset of H3F3A-mutated pediatric glioblastoma with additional alterations in the NF1 gene. Our findings thus identify new potential therapeutic targets in distinct subsets of pilocytic astrocytoma and childhood glioblastoma.

Integrative genomic analyses reveal an androgen-driven somatic alteration landscape in early-onset prostate cancer

Early-onset prostate cancer (EO-PCA) represents the earliest clinical manifestation of prostate cancer. To compare the genomic alteration landscapes of EO-PCA with "classical" (elderly-onset) PCA, we performed deep sequencing-based genomics analyses in 11 tumors diagnosed at young age, and pursued comparative assessments with seven elderly-onset PCA genomes. Remarkable age-related differences in structural rearrangement (SR) formation became evident, suggesting distinct disease pathomechanisms. Whereas EO-PCAs harbored a prevalence of balanced SRs, with a specific abundance of androgen-regulated ETS gene fusions including TMPRSS2:ERG, elderly-onset PCAs displayed primarily non-androgen-associated SRs. Data from a validation cohort of > 10,000 patients showed age-dependent androgen receptor levels and a prevalence of SRs affecting androgen-regulated genes, further substantiating the activity of a characteristic "androgen-type" pathomechanism in EO-PCA.

Impact of genomic structural variation in Drosophila melanogaster based on population-scale sequencing

Genomic structural variation (SV) is a major determinant for phenotypic variation. Although it has been extensively studied in humans, the nucleotide resolution structure of SVs within the widely used model organism Drosophila remains unknown. We report a highly accurate, densely validated map of unbalanced SVs comprising 8962 deletions and 916 tandem duplications in 39 lines derived from short-read DNA sequencing in a natural population (the “Drosophila melanogaster Genetic Reference Panel,” DGRP). Most SVs (>90%) were inferred at nucleotide resolution, and a large fraction was genotyped across all samples. Comprehensive analyses of SV formation mechanisms using the short-read data revealed an abundance of SVs formed by mobile element and nonhomologous end-joining-mediated rearrangements, and clustering of variants into SV hotspots. We further observed a strong depletion of SVs overlapping genes, which, along with population genetics analyses, suggests that these SVs are often deleterious. We inferred several gene fusion events also highlighting the potential role of SVs in the generation of novel protein products. Expression quantitative trait locus (eQTL) mapping revealed the functional impact of our high-resolution SV map, with quantifiable effects at >100 genic loci. Our map represents a resource for population-level studies of SVs in an important model organism.

DELLY: structural variant discovery by integrated paired-end and split-read analysis

MOTIVATION

The discovery of genomic structural variants (SVs) at high sensitivity and specificity is an essential requirement for characterizing naturally occurring variation and for understanding pathological somatic rearrangements in personal genome sequencing data. Of particular interest are integrated methods that accurately identify simple and complex rearrangements in heterogeneous sequencing datasets at single-nucleotide resolution, as an optimal basis for investigating the formation mechanisms and functional consequences of SVs.

RESULTS

We have developed an SV discovery method, called DELLY, that integrates short insert paired-ends, long-range mate-pairs and split-read alignments to accurately delineate genomic rearrangements at single-nucleotide resolution. DELLY is suitable for detecting copy-number variable deletion and tandem duplication events as well as balanced rearrangements such as inversions or reciprocal translocations. DELLY, thus, enables to ascertain the full spectrum of genomic rearrangements, including complex events. On simulated data, DELLY compares favorably to other SV prediction methods across a wide range of sequencing parameters. On real data, DELLY reliably uncovers SVs from the 1000 Genomes Project and cancer genomes, and validation experiments of randomly selected deletion loci show a high specificity.

AVAILABILITY

DELLY is available at www.korbel.embl.de/software.html

CONTACT

tobias.rausch@embl.de.

Subgroup-specific structural variation across 1,000 medulloblastoma genomes

Medulloblastoma, the most common malignant paediatric brain tumour, is currently treated with nonspecific cytotoxic therapies including surgery, whole-brain radiation, and aggressive chemotherapy. As medulloblastoma exhibits marked intertumoural heterogeneity, with at least four distinct molecular variants, previous attempts to identify targets for therapy have been underpowered because of small samples sizes. Here we report somatic copy number aberrations (SCNAs) in 1,087 unique medulloblastomas. SCNAs are common in medulloblastoma, and are predominantly subgroup-enriched. The most common region of focal copy number gain is a tandem duplication of SNCAIP, a gene associated with Parkinson's disease, which is exquisitely restricted to Group 4α. Recurrent translocations of PVT1, including PVT1-MYC and PVT1-NDRG1, that arise through chromothripsis are restricted to Group 3. Numerous targetable SCNAs, including recurrent events targeting TGF-β signalling in Group 3, and NF-κB signalling in Group 4, suggest future avenues for rational, targeted therapy.

Dissecting the genomic complexity underlying medulloblastoma

Medulloblastoma is an aggressively growing tumour, arising in the cerebellum or medulla/brain stem. It is the most common malignant brain tumour in children, and shows tremendous biological and clinical heterogeneity. Despite recent treatment advances, approximately 40% of children experience tumour recurrence, and 30% will die from their disease. Those who survive often have a significantly reduced quality of life. Four tumour subgroups with distinct clinical, biological and genetic profiles are currently identified. WNT tumours, showing activated wingless pathway signalling, carry a favourable prognosis under current treatment regimens. SHH tumours show hedgehog pathway activation, and have an intermediate prognosis. Group 3 and 4 tumours are molecularly less well characterized, and also present the greatest clinical challenges. The full repertoire of genetic events driving this distinction, however, remains unclear. Here we describe an integrative deep-sequencing analysis of 125 tumour-normal pairs, conducted as part of the International Cancer Genome Consortium (ICGC) PedBrain Tumor Project. Tetraploidy was identified as a frequent early event in Group 3 and 4 tumours, and a positive correlation between patient age and mutation rate was observed. Several recurrent mutations were identified, both in known medulloblastoma-related genes (CTNNB1, PTCH1, MLL2, SMARCA4) and in genes not previously linked to this tumour (DDX3X, CTDNEP1, KDM6A, TBR1), often in subgroup-specific patterns. RNA sequencing confirmed these alterations, and revealed the expression of what are, to our knowledge, the first medulloblastoma fusion genes identified. Chromatin modifiers were frequently altered across all subgroups. These findings enhance our understanding of the genomic complexity and heterogeneity underlying medulloblastoma, and provide several potential targets for new therapeutics, especially for Group 3 and 4 patients.

An improved protocol for sequencing of repetitive genomic regions and structural variations using mutagenesis and next generation sequencing

The rise of Next Generation Sequencing (NGS) technologies has transformed de novo genome sequencing into an accessible research tool, but obtaining high quality eukaryotic genome assemblies remains a challenge, mostly due to the abundance of repetitive elements. These also make it difficult to study nucleotide polymorphism in repetitive regions, including certain types of structural variations. One solution proposed for resolving such regions is Sequence Assembly aided by Mutagenesis (SAM), which relies on the fact that introducing enough random mutations breaks the repetitive structure, making assembly possible. Sequencing many different mutated copies permits the sequence of the repetitive region to be inferred by consensus methods. However, this approach relies on molecular cloning in order to isolate and amplify individual mutant copies, making it hard to scale-up the approach for use in conjunction with high-throughput sequencing technologies. To address this problem, we propose NG-SAM, a modified version of the SAM protocol that relies on PCR and dilution steps only, coupled to a NGS workflow. NG-SAM therefore has the potential to be scaled-up, e.g. using emerging microfluidics technologies. We built a realistic simulation pipeline to study the feasibility of NG-SAM, and our results suggest that under appropriate experimental conditions the approach might be successfully put into practice. Moreover, our simulations suggest that NG-SAM is capable of reconstructing robustly a wide range of potential target sequences of varying lengths and repetitive structures.

Genome sequencing of pediatric medulloblastoma links catastrophic DNA rearrangements with TP53 mutations

Genomic rearrangements are thought to occur progressively during tumor development. Recent findings, however, suggest an alternative mechanism, involving massive chromosome rearrangements in a one-step catastrophic event termed chromothripsis. We report the whole-genome sequencing-based analysis of a Sonic-Hedgehog medulloblastoma (SHH-MB) brain tumor from a patient with a germline TP53 mutation (Li-Fraumeni syndrome), uncovering massive, complex chromosome rearrangements. Integrating TP53 status with microarray and deep sequencing-based DNA rearrangement data in additional patients reveals a striking association between TP53 mutation and chromothripsis in SHH-MBs. Analysis of additional tumor entities substantiates a link between TP53 mutation and chromothripsis, and indicates a context-specific role for p53 in catastrophic DNA rearrangements. Among these, we observed a strong association between somatic TP53 mutations and chromothripsis in acute myeloid leukemia. These findings connect p53 status and chromothripsis in specific tumor types, providing a genetic basis for understanding particularly aggressive subtypes of cancer.

A comprehensive map of mobile element insertion polymorphisms in humans

As a consequence of the accumulation of insertion events over evolutionary time, mobile elements now comprise nearly half of the human genome. The Alu, L1, and SVA mobile element families are still duplicating, generating variation between individual genomes. Mobile element insertions (MEI) have been identified as causes for genetic diseases, including hemophilia, neurofibromatosis, and various cancers. Here we present a comprehensive map of 7,380 MEI polymorphisms from the 1000 Genomes Project whole-genome sequencing data of 185 samples in three major populations detected with two detection methods. This catalog enables us to systematically study mutation rates, population segregation, genomic distribution, and functional properties of MEI polymorphisms and to compare MEI to SNP variation from the same individuals. Population allele frequencies of MEI and SNPs are described, broadly, by the same neutral ancestral processes despite vastly different mutation mechanisms and rates, except in coding regions where MEI are virtually absent, presumably due to strong negative selection. A direct comparison of MEI and SNP diversity levels suggests a differential mobile element insertion rate among populations.

Molecular mechanisms for activation of the agouti-related protein and stimulation of appetite

OBJECTIVE

The agouti-related protein (Agrp) is a powerful orexigenic peptide, but little is known about its transcriptional regulation. The objective of this study was to determine molecular mechanisms for the activation of hypothalamic Agrp and identify compounds that stimulate appetite.

RESEARCH DESIGN AND METHODS

We used promoter analyses methods, hypothalamic cell culture and transfection, immunohistochemistry, luciferase-expressing transgenic mice, in vivo bioluminescence, anitisense RNA, mouse feeding studies, indirect calorimetry, real-time PCR, and Western blots.

RESULTS

We found that the Krüppel-like factor 4 (Klf4) is a potent activator of Agrp by binding to a specific CACCC-box in its minimal promoter. We also found that an extract of tarragon, termed PMI-5011, activated hypothalamic Klf4 and Agrp. In vivo, PMI-5011 increased Agrp promoter activity in luciferase-expressing transgenic mice, increased hypothalamic Klf4 and Agrp expression, increased hypothalamic Orexin and melanin-concentrating hormone, increased food intake, reduced circulating insulin and leptin levels, attenuated energy expenditure, and enhanced body weight but only when using a high-fat diet.

CONCLUSIONS

These data show that Klf4 augmented hypothalamic Agrp by binding to a specific CACCC-box onto its minimal promoter. In addition, the tarragon extract PMI-5011 activated Klf4 and orexigenic neuropeptides and reduced peripheral insulin and leptin levels leading to positive energy balance.

Systematic inference of copy-number genotypes from personal genome sequencing data reveals extensive olfactory receptor gene content diversity

Copy-number variations (CNVs) are widespread in the human genome, but comprehensive assignments of integer locus copy-numbers (i.e., copy-number genotypes) that, for example, enable discrimination of homozygous from heterozygous CNVs, have remained challenging. Here we present CopySeq, a novel computational approach with an underlying statistical framework that analyzes the depth-of-coverage of high-throughput DNA sequencing reads, and can incorporate paired-end and breakpoint junction analysis based CNV-analysis approaches, to infer locus copy-number genotypes. We benchmarked CopySeq by genotyping 500 chromosome 1 CNV regions in 150 personal genomes sequenced at low-coverage. The assessed copy-number genotypes were highly concordant with our performed qPCR experiments (Pearson correlation coefficient 0.94), and with the published results of two microarray platforms (95–99% concordance). We further demonstrated the utility of CopySeq for analyzing gene regions enriched for segmental duplications by comprehensively inferring copy-number genotypes in the CNV-enriched >800 olfactory receptor (OR) human gene and pseudogene loci. CopySeq revealed that OR loci display an extensive range of locus copy-numbers across individuals, with zero to two copies in some OR loci, and two to nine copies in others. Among genetic variants affecting OR loci we identified deleterious variants including CNVs and SNPs affecting ∼15% and ∼20% of the human OR gene repertoire, respectively, implying that genetic variants with a possible impact on smell perception are widespread. Finally, we found that for several OR loci the reference genome appears to represent a minor-frequency variant, implying a necessary revision of the OR repertoire for future functional studies. CopySeq can ascertain genomic structural variation in specific gene families as well as at a genome-wide scale, where it may enable the quantitative evaluation of CNVs in genome-wide association studies involving high-throughput sequencing.

Human individual genome sequencing has recently become affordable, enabling highly detailed genetic sequence comparisons. While the identification and genotyping of single-nucleotide polymorphisms has already been successfully established for different sequencing platforms, the detection, quantification and genotyping of large-scale copy-number variants (CNVs), i.e., losses or gains of long genomic segments, has remained challenging. We present a computational approach that enables detecting CNVs in sequencing data and accurately identifies the actual copy-number at which DNA segments of interest occur in an individual genome. This approach enabled us to obtain novel insights into the largest human gene family – the olfactory receptors (ORs) – involved in smell perception. While previous studies reported an abundance of CNVs in ORs, our approach enabled us to globally identify absolute differences in OR gene counts that exist between humans. While several OR genes have very high gene counts, other ORs are found only once or are missing entirely in some individuals. The latter have a particularly high probability of influencing individual differences in the perception of smell, a question that future experimental efforts can now address. Furthermore, we observed differences in OR gene counts between populations, pointing at ORs that might contribute to population-specific differences in smell.

Case-control association study of TGOLN2 in attempted suicide

Family, twin, and adoption studies provide convincing evidence for a genetic contribution to suicidal behavior. The heritability for suicidal behavior depends in part on the transmission of psychiatric disorders, such as mood disorders and substance use disorders, but is also partly independent of them. Three linkage studies using the attempted suicide phenotype in pedigrees with bipolar disorder, major depression, or alcoholism have provided consistent evidence that 2p11-12 harbors a susceptibility gene for attempted suicide. A microarray expression study using postmortem brain samples has implicated a gene from the 2p11-12 candidate region, the trans-Golgi network protein 2 (TGOLN2) gene, as being consistently up-regulated in suicide cases as compared to controls. Here, we present a TGOLN2 case-control association study using nine single nucleotide polymorphisms (SNPs). These nine SNPs, which include seven tag SNPs and two coding SNPs, have been genotyped in 517 mood disorder subjects with a history of attempted suicide and 515 normal controls. Allelic and genotypic analyses of the case-control sample did not provide evidence for association with the attempted suicide phenotype. Eight of the nine SNPs provided supportive evidence for association (P-values ranging from 0.008 to 0.03) when we compared the attempted suicide cases with a history of alcoholism to the attempted suicide cases without a history of alcoholism. However, this association finding was not replicated in an independent sample. Taken together, these analyses do not provide support for the hypothesis that common genetic variation in TGOLN2 contributes significantly to the risk for attempted suicide in subjects with major mood disorders.

Functional identification of the promoter of SLC4A5, a gene associated with cardiovascular and metabolic phenotypes in the HERITAGE Family Study

The sodium bicarbonate cotransporter gene SLC4A5, associated earlier with cardiovascular phenotypes, was tested for associations in the HERITAGE Family Study, and possible mechanisms were investigated. Twelve tag-single nucleotide polymorphisms (SNPs) covering the SLC4A5 gene were analyzed in 276 Black and 503 White healthy, sedentary subjects. Associations were tested using a variance components-based (QTDT) method with data adjusted for age, sex and body size. In Whites, rs6731545 and rs7571842 were significantly associated with resting and submaximal exercise pulse pressure (PP) (0.0004 <P<0.0007 and 0.002<P<0.003 respectively). Additionally, rs6731545 was associated with submaximal-exercise systolic blood pressure (SBP) and rate pressure product (P=0.002, both). New associations between rs6731545 and submaximal-exercise VO(2) (P=0.003), rs7587117 and rs7571842 and VCO(2) (0.0005<P<0.0009) and rs828853 and VE (P=0.002) were found. All these associations had a FDR<0.05. Single-marker associations were confirmed in haplotype analyses. Using in silico analysis, evidence was found for a main and an alternative promoter for SLC4A5. Specific promoter activity was experimentally confirmed using reporter constructs targeting both promoters in three physiologically relevant cell lines. Re-sequencing of 32 individuals having opposite homozygotes for rs7571842 and rs6731545 and exhibiting significantly different phenotypes showed no SNPs in the alternative promoter and no differences between the groups with SNPs in the main promoter. Also, of all known SLC4A5-coding SNPs, only one synonymous SNP was detected. Summarizing, the observed associations with resting and submaximal-exercise cardiovascular and metabolic traits in the HERITAGE Family Study are likely due to neither variation in the promoter nor known coding SNPs of SLC4A5.

KIF5B gene sequence variation and response of cardiac stroke volume to regular exercise

A genome-wide linkage scan for endurance training-induced changes in stroke volume detected a quantitative trait locus on chromosome 10p11 in white families of the HERITAGE Family Study. Dense microsatellite mapping narrowed down the linkage region to a 7 Mb area containing 16 known and 14 predicted genes. Association analyses with 90 single nucleotide polymorphisms (SNPs) provided suggestive evidence (P values from 0.03 to 0.06) for association in the kinesin heavy chain (KIF5B) gene locus in the whole cohort. The associations at the KIF5B locus were stronger (P values from 0.001 to 0.008) when the analyses were performed on linkage-informative families only (family-specific logarithm of the odds ratio scores >0.025 at peak linkage location). Resequencing the coding and regulatory regions of KIF5B revealed no new exonic SNPs. However, the putative promoter region was particularly polymorphic, containing eight SNPs with at least 5% minor allele frequency within 1850 bp upstream of the start codon. Functional analyses using promoter haplotype reporter constructs led to the identification of sequence variants that had significant effects on KIF5B promoter activity. Analogous inhibition and overexpression experiments showed that changes in KIF5B expression alter mitochondrial localization and biogenesis in a manner that could affect the ability of the heart to adjust to regular exercise. Our data suggest that KIF5B is a strong candidate gene for the response of stroke volume to regular exercise. Furthermore, training-induced changes in submaximal exercise stroke volume may be due to mitochondrial function and variation in KIF5B expression as determined by functional SNPs in its promoter.

Circadian expression of genes regulating food intake

OBJECTIVE

The Agouti-related protein (AgRP), neuropeptide Y (NPY), proopiomelanocortin (POMC), cocaine and amphetamine-regulated transcript (CART), Orexin, melanin concentrating hormone (MCH), leptin, and its hypothalamic receptor (LR) are key regulators of food intake and energy homeostasis. In the present study, we examined the circadian expression profiles of these genes.

RESEARCH METHODS AND PROCEDURES

We used quantitative reverse transcriptase-polymerase chain reaction (RT-PCR) to measure mRNA levels, spectral analysis to evaluate periodicity, and correlation analysis to examine for associations with diurnal food intake.

RESULTS

No gene in particular stood out as a strong candidate, but the overall circadian expression profiles of leptin and its hypothalamic receptor came close to statistically and graphically resembling the diurnal feeding behavior of mice. In mathematical terms, adrenal AgRP exhibited strong circadian expression and had the highest correlation with food intake, followed by leptin. Yet its highest point of expression occurred 8 hours after nocturnal food intake had peaked, suggesting that adrenal AgRP could not play a direct role in the initiation of nocturnal feeding; neither did hypothalamic AgRP, NPY, POMC, CART, Orexin, or MCH.

DISCUSSION

These data show that ad libitum feeding in mice is influenced by complex central and peripheral circuits involving orexigenic and anorectic agents of which leptin and its hypothalamic receptor could play more prevalent roles.

Selective tissue uptake of agouti-related protein(82-131) and its modulation by fasting

The blood concentration of agouti-related protein (AgRP), a protein related to hyperphagia and obesity, is increased in obese human and fasted lean subjects. Because there is no saturable transport system at the blood-brain barrier for circulating AgRP to reach its central nervous system target, uptake of AgRP by peripheral organs might be physiologically meaningful. Using the biologically active fragment AgRP(82-131), we determined the pharmacokinetics of its radioactively labeled tracer after iv bolus injection and compared it with that of the vascular marker albumin. AgRP enters peripheral organs at different influx rates, all of which were higher than into brain and spinal cord. At 10 min after iv injection, the radioactivity recovered in the liver, which had the fastest influx rate for AgRP, represented intact (125)I-AgRP. The adrenal gland had a moderately fast uptake (but the highest initial volume of distribution), followed by the heart, lungs, and skeletal muscle. By comparison, epididymal fat, testis, and pancreas had low permeability to AgRP. Saturation of influx was determined by coadministration of excess unlabeled AgRP and was shown to be present in the liver and adrenal gland. The influx rate and initial volume of distribution did not show a linear correlation with vascular permeability or regional blood flow. AgRP uptake by the liver and epididymal fat was significantly increased by overnight fasting, whereas that by the adrenal gland was significantly decreased in fasted mice. Thus, the differential uptake of AgRP by peripheral organs could be a regulated process that is modulated by food deprivation.

The agouti-related protein and its role in energy homeostasis

The melanocortin system plays an important role in the regulation of energy homeostasis. The Agouti-related protein (AGRP) is a natural antagonist of the action of alpha-melanocyte stimulating hormone (alpha-MSH) at the melanocortin receptors (MCR). AGRP is upregulated by fasting while intracerebroventricular injections of synthetic AGRP lead to increased appetite and food intake. Transgenic mice overexpressing AGRP are also hyperphagic and eventually become obese. AGRP is, therefore, a significant regulator of energy balance and a candidate gene for human fatness. Indeed, humans with common single nucleotide polymorphisms (SNPs) in the promoter or the coding region are leaner and resistant to late-onset obesity than wild-type individuals. AGRP is also expressed in the periphery. Recent studies show that AGRP in the adrenal gland is upregulated by fasting as much as it is in the hypothalamus. These data open up the possibility for a wider role by AGRP not only in food intake but also in the regulation of energy balance through its actions on peripheral tissues. This review summarizes recent advances in the biochemical and physiological properties of AGRP in an effort to enhance our understanding of the role this powerful neuropeptide plays in mammalian energy homeostasis.

The Th2 cell cytokines IL-4 and IL-13 regulate found in inflammatory zone 1/resistin-like molecule alpha gene expression by a STAT6 and CCAAT/enhancer-binding protein-dependent mechanism

The onset of allergic inflammation in the lung is driven by a complex genetic program. This study shows that found in inflammatory zone (FIZZ)1 and FIZZ2, but not FIZZ3, gene expression was up-regulated 6 h after Ag challenge in a mouse model of acute pulmonary inflammation. Induction of both genes was abolished in allergen-challenged STAT6-deficient mice. FIZZ1, but not FIZZ2, mRNA was up-regulated upon incubation of the myeloid cell line BMnot with IL-4. The promoter region of FIZZ1 contains functional binding sites for STAT6 and C/EBP. FIZZ1 promoter reporter gene constructs responded to IL-4 and IL-13 stimulation in transiently transfected cells. Point mutations in the STAT6 or the C/EBP site led to loss of cytokine responsiveness indicating that IL-4-mediated induction of murine FIZZ1 is orchestrated by the coordinate action of STAT6 and C/EBP. It is concluded that the expression of the genes encoding FIZZ1 and FIZZ2, but not FIZZ3, is induced in allergen-challenged lungs in a STAT6-dependent fashion. STAT6 directly regulates IL-4- and IL-13-triggered induction of FIZZ1 expression at the transcriptional level by cooperation with C/EBP. Induction of FIZZ2 gene expression most likely occurs independent of a direct effect by these cytokines and may be due to indirect STAT6-driven mechanisms.

Inhibition of interleukin-4- and CD40-induced IgE germline gene promoter activity by 2'-aminoethoxy-modified triplex-forming oligonucleotides

Elevated levels of IgE are intimately associated with a number of allergic diseases, such as allergic rhinitis or asthma. Therefore, prevention of IgE production in human B-cells represents an attractive therapeutic target. IL-4-induced IgE germline gene transcription represents a crucial early step during IgE isotype switch differentiation. Gene induction is orchestrated by the coordinated action of the transcription factors STAT6 (signal transducer and activator of transcription), NF-kappaB, PU.1, and C/EBP. This study shows that 2'-aminoethoxy-modified oligonucleotides, which partially overlap with the STAT6 and the adjacent PU.1/NF-kappaB binding site, inhibit DNA binding of all three proteins with high affinity in a dose- and time-dependent fashion in vitro. Loss of protein binding correlated strongly with increasing DNA triplex formation. Importantly, the oligomers also effectively displaced pre-bound recombinant NF-kappaB p50 from double-stranded DNA in vitro. Functionally, the oligonucleotides led to a selective inhibition of IL-4-induced reporter gene activity from a construct driven by the IgE germline gene promoter in human B-cells. These data confirm the critical role of this cytokine-responsive regulatory region in IgE germline gene induction and further support the concept of specific modulation of gene expression by DNA triplex formation induced with chemically modified oligonucleotides.

Functional synergism of STAT6 with either NF-kappa B or PU.1 to mediate IL-4-induced activation of IgE germline gene transcription

Ig heavy chain class switching to IgE is directed by IL-4 and IL-13 by inducing transcription from the IgE germline promoter. A crucial transcription factor in this process is STAT6, which binds to a specific DNA element upon cytokine activation. In this paper it is shown that the B cell- and monocyte-specific factor PU.1 interacts with a closely spaced sequence in the human IgE germline promoter that overlaps with a previously described binding site for NF kappa B/rel. The authenticity of PU.1 was demonstrated by specific competition and supershifts in EMSA experiments. In addition, in vitro translated PU.1 could interact with an oligonucleotide derived from the IgE germline promoter containing the PU.1 binding site and migrated with the same mobility compared with the complex formed with nuclear extracts. Transient transfection experiments using IgE germline promoter reporter gene constructs demonstrated that mutations affecting DNA binding of PU.1 or NF kappa B/rel had no or little effect on IL-4 inducibility of these plasmids. However, point mutations that abolished binding of both factors abrogated cytokine inducibility. No strict spacing of the STAT6 and the composite PU. 1/NF-kappa B elements is required for IL-4 induction. IL-4-induced STAT6 DNA binding was retained in PU.1-/NF kappa B/rel- double mutants. The data demonstrate that cooperation of STAT6 with at least PU.1 or NF kappa B/rel is necessary for IL-4-induced activation of IgE germline gene transcription.

Cooperation of binding sites for STAT6 and NF kappa B/rel in the IL-4-induced up-regulation of the human IgE germline promoter

Ig heavy chain class switching is directed by cytokines inducing transcription from unrearranged CH genes. Subsequently, such primed cells can undergo switch recombination to express the selected new isotype. In the case of IgE class switching, IL-4 activates the IgE germline promoter by inducing the interaction of the transcription factor STAT6 (IL-4STAT) with a responsive DNA element in the proximal region of the promoter. This study describes the characterization of two additional cis-acting elements that interact with members of the NF kappa B/rel transcription factor family in an IL-4-independent fashion. Electrophoretic mobility shift assays show that the nucleoprotein complex formed on the upstream site (NF kappa B1) contains the classical p50/p65 heterodimer. The complex on the proximal site (NF kappa B2) appears to be composed of p50 and relB. IgE germline promoter reporter gene constructs carrying point mutations in the NF kappa B2 site were largely unresponsive to IL-4 stimulation in transient transfection experiments, while plasmids with similar mutations in the NF kappa B1 site responded to cytokine stimulation better than the wild-type promoter. The NF kappa B2 effect was dependent on the presence of the STAT6 binding site, demonstrating that the NF kappa B2 motif is necessary but not sufficient for mediating cytokine up-regulation. In addition, the combination of a NF kappa B/rel binding site and the STAT6 response element conferred IL-4 inducibility to a heterologous minimal promoter, while the individual sites had no effect. The available data suggest that the NF kappa B2 nucleoprotein complex may cooperate with DNA-bound STAT6 to achieve IL-4-dependent activation of the human IgE germline gene.

Cooperation of binding sites for STAT6 and NF kappa B/rel in the IL-4-induced up-regulation of the human IgE germline promoter

Ig heavy chain class switching is directed by cytokines inducing transcription from unrearranged CH genes. Subsequently, such primed cells can undergo switch recombination to express the selected new isotype. In the case of IgE class switching, IL-4 activates the IgE germline promoter by inducing the interaction of the transcription factor STAT6 (IL-4STAT) with a responsive DNA element in the proximal region of the promoter. This study describes the characterization of two additional cis-acting elements that interact with members of the NF kappa B/rel transcription factor family in an IL-4-independent fashion. Electrophoretic mobility shift assays show that the nucleoprotein complex formed on the upstream site (NF kappa B1) contains the classical p50/p65 heterodimer. The complex on the proximal site (NF kappa B2) appears to be composed of p50 and relB. IgE germline promoter reporter gene constructs carrying point mutations in the NF kappa B2 site were largely unresponsive to IL-4 stimulation in transient transfection experiments, while plasmids with similar mutations in the NF kappa B1 site responded to cytokine stimulation better than the wild-type promoter. The NF kappa B2 effect was dependent on the presence of the STAT6 binding site, demonstrating that the NF kappa B2 motif is necessary but not sufficient for mediating cytokine up-regulation. In addition, the combination of a NF kappa B/rel binding site and the STAT6 response element conferred IL-4 inducibility to a heterologous minimal promoter, while the individual sites had no effect. The available data suggest that the NF kappa B2 nucleoprotein complex may cooperate with DNA-bound STAT6 to achieve IL-4-dependent activation of the human IgE germline gene.