Control-FREEC: a tool for assessing copy number and allelic content using next-generation sequencing data

Sequential antigen loss and branching evolution in lymphoma after CD19- and CD20-targeted T-cell-redirecting therapy

ABSTRACT

CD19 chimeric antigen receptor (CAR) T cells and CD20 targeting T-cell-engaging bispecific antibodies (bispecs) have been approved in B-cell non-Hodgkin lymphoma lately, heralding a new clinical setting in which patients are treated with both approaches, sequentially. The aim of our study was to investigate the selective pressure of CD19- and CD20-directed therapy on the clonal architecture in lymphoma. Using a broad analytical pipeline on 28 longitudinally collected specimen from 7 patients, we identified truncating mutations in the gene encoding CD20 conferring antigen loss in 80% of patients relapsing from CD20 bispecs. Pronounced T-cell exhaustion was identified in cases with progressive disease and retained CD20 expression. We also confirmed CD19 loss after CAR T-cell therapy and reported the case of sequential CD19 and CD20 loss. We observed branching evolution with re-emergence of CD20+ subclones at later time points and spatial heterogeneity for CD20 expression in response to targeted therapy. Our results highlight immunotherapy as not only an evolutionary bottleneck selecting for antigen loss variants but also complex evolutionary pathways underlying disease progression from these novel therapies.

Control-FREEC viewer: a tool for the visualization and exploration of copy number variation data

BACKGROUND

Copy number alterations (CNAs) are genetic changes commonly found in cancer that involve different regions of the genome and impact cancer progression by affecting gene expression and genomic stability. Computational techniques can analyze copy number data obtained from high-throughput sequencing platforms, and various tools visualize and analyze CNAs in cancer genomes, providing insights into genetic mechanisms driving cancer development and progression. However, tools for visualizing copy number data in cancer research have some limitations. In fact, they can be complex to use and require expertise in bioinformatics or computational biology. While copy number data analysis and visualization provide insights into cancer biology, interpreting results can be challenging, and there may be multiple explanations for observed patterns of copy number alterations.

RESULTS

We created Control-FREEC Viewer, a tool that facilitates effective visualization and exploration of copy number data. With Control-FREEC Viewer, experimental data can be easily loaded by the user. After choosing the reference genome, copy number data are displayed in whole genome or single chromosome view. Gain or loss on a specific gene can be found and visualized on each chromosome. Analysis parameters for subsequent sessions can be stored and images can be exported in raster and vector formats.

CONCLUSIONS

Control-FREEC Viewer enables users to import and visualize data analyzed by the Control-FREEC tool, as well as by other tools sharing a similar tabular output, providing a comprehensive and intuitive graphical user interface for data visualization.

The C250T Mutation of TERTp Might Grant a Better Prognosis to Glioblastoma by Exerting Less Biological Effect on Telomeres and Chromosomes Than the C228T Mutation

The aim of this study was to determine how TERTp mutations impact glioblastoma prognosis.

MATERIALS AND METHODS

TERTp mutations were assessed in a retrospective cohort of 258 uniformly treated glioblastoma patients. RNA-sequencing and whole exome sequencing results were available in a subset of patients.

RESULTS

Overall, there were no differences in outcomes between patients with mutated TERTp-wt or TERTp. However, we found significant differences according to the type of TERTp mutation. Progression-free survival (mPFS) was 9.1 months for those with the C250T mutation and 7 months for those with either the C228T mutation or TERTp-wt (p = 0.016). Overall survival (mOS) was 21.9 and 15 months, respectively (p = 0.026). This differential effect was more pronounced in patients with MGMTp methylation (mPFS: p = 0.008; mOS: p = 0.021). Multivariate analysis identified the C250T mutation as an independent prognostic factor for longer mOS (HR 0.69; p = 0.044). We found no differences according to TERTp mutation status in molecular alterations common in glioblastoma, nor in copy number variants in genes related to alternative lengthening of telomeres. Nevertheless, in the gene enrichment analysis adjusted for MGMTp methylation status, some Reactome gene sets were differentially enriched, suggesting that the C250T mutation may exert a lesser effect on telomeres or chromosomes.

CONCLUSIONS

In our series, patients exhibiting the C250T mutation had a more favorable prognosis compared to those with either TERPp-wt or TERTp C228T mutations. Additionally, our findings suggest a reduced involvement of the C250T mutation in the underlying biological mechanisms related to telomeres.

Comprehensive Genomic Analysis of Cemento-Ossifying Fibroma

Cemento-ossifying fibroma (COF) of the jaws is currently classified as a benign mesenchymal odontogenic tumor, and only targeted approaches have been used to assess its genetic alterations. A minimal proportion of COFs harbor CDC73 somatic mutations, and copy number alterations (CNAs) involving chromosomes 7 and 12 have recently been reported in a small proportion of cases. However, the genetic background of COFs remains obscure. We used a combination of whole-exome sequencing and RNA sequencing to assess somatic mutations, fusion transcripts, and CNAs in a cohort of 12 freshly collected COFs. No recurrent fusions have been identified among the 5 cases successfully analyzed by RNA sequencing, with in-frame fusions being detected in 2 cases (MARS1::GOLT1B and PARG::BMS1 in one case and NCLN::FZR1 and NFIC::SAMD1 in the other case) and no candidate fusions identified for the remaining 3 cases. No recurrent pathogenic mutations were detected in the 11 cases that had undergone whole-exome sequencing. A KRAS p.L19F missense variant was detected in one case, and 2 CDC73 deletions were detected in another case. The other variants were of uncertain significance and included variants in PC, ACTB, DOK6, HACE1, and COL1A2 and previously unreported variants in PTPN14, ATP5F1C, APOBEC1, HDAC5, ATF7IP, PARP2, and ACTR3B. The affected genes do not clearly converge on any signaling pathway. CNAs were detected in 5/11 cases (45%), with copy gains involving chromosome 12 occurring in 3/11 cases (27%). In conclusion, no recurrent fusions or pathogenic variants have been detected in the present COF cohort, with copy gains involving chromosome 12 occurring in 27% of cases.

OTSUCNV: an adaptive segmentation and OTSU-based anomaly classification method for CNV detection using NGS data

Copy-number variations (CNVs), which refer to deletions and duplications of chromosomal segments, represent a significant source of variation among individuals, contributing to human evolution and being implicated in various diseases ranging from mental illness and developmental disorders to cancer. Despite the development of several methods for detecting copy number variations based on next-generation sequencing (NGS) data, achieving robust detection performance for CNVs with arbitrary coverage and amplitude remains challenging due to the inherent complexity of sequencing samples. In this paper, we propose an alternative method called OTSUCNV for CNV detection on whole genome sequencing (WGS) data. This method utilizes a newly designed adaptive sequence segmentation algorithm and an OTSU-based CNV prediction algorithm, which does not rely on any distribution assumptions or involve complex outlier factor calculations. As a result, the effective detection of CNVs is achieved with lower computational complexity. The experimental results indicate that the proposed method demonstrates outstanding performance, and hence it may be used as an effective tool for CNV detection.

On the core segmentation algorithms of copy number variation detection tools

Shotgun sequencing is a high-throughput method used to detect copy number variants (CNVs). Although there are numerous CNV detection tools based on shotgun sequencing, their quality varies significantly, leading to performance discrepancies. Therefore, we conducted a comprehensive analysis of next-generation sequencing-based CNV detection tools over the past decade. Our findings revealed that the majority of mainstream tools employ similar detection rationale: calculates the so-called read depth signal from aligned sequencing reads and then segments the signal by utilizing either circular binary segmentation (CBS) or hidden Markov model (HMM). Hence, we compared the performance of those two core segmentation algorithms in CNV detection, considering varying sequencing depths, segment lengths and complex types of CNVs. To ensure a fair comparison, we designed a parametrical model using mainstream statistical distributions, which allows for pre-excluding bias correction such as guanine-cytosine (GC) content during the preprocessing step. The results indicate the following key points: (1) Under ideal conditions, CBS demonstrates high precision, while HMM exhibits a high recall rate. (2) For practical conditions, HMM is advantageous at lower sequencing depths, while CBS is more competitive in detecting small variant segments compared to HMM. (3) In case involving complex CNVs resembling real sequencing, HMM demonstrates more robustness compared with CBS. (4) When facing large-scale sequencing data, HMM costs less time compared with the CBS, while their memory usage is approximately equal. This can provide an important guidance and reference for researchers to develop new tools for CNV detection.

labelSeg: segment annotation for tumor copy number alteration profiles

Somatic copy number alterations (SCNAs) are a predominant type of oncogenomic alterations that affect a large proportion of the genome in the majority of cancer samples. Current technologies allow high-throughput measurement of such copy number aberrations, generating results consisting of frequently large sets of SCNA segments. However, the automated annotation and integration of such data are particularly challenging because the measured signals reflect biased, relative copy number ratios. In this study, we introduce labelSeg, an algorithm designed for rapid and accurate annotation of CNA segments, with the aim of enhancing the interpretation of tumor SCNA profiles. Leveraging density-based clustering and exploiting the length-amplitude relationships of SCNA, our algorithm proficiently identifies distinct relative copy number states from individual segment profiles. Its compatibility with most CNA measurement platforms makes it suitable for large-scale integrative data analysis. We confirmed its performance on both simulated and sample-derived data from The Cancer Genome Atlas reference dataset, and we demonstrated its utility in integrating heterogeneous segment profiles from different data sources and measurement platforms. Our comparative and integrative analysis revealed common SCNA patterns in cancer and protein-coding genes with a strong correlation between SCNA and messenger RNA expression, promoting the investigation into the role of SCNA in cancer development.

The origins of the Guinness stout yeast

Beer is made via the fermentation of an aqueous extract predominantly composed of malted barley flavoured with hops. The transforming microorganism is typically a single strain of Saccharomyces cerevisiae, and for the majority of major beer brands the yeast strain is a unique component. The present yeast used to make Guinness stout brewed in Dublin, Ireland, can be traced back to 1903, but its origins are unknown. To that end, we used Illumina and Nanopore sequencing to generate whole-genome sequencing data for a total of 22 S. cerevisiae yeast strains: 16 from the Guinness collection and 6 other historical Irish brewing. The origins of the Guinness yeast were determined with a SNP-based analysis, demonstrating that the Guinness strains occupy a distinct group separate from other historical Irish brewing yeasts. Assessment of chromosome number, copy number variation and phenotypic evaluation of key brewing attributes established Guinness yeast-specific SNPs but no specific chromosomal amplifications. Our analysis also demonstrated the effects of yeast storage on phylogeny. Altogether, our results suggest that the Guinness yeast used today is related to the first deposited Guinness yeast; the 1903 Watling Laboratory Guinness yeast.

BARD1 germline variants induce haploinsufficiency and DNA repair defects in neuroblastoma

BACKGROUND

High-risk neuroblastoma is a complex genetic disease that is lethal in more than 50% of patients despite intense multimodal therapy. Through genome-wide association studies (GWAS) and next-generation sequencing, we have identified common single nucleotide polymorphisms and rare, pathogenic or likely pathogenic germline loss-of-function variants in BARD1 enriched in neuroblastoma patients. The functional implications of these findings remain poorly understood.

METHODS

We correlated BARD1 genotype with expression in normal tissues and neuroblastomas, along with the burden of DNA damage in tumors. To validate the functional consequences of germline pathogenic or likely pathogenic BARD1 variants, we used CRISPR-Cas9 to generate isogenic neuroblastoma (IMR-5) and control (RPE1) cellular models harboring heterozygous BARD1 loss-of-function variants (R112*, R150*, E287fs, and Q564*) and quantified genomic instability in these cells via next-generation sequencing and with functional assays measuring the efficiency of DNA repair.

RESULTS

Both common and rare neuroblastoma-associated BARD1 germline variants were associated with lower levels of BARD1 mRNA and an increased burden of DNA damage. Using isogenic heterozygous BARD1 loss-of-function variant cellular models, we functionally validated this association with inefficient DNA repair. BARD1 loss-of-function variant isogenic cells exhibited reduced efficiency in repairing Cas9-induced DNA damage, ineffective RAD51 focus formation at DNA double-strand break sites, and enhanced sensitivity to cisplatin and poly (ADP-ribose) polymerase (PARP) inhibition both in vitro and in vivo.

CONCLUSIONS

Taken together, we demonstrate that germline BARD1 variants disrupt DNA repair fidelity. This is a fundamental molecular mechanism contributing to neuroblastoma initiation that may have important therapeutic implications.

Classification of Brain Tumors by Nanopore Sequencing of Cell-Free DNA from Cerebrospinal Fluid

BACKGROUND

Molecular brain tumor diagnosis is usually dependent on tissue biopsies or resections. This can pose several risks associated with anesthesia or neurosurgery, especially for lesions in the brain stem or other difficult-to-reach anatomical sites. Apart from initial diagnosis, tumor progression, recurrence, or the acquisition of novel genetic alterations can only be proven by re-biopsies.

METHODS

We employed Nanopore sequencing on cell-free DNA (cfDNA) from cerebrospinal fluid (CSF) and analyzed copy number variations (CNV) and global DNA methylation using a random forest classifier. We sequenced 129 samples with sufficient DNA. These samples came from 99 patients and encompassed 22 entities. Results were compared to clinical diagnosis and molecular analysis of tumor tissue, if available.

RESULTS

110/129 samples were technically successful, and 50 of these contained detectable circulating tumor DNA (ctDNA) by CNV or methylation profiling. ctDNA was detected in samples from patients with progressive disease but also from patients without known residual disease. CNV plots showed diagnostic and prognostic alterations, such as C19MC amplifications in embryonal tumors with multilayered rosettes or Chr.1q gains and Chr.6q losses in posterior fossa group A ependymoma, respectively. Most CNV profiles mirrored the profiles of the respective tumor tissue. DNA methylation allowed exact classification of the tumor in 22/110 cases and led to incorrect classification in 2/110 cases. Only 5/50 samples with detected ctDNA contained tumor cells detectable through microscopy.

CONCLUSIONS

Our results suggest that Nanopore sequencing data of cfDNA from CSF samples may be a promising approach for initial brain tumor diagnostics and an important tool for disease monitoring.

ECOLE: Learning to call copy number variants on whole exome sequencing data

Copy number variants (CNV) are shown to contribute to the etiology of several genetic disorders. Accurate detection of CNVs on whole exome sequencing (WES) data has been a long sought-after goal for use in clinics. This was not possible despite recent improvements in performance because algorithms mostly suffer from low precision and even lower recall on expert-curated gold standard call sets. Here, we present a deep learning-based somatic and germline CNV caller for WES data, named ECOLE. Based on a variant of the transformer architecture, the model learns to call CNVs per exon, using high-confidence calls made on matched WGS samples. We further train and fine-tune the model with a small set of expert calls via transfer learning. We show that ECOLE achieves high performance on human expert labelled data for the first time with 68.7% precision and 49.6% recall. This corresponds to precision and recall improvements of 18.7% and 30.8% over the next best-performing methods, respectively. We also show that the same fine-tuning strategy using tumor samples enables ECOLE to detect RT-qPCR-validated variations in bladder cancer samples without the need for a control sample. ECOLE is available at https://github.com/ciceklab/ECOLE .

Characterization of dabrafenib-induced drug insensitivity via cellular barcoding and collateral sensitivity to second-line therapeutics

Drug insensitivity is arguably one of the biggest challenges in cancer therapeutics. Although effective therapeutic solutions in cancer are limited due to the emergence of drug insensitivity, exploiting evolutionary understanding in this context can provide potential second-line therapeutics sensitizing the drug insensitive populations. Targeted therapeutic agent dabrafenib is used to treat CRC patients with BRAF V600E genotype and insensitivity to dabrafenib is often observed. Understanding underlying clonal architecture of dabrafenib-induced drug insensitivity and identification of potential second-line therapeutics that could sensitize dabrafenib insensitive populations remain to be elucidated. For this purpose, we utilized cellular barcoding technology to decipher dabrafenib-induced clonal evolution in BRAF V600E mutant HT-29 cells. This approach revealed the detection of both pre-existing and de novo barcodes with increased frequencies as a result of dabrafenib insensitivity. Furthermore, our longitudinal monitoring of drug insensitivity based on barcode detection from floating DNA within used medium enabled to identify temporal dynamics of pre-existing and de novo barcodes in relation to dabrafenib insensitivity in HT-29 cells. Moreover, whole-exome sequencing analysis exhibited possible somatic CNVs and SNVs contributing to dabrafenib insensitivity in HT-29 cells. Last, collateral drug sensitivity testing demonstrated oxaliplatin and capecitabine, alone or in combination, as successful second-like therapeutics in inducing collateral sensitivity in dabrafenib-insensitive HT-29 cells. Overall, our findings demonstrate clonal dynamics of dabrafenib-insensitivity in HT-29 cells. In addition, oxaliplatin and capecitabine, alone or in combination, were successful second-line therapeutics in inducing collateral sensitivity in dabrafenib-insensitive HT-29 cells.

Whole‑exome sequencing reveals Lewis lung carcinoma is a hypermutated Kras/Nras-mutant cancer with extensive regional mutation clusters in its genome

Lewis lung carcinoma (LLC), as a widely used preclinical cancer model, has still not been genetically and genomically characterized. Here, we performed a whole-exome sequencing analysis on the LLC cell line to elucidate its molecular characteristics and etiologies. Our data showed that LLC originated from a male mouse belonging to C57BL/6L (a transitional strain between C57BL/6J and C57BL/6N) and contains substantial somatic SNV and InDel mutations (> 20,000). Extensive regional mutation clusters are present in its genome, which were caused mainly by the mutational processes underlying the SBS1, SBS5, SBS15, SBS17a, and SBS21 signatures during frequent structural rearrangements. Thirty three deleterious mutations are present in 30 cancer genes including Kras, Nras, Trp53, Dcc, and Cacna1d. Cdkn2a and Cdkn2b are biallelically deleted from the genome. Five pathways (RTK/RAS, p53, cell cycle, TGFB, and Hippo) are oncogenically deregulated or affected. The major mutational processes in LLC include chromosomal instability, exposure to metabolic mutagens, spontaneous 5-methylcytosine deamination, defective DNA mismatch repair, and reactive oxygen species. Our data also suggest that LLC is a lung cancer similar to human lung adenocarcinoma. This study lays a molecular basis for the more targeted application of LLC in preclinical research.

Ovarian Cancer Patient-Derived Organoids Used as a Model for Replicating Genetic Characteristics and Testing Drug Responsiveness: A Preliminary Study

This study aimed to explore the role of ovarian cancer patient-derived organoids (PDOs) in their replicating genetic characteristics and testing drug responsiveness. Ovarian cancer PDOs were cultured in Matrigel with a specialized medium. The successful rate and proliferation rate were calculated. Morphology, histology, and immunohistochemistry (IHC) (PAX8, P53, and WT1) were used to identify the tumor characteristics. Gene sequencing, variant allele frequency (VAF), and copy number variation were used to explore the mutation profile. The sensitivity to chemodrugs (carboplatin, paclitaxel, gemcitabine, doxorubicin, and olaparib) was conducted. Successful generation of organoids occurred in 54% (7/13) of attempts, encompassing 4 high-grade serous carcinomas (HGSC), 1 mucinous carcinoma (MC), 1 clear cell carcinoma (CCC), and 1 carcinosarcoma. The experiments used six organoids (3 HGSC, 1 CCC, 1 MC, and 1 carcinosarcoma). The derived organoids exhibited spherical-like morphology, and the diameter ranged from 100 to 500 μm. The histology and IHC exhibited the same between organoids and primary tumors. After cryopreservation, the organoid's growth rate was slower than the primary culture (14 days vs 10 days, P < 0.01). Targeted sequencing revealed shared DNA variants, including mutations in key genes, such as BRCA1, PIK3CA, ARID1A, and TP53. VAF was similar between primary tumors and organoids. The organoids maintained inherited most copy number alterations. Drug sensitivity testing revealed varying responses, with carcinosarcoma organoids showing higher sensitivity to paclitaxel and gemcitabine than HGSC organoids. Our preliminary results showed that ovarian cancer PDOs could be successfully derived and histology, mutations, and diverse copy numbers of genotypes could be faithfully captured. Drug testing could reveal the individual PDO's responsiveness to drugs. PDOs might be as valuable resources for investigating genomic biomarkers for personalized treatment.

Downregulation of the expression of subgenomic chromosome A7 genes promotes plant height in resynthesized allopolyploid Brassica napus

KEY MESSAGE

Homoeolog expression bias and the gene dosage effect induce downregulation of genes on chromosome A7, causing a significant increase in the plant height of resynthesized allopolyploid Brassica napus. Gene expression levels in allopolyploid plants are not equivalent to the simple average of the expression levels in the parents and are associated with several non-additive expression phenomena, including homoeolog expression bias. However, hardly any information is available on the effect of homoeolog expression bias on traits. Here, we studied the effects of gene expression-related characteristics on agronomic traits using six isogenic resynthesized Brassica napus lines across the first ten generations. We found a group of genes located on chromosome A7 whose expression levels were significantly negatively correlated with plant height. They were expressed at significantly lower levels than their homoeologous genes, owing to allopolyploidy rather than inheritance from parents. Homoeolog expression bias resulted in resynthesized allopolyploids with a plant height similar to their female Brassica oleracea parent, but significantly higher than that of the male Brassica rapa parent. Notably, aneuploid lines carrying monosomic and trisomic chromosome A7 had the highest and lowest plant heights, respectively, due to changes in the expression bias of homoeologous genes because of alterations in the gene dosage. These findings suggest that the downregulation of the expression of homoeologous genes on a single chromosome can result in the partial improvement of traits to a significant extent in the nascent allopolyploid B. napus.

Whole-genome sequencing reveals the molecular implications of the stepwise progression of lung adenocarcinoma

The mechanism underlying the development of tumors, particularly at early stages, still remains mostly elusive. Here, we report whole-genome long and short read sequencing analysis of 76 lung cancers, focusing on very early-stage lung adenocarcinomas such as adenocarcinoma in situ (AIS) and minimally invasive adenocarcinoma. The obtained data is further integrated with bulk and spatial transcriptomic data and epigenomic data. These analyses reveal key events in lung carcinogenesis. Minimal somatic mutations in pivotal driver mutations and essential proliferative factors are the only detectable somatic mutations in the very early-stage of AIS. These initial events are followed by copy number changes and global DNA hypomethylation. Particularly, drastic changes are initiated at the later AIS stage, i.e., in Noguchi type B tumors, wherein cancer cells are exposed to the surrounding microenvironment. This study sheds light on the pathogenesis of lung adenocarcinoma from integrated pathological and molecular viewpoints.

Genomic determinants of response and resistance to inotuzumab ozogamicin in B-cell ALL

Inotuzumab ozogamicin (InO) is an antibody-drug conjugate that delivers calicheamicin to CD22-expressing cells. In a retrospective cohort of InO treated patients with B-cell acute lymphoblastic leukemia, we sought to understand the genomic determinants of response to InO. Acquired CD22 mutations were observed in 11% (3/27) of post-InO relapsed tumor samples. There were multiple CD22 mutations per sample and the mechanisms of CD22 escape included protein truncation, protein destabilization, and epitope alteration. Hypermutation by error-prone DNA damage repair (alternative end-joining, mismatch repair deficiency) drove CD22 escape. Acquired loss-of-function mutations in TP53 , ATM and CDKN2A were observed, suggesting compromise of the G1/S DNA damage checkpoint as a mechanism of evading InO-induced apoptosis. In conclusion, genetic alterations modulating CD22 expression and DNA damage response influence InO efficacy. The escape strategies within and beyond antigen loss to CD22-targeted therapy elucidated in this study provide insights into improving therapeutic approaches and overcoming resistance.

KEY POINTS

We identified multiple mechanisms of CD22 antigen escape from inotuzumab ozogamicin, including protein truncation, protein destabilization, and epitope alteration.Hypermutation caused by error-prone DNA damage repair was a driver of CD22 mutation and escape.

VISUAL ABSTRACT

Karyotyping of aneuploid and polyploid plants from low coverage whole-genome resequencing

BACKGROUND

Karyotype, as a basic characteristic of species, provides valuable information for fundamental theoretical research and germplasm resource innovation. However, traditional karyotyping techniques, including fluorescence in situ hybridization (FISH), are challenging and low in efficiency, especially when karyotyping aneuploid and polyploid plants. The use of low coverage whole-genome resequencing (lcWGR) data for karyotyping was explored, but existing methods are complicated and require control samples.

RESULTS

In this study, a new protocol for molecular karyotype analysis was provided, which proved to be a simpler, faster, and more accurate method, requiring no control. Notably, our method not only provided the copy number of each chromosome of an individual but also an accurate evaluation of the genomic contribution from its parents. Moreover, we verified the method through FISH and published resequencing data.

CONCLUSIONS

This method is of great significance for species evolution analysis, chromosome engineering, crop improvement, and breeding.

A Bioinformatics Toolkit for Next-Generation Sequencing in Clinical Oncology

Next-generation sequencing (NGS) has taken on major importance in clinical oncology practice. With the advent of targeted therapies capable of effectively targeting specific genomic alterations in cancer patients, the development of bioinformatics processes has become crucial. Thus, bioinformatics pipelines play an essential role not only in the detection and in identification of molecular alterations obtained from NGS data but also in the analysis and interpretation of variants, making it possible to transform raw sequencing data into meaningful and clinically useful information. In this review, we aim to examine the multiple steps of a bioinformatics pipeline as used in current clinical practice, and we also provide an updated list of the necessary bioinformatics tools. This resource is intended to assist researchers and clinicians in their genetic data analyses, improving the precision and efficiency of these processes in clinical research and patient care.

Low microsatellite instability: A distinct instability type in gastric cancer?

PURPOSE

We recently showed that low microsatellite instability (MSI-L) is associated with a good response to platinum/5-fluorouracil (5-FU) neoadjuvant chemotherapy (CTx) in gastric cancer. The purpose of this study was to characterize the instability pattern and to investigate an association of MSI-L tumors with mutations in genes of DNA repair pathways and with total tumor mutation burden (TMB).

METHODS

MSI patterns were compared between 67 MSI high (-H) and 35 MSI-L tumors. Whole-exome sequencing was performed in 34 microsatellite stable (MSS) and 20 MSI-L tumors after or without neoadjuvant CTx.

RESULTS

Of the 35 MSI-L tumors, 33 tumors had instability at a dinucleotide repeat marker. In the homologous recombination (HR) pathway, 10 of the 34 (29%) MSS and 10 of the 20 (50%) MSI-L tumors showed variants (p = 0.154). In the DNA damage tolerance pathway, 6 of the 34 (18%) MSS and 7 of the 20 (35%) MSI-L tumors had variants (p = 0.194). The HR deficiency score was similar in both tumor groups. TMB was significantly higher in MSI-L compared to MSS tumors after CTx (p = 0.046). In the MSS and MSI-L tumors without CTx no difference was observed (p = 1.00).

CONCLUSION

MSI-L due to instability at dinucleotide repeat markers was associated with increased TMB after neoadjuvant CTx treatment, indicating sensitivity to platinum/5-FU CTx. If confirmed in further studies, this could contribute to refined chemotherapeutic options including immune-based strategies for GC patients with MSI-L tumors.

Systematic evaluation of multiple NGS platforms for structural variants detection

Structural variations (SV) are critical genome changes affecting human diseases. Although many hybridization-based methods exist, evaluating SVs through next-generation sequencing (NGS) data is still necessary for broader research exploration. Here, we comprehensively compared the performance of 16 SV callers and multiple NGS platforms using NA12878 whole genome sequencing (WGS) datasets. The results indicated that several SV callers performed well relatively, such as Manta, GRIDSS, LUMPY, TARDIS, FermiKit, and Wham. Meanwhile, all NGS platforms have a similar performance using a single software. Additionally, we found that the source of undetected SVs was mostly from long reads datasets, therefore, the more appropriate strategy for accurate SV detection will be an integration of long and shorter reads in the future. At present, in the period of NGS as a mainstream method in bioinformatics, our study would provide helpful and comprehensive guidelines for specific categories of SV research.

CNV-FB: A Feature bagging strategy-based approach to detect copy number variants from NGS data

Copy number variation (CNV), as a type of genomic structural variation, accounts for a large proportion of structural variation and is related to the pathogenesis and susceptibility to some human diseases, playing an important role in the development and change of human diseases. The development of next-generation sequencing technology (NGS) provides strong support for the design of CNV detection algorithms. Although a large number of methods have been developed to detect CNVs using NGS data, it is still considered a difficult problem to detect CNVs with low purity and coverage. In this paper, a new calculation method CNV-FB is proposed to detect CNVs from NGS data. The core idea of CNV-FB is to randomly sample the read depth values of the genome fragment, and then each sample is individually detected for outliers, and finally combined into a final outlier score. The CNV-FB method was applied to simulation data and real data experiments and compared with the other five methods of the same type. The results show that the CNV-FB method has a better detection effect than other methods. Therefore, the CNV-FB method may be an effective algorithm for detecting genomic mutations.

Completing a genomic characterisation of microscopic tumour samples with copy number

BACKGROUND

Genomic insights in settings where tumour sample sizes are limited to just hundreds or even tens of cells hold great clinical potential, but also present significant technical challenges. We previously developed the DigiPico sequencing platform to accurately identify somatic mutations from such samples.

RESULTS

Here, we complete this genomic characterisation with copy number. We present a novel protocol, PicoCNV, to call allele-specific somatic copy number alterations from picogram quantities of tumour DNA. We find that PicoCNV provides exactly accurate copy number in 84% of the genome for even the smallest samples, and demonstrate its clinical potential in maintenance therapy.

CONCLUSIONS

PicoCNV complements our existing platform, allowing for accurate and comprehensive genomic characterisations of cancers in settings where only microscopic samples are available.

Therapeutic targeting of CPSF3-dependent transcriptional termination in ovarian cancer

Transcriptional dysregulation is a recurring pathogenic hallmark and an emerging therapeutic vulnerability in ovarian cancer. Here, we demonstrated that ovarian cancer exhibited a unique dependency on the regulatory machinery of transcriptional termination, particularly, cleavage and polyadenylation specificity factor (CPSF) complex. Genetic abrogation of multiple CPSF subunits substantially hampered neoplastic cell viability, and we presented evidence that their indispensable roles converged on the endonuclease CPSF3. Mechanistically, CPSF perturbation resulted in lengthened 3'-untranslated regions, diminished intronic polyadenylation and widespread transcriptional readthrough, and consequently suppressed oncogenic pathways. Furthermore, we reported the development of specific CPSF3 inhibitors building upon the benzoxaborole scaffold, which exerted potent antitumor activity. Notably, CPSF3 blockade effectively exacerbated genomic instability by down-regulating DNA damage repair genes and thus acted in synergy with poly(adenosine 5'-diphosphate-ribose) polymerase inhibition. These findings establish CPSF3-dependent transcriptional termination as an exploitable driving mechanism of ovarian cancer and provide a promising class of boron-containing compounds for targeting transcription-addicted human malignancies.

Genetic characteristics of common variable immunodeficiency patients with autoimmunity

Background: The pathogenesis of common variable immunodeficiency disorder (CVID) is complex, especially when combined with autoimmunity. Genetic factors may be potential explanations for this complex situation, and whole genome sequencing (WGS) provide the basis for this potential. Methods: Genetic information of patients with CVID with autoimmunity, together with their first-degree relatives, was collected through WGS. The association between genetic factors and clinical phenotypes was studied using genetic analysis strategies such as sporadic and pedigree. Results: We collected 42 blood samples for WGS (16 CVID patients and 26 first-degree relatives of healthy controls). Through pedigree, sporadic screening strategies and low-frequency deleterious screening of rare diseases, we obtained 9,148 mutation sites, including 8,171 single-nucleotide variants (SNVs) and 977 Insertion-deletions (InDels). Finally, we obtained a total of 28 candidate genes (32 loci), of which the most common mutant was LRBA. The most common autoimmunity in the 16 patients was systematic lupus erythematosis. Through KEGG pathway enrichment, we identified the top ten signaling pathways, including "primary immunodeficiency", "JAK-STAT signaling pathway", and "T-cell receptor signaling pathway". We used PyMOL to predict and analyse the three-dimensional protein structures of the NFKB1, RAG1, TIRAP, NCF2, and MYB genes. In addition, we constructed a PPI network by combining candidate mutants with genes associated with CVID in the OMIM database via the STRING database. Conclusion: The genetic background of CVID includes not only monogenic origins but also oligogenic effects. Our study showed that immunodeficiency and autoimmunity may overlap in genetic backgrounds. Clinical Trial Registration: identifier ChiCTR2100044035.

Translational practice of fluorescence in situ hybridisation to identify neuroblastic tumours with TERT rearrangements

Recently, telomerase reverse transcriptase (TERT) gene rearrangements have been identified in neuroblastoma (NB), the typical pathological type of neuroblastic tumours (NTs); however, the prevalence of TERT rearrangements in other types of NT remains unknown. This study aimed to develop a practical method for detecting TERT defects and to evaluate the clinical relevance of TERT rearrangements as a biomarker for NT prognosis. A TERT break-apart probe for fluorescence in situ hybridisation (FISH) was designed, optimised, and applied to assess the genomic status of TERT in Chinese children with NTs at the Beijing Children's Hospital from 2016 to 2019. Clinical, histological, and genetic characteristics of TERT-rearranged NTs were further addressed. Genomic TERT rearrangements could be effectively detected by FISH and were mutually exclusive with MYCN amplification. TERT rearrangements were identified in 6.0% (38/633) of NTs overall, but 12.4% (31/250) in high-risk patients. TERT rearrangements identified a subtype of aggressive NTs with the characteristics of Stage 3/4, high-risk category, over 18 months old, and presenting all histological subtypes of NB and ganglioneuroblastoma nodular. Moreover, TERT rearrangements were significantly associated with elevated TERT expression levels and decreased survival chances. Multivariable analysis confirmed that it was an independent prognostic marker for NTs. FISH is an easily applicable method for evaluating TERT defects, which define a subgroup of NTs with unfavourable prognosis. TERT rearrangements would contribute to characterising NT molecular signatures in clinical practice.

A novel subclonal rearrangement of the STRN3::PDGFRB gene in de novo acute myeloid leukemia with NPM1 mutation and its leukemogenic effects

Chromosome translocations in the 5q31-33 region are associated with a range of hematologic malignancies, some of which involve the platelet-derived growth factor receptor beta (PDGFRB) gene. We report a case of acute myeloid leukemia (AML) with a mutation in the NPM1 gene (NPM1-mut AML) and a subclonal gene rearrangement involving the PDGFRB gene. We identified a novel fusion gene, STRN3::PDGFRB, resulting from t(5;14) (q32;q12) chromosomal rearrangement. Sequential FISH confirmed that ~15% of leukemic cells carried the PDGFRB gene rearrangement, which suggests that STRN3::PDGFRB is a previously unreported fusion gene in a subclone. Reverse transcription PCR (RT-PCR) and Sanger sequencing confirmed that the fusion gene consisted of STRN3 exon 7 fused to PDGFRB exon 11, resulting in a chimeric protein containing the coiled-coil domain of striatin-3 and the transmembrane and intracellular tyrosine kinase domains of the PDGFRB. The new protein exhibited distinct cytoplasmic localization and had leukemogenic effects, as demonstrated by its ability to transform Ba/F3 cells to growth factor independence and cause a fatal myelodysplastic/myeloproliferative neoplasm (MDS/MPN)-like disease in mice, which then transformant to T-cell lymphoblastic lymphoma in secondary recipients. Ba/F3 cells expressing STRN3::PDGFRB or ETV6::PDGFRB were sensitive to tyrosine kinase inhibitors (TKIs) and selinexor, but in vitro experiments showed that the combination of imatinib and selinexor had a marked synergistic effect, although only the imatinib alone group could prolong the survival of T-cell blast transformation recipient mice. Our findings demonstrate the leukemogenic effects of the novel fusion gene and provide insights into the clone evolution of AML, which can be influenced by therapy selection. Furthermore, our results provide insight into the potential therapeutic options for patients with this type of mutation, as well as the need for careful consideration of treatment selection to prevent undesirable side effects.

Shallow whole genome sequencing approach to detect Homologous Recombination Deficiency in the PAOLA-1/ENGOT-OV25 phase-III trial

The bevacizumab (bev)/olaparib (ola) maintenance regimen was approved for BRCA1/2-mutated (BRCAmut) and Homologous Recombination Deficient (HRD) high-grade Advanced Ovarian Cancer (AOC) first line setting, based on a significantly improved progression-free survival (PFS) compared to bev alone in the PAOLA-1/ENGOT-ov25 trial (NCT02477644), where HRD was detected by MyChoice CDx PLUS test. The academic shallowHRDv2 test was developed based on shallow whole-genome sequencing as an alternative to MyChoice. Analytical and clinical validities of shallowHRDv2 as compared to MyChoice on 449 PAOLA-1 tumor samples are presented. The overall agreement between shallowHRDv2 and MyChoice was 94% (369/394). Less non-contributive tests were observed with shallowHRDv2 (15/449; 3%) than with MyChoice (51/449; 11%). Patients with HRD tumors according to shallowHRDv2 (including BRCAmut) showed a significantly prolonged PFS with bev+ola versus bev (median PFS: 65.7 versus 20.3 months, hazard ratio (HR): 0.36 [95% CI: 0.24-0.53]). This benefit was significant also for BRCA1/2 wild-type tumors (40.8 versus 19.5 months, HR: 0.45 [95% CI: 0.26-0.76]). ShallowHRDv2 is a performant, clinically validated, and cost-effective test for HRD detection.

A human mucosal melanoma organoid platform for modeling tumor heterogeneity and exploring immunotherapy combination options

Mucosal melanoma (MM), an aggressive rare subtype of melanoma, is distinct from cutaneous melanoma and has poor prognoses. We addressed the lack of cell models for MM by establishing 30 organoids of human oral MM (OMM), which retained major histopathological and functional features of parental tumors. Organoid groups derived from chronologically or intratumorally distinct lesions within the same individual displayed heterogeneous genetics, expression profiles, and drug responses, indicating rapid tumor evolution and poor clinical response. Furthermore, transcriptome analysis revealed receptor tyrosine kinases (RTKs) signaling, particularly NGFR, a nerve growth factor receptor, was significantly up-regulated in OMMs and organoids from patients resistant to anti-programmed cell death protein 1 (anti-PD-1) therapy. Combining anti-PD-1 with anlotinib (a phase 2 multitarget RTK inhibitor for OMM) or NGFR knockdown enhanced the effective activity of immune cells in organoid-immune cell coculture systems. Together, our study suggested that OMM organoids serve as faithful models for exploring tumor evolution and immunotherapy combination strategies.

Newly discovered genomic mutation patterns in radiation-induced small intestinal tumors of ApcMin/+ mice

Among the small intestinal tumors that occur in irradiated mice of the established mouse model B6/B6-Chr18MSM-F1 ApcMin/+, loss of heterozygosity analysis can be utilized to estimate whether a deletion in the wild-type allele containing the Adenomatous polyposis coli (Apc) region (hereafter referred to as Deletion), a duplication in the mutant allele with a nonsense mutation at codon 850 of Apc (Duplication), or no aberration (Unidentified) has occurred. Previous research has revealed that the number of Unidentified tumors tends to increase with the radiation dose. In the present study, we investigated the molecular mechanisms underlying the development of an Unidentified tumor type in response to radiation exposure. The mRNA expression levels of Apc were significantly lower in Unidentified tumors than in normal tissues. We focused on epigenetic suppression as the mechanism underlying this decreased expression; however, hypermethylation of the Apc promoter region was not observed. To investigate whether deletions occur that cannot be captured by loss of heterozygosity analysis, we analyzed chromosome 18 using a customized array comparative genomic hybridization approach designed to detect copy-number changes in chromosome 18. However, the copy number of the Apc region was not altered in Unidentified tumors. Finally, gene mutation analysis of the Apc region using next-generation sequencing suggested the existence of a small deletion (approximately 3.5 kbp) in an Unidentified tumor from a mouse in the irradiated group. Furthermore, nonsense and frameshift mutations in Apc were found in approximately 30% of the Unidentified tumors analyzed. These results suggest that radiation-induced Unidentified tumors arise mainly due to decreased Apc expression of an unknown regulatory mechanism that does not depend on promoter hypermethylation, and that some tumors may result from nonsense mutations which are as-yet undefined point mutations.

scNanoHi-C: a single-cell long-read concatemer sequencing method to reveal high-order chromatin structures within individual cells

The high-order three-dimensional (3D) organization of regulatory genomic elements provides a topological basis for gene regulation, but it remains unclear how multiple regulatory elements across the mammalian genome interact within an individual cell. To address this, herein, we developed scNanoHi-C, which applies Nanopore long-read sequencing to explore genome-wide proximal high-order chromatin contacts within individual cells. We show that scNanoHi-C can reliably and effectively profile 3D chromatin structures and distinguish structure subtypes among individual cells. This method could also be used to detect genomic variations, including copy-number variations and structural variations, as well as to scaffold the de novo assembly of single-cell genomes. Notably, our results suggest that extensive high-order chromatin structures exist in active chromatin regions across the genome, and multiway interactions between enhancers and their target promoters were systematically identified within individual cells. Altogether, scNanoHi-C offers new opportunities to investigate high-order 3D genome structures at the single-cell level.

Whole genome scanning of a Mediterranean basin hotspot collection provides new insights into olive tree biodiversity and biology

Olive tree (Olea europaea L. subsp. europaea var. europaea) is one of the most important species of the Mediterranean region and one of the most ancient species domesticated. The availability of whole genome assemblies and annotations of olive tree cultivars and oleaster (O. europaea subsp. europaea var. sylvestris) has contributed to a better understanding of genetic and genomic differences between olive tree cultivars. However, compared to other plant species there is still a lack of genomic resources for olive tree populations that span the entire Mediterranean region. In the present study we developed the most complete genomic variation map and the most comprehensive catalog/resource of molecular variation to date for 89 olive tree genotypes originating from the entire Mediterranean basin, revealing the genetic diversity of this commercially significant crop tree and explaining the divergence/similarity among different variants. Additionally, the monumental ancient tree 'Throuba Naxos' was studied to characterize the potential origin or routes of olive tree domestication. Several candidate genes known to be associated with key agronomic traits, including olive oil quality and fruit yield, were uncovered by a selective sweep scan to be under selection pressure on all olive tree chromosomes. To further exploit the genomic and phenotypic resources obtained from the current work, genome-wide association analyses were performed for 23 morphological and two agronomic traits. Significant associations were detected for eight traits that provide valuable candidates for fruit tree breeding and for deeper understanding of olive tree biology.

Investigation of evolutionary dynamics for drug resistance in 3D spheroid model system using cellular barcoding technology

Complex evolutionary dynamics governing the drug resistance is one of the major challenges in cancer treatment. Understanding these mechanisms requires a sequencing technology with higher resolution to delineate whether pre-existing or de novo drug mechanisms are behind the drug resistance. Combining this technology with clinically very relevant model system, namely 3D spheroids, better mimicking tumorigenesis and drug resistance have so far been lacking. Thus, we sought to establish dabrafenib and irinotecan resistant derivatives of barcoded 3D spheroids with the ultimate aim to quantify the selection-induced clonal dynamics and identify the genomic determinants in this model system. We found that dabrafenib and irinotecan induced drug resistance in 3D-HT-29 and 3D-HCT-116 spheroids are mediated by pre-existing and de novo resistant barcodes, indicating the presence of polyclonal drug resistance in this system. Moreover, whole-exome sequencing analysis found chromosomal gains and mutations associated with dabrafenib and irinotecan resistance in 3D-HT-29 and 3D-HCT-116 spheroids. Last, we show that dabrafenib and irinotecan resistance are also mediated by multiple drug resistance by detection of upregulation of the drug efflux pumps, ABCB1 and ABCG2, in our spheroid model system. Overall, we present the quantification of drug resistance and evolutionary dynamics in spheroids for the first time using cellular barcoding technology and the underlying genomic determinants of the drug resistance in our model system.

Copy Number Variation on ABCC2-DNMBP Loci Affects the Diversity and Composition of the Fecal Microbiota in Pigs

Genetic variation in the pig genome partially modulates the composition of porcine gut microbial communities. Previous studies have been focused on the association between single nucleotide polymorphisms (SNPs) and the gut microbiota, but little is known about the relationship between structural variants and fecal microbial traits. The main goal of this study was to explore the association between porcine genome copy number variants (CNVs) and the diversity and composition of pig fecal microbiota. For this purpose, we used whole-genome sequencing data to undertake a comprehensive identification of CNVs followed by a genome-wide association analysis between the estimated CNV status and the fecal bacterial diversity in a commercial Duroc pig population. A CNV predicted as gain (DUP) partially harboring ABCC2-DNMBP loci was associated with richness (P = 5.41 × 10-5, false discovery rate [FDR] = 0.022) and Shannon α-diversity (P = 1.42 × 10-4, FDR = 0.057). The in silico predicted gain of copies was validated by real-time quantitative PCR (qPCR), and its segregation, and positive association with the richness and Shannon α-diversity of the porcine fecal bacterial ecosystem was confirmed in an unrelated F1 (Duroc × Iberian) cross. Our results advise the relevance of considering the role of host-genome structural variants as potential modulators of microbial ecosystems and suggest the ABCC2-DNMBP CNV as a host-genetic factor for the modulation of the diversity and composition of the fecal microbiota in pigs. IMPORTANCE A better understanding of the environmental and host factors modulating gut microbiomes is a topic of greatest interest. Recent evidence suggests that genetic variation in the pig genome partially controls the composition of porcine gut microbiota. However, since previous studies have been focused on the association between single nucleotide polymorphisms and the fecal microbiota, little is known about the relationship between other sources of genetic variation, like the structural variants and microbial traits. Here, we identified, experimentally validated, and replicated in an independent population a positive link between the gain of copies of ABCC2-DNMBP loci and the diversity and composition of pig fecal microbiota. Our results advise the relevance of considering the role of host-genome structural variants as putative modulators of microbial ecosystems and open the possibility of implementing novel holobiont-based management strategies in breeding programs for the simultaneous improvement of microbial traits and host performance.

Whole exome sequencing of well-differentiated liposarcoma and dedifferentiated liposarcoma in older woman: a case report

Background

Common kinds of soft tissue sarcomas (STS) include well-differentiated liposarcoma (WDLPS) and dedifferentiated liposarcoma (DDLPS). In this case, we present a comprehensive clinical profile of a patient who underwent multiple recurrences during the progression from WDLPS to DDLPS.

Case presentation

A 62-year-old Asian female underwent retroperitoneal resection of a large tumor 11 years ago, the initial pathology revealed a fibrolipoma-like lesion. Over the next six years, the patient underwent three resections for recurrence of abdominal tumors. Postoperative histology shows mature adipose tissue with scattered "adipoblast"-like cells with moderate-to-severe heterogeneous spindle cells, pleomorphic cells, or tumor giant cells. Immunohistochemistry (IHC) demonstrated positive staining for MDM2 and CDK4, confirming that the abdominal tumor was WDLPS and gradually progressing to DDLPS. Post-operative targeted sequencing and IHC confirmed the POC1B::ROS1 fusion gene in DDLPS. Whole-exome sequencing (WES) revealed that WDLPS and DDLPS shared similar somatic mutations and copy number variations (CNVs), whereas DDLPS had more mutated genes and a higher and more concentrated amplification of the chromosome 12q region. Furthermore, somatic mutations in DDLPS were significantly reduced after treatment with CDK4 inhibitors, while CNVs remained elevated.

Conclusion

Due to the high likelihood of recurrence of liposarcoma, various effective treatments should be taken into consideration even if surgery is the primary treatment for recurrent liposarcoma. To effectively control the course of the disease following surgery, combination targeted therapy may be a viable alternative to chemotherapy and radiotherapy in the treatment of liposarcoma.

Non-homologous end joining shapes the genomic rearrangement landscape of chromothripsis from mitotic errors

Errors in mitosis can generate micronuclei that entrap mis-segregated chromosomes, which are susceptible to catastrophic fragmentation through a process termed chromothripsis. The reassembly of fragmented chromosomes by error-prone DNA double-strand break (DSB) repair generates a spectrum of simple and complex genomic rearrangements that are associated with human cancers and disorders. How specific DSB repair pathways recognize and process these lesions remains poorly understood. Here we used CRISPR/Cas9 to systematically inactivate distinct DSB processing or repair pathways and interrogated the rearrangement landscape of fragmented chromosomes from micronuclei. Deletion of canonical non-homologous end joining (NHEJ) components, including DNA-PKcs, LIG4, and XLF, substantially reduced the formation of complex rearrangements and shifted the rearrangement landscape toward simple alterations without the characteristic patterns of cancer-associated chromothripsis. Following reincorporation into the nucleus, fragmented chromosomes localize within micronuclei bodies (MN bodies) and undergo successful ligation by NHEJ within a single cell cycle. In the absence of NHEJ, chromosome fragments were rarely engaged by polymerase theta-mediated alternative end-joining or recombination-based mechanisms, resulting in delayed repair kinetics and persistent 53BP1-labeled MN bodies in the interphase nucleus. Prolonged DNA damage signaling from unrepaired fragments ultimately triggered cell cycle arrest. Thus, we provide evidence supporting NHEJ as the exclusive DSB repair pathway generating complex rearrangements following chromothripsis from mitotic errors.

Evolutionary histories of breast cancer and related clones

Recent studies have documented frequent evolution of clones carrying common cancer mutations in apparently normal tissues, which are implicated in cancer development1-3. However, our knowledge is still missing with regard to what additional driver events take place in what order, before one or more of these clones in normal tissues ultimately evolve to cancer. Here, using phylogenetic analyses of multiple microdissected samples from both cancer and non-cancer lesions, we show unique evolutionary histories of breast cancers harbouring der(1;16), a common driver alteration found in roughly 20% of breast cancers. The approximate timing of early evolutionary events was estimated from the mutation rate measured in normal epithelial cells. In der(1;16)(+) cancers, the derivative chromosome was acquired from early puberty to late adolescence, followed by the emergence of a common ancestor by the patient's early 30s, from which both cancer and non-cancer clones evolved. Replacing the pre-existing mammary epithelium in the following years, these clones occupied a large area within the premenopausal breast tissues by the time of cancer diagnosis. Evolution of multiple independent cancer founders from the non-cancer ancestors was common, contributing to intratumour heterogeneity. The number of driver events did not correlate with histology, suggesting the role of local microenvironments and/or epigenetic driver events. A similar evolutionary pattern was also observed in another case evolving from an AKT1-mutated founder. Taken together, our findings provide new insight into how breast cancer evolves.

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.

A complex structural variant near SOX3 causes X-linked split-hand/foot malformation

Split-hand/foot malformation (SHFM) is a congenital limb defect most typically presenting with median clefts in hands and/or feet, that can occur in a syndromic context as well as in isolated form. SHFM is caused by failure to maintain normal apical ectodermal ridge function during limb development. Although several genes and contiguous gene syndromes are implicated in the monogenic etiology of isolated SHFM, the disorder remains genetically unexplained for many families and associated genetic loci. We describe a family with isolated X-linked SHFM, for which the causative variant could be detected after a diagnostic journey of 20 years. We combined well-established approaches including microarray-based copy number variant analysis and fluorescence in situ hybridization coupled with optical genome mapping and whole genome sequencing. This strategy identified a complex structural variant (SV) comprising a 165-kb gain of 15q26.3 material ([GRCh37/hg19] chr15:99795320-99960362dup) inserted in inverted position at the site of a 38-kb deletion on Xq27.1 ([GRCh37/hg19] chrX:139481061-139518989del). In silico analysis suggested that the SV disrupts the regulatory framework on the X chromosome and may lead to SOX3 misexpression. We hypothesize that SOX3 dysregulation in the developing limb disturbed the fine balance between morphogens required for maintaining AER function, resulting in SHFM in this family.

Spatial and clonality-resolved 3D cancer genome alterations reveal enhancer-hijacking as a potential prognostic marker for colorectal cancer

The regulatory effect of non-coding large-scale structural variations (SVs) on proto-oncogene activation remains unclear. This study investigated SV-mediated gene dysregulation by profiling 3D cancer genome maps from 40 patients with colorectal cancer (CRC). We developed a machine learning-based method for spatial characterization of the altered 3D cancer genome. This revealed a frequent establishment of "de novo chromatin contacts" that can span multiple topologically associating domains (TADs) in addition to the canonical TAD fusion/shuffle model. Using this information, we precisely identified super-enhancer (SE)-hijacking and its clonal characteristics. Clonal SE-hijacking genes, such as TOP2B, are recurrently associated with cell-cycle/DNA-processing functions, which can potentially be used as CRC prognostic markers. Oncogene activation and increased drug resistance due to SE-hijacking were validated by reconstructing the patient's SV using CRISPR-Cas9. Collectively, the spatial and clonality-resolved analysis of the 3D cancer genome reveals regulatory principles of large-scale SVs in oncogene activation and their clinical implications.

OpenPBTA: The Open Pediatric Brain Tumor Atlas

Pediatric brain and spinal cancers are collectively the leading disease-related cause of death in children; thus, we urgently need curative therapeutic strategies for these tumors. To accelerate such discoveries, the Children's Brain Tumor Network (CBTN) and Pacific Pediatric Neuro-Oncology Consortium (PNOC) created a systematic process for tumor biobanking, model generation, and sequencing with immediate access to harmonized data. We leverage these data to establish OpenPBTA, an open collaborative project with over 40 scalable analysis modules that genomically characterize 1,074 pediatric brain tumors. Transcriptomic classification reveals universal TP53 dysregulation in mismatch repair-deficient hypermutant high-grade gliomas and TP53 loss as a significant marker for poor overall survival in ependymomas and H3 K28-mutant diffuse midline gliomas. Already being actively applied to other pediatric cancers and PNOC molecular tumor board decision-making, OpenPBTA is an invaluable resource to the pediatric oncology community.

Defining the landscape of circular RNAs in neuroblastoma unveils a global suppressive function of MYCN

Circular RNAs (circRNAs) are a regulatory RNA class. While cancer-driving functions have been identified for single circRNAs, how they modulate gene expression in cancer is not well understood. We investigate circRNA expression in the pediatric malignancy, neuroblastoma, through deep whole-transcriptome sequencing in 104 primary neuroblastomas covering all risk groups. We demonstrate that MYCN amplification, which defines a subset of high-risk cases, causes globally suppressed circRNA biogenesis directly dependent on the DHX9 RNA helicase. We detect similar mechanisms in shaping circRNA expression in the pediatric cancer medulloblastoma implying a general MYCN effect. Comparisons to other cancers identify 25 circRNAs that are specifically upregulated in neuroblastoma, including circARID1A. Transcribed from the ARID1A tumor suppressor gene, circARID1A promotes cell growth and survival, mediated by direct interaction with the KHSRP RNA-binding protein. Our study highlights the importance of MYCN regulating circRNAs in cancer and identifies molecular mechanisms, which explain their contribution to neuroblastoma pathogenesis.

MIRACUM-Pipe: An Adaptable Pipeline for Next-Generation Sequencing Analysis, Reporting, and Visualization for Clinical Decision Making

(1) Background: Next-generation sequencing (NGS) of patients with advanced tumors is becoming an established method in Molecular Tumor Boards. However, somatic variant detection, interpretation, and report generation, require in-depth knowledge of both bioinformatics and oncology. (2) Methods: MIRACUM-Pipe combines many individual tools into a seamless workflow for comprehensive analyses and annotation of NGS data including quality control, alignment, variant calling, copy number variation estimation, evaluation of complex biomarkers, and RNA fusion detection. (3) Results: MIRACUM-Pipe offers an easy-to-use, one-prompt standardized solution to analyze NGS data, including quality control, variant calling, copy number estimation, annotation, visualization, and report generation. (4) Conclusions: MIRACUM-Pipe, a versatile pipeline for NGS, can be customized according to bioinformatics and clinical needs and to support clinical decision-making with visual processing and interactive reporting.

Culture conditions of mouse ESCs impact the tumor appearance in vivo

Various culture conditions by small molecules have been explored to extend pluripotency of stem cells, but their impacts on cell fate in vivo remain elusive. We systematically compared the effects of various culture conditions on the pluripotency and cell fate in vivo of mouse embryonic stem cells (ESCs) by tetraploid embryo complementation assay. Conventional ESC cultures in serum/LIF-based condition produced complete ESC mice and also the survival to adulthood at the highest rates of all other chemical-based cultures. Moreover, long-term examination of the survived ESC mice demonstrated that conventional ESC cultures did not lead to visible abnormality for up to 1.5-2 years, whereas the prolonged chemical-based cultures developed retroperitoneal atypical teratomas or leiomyomas. The chemical-based cultures exhibited transcriptomes and epigenomes that typically differed from those of conventional ESC cultures. Our results warrant further refinement of culture conditions in promoting the pluripotency and safety of ESCs in future applications.

NOTCH4ΔL12_16 sensitizes lung adenocarcinomas to EGFR-TKIs through transcriptional down-regulation of HES1

Resistance to epidermal growth factor tyrosine kinase inhibitors (EGFR-TKI) remains one of the major challenges in lung adenocarcinoma (LUAD) therapy. Here, we find an increased frequency of the L12_16 amino acid deletion mutation in the signal peptide region of NOTCH4 (NOTCH4^ΔL12_16) in EGFR-TKI-sensitive patients. Functionally, exogenous induction of NOTCH4^ΔL12_16 in EGFR-TKI -resistant LUAD cells sensitizes them to EGFR-TKIs. This process is mainly mediated by the reduction of the intracellular domain of NOTCH4 (NICD4) caused by the NOTCH4^ΔL12_16 mutation, which results in a lower localization of NOTCH4 in the plasma membrane. Mechanistically, NICD4 transcriptionally upregulates the expression of HES1 by competitively binding to the gene promoter relative to p-STAT3. Because p-STAT3 can downregulate the expression of HES1 in EGFR-TKI-resistant LUAD cells, the reduction of NICD4 induced by NOTCH4^ΔL12_16 mutation leads to a decrease in HES1. Moreover, inhibition of the NOTCH4-HES1 pathway using inhibitors and siRNAs abolishes the resistance of EGFR-TKI. Overall, we report that the NOTCH4^ΔL12_16 mutation sensitizes LUAD patients to EGFR-TKIs through transcriptional down-regulation of HES1 and that targeted blockade of this signaling cohort could reverse EGFR-TKI -resistance in LUAD, providing a potential approach to overcome resistance to EGFR-TKI -therapy.

Loss of Heterozygosity and Its Importance in Evolution

Loss of heterozygosity (LOH) is a mitotic recombination event that converts heterozygous loci to homozygous loci. This mutation event is widespread in organisms that have asexual reproduction like budding yeasts, and is also an important and frequent mutation event in tumorigenesis. Mutation accumulation studies have demonstrated that LOH occurs at a rate higher than the point mutation rate, and can impact large portions of the genome. Laboratory evolution experiments of heterozygous yeasts have revealed that LOH often unmasks beneficial recessive alleles that can confer large fitness advantages. Here, I highlight advances in understanding dominance, fitness, and phenotypes in laboratory evolved heterozygous yeast strains. I discuss best practices for detecting LOH in intraspecific and interspecific evolved clones and populations. Utilizing heterozygous strain backgrounds in laboratory evolution experiments offers an opportunity to advance our understanding of this important mutation type in shaping adaptation and genome evolution in wild, domesticated, and clinical populations.

Identifying Targets of Selection in Laboratory Evolution Experiments

Adaptive evolution navigates a balance between chance and determinism. The stochastic processes of mutation and drift generate phenotypic variation; however, once mutations reach an appreciable frequency in the population, their fate is governed by the deterministic action of selection, enriching for favorable genotypes and purging the less-favorable ones. The net result is that replicate populations will traverse similar-but not identical-pathways to higher fitness. This parallelism in evolutionary outcomes can be leveraged to identify the genes and pathways under selection. However, distinguishing between beneficial and neutral mutations is challenging because many beneficial mutations will be lost due to drift and clonal interference, and many neutral (and even deleterious) mutations will fix by hitchhiking. Here, we review the best practices that our laboratory uses to identify genetic targets of selection from next-generation sequencing data of evolved yeast populations. The general principles for identifying the mutations driving adaptation will apply more broadly.

The Swedish childhood tumor biobank: systematic collection and molecular characterization of all pediatric CNS and other solid tumors in Sweden

The Swedish Childhood Tumor Biobank (BTB) is a nonprofit national infrastructure for collecting tissue samples and genomic data from pediatric patients diagnosed with central nervous system (CNS) and other solid tumors. The BTB is built on a multidisciplinary network established to provide the scientific community with standardized biospecimens and genomic data, thereby improving knowledge of the biology, treatment and outcome of childhood tumors. As of 2022, over 1100 fresh-frozen tumor samples are available for researchers. We present the workflow of the BTB from sample collection and processing to the generation of genomic data and services offered. To determine the research and clinical utility of the data, we performed bioinformatics analyses on next-generation sequencing (NGS) data obtained from a subset of 82 brain tumors and patient blood-derived DNA combined with methylation profiling to enhance the diagnostic accuracy and identified germline and somatic alterations with potential biological or clinical significance. The BTB procedures for collection, processing, sequencing, and bioinformatics deliver high-quality data. We observed that the findings could impact patient management by confirming or clarifying the diagnosis in 79 of the 82 tumors and detecting known or likely driver mutations in 68 of 79 patients. In addition to revealing known mutations in a broad spectrum of genes implicated in pediatric cancer, we discovered numerous alterations that may represent novel driver events and specific tumor entities. In summary, these examples reveal the power of NGS to identify a wide number of actionable gene alterations. Making the power of NGS available in healthcare is a challenging task requiring the integration of the work of clinical specialists and cancer biologists; this approach requires a dedicated infrastructure, as exemplified here by the BTB.

TRIM32 biallelic defects cause limb-girdle muscular dystrophy R8: identification of two novel mutations and investigation of genotype-phenotype correlation

BACKGROUND

Limb-girdle muscular dystrophy R8 (LGMD R8) is a rare autosomal recessive muscle disease caused by TRIM32 gene biallelic defects. The genotype-phenotype correlation of this disease has been reported poorly. Here, we report a Chinese family with two female LGMD R8 patients.

METHODS

We performed whole-genome sequencing (WGS) and Sanger sequencing on the proband. Meanwhile, the function of mutant TRIM32 protein was analyzed by bioinformatics and experimental analysis. In addition, a summary of the reported TRIM32 deletions and point mutations and an investigation of genotype-phenotype correlation were performed through a combined analysis of the two patients and other cases reported in previous literature.

RESULTS

The two patients displayed typical symptoms of LGMD R8, which worsened during pregnancy. Genetic analysis by whole-genome sequencing (WGS) and Sanger sequencing showed that the patients were compound heterozygotes of a novel deletion (chr9.hg19:g.119431290_119474250del) and a novel missense mutation (TRIM32:c.1700A > G, p.H567R). The deletion encompassed 43 kb and resulted in the removal of the entire TRIM32 gene. The missense mutation altered the structure and further affected function by interfering with the self-association of the TRIM32 protein. Females with LGMD R8 showed less severe symptoms than males, and patients carrying two mutations in NHL repeats of the TRIM32 protein had earlier disease onset and more severe symptoms than other patients.

CONCLUSIONS

This research extended the spectrum of TRIM32 mutations and firstly provided useful data on the genotype-phenotype correlation, which is valuable for the accurate diagnosis and genetic counseling of LGMD R8.