c-fos induction in the choroid plexus, tanycytes and pars tuberalis is an early indicator of spontaneous arousal from torpor in a deep hibernator.

Hibernation is an extreme state of seasonal energy conservation, reducing metabolic rate to as little as 1% of the active state. During the hibernation season, many species of hibernating mammals cycle repeatedly between the active (aroused) and hibernating (torpid) states (T-A cycling), using brown adipose tissue (BAT) to drive cyclical rewarming. The regulatory mechanisms controlling this process remain undefined but are presumed to involve thermoregulatory centres in the hypothalamus. Here, we used the golden hamster (Mesocricetus auratus), and high-resolution monitoring of BAT, core body temperature and ventilation rate, to sample at precisely defined phases of the T-A cycle. Using c-fos as a marker of cellular activity, we show that although the dorsomedial hypothalamus is active during torpor entry, neither it nor the pre-optic area shows any significant changes during the earliest stages of spontaneous arousal. Contrastingly, in three non-neuronal sites previously linked to control of metabolic physiology over seasonal and daily time scales - the choroid plexus, pars tuberalis and third ventricle tanycytes - peak c-fos expression is seen at arousal initiation. We suggest that through their sensitivity to factors in the blood or cerebrospinal fluid, these sites may mediate metabolic feedback-based initiation of the spontaneous arousal process.

Circadian coupling of mitochondria in a deep-diving mammal.

Regulation of mitochondrial oxidative phosphorylation is essential to match energy supply to changing cellular energy demands, and to cope with periods of hypoxia. Recent work implicates the circadian molecular clock in control of mitochondrial function and hypoxia sensing. Because diving mammals experience intermittent episodes of severe hypoxia, with diel patterning in dive depth and duration, it is interesting to consider circadian-mitochondrial interaction in this group. Here, we demonstrate that the hooded seal (Cystophora cristata), a deep-diving Arctic pinniped, shows strong daily patterning of diving behaviour in the wild. Cultures of hooded seal skin fibroblasts exhibit robust circadian oscillation of the core clock genes per2 and arntl. In liver tissue collected from captive hooded seals, expression of arntl was some 4-fold higher in the middle of the night than in the middle of the day. To explore the clock-mitochondria relationship, we measured the mitochondrial oxygen consumption in synchronized hooded seal skin fibroblasts and found a circadian variation in mitochondrial activity, with higher coupling efficiency of complex I coinciding with the trough of arntl expression. These results open the way for further studies of circadian-hypoxia interactions in pinnipeds during diving.

Legacy genetics of Arachis cardenasii in the peanut crop shows the profound benefits of international seed exchange.

The narrow genetics of most crops is a fundamental vulnerability to food security. This makes wild crop relatives a strategic resource of genetic diversity that can be used for crop improvement and adaptation to new agricultural challenges. Here, we uncover the contribution of one wild species accession, Arachis cardenasii GKP 10017, to the peanut crop (Arachis hypogaea) that was initiated by complex hybridizations in the 1960s and propagated by international seed exchange. However, until this study, the global scale of the dispersal of genetic contributions from this wild accession had been obscured by the multiple germplasm transfers, breeding cycles, and unrecorded genetic mixing between lineages that had occurred over the years. By genetic analysis and pedigree research, we identified A. cardenasii-enhanced, disease-resistant cultivars in Africa, Asia, Oceania, and the Americas. These cultivars provide widespread improved food security and environmental and economic benefits. This study emphasizes the importance of wild species and collaborative networks of international expertise for crop improvement. However, it also highlights the consequences of the implementation of a patchwork of restrictive national laws and sea changes in attitudes regarding germplasm that followed in the wake of the Convention on Biological Diversity. Today, the botanical collections and multiple seed exchanges which enable benefits such as those revealed by this study are drastically reduced. The research reported here underscores the vital importance of ready access to germplasm in ensuring long-term world food security.

Biological timekeeping in polar environments: lessons from terrestrial vertebrates.

The polar regions receive less solar energy than anywhere else on Earth, with the greatest year-round variation in daily light exposure; this produces highly seasonal environments, with short summers and long, cold winters. Polar environments are also characterised by a reduced daily amplitude of solar illumination. This is obvious around the solstices, when the Sun remains continuously above (polar 'day') or below (polar 'night') the horizon. Even at the solstices, however, light levels and spectral composition vary on a diel basis. These features raise interesting questions about polar biological timekeeping from the perspectives of function and causal mechanism. Functionally, to what extent are evolutionary drivers for circadian timekeeping maintained in polar environments, and how does this depend on physiology and life history? Mechanistically, how does polar solar illumination affect core daily or seasonal timekeeping and light entrainment? In birds and mammals, answers to these questions diverge widely between species, depending on physiology and bioenergetic constraints. In the high Arctic, photic cues can maintain circadian synchrony in some species, even in the polar summer. Under these conditions, timer systems may be refined to exploit polar cues. In other instances, temporal organisation may cease to be dominated by the circadian clock. Although the drive for seasonal synchronisation is strong in polar species, reliance on innate long-term (circannual) timer mechanisms varies. This variation reflects differing year-round access to photic cues. Polar chronobiology is a productive area for exploring the adaptive evolution of daily and seasonal timekeeping, with many outstanding areas for further investigation.

Diversified regulation of circadian clock gene expression following whole genome duplication.

Across taxa, circadian control of physiology and behavior arises from cell-autonomous oscillations in gene expression, governed by a networks of so-called 'clock genes', collectively forming transcription-translation feedback loops. In modern vertebrates, these networks contain multiple copies of clock gene family members, which arose through whole genome duplication (WGD) events during evolutionary history. It remains unclear to what extent multiple copies of clock gene family members are functionally redundant or have allowed for functional diversification. We addressed this problem through an analysis of clock gene expression in the Atlantic salmon, a representative of the salmonids, a group which has undergone at least 4 rounds of WGD since the base of the vertebrate lineage, giving an unusually large complement of clock genes. By comparing expression patterns across multiple tissues, and during development, we present evidence for gene- and tissue-specific divergence in expression patterns, consistent with functional diversification of clock gene duplicates. In contrast to mammals, we found no evidence for coupling between cortisol and circadian gene expression, but cortisol mediated non-circadian regulated expression of a subset of clock genes in the salmon gill was evident. This regulation is linked to changes in gill function necessary for the transition from fresh- to sea-water in anadromous fish. Overall, this analysis emphasises the potential for a richly diversified clock gene network to serve a mixture of circadian and non-circadian functions in vertebrate groups with complex genomes.

A refined method to monitor arousal from hibernation in the European hamster.

BACKGROUND: Hibernation is a physiological and behavioural adaptation that permits survival during periods of reduced food availability and extreme environmental temperatures. This is achieved through cycles of metabolic depression and reduced body temperature (torpor) and rewarming (arousal). Rewarming from torpor is achieved through the activation of brown adipose tissue (BAT) associated with a rapid increase in ventilation frequency. Here, we studied the rate of rewarming in the European hamster (Cricetus cricetus) by measuring both BAT temperature, core body temperature and ventilation frequency. RESULTS: Temperature was monitored in parallel in the BAT (IPTT tags) and peritoneal cavity (iButtons) during hibernation torpor-arousal cycling. We found that increases in brown fat temperature preceded core body temperature rises by approximately 48 min, with a maximum re-warming rate of 20.9℃*h-1. Re-warming was accompanied by a significant increase in ventilation frequency. The rate of rewarming was slowed by the presence of a spontaneous thoracic mass in one of our animals. Core body temperature re-warming was reduced by 6.2℃*h-1 and BAT rewarming by 12℃*h-1. Ventilation frequency was increased by 77% during re-warming in the affected animal compared to a healthy animal. Inspection of the position and size of the mass indicated it was obstructing the lungs and heart. CONCLUSIONS: We have used a minimally invasive method to monitor BAT temperature during arousal from hibernation illustrating BAT re-warming significantly precedes core body temperature re-warming, informing future study design on arousal from hibernation. We also showed compromised re-warming from hibernation in an animal with a mass obstructing the lungs and heart, likely leading to inefficient ventilation and circulation.

The pars tuberalis: The site of the circannual clock in mammals?

Accurate timing and physiological adaptation to anticipate seasonal changes are an essential requirement for an organism's survival. In contrast to all other environmental cues, photoperiod offers a highly predictive signal that can be reliably used to activate a seasonal adaptive programme at the correct time of year. Coupled to photoperiod sensing, it is apparent that many organisms have evolved innate long-term timekeeping systems, allowing reliable anticipation of forthcoming environmental changes. The fundamental biological processes giving rise to innate long-term timing, with which the photoperiod-sensing pathway engages, are not known for any organism. There is growing evidence that the pars tuberalis (PT) of the pituitary, which acts as a primary transducer of photoperiodic input, may be the site of the innate long-term timer or "circannual clock". Current research has led to the proposition that the PT-specific thyrotroph may act as a seasonal calendar cell, driving both hypothalamic and pituitary endocrine circuits. Based on this research we propose that the mechanistic basis for the circannual rhythm appears to be deeply conserved, driven by a binary switching cell based accumulator, analogous to that proposed for development. We review the apparent conservation of function and pathways to suggest that these broad principles may apply across the vertebrate lineage and even share characteristics with processes driving seasonal adaptation in plants.

The Digital Ageing Atlas: integrating the diversity of age-related changes into a unified resource.

Multiple studies characterizing the human ageing phenotype have been conducted for decades. However, there is no centralized resource in which data on multiple age-related changes are collated. Currently, researchers must consult several sources, including primary publications, in order to obtain age-related data at various levels. To address this and facilitate integrative, system-level studies of ageing we developed the Digital Ageing Atlas (DAA). The DAA is a one-stop collection of human age-related data covering different biological levels (molecular, cellular, physiological, psychological and pathological) that is freely available online (http://ageing-map.org/). Each of the >3000 age-related changes is associated with a specific tissue and has its own page displaying a variety of information, including at least one reference. Age-related changes can also be linked to each other in hierarchical trees to represent different types of relationships. In addition, we developed an intuitive and user-friendly interface that allows searching, browsing and retrieving information in an integrated and interactive fashion. Overall, the DAA offers a new approach to systemizing ageing resources, providing a manually-curated and readily accessible source of age-related changes.

Dissecting the gene network of dietary restriction to identify evolutionarily conserved pathways and new functional genes.

Dietary restriction (DR), limiting nutrient intake from diet without causing malnutrition, delays the aging process and extends lifespan in multiple organisms. The conserved life-extending effect of DR suggests the involvement of fundamental mechanisms, although these remain a subject of debate. To help decipher the life-extending mechanisms of DR, we first compiled a list of genes that if genetically altered disrupt or prevent the life-extending effects of DR. We called these DR-essential genes and identified more than 100 in model organisms such as yeast, worms, flies, and mice. In order for other researchers to benefit from this first curated list of genes essential for DR, we established an online database called GenDR (http://genomics.senescence.info/diet/). To dissect the interactions of DR-essential genes and discover the underlying lifespan-extending mechanisms, we then used a variety of network and systems biology approaches to analyze the gene network of DR. We show that DR-essential genes are more conserved at the molecular level and have more molecular interactions than expected by chance. Furthermore, we employed a guilt-by-association method to predict novel DR-essential genes. In budding yeast, we predicted nine genes related to vacuolar functions; we show experimentally that mutations deleting eight of those genes prevent the life-extending effects of DR. Three of these mutants (OPT2, FRE6, and RCR2) had extended lifespan under ad libitum, indicating that the lack of further longevity under DR is not caused by a general compromise of fitness. These results demonstrate how network analyses of DR using GenDR can be used to make phenotypically relevant predictions. Moreover, gene-regulatory circuits reveal that the DR-induced transcriptional signature in yeast involves nutrient-sensing, stress responses and meiotic transcription factors. Finally, comparing the influence of gene expression changes during DR on the interactomes of multiple organisms led us to suggest that DR commonly suppresses translation, while stimulating an ancient reproduction-related process.

Genome-environment interactions that modulate aging: powerful targets for drug discovery.

Aging is the major biomedical challenge of this century. The percentage of elderly people, and consequently the incidence of age-related diseases such as heart disease, cancer, and neurodegenerative diseases, is projected to increase considerably in the coming decades. Findings from model organisms have revealed that aging is a surprisingly plastic process that can be manipulated by both genetic and environmental factors. Here we review a broad range of findings in model organisms, from environmental to genetic manipulations of aging, with a focus on those with underlying gene-environment interactions with potential for drug discovery and development. One well-studied dietary manipulation of aging is caloric restriction, which consists of restricting the food intake of organisms without triggering malnutrition and has been shown to retard aging in model organisms. Caloric restriction is already being used as a paradigm for developing compounds that mimic its life-extension effects and might therefore have therapeutic value. The potential for further advances in this field is immense; hundreds of genes in several pathways have recently emerged as regulators of aging and caloric restriction in model organisms. Some of these genes, such as IGF1R and FOXO3, have also been associated with human longevity in genetic association studies. The parallel emergence of network approaches offers prospects to develop multitarget drugs and combinatorial therapies. Understanding how the environment modulates aging-related genes may lead to human applications and disease therapies through diet, lifestyle, or pharmacological interventions. Unlocking the capacity to manipulate human aging would result in unprecedented health benefits.

Hypothalamic tanycytes as mediators of maternally programmed seasonal plasticity.

In mammals, maternal photoperiodic programming (MPP) provides a means whereby juvenile development can be matched to forthcoming seasonal environmental conditions.1,2,3,4 This phenomenon is driven by in utero effects of maternal melatonin5,6,7 on the production of thyrotropin (TSH) in the fetal pars tuberalis (PT) and consequent TSH receptor-mediated effects on tanycytes lining the 3rd ventricle of the mediobasal hypothalamus (MBH).8,9,10 Here we use LASER capture microdissection and transcriptomic profiling to show that TSH-dependent MPP controls the attributes of the ependymal region of the MBH in juvenile animals. In Siberian hamster pups gestated and raised on a long photoperiod (LP) and thereby committed to a fast trajectory for growth and reproductive maturation, the ependymal region is enriched for tanycytes bearing sensory cilia and receptors implicated in metabolic sensing. Contrastingly, in pups gestated and raised on short photoperiod (SP) and therefore following an over-wintering developmental trajectory with delayed sexual maturation, the ependymal region has fewer sensory tanycytes. Post-weaning transfer of SP-gestated pups to an intermediate photoperiod (IP), which accelerates reproductive maturation, results in a pronounced shift toward a ciliated tanycytic profile and formation of tanycytic processes. We suggest that tanycytic plasticity constitutes a mechanism to tailor metabolic development for extended survival in variable overwintering environments.

Serotonin: from top to bottom.

Serotonin is a monoamine neurotransmitter, which is phylogenetically conserved in a wide range of species from nematodes to humans. In mammals, age-related changes in serotonin systems are known risk factors of age-related diseases, such as diabetes, faecal incontinence and cardiovascular diseases. A decline in serotonin function with aging would be consistent with observations of age-related changes in behaviours, such as sleep, sexual behaviour and mood all of which are linked to serotonergic function. Despite this little is known about serotonin in relation to aging. This review aims to give a comprehensive analysis of the distribution, function and interactions of serotonin in the brain; gastrointestinal tract; skeletal; vascular and immune systems. It also aims to demonstrate how the function of serotonin is linked to aging and disease pathology in these systems. The regulation of serotonin via microRNAs is also discussed, as are possible applications of serotonergic drugs in aging research and age-related diseases. Furthermore, this review demonstrates that serotonin is potentially involved in whole organism aging through its links with multiple organs, the immune system and microRNA regulation. Methods to investigate these links are discussed.

A-to-I RNA editing does not change with age in the healthy male rat brain.

RNA editing is a post-transcriptional process, which results in base substitution modifications to RNA. It is an important process in generating protein diversity through amino acid substitution and the modulation of splicing events. Previous studies have suggested a link between gene-specific reductions in adenosine to inosine RNA editing and aging in the human brain. Here we demonstrate that changes in RNA editing observed in humans with age are not observed during aging in healthy rats. Furthermore, we identify a conserved editing site in rats, in Cog3. We propose that either age-related changes in RNA editing are specific to primates or humans, or that they are the manifestation of disease pathology. Since rodents are often used as model organisms for studying aging, these findings demonstrate the importance of understanding species-specific differences in RNA biology during aging.

GeneFriends: an online co-expression analysis tool to identify novel gene targets for aging and complex diseases.

BACKGROUND: Although many diseases have been well characterized at the molecular level, the underlying mechanisms are often unknown. Nearly half of all human genes remain poorly studied, yet these genes may contribute to a number of disease processes. Genes involved in common biological processes and diseases are often co-expressed. Using known disease-associated genes in a co-expression analysis may help identify and prioritize novel candidate genes for further study. RESULTS: We have created an online tool, called GeneFriends, which identifies co-expressed genes in over 1,000 mouse microarray datasets. GeneFriends can be used to assign putative functions to poorly studied genes. Using a seed list of disease-associated genes and a guilt-by-association method, GeneFriends allows users to quickly identify novel genes and transcription factors associated with a disease or process. We tested GeneFriends using seed lists for aging, cancer, and mitochondrial complex I disease. We identified several candidate genes that have previously been predicted as relevant targets. Some of the genes identified are already being tested in clinical trials, indicating the effectiveness of this approach. Co-expressed transcription factors were investigated, identifying C/ebp genes as candidate regulators of aging. Furthermore, several novel candidate genes, that may be suitable for experimental or clinical follow-up, were identified. Two of the novel candidates of unknown function that were co-expressed with cancer-associated genes were selected for experimental validation. Knock-down of their human homologs (C1ORF112 and C12ORF48) in HeLa cells slowed growth, indicating that these genes of unknown function, identified by GeneFriends, may be involved in cancer. CONCLUSIONS: GeneFriends is a resource for biologists to identify and prioritize novel candidate genes involved in biological processes and complex diseases. It is an intuitive online resource that will help drive experimentation. GeneFriends is available online at: http://genefriends.org/.

Gathering insights on disease etiology from gene expression profiles of healthy tissues.

MOTIVATION: Gene expression profiles have been widely used to study disease states. It may be possible, however, to gather insights into human diseases by comparing gene expression profiles of healthy organs with different disease incidence or severity. We tested this hypothesis and developed an approach to identify candidate genes associated with disease development by focusing on cancer incidence since it varies greatly across human organs. RESULTS: We normalized organ-specific cancer incidence by organ weight and found that reproductive organs tend to have a higher mass-normalized cancer incidence, which could be due to evolutionary trade-offs. Next, we performed a genome-wide scan to identify genes whose expression across healthy organs correlates with organ-specific cancer incidence. We identified a large number of genes, including genes previously associated with tumorigenesis and new candidate genes. Most genes exhibiting a positive correlation with cancer incidence were related to ribosomal and transcriptional activity, translation and protein synthesis. Organs with enhanced transcriptional and translational activation may have higher cell proliferation and therefore be more likely to develop cancer. Furthermore, we found that organs with lower cancer incidence tend to express lower levels of known cancer-associated genes. Overall, these results demonstrate how genes and processes that predispose organs to specific diseases can be identified using gene expression profiles from healthy tissues. Our approach can be applied to other diseases and serve as foundation for further oncogenomic analyses. CONTACT: jp@senescence.info SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Immunologic Profiling of the Atlantic Salmon Gill by Single Nuclei Transcriptomics.

Anadromous salmonids begin life adapted to the freshwater environments of their natal streams before a developmental transition, known as smoltification, transforms them into marine-adapted fish. In the wild, smoltification is a photoperiod-regulated process, involving radical remodeling of gill function to cope with the profound osmotic and immunological challenges of seawater (SW) migration. While prior work has highlighted the role of specialized "mitochondrion-rich" cells (MRCs) and accessory cells (ACs) in delivering this phenotype, recent RNA profiling experiments suggest that remodeling is far more extensive than previously appreciated. Here, we use single-nuclei RNAseq to characterize the extent of cytological changes in the gill of Atlantic salmon during smoltification and SW transfer. We identify 20 distinct cell clusters, including known, but also novel gill cell types. These data allow us to isolate cluster-specific, smoltification-associated changes in gene expression and to describe how the cellular make-up of the gill changes through smoltification. As expected, we noted an increase in the proportion of seawater mitochondrion-rich cells, however, we also identify previously unknown reduction of several immune-related cell types. Overall, our results provide fresh detail of the cellular complexity in the gill and suggest that smoltification triggers unexpected immune reprogramming.

Genome-wide association analysis of canine atopic dermatitis and identification of disease related SNPs.

In humans, genome-wide association studies (GWAS) have been shown to be an effective and thorough approach for identifying polymorphisms associated with disease phenotypes. Here, we describe the first study to perform a genome-wide association study in canine atopic dermatitis (cAD) using the Illumina Canine SNP20 array, containing 22,362 single-nucleotide polymorphisms (SNPs). The aim of the study was to identify SNPs associated with cAD using affected and unaffected Golden Retrievers. Further validation studies were performed for potentially associated SNPs using Sequenom genotyping of larger numbers of cases and controls across eight breeds (Boxer, German Shepherd Dog, Labrador, Golden Retriever, Shiba Inu, Shih Tzu, Pit Bull, and West Highland White Terriers). Using meta-analysis, two SNPs were associated with cAD in all breeds tested. RS22114085 was identified as a susceptibility locus (p=0.00014, odds ratio=2) and RS23472497 as a protective locus (p=0.0015, odds ratio=0.6). Both of these SNPs were located in intergenic regions, and their effects have been demonstrated to be independent of each other, highlighting that further fine mapping and resequencing is required of these areas. Further, 12 SNPs were validated by Sequenom genotyping as associated with cAD, but these were not associated with all breeds. This study suggests that GWAS will be a useful approach for identifying genetic risk factors for cAD. Given the clinical heterogeneity within this condition and the likelihood that the relative genetic effect sizes are small, greater sample sizes and further studies will be required.

An integrative view of mammalian seasonal neuroendocrinology.

Seasonal neuroendocrine cycles that govern annual changes in reproductive activity, energy metabolism and hair growth are almost ubiquitous in mammals that have evolved at temperate and polar latitudes. Changes in nocturnal melatonin secretion regulating gene expression in the pars tuberalis (PT) of the pituitary stalk are a critical common feature in seasonal mammals. The PT sends signal(s) to the pars distalis of the pituitary to regulate prolactin secretion and thus the annual moult cycle. The PT also signals in a retrograde manner via thyroid-stimulating hormone to tanycytes, which line the ventral wall of the third ventricle in the hypothalamus. Tanycytes show seasonal plasticity in gene expression and play a pivotal role in regulating local thyroid hormone (TH) availability. Within the mediobasal hypothalamus, the cellular and molecular targets of TH remain elusive. However, two populations of hypothalamic neurones, which produce the RF-amide neuropeptides kisspeptin and RFRP3 (RF-amide related peptide 3), are plausible relays between TH and the gonadotrophin-releasing hormone-pituitary-gonadal axis. By contrast, the ways by which TH also impinges on hypothalamic systems regulating energy intake and expenditure remain unknown. Here, we review the neuroendocrine underpinnings of seasonality and identify several areas that warrant further research.

Maternal Photoperiodic Programming: Melatonin and Seasonal Synchronization Before Birth.

This mini-review considers the phenomenon of maternal photoperiodic programming (MPP). In order to match neonatal development to environmental conditions at the time of birth, mammals use melatonin produced by the maternal pineal gland as a transplacental signal representing ambient photoperiod. Melatonin acts via receptors in the fetal pituitary gland, exerting actions on the developing medio-basal hypothalamus. Within this structure, a central role for specialized ependymal cells known as tanycytes has emerged, linking melatonin to control of hypothalamic thyroid metabolism and in turn to pup development. This review summarizes current knowledge of this programming mechanism, and its relevance in an eco-evolutionary context. Maternal photoperiodic programming emerges as a useful paradigm for understanding how in utero programing of hypothalamic function leads to life-long effects on growth, reproduction, health and disease in mammals, including humans.

Reference genes for canine skin when using quantitative real-time PCR.

Quantitative real-time PCR (qPCR) facilitates the quantification of mRNA expression. Accurate qPCR analysis of gene expression requires the normalisation of data using a reference or housekeeping gene which is expressed at a similar level in all tissues tested. GAPDH is the most well known and most widely used reference gene but many papers have demonstrated that it is not stably expressed in different tissues. The aim of this study was to measure reference gene stability in canine skin using real-time qPCR. Skin samples from healthy control dogs (n=7) and dogs with atopic dermatitis (lesional skin n=7 and non-lesional skin n=7) were used to quantify seven reference genes (IMP, CG14980, S7, HIRA, GAPDH, RPL13A and SDHA) in canine whole skin. Three different statistical programs (Bestkeeper, GeNorm and Normfinder) were used to assess the stability of the reference genes. The results confirmed that GAPDH is not a stably expressed reference gene in canine skin; this finding may influence interpretation of previous qPCR studies on canine skin using this as a reference gene. RPL13A and CG14980 were found to be the most stably expressed genes in canine whole skin and would be more suitable as reference genes in future studies.