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

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 sub-categorized pathogenic DNA variants in ENG/endoglin, ACVRL1/ALK1, and SMAD4 if they generated premature termination codons (PTCs) subject to nonsense mediated decay. In three pre-phenotyped patient cohorts, a PTC-based classification system explained some previously puzzling hemorrhage variability. In blood outgrowth endothelial cells (BOECs) derived from ACVRL1+/PTC, ENG+/PTC, and SMAD4+/PTC patients, PTC-containing RNA transcripts persisted at low levels (8-23% expected, varying between replicate cultures); genes differentially expressed to Bonferroni p<0.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 '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 (ATF)4 which controls expression of stress-adaptive genes, and showed that ENG Q436X but not ENG R93X directly induced ATF4. AlphaFold accurately modelled 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.

Genome-wide assessment of genetic risk for systemic lupus erythematosus and disease severity

Using three European and two Chinese genome-wide association studies (GWAS), we investigated the performance of genetic risk scores (GRSs) for predicting the susceptibility and severity of systemic lupus erythematosus (SLE), using renal disease as a proxy for severity. We used four GWASs to test the performance of GRS both cross validating within the European population and between European and Chinese populations. The performance of GRS in SLE risk prediction was evaluated by receiver operating characteristic (ROC) curves. We then analyzed the polygenic nature of SLE statistically. We also partitioned patients according to their age-of-onset and evaluated the predictability of GRS in disease severity in each age group. We found consistently that the best GRS in the prediction of SLE used SNPs associated at the level of P < 1e−05 in all GWAS data sets and that SNPs with P-values above 0.2 were inflated for SLE true positive signals. The GRS results in an area under the ROC curve ranging between 0.64 and 0.72, within European and between the European and Chinese populations. We further showed a significant positive correlation between a GRS and renal disease in two independent European GWAS (P_cohort1 = 2.44e−08; P_cohort2 = 0.00205) and a significant negative correlation with age of SLE onset (P_cohort1 = 1.76e−12; P_cohort2 = 0.00384). We found that the GRS performed better in the prediction of renal disease in the ‘later onset’ compared with the ‘earlier onset’ group. The GRS predicts SLE in both European and Chinese populations and correlates with poorer prognostic factors: young age-of-onset and lupus nephritis.