Mate-pair Sequencing Research Articles

Use of Mate-Pair Sequencing (MPseq) to Elucidate a Complex BCR-ABL1 Rearrangement Observed in a Newly Diagnosed Case of Chronic Myeloid Leukemia

Use of Mate-Pair Sequencing (MPseq) to Elucidate a Complex BCR-ABL1 Rearrangement Observed in a Newly Diagnosed Case of Chronic Myeloid Leukemia

Effects of hydroxyurea on CNV induction in the mouse germline.

Copy number variants (CNVs) are important in genome variation and genetic disease, with new mutations arising frequently in the germline and somatic cells. Replication stress caused by aphidicolin and hydroxyurea (HU) is a potent inducer of de novo CNVs in cultured mammalian cells. HU is used extensively for long-term management of sickle cell disease. Here, we examined the effects of HU treatment on germline CNVs in vivo in male mice to explore whether replication stress can act as a CNV mutagen in germline mitotic divisions as in cultured cells and whether this would support a concern for increased CNV mutations in offspring of men treated with HU. Several trials of HU administration were performed by oral gavage and subcutaneous pump, with CNVs characterized in C57BL/6 x C3H/HeJ hybrid mouse offspring by microarray and mate-pair sequencing. HU had a short half-life of ~14 min and a narrow dose window over which studies could be performed while maintaining fertility. Tissue histopathology and reticulocyte micronucleus assays verified that doses had a substantial tissue and genetic toxicity. CNVs were readily detected in offspring that originated in both paternal and maternal mouse strains, as de novo and inherited events. However, HU did not increase CNV formation above baseline levels. These results reveal a high rate of CNV mutagenesis in the mouse germline but do not support the hypothesis that HU would increase CNV formation during mammalian spermatogenesis, perhaps due to highly toxic effects on sperm development or experimental variables related to HU pharmacology in mice. Environ. Mol. Mutagen. 59:698-714, 2018. © 2018 Wiley Periodicals, Inc.

The first complete genomic structure of Butyrivibrio fibrisolvens and its chromid

Butyrivibrio fibrisolvens forms part of the gastrointestinal microbiome of ruminants and other mammals, including humans. Indeed, it is one of the most common bacteria found in the rumen and plays an important role in ruminal fermentation of polysaccharides, yet, to date, there is no closed reference genome published for this species in any ruminant animal. We successfully assembled the nearly complete genome sequence of B. fibrisolvens strain INBov1 isolated from cow rumen using Illumina paired-end reads, 454 Roche single-end and mate pair sequencing technology. Additionally, we constructed an optical restriction map of this strain to aid in scaffold ordering and positioning, and completed the first genomic structure of this species. Moreover, we identified and assembled the first chromid of this species (pINBov266). The INBov1 genome encodes a large set of genes involved in the cellulolytic process but lacks key genes. This seems to indicate that B. fibrisolvens plays an important role in ruminal cellulolytic processes, but does not have autonomous cellulolytic capacity. When searching for genes involved in the biohydrogenation of unsaturated fatty acids, no linoleate isomerase gene was found in this strain. INBov1 does encode oleate hydratase genes known to participate in the hydrogenation of oleic acids. Furthermore, INBov1 contains an enolase gene, which has been recently determined to participate in the synthesis of conjugated linoleic acids. This work confirms the presence of a novel chromid in B. fibrisolvens and provides a new potential reference genome sequence for this species, providing new insight into its role in biohydrogenation and carbohydrate degradation.

Genome sequencing and conservation genomics in the Scandinavian wolverine population.

Genetic approaches have proved valuable to the study and conservation of endangered populations, especially for monitoring programs, and there is potential for further developments in this direction by extending analyses to the genomic level. We assembled the genome of the wolverine (Gulo gulo), a mustelid that in Scandinavia has recently recovered from a significant population decline, and obtained a 2.42Gb draft sequence representing >85% of the genome and including >21,000protein-coding genes. We then performed whole-genome resequencing of 10Scandinavian wolverines for population genomic and demographic analyses. Genetic diversity was among the lowest detected in a red-listed population (mean genome-wide nucleotide diversity of 0.05%). Results of the demographic analyses indicated a long-term decline of the effective population size (Ne ) from 10,000well before the last glaciation to <500 after this period. Current Ne appeared even lower. The genome-wide FIS level was 0.089 (possibly signaling inbreeding), but this effect was not observed when analyzing a set of highly variable SNP markers, illustrating that such markers can give a biased picture of the overall character of genetic diversity. We found significant population structure, which has implications for population connectivity and conservation. We used an integrated microfluidic circuit chip technology to develop an SNP-array consisting of 96 highly informative markers that, together with a multiplex pre-amplification step, was successfully applied to low-quality DNA from scat samples. Our findings will inform management, conservation, and genetic monitoring of wolverines and serve as a genomic roadmap that can be applied to other endangered species. The approach used here can be generally utilized in other systems, but we acknowledge the trade-off between investing in genomic resources and direct conservation actions.

Different whole genome sequencing strategies for the clinical management of patients with oropharyngeal squamous cell carcinoma

Oropharyngeal squamous cell carcinoma (OPSCC) is a cancer of the base of the tongue and tonsils. This cancer was previously found in individuals with a smoking and drinking history but is increasingly being promoted by the presence of high risk types of human papillomavirus (HPV). In an attempt to determine the exact role that HPV plays in the development of this cancer we have been using mate-pair next generation sequencing (MP-Seq). Using MP-Seq we have found that HPV is only integrated into the human genome in 30% of HPV-positive OPSCCs, which is dramatically different than what is found in cervical cancer. MP-Seq not only can determine the physical status of HPV in a given sample, but can also ascertain genome-wide changes in any cancer. We describe how MP-Seq could be used in the clinical management of HPV-positive OPSCC patients. We have also explored another whole genome sequencing strategy that has been developed by the Beijing Genomics Institute. They can currently do whole genome sequencing for a total cost (and this includes analysis) of $600 and are developing methodologies to reduce this to $100. I will describe how these types of whole genome sequencing strategies could be used to direct therapies for these patients including facilitating the use of the liquid biopsy to monitor patients during their clinical treatment.

Optical and physical mapping with local finishing enables megabase-scale resolution of agronomically important regions in the wheat genome

BackgroundNumerous scaffold-level sequences for wheat are now being released and, in this context, we report on a strategy for improving the overall assembly to a level comparable to that of the human genome.ResultsUsing chromosome 7A of wheat as a model, sequence-finished megabase-scale sections of this chromosome were established by combining a new independent assembly using a bacterial artificial chromosome (BAC)-based physical map, BAC pool paired-end sequencing, chromosome-arm-specific mate-pair sequencing and Bionano optical mapping with the International Wheat Genome Sequencing Consortium RefSeq v1.0 sequence and its underlying raw data. The combined assembly results in 18 super-scaffolds across the chromosome. The value of finished genome regions is demonstrated for two approximately 2.5 Mb regions associated with yield and the grain quality phenotype of fructan carbohydrate grain levels. In addition, the 50 Mb centromere region analysis incorporates cytological data highlighting the importance of non-sequence data in the assembly of this complex genome region.ConclusionsSufficient genome sequence information is shown to now be available for the wheat community to produce sequence-finished releases of each chromosome of the reference genome. The high-level completion identified that an array of seven fructosyl transferase genes underpins grain quality and that yield attributes are affected by five F-box-only-protein-ubiquitin ligase domain and four root-specific lipid transfer domain genes. The completed sequence also includes the centromere.

Copy number variant analysis using genome-wide mate-pair sequencing.

Copy number variation (CNV) is a common form of structural variation detected in human genomes, occurring as both constitutional and somatic events. Cytogenetic techniques like chromosomal microarray (CMA) are widely used in analyzing CNVs. However, CMA techniques cannot resolve the full nature of these structural variations (i.e. the orientation and location of associated breakpoint junctions) and must be combined with other cytogenetic techniques, such as karyotyping or FISH, to do so. This makes the development of a next-generation sequencing (NGS) approach capable of resolving both CNVs and breakpoint junctions desirable. Mate-pair sequencing (MPseq) is a NGS technology designed to find large structural rearrangements across the entire genome. Here we present an algorithm capable of performing copy number analysis from mate-pair sequencing data. The algorithm uses a step-wise procedure involving normalization, segmentation, and classification of the sequencing data. The segmentation technique combines both read depth and discordant mate-pair reads to increase the sensitivity and resolution of CNV calls. The method is particularly suited to MPseq, which is designed to detect breakpoint junctions at high resolution. This allows for the classification step to accurately calculate copy number levels at the relatively low read depth of MPseq. Here we compare results for a series of hematological cancer samples that were tested with CMA and MPseq. We demonstrate comparable sensitivity to the state-of-the-art CMA technology, with the benefit of improved breakpoint resolution. The algorithm provides a powerful analytical tool for the analysis of MPseq results in cancer.

1. Clinical utility of mate pair sequencing to detect diagnostic and prognostic chromosomal rearrangements and copy number changes in patients with acute myeloid leukemia

1. Clinical utility of mate pair sequencing to detect diagnostic and prognostic chromosomal rearrangements and copy number changes in patients with acute myeloid leukemia

Abstract 5726: Rearrangement-related peptides with neoantigenic potential in malignant pleural mesothelioma

Abstract Malignant pleural mesothelioma is a disease primarily associated with exposure to the carcinogen asbestos. Consistent with this carcinogenic exposure, cytogenetic analyses have identified multiple recurrent structural chromosomal abnormalities in this malignancy, but more recent high-throughput sequencing evaluations of point mutations suggest that there is a low mutational burden in mesothelioma. Since tumor mutational burden has been correlated with responses to treatment with immune checkpoint inhibitors such as nivolumab, it was not consistent that patients with mesothelioma and low mutation burdens would have similar response rates in clinical trials with immune checkpoint inhibitors as patients with non-small cell lung cancer which is associated with a high mutation burden. In order to reconcile these differences, and given the potential for an improved understanding of the molecular pathogenesis of mesothelioma to improve therapeutic options, we used mate-pair sequencing (MPseq) and RNA sequencing (RNAseq) to understand how structural variants affect the transcriptome. MPseq differs from standard next generation sequencing approaches by tiling the whole genome with larger fragments (2-5kb) to reliably detect structural variants such as insertions, deletions and rearrangements. Amongst 22 mesothelioma specimens there were 1535 chromosomal rearrangements (median 41, range 3-298 per specimen), that resulted in junctions or novel fusions of non-coding DNA or genes. Six-hundred thirty-seven of these rearrangements (median 22, range 5-103 range per specimen) resulted in novel fusions of genes. Many of these inter- or intra-chromosomal rearrangements were consistent with a pattern of chromoanagesis such as chromoplexy or chromothripsis. Chromosomal rearrangements detected by MPseq were used to guide analysis of RNAseq data and revealed that these chromosomal junctions resulted in the expression of 179 novel amino acid sequences (median 5, 0-51 range per specimen). To determine whether transcription of chromosomal rearrangement-related junctions have neoantigenic potential, we used in silico tools to determine whether any of the expressed junctions contained peptides that could be presented by patient-specific HLA molecules. The top candidate rearrangement-related peptides with neoantigenic potential bound patient-specific HLA molecules nearly as well or as well as a positive control in competitive binding assays. Our findings represent the discovery of potential neoantigen expression driven by structural chromosomal rearrangements. These results may have implications for the development of novel therapeutic strategies, the selection of patients to receive immunotherapy, and blood-based treatment monitoring strategies. Citation Format: Aaron S. Mansfield, Tobias Peikert, James B. Smadbeck, Julia B. Udell, Farhad Kosari, Stephen J. Murphy, Hongzheng Ren, Vishnu V. Serla, Janet L. Schaefer Klein, Giannoula Karagouga, Faye R. Harris, Carlos Sosa, Sarah H. Johnson, Wendy Nevala, Svetomir N. Markovic, Aaron O. Bungum, Eric S. Edell, Haidong Dong, John C. Cheville, Marie Christine Aubry, Jin Jen, George Vasmatzis. Rearrangement-related peptides with neoantigenic potential in malignant pleural mesothelioma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5726.

Cryptic breakpoint identified by whole-genome mate-pair sequencing in a rare paternally inherited complex chromosomal rearrangement

BackgroundPrecise characterization of apparently balanced complex chromosomal rearrangements in non-affected individuals is crucial as they may result in reproductive failure, recurrent miscarriages or affected offspring.Case presentationWe present a family, where the non-affected father and daughter were found, using FISH and karyotyping, to be carriers of a three-way complex chromosomal rearrangement [t(6;7;10)(q16.2;q34;q26.1), de novo in the father]. The family suffered from two stillbirths, one miscarriage, and has a son with severe intellectual disability. In the present study, the family was revisited using whole-genome mate-pair sequencing. Interestingly, whole-genome mate-pair sequencing revealed a cryptic breakpoint on derivative (der) chromosome 6 rendering the rearrangement even more complex. FISH using a chromosome (chr) 6 custom-designed probe and a chr10 control probe confirmed that the interstitial chr6 segment, created by the two chr6 breakpoints, was translocated onto der(10). Breakpoints were successfully validated with Sanger sequencing, and small imbalances as well as microhomology were identified. Finally, the complex chromosomal rearrangement breakpoints disrupted the SIM1, GRIK2, CNTNAP2, and PTPRE genes without causing any phenotype development.ConclusionsIn contrast to the majority of maternally transmitted complex chromosomal rearrangement cases, our study investigated a rare case where a complex chromosomal rearrangement, which most probably resulted from a Type IV hexavalent during the pachytene stage of meiosis I, was stably transmitted from a fertile father to his non-affected daughter. Whole-genome mate-pair sequencing proved highly successful in identifying cryptic complexity, which consequently provided further insight into the meiotic segregation of chromosomes and the increased reproductive risk in individuals carrying the specific complex chromosomal rearrangement. We propose that such complex rearrangements should be characterized in detail using a combination of conventional cytogenetic and NGS-based approaches to aid in better prenatal preimplantation genetic diagnosis and counseling in couples with reproductive problems.

Proposal for Modification of Cahan's Criteria Utilizing Molecular Genetic Analyses for Cases without Baseline Histopathology: A Unique Method Applicable to Primary Radiosurgery.

Background Cahan's criteria have been utilized since 1948 to establish causality between prior radiation treatment and the development of secondary malignancy. One major criterion specifies that histological and radiographic evidence collected before and after radiation treatment must confirm separate tumor types; however, pretreatment biopsy is rarely obtained prior to radiosurgery for vestibular schwannoma and many other skull base and cranial lesions. Therefore, in these cases Cahan's criteria cannot be validly applied. Objective This article proposes an update to Cahan's criteria using modern molecular genetic analysis for cases lacking baseline histopathology. Methods Mate-pair sequencing and whole exome sequencing of a cerebellopontine angle undifferentiated high-grade pleomorphic sarcoma (UHGPS) that developed after stereotactic radiosurgery of a presumed benign vestibular schwannoma. Results Mate-pair sequencing and whole exome sequencing of the sarcoma revealed complex chromosomal aberrations. Notably, the tumor contained a deletion in the NF2 gene at 22q12 and an in-frame deletion on exon 5 of the remaining copy of NF2 . Biallelic events impacting NF2 are atypical for UHGPS but are characteristic for vestibular schwannoma. These findings help support the conclusion that the UHGPS arose from a benign vestibular schwannoma all along. Conclusions Next-generation sequencing can be successfully applied to a radiation-induced sarcoma when both the original and malignant tumors harbor separate signature genetic markers. As our understanding of the genetic profile of various tumors expand, we believe that next-generation sequencing and other genomic tools will play an increasingly important role in establishing causality between radiation and the development of secondary malignancy.

Risks and Recommendations in Prenatally Detected De Novo Balanced Chromosomal Rearrangements from Assessment of Long-Term Outcomes

The 6%-9% risk of an untoward outcome previously established by Warburton for prenatally detected de novo balanced chromosomal rearrangements (BCRs) does not account for long-term morbidity. We performed long-term follow-up (mean 17 years) of a registry-based nationwide cohort of 41 individuals carrying a prenatally detected de novo BCR with normal first trimester screening/ultrasound scan. We observed a significantly higher frequency of neurodevelopmental and/or neuropsychiatric disorders than in a matched control group (19.5% versus 8.3%, p = 0.04), which was increased to 26.8% upon clinical follow-up. Chromosomal microarray of 32 carriers revealed no pathogenic imbalances, illustrating a low prognostic value when fetal ultrasound scan is normal. In contrast, mate-pair sequencing revealed disrupted genes (ARID1B, NPAS3, CELF4), regulatory domains of known developmental genes (ZEB2, HOXC), and complex BCRs associated with adverse outcomes. Seven unmappable autosomal-autosomal BCRs with breakpoints involving pericentromeric/heterochromatic regions may represent a low-risk group. We performed independent phenotype-aware and blinded interpretation, which accurately predicted benign outcomes (specificity = 100%) but demonstrated relatively low sensitivity for prediction of the clinical outcome in affected carriers (sensitivity = 45%-55%). This sensitivity emphasizes the challenges associated with prenatal risk prediction for long-term morbidity in the absence of phenotypic data given the still immature annotation of the morbidity genome and poorly understood long-range regulatory mechanisms. In conclusion, we upwardly revise the previous estimates of Warburton to a morbidity risk of 27% and recommend sequencing of the chromosomal breakpoints as the first-tier diagnostic test in pregnancies with a de novo BCR.

12. Mate pair sequencing: Unveiling underappreciated complexity and providing clarity to the previously unanswered questions of cytogenetics

12. Mate pair sequencing: Unveiling underappreciated complexity and providing clarity to the previously unanswered questions of cytogenetics

13. Clinical utility of mate pair sequencing to detect diagnostic and prognostic chromosomal rearrangements and copy number changes in patients with acute myeloid leukemia

Acute myeloid leukemia (AML) is the most common acute leukemia in adults and can be grouped based on recurrent genomic abnormalities with diagnostic and prognostic significance. Although cytogenetic characterization classically involves conventional chromosome and/or fluorescence in situ hybridization (FISH) assays, limitations to these labor-intensive techniques include poor resolution of chromosomes banding leading to the inability to precisely/completely identify breakpoints/gene content and difficulty visualizing cryptic rearrangements.

16. The evolutionary conserved developmental regulome define high risk regions for long range position effects

Balanced chromosomal rearrangements (BCRs) have been used for decades to identify disease loci, an approach which has been facilitated by global paired-end sequencing methods to detect chromosomal breakpoints at near sequence level. BCRs do not only cause disease by gene truncation, but also by removing cis-acting regulatory elements, e.g. tissue-specific enhancers, from specific target genes (long-range position effects – LRPE), and/or by bringing external regulatory elements into contact with a new target gene (enhancer adoption). In this study, we show that the majority of the known LRPE-associated breakpoints occur in Topologically Associating Domains (TADs) overlapping with clusters of evolutionary conserved non-genic elements (CNE), and where the dysregulated target genes frequently are marked by large transposon free promoters. By a genome-wide search we defined 316 TADs (cneTADs) with these conservation features, covering 18% of the genome, with 382 predicted target genes that are highly enriched in developmental genes including homeobox genes and other transcription factors. By systematic whole genome mate-pair sequencing, where we more than double the number of mapped two-way BCRs, and include the first large cohort of healthy BCR carriers, we show that >40% of all affected BCR-carriers, including all published ones, have one or two breakpoints in cneTADs, a significantly increase compared to the frequency in healthy BCR carriers. Our data support not only that cneTADs are candidate high-risk LRPE regions, but also that TAD-boundaries are LRPE-boundaries. Moreover, cneTADs are high-risk regions for enhancer adoption events, and a new extreme variant of this, enhancer swapping, where exchanges between two cneTADs may lead to unpredictable, novel phenotypes. Our concerted action will have the potential to establish a saturated map of genotype-phenotype links for this major part of the human developmental regulome, in addition to provide genotype-phenotype links for thousands of truncated protein-coding and non-coding genes.

9. Detection of fusion genes from complex rearrangements reported by genome-wide mate-pair sequencing (MPseq)

9. Detection of fusion genes from complex rearrangements reported by genome-wide mate-pair sequencing (MPseq)

Phenotypic consequences of gene disruption by a balanced de novo translocation involving SLC6A1 and NAA15

Mapping of de novo balanced chromosomal translocations (BCTs) in patients with sporadic poorly characterized disease(s) is an unbiased method of finding candidate gene(s) responsible for the observed symptoms. We present a paediatric patient suffering from epilepsy, developmental delay (DD) and atrial septal defect IIº (ASD) requiring surgery. Karyotyping indicated an apparently balanced de novo reciprocal translocation 46,XX,t(3;4)(p25.3;q31.1), whereas aCGH did not reveal any copy number changes. Using shallow mate-pair whole genome sequencing and direct Sanger sequencing of breakpoint regions we found that translocation disrupted SLC6A1 and NAA15 genes. Our results confirm two previous reports indicating that loss of function of a single allele of SLC6A1 causes epilepsy. In addition, we extend existing evidence that disruption of NAA15 is associated with DD and with congenital heart defects.

P_RNA_scaffolder: a fast and accurate genome scaffolder using paired-end RNA-sequencing reads

BackgroundObtaining complete gene structures is one major goal of genome assembly. Some gene regions are fragmented in low quality and high-quality assemblies. Therefore, new approaches are needed to recover gene regions. Genomes are widely transcribed, generating messenger and non-coding RNAs. These widespread transcripts can be used to scaffold genomes and complete transcribed regions.ResultsWe present P_RNA_scaffolder, a fast and accurate tool using paired-end RNA-sequencing reads to scaffold genomes. This tool aims to improve the completeness of both protein-coding and non-coding genes. After this tool was applied to scaffolding human contigs, the structures of both protein-coding genes and circular RNAs were almost completely recovered and equivalent to those in a complete genome, especially for long proteins and long circular RNAs. Tested in various species, P_RNA_scaffolder exhibited higher speed and efficiency than the existing state-of-the-art scaffolders. This tool also improved the contiguity of genome assemblies generated by current mate-pair scaffolding and third-generation single-molecule sequencing assembly.ConclusionsThe P_RNA_scaffolder can improve the contiguity of genome assembly and benefit gene prediction. This tool is available at http://www.fishbrowser.org/software/P_RNA_scaffolder.

Abstract LB-A19: Targeting ERBB2 pathway in ovarian cancer

Abstract Ovarian cancer (OC) is most lethal gynecologic malignancy. About 75% of OC patients relapse and/or develop chemo-resistant disease after initial response to standard-of-care treatment with platinum-based therapies. HER2 amplifications and overexpression in OC are reported to vary and responses to HER2 inhibitors have been feeble. Next generation sequencing technologies in conjunction with testing using patient-derived xenografts (PDX) allow validation of personalized treatments. Using whole-genome mate-pair next generation sequencing protocol we profiled a set of serous OC to identify copy number changes and gene fusions. Specific focus was on DNA alterations which were therapeutically targetable. The efficiency of anti-Her2 therapy was tested in three different PDX lines with alterations at genes of ERBB pathway and various levels of HER2 protein expression. Tumor in one line harbored a duplication of ERBB2 and high expression of HER2 protein. Second PDX model had a fusion at NRG1 gene, a ligand for HER4 and HER4 receptors, that was highly expressed on protein level. Third line harbored an amplification of NRG3 gene, a known ligand for HER4 receptor. Treatment responses to pertzumab/trastuzumab were compared in each PDX line to standard carboplatin and paclitaxel combination treatment. In all three PDX models HER2-targeted therapy resulted in significant inhibition of tumor growth compared to untreated control. However, the responses in each were inferior to chemotherapy, even in chemo-resistant case. There was an extra benefit when chemotherapy and HER-2 targeted therapy were administered together, as a significant regression of tumor was observed after 6 weeks of treatment using this regiment compared to chemotherapy alone. Post-treatment analysis of the tissues, nevertheless, revealed inhibition of ERBB2 pathway on the level of phosphorylation and expression of downstream targets. In conclusion, targeting of presumably activated ERBB2 pathway alone in OC results in a modest treatment benefit. A combination therapy including both chemo drugs and Her2 inhibitors provides far better response. Further studies are needed to address development of recurrence and sensitivity of recurrent disease to Her2 targeted therapy. Citation Format: Irina V. Kovtun, Zhang Piyan, Faye R. Harris, Xiaonan Hou, Saravut J. Weroha, George Vasmatzis. Targeting ERBB2 pathway in ovarian cancer [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2017 Oct 26-30; Philadelphia, PA. Philadelphia (PA): AACR; Mol Cancer Ther 2018;17(1 Suppl):Abstract nr LB-A19.

Analysis of Chromothripsis by Combined FISH and Microarray Analysis.

Fluorescence in situ hybridization (FISH) to metaphase chromosomes, in conjunction with SNP array, array CGH, or whole genome sequencing, can help determine the organization of abnormal genomes after chromothripsis and other types of complex genome rearrangement. DNA microarrays can identify the changes in copy number, but they do not give information on the organization of the abnormal chromosomes, balanced rearrangements, or abnormalities of the centromeres and other regions comprised of highly repetitive DNA. Many of these details can be determined by the strategic use of metaphase FISH. FISH is a single-cell technique, so it can identify low-frequency chromosome abnormalities, and it can determine which chromosome abnormalities occur in the same or different clonal populations. These are important considerations in cancer. Metaphase chromosomes are intact, so information about abnormalities of the chromosome homologues is preserved. Here we describe strategies for working out the organization of highly rearranged genomes by combining SNP array data with various metaphase FISH methods. This approach can also be used to address some of the uncertainties arising from whole genome or mate-pair sequencing data.