Ultra-high resolution array painting facilitates breakpoint sequencing

Accurate Breakpoint Mapping in Apparently Balanced Translocation Families with Discordant Phenotypes Using Whole Genome Mate-Pair Sequencing

Familial apparently balanced translocations (ABTs) segregating with discordant phenotypes are extremely challenging for interpretation and counseling due to the scarcity of publications and lack of routine techniques for quick investigation. Recently, next generation sequencing has emerged as an efficacious methodology for precise detection of translocation breakpoints. However, studies so far have mainly focused on de novo translocations. The present study focuses specifically on familial cases in order to shed some light to this diagnostic dilemma. Whole-genome mate-pair sequencing (WG-MPS) was applied to map the breakpoints in nine two-way ABT carriers from four families. Translocation breakpoints and patient-specific structural variants were validated by Sanger sequencing and quantitative Real Time PCR, respectively. Identical sequencing patterns and breakpoints were identified in affected and non-affected members carrying the same translocations. PTCD1, ATP5J2-PTCD1, CADPS2, and STPG1 were disrupted by the translocations in three families, rendering them initially as possible disease candidate genes. However, subsequent mutation screening and structural variant analysis did not reveal any pathogenic mutations or unique variants in the affected individuals that could explain the phenotypic differences between carriers of the same translocations. In conclusion, we suggest that NGS-based methods, such as WG-MPS, can be successfully used for detailed mapping of translocation breakpoints, which can also be used in routine clinical investigation of ABT cases. Unlike de novo translocations, no associations were determined here between familial two-way ABTs and the phenotype of the affected members, in which the presence of cryptic imbalances and complex chromosomal rearrangements has been excluded. Future whole-exome or whole-genome sequencing will potentially reveal unidentified mutations in the patients underlying the discordant phenotypes within each family. In addition, larger studies are needed to determine the exact percentage for phenotypic risk in families with ABTs.

Contrasting origin of B chromosomes in two cervids (Siberian roe deer and grey brocket deer) unravelled by chromosome-specific DNA sequencing

BACKGROUND

B chromosomes are dispensable and variable karyotypic elements found in some species of animals, plants and fungi. They often originate from duplications and translocations of host genomic regions or result from hybridization. In most species, little is known about their DNA content. Here we perform high-throughput sequencing and analysis of B chromosomes of roe deer and brocket deer, the only representatives of Cetartiodactyla known to have B chromosomes.

RESULTS

In this study we developed an approach to identify genomic regions present on chromosomes by high-throughput sequencing of DNA generated from flow-sorted chromosomes using degenerate-oligonucleotide-primed PCR. Application of this method on small cattle autosomes revealed a previously described KIT gene region translocation associated with colour sidedness. Implementing this approach to B chromosomes from two cervid species, Siberian roe deer (Capreolus pygargus) and grey brocket deer (Mazama gouazoubira), revealed dramatically different genetic content: roe deer B chromosomes consisted of two duplicated genomic regions (a total of 1.42-1.98 Mbp) involving three genes, while grey brocket deer B chromosomes contained 26 duplicated regions (a total of 8.28-9.31 Mbp) with 34 complete and 21 partial genes, including KIT and RET protooncogenes, previously found on supernumerary chromosomes in canids. Sequence variation analysis of roe deer B chromosomes revealed a high frequency of mutations and increased heterozygosity due to either amplification within B chromosomes or divergence between different Bs. In contrast, grey brocket deer B chromosomes were found to be more homogeneous and resembled autosomes in patterns of sequence variation. Similar tendencies were observed in repetitive DNA composition.

CONCLUSIONS

Our data demonstrate independent origins of B chromosomes in the grey brocket and roe deer. We hypothesize that the B chromosomes of these two cervid species represent different stages of B chromosome sequences evolution: probably nascent and similar to autosomal copies in brocket deer, highly derived in roe deer. Based on the presence of the same orthologous protooncogenes in canids and brocket deer Bs we argue that genomic regions involved in B chromosome formation are not random. In addition, our approach is also applicable to the characterization of other evolutionary and clinical rearrangements.

X-linked intellectual disability related genes disrupted by balanced X-autosome translocations

Detailed molecular characterization of chromosomal rearrangements involving X-chromosome has been a key strategy in identifying X-linked intellectual disability-causing genes. We fine-mapped the breakpoints in four women with balanced X-autosome translocations and variable phenotypes, in order to investigate the corresponding genetic contribution to intellectual disability. We addressed the impact of the gene interruptions in transcription and discussed the consequences of their functional impairment in neurodevelopment. Three patients presented with cognitive impairment, reinforcing the association between the disrupted genes (TSPAN7-MRX58, KIAA2022-MRX98, and IL1RAPL1-MRX21/34) and intellectual disability. While gene expression analysis showed absence of TSPAN7 and KIAA2022 expression in the patients, the unexpected expression of IL1RAPL1 suggested a fusion transcript ZNF611-IL1RAPL1 under the control of the ZNF611 promoter, gene disrupted at the autosomal breakpoint. The X-chromosomal breakpoint definition in the fourth patient, a woman with normal intellectual abilities, revealed disruption of the ZDHHC15 gene (MRX91). The expression assays did not detect ZDHHC15 gene expression in the patient, thus questioning its involvement in intellectual disability. Revealing the disruption of an X-linked intellectual disability-related gene in patients with balanced X-autosome translocation is a useful tool for a better characterization of critical genes in neurodevelopment. © 2015 Wiley Periodicals, Inc.

Genetic mechanisms leading to primary amenorrhea in balanced X-autosome translocations

OBJECTIVE

To map the X-chromosome and autosome breakpoints in women with balanced X-autosome translocations and primary amenorrhea, searching candidate genomic loci for female infertility.

DESIGN

Retrospective and case-control study.

SETTING

University-based research laboratory.

PATIENT(S)

Three women with balanced X-autosome translocation and primary amenorrhea.

INTERVENTION(S)

Conventional cytogenetic methods, genomic array, array painting, fluorescence in situ hybridization, and quantitative reverse transcription-polymerase chain reaction.

MAIN OUTCOME MEASURE(S)

Karyotype, copy number variation, breakpoint mapping, and gene expression levels.

RESULT(S)

All patients presented with breakpoints in the Xq13q21 region. In two patients, the X-chromosome breakpoint disrupted coding sequences (KIAA2022 and ZDHHC15 genes). Although both gene disruptions caused absence of transcription in peripheral blood, there is no evidence that supports the involvement of these genes with ovarian function. The ZDHHC15 gene belongs to a conserved syntenic region that encompasses the FGF16 gene, which plays a role in female germ line development. The break in the FGF16 syntenic block may have disrupted the interaction between the FGF16 promoter and its cis-regulatory element. In the third patient, although both breakpoints are intergenic, a gene that plays a role in the DAX1 pathway (FHL2 gene) flanks distally the autosome breakpoint. The FHL2 gene may be subject to position effect due to the attachment of an autosome segment in Xq21 region.

CONCLUSION(S)

The etiology of primary amenorrhea in balanced X-autosome translocation patients may underlie more complex mechanisms than interruption of specific X-linked candidate genes, such as position effect. The fine mapping of the rearrangement breakpoints may be a tool for identifying genetic pathogenic mechanisms for primary amenorrhea.

Nonallelic homologous recombination between retrotransposable elements is a driver of de novo unbalanced translocations

Large-scale analysis of balanced chromosomal translocation breakpoints has shown nonhomologous end joining and microhomology-mediated repair to be the main drivers of interchromosomal structural aberrations. Breakpoint sequences of de novo unbalanced translocations have not yet been investigated systematically. We analyzed 12 de novo unbalanced translocations and mapped the breakpoints in nine. Surprisingly, in contrast to balanced translocations, we identify nonallelic homologous recombination (NAHR) between (retro)transposable elements and especially long interspersed elements (LINEs) as the main mutational mechanism. This finding shows yet another involvement of (retro)transposons in genomic rearrangements and exposes a profoundly different mutational mechanism compared with balanced chromosomal translocations. Furthermore, we show the existence of compound maternal/paternal derivative chromosomes, reinforcing the hypothesis that human cleavage stage embryogenesis is a cradle of chromosomal rearrangements.

Characterising chromosome rearrangements: recent technical advances in molecular cytogenetics

Genomic rearrangements can result in losses, amplifications, translocations and inversions of DNA fragments thereby modifying genome architecture, and potentially having clinical consequences. Many genomic disorders caused by structural variation have initially been uncovered by early cytogenetic methods. The last decade has seen significant progression in molecular cytogenetic techniques, allowing rapid and precise detection of structural rearrangements on a whole-genome scale. The high resolution attainable with these recently developed techniques has also uncovered the role of structural variants in normal genetic variation alongside single-nucleotide polymorphisms (SNPs). We describe how array-based comparative genomic hybridisation, SNP arrays, array painting and next-generation sequencing analytical methods (read depth, read pair and split read) allow the extensive characterisation of chromosome rearrangements in human genomes.

Mapping of balanced chromosome translocation breakpoints to the basepair level from microdissected chromosomes

The analysis of structural variants associated with specific phenotypic features is promising for the elucidation of the function of involved genes. There is, however, at present no approach allowing the rapid mapping of chromosomal translocation breakpoints to the basepair level from a single chromosome. Here we demonstrate that we have advanced both the microdissection and the subsequent unbiased amplification to an extent that breakpoint mapping to the basepair level has become possible. As a case in point we analysed the two breakpoints of a t(7;13) translocation observed in a patient with split hand/foot malformation (SHFM1). The amplification products of the der(7) and of the der(13) were hybridized to custom-made arrays, enabling us to define primers at flanking breakpoint regions and thus to fine-map the breakpoints to the basepair level. Consequently, our results will also contribute to a further delineation of causative mechanisms underlying SHFM1 which are currently unknown.

Chromosome-specific segmentation revealed by structural analysis of individually isolated chromosomes

Analysis of structural rearrangements at the individual chromosomal level is still technologically challenging. Here we optimized a chromosome isolation method using fluorescent marker-assisted laser-capture and laser-beam microdissection and applied it to structural analysis of two aberrant chromosomes found in a lung cancer cell line. A high-density array-comparative genomic hybridization (array-CGH) analysis of DNA samples prepared from each of the chromosomes revealed that these two chromosomes contained 296 and 263 segments, respectively, ranging from 1.5 kb to 784.3 kb in size, derived from different portions of chromosome 8. Among these segments, 242 were common in both aberrant chromosomes, but 75 were found to be chromosome-specific. Sequences of 263 junction sites connecting the ends of segments were determined using a PCR/Sanger-sequencing procedure. Overlapping microhomologies were found at 169 junction sites. Junction partners came from various portions of chromosome 8 and no biased pattern in the positional distribution of junction partners was detected. These structural characteristics suggested the occurrence of random fragmentation of the entire chromosome 8 followed by random rejoining of these fragments. Based on that, we proposed a model to explain how these aberrant chromosomes are formed. Through these structural analyses, it was demonstrated that the optimized chromosome isolation method described here can provide high-quality chromosomal DNA for high resolution array-CGH analysis and probably for massively parallel sequencing analysis.

AID produces DNA double-strand breaks in non-Ig genes and mature B cell lymphomas with reciprocal chromosome translocations

Cancer-initiating translocations such as those associated with lymphomas require the formation of paired DNA double-strand breaks (DSBs). Activation-induced cytidine deaminase (AID) produces widespread somatic mutation in mature B cells; however, the extent of "off-target" DSB formation and its role in translocation-associated malignancy is unknown. Here, we show that deregulated expression of AID causes widespread genome instability, which alone is insufficient to induce B cell lymphoma; transformation requires concomitant loss of the tumor suppressor p53. Mature B cell lymphomas arising as a result of deregulated AID expression are phenotypically diverse and harbor clonal reciprocal translocations involving a group of Immunoglobulin (Ig) and non-Ig genes that are direct targets of AID. This group includes miR-142, a previously unknown micro-RNA target that is translocated in human B cell malignancy. We conclude that AID produces DSBs throughout the genome, which can lead to lymphoma-associated chromosome translocations in mature B cells.

Breakpoint analysis of balanced chromosome rearrangements by next-generation paired-end sequencing

Characterisation of breakpoints in disease-associated balanced chromosome rearrangements (DBCRs), which disrupt or inactivate specific genes, has facilitated the molecular elucidation of a wide variety of genetic disorders. However, conventional methods for mapping chromosome breakpoints, such as in situ hybridisation with fluorescent dye-labelled bacterial artificial chromosome clones (BAC-FISH), are laborious, time consuming and often with insufficient resolution to unequivocally identify the disrupted gene. By combining DNA array hybridisation with chromosome sorting, the efficiency of breakpoint mapping has dramatically improved. However, this can only be applied when the physical properties of the derivative chromosomes allow them to be flow sorted. To characterise the breakpoints in all types of balanced chromosome rearrangements more efficiently and more accurately, we performed massively parallel sequencing using Illumina 1G analyser and ABI SOLiD systems to generate short sequencing reads from both ends of DNA fragments. We applied this method to four different DBCRs, including two reciprocal translocations and two inversions. By identifying read pairs spanning the breakpoints, we were able to map the breakpoints to a region of a few hundred base pairs that could be confirmed by subsequent PCR amplification and Sanger sequencing of the junction fragments. Our results show the feasibility of paired-end sequencing of systematic breakpoint mapping and gene finding in patients with disease-associated chromosome rearrangements.

Array painting: a protocol for the rapid analysis of aberrant chromosomes using DNA microarrays

Array painting is a technique that uses microarray technology to rapidly map chromosome translocation breakpoints. Previous methods to map translocation breakpoints have used fluorescence in situ hybridization (FISH) and have consequently been labor-intensive, time-consuming and restricted to the low breakpoint resolution imposed by the use of metaphase chromosomes. Array painting combines the isolation of derivative chromosomes (chromosomes with translocations) and high-resolution microarray analysis to refine the genomic location of translocation breakpoints in a single experiment. In this protocol, we describe array painting by isolation of derivative chromosomes using a MoFlo flow sorter, amplification of these derivatives using whole-genome amplification and hybridization onto commercially available oligonucleotide microarrays. Although the sorting of derivative chromosomes is a specialized procedure requiring sophisticated equipment, the amplification, labeling and hybridization of DNA is straightforward, robust and can be completed within 1 week. The protocol described produces good quality data; however, array painting is equally achievable using any combination of the available alternative methodologies for chromosome isolation, amplification and hybridization.

A genome assembly-integrated dog 1 Mb BAC microarray: a cytogenetic resource for canine cancer studies and comparative genomic analysis

Molecular cytogenetic studies have been instrumental in defining the nature of numerical and structural chromosome changes in human cancers, but their significance remains to be fully understood. The emergence of high quality genome assemblies for several model organisms provides exciting opportunities to develop novel genome-integrated molecular cytogenetic resources that now permit a comparative approach to evaluating the relevance of tumor-associated chromosome aberrations, both within and between species. We have used the dog genome sequence assembly to identify a framework panel of 2,097 bacterial artificial chromosome (BAC) clones, selected at intervals of approximately one megabase. Each clone has been evaluated by multicolor fluorescence in situ hybridization (FISH) to confirm its unique cytogenetic location in concordance with its reported position in the genome assembly, providing new information on the organization of the dog genome. This panel of BAC clones also represents a powerful cytogenetic resource with numerous potential applications. We have used the clone set to develop a genome-wide microarray for comparative genomic hybridization (aCGH) analysis, and demonstrate its application in detection of tumor-associated DNA copy number aberrations (CNAs) including single copy deletions and amplifications, regional aneuploidy and whole chromosome aneuploidy. We also show how individual clones selected from the BAC panel can be used as FISH probes in direct evaluation of tumor karyotypes, to verify and explore CNAs detected using aCGH analysis. This cytogenetically validated, genome integrated BAC clone panel has enormous potential for aiding gene discovery through a comparative approach to molecular oncology.

Clonal evolution in t(14;18)-positive follicular lymphoma, evidence for multiple common pathways, and frequent parallel clonal evolution

PURPOSE

Follicular lymphoma typically has acquired a t(14;18) translocation, but subsequent additional cytogenetic abnormalities contribute to disease progression. The main aims of the study are to (a) identify the frequency and temporal sequence of cytogenetic events in t(14;18)-positive follicular lymphoma, (b) determine if there are specific pathways in the evolution of follicular lymphoma, (c) determine the clonal divergence in cases with sequential biopsies or multiple clones from a single biopsy, and (d) determine the association of genetic imbalances with clinical outcome.

EXPERIMENTAL DESIGN

All cases with a histologically confirmed diagnosis of follicular lymphoma and cytogenetic analysis showing t(14;18)(q32;q21) were included. The karyotypes were reviewed and cytogenetic data were entered into a relational database for further computational analysis; 418 biopsies from 360 follicular lymphoma patients including 43 sequential biopsies were analyzed.

RESULTS

Of the cases with only one or two genomic imbalances, the most frequent chromosomal imbalances were +7, del(6q), +der(18)t(14;18), +18, and +X. These abnormalities were also among the most frequent ones encountered when all karyotypes were analyzed. Cytogenetically abnormal clones in the same (26%) and sequential biopsies (63%) often showed divergence of genetic alterations. Balanced translocations other than the t(14;18) were uncommon events, but chromosomal breaks involving 14q32, 18q21, 1p36, 1q21, 10q22, 10q24, and a large cluster at 6q occurred relatively frequently. del(6q), +5, +19, and +20 were associated with poorer overall survival, and del(17p) was associated with poorer event-free survival. Lower-grade tumors (1 and 2) were associated with fewer imbalances.

CONCLUSION

Our analysis suggested that +der(18)t(14;18) may be an entry point to a distinct pathway of genetic evolution in follicular lymphoma. The other common early events appeared to provide multiple entry points, and they might cooperate in the pathogenesis and progression of the follicular lymphoma. Cytogenetically abnormal clones from same patients often showed divergence of genetic alterations, suggesting that parallel evolution from precursor clones are frequent events. This study provides the framework for further analysis of genetic pathways of tumor progression.

Estimating genome-wide copy number using allele-specific mixture models

Genomic changes such as copy number alterations are one of the major underlying causes of human phenotypic variation among normal and disease subjects. Array comparative genomic hybridization (CGH) technology was developed to detect copy number changes in a high-throughput fashion. However, this technology provides only a >30-kb resolution, which limits the ability to detect copy number alterations spanning small regions. Higher resolution technologies such as single nucleotide polymorphism (SNP) microarrays allow detection of copy number alterations at least as small as several thousand base pairs. Unfortunately, strong probe effects and variation introduced by sample preparation procedures have made single-point copy number estimates too imprecise to be useful. Various groups have proposed statistical procedures that pool data from neighboring locations to successfully improve precision. However, these procedure need to average across relatively large regions to work effectively, thus greatly reducing resolution. Recently, regression-type models that account for probe effects have been proposed and appear to improve accuracy as well as precision. In this paper, we propose a mixture model solution, specifically designed for single-point estimation, that provides various advantages over the existing methodology. We use a 314-sample database, to motivate and fit models for the conditional distribution of the observed intensities given allele-specific copy number. We can then compute posterior probabilities that provide a useful prediction rule as well as a confidence measure for each call. Software to implement this procedure will be available in the Bioconductor oligo package (www.bioconductor.org).

Array painting reveals a high frequency of balanced translocations in breast cancer cell lines that break in cancer-relevant genes

Chromosome translocations in the common epithelial cancers are abundant, yet little is known about them. They have been thought to be almost all unbalanced and therefore dismissed as mostly mediating tumour suppressor loss. We present a comprehensive analysis by array painting of the chromosome translocations of breast cancer cell lines HCC1806, HCC1187 and ZR-75-30. In array painting, chromosomes are isolated by flow cytometry, amplified and hybridized to DNA microarrays. A total of 200 breakpoints were identified and all were mapped to 1 Mb resolution on bacterial artificial chromosome (BAC) arrays, then 40 selected breakpoints, including all balanced breakpoints, were further mapped on tiling-path BAC arrays or to around 2 kb resolution using oligonucleotide arrays. Many more of the translocations were balanced at 1 Mb resolution than expected, either reciprocal (eight in total) or balanced for at least one participating chromosome (19 paired breakpoints). Second, many of the breakpoints were at genes that are plausible targets of oncogenic translocation, including balanced breaks at CTCF, EP300/p300 and FOXP4. Two gene fusions were demonstrated, TAX1BP1-AHCY and RIF1-PKD1L1. Our results support the idea that chromosome rearrangements may play an important role in common epithelial cancers such as breast cancer.

Mapping translocation breakpoints by next-generation sequencing

Balanced chromosome rearrangements (BCRs) can cause genetic diseases by disrupting or inactivating specific genes, and the characterization of breakpoints in disease-associated BCRs has been instrumental in the molecular elucidation of a wide variety of genetic disorders. However, mapping chromosome breakpoints using traditional methods, such as in situ hybridization with fluorescent dye-labeled bacterial artificial chromosome clones (BAC-FISH), is rather laborious and time-consuming. In addition, the resolution of BAC-FISH is often insufficient to unequivocally identify the disrupted gene. To overcome these limitations, we have performed shotgun sequencing of flow-sorted derivative chromosomes using "next-generation" (Illumina/Solexa) multiplex sequencing-by-synthesis technology. As shown here for three different disease-associated BCRs, the coverage attained by this platform is sufficient to bridge the breakpoints by PCR amplification, and this procedure allows the determination of their exact nucleotide positions within a few weeks. Its implementation will greatly facilitate large-scale breakpoint mapping and gene finding in patients with disease-associated balanced translocations.

Methods and strategies for analyzing copy number variation using DNA microarrays

The association of DNA copy-number variation (CNV) with specific gene function and human disease has been long known, but the wide scope and prevalence of this form of variation has only recently been fully appreciated. The latest studies using microarray technology have demonstrated that as much as 12% of the human genome and thousands of genes are variable in copy number, and this diversity is likely to be responsible for a significant proportion of normal phenotypic variation. Current challenges involve developing methods not only for detecting and cataloging CNVs in human populations at increasingly higher resolution but also for determining the association of CNVs with biological function, recent human evolution, and common and complex human disease.

Getting it right: designing microarray (and not 'microawry') comparative genomic hybridization studies for cancer research

The development of high-resolution microarray-based comparative genomic hybridization (aCGH), using cDNA, bacterial artificial chromosome (BAC) and oligonucleotide probes, is providing tremendous opportunities for translational research by facilitating detailed analysis of entire cancer genomes in a single experiment. However, this technology will only fulfil its promise if studies incorporating aCGH are designed with a full understanding of its current limitations and the strategies available to circumvent them. While there have been several excellent reviews on the current status of this technology, there is currently very little guidance available regarding the appropriate design of experiments incorporating aCGH (including the strengths and weaknesses of each platform), and how best to combine the results obtained from aCGH with other 'omic' technologies, including gene expression. In this review, we present the key design issues that need to be considered in order to optimize aCGH studies, including sample selection, the definition of appropriate experimental objectives, arguments for and against the various microarray platforms that are currently available, and methods for data validation and integration. It is envisaged that future well-designed aCGH studies will enhance our understanding of the genetic basis of cancer, and lead to the identification of novel predictive and prognostic cancer biomarkers, as well as molecular therapeutic targets in cancer.

Mapping translocation breakpoints using a wheat microarray

We report mapping of translocation breakpoints using a microarray. We used complex RNA to compare normal hexaploid wheat (17,000 Mb genome) to a ditelosomic stock missing the short arm of chromosome 1B (1BS) and wheat-rye translocations that replace portions of 1BS with rye 1RS. Transcripts detected by a probe set can come from all three Triticeae genomes in ABD hexaploid wheat, and sequences of homoeologous genes on 1AS, 1BS and 1DS often differ from each other. Absence or replacement of 1BS therefore must sometimes result in patterns within a probe set that deviate from hexaploid wheat. We termed these 'high variance probe sets' (HVPs) and examined the extent to which HVPs associated with 1BS aneuploidy are related to rice genes on syntenic rice chromosome 5 short arm (5S). We observed an enrichment of such probe sets to 15-20% of all HVPs, while 1BS represents approximately 2% of the total genome. In total 257 HVPs constitute wheat 1BS markers. Two wheat-rye translocations subdivided 1BS HVPs into three groups, allocating translocation breakpoints to narrow intervals defined by rice 5S coordinates. This approach could be extended to the entire wheat genome or any organism with suitable aneuploid or translocation stocks.