NGS was performed using the Ion GeneStudio™ S5 Plus platform (Thermo Fisher Scientific, Massachusetts, USA). Libraries were amplified by the Ion AmpliSeq™ Library kit Plus (Thermo Fisher Scientific) and quantified using the Qubit 4 Fluorometer and 2100 Bioanalyzer with dsDNA HS assay and High Sensitivity DNA kit (Thermo Fisher Scientific), respectively. Sequencing data were processed with the Ion Torrent Suite software v.5.10.
Positive calls were selected with a read depth >30X and allelic frequency higher than 0.3 (range 0–1.0). Reads were aligned to human genome sequence (build GRCh37/human genome 19). Variants were collected using Variant Caller plugin and systematically evaluated, filtered, and annotated using tailored R scripts. Variants’ annotations were finally verified using ANNOVAR.