Target enriched sequencing was designed to cover 356 IEI genes as previously described [35 (link)] and used to diagnose P1 and P6. WES was perfomed for P4, P5, and P12 as described [24 (link), 35 (link)]. For Sanger sequencing of IL21R in P2 and P3, genomic DNA was extracted from whole blood (QIamp DNA Blood Kit, Qiagen). Targeted Sanger sequencing of IL21R in P7–P10 has been reported [22 (link), 24 (link)]. Primer sequences used to identify IL21R variants in P2, P3, P11, and P13 and validate NGS findings in P1, P4, P5, and P6 are described in Supplementary Table I. Pathogenicity of genetic variants was evaluated based on the allele frequency in public databases including ExAC and gnomAD and prediction tools such as Polyphen-2, SIFT and the combined annotation dependent depletion (CADD) score.