ClinVar and HGMD genomic variant classification accuracy has improved over time, as measured by implied disease burden

BACKGROUND

Curated databases of genetic variants assist clinicians and researchers in interpreting genetic variation. Yet, these databases contain some misclassified variants. It is unclear whether variant misclassification is abating as these databases rapidly grow and implement new guidelines.

METHODS

Using archives of ClinVar and HGMD, we investigated how variant misclassification has changed over 6 years, across different ancestry groups. We considered inborn errors of metabolism (IEMs) screened in newborns as a model system because these disorders are often highly penetrant with neonatal phenotypes. We used samples from the 1000 Genomes Project (1KGP) to identify individuals with genotypes that were classified by the databases as pathogenic. Due to the rarity of IEMs, nearly all such classified pathogenic genotypes indicate likely variant misclassification in ClinVar or HGMD.

RESULTS

While the false-positive rates of both ClinVar and HGMD have improved over time, HGMD variants currently imply two orders of magnitude more affected individuals in 1KGP than ClinVar variants. We observed that African ancestry individuals have a significantly increased chance of being incorrectly indicated to be affected by a screened IEM when HGMD variants are used. However, this bias affecting genomes of African ancestry was no longer significant once common variants were removed in accordance with recent variant classification guidelines. We discovered that ClinVar variants classified as Pathogenic or Likely Pathogenic are reclassified sixfold more often than DM or DM? variants in HGMD, which has likely resulted in ClinVar's lower false-positive rate.

CONCLUSIONS

Considering misclassified variants that have since been reclassified reveals our increasing understanding of rare genetic variation. We found that variant classification guidelines and allele frequency databases comprising genetically diverse samples are important factors in reclassification. We also discovered that ClinVar variants common in European and South Asian individuals were more likely to be reclassified to a lower confidence category, perhaps due to an increased chance of these variants being classified by multiple submitters. We discuss features for variant classification databases that would support their continued improvement.

Distinct sequence features underlie microdeletions and gross deletions in the human genome

Microdeletions and gross deletions are important causes (~20%) of human inherited disease and their genomic locations are strongly influenced by the local DNA sequence environment. This notwithstanding, no study has systematically examined their underlying generative mechanisms. Here, we obtained 42,098 pathogenic microdeletions and gross deletions from the Human Gene Mutation Database (HGMD) that together form a continuum of germline deletions ranging in size from 1 to 28,394,429 bp. We analyzed the DNA sequence within 1 kb of the breakpoint junctions and found that the frequencies of non‐B DNA‐forming repeats, GC‐content, and the presence of seven of 78 specific sequence motifs in the vicinity of pathogenic deletions correlated with deletion length for deletions of length ≤30 bp. Further, we found that the presence of DR, GQ, and STR repeats is important for the formation of longer deletions (>30 bp) but not for the formation of shorter deletions (≤30 bp) while significantly (χ², p < 2E−16) more microhomologies were identified flanking short deletions than long deletions (length >30 bp). We provide evidence to support a functional distinction between microdeletions and gross deletions. Finally, we propose that a deletion length cut‐off of 25–30 bp may serve as an objective means to functionally distinguish microdeletions from gross deletions.

Down-regulation of SORL1 is associated with Alzheimer's disease through activating ABC transporter pathway

Alzheimer's disease (AD) is a common neurodegenerative disease with high morbidity among elderly people. A genetic attribution has been extensively proved. Here, we propose to further prioritize genes that harbor single nucleotide variation (SNV) or structural variation (SV) for AD and explore the underlying potential mechanisms through exploiting their expression and methylation spectra. A high-confidence AD-associated candidate gene list was obtained from the ClinVar and Human Gene Mutation Database (HGMD). Genome-wide methylation and expression profiles of AD and normal subjects were downloaded from the Gene Expression Omnibus (GEO). Through comprehensive comparison of expression and methylation levels between AD and normal samples, as well as different stages of AD samples, SORL1 was identified as the most plausible gene for AD incidence and progression. Gene Set Enrichment Analysis (GSEA) revealed significant activation of the ABC (ATP binding cassette) transporter with the aberrant up-regulation of SORL1 within AD samples. This study unfolds the expression and methylation spectra of previously probed genes with SNV or SV in AD for the first time, and reports an aberrant activation of the ABC transporter pathway that might contribute to AD progression. This should shed some light on AD diagnosis and precision treatment.

Patients with coronary heart disease, dilated cardiomyopathy and idiopathic ventricular tachycardia share overlapping patterns of pathogenic variation in cardiac risk genes

Background

Ventricular tachycardia (VT) is a major cause of sudden cardiac death (SCD). Clinical investigations can sometimes fail to identify the underlying cause of VT and the event is classified as idiopathic (iVT). VT contributes significantly to the morbidity and mortality in patients with coronary artery disease (CAD) and dilated cardiomyopathy (DCM). Since mutations in arrhythmia-associated genes frequently determine arrhythmia susceptibility screening for disease-predisposing variants could improve VT diagnostics and prevent SCD in patients.

Methods

Ninety-two patients diagnosed with coronary heart disease (CHD), DCM, or iVT were included in our study. We evaluated genetic profiles and variants in known cardiac risk genes by targeted next generation sequencing (NGS) using a newly designed custom panel of 96 genes. We hypothesized that shared morphological and phenotypical features among these subgroups may have an overlapping molecular base. To our knowledge, this was the first study of the deep sequencing of 96 targeted cardiac genes in Kazakhstan. The clinical significance of the sequence variants was interpreted according to the guidelines developed by the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) in 2015. The ClinVar and Varsome databases were used to determine the variant classifications.

Results

Targeted sequencing and stepwise filtering of the annotated variants identified a total of 307 unique variants in 74 genes, totally 456 variants in the overall study group. We found 168 mutations listed in the Human Genome Mutation Database (HGMD) and another 256 rare/unique variants with elevated pathogenic potential. There was a predominance of high- to intermediate pathogenicity variants in LAMA2, MYBPC3, MYH6, KCNQ1, GAA, and DSG2 in CHD VT patients. Similar frequencies were observed in DCM VT, and iVT patients, pointing to a common molecular disease association. TTN, GAA, LAMA2, and MYBPC3 contained the most variants in the three subgroups which confirm the impact of these genes in the complex pathogenesis of cardiomyopathies and VT. The classification of 307 variants according to ACMG guidelines showed that nine (2.9%) variants could be classified as pathogenic, nine (2.9%) were likely pathogenic, 98 (31.9%) were of uncertain significance, 73 (23.8%) were likely benign, and 118 (38.4%) were benign. CHD VT patients carry rare genetic variants with increased pathogenic potential at a comparable frequency to DCM VT and iVT patients in genes related to sarcomere function, nuclear function, ion flux, and metabolism.

Conclusions

In this study we showed that in patients with VT secondary to coronary artery disease, DCM, or idiopathic etiology multiple rare mutations and clinically significant sequence variants in classic cardiac risk genes associated with cardiac channelopathies and cardiomyopathies were found in a similar pattern and at a comparable frequency.

Exome data clouds the pathogenicity of genetic variants in Pulmonary Arterial Hypertension

BACKGROUND

We aimed to provide a set of previously reported PAH-associated missense and nonsense variants, and evaluate the pathogenicity of those variants.

METHODS

The Human Gene Mutation Database, PubMed, and Google Scholar were searched for previously reported PAH-associated genes and variants. Thereafter, both exome sequencing project and exome aggregation consortium as background population searched for previously reported PAH-associated missense and nonsense variants. The pathogenicity of previously reported PAH-associated missense variants evaluated by using four in silico prediction tools.

RESULTS

In total, 14 PAH-associated genes and 180 missense and nonsense variants were gathered. The BMPR2, the most frequent reported gene, encompasses 135 of 180 missense and nonsense variants. The exome sequencing project comprised 9, and the exome aggregation consortium counted 25 of 180 PAH-associated missense and nonsense variants. The TOPBP1 and ENG genes are unlikely to be the monogenic cause of PAH pathogenesis based on allele frequency in background population and prediction analysis.

CONCLUSION

This is the first evaluation of previously reported PAH-associated missense and nonsense variants. The BMPR2 identified as the major gene out of 14 PAH-associated genes. Based on findings, the ENG and TOPBP1 gene are not likely to be the monogenic cause of PAH.

Aggregation of population-based genetic variation over protein domain homologues and its potential use in genetic diagnostics

Whole exomes of patients with a genetic disorder are nowadays routinely sequenced but interpretation of the identified genetic variants remains a major challenge. The increased availability of population‐based human genetic variation has given rise to measures of genetic tolerance that have been used, for example, to predict disease‐causing genes in neurodevelopmental disorders. Here, we investigated whether combining variant information from homologous protein domains can improve variant interpretation. For this purpose, we developed a framework that maps population variation and known pathogenic mutations onto 2,750 “meta‐domains.” These meta‐domains consist of 30,853 homologous Pfam protein domain instances that cover 36% of all human protein coding sequences. We find that genetic tolerance is consistent across protein domain homologues, and that patterns of genetic tolerance faithfully mimic patterns of evolutionary conservation. Furthermore, for a significant fraction (68%) of the meta‐domains high‐frequency population variation re‐occurs at the same positions across domain homologues more often than expected. In addition, we observe that the presence of pathogenic missense variants at an aligned homologous domain position is often paired with the absence of population variation and vice versa. The use of these meta‐domains can improve the interpretation of genetic variation.

Integration of 60,000 exomes and ACMG guidelines question the role of Catecholaminergic Polymorphic Ventricular Tachycardia-associated variants

Catecholaminergic Polymorphic Ventricular Tachycardia (CPVT) is a highly lethal cardiac arrhythmia disease occurring during exercise or psychological stress. CPVT has an estimated prevalence of 1:10,000 and has mainly been associated with variants in calcium-regulating genes. Identification of potential false-positive pathogenic variants was conducted by searching the Exome Aggregation Consortium (ExAC) database (n = 60,706) for variants reported to be associated with CPVT. The pathogenicity of the interrogated variants was assessed using guidelines from the American College of Medical Genetics and Genomics (ACMG) and in silico prediction tools. Of 246 variants 38 (15%) variants previously associated with CPVT were identified in the ExAC database. We predicted the CPVT prevalence to be 1:132. The ACMG standards classified 29% of ExAC variants as pathogenic or likely pathogenic. The in silico predictions showed a reduced probability of disease-causing effect for the variants identified in the exome database (p < 0.001). We have observed a large overrepresentation of previously CPVT-associated variants in a large exome database. Based on the frequency of CPVT in the general population, it is less likely that the previously proposed variants are associated with a highly penetrant monogenic form of the disease.

How well do whole exome sequencing results correlate with medical findings? A study of 89 Mayo Clinic Biobank samples

Whole exome sequencing (WES) is increasingly being used for diagnosis without adequate information on predictive characteristics of reportable variants typically found on any given individual and correlation with clinical phenotype. In this study, we performed WES on 89 deceased individuals (mean age at death 74 years, range 28-93) from the Mayo Clinic Biobank. Significant clinical diagnoses were abstracted from electronic medical record via chart review. Variants [Single Nucleotide Variant (SNV) and insertion/deletion] were filtered based on quality (accuracy >99%, read-depth >20, alternate-allele read-depth >5, minor-allele-frequency <0.1) and available HGMD/OMIM phenotype information. Variants were defined as Tier-1 (nonsense, splice or frame-shifting) and Tier-2 (missense, predicted-damaging) and evaluated in 56 ACMG-reportable genes, 57 cancer-predisposition genes, along with examining overall genotype-phenotype correlations. Following variant filtering, 7046 total variants were identified (~79/person, 644 Tier-1, 6402 Tier-2), 161 among 56 ACMG-reportable genes (~1.8/person, 13 Tier-1, 148 Tier-2), and 115 among 57 cancer-predisposition genes (~1.3/person, 3 Tier-1, 112 Tier-2). The number of variants across 57 cancer-predisposition genes did not differentiate individuals with/without invasive cancer history (P > 0.19). Evaluating genotype-phenotype correlations across the exome, 202(3%) of 7046 filtered variants had some evidence for phenotypic correlation in medical records, while 3710(53%) variants had no phenotypic correlation. The phenotype associated with the remaining 44% could not be assessed from a typical medical record review. These data highlight significant continued challenges in the ability to extract medically meaningful predictive results from WES.

A novel mutation in MID1 in a patient with X-linked Opitz G/BBB syndrome

Opitz G/BBB syndrome (OS) is a genetically heterogeneous disease. We report on an OS patient with a novel inherited mutation in MID1. Metaphase analysis showed a normal male karyotype. Array CGH revealed a maternally inherited duplication at Xp22.31 (6,467,203-7,992,261, hg18), the size was estimated to 1.5Mb. Sequence analysis of the MID1 coding region revealed a novel missense mutation in exon 8 (c.1561C>T/p. R521C) which resulted in an ammonia acid substitution (R521C) in the PRX domain of the MID1 protein. The mutation was inherited from unaffected grandmother and mildly affected mother. Prenatal diagnosis was performed for the third pregnancy after identification of the causative mutation in the family. The third fetus was found to be a female carrier. Postnatal follow-up at 2-month-old showed normal phenotype. In conclusion, we reported a familial OS patient with a novel mutation in exon 8 which provided another evidence for that mutation clustered in C-terminal domain of MID1. The newly identified mutation in our patient expands mutation spectrum in MID1 gene.

Two novel mutations in the SLCO2A1 gene in a Chinese patient with primary hypertrophic osteoarthropathy

Primary hypertrophic osteoarthropathy (PHO) is a rare monogenetic disease characterized by digital clubbing, periostosis and pachydermia. Mutations in the 15-hydroxy-prostaglandin dehydrogenase (HPGD) gene and solute carrier organic anion transporter family member 2A1 (SLCO2A1) gene have been shown to be associated with PHO. Here, we described clinical characteristics in a Chinese patient with PHO, and identified two novel mutations in SLCO2A1: a heterozygous guanine-to-thymidine transition at the invariant -1 position of the acceptor site of intron 2 (c.235-1G>T) and a heterozygous missense mutation p.Pro219Leu (c.656C>T) in exon 5.

Identification of mutation in NPC2 by exome sequencing results in diagnosis of Niemann-Pick disease type C

We report identification of a homozygous mutation in NPC2 in two Iranian siblings with a neurologic dysfunction whose disease had not been diagnosed prior to our genetic analysis. The mutation was identified by exome sequencing. The finding resulted in diagnosis of Niemann-Pick disease type C (NPC) in the siblings, and initiation of treatment with Miglustat. The clinical features of the patients are presented. It has been suggested that NPC is under diagnosed, particularly when presentations are not very severe, as was the situation in the cases studied here. NPC is a fatal autosomal recessive disorder clinically characterized by hepatosplenomegaly and progressive neurological deterioration. At the cellular level, it causes aberrant cholesterol trafficking and accumulation of unesterified cholesterol in lysosomes. Mutations in NPC1 and NPC2 are cause of disease in respectively, 95% and 5% of NPC patients. The p.Pro120Ser causing mutation in NPC2 observed in the Iranian patients was earlier observed in the only other NPC2 patient reported from the Middle East. The study demonstrates that in addition to greatly facilitating gene discovery, exome sequencing has notable potentials for diagnosis, particularly for diagnosis of atypical cases.

The effects of non-synonymous single nucleotide polymorphisms (nsSNPs) on protein-protein interactions

Non-synonymous single nucleotide polymorphisms (nsSNPs) are single base changes leading to a change to the amino acid sequence of the encoded protein. Many of these variants are associated with disease, so nsSNPs have been well studied, with studies looking at the effects of nsSNPs on individual proteins, for example, on stability and enzyme active sites. In recent years, the impact of nsSNPs upon protein-protein interactions has also been investigated, giving a greater insight into the mechanisms by which nsSNPs can lead to disease. In this review, we summarize these studies, looking at the various mechanisms by which nsSNPs can affect protein-protein interactions. We focus on structural changes that can impair interaction, changes to disorder, gain of interaction, and post-translational modifications before looking at some examples of nsSNPs at human-pathogen protein-protein interfaces and the analysis of nsSNPs from a network perspective.

Mutation spectrum of primary hyperoxaluria type 1 in Tunisia: implication for diagnosis in North Africa

BACKGROUND

Primary hyperoxaluria type 1 (PH1) is an autosomal recessive inherited metabolic disease, characterized by progressive kidney failure due to renal deposition of calcium oxalate. Mutations in the AGXT gene, encoding the liver-specific enzyme alanine glyoxylate aminotransferase, are responsible for the disease. We aimed to determine the mutational spectrum causing PH1 and to provide an accurate tool for diagnosis as well as for prenatal diagnosis in the affected families.

METHODS

Direct sequencing was used to detect mutations in the AGXT gene in DNA samples from 13 patients belonging to 12 Tunisian families.

RESULTS

Molecular analysis revealed five mutations causing PH1 in Tunisia. The mutations were identified along exons 1, 2, 4, 5 and 7. The most predominant mutations were the Maghrebian "p.I244T" and the Arabic "p.G190R". Furthermore, three other mutations characteristic of different ethnic groups were found in our study population. These results confirm the mutational heterogeneity related to PH1 in Tunisian population. All the mutations are in a homozygous state, reflecting the high impact of endogamy in our population.

CONCLUSION

Mutation analysis through DNA sequencing can provide a useful first line investigation for PH1. This identification could provide an accurate tool for prenatal diagnosis, genetic counseling and screen for potential presymptomatic individuals.

A novel mutation in a large family causes a unique phenotype of Mucolipidosis IV

Mucolipidosis type IV is a rare autosomal recessive lysosomal storage disorder reported among Ashkenazi Jews and to a lesser extent in other ethnic groups. Several mutations have been reported in MCOLN1 which is the only known gene associated with the disorder. Here we report the first Saudi patient with Mucolipidosis type IV from a consanguineous family with two branches having a total of five patients carrying a novel transition mutation, c.1307A>G (p.Y436C) in exon 11. The clinical course of the patient was nonspecific and a lysosomal storage disorder was not highly suspected due to lack of coarse facial features, organomegaly and skeletal findings of dysostosis multiplex. The detailed bioinformatics analysis on the deleterious effects of the mutation is discussed. Emphasis is made on the importance of brain magnetic resonance imaging (MRI) findings and serum gastrin level as key clues to the diagnosis of this often subtle neurodevelopmental disorder.