Adaptive evolution of mammalian aromatases: lessons from Suiformes

Conceptus estrogen and prostaglandins provide the maternal recognition of pregnancy signal to prevent luteolysis during early pregnancy in the pig†

Conceptus estrogens and prostaglandins have long been considered the primary signals for maternal recognition of pregnancy (MRP) in the pig. However, loss-of-function studies targeting conceptus aromatase genes (CYP19A1 and CYP19A2) and prostaglandin-endoperoxide synthase 2 (PTGS2) indicated that conceptuses can not only signal MRP without estrogens or prostaglandins but can maintain early pregnancy. However, complete loss of estrogen production leads to abortion after day 25 of gestation. Although neither conceptus estrogens nor prostaglandins had a significant effect on early maintenance of corpora lutea (CL) function alone, the two conceptus factors have a biological relationship. To investigate the role that both conceptus estrogens and prostaglandins have on MRP and maintenance of pregnancy, a triple loss-of function model (TKO) was generated for conceptus CYP19A1, CYP19A2, and PTGS2. In addition, a conceptus CYP19A2-/- model (A2KO) was established to determine the role of placental estrogen during later pregnancy. Estrogen and prostaglandin synthesis were greatly reduced in TKO concept uses which resulted in a failure to inhibit luteolysis after day 15 of pregnancy despite the presence of conceptuses in the uterine lumen. However, A2KO placentae not only maintained functional CL but were able to maintain pregnancy to day 32 of gestation. Despite the loss of placental CYP19A2 expression, the allantois fluid content of estrogen was not affected as the placenta compensated by expressing CYP19A1 and CYP19A3, which are normally absent in controls. Results suggest conceptuses can signal MRP through production of conceptus PGE or stimulating PGE synthesis from the endometrium through conceptus estrogen. Failure of conceptuses to produce both factors results in failure of MRP and loss of pregnancy.

Maternal recognition of pregnancy in the pig: A servomechanism involving sex steroids, cytokines and prostaglandins

Maternal recognition of pregnancy (MRP) is a term utilized in mammals to describe pathways in which the conceptus alters the endometrial environment to prevent regression of corpora lutea to ensure continued production of progesterone (P4) required for establishment and maintenance of pregnancy. For nearly 40 years after publication of the endocrine/exocrine theory, conceptus estrogen (E2) was considered the primary maternal recognition signal in the pig. Conceptus production of prostaglandin E2 (PGE2) was also considered to be a major factor in preventing luteolysis. An addition to E2 and PGE2, pig conceptuses produce interleukin 1B2 (IL1B2) and interferons (IFN) delta (IFND) and gamma (IFNG). The present review provides brief history of the discovery of E2, PGs and IFNS which led to research investigating the role of these conceptus secreted factors in establishing and maintaining pregnancy in the pig. The recent utilization of gene editing technology allowed a more direct approach to investigate the in vivo roles of IL1B2, E2, PGE2, AND IFNG for establishment of pregnancy. These studies revealed unknown functions for IFNG and ILB2 in addition to PGE2 and E2. Thus, pregnancy recognition signal is via a servomechanism in requiring sequential effects of P4, E2, IL1B2, PGE2 and IFNG. Results indicate that the original established dogma for the role of conceptus E2 and PGs in MRP is a far too simplified model that involves the interplay of numerous mechanisms for inhibiting luteolysis, inducing critical elongation of the conceptuses and resolution of inflammation in pigs.

Identification, Expression and Evolutional Analysis of Two cyp19-like Genes in Amphioxus

The mechanism of sex determination and differentiation in animals remains a central focus of reproductive and developmental biology research, and the regulation of sex differentiation in amphioxus remains poorly understood. Cytochrome P450 Family 19 Subfamily A member 1 (CYP19A1) is a crucial sex differentiation gene that catalyzes the conversion of androgens into estrogens. In this study, we identified two aromatase-like genes in amphioxus: cyp19-like1 and cyp19-like2. The cyp19-like1 is more primitive and may represent the ancestral form of cyp19 in zebrafish and other vertebrates, while the cyp19-like2 is likely the result of gene duplication within amphioxus. To gain further insights into the expression level of these two aromatase-like, we examined their expression in different tissues and during different stages of gonad development. While the expression level of the two genes differs in tissues, both are highly expressed in the gonad primordium and are primarily localized to microsomal membrane systems. However, as development proceeds, their expression level decreases significantly. This study enhances our understanding of sex differentiation mechanisms in amphioxus and provides valuable insights into the formation and evolution of sex determination mechanisms in vertebrates.

Charged Amino Acids in the Transmembrane Helix Strongly Affect the Enzyme Activity of Aromatase

Estrogens play critical roles in embryonic development, gonadal sex differentiation, behavior, and reproduction in vertebrates and in several human cancers. Estrogens are synthesized from testosterone and androstenedione by the endoplasmic reticulum membrane-bound P450 aromatase/cytochrome P450 oxidoreductase complex (CYP19/CPR). Here, we report the characterization of novel mammalian CYP19 isoforms encoded by CYP19 gene copies. These CYP19 isoforms are all defined by a combination of mutations in the N-terminal transmembrane helix (E42K, D43N) and in helix C of the catalytic domain (P146T, F147Y). The mutant CYP19 isoforms show increased androgen conversion due to the KN transmembrane helix. In addition, the TY substitutions in helix C result in a substrate preference for androstenedione. Our structural models suggest that CYP19 mutants may interact differently with the membrane (affecting substrate uptake) and with CPR (affecting electron transfer), providing structural clues for the catalytic differences.

Phylogenomic and expression analysis of Colossoma macropomum cyp19a1a and cyp19a1b and their non-classical role in tambaqui sex differentiation

The genes coding for Cytochrome P450 aromatase (cyp19a1a and cyp19a1b) and estrogen (E₂) receptors (esr1, esr2a and esr2b) play a conserved role in ovarian differentiation and development among teleosts. Classically, the "gonad form" of aromatase, coded by the cyp19a1a, is responsible for the ovarian differentiation in genetic females via ligation and activation of the Esr, which mediates the endocrine and exocrine signaling to allow or block the establishment of the feminine phenotype. However, in neotropical species, studies on the molecular and endocrine processes involved in gonad differentiation as well as on the effects of sex modulators are recent and scarce. In this study, we combined in silico analysis, real-time quantitative PCR (qPCR) assay and quantification of E₂ plasma levels of differentiating tambaqui (Colossoma macropomum) to unveil the roles of the paralogs cypa19a1a and cyp19a1b during sex differentiation. Although the synteny of each gene is very conserved among characids, the genomic environment displays striking differences in comparison to model teleost species, with many rearrangements in cyp19a1a and cyp19a1b adjacencies and transposable element traces in both regulatory regions. The high dissimilarity (DI) of SF-1 binding motifs in cyp19a1a (DI = 10.06 to 14.90 %) and cyp19a1b (DI = 8.41 to 13.50 %) regulatory region, respectively, may reflect in an alternative pathway in tambaqui. Indeed, while low transcription of cyp19a1a was detected prior to sex differentiation, the expression of cyp19a1b and esr2a presented a large variation at this phase, which could be associated with sex-specific differential expression. Histological analysis revealed that anti-estradiol treatments did not affect gonadal sex ratios, although Fadrozole (50 mg kg^-1 of food) reduced E₂ plasma levels (p < 0,005) as well cyp19a1a transcription; and tamoxifen (200 mg kg^-1 of food) down regulated both cyp19a1a and cyp19a1b but did not influence E₂ levels. Altogether, our results bring into light new insights about the evolutionary fate of cyp19a1 paralogs in neotropical fish, which may have generated uncommon roles for the gonadal and brain forms of cyp19a1 genes and the unexpected lack of effect of endocrine disruptors on tambaqui sexual differentiation.

Toxic effects of octocrylene on zebrafish larvae and liver cell line (ZFL)

Octocrylene (OC) is a broad-spectrum ultraviolet-absorbing chemical used in sunscreen and other personal care products. Its health effects are a concern because it has been detected in water, fish, humans, and food chains. In vivo and in vitro investigations were performed in zebrafish (Danio rerio) larvae and a zebrafish liver cell line (ZFL), respectively, to understand the potential risks and molecular mechanisms of OC toxicity. The 96-h median lethal concentration (LC₅₀) of OC was determined to be 251.8 μM in larvae and 5.5 μM in ZFL cells. Quantitative real-time PCR (qRT-PCR) showed that OC induced the expression of genes for CYPs (CYP1A, CYP3A65), estrogen receptors (ERα, ERβ1, GPER), vitellogenin (VTG1), and sex determination (BRCA2, CYP19A, DMRT1, SOX9A), both in vitro and in vivo. A whole-transcriptome sequencing method was used to evaluate the gene expression profile of larvae exposed to OC. OC was found to mediate the biosynthesis of estrogens (such as estriol) and affect the antioxidant pathway (glutathione transferases and peroxisome). These findings clarify the toxic effects and molecular mechanisms of OC and support banning its use in cosmetics.

Genomic Structure of the Porcine CYP19 Locus and Expression of the CYP19A3 Paralog

Proper, tissue-specific regulation of CYP19, the gene encoding aromatase, the key enzyme of estrogen synthesis, is essential for reproductive processes. Here, we analyzed transcriptional regulation of the porcine CYP19 in female and male gonads and brain by 5'RACE and RT-PCR and comprehensively mapped the pig CYP19 locus by in silico analysis. Our data revealed that the complete locus, including three paralogous copies, CYP19A1, CYP19A2 and CYP19A3, spans approximately 330 kb of the porcine chromosome 1. The locus also harbors the first exon of the Gliomedin gene (GLDN) in reverse orientation. Only transcripts of the CYP19A3 paralog were substantially expressed in gonads and hypothalamus. We identified CYP19A3-associated untranslated exons approximately 160 kb and 50 kb distal from the first codon. The 5´ untranslated regions of transcripts were derived from either a proximal or from one of these distal untranslated exons. Transcripts including only untranslated exons could be amplified from testis, thus suggesting long non-coding transcripts. The data revealed an additional layer of complexity in the regulation of the porcine CYP19 locus. Tissue-specific expression is not only achieved by tissue- and stage-specific expression of the three different CYP19 paralogs, but also by directing the expression of CYP19A3 from different, proximal and distal promoter regions.

Molecular and Structural Evolution of Cytochrome P450 Aromatase

Aromatase is the cytochrome P450 enzyme converting androgens into estrogen in the last phase of steroidogenesis. As estrogens are crucial in reproductive biology, aromatase is found in vertebrates and the invertebrates of the genus Branchiostoma, where it carries out the aromatization reaction of the A-ring of androgens that produces estrogens. Here, we investigate the molecular evolution of this unique and highly substrate-selective enzyme by means of structural, sequence alignment, and homology modeling, shedding light on its key role in species conservation. The alignments led to the identification of a core structure that, together with key and unique amino acids located in the active site and the substrate recognition sites, has been well conserved during evolution. Structural analysis shows what their roles are and the reason why they have been preserved. Moreover, the residues involved in the interaction with the redox partner and some phosphorylation sites appeared late during evolution. These data reveal how highly substrate-selective cytochrome P450 has evolved, indicating that the driving forces for evolution have been the optimization of the interaction with the redox partner and the introduction of phosphorylation sites that give the possibility of modulating its activity in a rapid way.

Toxic effects and transcriptome analyses of zebrafish (Danio rerio) larvae exposed to benzophenones

Sunscreen chemicals, such as benzophenones (BPs), are common environmental contaminants that are posing a growing health concern due to their increasing presence in water, fish, and human systems. Benzoresorcinol (BP1), oxybenzone (BP3), and dioxybenzone (BP8) are the most commonly used BPs for their ability to protect from sunburn by absorbing a broad spectrum of ultraviolet radiation. In this study, zebrafish larvae were used as an in vivo model to investigate the potential risks and molecular mechanisms of the toxic effects of BPs. The effects of these BPs on the gene expression in the aryl hydrocarbon receptor pathway, estrogen receptor pathway, and sex differentiation were detected using quantitative real-time PCR. All BPs were found to function as agonists of the estrogen receptors α and β1, indicating that these BPs likely undergo similar molecular metabolism in vivo, whereby they can activate cytochrome P450 genes and promote the expression of CYP19A and DMRT1. Furthermore, the gene expression profile of larvae after BP3 exposure was evaluated using a whole transcriptome sequencing approach. BP3 affected estradiol biosynthesis and sex differentiation. It also regulated gonadotropin-releasing hormone, thus interfering with the endocrine system. As a xenobiotic toxicant, BP3 upregulated the expression of cytochrome P450 genes (CYP1A and CYP3A65) and glutathione metabolism-related genes (GSTA, GSTM, and GSTP). It also interfered with the nervous system by regulating the calcium signaling pathway. These findings will be useful for understanding the toxicity mechanisms and metabolism of BPs in aquatic organisms and promote the regulation of these chemicals in the environment.

Sequence comparisons of cytochrome P450 aromatases from Australian animals predict differences in enzymatic activity and/or efficiency†

Aromatase (P450arom, CYP19A1) is the terminal enzyme in the synthesis of the steroid hormone family of estrogens. Not surprisingly, this enzyme has structural similarities between the limited number of species studied thus far. This study examined the structure of aromatases from four diverse Australian species including a marsupial (tammar wallaby; Macropus eugenii), monotreme (platypus; Ornithorhynchus anatinus), ratite (emu; Dromaius novaehollandiae) and lizard (bearded dragon; Pogona vitticeps). We successfully built homology models for each species, using the only crystallographically determined structure available, human aromatase. The amino acid sequences showed high amino acid sequence identity to the human aromatase: wallaby 81%, platypus 73%, emu 75% and bearded dragon at 74%. The overall structure was highly conserved among the five species, although there were non-secondary structures (loops and bends) that were variable and flexible that may result in some differences in catalytic activity. At the N-terminal regions, there were deletions and variations that suggest that functional distinctions may be found. We found that the active sites of all these proteins were identical, except for a slight variation in the emu. The electrostatic potential across the surfaces of these aromatases highlighted likely variations to the protein-protein interactions of these enzymes with both redox partner cytochrome P450 reductase and possibly homodimerization in the case of the platypus, which has been postulated for the human aromatase enzyme. Given the high natural selection pressures on reproductive strategies, the relatively high degree of conservation of aromatase sequence and structure across species suggests that there is biochemically very little scope for changes to have evolved without the loss of enzyme activity.

Basal teleosts provide new insights into the evolutionary history of teleost-duplicated aromatase

Duplicated cyp19a1 genes (cyp19a1a encoding aromatase a and cyp19a1b encoding aromatase b) have been identified in an increasing number of teleost species. Cyp19a1a is mainly expressed in the gonads, while cyp19a1b is mainly expressed in the brain, specifically in radial glial cells, as largely investigated by Kah and collaborators. The third round of whole-genome duplication that specifically occurred in the teleost lineage (TWGD or 3R) is likely at the origin of the duplicated cyp19a1 paralogs. In contrast to the situation in other teleosts, our previous studies identified a single cyp19a1 in eels (Anguilla), which are representative species of a basal group of teleosts, Elopomorpha. In the present study, using genome data mining and phylogenetic and synteny analyses, we confirmed that the whole aromatase genomic region was duplicated in eels, with most aromatase-neighboring genes being conserved in duplicate in eels, as in other teleosts. These findings suggest that specific gene loss of one of the 3R-duplicated cyp19a1 paralogs occurred in Elopomorpha after TWGD. Similarly, a single cyp19a1 gene was found in the arowana, which is a representative species of another basal group of teleosts, Osteoglossomorpha. In eels, the single cyp19a1 is expressed in both the brain and the gonads, as observed for the single CYP19A1 gene present in other vertebrates. The results of phylogenetic, synteny, closest neighboring gene, and promoter structure analyses showed that the single cyp19a1 of the basal teleosts shared conserved properties with both teleost cyp19a1a and cyp19a1b paralogs, which did not allow us to conclude which of the 3R-duplicated paralogs (cyp19a1a or cyp19a1b) was lost in Elopomorpha. Elopomorpha and Osteoglossomorpha cyp19a1 genes exhibited preserved ancestral functions, including expression in both the gonad and brain. We propose that the subfunctionalization of the 3R-duplicated cyp19a1 paralogs expressed specifically in the gonad or brain occurred in Clupeocephala, after the split of Clupeocephala from Elopomorpha and Osteoglossomorpha, which represented a driving force for the conservation of both 3R-duplicated paralogs in all extant Clupeocephala. In contrast, the functional redundancy of the undifferentiated 3R-duplicated cyp19a1 paralogs in elopomorphs and osteoglossomorphs would have favored the loss of one 3R paralog in basal teleosts.

Immunolocalization of steroidogenic enzymes in the vaginal mucous of Galea spixii during the estrous cycle

BACKGROUND

The synthesis of sex steroids is controlled by several enzymes such as17α-hydroxylase cytochrome P450 (P450c17) catalyzing androgen synthesis and aromatase cytochrome P450 (P450arom) catalyzing estrogen synthesis, both of which must complex with the redox partner NADPH-cytochrome P450 oxidoreductase (CPR) for activity. Previous studies have identified expression of steroidogenic enzymes in vaginal tissue, suggesting local sex steroid synthesis. The current studies investigate P450c17, P450aromatase and CPR expression in vaginal mucosa of Galea spixii (Spix cavy) by immuno-histochemical and western immunoblot analyses.

METHODS

Stages of estrous cyclicity were monitored by vaginal exfoliative cytology. After euthanasia, vaginal tissues were retrieved, fixed and frozen at diestrus, proestrus, estrus and metestrus. The ovaries and testis were used as positive control tissues for immunohistochemistry.

RESULTS

Data from cytological study allowed identification of different estrous cycle phases. Immunohistochemical analysis showed different sites of expression of steroidogenic enzymes along with tissue response throughout different phases of the estrous cycle. However, further studies are needed to characterize the derived hormones synthesized by, and the enzymes activities associated with, vaginal tissues.

CONCLUSION

Current results not only support the expression of enzymes involved in sex steroid synthesis in the wall of the vagina, they also indicate that expression changes with the stage of the cycle, both the levels and types of cells exhibiting expression. Thus, changes in proliferation of vaginal epithelial cells and the differentiation of the mucosa may be influenced by local steroid synthesis as well as circulating androgens and estrogens.

Evolutionary comparisons predict that dimerization of human cytochrome P450 aromatase increases its enzymatic activity and efficiency

Estrogen is an essential vertebrate hormone synthesized from androgens involving multiple hydroxylations, catalyzed by cytochrome P450 aromatase (P450arom or CYP19) enzymes. Despite their importance, very few comparative studies have been conducted on vertebrate and/or mammalian P450arom enzymes, either structurally or functionally. Here we directly compared the human (h-) and porcine gonadal (p(g)-) P450arom, as p(g)-P450arom has very low catalytic efficiency, with a ten-fold higher affinity (Km) for a substrate (androstenedione) and ten-fold reduction in turnover (Vmax). We recombinantly expressed these proteins and compared their interactions on a membrane using a quartz crystal microbalance (QCM) and also with the electron donor protein cytochrome P450 oxidoreductase (CPR). Changes in frequency and dissipation in the QCM supported the h-P450arom forming a homodimer that agreed with the FRET data, but not p(g)-P450arom. Analysis of the X-ray crystal structure of the h-P450arom suggested a likely site of homo-dimerization and found that certain key interacting residues were not conserved in pg-P450arom. Molecular dynamics simulations provide support for the importance of these residues in homo-dimerization. Here we propose that the lower affinity and higher activity with reduced release of intermediate metabolites by the h-P450arom is as a consequence of its ability to form homodimers. The functional implications of dimerization provide an important mechanistic step in the requirement for efficient aromatization.

Pathways and genes involved in steroid hormone metabolism in male pigs: a review and update

This paper reviews state-of-the-art knowledge on steroid biosynthesis pathways in the pig and provides an updated characterization of the porcine genes involved in these pathways with particular focus on androgens, estrogens, and 16-androstenes. At least 21 different enzymes appear to be involved in these pathways in porcine tissues together with at least five cofactors. Until now, data on several porcine genes were scarce or confusing. We characterized the complete genomic and transcript sequences of the single porcine CYP11B gene. We analyzed the porcine AKR1 gene cluster and identified four AKR1C, one AKR1C like genes and one AKR1E2 gene. We provide evidence that porcine AKR1C genes are not orthologous to human AKR1C. A new nomenclature is thus needed for this gene family in the pig. Thirty-two genes are now described: transcript (30+2 characterized in this study) and genomic (complete: 18+1 and partial: 12+1) sequences are identified. However, despite increasing knowledge on steroid metabolism in the pig, there is still no explanation of why porcine testes can produce androstenone and epiandrosterone, but not dihydrotestosterone (DHT), which is also a reduced steroid.

Increased testicular Sertoli cell population induced by an estrogen receptor antagonist

Sertoli cell proliferation is prolonged in neonatal boars treated with the aromatase inhibitor letrozole, but porcine testicular aromatase synthesizes a potent, non-aromatizable androgen, 1-hydroxytestosterone, as well as estradiol. Therefore, experiments were conducted to determine whether the Sertoli cell proliferative response to letrozole is due to a loss of estrogen or a loss of androgen signaling. Littermate boars were treated with letrozole, the estrogen receptor blocker ICI 182,780, or vehicle, from 1 week of age and testes collected at 6.5 weeks. Sertoli cell number was increased 30% by letrozole or ICI 182,780 compared with vehicle. Neither treatment affected testosterone, gonadotropins or prolactin. We conclude that Sertoli cell proliferation in neonatal boars is restricted by the local activation of estrogen receptors. The response to letrozole is apparently not mediated by the novel capacity of the porcine gonadal aromatase for 1-hydroxytestosterone but by estradiol synthesis; therefore, aromatase inhibition may have similar effects on Sertoli cell proliferation in other species.

Evolutionary origins of the estrogen signaling system: insights from amphioxus

Classically, the estrogen signaling system has two core components: cytochrome P450 aromatase (CYP19), the enzyme complex that catalyzes the rate limiting step in estrogen biosynthesis; and estrogen receptors (ERs), ligand activated transcription factors that interact with the regulatory region of target genes to mediate the biological effects of estrogen. While the importance of estrogens for regulation of reproduction, development and physiology has been well-documented in gnathostome vertebrates, the evolutionary origins of estrogen as a hormone are still unclear. As invertebrates within the phylum Chordata, cephalochordates (e.g., the amphioxus of the genus Branchiostoma) are among the closest invertebrate relatives of the vertebrates and can provide critical insight into the evolution of vertebrate-specific molecules and pathways. To address this question, this paper briefly reviews relevant earlier studies that help to illuminate the history of the aromatase and ER genes, with a particular emphasis on insights from amphioxus and other invertebrates. We then present new analyses of amphioxus aromatase and ER sequence and function, including an in silico model of the amphioxus aromatase protein, and CYP19 gene analysis. CYP19 shares a conserved gene structure with vertebrates (9 coding exons) and moderate sequence conservation (40% amino acid identity with human CYP19). Modeling of the amphioxus aromatase substrate binding site and simulated docking of androstenedione in comparison to the human aromatase shows that the substrate binding site is conserved and predicts that androstenedione could be a substrate for amphioxus CYP19. The amphioxus ER is structurally similar to vertebrate ERs, but differs in sequence and key residues of the ligand binding domain. Consistent with results from other laboratories, amphioxus ER did not bind radiolabeled estradiol, nor did it modulate gene expression on an estrogen-responsive element (ERE) in the presence of estradiol, 4-hydroxytamoxifen, diethylstilbestrol, bisphenol A or genistein. Interestingly, it has been shown that a related gene, the amphioxus "steroid receptor" (SR), can be activated by estrogens and that amphioxus ER can repress this activation. CYP19, ER and SR are all primarily expressed in gonadal tissue, suggesting an ancient paracrine/autocrine signaling role, but it is not yet known how their expression is regulated and, if estrogen is actually synthesized in amphioxus, whether it has a role in mediating any biological effects. Functional studies are clearly needed to link emerging bioinformatics and in vitro molecular biology results with organismal physiology to develop an understanding of the evolution of estrogen signaling. This article is part of a Special Issue entitled 'Marine organisms'.

Correlated evolution of androgen receptor and aromatase revisited

Conserved interactions among proteins or other molecules can provide strong evidence for coevolution across their evolutionary history. Diverse phylogenetic methods have been applied to identify potential coevolutionary relationships. In most cases, these methods minimally require comparisons of orthologous sequences and appropriate controls to separate effects of selection from the overall evolutionary relationships. In vertebrates, androgen receptor (AR) and cytochrome p450 aromatase (CYP19) share an affinity for androgenic steroids, which serve as receptor ligands and enzyme substrates. In a recent study, Tiwary and Li (Tiwary BK, Li W-H. 2009. Parallel evolution between aromatase and androgen receptor in the animal kingdom. Mol Biol Evol. 26:123-129) reported that AR and CYP19 displayed a signature of ancient and conserved interactions throughout all the Eumetazoa (i.e., cnidarians, protostomes, and deuterostomes). Because these findings conflicted with a number of previous studies, we reanalyzed the data set used by Tiwary and Li. First, our analyses demonstrate that the invertebrate genes used in the previous analysis are not orthologous sequences but instead represent a diverse set of nuclear receptors and CYP enzymes with no confirmed or hypothesized relationships with androgens. Second, we show that 1) their analytical approach, which measures correlations in evolutionary distances between proteins, potentially led to spurious significant relationships due simply to conserved domains and 2) control comparisons provide positive evidence for a strong influence of evolutionary history. We discuss how corrections to this method and analysis of key taxa (e.g., duplications in the teleost fish and suiform lineages) can inform investigations of the coevolutionary relationships between AR and aromatase.