The genome of the white-throated sparrow (Zonotrichia albicollis) contains an inversion polymorphism on chromosome 2 that is linked to predictable variation in a suite of phenotypic traits including plumage color, aggression and parental behavior. Differences in gene expression between the two color morphs, which represent the two common inversion genotypes (ZAL2/ZAL2 and ZAL2/ZAL2m), may therefore advance our understanding of the molecular underpinnings of these phenotypes. To identify genes that are differentially expressed between the two morphs and correlated with behavior, we quantified gene expression and terrirorial aggression, including song, in a population of free-living white-throated sparrows. We analyzed gene expression in two brain regions, the medial amygdala (MeA) and hypothalamus. Both regions are part of a 'social behavior network', which is rich in steroid hormone receptors and previously linked with territorial behavior. Using weighted gene co-expression network analyses, we identified modules of genes that were correlated with both morph and singing behavior. The majority of these genes were located within the inversion, showing the profound effect of the inversion on the expression of genes captured by the rearrangement. These modules were enriched with genes related to retinoic acid signaling and basic cellular functioning. In the MeA, the most prominent pathways were those related to steroid hormone receptor activity. Within these pathways, the only gene encoding such a receptor was ESR1 (estrogen receptor 1), a gene previously shown to predict song rate in this species. The set of candidate genes we identified may mediate the effects of a chromosomal inversion on territorial behavior.

Quantitative real-time PCR (qPCR) is becoming a popular tool for the quantification of gene expression in the brain and endocrine tissues of songbirds. Accurate analysis of qPCR data relies on the selection of appropriate reference genes for normalization, yet few papers on songbirds contain evidence of reference gene validation. Here, we evaluated the expression of ten potential reference genes (18S, ACTB, GAPDH, HMBS, HPRT, PPIA, RPL4, RPL32, TFRC, and UBC) in brain, pituitary, ovary, and testis in two species of songbirds: zebra finch and white-throated sparrow. We used two algorithms, geNorm and NormFinder, to assess the stability of these reference genes in our samples. We found that the suitability of some of the most popular reference genes for target gene normalization in mammals, such as 18S, depended highly on tissue type. Thus, they are not the best choices for brain and gonad in these songbirds. In contrast, we identified alternative genes, such as HPRT, RPL4 and PPIA, that were highly stable in brain, pituitary, and gonad in these species. Our results suggest that the validation of reference genes in mammals does not necessarily extrapolate to other taxonomic groups. For researchers wishing to identify and evaluate suitable reference genes for qPCR in songbirds, our results should serve as a starting point and should help increase the power and utility of songbird models in behavioral neuroendocrinology.

Naïve CD4+ T cells experience weak T cell receptor (TCR) signals induced by self-peptides presented by MHC II. To investigate how these “basal” TCR signals influence responses to agonist TCR ligand stimulation, we analyzed naïve CD4+ cells expressing varying amounts of CD5, Ly6C, and Nur77-GFP, markers that reflect the strength of basal TCR signaling. Phenotypic analyses indicate that the broadest range of basal TCR signal strength can be visualized by a combination of Nur77-GFP and Ly6C. A range of basal TCR signaling is detectable even in populations that express identical TCRs. Whereas moderate basal TCR signal strength correlates with higher IL-2 secretion at early time points following TCR stimulation, weak basal TCR signaling correlated with higher IL-2 secretion at later time points. We identify a population of Nur77-GFPHI Ly6C− cells that could not be reliably marked by either of CD5, Ly6C, or Nur77-GFP alone. These cells experience the strongest basal TCR signaling, consistently produce less IL-2, and express PD-1 and markers associated with anergy, such as Grail and Cbl-b. We propose that adaptation to the strength of basal TCR signaling drives the phenotypic and functional heterogeneity of naïve CD4+ cells.

Major histocompatibility class II (MHC-II) expression is critical for immune responses and is controlled by the MHC-II transactivator CIITA. CIITA is primarily regulated at the transcriptional level and is expressed from three main promoters with myeloid, lymphoid and interferon (IFN)-γ-treated non-hematopoietic cells using promoters pI, pIII and pIV, respectively. Recent studies in non-hematopoietic cells suggest that a series of distal regulatory elements may be involved in regulating CIITA transcription. To identify distal elements in B cells, a DNase I hypersensitivity screen was performed, revealing a series of potential novel regulatory elements. These elements were analyzed computationally and biochemically. Several regions displayed active histone modifications and/or enhanced expression of a reporter gene. Four of the elements interacted with pIII in B cells. These same four regions were also found to interact with pI in splenic dendritic cells (spDC). Intriguingly, examination of the above interactions in pI-knockout-derived spDC showed a switch to the next available promoter, pIII. Extensive DNA methylation was found at the pI region in B cells, suggesting that this promoter is not accessible in B cells. Thus, CIITA expression is likely mediated in hematopoietic cells by common elements with promoter accessibility having a part in promoter choice.

The cumulative effects of T cell receptor (TCR) signal transduction over extended periods of time influences T cell biology, such as the positive selection of immature thymocytes or the proliferative responses of naive T cells. Naive T cells experience recurrent TCR signaling in response to self-antigens in the steady state. However, how these signals influence the responsiveness of naive CD8 + T cells to subsequent agonist TCR stimulation remains incompletely understood. We investigated how naive CD8 + T cells that experienced relatively low or high levels of TCR signaling in response to self-antigens respond to stimulation with foreign antigens. A transcriptional reporter of Nr4a1 (Nur77-GFP) revealed substantial heterogeneity of the amount of TCR signaling naive CD8 + T cells accumulate in the steady state. Nur77-GFP HI cells exhibited diminished T cell activation and secretion of IFNγ and IL-2 relative to Nur77-GFP LO cells in response to agonist TCR stimulation. Differential gene expression analyses revealed upregulation of genes associated with acutely stimulated T cells in Nur77-GFP HI cells but also increased expression of negative regulators such as the phosphatase Sts1. Responsiveness of Nur77-GFP HI cells to TCR stimulation was partially restored at the level of IFNγ secretion by deficiency of Sts1 or the ubiquitin ligase Cbl-b. Our data suggest that extensive accumulation of TCR signaling during steady state conditions induces a recalibration of the responsiveness of naive CD8 + T cells through gene expression changes and negative regulation, at least in part, dependent on Sts1 and Cbl-b. This cell-intrinsic negative feedback loop may allow the immune system to limit the autoreactive potential of highly self-reactive naive CD8 + T cells.

Three distinct promoters control the master regulator of MHC class II expression, CIITA, in a cell type specific manner. Promoter I (pI) CIITA, expressed primarily by dendritic cells and macrophages, expresses a unique isoform that contains a caspase recruitment domain. The activity and function of this isoform is not understood but has been thought to enhance the function of CIITA in antigen presenting cells. To determine if isoform I of CIITA has specific functions, CIITA mutant mice were created in which isoform I was replaced with isoform III sequences. Mice in which pI and the CARD encoding exon were deleted were also created. No defect in the formation of CD4 T cells, the ability to respond to a model antigen, or bacterial or viral challenge was observed in mice lacking CIITA isoform I. Although CIITA and MHC-II expression was decreased in splenic DC, the pI knockout animals expressed CIITA from downstream promoters, suggesting that control of pI activity is mediated by unknown s II distal elements that could act at the pIII, the B cell promoter. Thus, no critical function is linked to the CARD domain of CIITA isoform I with respect to basic immune system development, function and challenge.

In the white-throated sparrow (Zonotrichia albicollis), the second chromosome bears a striking resemblance to sex chromosomes. First, within each breeding pair of birds, one bird is homozygous for the standard arrangement of the chromosome (ZAL2/ZAL2) and its mate is heterozygous for a different version (ZAL2/ZAL2 m ). Second, recombination is profoundly suppressed between the two versions, leading to genetic differentiation between them. Third, the ZAL2 m version is linked with phenotypic traits, such as bright plumage, high aggression, and low parental behavior, which are usually associated with males. These similarities to sex chromosomes suggest that the evolutionary mechanisms that shape sex chromosomes, in particular genetic degeneration of the heterogametic version due to the suppression of recombination, are likely important in this system as well. Here, we investigated patterns of protein sequence evolution and gene expression evolution between the ZAL2 and ZAL2 m chromosomes by whole-genome sequencing and transcriptome analyses. Patterns of protein evolution exhibited only weak signals of genetic degeneration, and few genes harbored signatures of positive selection. We found substantial evidence of transcriptome evolution, such as significant expression divergence between ZAL2 and ZAL2 m alleles and signatures of dosage compensation for highly expressed genes. These results suggest that, early in the evolution of heteromorphic chromosomes, gene expression divergence and dosage compensation can prevail before large-scale genetic degeneration. Our results show further that suppression of recombination between heteromorphic chromosomes can lead to the evolution of alternative (sex-like) behavioral phenotypes before substantial genetic degeneration.

Behavioral evolution relies on genetic changes, yet few behaviors can be traced to specific genetic sequences in vertebrates. Here we provide experimental evidence showing that differentiation of a single gene has contributed to the evolution of divergent behavioral phenotypes in the white-throated sparrow, a common backyard songbird. In this species, a series of chromosomal inversions has formed a supergene that segregates with an aggressive phenotype. The supergene has captured ESR1, the gene that encodes estrogen receptor α (ERα); as a result, this gene is accumulating changes that now distinguish the supergene allele from the standard allele. Our results show that in birds of the more aggressive phenotype, ERα knockdown caused a phenotypic change to that of the less aggressive phenotype. We next showed that in a freeliving population, aggression is predicted by allelic imbalance favoring the supergene allele. Finally, we identified cis-regulatory features, both genetic and epigenetic, that explain the allelic imbalance. This work provides a rare illustration of how genotypic divergence has led to behavioral phenotypic divergence in a vertebrate.