The DiscovEHR collaboration between the Regeneron Genetics Center and Geisinger Health System couples high-throughput sequencing to an integrated health care system using longitudinal electronic health records (EHRs). We sequenced the exomes of 50,726 adult participants in the DiscovEHR study to identify ~4.2 million rare single-nucleotide variants and insertion/deletion events, of which ~176,000 are predicted to result in a loss of gene function. Linking these data to EHR-derived clinical phenotypes, we find clinical associations supporting therapeutic targets, including genes encoding drug targets for lipid lowering, and identify previously unidentified rare alleles associated with lipid levels and other blood level traits. About 3.5% of individuals harbor deleterious variants in 76 clinically actionable genes. The DiscovEHR data set provides a blueprint for large-scale precision medicine initiatives and genomics-guided therapeutic discovery.

Primary ciliary dyskinesia (PCD) is caused when defects of motile cilia lead to chronic airway infections, male infertility, and situs abnormalities. Multiple causative PCD mutations account for only 65% of cases, suggesting that many genes essential for cilia function remain to be discovered. By using zebrafish morpholino knockdown of PCD candidate genes as an in vivo screening platform, we identified c21orf59, ccdc65, and c15orf26 as critical for cilia motility. c21orf59 and c15orf26 knockdown in zebrafish and planaria blocked outer dynein arm assembly, and ccdc65 knockdown altered cilia beat pattern. Biochemical analysis in Chlamydomonas revealed that the C21orf59 ortholog FBB18 is a flagellar matrix protein that accumulates specifically when cilia motility is impaired. The Chlamydomonas ida6 mutant identifies CCDC65/FAP250 as an essential component of the nexin-dynein regulatory complex. Analysis of 295 individuals with PCD identified recessive truncating mutations of C21orf59 in four families and CCDC65 in two families. Similar to findings in zebrafish and planaria, mutations in C21orf59 caused loss of both outer and inner dynein arm components. Our results characterize two genes associated with PCD-causing mutations and elucidate two distinct mechanisms critical for motile cilia function: dynein arm assembly for C21orf59 and assembly of the nexin-dynein regulatory complex for CCDC65.

Cilia are essential for fertilization, respiratory clearance, cerebrospinal fluid circulation, and to establish laterality1. Cilia motility defects cause Primary Ciliary Dyskinesia (PCD, MIM 242650), a disorder affecting 1:15-30,000 births. Cilia motility requires the assembly of multisubunit dynein arms that drive cilia bending2. Despite progress in understanding the genetic basis of PCD, mutations remain to be identified for several PCD linked loci3. Here we show that the zebrafish cilia paralysis mutant schmalhanstn222 (smh) mutant encodes the coiled-coil domain containing 103 protein (Ccdc103), a foxj1a regulated gene. Screening 146 unrelated PCD families identified patients in six families with reduced outer dynein arms, carrying mutations in CCDC103. Dynein arm assembly in smh mutant zebrafish was rescued by wild-type but not mutant human CCDC103. Chlamydomonas Ccdc103 functions as a tightly bound, axoneme-associated protein. The results identify Ccdc103 as a novel dynein arm attachment factor that when mutated causes Primary Ciliary Dyskinesia.

Genetic mapping of mutations in model systems has facilitated the identification of genes contributing to fundamental biological processes including human diseases. However, this approach has historically required the prior characterization of informative markers. Here we report a fast and cost-effective method for genetic mapping using next-generation sequencing that combines single nucleotide polymorphism discovery, mutation localization, and potential identification of causal sequence variants. In contrast to prior approaches, we have developed a hidden Markov model to narrowly define the mutation area by inferring recombination breakpoints of chromosomes in the mutant pool. In addition, we created an interactive online software resource to facilitate automated analysis of sequencing data and demonstrate its utility in the zebrafish and mouse models. Our novel methodology and online tools will make next-generation sequencing an easily applicable resource for mutation mapping in all model systems.[Supplemental material is available for this article.]There can be little argument that genetic mapping has made a substantial contribution to our understanding of biology. For many years these studies used phenotypically defined markers, such as those used by Morgan in Drosophila and Haldane in mice (Morgan 1911;Haldane et al. 1915). The modern era of genetic analysis was heralded by the recognition that variation in genomic DNA sequence itself could be used as a facile assay for mapping (Botstein et al. 1980). This was initially accomplished using analysis of restriction fragmentlength polymorphisms, which were later replaced by microsatellites and subsequently by single nucleotide polymorphisms (SNPs). Despite the remarkable technological advances, these approaches hold in common with those of Morgan and Haldane the utilization of prespecified markers. Next-generation sequencing (NGS) technology enables simultaneous discovery of very dense sets of informative markers and actual gene mapping in the same experiment. Here, we present a strategy and computational tools to map genes in model organisms using sequencing of pooled samples. The approach can be applied to any model organism with a characterized genome and also to both spontaneous and induced mutants. We demonstrate the utility of the strategy and efficiency of the computational approach by mapping spontaneous and ethylnitrosourea (ENU)-induced developmental mutants in zebrafish and mouse.Large-scale forward mutagenesis screens in zebrafish have been used with success to investigate fundamental developmental processes. While the recent completion of the zebrafish genome has greatly aided in the identification of genes, mapping analyses continue to rely on the use of traditional microsatellite markers. However, the utilization of SNPs for mapping of zebrafish mutants was proposed almost a decade ago (Stickney et al. 2002), large numbers of SNPs have been identified (Guryev et al. 2006;Bradley et al. 2007), and the application of NGS for SNP discovery and mutat...

No abstract