BACKGROUND: Biological aging is an important factor leading to the development of pathologies associated with metabolic dysregulation, including type 2 diabetes, cancer, cardiovascular and neurodegenerative diseases. Telomere length, a central feature of aging, has additionally been identified as inversely associated with glucose tolerance and the development of type 2 diabetes. However, the effects of shortened telomeres on body weight and metabolism remain incompletely understood. Here, we studied the metabolic consequences of moderate telomere shortening using second generation loss of telomerase activity in mice. RESULTS: Aged male and female G2 Terc-/- mice and controls were characterized with respect to body weight and composition, glucose homeostasis, insulin sensitivity and metabolic activity. This was complemented with molecular and histological analysis of adipose tissue, liver and the intestine as well as microbiota analysis. We show that moderate telomere shortening leads to improved insulin sensitivity and glucose tolerance in aged male and female G2 Terc-/- mice. This is accompanied by reduced fat and lean mass in both sexes. Mechanistically, the metabolic improvement results from reduced dietary lipid uptake in the intestine, characterized by reduced gene expression of fatty acid transporters in enterocytes of the small intestine. Furthermore, G2-Terc-/- mice showed significant alterations in the composition of gut microbiota, potentially contributing to the improved glucose metabolism. CONCLUSIONS: Our study shows that moderate telomere shortening reduces intestinal lipid absorption, resulting in reduced adiposity and improved glucose metabolism in aged mice. These findings will guide future murine and human aging studies and provide important insights into the age associated development of type 2 diabetes and metabolic syndrome.

BACKGROUND: Smell abilities differ greatly among vertebrate species due to distinct sensory needs, with exceptional variability reported in the number of olfactory genes and the size of the odour-processing regions of the brain. However, key environmental factors shaping genomic and phenotypic changes linked to the olfactory system remain difficult to identify at macroevolutionary scales. Here, we investigate the association between diverse ecological traits and the number of olfactory chemoreceptors in approximately two hundred ray-finned fishes. RESULTS: We found independent expansions producing large gene repertoires in several lineages of nocturnal amphibious fishes, generally able to perform active terrestrial exploration. We reinforced this finding with on-purpose genomic and transcriptomic analysis of Channallabes apus, a catfish species from a clade with chemosensory-based aerial orientation. Furthermore, we also detected an augmented information-processing capacity in the olfactory bulb of nocturnal amphibious fishes by estimating the number of cells contained in this brain region in twenty-four actinopterygian species. CONCLUSIONS: Overall, we report a convergent genomic and phenotypic magnification of the olfactory system in nocturnal amphibious fishes. This finding suggests the possibility of an analogous evolutionary event in fish-like tetrapod ancestors during the first steps of the water-to-land transition, favouring terrestrial adaptation through enhanced aerial orientation.

BACKGROUND: Individual functional modifications shape the ability of wildlife populations to cope with anthropogenic environmental changes. But instead of adaptive response, human-altered environments can generate a succession of deleterious functional changes leading to the extinction of the population. To study how persistent anthropogenic changes impacted local species' population status, we characterised population structure, genetic diversity and individual response of gene expression in the tree frog Hyla orientalis along a gradient of radioactive contamination around the Chernobyl nuclear power plant. RESULTS: We detected lower effective population size in populations most exposed to ionizing radiation in the Chernobyl Exclusion Zone that is not compensated by migrations from surrounding areas. We also highlighted a decreased body condition of frogs living in the most contaminated area, a distinctive transcriptomics signature and stop-gained mutations in genes involved in energy metabolism. While the association with dose will remain correlational until further experiments, a body of evidence suggests the direct or indirect involvement of radiation exposure in these changes. CONCLUSIONS: Despite ongoing migration and lower total dose rates absorbed than at the time of the accident, our results demonstrate that Hyla orientalis specimens living in the Chernobyl Exclusion Zone are still undergoing deleterious changes, emphasizing the long-term impacts of the nuclear disaster.

BACKGROUND: The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development. RESULTS: We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content. CONCLUSIONS: Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

BACKGROUND: Deceptive alternative mating tactics are commonly maintained at low frequencies within populations because males using them are less competitive and acquire lower fitness than those using dominant tactics. However, the successful invasion of a male deceptive tactic is plausible if deception carries no fitness cost to females. Among populations of the gift-giving spider Paratrechalea ornata, males very often offer females a deceptive worthless gift, rather than a nutritive gift. We tested the degree to which deceptive worthless gifts can occur in natural populations living under divergent environmental conditions (moderate and stressful). We examined the plasticity of morphological and behavioral traits and analyzed the fitness of females in relation to the gift type, also examining the paternity acquired by males offering either gift type. RESULTS: We demonstrated that worthless gifts can become dominant under highly stressful environmental conditions (84-100%). Individuals in such environment reach smaller sizes than those in moderate conditions. We suggest that the size reduction probably favors low metabolic demands in both sexes and may reduce the costs associated with receiving deceptive worthless gifts for females. In contrast, males living under moderate conditions varied the use of the deceptive tactic (0-95%), and worthless gifts negatively influenced female fecundity. Furthermore, male size, rather than gift content, positively impacted paternity success in the moderate but not in the stressful environment. CONCLUSIONS: Overall, this is the first empirical evidence that a reversible deceptive tactic can become dominant when the environment becomes harsh and mate choice becomes limited.

BACKGROUND: Pennisetum giganteum (AABB, 2n = 4x = 28) is a C4 plant in the genus Pennisetum with origin in Africa but currently also grown in Asia and America. It is a crucial forage and potential energy grass with significant advantages in yield, stress resistance, and environmental adaptation. However, the mechanisms underlying these advantageous traits remain largely unexplored. Here, we present a high-quality genome assembly of the allotetraploid P. giganteum aiming at providing insights into biomass accumulation. RESULTS: Our assembly has a genome size 2.03 Gb and contig N50 of 88.47 Mb that was further divided into A and B subgenomes. Genome evolution analysis revealed the evolutionary relationships across the Panicoideae subfamily lineages and identified numerous genome rearrangements that had occurred in P. giganteum. Comparative genomic analysis showed functional differentiation between the subgenomes. Transcriptome analysis found no subgenome dominance at the overall gene expression level; however, differentially expressed homoeologous genes and homoeolog-specific expressed genes between the two subgenomes were identified, suggesting that complementary effects between the A and B subgenomes contributed to biomass accumulation of P. giganteum. Besides, C4 photosynthesis-related genes were significantly expanded in P. giganteum and their sequences and expression patterns were highly conserved between the two subgenomes, implying that both subgenomes contributed greatly and almost equally to the highly efficient C4 photosynthesis in P. giganteum. We also identified key candidate genes in the C4 photosynthesis pathway that showed sustained high expression across all developmental stages of P. giganteum. CONCLUSIONS: Our study provides important genomic resources for elucidating the genetic basis of advantageous traits in polyploid species, and facilitates further functional genomics research and genetic improvement of P. giganteum.

BACKGROUND: Single-cell RNA sequencing (scRNA-seq) has revolutionized the transcriptomics field by advancing analyses from tissue-level to cell-level resolution. Despite the great advances in the development of computational methods for various steps of scRNA-seq analyses, one major bottleneck of the existing technologies remains in identifying the molecular relationship between disease phenotype and cell subpopulations, where "disease phenotype" refers to the clinical characteristics of each patient sample, and subpopulation refer to groups of single cells, which often do not correspond to clusters identified by standard single-cell clustering analysis. Here, we present PACSI, a method aimed at distinguishing cell subpopulations associated with disease phenotypes at the single-cell level. RESULTS: PACSI takes advantage of the topological properties of biological networks to introduce a proximity-based measure that quantifies the correlation between each cell and the disease phenotype of interest. Applied to simulated data and four case studies, PACSI accurately identified cells associated with disease phenotypes such as diagnosis, prognosis, and response to immunotherapy. In addition, we demonstrated that PACSI can also be applied to spatial transcriptomics data and successfully label spots that are associated with poor survival of breast carcinoma. CONCLUSIONS: PACSI is an efficient method to identify cell subpopulations associated with disease phenotypes. Our research shows that it has a broad range of applications in revealing mechanistic and clinical insights of diseases.

BACKGROUND: Conversion or editing of adenosine (A) into inosine (I) catalyzed by specialized cellular enzymes represents one of the most common post-transcriptional RNA modifications with emerging connection to disease. A-to-I conversions can happen at specific sites and lead to increase in proteome diversity and changes in RNA stability, splicing, and regulation. Such sites can be detected as adenine-to-guanine sequence changes by next-generation RNA sequencing which resulted in millions reported sites from multiple genome-wide surveys. Nonetheless, the lack of extensive independent validation in such endeavors, which is critical considering the relatively high error rate of next-generation sequencing, leads to lingering questions about the validity of the current compendiums of the editing sites and conclusions based on them. RESULTS: Strikingly, we found that the current analytical methods suffer from very high false positive rates and that a significant fraction of sites in the public databases cannot be validated. In this work, we present potential solutions to these problems and provide a comprehensive and extensively validated list of A-to-I editing sites in a human cancer cell line. Our findings demonstrate that most of true A-to-I editing sites in a human cancer cell line are located in the non-coding transcripts, the so-called RNA 'dark matter'. On the other hand, many ADAR editing events occurring in exons of human protein-coding mRNAs, including those that can recode the transcriptome, represent false positives and need to be interpreted with caution. Nonetheless, yet undiscovered authentic ADAR sites that increase the diversity of human proteome exist and warrant further identification. CONCLUSIONS: Accurate identification of human ADAR sites remains a challenging problem, particularly for the sites in exons of protein-coding mRNAs. As a result, genome-wide surveys of ADAR editome must still be accompanied by extensive Sanger validation efforts. However, given the vast number of unknown human ADAR sites, there is a need for further developments of the analytical techniques, potentially those that are based on deep learning solutions, in order to provide a quick and reliable identification of the editome in any sample.

BACKGROUND: Aphids are common crop pests. These insects reproduce by facultative parthenogenesis involving several rounds of clonal reproduction interspersed with an occasional sexual cycle. Furthermore, clonal aphids give birth to live young that are already pregnant. These qualities enable rapid population growth and have facilitated the colonisation of crops globally. In several cases, so-called "super clones" have come to dominate agricultural systems. However, the extent to which the sexual stage of the aphid life cycle has shaped global pest populations has remained unclear, as have the origins of successful lineages. Here, we used chromosome-scale genome assemblies to disentangle the evolution of two global pests of cereals-the English (Sitobion avenae) and Indian (Sitobion miscanthi) grain aphids. RESULTS: Genome-wide divergence between S. avenae and S. miscanthi is low. Moreover, comparison of haplotype-resolved assemblies revealed that the S. miscanthi isolate used for genome sequencing is likely a hybrid, with one of its diploid genome copies closely related to S. avenae (~ 0.5% divergence) and the other substantially more divergent (> 1%). Population genomics analyses of UK and China grain aphids showed that S. avenae and S. miscanthi are part of a cryptic species complex with many highly differentiated lineages that predate the origins of agriculture. The complex consists of hybrid lineages that display a tangled history of hybridisation and genetic introgression. CONCLUSIONS: Our analyses reveal that hybridisation has substantially contributed to grain aphid diversity, and hence, to the evolutionary potential of this important pest species. Furthermore, we propose that aphids are particularly well placed to exploit hybridisation events via the rapid propagation of live-born "frozen hybrids" via asexual reproduction, increasing the likelihood of hybrid lineage formation.

BACKGROUND: Neurotransmitter release depends on the fusion of synaptic vesicles with the presynaptic membrane and is mainly mediated by SNARE complex assembly. During the transition of Munc18-1/Syntaxin-1 to the SNARE complex, the opening of the Syntaxin-1 linker region catalyzed by Munc13-1 leads to the extension of the domain 3a hinge loop, which enables domain 3a to bind SNARE motifs in Synaptobrevin-2 and Syntaxin-1 and template the SNARE complex assembly. However, the exact mechanism of domain 3a extension remains elusive. RESULTS: Here, we characterized residues on the domain 3a hinge loop that are crucial for the extension of domain 3a by using biophysical and biochemical approaches and electrophysiological recordings. We showed that the mutation of residues T323/M324/R325 disrupted Munc13-1-mediated SNARE complex assembly and membrane fusion starting from Munc18-1/Syntaxin-1 in vitro and caused severe defects in the synaptic exocytosis of mouse cortex neurons in vivo. Moreover, the mutation had no effect on the binding of Synaptobrevin-2 to isolated Munc18-1 or the conformational change of the Syntaxin-1 linker region catalyzed by the Munc13-1 MUN domain. However, the extension of the domain 3a hinge loop in Munc18-1/Syntaxin-1 was completely disrupted by the mutation, leading to the failure of Synaptobrevin-2 binding to Munc18-1/Syntaxin-1. CONCLUSIONS: Together with previous results, our data further support the model that the template function of Munc18-1 in SNARE complex assembly requires the extension of domain 3a, and particular residues in the domain 3a hinge loop are crucial for the autoinhibitory release of domain 3a after the MUN domain opens the Syntaxin-1 linker region.

BACKGROUND: Prime editing (PE) is the most recent gene editing technology able to introduce targeted alterations to the genome, including single base pair changes, small insertions, and deletions. Several improvements to the PE machinery have been made in the past few years, and these have been tested in a range of model systems including immortalized cell lines, stem cells, and animal models. While double nicking RNA (dncRNA) PE systems PE3 and PE5 currently show the highest editing rates, they come with reduced accuracy as undesired indels or SNVs arise at edited loci. Here, we aimed to improve single ncRNA (sncRNA) systems PE2 and PE4max by generating novel all-in-one (pAIO) plasmids driven by an EF-1α promoter, which is especially suitable for human-induced pluripotent stem cell (hiPSC) models. RESULTS: pAIO-EF1α-PE2 and pAIO-EF1α-PE4max were used to edit the voltage gated potassium channel gene KCNQ2 and voltage gated sodium channel gene SCN1A. Two clinically relevant mutations were corrected using pAIO-EF1α-PE2 including the homozygous truncating SCN1A R612* variant in HEK293T cells and the heterozygous gain-of-function KCNQ2 R201C variant in patient-derived hiPSC. We show that sncRNA PE yielded detectable editing rates in hiPSC ranging between 6.4% and 9.8%, which was further increased to 41% after a GFP-based fluorescence-activated cell sorting (FACS) cell sorting step. Furthermore, we show that selecting the high GFP expressing population improved editing efficiencies up to 3.2-fold compared to the low GFP expressing population, demonstrating that not only delivery but also the number of copies of the PE enzyme and/or pegRNA per cell are important for efficient editing. Edit rates were not improved when an additional silent protospacer-adjacent motif (PAM)-removing alteration was introduced in hiPSC at the target locus. Finally, there were no genome-wide off-target effects using pAIO-EF1α-PE2 and no off-target editing activity near the edit locus highlighting the accuracy of snc prime editors. CONCLUSION: Taken together, our study shows an improved efficacy of EF-1α driven sncRNA pAIO-PE plasmids in hiPSC reaching high editing rates, especially after FACS sorting. Optimizing these sncRNA PE systems is of high value when considering future therapeutic in vivo use, where accuracy will be extremely important.

BACKGROUND: Adenine base editors (ABEs) are promising therapeutic gene editing tools that can efficiently convert targeted A•T to G•C base pairs in the genome. However, the large size of commonly used ABEs based on SpCas9 hinders its delivery in vivo using certain vectors such as adeno-associated virus (AAV) during preclinical applications. Despite a number of approaches having previously been attempted to overcome that challenge, including split Cas9-derived and numerous domain-deleted versions of editors, whether base editor (BE) and prime editor (PE) systems can also allow deletion of those domains remains to be proven. In this study, we present a new small ABE (sABE) with significantly reduced size. RESULTS: We discovered that ABE8e can tolerate large single deletions in the REC2 (Δ174-296) and HNH (Δ786-855) domains of SpCas9, and these deletions can be stacked together to create a new sABE. The sABE showed higher precision than the original ABE8e, with proximally shifted protospacer adjacent motif (PAM) editing windows (A3- A15), and comparable editing efficiencies to 8e-SaCas9-KKH. The sABE system efficiently generated A-G mutations at disease-relevant loci (T1214C in GAA and A494G in MFN2) in HEK293T cells and several canonical Pcsk9 splice sites in N2a cells. Moreover, the sABE enabled in vivo delivery in a single adeno-associated virus (AAV) vector with slight efficiency. Furthermore, we also successfully edited the genome of mouse embryos by microinjecting mRNA and sgRNA of sABE system into zygotes. CONCLUSIONS: We have developed a substantially smaller sABE system that expands the targeting scope and offers higher precision of genome editing. Our findings suggest that the sABE system holds great therapeutic potential in preclinical applications.

BACKGROUND: Rodents utilize chemical cues to recognize and avoid other conspecifics infected with pathogens. Infection with pathogens and acute inflammation alter the repertoire and signature of olfactory stimuli emitted by a sick individual. These cues are recognized by healthy conspecifics via the vomeronasal or accessory olfactory system, triggering an innate form of avoidance behavior. However, the molecular identity of the sensory neurons and the higher neural circuits involved in the detection of sick conspecifics remain poorly understood. RESULTS: We employed mice that are in an acute state of inflammation induced by systemic administration of lipopolysaccharide (LPS). Through conditional knockout of the G-protein Gαi2 and deletion of other key sensory transduction molecules (Trpc2 and a cluster of 16 vomeronasal type 1 receptors), in combination with behavioral testing, subcellular Ca2+ imaging, and pS6 and c-Fos neuronal activity mapping in freely behaving mice, we show that the Gαi2+ vomeronasal subsystem is required for the detection and avoidance of LPS-treated mice. The active components underlying this avoidance are contained in urine whereas feces extract and two selected bile acids, although detected in a Gαi2-dependent manner, failed to evoke avoidance behavior. Our analyses of dendritic Ca2+ responses in vomeronasal sensory neurons provide insight into the discrimination capabilities of these neurons for urine fractions from LPS-treated mice, and how this discrimination depends on Gαi2. We observed Gαi2-dependent stimulation of multiple brain areas including medial amygdala, ventromedial hypothalamus, and periaqueductal grey. We also identified the lateral habenula, a brain region implicated in negative reward prediction in aversive learning, as a previously unknown target involved in these tasks. CONCLUSIONS: Our physiological and behavioral analyses indicate that the sensing and avoidance of LPS-treated sick conspecifics depend on the Gαi2 vomeronasal subsystem. Our observations point to a central role of brain circuits downstream of the olfactory periphery and in the lateral habenula in the detection and avoidance of sick conspecifics, providing new insights into the neural substrates and circuit logic of the sensing of inflammation in mice.

BACKGROUND: Chronic kidney disease (CKD) accelerates atherosclerosis, but the mechanisms remain unclear. Tyrosine sulfation has been recognized as a key post-translational modification (PTM) in regulation of various cellular processes, and the sulfated adhesion molecules and chemokine receptors have been shown to participate in the pathogenesis of atherosclerosis via enhancement of monocyte/macrophage function. The levels of inorganic sulfate, the essential substrate for the sulfation reaction, are dramatically increased in patients with CKD, which indicates a change of sulfation status in CKD patients. Thus, in the present study, we detected the sulfation status in CKD patients and probed into the impact of sulfation on CKD-related atherosclerosis by targeting tyrosine sulfation function. RESULTS: PBMCs from individuals with CKD showed higher amounts of total sulfotyrosine and tyrosylprotein sulfotransferase (TPST) type 1 and 2 protein levels. The plasma level of O-sulfotyrosine, the metabolic end product of tyrosine sulfation, increased significantly in CKD patients. Statistically, O-sulfotyrosine and the coronary atherosclerosis severity SYNTAX score positively correlated. Mechanically, more sulfate-positive nucleated cells in peripheral blood and more abundant infiltration of sulfated macrophages in deteriorated vascular plaques in CKD ApoE null mice were noted. Knockout of TPST1 and TPST2 decreased atherosclerosis and peritoneal macrophage adherence and migration in CKD condition. The sulfation of the chemokine receptors, CCR2 and CCR5, was increased in PBMCs from CKD patients. CONCLUSIONS: CKD is associated with increased sulfation status. Increased sulfation contributes to monocyte/macrophage activation and might be involved in CKD-related atherosclerosis. Inhibition of sulfation may suppress CKD-related atherosclerosis and is worthy of further study.

BACKGROUND: Plants can perceive and respond to mechanical signals. For instance, cortical microtubule (CMT) arrays usually reorganize following the predicted maximal tensile stress orientation at the cell and tissue level. While research in the last few years has started to uncover some of the mechanisms mediating these responses, much remains to be discovered, including in most cases the actual nature of the mechanosensors. Such discovery is hampered by the absence of adequate quantification tools that allow the accurate and sensitive detection of phenotypes, along with high throughput and automated handling of large datasets that can be generated with recent imaging devices. RESULTS: Here we describe an image processing workflow specifically designed to quantify CMT arrays response to tensile stress in time-lapse datasets following an ablation in the epidermis - a simple and robust method to change mechanical stress pattern. Our Fiji-based workflow puts together several plugins and algorithms under the form of user-friendly macros that automate the analysis process and remove user bias in the quantification. One of the key aspects is also the implementation of a simple geometry-based proxy to estimate stress patterns around the ablation site and compare it with the actual CMT arrays orientation. Testing our workflow on well-established reporter lines and mutants revealed subtle differences in the response over time, as well as the possibility to uncouple the anisotropic and orientational response. CONCLUSION: This new workflow opens the way to dissect with unprecedented detail the mechanisms controlling microtubule arrays re-organization, and potentially uncover the still largely elusive plant mechanosensors.

BACKGROUND: The standard evolutionary theory of ageing proposes that ageing occurs because of a trade-off between reproduction and longevity. Eusocial insect queens exhibit positive fecundity-longevity associations and so have been suggested to be counter-examples through not expressing costs of reproduction and through remodelling conserved genetic and endocrine networks regulating ageing and reproduction. If so, eusocial evolution from solitary ancestors with negative fecundity-longevity associations must have involved a stage at which costs of reproduction were suppressed and fecundity and longevity became positively associated. Using the bumblebee (Bombus terrestris), we experimentally tested whether queens in annual eusocial insects at an intermediate level of eusocial complexity experience costs of reproduction, and, using mRNA-seq, the extent to which they exhibit a remodelling of relevant genetic and endocrine networks. Specifically, we tested whether costs of reproduction are present but latent, or whether a remodelling of relevant genetic and endocrine networks has already occurred allowing queens to reproduce without costs. RESULTS: We experimentally increased queens' costs of reproduction by removing their eggs, which caused queens to increase their egg-laying rate. Treatment queens had significantly reduced longevity relative to control queens whose egg-laying rate was not increased. Reduced longevity in treatment queens was not caused by increased worker-to-queen aggression or by increased overall activity in queens. In addition, treatment and control queens differed in age-related gene expression based on mRNA-seq in both their overall expression profiles and the expression of ageing-related genes. Remarkably, these differences appeared to occur principally with respect to relative age, not chronological age. CONCLUSIONS: This study represents the first simultaneously phenotypic and transcriptomic experimental test for a longevity cost of reproduction in eusocial insect queens. The results support the occurrence of costs of reproduction in annual eusocial insects of intermediate social complexity and suggest that reproductive costs are present but latent in queens of such species, i.e. that these queens exhibit condition-dependent positive fecundity-longevity associations. They also raise the possibility that a partial remodelling of genetic and endocrine networks underpinning ageing may have occurred in intermediately eusocial species such that, in unmanipulated conditions, age-related gene expression depends more on chronological than relative age.