Functional filter for whole-genome sequencing data identifies HHT and stress-associated non-coding SMAD4 polyadenylation site variants >5 kb from coding DNA

Rare disease gene association discovery in the 100,000 Genomes Project

Up to 80% of rare disease patients remain undiagnosed after genomic sequencing1, with many probably involving pathogenic variants in yet to be discovered disease-gene associations. To search for such associations, we developed a rare variant gene burden analytical framework for Mendelian diseases, and applied it to protein-coding variants from whole-genome sequencing of 34,851 cases and their family members recruited to the 100,000 Genomes Project2. A total of 141 new associations were identified, including five for which independent disease-gene evidence was recently published. Following in silico triaging and clinical expert review, 69 associations were prioritized, of which 30 could be linked to existing experimental evidence. The five associations with strongest overall genetic and experimental evidence were monogenic diabetes with the known β cell regulator3,4 UNC13A, schizophrenia with GPR17, epilepsy with RBFOX3, Charcot-Marie-Tooth disease with ARPC3 and anterior segment ocular abnormalities with POMK. Further confirmation of these and other associations could lead to numerous diagnoses, highlighting the clinical impact of large-scale statistical approaches to rare disease-gene association discovery.

Hereditary haemorrhagic telangiectasia

Hereditary haemorrhagic telangiectasia (HHT) is a vascular dysplasia inherited as an autosomal dominant trait and caused by loss-of-function pathogenic variants in genes encoding proteins of the BMP signalling pathway. Up to 90% of disease-causal variants are observed in ENG and ACVRL1, with SMAD4 and GDF2 less frequently responsible for HHT. In adults, the most frequent HHT manifestations relate to iron deficiency and anaemia owing to recurrent epistaxis (nosebleeds) or bleeding from gastrointestinal telangiectases. Arteriovenous malformations (AVMs) in the lungs, liver and the central nervous system cause additional major complications and often complex symptoms, primarily due to vascular shunting, which is right-to-left through pulmonary AVMs (causing ischaemic stroke or cerebral abscess) and left-to-right through systemic AVMs (causing high cardiac output). Children usually experience isolated epistaxis; in rare cases, childhood complications occur from large AVMs in the lungs or central nervous system. Management goals encompass control of epistaxis and intestinal bleeding from telangiectases, screening for and treatment of iron deficiency (with or without anaemia) and AVMs, genetic counselling and evaluation of at-risk family members. Novel therapeutics, such as systemic antiangiogenic therapies, are actively being investigated. Although HHT is associated with increased morbidity, the appropriate screening and treatment of visceral AVMs, and the effective management of bleeding and anaemia, improves quality of life and overall survival.

Endogenous plasma resuspension of peripheral blood mononuclear cells prevents preparative-associated stress that modifies polyA-enriched RNA responses to subsequent acute stressors

Human peripheral blood mononuclear cells (PBMCs) are used to examine biological processes and disease, when basal variability in cellular activation and splicing is described and unexplained. Using isolation systems that maintained buffy coat cells (PBMCs, platelets) in their own plasma, poly-A enriched RNA-sequencing (RNASeq) detected 42,720 Ensembl gene IDs, including >95% of the top 100 Genotype Tissue Expression Project (GTEx)-expressed genes in lung, colon, heart, skeletal muscle and liver, and 10/17 clinically-actionable genes listed by the Pharmacogenomics Knowledgebase. Transcriptome changes were defined after 1h treatment with 32°C hypothermia (hsp70 family member change), 10 μmol/L ferric citrate that had no discernible effect, and 100 μg/mL cycloheximide leading to induction of primary response (immediate early) genes including IL1B and TNF. Same-donor PBMCs prepared conventionally using washes then resuspension in serum-supplemented media demonstrated basal upregulation of stress signalling pathway genes that masked and overlapped differential gene expression profiles after 100 µg/L cycloheximide. Plasma-resuspended PBMCs demonstrated minor transcriptome changes after 40 μmol/L ferric citrate, whereas consistent and greater magnitude changes were observed for washed/media-resuspended PBMCs. We conclude that endogenous plasma-maintained PBMCs provide a more robust platform to interrogate acute cellular perturbations triggering innate immunity, and that varying susceptibility of PBMCs to preparative stresses is an important component of experimental variability.

Mutations causing premature termination codons discriminate and generate cellular and clinical variability in HHT

ABSTRACT

For monogenic diseases caused by pathogenic loss-of-function DNA variants, attention focuses on dysregulated gene-specific pathways, usually considering molecular subtypes together within causal genes. To better understand phenotypic variability in hereditary hemorrhagic telangiectasia (HHT), we subcategorized pathogenic DNA variants in ENG/endoglin, ACVRL1/ALK1, and SMAD4 if they generated premature termination codons (PTCs) subject to nonsense-mediated decay. In 3 patient cohorts, a PTC-based classification system explained some previously puzzling hemorrhage variability. In blood outgrowth endothelial cells (BOECs) derived from patients with ACVRL1+/PTC, ENG+/PTC, and SMAD4+/PTC genotypes, PTC-containing RNA transcripts persisted at low levels (8%-23% expected, varying between replicate cultures); genes differentially expressed to Bonferroni P < .05 in HHT+/PTC BOECs clustered significantly only to generic protein terms (isopeptide-bond/ubiquitin-like conjugation) and pulse-chase experiments detected subtle protein maturation differences but no evidence for PTC-truncated protein. BOECs displaying highest PTC persistence were discriminated in unsupervised hierarchical clustering of near-invariant housekeeper genes, with patterns compatible with higher cellular stress in BOECs with >11% PTC persistence. To test directionality, we used a HeLa reporter system to detect induction of activating transcription factor 4 (ATF4), which controls expression of stress-adaptive genes, and showed that ENG Q436X but not ENG R93X directly induced ATF4. AlphaFold accurately modeled relevant ENG domains, with AlphaMissense suggesting that readthrough substitutions would be benign for ENG R93X and other less rare ENG nonsense variants but more damaging for Q436X. We conclude that PTCs should be distinguished from other loss-of-function variants, PTC transcript levels increase in stressed cells, and readthrough proteins and mechanisms provide promising research avenues.

Pathogenic Variant Frequencies in Hereditary Haemorrhagic Telangiectasia Support Clinical Evidence of Protection from Myocardial Infarction

Hereditary haemorrhagic telangiectasia (HHT) is a vascular dysplasia inherited as an autosomal dominant trait, due to a single heterozygous loss-of-function variant, usually in ACVRL1 (encoding activin receptor-like kinase 1 [ALK1]), ENG (encoding endoglin [CD105]), or SMAD4. In a consecutive single-centre series of 37 positive clinical genetic tests performed in 2021-2023, a skewed distribution pattern was noted, with 30 of 32 variants reported only once, but ACVRL1 c.1231C>T (p.Arg411Trp) identified as the disease-causal gene in five different HHT families. In the same centre's non-overlapping 1992-2020 series where 110/134 (82.1%) HHT-causal variants were reported only once, ACVRL1 c.1231C>T (p.Arg411Trp) was identified in nine further families. In a 14-country, four-continent HHT Mutation Database where 181/250 (72.4%) HHT-causal variants were reported only once, ACVRL1 c.1231C>T (p.Arg411Trp) was reported by 12 different laboratories, the adjacent ACVRL1 c.1232G>A (p.Arg411Gln) by 14, and ACVRL1 c.1120C>T (p.Arg374Trp) by 18. Unlike the majority of HHT-causal ACVRL1 variants, these encode ALK1 protein that reaches the endothelial cell surface but fails to signal. Six variants of this type were present in the three series and were reported 6.8-25.5 (mean 8.9) times more frequently than the other ACVRL1 missense variants (all p-values < 0.0039). Noting lower rates of myocardial infarction reported in HHT, we explore potential mechanisms, including a selective paradigm relevant to ALK1's role in the initiating event of atherosclerosis, where a plausible dominant negative effect of these specific variants can be proposed. In conclusion, there is an ~9-fold excess of kinase-inactive, cell surface-expressed ACVRL1/ALK1 pathogenic missense variants in HHT. The findings support further examination of differential clinical and cellular phenotypes by HHT causal gene molecular subtypes.

Rare disease gene association discovery from burden analysis of the 100,000 Genomes Project data

To discover rare disease-gene associations, we developed a gene burden analytical framework and applied it to rare, protein-coding variants from whole genome sequencing of 35,008 cases with rare diseases and their family members recruited to the 100,000 Genomes Project (100KGP). Following in silico triaging of the results, 88 novel associations were identified including 38 with existing experimental evidence. We have published the confirmation of one of these associations, hereditary ataxia with UCHL1 , and independent confirmatory evidence has recently been published for four more. We highlight a further seven compelling associations: hypertrophic cardiomyopathy with DYSF and SLC4A3 where both genes show high/specific heart expression and existing associations to skeletal dystrophies or short QT syndrome respectively; monogenic diabetes with UNC13A with a known role in the regulation of β cells and a mouse model with impaired glucose tolerance; epilepsy with KCNQ1 where a mouse model shows seizures and the existing long QT syndrome association may be linked; early onset Parkinson's disease with RYR1 with existing links to tremor pathophysiology and a mouse model with neurological phenotypes; anterior segment ocular abnormalities associated with POMK showing expression in corneal cells and with a zebrafish model with developmental ocular abnormalities; and cystic kidney disease with COL4A3 showing high renal expression and prior evidence for a digenic or modifying role in renal disease. Confirmation of all 88 associations would lead to potential diagnoses in 456 molecularly undiagnosed cases within the 100KGP, as well as other rare disease patients worldwide, highlighting the clinical impact of a large-scale statistical approach to rare disease gene discovery.

Pharmacogenomic Considerations for Anticoagulant Prescription in Patients with Hereditary Haemorrhagic Telangiectasia

Hereditary haemorrhagic telangiectasia (HHT) is a vascular dysplasia that commonly results in bleeding but with frequent indications for therapeutic anticoagulation. Our aims were to advance the understanding of drug-specific intolerance and evaluate if there was an indication for pharmacogenomic testing. Genes encoding proteins involved in the absorption, distribution, metabolism, and excretion of warfarin, heparin, and direct oral anticoagulants (DOACs) apixaban, rivaroxaban, edoxaban, and dabigatran were identified and examined. Linkage disequilibrium with HHT genes was excluded, before variants within these genes were examined following whole genome sequencing of general and HHT populations. The 44 genes identified included 5/17 actionable pharmacogenes with guidelines. The 76,156 participants in the Genome Aggregation Database v3.1.2 had 28,446 variants, including 9668 missense substitutions and 1076 predicted loss-of-function (frameshift, nonsense, and consensus splice site) variants, i.e., approximately 1 in 7.9 individuals had a missense substitution, and 1 in 71 had a loss-of-function variant. Focusing on the 17 genes relevant to usually preferred DOACs, similar variant profiles were identified in HHT patients. With HHT patients at particular risk of haemorrhage when undergoing anticoagulant treatment, we explore how pre-emptive pharmacogenomic testing, alongside HHT gene testing, may prove beneficial in reducing the risk of bleeding and conclude that HHT patients are well placed to be at the vanguard of personalised prescribing.