What makes a mimic? Orange, red, and black color production in the mimic poison frog (Ranitomeya imitator).

Aposematic organisms rely on their conspicuous appearance to signal that they are defended and unpalatable. Such phenotypes are strongly tied to survival and reproduction. Aposematic colors and patterns are highly variable; however, the genetic, biochemical and physiological mechanisms producing this conspicuous coloration remain largely unidentified. Here, we identify genes potentially affecting color variation in two color morphs of Ranitomeya imitator: the orange-banded Sauce and the redheaded Varadero morphs. We examine gene expression in black and orange skin patches from the Sauce morph and black and red skin patches from the Varadero morph. We identified genes differentially expressed between skin color patches, including those that are involved in melanin synthesis (e.g., mlana, pmel, tyrp1), iridophore development (e.g., paics, ppat, ak1), pteridine synthesis (e.g., gch1, recql4, xdh), and carotenoid metabolism (e.g., dgat2, rbp1, scarb2). In addition, using weighted gene network analysis, we identified the top 50 genes with high connectivity from the most significant network associated with gene expression differences between color morphs. Of these 50 genes, 14 were known to be related to color production (gch1, gmps, gpr143, impdh1, mc1r, pax3-a, pax7, ppat, rab27a, rlbp1, tfec, trpm1, xdh).

Impacts of dietary fat on multi tissue gene expression in the desert-adapted cactus mouse.

Understanding the relationship between dietary fat and physiological responses is crucial in species adapted to arid environments where water scarcity is common. In this study, we present a comprehensive exploration of gene expression across five tissues (kidney, liver, lung, gastrointestinal tract, and hypothalamus) and 19 phenotypic measurements, investigating the effects of dietary fat in the desert-adapted cactus mouse ( Peromyscus eremicus ). We show impacts on immune function, circadian gene regulation, and mitochondrial function for mice fed a lower-fat diet compared to mice fed a higher-fat diet. In arid environments with severe water scarcity, even subtle changes in organismal health and water balance can affect physical performance, potentially impacting survival and reproductive success. The study sheds light on the complex interplay between diet, physiological processes, and environmental adaptation, providing valuable insights into the multifaceted impacts of dietary choices on organismal well-being and adaptation strategies in arid habitats.

When the tap runs dry: The multi-tissue gene expression and physiological responses of water deprived Peromyscus eremicus.

The harsh and dry conditions of desert environments have resulted in genomic adaptations, allowing for desert organisms to withstand prolonged drought, extreme temperatures, and limited food resources. Here, we present a comprehensive exploration of gene expression across five tissues (kidney, liver, lung, gastrointestinal tract, and hypothalamus) and 19 phenotypic measurements to explore the whole-organism physiological and genomic response to water deprivation in the desert-adapted cactus mouse (Peromyscus eremicus). The findings encompass the identification of differentially expressed genes and correlative analysis between phenotypes and gene expression patterns across multiple tissues. Specifically, we found robust activation of the vasopressin renin-angiotensin-aldosterone system (RAAS) pathways, whose primary function is to manage water and solute balance. Animals reduce food intake during water deprivation, and upregulation of PCK1 highlights the adaptive response to reduced oral intake via its actions aimed at maintained serum glucose levels. Even with such responses to maintain water balance, hemoconcentration still occurred, prompting a protective downregulation of genes responsible for the production of clotting factors while simultaneously enhancing angiogenesis which is thought to maintains tissue perfusion. In this study, we elucidate the complex mechanisms involved in water balance in the desert-adapted cactus mouse, P. eremicus. By prioritizing a comprehensive analysis of whole-organism physiology and multi-tissue gene expression in a simulated desert environment, we describe the complex and successful response of regulatory processes.

When the tap runs dry: the physiological effects of acute experimental dehydration in Peromyscus eremicus.

Desert organisms have evolved physiological, biochemical and genomic mechanisms to survive the extreme aridity of desert environments. Studying desert-adapted species provides a unique opportunity to investigate the survival strategies employed by organisms in some of the harshest habitats on Earth. Two of the primary challenges faced in desert environments are maintaining water balance and thermoregulation. We collected data in a simulated desert environment and a captive colony of cactus mice (Peromyscus eremicus) and used lab-based experiments with real time physiological measurements; energy expenditure, water loss rate and respiratory exchange rate, to characterize the response to water deprivation. Mice without access to water had significantly lower energy expenditures and in turn, reduced water loss compared to mice with access to water after the first 24 h of the experiment. Additionally, we observed significant mass loss that is probably due to dehydration-associated anorexia a response to limit fluid loss by reducing waste and the solute load as well as allowing water reabsorption from the kidneys and gastrointestinal tract. Finally, we observed body temperature correlated with sex, with males without access to water maintaining body temperature when compared with hydrated males, whereas body temperature decreased for females without access to water, suggesting daily metabolic depression in females.

White-nose syndrome restructures bat skin microbiomes.

IMPORTANCE: Inherent complexities in the composition of microbiomes can often preclude investigations of microbe-associated diseases. Instead of single organisms being associated with disease, community characteristics may be more relevant. Longitudinal microbiome studies of the same individual bats as pathogens arrive and infect a population are the ideal experiment but remain logistically challenging; therefore, investigations like our approach that are able to correlate invasive pathogens to alterations within a microbiome may be the next best alternative. The results of this study potentially suggest that microbiome-host interactions may determine the likelihood of infection. However, the contrasting relationship between Pd and the bacterial microbiomes of Myotis lucifugus and Perimyotis subflavus indicate that we are just beginning to understand how the bat microbiome interacts with a fungal invader such as Pd.

When the tap runs dry: The physiological effects of acute experimental dehydration in Peromyscus eremicus.

Desert organisms have evolved physiological, biochemical, and genomic mechanisms to survive the extreme aridity of desert environments. Studying desert-adapted species provides a unique opportunity to investigate the survival strategies employed by organisms in some of the harshest habitats on Earth. Two of the primary challenges faced in desert environments are maintaining water balance and thermoregulation. We collected data in a simulated desert environment and a captive colony of cactus mice (Peromyscus eremicus) and used lab-based experiments with real time physiological measurements to characterize the response to water-deprivation. Mice without access to water had significantly lower energy expenditures and in turn, reduced water loss compared to mice with access to water after the first 24 hours of the experiment. Additionally, we observed significant weight loss likely related to dehydration-associated anorexia a response to limit fluid loss by reducing waste and the solute load as well as allowing water reabsorption from the kidneys and gastrointestinal tract. Finally, we observed body temperature correlated with sex, with males without access to water maintaining body temperature when compared to hydrated males while body temperature decreased for females without access to water compared to hydrated, suggesting daily torpor in females.

Chromosome size affects sequence divergence between species through the interplay of recombination and selection.

The structure of the genome shapes the distribution of genetic diversity and sequence divergence. To investigate how the relationship between chromosome size and recombination rate affects sequence divergence between species, we combined empirical analyses and evolutionary simulations. We estimated pairwise sequence divergence among 15 species from three different mammalian clades-Peromyscus rodents, Mus mice, and great apes-from chromosome-level genome assemblies. We found a strong significant negative correlation between chromosome size and sequence divergence in all species comparisons within the Peromyscus and great apes clades but not the Mus clade, suggesting that the dramatic chromosomal rearrangements among Mus species may have masked the ancestral genomic landscape of divergence in many comparisons. Our evolutionary simulations showed that the main factor determining differences in divergence among chromosomes of different sizes is the interplay of recombination rate and selection, with greater variation in larger populations than in smaller ones. In ancestral populations, shorter chromosomes harbor greater nucleotide diversity. As ancestral populations diverge, diversity present at the onset of the split contributes to greater sequence divergence in shorter chromosomes among daughter species. The combination of empirical data and evolutionary simulations revealed that chromosomal rearrangements, demography, and divergence times may also affect the relationship between chromosome size and divergence, thus deepening our understanding of the role of genome structure in the evolution of species divergence.

Prolactin and prolactin receptor expression in the HPG axis and crop during parental care in both sexes of a biparental bird (Columba livia).

During breeding, multiple circulating hormones, including prolactin, facilitate reproductive transitions in species that exhibit parental care. Prolactin underlies parental behaviors and related physiological changes across many vertebrates, including birds and mammals. While circulating prolactin levels often fluctuate across breeding, less is known about how relevant target tissues vary in their prolactin responsiveness via prolactin receptor (PRLR) expression. Recent studies have also investigated prolactin (PRL) gene expression outside of the pituitary (i.e., extra-pituitary PRL), but how PRL gene expression varies during parental care in non-pituitary tissue (e.g., hypothalamus, gonads) remains largely unknown. Further, it is unclear if and how tissue-specific PRL and PRLR vary between the sexes during biparental care. To address this, we measured PRL and PRLR gene expression in tissues relevant to parental care, the endocrine reproductive hypothalamic-pituitary- gonadal (HPG) axis and the crop (a tissue with a similar function as the mammalian mammary gland), across various reproductive stages in both sexes of a biparental bird, the rock dove (Columba livia). We also assessed how these genes responded to changes in offspring presence by adding chicks mid-incubation, simulating an early hatch when prolactin levels were still moderately low. We found that pituitary PRL expression showed similar increases as plasma prolactin levels, and detected extra-pituitary PRL in the hypothalamus, gonads and crop. Hypothalamic and gonadal PRLR expression also changed as birds began incubation. Crop PRLR expression correlated with plasma prolactin, peaking when chicks hatched. In response to replacing eggs with a novel chick mid-incubation, hypothalamic and gonadal PRL and PRLR gene expression differed significantly compared to mid-incubation controls, even when plasma prolactin levels did not differ. We also found sex differences in PRL and PRLR that suggest gene expression may allow males to compensate for lower levels in prolactin by upregulating PRLR in all tissues. Overall, this study advances our understanding of how tissue-specific changes in responsiveness to parental hormones may differ across key reproductive transitions, in response to offspring cues, and between the sexes.

Uncovering the Sex-Specific Endocrine Responses to Reproduction and Parental Care.

Hormones mediate physiological and behavioral changes in adults as they transition into reproduction. In this study, we characterize the circulating levels of five key hormones involved in reproduction in rock doves (Columba livia): corticosterone, progesterone, estradiol, testosterone, and prolactin using univariate and multivariate approaches. We show similar patterns as previous studies in the overall patterns in circulating levels of these hormones, i.e., testosterone (males) and estradiol (females) high during nest-building or egg-laying, prolactin increasing at mid-incubation and peaking at hatching (both sexes), and elevated corticosterone levels in later incubation and early nestling development. In our investigation of hormone co-variation, we find a strong correlation between prolactin and corticosterone across sampling stages and similarities in earlier (early to mid-incubation) compared to later (late incubation to nestling d9) sampling stages in males and females. Finally, we utilized experimental manipulations to simulate nest loss or altered caregiving lengths to test whether external cues, internal timing, or a combination of these factors contributed most to hormone variation. Following nest loss, we found that both males and females responded to the external cue. Males generally responded quickly following nest loss by increasing circulating testosterone, but this response was muted when nest loss occurred early in reproduction. Similar treatment type, e.g., removal of eggs, clustered similarly in hormone space. These results suggest internal drivers limited male response early in reproduction to nest loss. In contrast, circulating levels of these hormones in females either did not change or decreased following nest manipulation suggesting responsiveness to external drivers, but unlike males, this result suggests that reproductive processes were decreasing.

Disentangling environmental drivers of circadian metabolism in desert-adapted mice.

Metabolism is a complex phenotype shaped by natural environmental rhythms, as well as behavioral, morphological and physiological adaptations. Metabolism has been historically studied under constant environmental conditions, but new methods of continuous metabolic phenotyping now offer a window into organismal responses to dynamic environments, and enable identification of abiotic controls and the timing of physiological responses relative to environmental change. We used indirect calorimetry to characterize metabolic phenotypes of the desert-adapted cactus mouse (Peromyscus eremicus) in response to variable environmental conditions that mimic their native environment versus those recorded under constant warm and constant cool conditions, with a constant photoperiod and full access to resources. We found significant sexual dimorphism, with males being more prone to dehydration than females. Under circadian environmental variation, most metabolic shifts occurred prior to physical environmental change and the timing was disrupted under both constant treatments. The ratio of CO2 produced to O2 consumed (the respiratory quotient) reached greater than 1.0 only during the light phase under diurnally variable conditions, a pattern that strongly suggests that lipogenesis contributes to the production of energy and endogenous water. Our results are consistent with historical descriptions of circadian torpor in this species (torpid by day, active by night), but reject the hypothesis that torpor is initiated by food restriction or negative water balance.

Isolating the Role of Corticosterone in the Hypothalamic-Pituitary-Gonadal Transcriptomic Stress Response.

Investigation of the negative impacts of stress on reproduction has largely centered around the effects of the adrenal steroid hormone, corticosterone (CORT), and its influence on a system of tissues vital for reproduction-the hypothalamus of the brain, the pituitary gland, and the gonads (the HPG axis). Research on the action of CORT on the HPG axis has predominated the stress and reproductive biology literature, potentially overshadowing other influential mediators. To gain a more complete understanding of how elevated CORT affects transcriptomic activity of the HPG axis, we experimentally examined its role in male and female rock doves (Columba livia). We exogenously administrated CORT to mimic circulating levels during the stress response, specifically 30 min of restraint stress, an experimental paradigm known to increase circulating CORT in vertebrates. We examined all changes in transcription within each level of the HPG axis as compared to both restraint-stressed birds and vehicle-injected controls. We also investigated the differential transcriptomic response to CORT and restraint-stress in each sex. We report causal and sex-specific effects of CORT on the HPG transcriptomic stress response. Restraint stress caused 1567 genes to uniquely differentially express while elevated circulating CORT was responsible for the differential expression of 304 genes. Only 108 genes in females and 8 in males differentially expressed in subjects that underwent restraint stress and those who were given exogenous CORT. In response to elevated CORT and restraint-stress, both sexes shared the differential expression of 5 genes, KCNJ5, CISH, PTGER3, CEBPD, and ZBTB16, all located in the pituitary. The known functions of these genes suggest potential influence of elevated CORT on immune function and prolactin synthesis. Gene expression unique to each sex indicated that elevated CORT affected more gene transcription in females than males (78 genes versus 3 genes, respectively). To our knowledge, this is the first study to isolate the role of CORT in HPG genomic transcription during a stress response. We present an extensive and openly accessible view of the role corticosterone in the HPG transcriptomic stress response. Because the HPG system is well conserved across vertebrates, these data have the potential to inspire new therapeutic strategies for reproductive dysregulation in multiple vertebrate systems, including our own.

The genomics of mimicry: Gene expression throughout development provides insights into convergent and divergent phenotypes in a Müllerian mimicry system.

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which multiple species share the same conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which vivid colour and pattern are produced in a Müllerian mimicry complex of poison frogs. We did this by first assembling a high-quality de novo genome assembly for the mimic poison frog Ranitomeya imitator. This assembled genome is 6.8 Gbp in size, with a contig N50 of 300 Kbp R. imitator and two colour morphs from both Ranitomeya fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes that are differentially expressed throughout development, many of them related to melanocyte development, melanin synthesis, iridophore development and guanine synthesis. Polytypic differences within species may be the result of differences in expression and/or timing of expression, whereas convergence for colour pattern between species does not appear to be due to the same changes in gene expression. In addition, we identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species. Finally, we hypothesize that genes in the keratin family are important for producing different structural colours within these frogs.

Feather Gene Expression Elucidates the Developmental Basis of Plumage Iridescence in African Starlings.

Iridescence is widespread in the living world, occurring in organisms as diverse as bacteria, plants, and animals. Yet, compared to pigment-based forms of coloration, we know surprisingly little about the developmental and molecular bases of the structural colors that give rise to iridescence. Birds display a rich diversity of iridescent structural colors that are produced in feathers by the arrangement of melanin-containing organelles called melanosomes into nanoscale configurations, but how these often unusually shaped melanosomes form, or how they are arranged into highly organized nanostructures, remains largely unknown. Here, we use functional genomics to explore the developmental basis of iridescent plumage using superb starlings (Lamprotornis superbus), which produce both iridescent blue and non-iridescent red feathers. Through morphological and chemical analyses, we confirm that hollow, flattened melanosomes in iridescent feathers are eumelanin-based, whereas melanosomes in non-iridescent feathers are solid and amorphous, suggesting that high pheomelanin content underlies red coloration. Intriguingly, the nanoscale arrangement of melanosomes within the barbules was surprisingly similar between feather types. After creating a new genome assembly, we use transcriptomics to show that non-iridescent feather development is associated with genes related to pigmentation, metabolism, and mitochondrial function, suggesting non-iridescent feathers are more energetically expensive to produce than iridescent feathers. However, iridescent feather development is associated with genes related to structural and cellular organization, suggesting that, while nanostructures themselves may passively assemble, barbules and melanosomes may require active organization to give them their shape. Together, our analyses suggest that iridescent feathers form through a combination of passive self-assembly and active processes.

Signal, bias, and the role of transcriptome assembly quality in phylogenomic inference.

BACKGROUND: Phylogenomic approaches have great power to reconstruct evolutionary histories, however they rely on multi-step processes in which each stage has the potential to affect the accuracy of the final result. Many studies have empirically tested and established methodology for resolving robust phylogenies, including selecting appropriate evolutionary models, identifying orthologs, or isolating partitions with strong phylogenetic signal. However, few have investigated errors that may be initiated at earlier stages of the analysis. Biases introduced during the generation of the phylogenomic dataset itself could produce downstream effects on analyses of evolutionary history. Transcriptomes are widely used in phylogenomics studies, though there is little understanding of how a poor-quality assembly of these datasets could impact the accuracy of phylogenomic hypotheses. Here we examined how transcriptome assembly quality affects phylogenomic inferences by creating independent datasets from the same input data representing high-quality and low-quality transcriptome assembly outcomes. RESULTS: By studying the performance of phylogenomic datasets derived from alternative high- and low-quality assembly inputs in a controlled experiment, we show that high-quality transcriptomes produce richer phylogenomic datasets with a greater number of unique partitions than low-quality assemblies. High-quality assemblies also give rise to partitions that have lower alignment ambiguity and less compositional bias. In addition, high-quality partitions hold stronger phylogenetic signal than their low-quality transcriptome assembly counterparts in both concatenation- and coalescent-based analyses. CONCLUSIONS: Our findings demonstrate the importance of transcriptome assembly quality in phylogenomic analyses and suggest that a portion of the uncertainty observed in such studies could be alleviated at the assembly stage.

Limited Evidence for Parallel Evolution Among Desert-Adapted Peromyscus Deer Mice.

Warming climate and increasing desertification urge the identification of genes involved in heat and dehydration tolerance to better inform and target biodiversity conservation efforts. Comparisons among extant desert-adapted species can highlight parallel or convergent patterns of genome evolution through the identification of shared signatures of selection. We generate a chromosome-level genome assembly for the canyon mouse (Peromyscus crinitus) and test for a signature of parallel evolution by comparing signatures of selective sweeps across population-level genomic resequencing data from another congeneric desert specialist (Peromyscus eremicus) and a widely distributed habitat generalist (Peromyscus maniculatus), that may be locally adapted to arid conditions. We identify few shared candidate loci involved in desert adaptation and do not find support for a shared pattern of parallel evolution. Instead, we hypothesize divergent molecular mechanisms of desert adaptation among deer mice, potentially tied to species-specific historical demography, which may limit or enhance adaptation. We identify a number of candidate loci experiencing selective sweeps in the P. crinitus genome that are implicated in osmoregulation (Trypsin, Prostasin) and metabolic tuning (Kallikrein, eIF2-alpha kinase GCN2, APPL1/2), which may be important for accommodating hot and dry environmental conditions.

A total crapshoot? Evaluating bioinformatic decisions in animal diet metabarcoding analyses.

Metabarcoding studies provide a powerful approach to estimate the diversity and abundance of organisms in mixed communities in nature. While strategies exist for optimizing sample and sequence library preparation, best practices for bioinformatic processing of amplicon sequence data are lacking in animal diet studies. Here we evaluate how decisions made in core bioinformatic processes, including sequence filtering, database design, and classification, can influence animal metabarcoding results. We show that denoising methods have lower error rates compared to traditional clustering methods, although these differences are largely mitigated by removing low-abundance sequence variants. We also found that available reference datasets from GenBank and BOLD for the animal marker gene cytochrome oxidase I (COI) can be complementary, and we discuss methods to improve existing databases to include versioned releases. Taxonomic classification methods can dramatically affect results. For example, the commonly used Barcode of Life Database (BOLD) Classification API assigned fewer names to samples from order through species levels using both a mock community and bat guano samples compared to all other classifiers (vsearch-SINTAX and q2-feature-classifier's BLAST + LCA, VSEARCH + LCA, and Naive Bayes classifiers). The lack of consensus on bioinformatics best practices limits comparisons among studies and may introduce biases. Our work suggests that biological mock communities offer a useful standard to evaluate the myriad computational decisions impacting animal metabarcoding accuracy. Further, these comparisons highlight the need for continual evaluations as new tools are adopted to ensure that the inferences drawn reflect meaningful biology instead of digital artifacts.

The Genome of the Softshell Clam Mya arenaria and the Evolution of Apoptosis.

Apoptosis is a fundamental feature of multicellular animals and is best understood in mammals, flies, and nematodes, with the invertebrate models being thought to represent a condition of ancestral simplicity. However, the existence of a leukemia-like cancer in the softshell clam Mya arenaria provides an opportunity to re-evaluate the evolution of the genetic machinery of apoptosis. Here, we report the whole-genome sequence for M. arenaria which we leverage with existing data to test evolutionary hypotheses on the origins of apoptosis in animals. We show that the ancestral bilaterian p53 locus, a master regulator of apoptosis, possessed a complex domain structure, in contrast to that of extant ecdysozoan p53s. Further, ecdysozoan taxa, but not chordates or lophotrochozoans like M. arenaria, show a widespread reduction in apoptosis gene copy number. Finally, phylogenetic exploration of apoptosis gene copy number reveals a striking linkage with p53 domain complexity across species. Our results challenge the current understanding of the evolution of apoptosis and highlight the ancestral complexity of the bilaterian apoptotic tool kit and its subsequent dismantlement during the ecdysozoan radiation.

De novo assembly and functional annotation of the heart + hemolymph transcriptome in the Caribbean spiny lobster Panulirus argus.

The spiny lobster, Panulirus argus, is an ecologically relevant species in shallow water coral reefs and a target of the most lucrative fishery in the greater Caribbean region. This study reports, for the first time, the heart + hemolymph transcriptome of the Caribbean spiny lobster Panulirus argus assembled from short Illumina 150bp PE raw reads. A total of 80,152,094 raw reads were assembled using the Oyster River Protocol pipeline. The assembly resulted in a total of 254,773 transcripts. Functional gene annotation was conducted using the software package 'dammit'. Lastly, gene enrichment analyses were conducted using the Gene Ontology (GO) and KEGG pathway (Kaas) databases. This resource will be of utmost importance in future research aiming at exploring the effect of local and regional anthropogenic disturbances, as well as global climate change on the molecular physiology of this overexploited species.

Stress-mediated convergence of splicing landscapes in male and female rock doves.

BACKGROUND: The process of alternative splicing provides a unique mechanism by which eukaryotes are able to produce numerous protein products from the same gene. Heightened variability in the proteome has been thought to potentiate increased behavioral complexity and response flexibility to environmental stimuli, thus contributing to more refined traits on which natural and sexual selection can act. While it has been long known that various forms of environmental stress can negatively affect sexual behavior and reproduction, we know little of how stress can affect the alternative splicing associated with these events, and less still about how splicing may differ between sexes. Using the model of the rock dove (Columba livia), our team previously uncovered sexual dimorphism in the basal and stress-responsive gene transcription of a biological system necessary for facilitating sexual behavior and reproduction, the hypothalamic-pituitary-gonadal (HPG) axis. In this study, we delve further into understanding the mechanistic underpinnings of how changes in the environment can affect reproduction by testing the alternative splicing response of the HPG axis to an external stressor in both sexes. RESULTS: This study reveals dramatic baseline differences in HPG alternative splicing between males and females. However, after subjecting subjects to a restraint stress paradigm, we found a significant reduction in these differences between the sexes. In both stress and control treatments, we identified a higher incidence of splicing activity in the pituitary in both sexes as compared to other tissues. Of these splicing events, the core exon event is the most abundant form of splicing and more frequently occurs in the coding regions of the gene. Overall, we observed less splicing activity in the 3'UTR (untranslated region) end of transcripts than the 5'UTR or coding regions. CONCLUSIONS: Our results provide vital new insight into sex-specific aspects of the stress response on the HPG axis at an unprecedented proximate level. Males and females uniquely respond to stress, yet exhibit splicing patterns suggesting a convergent, optimal splicing landscape for stress response. This information has the potential to inform evolutionary theory as well as the development of highly-specific drug targets for stress-induced reproductive dysfunction.