Directional next-generation RNA sequencing and examination of premature termination codon mutations in endoglin/hereditary haemorrhagic telangiectasia

Hereditary haemorrhagic telangiectasia (HHT) is a disease characterised by abnormal vascular structures, and most commonly caused by mutations in ENG, ACVRL1 or SMAD4 encoding endothelial cell-expressed proteins involved in TGF-β superfamily signalling. The majority of mutations reported on the HHT mutation database are predicted to lead to stop codons, either due to frameshifts or direct nonsense substitutions. The proportion is higher for ENG (67%) and SMAD4 (65%) than for ACVRL1 (42%), p < 0.0001. Here, by focussing on ENG, we report why conventional views of these mutations may need to be revised. Of the 111 stop codon-generating ENG mutations, on ExPASy translation, all except one were premature termination codons (PTCs), sited at least 50-55 bp upstream of the final exon-exon boundary of the main endoglin isoform, L-endoglin. This strongly suggests that the mutated RNA species will undergo nonsense-mediated decay. We provide new in vitro expression data to support dominant negative activity of stable truncated endoglin proteins but suggest these will not generate HHT: the single natural stop codon mutation in L-endoglin (sited within 50-55 nucleotides of the final exon-exon boundary) is unlikely to generate functional protein since it replaces the entire transmembrane domain, as would 8 further natural stop codon mutations, if the minor S-endoglin isoform were implicated in HHT pathogenesis. Finally, next-generation RNA sequencing data of 7 different RNA libraries from primary human endothelial cells demonstrate that multiple intronic regions of ENG are transcribed. The potential consequences of heterozygous deletions or duplications of such regions are discussed. These data support the haploinsufficiency model for HHT pathogenesis, explain why final exon mutations have not been detected to date in HHT, emphasise the potential need for functional examination of non-PTC-generating mutations, and lead to proposals for an alternate stratification system of mutational types for HHT genotype-phenotype correlations.