Richard Lao

Next generation genome-wide association tool: Design and coverage of a high-throughput European-optimized SNP array

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.

Genotyping Informatics and Quality Control for 100,000 Subjects in the Genetic Epidemiology Research on Adult Health and Aging (GERA) Cohort

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.

Exome Sequencing in 32 Patients with Anophthalmia/Microphthalmia and Developmental Eye Defects

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 ExomeTM (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, type IV, is caused by mutations in CYP26C1

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.

Novel mutations in PXDN cause microphthalmia and anterior segment dysgenesis

We used exome sequencing to study a non-consanguineous family with two children who had anterior segment dysgenesis, sclerocornea, microphthalmia, hypotonia and developmental delays. Sanger sequencing verified two Peroxidasin (PXDN) mutations in both sibs—a maternally inherited, nonsense mutation, c.1021C>T predicting p.(Arg341*), and a paternally inherited, 23-basepair deletion causing a frameshift and premature protein truncation, c.2375_2397del23, predicting p.(Leu792Hisfs*67). We re-examined exome data from 20 other patients with structural eye defects and identified two additional PXDN mutations in a sporadic male with bilateral microphthalmia, cataracts and anterior segment dysgenesis—a maternally inherited, frameshift mutation, c.1192delT, predicting p.(Tyr398Thrfs*40) and a paternally inherited, missense substitution that was predicted to be deleterious, c.947 A>C, predicting p.(Gln316Pro). Mutations in PXDN were previously reported in three families with congenital cataracts, microcornea, sclerocornea and developmental glaucoma. The gene is expressed in corneal epithelium and is secreted into the extracellular matrix. Defective peroxidasin has been shown to impair sulfilimine bond formation in collagen IV, a constituent of the basement membrane, implying that the eye defects result because of loss of basement membrane integrity in the developing eye. Our finding of a broader phenotype than previously appreciated for PXDN mutations is typical for exome-sequencing studies, which have proven to be highly effective for mutation detection in patients with atypical presentations. We conclude that PXDN sequencing should be considered in microphthalmia with anterior segment dysgenesis.

DLX4 is associated with orofacial clefting and abnormal jaw development

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.

Whole exome and whole genome sequencing with dried blood spot DNA without whole genome amplification

Newborn screening (NBS) for rare conditions is performed in all 50 states in the USA. We have partnered with the California Department of Public Health Genetic Disease Laboratory to determine whether sufficient DNA can be extracted from archived dried blood spots (DBS) for next‐generation sequencing in the hopes that next‐generation sequencing can play a role in NBS. We optimized the DNA extraction and sequencing library preparation protocols for residual infant DBS archived over 20 years ago and successfully obtained acceptable whole exome and whole genome sequencing data. This sequencing study using DBS DNA without whole genome amplification prior to sequencing library preparation provides evidence that properly stored residual newborn DBS are a satisfactory source of DNA for genetic studies.

Two missense mutations in SALL4 in a patient with microphthalmia, coloboma, and optic nerve hypoplasia

ABSTRACT To investigate the genetic etiology of anophthalmia and microphthalmia, we used exome sequencing in a Caucasian female with unilateral microphthalmia and coloboma, bilateral optic nerve hypoplasia, ventricular and atrial septal defects, and growth delays. We found two sequence variants in SALL4 - c.[575C>A], predicting p.(Ala192Glu), that was paternally inherited, and c.[2053G>C], predicting p.(Asp685His), that was maternally inherited. Haploinsufficiency for SALL4 due to nonsense or frameshift mutations has been associated with acro-renal ocular syndrome that is characterized by eye defects including Duane anomaly and coloboma, in addition to radial ray malformations and renal abnormalities. Our report is the first description of structural eye defects associated with two missense variants in SALL4 inherited in trans; the absence of reported findings in both parents suggests that both sequence variants are hypomorphic mutations and that both are needed for the ocular phenotype. SALL4 is expressed in the developing lens and regulates BMP4, leading us to speculate that altered BMP4 expression was responsible for the eye defects, but we could not demonstrate altered BMP4 expression in vitro after using small interfering RNAs (siRNAs) to reduce SALL4 expression. We conclude that SALL4 hypomorphic variants may influence eye development.

A recurrent, non-penetrant sequence variant, p.Arg266Cys in Growth/Differentiation Factor 3 (GDF3) in a female with unilateral anophthalmia and skeletal anomalies

Purpose The genetic causes of anophthalmia, microphthalmia and coloboma remain poorly understood. Missense mutations in Growth/Differentiation Factor 3 (GDF3) gene have previously been reported in patients with microphthalmia, iridial and retinal colobomas, Klippel-Feil anomaly with vertebral fusion, scoliosis, rudimentary 12th ribs and an anomalous right temporal bone. We used whole exome sequencing with a trio approach to study a female with unilateral anophthalmia, kyphoscoliosis and additional skeletal anomalies. Observations Exome sequencing revealed that the proposita was heterozygous for c.796C > T, predicting p.Arg266Cys, in GDF3. Sanger sequencing confirmed the mutation and showed that the unaffected mother was heterozygous for the same missense substitution. Conclusions and importance Although transfection studies with the p.Arg266Cys mutation have shown that this amino acid substitution is likely to impair function, non-penetrance for the ocular defects was apparent in this family and has been observed in other families with sequence variants in GDF3. We conclude p.Arg266Cys and other GDF3 mutations can be non-penetrant, making pathogenicity more difficult to establish when sequence variants in this gene are present in patients with structural eye defects.