How to Identify Pathogenic Mutations among All Those Variations: Variant Annotation and Filtration in the Genome Sequencing Era

An Integrated Pipeline for Trio-Rapid Genome Sequencing in Critically Ill Infants

Trio-rapid genome sequencing (trio-rGS) can assist the genetic diagnosis of critically ill infants given its ability to detect a broad range of pathogenic variants, as well as microbes, simultaneously with high efficiency. To achieve more comprehensive clinical diagnoses, it is essential to propose a recommended protocol in clinical practice. Here, we introduced an integrated pipeline to detect germline variants and microorganisms simultaneously from trio-RGS in critically ill infants, which provides step-by-step criteria for the semi-automatic processing procedures. With this pipeline in clinical application, only 1 ml of peripheral blood is needed for clinicians to provide both genetic and infectious causal information to a patient. The establishment and clinical practice of the method is of great significance for further mining of high-throughput sequencing data and for assisting clinicians in promoting diagnosis efficiency and accuracy. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Experimental pipeline for rapid whole-genome sequencing for the simultaneous detection of germline variants and microorganisms Basic Protocol 2: Computational pipeline for rapid whole-genome sequencing for the simultaneous detection of germline variants and microorganisms.

Genetic association study detected misalignment in previous whole exome sequence: association study of ZNF806 and SART3 in tardive dystonia

Tardive dystonia is one of the most serious types of extrapyramidal symptoms that antipsychotics can cause. There is no established treatment to relieve this symptom, and its etiology is unclear. Recently, we identified very rare single-nucleotide polymorphisms (SNPs) on ZNF806 and SART3 by exome sequencing in three patients with profoundly severe tardive dystonia. Here, we conducted an association study (case, N = 16 vs. control, N = 96) on the rarest SNP selected from each gene. The results showed that rs2287546 on SART3 was not related to tardive dystonia and that rs4953961 on ZNF806 was a heterozygote in all the subjects, implying the absence of a rare SNP in this locus. We found three other genomic regions with high similarity to the relevant region on ZNF806 by BLAT searches. This strongly suggested a misalignment error in this region in our previous exome sequence. In conclusion, ZNF806 and SART3 are unlikely to be related to tardive dystonia.

Inborn errors of immunity-recent advances in research on the pathogenesis

Primary immunodeficiency (PID) is a genetic disorder with a defect of one of the important components of our immune system. Classical PID has been recognized as a disorder with loss of function of the immune system. Recent studies have unveiled disorders with immune dysfunction with autoimmunity, autoinflammation, allergy, or predisposition to malignancy. Some of them were caused by an augmented immune function or a defect in immune regulation. With this background, the term inborn errors of immunity (IEI) is now used to refer to PID in the International Union of Immunological Societies (IUIS) classification. More than 400 responsible genes have been identified in patients with IEI so far, and importantly, many of them identified lately were caused by a heterologous mutation. Moreover, the onset is not necessarily in childhood, and we started seeing more and more IEI patients diagnosed in adulthood in the clinical settings. Recent advances in genetic analysis, including whole-exome analysis, whole-genome analysis, and RNA-seq have contributed to the identification of the disease-causing gene mutation. We also started to find heterogeneity of phenotype even in the patients with the same mutation in the same family, leading us to wonder if modifier gene or epigenetic modification is involved in the pathogenesis. In contrast, we accumulated many cases suggesting genetic heterogeneity is associated with phenotypic homogeneity. It has thus become difficult to deduce a responsible gene only from the phenotype in a certain type of IEI. Current curative therapy for IEI includes hematopoietic cell transplantation and gene therapy. Other curative therapeutic modalities have been long waited and are to be introduced in the future. These include a small molecule that inhibits the gain-of-function of the molecule- and genome-editing technology. Research on IEI will surely lead to a better understanding of other immune-related disorders including rheumatic diseases and atopic disorders.

Genetic Counseling and Genome Sequencing in Pediatric Rare Disease

Both genome sequencing (GS) and exome sequencing (ES) have proven to be revolutionary in the diagnosis of pediatric rare disease. The diagnostic potential and increasing affordability make GS and ES more accessible as a routine clinical test in some centers. Herein, I review aspects of rare disease in pediatrics associated with the use of genomic technologies with an emphasis on the benefits and limitations of both ES and GS, complexities of variant classification, and the importance of genetic counseling. Indications for testing, the role of genetic counselors in genomic test selection, and the diagnostic potential of ES and GS in various pediatric multisystem disorders are discussed. The neonatal population represents an important cohort in pediatric rare disease. Rapid ES and GS in critically ill neonates can have an immediate impact on medical management and present unique genetic counseling challenges. This work includes reviews of recommendations for genetic counseling for families considering genome-wide sequencing, and issues of access to genetic counseling that affect clinical use and will necessitate implementation of innovative methods such as online decision aids. Finally, this work will also review the challenges of having a child with a rare disease, the impact of results from ES and GS on these families, and the role of various support agencies.

Genetic Modifiers and Rare Mendelian Disease

Despite advances in high-throughput sequencing that have revolutionized the discovery of gene defects in rare Mendelian diseases, there are still gaps in translating individual genome variation to observed phenotypic outcomes. While we continue to improve genomics approaches to identify primary disease-causing variants, it is evident that no genetic variant acts alone. In other words, some other variants in the genome (genetic modifiers) may alleviate (suppress) or exacerbate (enhance) the severity of the disease, resulting in the variability of phenotypic outcomes. Thus, to truly understand the disease, we need to consider how the disease-causing variants interact with the rest of the genome in an individual. Here, we review the current state-of-the-field in the identification of genetic modifiers in rare Mendelian diseases and discuss the potential for future approaches that could bridge the existing gap.

Epilepsy and developmental disorders: Next generation sequencing in the clinic

BACKGROUND

The advent of Next Generation Sequencing (NGS) has led to a redefining of the genetic landscape of the epilepsies. Hundreds of single gene epilepsies have been described. Genes associated with epilepsy involve diverse processes. Now a substantial proportion of individuals with epilepsy can receive a high definition molecular genetic diagnosis.

METHODS

In this review we update the current genetic landscape of the epilepsies and categorise the major functional groupings of epilepsy-associated genes. We describe currently available genetic testing approaches. We perform a literature review of NGS studies and review the factors which determine yield in cohorts undergoing testing. We identify factors associated with positive genetic diagnosis and consider the utility of genetic testing in terms of treatment selection as well as more qualitative aspects of care.

FINDINGS

Epilepsy-associated genes can be grouped into five broad functional categories: ion transport; cell growth and differentiation; regulation of synaptic processes; transport and metabolism of small molecules within and between cells; and regulation of gene transcription and translation. Early onset of seizures, drug-resistance, and developmental comorbidity are associated with higher diagnostic yield. The most commonly implicated genes in NGS studies to date, in order, are SCN1A, KCNQ2, CDKL5, SCN2A, and STXBP1. In unselected infantile cohorts PRRT2, a gene associated with self-limited epilepsy, is frequently implicated. Genetic diagnosis provides utility in terms of treatment choice closing the diagnostic odyssey, avoiding unnecessary further testing, and informing future reproductive decisions.

CONCLUSIONS

Genetic testing has become a first line test in epilepsy. As techniques improve and understanding advances, its utility is set to increase. Genetic diagnosis, particularly in early onset developmental and epileptic encephalopathies, influences treatment choice in a significant proportion of patients. The realistic prospect of gene therapy is a cause for optimism.

Genetic Disorders in Prenatal Onset Syndromic Short Stature Identified by Exome Sequencing

OBJECTIVE

To perform a prospective genetic investigation using whole exome sequencing of a group of patients with syndromic short stature born small for gestational age of unknown cause.

STUDY DESIGN

For whole exome sequencing analysis, we selected 44 children born small for gestational age with persistent short stature, and additional features, such as dysmorphic face, major malformation, developmental delay, and/or intellectual disability. Seven patients had negative candidate gene testing based on clinical suspicion and 37 patients had syndromic conditions of unknown etiology.

RESULTS

Of the 44 patients, 15 (34%) had pathogenic/likely pathogenic variants in genes already associated with growth disturbance: COL2A1 (n = 2), SRCAP (n = 2), AFF4, ACTG1, ANKRD11, BCL11B, BRCA1, CDKN1C, GINS1, INPP5K, KIF11, KMT2A, and POC1A (n = 1 each). Most of the genes found to be deleterious participate in fundamental cellular processes, such as cell replication and DNA repair.

CONCLUSIONS

The rarity and heterogeneity of syndromic short stature make the clinical diagnosis difficult. Whole exome sequencing allows the diagnosis of previously undiagnosed patients with syndromic short stature.

Clinical sequencing: From raw data to diagnosis with lifetime value

High-throughput sequencing (HTS) has revolutionized genetics by enabling the detection of sequence variants at hitherto unprecedented large scale. Despite these advances, however, there are still remaining challenges in the complete coverage of targeted regions (genes, exome or genome) as well as in HTS data analysis and interpretation. Moreover, it is easy to get overwhelmed by the plethora of available methods and tools for HTS. Here, we review the step-by-step process from the generation of sequence data to molecular diagnosis of Mendelian diseases. Highlighting advantages and limitations, this review addresses the current state of (1) HTS technologies, considering targeted, whole-exome, and whole-genome sequencing on short- and long-read platforms; (2) read alignment, variant calling and interpretation; as well as (3) regulatory issues related to genetic counseling, reimbursement, and data storage.

VarAFT: a variant annotation and filtration system for human next generation sequencing data

With the rapidly developing high-throughput sequencing technologies known as next generation sequencing or NGS, our approach to gene hunting and diagnosis has drastically changed. In <10 years, these technologies have moved from gene panel to whole genome sequencing and from an exclusively research context to clinical practice. Today, the limit is not the sequencing of one, many or all genes but rather the data analysis. Consequently, the challenge is to rapidly and efficiently identify disease-causing mutations within millions of variants. To do so, we developed the VarAFT software to annotate and pinpoint human disease-causing mutations through access to multiple layers of information. VarAFT was designed both for research and clinical contexts and is accessible to all scientists, regardless of bioinformatics training. Data from multiple samples may be combined to address all Mendelian inheritance modes, cancers or population genetics. Optimized filtration parameters can be stored and re-applied to large datasets. In addition to classical annotations from dbNSFP, VarAFT contains unique features at the disease (OMIM), phenotypic (HPO), gene (Gene Ontology, pathways) and variation levels (predictions from UMD-Predictor and Human Splicing Finder) that can be combined to optimally select candidate pathogenic mutations. VarAFT is freely available at: http://varaft.eu.

Whole exome sequencing of men with multiple morphological abnormalities of the sperm flagella reveals novel homozygous QRICH2 mutations

Multiple morphological anomalies of the sperm flagella (MMAF syndrome) is a severe male infertility phenotype which has so far been formally linked to the presence of biallelic mutations in nine genes mainly coding for axonemal proteins overexpressed in the sperm flagellum. Homozygous mutations in QRICH2, a gene coding for a protein known to be required for stabilizing proteins involved in sperm flagellum biogenesis, have recently been identified in MMAF patients from two Chinese consanguineous families. Here, in order to better assess the contribution of QRICH2 in the etiology of the MMAF phenotype, we analyzed all QRICH2 variants from whole exome sequencing data of a cohort of 167 MMAF-affected subjects originating from North Africa, Iran, and Europe. We identified a total of 14 potentially deleterious variants in 18 unrelated individuals. Two unrelated subjects, representing 1% of the cohort, carried a homozygous loss-of-function variant: c.3501C>G [p.Tyr1167Ter] and c.4614C>G [p.Tyr1538Ter], thus confirming the implication of QRICH2 in the MMAF phenotype and human male infertility. Sixteen MMAF patients (9.6%) carried a heterozygous QRICH2 potentially deleterious variant. This rate was comparable to what was observed in a control group (15.5%) suggesting that the presence of QRICH2 heterozygous variants is not associated with MMAF syndrome.

NQO1: A target for the treatment of cancer and neurological diseases, and a model to understand loss of function disease mechanisms

NAD(P)H quinone oxidoreductase 1 (NQO1) is a multi-functional protein that catalyses the reduction of quinones (and other molecules), thus playing roles in xenobiotic detoxification and redox balance, and also has roles in stabilising apoptosis regulators such as p53. The structure and enzymology of NQO1 is well-characterised, showing a substituted enzyme mechanism in which NAD(P)H binds first and reduces an FAD cofactor in the active site, assisted by a charge relay system involving Tyr-155 and His-161. Protein dynamics play important role in physio-pathological aspects of this protein. NQO1 is a good target to treat cancer due to its overexpression in cancer cells. A polymorphic form of NQO1 (p.P187S) is associated with increased cancer risk and certain neurological disorders (such as multiple sclerosis and Alzheimer´s disease), possibly due to its roles in the antioxidant defence. p.P187S has greatly reduced FAD affinity and stability, due to destabilization of the flavin binding site and the C-terminal domain, which leading to reduced activity and enhanced degradation. Suppressor mutations partially restore the activity of p.P187S by local stabilization of these regions, and showing long-range allosteric communication within the protein. Consequently, the correction of NQO1 misfolding by pharmacological chaperones is a viable strategy, which may be useful to treat cancer and some neurological conditions, targeting structural spots linked to specific disease-mechanisms. Thus, NQO1 emerges as a good model to investigate loss of function mechanisms in genetic diseases as well as to improve strategies to discriminate between neutral and pathogenic variants in genome-wide sequencing studies.

Insight into the specificity and severity of pathogenic mechanisms associated with missense mutations through experimental and structural perturbation analyses

Most pathogenic missense mutations cause specific molecular phenotypes through protein destabilization. However, how protein destabilization is manifested as a given molecular phenotype is not well understood. We develop here a structural and energetic approach to describe mutational effects on specific traits such as function, regulation, stability, subcellular targeting or aggregation propensity. This approach is tested using large-scale experimental and structural perturbation analyses in over thirty mutations in three different proteins (cancer-associated NQO1, transthyretin related with amyloidosis and AGT linked to primary hyperoxaluria type I) and comprising five very common pathogenic mechanisms (loss-of-function and gain-of-toxic function aggregation, enzyme inactivation, protein mistargeting and accelerated degradation). Our results revealed that the magnitude of destabilizing effects and, particularly, their propagation through the structure to promote disease-associated conformational states largely determine the severity and molecular mechanisms of disease-associated missense mutations. Modulation of the structural perturbation at a mutated site is also shown to cause switches between different molecular phenotypes. When very common disease-associated missense mutations were investigated, we also found that they were not among the most deleterious possible missense mutations at those sites, and required additional contributions from codon bias and effects of CpG sites to explain their high frequency in patients. Our work sheds light on the molecular basis of pathogenic mechanisms and genotype–phenotype relationships, with implications for discriminating between pathogenic and neutral changes within human genome variability from whole genome sequencing studies.

Biophysical and functional perturbation analyses at cancer-associated P187 and K240 sites of the multifunctional NADP(H):quinone oxidoreductase 1

Once whole-genome sequencing has reached the clinical practice, a main challenge ahead is the high-throughput and accurate prediction of the pathogenicity of genetic variants. However, current prediction tools do not consider explicitly a well-known property of disease-causing mutations: their ability to affect multiple functional sites distant in the protein structure. Here we carried out an extensive biophysical characterization of fourteen mutant variants at two cancer-associated sites of the enzyme NQO1, a paradigm of multi-functional protein. We showed that the magnitude of destabilizing effects, their molecular origins (structural vs. dynamic) and their efficient propagation through the protein structure gradually led to functional perturbations at different sites. Modulation of these structural perturbations also led to switches between molecular phenotypes. Our work supports that experimental and computational perturbation analyses would improve our understanding of the molecular basis of many loss-of-function genetic diseases as well as our ability to accurately predict the pathogenicity of genetic variants in a high-throughput fashion.

CRIMEtoYHU: a new web tool to develop yeast-based functional assays for characterizing cancer-associated missense variants

Evaluation of the functional impact of cancer-associated missense variants is more difficult than for protein-truncating mutations and consequently standard guidelines for the interpretation of sequence variants have been recently proposed. A number of algorithms and software products were developed to predict the impact of cancer-associated missense mutations on protein structure and function. Importantly, direct assessment of the variants using high-throughput functional assays using simple genetic systems can help in speeding up the functional evaluation of newly identified cancer-associated variants. We developed the web tool CRIMEtoYHU (CTY) to help geneticists in the evaluation of the functional impact of cancer-associated missense variants. Humans and the yeast Saccharomyces cerevisiae share thousands of protein-coding genes although they have diverged for a billion years. Therefore, yeast humanization can be helpful in deciphering the functional consequences of human genetic variants found in cancer and give information on the pathogenicity of missense variants. To humanize specific positions within yeast genes, human and yeast genes have to share functional homology. If a mutation in a specific residue is associated with a particular phenotype in humans, a similar substitution in the yeast counterpart may reveal its effect at the organism level. CTY simultaneously finds yeast homologous genes, identifies the corresponding variants and determines the transferability of human variants to yeast counterparts by assigning a reliability score (RS) that may be predictive for the validity of a functional assay. CTY analyzes newly identified mutations or retrieves mutations reported in the COSMIC database, provides information about the functional conservation between yeast and human and shows the mutation distribution in human genes. CTY analyzes also newly found mutations and aborts when no yeast homologue is found. Then, on the basis of the protein domain localization and functional conservation between yeast and human, the selected variants are ranked by the RS. The RS is assigned by an algorithm that computes functional data, type of mutation, chemistry of amino acid substitution and the degree of mutation transferability between human and yeast protein. Mutations giving a positive RS are highly transferable to yeast and, therefore, yeast functional assays will be more predictable. To validate the web application, we have analyzed 8078 cancer-associated variants located in 31 genes that have a yeast homologue. More than 50% of variants are transferable to yeast. Incidentally, 88% of all transferable mutations have a reliability score >0. Moreover, we analyzed by CTY 72 functionally validated missense variants located in yeast genes at positions corresponding to the human cancer-associated variants. All these variants gave a positive RS. To further validate CTY, we analyzed 3949 protein variants (with positive RS) by the predictive algorithm PROVEAN. This analysis shows that yeast-based functional assays will be more predictable for the variants with positive RS. We believe that CTY could be an important resource for the cancer research community by providing information concerning the functional impact of specific mutations, as well as for the design of functional assays useful for decision support in precision medicine.

Rare disease registries: a call to action

When registries collect accurate clinical data over time, they can act as fundamental support structures for patients and their families and powerful cost-effective instruments to support clinical trials and translational research to improve quality of care, quality of life and survival. Registries are critical for rare diseases (RD) with low prevalence and propensity for variation in treatment and outcomes. Rare Voices Australia is leading a call for action to the research and clinical community to prioritise RD data collection and develop an integrated RD Registry strategy for Australia. Financial, operational and governance challenges exist for establishing and maintaining RD registries. As a multidisciplinary team whose interests converge on RD, we highlight the need for the establishment of an Australian RD Registry Alliance. This 'umbrella' organisation will: (i) bring together existing RD registries across Australia; (ii) establish National RD Registry Standards to support interoperability and cohesion across registries; (iii) develop strategies to attract sustainable funding from government and other sources to maximise the utility of existing RD registries and support the development of new RD registries. The most important role for the Alliance would be to use the RD registries for translational research to address current knowledge gaps about RD and to improve the care for the over 1.4 million Australians estimated to live with RD.

Paediatric genomics: diagnosing rare disease in children

The majority of rare diseases affect children, most of whom have an underlying genetic cause for their condition. However, making a molecular diagnosis with current technologies and knowledge is often still a challenge. Paediatric genomics is an immature but rapidly evolving field that tackles this issue by incorporating next-generation sequencing technologies, especially whole-exome sequencing and whole-genome sequencing, into research and clinical workflows. This complex multidisciplinary approach, coupled with the increasing availability of population genetic variation data, has already resulted in an increased discovery rate of causative genes and in improved diagnosis of rare paediatric disease. Importantly, for affected families, a better understanding of the genetic basis of rare disease translates to more accurate prognosis, management, surveillance and genetic advice; stimulates research into new therapies; and enables provision of better support.

Scuba: scalable kernel-based gene prioritization

BACKGROUND

The uncovering of genes linked to human diseases is a pressing challenge in molecular biology and precision medicine. This task is often hindered by the large number of candidate genes and by the heterogeneity of the available information. Computational methods for the prioritization of candidate genes can help to cope with these problems. In particular, kernel-based methods are a powerful resource for the integration of heterogeneous biological knowledge, however, their practical implementation is often precluded by their limited scalability.

RESULTS

We propose Scuba, a scalable kernel-based method for gene prioritization. It implements a novel multiple kernel learning approach, based on a semi-supervised perspective and on the optimization of the margin distribution. Scuba is optimized to cope with strongly unbalanced settings where known disease genes are few and large scale predictions are required. Importantly, it is able to efficiently deal both with a large amount of candidate genes and with an arbitrary number of data sources. As a direct consequence of scalability, Scuba integrates also a new efficient strategy to select optimal kernel parameters for each data source. We performed cross-validation experiments and simulated a realistic usage setting, showing that Scuba outperforms a wide range of state-of-the-art methods.

CONCLUSIONS

Scuba achieves state-of-the-art performance and has enhanced scalability compared to existing kernel-based approaches for genomic data. This method can be useful to prioritize candidate genes, particularly when their number is large or when input data is highly heterogeneous. The code is freely available at https://github.com/gzampieri/Scuba .

Targeted Genetic Screen in Amyotrophic Lateral Sclerosis Reveals Novel Genetic Variants with Synergistic Effect on Clinical Phenotype

Amyotrophic lateral sclerosis (ALS) is underpinned by an oligogenic rare variant architecture. Identified genetic variants of ALS include RNA-binding proteins containing prion-like domains (PrLDs). We hypothesized that screening genes encoding additional similar proteins will yield novel genetic causes of ALS. The most common genetic variant of ALS patients is a G4C2-repeat expansion within C9ORF72. We have shown that G4C2-repeat RNA sequesters RNA-binding proteins. A logical consequence of this is that loss-of-function mutations in G4C2-binding partners might contribute to ALS pathogenesis independently of and/or synergistically with C9ORF72 expansions. Targeted sequencing of genomic DNA encoding either RNA-binding proteins or known ALS genes (n = 274 genes) was performed in ALS patients to identify rare deleterious genetic variants and explore genotype-phenotype relationships. Genomic DNA was extracted from 103 ALS patients including 42 familial ALS patients and 61 young-onset (average age of onset 41 years) sporadic ALS patients; patients were chosen to maximize the probability of identifying genetic causes of ALS. Thirteen patients carried a G4C2-repeat expansion of C9ORF72. We identified 42 patients with rare deleterious variants; 6 patients carried more than one variant. Twelve mutations were discovered in known ALS genes which served as a validation of our strategy. Rare deleterious variants in RNA-binding proteins were significantly enriched in ALS patients compared to control frequencies (p = 5.31E-18). Nineteen patients featured at least one variant in a RNA-binding protein containing a PrLD. The number of variants per patient correlated with rate of disease progression (t-test, p = 0.033). We identified eighteen patients with a single variant in a G4C2-repeat binding protein. Patients with a G4C2-binding protein variant in combination with a C9ORF72 expansion had a significantly faster disease course (t-test, p = 0.025). Our data are consistent with an oligogenic model of ALS. We provide evidence for a number of entirely novel genetic variants of ALS caused by mutations in RNA-binding proteins. Moreover we show that these mutations act synergistically with each other and with C9ORF72 expansions to modify the clinical phenotype of ALS. A key finding is that this synergy is present only between functionally interacting variants. This work has significant implications for ALS therapy development.

The Impact of Whole-Genome Sequencing on the Primary Care and Outcomes of Healthy Adult Patients: A Pilot Randomized Trial

Background

Whole-genome sequencing (WGS) in asymptomatic adults might prevent disease but increase health care use without clinical value.

Objective

To describe the effect on clinical care and outcomes of adding WGS to standardized family history assessment in primary care.

Design

Pilot randomized trial. (ClinicalTrials.gov: NCT01736566).

Setting

Academic primary care practices.

Participants

9 primary care physicians (PCPs) and 100 generally healthy patients recruited at ages 40 to 65 years.

Intervention

Patients were randomly assigned to receive a family history report alone (FH group) or in combination with an interpreted WGS report (FH + WGS group), which included monogenic disease risk (MDR) results (associated with Mendelian disorders), carrier variants, pharmacogenomic associations, and polygenic risk estimates for cardiometabolic traits. Each patient met with his or her PCP to discuss the report.

Measurements

Clinical outcomes and health care use through 6 months were obtained from medical records and audio-recorded discussions between PCPs and patients. Patients' health behavior changes were surveyed 6 months after receiving results. A panel of clinician-geneticists rated the appropriateness of how PCPs managed MDR results.

Results

Mean age was 55 years; 58% of patients were female. Eleven FH + WGS patients (22% [95% CI, 12% to 36%]) had new MDR results. Only 2 (4% [CI, 0.01% to 15%]) had evidence of the phenotypes predicted by an MDR result (fundus albipunctatus due to RDH5 and variegate porphyria due to PPOX). Primary care physicians recommended new clinical actions for 16% (CI, 8% to 30%) of FH patients and 34% (CI, 22% to 49%) of FH + WGS patients. Thirty percent (CI, 17% to 45%) and 41% (CI, 27% to 56%) of FH and FH + WGS patients, respectively, reported making a health behavior change after 6 months. Geneticists rated PCP management of 8 MDR results (73% [CI, 39% to 99%]) as appropriate and 2 results (18% [CI, 3% to 52%]) as inappropriate.

Limitation

Limited sample size and ancestral and socioeconomic diversity.

Conclusion

Adding WGS to primary care reveals new molecular findings of uncertain clinical utility. Nongeneticist providers may be able to manage WGS results appropriately, but WGS may prompt additional clinical actions of unclear value.

Primary Funding Source

National Institutes of Health.

Coupling between Protein Stability and Catalytic Activity Determines Pathogenicity of G6PD Variants

G6PD deficiency, an enzymopathy affecting 7% of the world population, is caused by over 160 identified amino acid variants in glucose-6-phosphate dehydrogenase (G6PD). The clinical presentation of G6PD deficiency is diverse, likely due to the broad distribution of variants across the protein and the potential for multidimensional biochemical effects. In this study, we use bioinformatic and biochemical analyses to interpret the relationship between G6PD variants and their clinical phenotype. Using structural information and statistical analyses of known G6PD variants, we predict the molecular phenotype of five uncharacterized variants from a reference population database. Through multidimensional analysis of biochemical data, we demonstrate that the clinical phenotypes of G6PD variants are largely determined by a trade-off between protein stability and catalytic activity. This work expands the current understanding of the biochemical underpinnings of G6PD variant pathogenicity and suggests a promising avenue for correcting G6PD deficiency by targeting essential structural features of G6PD.

Actionable Genes, Core Databases, and Locus-Specific Databases

Adoption of next-generation sequencing (NGS) in a diagnostic context raises numerous questions with regard to identification and reports of secondary variants (SVs) in actionable genes. To better understand the whys and wherefores of these questioning, it is necessary to understand how they are selected during the filtering process and how their proportion can be estimated. It is likely that SVs are underestimated and that our capacity to label all true SVs can be improved. In this context, Locus-specific databases (LSDBs) can be key by providing a wealth of information and enabling classifying variants. We illustrate this issue by analyzing 318 SVs in 23 actionable genes involved in cancer susceptibility syndromes identified through sequencing of 572 participants selected for a range of atherosclerosis phenotypes. Among these 318 SVs, only 43.4% are reported in Human Gene Mutation Database (HGMD) Professional versus 71.4% in LSDB. In addition, 23.9% of HGMD Professional variants are reported as pathogenic versus 4.8% for LSDB. These data underline the benefits of LSDBs to annotate SVs and minimize overinterpretation of mutations thanks to their efficient curation process and collection of unpublished data.