Abstract
Paternal dietary patterns and obesity can affect offspring through epigenetic signals in sperm RNA. By activating hypothalamic Agouti-related peptide (AgRP) neurons in mice to simulate obesity’s effects, we observed changes in the small noncoding RNA payload of sperm. These alterations, particularly in transfer RNA-derived small RNAs (tsRNAs), are similar to those induced by short-term high-fat diets, suggesting a common upstream regulatory mechanism involving AgRP neurons that affects metabolic epigenetic inheritance.
Similar content being viewed by others
Introduction
Paternal metabolic state and dietary patterns affect the progeny through a cascade of epigenetic signals. These modifications, transported by the sperm RNA payload, are now recognized as vectors transmitting nuanced information that shapes offspring phenotypes1,2,3,4. Recent studies have shown that transfer RNA-derived small RNAs (tsRNAs) in sperm are important mediators of epigenetic inheritance, particularly influencing the metabolic pathways that predispose offspring to obesity4,5,6,7,8,9. The mechanisms underlying sperm RNA modifications are complex. Sperm RNA content can be altered through several established pathways: direct environmental effects on developing sperm cells in the testes, RNA transfer from epididymal cells via epididymosomes, mitochondrial dysfunction, and RNA modifications and processing4,5,10,11,12.
This complexity of RNA modification pathways suggests multiple potential points where metabolic signals could influence the sperm RNA payload. In this context, hypothalamic Agouti-related peptide (AgRP) neurons emerge as particularly interesting candidates for investigation, as these neurons are central players in the orchestration of energy balance and are known to modify their activity in response to excessive weight gain, reflecting an intricate dialogue between nutritional status and neuroendocrine function13. The influence of AgRP neurons extends beyond metabolic regulation, as their activation triggers a complex network of downstream effects that impact various physiological processes such as reproductive function and cognitive function among others14,15,16. While the transformative power of the sperm RNA payload is acknowledged, particularly concerning metabolic traits, the precise signals that sculpt these epigenetic landscapes remain to be discovered.
Here, we aimed to mimic the hyperactive state of AgRP neurons one piRNA that occur in obesity17,18. To achieve this, we employed Designer Receptors Exclusively Activated by Designer Drugs (DREADD), a well-known chemogenetic tool for selectively activating specific neuronal populations in the mouse brain19 .We bilaterally injected an excitatory DREADD adeno-associated viral (AAV-hM3Dq-mCherry) vector into the arcuate nucleus of the hypothalamus (ARH) in both AgrpCre/+ and Agrp+/+ (Fig. 1a). Three weeks post-injection, we confirmed the restriction of AAV to AgRP neurons and validated its efficacy (Supplemental Fig. 1a-c). The DREADD ligand, Clozapine-N-Oxide (CNO), was administered to both AgrpCre/+ and control Agrp+/+ mice in the absence of food. Four hours post-administration, mature sperm was harvested from the caudae epididymides, and RNAs were extracted from the sperm. The purity of sperm samples was validated by qPCR analysis of cell-type specific markers. This analysis showed strong enrichment for mature sperm markers such as Protamine-2 (Prm2) and low enrichment for somatic and stem cell markers such as E-cadherin and C-kit (Supplemental Fig. 1d). These RNAs were then used to prepare a small non-coding RNA (sncRNA) library, which was subsequently sequenced (Fig. 1a).
Activation of hypothalamic AgRP neurons leads to alteration on the sperm RNA payload. (a) Experimental design of this study. (b) Distribution of sncRNA biotypes in cauda spermatozoa from control (Agrp+/+ + CNO) and activated (AgrpCre/+ + CNO) mice. n = 4 per group. microRNAs (miRNAs), PIWI-interacting RNAs (piRNAs), transfer RNA-derived small RNAs (tsRNAs) and ribosomal small RNAs (rsRNAs). (c) Volcano plot representation of differentially expressed small non coding RNAs ( sncRNAs) fragments (false discovery rate < 0.05 and log2[fold change (FC)] > |1|). n = 4 per group. (d) Heat map representation of differentially expressed piRNA, miRNA, rsRNA and tsRNA fragment levels in cauda spermatozoa from control and activated mice. n = 4 per group. (e) Chord diagram illustrating the overlap of significantly (false discovery rate < 0.05) altered sncRNAs between our study and previous studies examining the impact of high-fat diet (HFD) or Western diet (WD) feeding on sperm RNA payload composition. (f) Relative abundance of tsRNAs (Asp-GTC, Ser-TGA-5’, Glu-TTC-5’, Glu-CTC-5’) differentially expressed fragments in control (n = 4) and AgRP-activated (n = 4) groups (this study). (g) Relative abundance of tsRNAs (Asp-GTC, Glu-TTC-5’, Glu-CTC-5’) differentially expressed fragments in control (n = 4) and 2-weeks HFD (n = 4) groups (Tomar et al., 2024). (h) Relative abundance of tsRNAs (Asp-GTC, Ser-TGA-5’, Glu-CTC-5’) differentially expressed fragments in control (n = 3) and 3-months WD (n = 3) groups (Raad et al., 2021). (i) Gene ontology pathway analysis of the target genes of miR-486a and b (false discovery rate < 0.05). Data are presented as mean ± s.e.m. and were analyzed by unpaired two-tailed Student’s t-tests (f-h).
After mapping reads to the reference genome and sncRNA annotation databases, we confirmed that the distribution and biotypes between the two groups were very similar (Fig. 1b). Sperm of control mice had microRNAs (miRNAs) at 32%, PIWI-interacting RNAs (piRNAs) at 14%, transfer RNA-derived small RNAs (tsRNAs) at 15% and ribosomal small RNAs (rsRNAs) at 25%. In comparison, sperm of AgRP neuron-activated mice showed miRNAs at 33%, piRNAs at 16%, tsRNAs at 12% and rsRNAs at 22%. This similarity suggests that the global sncRNA profiles of the two groups are largely comparable, indicating that broad changes in RNA expression or biotype distribution in sperm are not the primary effects of AgRP neuron activation on the sperm RNA payload.
Nevertheless, a closer examination revealed specific alterations in tsRNA fragment sizes upon AgRP neuron activation. We observed significant changes in tsRNA sizes of 30, 31 bp (P value < 0.05), and 32 bp (P value = 0.05), with a reduction of approximately 40% (Supplemental Fig. 2a). The rapid changes observed in tsRNA composition following AgRP neuron activation are particularly intriguing and could be explained by the dynamic nature of tsRNA biogenesis. tsRNAs are produced through enzymatic processing of mature tRNAs by specific enzymes, notably DNMT2 and angiogenin10,20. Additionally, single-cell analysis of the epididymis has revealed an abundance of nucleases capable of such RNA processing, highlighting the epididymis as a key site for rapid RNA modifications21. We hypothesize that AgRP neuron activation may influence these enzymatic processing pathways, providing a mechanistic explanation for the acute changes we observe in sperm tsRNA content. This mechanism would allow for rapid modifications in RNA fragments without requiring de novo RNA synthesis, consistent with the timeframe of our observations.
In contrast, miRNAs, piRNAs and rsRNAs did not show differences in the abundances of specific fragment sizes, suggesting that their processing and stability are not significantly affected by AgRP neuron activation (Supplemental Fig. 2b-d).
Building on these observations, we carried out a differential expression analysis to pinpoint sncRNAs with significant expression differences between the control and AgRP neuron-activated groups. Unexpectedly, this analysis revealed 195 sncRNAs that were upregulated and 45 sncRNAs that were downregulated in the sperm following AgRP neuron activation (Padj < 0.05 and log2[fold change (FC)] > |1|), constituting approximately 2% of the total sncRNA pool (Fig. 1c). Although these changes represent only a small fraction of the sncRNA payload (2%), this is expected given the acute nature of the treatment. Among the affected biotypes, 26% were piRNA, 6% were rsRNA, 6% were tsRNA, and 2% were miRNA (Supplemental Fig. 2e). Additionally, a high fraction (60%) of the differentially expressed sncRNAs were without annotation, so our analysis focused on the annotated sequences to provide clearer insights into the specific sncRNAs influenced by AgRP neuron activation (Supplemental Fig. 2e). Among the differentially expressed sncRNAs that we detected, all tsRNAs, miRNAs, and rsRNAs were downregulated. Only piRNAs showed both upregulation and downregulation, indicating a more complex regulatory pattern for this biotype compared to the others (Fig. 1d).
To further contextualize our findings, we compared our data with three previous studies that analyzed the impact of short- and long-term obesogenic diets on sperm RNA payload composition4,5,9. Given the established role of AgRP neurons in regulating whole-body energy homeostasis, this comparison is particularly relevant. This intersection analysis revealed that the alterations in tsRNA-Ser-TGA-5’ and tsRNA-Glu-TTC-5’ are consistently observed across multiple studies (Fig. 1e). The recurrence of these specific tsRNAs across different experimental conditions underscores their potential role as key markers or mediators in the physiological responses to both AgRP neuron activation and dietary changes. This is particularly significant given the established role of tsRNAs in the epigenetic inheritance of metabolic traits, such as obesity4,5,10,12. Additionally, similar patterns of alterations for other tsRNAs, such as tsRNA-Asp-GTC and tsRNA-Glu-CTC-5’, reinforce the idea that specific tsRNAs are reliably affected by various physiological stimuli and dietary interventions, further linking neural activation and metabolic epigenetics (Fig. 1e). Importantly, the expression levels of tsRNAs Asp-GTC, Glu-CTC-5’, and Glu-TGA-5’ were significantly decreased in sperm following short-term high-fat diet (HFD)5, closely mirroring the downregulation observed after AgRP neuron activation (Fig. 1f-g). In contrast, long-term Western diet (WD)9 feeding resulted in a differential response: while tsRNAs Glu-CTC-5’ and Ser-TGA-5’ were decreased, tsRNA Asp-GTC was notably increased in sperm (Fig. 1h). We could not make this comparison with the long-term HFD study because the data from individual samples were not available4. These findings suggest that short-term dietary intervention and AgRP activation have a more similar impact on sperm tsRNA profiles, while long-term HFD produces distinct changes, indicating different underlying mechanisms at play. This unexpected regulation by AgRP neurons highlights a potential upstream regulatory link that influences the response of sperm tsRNAs to varying metabolic conditions.
For the other biotypes, we did not see such consistency. However, we identified miR-451a as significantly altered in both our data and the long-term HFD and WD studies, while miR-486 (a and b) was significantly altered in both our study and the short-term HFD and long-term WD studies. Given the more immediate relevance of short-term dietary changes to our AgRP activation model, we focused our Gene Ontology analysis on miR-486 (a and b) to understand the specific biological pathways they might influence. This analysis revealed that many pathways were related to biological regulation and developmental processes (Fig. 1i). This suggests that while the overlap in miRNA changes is less pronounced, the functional implications of these changes might still be significant, particularly in the context of early developmental stage. Finally, we only saw one piRNA that was significantly altered by a short period of HFD and AgRP activation (Fig. 1e). While our analysis examined individual piRNAs, future studies would benefit from analyzing piRNAs at the cluster level, as this approach better reflects their biological organization and function. Given the complex and still-emerging understanding of piRNA functions, it is difficult to draw definitive conclusions about their role and significance in the context of AgRP neuron activation and HFD feeding. Overall, this indicates that while tsRNAs show a clear and consistent response to both AgRP neuron activation and HFD feeding, the responses of other sncRNA biotypes may be more specific to the type of metabolic perturbation or involve different regulatory mechanisms.
Our findings should be considered in the context of the diverse pathways through which sperm RNA content can be modified. These range from direct effects on developing sperm cells in the testes to epigenetic modifications and RNA transfer via epididymosomes. The rapid changes we observe in sperm RNA content following AgRP neuron activation suggest that these neurons might influence these pathways through systemic signals or neural circuits that affect the reproductive tract.
Our findings should be considered in the context of the multiple pathways that modulate the sperm RNA payload. First, environmental and metabolic challenges exert direct effects on developing sperm cells in the testes, influencing spermatogenesis and the biogenesis of tRNA fragments1,4,5. Second, post-testicular processes critically shape the sperm RNA payload, as exemplified by the rapid modification mediated through the transfer of RNAs via epididymosomes, which facilitates active RNA exchange between epididymal cells and maturing sperm11,12. Finally, emerging evidence suggests that RNA stability and processing in the epididymis are tightly regulated by a variety of RNA-modifying enzymes. In addition to nucleases identified in single-cell analyses21, enzymes such as angiogenin20 and DNMT210 have been implicated in the biogenesis and modulation of tsRNAs. This coordinated enzymatic control may enable a rapid reshaping of the sperm RNA landscape in response to external metabolic cues. The rapid changes we observe in sperm RNA content following AgRP neuron activation imply that these neurons might influence these pathways through systemic signals or neural circuits that affect the reproductive tract. This modulation could involve hormonal or autonomic mechanisms similar to those described in previous works addressing the neural control of reproductive functions15,20. Such mechanisms might enable a fast adjustment in the RNA processing machinery without requiring de novo RNA synthesis, thereby aligning with the acute time frame of our observations.
In summary, our study demonstrates for the first time that AgRP neuron activation significantly alters the sperm sncRNA payload, with notable downregulation of specific tsRNAs, miRNAs, and rsRNAs, while piRNAs exhibit both upregulation and downregulation patterns. The consistent downregulation of tsRNAs parallels earlier findings from studies on HFD and WD feeding. Our results indicate that the detrimental effects of obesity on sperm RNA composition may be partially mediated via a common upstream regulatory mechanism involving AgRP neurons. These results reveal a new and unexpected perspective on the modulation of epigenetic information carried by sperm mediated by hypothalamic AgRP neurons.
While our findings provide novel evidence that AgRP neurons can influence sperm RNA composition, we acknowledge that additional validation approaches would further strengthen these observations. Future work will be valuable to deepen our understanding of the molecular mechanisms connecting hypothalamic AgRP neuron activity to sperm RNA modifications, particularly focusing on the precise pathways mediating these rapid changes in RNA content.
Methods
Animals
The University of Barcelona Ethics Committee approved all animal research (Approval No. CEEA-390/19), which were conducted in accordance with the current laws in Catalonia, Spain, and Europe. Mice were housed in a controlled environment, maintaining a temperature range of 20–24ºC and a 12-hour light-dark cycle. The mice’s health condition was constantly evaluated. Agrp+/+ and AgrpCre/+ were purchased from Jackson Laboratory (Strain #:012899). We crossed AgrpCre/+ with Agrp+/+ to generate more Agrp+/+ with AgrpCre/+. Agrp+/+ were used as control. All the methods and procedures were reported in accordance with ARRIVE guidelines.
Stereotaxic surgery
Eight-week-old AgrpCre/+ and Agrp+/+ male mice, were anesthetized using a Ketamine/Xylazine mixture (i.p., 100 mg/Kg and 10 mg/Kg) and secured in a stereotaxic frame (Kopf Instruments). For pain management and postoperative recovery, the mice were administered Buprenorphine (i.p., 0.3 mg/Kg) prior to surgery and every 12 h thereafter for 72 h. The surgical procedure began by exposing the skull and drilling a small hole for the precise injection of an adeno-associated virus (AAV) into the ARH. The AAV, containing excitatory Designer Receptors Exclusively Activated by Designer Drugs (DREADDs) encoding hSYN-DIO-hM3D(Gq)-mCherry (Addgene - # 44361-AAV9–1.10 × 1013 gc/mL), was delivered bilaterally (300 nL per side) using a 33-gauge needle attached to a 5 µL syringe (NeuroSyringe, Hamilton) at a rate of 50 nL/min. The coordinates used for the injections relative to Bregma were AP: -1.5 mm; ML: ± 0.3 mm; DV: -5.8 mm below the skull surface. After an 8-minute interval, the needle was raised 1 mm, held for an additional minute, and then fully retracted. The surgical site was sealed and the mice were kept on a heating pad at 37ºC until they fully recovered from the anesthesia, after which they were returned to their home cage.
Chemogenetic activation of AgRP neurons
On experimental days, cage food was removed for 1 h. Then mice were intraperitoneally injected with Clozapine-N-oxide (CNO; 1 mg/kg body weight for activation; Tocris Bioscience - #4936) dissolved in sterile saline. Mice were euthanized 4-hours later. Correct injection sites were confirmed in every mouse post-mortem by assessing mCherry signal under fluorescent microscope (Nikon Eclipse Ni-U).
Brain histology
Mice were perfused transcardially with ice-cold phosphate-buffered saline (PBS), followed by a fixative solution (4% PFA). Brains were dissected and post-fixed in 4% PFA at 4 °C for 24 h, then cryoprotected in 30% sucrose in 1x PBS (pH 7.4). Brain sections, 26 μm thick, were prepared using a cryostat (Leica Biosystems, IL, USA) at − 20 °C. For cFOS immunodetection, brain slices containing the ARH were extensively washed in PBS buffer and blocked in 3% BSA, 2% chicken serum in PBS plus 0.2% Triton X-100 for 1 h. Sections were incubated with a rabbit anti-cFOS antibody (1:500; Synaptic Systems - # 226008) in blocking buffer overnight at 4 °C. After washing with PBS 0.2% Triton X-100, slices were incubated with a chicken anti-rabbit Alexa Fluor 488 antibody (1:300) for 1 h at room temperature. After washing with PBS 0.2% Triton X-100, slices were mounted with Vectashield antifade mounting medium with DAPI (Vector Laboratories, USA). Images were taken using an Olympus BX41 fluorescence microscope equipped with a 20x objective.
Mouse mature sperm extraction
Mature sperm extraction from mice was conducted according to the established protocol for murine in vitro fertilization. Initially, a Petri dish was set up with 1mL of 1x PBS in the middle, covered with paraffin oil (#76235 – Sigma). After euthanizing the male mouse via cervical dislocation, the caudae epididymides were removed and placed on the paraffin oil in the prepared dish. Each cauda was then cut using scissors at the midpoint to release the sperm clusters, which were gently moved toward the PBS in the center of the dish. The sperm remained in the PBS for 10 min at room temperature. Subsequently, the solution with the sperm was collected and spun at 600 g at 4ºC for 5 min. Finally, the sperm underwent two additional washes with 1x PBS.
Small non-coding RNA sequencing (sncRNA-seq)
Mature sperm RNAs were extracted using Trizol (#15596026 – Thermo Fisher Scientific) following standard protocols. Quantification of each samples was performed using a Qubit fluorometer (Thermo Fisher Scientific). RNA integrity analysis was performed by analyzing each sample in a Agilent High Sensitivity RNA ScreenTape System (Agilent). Small non-coding RNA libraries were generated with 100ng of RNA with NEXTFLEX Small RNA-Seq kit v4 with UDIs (NOVA-5132-31 – Perkin Elmer). Libraries were sequenced on an Illumina NextSeq2000 (Illumina, Inc.) in single-end mode with a read length of 1 × 50 bp.
Sperm purity assessment by qPCR
RNA was reverse transcribed using the High-Capacity cDNA Reverse Transcription Kit (#4368813, Thermo Fisher) following manufacturer’s instructions. qPCR was performed using TaqMan Gene Expression Master Mix (Applied Biosystems, #4369016) with TaqMan probes for Prm1 (Mm01342731_g1, Thermo Fisher), Prm2 (Mm03048199_m1, Thermo Fisher), C-kit (Mm00445212_m1, Thermo Fisher) and E-cadherin (Mm01247357_m1, Thermo Fisher) on a QuantStudio 5 Real-Time PCR System (Applied Biosystems). Relative gene expression was calculated using the 2^-ΔΔCt method with Prm1 as the reference gene.
SncRNA-seq quality control and data analysis
After sequencing, we used the Small non-coding RNA annotation Pipeline Optimized for rRNA- and tRNA-Derived Small RNAs (SPORTS1.1)22 to trim, map and annotate fragments. We used SPORTS1.1 with default parameters in which reads were trimmed and filtered for low quality base-pairs. Sequencing duplicates were not removed. In addition, reads with a length shorter than 15 bp or greater than 45 bp were excluded from the analysis. Then, the remaining reads were mapped using bowtie23 to mouse reference genome and small RNA annotation databases, with no mismatches allowed. These were obtained from SPORTS1.1 as provided by the authors; including mm10 reference genome, miRNA from miRBase (release-21), rRNA curated sequences from National Center for Biotechnology Information (NCBI) Nucleotide, tRNA from GtRNAdb, piRNA from piRBase, other ncRNA from Ensembl (release-89), and Rfam (release-12.3).
SPORTS1.1 output files were used to generate normalized RPM counts for the different small RNA biotypes annotations per sample and fragment length. Small RNA biotype abundance per experimental group was calculated averaging the previous values.
Differentially expressed fragments analysis was carried out in R v4.3.1 with DESeq2 R package v1.40.224 using raw fragment counts. Only sequences with counts > = 50 in at least 4 biological replicates were considered. Thresholds to define differentially expressed fragments (DEFs) were set using Padj < 0.05 (using Benjamini-Hochberg FDR) and log2[fold change (FC)] > |1|. For a better estimation of log2(FC) we used shrunken log2(FC) with apeglm method25, which looks at the largest fold changes that are not due to low counts and uses these to inform a prior distribution, allowing to compare estimated log2(FC) across experiments. Heatmaps were made using VST-normalized expression values of DEFs using pheatmap R package v1.0.12.
Functional enrichment analysis with gene ontologies was conducted using g: GOSt from g: Profiler R package v 0.2.3 with predicted gene targets for miR-486; only significant functional pathways (Padj < 0.05) were considered. These targets were obtained from the TargetScanMouse database (v8)26, which contains information for predicted microRNA targets in mice. In this case, we extracted the predicted targets for the miR-486 conserved family.
We used the same approach and thresholds to reanalyze Tomar et al. mapped reads to obtain their set of DEFs. For Raad et al., we used DEFs as included in the study’s source data. Intersection between small RNA biotypes was restricted to annotations of each study’s DEFs, in which we evaluated whether their results included the same biotype annotations as in our analysis. Results were represented in a chord diagram made with the circlize R package v0.4.1527. We aggregated RPM values of each intersected tsRNA DEFs to compare their abundance between the respective experimental groups per study.
Statistical analysis
Data is presented as the mean value accompanied by their respective standard errors (SEM). Statistical evaluations were conducted with GraphPad Prism software (version 10) or R v4.3.1. Details on the specific statistical tests applied to each study are indicated in the legends accompanying the figures.
Data availability
The sncRNA-Seq datasets generated and analyzed in this research can be found in the NCBI GEO repository with the accession number GSE277621.
References
Sharma, U. et al. Biogenesis and function of tRNA fragments during sperm maturation and fertilization in mammals. Science 351, 391–396 (2016).
Gapp, K. et al. Implication of sperm RNAs in transgenerational inheritance of the effects of early trauma in mice. Nat. Neurosci. 17, 667–669 (2014).
Wang, Y. et al. Sperm MicroRNAs confer depression susceptibility to offspring. Sci. Adv. 7, eabd7605 (2021).
Chen, Q. et al. Sperm TsRNAs contribute to intergenerational inheritance of an acquired metabolic disorder. Science 351, 397–400 (2016).
Tomar, A. et al. Epigenetic inheritance of diet-induced and sperm-borne mitochondrial RNAs. Nature 630, 720–727 (2024).
Argaw-Denboba, A. et al. Paternal microbiome perturbations impact offspring fitness. Nature 629, 652–659 (2024).
de Barbosa, C. High-fat diet reprograms the epigenome of rat spermatozoa and transgenerationally affects metabolism of the offspring. Mol. Metab. 5, 184–197 (2016).
Grandjean, V. et al. RNA-mediated paternal heredity of diet-induced obesity and metabolic disorders. Sci. Rep. 5, 18193 (2015).
Raad, G. et al. Paternal multigenerational exposure to an obesogenic diet drives epigenetic predisposition to metabolic diseases in mice. Elife 10, e61736 (2021).
Zhang, Y. et al. Dnmt2 mediates intergenerational transmission of paternally acquired metabolic disorders through sperm small non-coding RNAs. Nat. Cell. Biol. 20, 535–540 (2018).
Sharma, U. et al. Small RNAs are trafficked from the epididymis to developing mammalian sperm. Dev. Cell. 46, 481–494e6 (2018).
Conine, C. C., Sun, F., Song, L., Rivera-Pérez, J. A. & Rando, O. J. Small RNAs gained during epididymal transit of sperm are essential for embryonic development in mice. Dev. Cell. 46, 470–480e3 (2018).
Brüning, J. C. & Fenselau, H. Integrative neurocircuits that control metabolism and food intake. Science 381, eabl7398 (2023).
Dietrich, M. O., Zimmer, M. R., Bober, J. & Horvath, T. L. Hypothalamic Agrp neurons drive stereotypic behaviors beyond feeding. Cell 160, 1222–1232 (2015).
Padilla, S. L. et al. AgRP to Kiss1 neuron signaling links nutritional state and fertility. Proc. Natl. Acad. Sci. USA. 114, 2413–2418 (2017).
Zimmer, M. R., Schmitz, A. E. & Dietrich, M. O. Activation of Agrp neurons modulates memory-related cognitive processes in mice. Pharmacol. Res. 141, 303–309 (2019).
Baver, S. B. et al. Leptin modulates the intrinsic excitability of AgRP/NPY neurons in the arcuate nucleus of the hypothalamus. J. Neurosci. 34, 5486–5496 (2014).
Beutler, L. R. et al. Obesity causes selective and long-lasting desensitization of AgRP neurons to dietary fat. Elife 9, e55909 (2020).
Zhang, S. et al. Molecular basis for selective activation of DREADD-based chemogenetics. Nature 612, 354–362 (2022).
Cai, J. et al. Angiogenin-mediated TsRNAs control inflammation and metabolic disorder by regulating NLRP3 inflammasome. Cell. Death Differ. 31, 1057–1069 (2024).
Rinaldi, V. D. et al. An atlas of cell types in the mouse epididymis and Vas deferens. Elife 9, e55474 (2020).
Shi, J., Ko, E. A., Sanders, K. M., Chen, Q. & Zhou, T. SPORTS1.0: A tool for annotating and profiling non-coding RNAs optimized for rRNA- and tRNA-derived small RNAs. Genomics Proteom. Bioinf. 16, 144–151 (2018).
Langmead, B., Trapnell, C., Pop, M. & Salzberg, S. L. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10, R25 (2009).
Love, M. I., Huber, W. & Anders, S. Moderated Estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Zhu, A., Ibrahim, J. G. & Love, M. I. Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences. Bioinformatics 35, 2084–2092 (2019).
Agarwal, V., Bell, G. W., Nam, J. W. & Bartel, D. P. Predicting effective microRNA target sites in mammalian mRNAs. Elife 4, e05005 (2015).
Gu, Z., Gu, L., Eils, R., Schlesner, M. & Brors, B. Circlize implements and enhances circular visualization in R. Bioinformatics 30, 2811–2812 (2014).
Acknowledgements
We thank Dr. Paula Mera and all the members of the Neuronal Control of Metabolism laboratory for invaluable discussions about this work. Additionally, we acknowledge Dr Anna Enjuanes and all the members of the Functional Genomics facility of IDIBAPS for help and technical support with small RNA sequencing. This work was carried out in part at the Esther Koplowitz Centre, Barcelona. The study received funding by the Instituto de Salud Carlos III under the project PI20/00732 (co-funded by the European Regional Development Fund “A Way to Make Europe”) (given to A.O). A.O. was awarded a Miguel Servet contract (CP19/00083) sponsored by the Instituto de Salud Carlos III and co-supported by the European Social Fund (FSE) under the “Investing in your future” initiative. I.T. is the recipient of the PRE2020-094896 fellowship, which is sponsored by the MCIN/AEI/https://doi.org/10.13039/501100011033/ and co-funded by the European Social Fund. S.R.L. is supported by the INVESTIGO program, grant 2022 INV-1 00028 (contract 100028TC3), which is funded by the European Union’s “NextGeneration EU” and the Generalitat de Catalonia. F.D-C. is a recipient of a fellowship from the Agencia Nacional de Investigación y Desarrollo (ANID), Doctorado Nacional 21,210,611, Chile and Scientific Exchange Grant 10,594 from the European Molecular Biology Organization (EMBO), Germany.
Author information
Authors and Affiliations
Contributions
I.T performed experiments, interpreted data, and edited the manuscript. S.R.L analysed data, interpreted data and edited the manuscript. F.C.D performed experiments. M.C edited the manuscript. A.O supervised the study, designed the study, performed experiments, interpreted data, secured funding and wrote the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing 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/.
About this article
Cite this article
Tahiri, I., Llana, S.R., Díaz-Castro, F. et al. AgRP neurons shape the sperm small RNA payload. Sci Rep 15, 7206 (2025). https://doi.org/10.1038/s41598-025-91391-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-025-91391-4