BTR: training asynchronous Boolean models using single-cell expression data

Abstract Background Rapid technological innovation for the generation of single-cell genomics data presents new challenges and opportunities for bioinformatics analysis. One such area lies in the development of new ways to train gene regulatory networks. The use of single-cell expression profiling technique allows the profiling of the expression states of hundreds of cells, but these expression states are typically noisier due to the presence of technical artefacts such as drop-outs. While many algorithms exist to infer a gene regulatory network, very few of them are able to harness the extra expression states present in single-cell expression data without getting adversely affected by the substantial technical noise present. Results Here we introduce BTR, an algorithm for training asynchronous Boolean models with single-cell expression data using a novel Boolean state space scoring function. BTR is capable of refining existing Boolean models and reconstructing new Boolean models by improving the match between model prediction and expression data. We demonstrate that the Boolean scoring function performed favourably against the BIC scoring function for Bayesian networks. In addition, we show that BTR outperforms many other network inference algorithms in both bulk and single-cell synthetic expression data. Lastly, we introduce two case studies, in which we use BTR to improve published Boolean models in order to generate potentially new biological insights. Conclusions BTR provides a novel way to refine or reconstruct Boolean models using single-cell expression data. Boolean model is particularly useful for network reconstruction using single-cell data because it is more robust to the effect of drop-outs. In addition, BTR does not assume any relationship in the expression states among cells, it is useful for reconstructing a gene regulatory network with as few assumptions as possible. Given the simplicity of Boolean models and the rapid adoption of single-cell genomics by biologists, BTR has the potential to make an impact across many fields of biomedical research

Botrytis cinerea BcCDI1 protein triggers both plant cell death and immune response

Cell death-inducing proteins (CDIPs) play important roles in the infection of Botrytis cinerea, a broad host-range necrotrophic phytopathogen. Here, we show that the secreted protein BcCDI1 (Cell Death Inducing 1) can cause necrosis in tobacco leaves and at the same time elicit plant defense. The transcription of Bccdi1 was induced at the infection stage. Deletion or overexpression of Bccdi1 resulted in no notable change in disease lesion on bean, tobacco, and Arabidopsis leaves, indicating that Bccdi1 has no effect on the final outcome of B. cinerea infection. Furthermore, the plant receptor-like kinases BAK1 and SOBIR1 are required to transduce the cell death-promoting signal induced by BcCDI1. These findings suggest that BcCDI1 is possibly recognized by plant receptors and then induces plant cell death

An experimentally validated network of nine haematopoietic transcription factors reveals mechanisms of cell state stability.

Transcription factor (TF) networks determine cell-type identity by establishing and maintaining lineage-specific expression profiles, yet reconstruction of mammalian regulatory network models has been hampered by a lack of comprehensive functional validation of regulatory interactions. Here, we report comprehensive ChIP-Seq, transgenic and reporter gene experimental data that have allowed us to construct an experimentally validated regulatory network model for haematopoietic stem/progenitor cells (HSPCs). Model simulation coupled with subsequent experimental validation using single cell expression profiling revealed potential mechanisms for cell state stabilisation, and also how a leukaemogenic TF fusion protein perturbs key HSPC regulators. The approach presented here should help to improve our understanding of both normal physiological and disease processes.Research in the authors’ laboratories was supported by Bloodwise, The Wellcome Trust, Cancer Research UK, the Biotechnology and Biological Sciences Research Council, the National Institute of Health Research, the Medical Research Council, the MRC Molecular Haematology Unit (Oxford) core award, a Weizmann-UK “Making Connections” grant (Oxford) and core support grants by the Wellcome Trust to the Cambridge Institute for Medical Research (100140) and Wellcome Trust–MRC Cambridge Stem Cell Institute (097922).This is the final version of the article. It first appeared from eLife via http://dx.doi.org/10.7554/eLife.1146

Study of e+e−→ppˉe^+e^- \rightarrow p\bar{p} in the vicinity of ψ(3770)\psi(3770)

Using 2917 $\rm{pb}^{-1}$ of data accumulated at 3.773~$\rm{GeV}$, 44.5~$\rm{pb}^{-1}$ of data accumulated at 3.65~$\rm{GeV}$ and data accumulated during a $\psi(3770)$ line-shape scan with the BESIII detector, the reaction $e^+e^-\rightarrow p\bar{p}$ is studied considering a possible interference between resonant and continuum amplitudes. The cross section of $e^+e^-\rightarrow\psi(3770)\rightarrow p\bar{p}$, $\sigma(e^+e^-\rightarrow\psi(3770)\rightarrow p\bar{p})$, is found to have two solutions, determined to be ($0.059\pm0.032\pm0.012$) pb with the phase angle $\phi = (255.8\pm37.9\pm4.8)^\circ$ ($<$0.11 pb at the 90% confidence level), or $\sigma(e^+e^-\rightarrow\psi(3770)\rightarrow p\bar{p}) = (2.57\pm0.12\pm0.12$) pb with $\phi = (266.9\pm6.1\pm0.9)^\circ$ both of which agree with a destructive interference. Using the obtained cross section of $\psi(3770)\rightarrow p\bar{p}$, the cross section of $p\bar{p}\rightarrow \psi(3770)$, which is useful information for the future PANDA experiment, is estimated to be either ($9.8\pm5.7$) nb ($<17.2$ nb at 90% C.L.) or $(425.6\pm42.9)$ nb

Using probabilistic graphical models to reconstruct biological networks and linkage maps

Probabilistic graphical models (PGMs) offer a conceptual architecture where biological and mathematical objects can be expressed with a common, intuitive formalism. This facilitates the joint development of statistical and computational tools for quantitative analysis of biological data. Over the last few decades, procedures based on well-understood principles for constructing PGMs from observational and experimental data have been studied extensively, and they thus form a model-based methodology for analysis and discovery. In this thesis, we further explore the potential of this methodology in systems biology and quantitative genetics, and illustrate the capabilities of our proposed approaches by several applications to both real and simulated omics data. In quantitative genetics, we partition phenotypic variation into heritable, genetic, and non-heritable, environmental, parts. In molecular genetics, we identify chromosomal regions that drive genetic variation: quantitative trait loci (QTLs). In systems genetics, we would like to answer the question of whether relations between multiple phenotypic traits can be organized within wholly or partially directed network structures. Directed edges in those networks can be interpreted as causal relationships, causality meaning that the consequences of interventions are predictable: phenotypic interventions in upstream traits, i.e. traits occurring early in causal chains, will produce changes in downstream traits. The effect of a QTL allele can be considered to represent a genetic intervention on the phenotypic network. Various methods have been proposed for statistical reconstruction of causal phenotypic networks exploiting previously identified QTLs. In chapter 2, we present a novel heuristic search algorithm, namely the QTL+phenotype supervised orientation (QPSO) algorithm, to infer causal relationships between phenotypic traits. Our algorithm shows good performance in the common, but so far uncovered case, where some traits come without QTLs. Therefore, our algorithm is especially attractive for applications involving expensive phenotypes, like metabolites, where relatively few genotypes can be measured and population size is limited. Standard QTL mapping typically models phenotypic variations observable in nature in relation to genetic variation in gene expression, regardless of multiple intermediate-level biological variations. In chapter 3, we present an approach integrating Gaussian graphical modeling (GGM) and causal inference for simultaneous modeling of multilevel biological responses to DNA variations. More specifically, for ripe tomato fruits, the dependencies of 24 sensory traits on 29 metabolites and the dependencies of all the sensory and metabolic traits further on 21 QTLs were investigated by three GGM approaches including: (i) lasso-based neighborhood selection in combination with a stability approach to regularization selection, (ii) the PC-skeleton algorithm and (iii) the Lasso in combination with stability selection, and then followed by the QPSO algorithm. The inferred dependency network which, though not essentially representing biological pathways, suggests how the effects of allele substitutions propagate through multilevel phenotypes. Such simultaneous study of the underlying genetic architecture and multifactorial interactions is expected to enhance the prediction and manipulation of complex traits. And it is applicable to a range of population structures, including offspring populations from crosses between inbred parents and outbred parents, association panels and natural populations. In chapter 4, we report a novel method for linkage map construction using probabilistic graphical models. It has been shown that linkage map construction can be hampered by the presence of genotyping errors and chromosomal rearrangements such as inversions and translocations. Our proposed method is proven, both theoretically and practically, to be effective in filtering out markers that contain genotyping errors. In particular, it carries out marker filtering and ordering simultaneously, and is therefore superior to the standard post-hoc filtering using nearest-neighbour stress. Furthermore, we demonstrate empirically that the proposed method offers a promising solution to genetic map construction in the case of a reciprocal translocation. In the domain of PGMs, Bayesian networks (BNs) have proven, both theoretically and practically, to be a promising tool for the reconstruction of causal networks. In particular, the PC algorithm and the Metropolis-Hastings algorithm, which are representatives of mainstream methods to BN structure learning, are reported to have been successfully applied to the field of biology. In view of the fact that most biological systems exist in the form of random network or scale-free network, in chapter 5 we compare the performance of the two algorithms in constructing both random and scale-free BNs. Our simulation study shows that for either type of BN, the PC algorithm is superior to the M-H algorithm in terms of timeliness; the M-H algorithm is preferable to the PC algorithm when the completeness of reconstruction is emphasized; but when the fidelity of reconstruction is taken into account, the better one of the two algorithms varies from case to case. Moreover, whichever algorithm is adopted, larger sample sizes generally permit more accurate reconstructions, especially in regard to the completeness of the resulting networks. Finally, chapter 6 presents a further elaboration and discussion of the key concepts and results involved in this thesis.</p

Squamous Cell Carcinoma Antigen: Clinical Application and Research Status

The squamous cell carcinoma antigen (SCCA) is a tumor marker that has gained increasing attention for its biological functions and significance in normal physiological and pathological processes. Not only SCCA but also circulating immune complexes of SCCA and immunoglobulin M (IgM) are involved in normal physiological and pathological processes, providing a background for numerous clinical studies aimed at assessing the potential role of SCCA, SCCA–IgM, and SCCA isoform complexes in clinical practice. Previous studies support the clinical value of SCCA as a tumor marker for either diagnosing squamous cancers or monitoring the response to radiotherapy or chemotherapy, tumor relapse, and treatment failure. However, these studies show contrasting results, making the diagnostic or prognostic value of SCCA controversial. To reduce clinical heterogeneity across studies and achieve a more accurate and reliable comparison of results, a standardized detection method, scoring system, and cutoff level need to be established. Moreover, despite the fact that performances of different methods are comparable, the dynamic observation of tumor marker kinetics should be conducted under the same method

miRNAs regulate immune response and signaling during hepatitis C virus infection

Abstract Hepatitis C is one of the most common types of viral hepatitis that impair human health. At present, there is still no effective specific therapy for hepatitis C virus infection. As host immunity is an important mechanism to defend against or clear infections, the interactions between the virus and the host immune response are crucial to the progress of the disease. Of note, hepatitis C virus infection has been reported to regulate cellular miRNAs, which play significant roles in many processes, including infection and immunity. In this review, we describe how miRNAs regulate the host immune response to hepatitis C virus via complex signaling pathways

An example of LODG.

<p><i>Y</i>₁ and <i>Y</i>₂ are two correlated traits; <i>C</i>₂ and <i>C</i>₃ are two traits that have been newly determined as parent nodes of <i>Y</i>₁; <i>Y</i>₁, <i>C</i>₃ and <i>C</i>₅ are three traits newly determined as parent nodes of <i>Y</i>₂; <i>Y</i>₁ is a newly determined parent node of traits <i>C</i>₁ and <i>C</i>₄; <i>Y</i>₂ is a newly determined parent node of traits <i>C</i>₄ and <i>C</i>₆.</p

The general representations of resolvable LGPNs.

<p><i>Y</i>₁ and <i>Y</i>₂ are two correlated traits; <b><i>P</i></b>₁ = {<i>P</i>₁₁,…,<i>P</i>_1<i>k</i>} and <b><i>P</i></b>₂ = {<i>P</i>₂₁,…,<i>P</i>_2<i>l</i>} are, respectively, the unique parent nodes of <i>Y</i>₁ and <i>Y</i>₂; <b><i>P</i></b>₁₂ = {<i>P</i>₁,…,<i>P_s</i>} are the common parent nodes of <i>Y</i>₁ and <i>Y</i>₂; <b><i>C</i></b>₁ = {<i>C</i>₁₁,…,<i>C</i>_1<i>u</i>} and <b><i>C</i></b>₂ = {<i>C</i>₂₁,…,<i>C</i>_2<i>v</i>} are the unique neighboring traits of <i>Y</i>₁ and <i>Y</i>₂; <b><i>C</i></b>₁₂ = {<i>C</i>₁,…,<i>C_t</i>} are the common neighboring traits of <i>Y</i>₁ and <i>Y</i>₂. Note that each of the neighboring traits of <i>Y</i>₁ is nonadjacent to at least one of the parent nodes of <i>Y</i>₁, and the same is true of <i>Y</i>₂. Also note that <b><i>P</i></b>₁, <b><i>P</i></b>₂ and <b><i>P</i></b>₁₂ are allowed to have three different compositions: (1) a pure set of QTLs, if only genetic factors have been identified for <i>Y</i>₁ and/or <i>Y</i>₂; (2) a mixed set of QTLs and traits, if some traits in addition to QTLs have been determined to have causal effects on <i>Y</i>₁ and/or <i>Y</i>₂; (3) a pure set of traits, if only some traits have been found as causal factors of <i>Y</i>₁ and/or <i>Y</i>₂; in contrast, <b><i>C</i></b>₁, <b><i>C</i></b>₂ and <b><i>C</i></b>₁₂ only refer to those traits that are directly connected to <i>Y</i>₁ and/or <i>Y</i>₂ by an undirected edge. (A) The general representation of LGPNs where both <i>Y</i>₁ and <i>Y</i>₂ have parent nodes, and at least one of them has unique parent nodes; (B) the general representation of LGPNs where only <i>Y</i>₁ has parent nodes.</p