No PAR genes were enriched in the analysis although p.Pro500Ala (X-1601004:C-G) in AKAP17A, linked to autoimmunity and showed non-adjusted significance in the variant level analysis (variant allele frequency 0.75 in autoimmune and 0.39 in non-autoimmune; effect size 0.36; gnomAD males 0.54; p-value <0.01, adj.p-value 1.0) (Figure 8; Supplementary Table 2). Here, AKAP17A is linked to Autoimmunity.