The evolutionary rate variation among genes of HOG-signaling pathway in yeast genomes

Transcription factors evolve faster than their structural gene targets in the flavonoid pigment pathway

Dissecting the relationship between gene function and substitution rates is key to understanding genome-wide patterns of molecular evolution. Biochemical pathways provide powerful systems for investigating this relationship because the functional role of each gene is often well characterized. Here, we investigate the evolution of the flavonoid pigment pathway in the colorful Petunieae clade of the tomato family (Solanaceae). This pathway is broadly conserved in plants, both in terms of its structural elements and its MYB, basic helix–loop–helix, and WD40 transcriptional regulators, and its function has been extensively studied, particularly in model species of petunia. We built a phylotranscriptomic data set for 69 species of Petunieae to infer patterns of molecular evolution across pathway genes and across lineages. We found that transcription factors exhibit faster rates of molecular evolution (dN/dS) than their targets, with the highly specialized MYB genes evolving fastest. Using the largest comparative data set to date, we recovered little support for the hypothesis that upstream enzymes evolve slower than those occupying more downstream positions, although expression levels do predict molecular evolutionary rates. Although shifts in floral pigmentation were only weakly related to changes affecting coding regions, we found a strong relationship with the presence/absence patterns of MYB transcripts. Intensely pigmented species express all three main MYB anthocyanin activators in petals, whereas pale or white species express few or none. Our findings reinforce the notion that pathway regulators have a dynamic history, involving higher rates of molecular evolution than structural components, along with frequent changes in expression during color transitions.

History and Domestication of Saccharomyces cerevisiae in Bread Baking

The yeast Saccharomyces cerevisiae has been instrumental in the fermentation of foods and beverages for millennia. In addition to fermentations like wine, beer, cider, sake, and bread, S. cerevisiae has been isolated from environments ranging from soil and trees, to human clinical isolates. Each of these environments has unique selection pressures that S. cerevisiae must adapt to. Bread dough, for example, requires S. cerevisiae to efficiently utilize the complex sugar maltose; tolerate osmotic stress due to the semi-solid state of dough, high salt, and high sugar content of some doughs; withstand various processing conditions, including freezing and drying; and produce desirable aromas and flavors. In this review, we explore the history of bread that gave rise to modern commercial baking yeast, and the genetic and genomic changes that accompanied this. We illustrate the genetic and phenotypic variation that has been documented in baking strains and wild strains, and how this variation might be used for baking strain improvement. While we continue to improve our understanding of how baking strains have adapted to bread dough, we conclude by highlighting some of the remaining open questions in the field.

Gene expression predictions and networks in natural populations supports the omnigenic theory

Background

Recent literature on the differential role of genes within networks distinguishes core from peripheral genes. If previous works have shown contrasting features between them, whether such categorization matters for phenotype prediction remains to be studied.

Results

We measured 17 phenotypic traits for 241 cloned genotypes from a Populus nigra collection, covering growth, phenology, chemical and physical properties. We also sequenced RNA for each genotype and built co-expression networks to define core and peripheral genes. We found that cores were more differentiated between populations than peripherals while being less variable, suggesting that they have been constrained through potentially divergent selection. We also showed that while cores were overrepresented in a subset of genes statistically selected for their capacity to predict the phenotypes (by Boruta algorithm), they did not systematically predict better than peripherals or even random genes.

Conclusion

Our work is the first attempt to assess the importance of co-expression network connectivity in phenotype prediction. While highly connected core genes appear to be important, they do not bear enough information to systematically predict better quantitative traits than other gene sets.

Network inference reveals novel connections in pathways regulating growth and defense in the yeast salt response

Cells respond to stressful conditions by coordinating a complex, multi-faceted response that spans many levels of physiology. Much of the response is coordinated by changes in protein phosphorylation. Although the regulators of transcriptome changes during stress are well characterized in Saccharomyces cerevisiae, the upstream regulatory network controlling protein phosphorylation is less well dissected. Here, we developed a computational approach to infer the signaling network that regulates phosphorylation changes in response to salt stress. We developed an approach to link predicted regulators to groups of likely co-regulated phospho-peptides responding to stress, thereby creating new edges in a background protein interaction network. We then use integer linear programming (ILP) to integrate wild type and mutant phospho-proteomic data and predict the network controlling stress-activated phospho-proteomic changes. The network we inferred predicted new regulatory connections between stress-activated and growth-regulating pathways and suggested mechanisms coordinating metabolism, cell-cycle progression, and growth during stress. We confirmed several network predictions with co-immunoprecipitations coupled with mass-spectrometry protein identification and mutant phospho-proteomic analysis. Results show that the cAMP-phosphodiesterase Pde2 physically interacts with many stress-regulated transcription factors targeted by PKA, and that reduced phosphorylation of those factors during stress requires the Rck2 kinase that we show physically interacts with Pde2. Together, our work shows how a high-quality computational network model can facilitate discovery of new pathway interactions during osmotic stress.

Cells sense and respond to stressful environments by utilizing complex signaling networks that integrate diverse signals to coordinate a multi-faceted physiological response. Much of this response is controlled by post-translational protein phosphorylation. Although many regulators that mediate changes in protein phosphorylation are known, how these regulators inter-connect in a single regulatory network that can transmit cellular signals is not known. It is also unclear how regulators that promote growth and regulators that activate the stress response interconnect to reorganize resource allocation during stress. Here, we developed an integrated experimental and computational workflow to infer the signaling network that regulates phosphorylation changes during osmotic stress in the budding yeast Saccharomyces cerevisiae. The workflow integrates data measuring protein phosphorylation changes in response to osmotic stress with known physical interactions between yeast proteins from large-scale datasets, along with other information about how regulators recognize their targets. The resulting network suggested new signaling connections between regulators and pathways, including those involved in regulating growth and defense, and predicted new regulators involved in stress defense. Our work highlights the power of using network inference to deliver new insight on how cells coordinate a diverse adaptive strategy to stress.

Stronger selective constraint on downstream genes in the oxidative phosphorylation pathway of cetaceans

The oxidative phosphorylation (OXPHOS) pathway is an efficient way to produce energy via adenosine triphosphate (ATP), which is critical for sustaining an energy supply for cetaceans in a hypoxic environment. Several studies have shown that natural selection may shape the evolution of the genes involved in OXPHOS. However, how network architecture drives OXPHOS protein sequence evolution remains poorly explored. Here, we investigated the evolutionary patterns of genes in the OXPHOS pathway across six cetacean genomes within the framework of a functional network. Our results show a negative correlation between the strength of purifying selection and pathway position. This result indicates that downstream genes were subjected to stronger evolutionary constraints than upstream genes, which may be due to the dual function of ATP synthase in the OXPHOS pathway. Additionally, there was a positive correlation between codon usage bias and omega (ω = dN/dS) and a negative correlation with synonymous substitution rate (dS), indicating that the stronger selective constraint on genes (with less biased codon usage) along the OXPHOS pathway is attributable to an increase in the rate of synonymous substitution. Surprisingly, there was no significant correlation between protein-protein interactions and the evolutionary estimates, implying that highly connected enzymes may not always show greater evolutionary constraints. Compared with that observed for terrestrial mammals, we found that the signature of positive selection detected in five genes (ATP5J, LHPP, PPA1, UQCRC1 and UQCRQ) was cetacean-specific, reflecting the importance of OXPHOS for survival in hypoxic, aquatic environments.

Evolutionary patterns of Toll-like receptor signaling pathway genes in the Suidae

BACKGROUND

The Toll-like receptor (TLR) signaling pathway constitutes an essential component of the innate immune system. Highly conserved proteins, indicative of their critical roles in host survival, characterize this pathway. Selective constraints could vary depending on the gene's position within the pathway as TLR signaling is a sequential process and that genes downstream of the TLRs may be more selectively constrained to ensure efficient immune responses given the important role of downstream genes in the signaling process. Thus, we investigated whether gene position influenced protein evolution in the TLR signaling pathway of the Suidae. The members of the Suidae examined included the European Sus scrofa (wild boar), Asian Sus scrofa (wild boar), Sus verrucosus, Sus celebensis, Sus scebifrons, Sus barbatus, Babyrousa babyrussa, Potamochoerus larvatus, Potamochoerus porcus and Phacochoerus africanus.

RESULTS

A total of 33 TLR signaling pathway genes in the Suidae were retrieved from resequencing data. The evolutionary parameter ω (dn/ds) had an overall mean of 0.1668 across genes, indicating high functional conservation within the TLR signaling pathway. A significant relationship was inferred for the network parameters gene position, number of protein-protein interactions, protein length and the evolutionary parameter dn (nonsynonymous substitutions) such that downstream genes had lower nonsynonymous substitution rates, more interactors and shorter protein length than upstream genes. Gene position was significantly correlated with the number of protein-protein interactions and protein length. Thus, the polarity in the selective constraint along the TLR signaling pathway was due to the number of molecules a protein interacted with and the protein's length.

CONCLUSION

Results indicate that the level of selective constraints on genes within the TLR signaling pathway of the Suidae is dependent on the gene's position and network parameters. In particular, downstream genes evolve more slowly as a result of being highly connected and having shorter protein lengths. These findings highlight the critical role of gene network parameters in gene evolution.

Gene Diversification of an Emerging Pathogen: A Decade of Mutation in a Novel Fish Viral Hemorrhagic Septicemia (VHS) Substrain since Its First Appearance in the Laurentian Great Lakes

Viral Hemorrhagic Septicemia virus (VHSv) is an RNA rhabdovirus, which causes one of the world's most serious fish diseases, infecting >80 freshwater and marine species across the Northern Hemisphere. A new, novel, and especially virulent substrain—VHSv-IVb—first appeared in the Laurentian Great Lakes about a decade ago, resulting in massive fish kills. It rapidly spread and has genetically diversified. This study analyzes temporal and spatial mutational patterns of VHSv-IVb across the Great Lakes for the novel non-virion (Nv) gene that is unique to this group of novirhabdoviruses, in relation to its glycoprotein (G), phosphoprotein (P), and matrix (M) genes. Results show that the Nv-gene has been evolving the fastest (k = 2.0x10^-3 substitutions/site/year), with the G-gene at ~1/7 that rate (k = 2.8x10^-4). Most (all but one) of the 12 unique Nv- haplotypes identified encode different amino acids, totaling 26 changes. Among the 12 corresponding G-gene haplotypes, seven vary in amino acids with eight total changes. The P- and M- genes are more evolutionarily conserved, evolving at just ~1/15 (k = 1.2x10^-4) of the Nv-gene’s rate. The 12 isolates contained four P-gene haplotypes with two amino acid changes, and six M-gene haplotypes with three amino acid differences. Patterns of evolutionary changes coincided among the genes for some of the isolates, but appeared independent in others. New viral variants were discovered following the large 2006 outbreak; such differentiation may have been in response to fish populations developing resistance, meriting further investigation. Two 2012 variants were isolated by us from central Lake Erie fish that lacked classic VHSv symptoms, having genetically distinctive Nv-, G-, and M-gene sequences (with one of them also differing in its P-gene); they differ from each other by a G-gene amino acid change and also differ from all other isolates by a shared Nv-gene amino acid change. Such rapid evolutionary differentiation may allow new viral variants to evade fish host recognition and immune responses, facilitating long-time persistence along with expansion to new geographic areas.

Genetic mapping of MAPK-mediated complex traits Across S. cerevisiae

Signaling pathways enable cells to sense and respond to their environment. Many cellular signaling strategies are conserved from fungi to humans, yet their activity and phenotypic consequences can vary extensively among individuals within a species. A systematic assessment of the impact of naturally occurring genetic variation on signaling pathways remains to be conducted. In S. cerevisiae, both response and resistance to stressors that activate signaling pathways differ between diverse isolates. Here, we present a quantitative trait locus (QTL) mapping approach that enables us to identify genetic variants underlying such phenotypic differences across the genetic and phenotypic diversity of S. cerevisiae. Using a Round-robin cross between twelve diverse strains, we identified QTL that influence phenotypes critically dependent on MAPK signaling cascades. Genetic variants under these QTL fall within MAPK signaling networks themselves as well as other interconnected signaling pathways. Finally, we demonstrate how the mapping results from multiple strain background can be leveraged to narrow the search space of causal genetic variants.

Wild yeast strains differ in phenotypes that are controlled by highly conserved signaling pathways. Yet it remains unknown how naturally occurring genetic variants influence signaling pathways in yeast. We have developed an approach to facilitate the mapping of genetic variants that underlie these phenotypic differences in a set of wild strain. Our mapping approach requires minimal strain engineering and enables the rapid isolation of mapping populations from any strain background. In particular, we have mapped genetic variants in twelve highly diverse yeast strains. Further, we demonstrate how the mapping results from these twelve strains can be used jointly to narrow the number of genetic variants identified to a set of putative causal variants. We identify genetic variants in genes with various roles in cell signaling. Our results illustrate the interplay of different signaling pathways and which signaling genes are more likely to contain variants of large phenotypic impact.

Investigation of protein secretion and secretion stress in Ashbya gossypii

Background

Ashbya gossypii is a filamentous Saccharomycete used for the industrial production of riboflavin that has been recently explored as a host system for recombinant protein production. To gain insight into the protein secretory pathway of this biotechnologically relevant fungus, we undertook genome-wide analyses to explore its secretome and its transcriptional responses to protein secretion stress.

Results

A computational pipeline was used to predict the inventory of proteins putatively secreted by A. gossypii via the general secretory pathway. The proteins actually secreted by this fungus into the supernatants of submerged cultures in minimal and rich medium were mapped by two-dimensional gel electrophoresis, revealing that most of the A. gossypii secreted proteins have an isoelectric point between 4 and 6, and a molecular mass above 25 kDa. These analyses together indicated that 1-4% of A. gossypii proteins are likely to be secreted, of which less than 33% are putative hydrolases. Furthermore, transcriptomic analyses carried out in A. gossypii cells under recombinant protein secretion conditions and dithiothreitol-induced secretion stress unexpectedly revealed that a conventional unfolded protein response (UPR) was not activated in any of the conditions, as the expression levels of several well-known UPR target genes (e.g. IRE1, KAR2, HAC1 and PDI1 homologs) remained unaffected. However, several other genes involved in protein unfolding, endoplasmatic reticulum-associated degradation, proteolysis, vesicle trafficking, vacuolar protein sorting, secretion and mRNA degradation were up-regulated by dithiothreitol-induced secretion stress. Conversely, the transcription of several genes encoding secretory proteins, such as components of the glycosylation pathway, was severely repressed by dithiothreitol

Conclusions

This study provides the first insights into the secretion stress response of A. gossypii, as well as a basic understanding of its protein secretion potential, which is more similar to that of yeast than to that of other filamentous fungi. Contrary to what has been widely described for yeast and fungi, a conventional UPR was not observed in A. gossypii, but alternative protein quality control mechanisms enabled it to cope with secretion stress. These data will help provide strategies for improving heterologous protein secretion in A. gossypii.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1137) contains supplementary material, which is available to authorized users.

Patterns of variation during adaptation in functionally linked loci

An understanding of the distribution of natural patterns of genetic variation is relevant to such fundamental biological fields as evolution and development. One recent approach to understanding such patterns has been to focus on the constraints that may arise as a function of the network or pathway context in which genes are embedded. Despite theoretical expectations of higher evolutionary constraint for genes encoding upstream versus downstream enzymes in metabolic pathways, empirical results have varied. Here we combine two complementary models from population genetics and enzyme kinetics to explore genetic variation as a function of pathway position when selection acts on whole-pathway flux. We are able to qualitatively reproduce empirically observed patterns of polymorphism and divergence and suggest that expectations should vary depending on the evolutionary trajectory of a population. Upstream genes are initially more polymorphic and diverge faster after an environmental change, while we see the opposite trend as the population approaches its fitness optimum.

Hog1: 20 years of discovery and impact

The protein kinase Hog1 (high osmolarity glycerol 1) was discovered 20 years ago, being revealed as a central signaling mediator during osmoregulation in the budding yeast Saccharomyces cerevisiae. Homologs of Hog1 exist in all evaluated eukaryotic organisms, and this kinase plays a central role in cellular responses to external stresses and stimuli. Here, we highlight the mechanism by which cells sense changes in extracellular osmolarity, the method by which Hog1 regulates cellular adaptation, and the impacts of the Hog1 pathway upon cellular growth and morphology. Studies that have addressed these issues reveal the influence of the Hog1 signaling pathway on diverse cellular processes.

Effects of deletion of different PP2C protein phosphatase genes on stress responses in Saccharomyces cerevisiae

A key mechanism of signal transduction in eukaryotes is reversible protein phosphorylation, mediated through protein kinases and protein phosphatases (PPases). Modulation of signal transduction by this means regulates many biological processes. Saccharomyces cerevisiae has 40 PPases, including seven protein phosphatase 2C (PP2C PPase) genes (PTC1-PTC7). However, their precise functions remain poorly understood. To elucidate their cellular functions and to identify those that are redundant, we constructed 127 strains with deletions of all possible combinations of the seven PP2C PPase genes. All 127 disruptants were viable under nutrient-rich conditions, demonstrating that none of the combinations induced synthetic lethality under these conditions. However, several combinations exhibited novel phenotypes, e.g. the Δptc5Δptc7 double disruptant and the Δptc2Δptc3Δptc5Δptc7 quadruple disruptant exhibited low (13°C) and high (37°C) temperature-sensitive growth, respectively. Interestingly, the septuple disruptant Δptc1Δptc2Δptc3Δptc4Δptc5Δptc6Δptc7 showed an essentially normal growth phenotype at 37°C. The Δptc2Δptc3Δptc5Δptc7 quadruple disruptant was sensitive to LiCl (0.4 m). Two double disruptants, Δptc1Δptc2 and Δptc1Δptc4, displayed slow growth and Δptc1Δptc2Δptc4 could not grow on medium containing 1.5 m NaCl. The Δptc1Δptc6 double disruptant showed increased sensitivity to caffeine, congo red and calcofluor white compared to each single deletion. Our observations indicate that S. cerevisiae PP2C PPases have a shared and important role in responses to environmental stresses. These disruptants also provide a means for exploring the molecular mechanisms of redundant PTC gene functions under defined conditions.

Differential evolutionary constraints in the evolution of chemoreceptors: a murine and human case study

Chemoreception is among the most important sensory modalities in animals. Organisms use the ability to perceive chemical compounds in all major ecological activities. Recent studies have allowed the characterization of chemoreceptor gene families. These genes present strikingly high variability in copy numbers and pseudogenization degrees among different species, but the mechanisms underlying their evolution are not fully understood. We have analyzed the functional networks of these genes, their orthologs distribution, and performed phylogenetic analyses in order to investigate their evolutionary dynamics. We have modeled the chemosensory networks and compared the evolutionary constraints of their genes in Mus musculus, Homo sapiens, and Rattus norvegicus. We have observed significant differences regarding the constraints on the orthologous groups and network topologies of chemoreceptors and signal transduction machinery. Our findings suggest that chemosensory receptor genes are less constrained than their signal transducing machinery, resulting in greater receptor diversity and conservation of information processing pathways. More importantly, we have observed significant differences among the receptors themselves, suggesting that olfactory and bitter taste receptors are more conserved than vomeronasal receptors.

Evolutionary rate patterns of genes involved in the Drosophila Toll and Imd signaling pathway

BACKGROUND

To survive in a hostile environment, insects have evolved an innate immune system to defend against infection. Studies have shown that natural selection may drive the evolution of immune system-related proteins. Yet, how network architecture influences protein sequence evolution remains unclear. Here, we analyzed the molecular evolutionary patterns of genes in the Toll and Imd innate immune signaling pathways across six Drosophila genomes within the context of a functional network.

RESULTS

Based on published literature, we identified 50 genes that are directly involved in the Drosophila Toll and Imd signaling pathways. Of those genes, only two (Sphinx1 and Dnr1) exhibited signals of positive selection. There existed a negative correlation between the strength of purifying selection and gene position within the pathway; the downstream genes were more conserved, indicating that they were subjected to stronger evolutionary constraints. Interestingly, there was also a significantly negative correlation between the rate of protein evolution and the number of regulatory microRNAs, implying that genes regulated by more miRNAs experience stronger functional constraints and therefore evolve more slowly.

CONCLUSION

Taken together, our results suggested that both network architecture and miRNA regulation affect protein sequence evolution. These findings improve our understanding of the evolutionary patterns of genes involved in Drosophila innate immune pathways.

Comment on "The Molecular Evolutionary Patterns of the Insulin/FOXO Signaling Pathway"

Letter to the Editor on Wang M, Wang Q, Wang Z, Zhang X, Pan Y. The molecular evolutionary patterns of the insulin/FOXO signaling pathway. Evol Bioinform. 2013;9:1–16. doi: 10.4137/EBO.S10539.

The HOG signal transduction pathway in the halophilic fungus Wallemia ichthyophaga: identification and characterisation of MAP kinases WiHog1A and WiHog1B

The high-osmolarity glycerol (HOG) pathway is one of the several MAP kinase cascades in fungi. It is the main signal transduction system that is responsible for cellular stress responses, and has primarily been studied in the context of osmotic stress. In the present study, we provide the first insights into the HOG pathway of the obligatory halophilic basidiomycetous fungus Wallemia ichthyophaga, with the characterisation of its two Hog1-like kinases: WiHog1A and WiHog1B. These share high similarity to Hog1 kinase from Saccharomyces cerevisiae (ScHog1) at the level of amino-acid sequence. While WiHog1A could not optimally complement the function of ScHog1, WiHog1B was a fully functional Hog1-like kinase and could improve the halotolerance of the yeast, compared to the wild-type or the ScHog1-expressing hog1Δ strain. In W. ichthyophaga cells, Hog1 was constitutively phosphorylated under optimal osmotic conditions and dephosphorylated when the cells were challenged with hypo-osmolar or hyperosmolar stress. This pattern of phosphorylation kinetics is opposite to that of yeast. Transcriptional analysis of these two kinases in W. ichthyophaga shows that WiHOG1B is more responsive to changes in NaCl concentrations than WiHOG1A. Our identification and characterisation of these Hog1-like kinases from W. ichthyophaga confirm the existence of the HOG signalling pathway and its role in osmosensing in this halophilic fungus.

Molecular cloning and evolutionary analysis of the HOG-signaling pathway genes from Saccharomyces cerevisiae rice wine isolates

The high osmolarity glycerol (HOG) signaling pathway is crucial for yeast to cope with high osmolarity. Here, we showed that Saccharomyces cerevisiae rice wine isolates exhibited higher tolerance to osmotic stress, which was associated with the evolution of HOG pathway genes. Phylogenetic analysis of HOG genes revealed that Chinese rice wine strains were closely related to sake strains, indicating a common origin of rice wine strains. The DNA sequence diversity analysis showed that higher levels of polymorphism tended to accumulate on the osmosensor genes (MSB2 and SLN1), suggesting that most changes in a signaling transduction pathway were concentrated in the receptors. Moreover, the rapid evolution of osmosensors (Sln1/Msb2) and transcription factor (Msn4) might experience positive selection. Our results imply that the evolution of HOG pathway genes in S. cerevisiae rice wine strains is associated with their adaptation to high osmotic environments.

The evolution and origin of animal Toll-like receptor signaling pathway revealed by network-level molecular evolutionary analyses

Genes carry out their biological functions through pathways in complex networks consisting of many interacting molecules. Studies on the effect of network architecture on the evolution of individual proteins will provide valuable information for understanding the origin and evolution as well as functional conservation of signaling pathways. However, the relationship between the network architecture and the individual protein sequence evolution is yet little known. In current study, we carried out network-level molecular evolution analysis on TLR (Toll-like receptor ) signaling pathway, which plays an important role in innate immunity in insects and mammals, and we found that: 1) The selection constraint of genes was negatively correlated with its position along TLR signaling pathway; 2) all genes in TLR signaling pathway were highly conserved and underwent strong purifying selection; 3) the distribution of selective pressure along the pathway was driven by differential nonsynonymous substitution levels; 4) The TLR signaling pathway might present in a common ancestor of sponges and eumetazoa, and evolve via the TLR, IKK, IκB and NF-κB genes underwent duplication events as well as adaptor molecular enlargement, and gene structure and conservation motif of NF-κB genes shifted in their evolutionary history. Our results will improve our understanding on the evolutionary history of animal TLR signaling pathway as well as the relationship between the network architecture and the sequences evolution of individual protein.

The relationship between the hierarchical position of proteins in the human signal transduction network and their rate of evolution

Background

Proteins evolve at disparate rates, as a result of the action of different types and strengths of evolutionary forces. An open question in evolutionary biology is what factors are responsible for this variability. In general, proteins whose function has a great impact on organisms’ fitness are expected to evolve under stronger selective pressures. In biosynthetic pathways, upstream genes usually evolve under higher levels of selective constraint than those acting at the downstream part, as a result of their higher hierarchical position. Similar observations have been made in transcriptional regulatory networks, whose upstream elements appear to be more essential and subject to selection. Less well understood is, however, how selective pressures distribute along signal transduction pathways.

Results

Here, I combine comparative genomics and directed protein interaction data to study the distribution of evolutionary forces across the human signal transduction network. Surprisingly, no evidence was found for higher levels of selective constraint at the upstream network genes (those occupying more hierarchical positions). On the contrary, purifying selection was found to act more strongly on genes acting at the downstream part of the network, which seems to be due to downstream genes being more highly and broadly expressed, performing certain functions and, in particular, encoding proteins that are more highly connected in the protein–protein interaction network. When the effect of these confounding factors is discounted, upstream and downstream genes evolve at similar rates. The trends found in the overall signaling network are exemplified by analysis of the distribution of purifying selection along the mammalian Ras signaling pathway, showing that upstream and downstream genes evolve at similar rates.

Conclusions

These results indicate that the upstream/downstream position of proteins in the signal transduction network has, in general, no direct effect on their rates of evolution, suggesting that upstream and downstream genes are similarly important for the function of the network. This implies that natural selection differently distributes across signal transduction networks and across biosynthetic and transcriptional regulatory networks, which might reflect fundamental differences in their function and organization.

Axl2 integrates polarity establishment, maintenance, and environmental stress response in the filamentous fungus Ashbya gossypii

In budding yeast, new sites of polarity are chosen with each cell cycle and polarization is transient. In filamentous fungi, sites of polarity persist for extended periods of growth and new polarity sites can be established while existing sites are maintained. How the polarity establishment machinery functions in these distinct growth forms found in fungi is still not well understood. We have examined the function of Axl2, a transmembrane bud site selection protein discovered in Saccharomyces cerevisiae, in the filamentous fungus Ashbya gossypii. A. gossypii does not divide by budding and instead exhibits persistent highly polarized growth, and multiple axes of polarity coexist in one cell. A. gossypii axl2Δ (Agaxl2Δ) cells have wavy hyphae, bulbous tips, and a high frequency of branch initiations that fail to elongate, indicative of a polarity maintenance defect. Mutant colonies also have significantly lower radial growth and hyphal tip elongation speeds than wild-type colonies, and Agaxl2Δ hyphae have depolarized actin patches. Consistent with a function in polarity, AgAxl2 localizes to hyphal tips, branches, and septin rings. Unlike S. cerevisiae Axl2, AgAxl2 contains a Mid2 homology domain and may function to sense or respond to environmental stress. In support of this idea, hyphae lacking AgAxl2 also display hypersensitivity to heat, osmotic, and cell wall stresses. Axl2 serves to integrate polarity establishment, polarity maintenance, and environmental stress response for optimal polarized growth in A. gossypii.

Expression level drives the pattern of selective constraints along the insulin/Tor signal transduction pathway in Caenorhabditis

Genes do not act in isolation but perform their biological functions within genetic pathways that are connected in larger networks. Investigation of nucleotide variation within genetic pathways and networks has shown that topology can affect the rate of protein evolution; however, it remains unclear whether a same pattern of nucleotide variation is expected within functionally similar networks and whether it may be due to similar or different biological mechanisms. We address these questions by investigating nucleotide variation in the context of the structure of the insulin/Tor-signaling pathway in Caenorhabditis, which is well characterized and is functionally conserved across phylogeny. In Drosophila and vertebrates, the rate of protein evolution is negatively correlated with the position of a gene within the insulin/Tor pathway. Similarly, we find that in Caenorhabditis, the rate of amino acid replacement is lower for downstream genes. However, in Caenorhabditis, the rate of synonymous substitution is also strongly affected by the position of a gene in the pathway, and we show that the distribution of selective pressure along the pathway is driven by differential expression level. A full understanding of the effect of pathway structure on selective constraints is therefore likely to require inclusion of specific biological function into more general network models.

Evolution of apoptosis-like programmed cell death in unicellular protozoan parasites

Apoptosis-like programmed cell death (PCD) has recently been described in multiple taxa of unicellular protists, including the protozoan parasites Plasmodium, Trypanosoma and Leishmania. Apoptosis-like PCD in protozoan parasites shares a number of morphological features with programmed cell death in multicellular organisms. However, both the evolutionary explanations and mechanisms involved in parasite PCD are poorly understood. Explaining why unicellular organisms appear to undergo 'suicide' is a challenge for evolutionary biology and uncovering death executors and pathways is a challenge for molecular and cell biology. Bioinformatics has the potential to integrate these approaches by revealing homologies in the PCD machinery of diverse taxa and evaluating their evolutionary trajectories. As the molecular mechanisms of apoptosis in model organisms are well characterised, and recent data suggest similar mechanisms operate in protozoan parasites, key questions can now be addressed. These questions include: which elements of apoptosis machinery appear to be shared between protozoan parasites and multicellular taxa and, have these mechanisms arisen through convergent or divergent evolution? We use bioinformatics to address these questions and our analyses suggest that apoptosis mechanisms in protozoan parasites and other taxa have diverged during their evolution, that some apoptosis factors are shared across taxa whilst others have been replaced by proteins with similar biochemical activities.

Identification and characterization of putative osmosensors, HwSho1A and HwSho1B, from the extremely halotolerant black yeast Hortaea werneckii

In Saccharomyces cerevisiae, the Sho1 protein is one of two potential osmosensors that can activate the kinase cascade of the HOG pathway in response to increased extracellular osmolarity. Two novel SHO1-like genes, HwSHO1A and HwSHO1B, have been cloned from the saltern-inhabiting, extremely halotolerant black yeast Hortaea werneckii. The HwSho1 protein isoforms are 93.8% identical in their amino-acid sequences, and have a conserved SH3 domain. When the HwSHO1 genes were transferred into S. cerevisae cells lacking the SHO1 gene, both of the HwSho1 isoforms fully complemented the function of the native S. cerevisiae Sho1 protein. Through microscopic and biochemical validation, we demonstrate that in S. cerevisiae, both of the HwSho1 proteins have characteristic subcellular localizations similar to the S. cerevisiae Sho1 protein, and they can both activate the HOG pathway under conditions of osmotic stress. To a lower extent, crosstalk to the mating pathway expressing HwSho1 proteins is conserved in the PBS2 deleted S. cerevisiae strain. These data show that the HwSho1 proteins from H. werneckii are true functional homologs of the Sho1 protein of S. cerevisiae.