Bayesian inference of the number of factors in gene-expression analysis: application to human virus challenge studies.

BACKGROUND: Nonparametric Bayesian techniques have been developed recently to extend the sophistication of factor models, allowing one to infer the number of appropriate factors from the observed data. We consider such techniques for sparse factor analysis, with application to gene-expression data from three virus challenge studies. Particular attention is placed on employing the Beta Process (BP), the Indian Buffet Process (IBP), and related sparseness-promoting techniques to infer a proper number of factors. The posterior density function on the model parameters is computed using Gibbs sampling and variational Bayesian (VB) analysis. RESULTS: Time-evolving gene-expression data are considered for respiratory syncytial virus (RSV), Rhino virus, and influenza, using blood samples from healthy human subjects. These data were acquired in three challenge studies, each executed after receiving institutional review board (IRB) approval from Duke University. Comparisons are made between several alternative means of per-forming nonparametric factor analysis on these data, with comparisons as well to sparse-PCA and Penalized Matrix Decomposition (PMD), closely related non-Bayesian approaches. CONCLUSIONS: Applying the Beta Process to the factor scores, or to the singular values of a pseudo-SVD construction, the proposed algorithms infer the number of factors in gene-expression data. For real data the "true" number of factors is unknown; in our simulations we consider a range of noise variances, and the proposed Bayesian models inferred the number of factors accurately relative to other methods in the literature, such as sparse-PCA and PMD. We have also identified a "pan-viral" factor of importance for each of the three viruses considered in this study. We have identified a set of genes associated with this pan-viral factor, of interest for early detection of such viruses based upon the host response, as quantified via gene-expression data.

A noisy linear map underlies oscillations in cell size and gene expression in bacteria.

During bacterial growth, a cell approximately doubles in size before division, after which it splits into two daughter cells. This process is subjected to the inherent perturbations of cellular noise and thus requires regulation for cell-size homeostasis. The mechanisms underlying the control and dynamics of cell size remain poorly understood owing to the difficulty in sizing individual bacteria over long periods of time in a high-throughput manner. Here we measure and analyse long-term, single-cell growth and division across different Escherichia coli strains and growth conditions. We show that a subset of cells in a population exhibit transient oscillations in cell size with periods that stretch across several (more than ten) generations. Our analysis reveals that a simple law governing cell-size control-a noisy linear map-explains the origins of these cell-size oscillations across all strains. This noisy linear map implements a negative feedback on cell-size control: a cell with a larger initial size tends to divide earlier, whereas one with a smaller initial size tends to divide later. Combining simulations of cell growth and division with experimental data, we demonstrate that this noisy linear map generates transient oscillations, not just in cell size, but also in constitutive gene expression. Our work provides new insights into the dynamics of bacterial cell-size regulation with implications for the physiological processes involved.

Transsynaptic Tracing from Peripheral Targets with Pseudorabies Virus Followed by Cholera Toxin and Biotinylated Dextran Amines Double Labeling.

Transsynaptic tracing has become a powerful tool used to analyze central efferents that regulate peripheral targets through multi-synaptic circuits. This approach has been most extensively used in the brain by utilizing the swine pathogen pseudorabies virus (PRV)(1). PRV does not infect great apes, including humans, so it is most commonly used in studies on small mammals, especially rodents. The pseudorabies strain PRV152 expresses the enhanced green fluorescent protein (eGFP) reporter gene and only crosses functional synapses retrogradely through the hierarchical sequence of synaptic connections away from the infection site(2,3). Other PRV strains have distinct microbiological properties and may be transported in both directions (PRV-Becker and PRV-Kaplan)(4,5). This protocol will deal exclusively with PRV152. By delivering the virus at a peripheral site, such as muscle, it is possible to limit the entry of the virus into the brain through a specific set of neurons. The resulting pattern of eGFP signal throughout the brain then resolves the neurons that are connected to the initially infected cells. As the distributed nature of transsynaptic tracing with pseudorabies virus makes interpreting specific connections within an identified network difficult, we present a sensitive and reliable method employing biotinylated dextran amines (BDA) and cholera toxin subunit b (CTb) for confirming the connections between cells identified using PRV152. Immunochemical detection of BDA and CTb with peroxidase and DAB (3, 3'-diaminobenzidine) was chosen because they are effective at revealing cellular processes including distal dendrites(6-11).