Bulk Segregant Analysis Using Single Nucleotide Polymorphism Microarrays

Bulk segregant analysis (BSA) using microarrays, and extreme array mapping (XAM) have recently been used to rapidly identify genomic regions associated with phenotypes in multiple species. These experiments, however, require the identification of single feature polymorphisms (SFP) between the cross parents for each new combination of genotypes, which raises the cost of experiments. The availability of the genomic polymorphism data in Arabidopsis thaliana, coupled with the efficient designs of Single Nucleotide Polymorphism (SNP) genotyping arrays removes the requirement for SFP detection and lowers the per array cost, thereby lowering the overall cost per experiment. To demonstrate that these approaches would be functional on SNP arrays and determine confidence intervals, we analyzed hybridizations of natural accessions to the Arabidopsis ATSNPTILE array and simulated BSA or XAM given a variety of gene models, populations, and bulk selection parameters. Our results show a striking degree of correlation between the genotyping output of both methods, which suggests that the benefit of SFP genotyping in context of BSA can be had with the cheaper, more efficient SNP arrays. As a final proof of concept, we hybridized the DNA from bulks of an F2 mapping population of a Sulfur and Selenium ionomics mutant to both the Arabidopsis ATTILE1R and ATSNPTILE arrays, which produced almost identical results. We have produced R scripts that prompt the user for the required parameters and perform the BSA analysis using the ATSNPTILE1 array and have provided them as supplemental data files

Interaction Networks in Yeast Define and Enumerate the Signaling Steps of the Vertebrate Aryl Hydrocarbon Receptor

The aryl hydrocarbon receptor (AHR) is a vertebrate protein that mediates the toxic and adaptive responses to dioxins and related environmental pollutants. In an effort to better understand the details of this signal transduction pathway, we employed the yeast S. cerevisiae as a model system. Through the use of arrayed yeast strains harboring ordered deletions of open reading frames, we determined that 54 out of the 4,507 yeast genes examined significantly influence AHR signal transduction. In an effort to describe the relationship between these modifying genes, we constructed a network map based upon their known protein and genetic interactions. Monte Carlo simulations demonstrated that this network represented a description of AHR signaling that was distinct from those generated by random chance. The network map was then explored with a number of computational and experimental annotations. These analyses revealed that the AHR signaling pathway is defined by at least five distinct signaling steps that are regulated by functional modules of interacting modifiers. These modules can be described as mediating receptor folding, nuclear translocation, transcriptional activation, receptor level, and a previously undescribed nuclear step related to the receptor's Per–Arnt–Sim domain

Discovery and Expansion of Gene Modules by Seeking Isolated Groups in a Random Graph Process

BACKGROUND: A central problem in systems biology research is the identification and extension of biological modules-groups of genes or proteins participating in a common cellular process or physical complex. As a result, there is a persistent need for practical, principled methods to infer the modular organization of genes from genome-scale data. RESULTS: We introduce a novel approach for the identification of modules based on the persistence of isolated gene groups within an evolving graph process. First, the underlying genomic data is summarized in the form of ranked gene-gene relationships, thereby accommodating studies that quantify the relevant biological relationship directly or indirectly. Then, the observed gene-gene relationship ranks are viewed as the outcome of a random graph process and candidate modules are given by the identifiable subgraphs that arise during this process. An isolation index is computed for each module, which quantifies the statistical significance of its survival time. CONCLUSIONS: The Miso (module isolation) method predicts gene modules from genomic data and the associated isolation index provides a module-specific measure of confidence. Improving on existing alternative, such as graph clustering and the global pruning of dendrograms, this index offers two intuitively appealing features: (1) the score is module-specific; and (2) different choices of threshold correlate logically with the resulting performance, i.e. a stringent cutoff yields high quality predictions, but low sensitivity. Through the analysis of yeast phenotype data, the Miso method is shown to outperform existing alternatives, in terms of the specificity and sensitivity of its predictions

How to identify essential genes from molecular networks?

<p>Abstract</p> <p>Background</p> <p>The prediction of essential genes from molecular networks is a way to test the understanding of essentiality in the context of what is known about the network. However, the current knowledge on molecular network structures is incomplete yet, and consequently the strategies aimed to predict essential genes are prone to uncertain predictions. We propose that simultaneously evaluating different network structures and different algorithms representing gene essentiality (centrality measures) may identify essential genes in networks in a reliable fashion.</p> <p>Results</p> <p>By simultaneously analyzing 16 different centrality measures on 18 different reconstructed metabolic networks for <it>Saccharomyces cerevisiae</it>, we show that no single centrality measure identifies essential genes from these networks in a statistically significant way; however, the combination of at least 2 centrality measures achieves a reliable prediction of most but not all of the essential genes. No improvement is achieved in the prediction of essential genes when 3 or 4 centrality measures were combined.</p> <p>Conclusion</p> <p>The method reported here describes a reliable procedure to predict essential genes from molecular networks. Our results show that essential genes may be predicted only by combining centrality measures, revealing the complex nature of the function of essential genes.</p

Mapping Genetically Compensatory Pathways from Synthetic Lethal Interactions in Yeast

Background: Synthetic lethal genetic interaction analysis has been successfully applied to predicting the functions of genes and their pathway identities. In the context of synthetic lethal interaction data alone, the global similarity of synthetic lethal interaction patterns between two genes is used to predict gene function. With physical interaction data, such as proteinprotein interactions, the enrichment of physical interactions within subsets of genes and the enrichment of synthetic lethal interactions between those subsets of genes are used as an indication of compensatory pathways. Result: In this paper, we propose a method of mapping genetically compensatory pathways from synthetic lethal interactions. Our method is designed to discover pairs of gene-sets in which synthetic lethal interactions are depleted among the genes in an individual set and where such gene-set pairs are connected by many synthetic lethal interactions. By its nature, our method could select compensatory pathway pairs that buffer the deleterious effect of the failure of either one, without the need of physical interaction data. By focusing on compensatory pathway pairs where genes in each individual pathway have a highly homogenous cellular function, we show that many cellular functions have genetically compensatory properties. Conclusion: We conclude that synthetic lethal interaction data are a powerful source to map genetically compensatory pathways, especially in systems lacking physical interaction information, and that the cellular function network contain

The role of the protein kinase A pathway in the response to alkaline pH stress in yeast

Exposure of Saccharomyces cerevisiae to alkaline pH provokes a stress condition that generates a compensatory reaction. In the present study we examined a possible role for the PKA (protein kinase A) pathway in this response. Phenotypic analysis revealed that mutations that activate the PKA pathway (ira1 ira2, bcy1) tend to cause sensitivity to alkaline pH, whereas its deactivation enhances tolerance to this stress. We observed that alkalinization causes a transient decrease in cAMP, the main regulator of the pathway. Alkaline pH causes rapid nuclear localization of the PKA-regulated Msn2 transcription factor which, together with Msn4, mediates a general stress response by binding with STRE (stress response element) sequences in many promoters. Consequently, a synthetic STRE–LacZ reporter shows a rapid induction in response to alkaline stress. A msn2 msn4 mutant is sensitive to alkaline pH, and transcriptomic analysis reveals that after 10 min of alkaline stress, the expression of many induced genes (47%) depends, at least in part, on the presence of Msn2 and Msn4. Taken together, these results demonstrate that inhibition of the PKA pathway by alkaline pH represents a substantial part of the adaptive response to this kind of stress and that this response involves Msn2/Msn4-mediated genome expression remodelling. However, the relevance of attenuation of PKA in high pH tolerance is probably not restricted to regulation of Msn2 function

High-Resolution Genotyping via Whole Genome Hybridizations to Microarrays Containing Long Oligonucleotide Probes

To date, microarray-based genotyping of large, complex plant genomes has been complicated by the need to perform genome complexity reduction to obtain sufficiently strong hybridization signals. Genome complexity reduction techniques are, however, tedious and can introduce unwanted variables into genotyping assays. Here, we report a microarray-based genotyping technology for complex genomes (such as the 2.3 GB maize genome) that does not require genome complexity reduction prior to hybridization. Approximately 200,000 long oligonucleotide probes were identified as being polymorphic between the inbred parents of a mapping population and used to genotype two recombinant inbred lines. While multiple hybridization replicates provided ∼97% accuracy, even a single replicate provided ∼95% accuracy. Genotyping accuracy was further increased to >99% by utilizing information from adjacent probes. This microarray-based method provides a simple, high-density genotyping approach for large, complex genomes

A Two-Hybrid Assay to Study Protein Interactions within the Secretory Pathway

Interactions of transcriptional activators are difficult to study using transcription-based two-hybrid assays due to potent activation resulting in false positives. Here we report the development of the Golgi two-hybrid (G2H), a method that interrogates protein interactions within the Golgi, where transcriptional activators can be assayed with negligible background. The G2H relies on cell surface glycosylation to report extracellularly on protein-protein interactions occurring within the secretory pathway. In the G2H, protein pairs are fused to modular domains of the reporter glycosyltransferase, Och1p, and proper cell wall formation due to Och1p activity is observed only when a pair of proteins interacts. Cells containing interacting protein pairs are identified by selectable phenotypes associated with Och1p activity and proper cell wall formation: cells that have interacting proteins grow under selective conditions and display weak wheat germ agglutinin (WGA) binding by flow cytometry, whereas cells that lack interacting proteins display stunted growth and strong WGA binding. Using this assay, we detected the interaction between transcription factor MyoD and its binding partner Id2. Interfering mutations along the MyoD:Id2 interaction interface ablated signal in the G2H assay. Furthermore, we used the G2H to detect interactions of the activation domain of Gal4p with a variety of binding partners. Finally, selective conditions were used to enrich for cells encoding interacting partners. The G2H detects protein-protein interactions that cannot be identified via traditional two-hybrid methods and should be broadly useful for probing previously inaccessible subsets of the interactome, including transcriptional activators and proteins that traffic through the secretory pathway