Genomic DNA was extracted from whole blood or saliva using a Qiagen QIAamp DNA blood kit or tissue kit (Qiagen, Hilden, Germany). Raw WES data were generated and processed to BAM files at three sites (Table S1): The Broad Institute (ASD = 256, HC = 299) [11 (link)] using Illumina HiSeq sequencers and an Illumina Nextera exome capture kit; Yokohama City University (ASD = 51) [14 (link)] using Illumina HiSeq with Agilent Sure Select v5; and Nagoya University (ASD = 2) using Illumina HiSeq with Agilent Sure Select v5. Detailed descriptions of each sequencing method are presented elsewhere [11 (link), 14 (link), 19 (link)]. Each sample’s sequencing reads were aligned onto the human genome build 37 (GRCh37/hg19) and then aggregated into a BAM file. Genomic variant call format (gVCF) files were generated using Haplotype Caller, version 4.1, in the Genome Analysis Toolkit (GATK) [20 (link)]. SNVs and insertions/deletions (indels) were jointly called across all samples (ASD = 309, HC = 299) using the GenotypeGVCF function.
Free full text: Click here