Complex reorganization and predominant non-homologous repair following chromosomal breakage in karyotypically balanced germline rearrangements and transgenic integration

We defined the genetic landscape of balanced chromosomal rearrangements at nucleotide resolution by sequencing 141 breakpoints from cytogenetically-interpreted translocations and inversions. We confirm that the recently described phenomenon of “chromothripsis” (massive chromosomal shattering and reorganization) is not unique to cancer cells but also occurs in the germline where it can resolve to a karyotypically balanced state with frequent inversions. We detected a high incidence of complex rearrangements (19.2%) and substantially less reliance on microhomology (31%) than previously observed in benign CNVs. We compared these results to experimentally-generated DNA breakage-repair by sequencing seven transgenic animals, and revealed extensive rearrangement of the transgene and host genome with similar complexity to human germline alterations. Inversion is the most common rearrangement, suggesting that a combined mechanism involving template switching and non-homologous repair mediates the formation of balanced complex rearrangements that are viable, stably replicated and transmitted unaltered to subsequent generations.


INTRODUCTION
Our understanding of the genetic architecture of human chromosomal rearrangements has expanded in recent years, as rapid improvements in genomics technology have spawned a growing number of mechanistic hypotheses, many involving some degree of homology between participant sequences [1][2][3] . For complex events, replication-based mechanisms have been proposed such as template switching from a stalled or disrupted replication fork (Fork-Stalling and Template Switching, FoSTeS) 4 and microhomology-mediated break-induced replication (MMBIR) 3,5 . While it has been established that chromosomal exchanges which appear to be balanced at lower-resolution can actually involve considerable complexity that may contribute to human disease in unexpected ways 6,7 , few studies have assessed these events at base-pair resolution. The first massively-parallel sequencing of cancer cells suggested a complex rearrangement landscape 8 , and more recently an astonishing phenomenon was uncovered in cancer cells that included massive chromosomal shattering and rearrangement, with frequent change in copy number state across the region 9 . The authors dubbed this "chromothripsis" and hypothesized it to be a unique feature of somatic mutation that might occur in 2-3% of all cancers 9 , with similar complexity seen in several subsequent sequencing studies [10][11][12] .
The mutational mechanism(s) underlying such complex genomic reorganization is unknown, but an intriguing feature of cancer-related chromothripsis was a frequent transition between two copy number states throughout the impacted region. As copy number variants (CNVs) have emerged as a major component of genetic variation in humans 13 , two recent population-based sequencing studies found a high degree of microhomology at CNV breakpoints and suggested two predominant CNV generating mechanisms: microhomologymediated end joining (MHEJ) and non-homologous end joining (NHEJ) 14,15 , similar events postulated to occur by distinct pathways 16,17 . These mechanistic hypotheses may be somewhat limited as they were restricted to events defined by DNA dosage changes in unphenotyped individuals. The rearrangement landscape of other forms of genomic rearrangements, namely constitutional balanced structural variation (SV) such as reciprocal translocation and inversion, have not been comprehensively assessed at the sequence level. Breakpoint resolution of these SVs is fundamental to the prediction of which rearrangement mechanisms could underlie their formation and to our understanding of the full range of mutational mechanisms involved in human SV.
Here, we provide the first high-throughput, sequence-based assessment of mutational mechanisms associated with breakpoints from 52 subjects with chromosomal abnormalities that were previously defined at cytogenetic resolution as balanced (45 reciprocal translocations and seven inversions) and clinically assessed as likely to be pathogenic (50 arose de novo and two were inherited from an affected parent). We performed a series of next-generation sequencing experiments that included either whole-genome sequencing or targeted capture of breakpoints (see Supplementary Methods and Talkowski et al. 18,19 ) and all breakpoints were confirmed at base-pair resolution by capillary sequencing. Our results show definitively that, as in cancer, remarkably complex genomic reorganization can also occur in the human germline, but that the repair process can resolve to a relatively balanced state rather than yielding extensive gains and losses of DNA. When we surveyed an experimental system of chromosomal rearrangement (transgenic animals), we found that similar genomic reorganization can result from experimentally-generated double-stranded DNA breaks (DSB), in the absence of environmental mutagenic factors and in the presence of an abundance of homologous template. In both human and transgenic animal breakpoints studied, the mechanism(s) that mediated rearrangements did not depend primarily on microhomology, nor were they frequently associated with large DNA dosage changes; both findings in stark contrast to previous studies of chromothripsis in cancer and benign CNV formation. Instead, these results reveal substantial chromosomal reorganization in the germline that can be compatible with viability, stably replicated, and transmitted to subsequent generations.

Complex genomic reorganization in the human germline
We found cytogenetically defined and apparently balanced SVs to be far more complex than originally thought, detecting 141 breakpoints, an average of 2.71 breakpoints per subject. Only two subjects had two derivative chromosomes with no DNA imbalance, but most subjects also did not suffer a substantial loss of genetic material (arbitrarily defined here as >1 kb total genomic imbalance, Table 1, Table S1). Instead, multiple breakpoints were often reassembled in a relatively balanced manner with microdeletions of one to several bases. We found that 19.2% of all karyotypically balanced SVs actually met conventional criterion of a complex chromosomal rearrangement (CCR; three or more breakpoints), a significantly higher frequency than the previous estimate of 2.8% from 246 de novo balanced anomalies assessed from more than 269,000 prenatal diagnoses (p = 2.2 × 10 −4 ) 7 . Nine of the CCRs involved at least one inverted segment; only two CCRs involved exchanges between more than two chromosomes.
Among the 10 CCR subjects, there was extensive complexity of rearrangement. In the two most intricate cases (BSID42 and BSID43), we observed shattering and reorganization of multiple chromosomal segments similar to chromothripsis seen in cancer, except that these germline events resolved to a largely dosage balanced state rather than exhibiting frequent alterations in copy number 9 (Fig. 1, Fig. 2; Supplementary Table 2). For BSID42, we observed 14 different junctions between chromosome (chr) X and chr5, with extensive shattering of localized regions (Fig. 2a, Supplementary Video 1). The 5q14.3 DNA 'shards' integrated into chrX with frequent oscillations between inverted and same-strand orientation (Fig. 2a). In both derivative chromosomes, large segments remained fully intact between reassembled DNA fragments. Prior to the most distal telomeric 5q-Xq junction in the der(x), an ~871 kb segment of chrX (designated Xc) was excised, inverted, and inserted into the strand-oscillating chr5 segment approximately 20 Mb from its chrX origin. Similarly, an 8 kb segment of chr5 (designated 5f) was also excised, inverted, and reinserted in der (5). Sequencing thus delineated a fascinating composite wherein shattered exchanges between chromosomes also resulted in extensive intrachromosomal reorganization. Notably, both the clinical karyotype and a 1 million feature aCGH experiment were completely blind to this series of rearrangements that extensively reorganized more than 23.5 Mb of DNA (0.76% of the entire genome for this SV; Table S2), as the sum of DNA imbalance from all junction fragments was just 6,357 bases.
The second sample, BSID43, had DNA available from two sources, one extracted from whole blood and another from an EBV-transformed lymphoblastoid cell line, both of which were sequenced (214.6X average physical coverage of mapped inserts, Table S2). The fully resolved genomic reorganization involved 11 breakpoints resulting from the apparent shattering of a 3 Mb region of 7q31.31-q31.32 that integrated into the long arms of both chr3 and chr5. The chr3 breakpoints spanned 3.5 Mb (3q25.2-q25.31), whereas chr5 had just a single breakpoint on each derivative (5q14.3) with the loss of only 16 bases between derivatives (Fig 2b, Supplementary Video 2). Despite the karyotypic interpretation of a reciprocal translocation of chr3 and chr5, there were no actual junctions connecting 3q25.2-25.31 to 5q14.3. Rather, all breakpoint junctions detected were joined to shattered fragments of 7q31.31-q31.32, and all events were again cryptic to both karyotype and 1M aCGH diagnostics (a total of only 1,551 bases was lost from all rearrangements) (Fig 2b). All breakpoints were identical between the two DNA sources, indicating that, once established, the rearranged chromosomal organization was stable through EBV transformation and subsequent cell propagation.
Eight additional subjects had CCRs ranging from 3-6 breakpoints each, many of which involved less than 1 kb of total genomic imbalance. All but one had inverted segments associated with a breakpoint, usually directly at the breakpoint junction, with excision, inversion, and insertion events as a common signature (Tables S1 and S2). The only CCR event occurring without an inversion associated with a breakpoint (BSID45) contained a balanced excision of 13.5 Mb of chr4 and reinsertion 5.3 Mb away from the der(13) translocation breakpoint (Tables S1 and S2). Of the seven karyotypically-defined inversions, six were simple and largely balanced events and only 1 event (BSID39) was complex, containing a smaller 4.1 kb inversion nested within the larger 43.7 Mb inversion (Tables S1 and S2).

Complex rearrangements from experimental transgenesis
Each of the previous examples of cancer chromothripsis postulated that exogenous factors such as ionizing radiation, reactive oxygen species, or chemotherapy could prime such rearrangements by inducing multiple DSBs (reviewed by Khanna & Jackson 20 ). One previous study assessed DSB repair between germline and experimentally generated, aphidicolin-induced CNV formation, finding non-homologous repair as the driving mechanism in both 21 . The proposition that the complex reorganization of human and cancer genomes operates through aberrant repair of multiple, apparently simultaneous DNA breaks prompted us to test whether a similar outcome is experimentally induced during transgene integration, where exogenous DNA fragments, which are themselves substrates for DSB repair, are incorporated into the genome. We hypothesized that the process of transgene integration might provide insight into the general mechanisms involved in complex rearrangements. Interestingly, few such transgene integration sites have been characterized at the sequence level 22 . Using two sequencing approaches (targeted capture of transgenes and whole-genome jumping libraries), we analyzed seven transgenic animals: two wellestablished lines of mice produced by injection of a 1.9 kb fragment of human 4p16.3 genomic DNA containing exon 1 of the Huntington's disease gene (HTT) (R6/1 and R6/2) 23 and five HTT sheep transgenic lines created by injection of an 11.6 kb full-length HTT cDNA (G 0 /1, G 0 /2, G 0 /5, G 0 /4 and G 0 /6) 24 (see Methods and Supplemental Information).
Transgene integration sites were identified to base-pair resolution in all seven animals, revealing abundant complex rearrangements in both the internal structure of the final transgene (i.e., deletion, duplication, inversion, and excision/inversion/insertion) and in the host genomic DNA. Inverted segments were a common feature, displaying frequent oscillations in DNA orientation similar to the human CCRs described above (Fig. 3). At least two inverted excision/insertion events were detected in each sheep line with many containing significantly more such events. In the R6/2 line, we also found insertion of a 168 bp segment of foreign DNA with high homology to Corynebacteria. Repair of the DNA breaks commonly resulted in rearrangement of the host genome at the site of integration as well as at remote loci. Genomic deletions were detected in four transgenic animals, the largest of which was a 5.4 kb deletion in R6/2. In G0/2, we also found an 81 kb inversion at the integration site, while in G0/6, an apparent shattering and complex reintegration of the injected DNA into six different chromosomes was observed, including an interchromosomal rearrangement between chromosomes OAR7, OAR8, and OAR15 (Fig. S1). This rearrangement involved a small interchromosomal excision/insertion of OAR15 into an apparent translocation junction between OAR7, transgene sequence, and OAR8, reminiscent of the complex rearrangement seen in subject BSID43.
These complex rearrangements mirrored the human germline CCRs, suggesting that shared or overlapping mechanisms of DSB repair mediate formation of CCRs. Arguably, the study of transgenic animals has the potential for two independent breakpoint classes: (1) the sites of integration and rearrangements of the host genome, and (2) rearrangements of the injected fragment itself. Having established the gross similarities of rearrangements in the human germline and in transgenic animals, we proceeded to analyze the sequence features of our balanced SVs and to compare these with the breakpoint sequences of human CNVs and the transgenic animals.

Mechanistic signatures differ between balanced structural variation and copy number variants
We discovered a dramatic disparity between breakpoint features of the de novo and presumably pathogenic translocations or inversions and the results obtained from population-based CNV studies, with important mechanistic consequences. Only 30.5% of these 141 breakpoints were flanked by microhomology, suggesting that most events were unlikely to have arisen by MHEJ, homology mediated repair, or microhomology-mediated replication-based mechanisms such as MMBIR 3,5 (Fig. 4a). This distribution resembled the profile of six tumors sequenced by Stephens et al. (2010) that had undergone chromothripsis (45% of breakpoints contained microhomology; Fig. 4b). By comparison, from recent studies 70%-80% of all CNVs in unphenotyped individuals contained microhomology 14,15 (Fig. 4c, 4d).
The previous CNV estimates were almost exclusively based upon deletions as these events were most readily detected by the previous CNV capture of  and the low-depth sequencing of Mills et al. (2011); neither study evaluated balanced events. Kidd et al. (2010) surveyed SVs at fosmid resolution and suggested that NAHR was the dominant mechanism mediating large inversion formation in population-based samples (69% had flanking homologous segments of > 200 bp) 25 . In our study of potentially pathogenic balanced SVs, we found that just one of the 43 total inverted segments contained greater than 200 bp of homology (654 bp) (Fig. S3). This one instance was not a karyotypically defined large inversion but rather a small 2.3 kb inversion associated with a CCR (Table  S1). Collectively, the microhomology profile of the karyotypically-defined inversions mirrored the breakpoints of the translocations, 69% of which contained < 2 bp of microhomology. The disparity between these presumably pathogenic inversions and the variants analyzed by Kidd et al. (2010) 25 was not biased by sequencing methodology (NAHR is a notoriously difficult mechanism to delineate by short-read sequencing) as we resolved each of the seven karyotypically-defined de novo inversions to nucleotide resolution, none of which contained long stretches of homology at the breakpoint.
We surveyed breakpoint profiles in the transgenic animals and found that the microhomology associated with this experimentally induced DSB repair provided compelling evidence for mechanisms more similar to de novo balanced SVs than to benign CNVs, as 78% of the host genomic DNA integration sites contained no microhomology. To test sensitivity of the paired-end sequencing results, we performed capillary sequencing validation of the entire transgene sequence in two animals (R6/2 and G0/1, excluding 121 bp of a hairpin formation in G0/1), which confirmed all of the rearrangement junctions detected by paired-end sequencing in Fig. 3. We also found that 76.9% of the internal transgene junctions had less than 2 bp of homology, despite the presence of abundant homologous template for homology mediated mechanisms.
We further explored this striking difference in the frequency of microhomology between de novo, balanced SVs tested here and previous CNV studies. We initially tested for methodological differences using three independent methods to assess microhomology: (1) a high-throughput alignment using BWA Smith-Watherman (BWA-SW; Fig. S2), resulting in 26.2% of breakpoints containing two or more bases of microhomology. (2) The EMBOSS needle program as described by Kidd et al. (2010), which allows concurrent inserted template sequences with microhomology, identifying microhomology in seven additional individuals (31.2%), which was not statistically different from the BWA-SW procedure (p = 0.43). (3) A revised version of BreakSeq 26 previously used in Mills et al. (2011), that we customized to analyze balanced SVs with nearly identical results to BWA-SW and EMBOSS (R. Mills, unpublished, see Supplemental Information). We therefore considered both the final breakpoint junctions and the sequence features of the initially intact chromosomes prior to their disruption in our analyses with consistent conclusions between methods. When we compared these data to the breakpoint homology from 16,783 CNV breakpoints assembled for the 1,000 Genomes Pilot 1 study using identical methods, we find a highly significant difference from either the BWA-SW analysis (χ 2 = 244.6, 1 d.f., p = 2.6×10−55) or the EMBOSS method (χ 2 = 201.6, 1 d.f., p = 1.26×10−45) (Fig. 3). The EMBOSS results are reported in the Tables and Figures as the most conservative interpretation of the distinction between previously published population-based studies and our results in karyotypically-defined balanced rearrangements.
This persistent deficiency in microhomology led us to question the extent to which any mechanisms other than random NHEJ contributed to formation of chromosomal rearrangements in these subjects. We generated 1,000 simulated chimeric sequences of random breakpoints in the genome and analyzed the microhomology at these random junctions relative to the breakpoints in our subjects. The overall distributions of microhomology were significantly different (χ 2 = 24.5, 2 d.f., p = 4.88×10 −6 ). When we scrutinized these data we found a close concordance between the experimental set and random breakpoints for exactly two bp of homology (12.1% compared to 10.7%, respectively; p = 0.66). Of the remaining breakpoints (< 2 bp or > 2 bp homology), there was a significant enrichment of microhomology in the experimental set compared to the random breakpoints (experimental vs. simulated for < 2 bp = 68.8% vs. 82.3%, p = 3.6×10 −4 ; > 2 bp = 19.1% vs. 7.0%; p = 1.38×10 −5 ) (Fig. 5a), indicating a non-random reliance on a microhomology-mediated repair mechanism in the formation of these balanced SVs. In comparison with simulations of 10,000 datasets of equal size to our experimental set, we found a marginal enrichment of LINEs at the breakpoints (empirical p = 0.085) but not SINEs (Fig. 5b, 5c), and with 1 million simulations of windows ranging from 2 bp -500 bp surrounding random genomic breakpoints, we found no significant enrichment in frequency of repetitive DNA motifs that could fold into non-B DNA structure 27,28 (Table  S3).

DISCUSSION
This sequence-level evaluation of the genetic architecture of karyotypically balanced, potentially pathogenic chromosomal rearrangements reveals several remarkable and underappreciated features of such events: (1) karyotypically balanced SVs are rarely truly balanced at the primary sequence level; (2) human germline rearrangements can be extraordinarily complex, including localized deconstruction of chromosomal segments into many small fragments that can be rejoined and resolved in an aberrant but largely dosagebalanced manner that is compatible with life; (3) local inversions and what we refer to as 'inverted translocations' (inversions associated with translocation breakpoints) are a commonly observed signature of CCRs; and (4) the breakpoint characteristics of potentially pathogenic, karyotypically balanced SVs are markedly different from those of apparently benign CNVs. These studies also provide the first high-throughput sequencing assessment of the complexity resulting from integration of linear exogenous DNA in transgenic animals introduced by pronuclear injection, revealing resolved breakpoints that in many ways mimic those of human germline and cancer chromothripsis 9,11,12 . These similarities suggest that the mechanisms involved in establishing transgene integrations in model organisms may overlap with those that produce complex genomic DNA rearrangement in humans.
Cytogenetic studies previously established that a small subset of translocations involve multiple chromosomal breakpoints rather than simple reciprocal exchanges between chromosomes 6,29 . By sequencing karyotypically balanced SVs, we find a substantially greater proportion CCRs than previously estimated. Sequence resolution enabled us to expand the definition of CCR to include local inversions at translocation breakpoints (dubbed 'inverted translocations' here) which appear to be an underappreciated feature of CCRs and among the most pervasive characteristics of germline chromothripsis, which displayed frequent switching of strand orientation between integrated fragments. This is consistent with cancer chromothripsis 9 , as well as the finding that 4.5% of sequenced CNVs contained inserted sequence from the opposite strand in close proximity to the breakpoint in one study 14 , and small inversions found at the breakpoints of constitutional CNVs in an independent study 21 . We further observed that inverted excision/insertions in CCRs are frequently accompanied by additional genomic rearrangements, and that these can be interspersed with large, fully intact segments, suggesting an association between local microinversion and aberrant repair of DNA breaks.
Our findings confirm and further illuminate the phenomenon of chromothripsis resulting from highly aberrant DSB repair, however they foster further uncertainty about the cause of the initial event and the mechanism of repair. Our results strengthen the conclusion of Stephens et al. (2011) 9 that these events occur in a one-off chromosomal catastrophe rather than by progressive accumulation, as we find the same phenomenon in the human germline and in the first generation transgenic sheep. We also confirm in humans and transgenic animals that, once resolved, the events can be viable, replicate during mitosis with high fidelity, and be stably transmitted to subsequent generations. The cause of chromothripsis thus does not appear to be specific to aberrant repair pathways in cancer, nor is it likely driven by other mechanisms considered by Stephens et al. (2011) 9 such as disruption of tumor suppressor genes, formation of double minutes, or telomere disruption from end-toend chromosome fusion. Although cancer and germline chromothripsis appear to differ in the presence of extensive copy number change, dosage alterations may be favored in cancer cells due to dysregulation of growth pathways via loss of tumor suppressor genes. Frequent dosage alterations likely also occur in the germline, as one recent study found a surprisingly high number of segmental imbalances in cleavage stage embryos (~70%) 30 . At least a portion of these complex events are viable, as evidenced by Lupski and colleagues in a study published during review of this manuscript that selected for human cases with multiple copy number imbalances on a single chromosome 31 . It is possible that the co-occurrence of dosage alteration with substantial inter-chromosomal exchanges result in a non-viable embryo or fetus, eliminating such cases before birth. Thus, the relatively balanced dosage of individuals with CCRs may simply reflect a selection for viability. With the increasing application of aCGH in research and clinical diagnostics, additional human cases with extensive copy number alterations will undoubtedly emerge and it will be interesting to determine whether they display the high frequency of concomitant chromosomal reorganization seen here and in cancer chromothripsis.
It is unknown if the NAHR-dominated mechanism of large inversions observed previously by fosmid sequencing 25 is a characteristic of recurrent events as inversions were not assessed in other population-based sequencing studies. One recent study of unbalanced, recurrent translocations found predictable microhomology at the recurrent breakpoints 32 . In contrast, accumulating evidence indicates that formation of both common and rare CNVs is predominantly driven by non-homologous repair mechanisms that are reliant on much smaller degrees of homology between segments than NAHR 14,15,25 . Our results align with those of Mills et al. (2011) in supporting non-homologous repair as the favored mechanism in formation of deletions and the non-recurrent balanced rearrangements sequenced here; however, the level of microhomology associated with breakpoint signatures of these nonrecurrent de novo balanced rearrangements was distinctly lower than those previous studies. There was also a dearth of sequence homology at the sites of integration in the transgenic animals and the chromothripsis seen in cancer cells based on data obtained from Stephens et al. (2011) 9 . The distinguishing feature of all three groups (balanced SVs, cancer cells, and transgenic animals) is that they presumably all contribute to an abnormal phenotype. In this regard, it may not be surprising that this extreme form of aberrant DNA breakage repair proposes significantly different mechanisms of formation than presumably benign CNVs in the general population. Our data support a version of NHEJ for most exchanges of genetic material between chromosomes in the absence of microhomology and without significant alteration in copy number. However, this mechanism as currently understood is unlikely to account for the frequent strand inversions observed, which instead imply a model of template switching, such as FoSTeS or MMBIR in the formation of CCRs including chromothripsis 3,5 . Spatially localized DNA breaks due to a single DNA insult could also drive the rearrangements seen here, as the complex CCRs may require multiple simultaneous breaks. These data might also implicate operation of multiple mechanisms, where reorganization of shattered chromosomes to a largely balanced state is driven by blunt end-joining with some degree of template switching and homology mediated repair, or by an alternative and yet undefined mechanism that has the capacity to derive inverted segments frequently, which have been seen in all instances of chromothripsis. Such a model could account for CCR formation in the germline aberrations as well as those seen in cancer cells.
Our simulations indicate this possibility as we saw a modest but significant enrichment of microhomology in comparison to randomly selected breakpoints, which increased with the complexity of the SV.
These analyses provide an initial sequence-based survey of karyotypically balanced chromosomal rearrangements that, like the previous findings of chromothripsis in cancer, set the stage for extensive further investigation. Yet, our results already have a potentially wide-ranging impact as they indicate that traditional and current clinical diagnostic methods, including chromosomal microarrays, can be insufficient to understand the true nature of genomic disruption in affected patients. They also raise questions about whether the mechanisms suggested in CNV studies of unphenotyped populations are fully generalizable to other structural genomic alterations, and they suggest that more complex mechanistic hypotheses are required to explain how the genome can dynamically remodel itself in a surprisingly extensive manner to produce a stable, viable, complex constitutional rearrangement that can contribute to altered human development and disease.

Supplementary Material
Refer to Web version on PubMed Central for supplementary material. Three plots are shown for the two most intricately rearranged balanced SV subjects (BSID42 and BSID43) and the multiple fragmented integrations associated a transgenic animal (G 0 /6). In each plot, lines connect each of the inter and intra-chromosomal junction fragments sequenced and the chromosomes are labeled outside of the circle. 1A. BSID42: Fourteen junctions were confirmed between chromosomes 5 and X. 1B. BSID43: Eleven junction fragments were confirmed between chromosomes 3, 5, and 7. 1C. G 0 /6: Sequencing in a transgenic sheep revealed that the transgene was apparently 'shattered' and integrated into six different chromosomes, including one interchromosomal excision of sheep chromosome 8 (OAR8) and insertion into the junction fragment between OAR7 and the transgene, mimicking a translocation junction. All junction fragments were confirmed by capillary sequencing. Additional validation was performed using FISH to further confirm these results (Supplementary Material). All positions are hg19. Two independent karyotypic analyses indicated a balanced reciprocal translocation between chromosomes 3q and 5q; however, sequencing revealed the shattering of chromatin from 7q and re-integration of 7q DNA shards into the junction fragments of both derivative chromosomes, resulting in no direct 3q-5q junctions. There were four different inverted excision/insertion events, including intrachromosomal excision, inversion, and re-insertion in 7q with a co-occurring inverted insertion of an intact 3.5 Mb segment of 3q, all of which involved only 1,551 bases of total DNA imbalance. All positions are hg19. See also Supplementary Movies 1 and 2.

Figure 3. Complex rearrangements in transgenic animals
Results from sequencing of two transgenic mice (A) and five transgenic sheep (B). Each transgene prior to pronuclear injection is provided in the first line and resultant transgene integration site and internal structure is shown below for each animal. Genomic integration sites are given by the gray lines flanking the transgene and the most parsimonious transgene structure is shown (R6/1,2 mice and G 0 /1,2,4,5,6 sheep). The fragment injected into the R6 lines (1,905 bp) comprised the human HTT promotor, exon 1 (~130 CAG repeats), and 264 bp of human intron 1. The fragment injected into the sheep (11,625 bp) included the human HTT promoter and cDNA (69 CAG repeats), followed by exon 4, intron 4, and exon 5 of the bovine growth hormone (BGH) gene. Arrows represent strand orientation. *This transgene was fragmented and inserted into multiple chromosomal locations; one of the small fragments and resultant complex rearrangements of the host genome involving three independent chromosomes is shown here (see also Fig. 1, Supplementary Fig. 1). **Headto-tail junctions and read depth analyses indicate multiple copies of the segment, however the precise number of duplications could not be determined. Mouse positions are mm9 and sheep positions are OAR2.   Percentages represent percentage of total events (either by subject, n=52, or by breakpoint, n=141) in each category.
Results presented from analyses using EMBOSS needle for direct comparison with Kidd et al. (2010). § Six of the seven karyotypically defined inversions involved just two breakpoints; one involved multiple inversion breakpoints. * A single inverted breakpoint junction contained 654 bp of homology, no other junctions had more than 20 bp of homology. Deletions restricted to 127 breakpoints for which breakpoint pairs could be unambiguously determined so the precise loss of DNA could be calculated.
Nat Genet. Author manuscript; available in PMC 2012 October 01.