From SNPs to pathways: integration of functional effect of sequence variations on models of cell signalling pathways

Background: Single nucleotide polymorphisms (SNPs) are the most frequent type of sequence variation between individuals, and represent a promising tool for finding genetic determinants of complex diseases and understanding the differences in drug response. In this regard, it is of particular interest to study the effect of non-synonymous SNPs in the context of biological networks such as cell signalling pathways. UniProt provides curated information about the functional and phenotypic effects of sequence variation, including SNPs, as well as on mutations of protein sequences. However, no strategy has been developed to integrate this information with biological networks, with the ultimate goal of studying the impact of the functional effect of SNPs in the structure and dynamics of biological networks. Results: First, we identified the different challenges posed by the integration of the phenotypic effect of sequence variants and mutations with biological networks. Second, we developed a strategy for the combination of data extracted from public resources, such as UniProt, NCBI dbSNP, Reactome and BioModels. We generated attribute files containing phenotypic and genotypic annotations to the nodes of biological networks, which can be imported into network visualization tools such as Cytoscape. These resources allow the mapping and visualization of mutations and natural variations of human proteins and their phenotypic effect on biological networks (e.g. signalling pathways, protein-protein interaction networks, dynamic models). Finally, an example on the use of the sequence variation data in the dynamics of a network model is presented. Conclusion: In this paper we present a general strategy for the integration of pathway and sequence variation data for visualization, analysis and modelling purposes, including the study of the functional impact of protein sequence variations on the dynamics of signalling pathways. This is of particular interest when the SNP or mutation is known to be associated to disease. We expect that this approach will help in the study of the functional impact of disease-associated SNPs on the behaviour of cell signalling pathways, which ultimately will lead to a better understanding of the mechanisms underlying complex diseases.

Genetic and genomic analysis modeling of germline c-MYC overexpression and cancer susceptibility

Background: Germline genetic variation is associated with the differential expression of many human genes. The phenotypic effects of this type of variation may be important when considering susceptibility to common genetic diseases. Three regions at 8q24 have recently been identified to independently confer risk of prostate cancer. Variation at 8q24 has also recently been associated with risk of breast and colorectal cancer. However, none of the risk variants map at or relatively close to known genes, with c-MYC mapping a few hundred kilobases distally. Results: This study identifies cis-regulators of germline c-MYC expression in immortalized lymphocytes of HapMap individuals. Quantitative analysis of c-MYC expression in normal prostate tissues suggests an association between overexpression and variants in Region 1 of prostate cancer risk. Somatic c-MYC overexpression correlates with prostate cancer progression and more aggressive tumor forms, which was also a pathological variable associated with Region 1. Expression profiling analysis and modeling of transcriptional regulatory networks predicts a functional association between MYC and the prostate tumor suppressor KLF6. Analysis of MYC/Myc-driven cell transformation and tumorigenesis substantiates a model in which MYC overexpression promotes transformation by down-regulating KLF6. In this model, a feedback loop through E-cadherin down-regulation causes further transactivation of c-MYC.Conclusion: This study proposes that variation at putative 8q24 cis-regulator(s) of transcription can significantly alter germline c-MYC expression levels and, thus, contribute to prostate cancer susceptibility by down-regulating the prostate tumor suppressor KLF6 gene.

OSIRISv1.2: a named entity recognition system for sequence variants of genes in biomedical literature

Background: Single Nucleotide Polymorphisms, among other type of sequence variants, constitute key elements in genetic epidemiology and pharmacogenomics. While sequence data about genetic variation is found at databases such as dbSNP, clues about the functional and phenotypic consequences of the variations are generally found in biomedical literature. The identification of the relevant documents and the extraction of the information from them are hampered by the large size of literature databases and the lack of widely accepted standard notation for biomedical entities. Thus, automatic systems for the identification of citations of allelic variants of genes in biomedical texts are required. Results: Our group has previously reported the development of OSIRIS, a system aimed at the retrieval of literature about allelic variants of genes http://ibi.imim.es/osirisform.html. Here we describe the development of a new version of OSIRIS (OSIRISv1.2, http://ibi.imim.es/OSIRISv1.2.html webcite) which incorporates a new entity recognition module and is built on top of a local mirror of the MEDLINE collection and HgenetInfoDB: a database that collects data on human gene sequence variations. The new entity recognition module is based on a pattern-based search algorithm for the identification of variation terms in the texts and their mapping to dbSNP identifiers. The performance of OSIRISv1.2 was evaluated on a manually annotated corpus, resulting in 99% precision, 82% recall, and an F-score of 0.89. As an example, the application of the system for collecting literature citations for the allelic variants of genes related to the diseases intracranial aneurysm and breast cancer is presented. Conclusion: OSIRISv1.2 can be used to link literature references to dbSNP database entries with high accuracy, and therefore is suitable for collecting current knowledge on gene sequence variations and supporting the functional annotation of variation databases. The application of OSIRISv1.2 in combination with controlled vocabularies like MeSH provides a way to identify associations of biomedical interest, such as those that relate SNPs with diseases.

The Drosophila cdc25 homolog twine is required for meiosis.

We have identified a second cdc25 homolog in Drosophila. In contrast to string (the first homolog identified in Drosophila) this second homolog, twine, does not function in the mitotic cell cycle, but is specialized for meiosis. Expression of twine was observed exclusively in male and female gonads. twine transcripts are present in germ cells during meiosis, and appear only late during gametogenesis, well after the end of the mitotic germ cell divisions. The sterile Drosophila mutant, mat(2)synHB5, which had previously been isolated and mapped to the same genomic region as twine (35F), was found to carry a missense mutation in the twine gene. This missense mutation in twine abolished its ability to complement a mutation in Schizosaccharomyces pombe cdc25. Phenotypic analysis of mat(2)synHB5 mutant flies revealed a complete block of meiosis in males and severe meiotic defects in females.

Pathogen-induced hatching and population-specific life-history response to water-borne cues in brown trout (Salmo trutta)

Hatching is an important niche shift, and embryos in a wide range of taxa can either accelerate or delay this life-history switch in order to avoid stage-specific risks. Such behavior can occur in response to stress itself and to chemical cues that allow anticipation of stress. We studied the genetic organization of this phenotypic plasticity and tested whether there are differences among populations and across environments in order to learn more about the evolutionary potential of stress-induced hatching. As a study species, we chose the brown trout (Salmo trutta; Salmonidae). Gametes were collected from five natural populations (within one river network) and used for full-factorial in vitro fertilizations. The resulting embryos were either directly infected with Pseudomonas fluorescens or were exposed to waterborne cues from P. fluorescens-infected conspecifics. We found that direct inoculation with P. fluorescens increased embryonic mortality and induced hatching in all host populations. Exposure to waterborne cues revealed population-specific responses. We found significant additive genetic variation for hatching time, and genetic variation in trait plasticity. In conclusion, hatching is induced in response to infection and can be affected by waterborne cues of infection, but populations and families differ in their reaction to the latter.

Heterozygosity-fitness correlations among wild populations of European tree frogs (Hyla arborea) detect fixation load.

Quantifying the impacts of inbreeding and genetic drift on fitness traits in fragmented populations is becoming a major goal in conservation biology. Such impacts occur at different levels and involve different sets of loci. Genetic drift randomly fixes slightly deleterious alleles leading to different fixation load among populations. By contrast, inbreeding depression arises from highly deleterious alleles in segregation within a population and creates variation among individuals. A popular approach is to measure correlations between molecular variation and phenotypic performances. This approach has been mainly used at the individual level to detect inbreeding depression within populations and sometimes at the population level but without consideration about the genetic processes measured. For the first time, we used in this study a molecular approach considering both the interpopulation and intrapopulation level to discriminate the relative importance of inbreeding depression vs. fixation load in isolated and non-fragmented populations of European tree frog (Hyla arborea), complemented with interpopulational crosses. We demonstrated that the positive correlations observed between genetic heterozygosity and larval performances on merged data were mainly caused by co-variations in genetic diversity and fixation load among populations rather than by inbreeding depression and segregating deleterious alleles within populations. Such a method is highly relevant in a conservation perspective because, depending on how populations lose fitness (inbreeding vs. fixation load), specific management actions may be designed to improve the persistence of populations.

Draft genome of the red harvester ant Pogonomyrmex barbatus.

We report the draft genome sequence of the red harvester ant, Pogonomyrmex barbatus. The genome was sequenced using 454 pyrosequencing, and the current assembly and annotation were completed in less than 1 y. Analyses of conserved gene groups (more than 1,200 manually annotated genes to date) suggest a high-quality assembly and annotation comparable to recently sequenced insect genomes using Sanger sequencing. The red harvester ant is a model for studying reproductive division of labor, phenotypic plasticity, and sociogenomics. Although the genome of P. barbatus is similar to other sequenced hymenopterans (Apis mellifera and Nasonia vitripennis) in GC content and compositional organization, and possesses a complete CpG methylation toolkit, its predicted genomic CpG content differs markedly from the other hymenopterans. Gene networks involved in generating key differences between the queen and worker castes (e.g., wings and ovaries) show signatures of increased methylation and suggest that ants and bees may have independently co-opted the same gene regulatory mechanisms for reproductive division of labor. Gene family expansions (e.g., 344 functional odorant receptors) and pseudogene accumulation in chemoreception and P450 genes compared with A. mellifera and N. vitripennis are consistent with major life-history changes during the adaptive radiation of Pogonomyrmex spp., perhaps in parallel with the development of the North American deserts.

Ring recoveries of dead birds confirm that darker pheomelanic Barn Owls disperse longer distances

Variation in melanin coloration is widespread and often associated with other phenotypic traits. A recent study showed that darker-reddish pheomelanic Barn Owls (Tyto alba) move longer distances between birth and breeding sites. Because this study considered only individuals recovered within a limited study area, it remains unclear whether the association between melanism and dispersal applies to a larger geographic scale. I analysed an independent dataset of birds ringed in the same study area but recovered dead along roads within and outside this area. As expected, dark pheomelanic owls dispersed further than lighter reddish conspecifics at a larger spatial scale.

Bidirectional shifts in colony queen number in a socially polymorphic ant population.

The breeding system of social organisms affects many important aspects of social life. Some species vary greatly in the number of breeders per group, but the mechanisms and selective pressures contributing to the maintenance of this polymorphism in social structure remain poorly understood. Here, we take advantage of a genetic dataset that spans 15 years to investigate the dynamics of colony queen number within a socially polymorphic ant species. Our study population of Formica selysi has single- and multiple-queen colonies. We found that the social structure of this species is somewhat flexible: on average, each year 3.2% of the single-queen colonies became polygynous, and conversely 1.4% of the multiple-queen colonies became monogynous. The annualized queen replacement rates were 10.3% and 11.9% for single- and multiple-queen colonies, respectively. New queens were often but not always related to previous colony members. At the population level, the social polymorphism appeared stable. There was no genetic differentiation between single- and multiple-queen colonies at eight microsatellite loci, suggesting ongoing gene flow between social forms. Overall, the regular and bidirectional changes in queen number indicate that social structure is a labile trait in F. selysi, with neither form being favored within a time-frame of 15 years.

The influence of learning on evolution: a mathematical framework.

The Baldwin effect can be observed if phenotypic learning influences the evolutionary fitness of individuals, which can in turn accelerate or decelerate evolutionary change. Evidence for both learning-induced acceleration and deceleration can be found in the literature. Although the results for both outcomes were supported by specific mathematical or simulation models, no general predictions have been achieved so far. Here we propose a general framework to predict whether evolution benefits from learning or not. It is formulated in terms of the gain function, which quantifies the proportional change of fitness due to learning depending on the genotype value. With an inductive proof we show that a positive gain-function derivative implies that learning accelerates evolution, and a negative one implies deceleration under the condition that the population is distributed on a monotonic part of the fitness landscape. We show that the gain-function framework explains the results of several specific simulation models. We also use the gain-function framework to shed some light on the results of a recent biological experiment with fruit flies.

Epidermal growth factor (EGF) stimulation of cultured brain cells. II. Increased production of extracellular soluble proteins.

The production of extracellular soluble proteins was studied in serum-free aggregating cell cultures of fetal rat telencephalon labeled on culture day 7 with a mixture of radioactive amino acid precursors. Cultures treated continuously with epidermal growth factor (EGF; 20 ng/ml) showed a generally increased protein secretion and a particularly enhanced production of a few distinct extracellular proteins. The time lag of this response after an initial dose of EGF (25 ng/ml) on day 7 was 48 h. The total macromolecular radioactivity that accumulated within 96 h of labeling in the media of EGF-treated cultures was 175% of untreated controls, whereas no difference was found in the proportions of intracellular amino acid incorporation. Cultures which received a single dose of EGF (25 ng/ml) on day 1 showed still a greatly increased protein secretion on day 7. Prevention of extracellular protein accumulation by reducing the initial cell number and increasing the rate of media changes did not affect the EGF-induced stimulation of the two glial enzymes, glutamine synthetase and 2',3'-cyclic nucleotide 3'-phosphohydrolase. The results suggest that both the increased production of extracellular proteins and the enhanced development of glial enzymatic activities reflect the stimulated phenotypic expression of EGF-sensitive brain cells.

Testing the predictive adaptive response in a host-parasite system

1. Harsh environmental conditions experienced during development can reduce the performance of the same individuals in adulthood. However, the 'predictive adaptive response' hypothesis postulates that if individuals adapt their phenotype during development to the environments where they are likely to live in the future, individuals exposed to harsh conditions in early life perform better when encountering the same harsh conditions in adulthood compared to those never exposed to these conditions before. 2. Using the common vole (Microtus arvalis) as study organism, we tested how exposure to flea parasitism during the juvenile stage affects the physiology (haematocrit, resistance to oxidative stress, resting metabolism, spleen mass, and testosterone), morphology (body mass, testis mass) and motor performance (open field activity and swimming speed) of the same individuals when infested with fleas in adulthood. According to the 'predictive adaptive response' hypothesis, we predicted that voles parasitized at the adult stage would perform better if they had already been parasitized with fleas at the juvenile stage. 3. We found that voles exposed to fleas in adulthood had a higher metabolic rate if already exposed to fleas when juvenile, compared to voles free of fleas when juvenile and voles free of fleas in adulthood. Independently of juvenile parasitism, adult parasitism impaired adult haematocrit and motor performances. Independently of adult parasitism, juvenile parasitism slowed down crawling speed in adult female voles. 4. Our results suggest that juvenile parasitism has long-term effects that do not protect from the detrimental effects of adult parasitism. On the contrary, experiencing parasitism in early-life incurs additional costs upon adult parasitism measured in terms of higher energy expenditure, rather than inducing an adaptive shift in the developmental trajectory. 5. Hence, our study provides experimental evidence for long term costs of parasitism. We found no support for a predictive adaptive response in this host-parasite system.

High genetic variability and low local diversity in a population of arbuscular mycorrhizal fungi.

Arbuscular mycorrhizal fungi (AMF) are ecologically important root symbionts of most terrestrial plants. Ecological studies of AMF have concentrated on differences between species; largely assuming little variability within AMF species. Although AMF are clonal, they have evolved to contain a surprisingly high within-species genetic variability, and genetically different nuclei can coexist within individual spores. These traits could potentially lead to within-population genetic variation, causing differences in physiology and symbiotic function in AMF populations, a consequence that has been largely neglected. We found highly significant genetic and phenotypic variation among isolates of a population of Glomus intraradices but relatively low total observed genetic diversity. Because we maintained the isolated population in a constant environment, phenotypic variation can be considered as variation in quantitative genetic traits. In view of the large genetic differences among isolates by randomly sampling two individual spores, <50% of the total observed population genetic diversity is represented. Adding an isolate from a distant population did not increase total observed genetic diversity. Genetic variation exceeded variation in quantitative genetic traits, indicating that selection acted on the population to retain similar traits, which might be because of the multigenomic nature of AMF, where considerable genetic redundancy could buffer the effects of changes in the genetic content of phenotypic traits. These results have direct implications for ecological research and for studying AMF genes, improving commercial AMF inoculum, and understanding evolutionary mechanisms in multigenomic organisms.

The adaptative function of melanin-based colour variants

Résumé : Les mécanismes de contrôle des couleurs mélaniques chez les vertébrés sont encore discutés parmi les biologistes de l'évolution. Une hypothèse récente affirme que les effets pléiotropies du système des mélanocortines expliquent l'association fréquente entre la coloration eumélanique noire (due à la déposition d'eumélanine) et de nombreux traits physiologiques et comportementaux. De nombreuses études suggèrent, en effet, que des niveaux plus élevés des mélanocortines induisent l'assombrissement des téguments eumélaniques et affectent d'autres traits phénotypiques simultanément. Cependant, il n'est pas encore établi si ce mécanisme de pléiotropie peut s'appliquer aux colorations dues à la déposition de phaeomélanine, une autre forme commune de mélanine. Les antagonistes des mélanocortines déclenchent le phaeomélanogenèse et bloquent l'effet des mélanocortines ou ont un effet pharmacologique opposé. Nous nous proposons donc d'évaluer l'hypothèse que les effets pléiotropes des antagonistes des mélanocortines génèrent des covariations entre la coloration phaeomélanique et des aspects de la qualité individuelle. Comme prédit par cette hypothèse, nous constatons chez la chouette effraie (Tyto alba) que les traits phénotypiques (résistance au stress oxydatif et aux parasites) corrèlent positivement au degré d'expression d'une couleur eumélanique mais négativement au degré d'expression d'une coloration phaeomélanique. Puis, nous montrons chez la chouette hulotte (Strix aluco) que les associations génétiques entre la coloration phaeomélanique et la physiologie (immunité et la régulation de l'homéostasie) confèrent des avantages aux individus de différentes couleurs dans différents environnements caractérisés par l'abondance de nourriture et le niveau d'exposition aux parasites. Ainsi, nos études soutiennent l'hypothèse que les effets pléiotropes des antagonistes des mélanocortines génèrent des covariations entre les traits mélaniques et divers aspects de la qualité individuelle. Finalement, nous montrons chez le faucon crécerelle (Falco Tinnunculus) que l'expression des ornements mélaniques est sensible à la qualité de l'environnement dans lequel les individus grandissent. Ceci suggère que les gènes codant pour les mélanocortines et leurs antagonistes pourraient induire une expression des traits mélaniques dépendante de la condition de l'individu, un pattern d'expression rarement observé pour des traits généralement sous fort contrôle génétique. Summary : The information content and control mechanisms of melanin-based colour signals in vertebrates are still debated among evolutionary biologists. A recent hypothesis contends that pleiotropic effects of the melanocortin system accounts for the frequent association between black eumelanic coloration and physiological and behavioural traits. Accordingly, empirical evidence suggests that higher levels of melanocortins concurrently promote darker eumelanic integuments and affect other phenotypic traits. However, whether this mechanism may apply to signals relying on phaeomelanin, another common form of melanin pigments, remains to be established. Melanocortin antagonists trigger phaeomelanogenesis and block the effect of melanocortins or result in the opposite pharmacological effect. Therefore, we tested the hypothesis that pleiotropic effects of melanocortin antagonists and inverse agonists account for covariations between phaeomelanin-based coloration and aspects of individual quality. As predicted, we found that phenotypic traits (resistance to oxidative stress and parasites) correlated positively with a eumelanic trait and negatively with a phaeomelanic trait in the barn owl (Tyto alba). Then, we showed in the tawny owl (Strix aluco) that genetic associations between phaeomelanin-based coloration and physiology (immunity and regulation of energy homeostasis) confer benefits to differently coloured individuals under different levels of food abundance and parasite exposure. Altogether, our studies support the hypothesis that pleiotropic effects of melanocortins antagonists can indeed account for covariations between phaeomelanin-based traits and aspects of individual quality. Eventually, we show in the Eurasian kestrel (Falco Tinnunculus) that expression of melanin-based ornaments is sensitive to the environment in which individuals grow. This suggests that genes coding for melanocortins and their antagonists can mediate the condition-dependent component of melanin-based traits.

Population density and offspring size in the cladoceran Simocephalus vetulus (Müller)

r/K theory classically predicts that offspring size should increase under density-dependent selection. However, this is questionable, being based on implicit rather than explicit assumption (the logistic model does not include offsring size as a parameter). From recent models of optimal offspring size (Sibly & Calow, 1983; Taylor & Williams, 1984) it can be shown that density should select for larger offspring if density-dependence in the per capita rate of increase is mainly due to a reduction of the juvenile growth rate or survivorship. In contrast, density should select for smaller offspring if such density-dependence is mainly due to a reduction of adult fecundity or survivorship. Therfore, the outcome of selection cannot be predicted without precise knowledge of the density-dependence of age-specific reproduction and mortality rates. To test the above models, genetically identical individuals of Simocephalus vetulus (Müller) were reared in a density gradient; density-dependence in the per capita rate of increase was shown to be mainly due to a reduction of the juvenile growth rate, thereby selecting for larger offspring; offspring size at birth appeared to be phenotypically plastic and to increase with density. Models were therefore qualitatively supported. However, a discrepancy occurred in quantitative predictions; offspring were produced larger than predicted. Field and laboratory studies are suggested to address this.