Haplotype-resolved karyotype construction from Hi-C data using refLinker.

Chromosomal aberrations are prevalent in cancer genomes, yet it remains challenging to resolve the long-range structure of rearranged chromosomes. A key problem is to determine the chromosomal origin of rearranged genomic segments, which requires chromosome-length haplotype information. Here we describe refLinker, a new computational method for whole-chromosome haplotype inference using external reference panels and Hi-C. We show that refLinker ensures consistent long-range phasing accuracy in both diploid human genomes and aneuploid cancers, including regions with loss-of-heterozygosity and high-level focal amplification. We further demonstrate the feasibility of complex genome reconstruction using haplotype-specific Hi-C contacts, revealing new karyotype features in two widely studied cancer cell lines. Together, these findings provide a new framework for the complete resolution of long-range chromosome structure in complex genomes and highlight the unique advantages of Hi-C data for reconstructing cancer genomes with chromosome-scale continuity.

Genomic signatures of past and present chromosomal instability in Barrett's esophagus and early esophageal adenocarcinoma.

The progression of precancerous lesions to malignancy is often accompanied by increasing complexity of chromosomal alterations but how these alterations arise is poorly understood. Here we perform haplotype-specific analysis of chromosomal copy-number evolution in the progression of Barrett's esophagus (BE) to esophageal adenocarcinoma (EAC) on multiregional whole-genome sequencing data of BE with dysplasia and microscopic EAC foci. We identify distinct patterns of copy-number evolution indicating multigenerational chromosomal instability that is initiated by cell division errors but propagated only after p53 loss. While abnormal mitosis, including whole-genome duplication, underlies chromosomal copy-number changes, segmental alterations display signatures of successive breakage-fusion-bridge cycles and chromothripsis of unstable dicentric chromosomes. Our analysis elucidates how multigenerational chromosomal instability generates copy-number variation in BE cells, precipitates complex alterations including DNA amplifications, and promotes their independent clonal expansion and transformation. In particular, we suggest sloping copy-number variation as a signature of ongoing chromosomal instability that precedes copy-number complexity. These findings suggest copy-number heterogeneity in advanced cancers originates from chromosomal instability in precancerous cells and such instability may be identified from the presence of sloping copy-number variation in bulk sequencing data.

A genetic basis for cancer sex differences revealed in Xp11 translocation renal cell carcinoma.

Xp11 translocation renal cell carcinoma (tRCC) is a female-predominant kidney cancer driven by translocations between the TFE3 gene on chromosome Xp11.2 and partner genes located on either chrX or on autosomes. The rearrangement processes that underlie TFE3 fusions, and whether they are linked to the female sex bias of this cancer, are largely unexplored. Moreover, whether oncogenic TFE3 fusions arise from both the active and inactive X chromosomes in females remains unknown. Here we address these questions by haplotype-specific analyses of whole-genome sequences of 29 tRCC samples from 15 patients and by re-analysis of 145 published tRCC whole-exome sequences. We show that TFE3 fusions universally arise as reciprocal translocations with minimal DNA loss or insertion at paired break ends. Strikingly, we observe a near exact 2:1 female:male ratio in TFE3 fusions arising via X:autosomal translocation (but not via X inversion), which accounts for the female predominance of tRCC. This 2:1 ratio is at least partially attributable to oncogenic fusions involving the inactive X chromosome and is accompanied by partial re-activation of silenced chrX genes on the rearranged chromosome. Our results highlight how somatic alterations involving the X chromosome place unique constraints on tumor initiation and exemplify how genetic rearrangements of the sex chromosomes can underlie cancer sex differences.

Heritable transcriptional defects from aberrations of nuclear architecture.

Transcriptional heterogeneity due to plasticity of the epigenetic state of chromatin contributes to tumour evolution, metastasis and drug resistance1-3. However, the mechanisms that cause this epigenetic variation are incompletely understood. Here we identify micronuclei and chromosome bridges, aberrations in the nucleus common in cancer4,5, as sources of heritable transcriptional suppression. Using a combination of approaches, including long-term live-cell imaging and same-cell single-cell RNA sequencing (Look-Seq2), we identified reductions in gene expression in chromosomes from micronuclei. With heterogeneous penetrance, these changes in gene expression can be heritable even after the chromosome from the micronucleus has been re-incorporated into a normal daughter cell nucleus. Concomitantly, micronuclear chromosomes acquire aberrant epigenetic chromatin marks. These defects may persist as variably reduced chromatin accessibility and reduced gene expression after clonal expansion from single cells. Persistent transcriptional repression is strongly associated with, and may be explained by, markedly long-lived DNA damage. Epigenetic alterations in transcription may therefore be inherently coupled to chromosomal instability and aberrations in nuclear architecture.

Determining Complete Chromosomal Haplotypes by mLinker.

Haplotype ("haploid genotype") phase is the combination of genotypes at sites of genetic variation along a chromosome [1]. We previously demonstrated that the complete chromosomal haplotype of diploid human genomes can be determined using molecular linkage from Hi-C sequencing and linked-reads sequencing [2]. In this chapter, we present a step-by-step guide to perform this analysis using mLinker, a software package for haplotype inference.

Examining the cooperation between extrachromosomal DNA circles.

In a departure from previous findings, new results suggest that free-floating pieces of DNA which carry additional copies of cancer-driving genes do not tend to cluster or have increased transcription.

Somatic XIST activation and features of X chromosome inactivation in male human cancers.

Expression of the non-coding RNA XIST is essential for initiating X chromosome inactivation (XCI) during early development in female mammals. As the main function of XCI is to enable dosage compensation of chromosome X genes between the sexes, XCI and XIST expression are generally absent in male normal tissues, except in germ cells and in individuals with supernumerary X chromosomes. Via a systematic analysis of public sequencing data of both cancerous and normal tissues, we report that XIST is somatically activated in a subset of male human cancers across diverse lineages. Some of these cancers display hallmarks of XCI, including silencing of gene expression, reduced chromatin accessibility, and increased DNA methylation across chromosome X, suggesting that the developmentally restricted, female-specific program of XCI can be somatically accessed in male cancers.

Cell-free DNA for the detection of emerging treatment failure in relapsed/ refractory multiple myeloma.

Interrogation of cell-free DNA (cfDNA) represents an emerging approach to non-invasively estimate disease burden in multiple myeloma (MM). Here, we examined low-pass whole genome sequencing (LPWGS) of cfDNA for its predictive value in relapsed/ refractory MM (RRMM). We observed that cfDNA positivity, defined as ≥10% tumor fraction by LPWGS, was associated with significantly shorter progression-free survival (PFS) in an exploratory test cohort of 16 patients who were actively treated on diverse regimens. We prospectively determined the predictive value of cfDNA in 86 samples from 45 RRMM patients treated with elotuzumab, pomalidomide, bortezomib, and dexamethasone in a phase II clinical trial (NCT02718833). PFS in patients with tumor-positive and -negative cfDNA after two cycles of treatment was 1.6 and 17.6 months, respectively (HR 7.6, P < 0.0001). Multivariate hazard modelling confirmed cfDNA as independent risk factor (HR 96.6, P = 6.92e-05). While correlating with serum-free light chains and bone marrow, cfDNA additionally discriminated patients with poor PFS among those with the same response by IMWG criteria. In summary, detectability of MM-derived cfDNA, as a measure of substantial tumor burden with therapy, independently predicts poor PFS and may provide refinement for standard-of-care response parameters to identify patients with poor response to treatment earlier than is currently feasible.

Whole chromosome loss and genomic instability in mouse embryos after CRISPR-Cas9 genome editing.

Karyotype alterations have emerged as on-target complications from CRISPR-Cas9 genome editing. However, the events that lead to these karyotypic changes in embryos after Cas9-treatment remain unknown. Here, using imaging and single-cell genome sequencing of 8-cell stage embryos, we track both spontaneous and Cas9-induced karyotype aberrations through the first three divisions of embryonic development. We observe the generation of abnormal structures of the nucleus that arise as a consequence of errors in mitosis, including micronuclei and chromosome bridges, and determine their contribution to common karyotype aberrations including whole chromosome loss that has been recently reported after editing in embryos. Together, these data demonstrate that Cas9-mediated germline genome editing can lead to unwanted on-target side effects, including major chromosome structural alterations that can be propagated over several divisions of embryonic development.

Determination of complete chromosomal haplotypes by bulk DNA sequencing.

Haplotype phase represents the collective genetic variation between homologous chromosomes and is an essential feature of non-haploid genomes. Here we describe a computational strategy to reliably determine complete whole-chromosome haplotypes using a combination of bulk long-range sequencing and Hi-C sequencing. We demonstrate that this strategy can resolve the haplotypes of parental chromosomes in diploid human genomes with high precision (>99%) and completeness (>98%) and assemble the syntenic structure of rearranged chromosomes in aneuploid cancer genomes at base pair level resolution. Our work enables direct interrogation of chromosome-specific alterations and chromatin reorganization using bulk DNA sequencing.

Chromothripsis as an on-target consequence of CRISPR-Cas9 genome editing.

Genome editing has therapeutic potential for treating genetic diseases and cancer. However, the currently most practicable approaches rely on the generation of DNA double-strand breaks (DSBs), which can give rise to a poorly characterized spectrum of chromosome structural abnormalities. Here, using model cells and single-cell whole-genome sequencing, as well as by editing at a clinically relevant locus in clinically relevant cells, we show that CRISPR-Cas9 editing generates structural defects of the nucleus, micronuclei and chromosome bridges, which initiate a mutational process called chromothripsis. Chromothripsis is extensive chromosome rearrangement restricted to one or a few chromosomes that can cause human congenital disease and cancer. These results demonstrate that chromothripsis is a previously unappreciated on-target consequence of CRISPR-Cas9-generated DSBs. As genome editing is implemented in the clinic, the potential for extensive chromosomal rearrangements should be considered and monitored.

No pains, no gains: how chromosome fragmentation promotes gene amplification.

Combining live-cell imaging, cytogenetics, genome sequencing, and in vitro evolution, Shoshani et al. (2020) revealed deep connections between chromothripsis, the catastrophic shattering of a chromosome in abnormal nuclear structures, and gene amplification, a frequent culprit of oncogenic activation.

Mechanisms generating cancer genome complexity from a single cell division error.

The chromosome breakage-fusion-bridge (BFB) cycle is a mutational process that produces gene amplification and genome instability. Signatures of BFB cycles can be observed in cancer genomes alongside chromothripsis, another catastrophic mutational phenomenon. We explain this association by elucidating a mutational cascade that is triggered by a single cell division error-chromosome bridge formation-that rapidly increases genomic complexity. We show that actomyosin forces are required for initial bridge breakage. Chromothripsis accumulates, beginning with aberrant interphase replication of bridge DNA. A subsequent burst of DNA replication in the next mitosis generates extensive DNA damage. During this second cell division, broken bridge chromosomes frequently missegregate and form micronuclei, promoting additional chromothripsis. We propose that iterations of this mutational cascade generate the continuing evolution and subclonal heterogeneity characteristic of many human cancers.

Comprehensive analysis of chromothripsis in 2,658 human cancers using whole-genome sequencing.

Chromothripsis is a mutational phenomenon characterized by massive, clustered genomic rearrangements that occurs in cancer and other diseases. Recent studies in selected cancer types have suggested that chromothripsis may be more common than initially inferred from low-resolution copy-number data. Here, as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA), we analyze patterns of chromothripsis across 2,658 tumors from 38 cancer types using whole-genome sequencing data. We find that chromothripsis events are pervasive across cancers, with a frequency of more than 50% in several cancer types. Whereas canonical chromothripsis profiles display oscillations between two copy-number states, a considerable fraction of events involve multiple chromosomes and additional structural alterations. In addition to non-homologous end joining, we detect signatures of replication-associated processes and templated insertions. Chromothripsis contributes to oncogene amplification and to inactivation of genes such as mismatch-repair-related genes. These findings show that chromothripsis is a major process that drives genome evolution in human cancer.

Cardioprotective mechanisms of salvianic acid A sodium in rats with myocardial infarction based on proteome and transcriptome analysis.

Ischemic heart diseases (IHDs) cause great morbidity and mortality worldwide, necessitating effective treatment. Salvianic acid A sodium (SAAS) is an active compound derived from the well-known herbal medicine Danshen, which has been widely used for clinical treatment of cardiovascular diseases in China. This study aimed to confirm the cardioprotective effects of SAAS in rats with myocardial infarction and to investigate the underlying molecular mechanisms based on proteome and transcriptome profiling of myocardial tissue. The results showed that SAAS effectively protected against myocardial injury and improved cardiac function. The differentially expressed proteins and genes included important structural molecules, receptors, transcription factors, and cofactors. Functional enrichment analysis indicated that SAAS participated in the regulation of actin cytoskeleton, phagosome, focal adhesion, tight junction, apoptosis, MAPK signaling, and Wnt signaling pathways, which are closely related to cardiovascular diseases. SAAS may exert its cardioprotective effect by targeting multiple pathways at both the proteome and transcriptome levels. This study has provided not only new insights into the pathogenesis of myocardial infarction but also a road map of the cardioprotective molecular mechanisms of SAAS, which may provide pharmacological evidence to aid in its clinical application.

Using whole genome scores to compare three clinical phenotyping methods in complex diseases.

Genome-wide association studies depend on accurate ascertainment of patient phenotype. However, phenotyping is difficult, and it is often treated as an afterthought in these studies because of the expense involved. Electronic health records (EHRs) may provide higher fidelity phenotypes for genomic research than other sources such as administrative data. We used whole genome association models to evaluate different EHR and administrative data-based phenotyping methods in a cohort of 16,858 Caucasian subjects for type 1 diabetes mellitus, type 2 diabetes mellitus, coronary artery disease and breast cancer. For each disease, we trained and evaluated polygenic models using three different phenotype definitions: phenotypes derived from billing data, the clinical problem list, or a curated phenotyping algorithm. We observed that for these diseases, the curated phenotype outperformed the problem list, and the problem list outperformed administrative billing data. This suggests that using advanced EHR-derived phenotypes can further increase the power of genome-wide association studies.