Defining the condensate landscape of fusion oncoproteins.

Fusion oncoproteins (FOs) arise from chromosomal translocations in ~17% of cancers and are often oncogenic drivers. Although some FOs can promote oncogenesis by undergoing liquid-liquid phase separation (LLPS) to form aberrant biomolecular condensates, the generality of this phenomenon is unknown. We explored this question by testing 166 FOs in HeLa cells and found that 58% formed condensates. The condensate-forming FOs displayed physicochemical features distinct from those of condensate-negative FOs and segregated into distinct feature-based groups that aligned with their sub-cellular localization and biological function. Using Machine Learning, we developed a predictor of FO condensation behavior, and discovered that 67% of ~3000 additional FOs likely form condensates, with 35% of those predicted to function by altering gene expression. 47% of the predicted condensate-negative FOs were associated with cell signaling functions, suggesting a functional dichotomy between condensate-positive and -negative FOs. Our Datasets and reagents are rich resources to interrogate FO condensation in the future.

The genomic landscape of pediatric acute lymphoblastic leukemia.

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Here, using whole-genome, exome and transcriptome sequencing of 2,754 childhood patients with ALL, we find that, despite a generally low mutation burden, ALL cases harbor a median of four putative somatic driver alterations per sample, with 376 putative driver genes identified varying in prevalence across ALL subtypes. Most samples harbor at least one rare gene alteration, including 70 putative cancer driver genes associated with ubiquitination, SUMOylation, noncoding transcripts and other functions. In hyperdiploid B-ALL, chromosomal gains are acquired early and synchronously before ultraviolet-induced mutation. By contrast, ultraviolet-induced mutations precede chromosomal gains in B-ALL cases with intrachromosomal amplification of chromosome 21. We also demonstrate the prognostic significance of genetic alterations within subtypes. Intriguingly, DUX4- and KMT2A-rearranged subtypes separate into CEBPA/FLT3- or NFATC4-expressing subgroups with potential clinical implications. Together, these results deepen understanding of the ALL genomic landscape and associated outcomes.

SequencErr: measuring and suppressing sequencer errors in next-generation sequencing data.

BACKGROUND: There is currently no method to precisely measure the errors that occur in the sequencing instrument/sequencer, which is critical for next-generation sequencing applications aimed at discovering the genetic makeup of heterogeneous cellular populations. RESULTS: We propose a novel computational method, SequencErr, to address this challenge by measuring the base correspondence between overlapping regions in forward and reverse reads. An analysis of 3777 public datasets from 75 research institutions in 18 countries revealed the sequencer error rate to be ~ 10 per million (pm) and 1.4% of sequencers and 2.7% of flow cells have error rates > 100 pm. At the flow cell level, error rates are elevated in the bottom surfaces and > 90% of HiSeq and NovaSeq flow cells have at least one outlier error-prone tile. By sequencing a common DNA library on different sequencers, we demonstrate that sequencers with high error rates have reduced overall sequencing accuracy, and removal of outlier error-prone tiles improves sequencing accuracy. We demonstrate that SequencErr can reveal novel insights relative to the popular quality control method FastQC and achieve a 10-fold lower error rate than popular error correction methods including Lighter and Musket. CONCLUSIONS: Our study reveals novel insights into the nature of DNA sequencing errors incurred on DNA sequencers. Our method can be used to assess, calibrate, and monitor sequencer accuracy, and to computationally suppress sequencer errors in existing datasets.

RNAIndel: discovering somatic coding indels from tumor RNA-Seq data.

MOTIVATION: Reliable identification of expressed somatic insertions/deletions (indels) is an unmet need due to artifacts generated in PCR-based RNA-Seq library preparation and the lack of normal RNA-Seq data, presenting analytical challenges for discovery of somatic indels in tumor transcriptome. RESULTS: We present RNAIndel, a tool for predicting somatic, germline and artifact indels from tumor RNA-Seq data. RNAIndel leverages features derived from indel sequence context and biological effect in a machine-learning framework. Except for tumor samples with microsatellite instability, RNAIndel robustly predicts 88-100% of somatic indels in five diverse test datasets of pediatric and adult cancers, even recovering subclonal (VAF range 0.01-0.15) driver indels missed by targeted deep-sequencing, outperforming the current best-practice for RNA-Seq variant calling which had 57% sensitivity but with 14 times more false positives. AVAILABILITY AND IMPLEMENTATION: RNAIndel is freely available at https://github.com/stjude/RNAIndel. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Molecular Mechanism of Telomere Length Dynamics and Its Prognostic Value in Pediatric Cancers.

BACKGROUND: We aimed to systematically evaluate telomere dynamics across a spectrum of pediatric cancers, search for underlying molecular mechanisms, and assess potential prognostic value. METHODS: The fraction of telomeric reads was determined from whole-genome sequencing data for paired tumor and normal samples from 653 patients with 23 cancer types from the Pediatric Cancer Genome Project. Telomere dynamics were characterized as the ratio of telomere fractions between tumor and normal samples. Somatic mutations were gathered, RNA sequencing data for 330 patients were analyzed for gene expression, and Cox regression was used to assess the telomere dynamics on patient survival. RESULTS: Telomere lengthening was observed in 28.7% of solid tumors, 10.5% of brain tumors, and 4.3% of hematological cancers. Among 81 samples with telomere lengthening, 26 had somatic mutations in alpha thalassemia/mental retardation syndrome X-linked gene, corroborated by a low level of the gene expression in the subset of tumors with RNA sequencing. Telomerase reverse transcriptase gene amplification and/or activation was observed in 10 tumors with telomere lengthening, including two leukemias of the E2A-PBX1 subtype. Among hematological cancers, pathway analysis for genes with expressions most negatively correlated with telomere fractions suggests the implication of a gene ontology process of antigen presentation by Major histocompatibility complex class II. A higher ratio of telomere fractions was statistically significantly associated with poorer survival for patients with brain tumors (hazard ratio = 2.18, 95% confidence interval = 1.37 to 3.46). CONCLUSION: Because telomerase inhibitors are currently being explored as potential agents to treat pediatric cancer, these data are valuable because they identify a subpopulation of patients with reactivation of telomerase who are most likely to benefit from this novel therapeutic option.

Pediatric Cancer Variant Pathogenicity Information Exchange (PeCanPIE): a cloud-based platform for curating and classifying germline variants.

Variant interpretation in the era of massively parallel sequencing is challenging. Although many resources and guidelines are available to assist with this task, few integrated end-to-end tools exist. Here, we present the Pediatric Cancer Variant Pathogenicity Information Exchange (PeCanPIE), a web- and cloud-based platform for annotation, identification, and classification of variations in known or putative disease genes. Starting from a set of variants in variant call format (VCF), variants are annotated, ranked by putative pathogenicity, and presented for formal classification using a decision-support interface based on published guidelines from the American College of Medical Genetics and Genomics (ACMG). The system can accept files containing millions of variants and handle single-nucleotide variants (SNVs), simple insertions/deletions (indels), multiple-nucleotide variants (MNVs), and complex substitutions. PeCanPIE has been applied to classify variant pathogenicity in cancer predisposition genes in two large-scale investigations involving >4000 pediatric cancer patients and serves as a repository for the expert-reviewed results. PeCanPIE was originally developed for pediatric cancer but can be easily extended for use for nonpediatric cancers and noncancer genetic diseases. Although PeCanPIE's web-based interface was designed to be accessible to non-bioinformaticians, its back-end pipelines may also be run independently on the cloud, facilitating direct integration and broader adoption. PeCanPIE is publicly available and free for research use.

VCF2CNA: A tool for efficiently detecting copy-number alterations in VCF genotype data and tumor purity.

VCF2CNA is a tool (Linux commandline or web-interface) for copy-number alteration (CNA) analysis and tumor purity estimation of paired tumor-normal VCF variant file formats. It operates on whole genome and whole exome datasets. To benchmark its performance, we applied it to 46 adult glioblastoma and 146 pediatric neuroblastoma samples sequenced by Illumina and Complete Genomics (CGI) platforms respectively. VCF2CNA was highly consistent with a state-of-the-art algorithm using raw sequencing data (mean F1-score = 0.994) in high-quality whole genome glioblastoma samples and was robust to uneven coverage introduced by library artifacts. In the whole genome neuroblastoma set, VCF2CNA identified MYCN high-level amplifications in 31 of 32 clinically validated samples compared to 15 found by CGI's HMM-based CNA model. Moreover, VCF2CNA achieved highly consistent CNA profiles between WGS and WXS platforms (mean F1 score 0.97 on a set of 15 rhabdomyosarcoma samples). In addition, VCF2CNA provides accurate tumor purity estimates for samples with sufficient CNAs. These results suggest that VCF2CNA is an accurate, efficient and platform-independent tool for CNA and tumor purity analyses without accessing raw sequence data.

Clinical genome sequencing uncovers potentially targetable truncations and fusions of MAP3K8 in spitzoid and other melanomas.

Spitzoid melanoma is a specific morphologic variant of melanoma that most commonly affects children and adolescents, and ranges on the spectrum of malignancy from low grade to overtly malignant. These tumors are generally driven by fusions of ALK, RET, NTRK1/3, MET, ROS1 and BRAF1,2. However, in approximately 50% of cases no genetic driver has been established2. Clinical whole-genome and transcriptome sequencing (RNA-Seq) of a spitzoid tumor from an adolescent revealed a novel gene fusion of MAP3K8, encoding a serine-threonine kinase that activates MEK3,4. The patient, who had exhausted all other therapeutic options, was treated with a MEK inhibitor and underwent a transient clinical response. We subsequently analyzed spitzoid tumors from 49 patients by RNA-Seq and found in-frame fusions or C-terminal truncations of MAP3K8 in 33% of cases. The fusion transcripts and truncated genes all contained MAP3K8 exons 1-8 but lacked the autoinhibitory final exon. Data mining of RNA-Seq from the Cancer Genome Atlas (TCGA) uncovered analogous MAP3K8 rearrangements in 1.5% of adult melanomas. Thus, MAP3K8 rearrangements-uncovered by comprehensive clinical sequencing of a single case-are the most common genetic event in spitzoid melanoma, are present in adult melanomas and could be amenable to MEK inhibition.

Polygenic Determinants for Subsequent Breast Cancer Risk in Survivors of Childhood Cancer: The St Jude Lifetime Cohort Study (SJLIFE).

PURPOSE: The risk of subsequent breast cancer among female childhood cancer survivors is markedly elevated. We aimed to determine genetic contributions to this risk, focusing on polygenic determinants implicated in breast cancer susceptibility in the general population. EXPERIMENTAL DESIGN: Whole-genome sequencing (30×) was performed on survivors in the St Jude Lifetime Cohort, and germline mutations in breast cancer predisposition genes were classified for pathogenicity. A polygenic risk score (PRS) was constructed for each survivor using 170 established common risk variants. Relative rate (RR) and 95% confidence interval (95% CI) of subsequent breast cancer incidence were estimated using multivariable piecewise exponential regression. RESULTS: The analysis included 1,133 female survivors of European ancestry (median age at last follow-up = 35.4 years; range, 8.4-67.4), of whom 47 were diagnosed with one or more subsequent breast cancers (median age at subsequent breast cancer = 40.3 years; range, 24.5-53.0). Adjusting for attained age, age at primary diagnosis, chest irradiation, doses of alkylating agents and anthracyclines, and genotype eigenvectors, RRs for survivors with PRS in the highest versus lowest quintiles were 2.7 (95% CI, 1.0-7.3), 3.0 (95% CI, 1.1-8.1), and 2.4 (95% CI, 0.1-81.1) for all survivors and survivors with and without chest irradiation, respectively. Similar associations were observed after excluding carriers of pathogenic/likely pathogenic mutations in breast cancer predisposition genes. Notably, the PRS was associated with the subsequent breast cancer rate under the age of 45 years (RR = 3.2; 95% CI, 1.2-8.3). CONCLUSIONS: Genetic profiles comprised of small-effect common variants and large-effect predisposing mutations can inform personalized breast cancer risk and surveillance/intervention in female childhood cancer survivors.

Genetic Risk for Subsequent Neoplasms Among Long-Term Survivors of Childhood Cancer.

Purpose Childhood cancer survivors are at increased risk of subsequent neoplasms (SNs), but the germline genetic contribution is largely unknown. We assessed the contribution of pathogenic/likely pathogenic (P/LP) mutations in cancer predisposition genes to their SN risk. Patients and Methods Whole-genome sequencing (30-fold) was performed on samples from childhood cancer survivors who were ≥ 5 years since initial cancer diagnosis and participants in the St Jude Lifetime Cohort Study, a retrospective hospital-based study with prospective clinical follow-up. Germline mutations in 60 genes known to be associated with autosomal dominant cancer predisposition syndromes with moderate to high penetrance were classified by their pathogenicity according to the American College of Medical Genetics and Genomics guidelines. Relative rates (RRs) and 95% CIs of SN occurrence by mutation status were estimated using multivariable piecewise exponential regression stratified by radiation exposure. Results Participants were 3,006 survivors (53% male; median age, 35.8 years [range, 7.1 to 69.8 years]; 56% received radiotherapy), 1,120 SNs were diagnosed among 439 survivors (14.6%), and 175 P/LP mutations were identified in 5.8% (95% CI, 5.0% to 6.7%) of survivors. Mutations were associated with significantly increased rates of breast cancer (RR, 13.9; 95% CI, 6.0 to 32.2) and sarcoma (RR, 10.6; 95% CI, 4.3 to 26.3) among irradiated survivors and with increased rates of developing any SN (RR, 4.7; 95% CI, 2.4 to 9.3), breast cancer (RR, 7.7; 95% CI, 2.4 to 24.4), nonmelanoma skin cancer (RR, 11.0; 95% CI, 2.9 to 41.4), and two or more histologically distinct SNs (RR, 18.6; 95% CI, 3.5 to 99.3) among nonirradiated survivors. Conclusion The findings support referral of all survivors for genetic counseling for potential clinical genetic testing, which should be prioritized for nonirradiated survivors with any SN and for those with breast cancer or sarcoma in the field of prior irradiation.