On the choice and number of microarrays for transcriptional regulatory network inference

BACKGROUND

Transcriptional regulatory network inference (TRNI) from large compendia of DNA microarrays has become a fundamental approach for discovering transcription factor (TF)-gene interactions at the genome-wide level. In correlation-based TRNI, network edges can in principle be evaluated using standard statistical tests. However, while such tests nominally assume independent microarray experiments, we expect dependency between the experiments in microarray compendia, due to both project-specific factors (e.g., microarray preparation, environmental effects) in the multi-project compendium setting and effective dependency induced by gene-gene correlations. Herein, we characterize the nature of dependency in an Escherichia coli microarray compendium and explore its consequences on the problem of determining which and how many arrays to use in correlation-based TRNI.

RESULTS

We present evidence of substantial effective dependency among microarrays in this compendium, and characterize that dependency with respect to experimental condition factors. We then introduce a measure neff of the effective number of experiments in a compendium, and find that corresponding to the dependency observed in this particular compendium there is a huge reduction in effective sample size i.e., neff = 14.7 versus n = 376. Furthermore, we found that the neff of select subsets of experiments actually exceeded neff of the full compendium, suggesting that the adage 'less is more' applies here. Consistent with this latter result, we observed improved performance in TRNI using subsets of the data compared to results using the full compendium. We identified experimental condition factors that trend with changes in TRNI performance and neff , including growth phase and media type. Finally, using the set of known E. coli genetic regulatory interactions from RegulonDB, we demonstrated that false discovery rates (FDR) derived from neff -adjusted p-values were well-matched to FDR based on the RegulonDB truth set.

CONCLUSIONS

These results support utilization of neff as a potent descriptor of microarray compendia. In addition, they highlight a straightforward correlation-based method for TRNI with demonstrated meaningful statistical testing for significant edges, readily applicable to compendia from any species, even when a truth set is not available. This work facilitates a more refined approach to construction and utilization of mRNA expression compendia in TRNI.

Predicting gene targets of perturbations via network-based filtering of mRNA expression compendia

MOTIVATION

DNA microarrays are routinely applied to study diseased or drug-treated cell populations. A critical challenge is distinguishing the genes directly affected by these perturbations from the hundreds of genes that are indirectly affected. Here, we developed a sparse simultaneous equation model (SSEM) of mRNA expression data and applied Lasso regression to estimate the model parameters, thus constructing a network model of gene interaction effects. This inferred network model was then used to filter data from a given experimental condition of interest and predict the genes directly targeted by that perturbation.

RESULTS

Our proposed SSEM-Lasso method demonstrated substantial improvement in sensitivity compared with other tested methods for predicting the targets of perturbations in both simulated datasets and microarray compendia. In simulated data, for two different network types, and over a wide range of signal-to-noise ratios, our algorithm demonstrated a 167% increase in sensitivity on average for the top 100 ranked genes, compared with the next best method. Our method also performed well in identifying targets of genetic perturbations in microarray compendia, with up to a 24% improvement in sensitivity on average for the top 100 ranked genes. The overall performance of our network-filtering method shows promise for identifying the direct targets of genetic dysregulation in cancer and disease from expression profiles.

AVAILABILITY

Microarray data are available at the Many Microbe Microarrays Database (M3D, http://m3d.bu.edu). Algorithm scripts are available at the Gardner Lab website (http://gardnerlab.bu.edu/SSEMLasso).

Towards environmental systems biology of Shewanella

Bacteria of the genus Shewanella are known for their versatile electron-accepting capacities, which allow them to couple the decomposition of organic matter to the reduction of the various terminal electron acceptors that they encounter in their stratified environments. Owing to their diverse metabolic capabilities, shewanellae are important for carbon cycling and have considerable potential for the remediation of contaminated environments and use in microbial fuel cells. Systems-level analysis of the model species Shewanella oneidensis MR-1 and other members of this genus has provided new insights into the signal-transduction proteins, regulators, and metabolic and respiratory subsystems that govern the remarkable versatility of the shewanellae.

Predicting sigma28 promoters in eleven Shewanella genomes

An iterative position-specific score matrix (PSSM)-based approach was used to predict sigma(28) promoters in 11 Shewanella genomes. The Shewanella Correlation Browser was used to distinguish true-positive predictions from false-positive predictions in Shewanella oneidensis MR-1 by generating a sigma(28)-regulated transcriptional network from transcriptional profiling data. This dual-pronged approach identified several genes that have sigma(28) promoters and that may be involved with motility or chemotaxis in Shewanella.

The mode-of-action by network identification (MNI) algorithm: a network biology approach for molecular target identification

This protocol details the use of the mode-of-action by network identification (MNI) algorithm to identify the gene targets of a drug treatment based on gene-expression data. Investigators might also use the MNI algorithm to identify the gene mediators of a disease or the physiological state of cells and tissues. The MNI algorithm uses a training data set of hundreds of expression profiles to construct a statistical model of gene-regulatory networks in a cell or tissue. The model describes combinatorial influences of genes on one another. The algorithm then uses the model to filter the expression profile of a particular experimental treatment and thereby distinguish the molecular targets or mediators of the treatment response from hundreds of additional genes that also exhibit expression changes. It takes approximately 1 h per run, although run time is significantly affected by the size of the genome and data set.

Many Microbe Microarrays Database: uniformly normalized Affymetrix compendia with structured experimental metadata

Many Microbe Microarrays Database (M^3D) is designed to facilitate the analysis and visualization of expression data in compendia compiled from multiple laboratories. M^3D contains over a thousand Affymetrix microarrays for Escherichia coli, Saccharomyces cerevisiae and Shewanella oneidensis. The expression data is uniformly normalized to make the data generated by different laboratories and researchers more comparable. To facilitate computational analyses, M^3D provides raw data (CEL file) and normalized data downloads of each compendium. In addition, web-based construction, visualization and download of custom datasets are provided to facilitate efficient interrogation of the compendium for more focused analyses. The experimental condition metadata in M^3D is human curated with each chemical and growth attribute stored as a structured and computable set of experimental features with consistent naming conventions and units. All versions of the normalized compendia constructed for each species are maintained and accessible in perpetuity to facilitate the future interpretation and comparison of results published on M^3D data. M^3D is accessible at http://m3d.bu.edu/.

Large-scale mapping and validation of Escherichia coli transcriptional regulation from a compendium of expression profiles

Machine learning approaches offer the potential to systematically identify transcriptional regulatory interactions from a compendium of microarray expression profiles. However, experimental validation of the performance of these methods at the genome scale has remained elusive. Here we assess the global performance of four existing classes of inference algorithms using 445 Escherichia coli Affymetrix arrays and 3,216 known E. coli regulatory interactions from RegulonDB. We also developed and applied the context likelihood of relatedness (CLR) algorithm, a novel extension of the relevance networks class of algorithms. CLR demonstrates an average precision gain of 36% relative to the next-best performing algorithm. At a 60% true positive rate, CLR identifies 1,079 regulatory interactions, of which 338 were in the previously known network and 741 were novel predictions. We tested the predicted interactions for three transcription factors with chromatin immunoprecipitation, confirming 21 novel interactions and verifying our RegulonDB-based performance estimates. CLR also identified a regulatory link providing central metabolic control of iron transport, which we confirmed with real-time quantitative PCR. The compendium of expression data compiled in this study, coupled with RegulonDB, provides a valuable model system for further improvement of network inference algorithms using experimental data.

Lightweight genome viewer: portable software for browsing genomics data in its chromosomal context

Background

Lightweight genome viewer (lwgv) is a web-based tool for visualization of sequence annotations in their chromosomal context. It performs most of the functions of larger genome browsers, while relying on standard flat-file formats and bypassing the database needs of most visualization tools. Visualization as an aide to discovery requires display of novel data in conjunction with static annotations in their chromosomal context. With database-based systems, displaying dynamic results requires temporary tables that need to be tracked for removal.

Results

lwgv simplifies the visualization of user-generated results on a local computer. The dynamic results of these analyses are written to transient files, which can import static content from a more permanent file. lwgv is currently used in many different applications, from whole genome browsers to single-gene RNAi design visualization, demonstrating its applicability in a large variety of contexts and scales.

Conclusion

lwgv provides a lightweight alternative to large genome browsers for visualizing biological annotations and dynamic analyses in their chromosomal context. It is particularly suited for applications ranging from short sequences to medium-sized genomes when the creation and maintenance of a large software and database infrastructure is not necessary or desired.

Identification of diverse carbon utilization pathways in Shewanella oneidensis MR-1 via expression profiling

To identify pathways of carbon utilization in the metal-reducing marine bacterium Shewanella oneidensis MR-1, we assayed the expression of cells grown with various carbon sources using a high-density oligonucleotide Affymetrix microarray. Our expression profiles reveal genes and regulatory mechanisms which govern the sensing, import, and utilization of the nucleoside inosine, the chitin monomer N-acetylglucosamine, and a casein-derived mixture of amino acids. Our analysis suggests a prominent role for the pentose-phosphate and Entner-Doudoroff pathways in energy metabolism, and regulatory coupling between carbon catabolism and electron acceptor pathways. In sum, these results indicate that S. oneidensis possesses a broader capacity for carbon utilization than previously reported, a view with implications for optimizing its role in microbial fuel cell and bioremediative applications.

Chemogenomic profiling on a genome-wide scale using reverse-engineered gene networks

A major challenge in drug discovery is to distinguish the molecular targets of a bioactive compound from the hundreds to thousands of additional gene products that respond indirectly to changes in the activity of the targets. Here, we present an integrated computational-experimental approach for computing the likelihood that gene products and associated pathways are targets of a compound. This is achieved by filtering the mRNA expression profile of compound-exposed cells using a reverse-engineered model of the cell's gene regulatory network. We apply the method to a set of 515 whole-genome yeast expression profiles resulting from a variety of treatments (compounds, knockouts and induced expression), and correctly enrich for the known targets and associated pathways in the majority of compounds examined. We demonstrate our approach with PTSB, a growth inhibitory compound with a previously unknown mode of action, by predicting and validating thioredoxin and thioredoxin reductase as its target.

Robust identification of large genetic networks

Temporal and spatial gene expression, together with the concentration of proteins and metabolites, is tightly controlled in the cell. This is possible thanks to complex regulatory networks between these different elements. The identification of these networks would be extremely valuable. We developed a novel algorithm to identify a large genetic network, as a set of linear differential equations, starting from measurements of gene expression at steady state following transcriptional perturbations. Experimentally, it is possible to overexpress each of the genes in the network using an episomal expression plasmid and measure the change in mRNA concentration of all the genes, following the perturbation. Computationally, we reduced the identification problem to a multiple linear regression, assuming that the network is sparse. We implemented a heuristic search method in order to apply the algorithm to large networks. The algorithm can correctly identify the network, even in the presence of large noise in the data, and can be used to predict the genes that directly mediate the action of a compound. Our novel approach is experimentally feasible and it is readily applicable to large genetic networks.

Programmable cells: interfacing natural and engineered gene networks

Novel cellular behaviors and characteristics can be obtained by coupling engineered gene networks to the cell's natural regulatory circuitry through appropriately designed input and output interfaces. Here, we demonstrate how an engineered genetic circuit can be used to construct cells that respond to biological signals in a predetermined and programmable fashion. We employ a modular design strategy to create Escherichia coli strains where a genetic toggle switch is interfaced with: (i) the SOS signaling pathway responding to DNA damage, and (ii) a transgenic quorum sensing signaling pathway from Vibrio fischeri. The genetic toggle switch endows these strains with binary response dynamics and an epigenetic inheritance that supports a persistent phenotypic alteration in response to transient signals. These features are exploited to engineer cells that form biofilms in response to DNA-damaging agents and cells that activate protein synthesis when the cell population reaches a critical density. Our work represents a step toward the development of "plug-and-play" genetic circuitry that can be used to create cells with programmable behaviors.

Inferring genetic networks and identifying compound mode of action via expression profiling

The complexity of cellular gene, protein, and metabolite networks can hinder attempts to elucidate their structure and function. To address this problem, we used systematic transcriptional perturbations to construct a first-order model of regulatory interactions in a nine-gene subnetwork of the SOS pathway in Escherichia coli. The model correctly identified the major regulatory genes and the transcriptional targets of mitomycin C activity in the subnetwork. This approach, which is experimentally and computationally scalable, provides a framework for elucidating the functional properties of genetic networks and identifying molecular targets of pharmacological compounds.

Construction of a genetic toggle switch in Escherichia coli

It has been proposed' that gene-regulatory circuits with virtually any desired property can be constructed from networks of simple regulatory elements. These properties, which include multistability and oscillations, have been found in specialized gene circuits such as the bacteriophage lambda switch and the Cyanobacteria circadian oscillator. However, these behaviours have not been demonstrated in networks of non-specialized regulatory components. Here we present the construction of a genetic toggle switch-a synthetic, bistable gene-regulatory network-in Escherichia coli and provide a simple theory that predicts the conditions necessary for bistability. The toggle is constructed from any two repressible promoters arranged in a mutually inhibitory network. It is flipped between stable states using transient chemical or thermal induction and exhibits a nearly ideal switching threshold. As a practical device, the toggle switch forms a synthetic, addressable cellular memory unit and has implications for biotechnology, biocomputing and gene therapy.

Problem of track offset in optical disk systems

In an optical disk drive, it is well known that a tilt of the disk causes an offset in the tracking-error signal (TES). One effect of disk tilt is the introduction of a dc component to the TES, which can be largely corrected by operation of the tracking system at the midpoint between the maximum and the minimum values of the open-loop TES. However, this method of correcting for the dc shift in the TES does not correct for the effect of coma in the focused spot, which leads to track offset. The track offset of a system is defined as the distance between the peak irradiance in the focused spot and the center of the groove when the tracking system is operating at the midpoint between the maximum and the minimum values of the open-loop TES in the presence of disk tilt. Calculations are performed that show the dependence of track offset on various system parameters, including track pitch, wavelength, and numerical aperture and rim intensity of the objective lens, and on the regions of the beam used to generate the TES. The track offsets for several beam-segmentation schemes are calculated for a digital versatile disk that uses push-pull and differential phase tracking. It is shown that for differential phase tracking the value of track offset depends on the mark length.

A theory for controlling cell cycle dynamics using a reversibly binding inhibitor

We demonstrate, by using mathematical modeling of cell division cycle (CDC) dynamics, a potential mechanism for precisely controlling the frequency of cell division and regulating the size of a dividing cell. Control of the cell cycle is achieved by artificially expressing a protein that reversibly binds and inactivates any one of the CDC proteins. In the simplest case, such as the checkpoint-free situation encountered in early amphibian embryos, the frequency of CDC oscillations can be increased or decreased by regulating the rate of synthesis, the binding rate, or the equilibrium constant of the binding protein. In a more complex model of cell division, where size-control checkpoints are included, we show that the same reversible binding reaction can alter the mean cell mass in a continuously dividing cell. Because this control scheme is general and requires only the expression of a single protein, it provides a practical means for tuning the characteristics of the cell cycle in vivo.