High-throughput integration site profiling has become a feasible tool to assess vector biosafety and to monitor the cell fate of the gene-corrected cell population in clinical gene therapy studies. Here we report a step-by-step protocol for universal genome-wide and comprehensive integrome analysis that can be performed on >10²–10³ samples of interest in parallel. This assay is composed of fast and cost-efficient non-restrictive linear amplification–mediated PCR; optimized sample preparation for pyrosequencing; and automated bioinformatic data mining, including sequence trimming, alignment to the cellular genome and further annotation. Moreover, the workflow of this large-scale assay can be adapted to any PCR-based method aiming to characterize unknown flanking DNA adjacent to a known DNA region. Thus, in combination with next-generation sequencing technologies, large-scale integrome analysis of >4 × 10⁵–1 × 10⁶ integration site sequences can be accomplished within a single week.