CNVkit: Genome-Wide Copy Number Detection and Visualization from Targeted DNA Sequencing

Haplotype-resolved de novo assembly revealed unique characteristics of alternative lengthening of telomeres in mouse embryonic stem cells

Telomeres protect chromosome ends from DNA damage responses, and their dysfunction triggers genomic alterations like chromosome fusion and rearrangement, which can lead to cellular death. Certain cells, including specific cancer cells, adopt alternative lengthening of telomere (ALT) to counteract dysfunctional telomeres and proliferate indefinitely. While telomere instability and ALT activity are likely major sources of genomic alteration, the patterns and consequences of such changes at the nucleotide level in ALT cells remain unexplored. Here we generated haplotype-resolved genome assemblies for type I ALT mouse embryonic stem cells, facilitated by highly accurate or ultra-long reads and Hi-C reads. High-quality genome revealed ALT-specific complex chromosome end structures and various genomic alterations including over 1000 structural variants (SVs). The unique sequence (mTALT) used as a template for type I ALT telomeres showed traces of being recruited into the genome, with mTALT being replicated with remarkably high accuracy. Subtelomeric regions exhibited distinct characteristics: resistance to the accumulation of SVs and small variants. We genotyped SVs at allele resolution, identifying genes (Rgs6, Dpf3 and Tacc2) crucial for maintaining ALT telomere stability. Our genome assembly-based approach elucidated the unique characteristics of ALT genome, offering insights into the genome evolution of cells surviving telomere-derived crisis.

Duchenne Muscular Dystrophy in Two Half-Brothers Due to Inherited 306 Kb Inverted Insertion of 10p15.1 into Intron 44 of the Dp427m Transcript of the DMD Gene

Duchenne muscular dystrophy (DMD) is a rare genetic disorder caused by the absence of a fully functional dystrophin protein in myocytes. In skeletal muscle, the lack of dystrophin ultimately results in muscle wasting and the replacement of myocytes with fatty or fibrous tissues. In the heart, cardiomyocytes eventually fail and cause fatal cardiomyopathy. We present a case of a male patient and his younger brother with a maternally inherited inverted insertion of approximately 306 kb of chromosome 10 in the deep intronic region between exons 44 and 45 of the DMD gene, leading to Duchenne muscular dystrophy. Chromosomal microarray, comprehensive muscular dystrophy genetic testing, and whole exome sequencing were negative. Targeted transcriptome RNA sequencing at an external lab showed no aberrant splicing. Research whole genome sequencing identified the copy number gain and insertion. Subsequent reanalysis of the RNA sequencing data showed possible aberrant splicing involving DMD exons 44-45, and research RNA sequencing revealed a fusion between the DMD gene on the minus strand of chromosome X and the PFKFB3 gene on the plus strand of chromosome 10. We demonstrate that whole genome sequencing can be valuable for identifying intronic events in the DMD gene previously undetected or not reported by traditional clinical testing.

Efficacy and safety of an oral combination therapy of niraparib and etoposide in platinum resistant/refractory ovarian cancer: a single arm, prospective, phase II study

OBJECTIVE

Non-platinum chemotherapy is used in platinum resistant/refractory ovarian cancer patients but offers limited efficacy, especially in those who develop platinum resistance after ≤2 lines of platinum based chemotherapy. This phase II study aimed to evaluate the efficacy and safety of oral niraparib plus etoposide in platinum resistant/refractory ovarian cancer.

METHODS

Platinum resistant/refractory ovarian cancer patients after ≤2 lines of platinum based chemotherapy, histologically confirmed as non-mucinous epithelial ovarian cancer, regardless of biomarker status, were eligible. Patients received niraparib with a starting dose of 200 mg/100 mg alternate once a day, and oral etoposide of 50 mg once a day, on days 1-20 of 30 days per cycle for a maximum of 6-8 cycles, followed by niraparib until disease progression or intolerable toxicity. The primary endpoint was investigator assessed progression free survival.

RESULTS

29 patients were enrolled from 22 May 2020 to 3 February 2023; 26 patients were included in the efficacy analysis set as per protocol. Median progression free survival was 4.2 months (95% confidence interval (CI) 3.9 to 4.4). Overall response rate was 26.9% (95% CI 8.7 to 45.2). Disease control rate was 57.7% (95% CI 37.3 to 78.0). Overall response rate in patients with a BRCA mutation and homologous recombination deficiency was 50% and 41.7%, respectively. Median progression free survival in patients with primary platinum resistance was 4.5 months (95% CI 3.6 to 5.3). 29 patients were included in the safety analysis set, and 8 (28%) patients experienced treatment related adverse events of grade ≥3. There was no treatment related discontinuation.

CONCLUSIONS

Niraparib combined with etoposide showed evidence of antitumor activity in platinum resistant/refractory ovarian cancer after ≤2 lines of platinum based chemotherapy, particularly in patients with a BRCA mutation, homologous recombination deficiency, or primary platinum resistance. This once-a-day oral combination was a convenient option.

TRIAL REGISTRATION NUMBER

NCT04217798.

Histopathologic, genomic, transcriptomic, and functional characteristics of eight melanocytic tumors with BRAF fusions showing stronger MAPK pathway activation compared to BRAF V600E tumors

BACKGROUND

Activating BRAF gene alterations are central to melanocytic tumor pathogenesis. A small, emerging subset of melanocytic tumors driven by BRAF fusions has distinct therapeutic implications and has been described to have Spitzoid morphology patterns. However, such morphological patterns do not encompass all cases, and little is known about the functional molecular events.

MATERIALS AND METHODS

We conducted a retrospective search through our molecular archives to identify melanocytic tumors with BRAF fusions. We reviewed clinical, histopathological, and genomic features. We further explored transcriptomic and protein-level findings.

RESULTS

Histopathologic patterns varied, with many cases without a distinctive pattern. We identified novel and diverse BRAF gene fusion partners. Differential transcriptomic analysis between low-risk BRAF fusion tumors and reference BRAF V600E tumors showed no differentially expressed genes. However, quantitatively stronger MAPK pathway activation of BRAF fusion tumors over BRAF V600E tumors was demonstrated by statistically significant stronger staining of p-ERK immunohistochemistry. Gene-specific RNA analysis shows comparable BRAF transcript levels between the two groups.

DISCUSSION AND CONCLUSION

The quantitatively stronger activation of the MAPK pathway of BRAF fusion tumors, instead of qualitatively different transcriptomes, may account for the morphology difference from conventional BRAF V600E tumors. BRAF fusions likely act through dysregulated protein function rather than RNA upregulation related to the characteristics of the fusion partners.

Evaluation and integration of cell-free DNA signatures for detection of lung cancer

Cell-free DNA (cfDNA) analysis has shown potential in detecting early-stage lung cancer based on non-genetic features. To distinguish patients with lung cancer from healthy individuals, peripheral blood were collected from 926 lung cancer patients and 611 healthy individuals followed by cfDNA extraction. Low-pass whole genome sequencing and targeted methylation sequencing were conducted and various features of cfDNA were evaluated. With our customized algorithm using the most optimal features, the ensemble stacked model was constructed, called ESim-seq (Early Screening tech with Integrated Model). In the independent validation cohort, the ESim-seq model achieved an area under the curve (AUC) of 0.948 (95 % CI: 0.915-0.981), with a sensitivity of 79.3 % (95 % CI: 71.5-87.0 %) across all stages at a specificity of 96.0 % (95 % CI: 90.6-100.0 %). Specifically, the sensitivity of the ESim-seq model was 76.5 % (95 % CI: 67.3-85.8 %) in stage I patients, 100 % (95 % CI: 100.0-100.0 %) in stage II patients, 100 % (95 % CI: 100.0-100.0 %) in stage III patients and 87.5 % (95 % CI: 64.6%-100.0 %) in stage IV patients in the independent validation cohort. Besides, we constructed LCSC model (Lung Cancer Subtype multiple Classification), which was able to accurately distinguish patients with small cell lung cancer from those with non-small cell lung cancer, achieving an AUC of 0.961 (95 % CI: 0.949-0.957). The present study has established a framework for assessing cfDNA features and demonstrated the benefits of integrating multiple features for early detection of lung cancer.

Coordinated inheritance of extrachromosomal DNAs in cancer cells

The chromosomal theory of inheritance dictates that genes on the same chromosome segregate together while genes on different chromosomes assort independently1. Extrachromosomal DNAs (ecDNAs) are common in cancer and drive oncogene amplification, dysregulated gene expression and intratumoural heterogeneity through random segregation during cell division2,3. Distinct ecDNA sequences, termed ecDNA species, can co-exist to facilitate intermolecular cooperation in cancer cells4. How multiple ecDNA species within a tumour cell are assorted and maintained across somatic cell generations is unclear. Here we show that cooperative ecDNA species are coordinately inherited through mitotic co-segregation. Imaging and single-cell analyses show that multiple ecDNAs encoding distinct oncogenes co-occur and are correlated in copy number in human cancer cells. ecDNA species are coordinately segregated asymmetrically during mitosis, resulting in daughter cells with simultaneous copy-number gains in multiple ecDNA species before any selection. Intermolecular proximity and active transcription at the start of mitosis facilitate the coordinated segregation of ecDNA species, and transcription inhibition reduces co-segregation. Computational modelling reveals the quantitative principles of ecDNA co-segregation and co-selection, predicting their observed distributions in cancer cells. Coordinated inheritance of ecDNAs enables co-amplification of specialized ecDNAs containing only enhancer elements and guides therapeutic strategies to jointly deplete cooperating ecDNA oncogenes. Coordinated inheritance of ecDNAs confers stability to oncogene cooperation and novel gene regulatory circuits, allowing winning combinations of epigenetic states to be transmitted across cell generations.

Deciphering the toxicity of polyhexamethylene guanidine phosphate in lung carcinogenesis: Mutational profiles and molecular mechanisms

Polyhexamethylene guanidine (PHMG) is widely utilized in personal hygiene products due to its bactericidal, non-volatile, and hydrophilic properties. However, the long-term toxic effects and underlying mechanisms associated with respiratory exposure to the commonly used form, PHMG phosphate (PHMG-p), are still insufficiently understood. This study aims to elucidate the types of pulmonary lesions and the incidence of lung cancer associated with varying concentrations of PHMG-p and observation periods, along with the molecular mechanisms underlying this relationship. To assess these effects, CT scans and pathological analyses were conducted for up to 54 weeks following initial exposure to PHMG-p. Furthermore, to investigate the underlying causes of pulmonary toxicity, TGF-beta-activated kinase 1 was identified as a PHMG-p-binding protein, and its associated signaling pathways, including necroptosis, apoptosis, and MKK7, were explored. Somatic mutational signature, and gene ontology (GO) analyses were performed to investigate the genetic characteristics of PHMG-p-induced lung carcinogenesis. PHMG-p exposure led to somatic mutations in lung cancer-related genes, including TP53, SOS1, KMT2D, MDM2, ERBB2, SETD2, MET, ARID1A, RBM10, and CDKN2A as well as in genes such as RAB31, WASHC1, DDX11, ECD, STAB2, MUC2, and MUC5AC. The mutated genes were primarily associated with impaired DNA repair mechanisms. GO analysis highlighted the activation of pathways related to cell cycle checkpoints, necroptosis, MAPK, and idiopathic pulmonary fibrosis, while also revealing the suppression of signaling pathways associated with natural killer cells, GADD45, LXR/RXR activation, and IL-15 production. Gain-of-function experiments confirmed the oncogenic roles of PLAU and HMGA2, as well as the tumor-suppressive functions of TBX4 and GPX3. These findings suggest that PHMG-p activates necroptosis and MAPK signaling, increases the frequency of somatic mutations, and inhibits apoptosis, thus fostering an environment conducive to carcinogenesis. This underscores the importance of understanding the potential health risks associated with PHMG-p exposure and provides insights for future research and regulatory considerations regarding the safety of personal hygiene products.

Stem cell factor and erythropoietin-independent production of cultured reticulocytes

Cultured reticulocytes can supplement transfusion needs and offer promise for drug delivery and immune tolerization. They can be produced from induced pluripotent stem cells (iPSC), but the 45-day culture time and cytokine costs make largescale production prohibitive. To overcome these limitations, we have generated iPSC that express constitutive stem cell factor (SCF) receptor and jak2 adaptor alleles. We show that iPSC lines carrying these alleles can differentiate into self-renewing erythroblasts that can proliferate for up to 70 cell-doubling in a cost-effective, chemically-defined, albumin- and cytokine-free medium. These kitjak2 self-renewing erythroblasts retain the ability to enucleate at a high rate up to senescence. Kitjak2-derived cultured reticulocytes should be safe for transfusion because they can be irradiated to eliminate residual nucleated cells. The kitjak2 cells express blood group 0 and test negative for RhD and other clinically significant red blood cell antigens and have sufficient proliferation capacity to meet global red blood cell needs.

Enhancing transcription-replication conflict targets ecDNA-positive cancers

Extrachromosomal DNA (ecDNA) presents a major challenge for cancer patients. ecDNA renders tumours treatment resistant by facilitating massive oncogene transcription and rapid genome evolution, contributing to poor patient survival1-7. At present, there are no ecDNA-specific treatments. Here we show that enhancing transcription-replication conflict enables targeted elimination of ecDNA-containing cancers. Stepwise analyses of ecDNA transcription reveal pervasive RNA transcription and associated single-stranded DNA, leading to excessive transcription-replication conflicts and replication stress compared with chromosomal loci. Nucleotide incorporation on ecDNA is markedly slower, and replication stress is significantly higher in ecDNA-containing tumours regardless of cancer type or oncogene cargo. pRPA2-S33, a mediator of DNA damage repair that binds single-stranded DNA, shows elevated localization on ecDNA in a transcription-dependent manner, along with increased DNA double strand breaks, and activation of the S-phase checkpoint kinase, CHK1. Genetic or pharmacological CHK1 inhibition causes extensive and preferential tumour cell death in ecDNA-containing tumours. We advance a highly selective, potent and bioavailable oral CHK1 inhibitor, BBI-2779, that preferentially kills ecDNA-containing tumour cells. In a gastric cancer model containing FGFR2 amplified on ecDNA, BBI-2779 suppresses tumour growth and prevents ecDNA-mediated acquired resistance to the pan-FGFR inhibitor infigratinib, resulting in potent and sustained tumour regression in mice. Transcription-replication conflict emerges as a target for ecDNA-directed therapy, exploiting a synthetic lethality of excess to treat cancer.

Cell-free DNA from ascites identifies clinically relevant variants and tumour evolution in patients with advanced ovarian cancer

The emergence of targeted therapies has transformed ovarian cancer treatment. However, biomarker profiling for precision medicine is limited by access to quality, tumour-enriched tissue samples. The use of cell-free DNA (cfDNA) in ascites presents a potential solution to this challenge. In this study, next-generation sequencing was performed on ascites-derived cfDNA samples (26 samples from 15 human participants with ovarian cancer), with matched DNA from ascites-derived tumour cells (n = 5) and archived formalin-fixed paraffin-embedded (FFPE) tissue (n = 5). Similar tumour purity and variant detection were achieved with cfDNA compared to FFPE and ascites cell DNA. Analysis of large-scale genomic alterations, loss of heterozygosity and tumour mutation burden identified six cases of high genomic instability (including four with pathogenic BRCA1 and BRCA2 mutations). Copy number profiles and subclone prevalence changed between sequential ascites samples, particularly in a case where deletions and chromothripsis in Chr17p13.1 and Chr8q resulted in changes in clinically relevant TP53 and MYC variants over time. Ascites cfDNA identified clinically actionable information, concordant to tissue biopsies, enabling opportunistic molecular profiling. This advocates for analysis of ascites cfDNA in lieu of accessing tumour tissue via biopsy.

Insight into mammary gland development and tumor progression in an E2F5 conditional knockout mouse model

Development of breast cancer is linked to altered regulation of mammary gland developmental processes. A better understanding of normal mammary gland development can thus reveal possible mechanisms of how normal cells are re-programmed to become malignant. E2Fs 1-4 are part of the E2F transcription factor family with varied roles in mammary development, but little is known about the role of E2F5. A combination of scRNAseq and predictive signature tools demonstrated the presence of E2F5 in the mammary gland and showed changes in predicted activity during the various phases of mammary gland development. Testing the hypothesis that E2F5 regulates mammary function, we generated a mammary-specific E2F5 knockout mouse model, resulting in modest mammary gland development changes. However, after a prolonged latency the E2F5 conditional knockout mice developed highly metastatic mammary tumors. Whole genome sequencing revealed significant intertumor heterogeneity. RNAseq and protein analysis identified altered levels of Cyclin D1, with similarities to MMTV-Neu tumors, suggesting that E2F5 conditional knockout mammary glands and tumors may be dependent on Cyclin D1. Transplantation of the tumors revealed metastases to lymph nodes that were enriched through serial transplantation in immune competent recipients. Based on these findings, we propose that loss of E2F5 leads to altered regulation of Cyclin D1, which facilitates the development of metastatic mammary tumors after long latency. More importantly, this study demonstrates that conditional loss of E2F5 in the mammary gland leads to tumor formation, revealing its role as a transcription factor regulating a network of genes that normally result in a tumor suppressor function.

Mapping extrachromosomal DNA amplifications during cancer progression

To understand the role of extrachromosomal DNA (ecDNA) amplifications in cancer progression, we detected and classified focal amplifications in 8,060 newly diagnosed primary cancers, untreated metastases and heavily pretreated tumors. The ecDNAs were detected at significantly higher frequency in untreated metastatic and pretreated tumors compared to newly diagnosed cancers. Tumors from chemotherapy-pretreated patients showed significantly higher ecDNA frequency compared to untreated cancers. In particular, tubulin inhibition associated with ecDNA increases, suggesting a role for ecDNA in treatment response. In longitudinally matched tumor samples, ecDNAs were more likely to be retained compared to chromosomal amplifications. EcDNAs shared between time points, and ecDNAs in advanced cancers were more likely to harbor localized hypermutation events compared to private ecDNAs and ecDNAs in newly diagnosed tumors. Relatively high variant allele fractions of ecDNA localized hypermutations implicated early ecDNA mutagenesis. Our findings nominate ecDNAs to provide tumors with competitive advantages during cancer progression and metastasis.

Molecular Profiling Defines Three Subtypes of Synovial Sarcoma

Synovial Sarcomas (SS) are characterized by the presence of the SS18::SSX fusion gene, which protein product induce chromatin changes through remodeling of the BAF complex. To elucidate the genomic events that drive phenotypic diversity in SS, we performed RNA and targeted DNA sequencing on 91 tumors from 55 patients. Our results were verified by proteomic analysis, public gene expression cohorts and single-cell RNA sequencing. Transcriptome profiling identified three distinct SS subtypes resembling the known histological subtypes: SS subtype I and was characterized by hyperproliferation, evasion of immune detection and a poor prognosis. SS subtype II and was dominated by a vascular-stromal component and had a significantly better outcome. SS Subtype III was characterized by biphasic differentiation, increased genomic complexity and immune suppression mediated by checkpoint inhibition, and poor prognosis despite good responses to neoadjuvant therapy. Chromosomal abnormalities were an independent significant risk factor for metastasis. KRT8 was identified as a key component for epithelial differentiation in biphasic tumors, potentially controlled by OVOL1 regulation. Our findings explain the histological grounds for SS classification and indicate that a significantly larger proportion of patients have high risk tumors (corresponding to SS subtype I) than previously believed.

Mechanisms of Response and Tolerance to Active RAS Inhibition in KRAS-Mutant Non-Small Cell Lung Cancer

Resistance to inactive state-selective RASG12C inhibitors frequently entails accumulation of RASGTP, rendering effective inhibition of active RAS potentially desirable. Here, we evaluated the antitumor activity of the RAS(ON) multiselective tricomplex inhibitor RMC-7977 and dissected mechanisms of response and tolerance in KRASG12C-mutant non-small cell lung cancer (NSCLC). Broad-spectrum reversible RASGTP inhibition with or without concurrent covalent targeting of active RASG12C yielded superior and differentiated antitumor activity across diverse comutational KRASG12C-mutant NSCLC mouse models of primary or acquired RASG12C(ON) or RASG12C(OFF) inhibitor resistance. Interrogation of time-resolved single-cell transcriptional responses established an in vivo atlas of multimodal acute and chronic RAS pathway inhibition in the NSCLC ecosystem and uncovered a regenerative mucinous transcriptional program that supports long-term tumor cell persistence. In patients with advanced KRASG12C-mutant NSCLC, the presence of mucinous histologic features portended poor response to sotorasib or adagrasib. Our results have potential implications for personalized medicine and the development of rational RAS inhibitor-anchored therapeutic strategies. Significance: Our work reveals robust and durable antitumor activity of the preclinical RAS(ON) multiselective inhibitor RMC-7977 against difficult-to-treat subsets of KRASG12C-mutant NSCLC with primary or acquired RASG12C inhibitor resistance and identifies a conserved mucinous transcriptional state that supports RAS inhibitor tolerance. See related commentary by Marasco and Misale, p. 2018.

Putative mechanisms of primary resistance to EGFR-targeted therapies: A retrospective study

BACKGROUNDS

Advanced lung adenocarcinoma (LUAD) patient with EGFR mutations often experience resistance to first-line epidermal growth factor tyrosine kinase inhibitors (EGFR-TKIs) therapy. Nonetheless, the mechanism and biomarkers of primary resistance remain unclear. Further exploration of independent prognostic factors will help clinicians identify patients who may not respond to EGFR-TKIs and select appropriate treatment strategies.

METHODS

A retrospective study involving 124 patients with stage IV LUAD harboring a common sensitizing EGFR mutation (exon 19 deletion or L858R mutation) who received EGFR-TKIs as first-line therapy was performed. All participants were tested by DNA-targeted sequencing in baseline samples, and there were 19 patients with progression-free survival (PFS) ≤ 3 months (cohort 1, C1, primary resistance), 22 patients with 3 < PFS < 8 months (cohort 2, C2, poor response) without known mutations associated with resistance, and 83 patients with PFS ≥ 8 months (cohort 3, C3, normal).

RESULTS

The most commonly mutated genes at baseline in patients prior to treatment within the entire study population. were TP53 (65 %), MYC (19 %), CDKN2A (12 %), MUC16 (12 %) and RBM10 (12 %). The baseline characteristics, except for the proportions of patients with EGFR L858R mutation and exon 19 deletion in C1 plus C2 compared to C3 (p = 0.036), were not significantly different among the cohorts. The frequencies of PIK3C2G, STK11, EPAS1, RARA and BTG2 variation were significantly higher in C1, the primary resistance group. Multivariate Cox analysis revealed that PIK3C2G (HR 15.70 95 % CI 3.24-76.05, p < 0.001), STK11 (HR 17.04, 95 % CI 3.68-78.92, p < 0.001), EPAS1 (HR 11.99, 95 % CI 2.57-56.03, p = 0.002), and BTG2 amplification (HR 9.53, 95 % CI 1.67-54.28, p = 0.011) were significantly associated with shorter PFS.

CONCLUSIONS

The genomic landscape varies significantly among patients with LUAD, which should be considered when making personalized treatment decisions. This information could provide insights into molecular changes and their effects on clinical treatment in diverse patients with LUAD harboring sensitizing EGFR mutations.

Molecular and Clinicopathologic Characterization of HER2-overexpressed Squamous Cell Carcinoma of the Cervix

HER2 amplification in cervical cancer has been associated with worse clinical prognosis and a potential favorable response to HER2 inhibitors. Immunohistochemistry for the HER2 receptor is a universally accepted surrogate test for HER2 amplification, but no standardized scoring system currently exists for cervical carcinomas. In this study, we investigated HER2 overexpression in cervical squamous cell carcinoma and correlated it with HER2 amplification using fluorescence in situ hybridization (FISH) and molecular methods. Seventy-two cases of human papillomavirus-associated cervical cancer were retrospectively reviewed, and at least 2 representative tumor sections were stained for HER2. HER2 scoring was performed using the 2018 American Society of Clinical Oncology/College of American Pathologist breast cancer criteria, and cases with equivocal (2+) to positive (3+) expression were analyzed for HER2 amplification using FISH and next-generation sequencing. The average patient age was 50 yrs (range: 27-85 yr), with most patients being African American (73.6%) and diagnosed at FIGO stage I (65.3%). Nineteen (26.4%) had equivocal HER2 expression and 4 (5.5%) showed positive expression. Three of the 4 cases with positive expression had enough tumors for FISH, and all 3 were amplified. Three cases with equivocal expression showed HER2 polysomy on FISH, and none showed HER2 amplification. Late clinical stage, high tumor grade, and regional lymph node metastasis were significantly correlated with HER2 overexpression and HER2 amplification. Next-generation sequencing of the 3 HER2-amplified tumors showed amplification of various genes, including CD274, JAK2, BIRC3, and ERBB2, and a PIK3CA missense mutation. In summary, HER2 immunohistochemistry is a reliable predictive marker of HER2 amplification in cervical cancer.

Tumor organoids improve mutation detection of pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) presents challenges in detecting somatic mutations due to its complex cellular composition. This study investigated the utility of patient-derived organoids (PDOs) to overcome these obstacles and enhance somatic mutation identification. Surgically resected PDAC tumors and their paired PDOs from 21 patients were examined. Whole-exome sequencing (WES) of tumor tissue, organoids, and peripheral blood mononuclear cells was performed to identify somatic mutations. Our findings demonstrate that PDOs retained about 80% of the somatic mutations from the original tumors, showing high concordance in mutation types. PDOs exhibited increased tumor purity and uncovered key driver mutations, aiding in identifying clinically relevant genomic alterations. Moreover, eight cycles of FOLFIRINOX treatment did not significantly alter the mutational landscape at the DNA level, indicating the stability of the mutational profile after therapeutic pressure in patients. In conclusion, PDOs are potentially important tools for exploring the somatic mutational landscape of PDAC. While they can reveal mutations that may be challenging to detect through traditional biopsy sequencing due to the inherently low tumor purity of PDAC, it is important to note that PDOs may not always fully recapitulate all mutations found in primary tumors. Despite this limitation, PDOs can still offer critical insights into the genomic complexities of PDAC, which is crucial for the development of personalized vaccines and therapies for this disease.

Multiomic and clinical analysis of multiply recurrent meningiomas reveals risk factors, underlying biology, and insights into evolution

An important subset of meningiomas behaves aggressively and is characterized by multiple recurrences. We identify clinical, genetic, and epigenetic predictors of multiply recurrent meningiomas (MRMs) and evaluate the evolution of these meningiomas in patient-matched samples. On multivariable binomial logistic regression, MRMs were significantly associated with male sex (P = 0.012), subtotal resection (P = 0.001), higher number of meningiomas on presentation (P = 0.017), and histopathological sheeting (P = 0.002). Multiomic analysis of primary meningiomas revealed that MRMs have greater copy number losses (P = 0.0313) and increased DNA methylation (P = 0.0155). In meningioma cells with knockdown of EDNRB, a locus with greater promoter methylation and decreased gene expression in MRMs had increased proliferation (P < 0.0001). MRM recurrences were found to be similar to primaries but have a greater burden of copy number gains (P < 0.0001) and increased methylation (P = 0.0045). This clinical and multiomic investigation of MRMs harbors implications for the future development of biomarkers and therapeutic agents for these challenging tumors.

WRAD core perturbation impairs DNA replication fidelity promoting immunoediting in pancreatic cancer

It is unclear how cells counteract the potentially harmful effects of uncoordinated DNA replication in the context of oncogenic stress. Here, we identify the WRAD (WDR5/RBBP5/ASH2L/DPY30) core as a modulator of DNA replication in pancreatic ductal adenocarcinoma (PDAC) models. Molecular analyses demonstrated that the WRAD core interacts with the replisome complex, with disruption of DPY30 resulting in DNA re-replication, DNA damage, and chromosomal instability (CIN) without affecting cancer cell proliferation. Consequently, in immunocompetent models, DPY30 loss induced T cell infiltration and immune-mediated clearance of highly proliferating cancer cells with complex karyotypes, thus improving anti-tumor efficacy upon anti-PD-1 treatment. In PDAC patients, DPY30 expression was associated with high tumor grade, worse prognosis, and limited response to immune checkpoint blockade. Together, our findings indicate that the WRAD core sustains genome stability and suggest that low intratumor DPY30 levels may identify PDAC patients who will benefit from immune checkpoint inhibitors.

Identification of osteoporosis genes using family studies

Osteoporosis is a multifactorial bone disease characterised by reduced bone mass and increased fracture risk. Family studies have made significant contribution in unravelling the genetics of osteoporosis. Yet, most of the underlying molecular and biological mechanisms remain unknown prompting the need for further studies. This review outlines the proper phenotyping and advanced genetic techniques in the form of high-throughput DNA sequencing used to identify genetic factors underlying monogenic osteoporosis in a family-based setting. The steps related to variant filtering prioritisation and curation are also described. From an evolutionary perspective, deleterious risk variants with higher penetrance tend to be rare as a result of negative selection. High-throughput sequencing (HTS) can identify rare variants with large effect sizes which are likely to be missed by candidate gene analysis or genome-wide association studies (GWAS) wherein common variants with small to moderate effect sizes are identified. We also describe the importance of replicating implicated genes, and possibly variants, identified following HTS to confirm their causality. Replication of the gene in other families, singletons or independent cohorts confirms that the shortlisted genes and/or variants are indeed causal. Furthermore, novel genes and/or variants implicated in monogenic osteoporosis require a thorough validation by means of in vitro and in vivo assessment. Therefore, analyses of families can continue to elucidate the genetic architecture of osteoporosis, paving the way for improved diagnostic and therapeutic strategies.

Modeling NK-cell lymphoma in mice reveals its cell-of-origin and microenvironmental changes and identifies therapeutic targets

Extranodal NK/T-cell lymphoma (ENKTCL) is an Epstein-Barr virus (EBV)-related neoplasm preferentially involving the upper aerodigestive tract. Here we show that NK-cell-specific Trp53 disruption in mice leads to the development of NK-cell lymphomas after long latency, which involve not only the hematopoietic system but also the salivary glands. Before tumor onset, Trp53 knockout causes extensive gene expression changes, resulting in immature NK-cell expansion, exclusively in the salivary glands. Both human and murine NK-cell lymphomas express tissue-resident markers, suggesting tissue-resident NK cells as their cell-of-origin. Murine NK-cell lymphomas show recurrent Myc amplifications and upregulation of MYC target gene signatures. EBV-encoded latent membrane protein 1 expression accelerates NK-cell lymphomagenesis and causes diverse microenvironmental changes, particularly myeloid propagation, through interferon-γ signaling. In turn, myeloid cells support tumor cells via CXCL16-CXCR6 signaling and its inhibition is effective against NK-cell tumors in vivo. Remarkably, KLRG1-expressing cells expand in the tumor and are capable of repopulating tumors in secondary recipients. Furthermore, targeting KLRG1 alone or combined with MYC inhibition using an eIF4 inhibitor is effective against NK-cell tumors. Therefore, our observations provide insights into the pathogenesis and highlight potential therapeutic targets, including CXCL16, KLRG1, and MYC, in ENKTCL, which can help improve its diagnostic and therapeutic strategies.

m6A-modified cenRNA stabilizes CENPA to ensure centromere integrity in cancer cells

m6A modification is best known for its critical role in controlling multiple post-transcriptional processes of the mRNAs. Here, we discovered elevated levels of m6A modification on centromeric RNA (cenRNA) in cancerous cells compared with non-cancerous cells. We then identified CENPA, an H3 variant, as an m6A reader of cenRNA. CENPA is localized at centromeres and is essential in preserving centromere integrity and function during mitosis. The m6A-modified cenRNA stabilizes centromeric localization of CENPA in cancer cells during the S phase of the cell cycle. Mutations of CENPA at the Leu61 and the Arg63 or removal of cenRNA m6A modification lead to loss of centromere-bound CENPA during S phase. This in turn results in compromised centromere integrity and abnormal chromosome separation and hinders cancer cell proliferation and tumor growth. Our findings unveil an m6A reading mechanism by CENPA that epigenetically governs centromere integrity in cancer cells, providing potential targets for cancer therapy.

Comprehensive evaluation and guidance of structural variation detection tools in chicken whole genome sequence data

BACKGROUND

Structural variations (SVs) are widespread across genome and have a great impact on evolution, disease, and phenotypic diversity. Despite the development of numerous bioinformatic tools, commonly referred to as SV callers, tailored for detecting SVs using whole genome sequence (WGS) data and employing diverse algorithms, their performance necessitates rigorous evaluation with real data and validated SVs. Moreover, a considerable proportion of these tools have been primarily designed and optimized using human genome data. Consequently, their applicability and performance in Avian species, characterized by smaller genomes and distinct genomic architectures, remain inadequately assessed.

RESULTS

We performed a comprehensive assessment of the performance of ten widely used SV callers using population-level real genomic data with the validated five common types of SVs. The performance of SV callers varies with the types and sizes of SVs. As compared with other tools, GRIDSS, Lumpy, Wham, and Manta present better detection accuracy. Pindel can detect more small SVs than others. CNVnator and CNVkit can detect more medium and large copy number variations. Given the poor consistency among different SV callers, the combination calling strategy is not recommended. All tools show poor ability in the detection of insertions (especially with size > 150 bp). At least 50× read depth is required to detect more than 80% of the SVs for most tools.

CONCLUSIONS

This study highlights the importance and necessity of using real sequencing data, rather than simulated data only, with validated SVs for SV caller evaluation. Some practical guidance and suggestions are provided for SV detection in future researches.

Whole Exome Sequencing of Intracranial Epidermoid Cysts Reveals Immune-Associated Mechanistic and Potential Targets

Background/Objectives: Intracranial Epidermoid Cysts (IECs) are rare intracranial tumors primarily treated through surgery. Cyst adherence complicates complete removal, leading to high rates of tumor progression after subtotal resection. The molecular drivers of IEC remain unknown. Consequently, advances in treatment have fallen short. Tumor genetic profiling has revealed potential targets for drug development, including FDA-approved options and reshaping treatment. The genetic landscape of IECs has not been explored. We applied Whole Exome Sequencing (WES) to IECs to gain insights into the mechanisms of oncogenesis and identify potential therapeutic targets. Methods: We performed WES on tumor tissue and matched blood samples, when available. Following GATK best practices, we conducted read processing, quality control, somatic variant calling, and copy-number inference. Data analyses and visualization were conducted in R. Results: Top altered genes are associated with the immune system and tumor microenvironment, suggesting a mechanism of immune evasion. Gene and pathway enrichment revealed a high mutation burden in genes associated with Extracellular Matrix (ECM) and PI3K-AKT-mTOR cascades. Recurrent and deleterious alterations in NOTCH2 and USP8 were identified in 50% and 30% of the cohort, respectively. Frequent amplifications in deubiquitinases and beta-defensins strengthened the involvement of immune mechanisms for oncogenic transformation. Conclusions: Top altered genes and recurrent mutations may play a role in shaping the microenvironment and modulating immune evasion in IECs. USP8 and NOTCH2 may serve as clinically relevant target for IECs. Finally, we present evidence that the crosstalk between the PI3K-Akt-mTOR and ECM signaling pathways may play a role in modulating the immune escape mechanism in IECs.

CoRAL accurately resolves extrachromosomal DNA genome structures with long-read sequencing

Extrachromosomal DNA (ecDNA) is a central mechanism for focal oncogene amplification in cancer, occurring in ∼15% of early-stage cancers and ∼30% of late-stage cancers. ecDNAs drive tumor formation, evolution, and drug resistance by dynamically modulating oncogene copy number and rewiring gene-regulatory networks. Elucidating the genomic architecture of ecDNA amplifications is critical for understanding tumor pathology and developing more effective therapies. Paired-end short-read (Illumina) sequencing and mapping have been utilized to represent ecDNA amplifications using a breakpoint graph, in which the inferred architecture of ecDNA is encoded as a cycle in the graph. Traversals of breakpoint graphs have been used to successfully predict ecDNA presence in cancer samples. However, short-read technologies are intrinsically limited in the identification of breakpoints, phasing together complex rearrangements and internal duplications, and deconvolution of cell-to-cell heterogeneity of ecDNA structures. Long-read technologies, such as from Oxford Nanopore Technologies, have the potential to improve inference as the longer reads are better at mapping structural variants and are more likely to span rearranged or duplicated regions. Here, we propose Complete Reconstruction of Amplifications with Long reads (CoRAL) for reconstructing ecDNA architectures using long-read data. CoRAL reconstructs likely cyclic architectures using quadratic programming that simultaneously optimizes parsimony of reconstruction, explained copy number, and consistency of long-read mapping. CoRAL substantially improves reconstructions in extensive simulations and 10 data sets from previously characterized cell lines compared with previous short- and long-read-based tools. As long-read usage becomes widespread, we anticipate that CoRAL will be a valuable tool for profiling the landscape and evolution of focal amplifications in tumors.

Reference-based genome assembly and comparative genomics of Calamus Brandisii Becc. for unveiling sex-specific genes for early gender detection

Calamus brandisii Becc. is an endangered rattan species indigenous to the Western Ghats of India and used in the furniture and handicraft industries. However, its dioecious nature and longer flowering time pose challenges for conservation efforts. Developing markers for early gender detection in seedlings is crucial for maintaining viable populations for in-situ and ex-situ conservation. Currently, no sex chromosomes or gender-specific genes have been reported in the species. We report the first comprehensive comparative genomics study between the male and female genomes of C. brandisii to identify polymorphisms and potential genes for gender determination. Reference-based assembly was conducted and the male and female genomes were predicted to contain 43,810 and 50,493 protein-coding genes respectively. The haploid genome size was ∼691 Mb and ∼884 Mb for male and female genomes respectively. Comparative analysis revealed significant genetic variation between the two genomes including 619,776 SNPs, 73,659 InDels, 212,123 Structural variants (SVs) and 305 copy number variations (CNVs). A total of 5 male-specific and 11 female-specific genes linked to the sex determining region was predicted. The genomic variants identified between the two genomes could be used in development of markers for early gender identification in C. brandisii for restoration programs. The gender-specific genes identified in this study also provide new insights into the mechanisms of sex determination and differentiation in rattans.

Analysis of actionable gene fusions in a large cohort of Chinese patients with colorectal cancer

Background

The prevalence of gene fusion is extremely low in unselected patients with colorectal cancer (CRC). Published data on gene fusions are limited by relatively small sample sizes, with a primary focus on Western populations. This study aimed to analyse actionable gene fusions in a large consecutive Chinese CRC population.

Methods

This study included 5,534 consecutive CRC patients from the Genecast database. Genomic profiling was performed using a panel of 769 cancer-related genes. Data for 34 CRC patients with actionable gene fusions were also collected from cBioPortal and ChimerSeq.

Results

Among 5,534 CRC patients, 54 (0.98%) had actionable gene fusions, with NTRK1/2/3 being the most common fusion (0.38%), accounting for 38.9% (21/54) of those with fusions. Actionable gene fusion enrichment was higher in patients with microsatellite instability-high (MSI-H) (6.7% vs. 0.5%, P < 0.001), RAS/BRAF wildtype (2.0% vs. 0.2%, P < 0.001) and RNF43 mutation (7.7% vs. 0.4%, P < 0.001) than in patients with microsatellite stability/MSI-low, RAS/BRAF mutation and RNF43 wildtype, respectively. When these markers were combined, the fusion detection rate increased. Among patients with RAS/BRAF wildtype and MSI-H, fusions were detected in 20.3% of patients. The fusion detection rate further increased to 37.5% when RNF43 mutation was added. The fusion detection rate was also higher in colon cancer than in rectal cancer. No significant differences in clinical or molecular features were found in patients with actionable gene fusions between the Genecast, cBioPortal, and ChimerSeq databases.

Conclusions

Approximately 1% of the unselected Chinese CRC population carries actionable gene fusions, mostly involving NTRK. Actionable gene fusions are more prevalent in MSI-H, RAS/BRAF wildtype, or RNF43-mutated CRC, as well as in colon cancer. Mapping of these molecular markers can markedly increase the fusion detection rate, which can help clinicians select candidates for fusion testing and targeted therapy.

Proof of Concept for Genome Profiling of the Neurofibroma/Sarcoma Sequence in Neurofibromatosis Type 1

Neurofibromatosis type 1 (NF1) is an autosomal dominant genetic disorder characterized by the predisposition to develop tumors such as malignant peripheral nerve sheath tumors (MPNSTs) which represents the primary cause of death for NF1-affected patients. Regardless of the high incidence and mortality, the molecular mechanisms underneath MPNST growth and metastatic progression remain poorly understood. In this proof-of-concept study, we performed somatic whole-exome sequencing (WES) to profile the genomic alterations in four samples from a patient with NF1-associated MPNST, consisting of a benign plexiform neurofibroma, a primary MPNST, and metastases from lung and skin tissues. By comparing genomic patterns, we identified a high level of variability across samples with distinctive genetic changes which allow for the definition of profiles of the early phase with respect to the late metastatic stages. Pathogenic and likely pathogenic variants were abundant in the primary tumor, whereas the metastatic samples exhibited a high level of copy-number variations (CNVs), highlighting a possible genomic instability in the late phases. The most known MPNST-related genes, such as TP53 and SUZ12, were identified in CNVs observed within the primary tumor. Pathway analysis of altered early genes in MPNST pointed to a potential role in cell motility, division and metabolism. Moreover, we employed survival analysis with the TCGA sarcoma genomic dataset on 262 affected patients, in order to corroborate the predictive significance of the identified early and metastatic MPNST driver genes. Specifically, the expression changes related to the mutated genes, such as in RBMX, PNPLA6 and AGAP2, were associated with reduced patient survival, distinguishing them as potential prognostic biomarkers. This study underlines the relevance of integrating genomic results with clinical information for early diagnosis and prognostic understanding of tumor aggressiveness.

Leveraging Patient-Derived Organoids for Personalized Liver Cancer Treatment

Primary liver cancer (PLC) is a primary cause of cancer-related death worldwide, and novel treatments are needed due to the limited options available for treatment and tumor heterogeneity. 66 surgically removed PLC samples were cultured using the self-developed 2:2 method, and the final success rate for organoid culture was 40.9%. Organoid performance has been evaluated using comprehensive molecular measurements, such as whole-exome and RNA sequencing, as well as anticancer drug testing. Multiple organoids and their corresponding tumor tissues contained several of the same mutations, with all pairs sharing conventional TP53 mutations. Regarding copy number variations and gene expression, significant correlations were observed between the organoids and their corresponding parental tumor tissues. Comparisons at the molecular level provided us with an assessment of organoid-to-tumor concordance, which, in combination with drug sensitivity testing provided direct guidance for treatment selection. Finally, we were able to determine an appropriate pharmacological regimen for a patient with ICC, demonstrating the clinical practicality in tailoring patient-specific drug regimens. Our study provides an organoid culture technology that can cultivate models that retain most of the molecular characteristics of tumors and can be used for drug sensitivity testing, demonstrating the broad potential application of organoid technology in precision medicine for liver cancer treatment.

Serial single-cell RNA sequencing unveils drug resistance and metastatic traits in stage IV breast cancer

Metastasis is a complex process that remains poorly understood at the molecular levels. We profiled single-cell transcriptomic, genomic, and epigenomic changes associated with cancer cell progression, chemotherapy resistance, and metastasis from a Stage IV breast cancer patient. Pretreatment- and posttreatment-specimens from the primary tumor and distant metastases were collected for single-cell RNA sequencing and subsequent cell clustering, copy number variation (CNV) estimation, transcriptomic factor estimation, and pseudotime analyses. CNV analysis revealed that a small population of pretreatment cancer cells resisted chemotherapy and expanded. New clones including Metastatic Precursor Cells (MPCs), emerged in the posttreatment primary tumors in CNV similar to metastatic cells. MPCs exhibited expression profiles indicative of epithelial-mesenchymal transition. Comparison of MPCs with metastatic cancer cells also revealed dynamic changes in transcription factors and calcitonin pathway gene expression. These findings demonstrate the utility of single-patient clinical sample analysis for understanding tumor drug resistance, regrowth, and metastasis.

Spatial single-cell protein landscape reveals vimentinhigh macrophages as immune-suppressive in the microenvironment of hepatocellular carcinoma

Tumor microenvironment heterogeneity in hepatocellular carcinoma (HCC) on a spatial single-cell resolution is unclear. Here, we conducted co-detection by indexing to profile the spatial heterogeneity of 401 HCC samples with 36 biomarkers. By parsing the spatial tumor ecosystem of liver cancer, we identified spatial patterns with distinct prognosis and genomic and molecular features, and unveiled the progressive role of vimentin (VIM)high macrophages. Integration analysis with eight independent cohorts demonstrated that the spatial co-occurrence of VIMhigh macrophages and regulatory T cells promotes tumor progression and favors immunotherapy. Functional studies further demonstrated that VIMhigh macrophages enhance the immune-suppressive activity of regulatory T cells by mechanistically increasing the secretion of interleukin-1β. Our data provide deep insights into the heterogeneity of tumor microenvironment architecture and unveil the critical role of VIMhigh macrophages during HCC progression, which holds potential for personalized cancer prevention and drug discovery and reinforces the need to resolve spatial-informed features for cancer treatment.

Cytomorphologic and molecular characterization of spindle cell carcinoid tumors of the lung

BACKGROUND

Spindle cell carcinoid tumor (SCCT) is a rare variant of lung carcinoid tumor consisting predominantly or exclusively of spindle cells. To the authors' knowledge, this is the first study to date investigating the molecular characteristics of SCCTs.

METHODS

Eighty-five carcinoid tumors initially diagnosed by fine-needle aspiration over a period of 10 years were reviewed. The final diagnostic classification was based on resection specimens. Six SCCTs were identified and characterized based on cytomorphology, and immunohistochemical and molecular features.

RESULTS

Most patients with SCCT were Caucasian (100.0%), women (83.3%), asymptomatic (66.7%), and nonsmokers (83.3%). The median age at diagnosis was 78.0 years (range, 58.2-80.3 years). A higher proportion of patients who had SCCT were diagnosed with distant metastasis. The smears were cellular and demonstrated clean backgrounds without necrosis or mitotic activity. SCCTs comprised of bipolar-to-elongated cells with finely granular chromatin, inconspicuous nucleoli, scant cytoplasm, and minimal atypia or pleomorphism. The tumor cells sometimes appeared boomerang-shaped and might mimic granulomas or blood vessels. SCCTs showed strong expression for pan-cytokeratin, synaptophysin, chromogranin, and CD56, with weak TTF-1 and a very low Ki-67 proliferation index. All SCCTs had low tumor mutational burden and were microsatellite-stable. One case showed multiple whole-gene losses in chromosome 11, whereas another harbored duplication in ARID1A. Two cases demonstrated gains in chromosomes 17, one of which also showed gains in chromosome 18. None had a single nucleotide mutation.

CONCLUSIONS

SCCT is a rare subset of lung carcinoid tumors. These tumors harbor unique cytologic, prognostic, and molecular features that may have significant diagnostic and clinical implications.

Clinicopathologic and Molecular Characterization of Inflammatory Bowel Disease-Associated Neuroendocrine Carcinomas and Mixed Neuroendocrine-Non-Neuroendocrine Neoplasms

The pathogenesis of neuroendocrine carcinomas (NECs) and mixed neuroendocrine-non-neuroendocrine neoplasms (MiNENs) in the gastrointestinal tract remains poorly understood. This study aims to characterize the clinicopathologic and molecular features of NEC/MiNEN in patients with inflammatory bowel disease (IBD). Eighteen surgically resected IBD-associated intestinal carcinomas with a minimum of 30% neuroendocrine component were collected from 6 academic centers and compared with a control group of 12 IBD-associated carcinomas lacking neuroendocrine differentiation. Both groups exhibited a male predominance and similar age distribution. The NEC/MiNEN group was more likely to have a higher percentage of Crohn disease (9/18 vs 1/12; P = .024), occur in the rectum (9/18 vs 3/12) and small intestine (4/18 vs 0/12) (P < .01), be diagnosed on resection without a preceding biopsy (6/18 vs 0/12; P = .057), and have unidentifiable precursor lesions (10/18 vs 1/12; P = .018) than the control group. Synchronous carcinoma, advanced tumor stage (pT3 and pT4), and lymph node metastasis occurred at similar rates; however, the NEC/MiNEN group had a higher incidence of angiovascular invasion (14/18 vs 4/12; P = .024), distant metastasis (8/18 vs 1/12; P = .049), mortality (8/18 vs 2/12; P = .058), and worse survival (Kaplan-Meier; P = .023) than the control group. All tested cases were mismatch repair proficient. A Ki-67 proliferation index ranged from 25% to 100%. Next-generation sequencing in 11 NEC/MiNEN cases revealed low tumor mutational burdens but complex genetic abnormalities commonly involving TP53 (9/11; 82%), FBXW7 (4/11; 36%), and APC (3/11; 27%) genes, with the other genetic alterations randomly occurring in 1 or 2 cases. The neuroendocrine component, which shared similar molecular alterations as the nonneuroendocrine component, was subcategorized into intermediate (G3a) and high grade (G3b); the higher grade correlated with more genetic alterations. In conclusion, IBD-associated NEC/MiNEN shows diverse histologic features, variable precursor lesions, intricate genetic abnormalities, and aggressive biologic behavior. The classification and grading of gastrointestinal NEC/MiNEN may be refined for better clinical management.

Aneuploidy as a driver of human cancer

Aneuploidy, an abnormal chromosome composition, is a major contributor to cancer development and progression and an important determinant of cancer therapeutic responses and clinical outcomes. Despite being recognized as a hallmark of human cancer, the exact role of aneuploidy as a 'driver' of cancer is still largely unknown. Identifying the specific genetic elements that underlie the recurrence of common aneuploidies remains a major challenge of cancer genetics. In this Review, we discuss recurrent aneuploidies and their function as drivers of tumor development. We then delve into the context-dependent identification and functional characterization of the driver genes underlying driver aneuploidies and examine emerging strategies to uncover these driver genes using cancer genomics data and cancer models. Lastly, we explore opportunities for targeting driver aneuploidies in cancer by leveraging the functional consequences of these common genetic alterations.

Characterization and integrated analysis of extrachromosomal DNA amplification in hematological malignancies

The study of extrachromosomal DNA (ecDNA), an element existing beyond classical chromosomes, contributes to creating a more comprehensive map of the cancer genome. In hematological malignancies, research on ecDNA has lacked comprehensive investigation into its frequency, structure, function, and mechanisms of formation. We re-analyzed WGS data from 208 hematological cancer samples across 11 types, focusing on ecDNA characteristics. Amplification of ecDNA was observed in 7 of these cancer types, with no instances found in normal blood cells. Patients with leukemia carrying ecDNA showed a low induction therapy remission rate (<30 %), a high relapse rate (75 %) among those who achieved complete remission, and a significantly lower survival rate compared to the general leukemia population, even those with complex chromosomal karyotypes. Among the 55 identified ecDNA amplicons, 268 genes were detected, of which 38 are known cancer-related genes exhibiting significantly increased copy numbers. By integrating RNA-Seq data, we discovered that the increased copy number, resulting in a higher amount of available DNA templates, indeed leads to the elevated expression of genes encoded on ecDNA. Additionally, through the integration of H3K4me3/H3K27ac chromatin immunoprecipitation sequencing, assay for transposase-accessible chromatin with sequencing, and high-throughput chromosome conformation capture data, we identified that ecDNA amplifications can also facilitate efficient, copy number-independent amplification of oncogenes. This process is linked to active histone modifications, improved chromatin accessibility, and enhancer hijacking, all of which are effects of ecDNA amplification. Mechanistically, chromothripsis and dysfunction of the DNA repair pathway can, to some extent, explain the origin of ecDNA.

Multiomic Characterization and Molecular Profiling of Nuclear Protein in Testis Carcinoma

PURPOSE

Nuclear protein in testis carcinoma (NC) is an underdiagnosed and aggressive squamous/poorly differentiated cancer characterized by rearrangement of the gene NUTM1 on chromosome 15q14. Co-occurring alternations have not been fully characterized.

METHODS

We analyzed the genomic and immune landscape of 54 cases of NC that underwent DNA- and RNA-based NGS sequencing (Caris).

RESULTS

While NC is driven by NUTM1 fusion oncoproteins, co-occurring DNA mutations in epigenetic or cell cycle pathways were observed in 26% of cases. There was no significant difference between the fusion partner of NUTM1 and co-occurring gene mutations. RNA sequencing analysis showed increased MYC pathway activity in NC compared with head and neck squamous cell carcinoma (HNSCC) and lung squamous cell carcinoma (LUSC), which is consistent with the known pathophysiology of NC. Characterization of the NC tumor microenvironment using RNA sequencing revealed significantly lower immune cell infiltration compared with HNSCC and LUSC. NC was 10× higher in patients with HNSCC and LUSC younger than 50 years than in those older than 70 years.

CONCLUSION

To our knowledge, this is the first series of NC profiled broadly at the DNA and RNA level. We observed fewer intratumoral immune cells by RNA sequencing, which may be associated with anecdotal data of lack of immunotherapy benefit in NC. High MYC pathway activity in NC supports ongoing trials targeting MYC suppression. The incidence of NC among patients younger than 50 years with LUSC/HNSCC supports testing for NC in these patients. The prognosis of NCs remains dismal, and future studies should focus on improving the response to immunotherapy and targeting MYC.

The estrogen response in fibroblasts promotes ovarian metastases of gastric cancer

Younger premenopausal women are more prone to developing ovarian metastases (OM) of gastric cancer (GC) than metastases of other organs; however, the molecular mechanisms remain unclear. Here we perform single-cell RNA sequencing on 45 tumor samples from 18 GC patients with OM. Interestingly, fibroblasts in OM of GC express high levels of estrogen receptor (ER) and midkine (MDK), interacting with tumor cells through activating ER-MDK-LRP1 (low-density lipoprotein receptor-related protein 1) signaling axis. Functional experiments demonstrate that estrogen stimulation induces MDK secretion by ovarian fibroblasts, and binding of MDK to LRP1 increases GC cell migration and invasion. Furthermore, in vivo, estrogen stimulation remarkably augments ovarian engraftment and metastasis of LRP1+ GC cells. Collectively, our findings reveal that ER+ ovarian fibroblasts secrete MDK under estrogen influence, driving OM of GC via the MDK-LRP1 axis. Our study holds the potential to catalyze innovative therapeutic strategies aimed at intercepting and managing OM in GC.

Clinical sequencing reveals diagnostic, therapeutic, and prognostic biomarkers for adult-type diffuse gliomas

Diffuse gliomas in adults are highly infiltrative and largely incurable. Whole exome sequencing (WES) has been demonstrated very useful in genetic analysis. Here WES was performed to characterize genomic landscape of adult-type diffuse gliomas to discover the diagnostic, therapeutic and prognostic biomarkers. Somatic and germline variants of 66 patients with adult-type diffuse gliomas were detected by WES based on the next-generation sequencing. TCGA and CGGA datasets were included to analyze the integrated diagnosis and prognosis. Among 66 patients, the diagnosis of 9 cases was changed, in which 8 cases of astrocytoma were corrected into IDH-wildtype glioblastoma (GBM), and 1 oligodendroglioma without 1p/19q co-deletion into astrocytoma. The distribution of mutations including ATRX/TP53 differed in three cohorts. The genetic mutations in GBM mainly concentrated on the cell cycle, PI3K and RTK pathways. The mutational landscape of astrocytoma was more similar to that of GBM, with the highest frequency in germline variants. Patients with IDH-mutant astrocytoma harboring SNVs of PIK3CA and PIK3R1 showed a significantly worse overall survival (OS) than wild-type patients. AEBP1 amplification was associated with shorter OS in GBM. Our study suggests that clinical sequencing can recapitulate previous findings, which may provide a powerful approach to discover diagnostic, therapeutic and prognostic markers for precision medicine in adult-type diffuse gliomas.

Age- and ethnic-driven molecular and clinical disparity of East Asian breast cancers

BACKGROUND

Breast cancer (BC) is a complex disease with profound genomic aberrations. However, the underlying molecular disparity influenced by age and ethnicity remains elusive.

METHODS

In this study, we aimed to investigate the molecular properties of 843 primary and metastatic BC patients enrolled in the K-MASTER program. By categorizing patients into two distinct age subgroups, we explored their unique molecular properties. Additionally, we leveraged large-scale genomic data from the TCGA and MSK-IMPACT studies to examine the ethnic-driven molecular and clinical disparities.

RESULTS

We observed a high prevalence of PI3KCA mutations in K-MASTER HER2 + tumors, particularly in older patients. Moreover, we identified increased mutation rates in DNA damage response molecules, including ARID1A, MSH6, and MLH1. The K-MASTER patients were mainly comprised of triple-negative breast cancer (TNBC) and HER2-positive tumors, while the TCGA and MSK-IMPACT cohorts exhibited a predominance of hormone receptor-positive (HR +) subtype tumors. Importantly, GATA3 mutations were less frequently observed in East Asian patients, which correlated with poor clinical outcomes. In addition to characterizing the molecular disparities, we developed a gradient-boosting multivariable model to identify a new molecular signature that could predict the therapeutic response to platinum-based chemotherapy.

CONCLUSIONS

Our findings collectively provide unprecedented insights into the significance of age and ethnicity on the molecular and clinical characteristics of BC patients.

Germline functional variants contribute to somatic mutation and outcomes in neuroblastoma

Germline genetic context may play a significant role in the development and evolution of cancer, particularly in childhood cancers such as neuroblastoma. This study investigates the role of putatively functional germline variants in neuroblastoma, even if they do not directly increase disease risk. Our whole-exome sequencing analysis of 125 patients with neuroblastoma reveals a positive correlation between germline variant burden and somatic mutations. Moreover, patients with higher germline variant burden exhibit worse outcomes. Similar findings are observed in the independent neuroblastoma cohort where a higher germline variant burden correlates with a higher somatic mutational burden and a worse overall survival outcome. However, contrasting results emerge in adult-onset cancer, emphasizing the importance of germline genetics in neuroblastoma. The enrichment of putatively functional germline variants in cancer predisposition genes is borderline significant when compared to healthy populations (P = 0.077; Odds Ratio, 1.45; 95% confidence intervals, 0.94-2.21) and significantly more pronounced against adult-onset cancers (P = 0.016; Odds Ratio, 2.13; 95% confidence intervals, 1.10-3.91). Additionally, the presence of these variants proves to have prognostic significance in neuroblastoma (log-rank P < 0.001), and combining germline with clinical risk factors notably improves survival predictions.

Endocrine persistence in ER+ breast cancer is accompanied by metabolic vulnerability in oxidative phosphorylation

Despite adjuvant treatment with endocrine therapies, estrogen receptor-positive (ER+) breast cancers recur in a significant proportion of patients. Recurrences are attributable to clinically undetectable endocrine-tolerant persister cancer cells that retain tumor-forming potential. Therefore, strategies targeting such persister cells may prevent recurrent disease. Using CRISPR-Cas9 genome-wide knockout screening in ER+ breast cancer cells, we identified a survival mechanism involving metabolic reprogramming with reliance upon mitochondrial respiration in endocrine-tolerant persister cells. Quantitative proteomic profiling showed reduced levels of glycolytic proteins in persisters. Metabolic tracing of glucose revealed an energy-depleted state in persisters where oxidative phosphorylation was required to generate ATP. A phase II clinical trial was conducted to evaluate changes in mitochondrial markers in primary ER+/HER2-breast tumors induced by neoadjuvant endocrine therapy ( NCT04568616 ). In an analysis of tumor specimens from 32 patients, tumors exhibiting residual cell proliferation after aromatase inhibitor-induced estrogen deprivation with letrozole showed increased mitochondrial content. Genetic profiling and barcode lineage tracing showed that endocrine-tolerant persistence occurred stochastically without genetic predisposition. Mice bearing cell line- and patient-derived xenografts were used to measure the anti-tumor effects of mitochondrial complex I inhibition in the context of endocrine therapy. Pharmacological inhibition of complex I suppressed the tumor-forming potential of persisters and synergized with the anti-estrogen fulvestrant to induce regression of patient-derived xenografts. These findings indicate that mitochondrial metabolism is essential in endocrine-tolerant persister ER+ breast cancer cells and warrant the development of treatment strategies to leverage this vulnerability in the context of endocrine-sensitive disease.

Statement of Significance

Endocrine-tolerant persister cancer cells that survive endocrine therapy can cause recurrent disease. Persister cells exhibit increased energetic dependence upon mitochondria for survival and tumor re-growth potential.

Treatment resistance to melanoma therapeutics on a single cell level

Therapy targeting the BRAF-MEK cascade created a treatment revolution for patients with BRAF mutant advanced melanoma. Unfortunately, 80% patients treated will progress by 5 years follow-up. Thus, it is imperative we study mechanisms of melanoma progression and therapeutic resistance. We created a scRNA (single cell RNA) atlas of 128,230 cells from 18 tumors across the treatment spectrum, discovering melanoma cells clustered strongly by transcriptome profiles of patients of origins. Our cell-level investigation revealed gains of 1q and 7q as likely early clonal events in metastatic melanomas. By comparing patient tumors and their derivative cell lines, we observed that PD1 responsive tumor fraction disappears when cells are propagated in vitro. We further established three anti-BRAF-MEK treatment resistant cell lines using three BRAF mutant tumors. ALDOA and PGK1 were found to be highly expressed in treatment resistant cell populations and metformin was effective in targeting the resistant cells. Our study suggests that the investigation of patient tumors and their derivative lines is essential for understanding disease progression, treatment response and resistance.

Clinical characteristics and genomic profiling of outpatients with endometrial cancer at a Chinese tertiary cancer center

OBJECTIVE

Endometrial cancer is stepping into the era of precision therapy. Genomic test is recommended for newly diagnostic patients. However, outpatients displayed more complex characteristics. Here, we elucidated the clinical characteristics and genomic profiling of outpatients with endometrial cancer at our institution.

METHODS

Between 2018 and 2023, 68 endometrial cancer received genomic tests at outpatient department of Fudan University Shanghai Cancer Center. Data, including age, pathological histology, FIGO stage and treatment strategy were collected. Germline mutations, molecular subtypes and other somatic mutations were also summarized.

RESULTS

Overall, 72.1% (49/68) of patients receive genomic tests at primary diagnosis, while 27.9% (19/68) of patients received tests at recurrence. Nine patients had deleterious germline mutations, including BRCA1(2), MLH1(1), MSH2(2, including one with co-mutation of RAD50), MSH6(2), FANCA(1), MUTYH(1). Molecular subtypes were recognized among 62 patients, as POLE super-mutation(4, 6.5%), MSI-H(7, 11.3%), CN-Low(36, 58.1%) and CN-High(15, 24.2%). Ten patients received anti-PD1 monotherapy or in combination with chemotherapy or anti-angiogenic therapy, with the duration of disease control of 1 to 35 months. The ORR rate was 30%, and six patients had stable disease. The median (range) follow-up time was 18(2-160) months. 23(33.8%) relapses were recorded, and CN-High subtype displayed worst PFS compared with other subtypes (P < 0.01). 6 deaths were reported including 2(5.6%) of CN-Low and 4(26.7%) of CN-High.

CONCLUSION

Outpatients department gathered a considerable proportion of recurrent patients with complex genomic features. Patients with worse prognosis could be well studied, and anti-PD1 therapy was a promising salvage therapy in the real world.

Mechanistic insights into lethal hyper progressive disease induced by PD-L1 inhibitor in metastatic urothelial carcinoma

Hyper progressive disease (HPD) is a paradoxical phenomenon characterized by accelerated tumor growth following treatment with immune checkpoint inhibitors. However, the pathogenic causality and its predictor remain unknown. We herein report a fatal case of HPD in a 50-year-old man with metastatic bladder cancer. He had achieved a complete response (CR) through chemoradiation therapy followed by twelve cycles of chemotherapy, maintaining CR for 24 months. Three weeks after initiating maintenance use of a PD-L1 inhibitor, avelumab, a massive amount of metastases developed, leading to the patient's demise. Omics analysis, utilizing metastatic tissues obtained from an immediate autopsy, implied the contribution of M2 macrophages, TGF-β signaling, and interleukin-8 to HPD pathogenesis.

Novel ERBB2 Variant Potentially Associated with Resistance against Anti-HER2 Monoclonal Antibody-Based Therapy in ERBB2-Amplified Metastatic Colorectal Cancer

PURPOSE

HER2-targeted therapies in ERBB2-amplified metastatic colorectal cancer (mCRC) are effective; however, a notable portion of patients do not respond to treatment, and secondary resistance occurs in most patients receiving these treatments. The purpose of this study was to investigate determinants of treatment efficacy and resistance in patients with ERBB2-amplified mCRC who received HER2-targeted therapy by analyzing multiomics data.

EXPERIMENTAL DESIGN

We investigated genomic data from a nationwide large cancer genomic screening project, the SCRUM-Japan project. We analyzed paired genome and transcriptome data of tissue and genomic data of ctDNA collected pre- and postprogression in patients enrolled in the related trial, TRIUMPH, in ERBB2-amplified mCRC.

RESULTS

In 155 patients with ERBB2-amplified solid tumors who received HER2-targeted therapy based on the SCRUM-Japan project, the objective response rate was 50%, 51%, and 35% in ERBB2 wild-type, variant of unknown significance, and pathogenic variant groups, respectively. In the paired genome and transcriptome data analyses in TRIUMPH, we identified the novel splicing-associated variant c.644-66_-2del in one of the 11 patients with paired whole-exome sequencing and whole-transcriptome sequencing data sets, which lacks the binding domain of pertuzumab, in progressed metastatic tumor as a variant with potential pathogenicity. The time-course ctDNA analysis detected c.644-66_-2del as an acquired variant.

CONCLUSIONS

This study highlighted the importance of ERBB2 genomic status when evaluating the efficacy of HER2-targeted therapies in ERBB2-amplified mCRC. The identification of a novel splicing-associated variant may provide insights into potential mechanisms of treatment resistance. Furthermore, we demonstrated the utility of ctDNA to follow the acquired genomic status of mCRC tumors.

SJPedPanel: A Pan-Cancer Gene Panel for Childhood Malignancies to Enhance Cancer Monitoring and Early Detection

PURPOSE

The purpose of the study was to design a pan-cancer gene panel for childhood malignancies and validate it using clinically characterized patient samples.

EXPERIMENTAL DESIGN

In addition to 5,275 coding exons, SJPedPanel also covers 297 introns for fusions/structural variations and 7,590 polymorphic sites for copy-number alterations. Capture uniformity and limit of detection are determined by targeted sequencing of cell lines using dilution experiment. We validate its coverage by in silico analysis of an established real-time clinical genomics (RTCG) cohort of 253 patients. We further validate its performance by targeted resequencing of 113 patient samples from the RTCG cohort. We demonstrate its power in analyzing low tumor burden specimens using morphologic remission and monitoring samples.

RESULTS

Among the 485 pathogenic variants reported in RTCG cohort, SJPedPanel covered 86% of variants, including 82% of 90 rearrangements responsible for fusion oncoproteins. In our targeted resequencing cohort, 91% of 389 pathogenic variants are detected. The gene panel enabled us to detect ∼95% of variants at allele fraction (AF) 0.5%, whereas the detection rate is ∼80% at AF 0.2%. The panel detected low-frequency driver alterations from morphologic leukemia remission samples and relapse-enriched alterations from monitoring samples, demonstrating its power for cancer monitoring and early detection.

CONCLUSIONS

SJPedPanel enables the cost-effective detection of clinically relevant genetic alterations including rearrangements responsible for subtype-defining fusions by targeted sequencing of ∼0.15% of human genome for childhood malignancies. It will enhance the analysis of specimens with low tumor burdens for cancer monitoring and early detection.

Mutational spectrum of breast cancer by shallow whole-genome sequencing of cfDNA and tumor gene panel analysis

Breast cancer (BC) has different molecular subgroups related to different risks and treatments. Tumor biopsies for BC detection are invasive and may not reflect tumor heterogeneity. Liquid biopsies have become relevant because they might overcome these limitations. We rationalize that liquid cfDNA biopsies through shallow whole genome sequencing (sWGS) could improve the detection of tumor alterations, complementing the genomic profiling. We evaluated the feasibility to detect somatic copy number alterations (SCNAs) in BC using shallow whole genome sequencing (sWGS) in cfDNA from archived samples from National Cancer Institute of Colombia patients. We sequenced tumor tissues from 38 BC patients with different molecular subtypes using a gene panel of 176 genes significantly mutated in cancer, and by liquid biopsies using sWGS on 20 paired samples to detect SCNAs and compare with the tumor samples. We identified an extensive intertumoral heterogeneity between the molecular subtypes of BC, with a mean tumor load of 602 mutations in the gene panel of tumor tissues. There was a 12.3% of concordance in deletions in the cfDNA-tumor pairs considering only the genes covered by the panel encompassing seven genes: BRCA1, CDK12, NF1, MAP2K4, NCOR1, TP53, and KEAP1 in three patients. This study shows the feasibility to complement the genomic analysis of tumor tissue biopsies to detect SCNA in BC using sWGS in cfDNA, providing a wider identification of potential therapeutic targets.

Functional pathogenicity of ESRRB variant of uncertain significance contributes to hearing loss (DFNB35)

Advances in next-generation sequencing technologies have led to elucidation of sensorineural hearing loss genetics and associated clinical impacts. However, studies on the functional pathogenicity of variants of uncertain significance (VUS), despite their close association with clinical phenotypes, are lacking. Here we identified compound heterozygous variants in ESRRB transcription factor gene linked to DFNB35, specifically a novel splicing variant (NM_004452.4(ESRRB): c.397 + 2T>G) in trans with a missense variant (NM_004452.4(ESRRB): c.1144C>T p.(Arg382Cys)) whose pathogenicity remains unclear. The splicing variant (c.397 + 2T>G) caused exon 4 skipping, leading to premature stop codon formation and nonsense-mediated decay. The p.(Arg382Cys) variant was classified as a VUS due to its particularly higher allele frequency among East Asian population despite disease-causing in-silico predictions. However, functional assays showed that p.(Arg382Cys) variant disrupted key intramolecular interactions, leading to protein instability. This variant also reduced transcriptional activity and altered expression of downstream target genes essential for inner ear function, suggesting genetic contribution to disease phenotype. This study expanded the phenotypic and genotypic spectrum of ESRRB in DFNB35 and revealed molecular mechanisms underlying ESRRB-associated DFNB35. These findings suggest that variants with high allele frequencies can also possess functional pathogenicity, providing a breakthrough for cases where VUS, previously unexplored, could be reinterpreted by elucidating their functional roles and disease-causing characteristics.

Prevalence of germline variants in Brazilian pancreatic carcinoma patients

We evaluated the prevalence of pathogenic/likely pathogenic germline variants (PGV) in Brazilian pancreatic adenocarcinoma (PC) patients, that represent a multiethnic population, in a cross-sectional study. We included 192 PC patients unselected for family history of cancer. We evaluated a panel of 113 cancer genes, through genomic DNA sequencing and 46 ancestry-informative markers, through multiplex PCR. The median age was 61 years; 63.5% of the patients presented disease clinical stages III or IV; 8.3% reported personal history of cancer; 4.7% and 16.1% reported first-degree relatives with PC or breast and/or prostate cancer, respectively. Although the main ancestry was European, there was considerable genetic composition admixture. Twelve patients (6.25%) were PGV carriers in PC predisposition genes (ATM, BRCA1, BRCA2, CDKN2A, MSH2, PALB2) and another 25 (13.0%) were PGV carriers in genes with a limited association or not previously associated with PC (ACD, BLM, BRIP1, CHEK2, ERCC4, FANCA, FANCE, FANCM, GALNT12, MITF, MRE11, MUTYH, POLE, RAD51B, RAD51C, RECQL4, SDHA, TERF2IP). The most frequently affected genes were CHEK2, ATM and FANC. In tumor samples from PGV carriers in ACD, BRIP1, MRE11, POLE, SDHA, TERF2IP, which were examined through exome sequencing, the main single base substitutions (SBS) mutational signature was SBS1+5+18, probably associated with age, tobacco smoking and reactive oxygen species. SBS3 associated with homologous repair deficiency was also represented, but on a lower scale. There was no difference in the frequency of PGV carriers between: (a) patients with or without first-degree relatives with cancer; and (b) patients with admixed ancestry versus those with predominantly European ancestry. Furthermore, there was no difference in overall survival between PGV carriers and non-carriers. Therefore, genetic testing should be offered to all Brazilian pancreatic cancer patients, regardless of their ancestry. Genes with limited or previously unrecognized associations with pancreatic cancer should be further investigated to clarify their role in cancer risk.

Human aneuploid cells depend on the RAF/MEK/ERK pathway for overcoming increased DNA damage

Aneuploidy is a hallmark of human cancer, yet the molecular mechanisms to cope with aneuploidy-induced cellular stresses remain largely unknown. Here, we induce chromosome mis-segregation in non-transformed RPE1-hTERT cells and derive multiple stable clones with various degrees of aneuploidy. We perform a systematic genomic, transcriptomic and proteomic profiling of 6 isogenic clones, using whole-exome DNA, mRNA and miRNA sequencing, as well as proteomics. Concomitantly, we functionally interrogate their cellular vulnerabilities, using genome-wide CRISPR/Cas9 and large-scale drug screens. Aneuploid clones activate the DNA damage response and are more resistant to further DNA damage induction. Aneuploid cells also exhibit elevated RAF/MEK/ERK pathway activity and are more sensitive to clinically-relevant drugs targeting this pathway, and in particular to CRAF inhibition. Importantly, CRAF and MEK inhibition sensitize aneuploid cells to DNA damage-inducing chemotherapies and to PARP inhibitors. We validate these results in human cancer cell lines. Moreover, resistance of cancer patients to olaparib is associated with high levels of RAF/MEK/ERK signaling, specifically in highly-aneuploid tumors. Overall, our study provides a comprehensive resource for genetically-matched karyotypically-stable cells of various aneuploidy states, and reveals a therapeutically-relevant cellular dependency of aneuploid cells.