@article {750, title = {Exome sequencing reveals a novel variant in causing intracranial aneurysm in a Chinese family.}, journal = {J Neurointerv Surg}, volume = {12}, year = {2020}, month = {2020 Feb}, pages = {221-226}, abstract = {

BACKGROUND: Genetic risk factors play an important role in the pathogenesis of familial intracranial aneurysms (FIAs); however, the molecular mechanisms remain largely unknown.

OBJECTIVE: To investigate potential FIA-causing genetic variants by rare variant interrogation and a family-based genomics approach in a large family with an extensive multigenerational pedigree with FIAs.

METHOD: Exome sequencing (ES) was performed in a dominant likely family with intracranial aneurysms (IAs). Variants were analyzed by an in-house developed pipeline and prioritized using various filtering strategies, including population frequency, variant type, and predicted variant pathogenicity. Sanger sequencing was also performed to evaluate the segregation of the variants with the phenotype.

RESULTS: Based on the ES data obtained from five individuals from a family with 7/21 living members affected with IAs, a total of 14 variants were prioritized as candidate variants. Familial segregation analysis revealed that c.2519T>C (p.Leu840Pro) segregated in accordance with Mendelian expectations with the phenotype within the family-that is, present in all IA-affected cases and absent from all unaffected members of the second generation. This missense variant is absent from public databases (1000genome, ExAC, gnomAD, ESP5400), and has damaging predictions by bioinformatics tools (Gerp ++ score = 5.88, CADD score = 16.43, MutationTaster score = 1, LRT score = 0). In addition, 840Leu in NFX1 is robustly conserved in mammals and maps in a region before the RING-type zinc finger domain.

CONCLUSION: c.2519T>C (p.Leu840Pro) may contribute to the pathogenetics of a subset of FIAs.

}, keywords = {Adult, Asian Continental Ancestry Group, Computational Biology, Exome, Female, Genetic Variation, Humans, Intracranial Aneurysm, Male, Middle Aged, Pedigree, Phenotype, Repressor Proteins}, issn = {1759-8486}, doi = {10.1136/neurintsurg-2019-014900}, author = {Ding, Xinghuan and Zhao, Sen and Zhang, Qianqian and Yan, Zihui and Wang, Yang and Wu, Yong and Li, Xiaoxin and Liu, Jian and Niu, Yuchen and Zhang, Yisen and Zhang, Mingqi and Wang, Huizi and Zhang, Ying and Chen, Weisheng and Yang, Xin-Zhuang and Liu, Pengfei and Posey, Jennifer E and Lupski, James R and Wu, Zhihong and Yang, Xinjian and Wu, Nan and Wang, Kun} } @article {872, title = {Rapid, Paralog-Sensitive CNV Analysis of 2457 Human Genomes Using QuicK-mer2.}, journal = {Genes (Basel)}, volume = {11}, year = {2020}, month = {2020 01 29}, abstract = {

Gene duplication is a major mechanism for the evolution of gene novelty, and copy-number variation makes a major contribution to inter-individual genetic diversity. However, most approaches for studying copy-number variation rely upon uniquely mapping reads to a genome reference and are unable to distinguish among duplicated sequences. Specialized approaches to interrogate specific paralogs are comparatively slow and have a high degree of computational complexity, limiting their effective application to emerging population-scale data sets. We present QuicK-mer2, a self-contained, mapping-free approach that enables the rapid construction of paralog-specific copy-number maps from short-read sequence data. This approach is based on the tabulation of unique k-mer sequences from short-read data sets, and is able to analyze a 20X coverage human genome in approximately 20 min. We applied our approach to newly released sequence data from the 1000 Genomes Project, constructed paralog-specific copy-number maps from 2457 unrelated individuals, and uncovered copy-number variation of paralogous genes. We identify nine genes where none of the analyzed samples have a copy number of two, 92 genes where the majority of samples have a copy number other than two, and describe rare copy number variation effecting multiple genes at the APOBEC3 locus.

}, keywords = {Algorithms, Computational Biology, DNA Copy Number Variations, Evolution, Molecular, Gene Duplication, Genome, Human, Humans, Sequence Analysis, DNA}, issn = {2073-4425}, doi = {10.3390/genes11020141}, author = {Shen, Feichen and Kidd, Jeffrey M} } @article {703, title = {COL4A1 mutations as a potential novel cause of autosomal dominant CAKUT in humans.}, journal = {Hum Genet}, volume = {138}, year = {2019}, month = {2019 Oct}, pages = {1105-1115}, abstract = {

Congenital anomalies of the kidney and urinary tract (CAKUT) are the most common cause of chronic kidney disease (~ 45\%) that manifests before 30~years of age. The genetic locus containing COL4A1 (13q33-34) has been implicated in vesicoureteral reflux (VUR), but mutations in COL4A1 have not been reported in CAKUT. We hypothesized that COL4A1 mutations cause CAKUT in humans. We performed whole exome sequencing (WES) in 550 families with CAKUT. As negative control cohorts we used WES sequencing data from patients with nephronophthisis (NPHP) with no genetic cause identified (n = 257) and with nephrotic syndrome (NS) due to monogenic causes (n = 100). We identified a not previously reported heterozygous missense variant in COL4A1 in three siblings with isolated VUR. When examining 549 families with CAKUT, we identified nine additional different heterozygous missense mutations in COL4A1 in 11 individuals from 11 unrelated families with CAKUT, while no COL4A1 mutations were identified in a control cohort with NPHP and only one in the cohort with NS. Most individuals (12/14) had isolated CAKUT with no extrarenal features. The predominant phenotype was VUR (9/14). There were no clinical features of the COL4A1-related disorders (e.g., HANAC syndrome, porencephaly, tortuosity of retinal arteries). Whereas COL4A1-related disorders are typically caused by glycine substitutions in the collagenous domain (84.4\% of variants), only one variant in our cohort is a glycine substitution within the collagenous domain (1/10). We identified heterozygous COL4A1 mutations as a potential novel autosomal dominant cause of CAKUT that is allelic to the established COL4A1-related disorders and predominantly caused by non-glycine substitutions.

}, keywords = {Alleles, Amino Acid Substitution, Collagen Type IV, Computational Biology, Congenital Abnormalities, Databases, Genetic, DNA Mutational Analysis, Evolution, Molecular, Female, Genetic Association Studies, Genetic Loci, Genomics, Heterozygote, Humans, Kidney, Kidney Diseases, Cystic, Male, Mutation, Nephrotic Syndrome, Phenotype, Urinary Tract, Web Browser, Whole Exome Sequencing}, issn = {1432-1203}, doi = {10.1007/s00439-019-02042-4}, author = {Kitzler, Thomas M and Schneider, Ronen and Kohl, Stefan and Kolvenbach, Caroline M and Connaughton, Dervla M and Dai, Rufeng and Mann, Nina and Nakayama, Makiko and Majmundar, Amar J and Wu, Chen-Han W and Kari, Jameela A and El Desoky, Sherif M and Senguttuvan, Prabha and Bogdanovic, Radovan and Staji{\'c}, Natasa and Valivullah, Zaheer and Lek, Monkol and Mane, Shrikant and Lifton, Richard P and Tasic, Velibor and Shril, Shirlee and Hildebrandt, Friedhelm} } @article {512, title = {Confirmation of the Role of DHX38 in the Etiology of Early-Onset Retinitis Pigmentosa.}, journal = {Invest Ophthalmol Vis Sci}, volume = {59}, year = {2018}, month = {2018 09 04}, pages = {4552-4557}, abstract = {

Purpose: Retinitis pigmentosa (RP) is a genetically heterogeneous trait with autosomal-recessive (ar) inheritance underlying 50\% of genetic disease cases. Sixty-one arRP genes have been identified, and recently, DHX38 has been reported as a potential candidate gene for arRP with only a single family reported with a variant of unknown significance. We identified a missense variant in DHX38 that co-segregates with the arRP phenotype in two Pakistani families confirming the involvement of DHX38 in the etiology of early-onset RP.

Methods: Exome sequencing was performed using two DNA samples from affected members of Pakistani families (MA88 and MA157) with early onset arRP. Sanger sequencing of DNA samples from all family members confirmed the segregation of candidate variant within both families.

Results: A novel missense DHX38 variant c.971G>A; p.(Arg324Gln) was identified which segregates with the arRP phenotype and yielded a logarithm of the odds (LOD) score of 5.0 and 4.3 for families MA88 and MA157, respectively. This variant is predicted to be conserved and deleterious by several bioinformatics tools.

Conclusions: We identified a second deleterious DHX38 variant that segregates with arRP in two families, providing additional evidence that DHX38 is involved in RP etiology. DHX38 encodes for pre-mRNA splicing factor PRP16, which is important in catalyzing pre-mRNA splicing.

}, keywords = {Adolescent, Adult, Cataract, Computational Biology, DEAD-box RNA Helicases, Female, Genes, Recessive, Genetic Association Studies, Genetic Linkage, Humans, Male, Mutation, Missense, Nucleotide Mapping, Ophthalmoscopy, Pedigree, Retinitis Pigmentosa, RNA Splicing Factors, Sequence Analysis, DNA, Whole Exome Sequencing, Young Adult}, issn = {1552-5783}, doi = {10.1167/iovs.18-23849}, author = {Latif, Zahid and Chakchouk, Imen and Schrauwen, Isabelle and Lee, Kwanghyuk and Santos-Cortez, Regie Lyn P and Abbe, Izoduwa and Acharya, Anushree and Jarral, Afeefa and Ali, Imran and Ullah, Ehsan and Khan, Muhammad Nasim and Ali, Ghazanfar and Tahir, Tufail Hussain and Bamshad, Michael J and Nickerson, Deborah A and Ahmad, Wasim and Ansar, Muhammad and Leal, Suzanne M} } @article {551, title = {Ximmer: a system for improving accuracy and consistency of CNV calling from exome data.}, journal = {Gigascience}, volume = {7}, year = {2018}, month = {2018 10 01}, abstract = {

Background: While exome and targeted next-generation DNA sequencing are primarily used for detecting single nucleotide changes and small indels, detection of copy number variants (CNVs) can provide highly valuable additional information from the data. Although there are dozens of exome CNV detection methods available, these are often difficult to use, and accuracy varies unpredictably between and within datasets.

Findings: We present Ximmer, a tool that supports an end-to-end process for evaluating, tuning, and running analysis methods for detection of CNVs in germline samples. Ximmer includes a simulation framework, implementations of several commonly used CNV detection methods, and a visualization and curation tool that together enable interactive exploration and quality control of CNV results. Using Ximmer, we comprehensively evaluate CNV detection on four datasets using five different detection methods. We show that application of Ximmer can improve accuracy and aid in quality control of CNV detection results. In addition, Ximmer can be used to run analyses and explore CNV results in exome data.

Conclusions: Ximmer offers a comprehensive tool and method for applying and improving accuracy of CNV detection methods for exome data.

}, keywords = {Algorithms, Computational Biology, Computer Simulation, DNA Copy Number Variations, Exome, High-Throughput Nucleotide Sequencing, Humans, Reproducibility of Results, ROC Curve, Sensitivity and Specificity, Software}, issn = {2047-217X}, doi = {10.1093/gigascience/giy112}, author = {Sadedin, Simon P and Ellis, Justine A and Masters, Seth L and Oshlack, Alicia} } @article {411, title = {denovo-db: a compendium of human de novo variants.}, journal = {Nucleic Acids Res}, volume = {45}, year = {2017}, month = {2017 01 04}, pages = {D804-D811}, abstract = {

Whole-exome and whole-genome sequencing have facilitated the large-scale discovery of de novo variants in human disease. To date, most de novo discovery through next-generation sequencing focused on congenital heart disease and neurodevelopmental disorders (NDDs). Currently, de novo variants are one of the most significant risk factors for NDDs with a substantial overlap of genes involved in more than one NDD. To facilitate better usage of published data, provide standardization of annotation, and improve accessibility, we created denovo-db (http://denovo-db.gs.washington.edu), a database for human de novo variants. As of July 2016, denovo-db contained 40 different studies and 32,991 de novo variants from 23,098 trios. Database features include basic variant information (chromosome location, change, type); detailed annotation at the transcript and protein levels; severity scores; frequency; validation status; and, most importantly, the phenotype of the individual with the variant. We included a feature on our browsable website to download any query result, including a downloadable file of the full database with additional variant details. denovo-db provides necessary information for researchers to compare their data to other individuals with the same phenotype and also to controls allowing for a better understanding of the biology of de novo variants and their contribution to disease.

}, keywords = {Computational Biology, Databases, Nucleic Acid, Genetic Association Studies, Genetic Variation, Germ-Line Mutation, Humans, Molecular Sequence Annotation, Polymorphism, Single Nucleotide, Web Browser}, issn = {1362-4962}, doi = {10.1093/nar/gkw865}, author = {Turner, Tychele N and Yi, Qian and Krumm, Niklas and Huddleston, John and Hoekzema, Kendra and F Stessman, Holly A and Doebley, Anna-Lisa and Bernier, Raphael A and Nickerson, Deborah A and Eichler, Evan E} } @article {372, title = {Homozygous and hemizygous CNV detection from exome sequencing data in a Mendelian disease cohort.}, journal = {Nucleic Acids Res}, volume = {45}, year = {2017}, month = {2017 02 28}, pages = {1633-1648}, abstract = {

We developed an algorithm, HMZDelFinder, that uses whole exome sequencing (WES) data to identify rare and intragenic homozygous and hemizygous (HMZ) deletions that may represent complete loss-of-function of the indicated gene. HMZDelFinder was applied to 4866 samples in the Baylor-Hopkins Center for Mendelian Genomics (BHCMG) cohort and detected 773 HMZ deletion calls (567 homozygous or 206 hemizygous) with an estimated sensitivity of 86.5\% (82\% for single-exonic and 88\% for multi-exonic calls) and precision of 78\% (53\% single-exonic and 96\% for multi-exonic calls). Out of 773 HMZDelFinder-detected deletion calls, 82 were subjected to array comparative genomic hybridization (aCGH) and/or breakpoint PCR and 64 were confirmed. These include 18 single-exon deletions out of which 8 were exclusively detected by HMZDelFinder and not by any of seven other CNV detection tools examined. Further investigation of the 64 validated deletion calls revealed at least 15 pathogenic HMZ deletions. Of those, 7 accounted for 17-50\% of pathogenic CNVs in different disease cohorts where 7.1-11\% of the molecular diagnosis solved rate was attributed to CNVs. In summary, we present an algorithm to detect rare, intragenic, single-exon deletion CNVs using WES data; this tool can be useful for disease gene discovery efforts and clinical WES analyses.

}, keywords = {Algorithms, Alternative Splicing, Cohort Studies, Computational Biology, Consanguinity, Datasets as Topic, DNA Copy Number Variations, Exome, Genetic Diseases, Inborn, Hemizygote, High-Throughput Nucleotide Sequencing, Homozygote, Humans, Inheritance Patterns, Models, Genetic, Pedigree, Reproducibility of Results, Sequence Deletion, Workflow}, issn = {1362-4962}, doi = {10.1093/nar/gkw1237}, author = {Gambin, Tomasz and Akdemir, Zeynep C and Yuan, Bo and Gu, Shen and Chiang, Theodore and Carvalho, Claudia M B and Shaw, Chad and Jhangiani, Shalini and Boone, Philip M and Eldomery, Mohammad K and Karaca, Ender and Bayram, Yavuz and Stray-Pedersen, Asbj{\o}rg and Muzny, Donna and Charng, Wu-Lin and Bahrambeigi, Vahid and Belmont, John W and Boerwinkle, Eric and Beaudet, Arthur L and Gibbs, Richard A and Lupski, James R} } @article {352, title = {Intensification: A Resource for Amplifying Population-Genetic Signals with Protein Repeats.}, journal = {J Mol Biol}, volume = {429}, year = {2017}, month = {2017 02 03}, pages = {435-445}, abstract = {

Large-scale genome sequencing holds great promise for the interpretation of protein structures through the discovery of many, rare functional variants in the human population. However, because protein-coding regions are under high selective constraints, these variants occur at low frequencies, such that there is often insufficient statistics for downstream calculations. To address this problem, we develop the Intensification approach, which uses the modular structure of repeat protein domains to amplify signals of selection from population genetics and traditional interspecies conservation. In particular, we are able to aggregate variants at the codon level to identify important positions in repeat domains that show strong conservation signals. This allows us to compare conservation over different evolutionary timescales. It also enables us to visualize population-genetic measures on protein structures. We make available the Intensification results as an online resource (http://intensification.gersteinlab.org) and illustrate the approach through a case study on the tetratricopeptide repeat.

}, keywords = {Amino Acid Sequence, Chromosome Mapping, Codon, Computational Biology, Databases, Genetic, Evolution, Molecular, Genetics, Population, Humans, Internet, Molecular Structure, Protein Conformation, Proteins, Signal Processing, Computer-Assisted, X-Rays}, issn = {1089-8638}, doi = {10.1016/j.jmb.2016.12.003}, author = {Chen, Jieming and Wang, Bo and Regan, Lynne and Gerstein, Mark} } @article {357, title = {Lessons learned from additional research analyses of unsolved clinical exome cases.}, journal = {Genome Med}, volume = {9}, year = {2017}, month = {2017 03 21}, pages = {26}, abstract = {

BACKGROUND: Given the rarity of most single-gene Mendelian disorders, concerted efforts of data exchange between clinical and scientific communities are critical to optimize molecular diagnosis and novel disease gene discovery.

METHODS: We designed and implemented protocols for the study of cases for which a plausible molecular diagnosis was not achieved in a clinical genomics diagnostic laboratory (i.e. unsolved clinical exomes). Such cases were recruited to a research laboratory for further analyses, in order to potentially: (1) accelerate novel disease gene discovery; (2) increase the molecular diagnostic yield of whole exome sequencing (WES); and (3) gain insight into the genetic mechanisms of disease. Pilot project data included 74 families, consisting mostly of parent-offspring trios. Analyses performed on a research basis employed both WES from additional family members and complementary bioinformatics approaches and protocols.

RESULTS: Analysis of all possible modes of Mendelian inheritance, focusing on both single nucleotide variants (SNV) and copy number variant (CNV) alleles, yielded a likely contributory variant in 36\% (27/74) of cases. If one includes candidate genes with variants identified within a single family, a potential contributory variant was identified in a total of ~51\% (38/74) of cases enrolled in this pilot study. The molecular diagnosis was achieved in 30/63 trios (47.6\%). Besides this, the analysis workflow yielded evidence for pathogenic variants in disease-associated genes in 4/6 singleton cases (66.6\%), 1/1 multiplex family involving three affected siblings, and 3/4 (75\%) quartet families. Both the analytical pipeline and the collaborative efforts between the diagnostic and research laboratories provided insights that allowed recent disease gene discoveries (PURA, TANGO2, EMC1, GNB5, ATAD3A, and MIPEP) and increased the number of novel genes, defined in this study as genes identified in more than one family (DHX30 and EBF3).

CONCLUSION: An efficient genomics pipeline in which clinical sequencing in a diagnostic laboratory is followed by the detailed reanalysis of unsolved cases in a research environment, supplemented with WES data from additional family members, and subject to adjuvant bioinformatics analyses including relaxed variant filtering parameters in informatics pipelines, can enhance the molecular diagnostic yield and provide mechanistic insights into Mendelian disorders. Implementing these approaches requires collaborative clinical molecular diagnostic and research efforts.

}, keywords = {Adenosine Triphosphatases, ATPases Associated with Diverse Cellular Activities, Computational Biology, DNA Copy Number Variations, DNA-Binding Proteins, Exome, Female, Genetic Diseases, Inborn, Genomics, GTP-Binding Protein beta Subunits, Humans, Male, Membrane Proteins, Metalloendopeptidases, Mitochondrial Proteins, Pilot Projects, Polymorphism, Single Nucleotide, Sequence Analysis, DNA, Transcription Factors}, issn = {1756-994X}, doi = {10.1186/s13073-017-0412-6}, author = {Eldomery, Mohammad K and Coban-Akdemir, Zeynep and Harel, Tamar and Rosenfeld, Jill A and Gambin, Tomasz and Stray-Pedersen, Asbj{\o}rg and K{\"u}ry, S{\'e}bastien and Mercier, Sandra and Lessel, Davor and Denecke, Jonas and Wiszniewski, Wojciech and Penney, Samantha and Liu, Pengfei and Bi, Weimin and Lalani, Seema R and Schaaf, Christian P and Wangler, Michael F and Bacino, Carlos A and Lewis, Richard Alan and Potocki, Lorraine and Graham, Brett H and Belmont, John W and Scaglia, Fernando and Orange, Jordan S and Jhangiani, Shalini N and Chiang, Theodore and Doddapaneni, Harsha and Hu, Jianhong and Muzny, Donna M and Xia, Fan and Beaudet, Arthur L and Boerwinkle, Eric and Eng, Christine M and Plon, Sharon E and Sutton, V Reid and Gibbs, Richard A and Posey, Jennifer E and Yang, Yaping and Lupski, James R} } @article {440, title = {Matchmaker Exchange.}, journal = {Curr Protoc Hum Genet}, volume = {95}, year = {2017}, month = {2017 10 18}, pages = {9.31.1-9.31.15}, abstract = {

In well over half of the individuals with rare disease who undergo clinical or research next-generation sequencing, the responsible gene cannot be determined. Some reasons for this relatively low yield include unappreciated phenotypic heterogeneity; locus heterogeneity; somatic and germline mosaicism; variants of uncertain functional significance; technically inaccessible areas of the genome; incorrect mode of inheritance investigated; and inadequate communication between clinicians and basic scientists with knowledge of particular genes, proteins, or biological systems. To facilitate such communication and improve the search for patients or model organisms with similar phenotypes and variants in specific candidate genes, we have developed the Matchmaker Exchange (MME). MME was created to establish a federated network connecting databases of genomic and phenotypic data using a common application programming interface (API). To date, seven databases can exchange data using the API (GeneMatcher, PhenomeCentral, DECIPHER, MyGene2, matchbox, Australian Genomics Health Alliance Patient Archive, and Monarch Initiative; the latter included for model organism matching). This article guides usage of the MME for rare disease gene discovery. {\textcopyright} 2017 by John Wiley \& Sons, Inc.

}, keywords = {Animals, Computational Biology, Databases, Genetic, Genetic Association Studies, Genetic Predisposition to Disease, Genomics, Humans, Rare Diseases, Software, Web Browser}, issn = {1934-8258}, doi = {10.1002/cphg.50}, author = {Sobreira, Nara L M and Arachchi, Harindra and Buske, Orion J and Chong, Jessica X and Hutton, Ben and Foreman, Julia and Schiettecatte, Fran{\c c}ois and Groza, Tudor and Jacobsen, Julius O B and Haendel, Melissa A and Boycott, Kym M and Hamosh, Ada and Rehm, Heidi L} } @article {22, title = {My46: a Web-based tool for self-guided management of genomic test results in research and clinical settings.}, journal = {Genet Med}, volume = {19}, year = {2017}, month = {2017 Apr}, pages = {467-475}, abstract = {

A major challenge to implementing precision medicine is the need for an efficient and cost-effective strategy for returning individual genomic test results that is easily scalable and can be incorporated into multiple models of clinical practice. My46 is a Web-based tool for managing the return of genetic results that was designed and developed to support a wide range of approaches to disclosing results, ranging from traditional face-to-face disclosure to self-guided models. My46 has five key functions: set and modify results-return preferences, return results, educate, manage the return of results, and assess the return of results. These key functions are supported by six distinct modules and a suite of features that enhance the user experience, ease site navigation, facilitate knowledge sharing, and enable results-return tracking. My46 is a potentially effective solution for returning results and supports current trends toward shared decision making between patients and providers and patient-driven health management.Genet Med 19 4, 467-475.

}, keywords = {Biomedical Research, Computational Biology, Decision Making, Genomics, Humans, Internet, Medical Informatics, Patient Access to Records, Precision Medicine}, issn = {1530-0366}, doi = {10.1038/gim.2016.133}, author = {Tabor, Holly K and Jamal, Seema M and Yu, Joon-Ho and Crouch, Julia M and Shankar, Aditi G and Dent, Karin M and Anderson, Nick and Miller, Damon A and Futral, Brett T and Bamshad, Michael J} } @article {535, title = {Localized structural frustration for evaluating the impact of sequence variants.}, journal = {Nucleic Acids Res}, volume = {44}, year = {2016}, month = {2016 12 01}, pages = {10062-10073}, abstract = {

Population-scale sequencing is increasingly uncovering large numbers of rare single-nucleotide variants (SNVs) in coding regions of the genome. The rarity of these variants makes it challenging to evaluate their deleteriousness with conventional phenotype-genotype associations. Protein structures provide a way of addressing this challenge. Previous efforts have focused on globally quantifying the impact of SNVs on protein stability. However, local perturbations may severely impact protein functionality without strongly disrupting global stability (e.g. in relation to catalysis or allostery). Here, we describe a workflow in which localized frustration, quantifying unfavorable local interactions, is employed as a metric to investigate such effects. Using this workflow on the Protein Databank, we find that frustration produces many immediately intuitive results: for instance, disease-related SNVs create stronger changes in localized frustration than non-disease related variants, and rare SNVs tend to disrupt local interactions to a larger extent than common variants. Less obviously, we observe that somatic SNVs associated with oncogenes and tumor suppressor genes (TSGs) induce very different changes in frustration. In particular, those associated with TSGs change the frustration more in the core than the surface (by introducing loss-of-function events), whereas those associated with oncogenes manifest the opposite pattern, creating gain-of-function events.

}, keywords = {Computational Biology, Databases, Protein, Evolution, Molecular, Genes, Tumor Suppressor, Genetic Variation, Humans, Neoplasms, Oncogenes, Polymorphism, Single Nucleotide, Proteins, Workflow}, issn = {1362-4962}, doi = {10.1093/nar/gkw927}, author = {Kumar, Sushant and Clarke, Declan and Gerstein, Mark} } @article {96, title = {A patient with a novel homozygous missense mutation in FTO and concomitant nonsense mutation in CETP.}, journal = {J Hum Genet}, volume = {61}, year = {2016}, month = {2016 May}, pages = {395-403}, abstract = {

The fat mass and obesity associated (FTO) gene has previously been associated with a variety of diseases and conditions, notably obesity, acute coronary syndrome and metabolic syndrome. Reports describing mutations in FTO as well as in FTO animal models have further demonstrated a role for FTO in the development of the brain and other organs. Here, we describe a patient born of consanguineous union who presented with microcephaly, developmental delay, behavioral abnormalities, dysmorphic facial features, hypotonia and other various phenotypic abnormalities. Whole-exome sequencing revealed a novel homozygous missense mutation in FTO and a nonsense mutation in the cholesteryl ester transfer protein (CETP). Exome copy number variation analysis revealed no disease-causing large duplications or deletions within coding regions. Patient{\textquoteright}s, her parents{\textquoteright} and non-related control{\textquoteright} fibroblasts were analyzed for morphologic defects, abnormal proliferation, apoptosis and transcriptome profile. We have shown that FTO is located in the nucleus of cells from each tested sample. Western blot analysis demonstrated no changes in patient FTO. Quantitative (qPCR) analysis revealed slightly decreased levels of FTO expression in patient cells compared with controls. No morphological or proliferation differences between the patient and control fibroblasts were observed. There is still much to be learned about the molecular mechanisms by which mutations in FTO contribute to such severe phenotypes.

}, keywords = {Alpha-Ketoglutarate-Dependent Dioxygenase FTO, Apoptosis, Biopsy, Child, Preschool, Cholesterol Ester Transfer Proteins, Computational Biology, Consanguinity, DNA Copy Number Variations, DNA Mutational Analysis, Exome, Female, Gene Expression, Gene Expression Profiling, Genetic Association Studies, Genotype, High-Throughput Nucleotide Sequencing, Homozygote, Humans, Mutation, Missense, Phenotype, Transcriptome}, issn = {1435-232X}, doi = {10.1038/jhg.2015.160}, author = {{\c C}a{\u g}layan, Ahmet O and Tuysuz, Beyhan and Co{\c s}kun, S{\"u}leyman and Quon, Jennifer and Harmanc{\i}, Akdes S and Baranoski, Jacob F and Baran, Bur{\c c}in and Erson-Omay, E Zeynep and Henegariu, Octavian and Mane, Shrikant M and Bilguvar, Kaya and Yasuno, Katsuhito and G{\"u}nel, Murat} } @article {88, title = {A visual and curatorial approach to clinical variant prioritization and disease gene discovery in genome-wide diagnostics.}, journal = {Genome Med}, volume = {8}, year = {2016}, month = {2016 Feb 02}, pages = {13}, abstract = {

BACKGROUND: Genome-wide data are increasingly important in the clinical evaluation of human disease. However, the large number of variants observed in individual patients challenges the efficiency and accuracy of diagnostic review. Recent work has shown that systematic integration of clinical phenotype data with genotype information can improve diagnostic workflows and prioritization of filtered rare variants. We have developed visually interactive, analytically transparent analysis software that leverages existing disease catalogs, such as the Online Mendelian Inheritance in Man database (OMIM) and the Human Phenotype Ontology (HPO), to integrate patient phenotype and variant data into ranked diagnostic alternatives.

METHODS: Our tool, "OMIM Explorer" ( http://www.omimexplorer.com ), extends the biomedical application of semantic similarity methods beyond those reported in previous studies. The tool also provides a simple interface for translating free-text clinical notes into HPO terms, enabling clinical providers and geneticists to contribute phenotypes to the diagnostic process. The visual approach uses semantic similarity with multidimensional scaling to collapse high-dimensional phenotype and genotype data from an individual into a graphical format that contextualizes the patient within a low-dimensional disease map. The map proposes a differential diagnosis and algorithmically suggests potential alternatives for phenotype queries--in essence, generating a computationally assisted differential diagnosis informed by the individual{\textquoteright}s personal genome. Visual interactivity allows the user to filter and update variant rankings by interacting with intermediate results. The tool also implements an adaptive approach for disease gene discovery based on patient phenotypes.

RESULTS: We retrospectively analyzed pilot cohort data from the Baylor Miraca Genetics Laboratory, demonstrating performance of the tool and workflow in the re-analysis of clinical exomes. Our tool assigned to clinically reported variants a median rank of 2, placing causal variants in the top 1 \% of filtered candidates across the 47 cohort cases with reported molecular diagnoses of exome variants in OMIM Morbidmap genes. Our tool outperformed Phen-Gen, eXtasy, PhenIX, PHIVE, and hiPHIVE in the prioritization of these clinically reported variants.

CONCLUSIONS: Our integrative paradigm can improve efficiency and, potentially, the quality of genomic medicine by more effectively utilizing available phenotype information, catalog data, and genomic knowledge.

}, keywords = {Computational Biology, Diagnosis, Differential, Disease, Genetic Variation, Genome-Wide Association Study, Humans, Phenotype, Pilot Projects, Retrospective Studies, Software}, issn = {1756-994X}, doi = {10.1186/s13073-016-0261-8}, author = {James, Regis A and Campbell, Ian M and Chen, Edward S and Boone, Philip M and Rao, Mitchell A and Bainbridge, Matthew N and Lupski, James R and Yang, Yaping and Eng, Christine M and Posey, Jennifer E and Shaw, Chad A} } @article {177, title = {Assessing structural variation in a personal genome-towards a human reference diploid genome.}, journal = {BMC Genomics}, volume = {16}, year = {2015}, month = {2015 Apr 11}, pages = {286}, abstract = {

BACKGROUND: Characterizing large genomic variants is essential to expanding the research and clinical applications of genome sequencing. While multiple data types and methods are available to detect these structural variants (SVs), they remain less characterized than smaller variants because of SV diversity, complexity, and size. These challenges are exacerbated by the experimental and computational demands of SV analysis. Here, we characterize the SV content of a personal genome with Parliament, a publicly available consensus SV-calling infrastructure that merges multiple data types and SV detection methods.

RESULTS: We demonstrate Parliament{\textquoteright}s efficacy via integrated analyses of data from whole-genome array comparative genomic hybridization, short-read next-generation sequencing, long-read (Pacific BioSciences RSII), long-insert (Illumina Nextera), and whole-genome architecture (BioNano Irys) data from the personal genome of a single subject (HS1011). From this genome, Parliament identified 31,007 genomic loci between 100 bp and 1 Mbp that are inconsistent with the hg19 reference assembly. Of these loci, 9,777 are supported as putative SVs by hybrid local assembly, long-read PacBio data, or multi-source heuristics. These SVs span 59 Mbp of the reference genome (1.8\%) and include 3,801 events identified only with long-read data. The HS1011 data and complete Parliament infrastructure, including a BAM-to-SV workflow, are available on the cloud-based service DNAnexus.

CONCLUSIONS: HS1011 SV analysis reveals the limits and advantages of multiple sequencing technologies, specifically the impact of long-read SV discovery. With the full Parliament infrastructure, the HS1011 data constitute a public resource for novel SV discovery, software calibration, and personal genome structural variation analysis.

}, keywords = {Computational Biology, Databases, Genetic, Diploidy, Genome, Human, Genomic Structural Variation, Humans, Sequence Analysis, DNA, Software}, issn = {1471-2164}, doi = {10.1186/s12864-015-1479-3}, author = {English, Adam C and Salerno, William J and Hampton, Oliver A and Gonzaga-Jauregui, Claudia and Ambreth, Shruthi and Ritter, Deborah I and Beck, Christine R and Davis, Caleb F and Dahdouli, Mahmoud and Ma, Singer and Carroll, Andrew and Veeraraghavan, Narayanan and Bruestle, Jeremy and Drees, Becky and Hastie, Alex and Lam, Ernest T and White, Simon and Mishra, Pamela and Wang, Min and Han, Yi and Zhang, Feng and Stankiewicz, Pawel and Wheeler, David A and Reid, Jeffrey G and Muzny, Donna M and Rogers, Jeffrey and Sabo, Aniko and Worley, Kim C and Lupski, James R and Boerwinkle, Eric and Gibbs, Richard A} } @article {139, title = {GeneMatcher: a matching tool for connecting investigators with an interest in the same gene.}, journal = {Hum Mutat}, volume = {36}, year = {2015}, month = {2015 Oct}, pages = {928-30}, abstract = {

Here, we describe an overview and update on GeneMatcher (http://www.genematcher.org), a freely accessible Web-based tool developed as part of the Baylor-Hopkins Center for Mendelian Genomics. We created GeneMatcher with the goal of identifying additional individuals with rare phenotypes who had variants in the same candidate disease gene. We also wanted to facilitate connections to basic scientists working on orthologous genes in model systems with the goal of connecting their work to human Mendelian phenotypes. Meeting these goals will enhance the identification of novel Mendelian genes. Launched in September, 2013, GeneMatcher now has 2,178 candidate genes from 486 submitters spread across 38 countries entered in the database (June 1, 2015). GeneMatcher is also part of the Matchmaker Exchange (http://matchmakerexchange.org/) with an Application Programing Interface enabling submitters to query other databases of genetic variants and phenotypes without having to create accounts and data entries in multiple systems.

}, keywords = {Computational Biology, Databases, Genetic, Genetic Predisposition to Disease, Genetic Variation, Humans, Information Dissemination, Phenotype, Rare Diseases, Software, User-Computer Interface, Web Browser}, issn = {1098-1004}, doi = {10.1002/humu.22844}, author = {Sobreira, Nara and Schiettecatte, Fran{\c c}ois and Valle, David and Hamosh, Ada} } @article {136, title = {The Matchmaker Exchange API: automating patient matching through the exchange of structured phenotypic and genotypic profiles.}, journal = {Hum Mutat}, volume = {36}, year = {2015}, month = {2015 Oct}, pages = {922-7}, abstract = {

Despite the increasing prevalence of clinical sequencing, the difficulty of identifying additional affected families is a key obstacle to solving many rare diseases. There may only be a handful of similar patients worldwide, and their data may be stored in diverse clinical and research databases. Computational methods are necessary to enable finding similar patients across the growing number of patient repositories and registries. We present the Matchmaker Exchange Application Programming Interface (MME API), a protocol and data format for exchanging phenotype and genotype profiles to enable matchmaking among patient databases, facilitate the identification of additional cohorts, and increase the rate with which rare diseases can be researched and diagnosed. We designed the API to be straightforward and flexible in order to simplify its adoption on a large number of data types and workflows. We also provide a public test data set, curated from the literature, to facilitate implementation of the API and development of new matching algorithms. The initial version of the API has been successfully implemented by three members of the Matchmaker Exchange and was immediately able to reproduce previously identified matches and generate several new leads currently being validated. The API is available at https://github.com/ga4gh/mme-apis.

}, keywords = {Algorithms, Computational Biology, Databases, Genetic, Genetic Predisposition to Disease, Genotype, Humans, Information Dissemination, Phenotype, Rare Diseases, Web Browser}, issn = {1098-1004}, doi = {10.1002/humu.22850}, author = {Buske, Orion J and Schiettecatte, Fran{\c c}ois and Hutton, Benjamin and Dumitriu, Sergiu and Misyura, Andriy and Huang, Lijia and Hartley, Taila and Girdea, Marta and Sobreira, Nara and Mungall, Chris and Brudno, Michael} } @article {193, title = {New tools for Mendelian disease gene identification: PhenoDB variant analysis module; and GeneMatcher, a web-based tool for linking investigators with an interest in the same gene.}, journal = {Hum Mutat}, volume = {36}, year = {2015}, month = {2015 Apr}, pages = {425-31}, abstract = {

Identifying the causative variant from among the thousands identified by whole-exome sequencing or whole-genome sequencing is a formidable challenge. To make this process as efficient and flexible as possible, we have developed a Variant Analysis Module coupled to our previously described Web-based phenotype intake tool, PhenoDB (http://researchphenodb.net and http://phenodb.org). When a small number of candidate-causative variants have been identified in a study of a particular patient or family, a second, more difficult challenge becomes proof of causality for any given variant. One approach to this problem is to find other cases with a similar phenotype and mutations in the same candidate gene. Alternatively, it may be possible to develop biological evidence for causality, an approach that is assisted by making connections to basic scientists studying the gene of interest, often in the setting of a model organism. Both of these strategies benefit from an open access, online site where individual clinicians and investigators could post genes of interest. To this end, we developed GeneMatcher (http://genematcher.org), a freely accessible Website that enables connections between clinicians and researchers across the world who share an interest in the same gene(s).

}, keywords = {Computational Biology, Databases, Genetic, Genetic Association Studies, Genetic Diseases, Inborn, Genetic Variation, Genomics, Internet, Phenotype, Software, Web Browser}, issn = {1098-1004}, doi = {10.1002/humu.22769}, author = {Sobreira, Nara and Schiettecatte, Fran{\c c}ois and Boehm, Corinne and Valle, David and Hamosh, Ada} } @article {277, title = {DIAMUND: direct comparison of genomes to detect mutations.}, journal = {Hum Mutat}, volume = {35}, year = {2014}, month = {2014 Mar}, pages = {283-8}, abstract = {

DNA sequencing has become a powerful method to discover the genetic basis of disease. Standard, widely used protocols for analysis usually begin by comparing each individual to the human reference genome. When applied to a set of related individuals, this approach reveals millions of differences, most of which are shared among the individuals and unrelated to the disease being investigated. We have developed a novel algorithm for variant detection, one that compares DNA sequences directly to one another, without aligning them to the reference genome. When used to find de novo mutations in exome sequences from family trios, or to compare normal and diseased samples from the same individual, the new method, direct alignment for mutation discovery (DIAMUND), produces a dramatically smaller list of candidate mutations than previous methods, without losing sensitivity to detect the true cause of a genetic disease. We demonstrate our results on several example cases, including two family trios in which it correctly found the disease-causing variant while excluding thousands of harmless variants that standard methods had identified.

}, keywords = {Algorithms, Computational Biology, Exome, Genetic Loci, Genetic Variation, Genome, Human, Humans, Mutation, Sequence Alignment, Sequence Analysis, DNA, Software}, issn = {1098-1004}, doi = {10.1002/humu.22503}, author = {Salzberg, Steven L and Pertea, Mihaela and Fahrner, Jill A and Sobreira, Nara} } @article {299, title = {Mutations in KARS, encoding lysyl-tRNA synthetase, cause autosomal-recessive nonsyndromic hearing impairment DFNB89.}, journal = {Am J Hum Genet}, volume = {93}, year = {2013}, month = {2013 Jul 11}, pages = {132-40}, abstract = {

Previously, DFNB89, a locus associated with autosomal-recessive nonsyndromic hearing impairment (ARNSHI), was mapped to chromosomal region 16q21-q23.2 in three unrelated, consanguineous Pakistani families. Through whole-exome sequencing of a hearing-impaired individual from each family, missense mutations were identified at highly conserved residues of lysyl-tRNA synthetase (KARS): the c.1129G>A (p.Asp377Asn) variant was found in one family, and the c.517T>C (p.Tyr173His) variant was found in the other two families. Both variants were predicted to be damaging by multiple bioinformatics tools. The two variants both segregated with the nonsyndromic-hearing-impairment phenotype within the three families, and neither mutation was identified in ethnically matched controls or within variant databases. Individuals homozygous for KARS mutations had symmetric, severe hearing impairment across all frequencies but did not show evidence of auditory or limb neuropathy. It has been demonstrated that KARS is expressed in hair cells of zebrafish, chickens, and mice. Moreover, KARS has strong localization to the spiral ligament region of the cochlea, as well as to Deiters{\textquoteright} cells, the sulcus epithelium, the basilar membrane, and the surface of the spiral limbus. It is hypothesized that KARS variants affect aminoacylation in inner-ear cells by interfering with binding activity to tRNA or p38 and with tetramer formation. The identification of rare KARS variants in ARNSHI-affected families defines a gene that is associated with ARNSHI.

}, keywords = {Animals, Case-Control Studies, Chickens, Cochlea, Computational Biology, Consanguinity, Female, Genetic Linkage, Hair Cells, Auditory, Inner, Haplotypes, Hearing Loss, Homozygote, Humans, Lysine-tRNA Ligase, Male, Mice, Mutation, Missense, Pedigree, Transfer RNA Aminoacylation, Zebrafish}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2013.05.018}, author = {Santos-Cortez, Regie Lyn P and Lee, Kwanghyuk and Azeem, Zahid and Antonellis, Patrick J and Pollock, Lana M and Khan, Saadullah and Andrade-Elizondo, Paula B and Chiu, Ilene and Adams, Mark D and Basit, Sulman and Smith, Joshua D and Nickerson, Deborah A and McDermott, Brian M and Ahmad, Wasim and Leal, Suzanne M} }