JAGuaR is an alignment protocol for RNA-seq reads that uses an extended reference to increase alignment sensitivity. It uses BWA to align reads to the genome and reference transcript models ...(including annotated exon-exon junctions) specifically allowing for the possibility of a single read spanning multiple exons. Reads aligned to the transcript models are then re-mapped on to genomic coordinates, transforming alignments that span multiple exons into large-gapped alignments on the genome. While JAGuaR does not detect novel junctions, we demonstrate how JAGuaR generates fast and accurate transcriptome alignments, which allows for both sensitive and specific SNV calling.
Colorectal cancer (CRC) displays a complex pattern of inheritance. It is postulated that much of the missing heritability of CRC is enshrined in high-impact rare alleles, which are mechanistically ...and clinically important. In this study, we assay the impact of rare germline mutations on CRC, analysing high-coverage exome sequencing data on 1,006 early-onset familial CRC cases and 1,609 healthy controls, with additional sequencing and array data on up to 5,552 cases and 6,792 controls. We identify highly penetrant rare mutations in 16% of familial CRC. Although the majority of these reside in known genes, we identify POT1, POLE2 and MRE11 as candidate CRC genes. We did not identify any coding low-frequency alleles (1-5%) with moderate effect. Our study clarifies the genetic architecture of CRC and probably discounts the existence of further major high-penetrance susceptibility genes, which individually account for >1% of the familial risk. Our results inform future study design and provide a resource for contextualizing the impact of new CRC genes.
The comprehensive multiplatform genomics data generated by The Cancer Genome Atlas (TCGA) Research Network is an enabling resource for cancer research. It includes an unprecedented amount of microRNA ...sequence data: ~11 000 libraries across 33 cancer types. Combined with initiatives like the National Cancer Institute Genomics Cloud Pilots, such data resources will make intensive analysis of large-scale cancer genomics data widely accessible. To support such initiatives, and to enable comparison of TCGA microRNA data to data from other projects, we describe the process that we developed and used to generate the microRNA sequence data, from library construction through to submission of data to repositories. In the context of this process, we describe the computational pipeline that we used to characterize microRNA expression across large patient cohorts.
Knowledge of the contribution of high-penetrance susceptibility to familial colorectal cancer (CRC) is relevant to the counseling, treatment, and surveillance of CRC patients and families.
To ...quantify the impact of germline mutation to familial CRC, we sequenced the mismatch repair genes (MMR) APC, MUTYH, and SMAD4/BMPR1A in 626 early-onset familial CRC cases ascertained through a population-based United Kingdom national registry. In addition, we evaluated the contribution of mutations in the exonuclease domain (exodom) of POLE and POLD1 genes that have recently been reported to confer CRC risk.
Overall mutations (pathogenic, likely pathogenic) in MMR genes make the highest contribution to familial CRC (10.9%). Mutations in the other established CRC genes account for 3.3% of cases. POLE/POLD1 exodom mutations were identified in three patients with family histories consistent with dominant transmission of CRC. Collectively, mutations in the known genes account for 14.2% of familial CRC (89 of 626 cases; 95% CI = 11.5, 17.2).
A genetic diagnosis is feasible in a high proportion of familial CRC. Mainstreaming such analysis in clinical practice should enable the medical management of patients and their families to be optimized. Findings suggest CRC screening of POLE and POLD1 mutation carriers should be comparable to that afforded to those at risk of HNPCC. Although the risk of CRC associated with unexplained familial CRC is in general moderate, in some families the risk is substantive and likely to be the consequence of unidentified genes, as exemplified by POLE and POLD1. Our findings have utility in the design of genetic analyses to identify such novel CRC risk genes.
Aerobic glycolysis supports proliferation through unresolved mechanisms. We have previously shown that aerobic glycolysis is required for the regulated proliferation of cerebellar granule neuron ...progenitors (CGNP) and for the growth of CGNP-derived medulloblastoma. Blocking the initiation of glycolysis via deletion of
(
) disrupts CGNP proliferation and restricts medulloblastoma growth. Here, we assessed whether disrupting
(
), an enzyme that acts in the terminal steps of glycolysis, would alter CGNP metabolism, proliferation, and tumorigenesis. We observed a dichotomous pattern of PKM expression, in which postmitotic neurons throughout the brain expressed the constitutively active PKM1 isoform, while neural progenitors and medulloblastomas exclusively expressed the less active PKM2. Isoform-specific
deletion in CGNPs blocked all
expression.
-deleted CGNPs showed reduced lactate production and increased SHH-driven proliferation.
C-flux analysis showed that
deletion reduced the flow of glucose carbons into lactate and glutamate without markedly increasing glucose-to-ribose flux.
deletion accelerated tumor formation in medulloblastoma-prone
mice, indicating the disrupting PKM releases CGNPs from a tumor-suppressive effect. These findings show that distal and proximal disruptions of glycolysis have opposite effects on proliferation, and that efforts to block the oncogenic effect of aerobic glycolysis must target reactions upstream of PKM.
.
Acute lymphoblastic leukemia (ALL) is the major pediatric cancer diagnosed in economically developed countries with B-cell precursor (BCP)-ALL, accounting for approximately 70% of ALL. Recent ...genome-wide association studies (GWAS) have provided the first unambiguous evidence for common inherited susceptibility to BCP-ALL, identifying susceptibility loci at 7p12.2, 9p21.3, 10q21.2, and 14q11.2. To identify additional BCP-ALL susceptibility loci, we conducted a GWAS and performed a meta-analysis with a published GWAS totaling 1658 cases and 4723 controls, with validation in 1449 cases and 1488 controls. Combined analysis identified novel loci mapping to 10p12.2 (rs10828317, odds ratio OR = 1.23; P = 2.30 × 10−9) and 10p14 marked by rs3824662 (OR = 1.31; P = 8.62 × 10−12). The single nucleotide polymorphism rs10828317 is responsible for the N215S polymorphism in exon 7 of PIP4K2A, and rs3824662 localizes to intron 3 of the transcription factor and putative tumor suppressor gene GATA3. The rs10828317 association was shown to be specifically associated with hyperdiploid ALL, whereas the rs3824662-associated risk was confined to nonhyperdiploid non–TEL-AML1 + ALL. The risk allele of rs3824662 was correlated with older age at diagnosis (P < .001) and significantly worse event-free survivorship (P < .0001). These findings provide further insights into the genetic and biological basis of inherited genetic susceptibility to BCP-ALL and the influence of constitutional genotype on disease development.
Key Points
Malignant rhabdoid tumors (MRTs) are rare lethal tumors of childhood that most commonly occur in the kidney and brain. MRTs are driven by SMARCB1 loss, but the molecular consequences of SMARCB1 loss ...in extra-cranial tumors have not been comprehensively described and genomic resources for analyses of extra-cranial MRT are limited. To provide such data, we used whole-genome sequencing, whole-genome bisulfite sequencing, whole transcriptome (RNA-seq) and microRNA sequencing (miRNA-seq), and histone modification profiling to characterize extra-cranial MRTs. Our analyses revealed gene expression and methylation subgroups and focused on dysregulated pathways, including those involved in neural crest development.
•Extra-cranial malignant rhabdoid tumors (MRTs) exhibit molecular heterogeneity•Evidence is presented for epigenetic reprogramming of HOX genes•MRTs exhibit dysregulated expression of genes involved in neural crest development•Dysregulation of oncogenes and tumor suppressor genes is reported
Chun et al. perform integrated molecular analyses of extra-cranial malignant rhabdoid tumors (MRTs) and show that, although SMARCB1 loss drives nearly all MRTs, there are two subgroups of MRTs that are associated with patient age and differentially expressed genes.
To identify risk variants for multiple myeloma, we conducted a genome-wide association study of 1,675 individuals with multiple myeloma and 5,903 control subjects. We identified risk loci for ...multiple myeloma at 3p22.1 (rs1052501 in ULK4; odds ratio (OR) = 1.32; P = 7.47 × 10(-9)) and 7p15.3 (rs4487645, OR = 1.38; P = 3.33 × 10(-15)). In addition, we observed a promising association at 2p23.3 (rs6746082, OR = 1.29; P = 1.22 × 10(-7)). Our study identifies new genomic regions associated with multiple myeloma risk that may lead to new etiological insights.
We investigated the role of 3D genome architecture in instructing functional properties of glioblastoma stem cells (GSCs) by generating sub-5-kb resolution 3D genome maps by in situ Hi-C. Contact ...maps at sub-5-kb resolution allow identification of individual DNA loops, domain organization, and large-scale genome compartmentalization. We observed differences in looping architectures among GSCs from different patients, suggesting that 3D genome architecture is a further layer of inter-patient heterogeneity for glioblastoma. Integration of DNA contact maps with chromatin and transcriptional profiles identified specific mechanisms of gene regulation, including the convergence of multiple super enhancers to individual stemness genes within individual cells. We show that the number of loops contacting a gene correlates with elevated transcription. These results indicate that stemness genes are hubs of interaction between multiple regulatory regions, likely to ensure their sustained expression. Regions of open chromatin common among the GSCs tested were poised for expression of immune-related genes, including
We demonstrate that this gene is co-expressed with stemness genes in GSCs and that CD276 can be targeted with an antibody-drug conjugate to eliminate self-renewing cells. Our results demonstrate that integrated structural genomics data sets can be employed to rationally identify therapeutic vulnerabilities in self-renewing cells.
Abstract
Summary
Reliably identifying genomic rearrangements and interpreting their impact is a key step in understanding their role in human cancers and inherited genetic diseases. Many short read ...algorithmic approaches exist but all have appreciable false negative rates. A common approach is to evaluate the union of multiple tools increasing sensitivity, followed by filtering to retain specificity. Here we describe an application framework for the rapid generation of structural variant consensus, unique in its ability to visualize the genetic impact and context as well as process both genome and transcriptome data.
Availability and implementation
http://mavis.bcgsc.ca
Supplementary information
Supplementary data are available at Bioinformatics online.