EMBO/SourceData huggingface dataset:
Figure 4Transmission electron microscopy image of cardiac muscle from scaf1Δ1/Δ1 (n=3) and scaf1+/+ fish (n=3). (A) Representative images showing mitochondria.Transmission electron microscopy image of cardiac muscle from scaf1Δ1/Δ1 (n=3) and scaf1+/+ fish (n=3). (B) Mitochondria size (100-150 mitochondria per biological sample). (C) Cristae lumen width (average of three cristae per mitochondria, 20 mitochondria per biological sample). Different biological replicates are represented with different color tones.D, E, Mitochondrial DNA copy number per nuclear copy number in muscle in females (D) and males (E) (Δ1 n=6, Δ2 n=6 and same number for their respective controls).F, Survival curve of 4 days post-fertilization embryos treated with different concentrations of the indicated OXPHOS inhibitors (three experimental replicates per biological replicate and three biological replicates).G, H, Oxygen consumption of 48 hours post-fertilization embryos using the XFe24 Seahorse analyzer, (G) Representative oxygen consumption rate (OCR) profile along time and (H) maximum OCR (Δ1 n=11, Δ2 n=11, n=10, n=11 respectively for their controls).I, Maximum uncoupled (FCCP) OCR in isolated mitochondria from adult fish (male and females Δ1 n=4 and Δ2 n=4, and same number for their respective controls) with the indicated site I [pyruvate (Pyr), glutamate (Glu), malate (Mal)] or site II [succinate (Succ)] substrates.J, K, Respiratory control ratio (RCR)(State 3/State4) (J) and P/O ratio (K) in isolated mitochondria from adult fish (male and females Δ1 n=4, and same number for their respective controls) with the indicated substrates.
Figure 5A, Representative images of scaf1-/- and scaf1+/+ fish fed with the indicated diets.B, C, Changes in (B) length and (C) weight over time (Δ1 +/+ n=10, Δ1 -/- n=10, Δ2 +/+ n=10, Δ2 -/- n=12-13.D-F, Adipose tissue measurements on hematoxylin-eosin (H&E)-stained adult zebrafish sagittal sections. (D) Adipose tissue area per total section area (average of three sections/biological replicate) and (E) adipocyte size (average of 20-30 adipocytes of ventral adipose tissue per biological replicate) (Δ1 +/+ n=3, Δ2 n=8). scaf1Δ1 and scaf1Δ2 are represented with circles and squares, respectively.Adipose tissue measurements on hematoxylin-eosin (H&E)-stained adult zebrafish sagittal sections. (F) Representative images of ventral fat deposits (dotted lines).G, Number of eggs per clutch (standard diet Δ1 +/+ n=8, Δ1 -/- n=8, Δ2 +/+ n=7, Δ2 -/- n=12, double diet Δ1 +/+ n=7, Δ1 -/- n=8, Δ2 +/+ n=8, Δ2 -/- n=9).H, Representative images of H&E-stained ovaries. Black dotted lines outline the ovaries, blue dotted line indicates a mature follicle. I, Quantification of mature ovary follicles per ovary section (average of three sections/biological sample) (Δ1 n=2-3, Δ2 n=6).
Figure 6A, B, Representative images of scaf1-/- and scaf1+/+ females (A) and males (B) fed with the indicated diets.C, D, Length of females (C) and males (D) after the indicated diets (Δ1 +/+ n=5, Δ1 -/- n=5, Δ2 +/+ n=5, Δ2 -/- n=5).E, F, Weight of females (E) and males (F) after the indicated diets (Δ1 +/+ n=5, Δ1 -/- n=5, Δ2 +/+ n=5, Δ2 -/- n=5).G, number of eggs per clutch (standard diet: Δ1 +/+ n=4, Δ1 -/- n=3, Δ2 +/+ n=3, Δ2 -/- n=2, double diet: Δ1 +/+ n=5, Δ1 -/-, n=5 Δ2 +/+ n=1, Δ2 -/- n=2). H, Quantification of mature ovary follicles per ovary section in haematoxylin-eosin (H&E) histological sections (average of three sections/biological sample) (Δ1 n=3, Δ2 n=3). adipose tissue quantification in H&E sections (Δ1 n=3, Δ2 n=3) in females and males : adipose tissue area (I, J)adipose tissue quantification in H&E sections (Δ1 n=3, Δ2 n=3) in females and males adipocyte size (K, L).
Figure 7A, Immunoblot of the indicated proteins of BNGE from female scaf1+/+ and scaf1-/- whole zebrafish mitochondria for the indicated diet (representative of n=2).B, Maximum uncoupled (FCCP) oxygen consumption rate in whole zebrafish mitochondria (females Δ1 n=4 and Δ2 n=4, and same number for their respective controls) with glutamate (Glu), malate (Mal) and succinate (Succ).RNAseq data from scaf1+/+ and scaf1-/- skeletal muscle for the indicated diet (standard diet scaf1+/+ and scaf1-/- n=4, double diet scaf1+/+ n=3, scaf1-/- n=4). (C-F) Volcano plots of differentially expressed genes (DEGs). (C) Comparison between scaf1-/- and scaf1+/+ zebrafish in standard diet. (D) Comparison between scaf1-/- and scaf1+/+ zebrafish in double diet. (E) Comparison of scaf1+/+ zebrafish in double diet and standard diet. (F) Comparison of scaf1-/- zebrafish in double diet and standard diet. In blue, significant DEGs (p-adj < 0.05, log2FC > |1|), in grey not significant DEGs, red circles represent non-significant differentially regulated OXPHOS genes, green circles represent significant differentially regulated OXPHOS genes, purple scaf1 (cox7a2l).
Figure 1HMRs express RNAs. Heat maps show RNA expression levels in 16.5 dpc prospermatogonia. Only uniquely mapped RNA-seq reads were used for the analyses.metaplots (B) show RNA expression levels in 16.5 dpc prospermatogonia. Only uniquely mapped RNA-seq reads were used for the analyses. Positions (1)-(4) in (A and B) (green arrows) indicate the same positions (±2 kbp from the centers). The lower chart in (B) represents mappability defined as the percentage of unique sequences in all possible 50-nt sequences (see Appendix Figure S2B). TSS isoforms are included in 55798 TSSs. Low mappability at the centers of HMR_-L1s is due to high occurrence of repeat (retrotransposon) sequences. RPM, reads per million mapped reads.
Figure 2Chromatin marks associated with promoters are observed in hypomethylated regions without young LINE1 retrotransposons (HMR_-L1s). Metaplots show the average ChIP-seq signals across HMR_-L1s, transcription start sites (TSSs), and randomly selected sites. The signals were normalized by either H3 or input signal. The results shown in heat maps were not normalized by H3 or input.The initiation form of Pol2 with an unphosphorylated carboxy terminal domain (Pol2 null-P CTD) is highly enriched in HMR_-L1s.
Figure 3RNA-immunoprecipitation (RIP) analyses of EGFP-MIWI2 and EGFP. EGFP-MIWI2 and EGFP were immunoprecipitated from 16.5 dpc prospermatogonia, and long RNAs bound to these proteins and input RNAs were sequenced. RNAs expressed from HMR_-L1s and piRNA clusters, but not Ensembl and prospermatogonium-expressed mRNAs, are overall enriched by immunoprecipitation (IP) of EGFP-MIWI2 (left). This enrichment is not observed when EGFP is immunoprecipitated (right). The numbers of mRNAs, piRNA clusters, and HMR_-L1s analyzed were shown in parentheses. Prospermatogonium-expressed mRNAs are those that are highly and specifically expressed in prospermatogonia. The bottom and top of boxes represent 25% percentile and 75% percentile, respectively. The line in the box represents the median. Circles represent outliers, and the bottom (top) of whiskers represents the minimum (maximum) values except for outliers.RIP analyses of EGFP-MIWI2 in newborn testes from WT and Mili KO mice. RNA bound to EGFP-MIWI2 and input RNA were sequenced. All RNAs derived from each class was counted for this analysis.HMR_-L1s preferentially generate piRNAs. X- and Y-axes represent expression levels of long RNAs (16.5 dpc prospermatogonia) as determined by RNA-seq and MIWI2-bound piRNAs (18.5 dpc testes) as determined by small RNA-seq, respectively. Circles in the charts represent individual mRNAs, piRNA clusters and HMR_-L1s. More than half of the HMR_-L1s (right) are located in the area of y > x (upper part of the diagonal), whereas most mRNAs (left) are located in the area of y < x (lower part of the diagonal). RPKM, reads per kilobase per million mapped reads.Frequent 10-nt overlap between sense and antisense piRNAs derived from HMR_-L1s. Top 100 HMR_-L1s that generate the largest number of unique piRNAs were analyzed. Frequencies of 1−18-nt overlaps were analyzed for each HMR_-L1.piRNAs derived from HMR_-L1s are enriched in adenine at their 10th nucleotide. HMR_-L1s generating more than 40 unique piRNAs were analyzed (678 HMR_-L1s). The proportions of piRNAs with adenine (left) and uridine (right) at the indicated positions were analyzed for each HMR_-L1. See A for the explanation of box plots.
Figure 4DNA methylation status of a portion of the Rasgrf1 differentially methylated region (DMR) in 0−2-dpp Mili+/− and Mili−/− prospermatogonia. Regions corresponding to B, C, D, E, and F are demethylated in Mili−/− prospermatogonia. Vertical bars represent DNA methylation levels at CpG sites determined by whole genome bisulfite sequencing (see below). The top dark gray bar represents a portion of the Rasgrf1 gametic DMR (Tomizawa et al, 2011). The CpG cluster is shown. The deleted RMER4B sequence is indicated.DNA methylation status of the Rasgrf1 DMR in Rasgrf1ΔRMER4B/ΔRMER4B and control spermatogonia are shown. The regions examined are shown in (A). Black (white) circles represent methylated (unmethylated) CpGs.Rasgrf1 RMER4B functions in cis. Wild-type and mutant alleles of the Rasgrf1 DMR in Rasgrf1+/ΔRMER4B spermatogonia were analyzed separately.
Figure 5Structure of the chr7 piRNA cluster. Top gray bars represent a putative piRNA precursor predicted by H3K4me3 ChIP-seq, RNA-seq, and small RNA-seq. Only uniquely mapped reads were used for displaying the sequencing result. For RNA-seq and small RNA-seq, only the results of minus strand are shown, because plus strand does not express appreciable level of RNA. A promoter region deleted is shown above the putative piRNA precursor. RPM, Read Per Million.Low expression level of piRNAs from chr7 cluster in chr7 piRNA clusterΔpromoter/Δpromoter testes. Unchanged expression of piRNAs from chr10 cluster (bottom) shows chr7 cluster is specifically affected. Only uniquely mapped piRNAs are shown.DNA methylation status of the Rasgrf1 differentially methylated region (DMR) in spermatogonia from chr7 piRNA clusterΔpromoter/Δpromoter. See Figures 4A and B for information about the region examined.
Figure 6Depletion of piRNAs from chr7 piRNA cluster RMER4B sequence in chr7 piRNA clusterΔRMER4B/ΔRMER4B mice. Top panel shows overall view of chr7 piRNA cluster, and lower panel shows magnified view of RMER4B region in chr7 piRNA cluster. Only uniquely mapped piRNAs are shown.DNA methylation status of the Rasgrf1 differentially methylated region (DMR) in spermatogonia from chr7 piRNA clusterΔRMER4B/ΔRMER4B mice.DNA methylation status of the Rasgrf1 differentially methylated region (DMR) in spermatocytes, spermatids, and sperms (D) from chr7 piRNA clusterΔRMER4B/ΔRMER4B mice.DNA methylation pattern in the progeny of male Rasgrf1ΔRMER4B/ΔRMER4B mice. DNA methylation status of the Rasgrf1 DMR in E12.5 whole embryos is shown. E12.5 whole embryos were obtained by the crosses shown in the panel. Two independent bisulfite-treated DNA sample were examined for Homo13 embryo.
Figure 1A Domain architecture of S. pombe bromodomain AAA-ATPases. The D1 and D2 ATPase domains (AAA) are shaded blue and orange, respectively, and the bromodomain (BD) is grey. The percentage similarity to each other and the indicated proteins was determined by FASTA sequence comparison [74] using the scoring matrix BLOSUM50.B Venn diagrams showing overlap between genes up-regulated (1.5 fold) in abo1Δ mutants with genes up-regulated under the indicated condition, along with the significance of the overlaps (based on hypergeometric distribution).
Figure 2A Normalised cumulative nucleosome (150 ± 30 bp size class) position frequency profiles for 4013 S. pombe genes aligned at the transcription start site (TSS) plotted from low MNase (biorep1) and high MNase (biorep 2) datasets. P-values (calculated using a two-tailed unpaired t-test) for the difference in means between random n = 100 subsets of the frequency values at the points indicated by the dotted lines ("-4" nucleosome and "+4" nucleosome) are shown.B Nucleosome position frequency values for the coding regions of 4013 S. pombe genes were k-means clustered (k = 9) using the abo1Δ data from biorep 1 (low MNase) and displayed with positive values coloured yellow and other values coloured blue (left-hand panel). The cluster order was then used to display the equivalent wild type frequency values in the right-hand panel.C Level of histone gene mRNAs was determined by qRT-PCR. Data are the mean of four independent biological repeats and error bars are ±SEM. Two-tailed unpaired t-tests showed no significant differences (P > 0.05) between wild type and abo1Δ cells.D Whole cell extracts were subjected to western blotting with histone H3 (Abcam), histone H2A (Abcam) and tubulin antibodies. Examples of the primary data are shown (left) along with a quantification of histone H2A and H3 levels normalized to tubulin (right). Data are the mean of at least three independent repeats and error bars represent ±SEM. P-values were calculated using a two-tailed unpaired t-test.
Figure 3A Whole cell extracts (WCE) were prepared from the indicated strains, immunoprecipitated (IP) with anti-GFP antibody and subjected to western blotting with anti-V5-Pk (Serotec) and anti-GFP (Life Technologies) antibodies. Data are representative of three independent biological repeats.B As for A.C RNA purified from wild type, abo1Δ and spt16-18 cells was analysed by northern blotting using a probe to the 3′ end of spbc19c7.11 (top panel). RNA (5 µg) used for northern blotting was analysed on an ethidium bromide stained 1% TAE agarose gel (bottom panel). Data are representative of two independent biological repeats.D RNA purified from wild type, and abo1Δ cells was analysed by strand specific RT-PCR. RNA from hip1Δ cells was analysed as a control. One primer, complementary to either the forward or reverse transcripts, was included during the reverse transcription step the second primer was then added during PCR amplification. Control reactions omitting the reverse transcription step (-RT) were included to demonstrate the absence of contaminating genomic DNA. Data are representative of two independent biological repeats.E The indicated strains were subjected to ChIP analysis with anti-GFP antibodies and the resulting DNA was analysed by qPCR for the indicated locus. Data are the mean of three independent biological repeats and error bars represent ±SEM. P-values were calculated using a two-tailed unpaired t-test. Significant (P < 0.05) values are indicated.F ChIP analysis of histone H3 levels at act1+ was determined by qPCR. Data are the mean of three independent biological repeats and error bars represent ±SEM. P-value was calculated using a two-tailed unpaired t-test.G The indicated strains were subjected to ChIP analysis as described for E.H Log phase cells were subjected to five-fold serial dilution and spotted onto rich (YE5S) agar or agar supplemented with MMS (0.006%). Plates were incubated for 3-5 days at 30°C or 7 days at 20°C. Images are representative of three independent biological repeats.
Figure 4A The position of ura4+ reporter alleles in centromere 1 is shown in the top panel. Strains containing the imr::ura4+ allele were grown to log phase in YE5S medium, subjected to 5-fold serial dilutions and spotted onto YE5S agar or YE5S agar supplemented with 5-FOA (1 mg/ml). Data are representative of three independent biological repeats.B Strains carrying the otr::ura4+ were analysed as described for A.C RNA was purified from the indicated strains and the level of centromeric transcripts determined by qRT-PCR. Data is the mean of three independent repeats and error bars indicate ± SEM. P-value was calculated using a two-tailed unpaired t-test.D RNA was purified from the indicated strains and subjected to qRT-PCR for the subtelomeric gene, tlh1+. Data is the mean of three biological repeats and error bars indicate ± SEM. P-value was calculated using a two-tailed unpaired Student's t-test.E The position of insertion of ade6+ reporter allele in the silent mat locus is shown in the top panel. Log phase cells containing the mat3-M::ade6+ allele were subjected to five-fold serial dilution and spotted onto YE5S agar plates lacking adenine (Low Ade) were incubated for 4 days at 30°C. Data are representative of at least three independent biological repeats. Note, microarray analysis indicates that deletion of abo1+ does not influence the expression of ade6+ (or ura4+) when these genes are present at their normal genomic loci.
Figure 5A ChIP analysis was performed on wild type (untagged) and abo1-GFP cells and the resulting DNA analysed by qPCR for centromeric (dh, dg and imr) repeat sequences. Data are the mean of four independent biological repeats and error bars represent ±SEM. P-values calculated using a two-tailed unpaired t-test indicates that all loci are significantly enriched (P < 0.05) relative the untagged control.B The indicated strains expressing GFP-Swi6 were subjected to ChIP analysis. The level of centromeric (dg) repeat sequences relative to the euchromatic control locus, adh1+ was determined by qPCR and scaled to a clr4Δ (- H3 K9me control) mutant. Data are the mean of two independent ChIP experiments and error bars represent the range of the data.C ChIP analysis of histone H3 levels at centromeric (dg) repeat sequences was determined by qPCR. Data are the mean of three independent ChIP experiments and error bars are ±SEM. P-values were calculated using a two-tailed unpaired t-test.D. ChIP analysis of Abo1-GFP at centromeric (dh and imr) repeat sequences regions in wild type and pob3Δ cells. Data are the mean of three independent experiments and error bars represent ±SEM. P-values, calculated using a two-tailed unpaired t-test, indicated no significant difference (P > 0.05) for wild type and pob3Δ cells.E ChIP analysis of Abo1-GFP at centromeric (dh and imr) repeat sequences in wild type and pob3Δ cells. Data are the mean of duplicate experiments and error bars represent the range of the data.F ChIP analysis of Abo1-GFP over the dh and imr regions in wild type and swi6Δ cells. Data are the mean of three independent experiments and error bars represent ±SEM. P-values, calculated using a two-tailed unpaired t-test, indicated no significant differences (P > 0.05) for wild type and swi6Δ cells.
Figure 6A RNA was extracted from mid log phase cells and Tf2 mRNA levels were determined by qRT-PCR, normalised to act1+ mRNA and scaled relative to the wild type level. Data are the mean of three independent biological repeats and error bars represent ±SEM.B Mid log phase cells with the indicated integrated lacZ reporter were subjected to quantitative β-galactosidase assays. Data are the mean of three independent biological repeats and error bars represent ±SEM.C As for BD ChIP DNA samples from wild type (untagged) and abo1-GFP cells were analysed by qPCR for Tf2 LTR. Data are the mean of four independent biological repeats and error bars represent ±SEM.E Normalised cumulative nucleosome (150 ± 30 bp) position frequency profiles for Tf2 LTR retrotransposons aligned at the ATG plotted from low MNase (biorep1) and high MNase (biorep 2) datasets.P-values were calculated using a two-tailed unpaired t-test.
Figure 7A Log phase cells were subjected to five-fold serial dilution and spotted onto rich (YE5S) agar or agar supplemented with TBZ (15 μg/ml) and incubated for 3-5 days 30°C. Data are representative of at three independent biological repeats.B The frequency of Ch16 minichromosome loss was assayed by plating cells onto adenine limiting agar and determining the frequency of half sectored (red-white) colonies. At least two independently derived strain isolates were used for each genotype and the total number of colonies counted for each genotype is indicated in parenthesis.C Heterothallic (h90) strains were cultured on nitrogen limiting (EMMG) medium for 3 days at 25°C to induce mating and meiosis. The resulting asci were stained with DAPI and visualised by fluorescence microscopy. Data are representative of three biological repeatsD Spore viability was measured by dissecting spores onto YE5S agar followed by incubation for 4-5 days at 30°C. Data are the mean of three independent biological repeats and error bars are ±SEM. P-value was calculated using a two-tailed unpaired t-test.
Fig 1H&amp;amp;E staining of adult head sections of (A) repo>GFP.nls and (B) repo>Shits1 at 29°C for 3 days. Lamina degeneration was identified as vacuoles in (B). Cryosection of adult (C) repo>H2B-RFP exhibited the expression of the nuclear red fluorescent protein (RFP) in the glia (epithelialglialnucleus: arrowhead) and (D) repo>Shits1 at 28°C for 12 days. Vacuoles were identified in the lamina neuropile (D). Epithelial (arrowhead) (E), marginal (arrow) and distal satellite (arrowhead) glianuclei (F) are labeled by HisCl-GAL4 and NP2109-GAL4, respectively. Note that HisCl-GAL4 is not expressed in all epithelialglia. (G, H) Weak lamina vacuolization was identified in HisCl>Shits1 (G, 0.92%) and NP2109>Shits1 (H, 1.46%) at 29°C for 14 days. (I, J) A single MARCMglia clone (GFP, green) expressed Shits1 at 21°C (I) and 29°C (J). (J) A vacuole occurred within a glial clone. DAPI (white) stains the nuclei (C, D, and G-J). (K) The percentage of the vacuole area in the lamina of the repo>Shits1flies at 28°C progressively increased. n indicated in each column. P-values were calculated using one-way ANOVA with Bonferroni's post-test. (L) When repo>Shits1flies were shifted to 28°C for 12 days and then shifted to 17°C for 9 additional days, the vacuolization was not alleviated. P-values were calculated using one-way ANOVA with Tukey's post-test. n indicated in each column.
Fig 2(A-E) Horizontal sections of the adult head lamina cartridge. (A) In repo>GFP.nls flies, two lamina neurons L1, L2 axons (L) and R1-6 axons (blue area) were surrounded by the electron-dense cytoplasm of the epithelial glia (G). The section was examined at three different depths, and the size of a single cartridge is not significantly different at different depths. (B-D) repo>Shits1 adults maintained at 29°C for 4 days. Axons were enlarged but intact. A large number of small and medium vacuoles were identified in the epithelial glia. (C) Small vacuoles in the glia that contained a double-membrane structure (arrowhead). (C') Higher magnification of the boxed area in (C). (D) A double-membrane autophagosome-like structure (arrowhead) in the glialcytoplasm. (E) In repots>DERDN adults maintained at 29°C for 2 days, large vacuoles were identified in between the lamina cartridges. N: glialnuclei. (E') Higher magnification of the boxed area in (E), which shows autophagosome-like vesicles (arrowhead) in the gliacytoplasm.
Fig 3(A) Rh1>lacZ+H2B-RFP exhibited nuclearRFP expression in the R1-6 photoreceptors (red). (B) Rh1>lacZ+Shits1 maintained at 29°C for 14 days resulted in lamina vacuolization. (C) Pan-Rh7>H2B-RFP exhibited expression in R7 (arrow) and R8 (arrowhead) (red). (D) Pan-Rh7>Shits1 incubated at 29°C for 14 days did not exhibit lamina degeneration. (E) Lamina monopolar (L) neurons (red) were selectively labeled by the Ln-GAL4 driven H2B-RFP with repo-GAL80. (F) Ln>Shits1, repoGAL80 shifted to 29°C for 14 days did not exhibit lamina degeneration. (G) R1-6 photoreceptors were killed in Rh1ts>lacZ+Hid shifted to 29°C for 14 days. Rh1ts indicates tubGAL80ts; Rh1-GAL4. Degeneration was induced in the lamina in addition to the retina. (H) rdgC306/+ mutant exposed to constant light for 14 days exhibited degeneration in the retina, as well as the lamina. DAPI: nuclei (white in A-H). (I, K) Rh1>lacZ+Shits1 at 29°C for 4 days exhibited vacuoles in the electron dense cytoplasm of epithelialglia or near the glianucleus. Arrow: vacuole with internal debris. Arrowhead: large vacuole. (J) Lamina cartridge of Rh1>lacZ+H2B-RFP. N indicates glialnuclei. (L) The percentage of the vacuole area in the lamina at 29°C for 14 days was examined. All P-values were calculated using one-way ANOVA with Tukey's post-test. Scale bar: 20 μm.