Sufficient amounts of DNA were available for exome-sequencing from 50 JMML patients (discovery cohort and validation cohort). Sequencing libraries were prepared at the DKFZ Genomics and Proteomics Core Facility using the Agilent “SureSelectXT Human all Exon V4” kit and subsequently sequenced on a HiSeq2000 instrument using the 100 bp paired-end mode. Alignment to the hs37d5 reference genome was performed using BWA-MEM63 (link). SAMtools/BCFtools (version 0.1.19) were used for single nucleotide variants (SNV) calling64 (link). SNV found with high frequency in other mutation databases (i.e., ‘common = 1’ tag in dbSNP or > 1% frequency in ExAC 0.3.1) were filtered out and only high-confidence mutations in the coding regions were kept. Calling of small insertions/deletions (INDEL) was performed by Platypus 0.8.1 using the same filtering process as with the SNV calling65 (link). Both SNV and INDEL calling were performed without using paired germline control samples. The data were plotted as an OncoPrint using the ComplexHeatmap package66 (link). Statistics between groups were performed using the non-parametric Wilcoxon's test.
Free full text: Click here