Animals and housing<\/p>\n
Long-term laboratory colonies of naked mole-rats (H. glaber) were maintained in Berlin, Germany and Pretoria, South Africa. Damaraland mole-rat (F. damarensis) colonies were maintained in Pretoria, and Micklem\u2019s mole-rat (Fukomys micklemi), Mashona mole-rat (Fukomys darlingi) and Mechow\u2019s mole-rat (Fukomys mechowii) were maintained in \u010cesk\u00e9 Bud\u011bjovice, Czech Republic. (The following species were kept in the laboratory (Pretoria, South Africa) for unrelated studies after being collected from the wild: Highveld mole-rat (Cryptomys hottentotus pretoriae), Natal mole-rat (Cryptomys hottentotus natalensis), common mole-rat (Cryptomys hottentotus hottentotus), Cape mole-rat (G. capensis), Cape dune mole-rat (B. suillus), and Emin\u2019s mole-rat (H. emini).<\/p>\n
In Berlin, naked mole-rats were kept under controlled environmental conditions, with ambient temperature kept at 30\u201332\u2009\u00b0C and humidity at 50\u201370%, under dim illumination. Eighteen colonies were housed in custom-designed, interconnected plastic chamber systems (Fr\u00e4ntzel Kunststoffe). Animals were provided with a daily ad libitum diet of tubers (primarily sweet potatoes, celery root, cucumbers, bananas and carrots) and supplemented weekly with ProNutro (Bokomo). All husbandry and experimental procedures were approved by the local governmental authorities in Berlin (Landesamt f\u00fcr Gesundheit und Soziales, licenses G 0196\/17 and G 0121\/23). When selecting animals for experiments, a non-invasive salivary swab assay was established to determine sex by genotyping. In brief, saliva was collected from live naked mole-rats using a swab (Geyer). The extraction of genomic DNA was performed with proteinase K. Precipitation of the DNA was performed to concentrate the DNA solution used as template for the PCR. To determine whether the animal is a male, we amplified a 163-bp region of the sex-determining gene on the Y chromosome, Sry using Hg-SRY-For 5\u2032-GAAGAACGGCCATTTTTCGG-3\u2032 Hg-SRY-Rev 5\u2032-GCATTCATGGTGTGGTCTCG-3\u2032.<\/p>\n
To check the DNA quality and PCR condition, we amplified a 446-bp region of the mitochondrial 16S rRNA gene using Hg-16S-For, 5\u2032-TGGTGATAGCTGGTTGTCCA-3\u2032 and Hg-16S-Rev, 5\u2032-TAGTCTTTCCTTGCGGCACT-3\u2032, amplicon detectable in both sexes as previously described57<\/a>.<\/p>\n In South Africa, eight colonies of naked mole-rats were maintained in tunnel systems consisting of multiple plastic chambers (designated for food storage, toileting and nesting) connected by acrylic glass tunnels. Animals were fed a varied diet of chopped vegetables (primarily sweet potatoes, cucumbers and carrots) with weekly ProNutro supplementation. Nesting material consisted of wood shavings. The holding rooms were kept at 29\u201332\u2009\u00b0C with 50\u201370% relative humidity.<\/p>\n In the Czech Republic, the animal rooms of Micklem\u2019s mole-rat, Mashona mole-rat and Mechow\u2019s mole-rat were kept at 25\u2009\u00b1\u20091\u2009\u00b0C, 50\u2009\u00b1\u200910% relative humidity and photoperiod of 12\u2009h light:12\u2009h dark. Animals were fed ad libitum with vegetables (such as carrots, potatoes, sweet potatoes, beetroot, apple and cucumbers) and a rodent dry food mix. Animals were given the opportunity to carry out their natural digging behaviours in peat and were provided extra enrichment such as tree branches and plastic tubes for gnawing.<\/p>\n All procedures were approved by the relevant institutional animal care and use committees and local governmental authorities in Berlin (Max Delbr\u00fcck Center; Landesamt f\u00fcr Gesundheit und Soziales, licenses G 0196\/17 and G 0121\/23) and Pretoria (University of Pretoria), in accordance with national and international guidelines for the ethical treatment of laboratory animals. Procedures in Pretoria were approved by the Animal Ethics Committee of the University of Pretoria (license no. NAS199\/2020, NAS071-2023, NAS313\/2022, NAS324\/2022, NAS209-2021, NAS011\/2025 and NAS130_2025) and received DALRRD Section 20 approval (SEpi-Bizhub24110620262, SDAH-Epi-22101309360, SDAH-Epi-23032315040, SDAH-Epi-23041710040, 12\/11\/1\/1\/8 (2002 LH), SEpi-Bizhub25021311032 and SEpi-Bizhub25061009490). All naked mole-rats used in this study were derived from colonies originally captured by J. Jarvis, primarily in Mtito Andei and Lerata, Kenya, and represent a mixed parentage1<\/a>.<\/p>\n Behavioural assaysSocial rank determination<\/p>\n To assess dominance hierarchies within the eight colonies from South Africa, we followed a previously established assay for naked mole-rat social rank22<\/a>,58<\/a>. In brief, two plastic chambers were connected by a transparent tube, with one animal placed in each chamber. When both entered the tube simultaneously, the individual climbing over the other was scored as dominant. Each dyad was tested in at least three trials, with pairings pseudo-randomized and repeated across multiple months. A ranking index (R.I.) was calculated as wins divided by total trials and normalized to the colony maximum (R.I.\u2009=\u20091 for the top-ranking individual). Ranks were assigned categorically (rank 1: R.I.\u2009>\u20090.8, rank 5: R.I.\u2009<\u20090.2). For experiments, we used high-ranking individuals with R.I.\u2009\u2265\u20090.7 and low-ranking individuals with R.I.\u2009\u2264\u20090.3.<\/p>\n Social recognition and sniffing discrimination<\/p>\n Pairs of naked mole-rats were introduced into transparent acrylic tubes and video-recorded with a high-speed camera (BASLER, boA1936-400cm) at 100 frames per second for 5\u2009min. Sniffing bouts, characterized by rhythmic nose movements and defined as nose-to-body contact, were scored by an observer blind to colony identity. Each dot represents one animal-pair interaction.<\/p>\n Aggression assays<\/p>\n Animals from seven colonies from South Africa were paired and introduced into a rectangular arena (60\u2009cm\u2009\u00d7\u20097\u2009cm). Aggressive interactions\u2014including head-to-head lunges, shoves and bites\u2014were video-recorded, and head-to-head positions (distances <4\u2009cm) were analysed using a custom-written MATLAB script in combination with DeepLabCut59<\/a>. All frames were subsequently checked by human observers, and any errors were corrected to ensure accurate position identification. Recordings lasted 5\u2009min or were terminated earlier if there was a risk of severe injury. An aggression index was calculated as the proportion of time an animal engaged in aggressive behaviours relative to the total assay duration. Spatial distributions of head-to-head encounters were visualized as heat maps. Higher-ranking individuals, previously identified to defend their colonies in observations across these seven colonies, were used.<\/p>\n T-maze odour preference assay<\/p>\n Odour-guided attraction was tested in a two-choice T-tube olfactometer. The apparatus comprised a start chamber (30\u2009cm\u2009\u00d7\u200920\u2009cm\u2009\u00d7\u200914\u2009cm) connected to 2 entrance tunnels (64\u2009cm\u2009\u00d7\u20097\u2009cm), which extended into side tunnels (80\u2009cm\u2009\u00d7\u20097\u2009cm) leading to either a control or a stimulus chamber (30\u2009cm\u2009\u00d7\u200920\u2009\u00d7\u200914\u2009cm). In the food odour assay, one arm was scented with ProNutro (cereal-based supplement) food odour (odour source placed outside the maze to prevent contact), whereas the opposite arm contained no odour. In the IPM preference assay, bedding was treated with IPM diluted in DMSO (in total 50\u2009\u00b5l; 8.5\u2009\u00b5g) in one arm and 50\u2009\u00b5l DMSO as a solvent control in the other. Animals were placed in the start chamber and allowed to explore freely for 10\u2009min. Position and movement trajectories were recorded using a custom-written MATLAB script in combination with DeepLabCut59<\/a>. Attraction indices were calculated as (time in odour arm\u2009\u2013\u2009time in control arm)\/total assay time. To control for side bias, the odour-stimulus arm was alternated between replicates.<\/p>\n Maze exploration<\/p>\n Individuals were placed in a novel multi-arm maze (80\u2009cm\u2009\u00d7\u20097\u2009cm, 4 arms) and video-recorded for 15\u2009min. Animal trajectories were quantified and total coverage was calculated.<\/p>\n Pregnancy-suppression assay<\/p>\n Opposite-sex pairs were formed from low-ranking non-breeders, matched for body size and drawn from the same colonies to avoid confounds related to isolation history, and were housed in three chambers (30\u2009cm\u2009\u00d7\u200920\u2009cm\u2009\u00d7\u200914\u2009cm) connected by two tunnels (64\u2009cm\u2009\u00d7\u20097\u2009cm). Pairs were randomly assigned to blank bedding (n\u2009=\u20096), bedding supplemented with toilet bedding from their original colony (n\u2009=\u20096), or bedding treated daily with IPM (n\u2009=\u20097; 500\u2009\u00b5l applied to bedding, equivalent to 425\u2009mg; see Extended Data Fig. 3c,d<\/a> for time-dependent decay of IPM). Because testes in naked mole-rats remain intra-abdominal, external genital morphology was used only as a qualitative indicator and not as a quantitative measure of gonadal activation. Animals were monitored for 18 weeks, during which pregnancy status, body weight, genital morphology, vaginal perforation and reproductive activation were assessed; reproductive activation was assessed using endocrine measures and qualitative external genital observations, including vaginal perforation, which were documented photographically. Pregnancy was defined by progressive body mass gain of ~20\u201330\u2009g beginning at weeks 3\u20139 and continuing until weeks 15\u201318. In the blank group (3 replicates in Berlin, 3 in Pretoria), 5 of 6 females gained body mass and gave birth, whereas none of the females in the colony bedding (3 Berlin, 3 Pretoria) or IPM groups (3 Berlin\u00a0with 500 \u00b5l\u00a0isopropyl myristate + 500 \u00b5l\u00a0DMSO, 4 Pretoria\u00a0with 500 \u00b5l\u00a0neat isopropyl myristate;\u00a0both delivering ~425 mg IPM) became pregnant and body masses remained stable. Pregnancy outcomes were interpreted as a downstream consequence of reproductive activation and were evaluated alongside longitudinal endocrine measures and body mass trajectories. Males showed no significant body mass changes in any group.<\/p>\n Queen removal and colony stability<\/p>\n Three weeks after blood collection from an intact colony, the queen was removed from an intact colony that consisted of 19 individuals (8 females and 11 males) to initiate destabilization, after which IPM (500\u2009\u00b5l equal to 425\u2009mg; see Extended Data Fig. 3d<\/a> for time-dependent decay of IPM) was applied daily to the bedding for 12 weeks to simulate persistent queen-derived olfactory cues. The colony was inspected daily, and behavioural interactions were scored. Blood samples were collected at weeks 0, 15, 19 and 34 by tail venipuncture, and plasma was stored at \u201380\u2009\u00b0C until analysis. Prolactin concentrations were measured as described below. Body masses of all females were measured weekly using a precision balance (Scout Pro SPU123, Ohaus).<\/p>\n Odour collection and chemical analyses<\/p>\n Volatile compounds were collected by placing individual naked mole-rats in a Plexiglas chamber (30\u2009cm\u2009\u00d7\u200920\u2009cm\u2009\u00d7\u200914\u2009cm) for 30\u2009min, during which headspace odours were adsorbed onto PDMS tubes suspended from the chamber lid. PDMS tubing (1\u2009mm internal diameter\u2009\u00d7\u20090.4\u2009mm wall thickness; VWR International) was cut into 1-cm segments, soaked in 100% methanol for up to 24\u2009h, and conditioned under a constant flow of purified nitrogen at 180\u2009\u00b0C for 1.5\u2009h in a modified heating oven. Odour samples were collected from non-breeding males, non-breeding females, breeding males, and queens across reproductive states (pre-mating, mating, pregnancy and lactation (n\u2009=\u20093 for each)). Ovulation was inferred retrospectively using established reproductive staging in naked mole-rats, based on longitudinal behavioural observations, and defined reproductive states (pre-mating, mating, pregnancy, and lactation). As direct determination of the ovulatory window was not feasible, this classification reflects inferred reproductive state rather than precise measurement of ovulation. Additional samples were non-invasively obtained from eight African mole-rat species, and in naked mole-rat queens, swabs were taken from vaginal, anal, oral, teat, and skin regions using PDMS tubes on the same day.<\/p>\n Breeding females of the four Fukomys spp. were sampled across two reproductive stages (ovulation and pregnancy), as were naked mole-rats. The three Cryptomys spp. females were likewise sampled during ovulation and pregnancy, whereas Bathyergus females were sampled during lactation. We could not identify the reproductive stage of G. capensis. Importantly, IPM was detectable across all sampled reproductive stages (Fig. 2g<\/a>), but its abundance was most strongly elevated during ovulation, indicating that the observed species differences are not attributable to mismatched reproductive staging. Reproductive state in all breeding females was retrospectively determined using established morphological and behavioural criteria.<\/p>\n All samples were analysed by TD-GC\u2013MS as previously described23<\/a> with minor modifications. In brief, analyses were performed on an Agilent 7890A GC system coupled to a 5975C inert XL MSD and fitted with an HP5-MS UI column (19091S-433UI; Agilent Technologies). Following desorption at 250\u2009\u00b0C for 3\u2009min, volatiles were cryo-trapped at \u221250\u2009\u00b0C with liquid nitrogen and transferred to the GC column via a vaporizer injector heated to 270\u2009\u00b0C (12\u2009\u00b0C\u2009s\u22121, 5-min hold). The GC oven program was: 40\u2009\u00b0C for 3\u2009min, ramped at 5\u2009\u00b0C\u2009min\u22121 to 260\u2009\u00b0C (10-min hold), then 5\u2009\u00b0C\u2009min\u22121 to 280\u2009\u00b0C (5-min hold). Mass spectrometry parameters were: transfer line 260\u2009\u00b0C, ion source 230\u2009\u00b0C, quadrupole 150\u2009\u00b0C. Compounds were ionized by electron impact (70\u2009eV), detected in positive ion mode (m\/z 33\u2013500), and identified using NIST spectral libraries. XCMS (v3.7.1) was used to analyse and compare odour profiles. The identity of IPM was confirmed by comparison of mass spectra and retention time with a commercial standard (Sigma-Aldrich, 172472).<\/p>\n Calibration of IPM on PDMS tubes<\/p>\n To generate a calibration curve, defined amounts of IPM (1\u2009\u00b5l of 8.5\u2009ng, 850\u2009ng, 85\u2009\u00b5g, and 850\u2009\u00b5g; Sigma-Aldrich) were applied directly onto conditioned PDMS tubes (1\u2009cm\u2009\u00d7\u20090.3\u2009cm). After solvent evaporation, tubes were analysed using TD-GC\u2013MS (Agilent 7890\u2009A\/5975\u2009C, equipped with an HP-5MS UI column). Peak areas of the IPM chromatographic signal were quantified and log10-transformed. A linear regression of log10(peak area) versus log10(mass) was fitted (Y\u2009=\u20090.6016\u2009\u00d7\u2009X\u2009+\u20098.656). The mean peak area measured from queen samples was converted to mass equivalents using this calibration, corresponding to 660\u2009\u00b1\u200980\u2009ng IPM (mean\u2009\u00b1\u2009s.d., n\u2009=\u20094\u20137).<\/p>\n Time-course decay of IPM in a chamber<\/p>\n For breeding and queen-removal experiments, we examined the temporal dynamics of volatilized IPM, 500\u2009\u00b5l (425\u2009mg) of neat IPM was placed on the floor of a sealed Plexiglas chamber (30\u2009cm\u2009\u00d7\u200920\u2009cm\u2009\u00d7\u200914\u2009cm). PDMS tubes were suspended from the chamber lid and collected volatiles at defined timepoints (0, 6, 12, 18 and 24\u2009h). Tubes were analysed by TD-GC\u2013MS as described above. Peak areas were log10-transformed and plotted against time, revealing a progressive decline in IPM signal. Linear regression yielded Y\u2009=\u2009\u20130.05602\u2009\u00d7\u2009X\u2009+\u20098.827. The average queen IPM level corresponded to the concentration present at 5.2\u2009\u00b1\u20090.75\u2009h after application (n\u2009=\u20094).<\/p>\n Methimazole-induced olfactory ablation<\/p>\n Methimazole (Sigma-Aldrich; M8506) was diluted in 0.9% NaCl and administered by intraperitoneal injection at 75\u2009mg\u2009kg\u22121. Following methimazole injury, the olfactory epithelium in mice regenerates over approximately 1 month, with newly generated OSNs beginning to establish synaptic contacts with the olfactory bulb within ~1\u20132 weeks (ref. 60<\/a>). Accordingly, all behavioural experiments were conducted within 3\u20137 days post-injection, the established window of maximal OSN loss, ensuring effective olfactory ablation. Control animals received sterile saline only. After 72\u2009h, olfactory epithelium was collected for histological analysis as described18<\/a>. In brief, animals were perfused, and tissue was fixed in 4% paraformaldehyde (PFA), embedded in OCT, cryosectioned at 60\u2009\u00b5m, and immunostained with antibodies against OMP (mature OSNs; goat polyclonal, 1:1,000; Wako Chemicals, 544-10001) and STMN1 (immature OSNs; rabbit polyclonal, 1:500; Abcam, ab24445). Images were acquired using an Airyscan confocal fluorescence microscope (Axio, Zeiss) with a 20\u00d7 objective. Exposure parameters were kept constant across sections, and identical contrast adjustments were applied to all images. Mean fluorescence intensity was quantified using ImageJ61<\/a>.<\/p>\n Electro-olfactogram recordings<\/p>\n Recordings were performed as previously described62<\/a> with minor modifications. Animals were deeply anaesthetized with an overdose of ketamine and xylazine (100 and 20\u2009mg\u2009kg\u22121, respectively) and subsequently decapitated. Heads were bisected along the sagittal midline and immediately transferred to a dissecting microscope, where the olfactory epithelium was surgically exposed. Extracellular activity from OSNs was recorded using a tungsten wire electrode inserted into the olfactory mucosa. For stimulus delivery, 10\u2009\u03bcl of freshly prepared odorant solution was applied to a 1\u2009cm2 piece of filter paper (Whatman) and inserted into a glass Pasteur pipette. Odour stimuli (1\u2009s air puffs) of IPM (1:10 dilution in DMSO; equal to 85\u2009\u00b5g\u2009\u00b5l\u22121) or DMSO alone were delivered at a controlled flow rate of 40\u2009cm\u2009s\u22121 (Stimulus Controller CS 55, Ockenfels-Syntech). Signals were amplified at a sampling rate of 10\u2009kHz, low-pass filtered at 50\u2009Hz using a NeuroLog Amplifier (Digitimer), and recorded with a PowerLab 4\/30 system (ADInstruments) running LabChart 8 software for subsequent analysis.<\/p>\n Whole-brain FOS imaging<\/p>\n Animals were exposed to IPM (500\u2009\u00b5l, equivalent to 425\u2009mg) or solvent for 90\u2009min in a Plexiglas chamber (30\u2009cm\u2009\u00d7\u200920\u2009cm\u2009\u00d7\u200914\u2009cm). At the end of exposure time, animals were deeply anaesthetized with an intraperitoneal single-dose injection of ketamine and xylazine (100 and 20\u2009mg\u2009kg\u22121, respectively), transcardially perfused with phosphate-buffered saline (PBS) followed by 4% PFA. Whole-brain tissue was collected, incubated with 4% PFA overnight at 4\u2009\u00b0C and then transferred to 0.02% sodium azide in PBS for storage until clearing. One brain was excluded from the FOS analysis because unsuccessful perfusion compromised tissue quality and prevented reliable quantification. All other samples were processed according to the LifeCanvas Technologies protocol (v5.03, based on63<\/a>) in the following order: SHIELD preservation, Delipidation with SmartClear II Pro, Immunolabelling with SmartLabel. Cleared brains were incubated with 3.5\u2009\u03bcg of rabbit anti-FOS monoclonal antibody (Abcam, ab214672) together with 24\u2009\u03bcl of propidium iodide (Thermo Fisher, P3566, 1.0\u2009mg\u2009ml\u22121 solution in water) for nuclear stain. Primary antibody against FOS was fluorescently conjugated (Alexa Fluor\u00ae 647) with a donkey anti-rabbit secondary antibody (Jackson ImmunoResearch, 711-605-152) in 1:1.5 primary:secondary molar ratio. After active labelling, samples were matched to a refractive index of n\u2009=\u20091.52 using EasyIndex (LifeCanvas Technologies). Two datasets were generated, each containing equal numbers of animals from both treatment conditions. The first dataset (n\u2009=\u20094) was processed as described above and imaged with a custom built mesoSPIM29<\/a> light-sheet fluorescence microscope using 4\u00d7 magnification. Pixel size was 1.6\u2009\u00b5m for nominal lateral spatial resolution with a 6\u2009\u00b5m step size in the z-direction. 488, 561 and 647\u2009nm excitation was used with the respective 520\/50, 590\/50 and quad band 405\/488\/561\/640\u2009nm emission filters to image the entire brain for autofluorescence background, nuclear stain and FOS signal. The second dataset (n\u2009=\u20092) was sent to LifeCanvas Technologies where samples were passively delipidated for 7 days, batch labelled using SmartBatch+ with donkey anti-rabbit SeTau647 secondary antibody (LifeCanvas Technologies, DkxRb-ST), and imaged with SmartSPIM using a 3.6\u00d7 objective. Pixel size was 1.8\u2009\u03bcm laterally with a 4\u2009\u03bcm z-step. Laser excitation at 561\u2009nm and 640\u2009nm was used with 600\/50 and 690\/50 emission filters to image nuclear stain and FOS signal. All samples in each dataset were imaged using the same imaging protocol and excitation intensities.<\/p>\n Whole-brain image processing and analysis<\/p>\n FOS+ neurons in the olfactory bulb were quantified using a difference-of-Gaussian filter and local maxima detection implemented in Matlab 2023b (Mathworks). Whole-brain images were stitched and reconstructed using BigStitcher and Fiji64<\/a>,65<\/a>. All analyses were performed blinded to treatment conditions.<\/p>\n To enable quantitative comparison of FOS-positive cell densities acquired using different light-sheet microscopes, we assessed and matched the effective spatial resolution of mesoSPIM and SmartSPIM datasets. Intensity profiles along the axial (z) dimension were extracted from individual FOS-positive cells and fitted with Gaussian functions to estimate the axial point spread function (PSF). SmartSPIM data exhibited substantially higher axial resolution than mesoSPIM. To harmonize resolution, SmartSPIM image stacks were filtered with a 3D Gaussian kernel (\u03c3\u2009=\u20091, 1 and 2 pixels in x, y and z, respectively), preserving lateral resolution while increasing axial spread to match that of the mesoSPIM. Filtered SmartSPIM datasets were then resampled to the mesoSPIM voxel size prior to downstream cell counting and density analyses. Importantly, filtering parameters were derived solely from single-cell PSF measurements and were independent of cell counts.<\/p>\n The total volume of the olfactory bulb from which FOS-positive cells were quantified was estimated by manually thresholding the light-sheet fluorescence imaging data. The threshold was selected to include all visible olfactory bulb tissue while excluding background signal. The number of voxels exceeding this threshold was multiplied by the voxel volume to obtain the total sampled volume. FOS+ cell density was calculated as the number of FOS-positive cells divided by the total sampled olfactory bulb volume (cells per mm3).<\/p>\n fUSI experiment<\/p>\n Plane wave fUSI was performed using an Iconeus One ultrafast ultrasound system (Iconeus) equipped with a 15-MHz IcoPrime 4D MultiArray probe. Whole-brain volumes were acquired by scanning 40 coronal planes with a slice spacing of 0.21\u2009mm, using the manufacturer\u2019s default ultrafast plane-wave sequence. Experiments were performed in naked mole-rats aged 1\u20133 years with body weights not exceeding 40\u2009g. Animals were anaesthetized with ketamine (100\u2009mg\u2009kg\u22121) and xylazine (2\u2009mg\u2009kg\u22121) and secured in a stereotaxic frame. Rectal body temperature was maintained at 29.5\u2009\u00b0C using a heating pad placed beneath the animal. Odour delivery was controlled by a stimulus controller (CS 55, Ockenfels-Syntech) connected to a glass tube, from which gas was released at the distal end (2\u2009mm diameter). The glass tube was positioned in front of the anterior nostrils to allow efficient delivery of odorants into the nasal cavity. Each imaging session lasted 6\u2009min and followed an off\u2013on stimulation paradigm consisting of alternating 30-s stimulation and 30-s rest periods.<\/p>\n Signal analysis of fUSI data<\/p>\n Although no anatomical atlas is currently available for the adult naked mole-rat, we found that overall brain size and major vascular architecture resemble those of the mouse. Therefore, imaging data were registered to the Allen Mouse Brain Atlas66<\/a>. Following an initial anatomical registration using Iconeus\u2019 proprietary Brain Positioning System (BPS)67<\/a>, a second, animal-specific registration step was performed. For each animal, a randomly selected imaging session was chosen as a reference, and all remaining acquisitions from that animal were registered to this reference session, thereby inheriting the transformation matrix linking the data to atlas space. Registration was conducted in two stages\u2014coarse and fine\u2014using custom Python code implemented with SimpleITK registration utilities and methods68<\/a>.<\/p>\n First-level event-based analysis with generalized linear modelling<\/p>\n As previously described, olfactory stimulation was performed using a single event-based experimental design. For statistical analysis, we used previously developed custom code69<\/a> that adapts a generalized linear model framework70<\/a> from fMRI Python packages71<\/a> for fUSI data. Power Doppler images were first spatially smoothed using a Gaussian kernel (full width at half maximum, 300\u2009\u00b5m). A design matrix was then constructed incorporating the stimulus paradigm, cosine drift regressors, and data-driven confounding regressors. The stimulus pattern was convolved with a previously described hemodynamic response function72<\/a>. Cosine drift terms were computed up to 0.01\u2009Hz based on the power Doppler acquisition frequency (1\/6\u2009Hz). Data-driven confounds were estimated using aCompCor73<\/a> by extracting the first three principal components of the Power Doppler signal from white matter and cerebrospinal fluid voxels. Using a first-order autoregressive model74<\/a>, a generalized linear model was fitted for all imaging sessions in a given pair of animal and experimental condition (for IPM, DMSO (negative control), ethanol (positive control) and, when available, air olfactory stimulation). z-score maps were derived from model effect sizes using stimulus-specific contrast matrices. Statistical thresholds were determined using the Benjamini\u2013Hochberg procedure for false discovery rate correction (\u03b1\u2009=\u20090.01), and voxels not forming 6-connected clusters of at least 30 voxels above threshold were excluded. Finally, prominent clusters in thresholded IPM z-maps were identified, and spherical regions of interest (r\u2009\u2248\u2009400\u2009\u00b5m) were defined around the voxel exhibiting the maximal z-score.<\/p>\n rCBV in the piriform cortex<\/p>\n Following cluster detection in the z-score maps, one spherical mask per animal (n\u2009=\u20094) was selected to encompass a prominent activation cluster located within or adjacent to the piriform cortex. For each animal, power Doppler time series were extracted from all voxels within the corresponding mask for each experimental condition (IPM, DMSO and air). For each imaging session within a given condition, rCBV signals were obtained by regressing out previously computed data-driven aCompCor confounds from the extracted power Doppler time series using QR decomposition, followed by standardization of the residual signals as percent change relative to baseline. The baseline was defined as the first OFF period of the stimulation paradigm and the second half of subsequent OFF periods. For each rCBV time series, the average stimulus-evoked response was computed by windowing the stimulus ON period, padded at both onset and offset by half of the OFF period, and then averaging these windows across imaging sessions for each experimental condition.<\/p>\n Bootstrap and permutation testing<\/p>\n Permutation testing75<\/a> was performed separately for each animal to compare rCBV responses between pairs of experimental conditions. For each voxel within the spherical mask, the mean trapezoidal AUC of the rCBV signal was calculated over the stimulus ON period. The difference between mean AUC values for condition A (for example, IPM) and condition B (for example, DMSO) was then computed. Effect sizes and 99% confidence intervals were estimated using 5,000 permutation resamples, in which voxels were randomly reassigned to conditions A or B with replacement. Permutation testing was implemented using the DABEST Python package76<\/a>.<\/p>\n Plasma collection and prolactin measurement<\/p>\n Adult naked mole-rats were removed from their natal tunnel systems and weighed. Venous blood samples, corresponding to ~1% of the animal\u2019s body mass were collected from the tail vein using heparinized micro-haematocrit tubes and transferred into Eppendorf tubes. Samples were centrifuged immediately at 500g, and the resulting plasma was decanted and stored at \u221280\u2009\u00b0C until analysis. Plasma prolactin concentrations were quantified using a commercial enzyme-linked immunosorbent assay (ELISA; Elabscience, guinea pig prolactin kit, E-EL-GP0358) following the manufacturer\u2019s protocol. This assay has been previously validated for naked mole-rats38<\/a>. The intra-assay coefficient of variation was <10%, and assay sensitivity was 0.09\u2009ng\u2009ml\u22121.<\/p>\n Faecal sample collection<\/p>\n The opposite-sex pairs experiments were conducted across two locations (as stated above). Faecal hormone analyses were performed only on samples collected in South Africa (n\u2009=\u20093 pairs). Faecal samples were collected opportunistically during routine animal management procedures, including scheduled weighing and blood sampling events. Samples were collected immediately after defecation to avoid environmental contamination, transferred into sterile containers using sterilized forceps, and stored frozen until further processing.<\/p>\n Faecal steroid extraction<\/p>\n Faecal samples were freeze-dried and homogenized to a fine powder following established protocols42<\/a>. Depending on sample mass (0.0150\u20130.0249\u2009g, 0.0250\u20130.0366\u2009g and 0.0370\u20130.0550\u2009g), steroids were extracted using 0.5, 1.0 or 1.5\u2009ml of 80% ethanol, respectively. Samples were vortexed for 15\u2009min, centrifuged at 1,500g for 10\u2009min, and the resulting supernatant was transferred into microcentrifuge tubes77<\/a>. Extracts were stored at \u221220\u2009\u00b0C until hormone analysis.<\/p>\n Quantification of faecal androgen and progestagen metabolites<\/p>\n Faecal androgen metabolites and faecal progestagen metabolites were quantified using enzyme immunoassays (EIAs) for epiandrosterone78<\/a> and 5\u03b1-tetrahydroprogesterone79<\/a>, respectively. These assays have been validated for use in naked mole-rats and shown to reliably reflect reproductive endocrine activity80<\/a>. It should be noted that faecal steroid metabolites reflect circulating hormone changes with a physiological delay, as they represent metabolized hormone excretion rather than real-time endocrine levels.<\/p>\n Hormone quantification followed established protocols81<\/a>. Assay sensitivities were 6.0\u2009ng\u2009g\u22121 faecal dry mass for the 5\u03b1-tetrahydroprogesterone EIA and 7.2\u2009ng\u2009g\u22121 faecal dry mass for the epiandrosterone EIA. Intra-assay coefficients of variation were 4.38% and 5.32% for the 5\u03b1-tetrahydroprogesterone EIA and 7.40% and 7.94% for the epiandrosterone EIA. Inter-assay coefficients of variation for the 5\u03b1-tetrahydroprogesterone EIA were 11.10% and 12.75%. Only a single plate was used for the epiandrosterone EIA; therefore, inter-assay variation was not calculated.<\/p>\n Human odorant receptor screening and luminescence assay<\/p>\n HEK-293 cells (RRID:CVCL_0045<\/a>), a human embryonic kidney cell line, were used as a heterologous system for recombinant olfactory receptor expression. Cells were maintained at 37\u2009\u00b0C in a humidified incubator (5% CO2) in DMEM supplemented with 4.5\u2009g\u2009l\u22121 d-glucose, 10% fetal bovine serum, 2\u2009mM l-glutamine, 100\u2009U\u2009ml\u22121 penicillin and 100\u2009\u03bcg\u2009ml\u22121 streptomycin. Luminescence-based functional assays were performed as previously described35<\/a>. For receptor screening, a cDNA expression plasmid library comprising 766 constructs, encoding 386 human odorant receptor reference sequences and 380 common variants, was transfected into HEK-293 cells.<\/p>\n cAMP-dependent luminescence signals were acquired using a GloMax Discover detection system (Promega) and analysed in Excel (Microsoft, v16.109.2). Raw traces were averaged, and basal activity was subtracted from odorant-induced responses. For concentration\u2013response relationships, baseline-corrected values were normalized to the maximal response of a reference odorant receptor\u2013odorant pair, and signals from mock-transfected controls were subtracted. Dose\u2013response curves and half-maximal effective concentration (EC50) values were determined by nonlinear regression (SigmaPlot 14.0, Systat Software) using the equation f(x)\u2009=\u2009(min\u2009\u2013\u2009max)\/(1\u2009+\u2009(x\/EC50)Hillslope)+max, with Hill slope constrained to \u22642.9 and data points weighted by the reciprocal of their standard deviation (0.1\/s.d.). Dose\u2013response curves were fitted using nonlinear regression (log[agonist] versus response, variable slope).<\/p>\n Phylogenomic analysis<\/p>\n Protein-coding transcripts from ten African mole-rat species (H.\u00a0glaber, H. emini, G. capensis, F. damarensis, F. micklemi, F. darlingi, F. mechowii, C. hottentotus, C. natalensis and C. pretoriae) were combined with transcripts obtained from published data82<\/a> and from a newly assembled transcriptome of B. suillus. Following previous methodology82<\/a>, orthologous datasets were cleaned, resulting in a set of 3,393 single-copy genes (SCGs) shared across all 11 species. These genes were aligned using the mafft-linsi algorithm implemented in MAFFT83<\/a>, and alignments were further trimmed to remove columns containing gaps in more than 50% of the species. Orthofinder84<\/a> was used to discover orthogroups and build a species tree from the SCGs. Maximum-likelihood gene trees and species trees were inferred using default settings (-m MSA) based on the individual gene trees. For divergence time estimation, the species tree was calibrated with the TimeTree R package85<\/a>, constraining the split between M. musculus and Bathyergidae to 56\u201377 million years ago.<\/p>\n Statistics and figure preparation<\/p>\n Sample sizes were chosen on the basis of feasibility, ethical considerations and consistency with prior studies in naked mole-rats and other mole-rats, as well as established practice for behavioural, endocrine, imaging and electrophysiological experiments in this field. No statistical methods were used to predetermine sample size.\u00a0Where possible, multiple independent animals, colonies or experimental repetitions were used. Opposite-sex pairs in the pregnancy-suppression assay were randomly assigned to blank bedding, natal colony bedding or IPM. Other groups were defined by social\/reproductive status or colony manipulation, so randomization was not applicable. Behavioural and FOS analyses were performed and scored by observers blinded to treatment. Blinding was not applicable where groups were defined by social or reproductive status or by chemical or colony treatment. Normality was first assessed on datasets using a Shapiro\u2013Wilk test. Statistical analyses (see the corresponding legends of each figure) and figures were generated using GraphPad Prism v8 (https:\/\/www.graphpad.com<\/a>) and further processed with Adobe Illustrator v25.2.1.<\/p>\n Reporting summary<\/p>\n