Tissue-specific expression profiles and positive selection analysis in the tree swallow (Tachycineta bicolor) using a de novo transcriptome assembly

Maternal social environment shapes yolk testosterone allocation and embryonic neural gene expression in tree swallows

Offspring from females breeding in competitive social environments are often exposed to more testosterone (T) during embryonic development, which can affect traits from growth to behavior in potentially adaptive ways. Despite the important role of maternally derived steroids in shaping offspring development, the molecular mechanisms driving these processes are currently unclear. Here, we use tree swallows (Tachycineta bicolor) to explore the effects of the maternal social environment on yolk T concentrations and genome-wide patterns of neural gene expression in embryos. We measured aggressive interactions among females breeding at variable densities and collected their eggs at two timepoints, including the day laid to measure yolk T concentrations and on embryonic day 11 to measure gene expression in whole brain samples. We found that females breeding in high-density sites experienced elevated rates of physical aggression and their eggs had higher yolk T concentrations. A differential gene expression and weighted gene co-expression network analysis indicated that embryos from high-density sites experienced an upregulation of genes involved in hormone, circulatory, and immune processes, and these gene expression patterns were correlated with yolk T levels and aggression. Genes implicated in neural development were additionally downregulated in embryos from high-density sites. These data highlight how early neurogenomic processes may be affected by the maternal social environment, giving rise to phenotypic plasticity in offspring.

Among-population variation in telomere regulatory proteins and their potential role as hidden drivers of intraspecific variation in life history

Biologists aim to explain patterns of growth, reproduction and ageing that characterize life histories, yet we are just beginning to understand the proximate mechanisms that generate this diversity. Existing research in this area has focused on telomeres but has generally overlooked the telomere's most direct mediator, the shelterin protein complex. Shelterin proteins physically interact with the telomere to shape its shortening and repair. They also regulate metabolism and immune function, suggesting a potential role in life history variation in the wild. However, research on shelterin proteins is uncommon outside of biomolecular work. Intraspecific analyses can play an important role in resolving these unknowns because they reveal subtle variation in life history within and among populations. Here, we assessed ecogeographic variation in shelterin protein abundance across eight populations of tree swallow (Tachycineta bicolor) with previously documented variation in environmental and life history traits. Using the blood gene expression of four shelterin proteins in 12-day-old nestlings, we tested the hypothesis that shelterin protein gene expression varies latitudinally and in relation to both telomere length and life history. Shelterin protein gene expression differed among populations and tracked non-linear variation in latitude: nestlings from mid-latitudes expressed nearly double the shelterin mRNA on average than those at more northern and southern sites. However, telomere length was not significantly related to latitude. We next assessed whether telomere length and shelterin protein gene expression correlate with 12-day-old body mass and wing length, two proxies of nestling growth linked to future fecundity and survival. We found that body mass and wing length correlated more strongly (and significantly) with shelterin protein gene expression than with telomere length. These results highlight telomere regulatory shelterin proteins as potential mediators of life history variation among populations. Together with existing research linking shelterin proteins and life history variation within populations, these ecogeographic patterns underscore the need for continued integration of ecology, evolution and telomere biology, which together will advance understanding of the drivers of life history variation in nature.

Genomic Variation, Population History, and Long-Term Genetic Adaptation to High Altitudes in Tibetan Partridge (Perdix hodgsoniae)

Species residing across elevational gradients display adaptations in response to environmental changes such as oxygen availability, ultraviolet radiation, and temperature. Here, we study genomic variation, gene expression, and long-term adaptation in Tibetan Partridge (Perdix hodgsoniae) populations residing across the elevational gradient of the Tibetan Plateau. We generated a high-quality draft genome and used it to carry out downstream population genomic and transcriptomic analysis. The P. hodgsoniae populations residing across various elevations were genetically distinct, and their phylogenetic clustering was consistent with their geographic distribution. We identified possible evidence of gene flow between populations residing in <3,000 and >4,200 m elevation that is consistent with known habitat expansion of high-altitude populations of P. hodgsoniae to a lower elevation. We identified a 60 kb haplotype encompassing the Estrogen Receptor 1 (ESR1) gene, showing strong genetic divergence between populations of P. hodgsoniae. We identified six single nucleotide polymorphisms within the ESR1 gene fixed for derived alleles in high-altitude populations that are strongly conserved across vertebrates. We also compared blood transcriptome profiles and identified differentially expressed genes (such as GAPDH, LDHA, and ALDOC) that correlated with differences in altitude among populations of P. hodgsoniae. These candidate genes from population genomics and transcriptomics analysis were enriched for neutrophil degranulation and glycolysis pathways, which are known to respond to hypoxia and hence may contribute to long-term adaptation to high altitudes in P. hodgsoniae. Our results highlight Tibetan Partridges as a useful model to study molecular mechanisms underlying long-term adaptation to high altitudes.

Genome assembly of the foot-flagging frog, Staurois parvus: a resource for understanding mechanisms of behavior

Elaborate and skilled movements of the body have been selected in a variety of species as courtship and rivalry signals. One roadblock in studying these behaviors has been a lack of resources for understanding how they evolved at the genetic level. The Bornean rock frog (Staurois parvus) is an ideal species in which to address this issue. Males wave their hindlimbs in a "foot-flagging" display when competing for mates. The evolution of foot flagging in S. parvus and other species is accompanied by increases in the expression of the androgen receptor gene within its neuromuscular system, but it remains unclear what genetic or transcriptional changes are associated with this behavioral phenotype. We have now assembled the genome of S. parvus, resulting in 3.98 Gbp of 22,402 contigs with an N50 of 611,229 bp. The genome will be a resource for finding genes related to the physiology underlying foot flagging and to adaptations of the neuromuscular system. As a first application of the genome, we also began work in comparative genomics and differential gene expression analysis. We show that the androgen receptor is diverged from other anuran species, and we identify unique expression patterns of genes in the spinal cord and leg muscle that are important for axial patterning, cell specification and morphology, or muscle contraction. This genome will continue to be an important tool for future -omics studies to understand the evolution of elaborate signaling behaviors in this and potentially related species.

Rapid adaptation to a novel pathogen through disease tolerance in a wild songbird

Animal hosts can adapt to emerging infectious disease through both disease resistance, which decreases pathogen numbers, and disease tolerance, which limits damage during infection without limiting pathogen replication. Both resistance and tolerance mechanisms can drive pathogen transmission dynamics. However, it is not well understood how quickly host tolerance evolves in response to novel pathogens or what physiological mechanisms underlie this defense. Using natural populations of house finches (Haemorhous mexicanus) across the temporal invasion gradient of a recently emerged bacterial pathogen (Mycoplasma gallisepticum), we find rapid evolution of tolerance (<25 years). In particular, populations with a longer history of MG endemism have less pathology but similar pathogen loads compared with populations with a shorter history of MG endemism. Further, gene expression data reveal that more-targeted immune responses early in infection are associated with tolerance. These results suggest an important role for tolerance in host adaptation to emerging infectious diseases, a phenomenon with broad implications for pathogen spread and evolution.

Gene expression in the female tree swallow brain is associated with inter- and intra-population variation in glucocorticoid levels

Studies of the evolutionary causes and consequences of variation in circulating glucocorticoids (GCs) have begun to reveal how they are shaped by selection. Yet the extent to which variation in circulating hormones reflects variation in other important regulators of the hypothalamic-pituitary-adrenal (HPA) axis, and whether these relationships vary among populations inhabiting different environments, remain poorly studied. Here, we compare gene expression in the brain of female tree swallows (Tachycineta bicolor) from populations that breed in environments that differ in their unpredictability. We find evidence of inter-population variation in the expression of glucocorticoid and mineralocorticoid receptors in the hypothalamus, with the highest gene expression in a population from an extreme environment, and lower expression in a population from a more consistent environment as well as in birds breeding at an environmentally variable high-altitude site that are part of a population that inhabits a mixture of high and low altitude habitats. Within some populations, variation in circulating GCs predicted differences in gene expression, particularly in the hypothalamus. However, some patterns were present in all populations, whereas others were not. These results are consistent with the idea that some combination of local adaptation and phenotypic plasticity may modify components of the HPA axis affecting stress resilience. Our results also underscore that a comprehensive understanding of the function and evolution of the stress response cannot be gained from measuring circulating hormones alone, and that future studies that apply a more explicitly evolutionary approach to important regulatory traits are likely to provide significant insights.

The telomere regulatory gene POT1 responds to stress and predicts performance in nature: Implications for telomeres and life history evolution

Telomeres are emerging as correlates of fitness-related traits and may be important mediators of ecologically relevant variation in life history strategies. Growing evidence suggests that telomere dynamics can be more predictive of performance than length itself, but very little work considers how telomere regulatory mechanisms respond to environmental challenges or influence performance in nature. Here, we combine observational and experimental data sets from free-living tree swallows (Tachycineta bicolor) to assess how performance is predicted by the telomere regulatory gene POT1, which encodes a shelterin protein that sterically blocks telomerase from repairing the telomere. First, we show that lower POT1 gene expression in the blood was associated with higher female quality, that is, earlier breeding and heavier body mass. We next challenged mothers with an immune stressor (lipopolysaccharide injection) that led to "sickness" in mothers and 24 h of food restriction in their offspring. While POT1 did not respond to maternal injection, females with lower constitutive POT1 gene expression were better able to maintain feeding rates following treatment. Maternal injection also generated a 1-day stressor for chicks, which responded with lower POT1 gene expression and elongated telomeres. Other putatively stress-responsive mechanisms (i.e., glucocorticoids, antioxidants) showed marginal responses in stress-exposed chicks. Model comparisons indicated that POT1 mRNA abundance was a largely better predictor of performance than telomere dynamics, indicating that telomere regulators may be powerful modulators of variation in life history strategies.

Mathematical model for the relationship between single-cell and bulk gene expression to clarify the interpretation of bulk gene expression data

BACKGROUND

Differential expression analysis is a standard approach in molecular biology. For example, genes whose expression levels differ between diseased and non-diseased samples are considered to be associated with that disease. On the other hand, differential variability analysis focuses on the differences of the variances of gene expression between sample groups. Although differential variability is also known to capture biological information, its interpretation remains unclear and controversial. Recent single-cell analyses have revealed that differences between sample groups can affect gene expression in a cellular subset-specific manner or by altering the proportion of a particular cellular subset. The aim of this study is to clarify the interpretation of mean and variance of bulk gene expression data.

METHOD

We developed a mathematical model in which the bulk gene expression value is proportional to the mean value of the single-cell gene expression profile. Based on this model, we performed theoretical, simulated and real single-cell RNA-seq data analyses.

RESULT AND CONCLUSION

We identified how differences in single-cell gene expression profiles affect the differences in the mean and the variance of bulk gene expression. It is shown that differential expression analysis of bulk expression data can overlook significant changes in gene expression at the single-cell level. Further, differential variability analysis capture the complex feature affected by different gene expression shifts for each subset, changes in the proportions of cellular subsets, and variation in single-cell distribution parameters among samples.

Evolution of hormone-phenotype couplings and hormone-genome interactions

When selection favors a new relationship between a cue and a hormonally mediated response, adaptation can proceed by altering the hormonal signal that is produced or by altering the phenotypic response to the hormonal signal. The field of evolutionary endocrinology has made considerable progress toward understanding the evolution of hormonal signals, but we know much less about the evolution of hormone-phenotype couplings, particularly at the hormone-genome interface. We briefly review and classify the mechanisms through which these hormone-phenotype couplings likely evolve, using androgens and their receptors and genomic response elements to illustrate our view. We then present two empirical studies of hormone-phenotype couplings, one rooted in evolutionary quantitative genetics and another in comparative transcriptomics, each focused on the regulation of sexually dimorphic phenotypes by testosterone (T) in the brown anole lizard (Anolis sagrei). First, we illustrate the potential for hormone-phenotype couplings to evolve by showing that coloration of the dewlap (an ornament used in behavioral displays) exhibits significant heritability in its responsiveness to T, implying that anoles harbor genetic variance in the architecture of hormonal pleiotropy. Second, we combine T manipulations with analyses of the liver transcriptome to ask whether and how statistical methods for characterizing modules of co-expressed genes and in silico techniques for identifying androgen response elements (AREs) can improve our understanding of hormone-genome interactions. We conclude by emphasizing important avenues for future work at the hormone-genome interface, particularly those conducted in a comparative evolutionary framework.

How experimental competition changes ovarian gene activity in free-living birds: Implications for steroidogenesis, maternal effects, and beyond

The ovary plays an important role in mediating both a female's response to her social environment and communicating it to her developing offspring via maternal effects. Past work has focused on how ovarian hormones respond to competition, but we know little about how the broader ovarian transcriptomic landscape changes, either during or after competition, giving us a narrow perspective on how socially induced phenotypes arise. Here, we experimentally generated social competition among wild, cavity-nesting female birds (tree swallows, Tachycineta bicolor), a species in which females lack a socially induced rise in circulating testosterone but they nevertheless increase allocation to eggs. After territory settlement, we reduced availability of nesting cavities, generating heightened competition; within 24 h we reversed the manipulation, causing aggressive interactions to subside. We measured ovarian transcriptomic responses at the peak of competition and 48 h later, along with date-matched controls. Network analyses indicated that competing females experienced an immediate and temporary decrease in the expression of genes involved in the early stages of steroidogenesis, and this was moderately correlated with plasma testosterone; however, two days after competition had ended, there was a marked increase in the expression of genes involved in the final stages of steroidogenesis, including HSD17B1. Gene networks related to the cell cycle, muscle performance, and extracellular matrix organization also displayed altered activity. Although the functional consequences of these findings are unclear, they shed light on socially responsive ovarian genomic mechanisms that could potentially exert lasting effects on behavior, reproduction, and maternal effects.

A chromosome-scale genome assembly of the false clownfish, Amphiprion ocellaris

The false clownfish Amphiprion ocellaris is a popular fish species and an emerging model organism for studying the ecology, evolution, adaptation, and developmental biology of reef fishes. Despite this, high-quality genomic resources for this species are scarce, hindering advanced genomic analyses. Leveraging the power of PacBio long-read sequencing and Hi-C chromosome conformation capture techniques, we constructed a high-quality chromosome-scale genome assembly for the clownfish A. ocellaris. The initial genome assembly comprised of 1,551 contigs of 861.42 Mb, with an N50 of 863.85 kb. Hi-C scaffolding of the genome resulted in 24 chromosomes containing 856.61 Mb. The genome was annotated with 26,797 protein-coding genes and had 96.62% completeness of conserved actinopterygian genes, making this genome the most complete and high quality among published anemonefish genomes. Transcriptomic analysis identified tissue-specific gene expression patterns, with the brain and optic lobe having the largest number of expressed genes. Further, comparative genomic analysis revealed 91 genome elements conserved only in A. ocellaris and its sister species Amphiprion percula, and not in other anemonefish species. These elements are close to genes that are involved in various nervous system functions and exhibited distinct expression patterns in brain tissue, potentially highlighting the genetic toolkits involved in lineage-specific divergence and behaviors of the clownfish branch. Overall, our study provides the highest quality A. ocellaris genome assembly and annotation to date, whilst also providing a valuable resource for understanding the ecology and evolution of reef fishes.

Computational Methods for Prediction of Human Protein-Phenotype Associations: A Review

Deciphering the relationship between human proteins (genes) and phenotypes is one of the fundamental tasks in phenomics research. The Human Phenotype Ontology (HPO) builds upon a standardized logical vocabulary to describe the abnormal phenotypes encountered in human diseases and paves the way towards the computational analysis of their genetic causes. To date, many computational methods have been proposed to predict the HPO annotations of proteins. In this paper, we conduct a comprehensive review of the existing approaches to predicting HPO annotations of novel proteins, identifying missing HPO annotations, and prioritizing candidate proteins with respect to a certain HPO term. For each topic, we first give the formalized description of the problem, and then systematically revisit the published literatures highlighting their advantages and disadvantages, followed by the discussion on the challenges and promising future directions. In addition, we point out several potential topics to be worthy of exploration including the selection of negative HPO annotations and detecting HPO misannotations. We believe that this review will provide insight to the researchers in the field of computational phenotype analyses in terms of comprehending and developing novel prediction algorithms.

Sex-Biased Gene Expression and Evolution in the Cerebrum and Syrinx of Chinese Hwamei (Garrulax canorus)

It is common that males and females display sexual dimorphisms, which usually result from sex-biased gene expression. Chinese hwamei (Garrulax canorus) is a good model for studying sex-biased gene expression because the song between the sexes is quite different. In this study, we analyze cerebrum and syrinx sex-biased gene expression and evolution using the de novo assembled Chinese hwamei transcriptome. In both the cerebrum and syrinx, our study revealed that most female-biased genes were actively expressed in females only, while most male-biased genes were actively expressed in both sexes. In addition, both male- and female-biased genes were enriched on the putative Z chromosome, suggesting the existence of sexually antagonistic genes and the insufficient dosage compensation of the Z-linked genes. We also identified a 9 Mb sex linkage region on the putative 4A chromosome which enriched more than 20% of female-biased genes. Resultantly, male-biased genes in both tissues had significantly higher Ka/Ks and effective number of codons (ENCs) than unbiased genes, and this suggested that male-biased genes which exhibit accelerated divergence may have resulted from positive selection. Taken together, our results initially revealed the reasons for the differences in singing behavior between males and females of Chinese hwamei.

Experimental competition induces immediate and lasting effects on the neurogenome in free-living female birds

Periods of social instability can elicit adaptive phenotypic plasticity to promote success in future competition. However, the underlying molecular mechanisms have primarily been studied in captive and laboratory-reared animals, leaving uncertainty as to how natural competition among free-living animals affects gene activity. Here, we experimentally generated social competition among wild, cavity-nesting female birds (tree swallows, Tachycineta bicolor). After territorial settlement, we reduced the availability of key breeding resources (i.e., nest boxes), generating heightened competition; within 24 h we reversed the manipulation, causing aggressive interactions to subside. We sampled females during the peak of competition and 48 h after it ended, along with date-matched controls. We measured transcriptomic and epigenomic responses to competition in two socially relevant brain regions (hypothalamus and ventromedial telencephalon). Gene network analyses suggest that processes related to energy mobilization and aggression (e.g., dopamine synthesis) were up-regulated during competition, the latter of which persisted 2 d after competition had ended. Cellular maintenance processes were also down-regulated after competition. Competition additionally altered methylation patterns, particularly in pathways related to hormonal signaling, suggesting those genes were transcriptionally poised to respond to future competition. Thus, experimental competition among free-living animals shifts gene expression in ways that may facilitate the demands of competition at the expense of self-maintenance. Further, some of these effects persisted after competition ended, demonstrating the potential for epigenetic biological embedding of the social environment in ways that may prime individuals for success in future social instability.

Transcriptome of the Southern Muriqui Brachyteles arachnoides (Primates:Platyrrhini), a Critically Endangered New World Monkey: Evidence of Adaptive Evolution

The southern muriqui (Brachyteles arachnoides) is the largest neotropical primate. This species is endemic to Brazil and is currently critically endangered due to its habitat destruction. The genetic basis underlying adaptive traits of New World monkeys has been a subject of interest to several investigators, with significant concern about genes related to the immune system. In the absence of a reference genome, RNA-seq and de novo transcriptome assembly have proved to be valuable genetic procedures for accessing gene sequences and testing evolutionary hypotheses. We present here a first report on the sequencing, assembly, annotation and adaptive selection analysis for thousands of transcripts of B. arachnoides from two different samples, corresponding to 13 different blood cells and fibroblasts. We assembled 284,283 transcripts with N50 of 2,940 bp, with a high rate of complete transcripts, with a median high scoring pair coverage of 88.2%, including low expressed transcripts, accounting for 72.3% of complete BUSCOs. We could predict and extract 81,400 coding sequences with 79.8% of significant BLAST hit against the Euarchontoglires SwissProt dataset. Of these 64,929 sequences, 34,084 were considered homologous to Supraprimate proteins, and of the remaining sequences (30,845), 94% were associated with a protein domain or a KEGG Orthology group, indicating potentially novel or specific protein-coding genes of B. arachnoides. We use the predicted protein sequences to perform a comparative analysis with 10 other primates. This analysis revealed, for the first time in an Atelid species, an expansion of APOBEC3G, extending this knowledge to all NWM families. Using a branch-site model, we searched for evidence of positive selection in 4,533 orthologous sets. This evolutionary analysis revealed 132 amino acid sites in 30 genes potentially evolving under positive selection, shedding light on primate genome evolution. These genes belonged to a wide variety of categories, including those encoding the innate immune system proteins (APOBEC3G, OAS2, and CEACAM1) among others related to the immune response. This work generated a set of thousands of complete sequences that can be used in other studies on molecular evolution and may help to unveil the evolution of primate genes. Still, further functional studies are required to provide an understanding of the underlying evolutionary forces modeling the primate genome.

Design, challenges, and the potential of transcriptomics to understand social behavior

Rapid advances in Ribonucleic Acid sequencing (or RNA-seq) technology for analyzing entire transcriptomes of desired tissue samples, or even of single cells at scale, have revolutionized biology in the past decade. Increasing accessibility and falling costs are making it possible to address many problems in biology that were once considered intractable, including the study of various social behaviors. RNA-seq is opening new avenues to understand long-standing questions on the molecular basis of behavioral plasticity and individual variation in the expression of a behavior. As whole transcriptomes are examined, it has become possible to make unbiased discoveries of underlying mechanisms with little or no necessity to predict genes involved in advance. However, researchers need to be aware of technical limitations and have to make specific decisions when applying RNA-seq to study social behavior. Here, we provide a perspective on the applications of RNA-seq and experimental design considerations for behavioral scientists who are unfamiliar with the technology but are considering using it in their research.

Influence of heat stress on reference genes stability in heart and liver of two chickens genotypes

Introduction

Real-time polymerase chain reaction (RT-qPCR) is an important tool for analyzing gene expression. However, before analyzing the expression of target genes, it is crucial to normalize the reference genes, in order to find the most stable gene to be used as an endogenous control. A gene that remains stable in all samples under different treatments is considered a suitable normalizer. In this sense, we aimed to identify stable reference genes for normalization of target genes in the heart and liver tissues from two genetically divergent groups of chickens (Cobb 500^® commercial line and Peloco backyard chickens) under comfort and acute heat stress environmental conditions. Eight reference genes (ACTB, HPRT1, RPL5, EEF1, MRPS27, MRPS30, TFRC and LDHA) were analyzed for expression stability. The samples were obtained from 24 chickens, 12 from the backyard Peloco and 12 from the Cobb 500® line, exposed to two environmental conditions (comfort and heat stress). Comfort temperature was 23°C and heat stress temperature was 39.5°C for one hour. Subsequently, the animals were euthanized, and heart and liver tissue fragments were collected for RNA extraction and amplification. To determine the stability rate of gene expression, three different statistical algorithms were applied: BestKeeper, geNorm and NormFinder, and to obtain an aggregated stability list, the RankAgregg package of R software was used.

Results

The most stable genes using BestKeeper tool, including the two factors (genetic group and environmental condition), were LDHA, RPL5 and MRPS27 for heart tissue, and TFRC, RPL5 and EEF1 for liver tissue. Applying geNorm algorithm, the best reference genes were RPL5, EEF1 and MRPS30 for heart tissue and LDHA, EEF1 and RPL5 for liver. Using the NormFinder algorithm, the best normalizer genes were EEF1, RPL5 and LDHA in heart, and EEF1, RPL5 and ACTB in liver tissue. In the overall ranking obtained by RankAggreg package, considering the three algorithms, the RPL5, EEF1 and LDHA genes were the most stable for heart tissue, whereas RPL5, EEF1 and ACTB were the most stable for liver tissue.

Conclusion

According to the RankAggreg tool classification based on the three different algorithms (BestKeeper, geNorm and NormFinder), the most stable genes were RPL5, EEF1 and LDHA for heart tissue and RPL5, EEF1 and ACTB for liver tissue of chickens subjected to comfort and acute heat stress environmental conditions. However, the best reference genes may vary depending on the experimental conditions of each study, such as different breeds, environmental stressors, and tissues analyzed. Therefore, the need to perform priori studies to assay the best reference genes at the outset of each study is emphasized.

The neurogenomic transition from territory establishment to parenting in a territorial female songbird

Background

The brain plays a critical role in upstream regulation of processes central to mating effort, parental effort, and self-maintenance. For seasonally breeding animals, the brain is likely mediating trade-offs among these processes within a short breeding season, yet research thus far has only explored neurogenomic changes from non-breeding to breeding states or select pathways (e.g., steroids) in male and/or lab-reared animals. Here, we use RNA-seq to explore neural plasticity in three behaviorally relevant neural tissues (ventromedial telencephalon [VmT], hypothalamus [HYPO], and hindbrain [HB]), comparing free-living female tree swallows (Tachycineta bicolor) as they shift from territory establishment to incubation. We additionally highlight changes in aggression-related genes to explore the potential for a neurogenomic shift in the mechanisms regulating aggression, a critical behavior both in establishing and maintaining a territory and in defense of offspring.

Results

HB had few differentially expressed genes, but VmT and HYPO had hundreds. In particular, VmT had higher expression of genes related to neuroplasticity and processes beneficial for competition during territory establishment, but down-regulated immune processes. HYPO showed signs of high neuroplasticity during incubation, and a decreased potential for glucocorticoid signaling. Expression of aggression-related genes also shifted from steroidal to non-steroidal pathways across the breeding season.

Conclusions

These patterns suggest trade-offs between enhanced activity and immunity in the VmT and between stress responsiveness and parental care in the HYPO, along with a potential shift in the mechanisms regulating aggression. Collectively, these data highlight important gene regulatory pathways that may underlie behavioral plasticity in females.