Males and females can be highly dimorphic in metabolism and physiology despite sharing nearly identical genomes, and both sexes respond phenotypically to elevated testosterone, a steroid hormone that alters gene expression. Only recently has it become possible to learn how a hormone such as testosterone affects global gene expression in non-model systems, and whether it affects the same genes in males and females. To investigate the transcriptional mechanisms by which testosterone exerts its metabolic and physiological effects on the periphery, we compared gene expression by sex and in response to experimentally elevated testosterone in a well-studied bird species, the dark-eyed junco (Junco hyemalis). We identified 291 genes in the liver and 658 in the pectoralis muscle that were differentially expressed between males and females. In addition, we identified 1727 genes that were differentially expressed between testosterone-treated and control individuals in at least one tissue and sex. Testosterone treatment altered the expression of only 128 genes in both males and females in the same tissue, and 847 genes were affected significantly differently by testosterone treatment in the two sexes. These substantial differences in transcriptional response to testosterone suggest that males and females may employ different pathways when responding to elevated testosterone, despite the fact that many phenotypic effects of experimentally elevated testosterone are similar in both sexes. In contrast, of the 121 genes that were affected by testosterone treatment in both sexes, 78% were regulated in the same direction (e.g. either higher or lower in testosterone-treated than control individuals) in both males and females. Thus, it appears that testosterone acts through both unique and shared transcriptional pathways in males and females, suggesting multiple mechanisms by which sexual conflict can be mediated.

Males and females often face divergent selective pressures because of inherent differences in reproductive strategy, and these differences can be reflected in life history traits, including reproductive effort, longevity, growth and metabolism (Cox and Calsbeek, 2009). In several species, some life history traits (e.g. longevity and basal metabolic rate) appear to be at sub-optimum levels for each sex, and selection related to these phenotypes acts in opposite directions on males and females (Berg and Maklakov, 2012; Boratynski et al., 2010). This suggests that sexually antagonistic selection has led to a phenotypic compromise (Bonduriansky and Chenoweth, 2009; Chapman et al., 2003). Sexually dimorphic gene expression is thought to provide a solution to sexual conflict (van Doorn, 2009), given that males and females share nearly identical genomes (reviewed in Ellegren and Parsch, 2007). Sexually dimorphic patterns of gene expression are thought to account for many of the physiological differences between the sexes (Xu et al., 2012). For example, sex differences in liver gene expression are substantial in rodents (Corton et al., 2012), and explain several known sex differences in liver metabolism (Gatti et al., 2010).

In many vertebrate species, androgens, such as testosterone, are one of the key regulators of sex differences in many aspects of adult phenotype, including growth and metabolism (Cox et al., 2009; Woodward, 1993; Arnold et al., 1997; Wikelski et al., 1999). Testosterone plays a major role in directing the balance of energy expenditure (Marler and Moore, 1988), generally shifting energy away from metabolic processes of self-maintenance, such as immune function (Folstad and Karter, 1992), toward short-term reproductive efforts, such as courtship (Arnold, 1975; Wiley and Goldizen, 2003) and territory defense (Marler et al., 1995). These phenotypic effects often occur in both males and females, but are likely to affect reproductive success differently in each sex (Ketterson et al., 2005). Thus, there is likely to be conflict between the sexes over the optimal level of circulating testosterone (Boratyński et al., 2010; Mokkonen et al., 2012).

Comparative studies have shown that endogenous levels of testosterone are highly correlated between males and females among species, including in birds (Ketterson et al., 2005; Møller et al., 2005) and fish (Mank, 2007), raising the possibility that selection on circulating testosterone levels in one sex may lead to a similar change in circulating testosterone in the opposite sex. If, however, the sexes differ in their phenotypic and transcriptional response to circulating testosterone, they may be able to reduce this conflict, and each sex may be better able to reach its own optimum phenotypic value (Rice, 1984). Behavior and physiology are known to be sensitive to experimentally elevated testosterone, sometimes in both sexes, sometimes in only one. For example, immune function is sensitive to experimental elevation of testosterone in males of some species (Roberts et al., 2004), but is responsive to testosterone in females of only a subset of these species (Ketterson et al., 2005). The fact that phenotypic sensitivity to testosterone varies between species and sexes strongly suggests evolutionary lability in the genes and phenotypes that respond to testosterone.

To address the role of sexual dimorphism and testosterone in mediating phenotypes via gene expression in a natural system, we studied gene expression in the liver and pectoralis of a wild songbird, the dark-eyed junco [Junco hyemalis (Linnaeus 1758)]. The dark-eyed junco is a mildly dimorphic North American sparrow (Nolan et al., 2002) that has been the focus of ecological research for nearly a century (Rowan, 1925; Miller, 1941; Ketterson et al., 2009), and recent genomic tools have expanded these studies (Peterson et al., 2012). Sex differences and the phenotypic effects of experimentally elevated testosterone have been studied extensively (Ketterson et al., 1991; Ketterson et al., 2009), providing a solid ecological foundation on which to interpret findings from genomic tools (Peterson et al., 2012).

In particular, past research on free-living male and female juncos has detailed many phenotypic consequences of experimental testosterone treatments that maintain levels of testosterone near the early breeding season peak for each sex (Ketterson et al., 1992; Ketterson et al., 1996; Ketterson et al., 2005). Both male and female juncos respond phenotypically to experimentally elevated testosterone by decreasing immune function (Casto et al., 2001; Zysling et al., 2006) and body mass (Ketterson et al., 1991; Clotfelter et al., 2004), along with a number of behavioral responses (reviewed in Ketterson et al., 2005; Ketterson et al., 2009). However, only males increase their activity and home-range size in response to experimental testosterone (Chandler et al., 1994; Lynn et al., 2000; Reichard and Ketterson, 2012). The net result of these and other phenotypic effects of testosterone treatment is an increase in reproductive fitness for males (Reed et al., 2006) but a decrease in fitness for females (Gerlach and Ketterson, 2013), providing direct experimental support for the hypothesis that there is sexual conflict over optimal testosterone levels in this species. As such, this is an ideal system in which to investigate the molecular mechanisms by which sexual conflict occurs and/or is resolved, by specifically asking whether the sexes diverge in the gene expression response to testosterone treatment.

Many sexually dimorphic and androgen-responsive phenotypes are mediated directly by changes in peripheral tissues such as liver and muscle. The liver plays a key role in whole-body metabolism, including gluconeogenesis, glycogenolysis, glycogen storage, amino acid synthesis, lipid synthesis and breakdown, and the production of insulin-like growth factor (Miura et al., 1992; Heubi, 1993). Further, the liver is a key regulator of sexually dimorphic immune function: male mice are more susceptible to liver infection than females (Diodato et al., 2001), and these differences are androgen-mediated (Mock and Nacy, 1988) through gene expression changes (Delić et al., 2010). Sex differences in gene expression in liver can be substantial (Corton et al., 2012), and are largely driven by activational effects of hormones (van Nas et al., 2009). The physiological demands of flight are thought to have resulted in a larger liver in birds compared with mammals (Proctor, 1993), making hormonal influences of this organ particularly important in birds.

Similarly, muscle tissues are also often sensitive to testosterone and play a primary role in mediating dimorphic behavior and physiology (Arnold et al., 1997; Baur et al., 2008; Fernando et al., 2010). Gene expression appears to account for many sexually dimorphic muscle features in humans (Maher et al., 2009; Welle et al., 2008) and mice (Yang et al., 2006). Androgen treatment leads to increases in strength and lean muscle mass (Hartgens and Kuipers, 2004), and these effects may be linked to testosterone-mediated changes in gene expression (Montano et al., 2007; Labrie et al., 2005). Further, the effects of exercise on gene expression in muscle are sex-specific in humans (Liu et al., 2010), suggesting that different transcriptional pathways may underlie some of the sex differences in muscle. The pectoralis muscle, which is the major avian flight muscle, accounts for ~20% of the mass of an individual bird (Marden, 1987). Androgen receptor is expressed in the pectoralis (Feng et al., 2010), and testosterone modifies the expression of at least two candidate genes related to muscle function in the pectoralis (Fuxjager et al., 2012). Thus, the pectoralis provides an important, androgen-sensitive tissue in which to investigate the sex-specific effects of hormones in the periphery.

We anticipated that many of the genes differentially expressed between sexes and in response to testosterone treatment in the liver and the pectoralis would have functions related to metabolism, muscle development and immune function. Because many of the metabolic effects of testosterone are similar in male and female juncos, we also predicted that many genes whose expression was altered in response to testosterone treatment in one sex would also be altered in the other sex. However, we also predicted that some genes would respond to testosterone treatment in one sex, but not the other, consistent with previous findings that (1) not all physiological effects of testosterone are present in both sexes (Ketterson et al., 2005) and (2) the sexes respond differently to testosterone treatment at the level of gene transcription in the brain (Peterson et al., 2013) providing a possible solution to the sexual conflict over testosterone levels observed in previous studies on free-living juncos (Gerlach and Ketterson, 2013; Reed et al., 2006).

Sex differences

We identified significant differences in expression between control males and control females in both the liver and the pectoralis. In the liver, 291 genes (of 12,206 expressed) were differentially expressed between control males and females (Fig. 1A, supplementary material Table S1), including 218 that were more highly expressed in males than in females and 73 that were more highly expressed in females than in males. Among these genes, nine gene ontology (GO) terms were significantly over-represented (Table 1).

Fig. 1.

Sex differences in gene expression. Differences in gene expression between the sexes are represented by heat maps that show scaled individual expression scores (Z-scores) for significantly differentially expressed genes in the liver (A) and pectoralis (B). Yellow represents high gene expression; blue represents low expression scaled to the levels of expression for each gene. (C) Venn diagram shows the overlap in significant genes between the two tissues. Each column represents an individual, and each row a gene.

Fig. 1.

Sex differences in gene expression. Differences in gene expression between the sexes are represented by heat maps that show scaled individual expression scores (Z-scores) for significantly differentially expressed genes in the liver (A) and pectoralis (B). Yellow represents high gene expression; blue represents low expression scaled to the levels of expression for each gene. (C) Venn diagram shows the overlap in significant genes between the two tissues. Each column represents an individual, and each row a gene.

Table 1.

Gene ontology (GO) terms over-represented among genes differentially expressed in the liver between male and female dark-eyed juncos

Gene ontology (GO) terms over-represented among genes differentially expressed in the liver between male and female dark-eyed juncos
Gene ontology (GO) terms over-represented among genes differentially expressed in the liver between male and female dark-eyed juncos

In the pectoralis, 658 genes (of 11,465 expressed) were differentially expressed between control males and females (Fig. 1B, supplementary material Table S1), including 450 that were more highly expressed in males than in females and 208 that were more highly expressed in females than in males. Among these genes, 18 GO terms were significantly over-represented (Table 2).

Among the genes differentially expressed between the sexes, 117 were significantly different in both liver and pectoralis (Fig. 1C, supplementary material Table S1). Of these genes, 91 were higher in control males than control females in both tissues, and 25 were higher in control females than control males in both tissues. Only one gene was differentially expressed by sex in opposite directions in the two tissues: protein tyrosine phosphatase, receptor type C was higher in control males than control females in the liver, but higher in control females than control males in the pectoralis. The general patterns of gene expression by sex were largely consistent between the two tissues. That is, genes that were more highly expressed in males than females in one tissue tended to be more highly expressed by males than females in the other tissue, and vice versa, more than expected by chance (Fisher's exact test, P<0.0001).

Effect of testosterone treatment in females

In both liver and pectoralis, we identified significant differences in expression between control females and testosterone-treated females. In the liver, 801 genes (of 12,064 expressed) were differentially expressed (Fig. 2A, supplementary material Table S1), including 645 that were expressed at a higher level in testosterone-treated females than controls and 156 that were expressed at a lower level in testosterone-treated females than controls. Among these genes, 26 GO terms were over-represented (Table 3).

In the pectoralis, 402 genes (of 11,413 expressed) were differentially expressed between control females and testosterone-treated females (Fig. 2B, supplementary material Table S1), including 226 that were expressed at a higher level in testosterone-treated females than controls and 174 that were expressed at a lower level in testosterone-treated females than controls. Among these genes, 17 GO terms were over-represented (Table 4).

Among the genes differentially expressed between the testosterone-treated and control females, 40 were significantly different in both liver and pectoralis (Fig. 2C, supplementary material Table S1). Of these genes, 21 were higher in testosterone-treated than control females in both tissues and seven were lower in testosterone-treated than control females in both tissues; 12 genes were differentially expressed in opposite directions in both tissues. More genes were affected in the same direction (i.e. either higher or lower in testosterone-treated than control females) in both tissues than expected by chance (Fisher's exact test, P<0.05, demonstrating significant similarity in the direction of gene expression change in response to testosterone treatment in the two tissues in females.

Table 2.

GO terms over-represented among genes differentially expressed in the pectoralis between male and female dark-eyed juncos

GO terms over-represented among genes differentially expressed in the pectoralis between male and female dark-eyed juncos
GO terms over-represented among genes differentially expressed in the pectoralis between male and female dark-eyed juncos
Fig. 2.

Gene expression in response to testosterone treatment in each sex. Differences in gene expression between testosterone-treated and control individuals in both the liver (left column) and the pectoralis (middle column) in females (A–C) and males (D–F). Heat maps show scaled individual expression scores (Z-scores) for genes that were significantly differentially expressed between testosterone treated and control individuals in each sex (A,B,D,E). Each column represents an individual, and each row a gene. Yellow represents high gene expression, and blue represents low expression scaled to the levels of expression for each gene. (C,F) Venn diagrams show the overlap of significant genes within each contrast between the tissues. See Results and supplementary material Table S1 for more information.

Fig. 2.

Gene expression in response to testosterone treatment in each sex. Differences in gene expression between testosterone-treated and control individuals in both the liver (left column) and the pectoralis (middle column) in females (A–C) and males (D–F). Heat maps show scaled individual expression scores (Z-scores) for genes that were significantly differentially expressed between testosterone treated and control individuals in each sex (A,B,D,E). Each column represents an individual, and each row a gene. Yellow represents high gene expression, and blue represents low expression scaled to the levels of expression for each gene. (C,F) Venn diagrams show the overlap of significant genes within each contrast between the tissues. See Results and supplementary material Table S1 for more information.

Effect of testosterone treatment in males

In the liver, 283 genes (of 12,229 expressed) were differentially expressed between testosterone-treated and control males (Fig. 2D, supplementary material Table S1), including 99 that were expressed at a higher level in testosterone-treated males than controls and 184 that were expressed at a lower level in testosterone-treated males than controls. Among these genes, one GO term was over-represented: acetylglucosaminyltransferase activity.

In the pectoralis, 450 genes (of 11,282 expressed) were differentially expressed between control males and testosterone-treated males (Fig. 2E, supplementary material Table S1), including 148 that were expressed at a higher level in testosterone-treated males than controls and 302 that were expressed at a lower level in testosterone-treated males than controls. Among these genes, eight GO terms were over-represented (Table 5).

Table 3.

GO terms over-represented among genes differentially expressed in the liver between testosterone-treated and control female dark-eyed juncos

GO terms over-represented among genes differentially expressed in the liver between testosterone-treated and control female dark-eyed juncos
GO terms over-represented among genes differentially expressed in the liver between testosterone-treated and control female dark-eyed juncos
Table 4.

GO terms over-represented among genes differentially expressed in the pectoralis between testosterone-treated and control female dark-eyed juncos

GO terms over-represented among genes differentially expressed in the pectoralis between testosterone-treated and control female dark-eyed juncos
GO terms over-represented among genes differentially expressed in the pectoralis between testosterone-treated and control female dark-eyed juncos

Among the genes differentially expressed between the testosterone-treated and control males, 21 were significantly different in both liver and pectoralis (Fig. 2F, supplementary material Table S1). Of these genes, six were higher in testosterone-treated than control males in both tissues and 10 were lower in testosterone-treated than control males in both tissues; five genes were differentially expressed in opposite directions in the tissues. More genes were affected by testosterone treatment in the same direction (i.e. either higher or lower in testosterone-treated than control males) in both tissues than expected by chance (Fisher's exact test; P<0.05), suggesting similar changes in response to testosterone treatment in the two tissues in males.

Effect of testosterone treatment in both sexes

In both liver and pectoralis, some genes were differentially expressed between testosterone-treated and control individuals of both sexes, though many genes were significantly differently affected in the two sexes (i.e. had a significant interaction effect). In the liver, 58 genes were differentially expressed in both sexes, representing only 5.6% of the 1026 genes differentially expressed in at least one sex. There was a significant interaction between sex and the effect of testosterone treatment in the liver for 550 genes, including 366 (38%) of the genes that were significantly affected by testosterone treatment in only one sex (Fig. 3A, supplementary material Table S1).

In the pectoralis, 68 genes were differentially expressed between testosterone-treated and control individuals of both sexes, representing only 8.7% of the 784 genes that were differentially expressed in at least one sex. There was a significant interaction between sex and the effect of testosterone treatment in the pectoralis for 297 genes, including 189 (26%) of the genes that were only significantly affect by testosterone treatment in one sex (Fig. 3B, supplementary material Table S1).

In the liver, the genes differentially expressed between testosterone-treated and control individuals in both sexes include 28 that were expressed at a higher level in testosterone-treated individuals than controls in both sexes, 11 that were expressed at a lower level in testosterone-treated individuals than controls in both sexes, and 19 that were differentially expressed in opposite directions in the two sexes (Table 6). That is, 67% of genes differentially expressed by testosterone treatment in both sexes were differentially expressed in the same direction (i.e. either higher or lower in testosterone-treated than control individuals in both sexes), more than expected by chance (Fisher's exact test, P<0.05).

In the pectoralis, the genes differentially expressed between testosterone-treated and control individuals in both sexes included 34 that were expressed at a higher level in testosterone-treated individuals than controls in both sexes, 27 that were expressed at a lower level in testosterone-treated individuals than controls in both sexes, and seven that were differentially expressed in opposite directions in the two sexes (Table 6). That is, 90% of genes differentially expressed by testosterone treatment in both sexes were differentially expressed in the same direction (i.e. either higher or lower in testosterone-treated than control individuals in both sexes), more than expected by chance (Fisher's exact test, P<0.0001).

Table 5.

GO terms over-represented among genes differentially expressed in the pectoralis between testosterone-treated and control male dark-eyed juncos

GO terms over-represented among genes differentially expressed in the pectoralis between testosterone-treated and control male dark-eyed juncos
GO terms over-represented among genes differentially expressed in the pectoralis between testosterone-treated and control male dark-eyed juncos
Fig. 3.

Comparing the effect of testosterone treatment in males and females. Venn diagrams for (A) liver and (B) pectoralis showing the number of genes significantly differentially expressed between testosterone-treated and control individuals in males and females, and those with a significant sex-by-treatment interaction effect.

Fig. 3.

Comparing the effect of testosterone treatment in males and females. Venn diagrams for (A) liver and (B) pectoralis showing the number of genes significantly differentially expressed between testosterone-treated and control individuals in males and females, and those with a significant sex-by-treatment interaction effect.

Using a microarray specific to the dark-eyed junco, we identified a large number of genes that were expressed differentially between males and females, and between testosterone-treated and control individuals of each sex, in the liver and pectoralis. As predicted, many of the differentially expressed genes were functionally related to previously described phenotypic effects of testosterone treatment as well as known sexual dimorphisms. Testosterone treatment tended to affect different genes in males and females; however, among the genes differentially expressed by testosterone treatment in both sexes, testosterone treatment affected most genes in the same direction in males and females. This suggests that sexually dimorphic transcriptional responses to testosterone may provide one solution to sexual conflict over circulating levels of testosterone. Not only do these results provide a detailed view of the molecular mechanisms by which sexual conflict may be resolved, but they also lay a strong foundation for ecologically relevant and evolutionarily significant advances in our understanding of the mechanisms underlying life-history trade-offs and behavioral evolution in natural systems, such as the junco. Furthermore, by focusing on the liver and muscle, our findings point to the mechanisms by which sexual dimorphic peripheral responses to circulating hormones may play a role in sexual conflict and dimorphism, in addition to previously identified effects in the brain of juncos (Peterson et al., 2013) and the sex-specific effects previously identified in rats (van Nas et al., 2009; Yang et al., 2006).

Sexually dimorphic gene expression

Similarly to previous studies on neural tissues in juncos (Peterson et al., 2013) and multiple tissues in other species (reviewed in Ellegren and Parsch, 2007), we identified many genes that were expressed differentially between males and females. In the pectoralis, we identified 658 genes that were sexually dimorphic, and as predicted, GO analysis revealed over-representation of terms related to muscle development, including muscle system process and both the I band and A band portions of the sarcomere. These genes were generally regulated in directions consistent with known sex differences in body mass in the junco (Nolan et al., 2002). For example, titin, a gene that regulates muscle elasticity (Itoh-Satoh et al., 2002), is expressed at a higher level in control males than females in the pectoralis. SMAD-related protein 2 was more highly expressed in the pectoralis of control females than control males, consistent with the known role of SMAD proteins in reducing cellular growth (Nakao et al., 1997). Consistent with other studies comparing transcriptional patterns in skeletal muscle of males and females (Yang et al., 2006; Roth et al., 2002; Welle et al., 2008), we found a large number of genes that differed in expression between the sexes, including several genes that were directly related to muscle development and growth.

We also identified 291 genes that showed significantly different expression between control males and control females in the liver, and several of the differentially expressed genes were related to known phenotypic differences between the sexes. For example, Lipid phosphate phosphohydrolase 1, a gene involved in glycerolipid synthesis and lipid uptake (Kai et al., 1997), was expressed at a higher level in control males than females, consistent with sex differences in metabolic activity (Fernando et al., 2010; Wikelski et al., 1999). Further, the GO term steroid binding was over-represented among these genes. For example, hydroxysteroid dehydrogenase like 2 (HSDL2) was expressed more highly in the liver of control males than females. HSDL2 plays a role in sterol binding (Dai et al., 2003), and is marginally more highly expressed in the liver of females than males in mice (Gatti et al., 2010). This suggests that the sexes might differ in their metabolism of sterol-based compounds in the liver, but that this difference may vary between taxa.

Among genes differentially expressed between control males and control females in the pectoralis were several transcription factors. In both liver and pectoralis, transcription factor III B 150 (TFIIIB150) was more highly expressed in control males than control females. TFIIIB150 mediates transcription via RNA polymerase III (Schramm et al., 2000), which is primarily involved in the expression of 5S rRNA, tRNA and other small RNAs (Dieci et al., 2007). In contrast, Basic Transcription Factor 3 (BTF3) was expressed more in control females than males in the liver, and activated RNA polymerase II transcriptional coactivator p15 (P15) was expressed more in control females than males in the pectoralis. Both BTF3 (Zheng et al., 1990) and P15 (Kretzschmar et al., 1994) activate expression via RNA polymerase II, which is the primary polymerase for the expression of protein-coding genes (Sims et al., 2004). These findings suggest that males and females may orchestrate gene expression differently, with males favoring expression of housekeeping-type genes that may increase translation rates, and females favoring expression of protein-coding genes. Transcription factors are among the genes that are differentially expressed by sex in human muscles (Roth et al., 2002), and BTF3 and P15 are differentially expressed by sex in the liver and muscle of mice (Yang et al., 2006). This suggests that transcription factors in general, and BTF3 and P15 in particular, may be involved in sexually dimorphic patterns of expression in many species.

Table 6.

Comparing gene expression in response to testosterone treatment in male and female dark-eyed juncos

Comparing gene expression in response to testosterone treatment in male and female dark-eyed juncos
Comparing gene expression in response to testosterone treatment in male and female dark-eyed juncos

Effect of testosterone treatment in females

Testosterone implants affected female gene expression in both tissues, and the effects were consistent with known phenotypic effects of testosterone treatment. For example, Immunoglobulin A (IgA) heavy chain had lower expression in the liver of testosterone-treated females than controls, and AF411388_1 basic, a gene containing a conserved immunoglobulin region (Yoder et al., 2002), was expressed at a lower level in the pectoralis of testosterone-treated females than controls. Immunoglobulins play a major role in immune function (Litman et al., 1993), so their lower expression in testosterone-treated females is consistent with the known suppressive effect of testosterone treatment on immune function in female juncos (Zysling et al., 2006).

Additionally, the GO term growth was over-represented among the genes differentially expressed in the pectoralis between testosterone-treated and control females. Eleven of the 17 genes annotated as growth were more highly expressed in testosterone-treated females than controls, and the other six were expressed at a lower level in testosterone-treated females than controls. Three of the genes that were expressed at a lower level are known repressors of growth [two representations of Ankyrin repeat domain-containing protein 26 (Bera et al., 2008)] or transcription [B-cell CLL/lymphoma 6 (zinc finger protein 51) (Lemercier et al., 2002)]. Both higher expression of growth promoting genes and lower expression of growth repressors are consistent with the role of elevated androgens in increasing muscular growth and maintenance (Woodward, 1993; Hartgens and Kuipers, 2004).

A number of the genes identified as differentially expressed between testosterone-treated and control females were similar to those identified in studies of other organisms. The GO term response to hormone stimulus consists of genes identified as mediators of phenotypic effects of hormones in other species (Ashburner et al., 2000), and was over-represented among differentially expressed genes in female pectoralis. Among these genes, Serotonin 1B receptor was expressed at a higher level in testosterone-treated than control females. Serotonin 1B receptor expression is upregulated by mineralcorticoids in the aorta of rats (Banes and Watts, 2002); thus, expression of this serotonin receptor may be mediated indirectly by testosterone treatment through changes in other signaling molecules. Serotonin 1B receptor has a range of effects on both behavior and physiology (Donaldson et al., 2013), though its role in skeletal muscle tissue is unclear. In addition, carbonic anhydrase II, a catalyst of the hydrolysis of carbon dioxide (Sterling et al., 2001), was expressed at lower levels in the pectoralis of testosterone-treated females than controls. Expression of a related gene, carbonic anhydrase III, is reduced by strength training in humans (Roth et al., 2002), suggesting that the action of testosterone treatment may be related to changes in muscle activity. Expression of carbonic anhydrase II is also reduced by exposure to estrogens in some tissues in rats (Caldarelli et al., 2005), consistent with the view that some of the effects seen in our study may be mediated by conversion of testosterone to estradiol. We anticipate that many of the genes we have identified play a role in mediating tissue-level responses to hormones in multiple species, patterns that will become clear in time.

Effect of testosterone treatment in males

We identified a large number of genes that were differentially expressed between testosterone-treated and control males in both liver and pectoralis, and several of them are related to known phenotypic effects of testosterone treatment. For example, heme oxygenase (decyclizing) 1 (HMOX1) was expressed at lower levels in the liver of testosterone-treated than control males. HMOX1 is a key enzyme in the breakdown of heme (Platt and Nath, 1998), and has been implicated in the disruption of human glucose regulation (Bao et al., 2012). Therefore, HMOX1's lower expression in testosterone-treated males is consistent with previous findings that testosterone increases metabolism (Oppliger et al., 2004; Fernando et al., 2010), as well as other studies that have linked heme-related enzymes with activational effects of androgens (van Nas et al., 2009). Aldehyde oxidase 1 (AOX1) was also expressed at lower levels in the liver of testosterone-treated males than controls. Aldehyde oxidases break down a number of metabolically active compounds (Hartmann et al., 2012). So the lower expression of AOX1 may indicate that testosterone treatment reduced catabolism in the liver, consistent with a previous study showing that castrated mice treated with androgens also showed significant changes in expression of a variety of metabolic genes in the liver (van Nas et al., 2009).

Previous studies in humans (Michael et al., 2005), rats (Wakley et al., 1991) and chickens (Pederson et al., 1999) have demonstrated that higher testosterone reduces bone resorption, though several of the effects may be related to the conversion of testosterone to estradiol (Oursler et al., 1991). Consistent with these findings, we observed that osteoclast inhibitory lectin, which blocks the formation of bone-resorption osteoclasts (Hu et al., 2004), was expressed at higher levels in testosterone-treated than control males in both the liver and the pectoralis.

Many of the effects of testosterone treatment that we have identified have the potential to play large, downstream roles, as evidenced by the over-representation of the GO term rRNA metabolic process among genes differentially expressed in the pectoralis. For example, Serine/arginine-rich splicing factor 5 modulates the splice variant selection of many genes (Sebbag-Sznajder et al., 2012) and thus plays a large role in cellular function. This gene was more highly expressed in both the liver and pectoralis of testosterone-treated than control males. In the pectoralis, MGC89063 was more highly expressed in testosterone-treated males than controls, which is similar to what was found in the hypothalamus and medial amygdala of the junco (Peterson et al., 2013). MGC89063 is a transcription factor (Ashburner et al., 2000; Hunter et al., 2009), and the fact that it was more highly expressed in the hypothalamus, medial amygdala and pectoralis (but not liver) of testosterone-treated males than controls, and not differentially expressed by testosterone-treatment in females, raises the possibility that MGC89063 may play a tissue- and sex-specific role in meditating the effects of testosterone treatment via downstream gene regulation. The specific downstream effects of this gene remain unclear, but given its role in multiple target tissues, further investigation into the pleiotropic roles of MGC89063 will likely provide novel insights into the integrated response to testosterone treatment. Continued focus on non-model organisms such as the junco in these future studies may provide greater insight into the fitness consequences of genes like these.

Effect of testosterone treatment in both sexes

Many genes were differentially expressed in the liver and pectoralis between testosterone-treated and control individuals of both sexes (63 genes in liver and 70 genes in pectoralis). However, this number represents only 5% of the genes differentially expressed by testosterone treatment in either sex, meaning that 95% of genes that were affected by testosterone were not significantly affected in both sexes. Further, in each tissue, over a quarter of these genes were affected significantly differently in each sex (i.e. had a significant interaction effect), suggesting that many, though not all, of the genes identified in only one sex are truly only affected in that sex. This result, especially when combined with similar findings in neural tissue in juncos (Peterson et al., 2013), lends some support to the hypothesis that testosterone treatment leads to transcriptional changes in largely different genes in the two sexes, and suggests a possible remedy to sexual conflict over testosterone levels. In contrast, among those genes that were differentially expressed in both males and females, most (78%) were differentially expressed in the same direction in both sexes. Collectively, these results suggest that there may be a core transcriptional response to testosterone treatment shared between the sexes, but this response is fine-tuned by sex-specific responses, which may reduce sexual conflict over circulating testosterone levels.

Among the genes that were significantly differentially expressed between testosterone-treated and control individuals in both sexes, several relate to the known effects of testosterone treatment on activity and metabolism (Wikelski et al., 1999; Lynn et al., 2000; Buchanan et al., 2001). For example, in the liver and pectoralis of both sexes, L-arginine:glycine amidinotransferase was more highly expressed in testosterone-treated individuals than in controls. This gene encodes the enzyme for the rate-limiting step in creatine biosynthesis (Humm et al., 1997), which in turn increases energy availability in muscle (Kraemer and Volek, 1999), and is also regulated by steroid hormones in rodents (Kriskó and Walker, 1966). Therefore, greater expression of L-arginine:glycine amidinotransferase in testosterone-treated individuals is consistent with steroid-induced increases in activity levels and metabolic rate. Similarly, 3-hydroxybutyrate dehydrogenase (3HBDH) was expressed at a higher level in testosterone-treated individuals in both sexes and both tissues. 3HBDH catalyzes the reversible reaction between beta-hydroxybutyric acid and acetoacetate, a key step in the breakdown of fatty acids for energy (Bergmeyer et al., 1967; Williamson et al., 1962). Both male and female rats respond to androgen treatment with changes in the expression of fatty acid metabolizing genes as well (van Nas et al., 2009). Together, the changes in the expression of these genes could be a major contributor to testosterone-induced shifts in metabolism and activity in juncos (Chandler et al., 1994; Lynn et al., 2000) and other species (Wikelski et al., 1999; Marler et al., 1995). However, it remains possible that these changes in gene expression are indirect effects of testosterone treatment, e.g. if testosterone affects metabolism or activity via other routes and these genes respond in kind to altered metabolism or activity.

Several genes related to insulin signaling were differentially expressed between testosterone-treated and control individuals of both sexes in the liver. Insulin receptor substrate 4 (IRS4), for example, was expressed more highly in the livers of testosterone-treated individuals than controls in both sexes. IRS4 mediates the activity of a number of growth factors (e.g. Hinsby et al., 2004), and the lack of IRS4 leads to a decrease in body size in knockout mice (Fantin et al., 2000). Therefore, higher expression of IRS4 in the liver of testosterone-treated individuals than controls may mediate some of the previously reported metabolic and growth effects of testosterone treatment (Cox et al., 2009; Wikelski et al., 1999; Lynn et al., 2000). In addition, insulin-like growth factor 2 receptor was more highly expressed in the liver of testosterone-treated females than controls, and insulin-like growth factor 1 was expressed at lower levels in the liver of testosterone-treated males than controls. However, neither gene had a significant sex-by-treatment interaction term, suggesting that both genes may also have been regulated in the opposite sex, but below our limits of detection. Insulin-like growth factors also mediate growth (Abuzzahab et al., 2003; Petry et al., 2005) and have been implicated in the expression of sexually selected traits (Emlen et al., 2012), some of which are also meditated by androgens (Folstad and Karter, 1992).

Several genes related to the regulation of growth were differentially expressed in the liver between testosterone-treated and control individuals of each sex, though some of the specific genes affected by testosterone treatment differed between males and females. Follistatin was more highly expressed in the liver of testosterone-treated individuals than controls in both sexes. Follistatin binds and inactivates members of the TGF-beta super family, including myostatin, such that increased follistatin is associated with increased muscle growth (Lee and McPherron, 2001). Further, epidermal growth factor receptor (EGFR) was more highly expressed in the liver of testosterone-treated than control females, and had a marginally significant sex-by-treatment interaction term (uncorrected P=0.04), suggesting some sex-specific hormone regulation. Likewise, opioid growth factor receptor (OGFr) was less expressed in the liver of testosterone-treated than control males and also had a significant sex-by-treatment interaction term. EGFR acts to increase cell proliferation and growth (Oda et al., 2005), but OGFr acts to reduce growth (Zagon et al., 2008). In men, similar changes in the expression of growth-related genes, including OGFr, are observed in response to hormone manipulation, and the changes are believed to be related to lean muscle mass growth (Montano et al., 2007). Thus, the expression changes seen in both male and female juncos are related to increased growth, but potentially via different transcriptional mechanisms.

Conclusions

In this study, we applied genomic tools to the dark-eyed junco in order to identify ecologically relevant sex differences in gene expression and transcriptional responses to experimentally elevated testosterone. As predicted, many of the specific genes affected were associated with known physiological and metabolic effects of testosterone treatment, but the expression response to testosterone treatment was different in the two sexes: only 5% of regulated genes overlap in the two sexes. Interestingly, among genes that were differentially expressed between testosterone-treated and control individuals in both sexes, most were differentially expressed in the same direction. Therefore, testosterone may be utilizing a shared core set of transcriptional paths in both sexes that are complemented and modified by sex-specific transcriptional responses. Characterizing these effects in the periphery is particularly notable in light of the prevailing view in behavioral neuroendocrinology that many sex differences are mediated at the level of the brain. Our results detail some of the molecular mechanisms by which hormones have sex-specific activational effects in two important peripheral tissues. Whether these sex-specific mechanisms represent adaptive mechanistic responses to testosterone is an open question that can be addressed by continued focus on ecological model species such as the junco.

Animal collection and treatment

Adult dark-eyed juncos (14 male, 12 female) from near Mountain Lake Biological Station (Pembroke, VA, USA; 37°22′31″N, 80°31′24″W) were captured, held in a semi-naturalist aviary, and treated as described in a previous study analyzing neural tissues (Peterson et al., 2013). Briefly, testosterone-treated individuals were implanted with silastic tubing filled with crystalline testosterone (males: two 10 mm implants; females: one 5 mm implant; Sigma-Aldrich, St Louis, MO, USA), and control individuals were implanted with one 10 mm empty implant. These testosterone implants result in levels of testosterone near the physiological maximum in each sex (Ketterson et al., 2005). Thus, while all animals had testosterone levels above any threshold necessary to maintain reproductive physiology and behavior, animals given testosterone implants had testosterone levels that were at the high end of natural variation. Notably, this implant regiment has repeatedly been shown to affect many different phenotypes in male and female juncos, and experimental treatment with these testosterone implants reveals that there is sexual conflict over testosterone levels in this system – above-average testosterone levels are selectively advantageous for males and disadvantageous for females (see Gerlach and Ketterson, 2013; Ketterson et al., 2009; Reed et al., 2006) (summarized in Introduction).

We note that direct and indirect mechanisms of action and interaction with natural hormones are important to consider when evaluating our results for two reasons. First, our implants used testosterone, which can be aromatized into estradiol, and thus several of the effects described here may be mediated directly by estradiol after local conversion via aromatase (Herbst and Bhasin, 2004). These sex steroids may act directly on muscle tissue, and they also may directly alter activity, metabolism or other aspects of behavior and physiology that lead to indirect effects on gene expression in liver or muscle (Park et al., 2012).

Second, we used intact animals in breeding condition to ensure that seasonally variable aspects of behavior and physiology were characteristic of the breeding season, mimicking previous studies that have demonstrated sex differences in the phenotypic and fitness consequences of testosterone in otherwise normal breeding birds. Importantly, the effects seen here likely reflect the mechanisms of action in previous implant studies (e.g. Ketterson et al., 1996; Gerlach and Ketterson, 2013; Reed et al., 2006; Clotfelter et al., 2004) as well as those that might occur in response to evolutionary increases in testosterone levels (Ketterson et al., 2009).

After 26 days of exposure to implants, individuals were euthanized by overdose of isoflurane, and tissues were collected rapidly (within 15 min) to ensure minimal RNA degradation (Cheviron et al., 2011). Approximately 2 cm3 from the tip of the right lobe of the liver and approximately 1 cm3 from near the midline of the pectoralis muscle were collected from all individuals. Other collected tissues remain available for future analyses. Sexes and treatments were balanced across day and time of euthanasia (between 07:00 and 12:30 h). All animal methods were reviewed and approved by the Institutional Animal Care and Use Committee at Indiana University Bloomington (protocol no. 09-037).

cDNA preparation and hybridization

Microarray experiments were conducted as described previously (Peterson et al., 2012; Peterson et al., 2013) following Lopez and Colbourne (Lopez and Colbourne, 2011). RNA from liver and pectoralis was extracted in TRIzol following the manufacturer's instructions (Invitrogen, Carlsbad, CA, USA). All extracted RNA was assessed on an Agilent Bioanalyzer (Santa Clara, CA, USA) and showed high quality: RNA integrity number (Schroeder et al., 2006) scores ranged from 6.7 to 9.2. We then performed double-stranded cDNA synthesis with the Invitrogen SuperScript Double-Stranded cDNA Synthesis Kit with labeled cDNA using 1 optical density CY-labeled random nonamer primer (either Cy3 or Cy5) and random hexamer primers and 100 U Klenow fragment per 1 μg double-stranded cDNA (following NimbleGen labeling protocols).

A full round robin design was used for each tissue (N=6 per treatment group for each tissue). Each sample was tested once, and each treatment group was hybridized against each other group twice (once with each dye direction; supplementary material Fig. S1). Thus, 15 μg of two labeled samples (one Cy3, one Cy5) were hybridized to each sub-array of a custom Nimblegen 12-plex microarray (Roche Nimblegen, Madison, WI, USA) for the dark-eyed junco containing 100,635 features representing 33,545 contigs (assembled sequencing reads) in triplicate covering 22,765 isogroups (putative genes) based on transcriptome sequencing (Peterson et al., 2012). Post-hybridization washing and scanning followed the manufacturer's directions (Roche NimbleGen). An Axon GenePix 4200A scanner (Molecular Devices, Sunnyvale, CA, USA) with GenePix 6.0 software captured array images and NimbleScan 2.4 (Roche NimbleGen) was used to extract data. We then used the limma package (Smyth, 2005) in R (R Development Core Team, 2010) to process and normalize raw microarray data. Microarray data are available in the NCBI Gene Expression Ominubus repository (Accession number GSE41076).

Microarray analysis

Three comparisons for each tissue were made using limma (Smyth, 2005) – (1) control males versus control females, (2) control males versus testosterone-treated males and (3) control females versus testosterone-treated females (N=6 per treatment group for each tissue) – as well as the interaction between testosterone and sex. Only contigs that were expressed in at least one of the compared treatment groups were analyzed [identified as described in Peterson et al. (Peterson et al., 2012)]. Briefly, a gene was considered expressed if at least half of the individuals in a treatment group had expression scores greater than 97.5% of the random probes on the array.

In most isogroups, the log2 fold changes between treatment groups, along with the modified t-statistic and P-value, calculated in the limma package were used for calculations, statistics and visualization. However, for isogroups represented by more than one contig (4288 of 22,765 isogroups), we calculated the mean t-value of all contigs, and calculated significance on degrees of freedom equal to the total number of probes scored for the isogroup minus two. The median fold change from contigs was assigned to each isogroup. We used the R package qvalue (Storey, 2002) to calculate q-values using a global (across all eight contrasts) false discovery threshold of 0.05 (Benjamini and Hochberg, 1995). To further assess similarity in the effects of testosterone treatment between males and females, the direction of gene expression difference between comparisons was examined using a Fisher's exact test on genes that were differentially expressed between testosterone-treated and control individuals in both sexes.

We then used topGO (Alexa and Rahnenfuhrer, 2010) with the weight algorithm (Alexa et al., 2006) to identify the gene ontology (GO) terms (Ashburner et al., 2000) that were significantly over-represented among the significantly differentially expressed genes in each comparison. Because we analyzed all three GO topologies, we used a Bonferroni-corrected P-value cut-off of 0.0125. GO terms with fewer than five annotations were excluded from the analysis, and only terms with at least three genes in the significant gene set are reported.

The authors thank Sarah Wanamaker for assistance with animal care, Baiju Parikh, at Roche NimbleGen, for contributions to the CGB Ecological Genomics Pipeline, and Mountain Lake Biological Station, University of Virginia, for facilities and access to property for conducting this research.

Funding

This material is based upon work supported by the National Science Foundation [ACI-0338618l, CNS-0521433, ABI-1062432, OCI-0451237, OCI-0535258, OCI-0504075, IOS-0820055 and DGE-0504627]; the National Institutes of Health [T32-HD049336]; and the Indiana METACyt Initiative. The Indiana METACyt Initiative of Indiana University is supported in part by Lilly Endowment, Inc. This work was supported in part by Shared University Research grants from IBM to Indiana University. Deposited in PMC for release after 12 months.

Abuzzahab
M. J.
,
Schneider
A.
,
Goddard
A.
,
Grigorescu
F.
,
Lautier
C.
,
Keller
E.
,
Kiess
W.
,
Klammt
J.
,
Kratzsch
J.
,
Osgood
D.
, et al. 
Intrauterine Growth Retardation (IUGR) Study Group
(
2003
).
IGF-I receptor mutations resulting in intrauterine and postnatal growth retardation
.
N. Engl. J. Med.
349
,
2211
-
2222
.
Alexa
A.
,
Rahnenfuhrer
J.
(
2010
).
topGO: Enrichment Analysis for Gene Ontology. R package version 2.8.0
. Available at http://cran.r-project.org.
Alexa
A.
,
Rahnenfuehrer
J.
,
Lengauer
T.
(
2006
).
Improved scoring of functional groups from gene expression data by decorrelating GO graph structure
.
Bioinformatics
22
,
1600
-
1607
.
Arnold
A. P.
(
1975
).
The effects of castration and androgen replacement on song, courtship, and aggression in zebra finches (Poephila guttata)
.
J. Exp. Zool.
191
,
309
-
325
.
Arnold
A. M.
,
Peralta
J. M.
,
Thonney
M. L.
(
1997
).
Effect of testosterone on differential muscle growth and on protein and nucleic acid concentrations in muscles of growing lambs
.
J. Anim. Sci.
75
,
1495
-
1503
.
Ashburner
M.
,
Ball
C. A.
,
Blake
J. A.
,
Botstein
D.
,
Butler
H.
,
Cherry
J. M.
,
Davis
A. P.
,
Dolinski
K.
,
Dwight
S. S.
,
Eppig
J. T.
, et al. 
The Gene Ontology Consortium
(
2000
).
Gene ontology: tool for the unification of biology
.
Nat. Genet.
25
,
25
-
29
.
Banes
A. K. L.
,
Watts
S. W.
(
2002
).
Upregulation of arterial serotonin 1B and 2B receptors in deoxycorticosterone acetate-salt hypertension
.
Hypertension
39
,
394
-
398
.
Bao
W.
,
Rong
S.
,
Zhang
M.
,
Yu
X.
,
Zhao
Y.
,
Xiao
X.
,
Yang
W.
,
Wang
D.
,
Yao
P.
,
Hu
F. B.
, et al. 
. (
2012
).
Plasma heme oxygenase-1 concentration in relation to impaired glucose regulation in a non-diabetic Chinese population
.
PLoS ONE
7
,
e32223
.
Baur
L. A.
,
Nasipak
B. T.
,
Kelley
D. B.
(
2008
).
Sexually differentiated, androgen-regulated, larynx-specific myosin heavy-chain isoforms in Xenopus tropicalis; comparison to Xenopus laevis
.
Dev. Genes Evol.
218
,
371
-
379
.
Benjamini
Y.
,
Hochberg
Y.
(
1995
).
Controlling the false discovery rate – a practical and powerful approach to multiple testing
.
J. R. Stat. Soc. B
57
,
289
-
300
.
Bera
T. K.
,
Liu
X. F.
,
Yamada
M.
,
Gavrilova
O.
,
Mezey
E.
,
Tessarollo
L.
,
Anver
M.
,
Hahn
Y.
,
Lee
B.
,
Pastan
I.
(
2008
).
A model for obesity and gigantism due to disruption of the Ankrd26 gene
.
Proc. Natl. Acad. Sci. USA
105
,
270
-
275
.
Berg
E.
,
Maklakov
A.
(
2012
).
Sexes suffer from suboptimal lifespan because of genetic conflict in a seed beetle
.
Proc. R. Soc. B
279
,
4296
-
4302
.
Bergmeyer
H. U.
,
Gawehn
K.
,
Klotzsch
H.
,
Krebs
H. A.
,
Williamson
D. H.
(
1967
).
Purification and properties of crystalline 3-hydroxybutyrate dehydrogenase from Rhodopseudomonas spheroides
.
Biochem. J.
102
,
423
-
431
.
Bonduriansky
R.
,
Chenoweth
S. F.
(
2009
).
Intralocus sexual conflict
.
Trends Ecol. Evol.
24
,
280
-
288
.
Boratyński
Z.
,
Koskela
E.
,
Mappes
T.
,
Oksanen
T. A.
(
2010
).
Sex-specific selection on energy metabolism – selection coefficients for winter survival
.
J. Evol. Biol.
23
,
1969
-
1978
.
Buchanan
K. L.
,
Evans
M. R.
,
Goldsmith
A. R.
,
Bryant
D. M.
,
Rowe
L. V.
(
2001
).
Testosterone influences basal metabolic rate in male house sparrows: a new cost of dominance signalling?
Proc. R. Soc. B
268
,
1337
-
1344
.
Caldarelli
A.
,
Diel
P.
,
Vollmer
G.
(
2005
).
Effect of phytoestrogens on gene expression of carbonic anhydrase II in rat uterus and liver
.
J. Steroid Biochem. Mol. Biol.
97
,
251
-
256
.
Casto
J. M.
,
Nolan
V.
 Jr
,
Ketterson
E. D.
(
2001
).
Steroid hormones and immune function: experimental studies in wild and captive dark-eyed juncos (Junco hyemalis)
.
Am. Nat.
157
,
408
-
420
.
Chandler
C. R.
,
Ketterson
E. D.
,
Nolan
V.
,
Ziegenfus
C.
(
1994
).
Effects of testosterone on spatial activity in free-ranging male dark-eyed juncso, Junco hyemalis
.
Anim. Behav.
47
,
1445
-
1455
.
Chapman
T.
,
Arnqvist
G.
,
Bangham
J.
,
Rowe
L.
(
2003
).
Sexual conflict
.
Trends Ecol. Evol.
18
,
41
-
47
.
Cheviron
Z. A.
,
Carling
M. D.
,
Brumfield
R. T.
(
2011
).
Effects of postmortem interval and preservation method on RNA isolated from field-preserved avian tissues
.
Condor
113
,
483
-
489
.
Clotfelter
E. D.
,
O'Neal
D. M.
,
Gaudioso
J. M.
,
Casto
J. M.
,
Parker-Renga
I. M.
,
Snajdr
E. A.
,
Duffy
D. L.
,
Nolan
V.
 Jr
,
Ketterson
E. D.
(
2004
).
Consequences of elevating plasma testosterone in females of a socially monogamous songbird: evidence of constraints on male evolution?
Horm. Behav.
46
,
171
-
178
.
Corton
J. C.
,
Bushel
P. R.
,
Fostel
J.
,
O'Lone
R. B.
(
2012
).
Sources of variance in baseline gene expression in the rodent liver
.
Mutat. Res.
746
,
104
-
112
.
Cox
R. M.
,
Calsbeek
R.
(
2009
).
Sexually antagonistic selection, sexual dimorphism, and the resolution of intralocus sexual conflict
.
Am. Nat.
173
,
176
-
187
.
Cox
R. M.
,
Stenquist
D. S.
,
Calsbeek
R.
(
2009
).
Testosterone, growth and the evolution of sexual size dimorphism
.
J. Evol. Biol.
22
,
1586
-
1598
.
Dai
J.
,
Xie
Y.
,
Wu
Q.
,
Wang
L.
,
Yin
G.
,
Ye
X.
,
Zeng
L.
,
Xu
J.
,
Ji
C.
,
Gu
S.
, et al. 
. (
2003
).
Molecular cloning and characterization of a novel human hydroxysteroid dehydrogenase-like 2 (HSDL2) cDNA from fetal brain
.
Biochem. Genet.
41
,
165
-
174
.
Delić
D.
,
Grosser
C.
,
Dkhil
M.
,
Al-Quraishy
S.
,
Wunderlich
F.
(
2010
).
Testosterone-induced upregulation of miRNAs in the female mouse liver
.
Steroids
75
,
998
-
1004
.
Dieci
G.
,
Fiorino
G.
,
Castelnuovo
M.
,
Teichmann
M.
,
Pagano
A.
(
2007
).
The expanding RNA polymerase III transcriptome
.
Trends Genet.
23
,
614
-
622
.
Diodato
M. D.
,
Knöferl
M. W.
,
Schwacha
M. G.
,
Bland
K. I.
,
Chaudry
I. H.
(
2001
).
Gender differences in the inflammatory response and survival following haemorrhage and subsequent sepsis
.
Cytokine
14
,
162
-
169
.
Donaldson
Z. R.
,
Nautiyal
K. M.
,
Ahmari
S. E.
,
Hen
R.
(
2013
).
Genetic approaches for understanding the role of serotonin receptors in mood and behavior
.
Curr. Opin. Neurobiol.
23
,
399
-
406
.
Ellegren
H.
,
Parsch
J.
(
2007
).
The evolution of sex-biased genes and sex-biased gene expression
.
Nat. Rev. Genet.
8
,
689
-
698
.
Emlen
D. J.
,
Warren
I. A.
,
Johns
A.
,
Dworkin
I.
,
Lavine
L. C.
(
2012
).
A mechanism of extreme growth and reliable signaling in sexually selected ornaments and weapons
.
Science
337
,
860
-
864
,
860-864
.
Fantin
V. R.
,
Wang
Q.
,
Lienhard
G. E.
,
Keller
S. R.
(
2000
).
Mice lacking insulin receptor substrate 4 exhibit mild defects in growth, reproduction, and glucose homeostasis
.
Am. J. Physiol. Endocrinol. Metab.
278
,
E127
-
E133
.
Feng
N. Y.
,
Katz
A.
,
Day
L. B.
,
Barske
J.
,
Schlinger
B. A.
(
2010
).
Limb muscles are androgen targets in an acrobatic tropical bird
.
Endocrinology
151
,
1042
-
1049
.
Fernando
S. M.
,
Rao
P.
,
Niel
L.
,
Chatterjee
D.
,
Stagljar
M.
,
Monks
D. A.
(
2010
).
Myocyte androgen receptors increase metabolic rate and improve body composition by reducing fat mass
.
Endocrinology
151
,
3125
-
3132
.
Folstad
I.
,
Karter
A. J.
(
1992
).
Parasites, bright males, and the immunocompetence handicap
.
Am. Nat.
139
,
603
-
622
.
Fuxjager
M. J.
,
Barske
J.
,
Du
S.
,
Day
L. B.
,
Schlinger
B. A.
(
2012
).
Androgens regulate gene expression in avian skeletal muscles
.
PLoS ONE
7
,
e51482
.
Gatti
D. M.
,
Zhao
N.
,
Chesler
E. J.
,
Bradford
B. U.
,
Shabalin
A. A.
,
Yordanova
R.
,
Lu
L.
,
Rusyn
I.
(
2010
).
Sex-specific gene expression in the BXD mouse liver
.
Physiol. Genomics
42
,
456
-
468
.
Gerlach
N. M.
,
Ketterson
E. D.
(
2013
).
Experimental elevation of testosterone lowers fitness in female dark-eyed juncos
.
Horm. Behav.
63
,
782
-
790
.
Hartgens
F.
,
Kuipers
H.
(
2004
).
Effects of androgenic-anabolic steroids in athletes
.
Sports Med.
34
,
513
-
554
.
Hartmann
T.
,
Terao
M.
,
Garattini
E.
,
Teutloff
C.
,
Alfaro
J. F.
,
Jones
J. P.
,
Leimkühler
S.
(
2012
).
The impact of single nucleotide polymorphisms on human aldehyde oxidase
.
Drug Metab. Dispos.
40
,
856
-
864
.
Herbst
K. L.
,
Bhasin
S.
(
2004
).
Testosterone action on skeletal muscle
.
Curr. Opin. Clin. Nutr. Metab. Care
7
,
271
-
277
.
Heubi
J. A.
(
1993
).
Liver and biliary systems
. In
Physiology
(ed.
Sperelakis
N.
,
Banks
R. O.
), pp.
631
-
644
.
New York, NY
:
Little, Brown and Company
.
Hinsby
A. M.
,
Olsen
J. V.
,
Mann
M.
(
2004
).
Tyrosine phosphoproteomics of fibroblast growth factor signaling: a role for insulin receptor substrate-4
.
J. Biol. Chem.
279
,
46438
-
46447
.
Hu
Y. S.
,
Zhou
H.
,
Myers
D.
,
Quinn
J. M. W.
,
Atkins
G. J.
,
Ly
C.
,
Gange
C.
,
Kartsogiannis
V.
,
Elliott
J.
,
Kostakis
P.
, et al. 
. (
2004
).
Isolation of a human homolog of osteoclast inhibitory lectin that inhibits the formation and function of osteoclasts
.
J. Bone Miner. Res.
19
,
89
-
99
.
Humm
A.
,
Fritsche
E.
,
Steinbacher
S.
,
Huber
R.
(
1997
).
Crystal structure and mechanism of human l-arginine:glycine amidinotransferase: a mitochondrial enzyme involved in creatine biosynthesis
.
EMBO J.
16
,
3373
-
3385
.
Hunter
S.
,
Apweiler
R.
,
Attwood
T. K.
,
Bairoch
A.
,
Bateman
A.
,
Binns
D.
,
Bork
P.
,
Das
U.
,
Daugherty
L.
,
Duquenne
L.
, et al. 
. (
2009
).
InterPro: the integrative protein signature database
.
Nucleic Acids Res.
37
,
D211
-
D215
.
Itoh-Satoh
M.
,
Hayashi
T.
,
Nishi
H.
,
Koga
Y.
,
Arimura
T.
,
Koyanagi
T.
,
Takahashi
M.
,
Hohda
S.
,
Ueda
K.
,
Nouchi
T.
, et al. 
. (
2002
).
Titin mutations as the molecular basis for dilated cardiomyopathy
.
Biochem. Biophys. Res. Commun.
291
,
385
-
393
.
Kai
M.
,
Wada
I.
,
Imai
S.
,
Sakane
F.
,
Kanoh
H.
(
1997
).
Cloning and characterization of two human isozymes of Mg2+-independent phosphatidic acid phosphatase
.
J. Biol. Chem.
272
,
24572
-
24578
.
Ketterson
E. D.
,
Nolan
V.
 Jr
,
Wolf
L.
,
Ziegenfus
C.
,
Dufty
A. M.
 Jr
,
Ball
G. F.
,
Johnsen
T. S.
(
1991
).
Testosterone and avian life histories: the effect of experimentally elevated testosterone on corticosterone and body mass in dark-eyed juncos
.
Horm. Behav.
25
,
489
-
503
.
Ketterson
E. D.
,
Nolan
V.
,
Wolf
L.
,
Ziegenfus
C.
(
1992
).
Testosterone and avian life histories – effects of experimentally elevated testosterone on behavior and correlates of fitness in the dark-eyed junco (Junco hyemalis)
.
Am. Nat.
140
,
980
-
999
.
Ketterson
E. D.
,
Nolan
V.
,
Cawthorn
M. J.
,
Parker
P. G.
,
Ziegenfus
C.
(
1996
).
Phenotypic engineering: using hormones to explore the mechanistic and functional bases of phenotypic variation in nature
.
Ibis
138
,
70
-
86
.
Ketterson
E. D.
,
Nolan
V.
 Jr
,
Sandell
M.
(
2005
).
Testosterone in females: mediator of adaptive traits, constraint on sexual dimorphism, or both?
Am. Nat.
166
Suppl. 4
,
S85
-
S98
.
Ketterson
E. D.
,
Atwell
J. W.
,
McGlothlin
J. W.
(
2009
).
Phenotypic integration and independence: hormones, performance, and response to environmental change
.
Integr. Comp. Biol.
49
,
365
-
379
.
Kraemer
W. J.
,
Volek
J. S.
(
1999
).
Creatine supplementation. Its role in human performance
.
Clin. Sports Med.
18
,
651
-
666
,
ix
.
Kretzschmar
M.
,
Kaiser
K.
,
Lottspeich
F.
,
Meisterernst
M.
(
1994
).
A novel mediator of class II gene transcription with homology to viral immediate-early transcriptional regulators
.
Cell
78
,
525
-
534
.
Kriskó
I.
,
Walker
J. B.
(
1966
).
Influence of sex hormones on amidinotransferase levels. Metabolic control of creatine biosynthesis
.
Acta Endocrinol. (Copenh.)
53
,
655
-
662
.
Labrie
F.
,
Luu-The
V.
,
Calvo
E.
,
Martel
C.
,
Cloutier
J.
,
Gauthier
S.
,
Belleau
P.
,
Morissette
J.
,
Lévesque
M. H.
,
Labrie
C.
(
2005
).
Tetrahydrogestrinone induces a genomic signature typical of a potent anabolic steroid
.
J. Endocrinol.
184
,
427
-
433
.
Lee
S. J.
,
McPherron
A. C.
(
2001
).
Regulation of myostatin activity and muscle growth
.
Proc. Natl. Acad. Sci. USA
98
,
9306
-
9311
.
Lemercier
C.
,
Brocard
M. P.
,
Puvion-Dutilleul
F.
,
Kao
H. Y.
,
Albagli
O.
,
Khochbin
S.
(
2002
).
Class II histone deacetylases are directly recruited by BCL6 transcriptional repressor
.
J. Biol. Chem.
277
,
22045
-
22052
.
Litman
G. W.
,
Rast
J. P.
,
Shamblott
M. J.
,
Haire
R. N.
,
Hulst
M.
,
Roess
W.
,
Litman
R. T.
,
Hinds-Frey
K. R.
,
Zilch
A.
,
Amemiya
C. T.
(
1993
).
Phylogenetic diversification of immunoglobulin genes and the antibody repertoire
.
Mol. Biol. Evol.
10
,
60
-
72
.
Liu
D.
,
Sartor
M. A.
,
Nader
G. A.
,
Gutmann
L.
,
Treutelaar
M. K.
,
Pistilli
E. E.
,
Iglayreger
H. B.
,
Burant
C. F.
,
Hoffman
E. P.
,
Gordon
P. M.
(
2010
).
Skeletal muscle gene expression in response to resistance exercise: sex specific regulation
.
BMC Genomics
11
,
659
.
Lopez
J.
,
Colbourne
J.
(
2011
).
Dual-labeled expression microarray protocol for high-throughput genomic investigations
.
CGB Technical Report 2011
201
,
2
.
Lynn
S. E.
,
Houtman
A. M.
,
Weathers
W. W.
,
Ketterson
E. D.
,
Nolan
V.
 Jr
(
2000
).
Testosterone increases activity but not daily energy expenditure in captive male dark-eyed juncos, Junco hyemalis
.
Anim. Behav.
60
,
581
-
587
.
Maher
A. C.
,
Fu
M. H.
,
Isfort
R. J.
,
Varbanov
A. R.
,
Qu
X. A.
,
Tarnopolsky
M. A.
(
2009
).
Sex differences in global mRNA content of human skeletal muscle
.
PLoS ONE
4
,
e6335
.
Mank
J. E.
(
2007
).
The evolution of sexually selected traits and antagonistic androgen expression in actinopterygiian fishes
.
Am. Nat.
169
,
142
-
149
.
Marden
J. H.
(
1987
).
Maximum lift production during takeoff in flying animals
.
J. Exp. Biol.
130
,
235
-
258
.
Marler
C. A.
,
Moore
M. C.
(
1988
).
Evolutionary costs of aggression revealed by testosterone manipulations in free-living male lizards
.
Behav. Ecol. Sociobiol.
23
,
21
-
26
.
Marler
C. A.
,
Walsberg
G.
,
White
M. L.
,
Moore
M.
(
1995
).
Increased energy-expenditure due to increased territorial defense in male lizards after phenotypic manipulation
.
Behav. Ecol. Sociobiol.
37
,
225
-
231
.
Michael
H.
,
Härkönen
P. L.
,
Väänänen
H. K.
,
Hentunen
T. A.
(
2005
).
Estrogen and testosterone use different cellular pathways to inhibit osteoclastogenesis and bone resorption
.
J. Bone Miner. Res.
20
,
2224
-
2232
.
Miller
A. H.
(
1941
).
Speciation in the avian species
Junco. Univ. Calif. Publ. Zool.
44
,
173
-
434
.
Miura
Y.
,
Kato
H.
,
Noguchi
T.
(
1992
).
Effect of dietary proteins on insulin-like growth factor-1 (IGF-1) messenger ribonucleic acid content in rat liver
.
Br. J. Nutr.
67
,
257
-
265
.
Mock
B. A.
,
Nacy
C. A.
(
1988
).
Hormonal modulation of sex differences in resistance to Leishmania major systemic infections
.
Infect. Immun.
56
,
3316
-
3319
.
Mokkonen
M.
,
Koskela
E.
,
Mappes
T.
,
Mills
S. C.
(
2012
).
Sexual antagonism for testosterone maintains multiple mating behaviour
.
J. Anim. Ecol.
81
,
277
-
283
.
Møller
A. P.
,
Garamszegi
L. Z.
,
Gil
D.
,
Hurtrez-Bousses
S.
,
Eens
M.
(
2005
).
Correlated evolution of male and female testosterone profiles in birds and its consequences
.
Behav. Ecol. Sociobiol.
58
,
534
-
544
.
Montano
M.
,
Flanagan
J. N.
,
Jiang
L.
,
Sebastiani
P.
,
Rarick
M.
,
LeBrasseur
N. K.
,
Morris
C. A.
,
Jasuja
R.
,
Bhasin
S.
(
2007
).
Transcriptional profiling of testosterone-regulated genes in the skeletal muscle of human immunodeficiency virus-infected men experiencing weight loss
.
J. Clin. Endocrinol. Metab.
92
,
2793
-
2802
.
Nakao
A.
,
Röijer
E.
,
Imamura
T.
,
Souchelnytskyi
S.
,
Stenman
G.
,
Heldin
C. H.
,
ten Dijke
P.
(
1997
).
Identification of Smad2, a human Mad-related protein in the transforming growth factor beta signaling pathway
.
J. Biol. Chem.
272
,
2896
-
2900
.
Nolan
V.
,
Ketterson
E. D.
,
Cristol
D. A.
,
Rogers
C. M.
,
Clotfelter
E. D.
,
Titus
R. C.
,
Schoech
S. J.
,
Snajdr
E.
(
2002
).
Dark-eyed junco: Junco hyemalis
. In
The Birds of North America
(ed.
Poole
A.
), pp.
1
-
44
.
Ithaca, NY
:
Cornell Laboratory of Ornithology
.
Oda
K.
,
Matsuoka
Y.
,
Funahashi
A.
,
Kitano
H.
(
2005
).
A comprehensive pathway map of epidermal growth factor receptor signaling
.
Mol Syst Biol
.
1
,
2005.0010
.
Oppliger
A.
,
Giorgi
M. S.
,
Conelli
A.
,
Nembrini
M.
,
John-Alder
H. B.
(
2004
).
Effect of testosterone on immunocompetence, parasite load, and metabolism in the common wall lizard (Podarcis muralis)
.
Can. J. Zool.
82
,
1713
-
1719
.
Oursler
M. J. J.
,
Osdoby
P.
,
Pyfferoen
J.
,
Riggs
B. L.
,
Spelsberg
T. C.
(
1991
).
Avian osteoclasts as estrogen target cells
.
Proc. Natl. Acad. Sci. USA
88
,
6613
-
6617
.
Park
K. D.
,
Park
J.
,
Ko
J.
,
Kim
B. C.
,
Kim
H. S.
,
Ahn
K.
,
Do
K. T.
,
Choi
H.
,
Kim
H. M.
,
Song
S.
, et al. 
. (
2012
).
Whole transcriptome analyses of six thoroughbred horses before and after exercise using RNA-Seq
.
BMC Genomics
13
,
473
.
Pederson
L.
,
Kremer
M.
,
Judd
J.
,
Pascoe
D.
,
Spelsberg
T. C.
,
Riggs
B. L.
,
Oursler
M. J.
(
1999
).
Androgens regulate bone resorption activity of isolated osteoclasts in vitro
.
Proc. Natl. Acad. Sci. USA
96
,
505
-
510
.
Peterson
M. P.
,
Whittaker
D. J.
,
Ambreth
S.
,
Sureshchandra
S.
,
Buechlein
A.
,
Podicheti
R.
,
Choi
J. H.
,
Lai
Z.
,
Mockatis
K.
,
Colbourne
J.
, et al. 
. (
2012
).
De novo transcriptome sequencing in a songbird, the dark-eyed junco (Junco hyemalis): genomic tools for an ecological model system
.
BMC Genomics
13
,
305
.
Peterson
M. P.
,
Rosvall
K. A.
,
Choi
J. H.
,
Ziegenfus
C.
,
Tang
H.
,
Colbourne
J. K.
,
Ketterson
E. D.
(
2013
).
Testosterone affects neural gene expression differently in male and female juncos: a role for hormones in mediating sexual dimorphism and conflict
.
PLoS ONE
8
,
e61784
.
Petry
C. J.
,
Ong
K. K.
,
Wingate
D. L.
,
Brown
J.
,
Scott
C. D.
,
Jones
E. Y.
,
Pembrey
M. E.
,
Dunger
D. B.
Alspac Study Team
(
2005
).
Genetic variation in the type 2 insulin-like growth factor receptor gene and disparity in childhood height
.
Growth Horm. IGF Res.
15
,
363
-
368
.
Platt
J. L.
,
Nath
K. A.
(
1998
).
Heme oxygenase: protective gene or Trojan horse
.
Nat. Med.
4
,
1364
-
1365
.
Proctor
N. S.
(
1993
).
Manual of Ornithology: Avian Structure and Function
.
North Haven, CT
:
Yale University Press
.
R Development Core Team
(
2010
).
R: A Language and Environment for Statistical Computing
.
R Foundation for Statistical Computing
.
Vienna, Austria
. Available at http://www.r-project.org.
Reed
W. L.
,
Clark
M. E.
,
Parker
P. G.
,
Raouf
S. A.
,
Arguedas
N.
,
Monk
D. S.
,
Snajdr
E.
,
Nolan
V.
 Jr
,
Ketterson
E. D.
(
2006
).
Physiological effects on demography: a long-term experimental study of testosterone's effects on fitness
.
Am. Nat.
167
,
667
-
683
.
Reichard
D.
,
Ketterson
E.
(
2012
).
Estimation of female home-range size during the nestling period of dark-eyed juncos
.
Wilson J. Ornithol.
124
,
614
-
620
.
Rice
W.
(
1984
).
Sex chromosomes and the evolution of sexual dimorphism
.
Evolution
38
,
735
-
742
.
Roberts
M. L.
,
Buchanan
K. L.
,
Evans
M. R.
(
2004
).
Testing the immunocompetence handicap hypothesis: a review of the evidence
.
Anim. Behav.
68
,
227
-
239
.
Roth
S. M.
,
Ferrell
R. E.
,
Peters
D. G.
,
Metter
E. J.
,
Hurley
B. F.
,
Rogers
M. A.
(
2002
).
Influence of age, sex, and strength training on human muscle gene expression determined by microarray
.
Physiol. Genomics
10
,
181
-
190
.
Rowan
W.
(
1925
).
Relation of light to bird migration and developmental changes
.
Nature
115
,
494
-
495
.
Schramm
L.
,
Pendergrast
P. S.
,
Sun
Y.
,
Hernandez
N.
(
2000
).
Different human TFIIIB activities direct RNA polymerase III transcription from TATA-containing and TATA-less promoters
.
Genes Dev.
14
,
2650
-
2663
.
Schroeder
A.
,
Mueller
O.
,
Stocker
S.
,
Salowsky
R.
,
Leiber
M.
,
Gassmann
M.
,
Lightfoot
S.
,
Menzel
W.
,
Granzow
M.
,
Ragg
T.
(
2006
).
The RIN: an RNA integrity number for assigning integrity values to RNA measurements
.
BMC Mol. Biol.
7
,
3
.
Sebbag-Sznajder
N.
,
Raitskin
O.
,
Angenitzki
M.
,
Sato
T. A.
,
Sperling
J.
,
Sperling
R.
(
2012
).
Regulation of alternative splicing within the supraspliceosome
.
J. Struct. Biol.
177
,
152
-
159
.
Sims
R. J.
 III
,
Mandal
S. S.
,
Reinberg
D.
(
2004
).
Recent highlights of RNA-polymerase-II-mediated transcription
.
Curr. Opin. Cell Biol.
16
,
263
-
271
.
Smyth
G. K.
(
2005
).
Limma: linear models for microarray data
.
Bioinformatics and Computational Biology Solution Using R and Bioconductor
2005
,
397
-
420
.
Sterling
D.
,
Reithmeier
R. A. F.
,
Casey
J. R.
(
2001
).
A transport metabolon. Functional interaction of carbonic anhydrase II and chloride/bicarbonate exchangers
.
J. Biol. Chem.
276
,
47886
-
47894
.
Storey
J.
(
2002
).
A direct approach to false discovery rates
.
J. R. Stat. Soc. B
64
,
479
-
498
.
Van Doorn
G. S.
(
2009
).
Intralocus sexual conflict
.
Ann. New York Acad. Sci.
1168
,
52
-
71
.
van Nas
A.
,
Guhathakurta
D.
,
Wang
S. S.
,
Yehya
N.
,
Horvath
S.
,
Zhang
B.
,
Ingram-Drake
L.
,
Chaudhuri
G.
,
Schadt
E. E.
,
Drake
T. A.
, et al. 
. (
2009
).
Elucidating the role of gonadal hormones in sexually dimorphic gene coexpression networks
.
Endocrinology
150
,
1235
-
1249
.
Wakley
G. K.
,
Schutte
H. D.
 Jr
,
Hannon
K. S.
,
Turner
R. T.
(
1991
).
Androgen treatment prevents loss of cancellous bone in the orchidectomized rat
.
J. Bone Miner. Res.
6
,
325
-
330
.
Welle
S.
,
Tawil
R.
,
Thornton
C. A.
(
2008
).
Sex-related differences in gene expression in human skeletal muscle
.
PLoS ONE
3
,
e1385
.
Wikelski
M.
,
Lynn
S.
,
Breuner
C.
,
Wingfield
J. C.
,
Kenagy
G. J.
(
1999
).
Energy metabolism, testosterone and corticosterone in white-crowned sparrows
.
J. Comp. Physiol. A
185
,
463
-
470
.
Wiley
C. J.
,
Goldizen
A. W.
(
2003
).
Testosterone is correlated with courtship but not aggression in the tropical buff-banded rail, Gallirallus philippensis
.
Horm. Behav.
43
,
554
-
560
.
Williamson
D. H.
,
Mellanby
J.
,
Krebs
H. A.
(
1962
).
Enzymic determination of d(–)-β-hydroxybutyric acid and acetoacetic acid in blood
.
Biochem. J.
82
,
90
-
96
.
Woodward
C. J. H.
(
1993
).
A re-evaluation of the anabolic effect of testosterone in rats: interactions with gonadectomy, adrenalectomy and hypophysectomy
.
Acta Endocrinol.
128
,
473
-
477
.
Xu
X.
,
Coats
J. K.
,
Yang
C. F.
,
Wang
A.
,
Ahmed
O. M.
,
Alvarado
M.
,
Izumi
T.
,
Shah
N. M.
(
2012
).
Modular genetic control of sexually dimorphic behaviors
.
Cell
148
,
596
-
607
.
Yang
X.
,
Schadt
E. E.
,
Wang
S.
,
Wang
H.
,
Arnold
A. P.
,
Ingram-Drake
L.
,
Drake
T. A.
,
Lusis
A. J.
(
2006
).
Tissue-specific expression and regulation of sexually dimorphic genes in mice
.
Genome Res.
16
,
995
-
1004
.
Yoder
J. A.
,
Hawke
N. A.
,
Eason
D. D.
,
Mueller
M. G.
,
Davids
B. J.
,
Gillin
F. D.
,
Litman
G. W.
(
2002
).
BIVM, a novel gene widely distributed among deuterostomes, shares a core sequence with an unusual gene in Giardia lamblia
.
Genomics
79
,
750
-
755
.
Zagon
I. S.
,
Donahue
R. N.
,
Rogosnitzky
M.
,
McLaughlin
P. J.
(
2008
).
Imiquimod upregulates the opioid growth factor receptor to inhibit cell proliferation independent of immune function
.
Exp. Biol. Med. (Maywood)
233
,
968
-
979
.
Zheng
X. M.
,
Black
D.
,
Chambon
P.
,
Egly
J. M.
(
1990
).
Sequencing and expression of complementary DNA for the general transcription factor BTF3
.
Nature
344
,
556
-
559
.
Zysling
D. A.
,
Greives
T. J.
,
Breuner
C. W.
,
Casto
J. M.
,
Demas
G. E.
,
Ketterson
E. D.
(
2006
).
Behavioral and physiological responses to experimentally elevated testosterone in female dark-eyed juncos (Junco hyemalis carolinensis)
.
Horm. Behav.
50
,
200
-
207
.

Competing interests

The authors declare no competing financial interests.

Supplementary information