PacMAGI: A pipeline including accurate indel detection for the analysis of PacBio sequencing data applied to RPE65

Abstract

Third generation sequencing methods, like PacBio, provide information about structural variants, introns, enhancers and promoters. We developed an automated pipeline, called PacMAGI, including quality control, alignment, SNV, INDELs, structural variant calling, phasing, annotation and variant interpretation, for the analysis of PacBio data for any target region. Bi-allelic mutations in the RPE65 gene are associated with different inherited retinal dystrophies, such as Leber congenital amaurosis (LCA) and retinitis pigmentosa (RP). Diagnostic panel-based NGS analysis is performed on coding regions and intron/exon junctions of genes. To obtain a more conclusive diagnosis, we applied PacMAGI to obtain a second hit on RPE65 in LCA or RP patients who showed a single heterozygous variant by NGS. We used PacBio to sequence the full gene and identify putative second-hits in intronic, problematic and promoter regions. All variants identified in the diagnostic setting with NGS were correctly detected by the pipeline, and thanks to our custom algorithm for INDELs, a previously undetected 'Pathogenic' frameshift variant was found in a RP patient already identified to carry a 'Likely Pathogenic' variant.

Keywords: Bioinformatics pipeline; Leber congenital amaurosis; PacBio sequencing; RPE65; Retinitis pigmentosa; pipeline PacBio.