Conserved missense variant pathogenicity and correlated phenotypes across paralogous genes

Background The majority of missense variants in clinical genetic tests are classified as variants of uncertain significance. Prior research has shown that the deleterious effects and the subsequent molecular consequence of variants are often conserved among paralogous protein sequences within a gene family. Here, we systematically quantified on an exome-wide scale if the existence of pathogenic variants in paralogous genes at a conserved position could serve as evidence for the pathogenicity of a new variant. For the gene family of voltage-gated sodium channels where variants and expert-curated clinical phenotypes were available, we also assessed whether phenotype patterns of multiple disorders for each gene were also conserved across variant positions within the gene family. Methods We developed a framework that assesses the presence of pathogenic missense variants located in conserved residues across paralogous genes. We systematically mapped 2.5 million pathogenic and general population variants from the ClinVar, HGMD, and gnomAD databases onto a total of 9,990 genes and aligned them by gene families. We evaluated the quantity of classifiable amino acids by utilizing pathogenic variants identified in databases alone and then compared this assessment to the inclusion of paralogous pathogenic variants. We validated and quantified the evidence of conserved pathogenic paralogous variants in variant pathogenicity classification. Results Considering conserved pathogenic variants in paralogous genes, increased the number of classifiable variants 2.8-fold across the exome, compared to pathogenic variants in the gene of interest alone. The presence of a pathogenic variant in a paralogous gene is associated with a positive likelihood ratio of 8.32 for variant pathogenicity. The likelihood ratio was gene family-specific. Across ten genes encoding voltage-gated sodium channels and 22 expert-curated disorders, we identified cross-paralog correlated phenotypes based on 3D structure spatial position. For example, the established loss-of-function disorders SCN1A-associated Dravet syndrome, SCN2A-associated autism, SCN5A-associated Brugarda Syndrome, and SCN8A-associated neurodevelopmental disorder without seizures were correlated in their spatial variant position on structure. Finally, we show that phenotype integration in paralog variant selection improves variant classification. Conclusion Our results show that paralogous variants, in particular with phenotype information can enhance our understanding of variant effects.

The population frequency of AD polycystic kidney disease (ADPKD) has not been known with certainty. We hypothesized that more accurate population frequencies for typical (PKD1, PKD2) and atypical (GANAB, ALG9, DNAJB11 ALG5, IFT140, NEK8) ADPKD were possible using a computational approach. Initially, this study calculated the number of predicted pathogenic structural, null, and missense variants in the ADPKD-associated genes in gnomAD v.2.1.1, and compared this with disease-causing changes according to the ClinVar, HGMD, or LOVD databases. However, because of the difficulty assessing missense variant pathogenicity, ClinVar assessments were used instead for gnomAD v.4.1, which includes a larger cohort, and more structural and copy number data. Predicted pathogenic variants were found in the typical ADPKD genes in 1 in 314 people, with PKD1 and PKD2 changes present in 1 in 417 and 916, respectively, using ClinVar assessments of gnomAD v.4.1 variants and loss-of-function structural and copy number changes. Predicted pathogenic null and missense variants in the atypical ADPKD genes were found in 1 in 283 people in gnomAD v.4.1 using ClinVar. No pathogenic missense changes in the kinase domain of NEK8 previously-reported in monoallelic disease were present. Predicted pathogenic variants in the genes associated with typical and atypical ADPKD are more common than suspected previously. Predicted pathogenic variants are more common for atypical than for typical ADPKD, especially when loss-of-function structural and copy number variants are included.

eP061: Genetic risk for breast and ovarian cancer in a diverse and unselected population

ACMG SF v3.0 list for reporting of secondary findings in clinical exome and genome sequencing: a policy statement of the American College of Medical Genetics and Genomics (ACMG)

BackgroundGermline variants of ten keratin genes (K1, K2, K5, K6A, K6B, K9, K10, K14, K16, and K17) have been reported for causing different types of genodermatoses with an autosomal dominant mode of inheritance. Among all the variants of these ten keratin genes, most of them are missense variants. Unlike pathogenic and likely pathogenic variants, understanding the clinical importance of novel missense variants or variants of uncertain significance (VUS) is the biggest challenge for clinicians or medical geneticists. Functional characterization is the only way to understand the clinical association of novel missense variants or VUS but it is time consuming, costly, and depends on the availability of patient’s samples. Existing databases report the pathogenic variants of the keratin genes, but never emphasize the systematic effects of these variants on keratin protein structure and genotype-phenotype correlation.ResultsTo address this need, we developed a comprehensive database KVarPredDB, which contains information of all ten keratin genes associated with genodermatoses. We integrated and curated 400 reported pathogenic missense variants as well as 4629 missense VUS. KVarPredDB predicts the pathogenicity of novel missense variants as well as to understand the severity of disease phenotype, based on four criteria; firstly, the difference in physico-chemical properties between the wild type and substituted amino acids; secondly, the loss of inter/intra-chain interactions; thirdly, evolutionary conservation of the wild type amino acids and lastly, the effect of the substituted amino acids in the heptad repeat. Molecular docking simulations based on resolved crystal structures were adopted to predict stability changes and get the binding energy to compare the wild type protein with the mutated one. We use this basic information to determine the structural and functional impact of novel missense variants on the keratin coiled-coil heterodimer. KVarPredDB was built under the integrative web application development framework SSM (SpringBoot, Spring MVC, MyBatis) and implemented in Java, Bootstrap, React-mutation-mapper, MySQL, Tomcat. The website can be accessed through http://bioinfo.zju.edu.cn/KVarPredDB. The genomic variants and analysis results are freely available under the Creative Commons license.ConclusionsKVarPredDB provides an intuitive and user-friendly interface with computational analytical investigation for each missense variant of the keratin genes associated with genodermatoses.

Simple SummaryHearing loss is the most common sensory impairment in humans. Globally, GJB2 is the most common responsible gene, with missense variants being the most frequent variant type. Pathogenic missense variants in the GJB2 gene cause isolated hearing loss and hearing loss combined with skin diseases. However, the mechanism by which different missense variants cause different hearing loss phenotypes is unknown. We have summarized 84 functionally studied missense variants and reviewed the clinical phenotypes and the molecular mechanisms that affected connexon and gap junction functions. We have also compiled a comprehensive dataset of GJB2 missense variants from public databases and published literature and found that over 2/3 of the GJB2 missense variants are currently classified as variants of uncertain significance. With the development of cutting-edge deep mutational scanning technology, we predict that the molecular mechanisms by which all missense variants of GJB2 lead to different clinical phenotypes will be elucidated. These gained insights will facilitate the use of genetic testing in the prevention and control of hearing loss.The GJB2 gene is the most common gene responsible for hearing loss (HL) worldwide, and missense variants are the most abundant type. GJB2 pathogenic missense variants cause nonsyndromic HL (autosomal recessive and dominant) and syndromic HL combined with skin diseases. However, the mechanism by which these different missense variants cause the different phenotypes is unknown. Over 2/3 of the GJB2 missense variants have yet to be functionally studied and are currently classified as variants of uncertain significance (VUS). Based on these functionally determined missense variants, we reviewed the clinical phenotypes and investigated the molecular mechanisms that affected hemichannel and gap junction functions, including connexin biosynthesis, trafficking, oligomerization into connexons, permeability, and interactions between other coexpressed connexins. We predict that all possible GJB2 missense variants will be described in the future by deep mutational scanning technology and optimizing computational models. Therefore, the mechanisms by which different missense variants cause different phenotypes will be fully elucidated.

Pathogenic missense and truncating variants in the GABRG2 gene cause a spectrum of epilepsies, from Dravet syndrome to milder simple febrile seizures. In most cases, pathogenic missense variants in the GABRG2 gene segregate with a febrile seizure phenotype. In this case series, we report a recurrent, de novo missense variant (c0.316 G > A; p.A106T) in the GABRG2 gene that was identified in five unrelated individuals. These patients were described to have a more severe phenotype than previously reported for GABRG2 missense variants. Common features include variable early-onset seizures, significant motor and speech delays, intellectual disability, hypotonia, movement disorder, dysmorphic features and vision/ocular issues. Our report further explores a recurrent pathogenic missense variant within the GABRG2 variant family and broadens the spectrum of associated phenotypes for GABRG2-associated disorders.

Functional Testing of Variants of Uncertain Significance in TERC, TERT,and RTEL1 from Adult Patients with Telomere Biology Disorders

Background and aims Alport syndrome (AS) and thin basement membrane nephropathy (TBMN) are genetically heterogeneous kidney diseases that manifest with kidney involvement, AS also with hearing and visual impairment. Electron microscopic examination shows characteristic ultrastructural changes in the glomerular basement membrane (GBM) in AS and TBMN, but cannot reliably distinguish between the two entities. Therefore, molecular genetic studies and the identification of individual pathogenic variants in the target genes are crucial to confirm the diagnosis. The aim of our study is to determine the significance of the genetic variants associated with the ultrastructural changes of GBM found in the renal biopsy. Method We retrospectively analyzed records of patients from the University Medical Centre Ljubljana in which ultrastructural GBM changes characteristic for AS and TBMN were found in the renal biopsy, including their clinical, laboratory and histological features. Molecular genetic investigations were made using the new generation sequencing method (NGS) to determine germline pathogenic variants in genes directly related to AS and TBMN (COL4A3, COL4A4 and COL4A5) and in genes that are important in syndromes with a related clinical picture (COL4A6, CD151, LMX1B, MYH9, FN1 and MYO1E). The pathogenicity of variants was determined according to the Clinical Molecular Genetic Society best practice guidelines. Results The study included 73 patients with ultrastructurally confirmed thin GBM. Indications for kidney biopsy were chronic nephritis syndrome in 33 patients (45%), acute worsening of renal function in 16 (22%), nephrotic syndrome in 9 (12%), nephrotic proteinuria and chronic kidney disease in 6 (8%), isolated proteinuria in 5 (7%) and suspected vasculitis in 4 patients (6%). At the time of kidney biopsy average age was 50 ± 16 years, average creatinine 112 ± 70 µmol/L and average estimated glomerular filtration rate (eGFR) 65 ± 24 ml/min/1.73 m2, 7 patients had eGFR below 30 ml/min/1.73 m2. Average daily proteinuria was 2.4 ± 3.3 gr and 88% of patients had hematuria. 38% of them have hearing impairment and 43% had positive familiar history of kidney disease. We recorded 33 patients (45%) with genetic variants associated with AS or TBMN: 22 patients (67%) with pathogenic variants, 2 patients (6%) with likely pathogenic and 9 patients (27%) with variants of uncertain significance. The presence of a positive genetic analysis was associated with a higher level of hematuria (p = 0.003), a lower level of proteinuria (p = 0.004) and a positive family history of kidney disease (p < 0.001). It was also associated with hearing loss, but the association was not statistically significant (p = 0.099). Patients with genetic variants had similar eGFR to those with negative genetics (63.9 ml/min/1.73 m2 (58.3-69.4) vs 66.2 ml/min/1.73 m2 (60.0-72.3)), p = 0.698. We also found no differences in GBM thickness in patients without compared to those with identified gene variants (213 nm (160-315) vs. 209 nm (140-410)), p = 0.761. 32 patients (44%) had only focal chronic changes and ultrastructural changes of thin GBM present on kidney biopsy and 41 patients (56%) had also been diagnosed with other renal pathology. When comparing these two groups, patients without associated pathology had lower daily proteinuria (p = 0.048), higher level of hematuria (p = 0.021), frequent hearing impairment (p = 0.008) and also positive genetics analysis (p = 0.016). Conclusion We have shown that in patients with proven ultrastructural alterations in GBM, the presence of so far recognized genetic variants associated with AS and TBMN is 45%, and corresponds to the frequency of other phenotypic characteristics of AS, including hearing impairment and positive family history. In the future, additional clinical tools will be required to accurately determine the role of coexistent thin GBM on the prognosis of patients with other relevant kidney pathology and other genetic tools beyond NGS should be employed to explore other pathologic genetic variants associated with the AS spectrum.

Infertility is a heterogeneous condition, with genetic causes thought to underlie a substantial fraction of cases. Genome sequencing is becoming increasingly important for genetic diagnosis of diseases including idiopathic infertility; however, most rare or minor alleles identified in patients are variants of uncertain significance (VUS). Interpreting the functional impacts of VUS is challenging but profoundly important for clinical management and genetic counseling. To determine the consequences of these variants in key fertility genes, we functionally evaluated 11 missense variants in the genes ANKRD31, BRDT, DMC1, EXO1, FKBP6, MCM9, M1AP, MEI1, MSH4 and SEPT12 by generating genome-edited mouse models. Nine variants were classified as deleterious by most functional prediction algorithms, and two disrupted a protein-protein interaction (PPI) in the yeast two hybrid (Y2H) assay. Though these genes are essential for normal meiosis or spermiogenesis in mice, only one variant, observed in the MCM9 gene of a male infertility patient, compromised fertility or gametogenesis in the mouse models. To explore the disconnect between predictions and outcomes, we compared pathogenicity calls of missense variants made by ten widely used algorithms to 1) those annotated in ClinVar and 2) those evaluated in mice. All the algorithms performed poorly in terms of predicting the effects of human missense variants modeled in mice. These studies emphasize caution in the genetic diagnoses of infertile patients based primarily on pathogenicity prediction algorithms and emphasize the need for alternative and efficient in vitro or in vivo functional validation models for more effective and accurate VUS description to either pathogenic or benign categories.

Co-observation of a gene variant with a pathogenic variant in another gene that explains the disease presentation has been designated as evidence against pathogenicity for commonly used variant classification guidelines. Multiple variant curation expert panels have specified, from consensus opinion, that this evidence type is not applicable for the classification of breast cancer predisposition gene variants. Statistical analysis of sequence data for 55,815 individuals diagnosed with breast cancer from the BRIDGES sequencing project was undertaken to formally assess the utility of co-observation data for germline variant classification. Our analysis included expected loss-of-function variants in 11 breast cancer predisposition genes and pathogenic missense variants in BRCA1, BRCA2, and TP53. We assessed whether co-observation of pathogenic variants in two different genes occurred more or less often than expected under the assumption of independence. Co-observation of pathogenic variants in each of BRCA1, BRCA2, and PALB2 with the remaining genes was less frequent than expected. This evidence for depletion remained after adjustment for age at diagnosis, study design (familial versus population-based), and country. Co-observation of a variant of uncertain significance in BRCA1, BRCA2, or PALB2 with a pathogenic variant in another breast cancer gene equated to supporting evidence against pathogenicity following criterion strength assignment based on the likelihood ratio and showed utility in reclassification of missense BRCA1 and BRCA2 variants identified in BRIDGES. Our approach has applicability for assessing the value of co-observation as a predictor of variant pathogenicity in other clinical contexts, including for gene-specific guidelines developed by ClinGen Variant Curation Expert Panels.

Pathogenic variants in genes related to channelopathy and cardiomyopathy are the most common cause of sudden unexplained cardiac death. However, few reports have investigated the frequency and/or spectrum of pathogenic variants in these genes in Korean sudden cardiac arrest survivors. This study aimed to investigate the causative genetic variants of cardiac-associated genes in Korean sudden cardiac arrest survivors. We performed exome sequencing followed by filtering and validation of variants in 100 genes related to channelopathy and cardiomyopathy in 19 Korean patients who survived sudden cardiac arrest. Five of the 19 patients (26.3%) had either a pathogenic variant or a likely pathogenic variant in MYBPC3 (n=1), MYH7 (n=1), RYR2 (n=2), or TNNT2 (n=1). All five variants were missense variants that have been reported previously in patients with channelopathies or cardiomyopathies. Furthermore, an additional 12 patients (63.2%) had more than one variant of uncertain significance. In conclusion, pathogenic or likely pathogenic variants in genes related to channelopathy and cardiomyopathy are not uncommon in Korean sudden cardiac arrest survivors and cardiomyopathy-related genes should be included in the molecular diagnosis of sudden cardiac arrest in Korea.

Autosomal dominant optic atrophy is one of the most common inherited optic neuropathies. This disease is genetically heterogeneous, but most cases are due to pathogenic variants in the OPA1 gene: depending on the population studied, 32–90% of cases harbor pathogenic variants in this gene. The aim of this study was to provide a comprehensive overview of the entire spectrum of likely pathogenic variants in the OPA1 gene in a large cohort of patients. Over a period of 20 years, 755 unrelated probands with a diagnosis of bilateral optic atrophy were referred to our laboratory for molecular genetic investigation. Genetic testing of the OPA1 gene was initially performed by a combined analysis using either single-strand conformation polymorphism or denaturing high performance liquid chromatography followed by Sanger sequencing to validate aberrant bands or melting profiles. The presence of copy number variations was assessed using multiplex ligation-dependent probe amplification. Since 2012, genetic testing was based on next-generation sequencing platforms. Genetic screening of the OPA1 gene revealed putatively pathogenic variants in 278 unrelated probands which represent 36.8% of the entire cohort. A total of 156 unique variants were identified, 78% of which can be considered null alleles. Variant c.2708_2711del/p.(V903Gfs*3) was found to constitute 14% of all disease-causing alleles. Special emphasis was placed on the validation of splice variants either by analyzing cDNA derived from patients´ blood samples or by heterologous splice assays using minigenes. Splicing analysis revealed different aberrant splicing events, including exon skipping, activation of exonic or intronic cryptic splice sites, and the inclusion of pseudoexons. Forty-eight variants that we identified were novel. Nine of them were classified as pathogenic, 34 as likely pathogenic and five as variant of uncertain significance. Our study adds a significant number of novel variants to the mutation spectrum of the OPA1 gene and will thereby facilitate genetic diagnostics of patients with suspected dominant optic atrophy.

To report the genetic etiologies of Emery-Dreifuss muscular dystrophy (EDMD), limb-girdle muscular dystrophy (LGMD), congenital muscular dystrophy (CMD), and distal muscular dystrophy (DD) in 6 geographically defined areas of the United States. This was a cross-sectional, population-based study in which we studied the genes and variants associated with muscular dystrophy in individuals who were diagnosed with and received care for EDMD, LGMD, CMD, and DD from January 1, 2008, through December 31, 2016, in the 6 areas of the United States covered by the Muscular Dystrophy Surveillance, Tracking, and Research Network (MD STARnet). Variants of unknown significance (VUSs) from the original genetic test reports were reanalyzed for changes in interpretation. Among 243 individuals with definite or probable muscular dystrophy, LGMD was the most common diagnosis (138 cases), followed by CMD (62 cases), DD (22 cases), and EDMD (21 cases). There was a higher proportion of male individuals compared with female individuals, which persisted after excluding X-linked genes (EMD) and autosomal genes reported to have skewed gender ratios (ANO5, CAV3, and LMNA). The most common associated genes were FKRP, CAPN3, ANO5, and DYSF. Reanalysis yielded more definitive variant interpretations for 60 of 144 VUSs, with a mean interval between the original clinical genetic test of 8.11 years for all 144 VUSs and 8.62 years for the 60 reclassified variants. Ten individuals were found to have monoallelic pathogenic variants in genes known to be primarily recessive. This study is distinct for being an examination of 4 types of muscular dystrophies in selected geographic areas of the United States. The striking proportion of resolved VUSs demonstrates the value of periodic re-examinations of these variants. Such re-examinations will resolve some genetic diagnostic ambiguities before initiating repeat testing or more invasive diagnostic procedures such as muscle biopsy. The presence of monoallelic pathogenic variants in recessive genes in our cohort indicates that some individuals with muscular dystrophy continue to face incomplete genetic diagnoses; further refinements in genetic knowledge and diagnostic approaches will optimize diagnostic information for these individuals.

BackgroundLynch syndrome is a hereditary cancer syndrome associated with high risks of colorectal and endometrial cancer that is caused by pathogenic variants in the mismatch repair genes (MLH1, MSH2, MSH6, PMS2, EPCAM). Accurate classification of variants identified in these genes as pathogenic or benign enables informed medical management decisions. Previously, we developed a clinical History Weighting Algorithm (HWA) for the classification of variants of uncertain significance (VUSs) in BRCA1 and BRCA2. The BRCA1/2 HWA is based on the premise that pathogenic variants in these genes will be identified more often in individuals with strong personal and/or family histories of breast and/or ovarian cancer, while the identification of benign variants should be independent of cancer history. Here we report the development of a similar HWA to allow for classification of VUSs in genes associated with Lynch syndrome using data collected through both syndrome-specific and pan-cancer panel testing.MethodsUpon completion of algorithm development, the HWA was tested using simulated variants constructed from 79,214 probands, as well as 379 true variants. Positive (PPV) and negative predictive values (NPV) were calculated on a per gene basis.Results25,500 pathogenic and 50,500 benign simulated variants were analyzed using the HWA and the PPVs and NPVs for each gene were greater than 0.997 and 0.999, respectively. The HWA was also evaluated using 100 trials for each of the 379 true variants. PPVs of >0.998 and NPVs of >0.999 were obtained for all genes.ConclusionsWe have developed and implemented a HWA to aid in the classification of VUSs in genes associated with Lynch syndrome. The work presented here demonstrates that this HWA is able to classify MLH1, MSH2, and MSH6 VUSs as either benign or pathogenic with high accuracy.