Population Variation and Genetic Control of Modular Chromatin Architecture in Humans.

Chromatin state variation at gene regulatory elements is abundant across individuals, yet we understand little about the genetic basis of this variability. Here, we profiled several histone modifications, the transcription factor (TF) PU.1, RNA polymerase II, and gene expression in lymphoblastoid cell lines from 47 whole-genome sequenced individuals. We observed that distinct cis-regulatory elements exhibit coordinated chromatin variation across individuals in the form of variable chromatin modules (VCMs) at sub-Mb scale. VCMs were associated with thousands of genes and preferentially cluster within chromosomal contact domains. We mapped strong proximal and weak, yet more ubiquitous, distal-acting chromatin quantitative trait loci (cQTL) that frequently explain this variation. cQTLs were associated with molecular activity at clusters of cis-regulatory elements and mapped preferentially within TF-bound regions. We propose that local, sequence-independent chromatin variation emerges as a result of genetic perturbations in cooperative interactions between cis-regulatory elements that are located within the same genomic domain.

A dual program for translation regulation in cellular proliferation and differentiation.

A dichotomous choice for metazoan cells is between proliferation and differentiation. Measuring tRNA pools in various cell types, we found two distinct subsets, one that is induced in proliferating cells, and repressed otherwise, and another with the opposite signature. Correspondingly, we found that genes serving cell-autonomous functions and genes involved in multicellularity obey distinct codon usage. Proliferation-induced and differentiation-induced tRNAs often carry anticodons that correspond to the codons enriched among the cell-autonomous and the multicellularity genes, respectively. Because mRNAs of cell-autonomous genes are induced in proliferation and cancer in particular, the concomitant induction of their codon-enriched tRNAs suggests coordination between transcription and translation. Histone modifications indeed change similarly in the vicinity of cell-autonomous genes and their corresponding tRNAs, and in multicellularity genes and their tRNAs, suggesting the existence of transcriptional programs coordinating tRNA supply and demand. Hence, we describe the existence of two distinct translation programs that operate during proliferation and differentiation.

Germline Elongator mutations in Sonic Hedgehog medulloblastoma.

Cancer genomics has revealed many genes and core molecular processes that contribute to human malignancies, but the genetic and molecular bases of many rare cancers remains unclear. Genetic predisposition accounts for 5 to 10% of cancer diagnoses in children1,2, and genetic events that cooperate with known somatic driver events are poorly understood. Pathogenic germline variants in established cancer predisposition genes have been recently identified in 5% of patients with the malignant brain tumour medulloblastoma3. Here, by analysing all protein-coding genes, we identify and replicate rare germline loss-of-function variants across ELP1 in 14% of paediatric patients with the medulloblastoma subgroup Sonic Hedgehog (MBSHH). ELP1 was the most common medulloblastoma predisposition gene and increased the prevalence of genetic predisposition to 40% among paediatric patients with MBSHH. Parent-offspring and pedigree analyses identified two families with a history of paediatric medulloblastoma. ELP1-associated medulloblastomas were restricted to the molecular SHHα subtype4 and characterized by universal biallelic inactivation of ELP1 owing to somatic loss of chromosome arm 9q. Most ELP1-associated medulloblastomas also exhibited somatic alterations in PTCH1, which suggests that germline ELP1 loss-of-function variants predispose individuals to tumour development in combination with constitutive activation of SHH signalling. ELP1 is the largest subunit of the evolutionarily conserved Elongator complex, which catalyses translational elongation through tRNA modifications at the wobble (U34) position5,6. Tumours from patients with ELP1-associated MBSHH were characterized by a destabilized Elongator complex, loss of Elongator-dependent tRNA modifications, codon-dependent translational reprogramming, and induction of the unfolded protein response, consistent with loss of protein homeostasis due to Elongator deficiency in model systems7-9. Thus, genetic predisposition to proteome instability may be a determinant in the pathogenesis of paediatric brain cancers. These results support investigation of the role of protein homeostasis in other cancer types and potential for therapeutic interference.

Genomic basis for RNA alterations in cancer.

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

The molecular landscape of ETMR at diagnosis and relapse.

Embryonal tumours with multilayered rosettes (ETMRs) are aggressive paediatric embryonal brain tumours with a universally poor prognosis1. Here we collected 193 primary ETMRs and 23 matched relapse samples to investigate the genomic landscape of this distinct tumour type. We found that patients with tumours in which the proposed driver C19MC2-4 was not amplified frequently had germline mutations in DICER1 or other microRNA-related aberrations such as somatic amplification of miR-17-92 (also known as MIR17HG). Whole-genome sequencing revealed that tumours had an overall low recurrence of single-nucleotide variants (SNVs), but showed prevalent genomic instability caused by widespread occurrence of R-loop structures. We show that R-loop-associated chromosomal instability can be induced by the loss of DICER1 function. Comparison of primary tumours and matched relapse samples showed a strong conservation of structural variants, but low conservation of SNVs. Moreover, many newly acquired SNVs are associated with a mutational signature related to cisplatin treatment. Finally, we show that targeting R-loops with topoisomerase and PARP inhibitors might be an effective treatment strategy for this deadly disease.

The whole-genome landscape of medulloblastoma subtypes.

Current therapies for medulloblastoma, a highly malignant childhood brain tumour, impose debilitating effects on the developing child, and highlight the need for molecularly targeted treatments with reduced toxicity. Previous studies have been unable to identify the full spectrum of driver genes and molecular processes that operate in medulloblastoma subgroups. Here we analyse the somatic landscape across 491 sequenced medulloblastoma samples and the molecular heterogeneity among 1,256 epigenetically analysed cases, and identify subgroup-specific driver alterations that include previously undiscovered actionable targets. Driver mutations were confidently assigned to most patients belonging to Group 3 and Group 4 medulloblastoma subgroups, greatly enhancing previous knowledge. New molecular subtypes were differentially enriched for specific driver events, including hotspot in-frame insertions that target KBTBD4 and 'enhancer hijacking' events that activate PRDM6. Thus, the application of integrative genomics to an extensive cohort of clinical samples derived from a single childhood cancer entity revealed a series of cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for the treatment of patients with medulloblastoma.

Cancer risk and tumour spectrum in 172 patients with a germline SUFU pathogenic variation: a collaborative study of the SIOPE Host Genome Working Group.

BACKGROUND: Little is known about risks associated with germline SUFU pathogenic variants (PVs) known as a cancer predisposition syndrome. METHODS: To study tumour risks, we have analysed data of a large cohort of 45 unpublished patients with a germline SUFU PV completed with 127 previously published patients. To reduce the ascertainment bias due to index patient selection, the risk of tumours was evaluated in relatives with SUFU PV (89 patients) using the Nelson-Aalen estimator. RESULTS: Overall, 117/172 (68%) SUFU PV carriers developed at least one tumour: medulloblastoma (MB) (86 patients), basal cell carcinoma (BCC) (25 patients), meningioma (20 patients) and gonadal tumours (11 patients). Thirty-three of them (28%) had multiple tumours. Median age at diagnosis of MB, gonadal tumour, first BCC and first meningioma were 1.5, 14, 40 and 44 years, respectively. Follow-up data were available for 160 patients (137 remained alive and 23 died). The cumulative incidence of tumours in relatives was 14.4% (95% CI 6.8 to 21.4), 18.2% (95% CI 9.7 to 25.9) and 44.1% (95% CI 29.7 to 55.5) at the age of 5, 20 and 50 years, respectively. The cumulative risk of an MB, gonadal tumour, BCC and meningioma at age 50 years was: 13.3% (95% CI 6 to 20.1), 4.6% (95% CI 0 to 9.7), 28.5% (95% CI 13.4 to 40.9) and 5.2% (95% CI 0 to 12), respectively. Sixty-four different PVs were reported across the entire SUFU gene and inherited in 73% of cases in which inheritance could be evaluated. CONCLUSION: Germline SUFU PV carriers have a life-long increased risk of tumours with a spectrum dominated by MB before the age of 5, gonadal tumours during adolescence and BCC and meningioma in adulthood, justifying fine-tuned surveillance programmes.

Active medulloblastoma enhancers reveal subgroup-specific cellular origins.

Medulloblastoma is a highly malignant paediatric brain tumour, often inflicting devastating consequences on the developing child. Genomic studies have revealed four distinct molecular subgroups with divergent biology and clinical behaviour. An understanding of the regulatory circuitry governing the transcriptional landscapes of medulloblastoma subgroups, and how this relates to their respective developmental origins, is lacking. Here, using H3K27ac and BRD4 chromatin immunoprecipitation followed by sequencing (ChIP-seq) coupled with tissue-matched DNA methylation and transcriptome data, we describe the active cis-regulatory landscape across 28 primary medulloblastoma specimens. Analysis of differentially regulated enhancers and super-enhancers reinforced inter-subgroup heterogeneity and revealed novel, clinically relevant insights into medulloblastoma biology. Computational reconstruction of core regulatory circuitry identified a master set of transcription factors, validated by ChIP-seq, that is responsible for subgroup divergence, and implicates candidate cells of origin for Group 4. Our integrated analysis of enhancer elements in a large series of primary tumour samples reveals insights into cis-regulatory architecture, unrecognized dependencies, and cellular origins.

Integrative genomics identifies the corepressor SMRT as a gatekeeper of adipogenesis through the transcription factors C/EBPß and KAISO.

The molecular role of corepressors is poorly understood. Here, we studied the transcriptional function of the corepressor SMRT during terminal adipogenesis. Genome-wide DNA-binding profiling revealed that this corepressor is predominantly located in active chromatin regions and that most distal SMRT binding events are lost after differentiation induction. Promoter-proximal tethering of SMRT in preadipocytes is primarily mediated by KAISO through the conserved TCTCGCGAGA motif. Further characterization revealed that KAISO, similar to SMRT, accelerates the cell cycle and increases fat accumulation upon knockdown, identifying KAISO as an adipogenic repressor that likely modulates the mitotic clonal expansion phase of this process. SMRT-bound promoter-distal sites tend to overlap with C/EBPß-bound regions, which become occupied by proadipogenic transcription factors after SMRT clearance. This reveals a role for SMRT in masking enhancers from proadipogenic factors in preadipocytes. Finally, we identified SMRT as an adipogenic gatekeeper as it directly fine-tunes transcription of pro- and antiadipogenic genes.

Antibiotics-induced monodominance of a novel gut bacterial order.

OBJECTIVE: The composition of the healthy human adult gut microbiome is relatively stable over prolonged periods, and representatives of the most highly abundant and prevalent species have been cultured and described. However, microbial abundances can change on perturbations, such as antibiotics intake, enabling the identification and characterisation of otherwise low abundant species. DESIGN: Analysing gut microbial time-series data, we used shotgun metagenomics to create strain level taxonomic and functional profiles. Community dynamics were modelled postintervention with a focus on conditionally rare taxa and previously unknown bacteria. RESULTS: In response to a commonly prescribed cephalosporin (ceftriaxone), we observe a strong compositional shift in one subject, in which a previously unknown species, UBorkfalki ceftriaxensis, was identified, blooming to 92% relative abundance. The genome assembly reveals that this species (1) belongs to a so far undescribed order of Firmicutes, (2) is ubiquitously present at low abundances in at least one third of adults, (3) is opportunistically growing, being ecologically similar to typical probiotic species and (4) is stably associated to healthy hosts as determined by single nucleotide variation analysis. It was the first coloniser after the antibiotic intervention that led to a long-lasting microbial community shift and likely permanent loss of nine commensals. CONCLUSION: The bloom of UB. ceftriaxensis and a subsequent one of Parabacteroides distasonis demonstrate the existence of monodominance community states in the gut. Our study points to an undiscovered wealth of low abundant but common taxa in the human gut and calls for more highly resolved longitudinal studies, in particular on ecosystem perturbations.

Risk-adapted therapy for young children with medulloblastoma (SJYC07): therapeutic and molecular outcomes from a multicentre, phase 2 trial.

BACKGROUND: Young children with medulloblastoma have a poor overall survival compared with older children, due to use of radiation-sparing therapy in young children. Radiotherapy is omitted or reduced in these young patients to spare them from debilitating long-term side-effects. We aimed to estimate event-free survival and define the molecular characteristics associated with progression-free survival in young patients with medulloblastoma using a risk-stratified treatment strategy designed to defer, reduce, or delay radiation exposure. METHODS: In this multicentre, phase 2 trial, we enrolled children younger than 3 years with newly diagnosed medulloblastoma at six centres in the USA and Australia. Children aged 3-5 years with newly diagnosed, non-metastatic medulloblastoma without any high-risk features were also eligible. Eligible patients were required to start therapy within 31 days from definitive surgery, had a Lansky performance score of at least 30, and did not receive previous radiotherapy or chemotherapy. Patients were stratified postoperatively by clinical and histological criteria into low-risk, intermediate-risk, and high-risk treatment groups. All patients received identical induction chemotherapy (methotrexate, vincristine, cisplatin, and cyclophosphamide), with high-risk patients also receiving an additional five doses of vinblastine. Induction was followed by risk-adapted consolidation therapy: low-risk patients received cyclophosphamide (1500 mg/m2 on day 1), etoposide (100 mg/m2 on days 1 and 2), and carboplatin (area under the curve 5 mg/mL per min on day 2) for two 4-week cycles; intermediate-risk patients received focal radiation therapy (54 Gy with a clinical target volume of 5 mm over 6 weeks) to the tumour bed; and high-risk patients received chemotherapy with targeted intravenous topotecan (area under the curve 120-160 ng-h/mL intravenously on days 1-5) and cyclophosphamide (600 mg/m2 intravenously on days 1-5). After consolidation, all patients received maintenance chemotherapy with cyclophosphamide, topotecan, and erlotinib. The coprimary endpoints were event-free survival and patterns of methylation profiling associated with progression-free survival. Outcome and safety analyses were per protocol (all patients who received at least one dose of induction chemotherapy); biological analyses included all patients with tissue available for methylation profiling. This trial is registered with ClinicalTrials.gov, number NCT00602667, and was closed to accrual on April 19, 2017. FINDINGS: Between Nov 27, 2007, and April 19, 2017, we enrolled 81 patients with histologically confirmed medulloblastoma. Accrual to the low-risk group was suspended after an interim analysis on Dec 2, 2015, when the 1-year event-free survival was estimated to be below the stopping rule boundary. After a median follow-up of 5·5 years (IQR 2·7-7·3), 5-year event-free survival was 31·3% (95% CI 19·3-43·3) for the whole cohort, 55·3% (95% CI 33·3-77·3) in the low-risk cohort (n=23) versus 24·6% (3·6-45·6) in the intermediate-risk cohort (n=32; hazard ratio 2·50, 95% CI 1·19-5·27; p=0·016) and 16·7% (3·4-30·0) in the high-risk cohort (n=26; 3·55, 1·66-7·59; p=0·0011; overall p=0·0021). 5-year progression-free survival by methylation subgroup was 51·1% (95% CI 34·6-67·6) in the sonic hedgehog (SHH) subgroup (n=42), 8·3% (95% CI 0·0-24·0%) in the group 3 subgroup (n=24), and 13·3% (95% CI 0·0-37·6%) in the group 4 subgroup (n=10). Within the SHH subgroup, two distinct methylation subtypes were identified and named iSHH-I and iSHH-II. 5-year progression-free survival was 27·8% (95% CI 9·0-46·6; n=21) for iSHH-I and 75·4% (55·0-95·8; n=21) for iSHH-II. The most common adverse events were grade 3-4 febrile neutropenia (48 patients [59%]), neutropenia (21 [26%]), infection with neutropenia (20 [25%]), leucopenia (15 [19%]), vomiting (15 [19%]), and anorexia (13 [16%]). No treatment-related deaths occurred. INTERPRETATION: The risk-adapted approach did not improve event-free survival in young children with medulloblastoma. However, the methylation subgroup analyses showed that the SHH subgroup had improved progression-free survival compared with the group 3 subgroup. Moreover, within the SHH subgroup, the iSHH-II subtype had improved progression-free survival in the absence of radiation, intraventricular chemotherapy, or high-dose chemotherapy compared with the iSHH-I subtype. These findings support the development of a molecularly driven, risk-adapted, treatment approach in future trials in young children with medulloblastoma. FUNDING: American Lebanese Syrian Associated Charities, St Jude Children's Research Hospital, NCI Cancer Center, Alexander and Margaret Stewart Trust, Sontag Foundation, and American Association for Cancer Research.

Spectrum and prevalence of genetic predisposition in medulloblastoma: a retrospective genetic study and prospective validation in a clinical trial cohort.

BACKGROUND: Medulloblastoma is associated with rare hereditary cancer predisposition syndromes; however, consensus medulloblastoma predisposition genes have not been defined and screening guidelines for genetic counselling and testing for paediatric patients are not available. We aimed to assess and define these genes to provide evidence for future screening guidelines. METHODS: In this international, multicentre study, we analysed patients with medulloblastoma from retrospective cohorts (International Cancer Genome Consortium [ICGC] PedBrain, Medulloblastoma Advanced Genomics International Consortium [MAGIC], and the CEFALO series) and from prospective cohorts from four clinical studies (SJMB03, SJMB12, SJYC07, and I-HIT-MED). Whole-genome sequences and exome sequences from blood and tumour samples were analysed for rare damaging germline mutations in cancer predisposition genes. DNA methylation profiling was done to determine consensus molecular subgroups: WNT (MBWNT), SHH (MBSHH), group 3 (MBGroup3), and group 4 (MBGroup4). Medulloblastoma predisposition genes were predicted on the basis of rare variant burden tests against controls without a cancer diagnosis from the Exome Aggregation Consortium (ExAC). Previously defined somatic mutational signatures were used to further classify medulloblastoma genomes into two groups, a clock-like group (signatures 1 and 5) and a homologous recombination repair deficiency-like group (signatures 3 and 8), and chromothripsis was investigated using previously established criteria. Progression-free survival and overall survival were modelled for patients with a genetic predisposition to medulloblastoma. FINDINGS: We included a total of 1022 patients with medulloblastoma from the retrospective cohorts (n=673) and the four prospective studies (n=349), from whom blood samples (n=1022) and tumour samples (n=800) were analysed for germline mutations in 110 cancer predisposition genes. In our rare variant burden analysis, we compared these against 53 105 sequenced controls from ExAC and identified APC, BRCA2, PALB2, PTCH1, SUFU, and TP53 as consensus medulloblastoma predisposition genes according to our rare variant burden analysis and estimated that germline mutations accounted for 6% of medulloblastoma diagnoses in the retrospective cohort. The prevalence of genetic predispositions differed between molecular subgroups in the retrospective cohort and was highest for patients in the MBSHH subgroup (20% in the retrospective cohort). These estimates were replicated in the prospective clinical cohort (germline mutations accounted for 5% of medulloblastoma diagnoses, with the highest prevalence [14%] in the MBSHH subgroup). Patients with germline APC mutations developed MBWNT and accounted for most (five [71%] of seven) cases of MBWNT that had no somatic CTNNB1 exon 3 mutations. Patients with germline mutations in SUFU and PTCH1 mostly developed infant MBSHH. Germline TP53 mutations presented only in childhood patients in the MBSHH subgroup and explained more than half (eight [57%] of 14) of all chromothripsis events in this subgroup. Germline mutations in PALB2 and BRCA2 were observed across the MBSHH, MBGroup3, and MBGroup4 molecular subgroups and were associated with mutational signatures typical of homologous recombination repair deficiency. In patients with a genetic predisposition to medulloblastoma, 5-year progression-free survival was 52% (95% CI 40-69) and 5-year overall survival was 65% (95% CI 52-81); these survival estimates differed significantly across patients with germline mutations in different medulloblastoma predisposition genes. INTERPRETATION: Genetic counselling and testing should be used as a standard-of-care procedure in patients with MBWNT and MBSHH because these patients have the highest prevalence of damaging germline mutations in known cancer predisposition genes. We propose criteria for routine genetic screening for patients with medulloblastoma based on clinical and molecular tumour characteristics. FUNDING: German Cancer Aid; German Federal Ministry of Education and Research; German Childhood Cancer Foundation (Deutsche Kinderkrebsstiftung); European Research Council; National Institutes of Health; Canadian Institutes for Health Research; German Cancer Research Center; St Jude Comprehensive Cancer Center; American Lebanese Syrian Associated Charities; Swiss National Science Foundation; European Molecular Biology Organization; Cancer Research UK; Hertie Foundation; Alexander and Margaret Stewart Trust; V Foundation for Cancer Research; Sontag Foundation; Musicians Against Childhood Cancer; BC Cancer Foundation; Swedish Council for Health, Working Life and Welfare; Swedish Research Council; Swedish Cancer Society; the Swedish Radiation Protection Authority; Danish Strategic Research Council; Swiss Federal Office of Public Health; Swiss Research Foundation on Mobile Communication; Masaryk University; Ministry of Health of the Czech Republic; Research Council of Norway; Genome Canada; Genome BC; Terry Fox Research Institute; Ontario Institute for Cancer Research; Pediatric Oncology Group of Ontario; The Family of Kathleen Lorette and the Clark H Smith Brain Tumour Centre; Montreal Children's Hospital Foundation; The Hospital for Sick Children: Sonia and Arthur Labatt Brain Tumour Research Centre, Chief of Research Fund, Cancer Genetics Program, Garron Family Cancer Centre, MDT's Garron Family Endowment; BC Childhood Cancer Parents Association; Cure Search Foundation; Pediatric Brain Tumor Foundation; Brainchild; and the Government of Ontario.

No correlation between NF1 mutation position and risk of optic pathway glioma in 77 unrelated NF1 patients.

Neurofibromatosis type 1 (NF1) is a common monogenic disorder whereby affected individuals are predisposed to developing CNS tumors, including optic pathway gliomas (OPGs, occurring in ~15 to 20 % of cases). So far, no definite genotype-phenotype correlation determining NF1 patients at risk for tumor formation has been described, although enrichment for mutations in the 5' region of the NF1 gene in OPG patients has been suggested. We used whole exome sequencing, targeted sequencing, and copy number analysis to screen 77 unrelated NF1 patients with (n = 41) or without (n = 36; age ≥10 years) optic pathway glioma for germline NF1 alterations. We identified germline NF1 mutations in 69 of 77 patients (90 %), but no genotype-phenotype correlation was observed. Our data using a larger patient cohort did not confirm the previously reported clustering of mutations in the 5' region of the NF1 gene in patients with OPG. Thus, NF1 mutation location should not currently be used as a clinical criterion to assess the risk of developing OPGs.

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.

Identification and removal of low-complexity sites in allele-specific analysis of ChIP-seq data.

MOTIVATION: High-throughput sequencing technologies enable the genome-wide analysis of the impact of genetic variation on molecular phenotypes at unprecedented resolution. However, although powerful, these technologies can also introduce unexpected artifacts. RESULTS: We investigated the impact of library amplification bias on the identification of allele-specific (AS) molecular events from high-throughput sequencing data derived from chromatin immunoprecipitation assays (ChIP-seq). Putative AS DNA binding activity for RNA polymerase II was determined using ChIP-seq data derived from lymphoblastoid cell lines of two parent-daughter trios. We found that, at high-sequencing depth, many significant AS binding sites suffered from an amplification bias, as evidenced by a larger number of clonal reads representing one of the two alleles. To alleviate this bias, we devised an amplification bias detection strategy, which filters out sites with low read complexity and sites featuring a significant excess of clonal reads. This method will be useful for AS analyses involving ChIP-seq and other functional sequencing assays. AVAILABILITY: The R package abs filter for library clonality simulations and detection of amplification-biased sites is available from http://updepla1srv1.epfl.ch/waszaks/absfilter

Genomic variation and its impact on gene expression in Drosophila melanogaster.

Understanding the relationship between genetic and phenotypic variation is one of the great outstanding challenges in biology. To meet this challenge, comprehensive genomic variation maps of human as well as of model organism populations are required. Here, we present a nucleotide resolution catalog of single-nucleotide, multi-nucleotide, and structural variants in 39 Drosophila melanogaster Genetic Reference Panel inbred lines. Using an integrative, local assembly-based approach for variant discovery, we identify more than 3.6 million distinct variants, among which were more than 800,000 unique insertions, deletions (indels), and complex variants (1 to 6,000 bp). While the SNP density is higher near other variants, we find that variants themselves are not mutagenic, nor are regions with high variant density particularly mutation-prone. Rather, our data suggest that the elevated SNP density around variants is mainly due to population-level processes. We also provide insights into the regulatory architecture of gene expression variation in adult flies by mapping cis-expression quantitative trait loci (cis-eQTLs) for more than 2,000 genes. Indels comprise around 10% of all cis-eQTLs and show larger effects than SNP cis-eQTLs. In addition, we identified two-fold more gene associations in males as compared to females and found that most cis-eQTLs are sex-specific, revealing a partial decoupling of the genomic architecture between the sexes as well as the importance of genetic factors in mediating sex-biased gene expression. Finally, we performed RNA-seq-based allelic expression imbalance analyses in the offspring of crosses between sequenced lines, which revealed that the majority of strong cis-eQTLs can be validated in heterozygous individuals.

Relating CNVs to transcriptome data at fine resolution: assessment of the effect of variant size, type, and overlap with functional regions.

Copy-number variants (CNVs) form an abundant class of genetic variation with a presumed widespread impact on individual traits. While recent advances, such as the population-scale sequencing of human genomes, facilitated the fine-scale mapping of CNVs, the phenotypic impact of most of these CNVs remains unclear. By relating copy-number genotypes to transcriptome sequencing data, we have evaluated the impact of CNVs, mapped at fine scale, on gene expression. Based on data from 129 individuals with ancestry from two populations, we identified CNVs associated with the expression of 110 genes, with 13% of the associations involving complex, multiallelic CNVs. Categorization of CNVs according to variant type, size, and gene overlap enabled us to examine the impact of different CNV classes on expression variation. While many small (<4 kb) CNVs were associated with expression variation, overall we observed an enrichment of large duplications and deletions, including large intergenic CNVs, relative to the entire set of expression-associated CNVs. Furthermore, the copy number of genes intersecting with CNVs typically correlated positively with the genes' expression, and also was more strongly correlated with expression than nearby single nucleotide polymorphisms, suggesting a frequent causal role of CNVs in expression quantitative trait loci (eQTLs). We also elucidated unexpected cases of negative correlations between copy number and expression by assessing the CNVs' effects on the structure and regulation of genes. Finally, we examined dosage compensation of transcript levels. Our results suggest that association studies can gain in resolution and power by including fine-scale CNV information, such as those obtained from population-scale sequencing.

A yeast one-hybrid and microfluidics-based pipeline to map mammalian gene regulatory networks.

The comprehensive mapping of gene promoters and enhancers has significantly improved our understanding of how the mammalian regulatory genome is organized. An important challenge is to elucidate how these regulatory elements contribute to gene expression by identifying their trans-regulatory inputs. Here, we present the generation of a mouse-specific transcription factor (TF) open-reading frame clone library and its implementation in yeast one-hybrid assays to enable large-scale protein-DNA interaction detection with mouse regulatory elements. Once specific interactions are identified, we then use a microfluidics-based method to validate and precisely map them within the respective DNA sequences. Using well-described regulatory elements as well as orphan enhancers, we show that this cross-platform pipeline characterizes known and uncovers many novel TF-DNA interactions. In addition, we provide evidence that several of these novel interactions are relevant in vivo and aid in elucidating the regulatory architecture of enhancers.