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.

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.

A next-generation-sequencing panel for mutational analysis of dominant acute hepatic porphyrias.

Molecular diagnosis of autosomal dominant acute hepatic porphyrias (AHPs) plays an important role in the management of these disorders. To introduce next generation sequencing (NGS) to the porphyria diagnosis, we designed a panel that contained four genes, ALAS1, HMBS, CPOX and PPOX for mutational analysis of acute intermittent porphyria (AIP), hereditary coproporphyria (HCP) and variegate porphyria (VP). To validate the AHP panel, 30 samples with known pathogenic variants as determined by Sanger sequencing, were analyzed using the Ion PGM™. Among them, nine have so far not been reported. The pathogenic variants were identified and annotated manually in IGV by three individuals who were blinded to the Sanger results. The AHP panel consists of 95 amplicons that covers 92% of the coding region of the four genes. Of the 95 amplicons, 93 had an average read-depth of >500 reads. In 29 of the 30 tested samples, pathogenic variants were correctly identified and annotated. The number of reads from the mutated alleles were approximately 50% of the total. The annotation of a 22-bp duplication with NGS differed from that of Sanger by one nucleotide. NGS showed an advantage in allelic discrimination over Sanger sequencing and was also able to detect a known somatic variant in the HMBS gene. The AHP panel will be applied in the initial diagnosis of new patients. Any sequence variations with a frequency of ≥10% will be confirmed by Sanger sequencing. The cost-effectiveness of a NGS approach for AHP in a diagnostic laboratory needs to be further assessed.

New insights into the performance of human whole-exome capture platforms.

Whole exome sequencing (WES) is increasingly used in research and diagnostics. WES users expect coverage of the entire coding region of known genes as well as sufficient read depth for the covered regions. It is, however, unknown which recent WES platform is most suitable to meet these expectations. We present insights into the performance of the most recent standard exome enrichment platforms from Agilent, NimbleGen and Illumina applied to six different DNA samples by two sequencing vendors per platform. Our results suggest that both Agilent and NimbleGen overall perform better than Illumina and that the high enrichment performance of Agilent is stable among samples and between vendors, whereas NimbleGen is only able to achieve vendor- and sample-specific best exome coverage. Moreover, the recent Agilent platform overall captures more coding exons with sufficient read depth than NimbleGen and Illumina. Due to considerable gaps in effective exome coverage, however, the three platforms cannot capture all known coding exons alone or in combination, requiring improvement. Our data emphasize the importance of evaluation of updated platform versions and suggest that enrichment-free whole genome sequencing can overcome the limitations of WES in sufficiently covering coding exons, especially GC-rich regions, and in characterizing structural variants.

Clinical sequencing: is WGS the better WES?

Current clinical next-generation sequencing is done by using gene panels and exome analysis, both of which involve selective capturing of target regions. However, capturing has limitations in sufficiently covering coding exons, especially GC-rich regions. We compared whole exome sequencing (WES) with the most recent PCR-free whole genome sequencing (WGS), showing that only the latter is able to provide hitherto unprecedented complete coverage of the coding region of the genome. Thus, from a clinical/technical point of view, WGS is the better WES so that capturing is no longer necessary for the most comprehensive genomic testing of Mendelian disorders.

Novel Insights into the Aortic Mechanical Properties of Mice Modeling Hereditary Aortic Diseases.

OBJECTIVE: Hereditary aortic diseases (hADs) increase the risk of aortic dissections and ruptures. Recently, we have established an objective approach to measure the rupture force of the murine aorta, thereby explaining the outcomes of clinical studies and assessing the added value of approved drugs in vascular Ehlers-Danlos syndrome (vEDS). Here, we applied our approach to six additional mouse hAD models. MATERIAL AND METHODS: We used two mouse models (Fbn1C1041G and Fbn1mgR ) of Marfan syndrome (MFS) as well as one smooth-muscle-cell-specific knockout (SMKO) of Efemp2 and three CRISPR/Cas9-engineered knock-in models (Ltbp1, Mfap4, and Timp1). One of the two MFS models was subjected to 4-week-long losartan treatment. Per mouse, three rings of the thoracic aorta were prepared, mounted on a tissue puller, and uniaxially stretched until rupture. RESULTS: The aortic rupture force of the SMKO and both MFS models was significantly lower compared with wild-type mice but in both MFS models higher than in mice modeling vEDS. In contrast, the Ltbp1, Mfap4, and Timp1 knock-in models presented no impaired aortic integrity. As expected, losartan treatment reduced aneurysm formation but surprisingly had no impact on the aortic rupture force of our MFS mice. CONCLUSION: Our read-out system can characterize the aortic biomechanical integrity of mice modeling not only vEDS but also related hADs, allowing the aortic-rupture-force-focused comparison of mouse models. Furthermore, aneurysm progression alone may not be a sufficient read-out for aortic rupture, as antihypertensive drugs reducing aortic dilatation might not strengthen the weakened aortic wall. Our results may enable identification of improved medical therapies of hADs.

Potential of whole-genome sequencing-based pharmacogenetic profiling.

Pharmacogenetics represents a major driver of precision medicine, promising individualized drug selection and dosing. Traditionally, pharmacogenetic profiling has been performed using targeted genotyping that focuses on common/known variants. Recently, whole-genome sequencing (WGS) is emerging as a more comprehensive short-read next-generation sequencing approach, enabling both gene diagnostics and pharmacogenetic profiling, including rare/novel variants, in a single assay. Using the example of the pharmacogene CYP2D6, we demonstrate the potential of WGS-based pharmacogenetic profiling as well as emphasize the limitations of short-read next-generation sequencing. In the near future, we envision a shift toward long-read sequencing as the predominant method for gene diagnostics and pharmacogenetic profiling, providing unprecedented data quality and improving patient care.

Celiprolol but not losartan improves the biomechanical integrity of the aorta in a mouse model of vascular Ehlers-Danlos syndrome.

AIMS: Antihypertensive drugs are included in the medical therapy of vascular Ehlers-Danlos syndrome (vEDS). The ß-blocker celiprolol has been suggested to prevent arterial damage in vEDS, but the underlying mechanism remains unclear. It is also unknown whether the widely used angiotensin II receptor type 1 antagonist losartan has a therapeutic effect in vEDS. Here, we evaluated the impact of celiprolol and losartan on the biomechanical integrity of the vEDS thoracic aorta. METHODS AND RESULTS: We established a new approach to measure the maximum tensile force at rupture of uniaxially stretched murine thoracic aortic rings. In a vEDS model, which we (re-)characterized here at molecular level, heterozygous mice showed a significant reduction in the rupture force compared to wild-type mice, reflecting the increased mortality due to aortic rupture. For the assessment of treatment effects, heterozygous mice at 4 weeks of age underwent a 4-week treatment with celiprolol, losartan, and, as a proof-of-concept drug, the matrix metalloproteinase inhibitor doxycycline. Compared to age- and sex-matched untreated heterozygous mice, treatment with doxycycline or celiprolol resulted in a significant increase of rupture force, whereas no significant change was detected upon losartan treatment. CONCLUSIONS: In a vEDS model, celiprolol or doxycycline, but not losartan, can improve the biomechanical integrity of the aortic wall, thereby potentially reducing the risk of dissection and rupture. As doxycycline is a broad-spectrum antibiotic with considerable side effects, celiprolol may be more suitable for a long-term therapy and thus rather indicated for the medication of patients with vEDS.

Hungarian Marfan family with large FBN1 deletion calls attention to copy number variation detection in the current NGS era.

Copy number variations (CNVs) comprise about 10% of reported disease-causing mutations in Mendelian disorders. Nevertheless, pathogenic CNVs may have been under-detected due to the lack or insufficient use of appropriate detection methods. In this report, on the example of the diagnostic odyssey of a patient with Marfan syndrome (MFS) harboring a hitherto unreported 32-kb FBN1 deletion, we highlight the need for and the feasibility of testing for CNVs (>1 kb) in Mendelian disorders in the current next-generation sequencing (NGS) era.

Precise breakpoint localization of large genomic deletions using PacBio and Illumina next-generation sequencers.

Herein we present the applicability of single-molecule (PacBio RS) and second-generation sequencing technology (Illumina) to the characterization of large genomic deletions. By testing samples previously characterized using a Sanger approach, our methods determined that both next-generation sequencing platforms were able to identify the position of deletion breakpoints. Our results point out various advantages of next-generation sequencing platforms when characterizing genomic deletions; however, special attention must be dedicated to identical sequences flanking the breakpoints, such as poly(N) motifs.

Hemizygous deletion of COL3A1, COL5A2, and MSTN causes a complex phenotype with aortic dissection: a lesson for and from true haploinsufficiency.

Aortic dilatation/dissection (AD) can occur spontaneously or in association with genetic syndromes, such as Marfan syndrome (MFS; caused by FBN1 mutations), MFS type 2 and Loeys-Dietz syndrome (associated with TGFBR1/TGFBR2 mutations), and Ehlers-Danlos syndrome (EDS) vascular type (caused by COL3A1 mutations). Although mutations in FBN1 and TGFBR1/TGFBR2 account for the majority of AD cases referred to us for molecular genetic testing, we have obtained negative results for these genes in a large cohort of AD patients, suggesting the involvement of additional genes or acquired factors. In this study we assessed the effect of COL3A1 deletions/duplications in this cohort. Multiplex ligation-dependent probe amplification (MLPA) analysis of 100 unrelated patients identified one hemizygous deletion of the entire COL3A1 gene. Subsequent microarray analyses and sequencing of breakpoints revealed the deletion size of 3,408,306 bp at 2q32.1q32.3. This deletion affects not only COL3A1 but also 21 other known genes (GULP1, DIRC1, COL5A2, WDR75, SLC40A1, ASNSD1, ANKAR, OSGEPL1, ORMDL1, LOC100129592, PMS1, MSTN, C2orf88, HIBCH, INPP1, MFSD6, TMEM194B, NAB1, GLS, STAT1, and STAT4), mutations in three of which (COL5A2, SLC40A1, and MSTN) have also been associated with an autosomal dominant disorder (EDS classical type, hemochromatosis type 4, and muscle hypertrophy). Physical and laboratory examinations revealed that true haploinsufficiency of COL3A1, COL5A2, and MSTN, but not that of SLC40A1, leads to a clinical phenotype. Our data not only emphasize the impact/role of COL3A1 in AD patients but also extend the molecular etiology of several disorders by providing hitherto unreported evidence for true haploinsufficiency of the underlying gene.