Lipid droplets are present in early-stage somatic and germ cells.

(A) Testis lipid droplets (LD) in w1118 animals visualized with neutral lipid dye BODIPY. (A,A’) Scale bar=50 μm; (A’’,A’’’) scale bar=15 μm. Asterisk indicates hub in all images. Arrows point to LD; arrowheads point to spermatocytes in A,B. (B) Testis LD visualized with BODIPY in newly-eclosed males from two wild-type genotypes. Scale bars: main image=50 μm; inset image=10 μm. (C) Testis LD from w1118 animals at different times post-eclosion. Scale bars=50 μm. (D) Testis LD visualized with LipidTox Red in animals with somatic cell overexpression of GFP-LD (Tj-GAL4>UAS-GFP-LD). GFP-and LipidTox Red-positive punctae are somatic LD (D–D’’ arrows); LipidTox punctae without GFP indicate germline LD (D–D’’ arrowheads). Scale bars=10 μm. (E) Histogram showing the spatial distribution of somatic cell LD; error bars represent standard error of the mean (SEM). (F) Cumulative frequency distributions of somatic LD (blue line, data reproduced from E), zfh-1-positive somatic cells (zfh-1+ cells, orange line), and Eya-positive somatic cells (Eya+ cells, grey line). (G) Testis LD visualized with LipidTox Red in males with germline overexpression of GFP-LD (nos-GAL4>UAS-GFP-LD). GFP-and LipidTox Red-positive punctae indicate germline LD (arrows); LipidTox punctae without GFP indicate non-germline LD (arrowheads). Scale bars=10 μm. (H) Histogram representing the spatial distribution of LD within the germline; error bars represent SEM. (I) Histogram representing the spatial distribution of LD and GFP fluorescence (green line) (arbitrary units, a.u.) in a representative testis of a bam-GFP animal (panel J). (J) Testis LD in a bam-GFP animal; arrows point to LD and arrowheads point to spermatocytes. Scale bar=50 μm. See also Supplemental Figure 1.

bmm regulates testis lipid droplets in a cell-autonomous manner.

(A) Testis lipid droplets (LD) indicated by LipidTox Red in bmm-GFP animals. Arrows point to LD in all images. Arrowheads point to spermatocytes. Scale bars=50 μm. Asterisks indicate the hub in all images. (B) Quantification of nuclear GFP intensity in testes isolated from bmm-GFP animals (n=3). Germline stem cell (GSC), spermatogonia (SG), spermatocyte (SC). (C) Spatial distribution of LD (grey histogram) and GFP expression (green line) in testes from bmm-GFP animals as a function of distance from the hub (n=3). (D,E) LD near the apical region of the testis in bmmrev (D) or bmm1 (E) animals. (F) LD further away from the apical tip in bmm1 animals. (D–F) Scale bars=50 μm. (G) Histogram representing testis LD size distribution in bmmrev (grey) and bmm1 (orange). (H) Apical tip of the testes is at the left of the graph; individual dots represent a single LD and its relative position to the hub marked by an asterisk. Cumulative frequency distribution of the distance between LDs and the apical tip of the testes are drawn as solid lines. (I) Number of testis LD in bmmrev (grey) or bmm1 (orange) in FLP-FRT clones 3 days post-clone induction; dots represent measurements from a single clone. The number of cells in each cyst (CC) counted is indicated. See also Supplemental Figure 2.

A cell-autonomous role for bmm in regulating spermatogenesis.

Testes isolated from bmmrev (A) and bmm1 (A’) animals raised at 29°C stained with phalloidin. Scale bars=100 μm. (B) Testis size in bmm1 and bmmrev animals raised at 29°C. (C) Spermatid bundle number in bmm1 and bmmrev testes from animals reared at 29°C. (D,E) Representative images of bmmrev (D) or bmm1 (E) testes stained with DAPI and anti-Vasa antibody. Arrows indicate germline stem cells (GSC). Scale bar=50 μm. The hub is marked by an asterisk in all images. (F) GSC number in bmm1 and bmmrev testes. (G) Proportion of GSCs that were either bmm1 or bmmrev clones at 3 and 14 days post-clone induction. (H) Representative images of bmmrev (H) and bmm1 (H’) testes carrying bam-GFP; data quantified in Figure S3J. Arrows indicate regions with high Bam-GFP. Scale bars=50 μm. (I) Representative images of bmmrev (I) or bmm1 (I’,I’’) testes stained with anti-Vasa antibody. Arrows indicate Vasa-positive cysts in bmm1 testis. Panel I’’ is magnified from the boxed region in I’. (I,I’) Scale bars=100 μm; (I’’) scale bar=50 μm. (J) Maximum projection of bmmrev (J) or bmm1 (J’) testes stained with anti-Boule antibody (green) and DAPI (blue). Scale bars=100 μm. Number of bmm1 and bmmrev spermatocyte clones (K) or post-meiotic clones (L) at 3 and 14 days post-clone induction. See also Supplemental Figure 3.

Loss of bmm disrupts triglyceride homeostasis and leads to spermatogenic defects.

(A) Hierarchical clustering of lipid species detected in bmmrev and bmm1 animals. (B) Histograms showing the proportion of species in each lipid class with different levels between bmm1 and bmmrev. Numbers on histograms indicate the number of species with differences in abundance. (C) Volcano plot showing fold change in abundance of triglycerides (green; 97 species) and non-triglyceride lipids (grey; 186 species) in our dataset. (D) Arrows indicate testis LD stained with LipidTox Red in bmmrev (D), bmm1 (D’), or mdyQX25/k03902; bmm1 (D’’) animals. (E) Whole testes isolated from bmmrev (E), bmm1 (E’), or mdyQX25/k03902;bmm1 (E’’) animals stained with anti-Vasa antibody (red) and DAPI (blue). Arrowheads indicate spermatid bundles. Scale bars=100 μm. (F) Testis size in bmmrev, bmm1, and mdyQX25/k03902;bmm1 animals. Spermatid bundles (G) and number of germline stem cells (H) in bmmrev, bmm1, and mdyQX25/k03902;bmm1 animals. (I) Testis size in animals with germline-specific mdy knockdown (nos-GAL4>mdy RNAi; bmm1) compared with controls (nos-GAL4>+; bmm1 and +>mdy RNAi; bmm1). See also Supplemental Figure 4.

Cholesterol is absent from testis lipid droplets.

(A) Testes stained with BODIPY (A) to detect neutral lipids and Filipin III (A’) to detect free cholesterol. Scale bars=50 μm.

bmm mRNA levels during spermatogenesis in germline and somatic lineages.

(A) Pseudotime trajectory of germline (black line) based on single-cell RNA sequencing data62. Individual cells are labeled according to the annotation within the data set. (B) Rolling average of normalized bmm transcript counts in the germline along the trajectory shown in panel A are plotted as a black line on the upper panel. Composition of cell types mapped on to the trajectory at each time point is shown at the lower half of panel B. (C) Pseudotime trajectory of the somatic cells (black line) based on publicly available single-cell RNA sequencing data62. Individual cells are labeled according to the annotation within the data set. (D) Rolling average of normalized bmm transcript counts in somatic cells plotted as a black line along the trajectory shown in C (upper panel). Composition of cell types mapped on to the trajectory at each time point (lower panel). (E-H) Representative images of bmmrev (E and F) and bmm1 (G and H) testes with somatic over-expression of GFP-LD (Tj-GAL4>UAS-GFP-LD). Panel F and H contain magnified images of the area indicated by the boxes in panel E and G, respectively. In bmmrev testes, LD were restricted to a region near the apical tip (E) of the testis in both somatic (F–F’’’ arrows) and germline cells (F–F’’’ arrowheads). In bmm1 testes, LD were present in both somatic (G–H arrows) and germline cells (G–H arrowheads), near the apical tip of the testis in a region corresponding to early-stage germ cells and in the region corresponding to spermatocytes. (E,G) Scale bars=50 μm; (F,H) scale bars=20 μm.

Additional characterization of testis development and spermatogenesis in animals lacking bmm.

(A) Testis size was smaller in bmm1 mutant animals compared with bmmrev controls at <24 hr post-eclosion when raised at 25°C (A; Welch two-sample t-test). (B) The number of spermatid bundles was significantly lower in bmm1 mutant animals compared with bmmrev controls at <24 hr post-eclosion when raised at 25°C (Kruskal-Wallis rank sum test). (C) Testis size was significantly smaller in bmm1 mutant males compared with bmmrev control males at 14-days post-eclosion (Welch two-sample t-test). (D) While the median number of spermatid bundles was not significantly different between bmm1 mutant males and bmmrev control males at 14 days post-eclosion (Welch two-sample t-test), 8/27 bmm1 testis had no spermatid bundles, a phenotype absent in age-matched bmmrev males (0/22) (p=0.0163, Pearson’s Chi-squared test), suggesting a subtle defect is present. (E) Food supplemented with 4% medium chain triglyceride (MCT), but not long chain triglyceride (LCT), significantly increased testis length in bmm1 animals but had no effect on this phenotype in bmmrev control animals (one-way ANOVA with Tukey multiple comparison test). (F) Food supplemented with 4% medium chain triglyceride significantly increased the number of spermatid bundles in bmm1 testes but had no effect on this phenotype in bmmrev control animals (one-way ANOVA with Tukey multiple comparison test). (G) Representative images of bmmrev (G–G’) or bmm1 (G’’– G’’’) testes stained for Fas3 (G and G’’) and Vas (G’ and G’’’). Scale bars=25 μm. (H) Quantification of hub area in bmmrev or bmm1 testes showed a significantly larger hub size in bmm1 testes (Welch two-sample t-test). (I) The number of germline stem cells (GSC) undergoing mitosis (phospho-histone H3+ GSC/total GSC) was not significantly different between bmm1 and bmmrev testes (Kruskal-Wallis rank sum test). (J) The distance between the hub and the first Bam-GFP positive cyst (Figure 3H) was significantly higher in bmm1 testes than in bmmrev testes (Welch two-sample t-test). (K) All bmmrev testes and most bmm1 testes contained spermatids when raised at 25°C; however, the most advanced stage of spermatogenesis observed in the majority of bmm1 testes isolated from animals reared at 29°C was the spermatocyte stage. (L) Testes isolated from bmm1 animals showed a significantly smaller Boule-positive area than control testes (Welch two-sample t-test). (M) Testes isolated from bmm1 animals contain fewer individualization complexes than bmmrev control testes (Kruskal-Wallis rank sum test). (N) Fewer waste bags were present in testes isolated from bmm1 animals compared with bmmrev control testes (Kruskal-Wallis rank sum test).

Lipidomic analysis of animals lacking bmm.

(A) Higher fold-changes of triglycerides in bmm1 animals were associated with less saturation in the acyl-groups (Kendall’s rank correlation test). (B) Higher fold-changes of triglycerides in bmm1 animals were associated with higher number of carbons in the acyl-groups (Kendall’s rank correlation test). Each dot represents a single triglyceride species for panel B and C. (C) Volcano plot of identified lipids; monoglycerides shown in blue and diglycerides shown in orange. Many monoglycerides and diglycerides show increase in fold-change in bmm1 males. (D) The number of carbon and the degree of saturation of monoglycerides (MAG) and diglycerides (DAG) with significant changes in abundance between bmm1 and bmmrev males. (E) Volcano plot of identified lipids; fatty acids shown in magenta and acyl-carnitine shown in green. Many fatty acids show an increase in fold-change while many acyl-carnitines show a decrease in fold-change in bmm1 males. (F) The number of carbon and the degree of saturation of fatty acids (FA) and acyl-carnitines (ACar) with significant changes in abundance between bmm1 and bmmrev males. (G) Volcano plot of identified lipids; membrane lipids shown in yellow. (H) The number of carbon and the degree of saturation of membrane lipids with significant changes in abundance between bmm1 and bmmrev males. For panel G and H, PC: phosphatidylcholine; PE: phosphatidylethanolamine; PI: phosphatidylinositol; LPC: lysophosphatidylcholine; LPE: lysophosphatidylethanolamine; SM: sphingomyelin; PG: phosphatidylglycerol. (I) Loss of mdy function rescued the elevated number of LD in bmm1 testes to control levels (one-way ANOVA with Tukey multiple comparison test). (J) Germline-specific loss of mdy in bmm1 animals did not reduce GSC numbers, but the variance in GSC number was significantly rescued (nos-GAL4>+; bmm1 vs nos-GAL4>mdy RNAi; bmm1: p=4.5 × 10-5; +>mdy RNAi; bmm1 vs nos-GAL4>mdy RNAi; bmm1: p=0.0082 by F-test).