GSE182463
Biallelic PAX5 mutations cause hypogammaglobulinemia,sensorimotor deficits and autism spectrum disorder
To investigate the molecular basis for the B cell developmental arrest in Pax5R31Q/– mice, we performed RNA-sequencing (RNA-seq) with ex vivo sorted Pax5+/+ and Pax5R31Q/– pro-B cells. We identified Pax5-activated and Pax5-repressed genes that were no longer properly regulated by the Pax5-R31Q protein in Pax5R31Q/– pro-B cells. Notably, these genes were only a subset of the activated and repressed genes identified by comparing Pax5+/+ and Pax5–/– pro-B cells, which raised the question whether binding of the Pax5-R31Q protein may be selectively lost at the subset of deregulated genes in Pax5R31Q/– pro-B cells. By using a Pax5 paired domain antibody for chromatin immunoprecipitation coupled with deep sequencing (ChIP-seq) of short-term cultured Pax5+/+ and Pax5R31Q/– pro-B cells, we identified all associated Pax5 binding sites. Common Pax5 peaks had a similar Pax5-binding density in contrast to the strong binding difference observed at the unique peaks present in Pax5+/+ pro-B cells. To investigate a correlation between the loss of Pax5 binding and gene expression, we focused our analysis on Pax5 peaks in the TSS region of activated genes. We systematically investigated the correlation between loss of Pax5 binding at the TSS and down-regulation of gene expression in Pax5R31Q/– pro-B cells. The ratio of Pax5 binding between Pax5R31Q/– and Pax5+/+ pro-B cells at the TSS of these activated genes was significantly reduced compared to that of expressed non-regulated genes. We next explored whether the Pax5-binding difference at the TSS also correlated with the magnitude of gene expression difference. The loss of Pax5 binding at the TSS also correlated with the degree of expression change in Pax5R31Q/– pro-B cells compared to Pax5+/+ pro-B cells. We conclude therefore that the selective DNA-binding of Pax5-R31Q is responsible for the observed gene expression differences in Pax5R31Q/– pro-B cells. Overall design: 8 RNA-Seq samples, 4 different conditions (A,B,C,D) two replicates each. A (wildtype): Pro-B cells Pax5(+/+)  111735 + 111736; B (R31Q mutant): Pro-B cells pax5(R31Q/-) 111739 + 111740; C (Pax5 knock out): Progenitor cells Pax5(fl/fl) VAV-cre 111726 + 111741; D (Pro-B cells Pax5 het): Pro-B cell Pax5(+/-) 111737 + 111738. 2 ChIP-Seq samples (65448 (wildtype) + 65449 (mutant)).