Publications

It remains unclear whether causal, rather than merely correlational, relationships in molecular networks can be inferred in complex biological settings. Here we describe the HPN-DREAM network inference challenge, which focused on learning causal influences in signaling networks. We used phosphoprotein data from cancer cell lines as well as in silico data from a nonlinear dynamical model. Using the phosphoprotein data, we scored more than 2,000 networks submitted by challenge participants. The networks spanned 32 biological contexts and were scored in terms of causal validity with respect to unseen interventional data. A number of approaches were effective, and incorporating known biology was generally advantageous. Additional sub-challenges considered time-course prediction and visualization. Our results suggest that learning causal relationships may be feasible in complex settings such as disease states. Furthermore, our scoring approach provides a practical way to empirically assess inferred molecular networks in a causal sense.

Interpreting the functional state of chromatin from the combinatorial binding patterns of chromatin factors, that is, the chromatin codes, is crucial for decoding the epigenetic state of the cell. Here we present a systematic map of Drosophila chromatin states derived from data-driven probabilistic modelling of dependencies between chromatin factors. Our model not only recapitulates enhancer-like chromatin states as indicated by widely used enhancer marks but also divides these states into three functionally distinct groups, of which only one specific group possesses active enhancer activity. Moreover, we discover a strong association between one specific enhancer state and RNA Polymerase II pausing, linking transcription regulatory potential and chromatin organization. We also observe that with the exception of long-intron genes, chromatin state transition positions in transcriptionally active genes align with an absolute distance to their corresponding transcription start site, regardless of gene length. Using our method, we provide a resource that helps elucidate the functional and spatial organization of the chromatin code landscape.

A number of problems in probability and statistics can be addressed using the multivariate normal (Gaussian) distribution. In the one-dimensional case, computing the probability for a given mean and variance simply requires the evaluation of the corresponding Gaussian density. In the $n$-dimensional setting, however, it requires the inversion of an $n \times n$ covariance matrix, $C$, as well as the evaluation of its determinant, $\det(C)$. In many cases, such as regression using Gaussian processes, the covariance matrix is of the form $C = \sigma^2 I + K$, where $K$ is computed using a specified covariance kernel which depends on the data and additional parameters (hyperparameters). The matrix $C$ is typically dense, causing standard direct methods for inversion and determinant evaluation to require $\mathcal O(n^3)$ work. This cost is prohibitive for large-scale modeling. Here, we show that for the most commonly used covariance functions, the matrix $C$ can be hierarchically factored into a product of block low-rank updates of the identity matrix, yielding an $\mathcal O (n\log^2 n) $ algorithm for inversion. More importantly, we show that this factorization enables the evaluation of the determinant $\det(C)$, permitting the direct calculation of probabilities in high dimensions under fairly broad assumptions on the kernel defining $K$. Our fast algorithm brings many problems in marginalization and the adaptation of hyperparameters within practical reach using a single CPU core. The combination of nearly optimal scaling in terms of problem size with high-performance computing resources will permit the modeling of previously intractable problems. We illustrate the performance of the scheme on standard covariance kernels.

Chirality or mirror asymmetry is a common theme in biology found in organismal body plans, tissue patterns and even in individual cells. In many cases the emergence of chirality is driven by actin cytoskeletal dynamics. Although it is well established that the actin cytoskeleton generates rotational forces at the molecular level, we are only beginning to understand how this can result in chiral behavior of the entire actin network in vivo. In this review, we will give an overview of actin driven chiralities across different length scales known until today. Moreover, we evaluate recent quantitative models demonstrating that chiral symmetry breaking of cells can be achieved by properly aligning molecular-scale torque generation processes in the actomyosin cytoskeleton.

The design of folded miniature proteins is predicated on establishing noncovalent interactions that direct the self-assembly of discrete thermostable tertiary structures. In this work, we describe how a network of cation−π interactions present in proteins containing “WSXWS motifs” can be emulated to stabilize the core of a miniature protein. This 19-residue protein sequence recapitulates a set of interdigitated arginine and tryptophan residues that stabilize a distinctive β-strand:loop:PPII-helix topology. Validation of the compact fold determined by NMR was carried out by mutagenesis of the cation−π network and by comparison to the corresponding disulfide-bridged structure. These results support the involvement of a coordinated set of cation−π interactions that stabilize the tertiary structure.

Rapid increases in e-cigarette use and potential exposure to harmful byproducts have shifted public health focus to e-cigarettes as a possible drug of abuse. Effective surveillance of use and prevalence would allow appropriate regulatory responses. An ideal surveillance system would collect usage data in real time, focus on populations of interest, include populations unable to take the survey, allow a breadth of questions to answer, and enable geo-location analysis. Social media streams may provide this ideal system. To realize this use case, a foundational question is whether we can detect ecigarette use at all. This work reports two pilot tasks using text classification to identify automatically Tweets that indicate e-cigarette use and/or e-cigarette use for smoking cessation. We build and define both datasets and compare performance of 4 state of the art classifiers and a keyword search for each task. Our results demonstrate excellent classifier performance of up to 0.90 and 0.94 area under the curve in each category. These promising initial results form the foundation for further studies to realize the ideal surveillance solution.

We present a modular approach for analyzing calcium imaging recordings of large neuronal ensembles. Our goal is to simultaneously identify the locations of the neurons, demix spatially overlapping components, and denoise and deconvolve the spiking activity from the slow dynamics of the calcium indicator. Our approach relies on a constrained nonnegative matrix factorization that expresses the spatiotemporal fluorescence activity as the product of a spatial matrix that encodes the spatial footprint of each neuron in the optical field and a temporal matrix that characterizes the calcium concentration of each neuron over time. This framework is combined with a novel constrained deconvolution approach that extracts estimates of neural activity from fluorescence traces, to create a spatiotemporal processing algorithm that requires minimal parameter tuning. We demonstrate the general applicability of our method by applying it to in vitro and in vivo multi-neuronal imaging data, whole-brain light-sheet imaging data, and dendritic imaging data.

Background
Obesity, type 2 diabetes, and non-alcoholic fatty liver disease (NAFLD) are serious health concerns, especially in Western populations. Antibiotic exposure and high-fat diet (HFD) are important and modifiable factors that may contribute to these diseases.

Methods
To investigate the relationship of antibiotic exposure with microbiome perturbations in a murine model of growth promotion, C57BL/6 mice received lifelong sub-therapeutic antibiotic treatment (STAT), or not (control), and were fed HFD starting at 13 weeks. To characterize microbiota changes caused by STAT, the V4 region of the 16S rRNA gene was examined from collected fecal samples and analyzed.

Results
In this model, which included HFD, STAT mice developed increased weight and fat mass compared to controls. Although results in males and females were not identical, insulin resistance and NAFLD were more severe in the STAT mice. Fecal microbiota from STAT mice were distinct from controls. Compared with controls, STAT exposure led to early conserved diet-independent microbiota changes indicative of an immature microbial community. Key taxa were identified as STAT-specific and several were found to be predictive of disease. Inferred network models showed topological shifts concurrent with growth promotion and suggest the presence of keystone species.

Conclusions
These studies form the basis for new models of type 2 diabetes and NAFLD that involve microbiome perturbation.

The 3D organization of chromosomes enables cells to balance the biophysical constraints of the crowded nucleus with the functional dynamics of gene regulation. Physical contacts between genes and their regulatory elements are essential for proper transcriptional control and maintenance of these interactions is critical for preventing aberrations in physiological processes that could manifest as disease states. The first insights into global nuclear organization came from imaging studies using FISH (fluorescent in-situ hybridization) analyses, which demonstrated that chromosomes occupy individual territories in the nucleus with minimal intermingling between them [1]. The development of chromosome conformation capture (3C) in which chromatin fragments in close physical proximity can be detected enabled the characterization of molecular interactions between different loci [2]. When 3C-based techniques incorporated massive parallel sequencing (such as in Hi-C) the description of molecular chromatin interactions at a genome-wide scale was finally possible [3]. Hi-C was the first unbiased approach aimed at capturing all interactions in the nucleus thereby providing a snapshot of nuclear organization at the global scale. The first Hi-C study revealed that each chromosomal territory is further divided into large domains of 5–10Mb that physically separate two compartments (A and B), which strongly correlate with active and inactive chromatin, respectively [3]. Furthermore, this study demonstrated that interactions between loci in the same compartment occur at a higher frequency than between loci in different compartments [3]. With the progressive decrease in sequencing costs, higher-resolution Hi-C revealed a new level of nuclear organization where compartments A and B can be further divided into “topologically associated domains” (TADs) [4–6]. In mammalian cells these domains range in size from a few 100kbs to 5Mbs in size (with an average of 1MB). Since they exhibit a high degree of conservation between cell types and species it was proposed that TADs represent the fundamental unit of physical organization of the genome [5].

Thanks to the confluence of genomic technology and computational developments, the possibility of network inference methods that automatically learn large comprehensive models of cellular regulation is closer than ever. This perspective focuses on enumerating the elements of computational strategies that, when coupled to appropriate experimental designs, can lead to accurate large-scale models of chromatin state and transcriptional regulatory structure and dynamics. We highlight 4 research questions that require further investigation in order to make progress in network inference: (1) using overall constraints on network structure such as sparsity, (2) use of informative priors and data integration to constrain individual model parameters, (3) estimation of latent regulatory factor activity under varying cell conditions, and (4) new methods for learning and modeling regulatory factor interactions. We conclude that methods combining advances in these 4 categories of required effort with new genomic technologies will result in biophysically motivated dynamic genome-wide regulatory network models for several of the best-studied organisms and cell types.