Introduction

Obesity and associated metabolic disorders, such as insulin resistance and Type 2 Diabetes (T2D), are major sources of morbidity and mortality that are reaching epidemic proportions (Cawley et al., 2021; Kelly, Yang, Chen, Reynolds, & He, 2008; Stokes & Preston, 2017; WHO, 2016, 2020; Zimmet, Alberti, Magliano, & Bennett, 2016). Lifestyle interventions, such as diet and exercise, can be effective in combating these afflictions in some patients. However, low adherence and counter-regulatory mechanisms that limit weight reduction and long-term weight maintenance impede the effectiveness of such interventions. (Brownell, 1998; Miller et al., 2002; Roberts & Barnard, 2005). The new anti-obesity drugs, semaglutide, a GLP-1 receptor agonist, and tirzepatide, a dual GIP and GLP-1 receptor agonist, have demonstrated remarkable efficacy for weight loss and hold great promise as anti-obesity treatments to alleviate metabolic dysfunction (Venniyoor, 2022; Wilding et al., 2021). However, even with their recent success, there are many patients who do not respond to GLP-1 receptor agonist treatment or have to discontinue the treatment due to side effects (Chao, Tronieri, Amaro, & Wadden, 2022). Thus, there remains a need for additional new therapies to combat obesity.

To identify new targets for the treatment of obesity and diabetes, we began a large-scale forward genetic phenotypic screen in mice (T. Wang et al., 2015). Random mutations in the mouse genome were generated using the chemical ENU (N-ethyl-N-nitrosourea). The genomes of the animals were sequenced, and various metabolic endpoints were evaluated to identify genes with mutations associated with obesity and insulin resistance phenotypes (Fig. 1A). Through this process, mutations in the gene spag7 (sperm-associated antigen 7) were identified. SPAG7 is well-conserved with 97% of the amino acid sequence being identical in humans and mice. Structurally, it contains a nuclear localization signal, making it likely a nuclear protein, and an R3H domain (Nagata, 2005). The function of the R3H domain is also not fully understood; however, it is predicted to bind polynucleotides, and other proteins that contain R3H domains are known to bind single-stranded DNA or RNA (Grishin, 1998; He, Zhang, Liu, & Yan, 2013).

Figure 1.

The biological function of SPAG7 has been largely unstudied. It was first identified in the inner acrosomal compartment of fox sperm (Beaton, Cleary, ten Have, & Bradley, 1994). However, the gene is expressed in every tissue and cell in humans and mice (Noguchi et al., 2017). It has been implicated in various disease states, including Asperger syndrome (Tentler et al., 2003), parvovirus B19 infection (Kerr et al., 2005), synovial sarcoma (Fernebro et al., 2006), azoospermia (Ayhan et al., 2014), periodic fever, aphthous stomatitis, pharyngitis, and adenopathy (PFAPA) syndrome (Bens et al., 2014), squamous cell carcinoma (Singh et al., 2020) and oligoasthenozoospermia (Abu-Halima et al., 2023). However, the role SPAG7 plays in these diseases remains elusive.

Mice carrying SPAG7 mutations displayed interesting metabolic phenotypes, suggesting that SPAG7 may play a role in regulating metabolism in vivo. The goals of this study are to determine whether SPAG7 loss of function affects systemic metabolism and to establish a mechanism of action for SPAG7.

Results

SPAG7-deficient animals are obese and insulin resistant

To investigate the role of SPAG7 in metabolism, we generated a germline SPAG7-deficient mouse model (SPAG7 KO). SPAG7 KO animals are obese in adulthood (Fig. 1B). But they have an unexpected growth curve; they are underweight early in life, catch up around six weeks of age, and then continue to outgrow WT littermates into and throughout adulthood (Fig. 1C). SPAG7-deficient mice demonstrate decreased lean mass early in life that remains below that of WT littermates throughout their lifespan (Fig. 1D). In contrast, fat mass in the animals increases continually, growing to three-to-five-times the fat mass of WTs (Fig. 1E). SPAG7-deficient animals display decreased body length, measured nose to base of tail (Fig. 1F). Bone mass density is also decreased in SPAG7 KO animals, which may contribute to the decreased lean mass (Fig. 1G).

The majority of the increased body weight in SPAG7-deficient animals stems from increased adipose tissue mass (Fig. S1A-D). All adipose tissue depots more than double in size (Fig. S1A-D). SPAG7-deficient adipose tissues display adipocyte hypertrophy and increased immune cell infiltration, consistent with chronic obesity (Fig. S1E-H). However, heart, skeletal muscle, kidney, and brain weights are all decreased in SPAG7-deficient animals (Fig 1H-L). Livers display no change in mass but do have increased triglyceride content, indicating hepatic steatosis (Fig. 1M, Fig. S1I-K). No changes in circulating liver enzymes were observed, indicating that the lipid accumulation was not severe enough to induce significant liver damage (Fig. S1L-O). Plasma cholesterol, triglycerides, and free fatty acids are all increased, consistent with an obesity phenotype (Fig. 1N-S).

SPAG7-deficient animals display severely impaired glucose uptake in response to oral glucose gavage (Fig. 1T). Glucose-stimulated insulin secretion is increased in these animals, consistent with an insulin resistance phenotype (Fig. 1U). The animals also display an insulin intolerance phenotype, with decreased glucose uptake following an insulin injection (Fig. 1V). This is, perhaps, unsurprising, given the rather extreme fat mass phenotype.

Obesity in SPAG7-deficient animals is caused by decreased energy expenditure

A change in food intake was not detectable in SPAG7-deficient animals (Fig. 2A). In fact, a trend towards decreased food intake was consistently observed. When placed in metabolic cages, SPAG7-deficient animals display a significant decrease in energy expenditure that is present during the light cycle, but particularly clear in the dark cycle, when we would expect the animals to be active and moving about the cage (Fig. 2B, 2C). Consistent with this finding, SPAG7-deficient animals display decreased locomotor activity, as measured by home-cage beam break (Fig. 2D).

Figure 2.

To determine if thermogenesis also contributes to the decreased energy expenditure observed in the SPAG7 KO mice, SPAG7-deficient animals were housed at thermoneutrality (TN). At TN the animals display the same phenotypes: increased body weight in adulthood, dramatically increased fat mass, no change in food intake, and a decrease in total energy expenditure (Fig. 2E-H). This indicates that thermogenesis is unlikely to play a significant role in the decreased energy expenditure and weight gain observed in the animals. SPAG7-deficient animals fed a 60% high-fat diet (HFD) also display the same phenotypes with no increase in food intake in SPAG7 KOs compared to WTs (Fig. 2I-L). The percent body weight gain in SPAG7-deficient animals vs WT was very similar at room temperature on normal chow, at thermoneutrality on normal chow, and at room temperature on HFD: an increase of 15-17% (Fig. 2M), indicating that decreased energy expenditure due to reduced locomotor activity is the key driver for the obesity phenotype in the SPAG7 KO mice.

SPAG7-deficient animals display decreased exercise capacity and muscle function

The observation of reduced locomotor activity in the SPAG7 KO mice prompted us to investigate if there are defects in exercise performance and muscle function in these mice. WT and SPAG7 KO animals were subjected to treadmill endurance tests; the KO mice ran a significantly shorter distance than the WT animals before becoming exhausted (Fig. 3A). Additionally, the maximum oxygen consumption reached during the treadmill endurance test was significantly reduced in SPAG7-deficient animals (Fig. 3B). The maximum force produced by gastrocnemius complex hindlimb muscles when stimulated in anaesthetized animals was measured as indicative of muscle function and is severely impaired in SPAG7-deficient animals (Fig. 3C). Skeletal muscle histology indicates a significant reduction in the cross-sectional area of Type IIa (fast-twitch) muscle fiber (Fig. 3D, S2A), indicating muscle atrophy. SPAG7-deficient skeletal muscle was also found to contain elevated levels of triglyceride content (Fig. 3E). Skeletal muscle vascularization appeared normal (Fig. 3F-H). We did, however, observe a marked reduction in oxidative capacity in skeletal muscle, as indicated by succinate dehydrogenase b activity staining and citrate synthase activity in gastrocnemius muscle (Fig. 3I, 3J, S2B, S2C).

Figure 3.

Transcriptomics analysis of gastrocnemius muscle revealed 1,848 genes were significantly upregulated and 2,291 genes were downregulated in SPAG7 KO muscle compared to WT (Fig. 3K, S3A, S3B). Gene ontology enrichment analysis revealed many gene expression pathways that were affected by SPAG7-deficiency; among the most significantly downregulated were pathways involved in mitochondrial function (Fig. 3L). The most significantly downregulated gene was hsd17b10. This is a ubiquitously expressed mitochondrial protein that acts as both a short-chain dehydrogenase reductase and as a part of the RNase P complex, which is responsible for the maturation of all mitochondrial tRNA (Holzmann et al., 2008; Vilardo et al., 2012; Yang, He, & Schulz, 2005). Missense mutations in hsd17b10 occur in humans and result in HSD10 disease (Chatfield et al., 2015). HSD10 disease varies significantly in severity, but the most severe forms are characterized by developmental delays, decreased muscle tone, and neonatal mortality (Ofman et al., 2003; Vilardo & Rossmanith, 2015; Zschocke, 2012). Hsd17b10 expression is reduced in nearly all SPAG7-deficient tissues (Fig. S3C-H). Overall, these studies imply significant impairment of mitochondrial function and muscle performance in SPAG7-deficient animals.

Induction of SPAG7-deficiency after early development does not cause obesity or skeletal muscle abnormalities

SPAG7 KO mice were smaller than WT mice during the neonatal stage and profoundly obese in adulthood (Fig. 1C). This raised the question whether SPAG7-deficiency affects embryonic development, which may potentially cause the metabolic perturbations later in life. To test this hypothesis, an inducible whole-body SPAG7 KO (iSPAG7) mouse was developed. Using a globally expressed Cre-ERT2 along with a globally floxed SPAG7 gene, tamoxifen injection would induce global recombination and deletion of the SPAG7 gene (Fig. 4A). To test the effects of SPAG7-deficiency in adulthood, iSPAG7 animals were aged to 7 weeks before being dosed with tamoxifen. In this manner, tamoxifen dosing was able to establish effective knockout of SPAG7 protein levels in all tissues examined, including liver, brain, kidney, muscle, and fat (Fig. 4B). In contrast to the germline SPAG7 KO, iSPAG7 KO animals displayed no differences in body weight, lean mass, fat mass, or food intake (Fig. 4C-F). There is no difference in glucose tolerance and total energy expenditure between the iSPAG7 KO and the WT mice (Fig. 4G-J). Home cage locomotor activity, hindlimb max force generation, and treadmill endurance were all normal, as well (Fig. 4K-M). These results indicate that the metabolic phenotypes are downstream of the effects of SPAG7-deficiency on the developing animal, suggesting that SPAG7 may play a role during embryonic development.

Figure 4.

SPAG7-deficiency causes placental insufficiency and intrauterine growth restriction

SPAG7-deficient mice are born significantly underweight (Fig. 5A. 5B). They are also born below the expected rate. In heterozygous x heterozygous breeding, Mendelian ratios would predict roughly 25% of pups born to be SPAG7-deficient. We observe roughly 10% of pups to be SPAG7-deficient, indicating significant loss of SPAG7-deficient animals prior to birth (Fig. 5C). Throughout mouse (and human) development there are key stages where significant loss of embryos is more likely to occur. These include fertilization (e0), gastrulation (e6.5-e8), organogenesis (e8-e11.5), and birth (e19-e20). Genotyping embryos before and after each of these developmental stages, we observe expected Mendelian ratios of SPAG7-deficient animals from e6.5 through e11.5, indicating no loss of embryos during fertilization, gastrulation, and organogenesis (Fig. 5D-F). We do see a significant decrease in SPAG7-deficient percentage at e18.5, indicating that SPAG7-deficient embryos are lost at some point after e11.5 and prior to birth (Fig. 5G). However, it should be noted that SPAG7-deficient embryos appear to be smaller and less well-developed than WT embryos within the same dam, at each of these stages.

Figure 5.

At e18.5, SPAG7-deficient animals display decreased fetal and placental weight (Fig. 5H-J), indicating intrauterine growth restriction (IUGR). Other than the smaller size, the SPAG7-deficient fetus lacks any specific morphologic abnormalities: all organs and limbs appear appropriately developed for their size. To further characterize the IUGR phenotype, we measured markers relevant to IUGR in the fetus. SPAG7-deficient fetuses have lower blood glucose and plasma insulin compared with WT controls (Fig. 5K, 5L). Insulin-like growth factor (IGF) signaling has been implicated as essential for fetal growth and is dysregulated in IUGR fetuses (Agrogiannis, Sifakis, Patsouris, & Konstantinidou, 2014; Reid, Flozak, & Simmons, 2002). SPAG7-deficient fetal liver displays decreased IGF1 and 2 protein and gene expression levels, indicating reduced IGF signaling (Fig. 5M-O). These findings are consistent with an intrauterine growth restriction phenotype (IUGR) caused by placental insufficiency and suggest SPAG7 may be necessary to maintain placental function.

SPAG7 is expressed highly in mouse and human placenta (Uhlen et al., 2015), it is also highly expressed in every mouse tissue during the earliest stages of development (Baldarelli et al., 2021; Brunskill et al., 2014; Schmitt et al., 2014). SPAG7-deficient placenta is smaller in cross-sectional area compared to the placenta of WT littermates (Fig. 6A-C). SPAG7-deficient placenta displayed no significant changes in vascularization as measured by CD34 IHC staining and gene expression of vascularization markers (Fig. 6B, 6D, 6E). There is a significant increase in mitochondrial DNA content in SPAG7 KO placenta (Fig. 6F). This finding is consistent with previous findings in IUGR and is generally thought to be a compensatory response (Mando et al., 2014; Naha et al., 2020). Histological analyses showed that SPAG7 KO placenta displayed abnormal morphology, exhibiting disruption of both the labyrinth and junctional zones. The junctional zone, however, was particularly severely affected, as demonstrated by decreased junctional:labyrinth zone ratios and percent junctional zone (Fig. 6G-I). The junctional zone of the placenta plays a key role in the production of hormones and the transport of metabolites from dam to pup (Woods, Perez-Garcia, & Hemberger, 2018). A disruption in this area is likely to reduce nutrient availability to the fetus.

Figure 6.

Discussion

We identified mutations in the gene spag7 that induce obesity and insulin resistance in mice from a mouse forward genetic screen. The biological function of SPAG7 is largely unknown. Utilizing constitutive and inducible SPAG7 knockout mouse models, we investigated the function of SPAG7 in vivo. The results show that SPAG7-deficient mice were born with lower body weight; however, they developed obesity and insulin resistance in adulthood. The metabolic disturbance in the SPAG7 KO animals is due to a decrease in energy expenditure, driven by decreased locomotor activity. Skeletal muscle function and exercise endurance in SPAG7-deficient animals are significantly impaired, displaying decreased force generation and oxidative capacity. In addition, these phenotypes are downstream of developmental effects, as knocking out SPAG7 in adulthood did not cause these metabolic disturbances. Furthermore, SPAG7 KO embryos display intrauterine growth restriction (IUGR), and aberrant development of the junctional zone of SPAG7-deficient placenta. These results demonstrate that SPAG7 is essential for fetal growth and placental function. Lack of SPAG7 results in IUGR, which leads to the development of obesity, insulin resistance, and other metabolic disturbances later in life.

Fetal nutritional status is well-known to drive metabolic health in adulthood. Perhaps the most famous example of this is the Dutch Hunger Winter. Pregnant women, exposed to acute and extreme decreases in nutritional intake, gave birth to low birthweight newborns, who were then at higher risk for developing obesity and diabetes in adulthood (Roseboom, de Rooij, & Painter, 2006). Intrauterine growth restriction (IUGR) has been linked to the development of metabolic syndrome later in life. However, the mechanisms that cause IUGR and influence its association with adult metabolic health remain poorly understood. Bilateral uterine artery ligation during gestation in rats leads to decreased birth weights, but obesity in adulthood (Simmons, Templeton, & Gertz, 2001). Placental insufficiency is a common cause of IUGR (Malhotra et al., 2019; Wood, 2020). In fact, junctional zone deficiency, in particular, has been shown to negatively affect fetal growth and birth weight (Mark, 2011; Woods et al., 2018).

IUGR has been shown to affect skeletal muscle health and performance. Humans that experience IUGR, for any given population, display decreased skeletal muscle mass that is not compensated for after birth and persists through adulthood (Nastase, Cretoiu, & Stoicescu, 2018). The skeletal muscle of IUGR fetuses in sheep display decreased intracellular ATP content and decreased amino acid uptake (Sarkar et al., 2014; Stremming, Jansson, Powell, Rozance, & Brown, 2020). Muscle mitochondrial oxidative capacity can also be lower in IUGR skeletal muscle (Stremming et al., 2022).

The deleterious effects of SPAG7-deficiency on the fetus clearly indicate an important function for SPAG7 during embryonic development. The most common cause for IUGR is placental insufficiency. In the rodent, the placenta can be roughly broken up into three zones: the decidua, the labyrinth zone, and the junctional zone. The decidua and labyrinth zones contain the vasculature of the mother and fetus, respectively. The junctional zone is the main site of nutrient and gas exchange and also the main endocrine compartment of the placenta, producing hormones, growth factors, and cytokines that are important for the normal progression of pregnancy and fetal growth (Woods et al., 2018). Histological examination of SPAG7 KO placenta revealed abnormal morphology of both the labyrinth and junctional zones. In addition, the junctional zone was markedly diminished in SPAG7 KO placenta, as demonstrated by decreased junctional:labyrinth zone ratios and decreased percent junctional zone (Fig. 7G-I). The mature mouse placenta is established mid-gestation (∼e10.5) and continues to grow, in size and complexity, throughout the pregnancy. Early in gestation, the yolk sac supports the nutritional needs of the embryo, however, from mid-gestation onward, the fetal nutrient supply requires the transport capacity of the placenta. The success of this transition is critical for fetal survival and continued growth. During our studies, we did not observe significant loss of SPAG7 KO embryos until after the switch to placental nutrition supply, which could be due to the disruption of labyrinth and junctional zone development in the placenta.

Figure 7.

We hypothesize the SPAG7 KO fetuses that survive do so by decreasing the mitochondrial metabolic capacity of their skeletal muscle and other tissues. Reducing muscle energy consumption in response to the decreased nutrient availability in the IUGR fetus is advantageous to ensure fetal survival. However, when these animals are born and age, this reduction predisposes them towards metabolic disease (Pendleton, Wesolowski, Regnault, Lynch, & Limesand, 2021). After birth, when nutrition is readily available, the decrease in energy expenditure, combined with food intake normalized to wild-type littermates, eventually leads to increased body weight and fat mass. It is worth noting that SPAG7 KO mice are not hyperphagic, a phenotype observed in some IUGR models, suggesting that other perturbations in SPAG7 KO animals may counter-regulate hyperphagia. The key driver for the development of obesity and metabolic syndrome in the SPAG7 KO mouse is reduced energy expenditure. Models that produce obesity, primarily driven by decreased locomotor activity, are rare. Further investigation of SPAG7-deficient mice may provide valuable information for drugs targeting energy expenditure in the treatment of obesity.

Overall, our results demonstrate that SPAG7 plays an important role in fetal development. SPAG7 is essential for the proper formation of the junctional and labyrinth zones of the placenta, which contributes to IUGR in the SPAG7 KO fetus (Fig. 7). The surviving SPAG7 KO mice develop a range of metabolic disorders including obesity, insulin resistance, glucose intolerance, exercise impairment, reduced energy expenditure, and reduced skeletal muscle mitochondrial oxidative capacity in adulthood. Our results uncover, for the first time, the biological function for the spag7 gene in vivo. Furthermore, these findings shed light into a new biological pathway regulating fetal development. The results provide insights into how skeletal muscle function contributes to systemic energy balance and how fetal nutritional status can have far-reaching effects, influencing the nutritional status of the adult. Deeper understanding of the mechanisms involved in placental insufficiency, IUGR, and SPAG7 biology may give us the opportunity to reverse aberrations in fetal growth rate and prevent the metabolic co-morbidities associated with low birth weight (Wood, 2020).

Materials and methods

Mice

All animal experiments were conducted following study protocols and procedures reviewed and approved by Pfizer Institutional Animal Care and Use Committee. The facilities that supported this work are fully accredited by AAALAC International.□

Adult female wild-type (WT), SPAG7 KO, (C57BL/6J;□various ages) were obtained from Jackson Laboratories (Farmington, CT). Mice were housed individually in Innovive cages (Innorack IVC Mouse 3.5) under a standard 12-h light/12-h dark cycle (06:00 h: 18:00 h) in a temperature- and humidity-controlled environment (22□± 1°C for room temperature, 28 ± 1°C for thermoneutrality). Mice were given□ad libitum□access to tap water and standard chow (Innovive mouse metal feeder; Purina rodent diet 5061; Purina Mills, St. Louis, MO) or 60% high-fat diet (Research Diets D12492, New Brunswick, NJ). Body weight was recorded weekly (Mettler Toledo ME4002TE, Mettler Toledo, Oakland, CA). Body composition was assessed using the EchoMRI 4-in-1 500 Body Composition Analyzer (EchoMRI, Houston, TX).

Generation of Spag7 cKO mouse model using the CRISPR/Cas9 technology

The Spag7 cKO model was created on the C57BL/6J strain (JAX #000664) by inserting two LoxP sites flanking the exon 2 of the gene in the genome utilizing the CRISPR/Cas9 technology (H. Wang et al., 2013). One LoxP is positioned at 316 bp upstream of the exon 2 and the other at 97 bp downstream of the exon 2. A donor plasmid used as a repair template consists of the floxed genomic sequence flanked by 0.7 kb 5’ homology arm and 0.4 kb 3’ homology arm. Protospacer sequences of the two sgRNAs targeting the upstream and downstream insertion sites are 5’-CCTGAAGCCAGTGTATATTG and 5’-GGGCTTTTAGTTCCACCATC, respectively. These sgRNAs were synthesized by in vitro transcription. For genome editing in mouse zygotes, Cas9 protein (IDT, #1081058), sgRNAs and the donor plasmid were mixed and microinjected into 1-cell embryos harvested from super-ovulated donors. Injected embryos were transferred to pseudo-pregnant females to generate live pups. For genotyping of founders and their offspring, genomic PCR was performed with two primers external to the 5’ and 3’ homology arms (5’-TCACATCACGGTCCATCATC and 5’-TCATGACATAGCGACAGTCA, respectively) and amplified products were sequenced to confirm correct editing. Selected founders were bred to wildtype C57BL/6J mice for germline transmission of the cKO allele.

Western Blotting

Animal tissues were snap-frozen in liquid nitrogen and stored at□−80□°C. Frozen tissues were homogenized in ice-cold RIPA buffer (Sigma Aldrich, St. Louis, MO), and protein concentrations were determined using a BCA protein assay (Thermo Fisher Scientific, Waltham, MA). Protein extracts were separated on NuPAGE 4–20% Bis-Tris gels (Bio-Rad, Hercules, CA) and blotted onto PVDF membranes (Bio-Rad, Hercules, CA). Membranes were blocked for 1□h at room temperature in TBST (0.1%) containing 5% milk. Membranes were then incubated overnight at 4□°C with primary antibodies. Following three washing steps with TBST (0.1%), membranes were incubated with HRP-conjugated secondary antibodies for 1□hour at room temperature. After thorough washing, proteins were visualized with SuperSignal West Dura Extended Duration Substrate (Thermo Fisher Scientific, Waltham, MA) on the Bio-Rad ChemiDoc MP. Immunoreactive bands were quantified using Image J Software (NIH). Primary antibodies used: Anti-SPAG7 (Sigma, HPA024032), Anti-IGF1 (ThermoFisher, MA5-35060), Anti-IGF2 (ThermoFisher, PA5-71494), and Anti-bACTIN (ThermoFisher, MA5-32540)

Quantitative RT-PCR

For mouse studies, animals were euthanized via CO2 and tissues were immediately dissected and snap frozen in liquid nitrogen. Ribonucleic acid (RNA) was extracted and purified using RNeasy RNA Isolation Kit (QIAGEN, Germantown, MD) then reverse transcribed into complementary deoxyribonucleic acid (cDNA) with a High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, Foster City, CA). Quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) was performed using TaqMan reagents and primer-probes (Applied Biosystems, Foster City, CA). Gene expression was normalized to the control gene TATA box binding protein (Tbp).

Bone Mineral Density and Body Length

Dual X-ray absorptiometry (DEXA) imaging was used to assess body composition, bone mineral density and body length of SPAG7 knock out and litter mate or age-matched control mice. The imaging was conducted on the Ultra-Focus system Faxitron (BioVision). The system was calibrated and warmed up prior to image acquisition. Animals were anesthetized using isoflurane anesthesia (1 to 4%), and imaged in the prone position, with paws stretched and taped to maintain position. Whole body images were captured using both energies (40 and 80 kVP). Body composition was measured from DEXA images of each animal using software embedded tools (Vison DEXA). Overall length of each was animal measured from tip of mouse nose to the distal end of the tail using in built DEXA image analysis tool. Screenshots of data were saved and later transcribed to an excel sheet.

Circulating Lipid and Liver Enzyme Analysis

At 13 weeks of age, SPAG7 KO animals were bled via tail nick and blood was collected in EDTA-coated tubes (Becton, Dickinson and Company, Franklin Lakes, NJ) and centrifuged to collect plasma (2,000g for 10 min at 4°C), which was stored at −80°C until analysis. For fasted measurements, animals were bled following a 16 hour overnight fast. For refed measurements, animals were bled following a 2-hour refeed following a 16 hour overnight fast. Plasma was analyzed using the ADVIA Chemistry Analyzer XPT/1800 (Siemens, Munich, Germany). Plasma triglycerides (ADVIA XPT, Siemens, Munch Germany), cholesterol (ADVIA XPT, Siemens, Munch Germany), and non-esterified fatty acids (ADVIA XPT, Siemens, Munch Germany) were analyzed in all mice.

In Vivo Metabolic Studies

All metabolic tests were performed using standard protocols on female littermates (aged: 6-32 weeks). Mice were acclimated to handling and injection (0.1□mL saline) or oral gavage (0.1 mL saline) stress for ∼5□days prior to the start of all studies. For glucose tolerance tests (GTTs), animals were fasted for 8 hours with ad libitum access to water. Fasting glucose measurements were taken before 2mg/kg oral glucose dose. Blood glucose measurements were taken via tail nick using the AlphaTRAK 2 Glucose Monitoring System (Zoetis, Parsippany-Troy Hills, NJ) at 15, 30, 60, and 120 minutes after dose. Blood was collected for insulin measurement in EDTA tubes (Becton, Dickinson and Company, Franklin Lakes, NJ) and centrifuged to collect plasma (2,000g for 10 min at 4°C), which was stored at −80°C until analysis at 0, 15, and 30 minutes after dose. Insulin levels were measured using an Ultrasensitive Mouse Insulin ELISA (ALPCO, Salem, NH), according to manufacturer’s instructions. For insulin tolerance tests (ITTs), animals were fasted for 8 hours and blood glucose was measured before (0 minutes) and at 15, 30, 60, and 120 minutes following IP injection of .75 IU/kg of insulin (Humulin R, Eli Lilly, Indianapolis, IN).

Adipose Tissue Analysis

At 32 weeks of age, SPAG7 KO animals were sacrificed, perigonadal/visceral, inguinal/subcutaneous, and intrascapular/brown adipose tissues were removed and weighed. They were then fixed in 10% formalin, processed to paraffin, sectioned at 5 µm, stained with H&E, and scanned at 40x magnification using the AxioScan7 (Zeiss, Jena, Germany). Two 1000 µm x 1000 µm samples from homogeneous adipocyte areas of each image were sent to ImageJ and further analyzed with the Adipocyte Tools plugin using the option for large cells and the settings of minimum size 500 µm2 and number of erosions rounds 3. Faulty detections were manually corrected. Cell diameter was calculated assuming circular shape of adipocytes. Crownlike structures (CLS) were identified in subcutaneous AT sections as single adipocytes surrounded by at least 4 F4/80-positive macrophages.

Triglyceride Content Assays

32-week-old WT and SPAG7 KO animals were taken down and gastrocnemius muscle and liver were dissected and snap-frozen in liquid nitrogen. 50-100mg of pulverized tissue was homogenized and analyzed for triglyceride content using the Roche Hitachi 912 Analyzer Triglyceride Reagent Set (GMI, Ramsey, MN).

Food Intake

Food intake measurements were made using in-cage hopper weights or in-cage ceramic bowls, measured weekly. Food intake was further monitored using the BioDAQ food and water intake monitoring system system (Research Diets, New Brunswick, NJ). Mice were acclimated for four days in the BioDAQ system before food intake experiments were run. Food intake was monitored continuously and recorded at time intervals indicated in the figures.

Metabolic Cage Studies

Energy expenditure (EE), and locomotor activity were measured with an Oxymax indirect calorimetry system and CLAMS monitoring system (Columbus Instruments, Columbus, OH) starting at 11 weeks of age. The animals were given free access to water and diet. Means were calculated for dark, light, and full day cycles. Animals were acclimatized to the cages for 48□hours prior to data collection.

Treadmill Endurance Test

Mice were progressively acclimatized to treadmill exposure by increasing intensity and duration for 3 days. For the incremental test, mice ran on an enclosed, single lane treadmill (molecular enclosed metabolic treadmill for mice; Columbus Instruments), and real-time measurements of oxygen consumption (VlJO₂) and carbon dioxide output (VlJCO₂) were performed using an Oxymax/Comprehensive Laboratory Animal Monitoring System (CLAMS; Columbus Instruments). Mice ran at 6, 9, 12, 18, 21, and 23□m·min^-1 for 3 min at each velocity at a 0°, 5°, 10°, and 15° inclination for 6□min at each inclination. Exhaustion was defined as the inability to return to treadmill running after 10 seconds.

In Vivo Muscle Max Force Generation

Mice were anaesthetized with isoflurane and placed supine on a platform heated via a circulating water bath at 37°C. The right leg was shaved up to the patella and right knee stabilized via knee clamp. Once stabilized, the right foot was affixed to a Dual Mode Foot Plate (300-C FP, Aurora Scientific Inc., Aurora, Canada), and two electrodes were placed subcutaneously near the mid-belly of the gastrocnemius to achieve plantar flexion. A 1 Hz electrical stimulation was delivered (0.2 s duration, 1 s between stimulations) via stimulator (701C, Aurora Scientific Inc.) while increasing amperes to generate a maximum twitch measurement. After a maximum twitch was established, a force frequency of isometric contractions was initiated (0.2 s duration, 120 s between stimulations). All data were collected and analysed using the manufacturer supplied software (DMC and DMA, Aurora Scientific Inc.). Mice were housed at thermoneutral temperature prior to analysis.

Histology and Immunohistochemistry

Formalin fixed adipose tissue, livers, and gastrocnemius muscles were dehydrated and embedded in paraffin. One paraffin-embedded cross-section (5 μm thickness) through the adipose tissue, the left lateral lobe of the liver, and the mid-body of the gastrocnemius muscle from each animal was stained with hematoxylin and eosin using an automatic slide stainer (Leica) and microscopic evaluation was performed by a veterinary pathologist.

For muscle vascularity staining, mouse gastrocnemius muscle was isolated, mounted on a cork and snap-frozen in cold isopentane (Sigma-Aldrich, 277258) for cryo-sectioning. Cross-sections of 10-μm thickness were collected from the belly of the samples and placed on Superfrost Plus microscope slides. For capillary density and fiber size evaluation, the sections were stained with CD31 antibody (Biocare, CM-303A) and Dylight 594 tomato lectin (Vector Laboratories, DL-1177) for 1 hr under room temperature followed by an 1-hr incubation of an alexa 488 goat anti-rat secondary antibody (Thermo Fisher, A-11006). The slides were mounted with VECTASHIELD Antifade Mounting Medium with 4′,6-diamidino-2-phenylindole (DAPI) (Vector Laboratories, H-1800) and were imaged using a Zeiss AxioScan Z.1 slide scanner (Carl Zeiss, Jena, Germany). The images were analyzed using Visiopharm (Version 2020.09.0.8195) and custom-designed applications. The capillary/fiber ratio was calculated by dividing the total count of capillary (CD31+ lectin+) by the total count of fiber in the cross-section. The fiber size was the median surface area of all fibers in the cross-section.

For placenta vascularity staining, paraffin-embedded sections (5 μm thickness) were deparaffinized and rehydrated. Sections were blocked in 5% goat serum for 1 hour at room temperature. Sections were incubated with anti-CD34 primary antibody (Abcam ab81289) at a 1:200 concentration for 1 hour at room temperature and then overnight at 4 degrees. Sections were washed with PBS before incubation with secondary antibody for 1 hour at room temperature. Sections were washed again and imaged using a Zeiss AxioScan Z.1 slide scanner (Carl Zeiss, Jena, Germany).

Skeletal Muscle SDH Activity Stain

Transverse sections (5 μm) were cut from the mid-body of the gastrocnemius muscle. The sections were dried at room temperature for 30 minutes before incubation in a solution made up of 0.2 M phosphate buffer (pH 7.4), 0.1 M MgCl₂, 0.2 M succinic acid (Sigma Chemical Company, St. Louis, MO, USA) and 2.4 mM nitroblue tetrazolium (NBT, Sigma) at 37°C in a humidity chamber for 45 minutes. The sections were then washed in deionized water for three minutes, dehydrated in 50% ethanol for two minutes, and mounted for viewing with DPX mount medium (Electron Microscopy Sciences, Hatfield, PA, USA). Digital photographs were taken from each section at 10X magnification under a Nikon Eclipse TE 2000-U microscope (Nikon, Melville, NY, USA) with a Nikon digital camera (Digital Sight DS-Fi1), and fibers were quantified with imaging software (Image J, NIH).

Skeletal Muscle Citrate Synthase Activity

Whole gastrocnemius muscle was removed from the animals and snap-frozen in liquid nitrogen. The tissue was homogenized in CelLytic MT Cell Lysis Reagent (Millipore-Sigma) with protease inhibitors (Sigma). Citrate Synthase activity was determined using a Citrate Synthase Assay Kit (Millipore-Sigma), in accordance with manufacturer’s instructions.

RNAseq Sample Preparation

28-week-old WT and SPAG7 KO animals were sacrificed and gastrocnemius muscle was snap-frozen in liquid nitrogen. Tissue was pulverized and total RNA was extracted using Qiagen Rneasy Plus Universal mini kit following manufacturer’s instructions (Qiagen, Hilden, Germany). The RNA samples were quantified using Qubit 2.0 Fluorometer (ThermoFisher Scientific, Waltham, MA, USA) and RNA integrity was checked using TapeStation (Agilent Technologies, Palo Alto, CA, USA). RNA sequencing libraries were prepared using the NEBNext Ultra Directional RNA Library Prep Kit for Illumina following manufacturer’s instructions (NEB, Ipswich, MA, USA). Briefly, mRNAs were first enriched with Oligo(dT) beads. Enriched mRNAs were fragmented for 15 minutes at 94 °C. First strand and second strand cDNAs were subsequently synthesized. cDNA fragments were end repaired and adenylated at 3’ends, and universal adapters were ligated to cDNA fragments, followed by index addition and library enrichment by limited-cycle PCR. The sequencing libraries were validated on the Agilent TapeStation (Agilent Technologies, Palo Alto, CA, USA), and quantified by using Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA) as well as by quantitative PCR (KAPA Biosystems, Wilmington, MA, USA).

RNA Sequencing

The sequencing libraries were multiplexed and clustered onto a Illumina flowcell. After clustering, the flow cell was loaded onto the Illumina NovaSeq 6000 instrument according to manufacturer’s instructions. The samples were sequenced using a 2×150bp Paired End (PE) configuration and targeting 30 million reads/sample. Image analysis and base calling were conducted by the Illumina Control Software (HCS). Raw sequence data (.bcl files) generated from Illumina were converted into fastq files and de-multiplexed using Illumina bcl2fastq 2.20 software. One mis-match was allowed for index sequence identification.

Transcriptomic Analysis

The raw FASTQ files were processed using an internal pipeline, Rnaseq Experiment Dashboard, utilizing STAR (v2.7.3a) (Dobin et al., 2013) and Salmon (v1.5.2) (Patro, Duggal, Love, Irizarry, & Kingsford, 2017) to map reads to the mouse genome and transcriptome (GRCm39, Ensembl release 106), respectively, with ERCC spike-in sequences included. Quality of the samples was assessed using FastQC (v0.11.9), Picard Tools CollectRnaSeqMetrics (v 2.26.0), and MultiQC (v1.11) (Ewels, Magnusson, Lundin, & Kaller, 2016) to explore base quality, rRNA content, intronic percentages, and mapping rates. The Salmon transcript abundance estimates were aggregated to the gene level using tximport (v1.22.0) (Soneson, Love, & Robinson, 2015). Additional quality control was performed on the gene-level counts using principal component analysis. Two samples were removed from the analysis due to a batch effect and a label mismatch. Differential expression analysis was performed with DESeq2 (v1.34.0) (Love, Huber, & Anders, 2014), while over-representation analysis and gene set enrichment analysis (GSEA) of the differential expression results was executed using clusterProfiler (v4.2.2) (Wu et al., 2021) with GO (Ashburner et al., 2000; Gene Ontology, 2021), KEGG (Kanehisa, 2019; Kanehisa, Furumichi, Sato, Kawashima, & Ishiguro-Watanabe, 2023; Kanehisa & Goto, 2000), and MsigDB (Liberzon et al., 2015; Liberzon et al., 2011; Subramanian et al., 2005) gene sets. GSEA was performed on the apeglm (v1.16.0) (Zhu, Ibrahim, & Love, 2019)-shrunken log2 fold changes output from DESeq2.

iSPAG7 KO Study Design

Cre/ERT2 heterozygous, SPAG7-floxed homozygous animal were evaluated alongside Cre/ERT2 WT, SPAG7-floxed homozygous littermates. At 7 weeks of age female animals were injected IP with 200mg/kg of tamoxifen (Millipore Sigma, Burlington, MA) dissolved in corn oil for two consecutive days.

Placenta Histological Analysis

Placenta from Het x Het breedings were collected at e18.5 and fixed in 10% formalin. Tissues were embedded in paraffin and cut into 5 um sections. A section from each animal was stained with hematoxylin and eosin using an automatic slide stainer (Leica) and microscopic evaluation was performed by a veterinary pathologist. The junctional zone and whole placenta area, excluding residual decidua, were measured using FIJI software (2.9.0) blinded to the experimental groups. Labyrinth area was calculated as the whole placenta area – junctional zone area. The percentage of junctional zone present was calculated as the junctional zone area / whole placenta area (excluding decidua) x 100. The labyrinth:junctional zone ratio was calculated as the labyrinth area / junctional zone area.

mtDNA Copy Number

Gastrocnemius muscle tissue was dissected from WT and SPAG7 KO animals. Tissues were pulverized and DNA was extracted using phenol:chloroform:isoamyl alcohol (25:24:1) followed by ethanol precipitation. mtDNA copy number was determined via RT-qPCR using a probe for Mitochondrial Cytochrome b (Fisher Scientific - Mm04225271_g1). Nuclear DNA copy number was determined using a probe for Ubiquitin C (Fisher Scientific - Mm02525934_g1). Data are expressed as a ratio of Mitochondrial Cytochrome b expression per Ubiquitin c expression.

Statistical analysis

Statistical analysis, including post-hoc tests and statistical significance, were determined using Welch’s two-way two sample t-test, one-way ANOVA, and Tukey HSD. Exact sample size (n) is included in figure legends. Data are expressed as mean +/- standard deviation. Data were estimated to be statistically significant when P < 0.05.

Data availability

The data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus (Edgar, Domrachev, & Lash, 2002) and are accessible through GEO Series accession number GSE246428 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE246428)