Variant filtering, digenic variants, and other challenges in clinical sequencing: a lesson from fibrillinopathies

XOR-Derived ROS in Tie2-Lineage Cells Including Endothelial Cells Promotes Aortic Aneurysm Progression in Marfan Syndrome

BACKGROUND

Marfan syndrome (MFS) is an inherited disorder caused by mutations in the FBN1 gene encoding fibrillin-1, a matrix component of extracellular microfibrils. The main cause of morbidity and mortality in MFS is thoracic aortic aneurysm and dissection, but the underlying mechanisms remain undetermined.

METHODS

To elucidate the role of endothelial XOR (xanthine oxidoreductase)-derived reactive oxygen species in aortic aneurysm progression, we inhibited in vivo function of XOR either by endothelial cell (EC)-specific disruption of the Xdh gene or by systemic administration of an XOR inhibitor febuxostat in MFS mice harboring the Fbn1 missense mutation p.(Cys1041Gly). We assessed the aberrant activation of mechanosensitive signaling in the ascending aorta of Fbn1C1041G/+ mice. Further analysis of human aortic ECs investigated the mechanisms by which mechanical stress upregulates XOR expression.

RESULTS

We found a significant increase in reactive oxygen species generation in the ascending aorta of patients with MFS and Fbn1C1041G/+ mice, which was associated with a significant increase in protein expression and enzymatic activity of XOR protein in aortic ECs. Genetic disruption of Xdh in ECs or treatment with febuxostat significantly suppressed aortic aneurysm progression and improved perivascular infiltration of macrophages. Mechanistically, mechanosensitive signaling involving FAK (focal adhesion kinase)-p38 MAPK (p38 mitogen-activated protein kinase) and Egr-1 (early growth response-1) was aberrantly activated in the ascending aorta of Fbn1C1041G/+ mice, and mechanical stress on human aortic ECs upregulated XOR expression through Egr-1 upregulation. Consistently, EC-specific knockout of XOR or systemic administration of febuxostat in Fbn1C1041G/+ mice suppressed reactive oxygen species generation, FAK-p38 MAPK activation, and Egr-1 upregulation.

CONCLUSIONS

Aberrant activation of mechanosensitive signaling in vascular ECs triggered endothelial XOR activation and reactive oxygen species generation, which contributes to the progression of aortic aneurysms in MFS. These findings highlight a drug repositioning approach using a uric acid-lowering drug febuxostat as a potential therapy for MFS.

Variants that get straight to your heart - Cardiogenetic secondary findings in exome sequencing

BACKGROUND

Exome sequencing has been established as a fundamental tool in genetic diagnostics. It may also provide information about variants in genes unrelated to the primary purpose, so-called secondary findings. Especially, diagnoses of unnoticed inborn cardiac diseases are of high clinical relevance due to therapeutic options in context of prevention of sudden cardiac death.

METHODS

Exome data of 9962 individuals was analysed for relevant cardiogenetic findings. Genes were selected according to ACMG recommendations for secondary findings (v.3.1). First, a filter for (likely) pathogenic variants, published in the ClinVar database, was used. Second, exome data was screened for loss of function (LoF) variants in genes in which LoF is a known disease pathomechanism. All variants were evaluated by geneticists regarding their pathogenicity.

RESULTS

Pathogenic or likely pathogenic variants were identified in 136 different individuals (136/9962, 1.4%), with the Low-Density Lipoprotein Receptor gene (LDLR, 24/136, 17.6%) and the Titin gene (TTN, 24/136, 17.6%), being the most frequently affected ones. 31.6% (43/136) of the identified variants had been reported beforehand, while 47.1% (64/136) had not been reported. The remaining cases (29/136, 21.3%) were part of research projects with no written reports. In 26.5% (36/136), the finding would have been missed, if only index patients and not their parents had been screened for secondary findings in case of trio ES.

CONCLUSION

As demonstrated in our study, at least one or two out of one hundred people are likely to carry a pathogenic cardiogenetic variant. Counselling geneticist and clinicians need to be aware of these findings in exome and genome sequencing. Informed consent of the patient regarding the report of secondary findings should absolutely be obtained beforehand.

Reviewing hereditary connective tissue disorders: Proposals of harmonic medicolegal assessments

Hereditary connective tissue disorders (HCTDs) are a heterogeneous group of inherited diseases. These disorders show genetic mutations with loss of function of primary components of connective tissue, such as collagen and elastic fibers. There are more than 200 conditions that involve hereditary connective tissue disorders, while the most known are Marfan syndrome, Osteogenesis Imperfecta, and Ehlers-Danlos syndromes. These disorders need continuous updates, multidisciplinary skills, and specific methodologic evaluations sharing many medicolegal issues. Marfan syndrome and Ehlers-Danlos syndromes show a high risk of early sudden death. As a consequence of this, postmortem genetic testing can identify novel genotype-phenotype correlations which help the clinicians to assess personalized cardiovascular screening programs among the ill subjects. Genetic testing is also essential to identify children suffering from Osteogenesis Imperfecta, especially when a physical abuse is clinically suspected. However, this is a well-known clinical problem even though there are still challenges to interpret genetic data and variants of unknown significance due to the current extensive use of new genetic/genomic techniques. Additionally, the more significant applications and complexities of genomic testing raise novel responsibilities on the clinicians, geneticists, and forensic practitioners as well, increasing potential liability and medical malpractice claims. This systematic review provides a detailed overview on how multidisciplinary skills belonging to clinicians, medicolegal consultants, radiologists, and geneticists can cooperate to manage HCTDs from autopsy or clinical findings to genetic testing. Thus, technical aspects need to be addressed to the medicolegal community since there is no consensus works or guidelines which specifically discuss these issues.

Smart variant filtering - A blueprint solution for massively parallel sequencing-based variant analysis

Massively parallel sequencing helps create new knowledge on genes, variants and their association with disease phenotype. This important technological advancement simultaneously makes clinical decision making, using genomic information for cancer patients, more complex. Currently, identifying actionable pathogenic variants with diagnostic, prognostic, or predictive impact requires substantial manual effort. Objective: The purpose is to design a solution for clinical diagnostics of lymphoma, specifically for systematic variant filtering and interpretation. Methods: A scoping review and demonstrations from specialists serve as a basis for a blueprint of a solution for massively parallel sequencing-based genetic diagnostics. Results: The solution uses machine learning methods to facilitate decision making in the diagnostic process. A validation round of interviews with specialists consolidated the blueprint and anchored it across all relevant expert disciplines. The scoping review identified four components of variant filtering solutions: algorithms and Artificial Intelligence (AI) applications, software, bioinformatics pipelines and variant filtering strategies. The blueprint describes the input, the AI model and the interface for dynamic browsing. Conclusion: An AI-augmented system is designed for predicting pathogenic variants. While such a system can be used to classify identified variants, diagnosticians should still evaluate the classification's accuracy, make corrections when necessary, and ultimately decide which variants are truly pathogenic.

Re-evaluation of a Fibrillin-1 Gene Variant of Uncertain Significance Using the ClinGen Guidelines

Background

Marfan syndrome (MFS) is caused by fibrillin-1 gene (FBN1) variants. Mutational hotspots and/or well-established critical functional domains of FBN1 include cysteine residues, calcium-binding consensus sequences, and amino acids related to interdomain packaging. Previous guidelines for variant interpretation do not reflect the features of genes or related diseases. Using the Clinical Genome Resource (ClinGen) FBN1 variant curation expert panel (VCEP), we re-evaluated FBN1 germline variants reported as variants of uncertain significance (VUSs).

Methods

We re-evaluated 26 VUSs in FBN1 reported in 161 patients with MFS. We checked the variants in the Human Genome Mutation Database, ClinVar, and VarSome databases and assessed their allele frequencies using the gnomAD database. Patients' clinical information was reviewed.

Results

Four missense variants affecting cysteines (c.460T>C, c.1006T>C, c.5330G>C, and c.8020T>C) were reclassified as likely pathogenic and were assigned PM1_strong or PM1. Two intronic variants were reclassified as benign by granting BA1 (stand-alone). Four missense variants were reclassified as likely benign. BP5 criteria were applied in cases with an alternate molecular basis for disease, one of which (c.7231G>A) was discovered alongside a pathogenic de novo COL3A1 variant (c.1988G>T, p.Gly633Val).

Conclusions

Considering the high penetrance of FBN1 variants and clinical variability of MFS, the detection of pathogenic variants is important. The ClinGen FBN1 VCEP encompasses mutational hotspots and/or well-established critical functional domains and adjusts the criteria specifically for MFS; therefore, it is beneficial not only for identifying pathogenic FBN1 variants but also for distinguishing these variants from those that cause other connective tissue disorders with overlapping clinical features.

A population-based survey of FBN1 variants in Iceland reveals underdiagnosis of Marfan syndrome

Marfan syndrome (MFS) is an autosomal dominant condition characterized by aortic aneurysm, skeletal abnormalities, and lens dislocation, and is caused by variants in the FBN1 gene. To explore causes of MFS and the prevalence of the disease in Iceland we collected information from all living individuals with a clinical diagnosis of MFS in Iceland (n = 32) and performed whole-genome sequencing of those who did not have a confirmed genetic diagnosis (27/32). Moreover, to assess a potential underdiagnosis of MFS in Iceland we attempted a genotype-based approach to identify individuals with MFS. We interrogated deCODE genetics' database of 35,712 whole-genome sequenced individuals to search for rare sequence variants in FBN1. Overall, we identified 15 pathogenic or likely pathogenic variants in FBN1 in 44 individuals, only 22 of whom were previously diagnosed with MFS. The most common of these variants, NM_000138.4:c.8038 C > T p.(Arg2680Cys), is present in a multi-generational pedigree, and was found to stem from a single forefather born around 1840. The p.(Arg2680Cys) variant associates with a form of MFS that seems to have an enrichment of abdominal aortic aneurysm, suggesting that this may be a particularly common feature of p.(Arg2680Cys)-associated MFS. Based on these combined genetic and clinical data, we show that MFS prevalence in Iceland could be as high as 1/6,600 in Iceland, compared to 1/10,000 based on clinical diagnosis alone, which indicates underdiagnosis of this actionable genetic disorder.

Case report: Identification of novel fibrillin-2 variants impacting disulfide bond and causing congenital contractural arachnodactyly

Background: Congenital contractural arachnodactyly (CCA) is an autosomal dominant connective tissue disorder with clinical features of arthrogryposis, arachnodactyly, crumpled ears, scoliosis, and muscular hypoplasia. The heterozygous pathogenic variants in FBN2 have been shown to cause CCA. Fibrillin-2 is related to the elasticity of the tissue and has been demonstrated to play an important role in the constitution of extracellular microfibrils in elastic fibers, providing strength and flexibility to the connective tissue that sustains the body's joints and organs. Methods: We recruited two Chinese families with arachnodactyly and bilateral arthrogryposis of the fingers. Whole-exome sequencing (WES) and co-segregation analysis were employed to identify their genetic etiologies. Three-dimensional protein models were used to analyze the pathogenic mechanism of the identified variants. Results: We have reported two CCA families and identified two novel missense variants in FBN2 (NM_001999.3: c.4093T>C, p.C1365R and c.2384G>T, p.C795F). The structural models of the mutant FBN2 protein in rats exhibited that both the variants could break disulfide bonds. Conclusion: We detected two FBN2 variants in two families with CCA. Our description expands the genetic profile of CCA and emphasizes the pathogenicity of disulfide bond disruption in FBN2.

Need for revision of the ACMG/AMP guidelines for interpretation of X-linked variants

The guidelines provided by American College of Medical Genetics and Genomics (ACMG) and the Association of Molecular Pathology (AMP) (ACMG/AMP guidelines) suggest a framework for the classification of clinical variants. However, the interpretations can be inconsistent, with each definition sometimes proving to be ambiguous. In particular, there can be difficulty with interpretation of variants related to the X-linked recessive trait. To confirm whether there are biases in the interpretation of inherited traits, we reanalyzed variants reported prior to the release of the ACMG/AMP guidelines. As expected, the interpretation ratio as pathogenic or likely pathogenic was significantly lower for variants related to the X-linked recessive trait. Evaluation of variants related to the X-linked recessive trait, hence, need to consider whether the variant is identified only in males in accordance with the X-linked recessive trait. The ACMG/AMP guidelines should be revised to eliminate the bias revealed in this study.

Surveillance and monitoring in vascular Ehlers-Danlos syndrome in European Reference Network For Rare Vascular Diseases (VASCERN)

Vascular Ehlers-Danlos syndrome (vEDS) is a rare genetic disorder clinically characterized by vascular, intestinal and uterine fragility and caused by heterozygous pathogenic variants in the COL3A1 gene. Management of patients with vEDS is difficult due to the unpredictability of the events and clear recommendations on the care of adults and children with vEDS are lacking. Therefore, we aimed to collect data on the current strategy of surveillance and monitoring of vEDS patients by expert centers in continental Europe and Great Britain, as a first step towards a consensus statement. A survey on the clinical management of vEDS was sent to all members of the Medium Sized Artery (MSA) Working Group of the European Reference Network for Rare Vascular Diseases (VASCERN) and other expert centers. All experts endorse the importance of monitoring patients with vEDS. Despite the absence of evidence based guidelines monitoring is considered in almost all countries, but screening intervals and modalities used for monitoring may differ among centers. There is a need for more prospective multicenter studies to define proper guidelines.

Potential predictors of severe cardiovascular involvement in Marfan syndrome: the emphasized role of genotype-phenotype correlations in improving risk stratification-a literature review

BACKGROUND

Marfan syndrome (MFS) is a genetically determined systemic connective tissue disorder, caused by a mutation in the FBN1 gene. In MFS mainly the cardiovascular, musculoskeletal and ocular systems are affected. The most dangerous manifestation of MFS is aortic dissection, which needs to be prevented by a prophylactic aortic root replacement.

MAIN BODY

The indication criteria for the prophylactic procedure is currently based on aortic diameter, however aortic dissections below the threshold defined in the guidelines have been reported, highlighting the need for a more accurate risk stratification system to predict the occurrence of aortic complications. The aim of this review is to present the current knowledge on the possible predictors of severe cardiovascular manifestations in MFS patients, demonstrating the wide range of molecular and radiological differences between people with MFS and healthy individuals, and more importantly between MFS patients with and without advanced aortic manifestations. These differences originating from the underlying common molecular pathological processes can be assessed by laboratory (e.g. genetic testing) and imaging techniques to serve as biomarkers of severe aortic involvement. In this review we paid special attention to the rapidly expanding field of genotype-phenotype correlations for aortic features as by collecting and presenting the ever growing number of correlations, future perspectives for risk stratification can be outlined.

CONCLUSIONS

Data on promising biomarkers of severe aortic complications of MFS have been accumulating steadily. However, more unifying studies are required to further evaluate the applicability of the discussed predictors with the aim of improving the risk stratification and therefore the life expectancy and quality of life of MFS patients.

Stepwise ABC system for classification of any type of genetic variant

The American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG-AMP) system for variant classification is score based with five classes: benign, likely benign, variant of unknown significance (VUS), likely pathogenic, and pathogenic. Here, we present a variant classification model that can be an add-on or alternative to ACMG classification: A stepwise system that can classify any type of genetic variant (e.g., hypomorphic alleles, imprinted alleles, copy number variants, runs of homozygosity, enhancer variants, and variants related to traits). We call it the ABC system because classification is first functional (A), then clinical (B), and optionally a standard comment that fits the clinical question is selected (C). Both steps A and B have 1–5 grading when knowledge is sufficient, if not, class “zero” is assigned. Functional grading (A) only concerns biological consequences with the stages normal function (1), likely normal function (2), hypothetical functional effect (3), likely functional effect (4), and proven functional effect (5). Clinical grading (B) is genotype–phenotype focused with the stages “right type of gene” (1), risk factor (2), and pathogenic (3–5, depending on penetrance). Both grades are listed for each variant and combined to generate a joint class ranging from A to F. Importantly, the A–F classes are linked to standard comments, reflecting laboratory or national policy. In step A, the VUS class is split into class 0 (true unknown) and class 3 (hypothetical functional effect based on molecular predictions or de novo occurrence), providing a rationale for variant-of-interest reporting when the clinical picture could fit the finding. The system gives clinicians a better guide to variant significance.

Low penetrance COL5A1 variants in a young patient with intracranial aneurysm and very mild signs of Ehlers-Danlos syndrome

Spontaneous cervical artery dissection (CeAD) is a major cause of ischemic stroke in young adults, whose genetic susceptibility factors are still largely unknown. Nevertheless, subtle ultrastructural connective tissue alterations (especially in the collagen fibril morphology) are recognized in a large proportion of CeAD patients, in which recent genetic investigations reported an enrichment of variants in genes associated with known connective tissue disorders. In this regard, COL5A1 variants have been reported in a small subset of CeAD patients, with or without classical Ehlers-Danlos syndrome (cEDS) features. We investigated a 22-year-old patient with intracranial aneurysm and mild connective tissue manifestations reminiscent of EDS. Whole-exome sequencing identified two COL5A1 missense variants in trans configuration: NM_000093.5:c.[1588G>A];[4135C>T], NP_000084.3:p.[(Gly530Ser)];[(Pro1379Ser)]. Functional assays demonstrated a significant decrease of collagen α1(V) chain expression in both heterozygous parents compared to control cells, and an additive effect of these two variants in the proband. Interestingly, both parents manifested very subtle EDS signs, such as atrophic scars, recurrent bone fractures, colonic diverticulosis, varicose veins, and osteoarthritis. Our findings emphasize the involvement of COL5A1 in the predisposition to vascular phenotypes and provide novel insights on the c.1588G>A variant, whose functional significance has not been definitely established. In fact, it was previously reported as both "disease modifying", and as a biallelic causative mutation (with heterozygous individuals showing subtle clinical signs of cEDS). We speculated that the c.1588G>A variant might lead to overt phenotype in combination with additional genetic "hits" lowering the collagen α1(V) chain expression below a hypothetical disease threshold.

Double heterozygous variants in FBN1 and FBN2 in a Thai woman with Marfan and Beals syndromes

A phenotype of an individual is resulted from an interaction among variants in several genes. Advanced molecular technologies allow us to identify more patients with mutations in more than one genes. Here, we studied a Thai woman with combined clinical features of Marfan (MFS) and Beals (BS) syndromes including frontal bossing, enophthalmos, myopia, the crumpled appearance to the top of the pinnae, midface hypoplasia, high arched palate, dermal stretch marks, aortic enlargement, mitral valve prolapse and regurgitation, aortic root dilatation, and progressive scoliosis. The aortic root enlargement was progressive to a diameter of 7.2 cm requiring an aortic root replacement at the age of 8 years. At her last visit when she was 19 years old, she had moderate aortic regurgitation. Exome sequencing revealed that she carried the c.3159C > G (p.Cys1053Trp) in exon 26 of FBN1 and c.2638G > A (p. Gly880Ser) in exon 20 of FBN2. The variant in FBN1 was de novo, while that in FBN2 was inherited from her unaffected mother. Both genes encode for fibrillins, which are essential for elastic fibers and can form the heterotypic microfibrils. Two defective fibrillins may synergistically worsen cardiovascular manifestations seen in our patient. In this study, we identified the fourth patient with both MFS and BS, carrying mutations in both FBN1 and FBN2.

An Improved Phenotype-Driven Tool for Rare Mendelian Variant Prioritization: Benchmarking Exomiser on Real Patient Whole-Exome Data

Next-generation sequencing has revolutionized rare disease diagnostics, but many patients remain without a molecular diagnosis, particularly because many candidate variants usually survive despite strict filtering. Exomiser was launched in 2014 as a Java tool that performs an integrative analysis of patients' sequencing data and their phenotypes encoded with Human Phenotype Ontology (HPO) terms. It prioritizes variants by leveraging information on variant frequency, predicted pathogenicity, and gene-phenotype associations derived from human diseases, model organisms, and protein-protein interactions. Early published releases of Exomiser were able to prioritize disease-causative variants as top candidates in up to 97% of simulated whole-exomes. The size of the tested real patient datasets published so far are very limited. Here, we present the latest Exomiser version 12.0.1 with many new features. We assessed the performance using a set of 134 whole-exomes from patients with a range of rare retinal diseases and known molecular diagnosis. Using default settings, Exomiser ranked the correct diagnosed variants as the top candidate in 74% of the dataset and top 5 in 94%; not using the patients' HPO profiles (i.e., variant-only analysis) decreased the performance to 3% and 27%, respectively. In conclusion, Exomiser is an effective support tool for rare Mendelian phenotype-driven variant prioritization.

Added Value of Clinical Sequencing: WGS-Based Profiling of Pharmacogenes

Although several pharmacogenetic (PGx) predispositions affecting drug efficacy and safety are well established, drug selection and dosing as well as clinical trials are often performed in a non-pharmacogenetically-stratified manner, ultimately burdening healthcare systems. Pre-emptive PGx testing offers a solution which is often performed using microarrays or targeted gene panels, testing for common/known PGx variants. However, as an added value, whole-genome sequencing (WGS) could detect not only disease-causing but also pharmacogenetically-relevant variants in a single assay. Here, we present our WGS-based pipeline that extends the genetic testing of Mendelian diseases with PGx profiling, enabling the detection of rare/novel PGx variants as well. From our in-house WGS (PCR-free 60× PE150) data of 547 individuals we extracted PGx variants with drug-dosing recommendations of the Dutch Pharmacogenetics Working Group (DPWG). Furthermore, we explored the landscape of DPWG pharmacogenes in gnomAD and our in-house cohort as well as compared bioinformatic tools for WGS-based structural variant detection in CYP2D6. We show that although common/known PGx variants comprise the vast majority of detected DPWG pharmacogene alleles, for better precision medicine, PGx testing should move towards WGS-based approaches. Indeed, WGS-based PGx profiling is not only feasible and future-oriented but also the most comprehensive all-in-one approach without generating significant additional costs.

Systematic dissection of biases in whole-exome and whole-genome sequencing reveals major determinants of coding sequence coverage

Advantages and diagnostic effectiveness of the two most widely used resequencing approaches, whole exome (WES) and whole genome (WGS) sequencing, are often debated. WES dominated large-scale resequencing projects because of lower cost and easier data storage and processing. Rapid development of 3^rd generation sequencing methods and novel exome sequencing kits predicate the need for a robust statistical framework allowing informative and easy performance comparison of the emerging methods. In our study we developed a set of statistical tools to systematically assess coverage of coding regions provided by several modern WES platforms, as well as PCR-free WGS. We identified a substantial problem in most previously published comparisons which did not account for mappability limitations of short reads. Using regression analysis and simple machine learning, as well as several novel metrics of coverage evenness, we analyzed the contribution from the major determinants of CDS coverage. Contrary to a common view, most of the observed bias in modern WES stems from mappability limitations of short reads and exome probe design rather than sequence composition. We also identified the ~ 500 kb region of human exome that could not be effectively characterized using short read technology and should receive special attention during variant analysis. Using our novel metrics of sequencing coverage, we identified main determinants of WES and WGS performance. Overall, our study points out avenues for improvement of enrichment-based methods and development of novel approaches that would maximize variant discovery at optimal cost.

Variant filtering, digenic variants, and other challenges in clinical sequencing: a lesson from fibrillinopathies

Genome-scale high-throughput sequencing enables the detection of unprecedented numbers of sequence variants. Variant filtering and interpretation are facilitated by mutation databases, in silico tools, and population-based reference datasets such as ExAC/gnomAD, while variants are classified using the ACMG/AMP guidelines. These methods, however, pose clinically relevant challenges. We queried the gnomAD dataset for (likely) pathogenic variants in genes causing autosomal-dominant disorders. Furthermore, focusing on the fibrillinopathies Marfan syndrome (MFS) and congenital contractural arachnodactyly (CCA), we screened 500 genomes of our patients for co-occurring variants in FBN1 and FBN2. In gnomAD, we detected 2653 (likely) pathogenic variants in 253 genes associated with autosomal-dominant disorders, enabling the estimation of variant-filtering thresholds and disease predisposition/prevalence rates. In our database, we discovered two families with hitherto unreported co-occurrence of FBN1/FBN2 variants causing phenotypes with mixed or modified MFS/CCA clinical features. We show that (likely) pathogenic gnomAD variants may be more frequent than expected and are challenging to classify according to the ACMG/AMP guidelines as well as that fibrillinopathies are likely underdiagnosed and may co-occur. Consequently, selection of appropriate frequency cutoffs, recognition of digenic variants, and variant classification represent considerable challenges in variant interpretation. Neglecting these challenges may lead to incomplete or missed diagnoses.