A custom NGS panel with 84 genes related to circadian rhythms and melanopsin (Table s1) was based on a commercial kit (RT2 Profiler PCR Array, Qiagen) and designed with the Nextera DNA Flex Library Prep (Illumina Inc., San Diego, CA). Libraries were prepared from total blood’s DNA and were sequenced as 151-bp paired-end reads on NextSeq 500 platform (Illumina Inc., San Diego, CA). BCL files were demultiplexed and converted to the FASTQ format with the Illumina standalone bcl2fastq program (v2.20.0.422). Generated reads were aligned with BWA [50 (link)] to the reference genome hg19, realignment and base quality score recalibration were performed with GATK [51 (link)] and duplicate removal with PicardTools (https://broadinstitute.github.io/picard/). Alignment and coverage statistics were collected with SAM tools [52 (link)] and GATK. Variants were called and filtered by quality with GATK UnifiedGenotyper and VariantFiltration, then annotated with RefSeq using SnpEff [53 (link)].
Free full text: Click here