- Research article
- Open access
- Published:
Unraveling the impact of hyperleptinemia on female reproduction: insights from transgenic pig model
Biological Research volume 57, Article number: 60 (2024)
Abstract
Background
Infertility is a growing global health concern affecting millions of couples worldwide. Among several factors, an extreme body weight adversely affects reproductive functions. Leptin is a well-known adipokine that serves as an endocrine signal between adiposity and fertility. However, the exact mechanisms underlying the effects of high leptin level on female reproduction remain unclear.
Methods
Transgenic pigs overexpressing leptin (♀) were produced by backcrossing and screened for leptin overexpression. The growth curve, fat deposition, reproductive performance, apoptosis, serum hormones and cholesterol production, RNA sequencing, and single-nucleus RNA sequencing (snRNA-seq) of the leptin-overexpressing pigs and wild-type group were evaluated.
Results
Transgenic pigs overexpressing leptin (♀) were obtained, which exhibited significantly reduced body weight, body size, and back fat thickness. These pigs manifested a late onset of puberty (330 ± 54.3 vs. 155 ± 14.7 days), irregular estrous behavior characterized by increased inter-estrous interval (29.2 ± 0 vs. 21.3 ± 0.7 days), and more number of matings until pregnancy (at least 3 times). This reproductive impairment in leptin pigs was related to hormonal imbalances characterized by increased levels of FSH, LH, prolactin, E2, P4, and TSH, altered steroidogenesis such as increased levels of serum cholesterol esters along with steroidogenic markers (StAR, CYP19A), and ovarian dysfunctions manifested by neutrophilic infiltration and low expression of caspase-3 positive cells in the ovaries. Moreover, bulk RNA sequencing of the ovaries also revealed neutrophilic infiltration followed by upregulation of inflammation-related genes. Furthermore, snRNA-seq reflected that leptin overexpression triggered immune response, suppressed follicle development and luteinization, resulting in metabolic dysfunction and hormone imbalance in the ovary.
Conclusions
Low body weight in leptin overexpressing pigs adversely affects the reproductive performance, causing delayed puberty, irregular estrous cycles, and reduced breeding efficiency. This is linked to metabolic imbalances, an increased immune response, and altered ovarian functions. This study provides a theoretical basis for the complex mechanisms underlying leptin, and infertility by employing leptin-overexpressing female pigs.
Background
Infertility is a complex disorder that affects approximately 10–15% of couples worldwide. Among several factors, optimum body weight is necessary for initiation and maintenance of reproductive cyclicity [1]. Extreme body weight adversely affects the reproductive functions starting from pubertal development, and menstrual cyclicity to pregnancy [2]. Underweight women exhibit delayed onset of puberty followed by disruption of the menstrual cycle that may lead to prolonged periods of amenorrhea [3,4,5]. Further, long-lasting underweight also disturbs the gonadotropin secretions that ultimately affects the female copulatory behavior leading to infertility [6]. Likewise, overweight females express a higher incidence of menstrual dysfunction, anovulation, and infertility [7]. This compromised reproduction at both extremes is due to impaired secretions of gonadotropins (FSH, LH) that in turn affects the production of steroids and disrupts the ovarian function resulting in infertility [8]. Thus, an adequate mass of adipose tissue is required for the onset of puberty and maintenance of fertility in females.
The adipose tissue is regarded as a passive organ for lipid storage. It provides free fatty acids in response to energy demands and works as an endocrine organ that releases several hormones and cytokines [9]. Leptin is an adipocyte-derived hormone that was initially known as a satiety factor controlling feed intake and energy metabolism but later, it was determined to be an essential regulator of several physiological processes including reproduction. In the context of reproduction, leptin sequentially induces the gonadotropin-releasing hormone (GnRH) expression; peripheral luteinizing hormone (LH) and follicle stimulating hormone (FSH) are stimulated to potentiate the follicular growth and steroid production [10]. Shreds of evidence from murine models has indicated that leptin deficiency leads to infertility, whereas exogenous administration of recombinant leptin restores fertility. Furthermore, transient elevation of plasma leptin concentrations before adolescence is related to change in neuroendocrine function, and it is postulated that the “leptin surge” in childhood may signal the initiation of puberty [11]. In brief, leptin serves as an endocrine signal between the degree of adiposity and fertility, providing a peripheral message to the central nervous system (CNS) on the adequacy of nutritional status for reproductive functions.
Gene-edited (GE) mouse models have contributed prodigiously to elucidate the role of leptin in reproduction and it was reported that GE mice overexpressing leptin were skinny with accelerated puberty and intact fertility at younger ages followed by successful delivery of healthy pups. However, at older ages, they develop hypothalamic hypogonadism characterized by prolonged menstrual cycles, atrophic ovaries, reduced hypothalamic GnRH contents, and poor pituitary luteinizing hormone secretion [7]. Nevertheless, the small size, short life, and other differences between humans and rodents limit their ability to model complex diseases like infertility. Contrarily, pigs being large in size and having anatomical, physiological, and genetic similarities are regarded as ideal models for humans. A large number of pig models for xenotransplantation, cancer, diabetes, cardiovascular, neurodegenerative, and immunological disorders have been developed in recent years [12]. Pig models have been utilized for evaluating reproductive disorders like polycystic syndrome, ovarian cancer, and preeclampsia [13]. Thus, in this study, we utilized a pig model to evaluate the effect of high leptin on female reproduction.
Indeed, transcriptome studies using bulk RNA samples have identified expression changes in expression in pig ovaries, but interpretation is limited due to the heterogeneity of the organ and the varying effects of leptin on different cell types, which are not captured by bulk sequencing or qPCR. However, recent advances in single-nucleus RNA sequencing (snRNA-seq) have enabled the transcriptomic analysis of large numbers of cells at single-cell resolution with high depth [14]. This method provides insights into the gene expression profiles of individual ovarian cells, allowing for a better understanding of the heterogeneity of the ovarian tissue and the specific effects of leptin on different cell types [15,16,17]. On the basis of leptin-overexpressing pigs previously generated in our laboratory [18], in this study; the growth curve, fat deposition, reproductive performance, serum hormones and cholesterol profile, steroidogenesis and apoptosis, and changes at individual ovarian cells of the leptin overexpressing pigs and wild-type (WT) group were evaluated systemically to further elucidate the effect of leptin on female reproduction and its mechanism, and to provide a certain theoretical basis for this area of research.
Methods
Experimental animals
Animals used in this study were regularly maintained in the Laboratory Animal Centre of Yunnan Agricultural University, Kunming, China. All animal experiments were performed with the approval of the Animal Care and Use Committee of the above institution. The female pigs overexpressing leptin were obtained by backcrossing the WT female with already generated leptin overexpressing transgenic boar [18] and offspring were screened for the presence of transgenic vector containing leptin and enhanced green fluorescent protein (EGFP) by polymerase chain reaction (PCR) using primers listed in additional file Table S1, and circulating serum leptin concentrations. Furthermore, the information regarding animal used in this study is provided in additional file Table S2.
DNA extraction and polymerase chain reaction (PCR)
The genomic DNA from blood of each newborn piglet was extracted using TIANamp genomic DNA kit (Tiangen, Beijing, China) and PCR was carried out using 2×Ex Taq Master Mix (TaKaRa CW0718). The PCR reaction was composed of 0.5 µg of template, 10 µl of Ex Taq Master Mix, 0.4 µM of each primer and distilled water to a final volume of 20 µl. The amplification of each primer pairs was performed in separate reaction with 35 cycles of 94 °C for 30 s, 55–60 °C for 30 s, and 72 °C for 30 s. PCR products were separated by 1% agarose gel electrophoresis and visualized by staining with DL 2000 DNA Marker (TSINGKE). Depending upon the presence of transgenic vector comprised of leptin gene with EGFP expression, the offspring were divided into two groups; Leptin pigs and WT pigs. Next, the circulating serum leptin concentrations were measured using Porcine Leptin Elisa Kit (SEKP-0278, Solarbio life sciences, China).
Body weight and size measurements
Body weight (BW) and size including length, height, heart girth, chest depth, chest width, and abdominal circumference of leptin (n = 5) and WT (n = 4) pigs were measured at 30d intervals from d180 ~ d390 of age. The live weight of each animal was determined by suspending the animal on an electronic scale balance and the weight of each animal was taken and recorded. The body size measurements were recorded. In brief, body length was measured using the tape rule as the distance from the occipital protuberance to the base of the tail and heart girth was determined by measuring the chest with a tape. Height was measured as a distance from the surface of the ground to the withers using a meter rule and chest depth was measured as the distance from the sternum at the region just in front of the forelimb to the withers, while chest width was recorded as distance between shoulders using meter ruler. Abdominal circumference was measured from the bottom of flank on one side to the bottom of flank on the other side of pig.
Back fat thickness
The back fat thickness of WT and leptin pigs was scanned ultrasonically. Briefly, pigs were immobilized and restrained by the head in a squeeze chute and the image site at the P2 position (at the level of head of last left rib) was determined by physical palpation. The pigs were held manually, avoiding any abnormal situation that could stress the animal, and were only scanned in a relaxed posture, permitting accurate measurements. A mix of Eco Gel and isopropyl alcohol was used as a sound-conducting material to allow a better acoustic contact surface between the probe and the skin. An Aplio 500 V real-time ultrasound machine (BSM34-0927, Toshiba, Japan) equipped with a 5.0–14.0 MHz linear array transducer (PLT-1005BT, Toshiba, Japan) was used for scanning the image site.
Reproductive performance
Reproductive parameters comprising of onset of puberty, pattern of estrous expression, and breeding performance were evaluated in leptin and WT pigs. The estrus observation in all pigs was started at d120 of age and gilts were visually observed two times (9:00 AM and 05:00 PM) a day for the expression of estrus signs in the presence of boar and stages of estrous i.e. proestrus, estrus, met-estrus, and di-estrus were ascertained according to previously described estrus scoring system [19] with slight modifications. The age at which gilts displayed the first estrus was recorded as the age of puberty and the pattern of estrous expression was recorded until 540 d. The inter-estrous interval was defined as the duration (days) from the start of one estrous to the start of the next estrous. Next, we evaluated the breeding performance of leptin and WT pigs by mating with a fertile male. Briefly, leptin females (n = 4) were used for mating, and gilts were mated 2 times in one estrous cycle using the same male. Usually, the first mating was performed 24 h after standing reflex followed by the second mating 18 h later. Pregnancy was monitored by non-return to estrous and confirmed by ultrasound scanner for around 33 d after service (HS-101 V, Honda Electronics Co., Ltd., Yamazuka, Japan).
Serum hormonal and lipid profile
Blood samples were collected at five different stages of the estrous cycle: proestrus, estrus, met-estrus, d6 post-met-estrus, and d12 post-met-estrus for hormonal analysis. Briefly, the pigs were fasted overnight and 5 ~ 10 mL blood was collected from the jugular vein using a vacutainer (Cat no. 367986, BD Vacutainer SST™, USA) and transported to the laboratory within 10 min. The serum was separated by centrifugation at 4˚C at 3000 rpm for 30 min and stored at -80˚C until assayed for hormonal analysis and cholesterol esters (CE) profile. FSH, LH, thyroid stimulating hormones (TSH), porcine prolactin (PRL), anti-müllarian hormone (AMH), and testosterone were measured by enzyme-linked immunosorbent assay (Shanghai Enzyme-linked Biotechnology Co., Ltd., China). The serum progesterone (P4) and estradiol (E2) concentrations were measured using a radioimmunoassay (Iodine [125I] Progesterone Radioimmunoassay Kit (B08B), and Iodine [125I] Estradiol Radioimmunoassay Kit (B05B), Beijing North Institute of Biotechnology, Co., Ltd. (Beijing, China). The serum leptin concentrations were measured using Porcine Leptin Elisa Kit (SEKP-0278, Solarbio life sciences, China).
Serum samples lipid extraction, separation of cholesterol esters and phospholipids via thin-layer chromatography, and determination of fatty acids via gas chromatography were conducted as described previously [20]. The fatty acid composition was analyzed with An Agilent 7890 A gas chromatographer equipped with 15 m × 0.25 mm × 0.25 μm DB-WAX columns (Agilent). While, C15:0 was used as a standard for quantitation. All the data were first normalized with the phospholipids profile and then presented as mean ± SD.
Luteal cell culture
Fresh ovaries of leptin and WT pigs were removed by laparotomy and immersed immediately in sterile phosphate buffer saline (PBS) with 5% penicillin-streptomycin solution and transported to laboratory. The luteal tissues were obtained from the ovaries and cut into small pieces using sterilized ophthalmic scissors, minced well and dispersed by pipetting in M199 (Sigma, M4530) medium containing collagenase (10 mg/mL) and the tissue was digested on a horizontal shaker in an incubator at 37 °C for 30 min. The supernatants containing luteal cells were then decanted through nylon mesh to remove debris. The cells were washed twice by centrifugation and resuspended in 10 ml Dulbecco’s modified Eagle’s medium with 10% fetal bovine serum and incubated at 37 °C, 5% CO2 in a 96-well plate for 24 h and then cells were collected and used for subsequent experiment.
Quantitative polymerase chain reaction (qPCR)
Total RNA from ovarian tissues and luteal cells were isolated using the Trizol reagent (Transgen Up, China) according to the manufacture’s instruction. cDNA was synthesized from the total RNA using a PrimeScript RT reagent Kit (TaKaRa, Japan) and was used as a template to perform qPCR in SYRB green-based qPCR instrument (CFX-96, Bio-Rad, USA). The reaction was performed in a 20 µL reaction mixtures comprising 10 µL of 2× SYBR (TaKaRa, Japan), 1 µL of cDNA, 1 µL of forward primer, 1 µL of reverse primer, and 7 µL of ddH2O. The reaction program as following: 95 °C for 30 s, followed by 40 cycles of 95 °C for 10 s and 62 °C for 45 s. The relative expression levels of target genes were quantified by 2-ΔΔCT. The primers are listed in additional file Table S3.
Ovarian histology and immunohistochemistry
For histological examination, ovaries were fixed in 4% paraformaldehyde for 48 ~ 72 h, processed by an automatic tissue processor (Yd-12p, Jinhua Yidi, medical appliance Co., Ltd, Jinhua, China) and embedded in a paraffin block (Yd-6D, Jinhua Yidi, medical appliance Co., Ltd, Jinhua, China). The paraffin blocks were cut into 5-um-thick sections using a Microm HM 325 microtome (Thermo Scientific, Waltham, MA, USA) and allowed to dry on glass slides overnight at 37 °C. Thereafter, the tissue sections were deparaffinized in xylene and rehydrated through graded ethanol dilutions. Sections were stained with hematoxylin–eosin (H&E) (G1120; Solarbio, China) according to manufacturer’s instruction.
For immunohistochemistry (IHC), after dewaxing and hydration, sections were incubated in 3% H2O2 solution for 30 min, and washed with PBS for thrice (each time 3 min). After that, the sections were blocked in PBS containing 5% BSA for 15 min at room temperature. Finally, the tissue sections were incubated with caspase-3 antibody (additional file Table S4) at 4 °C overnight. After washing with PBS for thrice, sections were incubated with specific secondary antibodies (R&D, USA) for 20 min. After washing thrice again (each time 3 min), sections were stained with fresh DAB (KIT-9901, Elivision TM plus Polyer HRP IHC Kit, Fuzhou, China) solution in dark for 5 min and then washed with PBS thrice (each time 3 min). Hematoxylin counter-staining, and neutral gum sealing of the slides was performed. Finally, the slides were imaged using OLYMPUS BX53 fluorescence microscope and analyzed using software of accessories.
Protein extraction and immunoblotting
The ovarian tissues from WT and leptin pigs were used to evaluate the expression of different protein levels (listed in additional file Table S4) using western blotting. In brief, the ovarian tissues were lysed in RIPA lysis buffer (Bestbio, China) with protease inhibitors at 4 °C. After lysis, supernatants were obtained by centrifugation at 13,000 × g for 15 min at 4 °C. Equal amounts of protein (100 µg) were run on SDS-PAGE gel, along with molecular weight marker. After electrophoresis, the proteins were transferred to PVDF membranes and reacted with primary antibodies against various antibodies and β-actin at 4 °C overnight. After incubation, the membranes were washed and incubated with specific secondary antibodies (R&D, USA). The membranes were then incubated with ECL (Easysee Western Blot Kit, China) and visualized with an Imaging System (Bio-Rad, Universal Hood II, USA).
RNA−seq analysis and quality assessment
Total RNA was extracted using TRIzol Reagent (Invitrogen, CA, United States) and purified using an RNeasy Mini Kit (Qiagen, CA, United States). The quality of RNA was assessed using an Agilent 2100 Bioanalyzer (Agilent, Palo Alto, CA, United States). The rRNA-depleted RNA samples were further processed in accordance with the Illumina protocol (New England Biolabs, Massachusetts, United States). After cDNA synthesis, the samples were sequenced with an Illumina Novaseq6000 by Gene Denovo Biotechnology Co. (Guangzhou, China). The raw data were recorded. The overall quality of the RNA−seq data was evaluated by fastp (version 0.18.0). Clean reads were aligned to the Sus scrofa reference version 11.1 using HISAT2. 2.4 [21] with the default parameters.
Screening and clustering analysis of differentially expressed genes
Data preprocessing and follow-up analysis were performed DESeq2 software. The lists of DEGs between leptin and WT pigs were generated using the edgeR package (version 3.32.0) (Robinson et al., 2010). To normalize the raw data, the genes with the parameter of false discovery rate (FDR) below 0.05 and absolute fold change ≥ 2 were considered differentially expressed genes (DEGs) between the leptin and WT groups. Hierarchical clustering analysis was performed based on the expression levels of all transcripts and significantly differentially expressed transcripts using the pheatmap R package based on Euclidean distance.
Functional enrichment analysis
The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis in each module and network was conducted using the Database for Annotation, Visualization and Integrated Discover (DAVID). DEGs and enriched pathways were mapped using KEGG pathway annotation with KOBAS3.02. The top 5 KEGG pathways were selected and ranked by the enrichment factor. Subsequently, the Venn diagram tool was used to identify the common genes that were the focus of this work. Finally, gene set enrichment analysis (GSEA) using software GSEA and MSigDB [22] was performed to identify a set of genes in specific GO terms.
Single nucleus transcriptome sequencing (snRNA-Seq) libraries construction
Cellular suspensions were loaded on a 10X Genomics GemCode Single-cell instrument that generates single-cell Gel Bead-In-EMlusion (GEMs). Libraries were generated and sequenced from the cDNAs with Chromium Next GEM Single Cell 3’ Reagent Kits v2. following the manufacturer’s instructions.
Single-cell RNA-seq data processing
FASTQ files of 2 samples were processed with the use of Cell Ranger (v.3.1.0) count pipeline coupled with Sus scrofa reference version 11.1 to generate feature-barcode matrices. Seurat object list was then generated by Seurat package (v.4.0.4) [23], with R software (v.4.1.0) following these criteria: (1) min.cells = 5; (2) 200 < nFeature_RNA < 3600; (3) percent.mt < 0.1. In other words, genes expressed in at least 5 cells and gene number detected in cells ranging from 800 to 10,000 were kept for further analysis, and low-quality cells were filtered if R 10% UMIs derived from the mitochondrial genome.
Dimensionality reduction and clustering
To remove batch effects across samples, canonical correlation analysis method was used for data integration [24]. In detail, we initially normalized the filtered gene expression data with Normalize Data function with parameter (normalization.method = “LogNormalize”, scale.factor = 10000). Then top 2,000 variable genes were identified using the ‘vst’ method with FindVariableFeatures function in each sample, respectively. At last, we used FindIntegrationAnchors with parameters (i.e. anchor.features = 2000, k.filter = 200, dims = 1:30) and IntegrateData functions with the top 30 dimensions among 2 samples. After filtering low quality cells with 200 < nFeature_RNA < 3600, percent.HB < 5 and nCount_RNA < 10,000, we ran ScaleData function. To perform dimensionality reduction, the RunPCA function was conducted on linear-transformation scaled data with 2000 variable features and we performed UMAP with the top 30 dimensions. Finally, we clustered cells by using the FindNeighbors and FindClusters (resolution = 0.5) functions, which got 16 clusters.
Cell-type annotation and cluster marker identification
Cell clusters were obtained, after dimensional reduction and projection of all cells into two-dimensional space by UMAP. The Seurat FindAllMarkers function was used to identify markers for each cell cluster with the default settings. Canonical markers of specific cell types were used for cluster annotation. Ovarian cell-type markers were selected from PanglaoDB database(https://panglaodb.se/), CellMarker database(http://biocc.hrbmu.edu.cn/CellMarker/) and literature review.
DEG identification and functional pathways enrichment
The FindMarkers function in Seurat was used to identify differential expressed genes (DEGs) between two groups of cells with default parameters (logfc.threshold = 0.25, test.use = ‘‘wilcox,’’ min.pct = 0.1), and the enrich GO function in the clusterProfiler package (v.3.15.2) [25], was used to perform functional analysis with differential gene sets annotated with GO database (http://geneontology.org/). Enrichment pathways were obtained with parameters (pvalueCutoff = 0.05, pAdjustMethod = “BH”, OrgDb = “org.Ss.eg.db”).
Cell trajectory construction by Monocle
The subpopulation of germ line or granulosa cell (GC) line was imported into Monocle (version 2.20.0) to dissect the cell differentiation fate, also termed “pseudotime analysis.” With the gene count matrix as input, the new dataset for Monocle object was created, and functions of “reduceDimension” and “orderCells” were carried out to generate the cell trajectory based on pseudotime. Particularly, the ordering genes were DEGs between clusters in each cell type calculated by “differentialGeneTest” function in Monocle. In addition, the root state (that is, a prebranch in the heatmap) was set and adjusted following consideration of the biological meanings of different cell branches.
The regulon activity of transcription factors with SCENIC
The SCENIC algorithm was utilized to analyze the activity of transcription factors (TFs) and identify regulons (TFs and their target genes) in individual cells. The gene expression matrix, with genes in rows and cells in columns, was input into SCENIC (version 0.9.1) [26]. Co-expressed genes for each TF were constructed using GENIE3 software, followed by Spearman’s correlation analysis between the TFs and their potential targets. The “runSCENIC” procedure was then employed to generate the Gene Regulatory Networks (GRNs), also known as regulons. Regulon activity was analyzed using the AUCell (Area Under the Curve) software, applying a default threshold to categorize specific regulons as “0” (off) or “1” (on). t-SNE parameters were set to visualize the data with 50 principal components and a perplexity of 50, determined through a consistency test across multiple perplexity values and number of PCs. Cell states were mapped using specific regulons, and the average binary regulon activity was calculated.
Cell–cell communication analysis
Cell–cell communication was analyzed using iTALK package. The input data were Seurat object, and the ligand–receptor pairs were detected from top 50% highly expressed genes. The communication types mainly included growth factors, cytokines, and checkpoint. The network plot was visualized using LRPlot function.
Statistical analysis
Analyses were performed in the SAS statistical software (SAS ver. 9.4, SAS/STAT, SAS Institute Inc., Cary, NC). Data were expressed as least square mean ± standard error of mean unless otherwise stated. Continuous variables and measured over time, such as body weight and size including length, height, heart girth, chest depth, chest width and abdominal circumference, and hormones were analyzed by using the MIXED procedure of SAS (SAS/STAT) with estimations carried out by the method of restricted estimated maximum likelihood. All the mixed models included the repeated statement as day or stage of the estrous cycle and random effect of the animal identification nested within treatment for the proper error term. Continuous variable but measured on single point, such as serum leptin levels, pubertal onset, inter-estrus interval, back fat thickness, serum circulating cholesterol esters profile, relative expression of mRNA and protein were analyzed by PROC TTEST in SAS. Statistical significance was defined as *P < 0.05, **P < 0.01.
Results
Screening of leptin pigs and growth performance
To evaluate the effect of leptin overexpression on female reproduction, leptin overexpressing transgenic sows were generated by breeding the WT sows with already generated leptin overexpressed male. In brief, leptin overexpressed male was re-cloned and the cloned offspring (♂) were mated with Duroc/Landrace/Yorkshire sow (DLY, ♀) to obtain the F1 generation, and then the F1 generation was further crossed to obtain the F2 generation and so on (Fig. 1. A). The overexpression of leptin in offspring was confirmed by PCR and circulating serum levels of leptin (Fig. 1. B and C). The offspring-harboring transgenic vector (leptin gene with EGFP) were designated as leptin pigs, while, the others lacking transgenic vector were defined as WT (Fig. 1. B). The leptin expression in leptin pigs was further elucidated by high levels of serum leptin (3045.9 ± 794.2 vs. 521.1 ± 98.0 pg/mL, P < 0.05; Fig. 1. C). To assess the growth characteristics of leptin overexpressing pigs, we documented the growth parameters, including body weight, length, height, heart girth, chest width, chest depth and abdominal circumference (Fig. 1. D-J). The results showed that leptin pigs displayed slow growth rate as their body weight was significantly lower than WT pigs as at 180 d of age (30.0 ± 2.4 vs. 101.8 ± 3.8 kg, P < 0.05, Fig. 1. D). Nevertheless, the leptin pigs gradually gained weight but they remained underweight as compared to WT pigs (P < 0.05). Similarly, the other growth parameters like length, height, heart girth, chest width, chest depth, and abdominal circumference were also significantly lower in leptin pigs as compared to WT (P < 0.05; Fig. 1. E-J). Additionally, the back fat thickness in leptin pigs was significantly lower than WT (5.6 ± 6.4 vs. 42.9 ± 1.9 mm, P < 0.05; Fig. 1. K). Collectively, these findings indicated that leptin pigs exhibited stunted growth with decreased back fat thickness.
Screening of leptin pigs and body development: (A) Schematic diagram of leptin pigs generation and breeding. (B-C) Screening of leptin pigs (B) Leptin and EGFP expression was confirmed by 1% AGE. Lane M indicates DNA marker (DL2000). β-actin was used as an internal control (C) Serum leptin concentration in WT (n = 5) and leptin pigs (n = 7). (D-J) Changes in body weight and size comprising of length, height, heart girth, chest depth, chest width and abdominal circumference between WT (n = 4) and leptin pigs (n = 5). (K) Measurement of back fat thickness by ultrasonography in WT and leptin pigs. The back fat thickness in leptin pigs was significantly low (P < 0.05), WT (n = 3) and leptin pigs (n = 3)
Reproductive performance of leptin pigs
Leptin has been involved in regulation of female reproduction and we evaluated the reproductive performance of leptin pigs comprising of pubertal onset, estrus expression pattern, and breeding performance (Fig. 2). It was observed that the onset of estrus in leptin pigs (n = 8) was at the age of 294d, 416d, 259d, 407d, 335d, 352d, 290d, and 287d respectively. Collectively, the leptin pigs showed delayed onset of puberty as compared to WT pigs (330.0 ± 54.3 vs. 155 ± 14.7 days, P < 0.05; Fig. 2. A). Next, we evaluated the pattern of estrus expression in both leptin and WT pigs and observed that leptin pigs manifested irregular estrous expression characterized by both long and short cycles as sometimes they displayed estrus as short as after 9 d and sometimes as long as after 63 d (Fig. 2. B). The average length of estrous cycle in leptin pigs (n = 8) was 35.8d, 40d, 28.9d, 33d, 23.4d, 23d, 24.5d and 25d respectively, while in WT pigs (n = 5) it was 21.4d, 22.4d, 21.3d, 20.6d and 21d respectively. Collectively, the inter-estrous interval in leptin pigs was significantly higher than WT (29.2 ± 6.0 vs. 21.3 ± 0.7 days, P < 0.05; Fig. 2. C). Meanwhile, we started the breeding of leptin pigs (n = 4) and it was observed that leptin pigs received more mating until they got pregnant. As two sows were bred for 3 times and at 3rd time breeding, they became pregnant (Fig. 2. D & E). While, the 3rd leptin sow was bred for 6 times (Fig. 2. F), and 4th sow was bred for 5 times (Fig. 2. G), but both heads remained non-pregnant (Fig. 2. F & G). Collectively, these data indicated that leptin pigs exhibited a greater number of mating (17 times) with a minimum of 3 times to a maximum of 6 times mating. Contrarily, WT pigs (n = 10) were mated for a total of 10 times, out of that, they exhibited successful pregnancy for 8 times, and remained non-pregnant for 2 times. Whilst, the leptin pigs were mated for a total of 17 times, out of which they remained non-pregnant for 15 times, while became pregnant for 2 times (Fig. 2. H). Thus, delayed onset of puberty and irregular estrous expression was followed by breeding insufficiency as a greater number of matings in leptin pigs were required until successful pregnancy.
Puberty, estrus expression, and breeding performance of leptin pigs. (A) Age of puberty in WT (n = 5) and leptin pigs (n = 8). (B) Pattern of estrus expression in WT (n = 5) and leptin pigs (n = 8). Each peak represents the estrous expression. (C). Comparison of inter-estrous interval between WT (n = 5) and leptin pigs (n = 8). (D) Breeding performance of individual leptin pigs (n = 4). Two heads (D, E) get pregnant at 3rd time of breeding with male, while other 2 heads (E, F) remain non-pregnant even after 5–6 times breeding. (H) Comparison of breeding efficiency between leptin (n = 4) and WT (n = 10) pigs
Hormonal and lipid profiles of leptin pigs and steroidogenesis
Gonadotropins and steroids regulate the reproductive axis and therefore, we evaluated the serum hormonal profile in leptin pigs at five different stages of the estrous cycle. Results showed that the concentrations of FSH, E2, TSH, and PRL in leptin pigs were higher (P < 0.05) as compared with WT regardless of the estrous cycle stage (Fig. 3. A & B, E & F). Furthermore, the concentrations of LH and P4 were also significantly higher in leptin pigs but interaction (P < 0.001) between treatment and the estrous stage was observed (Fig. 3. C & D) that indicated that the concentration of LH in WT pigs almost remained consistent across all stages of the estrous cycle, however, in leptin pigs it increased in early stages of the estrous cycle and then declined at the end of estrous cycle (d12 post met-estrus). Despite this, the concentrations of P4 were similar at earlier stages of the estrous cycle (proestrus, estrus, and met-estrus) in both leptin and WT pigs but it was significantly higher in leptin pigs than WT during d6 and d12 after met-estrus (P < 0.05; Fig. 3. D). However, the concentrations of testosterone and AMH hormones remain unchanged in leptin pigs (P > 0.05; Fig. 3. G & H). Collectively, these results indicated that the irregular estrous in leptin pigs followed an abnormal hormonal profile.
Serum hormonal changes in leptin pigs and steroidogenesis: Serum hormonal changes at five different stages of estrous cycle comprising of (A) follicle stimulating hormone (FSH) (B) estradiol (E2) (C) luteinizing hormone (LH) (D) progesterone (P4) (E) thyroid stimulating hormones (TSH) (F) porcine prolactin (PRL) (G) testosterone and (H) anti-müllerian hormone (AMH) levels in leptin (n = 4) and WT (n = 4) pigs. P: proestrus; E: estrus; M: metestrus; d6 post M: day 6 post metestrus and d6 post M: day 12 post metestrus. Leptin OE stands for leptin overexpression. (I) Serum levels of total cholesterol esters (CE) (J) Total saturated fatty acids (SFA) (K) Total monounsaturated fatty acids (MUFA), (L) Total polyunsaturated fatty acids (PUFA) in leptin and WT pigs. (M) Serum levels of α-linolenic acid (α-LA), (N) linolenic acid (LA), (O) γ-linolenic acids (γ-LA) and (P) arachidonic acids (AA) in WT (n = 3) and leptin pigs (n = 3). (Q) mRNA expression of steroidogenic markers at ovary of WT (n = 2) and leptin pigs (n = 2). GAPDH was used as an internal control. (R) mRNA expression of steroidogenic markers at luteal cells. GAPDH was used as an internal control. (S) Protein abundance of steroidogenic markers at ovary of WT (n = 2) and leptin pigs (n = 2). β-actin was used as an internal control. (T) Ovarian histology of WT and leptin pigs. (U) The expression of apoptotic markers at ovary of WT and leptin pigs. The number of cas-3 positive cells were calculated after IHC. *P < 0.05, **P < 0.01. (V) The number of differentially expressed genes in leptin and WT group. (W) Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis performed by Metascape. (X) Heatmap of genes enriched in different pathways
To clarify the essence of body fat composition on reproductive performance, we further assessed the blood cholesterol ester (CE) profile in WT and leptin pigs. It was observed that the total CE levels in the serum of leptin pigs were significantly higher than WT pigs (125.6 ± 17.11 vs. 100.0 ± 12.5, P < 0.05; Fig. 3. I). Levels of total saturated fatty acids (SFA) in leptin pigs remain unchanged (90.8 ± 16 vs. 100.0 ± 16.8, P > 0.05; Fig. 3. J), while the levels of monounsaturated fatty acids (MUFA) and polyunsaturated fatty acids (PUFA) were increased in leptin pigs as compared with WT (136.0 ± 20.1 vs. 100.0 ± 16.0, P < 0.05; Fig. 3. K), and (136.0 ± 19.5 vs. 100.0 ± 14.6, P < 0.05; Fig. 3. L), respectively. Specifically, α-linolenic acid (α-LA) levels in leptin pigs were significantly decreased (42.9 ± 8.9 vs. 100.0 ± 14.2, P < 0.05; Fig. 3. M). While, leptin pigs also exhibited (P < 0.05) high levels of linolenic acid (132.1 ± 23.2 vs. 100.0 ± 15.0), γ-linolenic acid (244.3 ± 80.8 vs. 100.0 ± 17.8), and arachidonic acid (195.5 ± 28.0 vs. 100.0 ± 14.7) than WT pigs (Fig. 3. N-P). Collectively, the levels of CE (precursors of steroids) were increased in leptin pigs.
Next, we evaluated the mRNA and protein expression patterns of steroidogenic markers. Results indicated that mRNA expression of leptin, steroidogenic acute regulatory protein (StAR), cytochrome P450 19A1 (CYP19A1), estrogen receptor 1 (ESR1), and follicle-stimulating hormone receptor (FSHR) were higher (P < 0.05) in leptin pig ovaries, while the mRNA expression of cytochrome P450 17A1 (CYP17A1) and luteinizing hormone/choriogonadotropin receptor (LHCGR) was comparable between leptin and WT pigs (P > 0.05, Fig. 3. Q). Likewise, the mRNA expression of leptin, StAR, CYP17A1, and CYP19A1 was significantly higher (P < 0.05) in leptin luteal cells, while CYP11A1 was comparable (P > 0.05) between leptin and WT luteal cells (Fig. 3. R). Furthermore, western blot analysis indicated that the protein abundance of StAR, CYP19A1, and ESR1 was significantly higher in leptin pigs ovaries (P < 0.05, Fig. 3. S), while CYP17A1, FSHR, and LHCGR abundance was comparable between leptin and WT pigs (P > 0.05, Fig. 3. S).
Although histological analysis revealed no significant difference in the number of early follicles and sinus follicles in leptin pigs relative to WT. However, leptin pigs exhibited more severe vascular inflammatory infiltration of the ovarian corpus luteum relative to WT pigs (Fig. 3. T). Furthermore, the IHC analysis of ovaries revealed that the number of Caspase-3 positive cells was significantly lower in leptin pig ovaries than that of WT group (P < 0.05, Fig. 3. U). Collectively, these results indicated that the enhanced expression of StAR and CYP19A1 at tissue and cell levels was followed by inflammatory changes and compromised apoptosis in leptin pig ovaries.
Transcriptome profile and GSEA of pig ovaries
Our RNA sequencing analysis uncovered significant differences in gene expression between the leptin-overexpressing pigs and WT group, with a total of 312 differentially expressed genes. Among these 97 were upregulated, while 205 were downregulated in leptin pigs (Fig. 3. V, Fig. S1. A). To gain insight into the functional implications of these differentially expressed genes, we conducted functional enrichment analysis, which revealed several significantly enriched processes. These included cytoskeleton organization, antibacterial humoral response, humoral immune response, positive regulation of lipid localization, and others (Fig. 3. W). Furthermore, GSEA revealed negative regulation of interleukin-2 production pathways was downregulated in leptin pigs (Fig. S1. B). To visualize the differentially expressed genes, a heatmap was generated, and some genes crucial for apoptosis (TP73, CLDN4, DTHD1, and KCNRG) and reproduction (APLN, OVGP1, FABP4, and CELSR1) were downregulated in leptin pigs. Whilst, the genes involved in metabolic changes (ACSL4, LEP, and ND4L) and inflammatory changes (PTX-3, TFRC, COX-2 and AQP9) were upregulated in leptin pigs (Fig. 3. X). Collectively, these data revealed that the ovaries of leptin pigs exhibited inflammatory changes and compromised apoptosis and reproduction related markers, accompanied by notable metabolic alterations. These insights contribute to our understanding of the molecular mechanisms underlying the effects of leptin overexpression on ovarian function.
Single-nucleus transcriptome profiling of pig ovaries identified different ovarian cell types and gene expression signatures
To gain a deeper understanding of the impact of high leptin levels on individual ovarian cells, we characterized the distinct cell lineage and gene expression dynamics in both WT and leptin-overexpressed pigs. This was achieved by performing snRNA-seq on the ovarian tissue of both wild-type sus scrofa and our transgenic leptin pig model, which were then analyzed using bioinformatics tools (Fig. S2. A). After stringent cell filtration, the high-quality transcriptomes of 17,108 single cells (8686 in the WT group and 8422 in the leptin group; Fig. S2. B) were retained for subsequent analyses. We performed unsupervised clustering and grouped the total ovarian cells into 10 clusters as visualized in the Uniform Manifold Approximation and Projection (UMAP) plot (Fig. 4. A). By using established marker genes of ovarian cells, we identified 10 clusters as large granulosa cells (LGC), theca-interstitial cells (TIC), luteal cells (LC), endothelial cells (EC), small granulosa cells (SGC), macrophage, lymphocytes (1 and 2), perivascular cells (PCs), and epithelial cells (EC). We found that both the leptin and WT group contained all kinds of cell types, while the cell fraction was different between groups (Fig. 4. B, Fig. S2. C).
Single-nucleus transcriptome profiling of pig ovaries identified different ovarian cell types and gene expression signatures. (A) UMAP plot of the merged dataset of leptin and WT pig ovaries showing 10 ovarian cell types. (B) Cell type fractions in WT and leptin group. (C) Representative GO terms of upregulated DEGs between ovarian cell types. Color in left is corresponding to Fig. 4.A to represent cell types. Gene ratio is showed by color from blue to yellow, Gene count is indicated by size of circle. (D) UMAP plot in four subtypes of small granulosa cells (E). UMAP plot shows WT and leptin cell distribution in SGC’s subtypes. (F) Single-cell trajectories of the 3 SGC states through the pseudotime (a). Single-cell trajectories along with cell cycle phase (b), WT and Leptin group distribution (c) and pseudotime (d). (G) Heatmap of two kinetic trends for state 2 and state 3. Middle line is the start of pseudotime
All clusters were identified according to the specific expression of markers (Fig. S2. E). Clusters corresponding to GCs were identified based on the expression of granulosa cell markers, such as CYP19A1 and FST. LGCs and SGCs were annotated according to their average RNA levels and INHA (an antagonist of granulosa cell proliferation) expression patterns. NR5A2, RNF150 and NLGN1 were specifically expressed at LGCs, while INHA, MYH11 and SERPINE2 were specifically expressed at SGCs. Theca-interstitial cell markers COL3A1 were specifically expressed in one cluster. Interestingly, some of cells in this cluster expressed high levels of the theca marker PDGFRA, implying that this cluster contained theca cells that are associated with the development and maturation of follicles and CDH11 and SLIT2 were specifically distributed at TIC (Fig. S2. E). We next analyzed biological functions for every cell cluster by performing GO analysis of DEGs across all the clusters, revealing unique characteristics of these ovarian cells (Fig. 4. D). For example, GO terms of LGCs were involved in “lipid metabolic process” and “steroid metabolic process”. GO terms including “extracellular matrix organization” and “cell-cell signaling by wnt” were enriched in TIC, suggesting that TICs provide structural support and mediate signaling between ovarian cells. And GO terms including “blood vessel development” and “immune response” were enriched in ECs and macrophages, respectively. Collectively, we verified 10 different ovarian cell types in pig ovaries and depicted gene expression signatures for each cell type.
Developmental patterns of small granulosa cells indicate inflammatory state in leptin pig
Granulosa cells play an important role in ovarian follicle development [27]. To investigate gene expression dynamics of SGCs, we performed unsupervised clustering and identified four subtypes of SGCs (Fig. 4. E). Notably, cluster 1 and 2 exhibit obvious overlap between WT and leptin groups, while cluster 3 and 4 were mainly consisted with cells from leptin group (Fig. 4. E). To further identify the clusters, we detected DEGs in comparisons between four SGC subpopulations and the highly variable expressed genes were visualized in UMAP plots (Fig. S2. A). For example, mitosis-associated genes, such as DIAPH3, HMMR, SGO2, and CENPE were mainly expressed in cluster 2 (green), while inflammation-related genes like SOX5, TOX, AKAP12, and SOD exhibited peak expression in cluster 4 (purple). Subsequently, GO analysis reveals biological variance between SGC subtypes (Fig. S2. B). For instance, GO terms including “androgen metabolic process” and “lipoprotein biosynthetic process” were enriched in cluster 1, genes involving “hormone biosynthetic process” and “steroid biosynthetic process” were enriched in cluster 3, GO terms including “mitotic nuclear division” and “regulation of cell cycle” were enriched in cluster 2, GO terms including “positive regulation of interleukin-2 production” and “cAMP catabolic process” were enriched in cluster 4. Therefore, we referred SGCs in subtypes as Metabolic cells, Mitotic cells, Synthetic cells, and Immunoreactive cells, respectively. In addition, the dominant population in Synthetic cells and Immunoreactive cells were mainly cells from the leptin group, indicating SGCs in leptin-overexpressed ovaries have distinct gene signatures regarding lipid metabolism and inflammation (Fig. 4. E).
To explore how SGC cell transition from one state to another, we constructed cell developmental trajectories by Monocle2. Cell pseudotime trajectories predicted 3 granulosa cell states and 1 branch point (Fig. 4. F). State1 was mainly consisted with cells from G2M phase, we assigned this state to Mitotic SGCs. In this stage leptin and WT cells were evenly distributed, showing similar transcriptional dynamics between two samples during early mitotic phase of SGCs. For further characterization of the small granulosa cell lineage, top 100 differentially expressed genes of the 3 states were divided in 4 clusters according to their distinct expression patterns, and a heatmap was generated (Fig. 4. G). Notably, state3 highly expressed SOX5, BCAS3, MITF and TOX that is essential for cell migration and immune response, showing state 3 is corresponding to immunoreactive cell. As state 3 is mainly consist of leptin cells, we confirmed the inflammation state in leptin pig occurs at later cell differentiation stage. Additionally, state2 exhibited sharply increased expression of metabolic genes, such as COX1, INHA, CYP19A1 and HSD17β1, while mitotic SGCs had higher expression of MEIOB, HMGB2 and SGO2. The gene plot (Fig. 4H) showed two kinetic trends for state 2 or state 3, from branch point to the end of trajectories. DIAPH3, a cell mitotic gene, downregulated as cells went along the trajectory in both states, indicating loss of mitotic potential during the development process. We also showed TOX, an immune cell differentiation gene, which is upregulated during the developmental process in state 3. HSD17β1 is a hormone metabolic enzyme; it is upregulated in state 2, indicating gain of metabolic function during normal development.
In order to describe the pre-transcriptomic profile of SGC subtypes, we further investigated regulons activity of transcriptional factors in SGCs using SCENIC (Fig. S2). Based on 125 regulons and 6048 filtered genes with default filter parameters, the cells were clustered by 4 subcluster identities of SGCs in t-distributed stochastic neighbor embedding (t-SNE) method calculated by SCENIC. The regulon activity was binarized and matched with cell identities and some representative regulons and their motifs were listed. As shown, SOX5, SOX6, STAT5A and RUNX2 were active mainly in the Immunoreactive cells, while IVD, MYBL1, E2F8 and RAD21 were mostly active in Mitotic cells, indicating higher cell proliferation and signaling activities in these sub-clusters. Moreover, several TFs, such as TBL1XR1 and TCF4, were significantly represented in all four sub-clusters, which supports the role of these TFs as likely candidates for sustaining SGC-specific transcription programs during all stages.
Differential expressed genes reveal altered hormone metabolic process in leptin pig
To further dissect the mechanism for high leptin-induced infertility, we compared all cells of the leptin group and WT group and chose the top 400 DEGs to perform GO analysis. Our results revealed that the overall differential expressed genes were enriched in the lipid metabolic process and steroid metabolic process (Fig. 5. A), which are important for hormone homeostasis. This suggests that the leptin group may experience changes in hormone metabolism compared to the WT group. The GO term analysis revealed that the DEGs in leptin cells were mainly involved in steroid biosynthesis and other hormone synthesis-related pathways like cortisol, aldosterone and thyroid hormone synthesis, and metabolism-related pathways. While the DEGs in WT cells were enriched in tissue development and biological adhesion and cell differentiation pathways (Fig. S2. D). Collectively, these findings indicated that the DEGs in leptin cells exhibited mainly hormone biosynthesis and lipid metabolism-related pathways.
Differential expressed genes revealed altered metabolic process in leptin pigs. (A) Representative GO terms of DEGs between leptin and WT cells. (B) Vln plots of the expression level of representative DEGs involved in reproductive and metabolic process across all cell population in leptin and WT pig ovaries. (C-D) Vln plots of the expression level of genes involved in cholesterol biosynthesis in leptin and WT cells in four clusters of SGCs. (E) NetView plot showing interactions between different cell types. The width of edges represents the strength of the communication. Labels on the edges show exactly how many interactions exist between two types of cells. (F-J) Ligand receptor interaction between different cell populations. All cell types with all ligand-receptor categories (F-G). Interaction of luteal cell and macrophage by top 20 differential expressed cytokines between leptin and WT groups (H). Interaction of macrophage and SGCs through top 20 differential expressed growth factors (I) and other factors (J)
To further evaluate the metabolic heterogeneity of SGCs, we analyzed DEGs of SGCs between WT and leptin groups. GO terms including “cholesterol metabolic process” and “steroid metabolic process” were enriched in upregulated genes of the leptin group (Fig. 5. A). Interestingly, genes related to cell proliferation and metabolism such as KIT, INHA, CYP19A1, HSD17B1, and HSD17B4 showed similar expression patterns across cell types. These were mainly expressed in GCs and downregulated in the leptin group, suggesting GCs are the main site regulating ovary development and function impairment in leptin pigs. While lipoprotein APOA was expressed in all cell types and slightly deceased in luteal cells and TICs, and APOE was downregulated in leptin macrophages (Fig. 5. B). Additionally, pro-apoptosis genes like RARRES1, EPB41L3, IRAK2, and reproductive hormone receptor FSHR were upregulated in the leptin group (Fig. 5. D). Furthermore, genes participating in the cholesterol biosynthesis pathway, such as SCARB1, STAR, HMGCS1, and SQLE have a higher level in the leptin group (Fig. 5. C), indicating increased cholesterol biosynthesis. Strikingly, these genes were equally expressed at all four states of SGC except SQLE which was specifically exhibited in the mitotic state, suggesting a varying metabolic pattern during the differentiation of SGCs in the leptin group (Fig. 5. C). The results were in line with the phenotype observed in leptin pigs, suggesting disrupted hormone homeostasis and providing insights into the mechanisms underlying infertility in leptin pigs, particularly concerning cholesterol biosynthesis, granulosa cell function, and apoptosis regulation.
Cell-cell communication indicates inflammation of corpus luteum and follicle growth arrest
To further explore the crosstalk among different clusters of pig ovaries, we established the cell-cell communication network of all cell clusters with the iTALK package (Fig. 5. E, F). Remarkably, the majority of ligand-receptor pairs were identified between large granulosa cells (LGCs) and theca-interstitial cells (TICs), indicating that the LGC and TIC complexes serve as the main hub for cell-cell communication. These cells were also found to exchange signals with other cell types such as macrophages and endothelial cells. Moreover, LGCs and TICs had crucial roles in sustaining steroid hormone metabolic homeostasis by APOE-SCARB1 ligand-receptor pairs and regulating follicle development via TIMP1-CD63, FST-BMPR1B, and SLIT2-ROBO1 ligand-receptor pairs. Additionally, when it comes to cytokines, nearly all the receptors were derived from luteal cells, and the associated ligands were enriched in immune cells including macrophages, lymphocyte cluster 1 (lym1), and lymphocyte cluster 2 (lym2), which suggests the important roles of luteal cells in ovarian immune response (Fig. 5. G).
Subsequently, we performed a comparative analysis of important cell-cell interactions such as those between macrophage and luteal cells, macrophage, and SGCs (Fig. 5. H-J). The analysis of cytokine communication between macrophages and luteal cells revealed upregulation of ligands such as CCL2, CCL4, and CCL8, along with receptors NRP1, TLR2, and DYSF in macrophages of leptin pigs (Fig. 5. H). This finding suggests that leptin may mediate inflammation in the corpus luteum by activating macrophages. Additionally, the analysis of growth factor family interactions between macrophages and luteal cells showed upregulation of follicle-stimulating hormone receptor (FSHR) in leptin-treated SGCs, indicating a potential reduction in the FSH sensitivity of granulosa cells (Fig. 5. I). Furthermore, growth-associated factors including KIT, RELN, and INHA were downregulated in SGCs, implying inhibition of follicle growth (Fig. 5. J). In a nutshell, the cell-cell communication analysis revealed that inflammation of the corpus luteum and follicle growth arrest may be attributed to altered interactions between macrophages and luteal cells, as well as changes in growth factor signaling within SGCs. These findings shed light on the complex interplay between different cell types in the ovary and provide insights into the potential mechanisms underlying the reproductive impairments observed in leptin-overexpressing pigs.
Discussion
Leptin overexpressing transgenic female pigs had reduced body weight, growth, and fat depots along with delayed puberty, irregular estrous cycles, and greater number of matings. This was linked to metabolic imbalances, increased immune response, and altered ovarian functions. This study provides a theoretical basis for the complex mechanisms underlying leptin, and reproduction by employing leptin-overexpressing female pigs.
The constitutional delay in growth is directly related to leptin concentrations [28]. The leptin pigs exhibited delayed body weight and growth that coincides with the previous findings in GE mouse models overexpressing leptin [29] and is attributed to the anorexigenic effect of leptin, as chronic administration of leptin has been shown to reduce body weight and growth [30, 31]. Previously, we reported that high leptin in pigs suppressed insulin and impacted glucose homeostasis. Therefore, this could be the possible mechanism behind the sluggish growth in leptin-overexpressing pigs [12]. Moreover, the transgenic leptin pigs also exhibited low back fat deposits, likely due to homeostasis imbalance in lipolysis and lipogenesis as evidenced by our previous study [18]. Therefore, the sluggish growth of these leptin pigs could be due to the anorexigenic and adipolytic activity of leptin.
Puberty is a fascinating developmental transition that gates the attainment of reproductive capacity and culminates in the somatic and sexual maturation of the organism [32]. Shreds of evidence from animal models and human studies suggested that leptin acts as a permissive factor for the onset of puberty and maintenance of reproductive cyclicity [11]. However, in our study, leptin overexpressing pigs displayed a late onset of puberty followed by irregular estrus comprising of short and long estrus and breeding insufficiency characterized by more mating. Contrarily, GE mice overexpressing leptin manifested accelerated puberty and maintained fertility at earlier ages, despite being lean like our leptin overexpressing female pigs [7]. This discrepancy may be attributed to species differences and threshold levels of leptin as several reports have indicated that the physiological doses of leptin influencing anorexigenic effect are not equally enough to boost reproductive functions as both functions of leptin are mediated via different hypothalamic pathways [11, 33, 34]. Therefore, it is also plausible to speculate that our leptin female pigs may not attain enough leptin to boost reproduction and hence exhibit reproductive impairment, however, such mechanisms need further exploration.
The initiation and maintenance of reproductive functions normally depend upon gonadotropins that in turn cause the production of steroids. However, endocrine hormonal imbalances due to inappropriate secretion of gonadotropins lead to ovulatory dysfunction that results in reproductive impairment [35]. In this study, we observed reproductive disorders in leptin pigs were characterized by elevated basal levels of gonadotropins (FSH & LH) and increased steroid production. This could be attributed to high leptin, as it is well-accepted that leptin augments the secretions of gonadotropins [36]. According to the two-cell–two-gonadotropin theory, LH stimulates thecal-interstitial cells to produce androgens, and FSH stimulates granulosa cells to produce estrogens from androgens [37]. Our study showed both FSH and FSHR, and its target granulosa cell fraction was elevated in leptin pigs, which aligns with the elevated E2 level. We also observed elevated TICs and endothelial cells in the ovaries of leptin pig corresponding to higher P4 and vascularization. These findings indicate the ovarian axis was overactivated in leptin pigs. It has been reported that elevated LH levels impair downstream ovarian folliculogenesis as the ‘FSH threshold’ required for follicle maturation is frequently not reached, causing follicles to arrest in a preovulatory stage, giving rise to cystic ovaries which results in menstrual irregularities [38]. This follicle arrest was also evident from the cell-cell communication network as KIT and INHA were downregulated in the SGC of the leptin group (Fig. 5). Enhanced production of KIT ligand has been shown to stimulate the oocyte growth through PI3K signaling pathways, while low expressions were insufficient to activate KIT on oocyte and thus maintained the oocyte in quiescent stage [39]. Furthermore, a reduced luteal cell population, as well as elevated granulosa cells in leptin ovaries, pretends the impaired granulosa cells’ transformation to luteal cells, which might impair ovulation. Nonetheless, apoptosis is one of the necessary pathways for luteolysis. While downregulated apoptosis pathways in leptin pigs’ ovaries-might be due to reduced luteal cell count, also seems to impede the luteolysis which was also evident from high levels of P4 at the end of the estrous cycle. Thus, it is intriguingly speculated that anovulation and incomplete luteolysis may lead to irregular estrous in leptin pigs.
Steroids are synthesized from CE, highlighting the importance of cholesterol homeostasis for the maintenance of fertility. Disturbances in serum cholesterol levels can affect steroidogenesis and reproductive functions. In this study, leptin pigs exhibited high serum levels of CE which was also evident from snRNA sequencing that genes involved in cholesterol synthesis i.e. HMGCS1 and SQLE were upregulated in leptin ovaries. These high levels of CE along with the upregulation of genes mediating the transport of cholesterol from the outer mitochondrial membrane to the inner mitochondrial membrane i.e. StAR and SCARB1 (Fig. 5C) pretend the increased steroidogenesis (P4) in leptin pigs. Furthermore, α-LA has been shown to have an inhibitory effect on cholesterol synthesis [40, 41], and arachidonic acid (AA) and its metabolites have long been implicated in steroidogenesis through direct effects on the steroidogenic machinery (e.g., acute steroid regulator [StAR] [42]. Thus, the low levels of α-LA, and high levels of AA along with enhanced expression of the regulatory protein (StAR) in leptin pigs indicate enhanced steroidogenesis which might be the cause of the estrous irregularity.
Macrophages play crucial and diverse roles in various intra-ovarian events, encompassing folliculogenesis and the formation and regression of the corpus luteum [43]. Our study reveals macrophages communicated with different ovarian cell types, including granulosa cells, theca-interstitial cells, and luteal cells, which are mediated by cytokine ligand-receptor interactions. On the one hand, the development of follicles relies on the presence of an appropriate cytokine and growth factor environment, and macrophages serve as a significant source of these factors. In leptin pigs, we observed that the SGC and macrophage interaction is biased in both the cytokine and growth factor categories. These findings coincide with previous findings in polycystic ovarian syndrome (PCOS) patients in which GC expressed elevated levels of transcripts encoding cytokines, chemokines, and immune cell markers [44]. On the other hand, macrophages play a pivotal role in various aspects of corpus luteum function. First, they contribute to angiogenesis within the corpus luteum by secreting important factors, our study showed increased VCAM1 and ITGB1 in macrophage-luteal cell communication of leptin pigs, indicating well-established vasculature in leptin pigs. This is also evidenced by the presence of rich vessels in the H&E section of leptin pig ovaries. Notably, the number of macrophages within the corpus luteum reaches its peak during regression, indicating their involvement in the process of luteolysis. The ovaries of leptin pig exhibited fewer macrophages than that of WT pigs did, indicating that fewer luteal cells need to be eliminated. Therefore, we hypothesized that alterations in the macrophage-mediated regulation of folliculogenesis and corpus luteum function contributed to the pathogenesis of ovarian disorders.
Conclusions
In this study, we found that leptin overexpression in pigs adversely affects reproductive performance, causing stunted growth, delayed puberty, irregular estrous cycles, and reduced breeding efficiency. This is linked to metabolic imbalances, increased immune response, and altered ovarian functions. This study provides insights into the complex mechanisms underlying ovarian function and suggests a potential role for leptin in ovarian pathophysiology. These findings provide some theoretical basis for the reproductive role of leptin.
Data availability
The dataset supporting the conclusions of this article is included within the article and its additional file. The raw data for snRNA sequencing are submitted as GSE237153.
Abbreviations
- CYP19A1:
-
Cytochrome P45019A1
- E2:
-
Estradiol
- FSH:
-
Follicle-stimulating hormone
- FSHR:
-
Follicle-stimulating hormone receptor
- GCs:
-
Granulosa cells
- GE:
-
Gene-edited
- GSEA:
-
Gene set enrichment analysis
- LC:
-
Luteal cells
- LGC:
-
Large granulosa cells
- LH:
-
Luteinizing hormone
- LHCGR:
-
Luteinizing hormone/choriogonadotropin receptor
- P4:
-
Progesterone
- PUFA:
-
Poly unsaturated fatty acids
- SFA:
-
Saturated fatty acids
- SGC:
-
Small granulosa cells
- snRNA-seq:
-
Single-nucleus RNA sequencing
- StAR:
-
Steroidogenic acute regulatory protein
- TIC:
-
Theca-interstitial cells
- WT:
-
Wild-type
- α-LA:
-
α-linolenic acid
References
Karoutsos P, Karoutsou E, Karoutsos D. Obesity and female fertility: the bridging role of leptin. J Data min Genom Proteom. 2017;8:1–4. https://doi.org/10.4172/2153-0602.1000209.
Tong Q, Xu Y. Central leptin regulation of obesity and fertility. Curr Obes Rep. 2012;1:236–44. https://doi.org/10.1007/s13679-012-0025-8.
He Q, Karlberg J. BMI in childhood and its association with height gain, timing of puberty, and final height. Pediatr Res. 2001;49:244–51. https://doi.org/10.1203/00006450-200102000-00019.
Laron Z. Is obesity associated with early sexual maturation? Pediatrics. 2004;113:171–2. https://doi.org/10.1542/peds.110.5.903.
Wronka I, Suliga E, Pawlińska-Chmara R. Socioeconomic determinants of underweight and overweight in female Polish students in 2009. Anthropol Anz. 2012;85–96. https://doi.org/10.1127/0003-5548/2011/0131.
Wade GN, Jones JE. Neuroendocrinology of nutritional infertility. Am J Physiol Regul Integr Comp Physiol. 2004;287:R1277–96. https://doi.org/10.1152/ajpregu.00475.2004.
Yura S, Ogawa Y, Sagawa N, Masuzaki H, Itoh H, Ebihara K, et al. Accelerated puberty and late-onset hypothalamic hypogonadism in female transgenic skinny mice overexpressing leptin. J Clin Investig. 2000;105:749–55. https://doi.org/10.1172/JCI8353.
Smith JT, Spencer SJ. Preweaning over- and underfeeding alters onset of puberty in the rat without affecting kisspeptin. Biol Reprod. 2012;86(145):141–8. https://doi.org/10.1095/biolreprod.111.097758.
AL-Thuwaini TM. Adiponectin and its physiological function in ruminant livestock. Rev Agri Sci. 2022;10:115–22. https://doi.org/10.7831/ras.10.0_115.
AL-Thuwaini TM, Rhadi FA, Kadhim AH. Adipokines and reproduction in ruminant livestock. Rev Agri Sci. 2024;12:13–23. https://doi.org/10.7831/ras.12.0_13.
Egan OK, Inglis MA, Anderson GM. Leptin signaling in AgRP neurons modulates puberty onset and adult fertility in mice. J Neurosci. 2017;37:3875–86. https://doi.org/10.1523/JNEUROSCI.3138-16.2017.
Chen J, Zeng W, Pan W, Peng C, Zhang J, Su J, et al. Symptoms of systemic lupus erythematosus are diagnosed in leptin transgenic pigs. PLoS Biol. 2018;16:e2005354. https://doi.org/10.1371/journal.pbio.2005354.
Mordhorst B, Prather R. Pig models of reproduction. Anim Models Hum Reprod. 2017;213–34. https://doi.org/10.1002/9781118881286.
Ding J, Adiconis X, Simmons SK, Kowalczyk MS, Hession CC, Marjanovic ND, et al. Systematic comparison of single-cell and single-nucleus RNA-sequencing methods. Nat Biotechnol. 2020;38:737–46. https://doi.org/10.1038/s41587-020-0465-8.
Fan X, Bialecka M, Moustakas I. Single-cell reconstruction of follicular remodeling in the human adult ovary. Nat Commun. 2019;10:3164. https://doi.org/10.1038/s41467-019-11036-9.
Wagner M, Yoshihara M. Single-cell analysis of human ovarian cortex identifies distinct cell populations but no oogonial stem cells. Nat Commun. 2020;11:1147. https://doi.org/10.1038/s41467-020-14936-3.
Wang S, Zheng Y, Li J, Yu Y, Zhang W, Song M, et al. Single-cell transcriptomic atlas of primate ovarian aging. Cell. 2020;180:585–e600519. https://doi.org/10.1016/j.cell.2020.01.009.
Qing Y, Jamal MA, Shi D, Zhao S, Xu K, Jiao D, et al. Delayed body development with reduced triglycerides levels in leptin transgenic pigs. Transgenic Res. 2021. https://doi.org/10.1007/s11248-021-00288-1.
Knauer M, Cassady J, Newcom D, See MT. Phenotypic and genetic correlations between gilt estrus, puberty, growth, composition, and structural conformation traits with first-litter reproductive measures. J Anim Sci. 2011;89:935–42. https://doi.org/10.2527/jas.2009-2673.
Wang J, Zhang L, Xiao R, Li Y, Liao S, Zhang Z, et al. Plasma lipidomic signatures of spontaneous obese rhesus monkeys. Lipids Health Dis. 2019;18:1–11. https://doi.org/10.1186/s12944-018-0952-9.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60. https://doi.org/10.1038/nmeth.3317.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102:15545–50. https://doi.org/10.1073/pnas.0506580102.
Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM, et al. Comprehensive integration of single-cell data. Cell. 2019;177:1888–902. https://doi.org/10.1016/j.cell.2019.05.031. e1821.
Butler A, Hoffman P, Smibert P, Papalexi E, Satija R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol. 2018;36:411–20. https://doi.org/10.1038/nbt.4096.
Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics J Integr Biol. 2012;16:284–7. https://doi.org/10.1089/omi.2011.0118.
Aibar S, González-Blas CB, Moerman T, Huynh-Thu VA, Imrichova H, Hulselmans G, et al. SCENIC: single-cell regulatory network inference and clustering. Nat Methods. 2017;14:1083–6. https://doi.org/10.1038/nmeth.4463.
Khristi V, Chakravarthi VP, Singh P, Ghosh S, Pramanik A, Ratri A, et al. ESR2 regulates granulosa cell genes essential for follicle maturation and ovulation. Mol Cell Endocrinol. 2018;474:214–26. https://doi.org/10.1016/j.mce.2018.03.012.
El-Eshmawy MM, Abdel Aal IA, Hawary AK. Association of ghrelin and leptin with reproductive hormones in constitutional delay of growth and puberty. Reprod Biol Endocrinol. 2010;8. https://doi.org/10.1186/1477-7827-8-153.
Ogawa Y, Masuzaki H, Hosoda K, Aizawa-Abe M, Suga J, Suda M, et al. Increased glucose metabolism and insulin sensitivity in transgenic skinny mice overexpressing leptin. Diabetes. 1999;48:1822–9. https://doi.org/10.1159/000438877.
Ravussin Y, LeDuc CA, Watanabe K, Mueller BR, Skowronski A, Rosenbaum M, et al. Effects of chronic leptin infusion on subsequent body weight and composition in mice: can body weight set point be reset? Mol Metab. 2014;3:432–40. https://doi.org/10.1016/j.molmet.2014.02.003.
31, George BT, Jhancy M, Dube R, Kar SS, Annamma LM. The molecular basis of male infertility in obesity: a literature review. Int J Mol Sci. 2023. https://doi.org/10.3390/ijms25010120. 25.120.
Castellano JM, Tena-Sempere M. Animal modeling of early programming and disruption of pubertal maturation. Puberty from bench to clinic. Volume 29. Karger; 2016. pp. 87–121. https://doi.org/10.1159/000438877.
Henry BA, Goding JW, Alexander WS, Tilbrook AJ, Canny BJ, Dunshea F, et al. Central administration of leptin to ovariectomized ewes inhibits food intake without affecting the secretion of hormones from the pituitary gland: evidence for a dissociation of effects on appetite and neuroendocrine function. Endocrinol. 1999;140:1175–82. https://doi.org/10.1210/endo.140.3.6604.
Plant TM, Durrant AR. Circulating leptin does not appear to provide a signal for triggering the initiation of puberty in the male rhesus monkey (Macaca mulatta). Endocrinol. 1997;138:4505–8. https://doi.org/10.1210/endo.138.10.5574.
Abedal-Majed MA, Cupp AS. Livestock animals to study infertility in women. Anim Front. 2019;9:28–33. https://doi.org/10.1093/af/vfz017.
Hausman GJ, Barb CR, Lents CA. Leptin and reproductive function. Biochimi. 2012;94:2075–81. https://doi.org/10.1016/j.biochi.2012.02.022.
Mlyczyńska E, Kieżun M, Kurowska P, Dawid M, Pich K, Respekta N, et al. New aspects of corpus luteum regulation in physiological and pathological conditions: involvement of adipokines and neuropeptides. Cells. 2022;11:957. https://doi.org/10.3390/cells11060957.
Sun BZ, Kangarloo T, Adams JM, Sluss P, Chandler DW, Zava DT, et al. The relationship between progesterone, sleep, and LH and FSH secretory dynamics in early postmenarchal girls. J Clin Endocrinol Metab. 2019;104:2184–94. https://doi.org/10.1210/jc.2018-02400.
Zhang H, Risal S, Gorre N, Busayavalasa K, Li X, Shen Y, et al. Somatic cells initiate primordial follicle activation and govern the development of dormant oocytes in mice. Curr Biol. 2014;24:2501–8. https://doi.org/10.1016/j.cub.2014.09.023.
Hu X, Weng X, Tian Y, Wang C, Yang Y, Xu K, et al. Effects of omega-3 polyunsaturated fatty acids on steroidogenesis and cellular development in PCOS rats. Food Func. 2019;10:2504–14. https://doi.org/10.1039/c8fo02319k.
Mumford S, Schisterman E, Dasharathy S, Pollack A, Zhang C, Wactawski-Wende J. Omega-3 fatty acids and ovulatory function. Fert Steril. 2011;96:S15. https://doi.org/10.1016/j.fertnstert.2011.07.065.
Wang X, Walsh LP, Reinhart AJ, Stocco DM. The role of arachidonic acid in steroidogenesis and steroidogenic acute regulatory (StAR) gene and protein expression. J Biol Chem. 2000;275:20204–9. https://doi.org/10.1074/jbc.M003113200.
Wu R, Van der Hoek KH, Ryan NK, Norman RJ, Robker RL. Macrophage contributions to ovarian function. Hum Reprod Update. 2004;10:119–33. https://doi.org/10.1093/humupd/dmh011.
Adams J, Liu Z, Ren YA, Wun WS, Zhou W, Kenigsberg S, et al. Enhanced inflammatory transcriptome in the granulosa cells of women with polycystic ovarian syndrome. J Clin Endocrinol Metab. 2016;101:3459–68. https://doi.org/10.1210/jc.2015-4275.
Acknowledgements
We thank the Department of Science and Technology of Yunnan Province, China for their financial support.
Funding
This work was supported by The Major Science and Technology Project of Yunnan Province (202102AA100054); Spring City Plan "The High-level Talent Promotion and Training Project of Kunming (2022SCP001)".
Author information
Authors and Affiliations
Contributions
HJW conceived the project. HJW, HYZ, and RZ designed the research experiment. MAJ, YC, DJ, JG, HZ, TW, HL, WC, DZ, WX, SP and YC conducted experiments. HYZ, KX, DJ, MAJ, YQ and BJ acquired data. MAJ, YC, and DJ analyzed data. HJW, HYZ, MAJ, YC, KX and RZ wrote the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
Animals used in this study were regularly maintained in the Laboratory Animal Centre of Yunnan Agricultural University. All animal experiments were performed with the approval of the Animal Care and Use Committee of Yunnan Agricultural University, Kunming, China.
Consent for publication
Not applicable.
Competing interests
The authors declared that they have no conflicts of interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Jamal, M.A., Cheng, Y., Jiao, D. et al. Unraveling the impact of hyperleptinemia on female reproduction: insights from transgenic pig model. Biol Res 57, 60 (2024). https://doi.org/10.1186/s40659-024-00545-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40659-024-00545-7