Abstract

Background

Structural variants (SVs) are less common than single nucleotide polymorphisms and indels in the population, but collectively account for a significant fraction of genetic polymorphism and diseases. Base pair differences arising from SVs are on a much higher order (>100 fold) than point mutations; however, none of the current detection methods are comprehensive, and currently available methodologies are incapable of providing sufficient resolution and unambiguous information across complex regions in the human genome. To address these challenges, we applied a high-throughput, cost-effective genome mapping technology to comprehensively discover genome-wide SVs and characterize complex regions of the YH genome using long single molecules (>150 kb) in a global fashion.

Results

Utilizing nanochannel-based genome mapping technology, we obtained 708 insertions/deletions and 17 inversions larger than 1 kb. Excluding the 59 SVs (54 insertions/deletions, 5 inversions) that overlap with N-base gaps in the reference assembly hg19, 666 non-gap SVs remained, and 396 of them (60%) were verified by paired-end data from whole-genome sequencing-based re-sequencing or de novo assembly sequence from fosmid data. Of the remaining 270 SVs, 260 are insertions and 213 overlap known SVs in the Database of Genomic Variants. Overall, 609 out of 666 (90%) variants were supported by experimental orthogonal methods or historical evidence in public databases. At the same time, genome mapping also provides valuable information for complex regions with haplotypes in a straightforward fashion. In addition, with long single-molecule labeling patterns, exogenous viral sequences were mapped on a whole-genome scale, and sample heterogeneity was analyzed at a new level.

Conclusion

Our study highlights genome mapping technology as a comprehensive and cost-effective method for detecting structural variation and studying complex regions in the human genome, as well as deciphering viral integration into the host genome.

Background

A structural variant (SV) is generally defined as a region of DNA 1 kb and larger in size that is different with respect to another DNA sample [1]; examples include inversions, translocations, deletions, duplications and insertions. Deletions and duplications are also referred to as copy number variants (CNVs). SVs have proven to be an important source of human genetic diversity and disease susceptibility [26]. Base pair differences arising from SVs occur on a significantly higher order (>100 fold) than point mutations [7,8], and data from the 1000 Genomes Project show population-specific patterns of SV prevalence [9,10]. Also, recent studies have firmly established that SVs are associated with a number of human diseases ranging from sporadic syndromes and Mendelian diseases to common complex traits, particularly neurodevelopmental disorders [1113]. Chromosomal aneuploidies, such as trisomy 21 and monosomy X have long been known to be the cause of Down's and Turner syndromes, respectively. A microdeletion at 15q11.2q12 has been shown causal for Prader-Willi syndrome [14], and many submicroscopic SV syndromes have been revealed since then [15]. In addition, rare, large de novo CNVs were identified to be enriched in autism spectrum disorder (ASD) cases [16], and other SVs were described as contributing factors for other complex traits including cancer, schizophrenia, epilepsy, Parkinson's disease and immune diseases, such as psoriasis (reviewed in [11] and [12]). With the increasing recognition of the important role of genomic aberrations in disease and the need for improved molecular diagnostics, comprehensive characterization of these genomic SVs is vital for, not only differentiating pathogenic events from benign ones, but also for rapid and full-scale clinical diagnosis.

While a variety of experimental and computational approaches exist for SV detection, each has its distinct biases and limitations. Hybridization-based approaches [1719] are subject to amplification, cloning and hybridization biases, incomplete coverage, and low dynamic range due to hybridization saturation. Moreover, detection of CNV events by these methods provides no positional context, which is critical to deciphering their functional significance. More recently, high-throughput next generation sequencing (NGS) technologies have been heavily applied to genome analysis based on alignment/mapping [2022] or de novo sequence assembly (SA) [23]. Mapping methods include paired-end mapping (PEM) [20], split-read mapping (SR) [21] and read depth analysis (RD) [22]. These techniques can be powerful, but are tedious and biased towards deletions owing to typical NGS short inserts and short reads [24,25]. De novo assembly methods are more versatile and can detect a larger range of SV types and sizes (0 ∼ 25 kb) by pair-wise genome comparison [2325]. All such NGS-based approaches lack power for comprehensiveness and are heavily biased against repeats and duplications because of short-read mapping ambiguity and assembly collapse [9,10,26]. David C. Schwartz's group promoted optical mapping [27] as an alternative to detect SVs along the genome with restriction mapping profiles of stretched DNA, highlighting the use of long single-molecule DNA maps in genome analysis. However, since the DNA is immobilized on glass surfaces and stretched, the technique suffers from low throughput and non-uniform DNA stretching, resulting in imprecise DNA length measurement and high error rate, hindering its utility and adoption [24,27–29]. Thus, an effective method to help detect comprehensive SVs and reveal complex genomic regions is needed.

The nanochannel-based genome mapping technology, commercialized as the “Irys” platform, automatically images fluorescently labeled DNA molecules in a massively parallel nanochannel array, and was introduced as an advanced technology [30] compared to other restriction mapping methods because of high-throughput data collection and its robust and highly uniform linearization of DNA in nanochannels. This technology has previously been described and used to map the 4.7-Mb highly variable human major histocompatibility complex (MHC) region [31], as well as for de novo assembly of a 2.1-Mb region in the highly complex Aegilops tauschii genome [32], lending great promise for use in complete genome sequence analysis. Here, we apply this rapid and high-throughput genome mapping method to discern genome wide SVs, as well as explore complex regions based on the YH (first Asian genome) [33] cell line. The workflow for mapping a human genome on Irys requires no library construction; instead, whole genomic DNA is labeled, stained and directly loaded into nanochannels for imaging. With the current throughput, one can collect enough data for de novo assembly of a human genome in less than three days. Additionally, comprehensive SV detection can be accomplished with genome mapping alone, without the addition of orthogonal technologies or multiple library preparations. Utilizing genome mapping, we identified 725 SVs including insertions/deletions, inversions, as well as SVs involved in N-base gap regions that are difficult to assess by current methods. For 50% of these SVs, we detected a signal of variation by re-sequencing and an additional 10% by fosmid sequence-based de novo assembly whereas the remainder had no signal by sequencing, hinting at the intractability of detection by sequencing. Detailed analyses showed most of the undetected SVs (80%, 213 out of 270) could be found overlapped in the Database of Genomic Variant (DGV) database indicating their reliability. Genome mapping also provides valuable haplotype information on complex regions, such as MHC, killer cell Immunoglobulin-like receptor (KIR), T cell receptor alpha/beta (TRA/TRB) and immunoglobulin light/heavy locus (IGH/IGL), which can help determine these hyper-variable regions’ sequences and downstream functional analyses. In addition, with long molecule labeling patterns, we were able to accurately map the exogenous virus’ sequence that integrated into the human genome, which is useful for the study of the mechanism of how virus sequence integration leads to serious diseases like cancer.

Data description

High-molecular weight DNA was extracted from the YH cell line, and high-quality DNA was labeled and run on the Irys system. After excluding DNA molecules smaller than 100 kb for analysis, we obtained 303 Gb of data giving 95× depth for the YH genome (Table 1). For subsequent analyses, only molecules larger than 150 kb (223 Gb, ∼70X) were used. De novo assembly resulted in a set of consensus maps with an N50 of 1.03 Mb. We performed “stitching” of neighboring genome maps that were fragmented by fragile sites associated with nick sites immediately adjacent to each other. After fragile site stitching, the N50 improved to 2.87 Mb, and the assembly covered 93.0% of the non-N base portion of the human genome reference assembly hg19. Structural variation was classified as a significant discrepancy between the consensus maps and the hg19 in silico map. Further analyses were performed for highly repetitive regions, complex regions and Epstein-Barr virus (EBV) integration. Supporting data is available from the GigaScience database, GigaDB [3436].

Analyses

Generation of single-molecule sequence motif maps

Genome maps were generated for the YH cell line by purifying high-molecular weight DNA in a gel plug and labeling at single-strand nicks created by the Nt.BspQI nicking endonuclease. Molecules were then linearized in nanochannel arrays etched in silicon wafers for imaging [31,32]. From these images, a set of label locations on each DNA molecule defined an individual single-molecule map. Single molecules have, on average, one label every 9 kb and were up to 1 Mb in length. A total of 932,855 molecules larger than 150 kb were collected for a total length of 223 Gb (∼70-fold average depth) (Table 1). Molecules can be aligned to a reference to estimate the error rates in the single molecules. Here, we estimated the missing label rate is 10%, and the extra label rate is 17%. Most of the error associated with these reference differences are averaged out in the consensus de novo assembly. Distinct genetic features intractable to sequencing technologies, such as long arrays of tandem repeats were observed in the raw single molecules (

: Figure S1).

Table 1

Molecule collection statistics under different length thresholds

Length cutoff (kb) No. Molecules Total length (Gb) Estimated depth (X)* 
100 1,568,969 303 95 
120 1,313,832 275 86 
150 932,855 223 70 
180 659,149 179 56 
200 523,146 153 48 
300 170,432 68 21 
500 22,944 14 
Length cutoff (kb) No. Molecules Total length (Gb) Estimated depth (X)* 
100 1,568,969 303 95 
120 1,313,832 275 86 
150 932,855 223 70 
180 659,149 179 56 
200 523,146 153 48 
300 170,432 68 21 
500 22,944 14 
*

Estimated depth based on 3.2 Gb genome size.

De novo assembly of genome maps from single-molecule data

Single molecules were assembled de novo into consensus genome maps using an implementation of the overlap-layout-consensus paradigm [37]. An overlap graph was constructed by an initial pairwise comparison of all molecules >150 kb, by pattern matching using commercial software from BioNano Genomics. Thresholds for the alignments were based on a p-value appropriate for the genome size (thresholds can be adjusted for different genome sizes and degrees of complexity) to prevent spurious edges. This graph was used to generate a draft consensus map set that was improved by alignment of single molecules and recalculation of the relative label positions. Next, the consensus maps were extended by aligning overhanging molecules to the consensus maps and calculating a consensus in the extended regions. Finally, the consensus maps were compared and merged where patterns matched (Figure 1). The result of this de novo assembly is a genome map set entirely independent of known reference or external data. In this case, YH was assembled with an N50 of 1.03 Mb in 3,565 maps and an N50 of 2.87 Mb in 1,634 maps after stitching fragile sites (

: Figure S2 and : Table S1). These genome maps define motif positions that occur on every 9 kb on average, and these label site positions have a resolution of 1.45 kb. The standard deviation for interval measurements between two labels varies with length. For example, for a 10 kb interval, the standard deviation (SD) is 502 bp, and for a 100 kb interval, it is 1.2 kb. Consensus genome maps were aligned to an in silico Nt.BspQI sequence motif map of hg19. Ninety-nine percent of the genome maps could align to hg19 and they overlap 93% of the non-gap portion of hg19.

Structural variation analysis

Using the genome map assembly as input, we performed structural variation detection (Figure 1), and the genome maps were compared to hg19. Strings of intervals between labels/nick motifs were compared and when they diverged, an outlier p-value was calculated and SVs were called at significant differences (See Methods for details), generating a list of 725 SVs including 59 that overlapped with N-base gaps in hg19 (

, Spreadsheet 3). Based on the standard deviation of interval measurements, 1.5 kb is the smallest insertion or deletion that can be confidently measured for an interval of about 10 kb if there is no pattern change. However, if label patterns deviate from the reference, SVs with a net size difference less than 1.5 kb can be detected. : Figure S1 shows three mapping examples (one deletion, one insertion, and one inversion) of gap region SVs. We present these 59 events separately although technically, in those cases, genome mapping detected structural differences between the genome maps and reference regions. For the remaining 666 SVs, 654 of them were insertions/deletions (Figure 2) while 12 were inversions (, Spreadsheet 1 & 2). Out of the 654 insertions/deletions, 503 were defined as insertions and 151 were deletions, demonstrating an enrichment of insertions for this individual with respect to the hg19 reference (Figure 2). Of the 59 SV events that span N-gap regions, 5 of them were inversions. Of the remaining 54 events, 51 were estimated to be shorter than indicated and 3 longer. These gap-region related SVs indicate a specific structure of gap regions of the YH genome compared to the hg19 reference.

Figure 1

Flowchart of consensus genome map assembly and structural variant discovery using genome mapping data.

Figure 1

Flowchart of consensus genome map assembly and structural variant discovery using genome mapping data.

In order to validate our SVs, we first cross-referenced them with the public SV database DGV (http://dgv.tcag.ca/dgv/app/home) [38]. For each query SV, we required 50% overlap with records in DGV. We found that the majority of the SVs (583 out of 666; 87.5%) could be found (

, Spreadsheet 1 & 2), confirming their reliability. Next, we applied the NGS discordant paired-end mapping and read depth-based methods, as well as fosmid-based de novo assembly (See Methods for detail), and as a result, detected an SV signal in 396 (60%, Figure 2) out of 666 SVs by at least one of the two methods (Figure 2, , Spreadsheet 1 & 2). For the remaining 270 SVs, 79% (213 out of 270, , Spreadsheet 1 & 2) were found in the DGV database. Overall, 91% (609 out of 666, , Spreadsheet 1 & 2) of SVs had supporting evidence by retrospectively applied sequencing-based methods or database entries.

Figure 2

Size distribution of total detected large insertions (green) and deletions (purple) using genome mapping. The comparative histogram bars in red and blue respectively represent deletions and insertions supported by NGS. NGS: next-generation sequencing.

Figure 2

Size distribution of total detected large insertions (green) and deletions (purple) using genome mapping. The comparative histogram bars in red and blue respectively represent deletions and insertions supported by NGS. NGS: next-generation sequencing.

We wanted to determine if SVs revealed by genome mapping, but without an NGS supported signal, had unique properties. We firstly investigated the distribution of NGS-supported SVs and NGS-unsupported SVs in repeat-rich and segmental duplication regions. However we did not find significant differences between them (data not shown) which was in concordance with previous findings [27]. We also compared the distribution of insertions and deletions of different SV categories and found that SV events that were not supported by sequencing evidence were 97% (260 out of 268) insertions; in contrast, the SVs that were supported by sequencing evidence were only 61% (243 out of 396, Figure 2,

, Spreadsheet 1) insertions showing insertion enrichment (p = 2.2e-16 Chi-squared test, Figure 2) in SVs without sequencing evidence. In addition, we further investigated the novel 57 SVs without either sequencing evidence or database supporting evidence. We found that the genes they covered had important functions, such as ion binding, enzyme activating and so forth, indicating their important role in cellular biochemical activities. Some of the genes like ELMO1, HECW1, SLC30A8, SLC16A12, JAM3 are reported to be associated with diseases like diabetic nephropathy, lateral sclerosis, diabetes mellitus and cataracts [39], providing valuable foundation for clinical application (, Spreadsheet 1 & 2).

Highly repetitive regions of the human genome

Highly repetitive regions of the human genome are known to be nearly intractable by NGS because short reads are often collapsed, and these regions are often refractory to cloning. We have searched for and analyzed one class of simple tandem repeats (unit size ranging from 2–13 kb) in long molecules derived from the genomes of YH (male) and CEPH-NA12878 (female). The frequencies of these repeating units from both genomes were plotted in comparison with hg19 (Figure 3). We found repeat units across the entire spectrum of sizes in YH and NA12878 while there were only sporadic peaks in hg19, implying an under representation of copy number variation as described in the current reference assembly. Furthermore, we have found a very large peak of approximately 2.5-kb repeats in YH (male, 691 copies) but not in NA19878 (female, 36 copies; Figure 3). This was further supported by additional genome mapping in other males and females demonstrating a consistent and significant quantity of male-specific repeats of 2.5 kb (unpublished). As an example,

: Figure S3 shows a raw image of an intact long molecule of 630 kb with two tracts of at least 53 copies and at least 21 copies of 2.5-kb tandem repeats (each 2.5-kb unit has one nick label site, creating the evenly spaced pattern) physically linked by another label-absent putative tandem repeat spanning over 435 kb, and : Figure S4 shows convincing mapping information. Unambiguously elucidating the absolute value and architecture of such complex repeat regions is not possible with other short fragment or hybridization-based methods.

Complex region analysis using genome mapping

Besides SV detection, genome mapping data also provide abundant information about other complex regions in the genome. For complex regions that are functionally important, an accurate reference map is critical for precise sequence assembly and integration for functional analysis [4043]. We analyzed the structure of some complex regions of the human genome. They include MHC also called Human leukocyte antigen (HLA), KIR, IGL/IGH, as well as TRA/TRB [4448]. In the highly variable HLA-A and −C loci, the YH genome shared one haplotype with the previously typed PGF genome (used in hg19) and also revealed an Asian/YH-specific variant on maps 209 and 153 (

: Figure S5), respectively. In the variant haplotype (Map ID 153), there is a large insertion at the HLA-A locus while at the HLA-D and RCCX loci, YH had an Asian/YH-specific insertion and a deletion. In addition to the MHC region, we also detected Asian/YH-specific structural differences in KIR (: Figure S6), IGH/IGL (: Figure S7), and TRA/TRB (: Figure S8), compared to the reference genome.

External sequence integration detection using genome mapping

External viral sequence integration detection is important for the study of diseases such as cancer, but current high-throughput methods are limited in discovering integration break points [4951]. Although fiber fluorescence in situ hybridization (FISH) was used to discriminate between integration and episomal forms of virus utilizing long dynamic DNA molecules [52], this method was laborious, low-resolution and low-throughput. Thus, long, intact high-resolution single-molecule data provided by genome mapping allows for rapid and effective analysis of which part of the virus sequence has been integrated into the host genome and its localization. We detected EBV integration into the genome of the cell line sample.

Figure 3

A plot of repeat units in two human genomes as seen in single molecules. A repeat unit is defined as five or more equidistant labels. Total units in bins are normalized to the average coverage depth in the genome.

Figure 3

A plot of repeat units in two human genomes as seen in single molecules. A repeat unit is defined as five or more equidistant labels. Total units in bins are normalized to the average coverage depth in the genome.

The EBV virus map was assembled de novo during the whole genome de novo assembly of the YH cell line genome. We mapped the de novo EBV map to in silico maps from public databases to determine the strain that was represented in the cell line. We found that the YH strain was most closely related, although not identical, to strain B95-8 (GenBank: V01555.2). To detect EBV integration, portions of the aligned molecules extending beyond the EBV map were extracted and aligned with hg19 to determine potential integration sites (

: Figure S9). There are 1,340 EBV integration events across the genome (Figure 4). We found that the frequency of EBV integration mapping was significantly lower than the average coverage depth (∼70X), implying the DNA sample derived from a clonal cell population is potentially more diverse than previously thought, and that this method could reveal the heterogeneity of a very complex sample population at the single-molecule level. Also, the integrated portion of the EBV genome sequence was detected with a larger fraction towards the tail (: Figure S10). Besides integration events, we also found EBV episome molecules whose single-molecule map could be mapped to the EBV genome, free of flanking human genomic regions.

Discussion

Structural variants are increasingly frequently shown to play important roles in human health. However, available technologies, such as array-CGH, SNP array and NGS are incapable of cataloguing them in a comprehensive and unbiased manner. Genome mapping, a technology successfully applied to the assembly of complex regions of a plant genome and characterization of structural variation and haplotype differences in the human MHC region, has been adopted to capture the genome-wide structure of a human individual in the current study. Evidence for over 600 SVs in this individual has been provided. Despite the difficulty of SV detection by sequencing methods, the majority of genome map-detected SVs were retrospectively found to have signals consistent with the presence of an SV, validating genome mapping for SV discovery. Approximately 75% of the SVs discovered by genome mapping were insertions; this interesting phenomenon may be a method bias or a genuine representation of the additional content in this genome of Asian descent that is not present in hg19, which was compiled based on genomic materials presumably derived from mostly non-Asians. Analysis of additional genomes is necessary for comparison. Insertion detection is refractory to many existing methodologies [24,25], so to some extent, genome mapping revealed its distinct potential to address this challenge. Furthermore, functional annotation results of the detected SVs show that 30% of them (

, Spreadsheet 1 & 2) affect exonic regions of relevant genes which may cause severe effects on gene function. Gene ontology (GO) analysis demonstrates that these SVs are associated with genes that contribute to important biological processes (, Spreadsheet 1 & 2 and : Figure S11), reflecting that the SVs detected here are likely to affect a large number of genes and may have a significant impact on human health. Genome mapping provides us with an effective way to study the impact of genome-wide SV on human conditions. Some N-base gaps are estimated to have longer or shorter length or more complex structurally compared to hg19, demonstrating that genome mapping is useful for improving the human and other large genome assemblies. We also present a genome-wide analysis of short tandem repeats in individual human genomes and structural information and differences for some of the most complex regions in the YH genome. Independent computational analysis has been performed to discern exogenous viral insertions, as well as exogenous episomes. All of these provide invaluable insights into the capacity of genome mapping as a promising new strategy for research and clinical application.

Figure 4

Circos plot of distribution of integration events throughout YH genome. The genome was divided into non-overlapping windows of 200 kb. The number of molecules with evidence of integration in each window is plotted with each concentric grey circle representing a two-fold increment in virus detection.

Figure 4

Circos plot of distribution of integration events throughout YH genome. The genome was divided into non-overlapping windows of 200 kb. The number of molecules with evidence of integration in each window is plotted with each concentric grey circle representing a two-fold increment in virus detection.

The basis for the genome mapping technology that enables us to effectively address shortcomings of existing methodologies is the use of motif maps derived from extremely long DNA molecules hundreds of kb in length. Using these motif maps, we are able to also access challenging loci where existing technologies fail. Firstly, global structural variations were easily and quickly detected. Secondly, evidence for a deletion bias which is commonly observed with both arrays and NGS technology, is absent in genome mapping. In fact, we observe more insertions than deletions in this study. Thirdly, for the first time, we are able to measure the length of regions of the YH genome that represent gaps in the human reference assembly. Fourthly, consensus maps could be assembled in highly variable regions in the YH genome which are important for subsequent functional analysis. Finally, both integrated and un-integrated EBV molecules are identified, and potential sub-strains differentiated, and the EBV genome sequence that integrated into the host genome was obtained directly. This information was previously inaccessible without additional PCR steps or NGS approaches [50]. All in all, we demonstrated advantages and strong potential of the genome mapping technology based on nanochannel arrays to help overcome problems that have severely limited our understanding of the human genome.

In addition to the advantages this study reveals about the genome mapping technology, aspects that need to be improved also are highlighted. As genome mapping technology generates sequence-specific motif-labeled DNA molecules and analyzes these motif maps using an overlap-layout-consensus algorithm, subsequent performance and resolution largely depends on motif density (any individual event endpoints can only be resolved to the nearest restriction sites). For example, the EBV integration analysis in this study was more powerful in the high-density regions (

: Figure S10). Hence, higher density labeling methods to increase the information density that may promote even higher accuracy and unbiased analysis of genomes are currently being further developed. When data from genome mapping is combined with another source of information, one can achieve even higher resolution for each event. In addition, reducing random errors like extra restriction sites, missing restriction sites and size measurement is important for subsequent analysis. Finally, improvements to the SV detection algorithm will provide further discovery potential, and balanced reciprocal translocations can be identified in genome maps generated from cancer model genomes (personal communication, Michael Rossi).

The throughput and speed of a technology remains one of the most important factors for routine use in clinical screening as well as scientific research. At the time of manuscript submission, genome mapping of a human individual could be accomplished with fewer than three nanochannel array chips in a few days. It is anticipated that a single nanochannel chip would cover a human size genome in less than one day within 6 months, facilitating new studies aimed at unlocking the inaccessible portions of the genome. In this way, genome mapping has an advantage over the use of multiple orthogonal methods that are often used to detect global SVs. Thus, it is now feasible to conduct large population-based comprehensive SV studies efficiently on a single platform.

Methods

High-molecular weight DNA extraction

High-molecular weight (HMW) DNA extraction was performed as recommended for the CHEF Mammalian Genomic DNA Plug Kit (BioRad #170-3591). Briefly, cells from the YH or NA12878 cell lines were washed with 2x with PBS and resuspended in cell resuspension buffer, after which 7.5 × 105 cells were embedded in each gel plug. Plugs were incubated with lysis buffer and proteinase K for four hours at 50°C. The plugs were washed and then solubilized with GELase (Epicentre). The purified DNA was subjected to four hours of drop dialysis (Millipore, #VCWP04700) and quantified using Nanodrop 1000 (Thermal Fisher Scientific) and/or the Quant-iT dsDNA Assay Kit (Invitrogen/Molecular Probes).

DNA labeling

DNA was labeled according to commercial protocols using the IrysPrep Reagent Kit (BioNano Genomics, Inc). Specifically, 300 ng of purified genomic DNA was nicked with 7 U nicking endonuclease Nt.BspQI (New England BioLabs, NEB) at 37°C for two hours in NEB Buffer 3. The nicked DNA was labeled with a fluorescent-dUTP nucleotide analog using Taq polymerase (NEB) for one hour at 72°C. After labeling, the nicks were ligated with Taq ligase (NEB) in the presence of dNTPs. The backbone of fluorescently labeled DNA was stained with YOYO-1 (Invitrogen).

Data collection

The DNA was loaded onto the nanochannel array of BioNano Genomics IrysChip by electrophoresis of DNA. Linearized DNA molecules were then imaged automatically followed by repeated cycles of DNA loading using the BioNano Genomics Irys system.

The DNA molecules backbones (YOYO-1 stained) and locations of fluorescent labels along each molecule were detected using the in-house software package, IrysView. The set of label locations of each DNA molecule defines an individual single-molecule map.

De novo genome map assembly

Single-molecule maps were assembled de novo into consensus maps using software tools developed at BioNano Genomics. Briefly, the assembler is a custom implementation of the overlap-layout-consensus paradigm with a maximum likelihood model. An overlap graph was generated based on pairwise comparison of all molecules as input. Redundant and spurious edges were removed. The assembler outputs the longest path in the graph and consensus maps were derived. Consensus maps are further refined by mapping single-molecule maps to the consensus maps and label positions are recalculated. Refined consensus maps are extended by mapping single molecules to the ends of the consensus and calculating label positions beyond the initial maps. After merging of overlapping maps, a final set of consensus maps was generated and used for subsequent analysis. Furthermore, we applied a “stitching” procedure to join neighboring genome maps. Two adjacent genome maps would be joined together if the junction a) was within 50 kb apart, b) contained at most 5 labels, c) contained, or was within 50 kb from, a fragile site, and d) also contained no more than 5 unaligned end labels. If these criteria were satisfied, the two genome maps would be joined together with the intervening label patterns taken from the reference in silico map.

Structural variation detection

Alignments between consensus genome maps and the hg19 in silico sequence motif map were obtained using a dynamic programming approach where the scoring function was the likelihood of a pair of intervals being similar [53]. Likelihood is calculated based on a noise model which takes into account fixed sizing error, sizing error which scales linearly with the interval size, misaligned sites (false positives and false negatives), and optical resolution. Within an alignment, an interval or range of intervals whose cumulative likelihood for matching the reference map is worse than 0.01 percent chance is classified as an outlier region. If such a region occurs between highly scoring regions (p-value of 10e−6), an insertion or deletion call is made in the outlier region, depending on the relative size of the region on the query and reference maps. Inversions are defined if adjacent match-groups between the genome map and reference are in reverse relative orientation.

Signals refined by re-sequencing and de novo assembly based methods

In order to demonstrate the capacity of genome mapping for the detection of large SVs, we tested the candidate SVs using whole-genome paired-end 100 bp sequencing (WGS) data with insert sizes of 500 bp and fosmid sequence based de novo assembly result. SVs were tested based on the expectation that authentic SVs would be supported by abnormally mapped read pairs, and that deletions with respect to the reference should have lower mapped read depth than average [20,22,23]. We performed single-end/(paired-end + single-end) reads ratio (sp ratio) calculations at the whole-genome level to assign an appropriate threshold for abnormal regions as well as depth coverage. We set sp ratio and depth cutoff thresholds based on the whole genome data to define SV signals. Insertions with aberrant sp ratio and deletions with either sp ratio or abnormal depth were defined to be a supported candidate.

We also utilized fosmid-based de novo assembly data to search for signals supporting candidate SVs. We used contigs and scaffolds assembled from short reads to check for linearity between a given assembly and hg19 using LASTZ [54]. WGS-based and fosmid-based SV validation showed inconsistency and/or lack of saturation as each supported unique variants (

: Figure S2) [24].

EBV integration detection

Single-molecule maps were aligned with a map generated in silico based on the EBV reference sequence (strain B95-8; GenBank: V01555.2). Portions of the aligned molecules extending beyond the EBV map were extracted and aligned with hg19 to determine potential integration sites.

Availability of supporting data

The data sets supporting the results of this article is available in the GigaScience GigaDB, repository [55]. See the individual GigaDB entries for the YH Bionano data [35] and YH fosmid validation data [36], which is also available in the SRA [PRJEB7886].

Abbreviations

  • Array-CGH

    Array-based comparative genomic hybridization

  • AS

    De novo sequence assembly

  • ASD

    Autism spectrum disorder

  • BCR

    B cell receptor

  • CNV

    Copy number variant

  • DGV

    Database of genomic variants

  • EBV

    Epstein-Barr virus

  • FISH

    Fluorescence in situ hybridization

  • GO

    Gene ontology

  • HLA

    Human leukocyte antigen

  • HMW

    High-molecular weight

  • IGH

    Immunoglobulin heavy locus

  • IGL

    Immunoglobulin light locus

  • KIR

    Killer cell immunoglobulin-like receptor

  • LRC

    Leukocyte Receptor Complex

  • MHC

    Major histocompatibility complex

  • NGS

    Next-generation sequencing

  • PCR

    Polymerase chain reaction

  • PEM

    Pair-end mapping

  • RD

    Read depth

  • SNP

    Single nucleotide polymorphism

  • SR

    Split read

  • SV

    Structural variation

  • TCR

    T cell receptor

  • TRA

    T cell receptor alpha locus

  • TRB

    T cell receptor beta locus

  • WGS

    Whole-genome sequencing

  • YH

    YanHuang

Competing interests

The authors have the following interests: Alex R. Hastie, Ernest T. Lam, Warren Andrews, Saki Chan, Michael Requa, Thomas Anantharaman and Han Cao were employees of BioNano Genomics at the time of the study, and they owned company stock options.

Authors’ contributions

XX, HC and HZC managed the project. HC, HZC, ARH, DC, ETL and SC designed the analyses. HC, ARH, ETL, W A, MR and TA produced genome mapping data. HZC, ARH, DC, ETL, YS, HH, XL, LY, WA, SC, SH, XT, MR, TA, AK and HY performed the data analyses. HZC, ARH and DC did most of the writing with contributions from all authors. All authors read and approved the final manuscript.

Acknowledgements

We wish to recognize BGI-Shenzhen's sequencing platform for generating the data in this study. We thank the faculty and staff at BGI-Shenzhen who contributed to this project. We also thank everyone at BioNano Genomics who contributed to the development of the lrys system and the team's helpful discussion and analysis. This work was supported by the State Key Development Program for Basic Research of China—973 Program (NO.2011CB809202); the National High Technology Research and Development Program of China — 863 Program (NO.2012AA02A201); the Shenzhen Municipal Government of China (NOJC201005260191A); Shenzhen Key Laboratory of Transomics Biotechnologies (NO.CXB201108250096A).

References

1.
Freeman
JL
Perry
GH
Feuk
L
Redon
R
McCarroll
SA
Altshuler
DM
Aburatani
H
Jones
KW
Tyler-Smith
C
Hurles
ME
Carter
NP
Scherer
SW
Lee
C
2006
Copy number variation: new insights in genome diversity
Genome Res
 , vol. 
16
 (pg. 
949
-
961
10.1101/gr.3677206
2.
Tuzun
E
Sharp
AJ
Bailey
JA
Kaul
R
Morrison
VA
Pertz
LM
Haugen
E
Hayden
H
Albertson
D
Pinkel
D
Olson
MV
Eichler
EE
2005
Fine-scale structural variation of the human genome
Nat Genet
 , vol. 
37
 (pg. 
727
-
732
10.1038/ng1562
3.
Feuk
L
Carson
AR
Scherer
SW
2006
Structural variation in the human genome
Nat Rev Genet
 , vol. 
7
 (pg. 
85
-
97
)
4.
Kidd
JM
Cooper
GM
Donahue
WF
Hayden
HS
Sampas
N
Graves
T
Hansen
N
Teague
B
Alkan
C
Antonacci
F
Haugen
E
Zerr
T
Yamada
NA
Tsang
P
Newman
TL
Tüzün
E
Cheng
Z
Ebling
HM
Tusneem
N
David
R
Gillett
W
Phelps
KA
Weaver
M
Saranga
D
Brand
A
Tao
W
Gustafson
E
McKernan
K
Chen
L
Malig
M
, et al.  . 
2008
Mapping and sequencing of structural variation from eight human genomes
Nature
 , vol. 
453
 (pg. 
56
-
64
10.1038/nature06862
5.
Sebat
J
Lakshmi
B
Troge
J
Alexander
J
Young
J
Lundin
P
Månér
S
Massa
H
Walker
M
Chi
M
Navin
N
Lucito
R
Healy
J
Hicks
J
Ye
K
Reiner
A
Gilliam
TC
Trask
B
Patterson
N
Zetterberg
A
Wigler
M
2004
Large-scale copy number polymorphism in the human genome
Science
 , vol. 
305
 (pg. 
525
-
528
10.1126/science.1098918
6.
Itsara
A
Cooper
GM
Baker
C
Girirajan
S
Li
J
Absher
D
Krauss
RM
Myers
RM
Ridker
PM
Chasman
DI
Mefford
H
Ying
P
Nickerson
DA
Eichler
EE
2009
Population analysis of large copy number variants and hotspots of human genetic disease
Am J Hum Genet
 , vol. 
84
 (pg. 
148
-
161
10.1016/j.ajhg.2008.12.014
7.
Cheng
Z
Ventura
M
She
X
Khaitovich
P
Graves
T
Osoegawa
K
Church
D
DeJong
P
Wilson
RK
Pääbo
S
Rocchi
M
Eichler
EE
2005
A genome-wide comparison of recent chimpanzee and human segmental duplications
Nature
 , vol. 
437
 (pg. 
88
-
93
10.1038/nature04000
8.
Lupski
JR
2007
Genomic rearrangements and sporadic disease
Nat Genet
 , vol. 
39
 (pg. 
S43
-
S47
10.1038/ng2084
9.
Abecasis
GR
Altshuler
D
Auton
A
Brooks
LD
Durbin
RM
Gibbs
RA
Hurles
ME
McVean
GA
2010
A map of human genome variation from population-scale sequencing
Nature
 , vol. 
467
 (pg. 
1061
-
1073
10.1038/nature09534
10.
Mills
RE
Walter
K
Stewart
C
Handsaker
RE
Chen
K
Alkan
C
Abyzov
A
Yoon
SC
Ye
K
Cheetham
RK
Chinwalla
A
Conrad
DF
Fu
Y
Grubert
F
Hajirasouliha
I
Hormozdiari
F
Iakoucheva
LM
Iqbal
Z
Kang
S
Kidd
JM
Konkel
MK
Korn
J
Khurana
E
Kural
D
Lam
HY
Leng
J
Li
R
Li
Y
Lin
CY
Luo
R
, et al.  . 
2011
Mapping copy number variation by population-scale genome sequencing
Nature
 , vol. 
470
 (pg. 
59
-
65
10.1038/nature09708
11.
Stankiewicz
P
Lupski
JR
2010
Structural variation in the human genome and its role in disease
Annu Rev Med
 , vol. 
61
 (pg. 
437
-
455
10.1146/annurev-med-100708-204735
12.
Girirajan
S
Campbell
CD
Eichler
EE
2011
Human copy number variation and complex genetic disease
Annu Rev Genet
 , vol. 
45
 (pg. 
203
-
226
10.1146/annurev-genet-102209-163544
13.
Weischenfeldt
J
Symmons
O
Spitz
F
Korbel
JO
2013
Phenotypic impact of genomic structural variation: insights from and for human disease
Nat Rev Genet
 , vol. 
14
 (pg. 
125
-
138
10.1038/nrg3373
14.
Ledbetter
DHRV
Airhart
SD
Strobel
RJ
Keenan
BS
Crawford
JD
1981
Deletions of chromosome 15 as a cause of the Prader-Willi syndrome
N Engl J Med
 , vol. 
304
 (pg. 
325
-
329
10.1056/NEJM198102053040604
15.
Vissers
LE
Stankiewicz
P
2012
Microdeletion and microduplication syndromes
Methods Mol Biol
 , vol. 
838
 (pg. 
29
-
75
10.1007/978-1-61779-507-7_2
16.
Sebat
J
Lakshmi
B
Malhotra
D
Troge
J
Lese-Martin
C
Walsh
T
Yamrom
B
Yoon
S
Krasnitz
A
Kendall
J
Leotta
A
Pai
D
Zhang
R
Lee
YH
Hicks
J
Spence
SJ
Lee
AT
Puura
K
Lehtimäki
T
Ledbetter
D
Gregersen
PK
Bregman
J
Sutcliffe
JS
Jobanputra
V
Chung
W
Warburton
D
King
MC
Skuse
D
Geschwind
DH
Gilliam
TC
, et al.  . 
2007
Strong association of de novo copy number mutations with autism
Science
 , vol. 
316
 (pg. 
445
-
449
10.1126/science.1138659
17.
Redon
R
Ishikawa
S
Fitch
KR
Feuk
L
Perry
GH
Andrews
TD
Fiegler
H
Shapero
MH
Carson
AR
Chen
W
Cho
EK
Dallaire
S
Freeman
JL
González
JR
Gratacòs
M
Huang
J
Kalaitzopoulos
D
Komura
D
MacDonald
JR
Marshall
CR
Mei
R
Montgomery
L
Nishimura
K
Okamura
K
Shen
F
Somerville
MJ
Tchinda
J
Valsesia
A
Woodwark
C
Yang
F
, et al.  . 
2006
Global variation in copy number in the human genome
Nature
 , vol. 
444
 (pg. 
444
-
454
10.1038/nature05329
18.
McCarroll
SA
Kuruvilla
FG
Korn
JM
Cawley
S
Nemesh
J
Wysoker
A
Shapero
MH
de Bakker
PI
Maller
JB
Kirby
A
Elliott
AL
Parkin
M
Hubbell
E
Webster
T
Mei
R
Veitch
J
Collins
PJ
Handsaker
R
Lincoln
S
Nizzari
M
Blume
J
Jones
KW
Rava
R
Daly
MJ
Gabriel
SB
Altshuler
D
2008
Integrated detection and population-genetic analysis of SNPs and copy number variation
Nat Genet
 , vol. 
40
 (pg. 
1166
-
1174
10.1038/ng.238
19.
Sharp
AJ
Locke
DP
McGrath
SD
Cheng
Z
Bailey
JA
Vallente
RU
Pertz
LM
Clark
RA
Schwartz
S
Segraves
R
Oseroff
VV
Albertson
DG
Pinkel
D
Eichler
EE
2005
Segmental duplications and copy-number variation in the human genome
Am J Hum Genet
 , vol. 
77
 (pg. 
78
-
88
10.1086/431652
20.
Chen
K
Wallis
JW
McLellan
MD
Larson
DE
Kalicki
JM
Pohl
CS
McGrath
SD
Wendl
MC
Zhang
Q
Locke
DP
Shi
X
Fulton
RS
Ley
TJ
Wilson
RK
Ding
L
Mardis
ER
2009
BreakDancer: an algorithm for high-resolution mapping of genomic structural variation
Nat Methods
 , vol. 
6
 (pg. 
677
-
681
10.1038/nmeth.1363
21.
Wang
J
Mullighan
CG
Easton
J
Roberts
S
Heatley
SL
Ma
J
Rusch
MC
Chen
K
Harris
CC
Ding
L
Holmfeldt
L
Payne-Turner
D
Fan
X
Wei
L
Zhao
D
Obenauer
JC
Naeve
C
Mardis
ER
Wilson
RK
Downing
JR
Zhang
J
2011
CREST maps somatic structural variation in cancer genomes with base-pair resolution
Nat Methods
 , vol. 
8
 (pg. 
652
-
654
10.1038/nmeth.1628
22.
Abyzov
A
Urban
AE
Snyder
M
Gerstein
M
2011
CNVnator: an approach to discover, genotype, and characterize typical and atypical CNVs from family and population genome sequencing
Genome Res
 , vol. 
21
 (pg. 
974
-
984
10.1101/gr.114876.110
23.
Li
Y
Zheng
H
Luo
R
Wu
H
Zhu
H
Li
R
Cao
H
Wu
B
Huang
S
Shao
H
Ma
H
Zhang
F
Feng
S
Zhang
W
Du
H
Tian
G
Li
J
Zhang
X
Li
S
Bolund
L
Kristiansen
K
de Smith
AJ
Blakemore
AI
Coin
LJ
Yang
H
Wang
J
Wang
J
2011
Structural variation in two human genomes mapped at single-nucleotide resolution by whole genome de novo assembly
Nature biotechnology
 , vol. 
29
 (pg. 
723
-
730
10.1038/nbt.1904
24.
Alkan
C
Coe
BP
Eichler
EE
2011
Genome structural variation discovery and genotyping
Nat Rev Genet
 , vol. 
12
 (pg. 
363
-
376
10.1038/nrg2958
25.
Zhao
M
Wang
Q
Jia
P
Zhao
Z
2013
Computational tools for copy number variation (CNV) detection using next-generation sequencing data: features and perspectives
BMC Bioinformatics
 , vol. 
14
 
11
pg. 
S1
 
26.
Alkan
C
Sajjadian
S
Eichler
EE
2011
Limitations of next-generation genome sequence assembly
Nat Methods
 , vol. 
8
 (pg. 
61
-
65
10.1038/nmeth.1527
27.
Teague
B
Waterman
MS
Goldstein
S
Potamousis
K
Zhou
S
Reslewic
S
Sarkar
D
Valouev
A
Churas
C
Kidd
JM
Kohn
S
Runnheim
R
Lamers
C
Forrest
D
Newton
MA
Eichler
EE
Kent-First
M
Surti
U
Livny
M
Schwartz
DC
2010
High-resolution human genome structure by single-molecule analysis
Proc Natl Acad Sci U S A
 , vol. 
107
 (pg. 
10848
-
10853
10.1073/pnas.0914638107
28.
Dong
Y
Xie
M
Jiang
Y
Xiao
N
Du
X
Zhang
W
Tosser-Klopp
G
Wang
J
Yang
S
Liang
J
Chen
W
Chen
J
Zeng
P
Hou
Y
Bian
C
Pan
S
Li
Y
Liu
X
Wang
W
Servin
B
Sayre
B
Zhu
B
Sweeney
D
Moore
R
Nie
W
Shen
Y
Zhao
R
Zhang
G
Li
J
Faraut
T
, et al.  . 
2013
Sequencing and automated whole-genome optical mapping of the genome of a domestic goat (Capra hircus)
Nat Biotechnol
 , vol. 
31
 (pg. 
135
-
141
)
29.
Levy-Sakin
M
Ebenstein
Y
2013
Beyond sequencing: optical mapping of DNA in the age of nanotechnology and nanoscopy
Curr Opin Biotechnol
 , vol. 
24
 (pg. 
690
-
698
10.1016/j.copbio.2013.01.009
30.
Das
SK
Austin
MD
Akana
MC
Deshpande
P
Cao
H
Xiao
M
2010
Single molecule linear analysis of DNA in nano-channel labeled with sequence specific fluorescent probes
Nucleic Acids Res
 , vol. 
38
 pg. 
e177
  
10.1093/nar/gkq673
31.
Lam
ET
Hastie
A
Lin
C
Ehrlich
D
Das
SK
Austin
MD
Deshpande
P
Cao
H
Nagarajan
N
Xiao
M
Kwok
PY
2012
Genome mapping on nanochannel arrays for structural variation analysis and sequence assembly
Nat Biotechnol
 , vol. 
30
 (pg. 
771
-
776
10.1038/nbt.2303
32.
Hastie
AR
Dong
L
Smith
A
Finklestein
J
Lam
ET
Huo
N
Cao
H
Kwok
PY
Deal
KR
Dvorak
J
Luo
MC
Gu
Y
Xiao
M
2013
Rapid genome mapping in nanochannel arrays for highly complete and accurate de novo sequence assembly of the complex Aegilops tauschii genome
PLoS One
 , vol. 
8
 pg. 
e55864
  
10.1371/journal.pone.0055864
33.
Wang
J
Wang
W
Li
R
Li
Y
Tian
G
Goodman
L
Fan
W
Zhang
J
Li
J
Zhang
J
Guo
Y
Feng
B
Li
H
Lu
Y
Fang
X
Liang
H
Du
Z
Li
D
Zhao
Y
Hu
Y
Yang
Z
Zheng
H
Hellmann
I
Inouye
M
Pool
J
Yi
X
Zhao
J
Duan
J
Zhou
Y
Qin
J
, et al.  . 
2008
The diploid genome sequence of an Asian individual
Nature
 , vol. 
456
 (pg. 
60
-
65
10.1038/nature07484
34.
Sneddon
TP
Li
P
Edmunds
SC
2012
GigaDB: announcing the GigaScience database
biotechnology
 , vol. 
1
 pg. 
11
 
35.
Cao
HZ
Hastie
AR
Cao
D
Lam
ET
Sun
Y
Huang
H
Liu
X
Lin
L
Andrews
W
Chan
S
Huang
S
Tong
X
Requa
M
Anantharaman
T
Krogh
A
Yang
H
Cao
H
Xu
X
2014
Supporting material for: "Rapid detection of structural variation in a Human genome using nanochannel-based genome mapping technology"
GigaScience Database
  
36.
Cao
HZ
Wu
H
Luo
R
Huang
S
Sun
Y
Tong
X
Xie
Y
Liu
B
Yang
H
Zheng
H
Li
J
Li
B
Wang
Y
Yang
F
Sun
P
Liu
S
Gao
P
Huang
H
Sun
J
Chan
D
Ha
G
Huang
W
Huang
Z
Li
Y
Tellier
LCAM
Liu
X
Feng
Q
Xu
X
Zhang
X
, et al.  . 
2014
Supporting material for: De novo assembly of a haplotype-resolved human genome
GigaScience Database
  
37.
Valouev
A
Schwartz
DC
Zhou
S
Waterman
MS
2006
An algorithm for assembly of ordered restriction maps from single DNA molecules
Proc Nat Acad Sci USA
 , vol. 
103
 (pg. 
15770
-
15775
10.1073/pnas.0604040103
38.
Macdonald
JR
Ziman
R
Yuen
RK
Feuk
L
Scherer
SW
2014
The Database of Genomic Variants: a curated collection of structural variation in the human genome
Nucleic Acids Res
 , vol. 
42
 (pg. 
D986
-
D992
10.1093/nar/gkt958
39.
McKusick
VA
2007
Mendelian Inheritance in Man and its online version, OMIM
Am J Hum Genet
 , vol. 
80
 (pg. 
588
-
604
10.1086/514346
40.
Horton
R
Gibson
R
Coggill
P
Miretti
M
Allcock
RJ
Almeida
J
Forbes
S
Gilbert
JG
Halls
K
Harrow
JL
Hart
E
Howe
K
Jackson
DK
Palmer
S
Roberts
AN
Sims
S
Stewart
CA
Traherne
JA
Trevanion
S
Wilming
L
Rogers
J
de Jong
PJ
Elliott
JF
Sawcer
S
Todd
JA
Trowsdale
J
Beck
S
2008
Variation analysis and gene annotation of eight MHC haplotypes: the MHC Haplotype Project
Immunogenetics
 , vol. 
60
 (pg. 
1
-
18
10.1007/s00251-007-0262-2
41.
Tesch
HMP
Krueger
GR
Fischer
R
Diehl
V
1990
Analysis of immunoglobulin, T cell receptor and bcr rearrangements in human malignant lymphoma and Hodgkin’s disease
Oncology
 , vol. 
47
 (pg. 
215
-
223
10.1159/000226819
42.
Rajagopalan
S
Long
EO
2005
Understanding how combinations of HLA and KIR genes influence disease
J Exp Med
 , vol. 
201
 (pg. 
1025
-
1029
10.1084/jem.20050499
43.
Gonzalez
D
van der Burg
M
Garcia-Sanz
R
Fenton
JA
Langerak
AW
Gonzalez
M
van Dongen
JJ
San Miguel
JF
Morgan
GJ
2007
Immunoglobulin gene rearrangements and the pathogenesis of multiple myeloma
Blood
 , vol. 
110
 (pg. 
3112
-
3121
10.1182/blood-2007-02-069625
44.
Beck
S
Geraghty
D
Inoko
H
Rowen
L
Aguado
B
Bahram
S
Campbell
RD
Forbes
SA
Guillaudeux
T
Hood
L
1999
Complete sequence and gene map of a human major histocompatibility complex
Nature
 , vol. 
401
 (pg. 
921
-
923
10.1038/44853
45.
Vilches
C
Parham
P
2002
KIR: diverse, rapidly evolving receptors of innate and adaptive immunity
Annu Rev Immunol
 , vol. 
20
 (pg. 
217
-
251
10.1146/annurev.immunol.20.092501.134942
46.
Katzmann
JA
Clark
RJ
Abraham
RS
Bryant
S
Lymp
JF
Bradwell
AR
Kyle
RA
2002
Serum reference intervals and diagnostic ranges for free kappa and free lambda immunoglobulin light chains: relative sensitivity for detection of monoclonal light chains
Clin Chem
 , vol. 
48
 (pg. 
1437
-
1444
)
47.
Tomlinson
IM
Cook
GP
Walter
G
Carter
NP
Riethman
H
Buluwela
L
Rabbitts
TH
Winter
G
1995
A complete map of the human immunoglobulin VH locus
Ann N Y Acad Sci
 , vol. 
764
 (pg. 
43
-
46
)
48.
Haynes
MR
Wu
GE
2007
Gene discovery at the human T-cell receptor alpha/delta locus
Immunogenetics
 , vol. 
59
 (pg. 
109
-
121
10.1007/s00251-006-0165-7
49.
Li
WZX
Lee
NP
Liu
X
Chen
S
Guo
B
Yi
S
Zhuang
X
Chen
F
Wang
G
Poon
RT
Fan
ST
Mao
M
Li
Y
Li
S
Wang
J
Jianwang
Xu
X
Jiang
H
Zhang
X
2013
HIVID: an efficient method to detect HBV integration using low coverage sequencing
Genomics
 , vol. 
102
 (pg. 
338
-
344
10.1016/j.ygeno.2013.07.002
50.
Kan
Z
Zheng
H
Liu
X
Li
S
Barber
TD
Gong
Z
Gao
H
Hao
K
Willard
MD
Xu
J
Hauptschein
R
Rejto
PA
Fernandez
J
Wang
G
Zhang
Q
Wang
B
Chen
R
Wang
J
Lee
NP
Zhou
W
Lin
Z
Peng
Z
Yi
K
Chen
S
Li
L
Fan
X
Yang
J
Ye
R
Ju
J
Wang
K
, et al.  . 
2013
Whole-genome sequencing identifies recurrent mutations in hepatocellular carcinoma
Genome Res
 , vol. 
23
 (pg. 
1422
-
1433
10.1101/gr.154492.113
51.
Zhao
G
Krishnamurthy
S
Cai
Z
Popov
VL
da Rosa Travassos
AP
Guzman
H
Cao
S
Virgin
HW
Tesh
RB
Wang
D
2013
Identification of Novel Viruses Using VirusHunter – an Automated Data Analysis Pipeline
PLoS One
 , vol. 
8
 pg. 
e78470
  
10.1371/journal.pone.0078470
52.
Reisinger
J
Rumpler
S
Lion
T
Ambros
PF
2006
Visualization of episomal and integrated Epstein-Barr virus DNA by fiber fluorescence in situ hybridization
Int J Cancer
 , vol. 
118
 (pg. 
1603
-
1608
10.1002/ijc.21498
53.
Anantharaman
T
Mishra
B
2001
A probabilistic analysis of false positives in optical map alignment and validation
Proc. of WABI
 (pg. 
27
-
40
)
54.
Harris
RS
2007
Improved pairwise alignment of genomic DNA
 
55.
Sneddon
TP
Zhe
XS
Edmunds
SC
Li
P
Goodman
L
Hunter
CI
2014
GigaDB: promoting data dissemination and reproducibility
Database
 , vol. 
2014
 pg. 
bau018
  
10.1093/database/bau018

Author notes

Equal contributors
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Supplementary data