Copy number analysis of whole-genome data using BIC-seq2 and its application to detection of cancer susceptibility variants

Detecting copy-number alterations from single-cell chromatin sequencing data by AtaCNA

Single-cell assay of transposase-accessible chromatin sequencing (scATAC-seq) unbiasedly profiles genome-wide chromatin accessibility in single cells. In single-cell tumor studies, identification of normal cells or tumor clonal structures often relies on copy-number alterations (CNAs). However, CNA detection from scATAC-seq is difficult due to the high noise, sparsity, and confounding factors. Here, we describe AtaCNA, a computational algorithm that accurately detects high-resolution CNAs from scATAC-seq data. We benchmark AtaCNA using simulation and real data and find AtaCNA's superior performance. Analyses of 10 scATAC-seq datasets show that AtaCNA could effectively distinguish malignant from non-malignant cells. In glioblastoma, endometrial, and ovarian cancer samples, AtaCNA identifies subclones at distinct cellular states, suggesting an important interplay between genetic and epigenetic plasticity. Some tumor subclones only differ in small-scale (10-20 Mb) CNAs, demonstrating the importance of high-resolution CNA detection. These data show that AtaCNA can aid in integrative analysis to understand the complex heterogeneity in cancer.

Prediction of the 3D cancer genome from whole-genome sequencing using InfoHiC

The 3D genome prediction in cancer is crucial for uncovering the impact of structural variations (SVs) on tumorigenesis, especially when they are present in noncoding regions. We present InfoHiC, a systemic framework for predicting the 3D cancer genome directly from whole-genome sequencing (WGS). InfoHiC utilizes contig-specific copy number encoding on the SV contig assembly, and performs a contig-to-total Hi-C conversion for the cancer Hi-C prediction from multiple SV contigs. We showed that InfoHiC can predict 3D genome folding from all types of SVs using breast cancer cell line data. We applied it to WGS data of patients with breast cancer and pediatric patients with medulloblastoma, and identified neo topologically associating domains. For breast cancer, we discovered super-enhancer hijacking events associated with oncogenic overexpression and poor survival outcomes. For medulloblastoma, we found SVs in noncoding regions that caused super-enhancer hijacking events of medulloblastoma driver genes (GFI1, GFI1B, and PRDM6). In addition, we provide trained models for cancer Hi-C prediction from WGS at https://github.com/dmcb-gist/InfoHiC , uncovering the impacts of SVs in cancer patients and revealing novel therapeutic targets.

Somatic mutation phasing and haplotype extension using linked-reads in multiple myeloma

Somatic mutation phasing informs our understanding of cancer-related events, like driver mutations. We generated linked-read whole genome sequencing data for 23 samples across disease stages from 14 multiple myeloma (MM) patients and systematically assigned somatic mutations to haplotypes using linked-reads. Here, we report the reconstructed cancer haplotypes and phase blocks from several MM samples and show how phase block length can be extended by integrating samples from the same individual. We also uncover phasing information in genes frequently mutated in MM, including DIS3, HIST1H1E, KRAS, NRAS, and TP53, phasing 79.4% of 20,705 high-confidence somatic mutations. In some cases, this enabled us to interpret clonal evolution models at higher resolution using pairs of phased somatic mutations. For example, our analysis of one patient suggested that two NRAS hotspot mutations occurred on the same haplotype but were independent events in different subclones. Given sufficient tumor purity and data quality, our framework illustrates how haplotype-aware analysis of somatic mutations in cancer can be beneficial for some cancer cases.

Genome-wide analysis and visualization of copy number with CNVpytor in igv.js

SUMMARY

Copy number variation (CNV) and alteration (CNA) analysis is a crucial component in many genomic studies and its applications span from basic research to clinic diagnostics and personalized medicine. CNVpytor is a tool featuring a read depth-based caller and combined read depth and B-allele frequency (BAF) based 2D caller to find CNVs and CNAs. The tool stores processed intermediate data and CNV/CNA calls in a compact HDF5 file-pytor file. Here, we describe a new track in igv.js that utilizes pytor and whole genome variant files as input for on-the-fly read depth and BAF visualization, CNV/CNA calling and analysis. Embedding into HTML pages and Jupiter Notebooks enables convenient remote data access and visualization simplifying interpretation and analysis of omics data.

AVAILABILITY AND IMPLEMENTATION

The CNVpytor track is integrated with igv.js and available at https://github.com/igvteam/igv.js. The documentation is available at https://github.com/igvteam/igv.js/wiki/cnvpytor. Usage can be tested in the IGV-Web app at https://igv.org/app and also on https://github.com/abyzovlab/CNVpytor.

Multi-scale signaling and tumor evolution in high-grade gliomas

Although genomic anomalies in glioblastoma (GBM) have been well studied for over a decade, its 5-year survival rate remains lower than 5%. We seek to expand the molecular landscape of high-grade glioma, composed of IDH-wildtype GBM and IDH-mutant grade 4 astrocytoma, by integrating proteomic, metabolomic, lipidomic, and post-translational modifications (PTMs) with genomic and transcriptomic measurements to uncover multi-scale regulatory interactions governing tumor development and evolution. Applying 14 proteogenomic and metabolomic platforms to 228 tumors (212 GBM and 16 grade 4 IDH-mutant astrocytoma), including 28 at recurrence, plus 18 normal brain samples and 14 brain metastases as comparators, reveals heterogeneous upstream alterations converging on common downstream events at the proteomic and metabolomic levels and changes in protein-protein interactions and glycosylation site occupancy at recurrence. Recurrent genetic alterations and phosphorylation events on PTPN11 map to important regulatory domains in three dimensions, suggesting a central role for PTPN11 signaling across high-grade gliomas.

Comprehensive proteogenomic characterization of rare kidney tumors

Non-clear cell renal cell carcinomas (non-ccRCCs) encompass diverse malignant and benign tumors. Refinement of differential diagnosis biomarkers, markers for early prognosis of aggressive disease, and therapeutic targets to complement immunotherapy are current clinical needs. Multi-omics analyses of 48 non-ccRCCs compared with 103 ccRCCs reveal proteogenomic, phosphorylation, glycosylation, and metabolic aberrations in RCC subtypes. RCCs with high genome instability display overexpression of IGF2BP3 and PYCR1. Integration of single-cell and bulk transcriptome data predicts diverse cell-of-origin and clarifies RCC subtype-specific proteogenomic signatures. Expression of biomarkers MAPRE3, ADGRF5, and GPNMB differentiates renal oncocytoma from chromophobe RCC, and PIGR and SOSTDC1 distinguish papillary RCC from MTSCC. This study expands our knowledge of proteogenomic signatures, biomarkers, and potential therapeutic targets in non-ccRCC.

HATCHet2: clone- and haplotype-specific copy number inference from bulk tumor sequencing data

Bulk DNA sequencing of multiple samples from the same tumor is becoming common, yet most methods to infer copy-number aberrations (CNAs) from this data analyze individual samples independently. We introduce HATCHet2, an algorithm to identify haplotype- and clone-specific CNAs simultaneously from multiple bulk samples. HATCHet2 extends the earlier HATCHet method by improving identification of focal CNAs and introducing a novel statistic, the minor haplotype B-allele frequency (mhBAF), that enables identification of mirrored-subclonal CNAs. We demonstrate HATCHet2's improved accuracy using simulations and a single-cell sequencing dataset. HATCHet2 analysis of 10 prostate cancer patients reveals previously unreported mirrored-subclonal CNAs affecting cancer genes.

Circulating Nucleic Acids in Colorectal Cancer: Diagnostic and Prognostic Value

Colorectal cancer (CRC) is the third most prevalent cancer in the world and the fourth leading cause of cancer-related mortality. DNA (cfDNA/ctDNA) and RNA (cfRNA/ctRNA) in the blood are promising noninvasive biomarkers for molecular profiling, screening, diagnosis, treatment management, and prognosis of CRC. Technological advancements that enable precise detection of both genetic and epigenetic abnormalities, even in minute quantities in circulation, can overcome some of these challenges. This review focuses on testing for circulating nucleic acids in the circulation as a noninvasive method for CRC detection, monitoring, detection of minimal residual disease, and patient management. In addition, the benefits and drawbacks of various diagnostic techniques and associated bioinformatics tools have been detailed.

Integrated drug response prediction models pinpoint repurposed drugs with effectiveness against rhabdomyosarcoma

Targeted therapies for inhibiting the growth of cancer cells or inducing apoptosis are urgently needed for effective rhabdomyosarcoma (RMS) treatment. However, identifying cancer-targeting compounds with few side effects, among the many potential compounds, is expensive and time-consuming. A computational approach to reduce the number of potential candidate drugs can facilitate the discovery of attractive lead compounds. To address this and obtain reliable predictions of novel cell-line-specific drugs, we apply prediction models that have the potential to improve drug discovery approaches for RMS treatment. The results of two prediction models were ensemble and validated via in vitro experiments. The computational models were trained using data extracted from the Genomics of Drug Sensitivity in Cancer database and tested on two RMS cell lines to select potential RMS drug candidates. Among 235 candidate drugs, 22 were selected following the result of the computational approach, and three candidate drugs were identified (NSC207895, vorinostat, and belinostat) that showed selective effectiveness in RMS cell lines in vitro via the induction of apoptosis. Our in vitro experiments have demonstrated that our proposed methods can effectively identify and repurpose drugs for treating RMS.

Tracing cancer evolution and heterogeneity using Hi-C

Chromosomal rearrangements can initiate and drive cancer progression, yet it has been challenging to evaluate their impact, especially in genetically heterogeneous solid cancers. To address this problem we developed HiDENSEC, a new computational framework for analyzing chromatin conformation capture in heterogeneous samples that can infer somatic copy number alterations, characterize large-scale chromosomal rearrangements, and estimate cancer cell fractions. After validating HiDENSEC with in silico and in vitro controls, we used it to characterize chromosome-scale evolution during melanoma progression in formalin-fixed tumor samples from three patients. The resulting comprehensive annotation of the genomic events includes copy number neutral translocations that disrupt tumor suppressor genes such as NF1, whole chromosome arm exchanges that result in loss of CDKN2A, and whole-arm copy-number neutral loss of homozygosity involving PTEN. These findings show that large-scale chromosomal rearrangements occur throughout cancer evolution and that characterizing these events yields insights into drivers of melanoma progression.

MGA-seq: robust identification of extrachromosomal DNA and genetic variants using multiple genetic abnormality sequencing

Genomic abnormalities are strongly associated with cancer and infertility. In this study, we develop a simple and efficient method - multiple genetic abnormality sequencing (MGA-Seq) - to simultaneously detect structural variation, copy number variation, single-nucleotide polymorphism, homogeneously staining regions, and extrachromosomal DNA (ecDNA) from a single tube. MGA-Seq directly sequences proximity-ligated genomic fragments, yielding a dataset with concurrent genome three-dimensional and whole-genome sequencing information, enabling approximate localization of genomic structural variations and facilitating breakpoint identification. Additionally, by utilizing MGA-Seq, we map focal amplification and oncogene coamplification, thus facilitating the exploration of ecDNA's transcriptional regulatory function.

Detection of copy number variations based on a local distance using next-generation sequencing data

As one of the main types of structural variation in the human genome, copy number variation (CNV) plays an important role in the occurrence and development of human cancers. Next-generation sequencing (NGS) technology can provide base-level resolution, which provides favorable conditions for the accurate detection of CNVs. However, it is still a very challenging task to accurately detect CNVs from cancer samples with different purity and low sequencing coverage. Local distance-based CNV detection (LDCNV), an innovative computational approach to predict CNVs using NGS data, is proposed in this work. LDCNV calculates the average distance between each read depth (RD) and its k nearest neighbors (KNNs) to define the distance of KNNs of each RD, and the average distance between the KNNs for each RD to define their internal distance. Based on the above definitions, a local distance score is constructed using the ratio between the distance of KNNs and the internal distance of KNNs for each RD. The local distance scores are used to fit a normal distribution to evaluate the significance level of each RDS, and then use the hypothesis test method to predict the CNVs. The performance of the proposed method is verified with simulated and real data and compared with several popular methods. The experimental results show that the proposed method is superior to various other techniques. Therefore, the proposed method can be helpful for cancer diagnosis and targeted drug development.

Proteogenomic insights suggest druggable pathways in endometrial carcinoma

We characterized a prospective endometrial carcinoma (EC) cohort containing 138 tumors and 20 enriched normal tissues using 10 different omics platforms. Targeted quantitation of two peptides can predict antigen processing and presentation machinery activity, and may inform patient selection for immunotherapy. Association analysis between MYC activity and metformin treatment in both patients and cell lines suggests a potential role for metformin treatment in non-diabetic patients with elevated MYC activity. PIK3R1 in-frame indels are associated with elevated AKT phosphorylation and increased sensitivity to AKT inhibitors. CTNNB1 hotspot mutations are concentrated near phosphorylation sites mediating pS45-induced degradation of β-catenin, which may render Wnt-FZD antagonists ineffective. Deep learning accurately predicts EC subtypes and mutations from histopathology images, which may be useful for rapid diagnosis. Overall, this study identified molecular and imaging markers that can be further investigated to guide patient stratification for more precise treatment of EC.

Integrative multi-omic cancer profiling reveals DNA methylation patterns associated with therapeutic vulnerability and cell-of-origin

DNA methylation plays a critical role in establishing and maintaining cellular identity. However, it is frequently dysregulated during tumor development and is closely intertwined with other genetic alterations. Here, we leveraged multi-omic profiling of 687 tumors and matched non-involved adjacent tissues from the kidney, brain, pancreas, lung, head and neck, and endometrium to identify aberrant methylation associated with RNA and protein abundance changes and build a Pan-Cancer catalog. We uncovered lineage-specific epigenetic drivers including hypomethylated FGFR2 in endometrial cancer. We showed that hypermethylated STAT5A is associated with pervasive regulon downregulation and immune cell depletion, suggesting that epigenetic regulation of STAT5A expression constitutes a molecular switch for immunosuppression in squamous tumors. We further demonstrated that methylation subtype-enrichment information can explain cell-of-origin, intra-tumor heterogeneity, and tumor phenotypes. Overall, we identified cis-acting DNA methylation events that drive transcriptional and translational changes, shedding light on the tumor's epigenetic landscape and the role of its cell-of-origin.

Pan-cancer proteogenomics connects oncogenic drivers to functional states

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

Acquired resistance to anti-PD1 therapy in patients with NSCLC associates with immunosuppressive T cell phenotype

Immune checkpoint inhibitor treatment has the potential to prolong survival in non-small cell lung cancer (NSCLC), however, some of the patients develop resistance following initial response. Here, we analyze the immune phenotype of matching tumor samples from a cohort of NSCLC patients showing good initial response to immune checkpoint inhibitors, followed by acquired resistance at later time points. By using imaging mass cytometry and whole exome and RNA sequencing, we detect two patterns of resistance¨: One group of patients is characterized by reduced numbers of tumor-infiltrating CD8+ T cells and reduced expression of PD-L1 after development of resistance, whereas the other group shows high CD8+ T cell infiltration and high expression of PD-L1 in addition to markedly elevated expression of other immune-inhibitory molecules. In two cases, we detect downregulation of type I and II IFN pathways following progression to resistance, which could lead to an impaired anti-tumor immune response. This study thus captures the development of immune checkpoint inhibitor resistance as it progresses and deepens our mechanistic understanding of immunotherapy response in NSCLC.

HATCHet2: clone- and haplotype-specific copy number inference from bulk tumor sequencing data

Multi-region DNA sequencing of primary tumors and metastases from individual patients helps identify somatic aberrations driving cancer development. However, most methods to infer copy-number aberrations (CNAs) analyze individual samples. We introduce HATCHet2 to identify haplotype- and clone-specific CNAs simultaneously from multiple bulk samples. HATCHet2 introduces a novel statistic, the mirrored haplotype B-allele frequency (mhBAF), to identify mirrored-subclonal CNAs having different numbers of copies of parental haplotypes in different tumor clones. HATCHet2 also has high accuracy in identifying focal CNAs and extends the earlier HATCHet method in several directions. We demonstrate HATCHet2's improved accuracy using simulations and a single-cell sequencing dataset. HATCHet2 analysis of 50 prostate cancer samples from 10 patients reveals previously-unreported mirrored-subclonal CNAs affecting cancer genes.

Detection of trait-associated structural variations using short-read sequencing

Genomic structural variation (SV) affects genetic and phenotypic characteristics in diverse organisms, but the lack of reliable methods to detect SV has hindered genetic analysis. We developed a computational algorithm (MOPline) that includes missing call recovery combined with high-confidence SV call selection and genotyping using short-read whole-genome sequencing (WGS) data. Using 3,672 high-coverage WGS datasets, MOPline stably detected ∼16,000 SVs per individual, which is over ∼1.7-3.3-fold higher than previous large-scale projects while exhibiting a comparable level of statistical quality metrics. We imputed SVs from 181,622 Japanese individuals for 42 diseases and 60 quantitative traits. A genome-wide association study with the imputed SVs revealed 41 top-ranked or nearly top-ranked genome-wide significant SVs, including 8 exonic SVs with 5 novel associations and enriched mobile element insertions. This study demonstrates that short-read WGS data can be used to identify rare and common SVs associated with a variety of traits.

Epigenetic and transcriptomic characterization reveals progression markers and essential pathways in clear cell renal cell carcinoma

Identifying tumor-cell-specific markers and elucidating their epigenetic regulation and spatial heterogeneity provides mechanistic insights into cancer etiology. Here, we perform snRNA-seq and snATAC-seq in 34 and 28 human clear cell renal cell carcinoma (ccRCC) specimens, respectively, with matched bulk proteogenomics data. By identifying 20 tumor-specific markers through a multi-omics tiered approach, we reveal an association between higher ceruloplasmin (CP) expression and reduced survival. CP knockdown, combined with spatial transcriptomics, suggests a role for CP in regulating hyalinized stroma and tumor-stroma interactions in ccRCC. Intratumoral heterogeneity analysis portrays tumor cell-intrinsic inflammation and epithelial-mesenchymal transition (EMT) as two distinguishing features of tumor subpopulations. Finally, BAP1 mutations are associated with widespread reduction of chromatin accessibility, while PBRM1 mutations generally increase accessibility, with the former affecting five times more accessible peaks than the latter. These integrated analyses reveal the cellular architecture of ccRCC, providing insights into key markers and pathways in ccRCC tumorigenesis.

Potent acyl-CoA synthetase 10 inhibitors kill Plasmodium falciparum by disrupting triglyceride formation

Identifying how small molecules act to kill malaria parasites can lead to new "chemically validated" targets. By pressuring Plasmodium falciparum asexual blood stage parasites with three novel structurally-unrelated antimalarial compounds (MMV665924, MMV019719 and MMV897615), and performing whole-genome sequence analysis on resistant parasite lines, we identify multiple mutations in the P. falciparum acyl-CoA synthetase (ACS) genes PfACS10 (PF3D7_0525100, M300I, A268D/V, F427L) and PfACS11 (PF3D7_1238800, F387V, D648Y, and E668K). Allelic replacement and thermal proteome profiling validates PfACS10 as a target of these compounds. We demonstrate that this protein is essential for parasite growth by conditional knockdown and observe increased compound susceptibility upon reduced expression. Inhibition of PfACS10 leads to a reduction in triacylglycerols and a buildup of its lipid precursors, providing key insights into its function. Analysis of the PfACS11 gene and its mutations point to a role in mediating resistance via decreased protein stability.

Contributions of rare and common variation to early-onset and atypical dementia risk

We collected and analyzed genomic sequencing data from individuals with clinician- diagnosed early-onset or atypical dementia. Thirty-two patients were previously described, with sixty-eight newly described in this report. Of those sixty-eight, sixty-two patients reported Caucasian, non-Hispanic ethnicity and six reported as African American, non-Hispanic. Fifty-three percent of patients had a returnable variant. Five patients harbored a pathogenic variant as defined by the American College of Medical Genetics criteria for pathogenicity. A polygenic risk score was calculated for Alzheimer's patients in the total cohort and compared to the scores of a late-onset Alzheimer's cohort and a control set. Patients with early-onset Alzheimer's had higher non- APOE polygenic risk scores than patients with late onset Alzheimer's, supporting the conclusion that both rare and common genetic variation associate with early-onset neurodegenerative disease risk.

A shortest path-based approach for copy number variation detection from next-generation sequencing data

Copy number variation (CNV) is one of the main structural variations in the human genome and accounts for a considerable proportion of variations. As CNVs can directly or indirectly cause cancer, mental illness, and genetic disease in humans, their effective detection in humans is of great interest in the fields of oncogene discovery, clinical decision-making, bioinformatics, and drug discovery. The advent of next-generation sequencing data makes CNV detection possible, and a large number of CNV detection tools are based on next-generation sequencing data. Due to the complexity (e.g., bias, noise, alignment errors) of next-generation sequencing data and CNV structures, the accuracy of existing methods in detecting CNVs remains low. In this work, we design a new CNV detection approach, called shortest path-based Copy number variation (SPCNV), to improve the detection accuracy of CNVs. SPCNV calculates the k nearest neighbors of each read depth and defines the shortest path, shortest path relation, and shortest path cost sets based on which further calculates the mean shortest path cost of each read depth and its k nearest neighbors. We utilize the ratio between the mean shortest path cost for each read depth and the mean of the mean shortest path cost of its k nearest neighbors to construct a relative shortest path score formula that is able to determine a score for each read depth. Based on the score profile, a boxplot is then applied to predict CNVs. The performance of the proposed method is verified by simulation data experiments and compared against several popular methods of the same type. Experimental results show that the proposed method achieves the best balance between recall and precision in each set of simulated samples. To further verify the performance of the proposed method in real application scenarios, we then select real sample data from the 1,000 Genomes Project to conduct experiments. The proposed method achieves the best F1-scores in almost all samples. Therefore, the proposed method can be used as a more reliable tool for the routine detection of CNVs.

Histopathologic and proteogenomic heterogeneity reveals features of clear cell renal cell carcinoma aggressiveness

Clear cell renal cell carcinomas (ccRCCs) represent ∼75% of RCC cases and account for most RCC-associated deaths. Inter- and intratumoral heterogeneity (ITH) results in varying prognosis and treatment outcomes. To obtain the most comprehensive profile of ccRCC, we perform integrative histopathologic, proteogenomic, and metabolomic analyses on 305 ccRCC tumor segments and 166 paired adjacent normal tissues from 213 cases. Combining histologic and molecular profiles reveals ITH in 90% of ccRCCs, with 50% demonstrating immune signature heterogeneity. High tumor grade, along with BAP1 mutation, genome instability, increased hypermethylation, and a specific protein glycosylation signature define a high-risk disease subset, where UCHL1 expression displays prognostic value. Single-nuclei RNA sequencing of the adverse sarcomatoid and rhabdoid phenotypes uncover gene signatures and potential insights into tumor evolution. In vitro cell line studies confirm the potential of inhibiting identified phosphoproteome targets. This study molecularly stratifies aggressive histopathologic subtypes that may inform more effective treatment strategies.

Integrated multi-omics profiling yields a clinically relevant molecular classification for esophageal squamous cell carcinoma

Integrated molecular analysis of human cancer has yielded molecular classification for precise management of cancer patients. Here, we analyzed the whole genomic, epigenomic, transcriptomic, and proteomic data of 155 esophageal squamous cell carcinomas (ESCCs). Multi-omics analysis led to the classification of ESCCs into four subtypes: cell cycle pathway activation, NRF2 oncogenic activation, immune suppression (IS), and immune modulation (IM). IS and IM cases were highly immune infiltrated but differed in the type and distribution of immune cells. IM cases showed better response to immune checkpoint blockade therapy than other subtypes in a clinical trial. We further developed a classifier with 28 features to identify the IM subtype, which predicted anti-PD-1 therapy response with 85.7% sensitivity and 90% specificity. These results emphasize the clinical value of unbiased molecular classification based on multi-omics data and have the potential to further improve the understanding and treatment of ESCC.

Read Count Patterns and Detection of Cancerous Copy Number Alterations in Plasma Cell-Free DNA Whole Exome Sequencing Data for Advanced Non-Small Cell Lung Cancer

Plasma cell-free DNA (cfDNA) sequencing data have been widely studied for early diagnosis and treatment response or recurrence monitoring of cancers because of the non-invasive benefits. In cancer studies, whole exome sequencing (WES) is mostly used for discovering single nucleotide variants (SNVs), but it also has the potential to detect copy number alterations (CNAs) that are mostly discovered by whole genome sequencing or microarray. In clinical settings where the quantity of the acquired blood from the patients is limited and where various sequencing experiments are not possible, providing various types of mutation information such as CNAs and SNVs using only WES will be helpful in the treatment decision. Here, we questioned whether the plasma cfDNA WES data for patients with advanced non-small cell lung cancer (NSCLC) could be exploited for CNA detection. When the read count (RC) signals of the WES data were investigated, a similar fluctuation pattern was observed among the signals of different samples, and it can be a major challenge hindering CNA detection. When these RC patterns among cfDNA were suppressed by the method we proposed, the cancerous CNAs were more distinguishable in some samples with higher cfDNA quantity. Although the potential to detect CNAs using the plasma cfDNA WES data for NSCLC patients was studied here, further studies with other cancer types, with more samples, and with more sophisticated techniques for bias correction are required to confirm our observation. In conclusion, the detection performance for cancerous CNAs can be improved by controlling RC bias, but it depends on the quantity of cfDNA in plasma.

A comprehensive benchmarking of WGS-based deletion structural variant callers

Advances in whole-genome sequencing (WGS) promise to enable the accurate and comprehensive structural variant (SV) discovery. Dissecting SVs from WGS data presents a substantial number of challenges and a plethora of SV detection methods have been developed. Currently, evidence that investigators can use to select appropriate SV detection tools is lacking. In this article, we have evaluated the performance of SV detection tools on mouse and human WGS data using a comprehensive polymerase chain reaction-confirmed gold standard set of SVs and the genome-in-a-bottle variant set, respectively. In contrast to the previous benchmarking studies, our gold standard dataset included a complete set of SVs allowing us to report both precision and sensitivity rates of the SV detection methods. Our study investigates the ability of the methods to detect deletions, thus providing an optimistic estimate of SV detection performance as the SV detection methods that fail to detect deletions are likely to miss more complex SVs. We found that SV detection tools varied widely in their performance, with several methods providing a good balance between sensitivity and precision. Additionally, we have determined the SV callers best suited for low- and ultralow-pass sequencing data as well as for different deletion length categories.

Molecular Landscape of ERBB2 Alterations in 14,956 Solid Tumors

ERBB2 abnormalities frequently occur and serve as rationale therapeutic targets in cancer. In this study, clinical and next-generation sequencing data from 14,956 patients across more than 20 tumor types were collected. A total of 406 (2.7%) patients were identified with ERBB2 amplifications, and 303 (2.0%) patients with pathogenic somatic ERBB2 mutations. ERBB2 amplifications fell most frequently in breast (15.9%) and stomach (8.3%) cancers. Somatic ERBB2 SNVs/indels occurred most common in bladder/urinary tract (7.3%) and intestine (6.1%) cancers. The top mutated ERBB2 SNVs/indels were p.Y772_A775dup (25.5%) and p.S310F/Y (19.9%). Significantly higher rates of ERBB2 SNV/indels were found in women compared to men (2.8% vs. 1.5%, p < 0.0001). CDK12 was the most common co-amplification gene with ERBB2 in cancers with a high frequency of ERBB2 amplifications. Patients with ERBB2 amplifications or mutations had higher TMB compared with patients with non-ERBB2 alterations. The study provided the landscape of ERBB2 alterations across a variety of solid tumors that may benefit from anti-HER2 agents.

SETD2 Haploinsufficiency Enhances Germinal Center-Associated AICDA Somatic Hypermutation to Drive B-cell Lymphomagenesis

SETD2 is the sole histone methyltransferase responsible for H3K36me3, with roles in splicing, transcription initiation, and DNA damage response. Homozygous disruption of SETD2 yields a tumor suppressor effect in various cancers. However, SETD2 mutation is typically heterozygous in diffuse large B-cell lymphomas. Here we show that heterozygous Setd2 deficiency results in germinal center (GC) hyperplasia and increased competitive fitness, with reduced DNA damage checkpoint activity and apoptosis, resulting in accelerated lymphomagenesis. Impaired DNA damage sensing in Setd2-haploinsufficient germinal center B (GCB) and lymphoma cells associated with increased AICDA-induced somatic hypermutation, complex structural variants, and increased translocations including those activating MYC. DNA damage was selectively increased on the nontemplate strand, and H3K36me3 loss was associated with greater RNAPII processivity and mutational burden, suggesting that SETD2-mediated H3K36me3 is required for proper sensing of cytosine deamination. Hence, Setd2 haploinsufficiency delineates a novel GCB context-specific oncogenic pathway involving defective epigenetic surveillance of AICDA-mediated effects on transcribed genes.

SIGNIFICANCE

Our findings define a B cell-specific oncogenic effect of SETD2 heterozygous mutation, which unleashes AICDA mutagenesis of nontemplate strand DNA in the GC reaction, resulting in lymphomas with heavy mutational burden. GC-derived lymphomas did not tolerate SETD2 homozygous deletion, pointing to a novel context-specific therapeutic vulnerability. This article is highlighted in the In This Issue feature, p. 1599.

Germline variants of ATG7 in familial cholangiocarcinoma alter autophagy and p62

Autophagy is a housekeeping mechanism tasked with eliminating misfolded proteins and damaged organelles to maintain cellular homeostasis. Autophagy deficiency results in increased oxidative stress, DNA damage and chronic cellular injury. Among the core genes in the autophagy machinery, ATG7 is required for autophagy initiation and autophagosome formation. Based on the analysis of an extended pedigree of familial cholangiocarcinoma, we determined that all affected family members had a novel germline mutation (c.2000C>T p.Arg659* (p.R659*)) in ATG7. Somatic deletions of ATG7 were identified in the tumors of affected individuals. We applied linked-read sequencing to one tumor sample and demonstrated that the ATG7 somatic deletion and germline mutation were located on distinct alleles, resulting in two hits to ATG7. From a parallel population genetic study, we identified a germline polymorphism of ATG7 (c.1591C>G p.Asp522Glu (p.D522E)) associated with increased risk of cholangiocarcinoma. To characterize the impact of these germline ATG7 variants on autophagy activity, we developed an ATG7-null cell line derived from the human bile duct. The mutant p.R659* ATG7 protein lacked the ability to lipidate its LC3 substrate, leading to complete loss of autophagy and increased p62 levels. Our findings indicate that germline ATG7 variants have the potential to impact autophagy function with implications for cholangiocarcinoma development.

A multi-platform reference for somatic structural variation detection

Accurate detection of somatic structural variation (SV) in cancer genomes remains a challenging problem. This is in part due to the lack of high-quality, gold-standard datasets that enable the benchmarking of experimental approaches and bioinformatic analysis pipelines. Here, we performed somatic SV analysis of the paired melanoma and normal lymphoblastoid COLO829 cell lines using four different sequencing technologies. Based on the evidence from multiple technologies combined with extensive experimental validation, we compiled a comprehensive set of carefully curated and validated somatic SVs, comprising all SV types. We demonstrate the utility of this resource by determining the SV detection performance as a function of tumor purity and sequence depth, highlighting the importance of assessing these parameters in cancer genomics projects. The truth somatic SV dataset as well as the underlying raw multi-platform sequencing data are freely available and are an important resource for community somatic benchmarking efforts.

Insights into Immune Escape During Tumor Evolution and Response to Immunotherapy Using a Rat Model of Breast Cancer

Animal models are critical for the preclinical validation of cancer immunotherapies. Unfortunately, mouse breast cancer models do not faithfully reproduce the molecular subtypes and immune environment of the human disease. In particular, there are no good murine models of estrogen receptor-positive (ER+) breast cancer, the predominant subtype in patients. Here, we show that Nitroso-N-methylurea-induced mammary tumors in outbred Sprague-Dawley rats recapitulate the heterogeneity for mutational profiles, ER expression, and immune evasive mechanisms observed in human breast cancer. We demonstrate the utility of this model for preclinical studies by dissecting mechanisms of response to immunotherapy using combination TGFBR inhibition and PD-L1 blockade. Short-term treatment of early-stage tumors induced durable responses. Gene expression profiling and spatial mapping classified tumors as inflammatory and noninflammatory, and identified IFNγ, T-cell receptor (TCR), and B-cell receptor (BCR) signaling, CD74/MHC II, and epithelium-interacting CD8+ T cells as markers of response, whereas the complement system, M2 macrophage phenotype, and translation in mitochondria were associated with resistance. We found that the expression of CD74 correlated with leukocyte fraction and TCR diversity in human breast cancer. We identified a subset of rat ER+ tumors marked by expression of antigen-processing genes that had an active immune environment and responded to treatment. A gene signature characteristic of these tumors predicted disease-free survival in patients with ER+ Luminal A breast cancer and overall survival in patients with metastatic breast cancer receiving anti-PD-L1 therapy. We demonstrate the usefulness of this preclinical model for immunotherapy and suggest examination to expand immunotherapy to a subset of patients with ER+ disease. See related Spotlight by Roussos Torres, p. 672.

Computational analysis of cancer genome sequencing data

Distilling biologically meaningful information from cancer genome sequencing data requires comprehensive identification of somatic alterations using rigorous computational methods. As the amount and complexity of sequencing data have increased, so has the number of tools for analysing them. Here, we describe the main steps involved in the bioinformatic analysis of cancer genomes, review key algorithmic developments and highlight popular tools and emerging technologies. These tools include those that identify point mutations, copy number alterations, structural variations and mutational signatures in cancer genomes. We also discuss issues in experimental design, the strengths and limitations of sequencing modalities and methodological challenges for the future.

Somatic whole genome dynamics of precancer in Barrett's esophagus reveals features associated with disease progression

While the genomes of normal tissues undergo dynamic changes over time, little is understood about the temporal-spatial dynamics of genomes in premalignant tissues that progress to cancer compared to those that remain cancer-free. Here we use whole genome sequencing to contrast genomic alterations in 427 longitudinal samples from 40 patients with stable Barrett’s esophagus compared to 40 Barrett’s patients who progressed to esophageal adenocarcinoma (ESAD). We show the same somatic mutational processes are active in Barrett’s tissue regardless of outcome, with high levels of mutation, ESAD gene and focal chromosomal alterations, and similar mutational signatures. The critical distinction between stable Barrett’s versus those who progress to cancer is acquisition and expansion of TP53−/− cell populations having complex structural variants and high-level amplifications, which are detectable up to six years prior to a cancer diagnosis. These findings reveal the timing of common somatic genome dynamics in stable Barrett’s esophagus and define key genomic features specific to progression to esophageal adenocarcinoma, both of which are critical for cancer prevention and early detection strategies.

Barrett’s esophagus is a pre-malignant condition that can progress to esophageal cancer. Here, the authors carry out whole genome sequencing of samples from patients who did or did not progress to cancer and find that mutations in many genes occur regardless of progression status, but also find features associated with progressive disease.

A neurodegenerative disease landscape of rare mutations in Colombia due to founder effects

BACKGROUND

The Colombian population, as well as those in other Latin American regions, arose from a recent tri-continental admixture among Native Americans, Spanish invaders, and enslaved Africans, all of whom passed through a population bottleneck due to widespread infectious diseases that left small isolated local settlements. As a result, the current population reflects multiple founder effects derived from diverse ancestries.

METHODS

We characterized the role of admixture and founder effects on the origination of the mutational landscape that led to neurodegenerative disorders under these historical circumstances. Genomes from 900 Colombian individuals with Alzheimer's disease (AD) [n = 376], frontotemporal lobar degeneration-motor neuron disease continuum (FTLD-MND) [n = 197], early-onset dementia not otherwise specified (EOD) [n = 73], and healthy participants [n = 254] were analyzed. We examined their global and local ancestry proportions and screened this cohort for deleterious variants in disease-causing and risk-conferring genes.

RESULTS

We identified 21 pathogenic variants in AD-FTLD related genes, and PSEN1 harbored the majority (11 pathogenic variants). Variants were identified from all three continental ancestries. TREM2 heterozygous and homozygous variants were the most common among AD risk genes (102 carriers), a point of interest because the disease risk conferred by these variants differed according to ancestry. Several gene variants that have a known association with MND in European populations had FTLD phenotypes on a Native American haplotype. Consistent with founder effects, identity by descent among carriers of the same variant was frequent.

CONCLUSIONS

Colombian demography with multiple mini-bottlenecks probably enhanced the detection of founder events and left a proportionally higher frequency of rare variants derived from the ancestral populations. These findings demonstrate the role of genomically defined ancestry in phenotypic disease expression, a phenotypic range of different rare mutations in the same gene, and further emphasize the importance of inclusiveness in genetic studies.

Progress in Methods for Copy Number Variation Profiling

Copy number variations (CNVs) are the predominant class of structural genomic variations involved in the processes of evolutionary adaptation, genomic disorders, and disease progression. Compared with single-nucleotide variants, there have been challenges associated with the detection of CNVs owing to their diverse sizes. However, the field has seen significant progress in the past 20–30 years. This has been made possible due to the rapid development of molecular diagnostic methods which ensure a more detailed view of the genome structure, further complemented by recent advances in computational methods. Here, we review the major approaches that have been used to routinely detect CNVs, ranging from cytogenetics to the latest sequencing technologies, and then cover their specific features.

Genomics landscape of 185 Streptococcus thermophilus and identification of fermentation biomarkers

Streptococcus (S.) thermophilus, an indispensable dairy starter, has been used in autochthonous as well as industrial milk fermentation. However, the genetic architecture underlying S. thermophilus traits and phenotypes is largely unknown. Here, we sequenced 185 S. thermophilus strains, isolated from natural fermented dairy products of China and Mongolia and used comparative genomic and genome wide association study to provide novel point for genetic architecture underlying its traits and phenotypes. Genome analysis of S. thermophilus showed association of phylogeny with environmental and phenotypic features and revealed clades with high acid production potential or with substantial genome decay. A few S. thermophilus isolated from areas with high chloramphenicol emissions had a chloramphenicol-resistant gene CatB8. Most importantly, we defined a growth score and identified a missense mutation G1118698T located at the gene AcnA that were both predictive of acidification capability of S. thermophilus. Our findings provide novel insight in S. thermophilus genetic traits, antibiotic resistant and predictive of acidification capability which both may had huge help in culture starter screening.

Pyrotinib in Patients with HER2-Amplified Advanced Non-Small Cell Lung Cancer: A Prospective, Multicenter, Single-Arm Trial

PURPOSE

In this study, we aimed to evaluate the efficacy and safety of pyrotinib, a pan-HER inhibitor, in HER2-amplified NSCLC patients.

PATIENTS AND METHODS

In this prospective, multicenter, single-arm trial (ChiCTR1800020262), advanced NSCLC patients with HER2 amplification, as determined by next-generation sequencing, were enrolled and administered pyrotinib orally at 400 mg per day. The primary endpoint was 6-month progression-free survival (PFS) rate. Other endpoints included objective response rate (ORR), disease control rate (DCR), PFS, overall survival (OS), and safety.

RESULTS

The enrolled cohort included twenty-seven patients with HER2 amplification. The 6-month PFS rate was 51.9% (95% confidence interval [CI] 34.0%-69.3%). The median PFS was 6.3 months (95%CI 3.0-9.6 months), and median OS was 12.5 months (95%CI 8.2-16.8 months). Pyrotinib elicited a confirmed ORR of 22.2% (95%CI 10.6%-40.8%). Patients administered pyrotinib as first-line treatment achieved an mPFS of 12.4 months. Moreover, 30.8% of the patients who had progressed on epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor (TKI) responded to pyrotinib. Patients with brain metastases had an ORR of 40%. Treatment-related adverse events (TRAEs) occurred in all patients (grade 3, 22.2%), but no grade 4 or higher TRAEs were documented. Diarrhea was the most frequent TRAE (all, 92.6%; grade 3, 7.4%). Loss of HER2 amplification was detected upon disease progression.

CONCLUSIONS

Pyrotinib provided antitumor efficacy with a manageable safety profile in HER2-amplified NSCLC patients.

Proteogenomic characterization of pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive cancer with poor patient survival. Toward understanding the underlying molecular alterations that drive PDAC oncogenesis, we conducted comprehensive proteogenomic analysis of 140 pancreatic cancers, 67 normal adjacent tissues, and 9 normal pancreatic ductal tissues. Proteomic, phosphoproteomic, and glycoproteomic analyses were used to characterize proteins and their modifications. In addition, whole-genome sequencing, whole-exome sequencing, methylation, RNA sequencing (RNA-seq), and microRNA sequencing (miRNA-seq) were performed on the same tissues to facilitate an integrated proteogenomic analysis and determine the impact of genomic alterations on protein expression, signaling pathways, and post-translational modifications. To ensure robust downstream analyses, tumor neoplastic cellularity was assessed via multiple orthogonal strategies using molecular features and verified via pathological estimation of tumor cellularity based on histological review. This integrated proteogenomic characterization of PDAC will serve as a valuable resource for the community, paving the way for early detection and identification of novel therapeutic targets.

NiPTUNE: an automated pipeline for noninvasive prenatal testing in an accurate, integrative and flexible framework

Noninvasive prenatal testing (NIPT) consists of determining fetal aneuploidies by quantifying copy number alteration from the sequencing of cell-free DNA (cfDNA) from maternal blood. Due to the presence of cfDNA of fetal origin in maternal blood, in silico approaches have been developed to accurately predict fetal aneuploidies. Although NIPT is becoming a new standard in prenatal screening of chromosomal abnormalities, there are no integrated pipelines available to allow rapid, accurate and standardized data analysis in any clinical setting. Several tools have been developed, however often optimized only for research purposes or requiring enormous amount of retrospective data, making hard their implementation in a clinical context. Furthermore, no guidelines have been provided on how to accomplish each step of the data analysis to achieve reliable results. Finally, there is no integrated pipeline to perform all steps of NIPT analysis. To address these needs, we tested several tools for performing NIPT data analysis. We provide extensive benchmark of tools performances but also guidelines for running them. We selected the best performing tools that we benchmarked and gathered them in a computational pipeline. NiPTUNE is an open source python package that includes methods for fetal fraction estimation, a novel method for accurate gender prediction, a principal component analysis based strategy for quality control and fetal aneuploidies prediction. NiPTUNE is constituted by seven modules allowing the user to run the entire pipeline or each module independently. Using two cohorts composed by 1439 samples with 31 confirmed aneuploidies, we demonstrated that NiPTUNE is a valuable resource for NIPT analysis.

Landscape of IDH1/2 mutations in Chinese patients with solid tumors: A pan-cancer analysis

BACKGROUND

Isocitrate dehydrogenase (IDH) is an enzyme family involved in cell aerobic metabolism of tricarboxylic acid cycle. However, the landscape of IDH mutations in pan-cancer has not been fully characterized.

METHODS

Tissue or blood samples were subjected to next-generation sequencing (NGS) for detection the IDH mutation.

RESULTS

A total of 28.868 patients from more than 20 solid tumor species were analyzed. A total of 374 cases (1.30%) with IDH mutations were identified. Among all the IDH mutations cases, 80 (21.4%) were biliary tract cancer (BTC), 80 (21.4%) were lung cancer, 57 (15.2%) were liver cancer, and 42 (11.2%) were colorectal cancer. The most common IDH variant were IDH1 and IDH2 which were discovered in 0.81% cases and 0.47% cases, respectively. However, there were significant differences in IDH1 and IDH2 mutation frequency among different tumor species (p = 0.0003). Of the patients with IDH1 mutations, about 53.0% of these mutations occur in codons 132. Codons 172 (25.4%) was high-frequency mutation subtypes in IDH2 mutation. TP53, PBRM1, and BAP1 were the most significantly mutated genes in BTC which were different from others cancer. Moreover, TMB were significantly higher in lung cancer, colorectal cancer, and gastric cancer than BTC (p = 0.0164, p < 0.0001, p = 0.0067, respectively) and BTC patients with IDH mutation had lower TMB compared with wild-type IDH.

CONCLUSION

Somatic IDH mutation was found in multiple solid tumors and IDH would be a driver gene in BTC.

A Cluster-Based Approach for the Discovery of Copy Number Variations From Next-Generation Sequencing Data

The next-generation sequencing technology offers a wealth of data resources for the detection of copy number variations (CNVs) at a high resolution. However, it is still challenging to correctly detect CNVs of different lengths. It is necessary to develop new CNV detection tools to meet this demand. In this work, we propose a new CNV detection method, called CBCNV, for the detection of CNVs of different lengths from whole genome sequencing data. CBCNV uses a clustering algorithm to divide the read depth segment profile, and assigns an abnormal score to each read depth segment. Based on the abnormal score profile, Tukey's fences method is adopted in CBCNV to forecast CNVs. The performance of the proposed method is evaluated on simulated data sets, and is compared with those of several existing methods. The experimental results prove that the performance of CBCNV is better than those of several existing methods. The proposed method is further tested and verified on real data sets, and the experimental results are found to be consistent with the simulation results. Therefore, the proposed method can be expected to become a routine tool in the analysis of CNVs from tumor-normal matched samples.

Patient-derived tumor-like cell clusters for drug testing in cancer therapy

Several patient-derived tumor models emerged recently as robust preclinical drug-testing platforms. However, their potential to guide clinical therapy remained unclear. Here, we report a model called patient-derived tumor-like cell clusters (PTCs). PTCs result from the self-assembly and proliferation of primary epithelial, fibroblast, and immune cells, which structurally and functionally recapitulate original tumors. PTCs enabled us to accomplish personalized drug testing within 2 weeks after obtaining the tumor samples. The defined culture conditions and drug concentrations in the PTC model facilitate its clinical application in precision oncology. PTC tests of 59 patients with gastric, colorectal, or breast cancers revealed an overall accuracy of 93% in predicting their clinical outcomes. We implemented PTC to guide chemotherapy selection for a patient with mucinous rectal adenocarcinoma who experienced recurrence with metastases after conventional therapy. After three cycles of a nonconventional therapy identified by the PTC, the patient showed a positive response. These findings need to be validated in larger clinical trials, but they suggest that the PTC model could be prospectively implemented in clinical decision-making for therapy selection.

Evaluation of tools for identifying large copy number variations from ultra-low-coverage whole-genome sequencing data

BACKGROUND

Detection of copy number variations (CNVs) from high-throughput next-generation whole-genome sequencing (WGS) data has become a widely used research method during the recent years. However, only a little is known about the applicability of the developed algorithms to ultra-low-coverage (0.0005-0.8×) data that is used in various research and clinical applications, such as digital karyotyping and single-cell CNV detection.

RESULT

Here, the performance of six popular read-depth based CNV detection algorithms (BIC-seq2, Canvas, CNVnator, FREEC, HMMcopy, and QDNAseq) was studied using ultra-low-coverage WGS data. Real-world array- and karyotyping kit-based validation were used as a benchmark in the evaluation. Additionally, ultra-low-coverage WGS data was simulated to investigate the ability of the algorithms to identify CNVs in the sex chromosomes and the theoretical minimum coverage at which these tools can accurately function. Our results suggest that while all the methods were able to detect large CNVs, many methods were susceptible to producing false positives when smaller CNVs (< 2 Mbp) were detected. There was also significant variability in their ability to identify CNVs in the sex chromosomes. Overall, BIC-seq2 was found to be the best method in terms of statistical performance. However, its significant drawback was by far the slowest runtime among the methods (> 3 h) compared with FREEC (~ 3 min), which we considered the second-best method.

CONCLUSIONS

Our comparative analysis demonstrates that CNV detection from ultra-low-coverage WGS data can be a highly accurate method for the detection of large copy number variations when their length is in millions of base pairs. These findings facilitate applications that utilize ultra-low-coverage CNV detection.

Whole genome sequencing reveals a frameshift mutation and a large deletion in YY1AP1 in a girl with a panvascular artery disease

Background

Rare diseases are pathologies that affect less than 1 in 2000 people. They are difficult to diagnose due to their low frequency and their often highly heterogeneous symptoms. Rare diseases have in general a high impact on the quality of life and life expectancy of patients, which are in general children or young people. The advent of high-throughput sequencing techniques has improved diagnosis in several different areas, from pediatrics, achieving a diagnostic rate of 41% with whole genome sequencing (WGS) and 36% with whole exome sequencing, to neurology, achieving a diagnostic rate between 47 and 48.5% with WGS. This evidence has encouraged our group to pursue a molecular diagnosis using WGS for this and several other patients with rare diseases.

Results

We used whole genome sequencing to achieve a molecular diagnosis of a 7-year-old girl with a severe panvascular artery disease that remained for several years undiagnosed. We found a frameshift variant in one copy and a large deletion involving two exons in the other copy of a gene called YY1AP1. This gene is related to Grange syndrome, a recessive rare disease, whose symptoms include stenosis or occlusion of multiple arteries, congenital heart defects, brachydactyly, syndactyly, bone fragility, and learning disabilities. Bioinformatic analyses propose these mutations as the most likely cause of the disease, according to its frequency, in silico predictors, conservation analyses, and effect on the protein product. Additionally, we confirmed one mutation in each parent, supporting a compound heterozygous status in the child.

Conclusions

In general, we think that this finding can contribute to the use of whole genome sequencing as a diagnosis tool of rare diseases, and in particular, it can enhance the set of known mutations associated with different diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-021-00328-1.

Co-evolution of tumor and immune cells during progression of multiple myeloma

Multiple myeloma (MM) is characterized by the uncontrolled proliferation of plasma cells. Despite recent treatment advances, it is still incurable as disease progression is not fully understood. To investigate MM and its immune environment, we apply single cell RNA and linked-read whole genome sequencing to profile 29 longitudinal samples at different disease stages from 14 patients. Here, we collect 17,267 plasma cells and 57,719 immune cells, discovering patient-specific plasma cell profiles and immune cell expression changes. Patients with the same genetic alterations tend to have both plasma cells and immune cells clustered together. By integrating bulk genomics and single cell mapping, we track plasma cell subpopulations across disease stages and find three patterns: stability (from precancer to diagnosis), and gain or loss (from diagnosis to relapse). In multiple patients, we detect “B cell-featured” plasma cell subpopulations that cluster closely with B cells, implicating their cell of origin. We validate AP-1 complex differential expression (JUN and FOS) in plasma cell subpopulations using CyTOF-based protein assays, and integrated analysis of single-cell RNA and CyTOF data reveals AP-1 downstream targets (IL6 and IL1B) potentially leading to inflammation regulation. Our work represents a longitudinal investigation for tumor and microenvironment during MM progression and paves the way for expanding treatment options.

Clonal evolution in multiple myeloma (MM) needs to be understood in both the tumor and its microenvironment. Here the authors perform single-cell multi-omics profiling of samples from MM patients at different stages, finding transitions in the immune cell composition throughout progression.

Loss of BubR1 acetylation provokes replication stress and leads to complex chromosomal rearrangements

Accurate chromosomal segregation during mitosis is regulated by the spindle assembly checkpoint (SAC). SAC failure results in aneuploidy, a hallmark of cancer. However, many studies have suggested that aneuploidy alone is not oncogenic. We have reported that BubR1 acetylation deficiency in mice (K243R/+) caused spontaneous tumorigenesis via weakened SAC signaling and unstable chromosome-spindle attachment, resulting in massive chromosomal mis-segregation. In addition to aneuploidy, cells derived from K243R/+ mice exhibited moderate genetic instability and chromosomal translocation. Here, we investigated how the loss of BubR1 acetylation led to genetic instability and chromosomal rearrangement. To rescue all chromosomal abnormalities generated by the loss of BubR1 acetylation during development, K243R/+ mice were crossed with p53-deficient mice. Genome-wide sequencing and spectral karyotyping of tumors derived from these double-mutant mice revealed that BubR1 acetylation deficiency was associated with complex chromosomal rearrangements, including Robertsonian-like whole-arm translocations. By analyzing the telomeres and centromeres in metaphase chromosome spreads, we found that BubR1 acetylation deficiency increased the collapse of stalled replication forks, commonly referred to as replication stress, and led to DNA damage and chromosomal rearrangements. BubR1 mutations that are critical in interacting with PCAF acetyltransferase and acetylating K250, L249F and A251P, were found from human cancers. Furthermore, a subset of human cancer cells exhibiting whole-arm translocation also displayed defects in BubR1 acetylation, supporting that defects in BubR1 acetylation in mitosis contributes to tumorigenesis. Collectively, loss of BubR1 acetylation provokes replication stress, particularly at the telomeres, leading to genetic instability and chromosomal rearrangement.

Integrative reconstruction of cancer genome karyotypes using InfoGenomeR

Annotation of structural variations (SVs) and base-level karyotyping in cancer cells remains challenging. Here, we present Integrative Framework for Genome Reconstruction (InfoGenomeR)-a graph-based framework that can reconstruct individual SVs into karyotypes based on whole-genome sequencing data, by integrating SVs, total copy number alterations, allele-specific copy numbers, and haplotype information. Using whole-genome sequencing data sets of patients with breast cancer, glioblastoma multiforme, and ovarian cancer, we demonstrate the analytical potential of InfoGenomeR. We identify recurrent derivative chromosomes derived from chromosomes 11 and 17 in breast cancer samples, with homogeneously staining regions for CCND1 and ERBB2, and double minutes and breakage-fusion-bridge cycles in glioblastoma multiforme and ovarian cancer samples, respectively. Moreover, we show that InfoGenomeR can discriminate private and shared SVs between primary and metastatic cancer sites that could contribute to tumour evolution. These findings indicate that InfoGenomeR can guide targeted therapies by unravelling cancer-specific SVs on a genome-wide scale.

Whole-genome sequencing of phenotypically distinct inflammatory breast cancers reveals similar genomic alterations to non-inflammatory breast cancers

BACKGROUND

Inflammatory breast cancer (IBC) has a highly invasive and metastatic phenotype. However, little is known about its genetic drivers. To address this, we report the largest cohort of whole-genome sequencing (WGS) of IBC cases.

METHODS

We performed WGS of 20 IBC samples and paired normal blood DNA to identify genomic alterations. For comparison, we used 23 matched non-IBC samples from the Cancer Genome Atlas Program (TCGA). We also validated our findings using WGS data from the International Cancer Genome Consortium (ICGC) and the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium. We examined a wide selection of genomic features to search for differences between IBC and conventional breast cancer. These include (i) somatic and germline single-nucleotide variants (SNVs), in both coding and non-coding regions; (ii) the mutational signature and the clonal architecture derived from these SNVs; (iii) copy number and structural variants (CNVs and SVs); and (iv) non-human sequence in the tumors (i.e., exogenous sequences of bacterial origin).

RESULTS

Overall, IBC has similar genomic characteristics to non-IBC, including specific alterations, overall mutational load and signature, and tumor heterogeneity. In particular, we observed similar mutation frequencies between IBC and non-IBC, for each gene and most cancer-related pathways. Moreover, we found no exogenous sequences of infectious agents specific to IBC samples. Even though we could not find any strongly statistically distinguishing genomic features between the two groups, we did find some suggestive differences in IBC: (i) The MAST2 gene was more frequently mutated (20% IBC vs. 0% non-IBC). (ii) The TGF β pathway was more frequently disrupted by germline SNVs (50% vs. 13%). (iii) Different copy number profiles were observed in several genomic regions harboring cancer genes. (iv) Complex SVs were more frequent. (v) The clonal architecture was simpler, suggesting more homogenous tumor-evolutionary lineages.

CONCLUSIONS

Whole-genome sequencing of IBC manifests a similar genomic architecture to non-IBC. We found no unique genomic alterations shared in just IBCs; however, subtle genomic differences were observed including germline alterations in TGFβ pathway genes and somatic mutations in the MAST2 kinase that could represent potential therapeutic targets.

Proteogenomic and metabolomic characterization of human glioblastoma

Glioblastoma (GBM) is the most aggressive nervous system cancer. Understanding its molecular pathogenesis is crucial to improving diagnosis and treatment. Integrated analysis of genomic, proteomic, post-translational modification and metabolomic data on 99 treatment-naive GBMs provides insights to GBM biology. We identify key phosphorylation events (e.g., phosphorylated PTPN11 and PLCG1) as potential switches mediating oncogenic pathway activation, as well as potential targets for EGFR-, TP53-, and RB1-altered tumors. Immune subtypes with distinct immune cell types are discovered using bulk omics methodologies, validated by snRNA-seq, and correlated with specific expression and histone acetylation patterns. Histone H2B acetylation in classical-like and immune-low GBM is driven largely by BRDs, CREBBP, and EP300. Integrated metabolomic and proteomic data identify specific lipid distributions across subtypes and distinct global metabolic changes in IDH-mutated tumors. This work highlights biological relationships that could contribute to stratification of GBM patients for more effective treatment.