Resolving the ortholog conjecture: orthologs tend to be weakly, but significantly, more similar in function than paralogs.

The function of most proteins is not determined experimentally, but is extrapolated from homologs. According to the "ortholog conjecture", or standard model of phylogenomics, protein function changes rapidly after duplication, leading to paralogs with different functions, while orthologs retain the ancestral function. We report here that a comparison of experimentally supported functional annotations among homologs from 13 genomes mostly supports this model. We show that to analyze GO annotation effectively, several confounding factors need to be controlled: authorship bias, variation of GO term frequency among species, variation of background similarity among species pairs, and propagated annotation bias. After controlling for these biases, we observe that orthologs have generally more similar functional annotations than paralogs. This is especially strong for sub-cellular localization. We observe only a weak decrease in functional similarity with increasing sequence divergence. These findings hold over a large diversity of species; notably orthologs from model organisms such as E. coli, yeast or mouse have conserved function with human proteins.

A protein interaction atlas for the nuclear receptors: properties and quality of a hub-based dimerisation network.

BACKGROUND: The nuclear receptors are a large family of eukaryotic transcription factors that constitute major pharmacological targets. They exert their combinatorial control through homotypic heterodimerisation. Elucidation of this dimerisation network is vital in order to understand the complex dynamics and potential cross-talk involved. RESULTS: Phylogeny, protein-protein interactions, protein-DNA interactions and gene expression data have been integrated to provide a comprehensive and up-to-date description of the topology and properties of the nuclear receptor interaction network in humans. We discriminate between DNA-binding and non-DNA-binding dimers, and provide a comprehensive interaction map, that identifies potential cross-talk between the various pathways of nuclear receptors. CONCLUSION: We infer that the topology of this network is hub-based, and much more connected than previously thought. The hub-based topology of the network and the wide tissue expression pattern of NRs create a highly competitive environment for the common heterodimerising partners. Furthermore, a significant number of negative feedback loops is present, with the hub protein SHP [NR0B2] playing a major role. We also compare the evolution, topology and properties of the nuclear receptor network with the hub-based dimerisation network of the bHLH transcription factors in order to identify both unique themes and ubiquitous properties in gene regulation. In terms of methodology, we conclude that such a comprehensive picture can only be assembled by semi-automated text-mining, manual curation and integration of data from various sources.

Meta-analysis of estrogen response in MCF-7 distinguishes early target genes involved in signaling and cell proliferation from later target genes involved in cell cycle and DNA repair.

ABSTRACT: BACKGROUND: Many studies have been published outlining the global effects of 17 beta-estradiol (E2) on gene expression in human epithelial breast cancer derived MCF-7 cells. These studies show large variation in results, reporting between ~100 and ~1500 genes regulated by E2, with poor overlap. RESULTS: We performed a meta-analysis of these expression studies, using the Rank product method to obtain a more accurate and stable list of the differentially expressed genes, and of pathways regulated by E2. We analyzed 9 time-series data sets, concentrating on response at 3-4 hrs (early) and at 24 hrs (late). We found >1000 statistically significant probe sets after correction for multiple testing at 3-4 hrs, and >2000 significant probe sets at 24 hrs. Differentially expressed genes were examined by pathway analysis. This revealed 15 early response pathways, mostly related to cell signaling and proliferation, and 20 late response pathways, mostly related to breast cancer, cell division, DNA repair and recombination. CONCLUSIONS: Our results show that meta-analysis identified more differentially expressed genes than the individual studies, and that these genes act together in networks. These results provide new insight into E2 regulated mechanisms, especially in the context of breast cancer.

Gene loss and evolutionary rates following whole-genome duplication in teleost fishes.

Teleost fishes provide the first unambiguous support for ancient whole-genome duplication in an animal lineage. Studies in yeast or plants have shown that the effects of such duplications can be mediated by a complex pattern of gene retention and changes in evolutionary pressure. To explore such patterns in fishes, we have determined by phylogenetic analysis the evolutionary origin of 675 Tetraodon duplicated genes assigned to chromosomes, using additional data from other species of actinopterygian fishes. The subset of genes, which was retained in double after the genome duplication, is enriched in development, signaling, behavior, and regulation functional categories. The evolutionary rate of duplicate fish genes appears to be determined by 3 forces: 1) fish proteins evolve faster than mammalian orthologs; 2) the genes kept in double after genome duplication represent the subset under strongest purifying selection; and 3) following duplication, there is an asymmetric acceleration of evolutionary rate in one of the paralogs. These results show that similar mechanisms are at work in fishes as in yeast or plants and provide a framework for future investigation of the consequences of duplication in fishes and other animals.

Serum androgen levels in elite female athletes.

OBJECTIVE: Prior to the implementation of the blood steroidal module of the Athlete Biological Passport, we measured the serum androgen levels among a large population of high-level female athletes as well as the prevalence of biochemical hyperandrogenism and some disorders of sex development (DSD). METHODS AND RESULTS: In 849 elite female athletes, serum T, dehydroepiandrosterone sulphate, androstenedione, SHBG, and gonadotrophins were measured by liquid chromatography-mass spectrometry high resolution or immunoassay. Free T was calculated. The sampling hour, age, and type of athletic event only had a small influence on T concentration, whereas ethnicity had not. Among the 85.5% that did not use oral contraceptives, 168 of 717 athletes were oligo- or amenorrhoic. The oral contraceptive users showed the lowest serum androgen and gonadotrophin and the highest SHBG concentrations. After having removed five doped athletes and five DSD women from our population, median T and free T values were close to those reported in sedentary young women. The 99th percentile for T concentration was calculated at 3.08 nmol/L, which is below the 10 nmol/L threshold used for competition eligibility of hyperandrogenic women with normal androgen sensitivity. Prevalence of hyperandrogenic 46 XY DSD in our athletic population is approximately 7 per 1000, which is 140 times higher than expected in the general population. CONCLUSION: This is the first study to establish normative serum androgens values in elite female athletes, while taking into account the possible influence of menstrual status, oral contraceptive use, type of athletic event, and ethnicity. These findings should help to develop the blood steroidal module of the Athlete Biological Passport and to refine more evidence-based fair policies and recommendations concerning hyperandrogenism in female athletes.

Are 10 min of seating enough to guarantee stable haemoglobin and haematocrit readings for the athlete's biological passport?

Haemoglobin (Hb) and haematocrit (Hct) are measured as indirect markers of doping in athletes. We studied the effect of posture on these parameters in a typical antidoping setting. Venous blood samples were obtained from nine endurance athletes (six males, three females) and nine control subjects (six males, three females) immediately and after 5, 10, 15, 20 and 30 min after having adopted a seated position from normal daily activity. Hb (CV 0.72%) and Hct (CV 0.87%) were determined using an automated cell counter, plasma volume changes were calculated. Differences between the time points, gender and groups were calculated using a mixed-model procedure. Significant changes were observed in the first 10 min after sitting down but no further changes were noted between 10 and 30 min. Mean directional change for Hb and Hct between 0 min and the average of the period from 10 to 30 min was -2.4% (-0.35 g/dl) for Hb and -2.7% (-1.2%) for Hct. Plasma volume increased accordingly. Neither group nor gender had significant effects. Under typical conditions encountered during blood testing in doping control, a period of 10 min in a seated position is sufficient for the vascular volumes to re-equilibrate and to adapt to the new posture.

Comparative modular analysis of gene expression in vertebrate organs.

ABSTRACT: BACKGROUND: The degree of conservation of gene expression between homologous organs largely remains an open question. Several recent studies reported some evidence in favor of such conservation. Most studies compute organs' similarity across all orthologous genes, whereas the expression level of many genes are not informative about organ specificity. RESULTS: Here, we use a modularization algorithm to overcome this limitation through the identification of inter-species co-modules of organs and genes. We identify such co-modules using mouse and human microarray expression data. They are functionally coherent both in terms of genes and of organs from both organisms. We show that a large proportion of genes belonging to the same co-module are orthologous between mouse and human. Moreover, their zebrafish orthologs also tend to be expressed in the corresponding homologous organs. Notable exceptions to the general pattern of conservation are the testis and the olfactory bulb. Interestingly, some co-modules consist of single organs, while others combine several functionally related organs. For instance, amygdala, cerebral cortex, hypothalamus and spinal cord form a clearly discernible unit of expression, both in mouse and human. CONCLUSIONS: Our study provides a new framework for comparative analysis which will be applicable also to other sets of large-scale phenotypic data collected across different species.

Patterns of positive selection in seven ant genomes.

The evolution of ants is marked by remarkable adaptations that allowed the development of very complex social systems. To identify how ant-specific adaptations are associated with patterns of molecular evolution, we searched for signs of positive selection on amino-acid changes in proteins. We identified 24 functional categories of genes which were enriched for positively selected genes in the ant lineage. We also reanalyzed genome-wide data sets in bees and flies with the same methodology to check whether positive selection was specific to ants or also present in other insects. Notably, genes implicated in immunity were enriched for positively selected genes in the three lineages, ruling out the hypothesis that the evolution of hygienic behaviors in social insects caused a major relaxation of selective pressure on immune genes. Our scan also indicated that genes implicated in neurogenesis and olfaction started to undergo increased positive selection before the evolution of sociality in Hymenoptera. Finally, the comparison between these three lineages allowed us to pinpoint molecular evolution patterns that were specific to the ant lineage. In particular, there was ant-specific recurrent positive selection on genes with mitochondrial functions, suggesting that mitochondrial activity was improved during the evolution of this lineage. This might have been an important step toward the evolution of extreme lifespan that is a hallmark of ants.

The amphioxus genome illuminates vertebrate origins and cephalochordate biology.

Cephalochordates, urochordates, and vertebrates evolved from a common ancestor over 520 million years ago. To improve our understanding of chordate evolution and the origin of vertebrates, we intensively searched for particular genes, gene families, and conserved noncoding elements in the sequenced genome of the cephalochordate Branchiostoma floridae, commonly called amphioxus or lancelets. Special attention was given to homeobox genes, opsin genes, genes involved in neural crest development, nuclear receptor genes, genes encoding components of the endocrine and immune systems, and conserved cis-regulatory enhancers. The amphioxus genome contains a basic set of chordate genes involved in development and cell signaling, including a fifteenth Hox gene. This set includes many genes that were co-opted in vertebrates for new roles in neural crest development and adaptive immunity. However, where amphioxus has a single gene, vertebrates often have two, three, or four paralogs derived from two whole-genome duplication events. In addition, several transcriptional enhancers are conserved between amphioxus and vertebrates--a very wide phylogenetic distance. In contrast, urochordate genomes have lost many genes, including a diversity of homeobox families and genes involved in steroid hormone function. The amphioxus genome also exhibits derived features, including duplications of opsins and genes proposed to function in innate immunity and endocrine systems. Our results indicate that the amphioxus genome is elemental to an understanding of the biology and evolution of nonchordate deuterostomes, invertebrate chordates, and vertebrates.

Interactome mapping of the phosphatidylinositol 3-kinase-mammalian target of rapamycin pathway identifies deformed epidermal autoregulatory factor-1 as a new glycogen synthase kinase-3 interactor.

The phosphatidylinositol 3-kinase-mammalian target of rapamycin (PI3K-mTOR) pathway plays pivotal roles in cell survival, growth, and proliferation downstream of growth factors. Its perturbations are associated with cancer progression, type 2 diabetes, and neurological disorders. To better understand the mechanisms of action and regulation of this pathway, we initiated a large scale yeast two-hybrid screen for 33 components of the PI3K-mTOR pathway. Identification of 67 new interactions was followed by validation by co-affinity purification and exhaustive literature curation of existing information. We provide a nearly complete, functionally annotated interactome of 802 interactions for the PI3K-mTOR pathway. Our screen revealed a predominant place for glycogen synthase kinase-3 (GSK3) A and B and the AMP-activated protein kinase. In particular, we identified the deformed epidermal autoregulatory factor-1 (DEAF1) transcription factor as an interactor and in vitro substrate of GSK3A and GSK3B. Moreover, GSK3 inhibitors increased DEAF1 transcriptional activity on the 5-HT1A serotonin receptor promoter. We propose that DEAF1 may represent a therapeutic target of lithium and other GSK3 inhibitors used in bipolar disease and depression.

Neofunctionalization in vertebrates: the example of retinoic acid receptors.

Understanding the role of gene duplications in establishing vertebrate innovations is one of the main challenges of Evo-Devo (evolution of development) studies. Data on evolutionary changes in gene expression (i.e., evolution of transcription factor-cis-regulatory elements relationships) tell only part of the story; protein function, best studied by biochemical and functional assays, can also change. In this study, we have investigated how gene duplication has affected both the expression and the ligand-binding specificity of retinoic acid receptors (RARs), which play a major role in chordate embryonic development. Mammals have three paralogous RAR genes--RAR alpha, beta, and gamma--which resulted from genome duplications at the origin of vertebrates. By using pharmacological ligands selective for specific paralogues, we have studied the ligand-binding capacities of RARs from diverse chordates species. We have found that RAR beta-like binding selectivity is a synapomorphy of all chordate RARs, including a reconstructed synthetic RAR representing the receptor present in the ancestor of chordates. Moreover, comparison of expression patterns of the cephalochordate amphioxus and the vertebrates suggests that, of all the RARs, RAR beta expression has remained most similar to that of the ancestral RAR. On the basis of these results together, we suggest that while RAR beta kept the ancestral RAR role, RAR alpha and RAR gamma diverged both in ligand-binding capacity and in expression patterns. We thus suggest that neofunctionalization occurred at both the expression and the functional levels to shape RAR roles during development in vertebrates.

The immune space: a concept and template for rationalizing vaccine development.

Abstract Empirical testing of candidate vaccines has led to the successful development of a number of lifesaving vaccines. The advent of new tools to manipulate antigens and new methods and vectors for vaccine delivery has led to a veritable explosion of potential vaccine designs. As a result, selection of candidate vaccines suitable for large-scale efficacy testing has become more challenging. This is especially true for diseases such as dengue, HIV, and tuberculosis where there is no validated animal model or correlate of immune protection. Establishing guidelines for the selection of vaccine candidates for advanced testing has become a necessity. A number of factors could be considered in making these decisions, including, for example, safety in animal and human studies, immune profile, protection in animal studies, production processes with product quality and stability, availability of resources, and estimated cost of goods. The "immune space template" proposed here provides a standardized approach by which the quality, level, and durability of immune responses elicited in early human trials by a candidate vaccine can be described. The immune response profile will demonstrate if and how the candidate is unique relative to other candidates, especially those that have preceded it into efficacy testing and, thus, what new information concerning potential immune correlates could be learned from an efficacy trial. A thorough characterization of immune responses should also provide insight into a developer's rationale for the vaccine's proposed mechanism of action. HIV vaccine researchers plan to include this general approach in up-selecting candidates for the next large efficacy trial. This "immune space" approach may also be applicable to other vaccine development endeavors where correlates of vaccine-induced immune protection remain unknown.