The role of the transcription factor Tcf-1 for the development and the function of NK cells

Natural Killer (NK) cells are innate immune cells that can eliminate malignant and foreign cells and that play an important role for the early control of viral and fungal infections. Further, they are important regulators of the adaptive and innate immune responses. During their development in the bone marrow (BM) NK cells undergo several maturation steps that directly establish an effector program. The transcriptional network that controls NK cell development and maturation is still incompletely understood. Based on earlier findings that NK cell numbers are reduced in the absence of the transcription factor T cell factor-1 (Tcf-1), my thesis has addressed the precise role of this transcription factor for NK cell development, maturation and function and whether Tcf-1 acts as a nuclear effector of the canonical Wnt signaling pathway to mediate its effects. It is shown that Tcf-1 is selectively required for the emergence of mature BM NK cells. Surprisingly, the emergence of BM NK cells depends on the repressor function of Tcf-1 and is independent of the Wnt pathway. In BM and peripheral NK cells Tcf-1 is found to suppress Granzyme B (GzmB) expression, a key cytotoxic effector molecule required to kill target cells. We provide evidence that GzmB over-expression in the absence of Tcf-1 results in accelerated spontaneous death of bone marrow NK cells and of cytokine stimulated peripheral NK cells. Moreover, Tcf-1 deficient NK cells show reduced target cell killing, which is due to enhanced GzmB-dependent NK cell death induced by the recognition of tumour target cells. Collectively, these data provide significant new insights into the transcriptional regulation of NK cell development and function and suggest a novel mechanism that protects NK cells from the deleterious effects of highly cytotoxic effector molecules. - Les cellules NK (de l'anglais Natural Killer) font partie du système immunitaire inné et sont capables d'éliminer à elles seules les cellules cancéreuses ou infectées. Ces cellules participent dans la régulation et la coordination des réponses innée et adaptative. Lors de leur développement dans la moelle osseuse, les cellules NK vont acquérir leurs fonctions effectrices, un processus contrôlé par des facteurs de transcription mais encore peu connu. Des précédentes travaux ont montré qu'une diminution du nombre de cellules NK corrélait avec l'absence du facteur de transcription Tcf-1 (T cell factor-1), suggérant un rôle important de Tcf-1 dans le développement de cellules NK. Cette thèse a pour but de mieux comprendre le rôle du facteur de transcription Tcf-1 lors du développement et la maturation des cellules NK, ainsi que son interaction avec la voie de signalisation Wnt. Nous avons montré que Tcf-1 est essentiel pour la transition des cellules immatures NK (iNK) à des cellules matures NK (mNK) dans la moelle osseuse, et cela de manière indépendamment de la voie de signalisation Wnt. De manière intéressante, nous avons observé qu'en absence du facteur de transcription Tcf-1, les cellules NK augmentaient l'expression de la protéine Granzyme B (GzmB), une protéine essentielle pour l'élimination des cellules cancéreuses ou infectées. Ceci a pour conséquence, une augmentation de la mort des cellules mNK dans la moelle osseuse ainsi qu'une diminution de leur fonction «tueuses». Ces résultats montrent pour la première fois, le rôle répresseur du facteur de transcription Tcf-1 dans l'expression de la protéine GzmB. L'ensemble de ces résultats apporte de nouveaux éléments concernant le rôle de Tcf-1 dans la régulation du développement et de la fonction des cellules NK et suggèrent un nouveau mécanisme cellulaire de protection contre les effets délétères d'une dérégulation de l'expression des molécules cytotoxique.

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.

Tax reforms and inequality: theoretical and empirical implications

In this paper we examine the effect of tax policy on the relationship between inequality and growth in a two-sector non-scale model. With non-scale models, the longrun equilibrium growth rate is determined by technological parameters and it is independent of macroeconomic policy instruments. However, this fact does not imply that fiscal policy is unimportant for long-run economic performance. It indeed has important effects on the different levels of key economic variables such as per capita stock of capital and output. Hence, although the economy grows at the same rate across steady states, the bases for economic growth may be different.The model has three essential features. First, we explicitly model skill accumulation, second, we introduce government finance into the production function, and we introduce an income tax to mirror the fiscal events of the 1980¿s and 1990¿s in the US. The fact that the non-scale model is associated with higher order dynamics enables it to replicate the distinctly non-linear nature of inequality in the US with relative ease. The results derived in this paper attract attention to the fact that the non-scale growth model does not only fit the US data well for the long-run (Jones, 1995b) but also that it possesses unique abilities in explaining short term fluctuations of the economy. It is shown that during transition the response of the relative simulated wage to changes in the tax code is rather non-monotonic, quite in accordance to the US inequality pattern in the 1980¿s and early 1990¿s.More specifically, we have analyzed in detail the dynamics following the simulation of an isolated tax decrease and an isolated tax increase. So, after a tax decrease the skill premium follows a lower trajectory than the one it would follow without a tax decrease. Hence we are able to reduce inequality for several periods after the fiscal shock. On the contrary, following a tax increase, the evolution of the skill premium remains above the trajectory carried on by the skill premium under a situation with no tax increase. Consequently, a tax increase would imply a higher level of inequality in the economy

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.

KCM, a general framework for codon-based model of evolution

With the advancement of high-throughput sequencing and dramatic increase of available genetic data, statistical modeling has become an essential part in the field of molecular evolution. Statistical modeling results in many interesting discoveries in the field, from detection of highly conserved or diverse regions in a genome to phylogenetic inference of species evolutionary history Among different types of genome sequences, protein coding regions are particularly interesting due to their impact on proteins. The building blocks of proteins, i.e. amino acids, are coded by triples of nucleotides, known as codons. Accordingly, studying the evolution of codons leads to fundamental understanding of how proteins function and evolve. The current codon models can be classified into three principal groups: mechanistic codon models, empirical codon models and hybrid ones. The mechanistic models grasp particular attention due to clarity of their underlying biological assumptions and parameters. However, they suffer from simplified assumptions that are required to overcome the burden of computational complexity. The main assumptions applied to the current mechanistic codon models are (a) double and triple substitutions of nucleotides within codons are negligible, (b) there is no mutation variation among nucleotides of a single codon and (c) assuming HKY nucleotide model is sufficient to capture essence of transition- transversion rates at nucleotide level. In this thesis, I develop a framework of mechanistic codon models, named KCM-based model family framework, based on holding or relaxing the mentioned assumptions. Accordingly, eight different models are proposed from eight combinations of holding or relaxing the assumptions from the simplest one that holds all the assumptions to the most general one that relaxes all of them. The models derived from the proposed framework allow me to investigate the biological plausibility of the three simplified assumptions on real data sets as well as finding the best model that is aligned with the underlying characteristics of the data sets. -- Avec l'avancement de séquençage à haut débit et l'augmentation dramatique des données géné¬tiques disponibles, la modélisation statistique est devenue un élément essentiel dans le domaine dé l'évolution moléculaire. Les résultats de la modélisation statistique dans de nombreuses découvertes intéressantes dans le domaine de la détection, de régions hautement conservées ou diverses dans un génome de l'inférence phylogénétique des espèces histoire évolutive. Parmi les différents types de séquences du génome, les régions codantes de protéines sont particulièrement intéressants en raison de leur impact sur les protéines. Les blocs de construction des protéines, à savoir les acides aminés, sont codés par des triplets de nucléotides, appelés codons. Par conséquent, l'étude de l'évolution des codons mène à la compréhension fondamentale de la façon dont les protéines fonctionnent et évoluent. Les modèles de codons actuels peuvent être classés en trois groupes principaux : les modèles de codons mécanistes, les modèles de codons empiriques et les hybrides. Les modèles mécanistes saisir une attention particulière en raison de la clarté de leurs hypothèses et les paramètres biologiques sous-jacents. Cependant, ils souffrent d'hypothèses simplificatrices qui permettent de surmonter le fardeau de la complexité des calculs. Les principales hypothèses retenues pour les modèles actuels de codons mécanistes sont : a) substitutions doubles et triples de nucleotides dans les codons sont négligeables, b) il n'y a pas de variation de la mutation chez les nucléotides d'un codon unique, et c) en supposant modèle nucléotidique HKY est suffisant pour capturer l'essence de taux de transition transversion au niveau nucléotidique. Dans cette thèse, je poursuis deux objectifs principaux. Le premier objectif est de développer un cadre de modèles de codons mécanistes, nommé cadre KCM-based model family, sur la base de la détention ou de l'assouplissement des hypothèses mentionnées. En conséquence, huit modèles différents sont proposés à partir de huit combinaisons de la détention ou l'assouplissement des hypothèses de la plus simple qui détient toutes les hypothèses à la plus générale qui détend tous. Les modèles dérivés du cadre proposé nous permettent d'enquêter sur la plausibilité biologique des trois hypothèses simplificatrices sur des données réelles ainsi que de trouver le meilleur modèle qui est aligné avec les caractéristiques sous-jacentes des jeux de données. Nos expériences montrent que, dans aucun des jeux de données réelles, tenant les trois hypothèses mentionnées est réaliste. Cela signifie en utilisant des modèles simples qui détiennent ces hypothèses peuvent être trompeuses et les résultats de l'estimation inexacte des paramètres. Le deuxième objectif est de développer un modèle mécaniste de codon généralisée qui détend les trois hypothèses simplificatrices, tandis que d'informatique efficace, en utilisant une opération de matrice appelée produit de Kronecker. Nos expériences montrent que sur un jeux de données choisis au hasard, le modèle proposé de codon mécaniste généralisée surpasse autre modèle de codon par rapport à AICc métrique dans environ la moitié des ensembles de données. En outre, je montre à travers plusieurs expériences que le modèle général proposé est biologiquement plausible.

Preoperative assessment of neurovesical function in children with anorectal malformation: association with vertebral and spinal malformations.

PURPOSE: We preoperatively assessed neurovesical function and spinal cord function in children with anorectal malformations. In cases of neurovesical dysfunction we looked for an association with vertebral malformation or myelodysplasia. MATERIALS AND METHODS: We prospectively evaluated 80 children with anorectal malformations via preoperative urodynamics and magnetic resonance imaging of the spine. Bladder compliance and volume, detrusor activity and vesicosphincteric synergy during voiding allowed urodynamic evaluation. Results were reported according to Wingspread and Krickenbeck classifications of anorectal malformations. RESULTS: Urodynamic findings were pathological in 14 children (18%). Pathological evaluations did not seem related to type of fistula or level of anorectal malformation. Vertebral anomalies were seen in 34 patients (43%) and myelodysplasia in 16 (20%). Neither vertebral anomaly nor myelodysplasia seemed associated with type of fistula or severity of anorectal malformation. Of 14 children with pathological urodynamics no vertebral anomaly or myelodysplasia was found in 7. Of 66 children with normal urodynamics 40 presented with vertebral or spinal malformation. CONCLUSIONS: Lower urinary tract dysfunction is common in patients with anorectal malformations. Normal spine or spinal cord does not exclude neurovesical dysfunction. Myelodysplasia or vertebral anomaly does not determine lower urinary tract dysfunction. Thus, we recommend preoperative urodynamic assessment of the bladder and magnetic resonance imaging of the spine in children with anorectal malformations.

Evolution of proteins after whole-genome duplication

Evolution of proteins after whole-genome duplicationGene and genome duplication are considered major mechanisms in the creation of newfunctions in genomes, or in the refinement of networks by the division of function amongmore genes. In animals, the best demonstrated whole genome duplication occurred at theorigin of Teleost fishes. This makes fishes an ideal model to study the consequences ofgenome duplication, particularly since we have a good sampling of genome sequences,abundant functional information, and a very well studied outgroup: the tetrapodes (includinghuman). More specifically, I studied the consequences of duplication on proteins usingevolutionary models to infer adaptive events. I analysed the influence of positive selection invertebrate genes, by contrasting singleton genes and duplicated genes. The conclusion of theanalyses was threefold: (i) positive selection affects diverse phylogenetic branches anddiverse gene categories during vertebrate evolution; (ii) it concerns only a small proportion ofsites (1%-5%); and (iii) whole genome duplication had no detectable impact on theprevalence of this positive selection.I also studied evolution at the amino acid level with different methods to detect functionalshifts (covarion process and constant-but-different process). As in my previous research, Ifound similar numbers of functional shifts between duplicates and between orthologs.The accepted framework for studies of molecular evolution is that orthologs share the samefunction, whereas the function of paralogs diverges. This framework gives a special place togene duplication in evolution, as the main mechanism for generating novelty. With myprevious results showing that duplication and speciation are not so different, we investigatedthe literature to question the evidence for similar or divergent evolution of gene function afterduplication relative to speciation genes. This led us to propose a more rigorous design offuture studies of gene duplication.Finally, based on my automated protocol, we built a database of positive selection invertebrates' genes, Selectome. This database is freely available on the web and will helpfuture evolutionary as well as biochemical studies.

Life-history evolution and the polyphenic regulation of somatic maintenance and survival.

Here we discuss life-history evolution from the perspective of adaptive phenotypic plasticity, with a focus on polyphenisms for somatic maintenance and survival. Polyphenisms are adaptive discrete alternative phenotypes that develop in response to changes in the environment. We suggest that dauer larval diapause and its associated adult phenotypes in the nematode (Caenorhabditis elegans), reproductive dormancy in the fruit fly (Drosophila melanogaster) and other insects, and the worker castes of the honey bee (Apis mellifera) are examples of what may be viewed as the polyphenic regulation of somatic maintenance and survival. In these and other cases, the same genotype can--depending upon its environment--express either of two alternative sets of life-history phenotypes that differ markedly with respect to somatic maintenance, survival ability, and thus life span. This plastic modulation of somatic maintenance and survival has traditionally been underappreciated by researchers working on aging and life history. We review the current evidence for such adaptive life-history switches and their molecular regulation and suggest that they are caused by temporally and/or spatially varying, stressful environments that impose diversifying selection, thereby favoring the evolution of plasticity of somatic maintenance and survival under strong regulatory control. By considering somatic maintenance and survivorship from the perspective of adaptive life-history switches, we may gain novel insights into the mechanisms and evolution of aging.

The expansion of amino-acid repeats is not associated to adaptive evolution in mammalian genes.

BACKGROUND: The expansion of amino acid repeats is determined by a high mutation rate and can be increased or limited by selection. It has been suggested that recent expansions could be associated with the potential of adaptation to new environments. In this work, we quantify the strength of this association, as well as the contribution of potential confounding factors. RESULTS: Mammalian positively selected genes have accumulated more recent amino acid repeats than other mammalian genes. However, we found little support for an accelerated evolutionary rate as the main driver for the expansion of amino acid repeats. The most significant predictors of amino acid repeats are gene function and GC content. There is no correlation with expression level. CONCLUSIONS: Our analyses show that amino acid repeat expansions are causally independent from protein adaptive evolution in mammalian genomes. Relaxed purifying selection or positive selection do not associate with more or more recent amino acid repeats. Their occurrence is slightly favoured by the sequence context but mainly determined by the molecular function of the gene.

Origins, evolution, and phenotypic impact of new genes.

Ever since the pre-molecular era, the birth of new genes with novel functions has been considered to be a major contributor to adaptive evolutionary innovation. Here, I review the origin and evolution of new genes and their functions in eukaryotes, an area of research that has made rapid progress in the past decade thanks to the genomics revolution. Indeed, recent work has provided initial whole-genome views of the different types of new genes for a large number of different organisms. The array of mechanisms underlying the origin of new genes is compelling, extending way beyond the traditionally well-studied source of gene duplication. Thus, it was shown that novel genes also regularly arose from messenger RNAs of ancestral genes, protein-coding genes metamorphosed into new RNA genes, genomic parasites were co-opted as new genes, and that both protein and RNA genes were composed from scratch (i.e., from previously nonfunctional sequences). These mechanisms then also contributed to the formation of numerous novel chimeric gene structures. Detailed functional investigations uncovered different evolutionary pathways that led to the emergence of novel functions from these newly minted sequences and, with respect to animals, attributed a potentially important role to one specific tissue--the testis--in the process of gene birth. Remarkably, these studies also demonstrated that novel genes of the various types significantly impacted the evolution of cellular, physiological, morphological, behavioral, and reproductive phenotypic traits. Consequently, it is now firmly established that new genes have indeed been major contributors to the origin of adaptive evolutionary novelties.

Monitoring the clinical introduction of a glutamine and antioxidant solution in critically ill trauma and burn patients

OBJECTIVE: Enteral glutamine supplementation and antioxidants have been shown to be beneficial in some categories of critically ill patients. This study investigated the impact on organ function and clinical outcome of an enteral solution enriched with glutamine and antioxidant micronutrients in patients with trauma and with burns. METHODS: This was a prospective study of a historical control group including critically ill, burned and major trauma patients (n = 86, 40 patients with burns and 46 with trauma, 43 in each group) on admission to an intensive care unit in a university hospital (matching for severity, age, and sex). The intervention aimed to deliver a 500-mL enteral solution containing 30 g of glutamine per day, selenium, zinc, and vitamin E (Gln-AOX) for a maximum of 10 d, in addition to control treatment consisting of enteral nutrition in all patients and intravenous trace elements in all burn patients. RESULTS: Patients were comparable at baseline, except for more inhalation injuries in the burn-Gln-AOX group (P = 0.10) and greater neurologic impairment in the trauma-Gln-AOX group (P = 0.022). Intestinal tolerance was good. The full 500-mL dose was rarely delivered, resulting in a low mean glutamine daily dose (22 g for burn patients and 16 g for trauma patients). In burn patients intravenous trace element delivery was superior to the enteral dose. The evolution of the Sequential Organ Failure Assessment score and other outcome variables did not differ significantly between groups. C-reactive protein decreased faster in the Gln-AOX group. CONCLUSION: The Gln-AOX supplement was well tolerated in critically ill, injured patients, but did not improve outcome significantly. The delivery of glutamine below the 0.5-g/kg recommended dose in association with high intravenous trace element substitution doses in burn patients are likely to have blunted the impact by not reaching an efficient treatment dose. Further trials testing higher doses of Gln are required.

MHC class I expression dependent on bacterial infection and parental factors in whitefish embryos (Salmonidae).

Ecological conditions can influence not only the expression of a phenotype, but also the heritability of a trait. As such, heritable variation for a trait needs to be studied across environments. We have investigated how pathogen challenge affects the expression of MHC genes in embryos of the lake whitefish Coregonus palaea. In order to experimentally separate paternal (i.e. genetic) from maternal and environmental effects, and determine whether and how stress affects the heritable variation for MHC expression, embryos were produced in full-factorial in vitro fertilizations, reared singly, and exposed at 208 degree days (late-eyed stage) to either one of two strains of Pseudomonas fluorescens that differ in their virulence characteristics (one increased mortality, while both delayed hatching time). Gene expression was assessed 48 h postinoculation, and virulence effects of the bacterial infection were monitored until hatching. We found no evidence of MHC class II expression at this stage of development. MHC class I expression was markedly down-regulated in reaction to both pseudomonads. While MHC expression could not be linked to embryo survival, the less the gene was expressed, the earlier the embryos hatched within each treatment group, possibly due to trade-offs between immune function and developmental rate or further factors that affect both hatching timing and MHC expression. We found significant additive genetic variance for MHC class I expression in some treatments. That is, changes in pathogen pressures could induce rapid evolution in MHC class I expression. However, we found no additive genetic variance in reaction norms in our study population.

Evolution of gene expression in fire ants: the effects of developmental stage, caste, and species.

Ants provide remarkable examples of equivalent genotypes developing into divergent and discrete phenotypes. Diploid eggs can develop either into queens, which specialize in reproduction, or workers, which participate in cooperative tasks such as building the nest, collecting food, and rearing the young. In contrast, the differentiation between males and females generally depends upon whether eggs are fertilized, with fertilized (diploid) eggs giving rise to females and unfertilized (haploid) eggs giving rise to males. To obtain a comprehensive picture of the relative contributions of gender (sex), caste, developmental stage, and species divergence to gene expression evolution, we investigated gene expression patterns in pupal and adult queens, workers, and males of two species of fire ants, Solenopsis invicta and S. richteri. Microarray hybridizations revealed that variation in gene expression profiles is influenced more by developmental stage than by caste membership, sex, or species identity. The second major contributor to variation in gene expression was the combination of sex and caste. Although workers and queens share equivalent diploid nuclear genomes, they have highly distinctive patterns of gene expression in both the pupal and the adult stages, as might be expected given their extraordinary level of phenotypic differentiation. Overall, the difference in the proportion of differentially expressed genes was greater between workers and males than between workers and queens or queens and males, consistent with the fact that workers and males share neither gender nor reproductive capability. Moreover, between-species comparisons revealed that the greatest difference in gene expression patterns occurred in adult workers, a finding consistent with the fact that adult workers most directly experience the distinct external environments characterizing the different habitats occupied by the two species. Thus, much of the evolution of gene expression in ants may occur in the worker caste, despite the fact that these individuals are largely or completely sterile. Analyses of gene expression evolution revealed a combination of positive selection and relaxation of stabilizing selection as important factors driving the evolution of such genes.

Pleiotropy in the melanocortin system: expression levels of this system are associated with melanogenesis and pigmentation in the tawny owl (Strix aluco).

The adaptive function of melanin-based coloration is a long-standing debate. A recent genetic model suggested that pleiotropy could account for covariations between pigmentation, behaviour, morphology, physiology and life history traits. We explored whether the expression levels of genes belonging to the melanocortin system (MC1R, POMC, PC1/3, PC2 and the antagonist ASIP), which have many pleiotropic effects, are associated with melanogenesis (through variation in the expression of the genes MITF, SLC7A11, TYR, TYRP1) and in turn melanin-based coloration. We considered the tawny owl (Strix aluco) because individuals vary continuously from light to dark reddish, and thus, colour variation is likely to stem from differences in the levels of gene expression. We measured gene expression in feather bases collected in nestlings at the time of melanin production. As expected, the melanocortin system was associated with the expression of melanogenic genes and pigmentation. Offspring of darker reddish fathers expressed PC1/3 to lower levels but tended to express PC2 to higher levels. The convertase enzyme PC1/3 cleaves the POMC prohormone to obtain ACTH, while the convertase enzyme PC2 cleaves ACTH to produce α-melanin-stimulating hormone (α-MSH). ACTH regulates glucocorticoids, hormones that modulate stress responses, while α-MSH induces eumelanogenesis. We therefore conclude that the melanocortin system, through the convertase enzymes PC1/3 and PC2, may account for part of the interindividual variation in melanin-based coloration in nestling tawny owls. Pleiotropy may thus account for the covariation between phenotypic traits involved in social interactions (here pigmentation) and life history, morphology, behaviour and physiology.

Gene Ontology Function prediction in Mollicutes using Protein-Protein Association Networks

Abstract Background: Many complex systems can be represented and analysed as networks. The recent availability of large-scale datasets, has made it possible to elucidate some of the organisational principles and rules that govern their function, robustness and evolution. However, one of the main limitations in using protein-protein interactions for function prediction is the availability of interaction data, especially for Mollicutes. If we could harness predicted interactions, such as those from a Protein-Protein Association Networks (PPAN), combining several protein-protein network function-inference methods with semantic similarity calculations, the use of protein-protein interactions for functional inference in this species would become more potentially useful. Results: In this work we show that using PPAN data combined with other approximations, such as functional module detection, orthology exploitation methods and Gene Ontology (GO)-based information measures helps to predict protein function in Mycoplasma genitalium. Conclusions: To our knowledge, the proposed method is the first that combines functional module detection among species, exploiting an orthology procedure and using information theory-based GO semantic similarity in PPAN of the Mycoplasma species. The results of an evaluation show a higher recall than previously reported methods that focused on only one organism network.