In addition, we used the model of y =  α + β1xIL2RA +β 2xPTPN2 +β 3xIL2RA*PTPN2 + ε to assess the interaction of each variant in the control, T1D and MS cohorts shown in Figure 1. Here, PTPN2 is linked to type 1 diabetes mellitus.