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. This evidence concerns the gene PTPN2 and myeloid sarcoma.