The success of genome-wide association studies has paralleled the development of efficient genotyping technologies. We describe the development of a next-generation microarray based on the new highly-efficient Affymetrix Axiom genotyping technology that we are using to genotype individuals of European ancestry from the Kaiser Permanente Research Program on Genes, Environment and Health (RPGEH). The array contains 674,517 SNPs, and provides excellent genome-wide as well as gene-based and candidate-SNP coverage. Coverage was calculated using an approach based on imputation and cross validation. Preliminary results for the first 80,301 saliva-derived DNA samples from the RPGEH demonstrate very high quality genotypes, with sample success rates above 94% and over 98% of successful samples having SNP call rates exceeding 98%. At steady state, we have produced 462 million genotypes per week for each Axiom system. The new array provides a valuable addition to the repertoire of tools for large scale genome-wide association studies.

The Kaiser Permanente (KP) Research Program on Genes, Environment and Health (RPGEH), in collaboration with the University of California-San Francisco, undertook genome-wide genotyping of .100,000 subjects that constitute the Genetic Epidemiology Research on Adult Health and Aging (GERA) cohort. The project, which generated .70 billion genotypes, represents the first large-scale use of the Affymetrix Axiom Genotyping Solution. Because genotyping took place over a short 14-month period, creating a near-real-time analysis pipeline for experimental assay quality control and final optimized analyses was critical. Because of the multi-ethnic nature of the cohort, four different ethnic-specific arrays were employed to enhance genome-wide coverage. All assays were performed on DNA extracted from saliva samples. To improve sample call rates and significantly increase genotype concordance, we partitioned the cohort into disjoint packages of plates with similar assay contexts. Using strict QC criteria, the overall genotyping success rate was 103,067 of 109,837 samples assayed (93.8%), with a range of 92.1-95.4% for the four different arrays. Similarly, the SNP genotyping success rate ranged from 98.1 to 99.4% across the four arrays, the variation depending mostly on how many SNPs were included as single copy vs. double copy on a particular array. The high quality and large scale of genotype data created on this cohort, in conjunction with comprehensive longitudinal data from the KP electronic health records of participants, will enable a broad range of highly powered genome-wide association studies on a diversity of traits and conditions. KEYWORDS genome-wide genotyping; GERA cohort; Affymetrix Axiom; saliva DNA; quality control T HE Genetic Epidemiology Research on Adult Health and Aging (GERA) resource is a cohort of .100,000 subjects who are participants in the Kaiser Permanente Medical Care Plan, Northern California Region (KPNC), Research Program on Genes, Environment and Health (RPGEH) (detailed description of the cohort and study design can be found in dbGaP, Study Accession: phs000674.v1.p1). Genome-wide genotyping was targeted for this cohort to enable large-scale genome-wide association studies by linkage to comprehensive longitudinal clinical data derived from extensive KPNC electronic health record databases. The cohort is multi-ethnic, with 20% minority representation (African American, East Asian, and Latino or mixed), and the remaining 80% nonHispanic white. For this project, four ethnic-specific arrays were designed based on the Affymetrix Axiom Genotyping System (Hoffmann et al. 2011a,b). The genotyping assay experiment took place over a 14-month period and to our knowledge, is the single largest genotyping experiment to date, producing .70 billion genotypes. The magnitude of the experiment, in conjunction with the long duration and simultaneous high throughput, required new protocols for assuring quality control (QC) during the assays and new genotyping strategies in postassay data analysis.Samp...

Anophthalmia/microphthalmia (A/M) is a genetically heterogeneous birth defect for which the etiology is unknown in more than 50% of patients. We used exome sequencing with the ACE Exome™ (Personalis, Inc; 18 cases) and UCSF Genomics Core (21 cases) to sequence 28 patients with A/M and four patients with varied developmental eye defects. In the 28 patients with A/M, we identified de novo mutations in three patients (OTX2, p.(Gln91His), RARB, p.Arg387Cys and GDF6, p.Ala249Glu) and inherited mutations in STRA6 in two patients. In patients with developmental eye defects, a female with cataracts and cardiomyopathy had a de novo COL4A1 mutation, p.(Gly773Arg), expanding the phenotype associated with COL4A1 to include cardiomyopathy. A male with a chorioretinal defect, microcephaly, seizures and sensorineural deafness had two PNPT1 mutations, p.(Ala507Ser) and c.401-1G>A, and we describe eye defects associated with this gene for the first time. Exome sequencing was efficient for identifying mutations in pathogenic genes for which there is no clinical testing available and for identifying cases that expand phenotypic spectra, such as the PNPT1 and COL4A1-associated disorders described here.

Focal facial dermal dysplasia (FFDD) Type IV is a rare syndrome characterized by facial lesions resembling aplasia cutis in a preauricular distribution along the line of fusion of the maxillary and mandibular prominences. To identify the causative gene(s), exome sequencing was performed in a family with two affected siblings. Assuming autosomal recessive inheritance, two novel sequence variants were identified in both siblings in CYP26C1-a duplication of seven base pairs, which was maternally inherited, c.844_851dupCCATGCA, predicting p.Glu284fsX128 and a missense mutation, c.1433G>A, predicting p.Arg478His, that was paternally inherited. The duplication predicted a frameshift mutation that led to a premature stop codon and premature chain termination, whereas the missense mutation was not functional based on its in vitro expression in mammalian cells. The FFDD skin lesions arise along the sites of fusion of the maxillary and mandibular prominences early in facial development, and Cyp26c1 was expressed exactly along the fusion line for these facial prominences in the first branchial arch in mice. Sequencing of four additional, unrelated Type IV FFDD patients and eight Type II or III TWIST2-negative FFDD patients revealed that three of the Type IV patients were homozygous for the duplication, whereas none of the Type II or III patients had CYP26C1 mutations. The seven base pairs duplication was present in 0.3% of healthy controls and 0.3% of patients with other birth defects. These findings suggest that the phenotypic manifestations of FFDD Type IV can be non-penetrant or underascertained. Thus, FFDD Type IV results from the loss of function mutations in CYP26C1.

Cleft lip and/or palate (CL/P) are common structural birth defects in humans. We used exome sequencing to study a patient with bilateral CL/P and identified a single nucleotide deletion in the patient and her similarly affected son—c.546_546delG, predicting p.Gln183Argfs*57 in the Distal-less 4 (DLX4) gene. The sequence variant was absent from databases, predicted to be deleterious and was verified by Sanger sequencing. In mammals, there are three Dlx homeobox clusters with closely located gene pairs (Dlx1/Dlx2, Dlx3/Dlx4, Dlx5/Dlx6). In situ hybridization showed that Dlx4 was expressed in the mesenchyme of the murine palatal shelves at E12.5, prior to palate closure. Wild-type human DLX4, but not mutant DLX4_c.546delG, could activate two murine Dlx conserved regulatory elements, implying that the mutation caused haploinsufficiency. We showed that reduced DLX4 expression after short interfering RNA treatment in a human cell line resulted in significant up-regulation of DLX3, DLX5 and DLX6, with reduced expression of DLX2 and significant up-regulation of BMP4, although the increased BMP4 expression was demonstrated only in HeLa cells. We used antisense morpholino oligonucleotides to target the orthologous Danio rerio gene, dlx4b, and found reduced cranial size and abnormal cartilaginous elements. We sequenced DLX4 in 155 patients with non-syndromic CL/P and CP, but observed no sequence variants. From the published literature, Dlx1/Dlx2 double homozygous null mice and Dlx5 homozygous null mice both have clefts of the secondary palate. This first finding of a DLX4 mutation in a family with CL/P establishes DLX4 as a potential cause of human clefts.