GRIDSS2: comprehensive characterisation of somatic structural variation using single breakend variants and structural variant phasing

Unbalanced MYC break-apart FISH patterns indicate the presence of a MYC rearrangement in HGBCL-DH-BCL2

ABSTRACT

Fluorescence in situ hybridization (FISH) using break-apart probes is recommended for identifying high-grade B-cell lymphoma with MYC and BCL2 rearrangements (HGBCL-DH-BCL2). Unbalanced MYC break-apart patterns, in which the red or green signal is lost, are commonly reported as an equivocal result by clinical laboratories. In a cohort of 297 HGBCL-DH-BCL2, 13% of tumors had unbalanced MYC break-apart patterns with loss of red (LR; 2%) or loss of green (LG; 11%) signal. To determine the significance of these patterns, MYC rearrangements were characterized by sequencing in 130 HGBCL-DH-BCL2, including 3 LR and 14 LG tumors. A MYC rearrangement was identified for 71% of tumors with LR or LG patterns, with the majority involving immunoglobulin loci or other recurrent MYC rearrangement partners. The architecture of these rearrangements consistently preserved the rearranged MYC allele, with the MYC gene predicted to be on the derivative chromosome containing the signal that is still present in nearly all cases. MYC protein expression, MYC messenger RNA expression, and the proportion of tumors expressing the dark-zone signature was not significantly different between balanced and unbalanced groups. These results support a recommendation that unbalanced MYC break-apart FISH patterns be reported as positive for MYC rearrangement in the context of diagnosing HGBCL-DH-BCL2.

Performance of somatic structural variant calling in lung cancer using Oxford Nanopore sequencing technology

BACKGROUND

Lung cancer is a heterogeneous disease and the primary cause of cancer-related mortality worldwide. Somatic mutations, including large structural variants, are important biomarkers in lung cancer for selecting targeted therapy. Genomic studies in lung cancer have been conducted using short-read sequencing. Emerging long-read sequencing technologies are a promising alternative to study somatic structural variants, however there is no current consensus on how to process data and call somatic events. In this study, we preformed whole genome sequencing of lung cancer and matched non-tumour samples using long and short read sequencing to comprehensively benchmark three sequence aligners and seven structural variant callers comprised of generic callers (SVIM, Sniffles2, DELLY in generic mode and cuteSV) and somatic callers (Severus, SAVANA, nanomonsv and DELLY in somatic modes).

RESULTS

Different combinations of aligners and variant callers influenced somatic structural variant detection. The choice of caller had a significant influence on somatic structural variant detection in terms of variant type, size, sensitivity, and accuracy. The performance of each variant caller was assessed by comparing to somatic structural variants identified by short-read sequencing. When compared to somatic structural variants detected with short-read sequencing, more events were detected with long-read sequencing. The mean recall of somatic variant events identified by long-read sequencing was higher for the somatic callers (72%) than generic callers (53%). Among the somatic callers when using the minimap2 aligner, SAVANA and Severus achieved the highest recall at 79.5% and 79.25% respectively, followed by nanomonsv with a recall of 72.5%.

CONCLUSION

Long-read sequencing can identify somatic structural variants in clincal samples. The longer reads have the potential to improve our understanding of cancer development and inform personalized cancer treatment.

CRISPR-mediated megabase-scale transgene de-duplication to generate a functional single-copy full-length humanized DMD mouse model

BACKGROUND

The development of sequence-specific precision treatments like CRISPR gene editing therapies for Duchenne muscular dystrophy (DMD) requires sequence humanized animal models to enable the direct clinical translation of tested strategies. The current available integrated transgenic mouse model containing the full-length human DMD gene, Tg(DMD)72Thoen/J (hDMDTg), has been found to have two copies of the transgene per locus in a tail-to-tail orientation, which does not accurately simulate the true (single) copy number of the DMD gene. This duplication also complicates analysis when testing CRISPR therapy editing outcomes, as large genetic alterations and rearrangements can occur between the cut sites on the two transgenes.

RESULTS

To address this, we performed long read nanopore sequencing on hDMDTg mice to better understand the structure of the duplicated transgenes. Following that, we performed a megabase-scale deletion of one of the transgenes by CRISPR zygotic microinjection to generate a single-copy, full-length, humanized DMD transgenic mouse model (hDMDTgSc). Functional, molecular, and histological characterisation shows that the single remaining human transgene retains its function and rescues the dystrophic phenotype caused by endogenous murine Dmd knockout.

CONCLUSIONS

Our unique hDMDTgSc mouse model simulates the true copy number of the DMD gene, and can potentially be used for the further generation of DMD disease models that would be better suited for the pre-clinical assessment and development of sequence specific CRISPR therapies.

ONCOLINER: A new solution for monitoring, improving, and harmonizing somatic variant calling across genomic oncology centers

The characterization of somatic genomic variation associated with the biology of tumors is fundamental for cancer research and personalized medicine, as it guides the reliability and impact of cancer studies and genomic-based decisions in clinical oncology. However, the quality and scope of tumor genome analysis across cancer research centers and hospitals are currently highly heterogeneous, limiting the consistency of tumor diagnoses across hospitals and the possibilities of data sharing and data integration across studies. With the aim of providing users with actionable and personalized recommendations for the overall enhancement and harmonization of somatic variant identification across research and clinical environments, we have developed ONCOLINER. Using specifically designed mosaic and tumorized genomes for the analysis of recall and precision across somatic SNVs, insertions or deletions (indels), and structural variants (SVs), we demonstrate that ONCOLINER is capable of improving and harmonizing genome analysis across three state-of-the-art variant discovery pipelines in genomic oncology.

Evolutionary adaptation under climate change: Aedes sp. demonstrates potential to adapt to warming

Climate warming is expected to shift the distributions of mosquitoes and mosquito-borne diseases, facilitating expansions at cool range edges and contractions at warm range edges. However, whether mosquito populations could maintain their warm edges through evolutionary adaptation remains unknown. Here, we investigate the potential for thermal adaptation in Aedes sierrensis, a congener of the major disease vector species that experiences large thermal gradients in its native range, by assaying tolerance to prolonged and acute heat exposure, and its genetic basis in a diverse, field-derived population. We found pervasive evidence of heritable genetic variation in acute heat tolerance, which phenotypically trades off with tolerance to prolonged heat exposure. A simple evolutionary model based on our data shows that the estimated maximum rate of evolutionary adaptation in mosquito heat tolerance typically exceeds that of projected climate warming under idealized conditions. Our findings indicate that natural mosquito populations may have the potential to track projected warming via genetic adaptation. Prior climate-based projections may thus underestimate the range of mosquito and mosquito-borne disease distributions under future climate conditions.

Peering into the gaps: Long-read sequencing illuminates structural variants and genomic evolution in the Australasian snapper

Even before genome sequencing, genetic resources have supported species management and breeding programs. Current technologies, such as long-read sequencing, resolve complex genomic regions, like those rich in repeats or high in GC content. Improved genome contiguity enhances accuracy in identifying structural variants (SVs) and transposable elements (TEs). We present an improved genome assembly and SV catalogue for the Australasian snapper (Chrysophrys auratus). The new assembly is more contiguous, allowing for putative identification of 14 centromeres and transfer of 26,115 gene annotations from yellowfin seabream. Compared to the previous assembly, 35,000 additional SVs, including larger and more complex rearrangements, were annotated. SVs and TEs exhibit a distribution pattern skewed towards chromosome ends, likely influenced by recombination. Some SVs overlap with growth-related genes, underscoring their significance. This upgraded genome serves as a foundation for studying natural and artificial selection, offers a reference for related species, and sheds light on genome dynamics shaped by evolution.

Implications of noncoding regulatory functions in the development of insulinomas

Insulinomas are rare neuroendocrine tumors arising from pancreatic β cells, characterized by aberrant proliferation and altered insulin secretion, leading to glucose homeostasis failure. With the aim of uncovering the role of noncoding regulatory regions and their aberrations in the development of these tumors, we coupled epigenetic and transcriptome profiling with whole-genome sequencing. As a result, we unraveled somatic mutations associated with changes in regulatory functions. Critically, these regions impact insulin secretion, tumor development, and epigenetic modifying genes, including polycomb complex components. Chromatin remodeling is apparent in insulinoma-selective domains shared across patients, containing a specific set of regulatory sequences dominated by the SOX17 binding motif. Moreover, many of these regions are H3K27me3 repressed in β cells, suggesting that tumoral transition involves derepression of polycomb-targeted domains. Our work provides a compendium of aberrant cis-regulatory elements affecting the function and fate of β cells in their progression to insulinomas and a framework to identify coding and noncoding driver mutations.

Genomic and transcriptomic features of androgen receptor signaling inhibitor resistance in metastatic castration-resistant prostate cancer

BACKGROUNDAndrogen receptor signaling inhibitors (ARSIs) have improved outcomes for patients with metastatic castration-resistant prostate cancer (mCRPC), but their clinical benefit is limited by treatment resistance.METHODSTo investigate the mechanisms of ARSI resistance, we analyzed the whole-genome (n = 45) and transcriptome (n = 31) sequencing data generated from paired metastatic biopsies obtained before initiation of first-line ARSI therapy for mCRPC and after radiographic disease progression. We investigated the effects of genetic and pharmacologic modulation of SSTR1 in 22Rv1 cells, a representative mCRPC cell line.RESULTSWe confirmed the predominant role of tumor genetic alterations converging on augmenting androgen receptor (AR) signaling and the increased transcriptional heterogeneity and lineage plasticity during the emergence of ARSI resistance. We further identified amplifications involving a putative enhancer downstream of the AR and transcriptional downregulation of SSTR1, encoding somatostatin receptor 1, in ARSI-resistant tumors. We found that patients with SSTR1-low mCRPC tumors derived less benefit from subsequent ARSI therapy in a retrospective cohort. We showed that SSTR1 was antiproliferative in 22Rv1 cells and that the FDA-approved drug pasireotide suppressed 22Rv1 cell proliferation.CONCLUSIONOur findings expand the knowledge of ARSI resistance and point out actionable next steps, exemplified by potentially targeting SSTR1, to improve patient outcomes.FUNDINGNational Cancer Institute (NCI), NIH; Prostate Cancer Foundation; Conquer Cancer, American Society of Clinical Oncology Foundation; UCSF Benioff Initiative for Prostate Cancer Research; Netherlands Cancer Institute.

3D genomic analysis reveals novel enhancer-hijacking caused by complex structural alterations that drive oncogene overexpression

Cancer genomes are composed of many complex structural alterations on chromosomes and extrachromosomal DNA (ecDNA), making it difficult to identify non-coding enhancer regions that are hijacked to activate oncogene expression. Here, we describe a 3D genomics-based analysis called HAPI (Highly Active Promoter Interactions) to characterize enhancer hijacking. HAPI analysis of HiChIP data from 34 cancer cell lines identified enhancer hijacking events that activate both known and potentially novel oncogenes such as MYC, CCND1, ETV1, CRKL, and ID4. Furthermore, we found enhancer hijacking among multiple oncogenes from different chromosomes, often including MYC, on the same complex amplicons such as ecDNA. We characterized a MYC-ERBB2 chimeric ecDNA, in which ERBB2 heavily hijacks MYC's enhancers. Notably, CRISPRi of the MYC promoter led to increased interaction of ERBB2 with MYC enhancers and elevated ERBB2 expression. Our HAPI analysis tool provides a robust strategy to detect enhancer hijacking and reveals novel insights into oncogene activation.

Identifying dysregulated regions in amyotrophic lateral sclerosis through chromatin accessibility outliers

The high heritability of amyotrophic lateral sclerosis (ALS) contrasts with its low molecular diagnosis rate post-genetic testing, pointing to potential undiscovered genetic factors. To aid the exploration of these factors, we introduced EpiOut, an algorithm to identify chromatin accessibility outliers that are regions exhibiting divergent accessibility from the population baseline in a single or few samples. Annotation of accessible regions with histone chromatin immunoprecipitation sequencing and Hi-C indicates that outliers are concentrated in functional loci, especially among promoters interacting with active enhancers. Across different omics levels, outliers are robustly replicated, and chromatin accessibility outliers are reliable predictors of gene expression outliers and aberrant protein levels. When promoter accessibility does not align with gene expression, our results indicate that molecular aberrations are more likely to be linked to post-transcriptional regulation rather than transcriptional regulation. Our findings demonstrate that the outlier detection paradigm can uncover dysregulated regions in rare diseases. EpiOut is available at github.com/uci-cbcl/EpiOut.

An assessment of the genomic structural variation landscape in Sub-Saharan African populations

Structural variants are responsible for a large part of genomic variation between individuals and play a role in both common and rare diseases. Databases cataloguing structural variants notably do not represent the full spectrum of global diversity, particularly missing information from most African populations. To address this representation gap, we analysed 1,091 high-coverage African genomes, 545 of which are public data sets, and 546 which have been analysed for structural variants for the first time. Variants were called using five different tools and datasets merged and jointly called using SURVIVOR. We identified 67,795 structural variants throughout the genome, with 10,421 genes having at least one variant. Using a conservative overlap in merged data, 6,414 of the structural variants (9.5%) are novel compared to the Database of Genomic Variants. This study contributes to knowledge of the landscape of structural variant diversity in Africa and presents a reliable dataset for potential applications in population genetics and health-related research.

3D genomic analysis reveals novel enhancer-hijacking caused by complex structural alterations that drive oncogene overexpression

Cancer genomes are composed of many complex structural alterations on chromosomes and extrachromosomal DNA (ecDNA), making it difficult to identify non-coding enhancer regions that are hijacked to activate oncogene expression. Here, we describe a 3D genomics-based analysis called HAPI (Highly Active Promoter Interactions) to characterize enhancer hijacking. HAPI analysis of HiChIP data from 34 cancer cell lines identified enhancer hijacking events that activate both known and potentially novel oncogenes such as MYC, CCND1 , ETV1 , CRKL , and ID4 . Furthermore, we found enhancer hijacking among multiple oncogenes from different chromosomes, often including MYC , on the same complex amplicons such as ecDNA. We characterized a MYC - ERBB2 chimeric ecDNA, in which ERBB2 heavily hijacks MYC 's enhancers. Notably, CRISPRi of the MYC promoter led to increased interaction of ERBB2 with MYC enhancers and elevated ERBB2 expression. Our HAPI analysis tool provides a robust strategy to detect enhancer hijacking and reveals novel insights into oncogene activation.

Rare pathogenic structural variants show potential to enhance prostate cancer germline testing for African men

Prostate cancer (PCa) is highly heritable, with men of African ancestry at greatest risk and associated lethality. Lack of representation in genomic data means germline testing guidelines exclude for African men. Established that structural variations (SVs) are major contributors to human disease and prostate tumourigenesis, their role is under-appreciated in familial and therapeutic testing. Utilising a clinico-methodologically matched African (n = 113) versus European (n = 57) deep-sequenced PCa resource, we interrogated 42,966 high-quality germline SVs using a best-fit pathogenicity prediction workflow. We identified 15 potentially pathogenic SVs representing 12.4% African and 7.0% European patients, of which 72% and 86% met germline testing standard-of-care recommendations, respectively. Notable African-specific loss-of-function gene candidates include DNA damage repair MLH1 and BARD1 and tumour suppressors FOXP1, WASF1 and RB1. Representing only a fraction of the vast African diaspora, this study raises considerations with respect to the contribution of kilo-to-mega-base rare variants to PCa pathogenicity and African associated disparity.

Cytochrome P450 CitCYP97B modulates carotenoid accumulation diversity by hydroxylating β-cryptoxanthin in Citrus

Carotenoids in plant foods provide health benefits by functioning as provitamin A. One of the vital provitamin A carotenoids, β-cryptoxanthin, is typically plentiful in citrus fruit. However, little is known about the genetic basis of β-cryptoxanthin accumulation in citrus. Here, we performed a widely targeted metabolomic analysis of 65 major carotenoids and carotenoid derivatives to characterize carotenoid accumulation in Citrus and determine the taxonomic profile of β-cryptoxanthin. We used data from 81 newly sequenced representative accessions and 69 previously sequenced Citrus cultivars to reveal the genetic basis of β-cryptoxanthin accumulation through a genome-wide association study. We identified a causal gene, CitCYP97B, which encodes a cytochrome P450 protein whose substrate and metabolic pathways in land plants were undetermined. We subsequently demonstrated that CitCYP97B functions as a novel monooxygenase that specifically hydroxylates the β-ring of β-cryptoxanthin in a heterologous expression system. In planta experiments provided further evidence that CitCYP97B negatively regulates β-cryptoxanthin content. Using the sequenced Citrus accessions, we found that two critical structural cis-element variations contribute to increased expression of CitCYP97B, thereby altering β-cryptoxanthin accumulation in fruit. Hybridization/introgression appear to have contributed to the prevalence of two cis-element variations in different Citrus types during citrus evolution. Overall, these findings extend our understanding of the regulation and diversity of carotenoid metabolism in fruit crops and provide a genetic target for production of β-cryptoxanthin-biofortified products.

Constitutional and acquired genetic variants in ARID5B in pediatric B-cell precursor acute lymphoblastic leukemia

Constitutional polymorphisms in ARID5B are associated with an increased risk of developing high hyperdiploid (HeH; 51-67 chromosomes) pediatric B-cell precursor acute lymphoblastic leukemia (BCP ALL). Here, we investigated constitutional and somatic ARID5B variants in 1335 BCP ALL cases from five different cohorts, with a particular focus on HeH cases. In 353 HeH ALL that were heterozygous for risk alleles and trisomic for chromosome 10, where ARID5B is located, a significantly higher proportion of risk allele duplication was seen for the SNPs rs7090445 (p = 0.009), rs7089424 (p = 0.005), rs7073837 (p = 0.03), and rs10740055 (p = 0.04). Somatic ARID5B deletions were seen in 16/1335 cases (1.2%), being more common in HeH than in other genetic subtypes (2.2% vs. 0.4%; p = 0.002). The expression of ARID5B in HeH cases with genomic deletions was reduced, consistent with a functional role in leukemogenesis. Whole-genome sequencing and RNA-sequencing in HeH revealed additional somatic events involving ARID5B, resulting in a total frequency of 3.6% of HeH cases displaying a somatic ARID5B aberration. Overall, our results show that both constitutional and somatic events in ARID5B are involved in the leukemogenesis of pediatric BCP ALL, particularly in the HeH subtype.

SNP-SVant: A Computational Workflow to Predict and Annotate Genomic Variants in Organisms Lacking Benchmarked Variants

Whole-genome sequencing is widely used to investigate population genomic variation in organisms of interest. Assorted tools have been independently developed to call variants from short-read sequencing data aligned to a reference genome, including single nucleotide polymorphisms (SNPs) and structural variations (SVs). We developed SNP-SVant, an integrated, flexible, and computationally efficient bioinformatic workflow that predicts high-confidence SNPs and SVs in organisms without benchmarked variants, which are traditionally used for distinguishing sequencing errors from real variants. In the absence of these benchmarked datasets, we leverage multiple rounds of statistical recalibration to increase the precision of variant prediction. The SNP-SVant workflow is flexible, with user options to tradeoff accuracy for sensitivity. The workflow predicts SNPs and small insertions and deletions using the Genome Analysis ToolKit (GATK) and predicts SVs using the Genome Rearrangement IDentification Software Suite (GRIDSS), and it culminates in variant annotation using custom scripts. A key utility of SNP-SVant is its scalability. Variant calling is a computationally expensive procedure, and thus, SNP-SVant uses a workflow management system with intermediary checkpoint steps to ensure efficient use of resources by minimizing redundant computations and omitting steps where dependent files are available. SNP-SVant also provides metrics to assess the quality of called variants and converts between VCF and aligned FASTA format outputs to ensure compatibility with downstream tools to calculate selection statistics, which are commonplace in population genomics studies. By accounting for both small and large structural variants, users of this workflow can obtain a wide-ranging view of genomic alterations in an organism of interest. Overall, this workflow advances our capabilities in assessing the functional consequences of different types of genomic alterations, ultimately improving our ability to associate genotypes with phenotypes. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Predicting single nucleotide polymorphisms and structural variations Support Protocol 1: Downloading publicly available sequencing data Support Protocol 2: Visualizing variant loci using Integrated Genome Viewer Support Protocol 3: Converting between VCF and aligned FASTA formats.

Comparison of structural variant callers for massive whole-genome sequence data

BACKGROUND

Detecting structural variations (SVs) at the population level using next-generation sequencing (NGS) requires substantial computational resources and processing time. Here, we compared the performances of 11 SV callers: Delly, Manta, GridSS, Wham, Sniffles, Lumpy, SvABA, Canvas, CNVnator, MELT, and INSurVeyor. These SV callers have been recently published and have been widely employed for processing massive whole-genome sequencing datasets. We evaluated the accuracy, sequence depth, running time, and memory usage of the SV callers.

RESULTS

Notably, several callers exhibited better calling performance for deletions than for duplications, inversions, and insertions. Among the SV callers, Manta identified deletion SVs with better performance and efficient computing resources, and both Manta and MELT demonstrated relatively good precision regarding calling insertions. We confirmed that the copy number variation callers, Canvas and CNVnator, exhibited better performance in identifying long duplications as they employ the read-depth approach. Finally, we also verified the genotypes inferred from each SV caller using a phased long-read assembly dataset, and Manta showed the highest concordance in terms of the deletions and insertions.

CONCLUSIONS

Our findings provide a comprehensive understanding of the accuracy and computational efficiency of SV callers, thereby facilitating integrative analysis of SV profiles in diverse large-scale genomic datasets.

Automated recognition of chromosome fusion using an alignment-free natural vector method

Chromosomal fusion is a significant form of structural variation, but research into algorithms for its identification has been limited. Most existing methods rely on synteny analysis, which necessitates manual annotations and always involves inefficient sequence alignments. In this paper, we present a novel alignment-free algorithm for chromosomal fusion recognition. Our method transforms the problem into a series of assignment problems using natural vectors and efficiently solves them with the Kuhn-Munkres algorithm. When applied to the human/gorilla and swamp buffalo/river buffalo datasets, our algorithm successfully and efficiently identifies chromosomal fusion events. Notably, our approach offers several advantages, including higher processing speeds by eliminating time-consuming alignments and removing the need for manual annotations. By an alignment-free perspective, our algorithm initially considers entire chromosomes instead of fragments to identify chromosomal structural variations, offering substantial potential to advance research in this field.

Chemotherapy induces myeloid-driven spatial T-cell exhaustion in ovarian cancer

To uncover the intricate, chemotherapy-induced spatiotemporal remodeling of the tumor microenvironment, we conducted integrative spatial and molecular characterization of 97 high-grade serous ovarian cancer (HGSC) samples collected before and after chemotherapy. Using single-cell and spatial analyses, we identify increasingly versatile immune cell states, which form spatiotemporally dynamic microcommunities at the tumor-stroma interface. We demonstrate that chemotherapy triggers spatial redistribution and exhaustion of CD8+ T cells due to prolonged antigen presentation by macrophages, both within interconnected myeloid networks termed "Myelonets" and at the tumor stroma interface. Single-cell and spatial transcriptomics identifies prominent TIGIT-NECTIN2 ligand-receptor interactions induced by chemotherapy. Using a functional patient-derived immuno-oncology platform, we show that CD8+T-cell activity can be boosted by combining immune checkpoint blockade with chemotherapy. Our discovery of chemotherapy-induced myeloid-driven spatial T-cell exhaustion paves the way for novel immunotherapeutic strategies to unleash CD8+ T-cell-mediated anti-tumor immunity in HGSC.

Detecting Somatic Mutations Without Matched Normal Samples Using Long Reads

DNA sequencing of tumours to identify somatic mutations has become a critical tool to guide the type of treatment given to cancer patients. The gold standard for mutation calling is comparing sequencing data from the tumour to a matched normal sample to avoid mis-classifying inherited SNPs as mutations. This procedure works extremely well, but in certain situations only a tumour sample is available. While approaches have been developed to find mutations without a matched normal, they have limited accuracy or require specific types of input data (e.g. ultra-deep sequencing). Here we explore the application of single molecule long read sequencing to calling somatic mutations without matched normal samples. We develop a simple theoretical framework to show how haplotype phasing is an important source of information for determining whether a variant is a somatic mutation. We then use simulations to assess the range of experimental parameters (tumour purity, sequencing depth) where this approach is effective. These ideas are developed into a prototype somatic mutation caller, smrest, and its use is demonstrated on two highly mutated cancer cell lines. Finally, we argue that this approach has potential to measure clinically important biomarkers that are based on the genome-wide distribution of mutations: tumour mutation burden and mutation signatures.

Osteosarcoma PDX-Derived Cell Line Models for Preclinical Drug Evaluation Demonstrate Metastasis Inhibition by Dinaciclib through a Genome-Targeted Approach

PURPOSE

Models to study metastatic disease in rare cancers are needed to advance preclinical therapeutics and to gain insight into disease biology. Osteosarcoma is a rare cancer with a complex genomic landscape in which outcomes for patients with metastatic disease are poor. As osteosarcoma genomes are highly heterogeneous, multiple models are needed to fully elucidate key aspects of disease biology and to recapitulate clinically relevant phenotypes.

EXPERIMENTAL DESIGN

Matched patient samples, patient-derived xenografts (PDX), and PDX-derived cell lines were comprehensively evaluated using whole-genome sequencing and RNA sequencing. The in vivo metastatic phenotype of the PDX-derived cell lines was characterized in both an intravenous and an orthotopic murine model. As a proof-of-concept study, we tested the preclinical effectiveness of a cyclin-dependent kinase inhibitor on the growth of metastatic tumors in an orthotopic amputation model.

RESULTS

PDXs and PDX-derived cell lines largely maintained the expression profiles of the patient from which they were derived despite the emergence of whole-genome duplication in a subset of cell lines. The cell lines were heterogeneous in their metastatic capacity, and heterogeneous tissue tropism was observed in both intravenous and orthotopic models. Single-agent dinaciclib was effective at dramatically reducing the metastatic burden.

CONCLUSIONS

The variation in metastasis predilection sites between osteosarcoma PDX-derived cell lines demonstrates their ability to recapitulate the spectrum of the disease observed in patients. We describe here a panel of new osteosarcoma PDX-derived cell lines that we believe will be of wide use to the osteosarcoma research community.

Haplotype-specific assembly of shattered chromosomes in esophageal adenocarcinomas

The epigenetic landscape of cancer is regulated by many factors, but primarily it derives from the underlying genome sequence. Chromothripsis is a catastrophic localized genome shattering event that drives, and often initiates, cancer evolution. We characterized five esophageal adenocarcinoma organoids with chromothripsis using long-read sequencing and transcriptome and epigenome profiling. Complex structural variation and subclonal variants meant that haplotype-aware de novo methods were required to generate contiguous cancer genome assemblies. Chromosomes were assembled separately and scaffolded using haplotype-resolved Hi-C reads, producing accurate assemblies even with up to 900 structural rearrangements. There were widespread differences between the chromothriptic and wild-type copies of chromosomes in topologically associated domains, chromatin accessibility, histone modifications, and gene expression. Differential epigenome peaks were most enriched within 10 kb of chromothriptic structural variants. Alterations in transcriptome and higher-order chromosome organization frequently occurred near differential epigenetic marks. Overall, chromothripsis reshapes gene regulation, causing coordinated changes in epigenetic landscape, transcription, and chromosome conformation.

Genomic and epigenomic basis of breast invasive lobular carcinomas lacking CDH1 genetic alterations

CDH1 (E-cadherin) bi-allelic inactivation is the hallmark alteration of breast invasive lobular carcinoma (ILC), resulting in its discohesive phenotype. A subset of ILCs, however, lack CDH1 genetic/epigenetic inactivation, and their genetic underpinning is unknown. Through clinical targeted sequencing data reanalysis of 364 primary ILCs, we identified 25 ILCs lacking CDH1 bi-allelic genetic alterations. CDH1 promoter methylation was frequent (63%) in these cases. Targeted sequencing reanalysis revealed 3 ILCs harboring AXIN2 deleterious fusions (n = 2) or loss-of-function mutation (n = 1). Whole-genome sequencing of 3 cases lacking bi-allelic CDH1 genetic/epigenetic inactivation confirmed the AXIN2 mutation and no other cell-cell adhesion genetic alterations but revealed a new CTNND1 (p120) deleterious fusion. AXIN2 knock-out in MCF7 cells resulted in lobular-like features, including increased cellular migration and resistance to anoikis. Taken together, ILCs lacking CDH1 genetic/epigenetic alterations are driven by inactivating alterations in other cell adhesion genes (CTNND1 or AXIN2), endorsing a convergent phenotype in ILC.

Genomic evolution and insights into agronomic trait innovations of Sesamum species

Sesame is an ancient oilseed crop with high oil content and quality. However, the evolutionary history and genetic mechanisms of its valuable agronomic traits remain unclear. Here, we report chromosome-scale genomes of cultivated sesame (Sesamum indicum L.) and six wild Sesamum species, representing all three karyotypes within this genus. Karyotyping and genome-based phylogenic analysis revealed the evolutionary route of Sesamum species from n = 13 to n = 16 and revealed that allotetraploidization occurred in the wild species Sesamum radiatum. Early divergence of the Sesamum genus (48.5-19.7 million years ago) during the Tertiary period and its ancient phylogenic position within eudicots were observed. Pan-genome analysis revealed 9164 core gene families in the 7 Sesamum species. These families are significantly enriched in various metabolic pathways, including fatty acid (FA) metabolism and FA biosynthesis. Structural variations in SiPT1 and SiDT1 within the phosphatidyl ethanolamine-binding protein gene family lead to the genomic evolution of plant-architecture and inflorescence-development phenotypes in Sesamum. A genome-wide association study (GWAS) of an interspecific population and genome comparisons revealed a long terminal repeat insertion and a sequence deletion in DIR genes of wild Sesamum angustifolium and cultivated sesame, respectively; both variations independently cause high susceptibility to Fusarium wilt disease. A GWAS of 560 sesame accessions combined with an overexpression study confirmed that the NAC1 and PPO genes play an important role in upregulating oil content of sesame. Our study provides high-quality genomic resources for cultivated and wild Sesamum species and insights that can improve molecular breeding strategies for sesame and other oilseed crops.

Deciphering cancer genomes with GenomeSpy: a grammar-based visualization toolkit

BACKGROUND

Visualization is an indispensable facet of genomic data analysis. Despite the abundance of specialized visualization tools, there remains a distinct need for tailored solutions. However, their implementation typically requires extensive programming expertise from bioinformaticians and software developers, especially when building interactive applications. Toolkits based on visualization grammars offer a more accessible, declarative way to author new visualizations. Yet, current grammar-based solutions fall short in adequately supporting the interactive analysis of large datasets with extensive sample collections, a pivotal task often encountered in cancer research.

FINDINGS

We present GenomeSpy, a grammar-based toolkit for authoring tailored, interactive visualizations for genomic data analysis. By using combinatorial building blocks and a declarative language, users can implement new visualization designs easily and embed them in web pages or end-user-oriented applications. A distinctive element of GenomeSpy's architecture is its effective use of the graphics processing unit in all rendering, enabling a high frame rate and smoothly animated interactions, such as navigation within a genome. We demonstrate the utility of GenomeSpy by characterizing the genomic landscape of 753 ovarian cancer samples from patients in the DECIDER clinical trial. Our results expand the understanding of the genomic architecture in ovarian cancer, particularly the diversity of chromosomal instability.

CONCLUSIONS

GenomeSpy is a visualization toolkit applicable to a wide range of tasks pertinent to genome analysis. It offers high flexibility and exceptional performance in interactive analysis. The toolkit is open source with an MIT license, implemented in JavaScript, and available at https://genomespy.app/.

Recent gene selection and drug resistance underscore clinical adaptation across Candida species

Understanding how microbial pathogens adapt to treatments, humans and clinical environments is key to infer mechanisms of virulence, transmission and drug resistance. This may help improve therapies and diagnostics for infections with a poor prognosis, such as those caused by fungal pathogens, including Candida. Here we analysed genomic variants across approximately 2,000 isolates from six Candida species (C. glabrata, C. auris, C. albicans, C. tropicalis, C. parapsilosis and C. orthopsilosis) and identified genes under recent selection, suggesting a highly complex clinical adaptation. These involve species-specific and convergently affected adaptive mechanisms, such as adhesion. Using convergence-based genome-wide association studies we identified known drivers of drug resistance alongside potentially novel players. Finally, our analyses reveal an important role of structural variants and suggest an unexpected involvement of (para)sexual recombination in the spread of resistance. Our results provide insights on how opportunistic pathogens adapt to human-related environments and unearth candidate genes that deserve future attention.

HRD related signature 3 predicts clinical outcome in advanced tubo-ovarian high-grade serous carcinoma

OBJECTIVES

We evaluated usability of single base substitution signature 3 (Sig3) as a biomarker for homologous recombination deficiency (HRD) in tubo-ovarian high-grade serous carcinoma (HGSC).

MATERIALS AND METHODS

This prospective observational trial includes 165 patients with advanced HGSC. Fresh tissue samples (n = 456) from multiple intra-abdominal areas at diagnosis and after neoadjuvant chemotherapy (NACT) were collected for whole-genome sequencing. Sig3 was assessed by fitting samples independently with COSMIC v3.2 reference signatures. An HR scar assay was applied for comparison. Progression-free survival (PFS) and overall survival (OS) were studied using Kaplan-Meier and Cox regression analysis.

RESULTS

Sig3 has a bimodal distribution, eliminating the need for an arbitrary cutoff typical in HR scar tests. Sig3 could be assessed from samples with low (10%) cancer cell proportion and was consistent between multiple samples and stable during NACT. At diagnosis, 74 (45%) patients were HRD (Sig3+), while 91 (55%) were HR proficient (HRP, Sig3-). Sig3+ patients had longer PFS and OS than Sig3- patients (22 vs. 13 months and 51 vs. 34 months respectively, both p < 0.001). Sig3 successfully distinguished the poor prognostic HRP group among BRCAwt patients (PFS 19 months for Sig3+ and 13 months for Sig3- patients, p < 0.001). However, Sig3 at diagnosis did not predict chemoresponse anymore in the first relapse. The patient-level concordance between Sig3 and HR scar assay was 87%, and patients with HRD according to both tests had the longest median PFS.

CONCLUSIONS

Sig3 is a prognostic marker in advanced HGSC and useful tool in patient stratification for HRD.

Unified framework for patient-derived, tumor-organoid-based predictive testing of standard-of-care therapies in metastatic colorectal cancer

Predictive drug testing of patient-derived tumor organoids (PDTOs) holds promise for personalizing treatment of metastatic colorectal cancer (mCRC), but prospective data are limited to chemotherapy regimens with conflicting results. We describe a unified framework for PDTO-based predictive testing across standard-of-care chemotherapy and biologic and targeted therapy options. In an Australian community cohort, PDTO predictions based on treatment-naive patients (n = 56) and response rates from first-line mCRC clinical trials achieve 83% accuracy for forecasting responses in patients receiving palliative treatments (18 patients, 29 treatments). Similar assay accuracy is achieved in a prospective study of third-line or later mCRC treatment, AGITG FORECAST-1 (n = 30 patients). "Resistant" predictions are associated with inferior progression-free survival; misclassification rates are similar by regimen. Liver metastases are the optimal site for sampling, with testing achievable within 7 weeks for 68.8% cases. Our findings indicate that PDTO drug panel testing can provide predictive information for multifarious standard-of-care therapies for mCRC.

Manipulating the 3D organization of the largest synthetic yeast chromosome

Whether synthetic genomes can power life has attracted broad interest in the synthetic biology field. Here, we report de novo synthesis of the largest eukaryotic chromosome thus far, synIV, a 1,454,621-bp yeast chromosome resulting from extensive genome streamlining and modification. We developed megachunk assembly combined with a hierarchical integration strategy, which significantly increased the accuracy and flexibility of synthetic chromosome construction. Besides the drastic sequence changes, we further manipulated the 3D structure of synIV to explore spatial gene regulation. Surprisingly, we found few gene expression changes, suggesting that positioning inside the yeast nucleoplasm plays a minor role in gene regulation. Lastly, we tethered synIV to the inner nuclear membrane via its hundreds of loxPsym sites and observed transcriptional repression of the entire chromosome, demonstrating chromosome-wide transcription manipulation without changing the DNA sequences. Our manipulation of the spatial structure of synIV sheds light on higher-order architectural design of the synthetic genomes.

The microtubule inhibitor eribulin demonstrates efficacy in platinum-resistant and refractory high-grade serous ovarian cancer patient-derived xenograft models

Background

Despite initial response to platinum-based chemotherapy and PARP inhibitor therapy (PARPi), nearly all recurrent high-grade serous ovarian cancer (HGSC) will acquire lethal drug resistance; indeed, ~15% of individuals have de novo platinum-refractory disease.

Objectives

To determine the potential of anti-microtubule agent (AMA) therapy (paclitaxel, vinorelbine and eribulin) in platinum-resistant or refractory (PRR) HGSC by assessing response in patient-derived xenograft (PDX) models of HGSC.

Design and methods

Of 13 PRR HGSC PDX, six were primary PRR, derived from chemotherapy-naïve samples (one was BRCA2 mutant) and seven were from samples obtained following chemotherapy treatment in the clinic (five were mutant for either BRCA1 or BRCA2 (BRCA1/2), four with prior PARPi exposure), recapitulating the population of individuals with aggressive treatment-resistant HGSC in the clinic. Molecular analyses and in vivo treatment studies were undertaken.

Results

Seven out of thirteen PRR PDX (54%) were sensitive to treatment with the AMA, eribulin (time to progressive disease (PD) ⩾100 days from the start of treatment) and 11 out of 13 PDX (85%) derived significant benefit from eribulin [time to harvest (TTH) for each PDX with p < 0.002]. In 5 out of 10 platinum-refractory HGSC PDX (50%) and one out of three platinum-resistant PDX (33%), eribulin was more efficacious than was cisplatin, with longer time to PD and significantly extended TTH (each PDX p < 0.02). Furthermore, four of these models were extremely sensitive to all three AMA tested, maintaining response until the end of the experiment (120d post-treatment start). Despite harbouring secondary BRCA2 mutations, two BRCA2-mutant PDX models derived from heavily pre-treated individuals were sensitive to AMA. PRR HGSC PDX models showing greater sensitivity to AMA had high proliferative indices and oncogene expression. Two PDX models, both with prior chemotherapy and/or PARPi exposure, were refractory to all AMA, one of which harboured the SLC25A40-ABCB1 fusion, known to upregulate drug efflux via MDR1.

Conclusion

The efficacy observed for eribulin in PRR HGSC PDX was similar to that observed for paclitaxel, which transformed ovarian cancer clinical practice. Eribulin is therefore worthy of further consideration in clinical trials, particularly in ovarian carcinoma with early failure of carboplatin/paclitaxel chemotherapy.

Context-dependent neocentromere activity in synthetic yeast chromosome VIII

Pioneering advances in genome engineering, and specifically in genome writing, have revolutionized the field of synthetic biology, propelling us toward the creation of synthetic genomes. The Sc2.0 project aims to build the first fully synthetic eukaryotic organism by assembling the genome of Saccharomyces cerevisiae. With the completion of synthetic chromosome VIII (synVIII) described here, this goal is within reach. In addition to writing the yeast genome, we sought to manipulate an essential functional element: the point centromere. By relocating the native centromere sequence to various positions along chromosome VIII, we discovered that the minimal 118-bp CEN8 sequence is insufficient for conferring chromosomal stability at ectopic locations. Expanding the transplanted sequence to include a small segment (∼500 bp) of the CDEIII-proximal pericentromere improved chromosome stability, demonstrating that minimal centromeres display context-dependent functionality.

Prolonged survival of a patient with active MDR-TB HIV co-morbidity: insights from a Mycobacterium tuberculosis strain with a unique genomic deletion

Coinfection of HIV and multidrug-resistant tuberculosis (MDR-TB) presents significant challenges in terms of the treatment and prognosis of tuberculosis, leading to complexities in managing the disease and impacting the overall outcome for TB patients. This study presents a remarkable case of a patient with MDR-TB and HIV coinfection who survived for over 8 years, despite poor treatment adherence and comorbidities. Whole genome sequencing (WGS) of the infecting Mycobacterium tuberculosis (Mtb) strain revealed a unique genomic deletion, spanning 18 genes, including key genes involved in hypoxia response, intracellular survival, immunodominant antigens, and dormancy. This deletion, that we have called "Del-X," potentially exerts a profound influence on the bacterial physiology and its virulence. Only few similar deletions were detected in other non-related Mtb genomes worldwide. In vivo evolution analysis identified drug resistance and metabolic adaptation mutations and their temporal dynamics during the patient's treatment course.

Transcriptomic and genomic effects of gamma-radiation exposure on strains of the black yeast Exophiala dermatitidis evolved to display increased ionizing radiation resistance

IMPORTANCE

Ionizing radiation poses a significant threat to living organisms and human health, given its destructive nature and widespread use in fields such as medicine and the potential for nuclear disasters. Melanized fungi exhibit remarkable survival capabilities, enduring doses up to 1,000-fold higher than mammals. Through adaptive laboratory evolution, we validated the protective role of constitutive upregulation of DNA repair genes in the black yeast Exophiala dermatitidis, enhancing survival after radiation exposure. Surprisingly, we found that evolved strains lacking melanin still achieved high levels of radioresistance. Our study unveiled the significance of robust activation and enhancement of redox homeostasis, as evidenced by the profound transcriptional changes and increased accumulation of mutations, in substantially improving ionizing radiation resistance in the absence of melanin. These findings underscore the delicate balance between DNA repair and redox homeostasis for an organism's ability to endure and recover from radiation exposure.

A 4.43-Kb deletion of chromosomal segment containing an ovate family protein confers long capsule in sesame (Sesamum indicum L.)

KEY MESSAGE

A 4.43-Kb structural variation in the sesame genome results in the deletion of the Siofp1 gene and induces the long capsule length trait. Capsule length (CL) has a positive effect on seed weight and yield in various agronomically important species; however, the molecular mechanism underlying long capsule trait regulation in sesame remains unknown. The inheritance analysis showed that long capsule traits (CL > 4.0 cm) were dominant over normal length (average CL = 3.0 cm) and were controlled by a single gene pair. Association mapping with a RIL population and 259 natural sesame germplasm accessions indicated that the target interval was 52,830-730,961 bp of SiChr.10 in sesame. Meanwhile, the structural variation (SV) of the association mapping revealed that only SV_414325 on chromosome 10 was significantly associated with the CL trait, with a P value of 1.1135E-19. SV_414325 represents a 4430-bp deletion from 414,325 to 418,756 bp on SiChr.10, covering Sindi_2155000 (named SiOFP1). In the normal length type, Siofp1 encodes 411 amino acids of the ovate family proteins and is highly expressed in the leaf, stem, bud, and capsule tissues of sesame. In accordance with the transcriptional repressor character, Siofp1 overexpression in transgenic Arabidopsis (T0 and T1 generations) induced a 25-39% greater shortening of silique length than the wild type (P < 0.05), as well as round cauline leaves and short carpels. These results confirm that SiOFP1 plays a key role in regulating CL trait in sesame and other flowering plants. These findings provide a theoretical and material basis for sesame capsule development and high-yield breeding research.

H&E image analysis pipeline for quantifying morphological features

Detecting cell types from histopathological images is essential for various digital pathology applications. However, large number of cells in whole-slide images (WSIs) necessitates automated analysis pipelines for efficient cell type detection. Herein, we present hematoxylin and eosin (H&E) Image Processing pipeline (HEIP) for automatied analysis of scanned H&E-stained slides. HEIP is a flexible and modular open-source software that performs preprocessing, instance segmentation, and nuclei feature extraction. To evaluate the performance of HEIP, we applied it to extract cell types from ovarian high-grade serous carcinoma (HGSC) patient WSIs. HEIP showed high precision in instance segmentation, particularly for neoplastic and epithelial cells. We also show that there is a significant correlation between genomic ploidy values and morphological features, such as major axis of the nucleus.

Comprehensive single-cell genome analysis at nucleotide resolution using the PTA Analysis Toolbox

Detection of somatic mutations in single cells has been severely hampered by technical limitations of whole-genome amplification. Novel technologies including primary template-directed amplification (PTA) significantly improved the accuracy of single-cell whole-genome sequencing (WGS) but still generate hundreds of artifacts per amplification reaction. We developed a comprehensive bioinformatic workflow, called the PTA Analysis Toolbox (PTATO), to accurately detect single base substitutions, insertions-deletions (indels), and structural variants in PTA-based WGS data. PTATO includes a machine learning approach and filtering based on recurrence to distinguish PTA artifacts from true mutations with high sensitivity (up to 90%), outperforming existing bioinformatic approaches. Using PTATO, we demonstrate that hematopoietic stem cells of patients with Fanconi anemia, which cannot be analyzed using regular WGS, have normal somatic single base substitution burdens but increased numbers of deletions. Our results show that PTATO enables studying somatic mutagenesis in the genomes of single cells with unprecedented sensitivity and accuracy.

A Matched Molecular and Clinical Analysis of the Epithelioid Haemangioendothelioma Cohort in the Stafford Fox Rare Cancer Program and Contextual Literature Review

BACKGROUND

Epithelioid haemangioendothelioma (EHE) is an ultra-rare malignant vascular tumour with a prevalence of 1 per 1,000,000. It is typically molecularly characterised by a WWTR1::CAMTA1 gene fusion in approximately 90% of cases, or a YAP1::TFE3 gene fusion in approximately 10% of cases. EHE cases are typically refractory to therapies, and no anticancer agents are reimbursed for EHE in Australia.

METHODS

We report a cohort of nine EHE cases with comprehensive histologic and molecular profiling from the Walter and Eliza Hall Institute of Medical Research Stafford Fox Rare Cancer Program (WEHI-SFRCP) collated via nation-wide referral to the Australian Rare Cancer (ARC) Portal. The diagnoses of EHE were confirmed by histopathological and immunohistochemical (IHC) examination. Molecular profiling was performed using the TruSight Oncology 500 assay, the TruSight RNA fusion panel, whole genome sequencing (WGS), or whole exome sequencing (WES).

RESULTS

Molecular analysis of RNA, DNA or both was possible in seven of nine cases. The WWTR1::CAMTA1 fusion was identified in five cases. The YAP1::TFE3 fusion was identified in one case, demonstrating unique morphology compared to cases with the more common WWTR1::CAMTA1 fusion. All tumours expressed typical endothelial markers CD31, ERG, and CD34 and were negative for pan-cytokeratin. Cases with a WWTR1::CAMTA1 fusion displayed high expression of CAMTA1 and the single case with a YAP1::TFE3 fusion displayed high expression of TFE3. Survival was highly variable and unrelated to molecular profile.

CONCLUSIONS

This cohort of EHE cases provides molecular and histopathological characterisation and matching clinical information that emphasises the molecular patterns and variable clinical outcomes and adds to our knowledge of this ultra-rare cancer. Such information from multiple studies will advance our understanding, potentially improving treatment options.

Barcoded multiple displacement amplification for high coverage sequencing in spatial genomics

Determining mutational landscapes in a spatial context is essential for understanding genetically heterogeneous cell microniches. Current approaches, such as Multiple Displacement Amplification (MDA), offer high genome coverage but limited multiplexing, which hinders large-scale spatial genomic studies. Here, we introduce barcoded MDA (bMDA), a technique that achieves high-coverage genomic analysis of low-input DNA while enhancing the multiplexing capabilities. By incorporating cell barcodes during MDA, bMDA streamlines library preparation in one pot, thereby overcoming a key bottleneck in spatial genomics. We apply bMDA to the integrative spatial analysis of triple-negative breast cancer tissues by examining copy number alterations, single nucleotide variations, structural variations, and kataegis signatures for each spatial microniche. This enables the assessment of subclonal evolutionary relationships within a spatial context. Therefore, bMDA has emerged as a scalable technology with the potential to advance the field of spatial genomics significantly.

Spontaneous mutagenesis in human cells is controlled by REV1-Polymerase ζ and PRIMPOL

Translesion DNA synthesis (TLS) facilitates replication over damaged or difficult-to-replicate templates by employing specialized DNA polymerases. We investigate the effect on spontaneous mutagenesis of three main TLS control mechanisms: REV1 and PCNA ubiquitylation that recruit TLS polymerases and PRIMPOL that creates post-replicative gaps. Using whole-genome sequencing of cultured human RPE-1 cell clones, we find that REV1 and Polymerase ζ are wholly responsible for one component of base substitution mutagenesis that resembles homologous recombination deficiency, whereas the remaining component that approximates oxidative mutagenesis is reduced in PRIMPOL-/- cells. Small deletions in short repeats appear in REV1-/-PCNAK164R/K164R double mutants, revealing an alternative TLS mechanism. Also, 500-5,000 bp deletions appear in REV1-/- and REV3L-/- mutants, and chromosomal instability is detectable in REV1-/-PRIMPOL-/- cells. Our results indicate that TLS protects the genome from deletions and large rearrangements at the expense of being responsible for the majority of spontaneous base substitutions.

Genomic landscape of TP53-mutated myeloid malignancies

TP53-mutated myeloid malignancies are associated with complex cytogenetics and extensive structural variants, which complicates detailed genomic analysis by conventional clinical techniques. We performed whole-genome sequencing (WGS) of 42 acute myeloid leukemia (AML)/myelodysplastic syndromes (MDS) cases with paired normal tissue to better characterize the genomic landscape of TP53-mutated AML/MDS. WGS accurately determines TP53 allele status, a key prognostic factor, resulting in the reclassification of 12% of cases from monoallelic to multihit. Although aneuploidy and chromothripsis are shared with most TP53-mutated cancers, the specific chromosome abnormalities are distinct to each cancer type, suggesting a dependence on the tissue of origin. ETV6 expression is reduced in nearly all cases of TP53-mutated AML/MDS, either through gene deletion or presumed epigenetic silencing. Within the AML cohort, mutations of NF1 are highly enriched, with deletions of 1 copy of NF1 present in 45% of cases and biallelic mutations in 17%. Telomere content is increased in TP53-mutated AMLs compared with other AML subtypes, and abnormal telomeric sequences were detected in the interstitial regions of chromosomes. These data highlight the unique features of TP53-mutated myeloid malignancies, including the high frequency of chromothripsis and structural variation, the frequent involvement of unique genes (including NF1 and ETV6) as cooperating events, and evidence for altered telomere maintenance.

Precise characterization of somatic complex structural variations from tumor/control paired long-read sequencing data with nanomonsv

We present our novel software, nanomonsv, for detecting somatic structural variations (SVs) using tumor and matched control long-read sequencing data with a single-base resolution. The current version of nanomonsv includes two detection modules, Canonical SV module, and Single breakend SV module. Using tumor/control paired long-read sequencing data from three cancer and their matched lymphoblastoid lines, we demonstrate that Canonical SV module can identify somatic SVs that can be captured by short-read technologies with higher precision and recall than existing methods. In addition, we have developed a workflow to classify mobile element insertions while elucidating their in-depth properties, such as 5' truncations, internal inversions, as well as source sites for 3' transductions. Furthermore, Single breakend SV module enables the detection of complex SVs that can only be identified by long-reads, such as SVs involving highly-repetitive centromeric sequences, and LINE1- and virus-mediated rearrangements. In summary, our approaches applied to cancer long-read sequencing data can reveal various features of somatic SVs and will lead to a better understanding of mutational processes and functional consequences of somatic SVs.

Whole Animal Genome Sequencing: user-friendly, rapid, containerized pipelines for processing, variant discovery, and annotation of short-read whole genome sequencing data

Advancements in massively parallel short-read sequencing technologies and the associated decreasing costs have led to large and diverse variant discovery efforts across species. However, processing high-throughput short-read sequencing data can be challenging with potential pitfalls and bioinformatics bottlenecks in generating reproducible results. Although a number of pipelines exist that address these challenges, these are often geared toward human or traditional model organism species and can be difficult to configure across institutions. Whole Animal Genome Sequencing (WAGS) is an open-source set of user-friendly, containerized pipelines designed to simplify the process of identifying germline short (SNP and indel) and structural variants (SVs) geared toward the veterinary community but adaptable to any species with a suitable reference genome. We present a description of the pipelines [adapted from the best practices of the Genome Analysis Toolkit (GATK)], along with benchmarking data from both the preprocessing and joint genotyping steps, consistent with a typical user workflow.

Characteristics of DNA macro-alterations in breast cancer with liver metastasis before treatment

BACKGROUND

Whole-genome doubling (WGD) has been observed in 30% of cancers, followed by a highly complex rearranged karyotype unfavourable to breast cancer's outcome. However, the macro-alterations that characterise liver metastasis in breast cancer(BC) are poorly understood. Here, we conducted a whole-genome sequencing analysis of liver metastases to explore the status and the time frame model of these macro-alterations in pre-treatment patients with metastatic breast cancer.

RESULTS

Whole-genome sequencing was conducted in 11 paired primary tumours, lymph node metastasis, and liver metastasis fresh samples from four patients with late-stage breast cancer. We also chose five postoperative frozen specimens from patients with early-stage breast cancer before any treatment as control. Surprisingly, all four liver metastasis samples were classified as WGD + . However, the previous study reported that WGD happened in 30% of cancers and 2/5 in our early-stage samples. WGD was not observed in the two separate primary tumours and one lymph node metastasis of one patient with metastatic BC, but her liver metastasis showed an early burst of bi-allelic copy number gain. The phylogenetic tree proves her 4 tumour samples were the polyclonal origin and only one WGD + clone metastasis to the liver. Another 3 metastatic BC patients' primary tumour and lymph node metastasis experienced WGD as well as liver metastasis, and they all showed similar molecular time-frame of copy number(CN) gain across locations within the same patient. These patients' tumours were of monoclonal origin, and WGD happened in a founding clone before metastasis, explaining that all samples share the CN-gain time frame. After WGD, the genomes usually face instability to evolve other macro-alterations. For example, a greater quantity and variety of complex structural variations (SVs) were detected in WGD + samples. The breakpoints were enriched in the chr17: 39 Mb-40 Mb tile, which contained the HER2 gene, resulting in the formation of tyfonas, breakage-fusion-bridge cycles, and double minutes. These complex SVs may be involved in the evolutionary mechanisms of the dramatic increase of HER2 copy number.

CONCLUSION

Our work revealed that the WGD + clone might be a critical evolution step for liver metastasis and favoured following complex SV of breast cancer.

Re-analysis of hepatitis B virus integration sites reveals potential new loci associated with oncogenesis in hepatocellular carcinoma

BACKGROUND

Hepatitis B virus (HBV) is a major cause of hepatocellular carcinoma (HCC). HBV DNA can get integrated into the hepatocyte genome to promote carcinogenesis. However, the precise mechanism by which the integrated HBV genome promotes HCC has not been elucidated.

AIM

To analyze the features of HBV integration in HCC using a new reference database and integration detection method.

METHODS

Published data, consisting of 426 Liver tumor samples and 426 paired adjacent non-tumor samples, were re-analyzed to identify the integration sites. Genome Reference Consortium Human Build 38 (GRCh38) and Telomere-to-Telomere Consortium CHM13 (T2T-CHM13 (v2.0)) were used as the human reference genomes. In contrast, human genome 19 (hg19) was used in the original study. In addition, GRIDSS VIRUSBreakend was used to detect HBV integration sites, whereas high-throughput viral integration detection (HIVID) was applied in the original study (HIVID-hg19).

RESULTS

A total of 5361 integration sites were detected using T2T-CHM13. In the tumor samples, integration hotspots in the cancer driver genes, such as TERT and KMT2B, were consistent with those in the original study. GRIDSS VIRUSBreakend detected integrations in more samples than by HIVID-hg19. Enrichment of integration was observed at chromosome 11q13.3, including the CCND1 pro-moter, in tumor samples. Recurrent integration sites were observed in mitochondrial genes.

CONCLUSION

GRIDSS VIRUSBreakend using T2T-CHM13 is accurate and sensitive in detecting HBV integration. Re-analysis provides new insights into the regions of HBV integration and their potential roles in HCC development.

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

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

Pan-cancer whole-genome comparison of primary and metastatic solid tumours

Metastatic cancer remains an almost inevitably lethal disease1-3. A better understanding of disease progression and response to therapies therefore remains of utmost importance. Here we characterize the genomic differences between early-stage untreated primary tumours and late-stage treated metastatic tumours using a harmonized pan-cancer analysis (or reanalysis) of two unpaired primary4 and metastatic5 cohorts of 7,108 whole-genome-sequenced tumours. Metastatic tumours in general have a lower intratumour heterogeneity and a conserved karyotype, displaying only a modest increase in mutations, although frequencies of structural variants are elevated overall. Furthermore, highly variable tumour-specific contributions of mutational footprints of endogenous (for example, SBS1 and APOBEC) and exogenous mutational processes (for example, platinum treatment) are present. The majority of cancer types had either moderate genomic differences (for example, lung adenocarcinoma) or highly consistent genomic portraits (for example, ovarian serous carcinoma) when comparing early-stage and late-stage disease. Breast, prostate, thyroid and kidney renal clear cell carcinomas and pancreatic neuroendocrine tumours are clear exceptions to the rule, displaying an extensive transformation of their genomic landscape in advanced stages. Exposure to treatment further scars the tumour genome and introduces an evolutionary bottleneck that selects for known therapy-resistant drivers in approximately half of treated patients. Our data showcase the potential of pan-cancer whole-genome analysis to identify distinctive features of late-stage tumours and provide a valuable resource to further investigate the biological basis of cancer and resistance to therapies.

ERα-associated translocations underlie oncogene amplifications in breast cancer

Focal copy-number amplification is an oncogenic event. Although recent studies have revealed the complex structure1-3 and the evolutionary trajectories4 of oncogene amplicons, their origin remains poorly understood. Here we show that focal amplifications in breast cancer frequently derive from a mechanism-which we term translocation-bridge amplification-involving inter-chromosomal translocations that lead to dicentric chromosome bridge formation and breakage. In 780 breast cancer genomes, we observe that focal amplifications are frequently connected to each other by inter-chromosomal translocations at their boundaries. Subsequent analysis indicates the following model: the oncogene neighbourhood is translocated in G1 creating a dicentric chromosome, the dicentric chromosome is replicated, and as dicentric sister chromosomes segregate during mitosis, a chromosome bridge is formed and then broken, with fragments often being circularized in extrachromosomal DNAs. This model explains the amplifications of key oncogenes, including ERBB2 and CCND1. Recurrent amplification boundaries and rearrangement hotspots correlate with oestrogen receptor binding in breast cancer cells. Experimentally, oestrogen treatment induces DNA double-strand breaks in the oestrogen receptor target regions that are repaired by translocations, suggesting a role of oestrogen in generating the initial translocations. A pan-cancer analysis reveals tissue-specific biases in mechanisms initiating focal amplifications, with the breakage-fusion-bridge cycle prevalent in some and the translocation-bridge amplification in others, probably owing to the different timing of DNA break repair. Our results identify a common mode of oncogene amplification and propose oestrogen as its mechanistic origin in breast cancer.

Comparative analyses define differences between BHD-associated renal tumour and sporadic chromophobe renal cell carcinoma

BACKGROUND

Birt-Hogg-Dubé (BHD) syndrome, caused by germline alteration of folliculin (FLCN) gene, develops hybrid oncocytic/chromophobe tumour (HOCT) and chromophobe renal cell carcinoma (ChRCC), whereas sporadic ChRCC does not harbor FLCN alteration. To date, molecular characteristics of these similar histological types of tumours have been incompletely elucidated.

METHODS

To elucidate renal tumourigenesis of BHD-associated renal tumours and sporadic renal tumours, we conducted whole genome sequencing (WGS) and RNA-sequencing (RNA-seq) of sixteen BHD-associated renal tumours from nine unrelated BHD patients, twenty-one sporadic ChRCCs and seven sporadic oncocytomas. We then compared somatic mutation profiles with FLCN variants and RNA expression profiles between BHD-associated renal tumours and sporadic renal tumours.

FINDINGS

RNA-seq analysis revealed that BHD-associated renal tumours and sporadic renal tumours have totally different expression profiles. Sporadic ChRCCs were clustered into two distinct clusters characterized by L1CAM and FOXI1 expressions, molecular markers for renal tubule subclasses. Increased mitochondrial DNA (mtDNA) copy number with fewer variants was observed in BHD-associated renal tumours compared to sporadic ChRCCs. Cell-of-origin analysis using WGS data demonstrated that BHD-associated renal tumours and sporadic ChRCCs may arise from different cells of origin and second hit FLCN alterations may occur in early third decade of life in BHD patients.

INTERPRETATION

These data further our understanding of renal tumourigenesis of these two different types of renal tumours with similar histology.

FUNDING

This study was supported by JSPS KAKENHI Grants, RIKEN internal grant, and the Intramural Research Program of the National Institutes of Health (NIH), National Cancer Institute (NCI), Center for Cancer Research.

Targeting homologous recombination deficiency in uterine leiomyosarcoma

BACKGROUND

Uterine leiomyosarcoma (uLMS) is a rare and aggressive gynaecological malignancy, with individuals with advanced uLMS having a five-year survival of < 10%. Mutations in the homologous recombination (HR) DNA repair pathway have been observed in ~ 10% of uLMS cases, with reports of some individuals benefiting from poly (ADP-ribose) polymerase (PARP) inhibitor (PARPi) therapy, which targets this DNA repair defect. In this report, we screened individuals with uLMS, accrued nationally, for mutations in the HR repair pathway and explored new approaches to therapeutic targeting.

METHODS

A cohort of 58 individuals with uLMS were screened for HR Deficiency (HRD) using whole genome sequencing (WGS), whole exome sequencing (WES) or NGS panel testing. Individuals identified to have HRD uLMS were offered PARPi therapy and clinical outcome details collected. Patient-derived xenografts (PDX) were generated for therapeutic targeting.

RESULTS

All 13 uLMS samples analysed by WGS had a dominant COSMIC mutational signature 3; 11 of these had high genome-wide loss of heterozygosity (LOH) (> 0.2) but only two samples had a CHORD score > 50%, one of which had a homozygous pathogenic alteration in an HR gene (deletion in BRCA2). A further three samples harboured homozygous HRD alterations (all deletions in BRCA2), detected by WES or panel sequencing, with 5/58 (9%) individuals having HRD uLMS. All five individuals gained access to PARPi therapy. Two of three individuals with mature clinical follow up achieved a complete response or durable partial response (PR) with the subsequent addition of platinum to PARPi upon minor progression during initial PR on PARPi. Corresponding PDX responses were most rapid, complete and sustained with the PARP1-specific PARPi, AZD5305, compared with either olaparib alone or olaparib plus cisplatin, even in a paired sample of a BRCA2-deleted PDX, derived following PARPi therapy in the patient, which had developed PARPi-resistance mutations in PRKDC, encoding DNA-PKcs.

CONCLUSIONS

Our work demonstrates the value of identifying HRD for therapeutic targeting by PARPi and platinum in individuals with the aggressive rare malignancy, uLMS and suggests that individuals with HRD uLMS should be included in trials of PARP1-specific PARPi.

The peptide woods are lovely, dark and deep: Hunting for novel cancer antigens

Harnessing the patient's immune system to control a tumor is a proven avenue for cancer therapy. T cell therapies as well as therapeutic vaccines, which target specific antigens of interest, are being explored as treatments in conjunction with immune checkpoint blockade. For these therapies, selecting the best suited antigens is crucial. Most of the focus has thus far been on neoantigens that arise from tumor-specific somatic mutations. Although there is clear evidence that T-cell responses against mutated neoantigens are protective, the large majority of these mutations are not immunogenic. In addition, most somatic mutations are unique to each individual patient and their targeting requires the development of individualized approaches. Therefore, novel antigen types are needed to broaden the scope of such treatments. We review high throughput approaches for discovering novel tumor antigens and some of the key challenges associated with their detection, and discuss considerations when selecting tumor antigens to target in the clinic.