Detection of Pathways Affected by Positive Selection in Primate Lineages Ancestral to Humans.

Gene set enrichment approaches have been increasingly successful in finding signals of recent polygenic selection in the human genome. In this study, we aim at detecting biological pathways affected by positive selection in more ancient human evolutionary history. Focusing on four branches of the primate tree that lead to modern humans, we tested all available protein coding gene trees of the Primates clade for signals of adaptation in these branches, using the likelihood-based branch site test of positive selection. The results of these locus-specific tests were then used as input for a gene set enrichment test, where whole pathways are globally scored for a signal of positive selection, instead of focusing only on outlier "significant" genes. We identified signals of positive selection in several pathways that are mainly involved in immune response, sensory perception, metabolism, and energy production. These pathway-level results are highly significant, even though there is no functional enrichment when only focusing on top scoring genes. Interestingly, several gene sets are found significant at multiple levels in the phylogeny, but different genes are responsible for the selection signal in the different branches. This suggests that the same function has been optimized in different ways at different times in primate evolution.

How many nuclear hormone receptors are there in the human genome?

The sequence of the human genome now allows the definition of the complete set of genes for specific protein families in humans. Because of their involvement in many physiological and pathological processes, the nuclear hormone receptors are a superfamily of crucial medical significance. Although 48 human nuclear receptor genes were identified previously, their total number is unclear from early human genome reports. Here, we report the identification and classification of all nuclear receptor genes in the human genome, and we discuss corresponding transcriptome and proteome diversity.

A conserved retinoid X receptor (RXR) from the mollusk Biomphalaria glabrata transactivates transcription in the presence of retinoids.

Retinoid X receptors (RXR) are members of the nuclear receptor superfamily of ligand-activated transcription factors that have been characterized in a wide variety of metazoan phyla. They act as heterodimer partners of other nuclear receptors, and in vertebrates also activate transcription as homodimers in the presence of a ligand, 9-cis retinoic acid. In order to test the hypothesis that retinoic acid signaling pathways involving RXRs are present in the Lophotrochozoa, we have sought to isolate conserved members of this family from the platyhelminth parasite Schistosoma mansoni and its intermediate host, the mollusk Biomphalaria glabrata. Here we report that an RXR ortholog from B. glabrata (BgRXR) is better conserved, compared with mouse RXRalpha, both in the DNA-binding domain (89% identity) and in the ligand-binding domain (LBD) (81% identity), than are arthropod homologs. In EMSA, BgRXR binds to the direct repeat response element DR1 as a homodimer or as a heterodimer with mammalian RARalpha, LXR, FXR or PPARalpha. When transfected alone into mammalian cell lines, BgRXR transactivated transcription of a reporter gene from the Apo-A1 promoter in the presence of 9-cis retinoic acid or DHA. Constructs with the Gal4 DNA binding domain fused to the hinge and LBDs of BgRXR were used to show that ligand-dependent activation of transcription by BgRXR required its intact AF-2 activation domain, and that the LBD can form homodimers. Finally, the binding of 9-cis retinoic acid preferentially protected the LBD of BgRXR from degradation by trypsin in a proteolysis protection assay. Our results show that BgRXR binds and is activated by retinoids and suggest that retinoid signaling pathways are conserved in the Lophotrochozoa. The nucleotide sequence reported in this paper has been submitted to the GenBank/EBI Data Bank with accession no. AY048663.

Characterization of oestrogen receptors in zebrafish (Danio rerio).

We cloned the cDNAs corresponding to three oestrogen receptors (ERs) in zebrafish (Danio rerio). Sequence analysis and phylogenetic studies demonstrated that two of these genes, ER beta.1 and ER beta.2, arose from duplication of the original ER beta in many species of the fish phylum, whereas ER alpha is unique. Zebrafish ERs behaved as oestrogen-dependent transcription factors in transactivation assays. However, their reactivity to various oestrogen modulators was different compared with that of mouse ERs. ER mRNA expression during zebrafish development is restricted to distinct time periods, as observed by RNase protection assays. ER beta.2 is initially expressed as maternally transmitted RNA, until 6 h after fertilization, when expression disappears. Between 6 and 48 h after fertilization, no ER expression could be observed. After 48 h after fertilization, all ERs, but predominantly ER alpha, began to be expressed. We conclude that oestrogen signal transduction can operate during zebrafish development only within discrete time windows.

Evidence for two distinct functional glucocorticoid receptors in teleost fish.

Using RT-PCR with degenerated primers followed by screening of a rainbow trout (Oncorhynchus mykiss) intestinal cDNA library, we have isolated from the rainbow trout a new corticosteroid receptor which shows high sequence homology with other glucocorticoid receptors (GRs), but is clearly different from the previous trout GR (named rtGR1). Phylogenetic analysis of these two sequences and other GRs known in mammals, amphibians and fishes indicate that the GR duplication is probably common to most teleost fish. The open reading frame of this new trout GR (named rtGR2) encodes a protein of 669 amino acids and in vitro translation produces a protein of 80 kDa that appears clearly different from rtGR1 protein (88 kDa). Using rtGR2 cDNA as a probe, a 7.3 kb transcript was observed in various tIssues suggesting that this gene would lead to expression of a steroid receptor. In vitro studies were used to further characterize this new corticosteroid receptor. Binding studies with recombinant rtGR1 and rtGR2 proteins show that the two receptors have a similar affinity for dexamethasone (GR1 K(d)=5.05+/-0.45 nM; GR2 K(d)=3.04+/-0.79 nM). Co-transfection of an rtGR1 or rtGR2 expression vector into CHO-K1 or COS-7 cells, along with a reporter plasmid containing multiple consensus glucocorticoid response elements, shows that both clones are able to induce transcriptional activity in the presence of cortisol and dexamethasone. Moreover, at 10(-)(6 )M 11-deoxycortisol and corticosterone partially induced rtGR2 transactivation activity but were without effect on rtGR1. The other major teleost reproductive hormones, as well as a number of their precursors or breakdown products of these and corticosteroid hormones, were without major effects on either receptor. Interestingly, rtGR2 transactivational activity was induced at far lower concentrations of dexamethasone or cortisol (cortisol EC(50)=0.72+/-0.87 nM) compared with rtGR1 (cortisol EC(50)=46+/-12 nM). Similarly, even though RU486 inhibited transactivation activity in both rtGR1 and rtGR2, rtGR1 was more sensitive to this GR antagonist. Altogether, these results indicate that these two GR sequences encode for two functionally distinct GRs acting as ligand-inducible transcription factors in rainbow trout.

One billion years of bZIP transcription factor evolution: conservation and change in dimerization and DNA-binding site specificity.

The genomic era has revealed that the large repertoire of observed animal phenotypes is dependent on changes in the expression patterns of a finite number of genes, which are mediated by a plethora of transcription factors (TFs) with distinct specificities. The dimerization of TFs can also increase the complexity of a genetic regulatory network manifold, by combining a small number of monomers into dimers with distinct functions. Therefore, studying the evolution of these dimerizing TFs is vital for understanding how complexity increased during animal evolution. We focus on the second largest family of dimerizing TFs, the basic-region leucine zipper (bZIP), and infer when it expanded and how bZIP DNA-binding and dimerization functions evolved during the major phases of animal evolution. Specifically, we classify the metazoan bZIPs into 19 families and confirm the ancient nature of at least 13 of these families, predating the split of the cnidaria. We observe fixation of a core dimerization network in the last common ancestor of protostomes-deuterostomes. This was followed by an expansion of the number of proteins in the network, but no major dimerization changes in interaction partners, during the emergence of vertebrates. In conclusion, the bZIPs are an excellent model with which to understand how DNA binding and protein interactions of TFs evolved during animal evolution.

RRTree: relative-rate tests between groups of sequences on a phylogenetic tree.

UNLABELLED: RRTree is a user-friendly program for comparing substitution rates between lineages of protein or DNA sequences, relative to an outgroup, through relative rate tests. Genetic diversity is taken into account through use of several sequences, and phylogenetic relations are integrated by topological weighting. AVAILABILITY: The ANSI C source code of RRTree, and compiled versions for Macintosh, MS-DOS/Windows, SUN Solaris, and CGI, are freely available at http://pbil.univ-lyon1.fr/software/rrtree.html CONTACT: marc.robinson@ens-lyon.fr

Fast evolution of interleukin-2 in mammals and positive selection in ruminants.

Interleukin-2 (IL-2) is a cytokine involved in induction and regulation of the immune response in mammals. There have been numerous reports about the search for IL-2 in species other than mammals, and recently an IL-2-like gene has been isolated in chicken. Using PCR, we searched for IL-2 gene sequences in a wide variety of mammals, including marsupials and monotremes, as well as in birds. Although we can readily amplify IL-2 gene fragments in placental mammals, no amplification was obtained in other species. This is best explained by very high substitution rates. This suggest that strategies to isolate IL-2 homologous genes outside mammals should involve functional assays, as for the chicken gene, and not hybridization-based techniques. Nonsynonymous substitution rates are especially high in ruminants, due to positive selection acting on regions important in term of structure-function. We suggest that, although globally similar, the immune response of various mammals is not identical, mainly at the level of cytokine-mediated regulations.

Euteleost fish genomes are characterized by expansion of gene families.

The presence of additional hox clusters in the zebrafish has led to the hypothesis that there was a whole genome duplication at the origin of modern fish. To investigate the generality of this assumption, we analyzed all available actinopterygian fish gene families, and sequenced nuclear receptors from diverse teleost fish. The origin and timing of duplications was systematically determined by phylogenetic analysis. More genes are indeed found in zebrafish than in mouse. This abundance is shared by all major groups of euteleost fish, but not by eels. Phylogenetic analysis shows that it may result from frequent independent duplications, rather than from an ancestral genome duplication. We predict two zebrafish paralogs for most mouse or human genes, thus expressing a note of caution in functional comparison of fish and mammalian genomes. Redundancy appears to be the rule in fish developmental genetics. Finally, our results imply that the outcome of genome projects cannot be extrapolated easily between fish species.