Extended Data Figure 3. Repertoire changes with age and isotpe usages with disease.
a) Correlation of p-values obtained using age as a covariate versus excluding age in the analysis across 178 BCR features (calculated by two-sided by ANOVA), and b) analyses where statistical significance was discordant (i.e. below significance threshold without accounting for age and above significance threshold while using age as a covariate, or vice versa, purple points from (a)). c) Analyses where statistically significant p-values were decreased by >1.5 fold after using age as a covariate. Grey dotted lines in (a) indicate the threshold of significance after account for multiple testing (FDR<0.05, determined by Šidák method). d) The percentages of normalized isotype usages (unique VDJ sequences per isotype) for PBMC BCR repertoires per disease. e) Normalized IgE immunoglobulin constant region transcript levels between disease groups from transcriptomic data. n=58, 23, 33, 13, 10, 8, 11 and 37 for healthy, AAV MP0+, AAV PR3+, Behçhet’s, CD, EGPA, IgAV and SLE patients respectively. f) IgE titre between healthy individuals (n=4) and EGPA patients (n=5). P-values calculated by two-sided by ANOVA for (d)-(e) and * denotes FDR <0.05, ** <0.005, *** <0.0005, where FDR determined by Šidák method. Boxplots show the 25th, 50th and 75th percentiles; whiskers show upper and lower quartiles.