VariantValidator: Accurate validation, mapping, and formatting of sequence variation descriptions

Improving access to exome sequencing in a medically underserved population through the Texome Project

PURPOSE

Genomic medicine can end diagnostic odysseys for patients with complex phenotypes; however, limitations in insurance coverage and other systemic barriers preclude individuals from accessing comprehensive genetics evaluation and testing.

METHODS

The Texome Project is a 4-year study that reduces barriers to genomic testing for individuals from underserved and underrepresented populations. Participants with undiagnosed, rare diseases who have financial barriers to obtaining exome sequencing (ES) clinically are enrolled in the Texome Project.

RESULTS

We highlight the Texome Project process and describe the outcomes of the first 60 ES results for study participants. Participants received a genetic evaluation, ES, and return of results at no cost. We summarize the psychosocial or medical implications of these genetic diagnoses. Thus far, ES provided molecular diagnoses for 18 out of 60 (30%) of Texome participants. Plus, in 11 out of 60 (18%) participants, a partial or probable diagnosis was identified. Overall, 5 participants had a change in medical management.

CONCLUSION

To date, the Texome Project has recruited a racially, ethnically, and socioeconomically diverse cohort. The diagnostic rate and medical impact in this cohort support the need for expanded access to genetic testing and services. The Texome Project will continue reducing barriers to genomic care throughout the future study years.

Genotype and phenotype correlation of PHACTR1-related neurological disorders

BACKGROUND

PHACTR1 (phosphatase and actin regulators) plays a key role in cortical migration and synaptic activity by binding and regulating G-actin and PPP1CA. This study aimed to expand the genotype and phenotype of patients with de novo variants in PHACTR1 and analyse the impact of variants on protein-protein interaction.

METHODS

We identified seven patients with PHACTR1 variants by trio-based whole-exome sequencing. Additional two subjects were ascertained from two centres through GeneMatcher. The genotype-phenotype correlation was determined, and AlphaFold-Multimer was used to predict protein-protein interactions and interfaces.

RESULTS

Eight individuals carried missense variants and one had CNV in the PHACTR1. Infantile epileptic spasms syndrome (IESS) was the unifying phenotype in eight patients with missense variants of PHACTR1. They could present with other types of seizures and often exhibit drug-resistant epilepsy with a poor prognosis. One patient with CNV displayed a developmental encephalopathy phenotype. Using AlphaFold-Multimer, our findings indicate that PHACTR1 and G-actin-binding sequences overlap with PPP1CA at the RPEL3 domain, which suggests possible competition between PPP1CA and G-actin for binding to PHACTR1 through a similar polymerisation interface. In addition, patients carrying missense variants located at the PHACTR1-PPP1CA or PHACTR1-G-actin interfaces consistently exhibit the IESS phenotype. These missense variants are mostly concentrated in the overlapping sequence (RPEL3 domain).

CONCLUSIONS

Patients with variants in PHACTR1 can have a phenotype of developmental encephalopathy in addition to IESS. Moreover, our study confirmed that the variants affect the binding of PHACTR1 to G-actin or PPP1CA, resulting in neurological disorders in patients.

REEV: review, evaluate and explain variants

In the era of high throughput sequencing, special software is required for the clinical evaluation of genetic variants. We developed REEV (Review, Evaluate and Explain Variants), a user-friendly platform for clinicians and researchers in the field of rare disease genetics. Supporting data was aggregated from public data sources. We compared REEV with seven other tools for clinical variant evaluation. REEV (semi-)automatically fills individual ACMG criteria facilitating variant interpretation. REEV can store disease and phenotype data related to a case to use these for phenotype similarity measures. Users can create public permanent links for individual variants that can be saved as browser bookmarks and shared. REEV may help in the fast diagnostic assessment of genetic variants in a clinical as well as in a research context. REEV (https://reev.bihealth.org/) is free and open to all users and there is no login requirement.

Comparative analysis of in-silico tools in identifying pathogenic variants in dominant inherited retinal diseases

Inherited retinal diseases (IRDs) are a group of rare genetic eye conditions that cause blindness. Despite progress in identifying genes associated with IRDs, improvements are necessary for classifying rare autosomal dominant (AD) disorders. AD diseases are highly heterogenous, with causal variants being restricted to specific amino acid changes within certain protein domains, making AD conditions difficult to classify. Here, we aim to determine the top-performing in-silico tools for predicting the pathogenicity of AD IRD variants. We annotated variants from ClinVar and benchmarked 39 variant classifier tools on IRD genes, split by inheritance pattern. Using area-under-the-curve (AUC) analysis, we determined the top-performing tools and defined thresholds for variant pathogenicity. Top-performing tools were assessed using genome sequencing on a cohort of participants with IRDs of unknown etiology. MutScore achieved the highest accuracy within AD genes, yielding an AUC of 0.969. When filtering for AD gain-of-function and dominant negative variants, BayesDel had the highest accuracy with an AUC of 0.997. Five participants with variants in NR2E3, RHO, GUCA1A, and GUCY2D were confirmed to have dominantly inherited disease based on pedigree, phenotype, and segregation analysis. We identified two uncharacterized variants in GUCA1A (c.428T>A, p.Ile143Thr) and RHO (c.631C>G, p.His211Asp) in three participants. Our findings support using a multi-classifier approach comprised of new missense classifier tools to identify pathogenic variants in participants with AD IRDs. Our results provide a foundation for improved genetic diagnosis for people with IRDs.

RASopathies are the most common set of monogenic syndromes identified by exome sequencing for nonimmune hydrops fetalis: A systematic review and meta-analysis

RASopathies are a group of malformation syndromes known to lead to nonimmune hydrops fetalis (NIHF) in severe presentations. Pathogenic variants can be de novo or parentally inherited. Despite being a known frequent presentation, the fraction of monogenic NIHF cases due to RASopathies is limited in the literature. Also, the specific parental contribution of RASopathies to NIHF is not well described. Our objective was to review pooled exome sequencing (ES) diagnostic yield of RASopathies for NIHF and to determine the parental contribution of RASopathy to NIHF. We performed a systematic review of prenatal ES studies from January 1, 2000 to August 1, 2022. Thirty-six studies met inclusion criteria. Cases with RASopathy gene variants were reviewed. NIHF cases were further classified as isolated or non-isolated. Thirty-six ES studies including 46 pregnancies with NIHF and a diagnosed RASopathy were reviewed. Forty-four diagnostic variants and 2 variants of uncertain significance in 12 RASopathy genes were identified. Expanding on what was previously published, a total of 506 NIHF cases were extracted with 191 cases yielding a positive diagnosis by ES. The overall rate of RASopathy diagnosis in clinically diagnosed NIHF cases was 9% (44/506). The rate of RASopathy diagnosis among NIHF cases with positive genetic diagnosis by ES was 23% (44/191). Of the 46 cases identified, 13 (28%) variants were parentally inherited; specifically, 5/13 (38%) maternal, 3/13 (23%) paternal, 2/13 (15%) biparental, and 3/13 (23%) unspecified. Majority of NIHF cases 29/46 (63%) were isolated. Among NIHF cases with positive ES diagnoses, RASopathy diagnostic yield by ES was 23%. NIHF secondary to RASopathies was parentally inherited in 28% of cases. Most cases of NIHF due to RASopathy were isolated, with no prenatal detection of associated anomalies.

Mosaic RASopathies concept: different skin lesions, same systemic manifestations?

BACKGROUND

Cutaneous epidermal nevi are genotypically diverse mosaic disorders. Pathogenic hotspot variants in HRAS, KRAS, and less frequently, NRAS and BRAF may cause isolated keratinocytic epidermal nevi and sebaceous nevi or several different syndromes when associated with extracutaneous anomalies. Therefore, some authors suggest the concept of mosaic RASopathies to group these different disorders.

METHODS

In this paper, we describe three new cases of syndromic epidermal nevi caused by mosaic HRAS variants: one associating an extensive keratinocytic epidermal nevus with hypomastia, another with extensive mucosal involvement and a third combining a small sebaceous nevus with seizures and intellectual deficiency. Moreover, we performed extensive literature of all cases of syndromic epidermal nevi and related disorders with confirmed pathogenic postzygotic variants in HRAS, KRAS, NRAS or BRAF.

RESULTS

Most patients presented with bone, ophthalmological or neurological anomalies. Rhabdomyosarcoma, urothelial cell carcinoma and pubertas praecox are also repeatedly reported. KRAS pathogenic variants are involved in 50% of the cases, especially in sebaceous nevi, oculoectodermal syndrome and encephalocraniocutaneous lipomatosis. They are frequently associated with eye and brain anomalies. Pathogenic variants in HRAS are rather present in syndromic keratinocytic epidermal nevi and phacomatosis pigmentokeratotica.

CONCLUSION

This review delineates genotype/phenotype correlations of syndromic epidermal nevi with somatic RAS and BRAF pathogenic variants and may help improve their follow-up.

Phenotypic characterisation of SMAD4 variant carriers

BACKGROUND

Both hereditary haemorrhagic telangiectasia (HHT) and juvenile polyposis syndrome (JPS) are known to be caused by SMAD4 pathogenic variants, with overlapping symptoms for both disorders in some patients. Additional connective tissue disorders have also been reported. Here, we describe carriers of SMAD4 variants followed in an HHT reference centre to further delineate the phenotype.

METHODS

Observational study based on data collected from the Clinical Investigation for the Rendu-Osler Cohort database.

RESULTS

Thirty-three participants from 15 families, out of 1114 patients with HHT, had an SMAD4 variant (3%).Regarding HHT, 26 out of 33 participants (88%) had a definite clinical diagnosis based on Curaçao criteria. Complication frequencies were as follows: epistaxis (n=27/33, 82%), cutaneous telangiectases (n=19/33, 58%), pulmonary arteriovenous malformations (n=17/32, 53%), hepatic arteriovenous malformations (AVMs) (n=7/18, 39%), digestive angiodysplasia (n=13/22, 59%). No cerebral AVMs were diagnosed.Regarding juvenile polyposis, 25 out of 31 participants (81%) met the criteria defined by Jass et al for juvenile polyposis syndrome. Seven patients (21%) had a prophylactic gastrectomy due to an extensive gastric polyposis incompatible with endoscopic follow-up, and four patients (13%) developed a digestive cancer.Regarding connective tissue disorders, 20 (61%) had at least one symptom, and 4 (15%) participants who underwent echocardiography had an aortic dilation.

CONCLUSION

We describe a large cohort of SMAD4 variant carriers in the context of HHT. Digestive complications are frequent, early and diffuse, justifying endoscopy every 2 years. The HHT phenotype, associating pulmonary and hepatic AVMs, warrants systematic screening. Connective tissue disorders broaden the phenotype associated with SMAD4 gene variants and justify systematic cardiac ultrasound and skeletal complications screening.

Molecular testing in newborn screening: VUS burden among true positives and secondary reproductive limitations via expanded carrier screening panels

PURPOSE

Expanded carrier screening (ECS) gene panels have several limitations, including variable content, current knowledge of disease-causing variants, and differing reporting policies. This study evaluated if the disease-associated variants identified in affected neonates who screened positive by California newborn screening (NBS) for an inherited metabolic disorder (IMD) by tandem mass spectrometry (MS/MS) would likely be reported by ECS gene panels.

METHODS

Retrospective review of neonates referred by the California Department of Public Health for a positive NBS by multianalyte MS/MS from January 1, 2020 through June 30, 2021.

RESULTS

One hundred thirty-six neonates screened positive for ≥1 NBS MS/MS indication. Nineteen neonates (14%) were ultimately diagnosed with an IMD, all of whom had abnormal biochemical testing. Eighteen of the 19 underwent molecular testing; 10 (56%) neonates had ≥1 variants of uncertain significance, 9 of whom were of non-White ancestry. ECS panels would have been negative for 56% (20/36) of parents with an affected neonate, 85% (17/20) of whom were of non-White ancestry.

CONCLUSION

The number of variants of uncertain significance identified in this cohort highlights the need for more diversified variant databases. Due in part to the lack of diversity in currently sequenced populations, genomic sequencing cannot replace biochemical testing for the diagnosis of an IMD.

A novel mutation in the NR3C1 gene associated with reversible glucocorticoid resistance

OBJECTIVE

Glucocorticoid resistance is a rare endocrine disease caused by variants of the NR3C1 gene encoding the glucocorticoid receptor (GR). We identified a novel heterozygous variant (GRR569Q) in a patient with uncommon reversible glucocorticoid resistance syndrome.

METHODS

We performed ex vivo functional characterization of the variant in patient fibroblasts and in vitro through transient transfection in undifferentiated HEK 293T cells to assess transcriptional activity, affinity, and nuclear translocation. We studied the impact of the variant on the tertiary structure of the ligand-binding domain through 3D modeling.

RESULTS

The patient presented initially with an adrenal adenoma with mild autonomous cortisol secretion and undetectable adrenocorticotropin hormone (ACTH) levels. Six months after surgery, biological investigations showed elevated cortisol and ACTH (urinary free cortisol 114 µg/24 h, ACTH 10.9 pmol/L) without clinical symptoms, evoking glucocorticoid resistance syndrome. Functional characterization of the GRR569Q showed decreased expression of target genes (in response to 100 nM cortisol: SGK1 control +97% vs patient +20%, P < .0001) and impaired nuclear translocation in patient fibroblasts compared to control. Similar observations were made in transiently transfected cells, but higher cortisol concentrations overcame glucocorticoid resistance. GRR569Q showed lower ligand affinity (Kd GRWT: 1.73 nM vs GRR569Q: 4.61 nM). Tertiary structure modeling suggested a loss of hydrogen bonds between H3 and the H1-H3 loop.

CONCLUSION

This is the first description of a reversible glucocorticoid resistance syndrome with effective negative feedback on corticotroph cells regarding increased plasma cortisol concentrations due to the development of mild autonomous cortisol secretion.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders

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 pedigrees 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.

Short-read whole genome sequencing identifies causative variants in most individuals with previously unexplained aniridia

BACKGROUND

Classic aniridia is a highly penetrant autosomal dominant disorder characterised by congenital absence of the iris, foveal hypoplasia, optic disc anomalies and progressive opacification of the cornea. >90% of cases of classic aniridia are caused by heterozygous, loss-of-function variants affecting the PAX6 locus.

METHODS

Short-read whole genome sequencing was performed on 51 (39 affected) individuals from 37 different families who had screened negative for mutations in the PAX6 coding region.

RESULTS

Likely causative mutations were identified in 22 out of 37 (59%) families. In 19 out of 22 families, the causative genomic changes have an interpretable deleterious impact on the PAX6 locus. Of these 19 families, 1 has a novel heterozygous PAX6 frameshift variant missed on previous screens, 4 have single nucleotide variants (SNVs) (one novel) affecting essential splice sites of PAX6 5' non-coding exons and 2 have deep intronic SNV (one novel) resulting in gain of a donor splice site. In 12 out of 19, the causative variants are large-scale structural variants; 5 have partial or whole gene deletions of PAX6, 3 have deletions encompassing critical PAX6 cis-regulatory elements, 2 have balanced inversions with disruptive breakpoints within the PAX6 locus and 2 have complex rearrangements disrupting PAX6. The remaining 3 of 22 families have deletions encompassing FOXC1 (a known cause of atypical aniridia). Seven of the causative variants occurred de novo and one cosegregated with familial aniridia. We were unable to establish inheritance status in the remaining probands. No plausibly causative SNVs were identified in PAX6 cis-regulatory elements.

CONCLUSION

Whole genome sequencing proves to be an effective diagnostic test in most individuals with previously unexplained aniridia.

Heterozygosity in factor XIII genes and the manifestation of mild inherited factor XIII deficiency

BACKGROUND

The characterization of inherited mild factor XIII deficiency is more imprecise than its rare, inherited severe forms. It is known that heterozygosity at FXIII genetic loci results in mild FXIII deficiency, characterized by circulating FXIII activity levels ranging from 20% to 60%. There exists a gap in information on 1) how genetic heterozygosity renders clinical bleeding manifestations among these individuals and 2) the reversal of unexplained bleeding upon FXIII administration in mild FXIII-deficient individuals.

OBJECTIVES

To assess the prevalence and burden of mild FXIII deficiency among the apparently healthy German-Caucasian population and correlate it with genetic heterozygosity at FXIII and fibrinogen gene loci.

METHODS

Peripheral blood was collected from 752 donors selected from the general population with essentially no bleeding complications to ensure asymptomatic predisposition. These were assessed for FXIII and fibrinogen activity, and FXIII and fibrinogen genes were resequenced using next-generation sequencing. For comparison, a retrospective analysis was performed on a cohort of mild inherited FXIII deficiency patients referred to us.

RESULTS

The prevalence of mild FXIII deficiency was high (∼0.8%) among the screened German-Caucasian population compared with its rare-severe forms. Although no new heterozygous missense variants were found, certain combinations were relatively dominant/prevalent among the mild FXIII-deficient individuals.

CONCLUSION

This extensive, population-based quasi-experimental approach revealed that the burden of heterozygosity in FXIII and fibrinogen gene loci causes the clinical manifestation of inherited mild FXIII deficiency, resulting in ''unexplained bleeding'' upon provocation.

Lessons learned from the first national population-based genetic carrier-screening program for Duchenne muscular dystrophy

PURPOSE

To summarize the results of first year implementation of pan-ethnic screening testing for Duchenne muscular dystrophy (DMD) and present the ensuing challenges.

METHODS

Data acquisition for this study was performed by retrospective search of Ministry of Health registry for reports of all laboratories performing genetic screening tests. DMD testing was performed by multiplex ligation-dependent probe amplification technology. In case of single-exon deletion, sequencing of the specific exon was performed to rule out underlying single-nucleotide variant.

RESULTS

Of overall 85,737 DMD tests, 82 clinically significant findings were noted (0.095%, or 1:1,046 women). In addition, 80 findings with uncertain clinical significance were detected (0.093%, or 1:1072), as well as 373 cases (0.4%, or 1:230) of single-exon deletions subsequently identified as false positives because of underlying single-nucleotide variant, mostly variants in exon 8 in North African Jewish population, and in exon 48 in Arab Muslim population.

CONCLUSION

Interpretation of population-based DMD carrier screening is complex, occasionally requiring additional genetic testing methods and ethical considerations. Multicenter data registry, including ethnic origin and familial segregation in selected cases, is crucial for optimal definition of the results during genetic counseling and informed decisions regarding prenatal testing.

Clinical outcomes and response to chemotherapy in a cohort of pancreatic cancer patients with germline variants of unknown significance (VUS) in BRCA1 and BRCA2 genes

PURPOSE

The clinical outcome and the efficacy of chemotherapy in pancreatic cancer patients with BRCA1/2 Variants of Unknown Significance (VUS) is unknown. We explored the effects of chemotherapy with or without Platinum in non metastatic and metastatic pancreatic cancer patients with BRCA1/2 VUS.

METHODS

A retrospective analysis of non-metastatic or metastatic pancreatic cancer patients with gBRCA1/2 VUS treated in 13 Italian centers between November 2015 and December 2020 was performed. All patients were assessed for toxicity and RECIST 1.1 response. Metastatic patients were evaluated for survival outcome.

RESULTS

30 pancreatic cancer patients with gBRCA1/2 VUS were considered: 20 were M+ and 10 were non-M+. Pl-CT was recommended to 16 patients: 10 M+ (6 FOLFIRINOX and 4 PAXG) and 6 non-M+ (3 FOLFIRINOX and 3 PAXG); 11 patients received Nabpaclitaxel-Gemcitabine (AG; 8 M+) and 3 patients (2 M+) were treated with Gemcitabine (G). The RECIST 1.1 response rate was 27% for AG and 44% for Pl-CT (22% for (m) FOLFIRINOX and 71% PAXG). 1 year Progression-Free Survival was 37.5% for patients treated with AG and 33% in the Pl-CT subgroup. Median Overall Survival (OS) was 23.5 months for patients treated with AG and 14 months for the Pl-CT subgroup. 1 Year and 2 Year OS were numerically better for AG (1 Year OS: 75% vs 60% and 2 Year OS: 50% and 20% in AG and Pl-CT subgroups, respectively) as well.

CONCLUSIONS

Pl-CT does not seem to be associated with a better outcome compared to AG chemotherapy in PDAC patients with BRCA 1/2 VUS.

Exome sequencing (ES) of a pediatric cohort with chronic endocrine diseases: a single-center study (within the framework of the TRANSLATE-NAMSE project)

BACKGROUND

Endocrine disorders are heterogeneous and include a significant number of rare monogenic diseases.

METHODS

We performed exome sequencing (ES) in 106 children recruited from a single center within the TRANSLATE‑NAMSE project. They were categorized into subgroups: proportionate short stature (PSS), disproportionate short stature (DSS), hypopituitarism (H), differences in sexual development (DSD), syndromic diseases (SD) and others.

RESULTS

The overall diagnostic yield was 34.9% (n = 37/106), including 5 patients with variants in candidate genes, which have contributed to collaborations to identify gene-disease associations. The diagnostic yield varied significantly between subgroups: PSS: 16.6% (1/6); DSS: 18.8% (3/16); H: 17.1% (6/35); DSD: 37.5% (3/8); SD: 66.6% (22/33); others: 25% (2/8). Confirmed diagnoses included 75% ultrarare diseases. Three patients harbored more than one disease-causing variant, resulting in dual diagnoses.

CONCLUSIONS

ES is an effective tool for genetic diagnosis in pediatric patients with complex endocrine diseases. An accurate phenotypic description, including comprehensive endocrine diagnostics, as well as the evaluation of variants in multidisciplinary case conferences involving geneticists, are necessary for personalized diagnostic care. Here, we illustrate the broad spectrum of genetic endocrinopathies that have led to the initiation of specific treatment, surveillance, and family counseling.

Maleic acid is a biomarker for maleylacetoacetate isomerase deficiency; implications for newborn screening of tyrosinemia type 1

Dried blood spot succinylacetone (SA) is often used as a biomarker for newborn screening (NBS) for tyrosinemia type 1 (TT1). However, false-positive SA results are often observed. Elevated SA may also be due to maleylacetoacetate isomerase deficiency (MAAI-D), which appears to be clinically insignificant. This study investigated whether urine organic acid (uOA) and quantitative urine maleic acid (Q-uMA) analyses can distinguish between TT1 and MAAI-D. We reevaluated/measured uOA (GC-MS) and/or Q-uMA (LC-MS/MS) in available urine samples of nine referred newborns (2 TT1, 7 false-positive), eight genetically confirmed MAAI-D children, and 66 controls. Maleic acid was elevated in uOA of 5/7 false-positive newborns and in the three available samples of confirmed MAAI-D children, but not in TT1 patients. Q-uMA ranged from not detectable to 1.16 mmol/mol creatinine in controls (n = 66) and from 0.95 to 192.06 mmol/mol creatinine in false-positive newborns and MAAI-D children (n = 10). MAAI-D was genetically confirmed in 4/7 false-positive newborns, all with elevated Q-uMA, and rejected in the two newborns with normal Q-uMA. No sample was available for genetic analysis of the last false-positive infant with elevated Q-uMA. Our study shows that MAAI-D is a recognizable cause of false-positive TT1 NBS results. Elevated urine maleic acid excretion seems highly effective in discriminating MAAI-D from TT1.

VarSCAT: A computational tool for sequence context annotations of genomic variants

The sequence contexts of genomic variants play important roles in understanding biological significances of variants and potential sequencing related variant calling issues. However, methods for assessing the diverse sequence contexts of genomic variants such as tandem repeats and unambiguous annotations have been limited. Herein, we describe the Variant Sequence Context Annotation Tool (VarSCAT) for annotating the sequence contexts of genomic variants, including breakpoint ambiguities, flanking bases of variants, wildtype/mutated DNA sequences, variant nomenclatures, distances between adjacent variants, tandem repeat regions, and custom annotation with user customizable options. Our analyses demonstrate that VarSCAT is more versatile and customizable than the currently available methods or strategies for annotating variants in short tandem repeat (STR) regions or insertions and deletions (indels) with breakpoint ambiguity. Variant sequence context annotations of high-confidence human variant sets with VarSCAT revealed that more than 75% of all human individual germline and clinically relevant indels have breakpoint ambiguities. Moreover, we illustrate that more than 80% of human individual germline small variants in STR regions are indels and that the sizes of these indels correlated with STR motif sizes. VarSCAT is available from https://github.com/elolab/VarSCAT.

Poison exon annotations improve the yield of clinically relevant variants in genomic diagnostic testing

PURPOSE

Neurodevelopmental disorders (NDDs) often result from rare genetic variation, but genomic testing yield for NDDs remains below 50%, suggesting that clinically relevant variants may be missed by standard analyses. Here, we analyze "poison exons" (PEs), which are evolutionarily conserved alternative exons often absent from standard gene annotations. Variants that alter PE inclusion can lead to loss of function and may be highly penetrant contributors to disease.

METHODS

We curated published RNA sequencing data from developing mouse cortex to define 1937 conserved PE regions potentially relevant to NDDs, and we analyzed variants found by genome sequencing in multiple NDD cohorts.

RESULTS

Across 2999 probands, we found 6 novel clinically relevant variants in PE regions. Five of these variants are in genes that are part of the sodium voltage-gated channel alpha subunit family (SCN1A, SCN2A, and SCN8A), which is associated with epilepsies. One variant is in SNRPB, associated with cerebrocostomandibular syndrome. These variants have moderate to high computational impact assessments, are absent from population variant databases, and in genes with gene-phenotype associations consistent with each probands reported features.

CONCLUSION

With a very minimal increase in variant analysis burden (average of 0.77 variants per proband), annotation of PEs can improve diagnostic yield for NDDs and likely other congenital conditions.

Pharmacovariome scanning using whole pharmacogene resequencing coupled with deep computational analysis and machine learning for clinical pharmacogenomics

BACKGROUND

This pilot study aims to identify and functionally assess pharmacovariants in whole exome sequencing data. While detection of known variants has benefited from pharmacogenomic-dedicated bioinformatics tools before, in this paper we have tested novel deep computational analysis in addition to artificial intelligence as possible approaches for functional analysis of unknown markers within less studied drug-related genes.

METHODS

Pharmacovariants from 1800 drug-related genes from 100 WES data files underwent (a) deep computational analysis by eight bioinformatic algorithms (overall containing 23 tools) and (b) random forest (RF) classifier as the machine learning (ML) approach separately. ML model efficiency was calculated by internal and external cross-validation during recursive feature elimination. Protein modelling was also performed for predicted highly damaging variants with lower frequencies. Genotype-phenotype correlations were implemented for top selected variants in terms of highest possibility of being damaging.

RESULTS

Five deleterious pharmacovariants in the RYR1, POLG, ANXA11, CCNH, and CDH23 genes identified in step (a) and subsequent analysis displayed high impact on drug-related phenotypes. Also, the utilization of recursive feature elimination achieved a subset of 175 malfunction pharmacovariants in 135 drug-related genes that were used by the RF model with fivefold internal cross-validation, resulting in an area under the curve of 0.9736842 with an average accuracy of 0.9818 (95% CI: 0.89, 0.99) on predicting whether a carrying individuals will develop adverse drug reactions or not. However, the external cross-validation of the same model indicated a possible false positive result when dealing with a low number of observations, as only 60 important variants in 49 genes were displayed, giving an AUC of 0.5384848 with an average accuracy of 0.9512 (95% CI: 0.83, 0.99).

CONCLUSION

While there are some technologies for functionally assess not-interpreted pharmacovariants, there is still an essential need for the development of tools, methods, and algorithms which are able to provide a functional prediction for every single pharmacovariant in both large-scale datasets and small cohorts. Our approaches may bring new insights for choosing the right computational assessment algorithms out of high throughput DNA sequencing data from small cohorts to be used for personalized drug therapy implementation.

Clinical implications of conflicting variant interpretations in the cancer genetics clinic

PURPOSE

The aim of this study was to describe the clinical impact of commercial laboratories issuing conflicting classifications of genetic variants.

METHODS

Results from 2000 patients undergoing a multigene hereditary cancer panel by a single laboratory were analyzed. Clinically significant discrepancies between the laboratory-provided test reports and other major commercial laboratories were identified, including differences between pathogenic/likely pathogenic and variant of uncertain significance (VUS) classifications, via review of ClinVar archives. For patients carrying a VUS, clinical documentation was assessed for evidence of provider awareness of the conflict.

RESULTS

Fifty of 975 (5.1%) patients with non-negative results carried a variant with a clinically significant conflict, 19 with a pathogenic/likely pathogenic variant reported in APC or MUTYH, and 31 with a VUS reported in CDKN2A, CHEK2, MLH1, MSH2, MUTYH, RAD51C, or TP53. Only 10 of 28 (36%) patients with a VUS with a clinically significant conflict had a documented discussion by a provider about the conflict. Discrepant counseling strategies were used for different patients with the same variant. Among patients with a CDKN2A variant or a monoallelic MUTYH variant, providers were significantly more likely to make recommendations based on the laboratory-reported classification.

CONCLUSION

Our findings highlight the frequency of variant interpretation discrepancies and importance of clinician awareness. Guidance is needed on managing patients with discrepant variants to support accurate risk assessment.

Non-oral manifestations in adults with a clinical and molecularly confirmed diagnosis of periodontal Ehlers-Danlos syndrome

Introduction: Periodontal Ehlers-Danlos Syndrome (pEDS) is a rare autosomal dominant type of EDS characterised by severe early-onset periodontitis, lack of attached gingiva, pretibial plaques, joint hypermobility and skin hyperextensibility as per the 2017 International EDS Classification. In 2016, deleterious pathogenic heterozygous variants were identified in C1R and C1S, which encode components of the complement system. Materials and Methods: Individuals with a clinical suspicion of pEDS were clinically and molecularly assessed through the National EDS Service in London and Sheffield and in genetic services in Austria, Sweden and Australia. Transmission electron microscopy and fibroblast studies were performed in a small subset of patients. Results: A total of 21 adults from 12 families were clinically and molecularly diagnosed with pEDS, with C1R variants in all families. The age at molecular diagnosis ranged from 21-73 years (mean 45 years), male: female ratio 5:16. Features of easy bruising (90%), pretibial plaques (81%), skin fragility (71%), joint hypermobility (24%) and vocal changes (38%) were identified as well as leukodystrophy in 89% of those imaged. Discussion: This cohort highlights the clinical features of pEDS in adults and contributes several important additional clinical features as well as novel deleterious variants to current knowledge. Hypothetical pathogenic mechanisms which may help to progress understanding and management of pEDS are also discussed.

Investigation on the high recurrence of the ATTRv-causing transthyretin variant Val142Ile in central Italy

The p.Val142Ile variant in transthyretin (encoded by the TTR gene) is the most common genetic cause of transthyretin-related amyloidosis. This allele is particularly prevalent in communities ofAfrican descent compared with populations of different ancestries, where its frequency is two orders of magnitude lower. For this reason, p.Val142Ile has always been considered an "African" variant, with limited studies performed on individuals of European descent. However, recent reports of higher-than-expected prevalence in European-ancestry populations question the African specificity of this allele. Here we show that the high recurrence of p.Val142Ile in central Italy is due to a founder effect and not to recent admixture from African populations, highlighting how this may be the case in other communities. This suggests a probable underestimate of the global prevalence of p.Val142Ile, and further emphasizes the importance of routine inclusion of TTR in gene panels used for clinical genetic testing in hypertrophic cardiomyopathy (independently of the patient's geographical origin), that transthyretin-related amyloidosis can mimic.

NAGS, CPS1, and SLC25A13 (Citrin) at the Crossroads of Arginine and Pyrimidines Metabolism in Tumor Cells

Urea cycle enzymes and transporters collectively convert ammonia into urea in the liver. Aberrant overexpression of carbamylphosphate synthetase 1 (CPS1) and SLC25A13 (citrin) genes has been associated with faster proliferation of tumor cells due to metabolic reprogramming that increases the activity of the CAD complex and pyrimidine biosynthesis. N-acetylglutamate (NAG), produced by NAG synthase (NAGS), is an essential activator of CPS1. Although NAGS is expressed in lung cancer derived cell lines, expression of the NAGS gene and its product was not evaluated in tumors with aberrant expression of CPS1 and citrin. We used data mining approaches to identify tumor types that exhibit aberrant overexpression of NAGS, CPS1, and citrin genes, and evaluated factors that may contribute to increased expression of the three genes and their products in tumors. Median expression of NAGS, CPS1, and citrin mRNA was higher in glioblastoma multiforme (GBM), glioma, and stomach adenocarcinoma (STAD) samples compared to the matched normal tissue. Median expression of CPS1 and citrin mRNA was higher in the lung adenocarcinoma (LUAD) sample while expression of NAGS mRNA did not differ. High NAGS expression was associated with an unfavorable outcome in patients with glioblastoma and GBM. Low NAGS expression was associated with an unfavorable outcome in patients with LUAD. Patterns of DNase hypersensitive sites and histone modifications in the upstream regulatory regions of NAGS, CPS1, and citrin genes were similar in liver tissue, lung tissue, and A549 lung adenocarcinoma cells despite different expression levels of the three genes in the liver and lung. Citrin gene copy numbers correlated with its mRNA expression in glioblastoma, GBM, LUAD, and STAD samples. There was little overlap between NAGS, CPS1, and citrin sequence variants found in patients with respective deficiencies, tumor samples, and individuals without known rare genetic diseases. The correlation between NAGS, CPS1, and citrin mRNA expression in the individual glioblastoma, GBM, LUAD, and STAD samples was very weak. These results suggest that the increased cytoplasmic supply of either carbamylphosphate, produced by CPS1, or aspartate may be sufficient to promote tumorigenesis, as well as the need for an alternative explanation of CPS1 activity in the absence of NAGS expression and NAG.

The first genetic landscape of inherited retinal dystrophies in Portuguese patients identifies recurrent homozygous mutations as a frequent cause of pathogenesis

Inherited retinal diseases (IRDs) are a group of ocular conditions characterized by an elevated genetic and clinical heterogeneity. They are transmitted almost invariantly as monogenic traits. However, with more than 280 disease genes identified so far, association of clinical phenotypes with genotypes can be very challenging, and molecular diagnosis is essential for genetic counseling and correct management of the disease. In addition, the prevalence and the assortment of IRD mutations are often population-specific. In this work, we examined 230 families from Portugal, with individuals suffering from a variety of IRD diagnostic classes (270 subjects in total). Overall, we identified 157 unique mutations (34 previously unreported) in 57 distinct genes, with a diagnostic rate of 76%. The IRD mutational landscape was, to some extent, different from those reported in other European populations, including Spanish cohorts. For instance, the EYS gene appeared to be the most frequently mutated, with a prevalence of 10% among all IRD cases. This was, in part, due to the presence of a recurrent and seemingly founder mutation involving the deletion of exons 13 and 14 of this gene. Moreover, our analysis highlighted that as many as 51% of our cases had mutations in a homozygous state. To our knowledge, this is the first study assessing a cross-sectional genotype-phenotype landscape of IRDs in Portugal. Our data reveal a rather unique distribution of mutations, possibly shaped by a small number of rare ancestral events that have now become prevalent alleles in patients.

Diagnostic yield from prenatal exome sequencing for non-immune hydrops fetalis: A systematic review and meta-analysis

Non-immune hydrops fetalis (NIHF) has multiple genetic etiologies diagnosable by exome sequencing (ES). We evaluated the yield of prenatal ES for NIHF, and the contribution of additional clinical findings and history. Systematic review was performed with PROSPERO tag 232951 using CINAHL, PubMed, and Ovid MEDLINE from January 1, 2000 through December 1, 2021. Selected studies performed ES to augment standard prenatal diagnostic approaches. Cases meeting a strict NIHF phenotype were tabulated with structured data imputed from papers or requested from authors. Genetic variants and diagnostic outcomes were harmonized across studies using current ACMG and ClinGen variant classification guidelines. Thirty-one studies reporting 445 NIHF cases had a 37% (95% CI: 32%-41%) diagnostic rate. There was no significant difference between isolated NIHF and NIHF with fetal malformations or between recurrent and simplex cases. Diagnostic rate was higher for consanguineous than non-consanguineous cases. Disease categories included RASopathies (24%), neuromuscular (21%), metabolic (17%), lymphatic (13%), other syndromes (9%), cardiovascular (5%), hematologic (2%), skeletal (2%), and other categories (7%). Inheritance patterns included recessive (55%), dominant (41%), and X-linked (4%). ES should be considered in the diagnostic workup of NIHF with and without associated ultrasound findings regardless of history of recurrence or consanguinity.

A systematic review on Treacher Collins syndrome: Correlation between molecular genetic findings and clinical severity

Treacher Collins syndrome (TCS, OMIM: 154500) is a rare congenital craniofacial disorder that is caused by variants in the genes TCOF1, POLR1D, POLR1C, and POLR1B. Studies on the association between phenotypic variability and their relative variants are very limited. This systematic review summarized the 53 literatures from PubMed and Scopus to explore the potential TCS genotype-phenotype correlations with statistical analysis. Studies reporting both complete molecular genetics and clinical data were included. We identified that the molecular anomaly within TCOF1 (88.71%) accounted for most TCS cases. The only true hot spot for TCOF1 was detected in exon 24, with recurrent c.4369_4373delAAGAA variant is identified. While the hot spot for POLR1D, POLR1C, and POLR1B were identified in exons 3, 8, and 15, respectively. Our result suggested that the higher severity level was likely to be observed in Asian patients harboring TCOF1 variants rather than POLR1. Moreover, common 5-bp deletions tended to have a higher severity degree in comparison to any variants within exon 24 of TCOF1. In summary, this report suggested the relationship between genetic and clinical data in TCS. Our findings could be used as a reference for clinical diagnosis and further biological studies.

Systematic errors in annotations of truncations, loss-of-function and synonymous variants

Description of genetic phenomena and variations requires exact language and concepts. Vast amounts of variation data are produced with next-generation sequencing pipelines. The obtained variations are automatically annotated, e.g., for their functional consequences. These tools and pipelines, along with systematic nomenclature, mainly work well, but there are still some problems in nomenclature, organization of some databases, misuse of concepts and certain practices. Therefore, systematic errors prevent correct annotation and often preclude further analysis of certain variation types. Problems and solutions are described for presumed protein truncations, variants that are claimed to be of loss-of-function based on the type of variation, and synonymous variants that are not synonymous and lead to sequence changes or to missing protein.

Poison exon annotations improve the yield of clinically relevant variants in genomic diagnostic testing

Purpose

Neurodevelopmental disorders (NDDs) often result from rare genetic variation, but genomic testing yield for NDDs remains around 50%, suggesting some clinically relevant rare variants may be missed by standard analyses. Here we analyze "poison exons" (PEs) which, while often absent from standard gene annotations, are alternative exons whose inclusion results in a premature termination codon. Variants that alter PE inclusion can lead to loss-of-function and may be highly penetrant contributors to disease.

Methods

We curated published RNA-seq data from developing mouse cortex to define 1,937 PE regions conserved between humans and mice and potentially relevant to NDDs. We then analyzed variants found by genome sequencing in multiple NDD cohorts.

Results

Across 2,999 probands, we found six clinically relevant variants in PE regions that were previously overlooked. Five of these variants are in genes that are part of the sodium voltage-gated channel alpha subunit family ( SCN1A, SCN2A , and SCN8A ), associated with epilepsies. One variant is in SNRPB , associated with Cerebrocostomandibular Syndrome. These variants have moderate to high computational impact assessments, are absent from population variant databases, and were observed in probands with features consistent with those reported for the associated gene.

Conclusion

With only a minimal increase in variant analysis burden (most probands had zero or one candidate PE variants in a known NDD gene, with an average of 0.77 per proband), annotation of PEs can improve diagnostic yield for NDDs and likely other congenital conditions.

The global prevalence and ethnic heterogeneity of iron-refractory iron deficiency anaemia

BACKGROUND

Iron-refractory iron deficiency anaemia (IRIDA) is an autosomal recessive iron deficiency anaemia caused by mutations in the TMPRSS6 gene. Iron deficiency anaemia is common, whereas IRIDA is rare. The prevalence of IRIDA is unclear. This study aimed to estimate the carrier frequency and genetic prevalence of IRIDA using Genome Aggregation Database (gnomAD) data.

METHODS

The pathogenicity of TMPRSS6 variants was interpreted according to the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) standards and guidelines. The minor allele frequency (MAF) of TMPRSS6 gene disease-causing variants in 141,456 unique individuals was examined to estimate the global prevalence of IRIDA in seven ethnicities: African/African American (afr), American Admixed/Latino (amr), Ashkenazi Jewish (asj), East Asian (eas), Finnish (fin), Non-Finnish European (nfe) and South Asian (sas). The global and population-specific carrier frequencies and genetic prevalence of IRIDA were calculated using the Hardy-Weinberg equation.

RESULTS

In total, 86 pathogenic/likely pathogenic variants (PV/LPV) were identified according to ACMG/AMP guideline. The global carrier frequency and genetic prevalence of IRIDA were 2.02 per thousand and 1.02 per million, respectively.

CONCLUSIONS

The prevalence of IRIDA is greater than previous estimates.

Enrichment of titin-truncating variants in exon 327 in dilated cardiomyopathy and its relevance to reduced nonsense-mediated mRNA decay efficiency

Titin truncating variants (TTNtvs) are the most common genetic cause of dilated cardiomyopathy (DCM). Among four regions of titin, A-band enrichment of DCM-causing TTNtvs is widely accepted but the underlying mechanism is still unknown. Meanwhile, few reports have identified exon 327 as a highly mutated A-band exon but the degree of exon 327 enrichment has not been quantitatively investigated. To find the real hotspot of DCM-causing TTNtvs, we aimed to reassess the degree of TTNtv enrichment in known titin regions and in exon 327, separately. In addition, we tried to explain exon 327 clustering in terms of nonsense-mediated mRNA decay (NMD) efficiency and a dominant negative mechanism recently proposed. Research papers focusing on TTNtvs found in patients with DCM were collected. A total of 612 patients with TTNtv-realated DCM were obtained from 10 studies. In the four regions of TTN and exon 327, the degree of TTNtvs enrichment was calculated in a way that the effect of distribution of highly expressed exons was normalized. As a result, exon 327 was the only region that showed significant enrichment for DCM-related TTNtv (p < .001). On the other hand, other A-band exons had almost the same number of TTNtv of random distribution. A review of RNAseq data revealed that the median allelic imbalance deviation of exon 327 TTNtvs was .04, indicating almost zero NMD. From these findings, we propose that the widely accepted A-band enrichment of DCM-related TTNtv is mostly attributable to exon 327 enrichment. In addition, based on the recently demonstrated dominant negative mechanism, the extremely low NMD efficiency seems to contribute to exon 327 enrichment.

A Boolean algebra for genetic variants

MOTIVATION

Beyond identifying genetic variants, we introduce a set of Boolean relations, which allows for a comprehensive classification of the relations of every pair of variants by taking all minimal alignments into account. We present an efficient algorithm to compute these relations, including a novel way of efficiently computing all minimal alignments within the best theoretical complexity bounds.

RESULTS

We show that these relations are common, and many non-trivial, for variants of the CFTR gene in dbSNP. Ultimately, we present an approach for the storing and indexing of variants in the context of a database that enables efficient querying for all these relations.

AVAILABILITY AND IMPLEMENTATION

A Python implementation is available at https://github.com/mutalyzer/algebra/tree/v0.2.0 as well as an interface at https://mutalyzer.nl/algebra.

A spectrum of BRCA1 and BRCA2 germline deleterious variants in ovarian cancer in Russia

PURPOSE

Pathogenic variants (PVs) in BRCA1 and BRCA2 genes are essential biomarkers of an increased breast and ovarian cancer risk and tumor sensitivity to poly ADP ribose polymerase inhibitors. In Russia, eight PVs were thought to be the most common, among which BRCA1 c.5266dup is the most frequently identified one.

METHODS

We show the distribution of BRCA1/2 PVs identified with quantitative PCR and targeted next-generation sequencing in 1399 ovarian cancer patients recruited into the study from 72 Russian regions in 2015-2021.

RESULTS

The most abundant PVs were c.5266dup (41.0%), c.4035del (7.0%), c.1961del (6.3%), c.181 T > G (5.2%), c.3756_3759del (1.8%), c.3700_3704del (1.5%), and c.68_69del (1.5%), all found in BRCA1 and known to be recurrent in Russia. Several other frequent PVs were identified: c.5152 + 1G > T (1.2%), c.1687C > T (1.0%), c.4689C > G (0.9%), c.1510del (0.6%), c.2285_2286del (0.6%) in the BRCA1 gene; and c.5286 T > G (1.2%), c.2808_2811del (0.8%), c.3847_3848del (0.8%), c.658_659del (0.7%), c.7879A > T (0.6%), in the BRCA2 gene. For the most common PV in the BRCA2 gene c.5286 T > G, we suggested that it arose about 700 years ago and is a new founder mutation.

CONCLUSION

This study extends our knowledge about the BRCA1 and BRCA2 pathogenic variants variability.

Biallelic variants in the calpain regulatory subunit CAPNS1 cause pulmonary arterial hypertension

Purpose

The aim of this study was to identify the monogenic cause of pulmonary arterial hypertension (PAH), a multifactorial and often fatal disease, in 2 unrelated consanguine families.

Methods

We performed exome sequencing and validated variant pathogenicity by whole-blood RNA and protein expression analysis in both families. Further RNA sequencing of preserved lung tissue was performed to investigate the consequences on selected genes that are involved in angiogenesis, proliferation, and apoptosis.

Results

We identified 2 rare biallelic variants in CAPNS1, encoding the regulatory subunit of calpain. The variants cosegregated with PAH in the families. Both variants lead to loss of function (LoF), which is demonstrated by aberrant splicing resulting in the complete absence of the CAPNS1 protein in affected patients. No other LoF CAPNS1 variant was identified in the genome data of more than 1000 patients with unresolved PAH.

Conclusion

The calpain holoenzyme was previously linked to pulmonary vascular development and progression of PAH in patients. We demonstrated that biallelic LoF variants in CAPNS1 can cause idiopathic PAH by the complete absence of CAPNS1 protein. Screening of this gene in patients who are affected by PAH, especially with suspected autosomal recessive inheritance, should be considered.

Large scale genotype- and phenotype-driven machine learning in Von Hippel-Lindau disease

Von Hippel-Lindau (VHL) disease is a hereditary cancer syndrome where individuals are predisposed to tumor development in the brain, adrenal gland, kidney, and other organs. It is caused by pathogenic variants in the VHL tumor suppressor gene. Standardized disease information has been difficult to collect due to the rarity and diversity of VHL patients. Over 4100 unique articles published until October 2019 were screened for germline genotype-phenotype data. Patient data were translated into standardized descriptions using Human Genome Variation Society gene variant nomenclature and Human Phenotype Ontology terms and has been manually curated into an open-access knowledgebase called Clinical Interpretation of Variants in Cancer. In total, 634 unique VHL variants, 2882 patients, and 1991 families from 427 papers were captured. We identified relationship trends between phenotype and genotype data using classic statistical methods and spectral clustering unsupervised learning. Our analyses reveal earlier onset of pheochromocytoma/paraganglioma and retinal angiomas, phenotype co-occurrences and genotype-phenotype correlations including hotspots. It confirms existing VHL associations and can be used to identify new patterns and associations in VHL disease. Our database serves as an aggregate knowledge translation tool to facilitate sharing information about the pathogenicity of VHL variants.

When a Synonymous Variant Is Nonsynonymous

Term synonymous variation is widely used, but frequently in a wrong or misleading meaning and context. Twenty three point eight % of possible nucleotide substitution types in the universal genetic code are for synonymous amino acid changes, but when these variants have a phenotype and functional effect, they are very seldom synonymous. Such variants may manifest changes at DNA, RNA and/or protein levels. Large numbers of variations are erroneously annotated as synonymous, which causes problems e.g., in clinical genetics and diagnosis of diseases. To facilitate precise communication, novel systematics and nomenclature are introduced for variants that when looking only at the genetic code seem like synonymous, but which have phenotypes. A new term, unsense variant is defined as a substitution in the mRNA coding region that affects gene expression and protein production without introducing a stop codon in the variation site. Such variants are common and need to be correctly annotated. Proper naming and annotation are important also to increase awareness of these variants and their consequences.

Biallelic variants in WARS1 cause a highly variable neurodevelopmental syndrome and implicate a critical exon for normal auditory function

Aminoacyl-tRNA synthetases (ARSs) are essential enzymes for faithful assignment of amino acids to their cognate tRNA. Variants in ARS genes are frequently associated with clinically heterogeneous phenotypes in humans and follow both autosomal dominant or recessive inheritance patterns in many instances. Variants in tryptophanyl-tRNA synthetase 1 (WARS1) cause autosomal dominantly inherited distal hereditary motor neuropathy and Charcot-Marie-Tooth disease. Presently, only one family with biallelic WARS1 variants has been described. We present three affected individuals from two families with biallelic variants (p.Met1? and p.(Asp419Asn)) in WARS1, showing varying severities of developmental delay and intellectual disability. Hearing impairment and microcephaly, as well as abnormalities of the brain, skeletal system, movement/gait, and behavior were variable features. Phenotyping of knocked down wars-1 in a C. elegans model showed depletion is associated with defects in germ cell development. A wars1 knockout vertebrate model recapitulates the human clinical phenotypes, confirms variant pathogenicity and uncovers evidence implicating the p.Met1? variant as potentially impacting an exon critical for normal hearing. Together, our findings provide consolidating evidence for biallelic disruption of WARS1 as causal for an autosomal recessive neurodevelopmental syndrome and present a vertebrate model that recapitulates key phenotypes observed in patients. This article is protected by copyright. All rights reserved.

Evaluating breast cancer predisposition genes in women of African ancestry

PURPOSE

Studies conducted primarily among European ancestry women reported 12 breast cancer predisposition genes. However, etiologic roles of these genes in breast cancer among African ancestry women have been less well-investigated.

METHODS

We conducted a case-control study in African American women, which included 1117 breast cancer cases and 2169 cancer-free controls, and a pooled analysis, which included 7096 cases and 8040 controls of African descent. Odds ratios of associations with breast cancer risk were estimated.

RESULTS

Using sequence data, we identified 61 pathogenic variants in 12 breast cancer predisposition genes, including 11 pathogenic variants not yet reported in previous studies. Pooled analysis showed statistically significant associations of breast cancer risk with pathogenic variants in BRCA1, BRCA2, PALB2, ATM, CHEK2, TP53, NF1, RAD51C, and RAD51D (all P < .05). The associations with BRCA1, PALB2, and RAD51D were stronger for estrogen receptor (ER)-negative than for ER-positive breast cancer (P heterogeneity < .05), whereas the association with CHEK2 was stronger for ER-positive than for ER-negative breast cancer.

CONCLUSION

Our study confirmed previously identified associations of breast cancer risk with BRCA1, BRCA2, PALB2, ATM, TP53, NF1, and CHEK2 and provided new evidence to extend the associations of breast cancer risk with RAD51C and RAD51D, which was identified previously in European ancestry populations, to African ancestry women.

A novel phenotype associated with the R162W variant in the KCNJ13 gene

BACKGROUND

Pathogenic variants in KCNJ13 have been associated with both autosomal dominant Snowflake vitreoretinal degeneration (SVD) and autosomal recessive Leber congenital amaurosis. SVD is characterized by aberrant vitreoretinal interface leading to increased risk of retinal detachment, crystalline retinal snowflake deposits, optic disc abnormalities, early-onset cataract, and cornea guttae. Reduced dark adaptation and reduced scotopic rod b-waves have also been described. We report a novel phenotype associated with the R162W variant in KCNJ13.

METHODS

Four affected members of a Swedish family were included. Three of them were examined with best corrected visual acuity, Goldmann perimetry, full-field-and multifocal electroretinography, optical coherence tomography, fundus color photographs, fundus autofluorescence images, slit lamp inspection, and genetic testing. The fourth subject only managed genetic testing.

RESULTS

All subjects carry the pathogenic missense variant; c.484C>T (NM_002242.4), R162W, in KCNJ13. ERG measurements revealed reduced macular-as well as general retinal function. Two of the subjects had a history of retinal detachment and the two younger subjects demonstrated early onset cataract. They all had structural macular changes and slightly gliotic optic discs.

CONCLUSION

In this family, the R162W variant in KCNJ13, previously described in association with SVD, causes a somewhat novel phenotype including macular dystrophy and moderate reduction of general retinal function as the main features combined with disc abnormalities, retinal detachment, and presenile cataract that has been described before. In times of up-coming gene-based therapies, it is important to report new genotype-phenotype associations to improve the possibilities to identify future treatment candidates.

The clinical features of polymerase proof-reading associated polyposis (PPAP) and recommendations for patient management

Pathogenic germline exonuclease domain (ED) variants of POLE and POLD1 cause the Mendelian dominant condition polymerase proof-reading associated polyposis (PPAP). We aimed to describe the clinical features of all PPAP patients with probably pathogenic variants. We identified patients with a variants mapping to the EDs of POLE or POLD1 from cancer genetics clinics, a colorectal cancer (CRC) clinical trial, and systematic review of the literature. We used multiple evidence sources to separate ED variants into those with strong evidence of pathogenicity and those of uncertain importance. We performed quantitative analysis of the risk of CRC, colorectal adenomas, endometrial cancer or any cancer in the former group. 132 individuals carried a probably pathogenic ED variant (105 POLE, 27 POLD1). The earliest malignancy was colorectal cancer at 14. The most common tumour types were colorectal, followed by endometrial in POLD1 heterozygotes and duodenal in POLE heterozygotes. POLD1-mutant cases were at a significantly higher risk of endometrial cancer than POLE heterozygotes. Five individuals with a POLE pathogenic variant, but none with a POLD1 pathogenic variant, developed ovarian cancer. Nine patients with POLE pathogenic variants and one with a POLD1 pathogenic variant developed brain tumours. Our data provide important evidence for PPAP management. Colonoscopic surveillance is recommended from age 14 and upper-gastrointestinal surveillance from age 25. The management of other tumour risks remains uncertain, but surveillance should be considered. In the absence of strong genotype-phenotype associations, these recommendations should apply to all PPAP patients.

Evaluation of in silico predictors on short nucleotide variants in HBA1, HBA2, and HBB associated with haemoglobinopathies

Haemoglobinopathies are the commonest monogenic diseases worldwide and are caused by variants in the globin gene clusters. With over 2400 variants detected to date, their interpretation using the American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) guidelines is challenging and computational evidence can provide valuable input about their functional annotation. While many in silico predictors have already been developed, their performance varies for different genes and diseases. In this study, we evaluate 31 in silico predictors using a dataset of 1627 variants in HBA1, HBA2, and HBB. By varying the decision threshold for each tool, we analyse their performance (a) as binary classifiers of pathogenicity and (b) by using different non-overlapping pathogenic and benign thresholds for their optimal use in the ACMG/AMP framework. Our results show that CADD, Eigen-PC, and REVEL are the overall top performers, with the former reaching moderate strength level for pathogenic prediction. Eigen-PC and REVEL achieve the highest accuracies for missense variants, while CADD is also a reliable predictor of non-missense variants. Moreover, SpliceAI is the top performing splicing predictor, reaching strong level of evidence, while GERP++ and phyloP are the most accurate conservation tools. This study provides evidence about the optimal use of computational tools in globin gene clusters under the ACMG/AMP framework.

Genome sequencing as a first-line diagnostic test for hospitalized infants

PURPOSE

SouthSeq is a translational research study that undertook genome sequencing (GS) for infants with symptoms suggestive of a genetic disorder. Recruitment targeted racial/ethnic minorities and rural, medically underserved areas in the Southeastern United States, which are historically underrepresented in genomic medicine research.

METHODS

GS and analysis were performed for 367 infants to detect disease-causal variation concurrent with standard of care evaluation and testing.

RESULTS

Definitive diagnostic (DD) or likely diagnostic (LD) genetic findings were identified in 30% of infants, and 14% of infants harbored an uncertain result. Only 43% of DD/LD findings were identified via concurrent clinical genetic testing, suggesting that GS testing is better for obtaining early genetic diagnosis. We also identified phenotypes that correlate with the likelihood of receiving a DD/LD finding, such as craniofacial, ophthalmologic, auditory, skin, and hair abnormalities. We did not observe any differences in diagnostic rates between racial/ethnic groups.

CONCLUSION

We describe one of the largest-to-date GS cohorts of ill infants, enriched for African American and rural patients. Our results show the utility of GS because it provides early-in-life detection of clinically relevant genetic variations not detected by current clinical genetic testing, particularly for infants exhibiting certain phenotypic features.

BDV Syndrome: An Emerging Syndrome With Profound Obesity and Neurodevelopmental Delay Resembling Prader-Willi Syndrome

CONTEXT

CPE encodes carboxypeptidase E, an enzyme that converts proneuropeptides and propeptide hormones to bioactive forms. It is widely expressed in the endocrine and central nervous system. To date, 4 individuals from 2 families with core clinical features including morbid obesity, neurodevelopmental delay, and hypogonadotropic hypogonadism, harboring biallelic loss-of-function (LoF) CPE variants, have been reported.

OBJECTIVE

We describe 4 affected individuals from 3 unrelated consanguineous families, 2 siblings of Syrian, 1 of Egyptian, and 1 of Pakistani descent, all harboring novel homozygous CPE LoF variants.

METHODS

After excluding Prader-Willi syndrome (PWS), exome sequencing was performed in both Syrian siblings. The variants identified in the other 2 individuals were reported as research variants in a large-scale exome study and in the ClinVar database. Computational modeling of all possible missense alterations allowed assessing CPE tolerance to missense variants.

RESULTS

All affected individuals were severely obese with neurodevelopmental delay and other endocrine anomalies. Three individuals from 2 families shared the same CPE homozygous truncating variant c.361C > T, p.(Arg121*), while the fourth carried the c.994del, p.(Ser333Alafs*22) variant. Comparison of clinical features with previously described cases and standardization according to the Human Phenotype Ontology terms indicated a recognizable clinical phenotype, which we termed Blakemore-Durmaz-Vasileiou (BDV) syndrome. Computational analysis indicated high conservation of CPE domains and intolerance to missense changes.

CONCLUSION

Biallelic truncating CPE variants are associated with BDV syndrome, a clinically recognizable monogenic recessive syndrome with childhood-onset obesity, neurodevelopmental delay, hypogonadotropic hypogonadism, and hypothyroidism. BDV syndrome resembles PWS. Our findings suggest missense variants may also be clinically relevant.

Comparison of in silico strategies to prioritize rare genomic variants impacting RNA splicing for the diagnosis of genomic disorders

The development of computational methods to assess pathogenicity of pre-messenger RNA splicing variants is critical for diagnosis of human disease. We assessed the capability of eight algorithms, and a consensus approach, to prioritize 249 variants of uncertain significance (VUSs) that underwent splicing functional analyses. The capability of algorithms to differentiate VUSs away from the immediate splice site as being 'pathogenic' or 'benign' is likely to have substantial impact on diagnostic testing. We show that SpliceAI is the best single strategy in this regard, but that combined usage of tools using a weighted approach can increase accuracy further. We incorporated prioritization strategies alongside diagnostic testing for rare disorders. We show that 15% of 2783 referred individuals carry rare variants expected to impact splicing that were not initially identified as 'pathogenic' or 'likely pathogenic'; one in five of these cases could lead to new or refined diagnoses.

Variations in Nomenclature of Clinical Variants between Annotation Tools

BACKGROUND

Accurate nomenclature of variants is an essential element for genetic diagnosis and patient care.

OBJECTIVE

To investigate annotation differences of clinical variants between annotation tools.

METHODS

We analyzed 218,156 clinical variants from the Human Gene Mutation Database. Multiple nomenclatures based on RefSeq transcripts were provided using ANNOVAR and snpEff.

RESULTS

The concordance rate between ANNOVAR and snpEff was approximately 85%. Based on the Human Genome Variation Society (HGVS) nomenclature, snpEff was more accurate than ANNOVAR (coding variants, 99.3% vs 84.9%; protein variants, 94.3% vs 79.8%). When annotating each variant with ANNOVAR and snpEff, the accuracy of nomenclature was 99.5%.

CONCLUSIONS

There were substantial differences between ANNOVAR and snpEff annotations. The findings of this study suggest that simultaneous use of multiple annotation tools could decrease nomenclature errors and contribute to providing standardized clinical reporting.

Whole exome sequencing in 17 consanguineous Iranian pedigrees expands the mutational spectrum of inherited retinal dystrophies

Inherited retinal dystrophies (IRDs) constitute one of the most heterogeneous groups of Mendelian human disorders. Using autozygome-guided next-generation sequencing methods in 17 consanguineous pedigrees of Iranian descent with isolated or syndromic IRD, we identified 17 distinct genomic variants in 11 previously-reported disease genes. Consistent with a recessive inheritance pattern, as suggested by pedigrees, variants discovered in our study were exclusively bi-allelic and mostly in a homozygous state (in 15 families out of 17, or 88%). Out of the 17 variants identified, 5 (29%) were never reported before. Interestingly, two mutations (GUCY2D:c.564dup, p.Ala189ArgfsTer130 and TULP1:c.1199G > A, p.Arg400Gln) were also identified in four separate pedigrees (two pedigrees each). In addition to expanding the mutational spectrum of IRDs, our findings confirm that the traditional practice of endogamy in the Iranian population is a prime cause for the appearance of IRDs.

A monoallelic SEC23A variant E599K associated with cranio-lenticulo-sutural dysplasia

Cranio-lenticulo-sutural dysplasia (CLSD; MIM 607812) is a rare or underdiagnosed condition, as only two families have been reported. The original family (Boyadjiev et al., Human Genetics, 2003, 113, 1-9 and Boyadjiev et al., Nature Genetics, 2006, 38, 1192-1197) showed recessive inheritance of the condition with a biallelic SEC23A missense variant in affected individuals. In contrast, another child with sporadic CLSD had a monoallelic SEC23A variant inherited from the reportedly unaffected father (Boyadjiev et al., Clinical Genetics, 2011, 80, 169-176), raising questions on possible digenism. Here, we report a 2-month-old boy seen because of large fontanels with wide cranial sutures, a large forehead, hypertelorism, a thin nose, a high arched palate, and micrognathia. His mother was clinically unremarkable, while his father had a history of large fontanels in infancy who had closed only around age 10 years; he also had a large forehead, hypertelorism, a thin, beaked nose and was operated for bilateral glaucoma with exfoliation of the lens capsule. Trio genome sequencing and familial segregation revealed a monoallelic c.1795G > A transition in SEC23A that was de novo in the father and transmitted to the proband. The variant predicts a nonconservative substitution (p.E599K) in an ultra-conserved residue that is seen in 3D models of yeast SEC23 to be involved in direct binding between SEC23 and SAR1 subunits of the coat protein complex II coat. This observation confirms the link between SEC23A variants and CLSD but suggests that in addition to the recessive inheritance described in the original family, SEC23A variants may result in dominant inheritance of CLSD, possibly by a dominant-negative disruptive effect on the SEC23 multimer.

The LOVD3 platform: efficient genome-wide sharing of genetic variants

Gene variant databases are the backbone of DNA-based diagnostics. These databases, also called Locus-Specific DataBases (LSDBs), store information on variants in the human genome and the observed phenotypic consequences. The largest collection of public databases uses the free, open-source LOVD software platform. To cope with the current demand for online databases, we have entirely redesigned the LOVD software. LOVD3 is genome-centered and can be used to store summary variant data, as well as full case-level data with information on individuals, phenotypes, screenings, and variants. While built on a standard core, the software is highly flexible and allows personalization to cope with the largely different demands of gene/disease database curators. LOVD3 follows current standards and includes tools to check variant descriptions, generate HTML files of reference sequences, predict the consequences of exon deletions/duplications on the reading frame, and link to genomic views in the different genomes browsers. It includes APIs to collect and submit data. The software is used by about 100 databases, of which 56 public LOVD instances are registered on our website and together contain 1,000,000,000 variant observations in 1,500,000 individuals. 42 LOVD instances share data with the federated LOVD data network containing 3,000,000 unique variants in 23,000 genes. This network can be queried directly, quickly identifying LOVD instances containing relevant information on a searched variant.

An Alu-mediated duplication in NMNAT1, involved in NAD biosynthesis, causes a novel syndrome, SHILCA, affecting multiple tissues and organs

We investigated the genetic origin of the phenotype displayed by three children from two unrelated Italian families, presenting with a previously unrecognized autosomal recessive disorder that included a severe form of spondylo-epiphyseal dysplasia, sensorineural hearing loss, intellectual disability and Leber congenital amaurosis (SHILCA), as well as some brain anomalies that were visible at the MRI. Autozygome-based analysis showed that these children shared a 4.76 Mb region of homozygosity on chromosome 1, with an identical haplotype. Nonetheless, whole-exome sequencing failed to identify any shared rare coding variants, in this region or elsewhere. We then determined the transcriptome of patients' fibroblasts by RNA sequencing, followed by additional whole-genome sequencing experiments. Gene expression analysis revealed a 4-fold downregulation of the gene NMNAT1, residing indeed in the shared autozygous interval. Short- and long-read whole-genome sequencing highlighted a duplication involving 2 out of the 5 exons of NMNAT1 main isoform (NM_022787.3), leading to the production of aberrant mRNAs. Pathogenic variants in NMNAT1 have been previously shown to cause non-syndromic Leber congenital amaurosis (LCA). However, no patient with null biallelic mutations has ever been described, and murine Nmnat1 knockouts show embryonic lethality, indicating that complete absence of NMNAT1 activity is probably not compatible with life. The rearrangement found in our cases, presumably causing a strong but not complete reduction of enzymatic activity, may therefore result in an intermediate syndromic phenotype with respect to LCA and lethality.