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.

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.

MCM8IP activates the MCM8-9 helicase to promote DNA synthesis and homologous recombination upon DNA damage.

Homologous recombination (HR) mediates the error-free repair of DNA double-strand breaks to maintain genomic stability. Here we characterize C17orf53/MCM8IP, an OB-fold containing protein that binds ssDNA, as a DNA repair factor involved in HR. MCM8IP-deficient cells exhibit HR defects, especially in long-tract gene conversion, occurring downstream of RAD51 loading, consistent with a role for MCM8IP in HR-dependent DNA synthesis. Moreover, loss of MCM8IP confers cellular sensitivity to crosslinking agents and PARP inhibition. Importantly, we report that MCM8IP directly associates with MCM8-9, a helicase complex mutated in primary ovarian insufficiency, and RPA1. We additionally show that the interactions of MCM8IP with MCM8-9 and RPA facilitate HR and promote replication fork progression and cellular viability in response to treatment with crosslinking agents. Mechanistically, MCM8IP stimulates the helicase activity of MCM8-9. Collectively, our work identifies MCM8IP as a key regulator of MCM8-9-dependent DNA synthesis during DNA recombination and replication.

The human Shu complex functions with PDS5B and SPIDR to promote homologous recombination.

RAD51 plays a central role in homologous recombination during double-strand break repair and in replication fork dynamics. Misregulation of RAD51 is associated with genetic instability and cancer. RAD51 is regulated by many accessory proteins including the highly conserved Shu complex. Here, we report the function of the human Shu complex during replication to regulate RAD51 recruitment to DNA repair foci and, secondly, during replication fork restart following replication fork stalling. Deletion of the Shu complex members, SWS1 and SWSAP1, using CRISPR/Cas9, renders cells specifically sensitive to the replication fork stalling and collapse caused by methyl methanesulfonate and mitomycin C exposure, a delayed and reduced RAD51 response, and fewer sister chromatid exchanges. Our additional analysis identified SPIDR and PDS5B as novel Shu complex interacting partners and genetically function in the same pathway upon DNA damage. Collectively, our study uncovers a protein complex, which consists of SWS1, SWSAP1, SPIDR and PDS5B, involved in DNA repair and provides insight into Shu complex function and composition.

Differential Requirements for the RAD51 Paralogs in Genome Repair and Maintenance in Human Cells.

Deficiency in several of the classical human RAD51 paralogs [RAD51B, RAD51C, RAD51D, XRCC2 and XRCC3] is associated with cancer predisposition and Fanconi anemia. To investigate their functions, isogenic disruption mutants for each were generated in non-transformed MCF10A mammary epithelial cells and in transformed U2OS and HEK293 cells. In U2OS and HEK293 cells, viable ablated clones were readily isolated for each RAD51 paralog; in contrast, with the exception of RAD51B, RAD51 paralogs are cell-essential in MCF10A cells. Underlining their importance for genomic stability, mutant cell lines display variable growth defects, impaired sister chromatid recombination, reduced levels of stable RAD51 nuclear foci, and hyper-sensitivity to mitomycin C and olaparib, with the weakest phenotypes observed in RAD51B-deficient cells. Altogether these observations underscore the contributions of RAD51 paralogs in diverse DNA repair processes, and demonstrate essential differences in different cell types. Finally, this study will provide useful reagents to analyze patient-derived mutations and to investigate mechanisms of chemotherapeutic resistance deployed by cancers.

Secondary Somatic Mutations Restoring RAD51C and RAD51D Associated with Acquired Resistance to the PARP Inhibitor Rucaparib in High-Grade Ovarian Carcinoma.

High-grade epithelial ovarian carcinomas containing mutated BRCA1 or BRCA2 (BRCA1/2) homologous recombination (HR) genes are sensitive to platinum-based chemotherapy and PARP inhibitors (PARPi), while restoration of HR function due to secondary mutations in BRCA1/2 has been recognized as an important resistance mechanism. We sequenced core HR pathway genes in 12 pairs of pretreatment and postprogression tumor biopsy samples collected from patients in ARIEL2 Part 1, a phase II study of the PARPi rucaparib as treatment for platinum-sensitive, relapsed ovarian carcinoma. In 6 of 12 pretreatment biopsies, a truncation mutation in BRCA1, RAD51C, or RAD51D was identified. In five of six paired postprogression biopsies, one or more secondary mutations restored the open reading frame. Four distinct secondary mutations and spatial heterogeneity were observed for RAD51CIn vitro complementation assays and a patient-derived xenograft, as well as predictive molecular modeling, confirmed that resistance to rucaparib was associated with secondary mutations.Significance: Analyses of primary and secondary mutations in RAD51C and RAD51D provide evidence for these primary mutations in conferring PARPi sensitivity and secondary mutations as a mechanism of acquired PARPi resistance. PARPi resistance due to secondary mutations underpins the need for early delivery of PARPi therapy and for combination strategies. Cancer Discov; 7(9); 984-98. ©2017 AACR.See related commentary by Domchek, p. 937See related article by Quigley et al., p. 999See related article by Goodall et al., p. 1006This article is highlighted in the In This Issue feature, p. 920.