EMBO/soda-nel huggingface dataset:
USP34 controls osteogenic differentiation and bone formation by regulating BMP2 signaling (2018)
Figure 4

Figure 4. Loss of USP34 blunts BMP2-induced responses

A Scatter plot of RNA-seq expression analysis. 508 genes were upregulated and 391 genes were downregulated. Human MSCs were transfected with USP34 or control siRNAs and cultured in osteogenic medium for 2 days. Two biological replicates per group.

B KEGG pathway analysis indicated the altered function of TGF-β/BMP signaling pathway.

C Heatmap of representative genes associated with BMP2 signaling pathway.

D Relative activity of BRE luciferase assay. 293T cells were depleted with serum overnight followed by a treatment with BMP2 for 6 hours. Results are shown as mean ± SEM; n=3. ***: p<0.001 by t test.

E Immunoblot analysis of BMP2-induce Sp7 expression. Prx1-Cre;Usp34f/f and Usp34f/f control MSCs were starved overnight and then stimulated with 100 ng/ml BMP2 for 24 hours.

F Quantitative RT-PCR analysis of BMP2-induce Sp7 transcription. Cells were depleted from serum overnight before BMP2 stimulation. Results are shown as mean ± SEM; n=3; *: p<0.05 and **: p<0.01 by t test.

G Representative images of ALP and ARS staining. MSCs were cultured with normal medium supplemented with 100 ng/ml BMP2 for 10 and 21 days, respectively.

H, I Quantitative analyses of the ALP activity (n=4) and calcium mineralization (n=5). Results are shown as mean ± SEM; *: p<0.05, **: p<0.01 and ***: p<0.001 by t test.

J, K Representative microCT images and quantification of newly formed bone in response to BMP2. Results are shown as mean ± SEM; n=8. *: p<0.05 by t test. Scale bar, 2500 μm.

L H&E staining of BMP2-induced newly formed bone. Scale bar, 50 μm.

M Quantitative analysis of bone volume versus total tissue volume (BV/TV). Results are shown as mean ± SEM; n=8; **: p<0.01 by t test.

N H&E staining of calvarial organ cultures stimulated with or without 100 ng/ml BMP2. Red dotted lines indicate the interface between original and newly formed bone. Scale bar, 50 μm.

O Quantitative analysis of newly formed bone volume versus original bone volume (%). Results are shown as mean ± SEM; n=8; **: p<0.01 by t test.
USP34 controls osteogenic differentiation and bone formation by regulating BMP2 signaling (2018)
Figure 6

Figure 6. USP34 deubiquitinates and stabilizes Smad1 and RUNX2

A Co-immunoprecipitation of USP34 with ectopically expressed FLAG-tagged Smad1 in 293T cells.

B Immunoblot analysis showing decreased FLAG-Smad1 after USP34 knockdown.

C Immunoblot of Smad1-linked polyUb. 293T cells were co-transfected with FLAG-Smad1 overexpression plasmid and USP34 siRNA, and treated with 10 µM MG132 for 4 hours before collection.

D Immunoblot analysis of the degradation of FLAG-Smad1 protein. 293T cells were co-transfected with FLAG-Smad1 overexpression plasmid and USP34 siRNA, and treated with 50 µg/ml cycloheximide (CHX) before collection at designated time points. Right: relative quantification of FLAG-Smad1 protein levels at different time points.

E Co-immunoprecipitation of USP34 with ectopically expressed FLAG-tagged RUNX2 in 293T cells.

F Co-immunoprecipitation of USP34 with endogenous RUNX2 in mouse MSCs.

G Immunoblot of RUNX2-linked polyUb. Differentiated mouse MSCs were treated with 10 µM MG132 for 4 hours before collection. Prx1-Cre(+) indicates MSCs isolated from Prx1-Cre;Usp34f/f mice, and Prx1-Cre(-) indicates those from Usp34f/f control mice.

H Immunoblot analysis of the degradation of endogenous RUNX2 protein in MSCs. Prx1-Cre;Usp34f/f and Usp34f/f control MSCs were stimulated with BMP2, treated with 50 µg/ml CHX, and collected at designated time points. Right: relative quantification of RUNX2 protein levels at different time points.

I Representative images of ALP and ARS staining. Usp34f/f and Prx1-Cre;Usp34f/f MSCs were transduced with GFP, Smad1, Runx2 or Smad1+Runx2 adenoviruses, and then cultured in osteogenic medium for 5 and 10 days, respectively.

J, K Quantitative analyses of the ALP activity and calcium mineralization. Results are shown as mean ± SEM; n=8; *: p<0.05, **: p<0.01 and ***: p<0.001 by ANOVA with Tukey's post hoc test.

L, M H&E staining and quantitative analysis of MSCs-mediated ectopic bone formation. Scale bar, 50 μm. Results are shown as mean ± SEM; n=8; *: p<0.05 and **: p<0.01 by ANOVA with Tukey's post hoc test.
USP34 controls osteogenic differentiation and bone formation by regulating BMP2 signaling (2018)
Figure 7

Figure 7. Depletion of Smurf1 restores the BMP2-induced responses and osteogenic potential of Usp34-deficient MSCs

A Relative activity of BRE luciferase assay. 293T cells were depleted with serum overnight followed by a treatment with BMP2 for 6 hours. Results are shown as mean ± SEM; n=3. *: p<0.05 by ANOVA with Tukey's post hoc test.

B Immunoblot analysis of BMP2-induce Sp7 expression. Prx1-Cre;Usp34f/f MSCs were treated with Smurf1 or control siRNAs, starved overnight and then stimulated with 100 ng/ml BMP2 for 24 hours.

C Quantitative RT-PCR analysis of BMP2-induce Sp7 transcription. Prx1-Cre;Usp34f/f MSCs were treated with Smurf1 or control siRNAs, and depleted from serum overnight before BMP2 stimulation. Results are shown as mean ± SEM; n=3; *: p<0.05 and **: p<0.01 by ANOVA with Tukey's post hoc test.

D ALP staining and quantitative measurement of ALP activity. MSCs were cultured with normal medium supplemented with 100 ng/ml BMP2 for 5 days. Results are shown as mean ± SEM; n=8; **: p<0.01 by ANOVA with Tukey's post hoc test.

E Representative images of ARS staining and quantification. Cells were cultured with normal medium supplemented with 100 ng/ml BMP2 for 14 days. Results are shown as mean ± SEM; n=8; *: p<0.05, and **: p<0.01 by ANOVA with Tukey's post hoc test.

F Representative images of ALP and ARS staining. Mouse MSCs were treated with Smurf1 or control siRNAs, and cultured in osteogenic medium for 5 and 10 days, respectively.

G, H Quantitative analyses of the ALP activity and calcium mineralization. Results are shown as mean ± SEM; n=8; *: p<0.05, **: p<0.01 and ***: p<0.001 by ANOVA with Tukey's post hoc test.

I, J Quantitative RT-PCR and western blot analyses of the expression of osteogenic markers. Cells were treated with Smurf1 or control siRNAs, and cultured in osteogenic medium for 5 days. Results are shown as mean ± SEM; n=3; *: p<0.05, **: p<0.01 and ***: p<0.001 by ANOVA with Tukey's post hoc test.
A novel thermoregulatory role for PDE10A in mouse and human adipocytes (2016)
Figure 1

Figure 1. PDE10A is expressed in brown adipose tissue of mice

Representative PET, MRI and fused PET/MRI images of the thoracic region of a mouse that received the PDE10A radioligand [18F]-AQ28A. Interscapular brown adipose tissue (BAT) in each image is highlighted by the crosshairs. The mean standardized uptake (SUV) value of [18F]-AQ28A in BAT and skeletal muscle throughout the dynamic scan was calculated and compared for each animal (n = 5). * P = 0.0488 (20 mins), * P = 0.0236 (25 mins), * P = 0.0167 (30 mins), ** P = 0.00301 (35 mins), ** P = 0.00269 (40 mins), ** P = 0.00215 (45 mins), ** P = 0.00129 (50mins) and ** P = 0.00141 (45mins) using unpaired 2-tailed Student’s t tests.

Data information: Data are represented as mean +/- SEM.
A novel thermoregulatory role for PDE10A in mouse and human adipocytes (2016)
Figure 2

Figure 2. Acute pharmacological inhibition of PDE10A with MP-10 stimulates in vivo glucose uptake by brown adipose tissue and thermogenesis in lean mice

A. Representative PET images of a sagittal view of the thoracic region of fasted lean mice treated with either MP-10 (30mg/kg) (n = 3) or vehicle (n = 4) prior to receiving [18F]-FDG. Interscapular brown adipose tissue (BAT) in each image is highlighted by the crosshairs. The mean standardized uptake value (SUV) of [18F]-FDG in BAT after both treatments was calculated for each group. ** P = 0.0013 using unpaired 2-tailed Student’s t test.

B. The percentage (%) change in core body temperature and associated area under the curve (AUC) of mice treated with MP-10 (30mg/kg) (n = 4) or vehicle (n=4) prior to a cold challenge study for 4 hours at 8oC. * P = 0.0108 (1st hour), * P = 0.0218 (2nd hour), * P = 0.0121 (3rd hour), * P = 0.0168 (4th hour) and * P = 0.0264 (AUC) using unpaired 2-tailed Student’s t tests.

C. RT-qPCR analysis of relative mRNA expression of Pde10a in BAT, peri-ovarian visceral white adipose tissue (VAT) and inguinal subcutaneous white adipose tissue (SAT) of lean mice (n = 5). ** P = 0.0063 (SAT vs VAT), ** P = 0.0011 (SAT vs BAT) and ** P = 0.0037 (VAT vs BAT) using unpaired 2-tailed Student’s t test.

D Cyclic GMP (cGMP) and cyclic AMP (cGMP) concentrations normalized to protein concentration in SAT, VAT and BAT 30 minutes after treatment of lean mice with MP-10 (30 mg/kg) (n=6) or vehicle control (n = 6). ** P = 0.0085 (cGMP in VAT), ** P = 0.0011 (cGMP in BAT) and * P = 0.04 (cAMP in VAT) using unpaired 2-tailed Student’s t tests.

E RT-qPCR analysis of relative mRNA expression of thermogenic genes in SAT, VAT and BAT after treatment of lean mice treated with MP-10 (30 mg/kg) (n=7) or vehicle control (n = 6). **** P < 0.0001 (Pgc1alpha in VAT), ** P = 0.0012 (Ucp1 in VAT), ** P = 0.0071 (Pgc1alpha in BAT) and * P = 0.0486 (Ucp1 in BAT) using unpaired 2-tailed Student’s t tests.

Data information: Data are represented as mean +/- SEM. Graphs in A. (right panel) and B. (right panel) show individual data where bars represent mean +/- SEM.
A novel thermoregulatory role for PDE10A in mouse and human adipocytes (2016)
Figure 3

Figure 3. Expression of PDE10A in brown adipose tissue is increased in different mouse models of obesity

Representative fused PET/MRI images of the transverse view of the thoracic region of lean (n = 5), diet induced obese (DIO) (n = 4) and leptin deficient (ob/ob) (n = 5) mice that received the PDE10A radioligand [18F]-AQ28A. Interscapular brown adipose tissue (BAT) in each image is highlighted by the crosshairs. The mean standardized uptake (SUV) value of [18F]-AQ28A in BAT throughout the dynamic scan was calculated with associated area under the curve (AUC). * P = 0.0128 (lean vs DIO) and ** P = 0.0031 (lean vs ob/ob) using unpaired 2-tailed Student’s t test.

Data information: Data are represented as mean +/- SEM. Graph in right panel shows individual data where bars represent mean +/- SEM
A novel thermoregulatory role for PDE10A in mouse and human adipocytes (2016)
Figure 4

Figure 4. Weight loss from chronic inhibition of PDE10A with MP-10 is associated with browning of visceral white adipose tissue and improved insulin sensitivity in diet induced obese mice

A, B The percentage (%) change in body weight and cumulative food intake of diet induced obese (DIO) mice in response to daily treatment with MP-10 (10mg/kg) or vehicle for a week (n= 4 per group). * P = 0.0264 using (day 4), *** P = 0.0002 (day 5), *** P < 0.0001 (day 6) and *** P < 0.0001 (day 7) using Two-Way analysis of variance (ANOVA) with Sidak’s post-hoc test.

C. Upon completion of the chronic feeding study, an insulin tolerance test was performed on mice with single treatment of 0.75 U/kg insulin (n = 4 per group). * P = 0.0246 and * P = 0.0305 using unpaired 2-tailed Student’s t tests.

D, E, F Upon completion of the chronic feeding study, resting energy expenditure (REE), respiratory quotient (RQ) and locomotor activity of DIO mice treated with MP-10 (30mg/kg) compared to vehicle control were measured in metabolic cage experiments (n = 4 per group). * P = 0.0239 (02:00, RQ), * P = 0.0424 (02:30, RQ), * P = 0.0420 (03:00, RQ), * P = 0.0230 (03:30, RQ), * P = 0.0213 (04:00, RQ), ** P = 0.0116 (00:00, REE), ** P = 0.00823 (00:30, REE), ** P = 0.00969 (01:00, REE), ** P = 0.00274 (01:30, REE), ** P = 0.00323 (02:00, REE), ** P = 0.000441 (02:30, RE), ** P = 0.00390 (03:00, REE), ** P = 0.000211 (03:30, REE), ** P = 0.00216 (04:00, REE) and ** P = 0.0185 (04:30, REE) using unpaired 2-tailed Student’s t tests.

G. RT-qPCR analysis of relative mRNA expression of candidate genes in inguinal subcutaneous white adipose tissue (SAT), periovarian visceral white adipose tissue (VAT) and interscapular brown adipose tissue (BAT) (n = 8 per group) and immunostaining for UCP1 protein in *** P = 0.0004 (Pgc1alpha VAT), ** P = 0.0036 (Ucp1 VAT), ** P = 0.0068 (Cidea VAT), * P = 0.0411 (Tmem26 VAT) and ** P = 0.0070 (Pgc1alpha BAT) using unpaired 2-tailed Student’s t tests.

H. Representative sections of inguinal SAT, peri-ovarian VAT and interscapular BAT immuno-stained for UCP1 of mice chronically treated with MP-10 (10 mg/kg) or vehicle control. Scale bars, 100µm.

Data information: Data are represented as mean +/- SEM. Graph in B. shows individual data where bars represent mean +/- SEM
A novel thermoregulatory role for PDE10A in mouse and human adipocytes (2016)
Figure 5

Figure 5. PDE10A is expressed in supraclavicular brown tissue of humans

Representative PET image of the upper body of a human subject who received the PDE10A radioligand [18F]-MNI-1869. Supraclavicular brown adipose tissue (BAT) is highlighted by the crosshairs. The mean standardized uptake value (SUV) of [18F]-MNI-1869 in supraclavicular BAT and skeletal muscle was calculated in human subjects (n = 5). ** P = 0.00210 (5 mins), **** P < 0.0001 (15 mins), **** P < 0.0001 (30 mins), **** P < 0.0001 (50 mins) and ** P = 0.00462 (75 mins) using unpaired 2-tailed Student’s t tests.

Data information: Data are represented as mean +/- SEM.
A novel thermoregulatory role for PDE10A in mouse and human adipocytes (2016)
Figure 6

Figure 6. Pharmacological inhibition of PDE10A with MP-10 activates human brown adipocytes and causes browning of human white adipocytes

A. RT-qPCR analysis of relative mRNA expression of thermogenic genes in cultured primary human brown adipocytes in response to chronic treatments (8 hours) with DMSO (control), MP-10 (1nM and 100nM), cyclic AMP (cAMP – 10µM), cyclic GMP (cGMP - 10µM) and cAMP + cGMP (cA + cG – 10µM (n = 4 separate cultures). * P = 0.0278 vs (100nM MP-10 vs control; Pgc1alpha), **P = 0.0097 (10µM cA + cG vs control; Pgc1alpha), ** P = 0.0063 (1nM MP-10 vs control; Ucp1), **** P < 0.0001 (100nM MP-10 vs control; Ucp1), ## P = 0.0017 (1nM MP-10 vs 100nM MP-10; Ucp1), §§§§ P < 0.0001 (10µM cAMP vs 100nM MP-10; Ucp1), §§§§ P < 0.0001 (10µM cGMP vs 100nM MP-10; Ucp1), §§§§ P < 0.0001 (10µM cA + cG vs 100nM MP-10; Ucp1), * P = 0.0112 vs (100nM MP-10 vs control; Cidea), *** P = 0.0007 (1nM MP-10 vs control; Dio2), **** P < 0.0001 (100nM MP-10 vs control; Dio2), # P = 0.0267 (1nM MP-10 vs 10µM cAMP; Dio2), §§§ P = 0.0001 (10µM cAMP vs 100nM MP-10; Dio2), §§§ P = 0.0007 (10µM cGMP vs 100nM MP-10; Dio2), * P = 0.0232 (10µM cA + cG vs control; Dio2) §§ P = 0.0058 (10µM cA + cG vs 100nM MP-10; Dio2), * P = 0.0416 (10µM cGMP vs control; Prdm16) and * P = 0.0105 (10µM cA + cG vs control; Prdm16) using one-way analysis of variance (ANOVA) with Tukey’s post-hoc test.

B. RT-qPCR analysis of relative mRNA expression of thermogenic genes in cultured primary human white adipocytes in response to chronic treatments (8 hours) with DMSO (control), MP-10 (1nM and 100nM), cyclic AMP (cAMP – 10µM), cyclic GMP (cGMP - 10µM) and cAMP + cGMP (cA + cG – 10µM (n = 4 separate cultures). * P = 0.0186 vs (1nM MP-10 vs control; Ucp1), **** P < 0.0001 (100nM MP-10 vs control; Ucp1), ## P = 0.0016 (1nM MP-10 vs 100nM MP-10; Ucp1), *** P = 0.0003 (10µM cAMP vs control; Ucp1), §§§ P = 0.0002 (10µM cGMP vs 100nM MP-10; Ucp1) and *** P = 0.0003 (cA + cG vs control; Ucp1) using one-way analysis of variance (ANOVA) with Tukey’s post-hoc test.

Data information: Individual data from separate experiments are shown where bars represent mean +/- SEM.
Histone H3.3 promotes IgV gene diversification by enhancing formation of AID-accessible single stranded DNA (2016)
Figure 1

Figure 1. Reduced diversification of IGVL in cells lacking H3.3.

A. Ig sequence diversification in wild type and h3.3. Each line represents an IGVLR gene from a population cells that have expanded for 36 divisions after cloning. Blue bars = gene conversions; green lollipops = point mutations; green bars = ambiguous mutations (Sale et al, 2001).

B. Quantification of sequence diversification in wild type and two h3.3 clones expressed as unique events (gene conversions or point mutations) per kb per generation. (Sequences analysed: WT n = 142, for h3.3 c20 n = 140, for h3.3 c32 n = 190).

C. Quantification of endogenous AID mRNA in wild type and h3.3 cells. The graph represents the mean and SD of four independent RNA extractions of wild type and h3.3c20. p value calculated with the unpaired T test.

D. Fluctuation analysis for the formation of sIg-loss variants in wild type (black circles) and h3.3 clones (black diamonds) and in two wild type clones (blue circles), two h3.3 clones (blue diamonds) expressing hAIDup (Wang et al, 2009) at 15 - 20 times the level of wild type AID. The final column shows h3.3 cells with hAIDup complemented with H3.3 (blue triangles). Clones were expanded for 36 cell divisions at which point the percentage of sIgM negative cells in each population was determined by flow cytometry. Error bars represent mean and 1 S.D. P values calculated with the Mann-Whitney test.

E. Ig sequence diversification in wild type, h3.3 and complemented h3.3 cells overexpressing hAIDup. Key as in (A) plus * = deletion or duplication.

F. Mutation rate as a function of hAIDup expression in wild type and h3.3. Filled symbols represent clones overexpressing hAIDup, open symbols clones with natural AID expression. Linear regression lines are shown.

G. Spectrum of diversification in wild type and h3.3. GC = gene conversion; PM = point mutation; Amb = ambiguous; Other = deletions or duplications.
Histone H3.3 promotes IgV gene diversification by enhancing formation of AID-accessible single stranded DNA (2016)
Figure 2

Figure 2. The limited IgV diversification in h3.3 cells results from reduced activity of AID.

A. Scheme outlining the mechanism of Ig diversification in DT40 cells, indicating the action of uracil glycosylase inhibitor (UGI). UNG = Uracil DNA glycosylase.

B. Diversification in wild type and h3.3 cells expressing similar levels of hAIDup and UGI. Red lollipops = transitions at C/T; Black lollipops = other point mutations; blue bars = gene conversions.

C. Mutation rate as a function of hAIDup expression in UGI-expressing wild type and h3.3 cells. Each clone was expanded for 36 cell divisions and hAIDup expression determined by qPCR during the last week of culture.
Histone H3.3 promotes IgV gene diversification by enhancing formation of AID-accessible single stranded DNA (2016)
Figure 3

Figure 3. Loss of H3.3 does not affect transcription of IGVL.

A. Enrichment of H3.3 in the rearranged allele of the Ig light chain (IGVLR) compared with the unrearranged allele (IGVLUR). Each bar represents the mean of 3 independent ChIPs each assayed with 3 repeat qPCRs (**** = p<0.0001, unpaired T-test).

B. Igλ mRNA levels in wild type and h3.3 cells determined by qPCR, normalised to EF-1α expression and to WT. n = 3 (ns = not significant, unpaired T-test).

C. Surface Ig (sIg) in multiple clones of wild type (blue) and h3.3 (red) cells. Green = unstained control.

D. ChIP for H3K4me3 and H3K36me3 in IGVLR and IGVLUR. Each bar represents the mean of 3 ChIPs for H3K4me3 and IgG and 2 for H3K36me3, each assayed with 3 qPCR repeats (statistical analysis as in Panel A).

E. 4sUDRB mapping of transcription kinetics through IGVL. Representative locations of qPCR amplicons (blue) are shown relative to the TSS of IGVL (distance is shown in parenthesis). qPCR signals were normalised to the enrichment at time zero (before release). Error bars represent 1SD of qPCR triplicates.
Histone H3.3 promotes IgV gene diversification by enhancing formation of AID-accessible single stranded DNA (2016)
Figure 4

Figure 4. H3.3 promotes formation of R-loops in IGVL.

A. DNA-RNA hybrids in IGVL and IGCL in wild type and h3.3. The DRIP signal is normalised to total input and to wild type IgVλ (*** = p<0.001, ns = not significant, unpaired T-test).

B. MNase mapping reveals that loss of H3.3 does not affect average Igλ nucleosome positioning. The red bars indicate the qPCR amplicons for the MNase mapping. The qPCR signals are normalised to the first amplicon at the TSS. Error bars represent 1SD of 2 experimental repeats.

C. DNA-RNA hybrids in IGVL of wild type clones either without overexpression of RNAseHI (red, clones 1 and 2, each assayed with three qPCR repeats) or with overexpression of RNAseHI (blue, clones 3 and 4, each assayed with three qPCR repeats) (*** = p<0.001, ** = p < 0.01; unpaired T-test). Expression of hAIDup-FLAG and chicken RNAseHI-YFP in these clones is shown in Figure EV6.

D. Rate of Ig diversification in wild type cells with and without overexpression of RNAseHI.
Histone H3.3 promotes IgV gene diversification by enhancing formation of AID-accessible single stranded DNA (2016)
Figure 5

Figure 5. H3.3 promotes the formation of ssDNA stretches along the IgV locus

A. Location of ssDNA tracts ≥ 8nt in IGVL in wild type (WT) and h3.3 cells detected by bisulphite sequencing. Data is from three independent experiments, the results from each being shown as black, red and green bars respectively. ssDNA on the coding strand is shown above the marker line (numbered relative to the TSS at 0) and on the template strand below. In the centre, the percentage of G and C bases on the coding strand in a 30 bp rolling average window is shown.

B. Quantitation of the bisulphite detection of single stranded DNA in wild type, h3.3 and h3.3 complemented cells, normalised to wild type. p-values calculated with a chi‐square test. Sequences analysed: Wild type n = 322, h3.3 n = 320 and h3.3 + H3.3 n = 95.

C. Bisulphite detection of single stranded DNA in wild type cells with and without RNAseHI overexpression. p-value calculated with a chi‐square test. Sequences analysed: RNAseH1+ve n = 94, RNAseH1-ve n = 90.