visPIG--a web tool for producing multi-region, multi-track, multi-scale plots of genetic data

A genome-wide association study identifies susceptibility loci for primary central nervous system lymphoma at 6p25.3 and 3p22.1: a LOC Network study

BACKGROUND

Primary central nervous system lymphoma (PCNSL) is a rare form of extra-nodal non-Hodgkin lymphoma. PCNSL is a distinct subtype of non-Hodgkin lymphoma, with over 95% of tumors belonging to the diffuse large B-cell lymphoma (DLBCL) group. We have conducted a genome-wide association study (GWAS) on immunocompetent patients to address the possibility that common genetic variants influence the risk of developing PCNSL.

METHODS

We performed a meta-analysis of 2 new GWASs of PCNSL totaling 475 cases and 1134 controls of European ancestry. To increase genomic resolution, we imputed >10 million single nucleotide polymorphisms using the 1000 Genomes Project combined with UK10K as reference. In addition we performed a transcription factor binding disruption analysis and investigated the patterns of local chromatin by Capture Hi-C data.

RESULTS

We identified independent risk loci at 3p22.1 (rs41289586, ANO10, P = 2.17 × 10-8) and 6p25.3 near EXOC2 (rs116446171, P = 1.95 x 10-13). In contrast, the lack of an association between rs41289586 and DLBCL suggests distinct germline predisposition to PCNSL and DLBCL. We found looping chromatin interactions between noncoding regions at 6p25.3 (rs11646171) with the IRF4 promoter and at 8q24.21 (rs13254990) with the MYC promoter, both genes with strong relevance to B-cell tumorigenesis.

CONCLUSION

To our knowledge this is the first study providing insight into the genetic predisposition to PCNSL. Our findings represent an important step in defining the contribution of common genetic variation to the risk of developing PCNSL.

Large-Scale Profiling of Signaling Pathways Reveals a Distinct Demarcation between Normal and Extended Liver Resection

Despite numerous studies addressing normal liver regeneration, we still lack comprehensive understanding of the biological processes underlying failed liver regeneration. Therefore, we analyzed the activity of 271 intracellular signaling pathways (ISPs) by genome wide profiling of differentially expressed RNAs in murine liver tissue biopsies after normal hepatectomy (nHx; 68% of liver removed) and extended hepatectomy (eHx; 86% of liver removed). Comprehensive, genome-wide transcriptome profiling using RNAseq was performed in liver tissue obtained from mice (sham, nHx, and eHx) harvested 1, 8, 16, 32, and 48 h after operation (n = 3 per group) and the OncoFinder toolkit was used for an unsupervised, unbiased identification of intracellular signaling pathways (ISP) activity. We observed that the normal regenerative process requires a transient activation and silencing of approximately two dozen of ISPs. After nHx, the Akt Pathway represented with 13 branches, the Chromatin Pathway and the DDR Pathways dominated. After eHx, the ATM main pathway and two of its branches (Cell Survival; G2_M Checkpoint Arrest) dominated, as well as the Hypoxia Pathways. Further, 14 ISPs demonstrated a strong inverse regulation, with the Hedgehog and the Brca1 Main Pathways as chief activators after nHx, and the ATM Pathway(G2_M Checkpoint Arrest) as the dominating constraining response after eHx.

Exploration of the Transcriptional Landscape of ALPPS Reveals the Pathways of Accelerated Liver Regeneration

Background and Aims: ALPPS (associating liver partition and portal vein ligation for staged hepatectomy), a novel 2-staged hepatectomy, dramatically accelerates liver regeneration and thus enables extensive liver tumor resection. The signaling networks underlying the ALPPS-induced accelerated regeneration process are largely unknown.

Methods: We performed transcriptome profiling (TP) of liver tissue obtained from a mouse model of ALPPS, standard hepatectomy (68% model), and additional control surgeries (sham, PVL and Tx). We also performed TP using human liver biopsies (n = 5) taken from the occluded lobe and the future liver remnant (FLR) during the first step of ALPPS surgery (4–5 h apart). We used Oncofinder computational tools, which covers 378 ISPs, for unsupervised, unbiased quantification of ISP activity.

Results: Gene expression cluster analysis revealed an ALPPS specific signature: the IGF1R Signaling Pathway (Cell survival), the ILK Pathway (Induced cell proliferation), and the IL-10 Pathway (Stability determination) were significantly enriched, whereas the activity of the Interferon Pathway (Transcription) was reduced (p < 0.05). Further, the PAK- and ILK-associated ISPs were activated at an earlier time point, reflecting significant acceleration of liver regeneration (p < 0.001). These pathways, which were also recovered in human liver biopsies, control cell growth and proliferation, inflammatory response, and hypoxia-related processes.

Conclusions: ALPPS is not a straightforward addition of portal vein ligation (PVL) plus transection—it is more. The early stages of normal and accelerated liver regeneration are clearly discernible by a significantly increased and earlier activation of a small number of signaling pathways. Compounds mimicking these responses may help to improve the ALPPS method and further reduce the hospitalization time of the patient.

Genome-wide association analysis identifies a meningioma risk locus at 11p15.5

Background

Meningiomas are adult brain tumors originating in the meningeal coverings of the brain and spinal cord, with significant heritable basis. Genome-wide association studies (GWAS) have previously identified only a single risk locus for meningioma, at 10p12.31.

Methods

To identify a susceptibility locus for meningioma, we conducted a meta-analysis of 2 GWAS, imputed using a merged reference panel from the 1000 Genomes Project and UK10K data, with validation in 2 independent sample series totaling 2138 cases and 12081 controls.

Results

We identified a new susceptibility locus for meningioma at 11p15.5 (rs2686876, odds ratio = 1.44, P = 9.86 × 10^–9). A number of genes localize to the region of linkage disequilibrium encompassing rs2686876, including RIC8A, which plays a central role in the development of neural crest-derived structures, such as the meninges.

Conclusions

This finding advances our understanding of the genetic basis of meningioma development and provides additional support for a polygenic model of meningioma.

In Vitro Digestion and Fermentation of Three Polysaccharide Fractions from Laminaria japonica and Their Impact on Lipid Metabolism-Associated Human Gut Microbiota

Our previous study has proved that the three polysaccharide fractions from L. japonica (LP-A4, LP-A6, and LP-A8) had significantly different structure characterization. Herein, we conducted in vitro simulated digestion and fermentation to study the digestive mechanism of LP-As. The results of gastrointestinal digestion indicated that LP-A6 and LP-A8 would be easier to trap the enzyme molecules for their denser interconnected macromolecule network compared with LP-A4. Fermentation of LP-As by human gut microbiota, especially for LP-A8, generated a large amount of short-chain fatty acids (SCFAs), which could upregulate the abundance of Firmicutes ( Lachnoclostridium and Eubacterium). The high content of sulfate and highly branched sugar residue of LP-A8 might help it be easily used by Firmicutes in gut microbiota of hyperlipidemic patients. Functional analysis revealed that the increased metabolic activities of glycerophospholipid metabolism, ether lipid metabolism, and fatty acid metabolism induced by LP-A8 treatment were closely associated with metabolic syndromes and hyperlipidemia.

Association analyses identify 31 new risk loci for colorectal cancer susceptibility

Colorectal cancer (CRC) is a leading cause of cancer-related death worldwide, and has a strong heritable basis. We report a genome-wide association analysis of 34,627 CRC cases and 71,379 controls of European ancestry that identifies SNPs at 31 new CRC risk loci. We also identify eight independent risk SNPs at the new and previously reported European CRC loci, and a further nine CRC SNPs at loci previously only identified in Asian populations. We use in situ promoter capture Hi-C (CHi-C), gene expression, and in silico annotation methods to identify likely target genes of CRC SNPs. Whilst these new SNP associations implicate target genes that are enriched for known CRC pathways such as Wnt and BMP, they also highlight novel pathways with no prior links to colorectal tumourigenesis. These findings provide further insight into CRC susceptibility and enhance the prospects of applying genetic risk scores to personalised screening and prevention.

CBX6 is negatively regulated by EZH2 and plays a potential tumor suppressor role in breast cancer

Chromobox 6 (CBX6) is a subunit of Polycomb Repressive Complex 1 (PRC1) that mediates epigenetic gene repression and acts as an oncogene or tumor suppressor in a cancer type-dependent manner. The specific function of CBX6 in breast cancer is currently undefined. In this study, a comprehensive analysis of The Cancer Genome Atlas (TCGA) dataset led to the identification of CBX6 as a consistently downregulated gene in breast cancer. We provided evidence showing enhancer of zeste homolog 2 (EZH2) negatively regulated CBX6 expression in a Polycomb Repressive Complex 2 (PRC2)-dependent manner. Exogenous overexpression of CBX6 inhibited cell proliferation and colony formation, and induced cell cycle arrest along with suppression of migration and invasion of breast cancer cells in vitro. Microarray analyses revealed that CBX6 governs a complex gene expression program. Moreover, CBX6 induced significant downregulation of bone marrow stromal cell antigen-2 (BST2), a potential therapeutic target, via interactions with its promoter region. Our collective findings support a tumor suppressor role of CBX6 in breast cancer.

Genome-wide association study implicates immune dysfunction in the development of Hodgkin lymphoma

To further our understanding of inherited susceptibility to Hodgkin lymphoma (HL), we performed a meta-analysis of 7 genome-wide association studies totaling 5325 HL cases and 22 423 control patients. We identify 5 new HL risk loci at 6p21.31 (rs649775; P = 2.11 × 10-10), 6q23.3 (rs1002658; P = 2.97 × 10-8), 11q23.1 (rs7111520; P = 1.44 × 10-11), 16p11.2 (rs6565176; P = 4.00 × 10-8), and 20q13.12 (rs2425752; P = 2.01 × 10-8). Integration of gene expression, histone modification, and in situ promoter capture Hi-C data at the 5 new and 13 known risk loci implicates dysfunction of the germinal center reaction, disrupted T-cell differentiation and function, and constitutive NF-κB activation as mechanisms of predisposition. These data provide further insights into the genetic susceptibility and biology of HL.

Omics approaches to study gene regulatory networks for development in echinoderms

Gene regulatory networks (GRNs) describe the interactions for a developmental process at a given time and space. Historically, perturbation experiments represent one of the key methods for analyzing and reconstructing a GRN, and the GRN governing early development in the sea urchin embryo stands as one of the more deeply dissected so far. As technology progresses, so do the methods used to address different biological questions. Next-generation sequencing (NGS) has become a standard experimental technique for genome and transcriptome sequencing and studies of protein-DNA interactions and DNA accessibility. While several efforts have been made toward the integration of different omics approaches for the study of the regulatory genome in many animals, in a few cases, these are applied with the purpose of reconstructing and experimentally testing developmental GRNs. Here, we review emerging approaches integrating multiple NGS technologies for the prediction and validation of gene interactions within echinoderm GRNs. These approaches can be applied to both 'model' and 'non-model' organisms. Although a number of issues still need to be addressed, advances in NGS applications, such as assay for transposase-accessible chromatin sequencing, combined with the availability of embryos belonging to different species, all separated by various evolutionary distances and accessible to experimental regulatory biology, place echinoderms in an unprecedented position for the reconstruction and evolutionary comparison of developmental GRNs. We conclude that sequencing technologies and integrated omics approaches allow the examination of GRNs on a genome-wide scale only if biological perturbation and cis-regulatory analyses are experimentally accessible, as in the case of echinoderm embryos.

Genome-wide association study identifies susceptibility loci for B-cell childhood acute lymphoblastic leukemia

Genome-wide association studies (GWAS) have advanced our understanding of susceptibility to B-cell precursor acute lymphoblastic leukemia (BCP-ALL); however, much of the heritable risk remains unidentified. Here, we perform a GWAS and conduct a meta-analysis with two existing GWAS, totaling 2442 cases and 14,609 controls. We identify risk loci for BCP-ALL at 8q24.21 (rs28665337, P = 3.86 × 10-9, odds ratio (OR) = 1.34) and for ETV6-RUNX1 fusion-positive BCP-ALL at 2q22.3 (rs17481869, P = 3.20 × 10-8, OR = 2.14). Our findings provide further insights into genetic susceptibility to ALL and its biology.

Validation of loci at 2q14.2 and 15q21.3 as risk factors for testicular cancer

Testicular germ cell tumor (TGCT), the most common cancer in men aged 18 to 45 years, has a strong heritable basis. Genome-wide association studies (GWAS) have proposed single nucleotide polymorphisms (SNPs) at a number of loci influencing TGCT risk. To further evaluate the association of recently proposed risk SNPs with TGCT at 2q14.2, 3q26.2, 7q36.3, 10q26.13 and 15q21.3, we analyzed genotype data on 3,206 cases and 7,422 controls. Our analysis provides independent replication of the associations for risk SNPs at 2q14.2 (rs2713206 at P = 3.03 × 10-2; P-meta = 3.92 × 10-8; nearest gene, TFCP2L1) and rs12912292 at 15q21.3 (P = 7.96 × 10-11; P-meta = 1.55 × 10-19; nearest gene PRTG). Case-only analyses did not reveal specific associations with TGCT histology. TFCP2L1 joins the growing list of genes located within TGCT risk loci with biologically plausible roles in developmental transcriptional regulation, further highlighting the importance of this phenomenon in TGCT oncogenesis.

Large-scale profiling of signalling pathways reveals an asthma specific signature in bronchial smooth muscle cells

Background

Bronchial smooth muscle (BSM) cells from asthmatic patients maintain in vitro a distinct hyper-reactive (“primed”) phenotype, characterized by increased release of pro-inflammatory factors and mediators, as well as hyperplasia and/or hypertrophy. This “primed” phenotype helps to understand pathogenesis of asthma, as changes in BSM function are essential for manifestation of allergic and inflammatory responses and airway wall remodelling.

Objective

To identify signalling pathways in cultured primary BSMs of asthma patients and non-asthmatic subjects by genome wide profiling of differentially expressed mRNAs and activated intracellular signalling pathways (ISPs).

Methods

Transcriptome profiling by cap-analysis-of-gene-expression (CAGE), which permits selection of preferentially capped mRNAs most likely to be translated into proteins, was performed in human BSM cells from asthmatic (n=8) and non-asthmatic (n=6) subjects and OncoFinder tool were then exploited for identification of ISP deregulations.

Results

CAGE revealed >600 RNAs differentially expressed in asthma vs control cells (p≤0.005), with asthma samples showing a high degree of similarity among them. Comprehensive ISP activation analysis revealed that among 269 pathways analysed, 145 (p<0.05) or 103 (p<0.01) are differentially active in asthma, with profiles that clearly characterize BSM cells of asthmatic individuals. Notably, we identified 7 clusters of coherently acting pathways functionally related to the disease, with ISPs down-regulated in asthma mostly targeting cell death-promoting pathways and up-regulated ones affecting cell growth and proliferation, inflammatory response, control of smooth muscle contraction and hypoxia-related signalization.

Conclusions

These first-time results can now be exploited toward development of novel therapeutic strategies targeting ISP signatures linked to asthma pathophysiology.

Genome-wide association study of classical Hodgkin lymphoma identifies key regulators of disease susceptibility

Several susceptibility loci for classical Hodgkin lymphoma have been reported. However, much of the heritable risk is unknown. Here, we perform a meta-analysis of two existing genome-wide association studies, a new genome-wide association study, and replication totalling 5,314 cases and 16,749 controls. We identify risk loci for all classical Hodgkin lymphoma at 6q22.33 (rs9482849, P = 1.52 × 10-8) and for nodular sclerosis Hodgkin lymphoma at 3q28 (rs4459895, P = 9.43 × 10-17), 6q23.3 (rs6928977, P = 4.62 × 10-11), 10p14 (rs3781093, P = 9.49 × 10-13), 13q34 (rs112998813, P = 4.58 × 10-8) and 16p13.13 (rs34972832, P = 2.12 × 10-8). Additionally, independent loci within the HLA region are observed for nodular sclerosis Hodgkin lymphoma (rs9269081, HLA-DPB1*03:01, Val86 in HLA-DRB1) and mixed cellularity Hodgkin lymphoma (rs1633096, rs13196329, Val86 in HLA-DRB1). The new and established risk loci localise to areas of active chromatin and show an over-representation of transcription factor binding for determinants of B-cell development and immune response.

C-State: an interactive web app for simultaneous multi-gene visualization and comparative epigenetic pattern search

Background

Comparative epigenomic analysis across multiple genes presents a bottleneck for bench biologists working with NGS data. Despite the development of standardized peak analysis algorithms, the identification of novel epigenetic patterns and their visualization across gene subsets remains a challenge.

Results

We developed a fast and interactive web app, C-State (Chromatin-State), to query and plot chromatin landscapes across multiple loci and cell types. C-State has an interactive, JavaScript-based graphical user interface and runs locally in modern web browsers that are pre-installed on all computers, thus eliminating the need for cumbersome data transfer, pre-processing and prior programming knowledge.

Conclusions

C-State is unique in its ability to extract and analyze multi-gene epigenetic information. It allows for powerful GUI-based pattern searching and visualization. We include a case study to demonstrate its potential for identifying user-defined epigenetic trends in context of gene expression profiles.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1786-6) contains supplementary material, which is available to authorized users.

Challenges for visualizing three-dimensional data in genomic browsers

Genomic interactions reveal the spatial organization of genomes and genomic domains, which is known to play key roles in cell function. Physical proximity can be represented as two-dimensional heat maps or matrices. From these, three-dimensional (3D) conformations of chromatin can be computed revealing coherent structures that highlight the importance of nonsequential relationships across genomic features. Mainstream genomic browsers have been classically developed to display compact, stacked tracks based on a linear, sequential, per-chromosome coordinate system. Genome-wide comparative analysis demands new approaches to data access and new layouts for analysis. The legibility can be compromised when displaying track-aligned second dimension matrices, which require greater screen space. Moreover, 3D representations of genomes defy vertical alignment in track-based genome browsers. Furthermore, investigation at previously unattainable levels of detail is revealing multiscale, multistate, time-dependent complexity. This article outlines how these challenges are currently handled in mainstream browsers as well as how novel techniques in visualization are being explored to address them. A set of requirements for coherent visualization of novel spatial genomic data is defined and the resulting potential for whole genome visualization is described.

Data aggregation at the level of molecular pathways improves stability of experimental transcriptomic and proteomic data

High throughput technologies opened a new era in biomedicine by enabling massive analysis of gene expression at both RNA and protein levels. Unfortunately, expression data obtained in different experiments are often poorly compatible, even for the same biologic samples. Here, using experimental and bioinformatic investigation of major experimental platforms, we show that aggregation of gene expression data at the level of molecular pathways helps to diminish cross- and intra-platform bias otherwise clearly seen at the level of individual genes. We created a mathematical model of cumulative suppression of data variation that predicts the ideal parameters and the optimal size of a molecular pathway. We compared the abilities to aggregate experimental molecular data for the 5 alternative methods, also evaluated by their capacity to retain meaningful features of biologic samples. The bioinformatic method OncoFinder showed optimal performance in both tests and should be very useful for future cross-platform data analyses.

Bioinformatic analysis of microRNA-mRNA expression profiles of bladder tissue induced by bladder outlet obstruction in a rat model

Various microRNAs (miRNAs) have previously been demonstrated to exhibit an association with the process of bladder remodeling, induced by bladder outlet obstruction (BOO). However, little is known about miRNA and gene expression profiles and the molecular mechanism underlying bladder pathophysiological alterations. The present study used bioinformatic analysis technology to examine the altered miRNA and mRNA expression profiles of bladder tissue in a rat model of BOO and validate the involved signaling pathways. The gene expression profile data was downloaded from Gene Expression Omnibus (GEO), and the differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) were screened. Potential target genes of DEMs were predicted. The target genes and DEGs were used for further gene ontology (GO) analysis followed by Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using the Database for Annotation, Visualization and Integrated Discovery. The present study additionally constructed a DEM‑DEG interaction network. A total of 9 DEMs (3 upregulated and 6 downregulated) were identified; 664 DEGs were screened. KEGG analysis revealed that the DEGs were involved in the regulation of the actin cytoske-leton, extracellular matrix (ECM) remodeling, cell adhesion and the cell cycle. Additionally, KEGG classification indicated that these genes were important in angiogenesis, and in the p53 and transforming growth factor‑β signaling pathways. Notably, rno‑miRNA (miR)‑26b and rno‑miR‑101b were the two larger nodes of the 7 obstruction‑associated DEMs and interacted with 32 and 27 DEGs, respectively. On removal of obstruction, few DEMs were present; however, 370 genes exhibited the opposite expression trend. The majority of pathways enriched for the DEGs were identified and were associated with ECM‑receptor interaction and focal adhesion. In the DEM‑DEG regulatory network, miR‑495, miR‑494 and their target genes were significantly differentially expressed. The present study demonstrated that miRNAs and genes may be potential biomarkers of bladder remodeling induced by BOO, and additionally provided novel insights into the molecular mechanisms underlying this disorder.

Gene expression and molecular pathway activation signatures of MYCN-amplified neuroblastomas

Neuroblastoma is a pediatric cancer arising from sympathetic nervous system. Remarkable heterogeneity in outcomes is one of its widely known features. One of the traits strongly associated with the unfavorable subtype is the amplification of oncogene MYCN. Here, we performed cross-platform biomarker detection by comparing gene expression and pathway activation patterns from the two literature reports and from our experimental dataset, combining profiles for the 761 neuroblastoma patients with known MYCN amplification status. We identified 109 / 25 gene expression / pathway activation biomarkers strongly linked with the MYCN amplification. The marker genes/pathways are involved in the processes of purine nucleotide biosynthesis, ATP-binding, tetrahydrofolate metabolism, building mitochondrial matrix, biosynthesis of amino acids, tRNA aminoacylation and NADP-linked oxidation-reduction processes, as well as in the tyrosine phosphatase activity, p53 signaling, cell cycle progression and the G1/S and G2/M checkpoints. To connect molecular functions of the genes involved in MYCN-amplified phenotype, we built a new molecular pathway using known intracellular protein interaction networks. The activation of this pathway was highly selective in discriminating MYCN-amplified neuroblastomas in all three datasets. Our data also suggest that the phosphoinositide 3-kinase (PI3K) inhibitors may provide new opportunities for the treatment of the MYCN-amplified neuroblastoma subtype.

Identification of 19 new risk loci and potential regulatory mechanisms influencing susceptibility to testicular germ cell tumor

Genome-wide association studies (GWAS) have transformed understanding of susceptibility to testicular germ cell tumors (TGCTs), but much of the heritability remains unexplained. Here we report a new GWAS, a meta-analysis with previous GWAS and a replication series, totaling 7,319 TGCT cases and 23,082 controls. We identify 19 new TGCT risk loci, roughly doubling the number of known TGCT risk loci to 44. By performing in situ Hi-C in TGCT cells, we provide evidence for a network of physical interactions among all 44 TGCT risk SNPs and candidate causal genes. Our findings implicate widespread disruption of developmental transcriptional regulators as a basis of TGCT susceptibility, consistent with failed primordial germ cell differentiation as an initiating step in oncogenesis. Defective microtubule assembly and dysregulation of KIT-MAPK signaling also feature as recurrently disrupted pathways. Our findings support a polygenic model of risk and provide insight into the biological basis of TGCT.

Genome-wide association study of glioma subtypes identifies specific differences in genetic susceptibility to glioblastoma and non-glioblastoma tumors

Genome-wide association studies (GWAS) have transformed our understanding of glioma susceptibility, but individual studies have had limited power to identify risk loci. We performed a meta-analysis of existing GWAS and two new GWAS, which totaled 12,496 cases and 18,190 controls. We identified five new loci for glioblastoma (GBM) at 1p31.3 (rs12752552; P = 2.04 × 10-9, odds ratio (OR) = 1.22), 11q14.1 (rs11233250; P = 9.95 × 10-10, OR = 1.24), 16p13.3 (rs2562152; P = 1.93 × 10-8, OR = 1.21), 16q12.1 (rs10852606; P = 1.29 × 10-11, OR = 1.18) and 22q13.1 (rs2235573; P = 1.76 × 10-10, OR = 1.15), as well as eight loci for non-GBM tumors at 1q32.1 (rs4252707; P = 3.34 × 10-9, OR = 1.19), 1q44 (rs12076373; P = 2.63 × 10-10, OR = 1.23), 2q33.3 (rs7572263; P = 2.18 × 10-10, OR = 1.20), 3p14.1 (rs11706832; P = 7.66 × 10-9, OR = 1.15), 10q24.33 (rs11598018; P = 3.39 × 10-8, OR = 1.14), 11q21 (rs7107785; P = 3.87 × 10-10, OR = 1.16), 14q12 (rs10131032; P = 5.07 × 10-11, OR = 1.33) and 16p13.3 (rs3751667; P = 2.61 × 10-9, OR = 1.18). These data substantiate that genetic susceptibility to GBM and non-GBM tumors are highly distinct, which likely reflects different etiology.

Genetic and regulatory mechanism of susceptibility to high-hyperdiploid acute lymphoblastic leukaemia at 10p21.2

Despite high-hyperdiploid acute lymphoblastic leukaemia (HD-ALL) being the most common subgroup of paediatric ALL, its aetiology remains unknown. Genome-wide association studies have demonstrated association at 10q21.2. Here, we sought to determine how this region influences HD-ALL risk. We impute genotypes across the locus, finding the single nucleotide polymorphism rs7090445 highly associated with HD-ALL (P=1.54 × 10⁻³⁸), and residing in a predicted enhancer element. We show this region physically interacts with the transcription start site of ARID5B, that alleles of rs7090445 have differential enhancer activity and influence RUNX3 binding. RUNX3 knock-down reduces ARID5B expression and rs7090445 enhancer activity. Individuals carrying the rs7090445-C risk allele also have reduced ARID5B expression. Finally, the rs7090445-C risk allele is preferentially retained in HD-ALL blasts consistent with inherited genetic variation contributing to arrest of normal lymphocyte development, facilitating leukaemic clonal expansion. These data provide evidence for a biological mechanism underlying hereditary risk of HD-ALL at 10q21.2.

Risk for the paediatric cancer high-hyperdiploid acute lymphoblastic leukaemia (HD-ALL) has been associated with genetic variants at 10q21.2. Here, the authors characterize this region, establishing a single risk variant and showing its role in dysregulated expression of ARID5B.

A genome-wide association study identifies risk loci for childhood acute lymphoblastic leukemia at 10q26.13 and 12q23.1

Genome-wide association studies (GWASs) have shown that common genetic variation contributes to the heritable risk of childhood acute lymphoblastic leukemia (ALL). To identify new susceptibility loci for the largest subtype of ALL, B-cell precursor ALL (BCP-ALL), we conducted a meta-analysis of two GWASs with imputation using 1000 Genomes and UK10K Project data as reference (totaling 1658 cases and 7224 controls). After genotyping an additional 2525 cases and 3575 controls, we identify new susceptibility loci for BCP-ALL mapping to 10q26.13 (rs35837782, LHPP, P=1.38 × 10-11) and 12q23.1 (rs4762284, ELK3, P=8.41 × 10-9). We also provide confirmatory evidence for the existence of independent risk loci at 9p21.3, but show that the association marked by rs77728904 can be accounted for by linkage disequilibrium with the rare high-impact CDKN2A p.Ala148Thr variant rs3731249. Our data provide further insights into genetic susceptibility to ALL and its biology.

Genome-wide association analysis implicates dysregulation of immunity genes in chronic lymphocytic leukaemia

Several chronic lymphocytic leukaemia (CLL) susceptibility loci have been reported; however, much of the heritable risk remains unidentified. Here we perform a meta-analysis of six genome-wide association studies, imputed using a merged reference panel of 1,000 Genomes and UK10K data, totalling 6,200 cases and 17,598 controls after replication. We identify nine risk loci at 1p36.11 (rs34676223, P=5.04 × 10-13), 1q42.13 (rs41271473, P=1.06 × 10-10), 4q24 (rs71597109, P=1.37 × 10-10), 4q35.1 (rs57214277, P=3.69 × 10-8), 6p21.31 (rs3800461, P=1.97 × 10-8), 11q23.2 (rs61904987, P=2.64 × 10-11), 18q21.1 (rs1036935, P=3.27 × 10-8), 19p13.3 (rs7254272, P=4.67 × 10-8) and 22q13.33 (rs140522, P=2.70 × 10-9). These new and established risk loci map to areas of active chromatin and show an over-representation of transcription factor binding for the key determinants of B-cell development and immune response.

Genome-wide association analysis of chronic lymphocytic leukaemia, Hodgkin lymphoma and multiple myeloma identifies pleiotropic risk loci

B-cell malignancies (BCM) originate from the same cell of origin, but at different maturation stages and have distinct clinical phenotypes. Although genetic risk variants for individual BCMs have been identified, an agnostic, genome-wide search for shared genetic susceptibility has not been performed. We explored genome-wide association studies of chronic lymphocytic leukaemia (CLL, N = 1,842), Hodgkin lymphoma (HL, N = 1,465) and multiple myeloma (MM, N = 3,790). We identified a novel pleiotropic risk locus at 3q22.2 (NCK1, rs11715604, P = 1.60 × 10-9) with opposing effects between CLL (P = 1.97 × 10-8) and HL (P = 3.31 × 10-3). Eight established non-HLA risk loci showed pleiotropic associations. Within the HLA region, Ser37 + Phe37 in HLA-DRB1 (P = 1.84 × 10-12) was associated with increased CLL and HL risk (P = 4.68 × 10-12), and reduced MM risk (P = 1.12 × 10-2), and Gly70 in HLA-DQB1 (P = 3.15 × 10-10) showed opposing effects between CLL (P = 3.52 × 10-3) and HL (P = 3.41 × 10-9). By integrating eQTL, Hi-C and ChIP-seq data, we show that the pleiotropic risk loci are enriched for B-cell regulatory elements, as well as an over-representation of binding of key B-cell transcription factors. These data identify shared biological pathways influencing the development of CLL, HL and MM. The identification of these risk loci furthers our understanding of the aetiological basis of BCMs.

Molecular pathway activation features of pediatric acute myeloid leukemia (AML) and acute lymphoblast leukemia (ALL) cells

Acute lymphoblast leukemia (ALL) is characterized by overproduction of immature white blood cells in the bone marrow. ALL is most common in the childhood and has high (>80%) cure rate. In contrast, acute myeloid leukemia (AML) has far greater mortality rate than the ALL and is most commonly affecting older adults. However, AML is a leading cause of childhood cancer mortality. In this study, we compare gene expression and molecular pathway activation patterns in three normal blood, seven pediatric ALL and seven pediatric AML bone marrow samples. We identified 172/94 and 148/31 characteristic gene expression/pathway activation signatures, clearly distinguishing pediatric ALL and AML cells, respectively, from the normal blood. The pediatric AML and ALL cells differed by 139/34 gene expression/pathway activation biomarkers. For the adult 30 AML and 17 normal blood samples, we found 132/33 gene expression/pathway AML-specific features, of which only 7/2 were common for the adult and pediatric AML and, therefore, age-independent. At the pathway level, we found more differences than similarities between the adult and pediatric forms. These findings suggest that the adult and pediatric AMLs may require different treatment strategies.

The epichaperome is an integrated chaperome network that facilitates tumour survival

Transient, multi-protein complexes are important facilitators of cellular functions. This includes the chaperome, an abundant protein family comprising chaperones, co-chaperones, adaptors, and folding enzymes-dynamic complexes of which regulate cellular homeostasis together with the protein degradation machinery. Numerous studies have addressed the role of chaperome members in isolation, yet little is known about their relationships regarding how they interact and function together in malignancy. As function is probably highly dependent on endogenous conditions found in native tumours, chaperomes have resisted investigation, mainly due to the limitations of methods needed to disrupt or engineer the cellular environment to facilitate analysis. Such limitations have led to a bottleneck in our understanding of chaperome-related disease biology and in the development of chaperome-targeted cancer treatment. Here we examined the chaperome complexes in a large set of tumour specimens. The methods used maintained the endogenous native state of tumours and we exploited this to investigate the molecular characteristics and composition of the chaperome in cancer, the molecular factors that drive chaperome networks to crosstalk in tumours, the distinguishing factors of the chaperome in tumours sensitive to pharmacologic inhibition, and the characteristics of tumours that may benefit from chaperome therapy. We find that under conditions of stress, such as malignant transformation fuelled by MYC, the chaperome becomes biochemically 'rewired' to form a network of stable, survival-facilitating, high-molecular-weight complexes. The chaperones heat shock protein 90 (HSP90) and heat shock cognate protein 70 (HSC70) are nucleating sites for these physically and functionally integrated complexes. The results indicate that these tightly integrated chaperome units, here termed the epichaperome, can function as a network to enhance cellular survival, irrespective of tissue of origin or genetic background. The epichaperome, present in over half of all cancers tested, has implications for diagnostics and also provides potential vulnerability as a target for drug intervention.

Genetic Predisposition to Chronic Lymphocytic Leukemia Is Mediated by a BMF Super-Enhancer Polymorphism

Chronic lymphocytic leukemia (CLL) is an adult B cell malignancy. Genome-wide association studies show that variation at 15q15.1 influences CLL risk. We deciphered the causal variant at 15q15.1 and the mechanism by which it influences tumorigenesis. We imputed all possible genotypes across the locus and then mapped highly associated SNPs to areas of chromatin accessibility, evolutionary conservation, and transcription factor binding. SNP rs539846 C>A, the most highly associated variant (p = 1.42 × 10(-13), odds ratio = 1.35), localizes to a super-enhancer defined by extensive histone H3 lysine 27 acetylation in intron 3 of B cell lymphoma 2 (BCL2)-modifying factor (BMF). The rs539846-A risk allele alters a conserved RELA-binding motif, disrupts RELA binding, and is associated with decreased BMF expression in CLL. These findings are consistent with rs539846 influencing CLL susceptibility through differential RELA binding, with direct modulation of BMF expression impacting on anti-apoptotic BCL2, a hallmark of oncogenic dependency in CLL.

Genome-wide association study identifies multiple susceptibility loci for multiple myeloma

Multiple myeloma (MM) is a plasma cell malignancy with a significant heritable basis. Genome-wide association studies have transformed our understanding of MM predisposition, but individual studies have had limited power to discover risk loci. Here we perform a meta-analysis of these GWAS, add a new GWAS and perform replication analyses resulting in 9,866 cases and 239,188 controls. We confirm all nine known risk loci and discover eight new loci at 6p22.3 (rs34229995, P=1.31 × 10(-8)), 6q21 (rs9372120, P=9.09 × 10(-15)), 7q36.1 (rs7781265, P=9.71 × 10(-9)), 8q24.21 (rs1948915, P=4.20 × 10(-11)), 9p21.3 (rs2811710, P=1.72 × 10(-13)), 10p12.1 (rs2790457, P=1.77 × 10(-8)), 16q23.1 (rs7193541, P=5.00 × 10(-12)) and 20q13.13 (rs6066835, P=1.36 × 10(-13)), which localize in or near to JARID2, ATG5, SMARCD3, CCAT1, CDKN2A, WAC, RFWD3 and PREX1. These findings provide additional support for a polygenic model of MM and insight into the biological basis of tumour development.

Variation at 2q35 (PNKD and TMBIM1) influences colorectal cancer risk and identifies a pleiotropic effect with inflammatory bowel disease

To identify new risk loci for colorectal cancer (CRC), we conducted a meta-analysis of seven genome-wide association studies (GWAS) with independent replication, totalling 13 656 CRC cases and 21 667 controls of European ancestry. The combined analysis identified a new risk association for CRC at 2q35 marked by rs992157 (P = 3.15 × 10-8, odds ratio = 1.10, 95% confidence interval = 1.06-1.13), which is intronic to PNKD (paroxysmal non-kinesigenic dyskinesia) and TMBIM1 (transmembrane BAX inhibitor motif containing 1). Intriguingly this susceptibility single-nucleotide polymorphism (SNP) is in strong linkage disequilibrium (r2 = 0.90, D' = 0.96) with the previously discovered GWAS SNP rs2382817 for inflammatory bowel disease (IBD). Following on from this observation we examined for pleiotropy, or shared genetic susceptibility, between CRC and the 200 established IBD risk loci, identifying an additional 11 significant associations (false discovery rate [FDR]) < 0.05). Our findings provide further insight into the biological basis of inherited genetic susceptibility to CRC, and identify risk factors that may influence the development of both CRC and IBD.

GenomicInteractions: An R/Bioconductor package for manipulating and investigating chromatin interaction data

Background

Precise quantitative and spatiotemporal control of gene expression is necessary to ensure proper cellular differentiation and the maintenance of homeostasis. The relationship between gene expression and the spatial organisation of chromatin is highly complex, interdependent and not completely understood. The development of experimental techniques to interrogate both the higher-order structure of chromatin and the interactions between regulatory elements has recently lead to important insights on how gene expression is controlled. The ability to gain these and future insights is critically dependent on computational tools for the analysis and visualisation of data produced by these techniques.

Results and conclusion

We have developed GenomicInteractions, a freely available R/Bioconductor package designed for processing, analysis and visualisation of data generated from various types of chromosome conformation capture experiments. The package allows the easy annotation and summarisation of large genome-wide datasets at both the level of individual interactions and sets of genomic features, and provides several different methods for interrogating and visualising this type of data. We demonstrate this package’s utility by showing example analyses performed on interaction datasets generated using Hi-C and ChIA-PET.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2140-x) contains supplementary material, which is available to authorized users.

Identification of four new susceptibility loci for testicular germ cell tumour

Genome-wide association studies (GWAS) have identified multiple risk loci for testicular germ cell tumour (TGCT), revealing a polygenic model of disease susceptibility strongly influenced by common variation. To identify additional single-nucleotide polymorphisms (SNPs) associated with TGCT, we conducted a multistage GWAS with a combined data set of >25,000 individuals (6,059 cases and 19,094 controls). We identified new risk loci for TGCT at 3q23 (rs11705932, TFDP2, P=1.5 × 10(-9)), 11q14.1 (rs7107174, GAB2, P=9.7 × 10(-11)), 16p13.13 (rs4561483, GSPT1, P=1.6 × 10(-8)) and 16q24.2 (rs55637647, ZFPM1, P=3.4 × 10(-9)). We additionally present detailed functional analysis of these loci, identifying a statistically significant relationship between rs4561483 risk genotype and increased GSPT1 expression in TGCT patient samples. These findings provide additional support for a polygenic model of TGCT risk and further insight into the biological basis of disease development.

The 9p21.3 risk of childhood acute lymphoblastic leukaemia is explained by a rare high-impact variant in CDKN2A

Genome-wide association studies (GWAS) have provided strong evidence for inherited predisposition to childhood acute lymphoblastic leukaemia (ALL) identifying a number of risk loci. We have previously shown common SNPs at 9p21.3 influence ALL risk. These SNP associations are generally not themselves candidates for causality, but simply act as markers for functional variants. By means of imputation of GWAS data and subsequent validation SNP genotyping totalling 2,177 ALL cases and 8,240 controls, we have shown that the 9p21.3 association can be ascribed to the rare high-impact CDKN2A p.Ala148Thr variant (rs3731249; Odds ratio = 2.42, P = 3.45 × 10(-19)). The association between rs3731249 genotype and risk was not specific to particular subtype of B-cell ALL. The rs3731249 variant is associated with predominant nuclear localisation of the CDKN2A transcript suggesting the functional effect of p.Ala148Thr on ALL risk may be through compromised ability to inhibit cyclin D within the cytoplasm.

Genome-wide association study identifies multiple susceptibility loci for glioma

Previous genome-wide association studies (GWASs) have shown that common genetic variation contributes to the heritable risk of glioma. To identify new glioma susceptibility loci, we conducted a meta-analysis of four GWAS (totalling 4,147 cases and 7,435 controls), with imputation using 1000 Genomes and UK10K Project data as reference. After genotyping an additional 1,490 cases and 1,723 controls we identify new risk loci for glioblastoma (GBM) at 12q23.33 (rs3851634, near POLR3B, P=3.02 × 10(-9)) and non-GBM at 10q25.2 (rs11196067, near VTI1A, P=4.32 × 10(-8)), 11q23.2 (rs648044, near ZBTB16, P=6.26 × 10(-11)), 12q21.2 (rs12230172, P=7.53 × 10(-11)) and 15q24.2 (rs1801591, near ETFA, P=5.71 × 10(-9)). Our findings provide further insights into the genetic basis of the different glioma subtypes.

Common variation at 1q24.1 (ALDH9A1) is a potential risk factor for renal cancer

So far six susceptibility loci for renal cell carcinoma (RCC) have been discovered by genome-wide association studies (GWAS). To identify additional RCC common risk loci, we performed a meta-analysis of published GWAS (totalling 2,215 cases and 8,566 controls of Western-European background) with imputation using 1000 Genomes Project and UK10K Project data as reference panels and followed up the most significant association signals [22 single nucleotide polymorphisms (SNPs) and 3 indels in eight genomic regions] in 383 cases and 2,189 controls from The Cancer Genome Atlas (TCGA). A combined analysis identified a promising susceptibility locus mapping to 1q24.1 marked by the imputed SNP rs3845536 (Pcombined =2.30x10-8). Specifically, the signal maps to intron 4 of the ALDH9A1 gene (aldehyde dehydrogenase 9 family, member A1). We further evaluated this potential signal in 2,461 cases and 5,081 controls from the International Agency for Research on Cancer (IARC) GWAS of RCC cases and controls from multiple European regions. In contrast to earlier findings no association was shown in the IARC series (P=0.94; Pcombined =2.73x10-5). While variation at 1q24.1 represents a potential risk locus for RCC, future replication analyses are required to substantiate our observation.

Capture Hi-C identifies the chromatin interactome of colorectal cancer risk loci

Multiple regulatory elements distant from their targets on the linear genome can influence the expression of a single gene through chromatin looping. Chromosome conformation capture implemented in Hi-C allows for genome-wide agnostic characterization of chromatin contacts. However, detection of functional enhancer-promoter interactions is precluded by its effective resolution that is determined by both restriction fragmentation and sensitivity of the experiment. Here we develop a capture Hi-C (cHi-C) approach to allow an agnostic characterization of these physical interactions on a genome-wide scale. Single-nucleotide polymorphisms associated with complex diseases often reside within regulatory elements and exert effects through long-range regulation of gene expression. Applying this cHi-C approach to 14 colorectal cancer risk loci allows us to identify key long-range chromatin interactions in cis and trans involving these loci.