NEMix: Single-cell Nested Effects Models for Probabilistic Pathway Stimulation

Nested effects models have been used successfully for learning subcellular networks from high-dimensional perturbation effects that result from RNA interference (RNAi) experiments. Here, we further develop the basic nested effects model using high-content single-cell imaging data from RNAi screens of cultured cells infected with human rhinovirus. RNAi screens with single-cell readouts are becoming increasingly common, and they often reveal high cell-to-cell variation. As a consequence of this cellular heterogeneity, knock-downs result in variable effects among cells and lead to weak average phenotypes on the cell population level. To address this confounding factor in network inference, we explicitly model the stimulation status of a signaling pathway in individual cells. We extend the framework of nested effects models to probabilistic combinatorial knock-downs and propose NEMix, a nested effects mixture model that accounts for unobserved pathway activation. We analyzed the identifiability of NEMix and developed a parameter inference scheme based on the Expectation Maximization algorithm. In an extensive simulation study, we show that NEMix improves learning of pathway structures over classical NEMs significantly in the presence of hidden pathway stimulation. We applied our model to single-cell imaging data from RNAi screens monitoring human rhinovirus infection, where limited infection efficiency of the assay results in uncertain pathway stimulation. Using a subset of genes with known interactions, we show that the inferred NEMix network has high accuracy and outperforms the classical nested effects model without hidden pathway activity. NEMix is implemented as part of the R/Bioconductor package ‘nem’ and available at www.cbg.ethz.ch/software/NEMix.     

Scrutinizing MHC-I Binding Peptides and Their Limits of Variation

Designed peptides that bind to major histocompatibility protein I (MHC-I) allomorphs bear the promise of representing epitopes that stimulate a desired immune response. A rigorous bioinformatical exploration of sequence patterns hidden in peptides that bind to the mouse MHC-I allomorph H-2K^b is presented. We exemplify and validate these motif findings by systematically dissecting the epitope SIINFEKL and analyzing the resulting fragments for their binding potential to H-2K^b in a thermal denaturation assay. The results demonstrate that only fragments exclusively retaining the carboxy- or amino-terminus of the reference peptide exhibit significant binding potential, with the N-terminal pentapeptide SIINF as shortest ligand. This study demonstrates that sophisticated machine-learning algorithms excel at extracting fine-grained patterns from peptide sequence data and predicting MHC-I binding peptides, thereby considerably extending existing linear prediction models and providing a fresh view on the computer-based molecular design of future synthetic vaccines. The server for prediction is available at http://modlab-cadd.ethz.ch (SLiDER tool, MHC-I version 2012).     

Malfunctioning of adipocytes in obesity is linked to quantitative surfaceome changes

Increased triglyceride accumulation in adipocytes caused by a misbalance between energy intake and energy consumption, results in increased adipocyte size, excess adipose tissue, increased body weight and ultimately, obesity. It is well established that enlarged adipocytes exhibit malfunctions that contribute to whole body insulin resistance, a key factor for the development of type 2 diabetes. However, the underlying molecular cause for dysfunctional adipocyte behavior and signaling is poorly understood. Since the adipocyte cell surface proteome, or surfaceome, represents the cellular signaling gateway to the microenvironment, we studied the contribution of this subproteome to adipocyte malfunctions in obesity. By using the chemoproteomic Cell Surface Capture (CSC) technology, we established surfaceome maps of primary adipocytes derived from different mouse models for metabolic disorders. Relative quantitative comparison between these surfaceome maps revealed a set of cell surface glycoproteins with modulated location-specific abundance levels. RNAi mediated targeting of a subset of the detected obesity modulated cell surface glycoproteins in an in vitro model system provided functional evidence for their role in adiponectin secretion and the lipolytic activity of adipocytes. Thus, we conclude that the identified cell surface glycoproteins which exhibit obesity induced abundance changes and impact adipocyte function at the same time contribute to adipocyte malfunction in obesity. The regulation of their concerted activities could improve adipocyte function in obesity.     

Large-Scale Proteomics of the Cassava Storage Root and Identification of a Target Gene to Reduce Postharvest Deterioration

Cassava (Manihot esculenta) is the most important root crop in the tropics, but rapid postharvest physiological deterioration (PPD) of the root is a major constraint to commercial cassava production. We established a reliable method for image-based PPD symptom quantification and used label-free quantitative proteomics to generate an extensive cassava root and PPD proteome. Over 2600 unique proteins were identified in the cassava root, and nearly 300 proteins showed significant abundance regulation during PPD. We identified protein abundance modulation in pathways associated with oxidative stress, phenylpropanoid biosynthesis (including scopoletin), the glutathione cycle, fatty acid α-oxidation, folate transformation, and the sulfate reduction II pathway. Increasing protein abundances and enzymatic activities of glutathione-associated enzymes, including glutathione reductases, glutaredoxins, and glutathione S-transferases, indicated a key role for ascorbate/glutathione cycles. Based on combined proteomics data, enzymatic activities, and lipid peroxidation assays, we identified glutathione peroxidase as a candidate for reducing PPD. Transgenic cassava overexpressing a cytosolic glutathione peroxidase in storage roots showed delayed PPD and reduced lipid peroxidation as well as decreased H₂O₂ accumulation. Quantitative proteomics data from ethene and phenylpropanoid pathways indicate additional gene candidates to further delay PPD. Cassava root proteomics data are available at www.pep2pro.ethz.ch for easy access and comparison with other proteomics data.     

A Conserved GEF for Rho-Family GTPases Acts in an EGF Signaling Pathway to Promote Sleep-like Quiescence in Caenorhabditis elegans

Sleep is evolutionarily conserved and required for organism homeostasis and survival. Despite this importance, the molecular and cellular mechanisms underlying sleep are not well understood. Caenorhabditis elegans exhibits sleep-like behavioral quiescence and thus provides a valuable, simple model system for the study of cellular and molecular regulators of this process. In C. elegans, epidermal growth factor receptor (EGFR) signaling is required in the neurosecretory neuron ALA to promote sleep-like behavioral quiescence after cellular stress. We describe a novel role for VAV-1, a conserved guanine nucleotide exchange factor (GEF) for Rho-family GTPases, in regulation of sleep-like behavioral quiescence. VAV-1, in a GEF-dependent manner, acts in ALA to suppress locomotion and feeding during sleep-like behavioral quiescence in response to cellular stress. Additionally, VAV-1 activity is required for EGF-induced sleep-like quiescence and normal levels of EGFR and secretory dense core vesicles in ALA. Importantly, the role of VAV-1 in promoting cellular stress–induced behavioral quiescence is vital for organism health because VAV-1 is required for normal survival after cellular stress.     

A Novel Cell Death Gene Acts to Repair Patterning Defects in Drosophila melanogaster

Cell death is a mechanism utilized by organisms to eliminate excess cells during development. Here, we describe a novel regulator of caspase-independent cell death, Mabiki (Mabi), that is involved in the repair of the head patterning defects caused by extra copies of bicoid in Drosophila melanogaster. Mabiki functions together with caspase-dependent cell death mechanisms to provide robustness during development.     

LifeMap Discovery?: The Embryonic Development,Stem Cells,and Regenerative Medicine Research Portal

LifeMap Discovery™ provides investigators with an integrated database of embryonic development, stem cell biology and regenerative medicine. The hand-curated reconstruction of cell ontology with stem cell biology; including molecular, cellular, anatomical and disease-related information, provides efficient and easy-to-use, searchable research tools. The database collates in vivo and in vitro gene expression and guides translation from in vitro data to the clinical utility, and thus can be utilized as a powerful tool for research and discovery in stem cell biology, developmental biology, disease mechanisms and therapeutic discovery. LifeMap Discovery is freely available to academic nonprofit institutions at http://discovery.lifemapsc.com     

Functional Genomic Identification of Genes Required for Male Gonadal Differentiation in Caenorhabditis elegans

The Caenorhabditis elegans somatic gonad develops from a four-cell primordium into a mature organ that differs dramatically between the sexes in overall morphology (two arms in hermaphrodites and one in males) and in the cell types comprising it. Gonadal development in C. elegans is well studied, but regulation of sexual differentiation, especially later in gonadal development, remains poorly elucidated. To identify genes involved in this process, we performed a genome-wide RNAi screen using sex-specifically expressed gonadal GFP reporters. This screen identified several phenotypic classes, including ∼70 genes whose depletion feminized male gonadal cells. Among the genes required for male cell fate specification are Wnt/β-catenin pathway members, cell cycle regulators, and genes required for mitotic spindle function and cytokinesis. We find that a Wnt/β-catenin pathway independent of extracellular Wnt ligand is essential for asymmetric cell divisions and male differentiation during gonadal development in larvae. We also find that the cell cycle regulators cdk-1 and cyb-3 and the spindle/cytokinesis regulator zen-4 are required for Wnt/β-catenin pathway activity in the developing gonad. After sex is determined in the gonadal primordium the global sex determination pathway is dispensable for gonadal sexual fate, suggesting that male cell fates are promoted and maintained independently of the global pathway during this period.THE Caenorhabditis elegans gonad derives from a simple primordium of four cells that coalesces during embryogenesis and contains two somatic gonad precursors (SGPs), Z1 and Z4, flanking two germline precursors, Z2 and Z3 (Kimble and Hirsh 1979). The SGPs undergo very different developmental programs in each sex, involving sexually dimorphic cell lineages and migrations and sex-specific cellular differentiation. The result is a two-armed bilaterally symmetrical gonad in the adult hermaphrodite or a single-armed asymmetric gonad in the adult male. The high degree of sexual dimorphism of the mature organ and variety of cellular events that occur sex specifically during its development make the C. elegans gonad an outstanding model for sex-specific organogenesis.Development of the somatic gonad occurs in two phases. The early phase defines the gonadal axes and establishes the precursors of the major gonadal cell types. This takes place during the first larval stage (L1), beginning shortly after hatching with the first division of the SGPs. In both sexes SGP division is asymmetric in terms of both the sizes and the fates of the daughter cells, and establishes the proximal/distal axis of the gonad (Hirsh et al. 1976; Kimble and Hirsh 1979). The global sex determination pathway establishes the future sex of the gonad around the time of hatching (Klass et al. 1976; Nelson et al. 1978), and sexual dimorphism is already apparent when the SGPs divide: the size asymmetry of the SGP daughters is much more pronounced in males than hermaphrodites. In both sexes the asymmetry of the first SGP division requires a Wnt/β-catenin pathway. Mutations compromising this pathway cause a “symmetrical sisters” phenotype in which both daughters adopt the same fate (Miskowski et al. 2001; Siegfried and Kimble 2002; Phillips and Kimble 2009). Sex specificity is imposed on the SGPs by the global sex determining gene tra-1 (Hodgkin 1987) and the gonad-specific sex determining gene fkh-6 (Chang et al. 2004). These genes play opposing roles in SGP sex determination, with tra-1 feminizing and fkh-6 masculinizing the somatic gonad, and they also act redundantly to promote mitotic proliferation of the SGP lineage (Chang et al. 2004). SGP sex determination is linked to cell cycle progression by cyclin D, which is required to overcome repression of fkh-6 expression in the SGPs by E2F (Tilmann and Kimble 2005).The later phase of gonadal development involves the elongation of the gonad, together with cellular proliferation and differentiation, and lasts from L2 to adulthood. During L2 the somatic cells enlarge and leader cells (distal tip cells in the hermaphrodite, linker cell in the male) begin long-range migrations that extend the gonad. During L3, somatic gonad cell division resumes in both sexes, leading to the formation of differentiated somatic cell types by the end of L3 or beginning of L4. Gonadal morphogenesis is completed and gametogenesis begins during L4 (Kimble and Hirsh 1979).Although SGP division and much of hermaphrodite gonadal development have been well studied (Hubbard and Greenstein 2000), sexual cell fate specification in the somatic gonad is more poorly understood, particularly after the L1 stage. Despite the importance of fkh-6 in promoting male differentiation, it is expressed in males only during early L1 and null mutants have incomplete gonadal sex reversal. We have therefore performed a genome-wide RNAi screen to identify additional genes required after hatching for gonadal development in each sex. Among the advantages of this approach is the ability to identify gonadal regulators that also are essential for embryonic development. To our knowledge this is the first functional genomic study of gonadal sex differentiation.The screen identified many genes whose depletion disrupts gonadogenesis in each sex and nearly 70 genes whose depletion causes gonadal feminization in males. Prominent among this latter class were components of a Wnt/β-catenin pathway, cell cycle regulators, and genes involved in mitotic spindle function and cytokinesis. We find that Wnt/β-catenin activity continues in both sexes after SGP division and is required for male cell fate commitment in the gonad. We also find that the cyclin-dependent kinase cdk-1 and its cognate cyclin cyb-3 as well as the mitotic spindle regulator zen-4 are required for gonadal Wnt/β-catenin pathway activity, providing a potential new link between the cell cycle, asymmetric division, and sexual differentiation. The feminization caused by depletion of Wnt/β-catenin pathway components or cdk-1 is independent of the global sex determination pathway, suggesting that sexual fates in the male gonad remain plastic after the primary sex determination decision.     

Host Cell Factors in HIV Replication: Meta-Analysis of Genome-Wide Studies

We have analyzed host cell genes linked to HIV replication that were identified in nine genome-wide studies, including three independent siRNA screens. Overlaps among the siRNA screens were very modest (<7% for any pairwise combination), and similarly, only modest overlaps were seen in pairwise comparisons with other types of genome-wide studies. Combining all genes from the genome-wide studies together with genes reported in the literature to affect HIV yields 2,410 protein-coding genes, or fully 9.5% of all human genes (though of course some of these are false positive calls). Here we report an “encyclopedia” of all overlaps between studies (available at http://www.hostpathogen.org), which yielded a more extensively corroborated set of host factors assisting HIV replication. We used these genes to calculate refined networks that specify cellular subsystems recruited by HIV to assist in replication, and present additional analysis specifying host cell genes that are attractive as potential therapeutic targets.     

Asymmetric Wnt Pathway Signaling Facilitates Stem Cell-Like Divisions via the Nonreceptor Tyrosine Kinase FRK-1 in Caenorhabditis elegans

Asymmetric cell division is critical during development, as it influences processes such as cell fate specification and cell migration. We have characterized FRK-1, a homolog of the mammalian Fer nonreceptor tyrosine kinase, and found it to be required for differentiation and maintenance of epithelial cell types, including the stem cell-like seam cells of the hypodermis. A genomic knockout of frk-1, allele ok760, results in severely uncoordinated larvae that arrest at the L1 stage and have an excess number of lateral hypodermal cells that appear to have lost asymmetry in the stem cell-like divisions of the seam cell lineage. frk-1(ok760) mutants show that there are excess lateral hypodermal cells that are abnormally shaped and smaller in size compared to wild type, a defect that could be rescued only in a manner dependent on the kinase activity of FRK-1. Additionally, we observed a significant change in the expression of heterochronic regulators in frk-1(ok760) mutants. However, frk-1(ok760) mutants do not express late, nonseam hypodermal GFP markers, suggesting the seam cells do not precociously differentiate as adult-hyp7 cells. Finally, our data also demonstrate a clear role for FRK-1 in seam cell proliferation, as eliminating FRK-1 during the L3–L4 transition results in supernumerary seam cell nuclei that are dependent on asymmetric Wnt signaling. Specifically, we observe aberrant POP-1 and WRM-1 localization that is dependent on the presence of FRK-1 and APR-1. Overall, our data suggest a requirement for FRK-1 in maintaining the identity and proliferation of seam cells primarily through an interaction with the asymmetric Wnt pathway.     

ReCount: A multi-experiment resource of analysis-ready RNA-seq gene count datasets

Background

With next-generation sequencing technologies, experiments that were considered prohibitive only a few years ago are now possible. However, while these technologies have the ability to produce enormous volumes of data, the sequence reads are prone to error. This poses fundamental hurdles when genetic diversity is investigated.

Results

We developed ShoRAH, a computational method for quantifying genetic diversity in a mixed sample and for identifying the individual clones in the population, while accounting for sequencing errors. The software was run on simulated data and on real data obtained in wet lab experiments to assess its reliability.

Conclusions

ShoRAH is implemented in C++, Python, and Perl and has been tested under Linux and Mac OS X. Source code is available under the GNU General Public License at http://www.cbg.ethz.ch/software/shorah.     

Flow Cytometry Bioinformatics

Flow cytometry bioinformatics is the application of bioinformatics to flow cytometry data, which involves storing, retrieving, organizing, and analyzing flow cytometry data using extensive computational resources and tools. Flow cytometry bioinformatics requires extensive use of and contributes to the development of techniques from computational statistics and machine learning. Flow cytometry and related methods allow the quantification of multiple independent biomarkers on large numbers of single cells. The rapid growth in the multidimensionality and throughput of flow cytometry data, particularly in the 2000s, has led to the creation of a variety of computational analysis methods, data standards, and public databases for the sharing of results. Computational methods exist to assist in the preprocessing of flow cytometry data, identifying cell populations within it, matching those cell populations across samples, and performing diagnosis and discovery using the results of previous steps. For preprocessing, this includes compensating for spectral overlap, transforming data onto scales conducive to visualization and analysis, assessing data for quality, and normalizing data across samples and experiments. For population identification, tools are available to aid traditional manual identification of populations in two-dimensional scatter plots (gating), to use dimensionality reduction to aid gating, and to find populations automatically in higher dimensional space in a variety of ways. It is also possible to characterize data in more comprehensive ways, such as the density-guided binary space partitioning technique known as probability binning, or by combinatorial gating. Finally, diagnosis using flow cytometry data can be aided by supervised learning techniques, and discovery of new cell types of biological importance by high-throughput statistical methods, as part of pipelines incorporating all of the aforementioned methods. Open standards, data, and software are also key parts of flow cytometry bioinformatics. Data standards include the widely adopted Flow Cytometry Standard (FCS) defining how data from cytometers should be stored, but also several new standards under development by the International Society for Advancement of Cytometry (ISAC) to aid in storing more detailed information about experimental design and analytical steps. Open data is slowly growing with the opening of the CytoBank database in 2010 and FlowRepository in 2012, both of which allow users to freely distribute their data, and the latter of which has been recommended as the preferred repository for MIFlowCyt-compliant data by ISAC. Open software is most widely available in the form of a suite of Bioconductor packages, but is also available for web execution on the GenePattern platform.

This is a “Topic Page” article for PLOS Computational Biology.

     

Caenorhabditis elegans PIG-1/MELK Acts in a Conserved PAR-4/LKB1 Polarity Pathway to Promote Asymmetric Neuroblast Divisions

Asymmetric cell divisions produce daughter cells with distinct sizes and fates, a process important for generating cell diversity during development. Many Caenorhabditis elegans neuroblasts, including the posterior daughter of the Q cell (Q.p), divide to produce a larger neuron or neuronal precursor and a smaller cell that dies. These size and fate asymmetries require the gene pig-1, which encodes a protein orthologous to vertebrate MELK and belongs to the AMPK-related family of kinases. Members of this family can be phosphorylated and activated by the tumor suppressor kinase LKB1, a conserved polarity regulator of epithelial cells and neurons. In this study, we present evidence that the C. elegans orthologs of LKB1 (PAR-4) and its partners STRAD (STRD-1) and MO25 (MOP-25.2) regulate the asymmetry of the Q.p neuroblast division. We show that PAR-4 and STRD-1 act in the Q lineage and function genetically in the same pathway as PIG-1. A conserved threonine residue (T169) in the PIG-1 activation loop is essential for PIG-1 activity, consistent with the model that PAR-4 (or another PAR-4-regulated kinase) phosphorylates and activates PIG-1. We also demonstrate that PIG-1 localizes to centrosomes during cell divisions of the Q lineage, but this localization does not depend on T169 or PAR-4. We propose that a PAR-4-STRD-1 complex stimulates PIG-1 kinase activity to promote asymmetric neuroblast divisions and the generation of daughter cells with distinct fates. Changes in cell fate may underlie many of the abnormal behaviors exhibited by cells after loss of PAR-4 or LKB1.