Search

Summary: Background: Belgian horses are commonly affected with ocular squamous cell carcinoma (SCC), the most common cancer of the equine eye. A missense mutation in damage‐specific DNA binding protein 2 (DDB2 c.1013C>T, p.Thr338Met) has been established as a recessive genetic risk factor for ocular SCC in the Haflinger breed. A sample of Belgian horses with unknown SCC phenotype was shown to possess this variant at a similar frequency to the Haflinger breed. Retrospective studies indicate that chestnut coat colour may predispose to the development of SCC. Objectives: To determine if DDB2 c.1013C>T is a risk factor for ocular SCC in a strictly phenotyped sample of Belgian horses. To investigate associations between coat colour loci genotypes and ocular SCC. Study design: Retrospective and prospective case identification, genetic investigation. Methods: Genomic DNA was isolated from blood, hair or formalin‐fixed paraffin‐embedded tissue from 25 Belgian horses with histologically confirmed ocular SCC and 18 unaffected Belgian horses. Association testing of 34 single nucleotide variants from 11 genomic loci and genotyping for DDB2 c.1013C>T and coat colour alleles were performed. Exons of DDB2 were sequenced in four cases and two controls. Associations were analysed by Chi‐square or Fisher's exact tests and relative risk was calculated. Results: Homozygosity for DDB2 c.1013C>T was significantly associated with ocular SCC (P = 7.4 × 10−7). Seventy‐six per cent of affected horses were homozygous for the variant. Relative risk for homozygous horses developing SCC was 4.0 (P = 1.0 × 10−4). Sequencing DDB2 did not identify a variant more concordant with disease phenotype. An association between disease and coat colour loci was not identified. Main limitations: Phenotyping was determined at a single timepoint. Each included horse genotyped as chestnut, so association with this MC1R variant could not be investigated. Conclusions: A missense variant, DDB2 c.1013C>T, p.Thr338Met, is a risk factor for ocular SCC in Belgian horses. A genetic risk test is commercially available. [ABSTRACT FROM AUTHOR]

BACKGROUND: Genetic variants that cause rare disorders may remain elusive even after expansive testing, such as exome sequencing. The diagnostic yield of genome sequencing, particularly after a negative evaluation, remains poorly defined. METHODS: We sequenced and analyzed the genomes of families with diverse phenotypes who were suspected to have a rare monogenic disease and for whom genetic testing had not revealed a diagnosis, as well as the genomes of a replication cohort at an independent clinical center. RESULTS: We sequenced the genomes of 822 families (744 in the initial cohort and 78 in the replication cohort) and made a molecular diagnosis in 218 of 744 families (29.3%). Of the 218 families, 61 (28.0%) -- 8.2% of families in the initial cohort -- had variants that required genome sequencing for identification, including coding variants, intronic variants, small structural variants, copy-neutral inversions, complex rearrangements, and tandem repeat expansions. Most families in which a molecular diagnosis was made after previous nondiagnostic exome sequencing (63.5%) had variants that could be detected by reanalysis of the exome-sequence data (53.4%) or by additional analytic methods, such as copy-number variant calling, to exome-sequence data (10.8%). We obtained similar results in the replication cohort: in 33% of the families in which a molecular diagnosis was made, or 8% of the cohort, genome sequencing was required, which showed the applicability of these findings to both research and clinical environments. CONCLUSIONS: The diagnostic yield of genome sequencing in a large, diverse research cohort and in a small clinical cohort of persons who had previously undergone genetic testing was approximately 8% and included several types of pathogenic variation that had not previously been detected by means of exome sequencing or other techniques. (Funded by the National Human Genome Research Institute and others.). [ABSTRACT FROM AUTHOR]

Background: PLA2G6-associated neurodegeneration (PLAN) comprises three diseases with overlapping features: infantile neuroaxonal dystrophy (INAD), atypical neuroaxonal dystrophy (atypical NAD), and PLA2G6-related dystonia-parkinsonism. INAD is an early onset disease characterized by progressive loss of vision, muscular control, and mental skills. The prevalence of PLA2G6-associated diseases has not been previously calculated., Methods: To provide the most accurate prevalence estimate, we utilized two independent approaches: database-based approach which included collecting variants from ClinVar, Human Gene Mutation Database (HGMD) and high confidence predicted loss-of-function (pLoF) from gnomAD (Rare Genomes Project Genetic Prevalence Estimator; GeniE), and literature-based approach which gathered variants through Mastermind Genomic Search Engine (Genomenon, Inc). Genetic prevalence of PLAN was calculated based on allele frequencies from gnomAD, assuming Hardy-Weinberg equilibrium., Results: In the PLA2G6 gene, our analysis found 122 pathogenic, 82 VUS, and 15 variants with conflicting interpretations (pathogenic vs VUS) between two approaches. Allele frequency was available for 58 pathogenic, 42 VUS, and 15 conflicting variants in gnomAD database. If pathogenic and/or conflicting variants are included, the overall genetic prevalence was estimated to be between 1 in 987,267 to 1 in 1,570,079 pregnancies, with the highest genetic prevalence in African/African-American (1 in 421,960 to 1 in 365,197) and East-Asian (1 in 683,978 to 1 in 190,771) populations., Conclusion: Our estimates highlight the significant underdiagnosis of PLA2G6-associated neurodegeneration and underscores the need for increased awareness and diagnostic efforts. Furthermore, our study revealed a higher carrier frequency of PLA2G6 variants in African and Asian populations, stressing the importance of expanded genetic sequencing in non-European populations to ensure accurate and comprehensive diagnosis. Future research should focus on confirming our findings and implementing expanded sequencing strategies to facilitate maximal and accurate diagnosis, particularly in non-European populations., (© 2024. The Author(s).)

Purpose: To identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs)., Methods: We coupled phenotyping with exome or genome sequencing of 467 probands (550 affected and 1108 total individuals) with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations., Results: Analyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43/467 probands (9.2%), and prioritized variants of uncertain significance in 70/467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (e.g., MYH10, KIF21B, TGFBR2, TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (e.g., CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (e.g., TUBA4A, KIF5C, CTNNA1, KLB, FGF21), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses., Conclusion: This study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies., (Copyright © 2024. Published by Elsevier Inc.)

Incomplete penetrance, or absence of disease phenotype in an individual with a disease-associated variant, is a major challenge in variant interpretation. Studying individuals with apparent incomplete penetrance can shed light on underlying drivers of altered phenotype penetrance. Here, we investigate clinically relevant variants from ClinVar in 807,162 individuals from the Genome Aggregation Database (gnomAD), demonstrating improved representation in gnomAD version 4. We then conduct a comprehensive case-by-case assessment of 734 predicted loss of function variants (pLoF) in 77 genes associated with severe, early-onset, highly penetrant haploinsufficient disease. We identified explanations for the presumed lack of disease manifestation in 701 of the variants (95%). Individuals with unexplained lack of disease manifestation in this set of disorders rarely occur, underscoring the need and power of deep case-by-case assessment presented here to minimize false assignments of disease risk, particularly in unaffected individuals with higher rates of secondary properties that result in rescue., Competing Interests: COMPETING INTERESTS A.O-D.L. is on the scientific advisory board for Congenica, receives research funding in the form of reagents from Pacific Biosciences, and is a paid advisor to Addition Therapeutics and former advisor to Tome Biosciences and Ono Pharma USA. D.G.M. is a paid adviser to GlaxoSmithKline, Insitro, and Overtone Therapeutics, and receives research funding from Microsoft Corporation. H.L.R. receives research funding from Microsoft. T.L. is an advisor and has equity in Variant Bio.

Copy number variants (CNVs) are significant contributors to the pathogenicity of rare genetic diseases and, with new innovative methods, can now reliably be identified from exome sequencing. Challenges still remain in accurate classification of CNV pathogenicity. CNV calling using GATK-gCNV was performed on exomes from a cohort of 6,633 families (15,759 individuals) with heterogeneous phenotypes and variable prior genetic testing collected at the Broad Institute Center for Mendelian Genomics of the Genomics Research to Elucidate the Genetics of Rare Diseases consortium and analyzed using the seqr platform. The addition of CNV detection to exome analysis identified causal CNVs for 171 families (2.6%). The estimated sizes of CNVs ranged from 293 bp to 80 Mb. The causal CNVs consisted of 140 deletions, 15 duplications, 3 suspected complex structural variants (SVs), 3 insertions, and 10 complex SVs, the latter two groups being identified by orthogonal confirmation methods. To classify CNV variant pathogenicity, we used the 2020 American College of Medical Genetics and Genomics/ClinGen CNV interpretation standards and developed additional criteria to evaluate allelic and functional data as well as variants on the X chromosome to further advance the framework. We interpreted 151 CNVs as likely pathogenic/pathogenic and 20 CNVs as high-interest variants of uncertain significance. Calling CNVs from existing exome data increases the diagnostic yield for individuals undiagnosed after standard testing approaches, providing a higher-resolution alternative to arrays at a fraction of the cost of genome sequencing. Our improvements to the classification approach advances the systematic framework to assess the pathogenicity of CNVs., Competing Interests: Declaration of interests H.L.R. and M.E.T. have received support from Illumina and Microsoft to support rare disease gene discovery and diagnosis. A.O’D.-L. is on the scientific advisory board for Congenica Inc. D.G.M. is a paid advisor to GlaxoSmithKline, Insitro, Variant Bio, and Overtone Therapeutics and has received research support from AbbVie, Astellas, Biogen, BioMarin, Eisai, Google, Merck, Microsoft, Pfizer, and Sanofi-Genzyme. C.A.W. is a paid advisor to Maze Therapeutics. M.E.T. has also received reagents and/or research support from Levo Therapeutics, Pacific Biosciences, Oxford Nanopore, and Ionis Pharmaceuticals., (Copyright © 2024 American Society of Human Genetics. All rights reserved.)

Background: A major obstacle faced by families with rare diseases is obtaining a genetic diagnosis. The average "diagnostic odyssey" lasts over five years and causal variants are identified in under 50%, even when capturing variants genome-wide. To aid in the interpretation and prioritization of the vast number of variants detected, computational methods are proliferating. Knowing which tools are most effective remains unclear. To evaluate the performance of computational methods, and to encourage innovation in method development, we designed a Critical Assessment of Genome Interpretation (CAGI) community challenge to place variant prioritization models head-to-head in a real-life clinical diagnostic setting., Methods: We utilized genome sequencing (GS) data from families sequenced in the Rare Genomes Project (RGP), a direct-to-participant research study on the utility of GS for rare disease diagnosis and gene discovery. Challenge predictors were provided with a dataset of variant calls and phenotype terms from 175 RGP individuals (65 families), including 35 solved training set families with causal variants specified, and 30 unlabeled test set families (14 solved, 16 unsolved). We tasked teams to identify causal variants in as many families as possible. Predictors submitted variant predictions with estimated probability of causal relationship (EPCR) values. Model performance was determined by two metrics, a weighted score based on the rank position of causal variants, and the maximum F-measure, based on precision and recall of causal variants across all EPCR values., Results: Sixteen teams submitted predictions from 52 models, some with manual review incorporated. Top performers recalled causal variants in up to 13 of 14 solved families within the top 5 ranked variants. Newly discovered diagnostic variants were returned to two previously unsolved families following confirmatory RNA sequencing, and two novel disease gene candidates were entered into Matchmaker Exchange. In one example, RNA sequencing demonstrated aberrant splicing due to a deep intronic indel in ASNS, identified in trans with a frameshift variant in an unsolved proband with phenotypes consistent with asparagine synthetase deficiency., Conclusions: Model methodology and performance was highly variable. Models weighing call quality, allele frequency, predicted deleteriousness, segregation, and phenotype were effective in identifying causal variants, and models open to phenotype expansion and non-coding variants were able to capture more difficult diagnoses and discover new diagnoses. Overall, computational models can significantly aid variant prioritization. For use in diagnostics, detailed review and conservative assessment of prioritized variants against established criteria is needed., (© 2024. The Author(s).)

Heterozygous missense variants and in-frame indels in SMC3 are a cause of Cornelia de Lange syndrome (CdLS), marked by intellectual disability, growth deficiency, and dysmorphism, via an apparent dominant-negative mechanism. However, the spectrum of manifestations associated with SMC3 loss-of-function variants has not been reported, leading to hypotheses of alternative phenotypes or even developmental lethality. We used matchmaking servers, patient registries, and other resources to identify individuals with heterozygous, predicted loss-of-function (pLoF) variants in SMC3, and analyzed population databases to characterize mutational intolerance in this gene. Here, we show that SMC3 behaves as an archetypal haploinsufficient gene: it is highly constrained against pLoF variants, strongly depleted for missense variants, and pLoF variants are associated with a range of developmental phenotypes. Among 14 individuals with SMC3 pLoF variants, phenotypes were variable but coalesced on low growth parameters, developmental delay/intellectual disability, and dysmorphism, reminiscent of atypical CdLS. Comparisons to individuals with SMC3 missense/in-frame indel variants demonstrated an overall milder presentation in pLoF carriers. Furthermore, several individuals harboring pLoF variants in SMC3 were nonpenetrant for growth, developmental, and/or dysmorphic features, and some had alternative symptomatologies with rational biological links to SMC3. Analyses of tumor and model system transcriptomic data and epigenetic data in a subset of cases suggest that SMC3 pLoF variants reduce SMC3 expression but do not strongly support clustering with functional genomic signatures of typical CdLS. Our finding of substantial population-scale LoF intolerance in concert with variable growth and developmental features in subjects with SMC3 pLoF variants expands the scope of cohesinopathies, informs on their allelic architecture, and suggests the existence of additional clearly LoF-constrained genes whose disease links will be confirmed only by multilayered genomic data paired with careful phenotyping., Competing Interests: Declaration of interests M.E.T. is supported by research funding and/or reagents from Illumina, Microsoft, Ionis Therapeutics, and Levo Therapeutics. M.P.N., J.H., and J.J. are employees of GeneDx., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Recessive diseases arise when both copies of a gene are impacted by a damaging genetic variant. When a patient carries two potentially causal variants in a gene, accurate diagnosis requires determining that these variants occur on different copies of the chromosome (that is, are in trans) rather than on the same copy (that is, in cis). However, current approaches for determining phase, beyond parental testing, are limited in clinical settings. Here we developed a strategy for inferring phase for rare variant pairs within genes, leveraging genotypes observed in the Genome Aggregation Database (v2, n = 125,748 exomes). Our approach estimates phase with 96% accuracy, both in trio data and in patients with Mendelian conditions and presumed causal compound heterozygous variants. We provide a public resource of phasing estimates for coding variants and counts per gene of rare variants in trans that can aid interpretation of rare co-occurring variants in the context of recessive disease., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Genetic association studies have made significant contributions to our understanding of the etiology of neurodevelopmental disorders (NDDs). However, these studies rarely focused on the African continent. The NeuroDev Project aims to address this diversity gap through detailed phenotypic and genetic characterization of children with NDDs from Kenya and South Africa. We present results from NeuroDev's first year of data collection, including phenotype data from 206 cases and clinical genetic analyses of 99 parent-child trios. Most cases met criteria for global developmental delay/intellectual disability (GDD/ID, 80.3%). Approximately half of the children with GDD/ID also met criteria for autism. Analysis of exome-sequencing data identified a pathogenic or likely pathogenic variant in 13 (17%) of the 75 cases from South Africa and 9 (38%) of the 24 cases from Kenya. Data from the trio pilot are publicly available, and the NeuroDev Project will continue to develop resources for the global genetics community., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Predicted loss of function (pLoF) variants are often highly deleterious and play an important role in disease biology, but many pLoF variants may not result in loss of function (LoF). Here we present a framework that advances interpretation of pLoF variants in research and clinical settings by considering three categories of LoF evasion: (1) predicted rescue by secondary sequence properties, (2) uncertain biological relevance, and (3) potential technical artifacts. We also provide recommendations on adjustments to ACMG/AMP guidelines' PVS1 criterion. Applying this framework to all high-confidence pLoF variants in 22 genes associated with autosomal-recessive disease from the Genome Aggregation Database (gnomAD v.2.1.1) revealed predicted LoF evasion or potential artifacts in 27.3% (304/1,113) of variants. The major reasons were location in the last exon, in a homopolymer repeat, in a low proportion expressed across transcripts (pext) scored region, or the presence of cryptic in-frame splice rescues. Variants predicted to evade LoF or to be potential artifacts were enriched for ClinVar benign variants. PVS1 was downgraded in 99.4% (162/163) of pLoF variants predicted as likely not LoF/not LoF, with 17.2% (28/163) downgraded as a result of our framework, adding to previous guidelines. Variant pathogenicity was affected (mostly from likely pathogenic to VUS) in 20 (71.4%) of these 28 variants. This framework guides assessment of pLoF variants beyond standard annotation pipelines and substantially reduces false positive rates, which is key to ensure accurate LoF variant prediction in both a research and clinical setting., Competing Interests: Declaration of interests E.E. is a current employee of Clalit Research Institute and a graduate student at the Tel Aviv University Faculty of Medicine; her partner is an employee of GSK. K.J.K. is a paid advisor of Tome Biosciences and Vor Biosciences and is a member of the scientific advisory board of Nurture Genomics. S.M.H. is an employee of Ambry Genetics. D.G.M. is a paid advisor to GlaxoSmithKline (GSK), Insitro, Variant Bio, and Overtone Therapeutics and has received research support from AbbVie, Astellas, Biogen, BioMarin, Eisai, Google, Merck, Microsoft, Pfizer, and Sanofi-Genzyme. H.L.R. has received support from Illumina and Microsoft to support rare disease gene discovery and diagnosis. A.O’D.L. is a paid advisor of Tome Biosciences and member of the scientific advisory board for Congenica and the Simons Foundation SPARK for Autism study., (Copyright © 2023 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Personalized genome sequencing has revealed millions of genetic differences between individuals, but our understanding of their clinical relevance remains largely incomplete. To systematically decipher the effects of human genetic variants, we obtained whole-genome sequencing data for 809 individuals from 233 primate species and identified 4.3 million common protein-altering variants with orthologs in humans. We show that these variants can be inferred to have nondeleterious effects in humans based on their presence at high allele frequencies in other primate populations. We use this resource to classify 6% of all possible human protein-altering variants as likely benign and impute the pathogenicity of the remaining 94% of variants with deep learning, achieving state-of-the-art accuracy for diagnosing pathogenic variants in patients with genetic diseases.

Purpose: TANGO2 deficiency disorder (TDD), an autosomal recessive disease first reported in 2016, is characterized by neurodevelopmental delay, seizures, intermittent ataxia, hypothyroidism, and life-threatening metabolic and cardiac crises. The purpose of this study was to define the natural history of TDD., Methods: Data were collected from an ongoing natural history study of patients with TDD enrolled between February 2019 and May 2022. Data were obtained through phone or video based parent interviews and medical record review., Results: Data were collected from 73 patients (59% male) from 57 unrelated families living in 16 different countries. The median age of participants at the time of data collection was 9.0 years (interquartile range = 5.3-15.9 years, range = fetal to 31.8 years). A total of 24 different TANGO2 alleles were observed. Patients showed normal development in early infancy, with progressive delay in developmental milestones thereafter. Symptoms included ataxia, dystonia, and speech difficulties, typically starting between the ages of 1 to 3 years. A total of 46/71 (65%) patients suffered metabolic crises, and of those, 30 (65%) developed cardiac crises. Metabolic crises were significantly decreased after the initiation of B-complex or multivitamin supplementation., Conclusion: We provide the most comprehensive review of natural history of TDD and important observational data suggesting that B-complex or multivitamins may prevent metabolic crises., Competing Interests: Conflict of Interest The Department of Molecular and Human Genetics at Baylor College of Medicine derives revenue from clinical laboratory testing conducted at Baylor Genetics. S.A.M. serves as Scientific Advisor and consultant for Aytu BioPharma, Inc for the AR101 Enzastaurin Clinical Trial. All other authors declare no conflicts of interest., (Copyright © 2022 American College of Medical Genetics and Genomics. Published by Elsevier Inc. All rights reserved.)

Reference population databases are an essential tool in variant and gene interpretation. Their use guides the identification of pathogenic variants amidst the sea of benign variation present in every human genome, and supports the discovery of new disease-gene relationships. The Genome Aggregation Database (gnomAD) is currently the largest and most widely used publicly available collection of population variation from harmonized sequencing data. The data is available through the online gnomAD browser (https://gnomad.broadinstitute.org/) that enables rapid and intuitive variant analysis. This review provides guidance on the content of the gnomAD browser, and its usage for variant and gene interpretation. We introduce key features including allele frequency, per-base expression levels, constraint scores, and variant co-occurrence, alongside guidance on how to use these in analysis, with a focus on the interpretation of candidate variants and novel genes in rare disease., (© 2022 The Authors. Human Mutation Published by Wiley Periodicals LLC.)

Hundreds of thousands of genetic variants have been reported to cause severe monogenic diseases, but the probability that a variant carrier develops the disease (termed penetrance) is unknown for virtually all of them. Additionally, the clinical utility of common polygenetic variation remains uncertain. Using exome sequencing from 77,184 adult individuals (38,618 multi-ancestral individuals from a type 2 diabetes case-control study and 38,566 participants from the UK Biobank, for whom genotype array data were also available), we apply clinical standard-of-care gene variant curation for eight monogenic metabolic conditions. Rare variants causing monogenic diabetes and dyslipidemias display effect sizes significantly larger than the top 1% of the corresponding polygenic scores. Nevertheless, penetrance estimates for monogenic variant carriers average 60% or lower for most conditions. We assess epidemiologic and genetic factors contributing to risk prediction in monogenic variant carriers, demonstrating that inclusion of polygenic variation significantly improves biomarker estimation for two monogenic dyslipidemias.

Squamous cell carcinoma (SCC) occurs frequently in the human Xeroderma Pigmentosum (XP) syndrome and is characterized by deficient UV-damage repair. SCC is the most common equine ocular cancer and the only associated genetic risk factor is a UV-damage repair protein. Specifically, a missense mutation in horse DDB2 (T338M) was strongly associated with both limbal SCC and third eyelid SCC in three breeds of horses (Halflinger, Belgian, and Rocky Mountain Horses) and was hypothesized to impair binding to UV-damaged DNA. Here, we investigate DDB2-T338M mutant's capacity to recognize UV lesions in vitro and in vivo, together with human XP mutants DDB2-R273H and -K244E. We show that the recombinant DDB2-T338M assembles with DDB1, but fails to show any detectable binding to DNA substrates with or without UV lesions, due to a potential structural disruption of the rigid DNA recognition β-loop. Consistently, we demonstrate that the cellular DDB2-T338M is defective in its recruitment to focally radiated DNA damages, and in its access to chromatin. Thus, we provide direct functional evidence indicating the DDB2-T338M recapitulates molecular defects of human XP mutants, and is the causal loss-of-function allele that gives rise to equine ocular SCCs. Our findings shed new light on the mechanism of DNA recognition by UV-DDB and on the initiation of ocular malignancy., (Copyright © 2020 Elsevier B.V. All rights reserved.)

The acceleration of DNA sequencing in samples from patients and population studies has resulted in extensive catalogues of human genetic variation, but the interpretation of rare genetic variants remains problematic. A notable example of this challenge is the existence of disruptive variants in dosage-sensitive disease genes, even in apparently healthy individuals. Here, by manual curation of putative loss-of-function (pLoF) variants in haploinsufficient disease genes in the Genome Aggregation Database (gnomAD) 1 , we show that one explanation for this paradox involves alternative splicing of mRNA, which allows exons of a gene to be expressed at varying levels across different cell types. Currently, no existing annotation tool systematically incorporates information about exon expression into the interpretation of variants. We develop a transcript-level annotation metric known as the 'proportion expressed across transcripts', which quantifies isoform expression for variants. We calculate this metric using 11,706 tissue samples from the Genotype Tissue Expression (GTEx) project 2 and show that it can differentiate between weakly and highly evolutionarily conserved exons, a proxy for functional importance. We demonstrate that expression-based annotation selectively filters 22.8% of falsely annotated pLoF variants found in haploinsufficient disease genes in gnomAD, while removing less than 4% of high-confidence pathogenic variants in the same genes. Finally, we apply our expression filter to the analysis of de novo variants in patients with autism spectrum disorder and intellectual disability or developmental disorders to show that pLoF variants in weakly expressed regions have similar effect sizes to those of synonymous variants, whereas pLoF variants in highly expressed exons are most strongly enriched among cases. Our annotation is fast, flexible and generalizable, making it possible for any variant file to be annotated with any isoform expression dataset, and will be valuable for the genetic diagnosis of rare diseases, the analysis of rare variant burden in complex disorders, and the curation and prioritization of variants in recall-by-genotype studies.

Genetic variants that inactivate protein-coding genes are a powerful source of information about the phenotypic consequences of gene disruption: genes that are crucial for the function of an organism will be depleted of such variants in natural populations, whereas non-essential genes will tolerate their accumulation. However, predicted loss-of-function variants are enriched for annotation errors, and tend to be found at extremely low frequencies, so their analysis requires careful variant annotation and very large sample sizes 1 . Here we describe the aggregation of 125,748 exomes and 15,708 genomes from human sequencing studies into the Genome Aggregation Database (gnomAD). We identify 443,769 high-confidence predicted loss-of-function variants in this cohort after filtering for artefacts caused by sequencing and annotation errors. Using an improved model of human mutation rates, we classify human protein-coding genes along a spectrum that represents tolerance to inactivation, validate this classification using data from model organisms and engineered human cells, and show that it can be used to improve the power of gene discovery for both common and rare diseases.

Squamous cell carcinoma (SCC) is the most common cancer of the equine eye, frequently originating at the limbus, with the potential to invade the cornea, cause visual impairment, and result in loss of the eye. Several breeds of horses have a high occurrence of limbal SCC implicating a genetic basis for limbal SCC predisposition. Pedigree analysis in the Haflinger breed supports a simple recessive mode of inheritance and a genome-wide association study (N = 23) identified a 1.5 Mb locus on ECA12 significantly associated with limbal SCC (P corrected = 0.04). Sequencing the most physiologically relevant gene from this locus, damage specific DNA binding protein 2 (DDB2), identified a missense mutation (c.1013 C > T p.Thr338Met) that was strongly associated with limbal SCC (P = 3.41 × 10 -10 ). Genotyping 42 polymorphisms narrowed the ECA12 candidate interval to 483 kb but did not identify another variant that was more strongly associated. DDB2 binds to ultraviolet light damaged DNA and recruits other proteins to perform global genome nucleotide excision repair. Computational modeling predicts this mutation to be deleterious by altering conformation of the β loop involved in photolesion recognition. This DDB2 variant was also detected in two other closely related breeds with reported cases of ocular SCC, the Belgian and the Percheron, suggesting it may also be a SCC risk factor in these breeds. Furthermore, in humans xeroderma pigmentosum complementation group E, a disease characterized by sun sensitivity and increased risk of cutaneous SCC and melanomas, is explained by mutations in DDB2. Cross-species comparison remains to be further evaluated., (© 2017 UICC.)