A TOPBP1 allele causing male infertility uncouples XY silencing dynamics from sex body formation

  1. Carolline Ascenção
  2. Jennie R Sims
  3. Alexis Dziubek
  4. William Comstock
  5. Elizabeth A Fogarty
  6. Jumana Badar
  7. Raimundo Freire
  8. Andrew Grimson
  9. Robert S Weiss
  10. Paula E Cohen  Is a corresponding author
  11. Marcus B Smolka  Is a corresponding author
  1. Department of Molecular Biology and Genetics, Weill Institute for Cell and Molecular Biology, Cornell University, United States
  2. Fundación Canaria del Instituto de Investigación Sanitaria de Canarias (FIISC), Unidad de Investigación, Hospital Universitario de Canarias, Spain
  3. Instituto de Tecnologías Biomédicas, Universidad de La Laguna, Spain
  4. Universidad Fernando Pessoa Canarias, Spain
  5. Department of Biomedical Sciences, Cornell University, United States
11 figures, 2 tables and 1 additional file

Figures

Figure 1 with 1 supplement
A new TOPBP1 mutant separating its role in fertility from organismal viability.

(A) Schematic showing mutations in the Topbp1 B5 allele. (B) Body mass (Topbp1+/+ mean = 25.26, SD = 2.38; Topbp1B5/B5 mean = 26.43, SD = 2.28, n = 9) and (C) appearance of Topbp1B5/B5 and Topbp1+/+ littermate mice. (D) Effect of full body ionizing radiation (IR) (7 Gy) on changes in body mass of Topbp1B5/B5 and Topbp1+/+ littermate mice (E) Breeding scheme and resulting litters. (F, G) Comparison of testes size (Topbp1+/+ mean = 0.038, SD = 0.006; Topbp1B5/B5 mean = 0.011, SD = 0.004, n = 9), and (H) sperm count, of Topbp1B5/B5 and Topbp1+/+ littermate mice (Topbp1+/+ mean = 2.1 × 107, SD = 6 × 106; Topbp1B5/B5 mean = 0.0, SD = 0.0, n = 9). (I) H&E-stained histological testes sections displaying loss of cellularity in Topbp1B5/B5. Green arrow = spermatogonia, red arrow = healthy spermatocyte, blue arrow = dying spermatocyte. (J, K) TUNEL assay performed on histological testes sections (Topbp1+/+ mean = 0.36, SD = 0.15; Topbp1B5/B5 mean = 4.80, SD = 0.43, n = 3). (L) Meiotic spreads stained for SYCP3, SYCP1, and H1t. ****p<0.0001, n = number of mice. p-Values were calculated using unpaired t-test.

Figure 1—figure supplement 1
Genotyping of Topbp1B5/B5 mice.

Sanger sequencing confirms the eight-charge reversal point mutations on the BRCT 5 domain of TOPBP1 from homozygous Topbp1B5/B5 mice.

Figure 2 with 4 supplements
Topbp1B5/B5 mouse embryonic fibroblasts (MEFs) display no detectable DNA damage response (DDR) defects.

(A) MEFs obtained from littermate mice of indicated genotypes were subjected to clonogenic survival assay in the indicated concentrations of hydroxyurea or (C) phleomycin. (B, D) Quantification of clonogenic survival assays from (A) and (C), respectively, performed in biological and experimental triplicates (n = 3), error bars are displayed as standard deviation. (E) Western blot for indicated DDR markers in MEFs obtained from Topbp1B5/B5 and Topbp1+/+ littermate mice. The data from MEFs were performed using littermate pairs and validated using a second pair of Topbp1B5/B5 and Topbp1+/+ littermate mice. (F) Immunofluorescence of MDC1 and phosphoH2AX_S139-stained nuclei from Topbp1B5/B5 and Topbp1+/+ MEFs treated with ionizing radiation (IR) (7 Gy).

Figure 2—source data 1

Original file for the western blot analysis in Figure 2E (anti-CHK1, anti-CHK2, and anti-GAPDH).

https://cdn.elifesciences.org/articles/90887/elife-90887-fig2-data1-v1.zip
Figure 2—figure supplement 1
Survival assays of Topbp1B5/B5 mouse embryonic fibroblasts (MEFs).

(A) MEFs obtained from littermate mice (second independent litter) of indicated genotypes were subjected to clonogenic survival assay in the indicated concentrations of hydroxyurea or (B) phleomycin. The clonogenic survival assays were performed in biological and experimental triplicates (n = 3), error bars are displayed as standard deviation.

Figure 2—figure supplement 2
Absence of an increased number of micronuclei in Topbp1B5/B5 mouse embryonic fibroblasts (MEFs).

Immunofluorescence showing DAPI-stained nuclei and micronuclei quantification of Topbp1+/+ and Topbp1B5/B5 (A) MEFs #1 and (B) MEFs #2. (MEFs #2 is a second independent litter.) Ns = p>0.05. p-Values were calculated using unpaired t-test.

Figure 2—figure supplement 3
TOPBP1-B5 shows disruption in binding 53BP1 and BLM.

(A) HEK-293T cells were transfected with plasmids expressing wild-type TOPBP1 or indicated BRCT mutants (K155A – BRCT 1, K250A – BRCT 2, K704A – BRCT 5, K1713A – BRCT 7, and B5 – 8 mutations in BRCT 5), with a 3xFlag tag. After cell lysis and anti-Flag IP, lysate and elution were analyzed by western blot for the indicated proteins. (B) IP-mass spectrometry of flagTOPBP1 wild type versus flag-TOPBP1-B5 was performed in silac-labeled HEK-293T cells. Cell transfections, lysis, protein IP, and elution were performed as described in (A). From left to right, the TOPBP1 interactors were plotted from highest to lowest binding. The size of circles represents the number of peptide hits.

Figure 2—figure supplement 4
Topbp1B5/B5 mouse embryonic fibroblasts (MEFs) show lower TOPBP1 protein abundance than Topbp1+/+ MEFs.

Western blot probing TOPBP1 in MEFs, Topbp1 +/+, Topbp1B5/+, and Topbp1B5/B5 in two biological replicates.

Figure 2—figure supplement 4—source data 1

Original file for the western blot analysis in Figure 2—figure supplement 4 (anti-TOPBP1 and anti-TUBULIN).

https://cdn.elifesciences.org/articles/90887/elife-90887-fig2-figsupp4-data1-v1.zip
Figure 2—figure supplement 4—source data 2

PDF containing Figure 2—figure supplement 4 and original scans of the relevant western blot analysis (anti-TOPBP1 and anti-TUBULIN) with highlighted bands and labels.

https://cdn.elifesciences.org/articles/90887/elife-90887-fig2-figsupp4-data2-v1.pdf
Figure 3 with 5 supplements
Markers of DNA repair, synapsis, sex body formation, and TOPBP1 localization are mostly normal in Topbp1B5/B5 spermatocytes.

(A) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 spermatocytes stained with SYCP3 and γH2AX, prepared as described in ‘Materials and methods’. (B) Quantification of γH2AX-stained pachytene spreads, upper graph XY body (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 386, number of mice = 8; Topbp1B5/B5 number of cells = 410, number of mice = 8; p-value=0.3063), bottom graph autosomal chromosomes (each dot represents the average of signal from all autosomes in each mouse, Topbp1+/+ number of cells = 160, number of mice = 8; Topbp1B5/B5 number of cells = 161, number of mice = 8; p-value=0.5081). (C) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and RAD51. (D) Quantification of RAD51 foci/cell of mid-pachytene meiotic spreads (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 149, number of mice = 3; Topbp1B5/B5 number of cells = 183, number of mice = 3; p-value=0.2174). (E) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and SYCP1. (F) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and TOPBP1. (G) Quantification of TOPBP1 on X and Y chromosome cores from (F) (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 246, number of mice = 3; Topbp1B5/B5 number of cells = 233, number of mice = 3; p-value=0.8546). (H) Quantification of TOPBP1 on X and Y chromatin loops from (F) (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 246, number of mice = 3; Topbp1B5/B5 number of cells = 233, number of mice = 3; p-value=0.6755). n = number of mice. p-Values were calculated using a linear mixed effect model (see ‘Materials and methods’ for details).

Figure 3—figure supplement 1
Topbp1B5/B5 mid-pachytene spermatocytes show normal HORMAD1, HORMAD2, pHORMAD2, and pCHK1_S345 localization and intensities.

(A) Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with HORMAD1 and SYCP3, prepared as described in ‘Materials and methods’. (B) Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with HORMAD2 and SYCP3. (C) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and pMDC1-T4 prepared as described in ‘Materials and methods’. (D) Quantification of pMDC1-T4-stained pachytene spreads from (C) (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 60, number of mice = 3; Topbp1B5/B5 number of cells = 159, number of mice = 3; p-value=0.8193). (E) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and pCHK1-S345. (F) Quantification of pCHK1-S345-stained pachytene spreads from (E) (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 168, number of mice = 2; Topbp1B5/B5 number of cells = 125, number of mice = 2; p-value=0.6536). (G) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and pHORMAD2-S271. (H) Quantification of pHORMAD2-S271-stained pachytene spreads from (G) (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 224, number of mice = 3; Topbp1B5/B5 number of cells = 254, number of mice = 3; p-value=0.2811). p-Values were calculated using a linear mixed effect model (see ‘Materials and methods’ for details).

Figure 3—figure supplement 2
Topbp1B5/B5 mid-pachytene spermatocytes show an increased number of MLH1 foci/cell but no difference in the number of MHL3 foci/cell.

(A) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and MLH1, prepared as described in ‘Materials and methods’ on the left, and on the right, the foci count of MLH1-stained pachytene spreads (each dot represents the number of foci/cell; Topbp1+/+ number of cells = 84, number of mice = 4; Topbp1B5/B5 number of cells = 94, number of mice = 4; p-value<0.0001). (B) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and MLH3 prepared as described in ‘Materials and methods’ on the left, and on the right, the Foci count of MLH3-stained pachytene spreads (each dot represents the number of foci/cell; Topbp1+/+ number of cells = 48, number of mice = 3; Topbp1B5/B5 number of cells = 91, number of mice = 3; p-value=0.4088). p-Values were calculated using a linear mixed effect model (see ‘Materials and methods’ for details).

Figure 3—figure supplement 3
Topbp1B5/B5 mid-pachytene spermatocytes form a normal sex body.

Topbp1+/+ and Topbp1B5/B5 mid-pachytene spermatocytes stained with DAPI and SYCP3, prepared as described in ‘Materials and methods’.

Figure 3—figure supplement 4
Topbp1B5/B5 mid-pachytene spermatocytes show subtle defects/delays in the spreading of γH2AX on XY chromosomes.

Quantification and representative images of Topbp1B5/B5 mid-pachytene spermatocytes showing: (A) normal spreading of γH2AX on XY chromosomes; (B) defect/delay in the spreading of γH2AX at the PAR; (C) defect/delay in the spreading of γH2AX at the pericentromeric region of the Y chromosome; (D) spreading of γH2AX into an autosomal chromosome fused to the X chromosome; and (E) other types of defects/delays in the spreading of γH2AX onto the XY chromosomes. Each dot in the graph represents the sum of the cells counted per mouse. (F) Percentage of the previously described defects found in Topbp1+/+ or (G) Topbp1B5/B5 mid-pachytene spermatocytes. *p<0.05, **p<0.01, ***p<0.001, and ****p<0.0001. p-Values were calculated using unpaired t-test in GraphPad. n = number of mice.

Figure 3—figure supplement 5
TOPBP1 is normally localized in all stages of prophase I in Topbp1B5/B5.

Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 spermatocytes stained with SYCP3 and TOPBP1, prepared as described in ‘Materials and methods’.

Figure 4 with 1 supplement
Topbp1B5/B5 spermatocytes display normal localization of ATR, BRCA1, and MDC1.

(A) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and ATR. (B) Quantification of ATR in pachytene spreads from (A). Left: ATR on X and Y chromosome cores (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 127, number of mice = 3; Topbp1B5/B5 number of cells = 127, number of mice = 3; p-value=0.4068). Right: ATR on X and Y chromatin loops (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 127, number of mice = 3; Topbp1B5/B5 number of cells = 127, number of mice = 3; p-value=0.9396). (C) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and BRCA1. (D) Quantification of BRCA1 in pachytene spreads from (C) (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 152, number of mice = 3; Topbp1B5/B5 number of cells = 140, number of mice = 3; p-value=0.6509). (E) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and MDC1. (F) Quantification of MDC1 in pachytene spreads from (E) (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 696, number of mice = 8; Topbp1B5/B5 number of cells = 988, number of mice = 8; p-value=0.3603). (G) Meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and pCHK1-S317. (H) Quantification of pCHK1-S317 in pachytene spreads from (G) (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 223, number of mice = 3; Topbp1B5/B5 number of cells = 254, number of mice = 3; **p-value=0.0023). p-Values were calculated using a linear mixed effect model (see ‘Materials and methods’ for details). (I) Table summarizing the normal or disrupted ATR and TOPBP1-dependent events during male fertility accessed in Topbp1B5/B5. n = number of mice.

Figure 4—figure supplement 1
53BP1 localization in Topbp1B5/B5 mid-pachytene spermatocytes.

Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and 53BP1 on the left, and quantification of 53BP1-stained pachytene spreads on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 95, number of mice = 2; Topbp1B5/B5 number of cells = 171, number of mice = 2; p-value=0.2384). p-Values were calculated using a linear mixed effect model (see ‘Materials and methods’ for details).

Defective Senataxin (SETX) phosphorylation and localization in Topbp1B5/B5 spermatocytes.

(A) Scatter plot of phosphoproteomic datasets corresponding to Topbp1+/+/Topbp1B5/B (Y axis) and Topbp1+/+(vehicle)/Topbp1+/+(AZ20) (X axis) from whole testes of mice. (B) SETX phosphopeptides identified in the Topbp1+/+/Topbp1B5/B phosphoproteomic experiment shown in (A). Red: reduced in Topbp1B5/B mutant; blue: unchanged. (C) Meiotic spreads showing pachytene spermatocytes from Topbp1+/+ and Topbp1B5/B5 mice stained with SYCP3 and SETX in regular immunofluorescence. (D) 3D-SIM analysis of meiotic spreads described in (C). (E) Quantification of SETX on X and Y chromatin loops in pachytene spreads from (C) (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 152, number of mice = 4; Topbp1B5/B5 number of cells = 174, number of mice = 4; *p-value=0.0452). (F) Quantification of SETX on X and Y chromosome cores in pachytene spreads from (C) (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 152, number of mice = 4; Topbp1B5/B5 number of cells = 174, number of mice = 4; p-value=0.5987). p-Values were calculated using a linear mixed effect model (see ‘Materials and methods’ for details).

Figure 6 with 16 supplements
Single-cell RNAseq reveals that Topbp1B5/B5 spermatocytes initiate meiotic sex chromosome inactivation (MSCI) but fail to promote full silencing.

(A) scRNAseq workflow for isolation and purification of single cells for RNAseq. (B) Uniform Manifold Approximation and Projection (UMAP) analysis of sub-clusters captured in the scRNAseq, representing all cells captured for both Topbp1+/+ and Topbp1B5/B5. (C) UMAP analysis of sub-clusters captured in the scRNAseq of Topbp1+/+. (D) UMAP analysis of sub-clusters captured in the scRNAseq of Topbp1+/+ (gray) and Topbp1B5/B5 (blue). (E) Violin plots displaying the ratio of the average expression of X chromosome genes by the average expression of chromosome 9 genes at different stages of spermatogenesis for Topbp1+/+ and Topbp1B5/B5 cells. The level of X-genes expression in spermatocyte 3 is significantly higher in Topbp1B5/B5 cells when compared to Topbp1+/+ cells, with a p-value of 1.5e-178 using a two-sided Wilcoxon rank-sum test.

Figure 6—figure supplement 1
Topbp1B5/B5 pachytene spermatocytes do not exhibit any difference in CHD4, Sumo2-3, and USP7 intensities and localization compared to Topbp1+/+.

(A) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and CHD4 on the left, and quantification of CHD4-stained pachytene spreads on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 176, number of mice = 4; Topbp1B5/B5 number of cells = 215, number of mice = 4; p-value=0.2125). (B) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and SUMO2.3 on the left, and quantification of SUMO2.3-stained pachytene spreads on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 155, number of mice = 2; Topbp1B5/B5 number of cells = 126, number of mice = 2; p-value=0.8951). (C) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and USP7 on the left, and quantification of USP7-stained pachytene spreads on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 73, number of mice = 2; Topbp1B5/B5 number of cells = 48, number of mice = 2; p-value=0.8610). p-Values were calculated using a linear mixed effect model (see ‘Materials and methods’ for details).

Figure 6—figure supplement 2
Topbp1B5/B5 pachytene spermatocytes show no difference in H3K9ac or K3K9me3 staining compared to Topbp1+/+.

(A) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and H3K9ac on the left, and quantification of H3K9ac-stained pachytene spreads on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 214, number of mice = 3; Topbp1B5/B5 number of cells = 203, number of mice = 3; p-value=0.4948). (B) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and H3K9me3 on the left, and quantification of H3K9me3-stained pachytene spreads on the right. (Early-pachytene spreads: each dot represents one cell measured; Topbp1+/+ number of cells = 31, number of mice = 2; Topbp1B5/B5 number of cells = 48, number of mice = 2; p-value=0.6374; mid-pachytene spreads: each dot represents one cell measured; Topbp1+/+ number of cells = 35, number of mice = 2; Topbp1B5/B5 number of cells = 86, number of mice = 2; p-value=0.3990; late-pachytene: each dot represents one cell measured; Topbp1+/+ number of cells = 49, number of mice = 2; Topbp1B5/B5 number of cells = 14, number of mice = 2; p-value=0.3474.) p-Values were calculated using a linear mixed effect model (see ‘Materials and methods’ for details).

Figure 6—figure supplement 3
Topbp1B5/B5 pachytene spermatocytes display no difference in H3K27ac, H3K36me3, or K3K4me3 staining compared to Topbp1+/+.

(A) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and H3K27ac on the left, and quantification of H3K27ac-stained pachytene spreads on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 143, number of mice = 2; Topbp1B5/B5 number of cells = 56, number of mice = 2; p-value=0.7881). (B) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and H3K36me3 on the left, and quantification of H3K36me3-stained pachytene spreads on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 79, number of mice = 2; Topbp1B5/B5 number of cells = 52, number of mice = 2; p-value=0.2778). (C) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and H3K4me3 on the left, and quantification of H3K4me3-stained pachytene spreads on the right. (Early-pachytene spreads: each dot represents one cell measured Topbp1+/+ number of cells = 7, number of mice = 2; Topbp1B5/B5 number of cells = 46, number of mice = 2; p-value=0.9895; mid-pachytene spreads: each dot represents one cell measured; Topbp1+/+ number of cells = 26, number of mice = 2; Topbp1B5/B5 number of cells = 58, number of mice = 2; p-value = 0.3292.). p-Values were calculated using a linear mixed effect model (see ‘Materials and methods’ for details).

Figure 6—figure supplement 4
Topbp1B5/B5 pachytene spermatocytes display no difference in H4K16ac, H4ac, or H2AK166ub staining compared to Topbp1+/+.

(A) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and H4K16ac on the left, and quantification of H4K16ac-stained pachytene spreads, on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 19, number of mice = 1; Topbp1B5/B5 number of cells = 53, number of mice = 2; p-value=0.7314). (B) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and H4ac on the left, and quantification of H4ac-stained pachytene spreads on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 26, number of mice = 2; Topbp1B5/B5 number of cells = 87, number of mice = 2; p-value=0.6795). (C) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and H2A-K199-ub on the left, and quantification of H2A-K199-ub-stained pachytene spreads on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 35, number of mice = 2; Topbp1B5/B5 number of cells = 20, number of mice = 2; p-value=0.9941). p-Values were calculated using a linear mixed effect model (see ‘Materials and methods’ for details).

Figure 6—figure supplement 5
From early pachytene to beginning of late pachytene, Topbp1B5/B5 spermatocytes exhibit the same localization pattern of H3K4me1 compared to Topbp1+/+.

Representative images of Topbp1+/+ and Topbp1B5/B5 early, mid, and late pachytene spermatocytes stained with H3K4me1 and SYCP3.

Figure 6—figure supplement 6
Topbp1B5/B5 pachytene spermatocytes exhibit no differences in the intensity or localization patterns of RNApol2, pRNApol2_Ser2, or pRNApol2_ser5 compared to Topbp1+/+.

(A) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and RNApol2 on the left, and quantification of RNApol2-stained pachytene spreads on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 242, number of mice = 6; Topbp1B5/B5 number of cells = 289, number of mice = 6; p-value for X and Y cores = 0.2026, p-value for X and Y chromatin loops = 0.4407). (B) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and pRNApol2-Ser2 on the left, and quantification of pRNApol2-Ser2-stained pachytene spreads on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 66, number of mice = 5; Topbp1B5/B5 number of cells = 66, number of mice = 5; p-value=0.3862). (C) Representative images of meiotic spreads showing Topbp1+/+ and Topbp1B5/B5 pachytene spermatocytes stained with SYCP3 and pRNApol2-Ser5 on the left, and quantification of pRNApol2-Ser5-stained pachytene spreads on the right (each dot represents one pachytene cell measured; Topbp1+/+ number of cells = 47, number of mice = 3; Topbp1B5/B5 number of cells = 21, number of mice = 2; p-value=0.5104). p-Values were calculated in GraphPad using a linear mixed effect model. p-Values were calculated using a linear mixed effect model (see ‘Materials and methods’ for details).

Figure 6—figure supplement 7
Genes used as markers to define the stages of the germ cells in the scRNAseq analysis.

Uniform Manifold Approximation and Projection (UMAP) plots showing markers for spermatogonia (Sal1 and Dmrt1), spermatocytes (Dazl, Id4, Sycp3, and Shcbp1l), and spermatids (Acrv1, Oaz3, and Prm2).

Figure 6—figure supplement 8
Genes used as markers to track spermatogenesis progression and form the 47 sub-cluster for all cells captured in the 10X platform in the scRNAseq analysis.

Bubble plot displaying gene markers and their respective pattern of expression during spermatogenesis.

Figure 6—figure supplement 9
Genes used as markers of somatic cells captured in the scRNAseq analysis.

Uniform Manifold Approximation and Projection (UMAP) plots show markers for the somatic cells captured on the germ cell-enriched population, endothelial and hematopoietic (Col1a2, Vcam1, Laptm5, and Insl3), smooth muscle (Acta2), and Sertoli (Ptgds and Wt1).

Figure 6—figure supplement 10
Analysis of sub-clusters of all cells captured in the 10X platform.

Uniform Manifold Approximation and Projection (UMAP) plot displaying the 47 sub-clusters of cells captured in the scRNAseq analysis.

Figure 6—figure supplement 11
Analysis of the main 24 sub-clusters captured in the germ cells-enriched scRNAseq analysis.

Uniform Manifold Approximation and Projection (UMAP) plot showing 24 sub-clusters for spermatogonia, spermatocytes, round spermatids, spermatids, Sertoli, smooth muscle, and endothelial and hematopoietic cells captured in the 10X platform.

Figure 6—figure supplement 12
Spermatocytes from Topbp1+/+ and Topbp1B5/B5 show correlation greater than 0.9.

Plot of Pearson correlation values between cell groups separated by cluster and genotype using log-normalized average expression of all genes.

Figure 6—figure supplement 13
Topbp1B5/B5 pachytene spermatocytes show increased expression of X-linked genes.

Heatmap displaying the expression levels of X-chromosome genes in spermatocyte 3 stage (pachytene) of Topbp1B5/B5 on the right and Topbp1+/+ on the left. Each line represents a gene, and each row is a cell.

Figure 6—figure supplement 14
Topbp1B5/B5 pachytene spermatocytes show increased expression of XY ‘killer genes’ and other X and Y genes typically used to illustrate meiotic sex chromosome inactivation (MSCI) defects.

Violin plots showing the level of expression and number of cells expressing the XY chromosome genes in SP3 (pachytene): Zfy1, Zfy2, Kdm6a, Lamp2, Scml2, Zfx, Uba1y, and Rhox13 in Topbp1B5/B5 and in Topbp1+/+. The Wilcoxon test was used to calculate the p-values. **p<0.01, ***p<0.001, and **** p<0.0001.

Figure 6—figure supplement 15
Topbp1B5/B5 pachytene spermatocytes show increased expression of XY-linked genes.

Bar plots showing the number of cells expressing the XY chromosome genes in SP3 (pachytene): Zfy1, Zfy2, Kdm6a, Lamp2, Scml2, Zfx, Uba1y, and Rhox13 in Topbp1B5/B5 and in Topbp1+/+. Fisher’s exact test was used to calculate the p-values. **p<0.01, ***p<0.001, and ****p<0.0001.

Figure 6—figure supplement 16
Quality control of the scRNAseq data.

Violin plots show the quality control for the scRNAseq data, displaying the distribution of gene counts, UMI counts, and mitochondrial gene percentage per library before filtering data.

Figure 7 with 2 supplements
TOPBP1 regulates silencing dynamics of X genes at the spermatocyte 3 stage.

(A) Illustration of the gene markers used to define spermatocyte 1 as leptotene, spermatocyte 2 as zygotene, and spermatocyte 3 as pachytene; hypothetical examples illustrating the categorization of transitions in silencing dynamics between the stages of pre-leptotene (PL), spermatocyte 1 (SP1), spermatocyte 2 (SP2), and spermatocyte 3 (SP3). (B) Number of genes in each of the categories described in (A), during the different stage transitions and respective p-values above each graph (the p-values were calculated using the Fisher’s exact test). (C) Examples of genes with altered silencing dynamics in the Topbp1B5/B5, red = reduced silencing, blue = unaltered silencing and green = increased silencing (D) Scatter plot showing the difference in RNA level between Topbp1+/+ and Topbp1B5/B5 for each of the indicated stage transitions. (E) Scatter plot showing expression level of X-chromosome genes, normalized to pre-leptotene levels, in Topbp1+/+ (gray) and Topbp1B5/B5 (blue) at SP3. (F) Graph plotting expression levels of X-chromosome genes, normalized to pre-leptotene levels, in Topbp1+/+ (Y axis) and Topbp1B5/B5 (X axis) and split in three categories based on the severity of silencing defect. (G) Box plot showing PL expression levels of X-chromosome genes in each of the categories of silencing defect severity shown in (F).

Figure 7—figure supplement 1
Analysis of the gene marker used to define the zygotene stage.

Uniform Manifold Approximation and Projection (UMAP) analysis of the gene Gm960 was used as a gene marker to define the zygotene stage.

Figure 7—figure supplement 2
Expression levels of genes with no defects, mild defects, and strong defects in Topbp1B5/B5 and Topbp1+/+.

Boxplot showing the expression levels of the clusters of genes with no defects, mild defects, and strong defects from pre-leptotene to Spermatocyte 3 (pachytene) in Topbp1B5/B5 and Topbp1+/+. The Wilcoxon rank-sum test was used to calculate the p-values. *p<0.05, **p<0.01, ***p<0.001, and ****p<0.0001.

A new TOPBP1 mutant separates XY silencing from sex body formation.

Schematic of sub-stages of meiotic prophase I. In wild-type mice, meiotic sex chromosome inactivation (MSCI) initiates following the accumulation of the DNA damage response (DDR) proteins at the XY chromosomes. During mid-pachytene, the XY body is fully formed, and transcription is restricted to the autosomes. In Atr or Topbp1 CKOs, the sex body is not formed, and the DDR proteins are not sequestered to the XY. Asynapsis events and transcription of toxic genes at the sex chromosomes are observed, triggering mid-pachytene arrest. In Topbp1B5/B5, MSCI initiates, the sex body is normally formed with normal recruitment of DDR proteins to the X and Y chromosomes, yet cells fail in the reinforcement/maintenance of silencing. Cells progress through mid-pachytene but not into diplonema.

Author response image 1
Topbp1B5/B5 displays normal condensation of the XY-pair.

(A) Immunostaining of XY condensation in Topbp1+/+ and Topbp1B5/B5 testes sections (γH2AX: green and DAPI: gray). (B) Quantification of all sex body-positive cells per tubule (Topbp1+/+ number of cells counted = 781, number of tubules counted = 28, number of mice = 3; Topbp1B5/B5 number of cells counted = 967, number of tubules counted = 28, number of mice = 3). (C) Quantification of elongated-sex body cells per tubule (Topbp1+/+ number of cells counted = 19 and 762 normal round/oval-sex bodies cells, number of tubules counted = 28, number of mice = 3; Topbp1B5/B5 number of cells counted = 45 and 922 normal round/oval-sex bodies cells, number of tubules counted = 28, number of mice = 3).

Author response image 2
Genes used as markers of spermatocytes captured in the scRNAseq analysis.

Violin plots display the distribution of cells expressing Gm960 (Leptotene marker), Meiob(Leptotene/Zygotene marker), Psma8 (Pachytene marker), Pwill1 (Pachytene marker), Pou5f2 (Diplotene marker), and Ccna1 (Secondary Spermatocytes marker).

Author response image 3
Single cell RNAseq reveals that Topbp1B5/B5 spermatocytes initiate MSCI but fail to promote full silencing of Y chromosome-linked genes.

Violin plot displaying the ratio of the average expression of Y chromosome genes by the average expression of chromosome 9 genes at different stages of spermatogenesis for Topbp1+/+ and Topbp1B5/B5 cells.

Tables

Appendix 1—key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Sequence-based reagentp3xflag-TOPBP1-Flag-N-terminal-FIDTPCR primerattcatcgatagatctgataAT
GTCCAGAAATGACAAAGA
Sequence-based reagentp3xflag-TOPBP1-Flag-N-terminal-RIDTPCR primertagagtcgactggtaccgatttag
TGTACTCTAGGTCGTT
Sequence-based reagentTOPBP1_K1317A_RIDTSite-directed mutagenesistgccactgaggctaaatacg
cctcgtttcgaagtggatgt
Sequence-based reagentTOPBP1_K1317A_FIDTSite-directed mutagenesisacatccacttcgaaacgagg
cgtatttagcctcagtggca
Sequence-based reagentTOPBP1_K155A_FIDTSite-directed mutagenesiscaggtttgcagcaactaaat
atgctttgctaccaacttctcctgca
Sequence-based reagentTOPBP1_K155A_RIDTSite-directed mutagenesistgcaggagaagttggtagcaaa
gcatatttagttgctgcaaacctg
Sequence-based reagentTOPBP1_K250A_RIDTSite-directed mutagenesistctcttggcacactcatacg
cctgaccttttggttcttgc
Sequence-based reagentTOPBP1_K250A_FIDTSite-directed mutagenesisgcaagaaccaaaaggtcagg
cgtatgagtgtgccaagaga
Sequence-based reagentTopbp1_BRCT5-Fwd genotypingIDTPCR primertgcatttccattaaccaacctc
Sequence-based reagentTopbp1_BRCT5 Rev genotypingIDTPCR primerggtagagttcaaatgtgtgtcatg
AntibodyAnti-phospho H2A.X (Ser139)
(mouse monoclonal)
MilliporeCat# 05-636; RRID:AB_309864IF (meiotic spreads) 1:50,000
IF 1:2000
AntibodyAnti-SCP3 antibody
(mouse monoclonal)
AbcamCat# ab97672; RRID:AB_10678841IF (meiotic spreads) 1:1000
AntibodyAnti-SYCP3
(rabbit polyclonal)
Lenzi et al., 2005CustomIF (meiotic spreads) 1:10,000
AntibodyAnti-SCP1
(rabbit polyclonal)
Abcamab15090IF (meiotic spreads) 1:1000
AntibodyAnti-Rad51
(rabbit polyclonal)
MilliporePC130IF (meiotic spreads) 1:1000
AntibodyAnti-ATR
(rabbit polyclonal)
Cell signaling2790IF (meiotic spreads) 1:1000
AntibodyAnti-TOPBP1
(rabbit polyclonal)
Danielsen et al., 2009CustomIF (meiotic spreads)
1:500
Western blot 1:1000
AntibodyAnti-phospho-Chk1 (ser317)
(rabbit monoclonal)
Cell Signaling12302IF (meiotic spreads)
1:100
Western blot 1:1000
AntibodyAnti-H1T
(guinea pig polyclonal)
A gift from Dr. Mary Ann Handel; Inselman et al., 2003CustomIF (meiotic spreads)
1:500
AntibodyAnti-HORMAD2
(rabbit polyclonal)
Wojtasz et al., 2012CustomIF (meiotic spreads)
1:500
AntibodyAnti-HORMAD1
(rabbit polyclonal)
Wojtasz et al., 2012CustomIF (meiotic spreads)
1:500
AntibodyPhospho HORMAD2 (S271)
(rabbit polyclonal)
Wojtasz et al., 2009CustomIF (meiotic spreads)
1:500
AntibodyAnti-GAPDH
(mouse monoclonal)
Thermo Fisher ScientificAM4300Western blot 1:5000
AntibodyAnti-β-actin
(rabbit polyclonal)
Cell Signaling4967Western blot 1:5000
AntibodyAnti-phospho MDC1 (T4)
(rabbit polyclonal)
AbcamAb35967IF (meiotic spreads)
1:500
AntibodyAnti-SETX
(rabbit polyclonal)
AbcamAb220827IF (meiotic spreads)
1:100
AntibodyAnti-MLH1
(mouse monoclonal)
BD Biosciences550838IF (meiotic spreads)
1:200
AntibodyAnti-MLH3
(guinea pig polyclonal)
Holloway et al., 2014CustomIF (meiotic spreads)
1:200
AntibodyAnti-CHK1
(mouse monoclonal)
Santa Cruzsc-8408Western blot 1:1000
AntibodyAnti-CHK2
(mouse monoclonal)
Millipore05-649 (clone7)Western blot 1:1000
AntibodyAnti-phospho CHK1 (S345)
(rabbit polyclonal)
 Cell Signaling2341IF (meiotic spreads)
1:200
Western blot 1:1000
AntibodyAnti-phospho RPA (S4/S8)
(rabbit polyclonal)
BethylA300-245AWestern blot 1:1000
AntibodyRPA (made against full length RPA2 expressed and purified in E. coli, and injected to rabbit)
(rabbit polyclonal)
N/ACustomWestern blot 1:1000
AntibodyAnti-KAP1
(rabbit monoclonal)
Bethyla700-014-TWestern blot 1:1000
AntibodyAnti-phospho KAP1 (S824)
(rabbit polyclonal)
BethylA300-767A-TWestern blot 1:1000
AntibodyAnti-53BP1
(rabbit polyclonal)
Cell Signaling4937Western blot 1:1000
AntibodyAnti-53BP1
(rabbit polyclonal)
Novus BiologicalsNB100-304IF (meiotic spreads)
1:200
AntibodyAnti-BLM
(rabbit polyclonal)
Abcamab2179Western blot 1:500
AntibodyAnti-Rad9
(rabbit polyclonal)
BethylA300-890A-TWestern blot 1:1000
AntibodyAnti-BRCA1
(rabbit polyclonal)
Kakarougkas et al., 2013CustomIF (meiotic spreads)
1:200
AntibodyAnti-BRCA1
(rabbit polyclonal)
Proteintech22362-1-APWestern blot 1:1000
AntibodyAnti-MDC1
(rabbit polyclonal)
Modzelewski et al., 2015CustomIF (meiotic spreads)
1:200
IF
1:200
AntibodyAnti-RNA pol 2
(mouse monoclonal)
Millipore05-623IF (meiotic spreads)
1:2000
AntibodyAnti-phospho RNA pol 2 (Ser2)
(rat monoclonal)
Millipore04-1571IF (meiotic spreads)
1:400
AntibodyAnti-phospho RNA pol 2 (Ser5)
(rat monoclonal)
Millipore04-1572IF (meiotic spreads)
1:400
AntibodyAnti-H3K9ac
(rabbit polyclonal)
AbclonalA7255IF (meiotic spreads)
1:200
AntibodyAnti-H3K9me3
(rabbit polyclonal)
Active Motif39766IF (meiotic spreads)
1:200
AntibodyAnti-H3K4me1
(rabbit polyclonal)
AbclonalA2355IF (meiotic spreads)
1:200
AntibodyAnti-CHD4
(rabbit polyclonal)
AbclonalA10557IF (meiotic spreads)
1:200
AntibodyAnti-Sumo_2/3
(rabbit polyclonal)
Proteintech11251-1-APIF (meiotic spreads)
1:200
AntibodyAnti-USP7
(mouse monoclonal)
Proteintech66514-1-IgIF (meiotic spreads)
1:200
AntibodyAnti-H3K27ac
(rabbit polyclonal)
Active Motif39134IF (meiotic spreads)
1:200
AntibodyAnti-H3K36ac
(rabbit polyclonal)
Active Motif61102IF (meiotic spreads)
1:200
AntibodyAnti-H3K4me3
(rabbit polyclonal)
Active Motif39160IF (meiotic spreads)
1:200
AntibodyAnti-H4K16ac
(rabbit polyclonal)
AbclonalA5280IF (meiotic spreads)
1:200
AntibodyAnti-H4c
(rabbit polyclonal)
Millipore06-598IF (meiotic spreads)
1:200
AntibodyAnti-H2AK199ub
(rabbit monoclonal)
AbcamAb193203IF (meiotic spreads)
1:200
AntibodyGoat anti-rabbit IgG (H+L) highly cross-adsorbed secondary antibody, Alexa Fluor 488
(goat polyclonal)
Thermo Fisher ScientificA-11034IF (meiotic spreads)
1:1000
AntibodyGoat anti-mouse IgG (H+L) antibody, Alexa Fluor 488 conjugated
(goat polyclonal)
Thermo Fisher ScientificA-11017IF (meiotic spreads)
1:1000
AntibodyGoat anti-rabbit IgG (H+L) antibody, Alexa Fluor 594 conjugated
(goat polyclonal)
Thermo Fisher ScientificA-11012IF (meiotic spreads)
1:1000
AntibodyGoat anti-mouse IgG (H+L) highly cross-adsorbed secondary antibody, Alexa Fluor Plus 594
(goat polyclonal)
Thermo Fisher ScientificA32742IF (meiotic spreads)
1:1000
AntibodyGoat anti-guinea pig IgG (H+L) highly cross-adsorbed secondary antibody, Alexa Fluor 647
(goat polyclonal)
Thermo Fisher ScientificA-21450IF (meiotic spreads)
1:1000
Commercial assay or kitApopTag Plus Peroxidase In Situ Apoptosis KitMilliporeS7101
Software, algorithmGraphPad Prism 9GraphPad PrismRRID:SCR_002798
Software, algorithmCellranger v7.0.010X Genomicshttps://support.10xgenomics.com/single-cell-gene-expression/software/downloads/3.1
Software, algorithmR v4.2.1R FoundationRRID:SCR_001905
Software, algorithmSeurat v4.3.0Hao et al., 2021RRID:SCR_007322
Software, algorithmsva v3.44.0Leek et al., 2023RRID:SCR_012836
Software, algorithmSingleR v1.10.0Aran et al., 2019RRID:SCR_023120
Software, algorithmReshape2 v1.4.4Wickham, 2007RRID:SCR_022679
Software, algorithmbiomaRt v2.52.0Durinck et al., 2009RRID:SCR_019214
Software, algorithmggpubr v0.5.0Kassambara, 2019RRID:SCR_021139
Software, algorithmTidyverse v1.3.2Wickham et al., 2019RRID:SCR_019186
Software, algorithmggnewscale v0.4.8Campitellihttps://eliocamp.github.io/ggnewscale/index.html
Software, algorithmHarmony v0.1.1Korsunsky et al., 2019RRID:SCR_022206
Software, algorithmOrg.Mm.eg.db v3.15.0Carlson et al., 2019RRID:SCR_023488
Software, algorithmCorrplot v0.92Wei and Simko, 2017RRID:SCR_023081
Software, algorithmSeuratWrappers v0.3.1Satija LabRRID:SCR_022555
Software, algorithmgghalves v0.1.4Tiedemann, 2022https://github.com/erocoar/gghalves
Author response table 1
MouseTopbplTopbpl
12643
24672
34384
47964
53729
64852
77467
8528

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Carolline Ascenção
  2. Jennie R Sims
  3. Alexis Dziubek
  4. William Comstock
  5. Elizabeth A Fogarty
  6. Jumana Badar
  7. Raimundo Freire
  8. Andrew Grimson
  9. Robert S Weiss
  10. Paula E Cohen
  11. Marcus B Smolka
(2024)
A TOPBP1 allele causing male infertility uncouples XY silencing dynamics from sex body formation
eLife 12:RP90887.
https://doi.org/10.7554/eLife.90887.3