Whole-genome sequencing identifies rare missense variants of WNT16 and ERVW-1 causing the systemic lupus erythematosus

Systemic lupus erythematosus (SLE, OMIM 152700) is a rare autoimmune disease with high heritability that affects ~0.1% of the population. Previous studies have revealed several common variants with small effects in European and East Asian SLE patients. However, there is still no rare variant study on Chinese SLE patients using the whole-genome sequencing technology (WGS). Here, we designed a family based WGS study to identify novel rare variants with large effects. Based on large-scale allele frequency data from the gnomAD database, we identified rare protein-coding gene variants with disruptive and sequence-altering impacts in SLE patients. We found that the burden of rare variants was significantly higher than that of common variants in patients, suggesting a larger effect of rare variants on the SLE pathogenesis. We identified the pathogenic risk of rare missense variants with significant odds ratios (p < 0.05) in two genes, including WNT16 (NC_000007.14:g.121329757G > C, NP_057171.2:p.(Ala86Pro) and 7 g.121329760G > C, NP_057171.2:p.(Ala87Pro)), which explains five out of seven patients covering all three families but are absent from all controls, and ERVW-1 (NC_000007.14:g.92469882A > G, NP_001124397.1:p.(Leu167Pro), rs74545114; NC_000007.14:g.92469907G > A, NP_001124397.1:p.(Arg159Cys), rs201142302; NC_000007.14:g.92469919G > A, NP_001124397.1:p.(His155Tyr), rs199552228), which explains the other two patients. None of these variants were identified in any of the controls. These associations are supported by known gene expression studies in SLE patients based on literature review. We further tested the wild and mutant types using the luciferase assays and qPCR in cells. We found that WNT16 can activate the canonical Wnt/β-catenin pathway while the mutant cannot. Additionally, the wild ERVW-1 expression can be significantly up-regulated by cAMP while the mutant cannot. Our study provides the first direct genetic and in vitro evidence for the pathogenic risk of mutant WNT16 and ERVW-1, which may facilitate the design of precision therapy for SLE.