Genomic DNA was extracted from peripheral blood samples using the blood genomic extraction kit (Qiagen, Germany) and subjected to targeted next-generation sequencing. The coding exons and flanking sequences of over 20,000 genes were captured for paired-end sequencing (2× 100 bp) on Illumina HiSeq2500 (Illumina, USA). The average sequencing depth of the target region was higher than 170×, with over 95% of bases covered by at least 30×. An in-house genome analysis pipeline composed of BWA 0.7.12 (Li and Durbin, 2009 (link)), GATK 3.5 (McKenna et al., 2010 (link); DePristo et al., 2011 (link)), snpEff 4.0 (Cingolani et al., 2012 (link)), BED Tools 2.25.0 (Quinlan, 2014 ), ANNOVA (Yang and Wang, 2015 (link)) were used to process demultiplexed fastq data. BAM files were generated to visualize read pairs and variant calling in IGV (Broad Institute, USA). Candidate variants were filtered based on the population frequency of 2%, predicated variant pathogenicity, inheritance mode, and protein interaction network.
Free full text: Click here