Loss of UBL7 leads to sperm malformation and male sterility
Previous studies have reported that the mRNA of Ubl7 is highly expressed in testis tissue33. Western blot (WB) in multiple mouse tissues showed that UBL7 protein is enriched in testis (Fig. 1a). UBL7 can be detected in testis at day postpartum 8 (dpp 8) and its level gradually increases till dpp 35 and adulthood (Fig. 1b). Immunofluorescence analysis of adult testis sections demonstrated that UBL7 is expressed in all germ cell types (Fig. 1c). Immunofluorescence of single germ cell showed that UBL7 localizes to both cytoplasm and nucleus of germ cells with highest density in cytoplasm of round spermatids (Supplementary Fig. 1). Additionally, UBL7 was enriched in the residual bodies that are attached to or shed off from elongating spermatids (Supplementary Fig. 1e, f).
To investigate the physiological function of UBL7 during spermatogenesis, we generated Ubl7 knockout (Ubl7−/−) mice with a deletion of exon 5 that causes a premature termination at exon 6 (Supplementary Fig. 2a). WB and immunofluorescence analysis confirmed the absence of UBL7 protein in Ubl7−/− testes (Supplementary Fig. 2b, c). Fertility tests showed that all Ubl7−/− male mice failed to produce offspring (Supplementary Fig. 2d), underscoring the essential role of UBL7 in spermatogenesis.
Further investigations demonstrated that the testes of Ubl7−/− mice were significantly smaller than those of Ubl7+/− mice (Fig. 1d, e). The stage distribution of seminiferous epithelium cycle in Ubl7−/− mice appeared disorganized (Supplementary Fig. 2g). Despite the disorganization, the seminiferous tubules at each stage are uniformly shorter in Ubl7−/− testes, resulting in an unchanged distribution of germ cell subpopulations (Fig. 1f, g, Supplementary Fig. 2e). Notably, we observed abnormal retention of elongated spermatids in stage IX-X seminiferous tubules of Ubl7−/− mice (Fig. 1f), indicating defective spermiation. Sperm isolated from cauda epididymis of Ubl7−/− mice exhibited a reduced count and abnormal morphology (Fig. 1h-j, Supplementary Fig. 2f). The sperm tails showed abnormal bending at the head-neck connecting region and the midpiece-principal piece connecting region (Fig. 1i). Additionally, the head shape was smaller in size and had a reduced diameter, particularly around the equator region (Fig. 1j). Collectively, our findings highlight that UBL7 is enriched in germ cells and loss of UBL7 leads to malformation in both sperm tail and head, leading to male sterility in mice.
Ubl7 −/− spermatids exhibit abnormal elongated manchette structures
Given that Ubl7−/− sperm exhibited a smaller head shape, we sought to investigate the underlying cause of this abnormality. The shaping of the mammalian sperm head involves several key processes, including histone replacement, acrosome development, and the formation of the manchette structure35. Our immunofluorescent staining of Histone H3, transition protein 1 (TNP1) and protamine 1 (PRM1) revealed that the expression patterns of these proteins were not altered in Ubl7−/− spermatids (Supplementary Fig. 3a-c), suggesting that histone replacement was not affected by UBL7 deficiency. Additionally, PNA staining indicated no abnormalities in acrosome morphology during spermatid development in Ubl7−/− testes (Supplementary Fig. 3a-c).
During spermiogenesis, a transient skirt-like structure named manchette progressively moves down the spermatid head, concomitantly constricts to sculpt the base of the nucleus4. We performed immunofluorescent staining using an anti-alpha tubulin (α-tubulin) antibody to label manchette structures in testis sections. A stage-by-stage comparison of sections from Ubl7+/− and Ubl7−/− testes revealed that manchettes in Ubl7−/− spermatids initiated formation at the correct time, but their removal was delayed (Fig. 2a). Morphological analysis of manchettes in isolated spermatids revealed that the formation of manchettes in Ubl7−/− spermatids proceeded normally at step 9, but their movement down the sperm head stalled around step 12–13 and the disassembly was also delayed around step 14 (Fig. 2b). Intriguingly, we occasionally observed a subset of spermatids with narrow, “rod-like” heads and conspicuously elongated manchettes in Ubl7−/− mice, which were rarely seen in Ubl7+/− mice (Fig. 2b). Furthermore, we isolated stage II seminiferous tubules under a transilluminating dissection microscope and extracted spermatids around step 14. We found that manchettes in Ubl7−/− spermatids around step 14 had increased length and wider coverage of the sperm nucleus (Fig. 2c, d). Malformation of manchettes leads to excessive constriction of the sperm head, ultimately resulting in the abnormal head shape of Ubl7−/− sperm.
Deficiency in ultrastructure of HTCA and annulus in Ubl7−/− sperm
Sperm Isolated from Ubl7−/− epididymis exhibited severe abnormalities in tail shape. We classified these abnormal tail shapes into four clusters: 1. Headless (Hless), sperm with separated head and tail; 2. Head-neck bending (HNB), sperm with a bending tail at the head-neck connecting region; 3. Midpiece-principal-piece bending (MPB), sperm with a bending tail at midpiece-principal-piece connecting region; 4. Midpiece-principal-piece gap (MPG), sperm with separation between midpiece and principal piece at connecting region (Fig. 3a). Quantification showed that the ratios of these four types of abnormal sperm were all significantly higher in Ubl7−/− mice (Fig. 3b). Interestingly, we also found that the abnormalities were exacerbated during sperm transiting through epididymis (Supplementary Fig. 4a, b).
To gain further insights into the morphological abnormalities of Ubl7−/− sperm tails, we compared the ultrastructure of cauda epididymal spermatozoa from Ubl7+/− and Ubl7−/− mice by transmission electron microscopy (TEM). First of all, Ubl7−/− sperm exhibited an integrated and regular “9 + 2” axonemal arrangement of the flagella (Supplementary Fig. 4c). The head-neck connecting region is associated with the HTCA. TEM analysis of longitudinal sperm sections revealed that in Ubl7−/− sperm, segmented columns (Sc) appeared blurry and disorganized, and capitulum (Cp) was incomplete and even completely separated from basal plate (Bp) (Fig. 3c). Notably, in the headless sperm from Ubl7−/− mice, the HTCA was entirely detached from sperm head (Fig. 3c). In summary, UBL7 deficiency leads to fragile HTCA and predisposes sperm to bend or even fracture at the head-neck connecting region, thus resulting in acephalic spermatozoa.
The midpiece-principal-piece connecting region is associated with the annulus structure. Disruption of the annulus typically results in a sharp bend at the midpiece-principal-piece connecting region36–38. Septin-4 (SEPT4) is one of the main components of annulus. We performed immunofluorescent staining of SEPT4 in Ubl7+/− and Ubl7−/− cauda sperm to label annulus. Results showed that SEPT4 localized to the boundary between the midpiece and principal piece, appearing as short cross-lines or two bright punctate signals (Supplementary Fig. 4d). Virtually, SEPT4 signal appeared normal in Ubl7−/− sperm, regardless of whether the tail exhibited bending abnormalities or not (Supplementary Fig. 4d). Further examination using TEM demonstrated that an electron-dense, wedge-shaped annulus was present in both Ubl7+/− and Ubl7−/− cauda sperm (Fig. 3d). However, in Ubl7−/− sperm with a gap between the midpiece and principal piece, we observed a membrane sink between the mitochondrial sheath and the fibrous sheath. The annulus shifted distally along with the fibrous sheath and lost its interaction with the mitochondrial sheath (Fig. 3d). In conclusion, the annulus exists but its localization is altered in Ubl7−/− sperm. Based on these results, UBL7 may not participate directly in the assembly of the annulus, but may play a role in ensuring the proper attachment of the annulus to both the mitochondrial sheath and the fibrous sheath. However, the factors influencing these connections remain obscure yet.
The proportion of sperm with cytoplasmic droplets is reduced in Ubl7−/− mice
The CD is a crucial cytoplasm-like structure attached to the sperm tail, associated with sperm maturation, energy supply and motility39–42. Deformations of CDs can also lead to abnormal bending of the sperm tail18,20,21,43–45. Using ubiquitin (Ub) to label CDs, we found that both caput and cauda sperm from Ubl7−/− mice exhibited a reduced Ub signal at the CD region (Fig. 3e, f). Moreover, we observed numerous free vesicles labeled by Ub in the cauda sperm smear of Ubl7+/− mice (Fig. 3e). We speculated that these vesicles were CDs detached from sperm tail during epididymal transit or sample preparation. Quantification revealed that the count of these CD-like vesicles was also lower in the cauda sperm smear of Ubl7−/− mice (Fig. 3g).
TEM analysis of CDs showed that numerous saccular elements and vesicles were present in most of the CDs (39 out of 52) from Ubl7+/− sperm (Fig. 3h), but only around one-third of CDs remaining in Ubl7−/− sperm (13 out of 34) appeared normal. The majority of CDs in Ubl7−/− sperm had fewer (9 out of 34) or no saccular elements (12 out of 34) and exhibited a reduced volume (Fig. 3h). To explore whether CD formation failed in Ubl7−/− spermatids, we performed whole-mount staining of stage VIII seminiferous tubules using a marker of CDs, phosphoglycerate kinase 2 (PGK2)40,42. Contrary to our expectations, CDs appeared to form normally in Ubl7−/− spermatids at step 16, albeit slightly smaller in size (Supplementary Fig. 4e). This result indicates that loss of UBL7 does not impede the formation of CDs in spermatids but impacts the morphology and alters the contents of CDs.
UBL7 binds to VCP/p97 complex through direct interaction with UBE4B
To understand the molecular mechanisms underlying the defects caused by UBL7 deficiency in sperm, we conducted a screen to identify interacting proteins of UBL7. We generated flag-Ubl7 knock-in mice using CRISPR-Cas9 (Supplementary Fig. 5a). The protein localization remained unchanged after insertion of the FLAG tag (Supplementary Fig. 5b, Fig. 1d). Furthermore, the sperm morphology was normal in Ubl7 flag/flag mice (Supplementary Fig. 5c), suggesting that the insertion of the FLAG tag in both alleles did not impair UBL7 function or spermatogenesis.
Next, we performed FLAG-UBL7 immunoprecipitation (IP) followed by mass spectrum using testes or spermatids (elongating and round spermatids, ES + RS) isolated from Ubl7 flag/flag and WT mice. Gene Ontology (GO) analysis of proteins identified in three biological replicates revealed that UBL7 predominantly interacted with proteins involved in endoplasmic reticulum associated degradation (ERAD) pathway, as well as processes related to protein folding and stabilization (Fig. 4a, Supplementary Tabel. 1). We validated several interacting proteins by WB analysis (Fig. 4b). Intriguingly, we observed extremely slowly migrating ubiquitinated proteins pulled down by FLAG-UBL7 IP (Fig. 4b), which are believed to be proteins with branched Ub chains that are prone to rapid degradation46.
The AAA ATPase valosin-containing protein (VCP) (also known as p97) plays a crucial role in protein degradation, especially in various cellular compartments, including endoplasmic reticulum, mitochondria, ribosome, and so on47–49. We identified VCP and several cofactors of it in FLAG-UBL7 IP. We wondered which components of the VCP complex directly interacted with UBL7. Pull down assays were carried out in 293T cells and results revealed that UBL7 specifically bound to the E4 ligase UBE4B, but not to VCP or its cofactors UFD1 and NPL4 (Fig. 4d, Supplementary Fig. 5d). To dissect the critical domains involved in this interaction, we generated truncated forms of UBL7 targeting its UBL and UBA domains (Fig. 4c). While UBL7-ΔUBL was unstable when expressed in cells (Fig. 4d), we constructed two other truncations, UBL7-ΔUBL-1 and UBL7-ΔUBL-2, by preserving specific regions of the UBL domain, and a site mutant (I44V70A) that disrupts the interacting surface of the UBL domain and the Ub binding domains22 (Fig. 4c). The results demonstrated that the affinity between UBL7-ΔUBA and UBE4B was similar to that of whole-length UBL7 (UBL7-WL), while loss of or disruption of the UBL domain disturbed it. Furthermore, UBE4B directly interacts with VCP (Fig. 4e). These data suggest that UBL7 was involved in VCP complex through direct interaction with UBE4B depending on its UBL domain.
Analogously, we conducted IP experiments of UBL7 or its truncations and Ub-conjugates in 293T cells. UBL7 effectively pulled down Ub-conjugates in cells. Notably, in cells treated with proteasome inhibitor MG132, more Ub-conjugates were co-precipitated, while in cells treated with Ub-activating enzyme E1 inhibitor TAK243, little Ub-conjugates were pulled down (Fig. 4f). Furthermore, the binding of UBL7 and Ub-conjugates was mediated by the UBA domain and deletion of the UBL domain released the UBA domain and enhanced the affinity (Fig. 4f). Based on these findings, we propose a schematic model in which UBL7 interacts with UBE4B and binds to ubiquitinated substrates translocated by the VCP complex (Fig. 4g).
UBL7 binds to testis 26S proteasome
Next, we investigated whether UBL7 binds to proteasome Ub receptors and act as a shuttling factor like other well-studied UBL-UBA proteins, such as Rad23 and Dsk250–52. In mammalian cells, the major intrinsic proteasome Ub receptors include PSMD2, PSMD4 and ADRM16. We performed IP in 293T cells to analyze the interaction between UBL7 and these receptors. The results showed that UBL7 specifically interacted with PSMD2 and PSMD4, but not with ADRM1 (Fig. 5a). Interestingly, in cells treated with MG132, the amount of PSMD2 interacting with UBL7 is unchanged, while in cells treated with TAK243, UBL7 failed to pull down PSMD2 efficiently (Fig. 5b). Consistent with this, the UBL7-ΔUBA, which lost its ability to bind to Ub-conjugates, also showed reduced interaction with PSMD2, and both UBL7-ΔUBL and UBL7-I44V70A pulled down more PSMD2 (Fig. 5b). Together, these results suggest that in cells, UBL7 interacts with the proteasome Ub receptors PSMD2 and PSMD4 and this interaction is dependent on the presence of Ub-conjugates.
To investigate the association between UBL7 and proteasomes in germ cells, we purified 26S proteasomes from testes using the UBL-affinity isolation method53 (Supplementary Fig. 6a). We confirmed the successful isolation of testis 26S proteasomes by WB detection of three proteasome subunits-PSMD4, PSMD2 and PSMB5 (Supplementary Fig. 6b, c). Silver staining further demonstrated that the affinity-purified products from mouse testes exhibited the classical distribution of proteasome subunits53 (Supplementary Fig. 6d). Moreover, we measured the peptidase activity of the products using N-Succinyl-Leu-Leu-Val-Tyr-7-Amido-4-Methylcoumarin (LLVY-amc) as a substrate. Results showed that affinity-isolated testis proteasomes efficiently hydrolyzed LLVY-amc and released the fluorescent amc, while in the presence of proteasome inhibitor MG132 the substrate remained undigested (Supplementary Fig. 6e). Thus, we successfully isolated functional 26S proteasomes from mouse testes under non-reducing conditions.
Using an in vitro pull-down assay, we discovered that purified UBL7 protein directly bound to the isolated testis 26S proteasomes in a concentration-dependent manner in absence of Ub-conjugates (Fig. 5c). Additionally, we added different concentrations of purified UBL7 protein to the proteasome degradation reaction using Boc-Leu-Gly-Arg-AMC (LGR-AMC) as a substrate. The peptidase activation curve indicated that UBL7 slightly stimulated the peptidase activity of testis proteasomes, resulting in an approximate 1.3-fold increase over the proteasome alone (Fig. 5d). However, when compared to a linear Ub chain (6Ub), which activated activity of testis proteasomes to a 2.4-fold increase, the impact from UBL7was negligible (Fig. 5e). Thus, we elucidate that UBL7 directly interacts with testis 26S proteasomes and induces a much smaller increase in proteasome activity than Ub chains.
UBL7 prevents substrates degradation in cells
Next, we investigated the impact of UBL7 on the degradation of proteasome substrates. Given the vital role of VCP complex in ERAD, we selected two typical ERAD substrates—Tyrosinase-C89R (TYR-C89R) and CD3δ—and measured their degradation rates in the presence or absence of UBL7 protein in 293T cells. We monitored their protein levels over time in cells under the treatment of cycloheximide (CHX) to assess their degradation rates. The results showed that the degradation rates of both TYR-C89R and CD3δ were significantly reduced in the presence of UBL7 (Fig. 6a-d). Additionally, the VCP complex has been implicated in the degradation of cytoplasmic substrates via the ubiquitin fusion degradation (UFD) pathway54–56. We constructed an artificial substrate for UFD pathway, UbV76-V-GFP, and found its degradation rate was not significantly changed by UBL7 (Fig. 6e, f). In conclusion, these CHX chasing assays demonstrate that UBL7 slow down the degradation rates of substrates associated with the VCP complex in cells and shows preference to substrates in organelles, such as in ER.
UBL7 safeguards critical factors from excessive degradation during spermiogenesis
To explore whether UBL7 could also protect substrates in spermtids, we performed tandem mass tag (TMT)-labeled quantitative proteomics of sperm from Ubl7+/− and Ubl7−/− epididymis. Overall, 31,121 unique peptides were detected and 5,034 proteins were identified in our samples. Among these, 3,886 proteins contained quantitative values in all three replicates. We observed that 405 proteins showed significantly reduced abundance in Ubl7−/− sperm (with a 1.5-fold cutoff), while 22 proteins showed increased abundance (with a 1.5-fold cutoff) (Fig. 6g, Supplementary Table 2). GO analysis of the reduced proteins in Ubl7−/− sperm identified functional enrichment of terms such as “glycolytic process”, “ATP metabolic process”, “glutathione metabolic process”, and “flagellated sperm motility”, among others (Supplementary Fig. 7a). These processes shared a large overlap with the processes in which proteins enriched in CDs participate40. Then we noticed that most of the CD-enriched proteins reported previously were reduced in Ubl7−/− sperm (Supplementary Fig. 7b). We validated the levels of several proteins by WB, including SCCPDH and DNPEP (amino acid metabolism), HK1, SDHA and PDHB (ATP metabolism), MAN2C1 (post-translational modification of proteins), PSMD2 and SYPL142 (CD formation). All of the proteins showed significantly lower levels in Ubl7−/− sperm compared to Ubl7+/− sperm (Fig. 6h, i). Immunofluorescent staining showed that many of the reduced proteins localized to CDs, such as DNPEP, MAN2C1 and SYPL1 (Fig. 6j and Supplementary Fig. 7c). Some others were enriched in CDs and were also detectable in sperm tails, such as TCP1 and HK1 (Fig. 6j and Supplementary Fig. 7c). Proteins exclusively localized to sperm tails also existed, such as SCCPDH (Fig. 6j). All of the proteins we validated exhibited lower levels in Ubl7−/− sperm, no matter in CDs or in sperm tails (Fig. 6j and Supplementary Fig. 7c).
The quantitative proteome analysis indicated that UBL7 deficiency led to loss of numerous functional proteins in sperm, suggesting UBL7 probably protected critical factors from degradation. To investigate the direct substrates targeted by UBL7, we set up a two-step IP method to enrich UBL7-interacting Ub-conjugates. For short, we performed FLAG-UBL7 IP in Ubl7 flag/flag testes after crosslinking with dithiobis(succinimidyl propionate) (DSP). Then the elutions were subjected to Tandem Ubiquitin Binding Entities (TUBE)-IP to enrich Ub-conjugates in the UBL7 interactome, as previously described 57. WB analysis confirmed that Ub-conjugates binding to UBL7 were successfully isolated (Fig. 6k). In total, we identified 142 proteins that were enriched by the two-step IP in Ubl7 flag/flag testes compared to WT testes. Excluding proteins nonspecifically isolated due to their extremely high abundance in FLAG-IP elution, many other proteins in the list have known functions in spermiogenesis (Supplementary Table 3). We classified these proteins into nine categories based on their functions or subcellular localization (Fig. 6l). For example, IFT88, IFT140, IFT122, CCT3, CCT5 and GCAP14 are known to be involved in manchette development58–65. SPATA20 and IFT88 contribute to the integrity of HTCA66 67. Additionally, plenty of proteins enriched in CDs were also detected by the two-step IP, including HK1, PEBP1, GAPDHS, IFT88, ATP5a1, SUCLA2, SLC2a3, LDHAL6b, PSMD3 and PSMD1. Apart from proteins related to the morphological defects of Ubl7−/− sperm, proteins involved in anit-oxidative stress (GSTM5, ALDH2 and LANCL1) and glycolysis pathway (PFKP, LDHAL6b, HK1, TKTL1, PKM2, GAPDHS, LDHA, ACYP1, ENO1, PGAM2 and PFKM) also stood out. Their loss might impact energy metabolism and motility of sperm. In addition, based on subcellular localization, many proteins detected by the two-step IP were mitochondrial proteins (ERCC6L2, DLAT, LDHAL6b, ACO2, ATP5a1, AKR1b3, SUCLA2), Golgi-resident proteins (VPS26b, VPS16, TMEM87a), ER-resident proteins (POMT1, ATP2a1, ATP2a2, BSCL2, SPPL2c) or plasma membrane proteins (SLC2a3, SLC7a5, ATP1a3, ATP1a4, SLMAP, DGKK). Degradation of these proteins may rely on translocation by the VCP complex, thereby may be affected by UBL7. Furthermore, we compared protein levels of the targets identified by the two-step IP in Ubl7−/− sperm with that in Ubl7+/− sperm. 81 out of the 142 targets were detectable in the TMT quantitative proteome, and the heat map showed that the majority of the putative UBL7 targets tended to decrease in Ubl7−/− sperm (Fig. 6n). Bolstered by these results, we speculate that UBL7 protects functional proteins in sperm from undue degradation, thereby ensuring successful spermiogenesis.