Long insert whole genome sequencing for copy number variant and translocation detection

Microfluidic chip immunoassay based on rolling circle amplification and G-quadruplex/Thioflavin T for multiplex detection of CTX I

A label-free immunoassay based on rolling circle amplification (RCA) and G-quadruplex/Thioflavin T (G4/ThT) is proposed to realize the sensitive detection of carboxy-terminal cross-linked fragment of type I collagen (CTX I) for bone loss. Under the optimal conditions, as low as 38.02 pg/mL of CTX I can be detected. To improve the detecting throughput and simplify the operation, a microfluidic chip was designed, fabricated, and used for CTX I detection based on the proposed assay. The detection can be completed with only a single on-chip magnetic separation step, which was easy to operate, less time-consuming, and has only low reagent consumption. The limit of detection was 131.83 pg/mL by observing with fluorescence microscope. With further improvement of detection equipment, the sensitivity of on-chip detection can be improved. It can be expected that the proposed RCA/G4/ThT immunoassay for sensitive and high-throughput automated detection of CTX I might be chosen as a potential analytical tool for clinical osteoporosis diagnosis and in-orbit bone loss detection.

Genomic Staging of Multifocal Lung Squamous Cell Carcinomas Is Independent of the Comprehensive Morphologic Assessment

INTRODUCTION

Morphologic and molecular data for staging of multifocal lung squamous cell carcinomas (LSCCs) are limited. In this study, whole exome sequencing (WES) was used as the gold standard to determine whether multifocal LSCC represented separate primary lung cancers (SPLCs) or intrapulmonary metastases (IPMs). Genomic profiles were compared with the comprehensive morphologic assessment.

METHODS

WES was performed on 20 tumor pairs of multifocal LSCC and matched normal lymph nodes using the Illumina NovaSeq6000 S4-Xp (Illumina, San Diego, CA). WES clonal and subclonal analysis data were compared with histologic assessment by 16 thoracic pathologists. In addition, the immune gene profiling of the study cases was characterized by the HTG EdgeSeq Precision Immuno-Oncology Panel.

RESULTS

By WES data, 11 cases were classified as SPLC and seven cases as IPM. Two cases were technically suboptimal. Analysis revealed marked genomic and immunogenic heterogeneity, but immune gene expression profiles highly correlated with mutation profiles. Tumors classified as IPM have a large number of shared mutations (ranging from 33.5% to 80.7%). The agreement between individual morphologic assessments for each case and WES was 58.3%. One case was unanimously interpreted morphologically as IPM and was in agreement with WES. In a further 17 cases, the number of pathologists whose morphologic interpretation was in agreement with WES ranged from two (one case) to 15 pathologists (one case) per case. Pathologists showed a fair interobserver agreement in the morphologic staging of multiple LSCCs, with an overall kappa of 0.232.

CONCLUSIONS

Staging of multifocal LSCC based on morphologic assessment is unreliable. Comprehensive genomic analyses should be adopted for the staging of multifocal LSCC.

Resolving Breakpoints of Chromosomal Rearrangements at the Nucleotide Level Using Sanger Sequencing

Novel cytogenetic tools are increasingly based on genome sequencing for detecting chromosomal abnormalities. Different sequence-based techniques optimized for diagnosis of structural variants can be useful for narrowing down the localization of breakpoints of chromosomal abnormalities, but do not offer nucleotide resolution of breakpoints for proper interpretation of gene disruption. This protocol presents the characterization of structural variants at nucleotide resolution using Sanger sequencing after low-pass large-insert genome sequencing or other long-molecule methods. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Primer design for junction amplification at translocations and inversions Basic Protocol 2: Amplification of derivative chromosomes using a long-range polymerase Alternate Protocol: Amplification of derivative chromosomes using a hot-start polymerase Basic Protocol 3: Preparation of DNA for Sanger sequencing Basic Protocol 4: Interpretation and reporting of breakpoints based on Sanger sequencing.

Single-cell sequencing of genomic DNA resolves sub-clonal heterogeneity in a melanoma cell line

We performed shallow single-cell sequencing of genomic DNA across 1475 cells from a cell-line, COLO829, to resolve overall complexity and clonality. This melanoma tumor-line has been previously characterized by multiple technologies and is a benchmark for evaluating somatic alterations. In some of these studies, COLO829 has shown conflicting and/or indeterminate copy number and, thus, single-cell sequencing provides a tool for gaining insight. Following shallow single-cell sequencing, we first identified at least four major sub-clones by discriminant analysis of principal components of single-cell copy number data. Based on clustering, break-point and loss of heterozygosity analysis of aggregated data from sub-clones, we identified distinct hallmark events that were validated within bulk sequencing and spectral karyotyping. In summary, COLO829 exhibits a classical Dutrillaux’s monosomic/trisomic pattern of karyotype evolution with endoreduplication, where consistent sub-clones emerge from the loss/gain of abnormal chromosomes. Overall, our results demonstrate how shallow copy number profiling can uncover hidden biological insights.

Through shallow single-cell sequencing of genomic DNA followed by clustering analysis, Velazquez-Villarreal et al. reveal sub-clones of the melanoma cell line COLO829 and further identify and validate chromosome translocations and copy number changes. This study illustrates how copy number variation analysis can provide insights into cancer cell heterogeneity.

Clonality analysis of pulmonary tumors by genome-wide copy number profiling

Multiple tumors in patients are frequently diagnosed, either synchronous or metachronous. The distinction between a second primary and a metastasis is important for treatment. Chromosomal DNA copy number aberrations (CNA) patterns are highly unique to specific tumors. The aim of this study was to assess genome-wide CNA-patterns as method to identify clonally related tumors in a prospective cohort of patients with synchronous or metachronous tumors, with at least one intrapulmonary tumor. In total, 139 tumor pairs from 90 patients were examined: 35 synchronous and 104 metachronous pairs. Results of CNA were compared to histological type, clinicopathological methods (Martini-Melamed-classification (MM) and ACCP-2013-criteria), and, if available, EGFR- and KRAS-mutation analysis. CNA-results were clonal in 74 pairs (53%), non-clonal in 33 pairs (24%), and inconclusive in 32 pairs (23%). Histological similarity was found in 130 pairs (94%). Concordance between histology and conclusive CNA-results was 69% (74 of 107 pairs: 72 clonal and two non-clonal). In 31 of 103 pairs with similar histology, genetics revealed non-clonality. In two out of four pairs with non-matching histology, genetics revealed clonality. The subgroups of synchronous and metachronous pairs showed similar outcome for the comparison of histological versus CNA-results. MM-classification and ACCP-2013-criteria, applicable on 34 pairs, and CNA-results were concordant in 50% and 62% respectively. Concordance between mutation matching and conclusive CNA-results was 89% (8 of 9 pairs: six clonal and two non-clonal). Interestingly, in one patient both tumors had the same KRAS mutation, but the CNA result was non-clonal. In conclusion, although some concordance between histological comparison and CNA profiling is present, arguments exist to prefer extensive molecular testing to determine whether a second tumor is a metastasis or a second primary.

Whole genome paired-end sequencing elucidates functional and phenotypic consequences of balanced chromosomal rearrangement in patients with developmental disorders

BACKGROUND

Balanced chromosomal rearrangements associated with abnormal phenotype are rare events, but may be challenging for genetic counselling, since molecular characterisation of breakpoints is not performed routinely. We used next-generation sequencing to characterise breakpoints of balanced chromosomal rearrangements at the molecular level in patients with intellectual disability and/or congenital anomalies.

METHODS

Breakpoints were characterised by a paired-end low depth whole genome sequencing (WGS) strategy and validated by Sanger sequencing. Expression study of disrupted and neighbouring genes was performed by RT-qPCR from blood or lymphoblastoid cell line RNA.

RESULTS

Among the 55 patients included (41 reciprocal translocations, 4 inversions, 2 insertions and 8 complex chromosomal rearrangements), we were able to detect 89% of chromosomal rearrangements (49/55). Molecular signatures at the breakpoints suggested that DNA breaks arose randomly and that there was no major influence of repeated elements. Non-homologous end-joining appeared as the main mechanism of repair (55% of rearrangements). A diagnosis could be established in 22/49 patients (44.8%), 15 by gene disruption (KANSL1, FOXP1, SPRED1, TLK2, MBD5, DMD, AUTS2, MEIS2, MEF2C, NRXN1, NFIX, SYNGAP1, GHR, ZMIZ1) and 7 by position effect (DLX5, MEF2C, BCL11B, SATB2, ZMIZ1). In addition, 16 new candidate genes were identified. Systematic gene expression studies further supported these results. We also showed the contribution of topologically associated domain maps to WGS data interpretation.

CONCLUSION

Paired-end WGS is a valid strategy and may be used for structural variation characterisation in a clinical setting.

Identification of therapeutic targets in chordoma through comprehensive genomic and transcriptomic analyses

Chordoma is a rare, orphan cancer arising from embryonal precursors of bone. Surgery and radiotherapy (RT) provide excellent local control, often at the price of significant morbidity because of the structures involved and the need for relatively high doses of RT; however, recurrence remains high. Although our understanding of the genetic changes that occur in chordoma is evolving rapidly, this knowledge has yet to translate into treatments. We performed comprehensive DNA (paired tumor/normal whole-exome and shallow whole-genome) and RNA (tumor whole-transcriptome) next-generation sequencing analyses of archival sacral and clivus chordoma specimens. Incorporation of transcriptomic data enabled the identification of gene overexpression and expressed DNA alterations, thus providing additional support for potential therapeutic targets. In three patients, we identified alterations that may be amenable to off-label FDA-approved treatments for other tumor types. These alterations include FGFR1 overexpression (ponatinib, pazopanib) and copy-number duplication of CDK4 (palbociclib) and ERBB3 (gefitinib). In a third patient, germline DNA demonstrated predicted pathogenic changes in CHEK2 and ATM, which may have predisposed the patient to developing chordoma at a young age and may also be associated with potential sensitivity to PARP inhibitors because of homologous recombination repair deficiency. Last, in the fourth patient, a missense mutation in IGF1R was identified, suggesting potential activity for investigational anti-IGF1R strategies. Our findings demonstrate that chordoma patients present with aberrations in overlapping pathways. These results provide support for targeting the IGF1R/FGFR/EGFR and CDK4/6 pathways as treatment strategies for chordoma patients. This study underscores the value of comprehensive genomic and transcriptomic analysis in the development of rational, individualized treatment plans for chordoma.

Somatic inactivating PTPRJ mutations and dysregulated pathways identified in canine malignant melanoma by integrated comparative genomic analysis

Canine malignant melanoma, a significant cause of mortality in domestic dogs, is a powerful comparative model for human melanoma, but little is known about its genetic etiology. We mapped the genomic landscape of canine melanoma through multi-platform analysis of 37 tumors (31 mucosal, 3 acral, 2 cutaneous, and 1 uveal) and 17 matching constitutional samples including long- and short-insert whole genome sequencing, RNA sequencing, array comparative genomic hybridization, single nucleotide polymorphism array, and targeted Sanger sequencing analyses. We identified novel predominantly truncating mutations in the putative tumor suppressor gene PTPRJ in 19% of cases. No BRAF mutations were detected, but activating RAS mutations (24% of cases) occurred in conserved hotspots in all cutaneous and acral and 13% of mucosal subtypes. MDM2 amplifications (24%) and TP53 mutations (19%) were mutually exclusive. Additional low-frequency recurrent alterations were observed amidst low point mutation rates, an absence of ultraviolet light mutational signatures, and an abundance of copy number and structural alterations. Mutations that modulate cell proliferation and cell cycle control were common and highlight therapeutic axes such as MEK and MDM2 inhibition. This mutational landscape resembles that seen in BRAF wild-type and sun-shielded human melanoma subtypes. Overall, these data inform biological comparisons between canine and human melanoma while suggesting actionable targets in both species.

Whole Genome Library Construction for Next Generation Sequencing

With the rapid evolution of genomics technologies over the past decade, whole genome sequencing (WGS) has become an increasingly accessible tool in biomedical research. WGS applications include analysis of genomic DNA from single individuals, multiple related family members, and tumor/normal samples from the same patient in the context of oncology. A number of different modalities are available for performing WGS; this chapter focuses on wet lab library construction procedures for complex short insert WGS libraries using the KAPA Hyper Prep Kit (Kapa Biosystems), and includes a discussion of appropriate quality control measures for sequencing on the Illumina HiSeq2000 platform. Additional modifications to the protocol for long insert WGS library construction, to assess structural alterations and copy number changes, are also described.

Cat-D: a targeted sequencing method for the simultaneous detection of small DNA mutations and large DNA deletions with flexible boundaries

We developed a targeted DNA sequencing method that is capable of detecting a comprehensive panel of DNA mutations including small DNA mutations and large DNA deletions with unknown/flexible boundaries. The method directly identifies the large DNA deletions (Cat-D) without relying on sequencing coverage to make the genotype calls. We performed the method to simultaneously detect 10 small DNA mutations in β-thalassemia and 2 large genomic deletions in α-thalassemia from 10 genomic DNA samples. Cat-D was performed on 8 genomic DNA samples in duplicate. The 18 Cat-D samples were combined in one sequencing run. In total, 216 genotype calls were made, and 215 of the genotype calls were accurate. No false negative genotype calls were made. One false positive genotype call was made on one target mutation in one experimental duplicate from a genomic DNA sample. In summary, Cat-D can be developed into a robust, high-throughput and cost-effective method suitable for population-based carrier screens.

Prospective Feasibility Trial for Genomics-Informed Treatment in Recurrent and Progressive Glioblastoma

Purpose: Glioblastoma is an aggressive and molecularly heterogeneous cancer with few effective treatment options. We hypothesized that next-generation sequencing can be used to guide treatment recommendations within a clinically acceptable time frame following surgery for patients with recurrent glioblastoma.Experimental Design: We conducted a prospective genomics-informed feasibility trial in adults with recurrent and progressive glioblastoma. Following surgical resection, genome-wide tumor/normal exome sequencing and tumor RNA sequencing were performed to identify molecular targets for potential matched therapy. A multidisciplinary molecular tumor board issued treatment recommendations based on the genomic results, blood-brain barrier penetration of the indicated therapies, drug-drug interactions, and drug safety profiles. Feasibility of generating genomics-informed treatment recommendations within 35 days of surgery was assessed.Results: Of the 20 patients enrolled in the study, 16 patients had sufficient tumor tissue for analysis. Exome sequencing was completed for all patients, and RNA sequencing was completed for 14 patients. Treatment recommendations were provided within the study's feasibility time frame for 15 of 16 (94%) patients. Seven patients received treatment based on the tumor board recommendations. Two patients reached 12-month progression-free survival, both adhering to treatments based on the molecular profiling results. One patient remained on treatment and progression free 21 months after surgery, 3 times longer than the patient's previous time to progression. Analysis of matched nonenhancing tissue from 12 patients revealed overlapping as well as novel putatively actionable genomic alterations.Conclusions: Use of genome-wide molecular profiling is feasible and can be informative for guiding real-time, central nervous system-penetrant, genomics-informed treatment recommendations for patients with recurrent glioblastoma. Clin Cancer Res; 24(2); 295-305. ©2017 AACRSee related commentary by Wick and Kessler, p. 256.

Integrated genomic analyses reveal frequent TERT aberrations in acral melanoma

Genomic analyses of cutaneous melanoma (CM) have yielded biological and therapeutic insights, but understanding of non-ultraviolet (UV)-derived CMs remains limited. Deeper analysis of acral lentiginous melanoma (ALM), a rare sun-shielded melanoma subtype associated with worse survival than CM, is needed to delineate non-UV oncogenic mechanisms. We thus performed comprehensive genomic and transcriptomic analysis of 34 ALM patients. Unlike CM, somatic alterations were dominated by structural variation and absence of UV-derived mutation signatures. Only 38% of patients demonstrated driver BRAF/NRAS/NF1 mutations. In contrast with CM, we observed PAK1 copy gains in 15% of patients, and somatic TERT translocations, copy gains, and missense and promoter mutations, or germline events, in 41% of patients. We further show that in vitro TERT inhibition has cytotoxic effects on primary ALM cells. These findings provide insight into the role of TERT in ALM tumorigenesis and reveal preliminary evidence that TERT inhibition represents a potential therapeutic strategy in ALM.

Case report: whole exome sequencing of primary cardiac angiosarcoma highlights potential for targeted therapies

Background

Primary cardiac angiosarcomas are rare, but they are the most aggressive type of primary cardiac neoplasms. When patients do present, it is with advanced pulmonary and/or cardiac symptoms. Therefore, many times the correct diagnosis is not made at the time of initial presentation. These patients have metastatic disease and the vast majority of these patients die within a few months after diagnosis. Currently the treatment choices are limited and there are no targeted therapies available.

Case presentation

A 56-year-old male presented with shortness of breath, night sweats, and productive cough for a month. Workup revealed pericardial effusion and multiple bilateral pulmonary nodules suspicious for metastatic disease. Transthoracic echocardiogram showed a large pericardial effusion and a large mass in the base of the right atrium. Results of biopsy of bilateral lung nodules established a diagnosis of primary cardiac angiosarcoma. Aggressive pulmonary disease caused rapid deterioration; the patient went on hospice and subsequently died. Whole exome sequencing of the patient’s postmortem tumor revealed a novel KDR (G681R) mutation, and focal high-level amplification at chromosome 1q encompassing MDM4, a negative regulator of TP53.

Conclusion

Mutations in KDR have been reported previously in angiosarcomas. Previous studies also demonstrated that KDR mutants with constitutive KDR activation could be inhibited with specific KDR inhibitors in vitro. Thus, patients harboring activating KDR mutations could be candidates for treatment with KDR-specific inhibitors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-3000-z) contains supplementary material, which is available to authorized users.

Clinical Implementation of Integrated Genomic Profiling in Patients with Advanced Cancers

DNA focused panel sequencing has been rapidly adopted to assess therapeutic targets in advanced/refractory cancer. Integrated Genomic Profiling (IGP) utilising DNA/RNA with tumour/normal comparisons in a Clinical Laboratory Improvement Amendments (CLIA) compliant setting enables a single assay to provide: therapeutic target prioritisation, novel target discovery/application and comprehensive germline assessment. A prospective study in 35 advanced/refractory cancer patients was conducted using CLIA-compliant IGP. Feasibility was assessed by estimating time to results (TTR), prioritising/assigning putative therapeutic targets, assessing drug access, ascertaining germline alterations, and assessing patient preferences/perspectives on data use/reporting. Therapeutic targets were identified using biointelligence/pathway analyses and interpreted by a Genomic Tumour Board. Seventy-five percent of cases harboured 1–3 therapeutically targetable mutations/case (median 79 mutations of potential functional significance/case). Median time to CLIA-validated results was 116 days with CLIA-validation of targets achieved in 21/22 patients. IGP directed treatment was instituted in 13 patients utilising on/off label FDA approved drugs (n = 9), clinical trials (n = 3) and single patient IND (n = 1). Preliminary clinical efficacy was noted in five patients (two partial response, three stable disease). Although barriers to broader application exist, including the need for wider availability of therapies, IGP in a CLIA-framework is feasible and valuable in selection/prioritisation of anti-cancer therapeutic targets.

Sequencing Structural Variants in Cancer for Precision Therapeutics

The identification of mutations that guide therapy selection for patients with cancer is now routine in many clinical centres. The majority of assays used for solid tumour profiling use DNA sequencing to interrogate somatic point mutations because they are relatively easy to identify and interpret. Many cancers, however, including high-grade serous ovarian, oesophageal, and small-cell lung cancer, are driven by somatic structural variants that are not measured by these assays. Therefore, there is currently an unmet need for clinical assays that can cheaply and rapidly profile structural variants in solid tumours. In this review we survey the landscape of 'actionable' structural variants in cancer and identify promising detection strategies based on massively-parallel sequencing.

A somatic reference standard for cancer genome sequencing

Large-scale multiplexed identification of somatic alterations in cancer has become feasible with next generation sequencing (NGS). However, calibration of NGS somatic analysis tools has been hampered by a lack of tumor/normal reference standards. We thus performed paired PCR-free whole genome sequencing of a matched metastatic melanoma cell line (COLO829) and normal across three lineages and across separate institutions, with independent library preparations, sequencing, and analysis. We generated mean mapped coverages of 99X for COLO829 and 103X for the paired normal across three institutions. Results were combined with previously generated data allowing for comparison to a fourth lineage on earlier NGS technology. Aggregate variant detection led to the identification of consensus variants, including key events that represent hallmark mutation types including amplified BRAF V600E, a CDK2NA small deletion, a 12 kb PTEN deletion, and a dinucleotide TERT promoter substitution. Overall, common events include >35,000 point mutations, 446 small insertion/deletions, and >6,000 genes affected by copy number changes. We present this reference to the community as an initial standard for enabling quantitative evaluation of somatic mutation pipelines across institutions.

Rare variant discovery by deep whole-genome sequencing of 1,070 Japanese individuals

The Tohoku Medical Megabank Organization reports the whole-genome sequences of 1,070 healthy Japanese individuals and construction of a Japanese population reference panel (1KJPN). Here we identify through this high-coverage sequencing (32.4 × on average), 21.2 million, including 12 million novel, single-nucleotide variants (SNVs) at an estimated false discovery rate of <1.0%. This detailed analysis detected signatures for purifying selection on regulatory elements as well as coding regions. We also catalogue structural variants, including 3.4 million insertions and deletions, and 25,923 genic copy-number variants. The 1KJPN was effective for imputing genotypes of the Japanese population genome wide. These data demonstrate the value of high-coverage sequencing for constructing population-specific variant panels, which covers 99.0% SNVs of minor allele frequency ≥0.1%, and its value for identifying causal rare variants of complex human disease phenotypes in genetic association studies.

The Tohoku Medical Megabank Organization establishes a biobank with detailed patient health care and genome information. Here the authors analyse whole-genome sequences of 1,070 Japanese individuals, allowing them to catalogue 21 million single-nucleotide variants including 12 million novel ones.

Evidence Suggesting That Discontinuous Dosing of ALK Kinase Inhibitors May Prolong Control of ALK+ Tumors

The anaplastic lymphoma kinase (ALK) is chromosomally rearranged in a subset of certain cancers, including 2% to 7% of non-small cell lung cancers (NSCLC) and ∼70% of anaplastic large cell lymphomas (ALCL). The ALK kinase inhibitors crizotinib and ceritinib are approved for relapsed ALK(+) NSCLC, but acquired resistance to these drugs limits median progression-free survival on average to ∼10 months. Kinase domain mutations are detectable in 25% to 37% of resistant NSCLC samples, with activation of bypass signaling pathways detected frequently with or without concurrent ALK mutations. Here we report that, in contrast to NSCLC cells, drug-resistant ALCL cells show no evidence of bypassing ALK by activating alternate signaling pathways. Instead, drug resistance selected in this setting reflects upregulation of ALK itself. Notably, in the absence of crizotinib or ceritinib, we found that increased ALK signaling rapidly arrested or killed cells, allowing a prolonged control of drug-resistant tumors in vivo with the administration of discontinuous rather than continuous regimens of drug dosing. Furthermore, even when drug resistance mutations were detected in the kinase domain, overexpression of the mutant ALK was toxic to tumor cells. We confirmed these findings derived from human ALCL cells in murine pro-B cells that were transformed to cytokine independence by ectopic expression of an activated NPM-ALK fusion oncoprotein. In summary, our results show how ALK activation functions as a double-edged sword for tumor cell viability, with potential therapeutic implications.

Translating next-generation sequencing from clinical trials to clinical practice for the treatment of advanced cancers

Next-generation sequencing (NGS) is being applied in oncology care to identify specific molecular aberrations of patient's tumors. The use of NGS now allows for sequencing entire human genomes within a reasonable cost and practical time frames for treatment decision making. Further delineation of epigenetics, transcriptomics, metagenomics and NGS at the level of circulating tumor DNA reveal ever increasing complexity to understand these interactions and the roles they play in cancer. With the improvement in understanding the study of proteomics, it has become clear that NGS has room for innovation to someday include sequencing of proteins. Early embarkation of NGS incorporated into clinical trials has begun. Here, we review the feasibility and practicality of translating NGS from clinical trials to clinical practice.

Non-random DNA fragmentation in next-generation sequencing

Next Generation Sequencing (NGS) technology is based on cutting DNA into small fragments, and their massive parallel sequencing. The multiple overlapping segments termed “reads” are assembled into a contiguous sequence. To reduce sequencing errors, every genome region should be sequenced several dozen times. This sequencing approach is based on the assumption that genomic DNA breaks are random and sequence-independent. However, previously we showed that for the sonicated restriction DNA fragments the rates of double-stranded breaks depend on the nucleotide sequence. In this work we analyzed genomic reads from NGS data and discovered that fragmentation methods based on the action of the hydrodynamic forces on DNA, produce similar bias. Consideration of this non-random DNA fragmentation may allow one to unravel what factors and to what extent influence the non-uniform coverage of various genomic regions.