Estimating genomic coexpression networks using first-order conditional independence

We describe a computationally efficient statistical framework for estimating networks of coexpressed genes. This framework exploits first-order conditional independence relationships among gene-expression measurements to estimate patterns of association. We use this approach to estimate a coexpression network from microarray gene-expression measurements from Saccharomyces cerevisiae. We demonstrate the biological utility of this approach by showing that a large number of metabolic pathways are coherently represented in the estimated network. We describe a complementary unsupervised graph search algorithm for discovering locally distinct subgraphs of a large weighted graph. We apply this algorithm to our coexpression network model and show that subgraphs found using this approach correspond to particular biological processes or contain representatives of distinct gene families

Modeling mutant phenotypes and oscillatory dynamics in the Saccharomyces cerevisiae cAMP-PKA pathway

Background The cyclic AMP-Protein Kinase A (cAMP-PKA) pathway is an evolutionarily conserved signal transduction mechanism that regulates cellular growth and differentiation in animals and fungi. We present a mathematical model that recapitulates the short-term and long-term dynamics of this pathway in the budding yeast, Saccharomyces cerevisiae. Our model is aimed at recapitulating the dynamics of cAMP signaling for wild-type cells as well as single (pde1Δ and pde2Δ) and double (pde1Δpde2Δ) phosphodiesterase mutants. Results Our model focuses on PKA-mediated negative feedback on the activity of phosphodiesterases and the Ras branch of the cAMP-PKA pathway. We show that both of these types of negative feedback are required to reproduce the wild-type signaling behavior that occurs on both short and long time scales, as well as the the observed responses of phosphodiesterase mutants. A novel feature of our model is that, for a wide range of parameters, it predicts that intracellular cAMP concentrations should exhibit decaying oscillatory dynamics in their approach to steady state following glucose stimulation. Experimental measurements of cAMP levels in two genetic backgrounds of S. cerevisiae confirmed the presence of decaying cAMP oscillations as predicted by the model. Conclusions Our model of the cAMP-PKA pathway provides new insights into how yeast respond to alterations in their nutrient environment. Because the model has both predictive and explanatory power it will serve as a foundation for future mathematical and experimental studies of this important signaling network

Elucidation of Directionality for Co-Expressed Genes: Predicting Intra-Operon Termination Sites

We present a novel framework for inferring regulatory and sequence-level information from gene co-expression networks. The key idea of our methodology is the systematic integration of network inference and network topological analysis approaches for uncovering biological insights. We determine the gene co-expression network of Bacillus subtilis using Affymetrix GeneChip time series data and show how the inferred network topology can be linked to sequence-level information hard-wired in the organism's genome. We propose a systematic way for determining the correlation threshold at which two genes are assessed to be co-expressed by using the clustering coefficient and we expand the scope of the gene co-expression network by proposing the slope ratio metric as a means for incorporating directionality on the edges. We show through specific examples for B. subtilis that by incorporating expression level information in addition to the temporal expression patterns, we can uncover sequence-level biological insights. In particular, we are able to identify a number of cases where (i) the co-expressed genes are part of a single transcriptional unit or operon and (ii) the inferred directionality arises due to the presence of intra-operon transcription termination sites.Comment: 7 pages, 8 figures, accepted in Bioinformatic

On excitable beta-skeletons

A beta-skeleton is a planar proximity undirected graph of an Euclidean point set where nodes are connected by an edge if their lune-based neighborhood contains no other points of the given set. Parameter $\beta$ determines size and shape of the nodes' neighborhoods. In an excitable beta-skeleton every node takes three states --- resting, excited and refractory, and updates its state in discrete time depending on states of its neighbors. We design families of beta-skeletons with absolute and relative thresholds of excitability and demonstrate that several distinct classes of space-time excitation dynamics can be selected using beta. The classes include spiral and target waves of excitation, branching domains of excitation and oscillating localizations

Amoeba predation of <i>Cryptococcus</i>:A quantitative and population genomic evaluation of the accidental pathogen hypothesis

The “Amoeboid Predator-Fungal Animal Virulence Hypothesis” posits that interactions with environmental phagocytes shape the evolution of virulence traits in fungal pathogens. In this hypothesis, selection to avoid predation by amoeba inadvertently selects for traits that contribute to fungal escape from phagocytic immune cells. Here, we investigate this hypothesis in the human fungal pathogens Cryptococcus neoformans and Cryptococcus deneoformans. Applying quantitative trait locus (QTL) mapping and comparative genomics, we discovered a cross-species QTL region that is responsible for variation in resistance to amoeba predation. In C. neoformans, this same QTL was found to have pleiotropic effects on melanization, an established virulence factor. Through fine mapping and population genomic comparisons, we identified the gene encoding the transcription factor Bzp4 that underlies this pleiotropic QTL and we show that decreased expression of this gene reduces melanization and increases susceptibility to amoeba predation. Despite the joint effects of BZP4 on amoeba resistance and melanin production, we find no relationship between BZP4 genotype and escape from macrophages or virulence in murine models of disease. Our findings provide new perspectives on how microbial ecology shapes the genetic architecture of fungal virulence, and suggests the need for more nuanced models for the evolution of pathogenesis that account for the complexities of both microbe-microbe and microbe-host interactions

The geographic distribution of saccharomyces cerevisiae isolates within three Italian neighboring winemaking regions reveals strong differences in yeast abundance, genetic diversity and industrial strain dissemination

In recent years the interest for natural fermentations has been re-evaluated in terms of increasing the wine terroir and managing more sustainable winemaking practices. Therefore, the level of yeast genetic variability and the abundance of Saccharomyces cerevisiae native populations in vineyard are becoming more and more crucial at both ecological and technological level. Among the factors that can influence the strain diversity, the commercial starter release that accidentally occur in the environment around the winery, has to be considered. In this study we led a wide scale investigation of S. cerevisiae genetic diversity and population structure in the vineyards of three neighboring winemaking regions of Protected Appellation of Origin, in North-East of Italy. Combining mtDNA RFLP and microsatellite markers analyses we evaluated 634 grape samples collected over 3 years. We could detect major differences in the presence of S. cerevisiae yeasts, according to the winemaking region. The population structures revealed specificities of yeast microbiota at vineyard scale, with a relative Appellation of Origin area homogeneity, and transition zones suggesting a geographic differentiation. Surprisingly, we found a widespread industrial yeast dissemination that was very high in the areas where the native yeast abundance was low. Although geographical distance is a key element involved in strain distribution, the high presence of industrial strains in vineyard reduced the differences between populations. This finding indicates that industrial yeast diffusion it is a real emergency and their presence strongly interferes with the natural yeast microbiota

The Statistics of Bulk Segregant Analysis Using Next Generation Sequencing

We describe a statistical framework for QTL mapping using bulk segregant analysis (BSA) based on high throughput, short-read sequencing. Our proposed approach is based on a smoothed version of the standard statistic, and takes into account variation in allele frequency estimates due to sampling of segregants to form bulks as well as variation introduced during the sequencing of bulks. Using simulation, we explore the impact of key experimental variables such as bulk size and sequencing coverage on the ability to detect QTLs. Counterintuitively, we find that relatively large bulks maximize the power to detect QTLs even though this implies weaker selection and less extreme allele frequency differences. Our simulation studies suggest that with large bulks and sufficient sequencing depth, the methods we propose can be used to detect even weak effect QTLs and we demonstrate the utility of this framework by application to a BSA experiment in the budding yeast Saccharomyces cerevisiae

A computational pipeline to discover highly phylogenetically informative genes in sequenced genomes: application to Saccharomyces cerevisiae natural strains

The quest for genes representing genetic relationships of strains or individuals within populations and their evolutionary history is acquiring a novel dimension of complexity with the advancement of next-generation sequencing (NGS) technologies. In fact, sequencing an entire genome uncovers genetic variation in coding and non-coding regions and offers the possibility of studying Saccharomyces cerevisiae populations at the strain level. Nevertheless, the disadvantageous cost-benefit ratio (the amount of details disclosed by NGS against the time-expensive and expertise-demanding data assembly process) still precludes the application of these techniques to the routinely assignment of yeast strains, making the selection of the most reliable molecular markers greatly desirable. In this work we propose an original computational approach to discover genes that can be used as a descriptor of the population structure. We found 13 genes whose variability can be used to recapitulate the phylogeny obtained from genome-wide sequences. The same approach that we prove to be successful in yeasts can be generalized to any other population of individuals given the availability of high-quality genomic sequences and of a clear population structure to be targeted

The search for Pleiades in trait constellations: functional integration and phenotypic selection in the complex flowers of Morrenia brachystephana (Apocynaceae)

Pollinator‐mediated natural selection on single traits, such as corolla tube or spur length, has been well documented. However, flower phenotypes are usually complex, and selection is expected to act on several traits that functionally interact rather than on a single isolated trait. Despite the fact that selection on complex phenotypes is expectedly widespread, multivariate selection modelling on such phenotypes still remains under‐explored in plants. Species of the subfamily Asclepiadoideae (Apocynaceae) provide an opportunity to study such complex flower contrivances integrated by fine‐scaled organs from disparate developmental origin. We studied the correlation structure among linear floral traits (i) by testing a priori morphological, functional or developmental hypotheses among traits and (ii) by exploring the organization of flower covariation, considering alternative expectations of modular organization or whole flower integration through conditional dependence analysis (CDA) and integration matrices. The phenotypic selection approach was applied to determine whether floral traits involved in the functioning of the pollination mechanism were affected by natural selection. Floral integration was low, suggesting that flowers are organized in more than just one correlation pleiad; our hypothetical functional correlation matrix was significantly correlated with the empirical matrix, and the CDA revealed three putative modules. Analyses of phenotypic selection showed significant linear and correlational gradients, lending support to expectations of functional interactions between floral traits. Significant correlational selection gradients found involved traits of different floral whorls, providing evidence for the existence of functional integration across developmental domains.publishedVersionFil: Baranzelli, Matías Cristian. Universidad Nacional de Córdoba. Facultad de Ciencias Exactas, Físicas y Naturales; Argentina.Fil: Baranzelli, Matías Cristian. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Investigaciones Biológicas y Tecnológicas; Argentina.Fil: Sérsic, A. N. Universidad Nacional de Córdoba. Facultad de Ciencias Exactas, Físicas y Naturales; Argentina.Fil: Sérsic, A. N. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Investigaciones Biológicas y Tecnológicas; Argentina.Fil: Cocucci, A. A. Universidad Nacional de Córdoba. Facultad de Ciencias Exactas, Físicas y Naturales; Argentina.Fil: Cocucci, A. A. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Investigaciones Biológicas y Tecnológicas; Argentina

HAMSTER: visualizing microarray experiments as a set of minimum spanning trees

<p>Abstract</p> <p>Background</p> <p>Visualization tools allow researchers to obtain a global view of the interrelationships between the probes or experiments of a gene expression (<it>e.g. microarray</it>) data set. Some existing methods include hierarchical clustering and k-means. In recent years, others have proposed applying minimum spanning trees (MST) for microarray clustering. Although MST-based clustering is formally equivalent to the dendrograms produced by hierarchical clustering under certain conditions; visually they can be quite different.</p> <p>Methods</p> <p>HAMSTER (Helpful Abstraction using Minimum Spanning Trees for Expression Relations) is an open source system for generating a <b>set </b>of MSTs from the experiments of a microarray data set. While previous works have generated a single MST from a data set for data clustering, we recursively merge experiments and repeat this process to obtain a set of MSTs for data visualization. Depending on the parameters chosen, each tree is analogous to a snapshot of one step of the hierarchical clustering process. We scored and ranked these trees using one of three proposed schemes. HAMSTER is implemented in C++ and makes use of Graphviz for laying out each MST.</p> <p>Results</p> <p>We report on the running time of HAMSTER and demonstrate using data sets from the NCBI Gene Expression Omnibus (GEO) that the images created by HAMSTER offer insights that differ from the dendrograms of hierarchical clustering. In addition to the C++ program which is available as open source, we also provided a web-based version (HAMSTER⁺) which allows users to apply our system through a web browser without any computer programming knowledge.</p> <p>Conclusion</p> <p>Researchers may find it helpful to include HAMSTER in their microarray analysis workflow as it can offer insights that differ from hierarchical clustering. We believe that HAMSTER would be useful for certain types of gradient data sets (e.g time-series data) and data that indicate relationships between cells/tissues. Both the source and the web server variant of HAMSTER are available from <url>http://hamster.cbrc.jp/</url>.</p