The Computational Modeling in Biology Network (COMBINE) is a consortium of groups involved in the development of open community standards and formats used in computational modeling in biology. COMBINE's aim is to act as a coordinator, facilitator, and resource for different standardization efforts whose domains of use cover related areas of the computational biology space. In this perspective article, we summarize COMBINE, its general organization, and the community standards and other efforts involved in it. Our goals are to help guide readers toward standards that may be suitable for their research activities, as well as to direct interested readers to relevant communities where they can best expect to receive assistance in how to develop interoperable computational models.

Smoothened (SMO) inhibitors are under clinical investigation for the treatment of several cancers. Vismodegib is approved for the treatment of locally advanced and metastatic basal cell carcinoma (BCC). Most BCC patients experience significant clinical benefit on vismodegib, but some develop resistance. Genomic analysis of tumor biopsies revealed that vismodegib resistance is associated with Hedgehog (Hh) pathway reactivation, predominantly through mutation of the drug target SMO and to a lesser extent through concurrent copy number changes in SUFU and GLI2. SMO mutations either directly impaired drug binding or activated SMO to varying levels. Furthermore, we found evidence for intra-tumor heterogeneity, suggesting that a combination of therapies targeting components at multiple levels of the Hh pathway is required to overcome resistance.

The mammalian genome harbors up to one million regulatory elements often located at great distances from their target genes. Long-range elements control genes through physical contact with promoters and can be recognized by the presence of specific histone modifications and transcription factor binding. Linking regulatory elements to specific promoters genome-wide is currently impeded by the limited resolution of high-throughput chromatin interaction assays. Here we apply a sequence capture approach to enrich Hi-C libraries for >22,000 annotated mouse promoters to identify statistically significant, long-range interactions at restriction fragment resolution, assigning long-range interacting elements to their target genes genome-wide in embryonic stem cells and fetal liver cells. The distal sites contacting active genes are enriched in active histone modifications and transcription factor occupancy, whereas inactive genes contact distal sites with repressive histone marks, demonstrating the regulatory potential of the distal elements identified. Furthermore, we find that coregulated genes cluster nonrandomly in spatial interaction networks correlated with their biological function and expression level. Interestingly, we find the strongest gene clustering in ES cells between transcription factor genes that control key developmental processes in embryogenesis. The results provide the first genome-wide catalog linking gene promoters to their long-range interacting elements and highlight the complex spatial regulatory circuitry controlling mammalian gene expression.

Contrasting phenotypes arise from similar genomes through a combination of losses, gains, co-option and modifications of inherited genomic material. Understanding the molecular basis of this phenotypic diversity is a fundamental challenge in modern evolutionary biology. Comparisons of the genes and their expression patterns underlying traits in closely related species offer an unrivaled opportunity to evaluate the extent to which genomic material is reorganized to produce novel traits. Advances in molecular methods now allow us to dissect the molecular machinery underlying phenotypic diversity in almost any organism, from single-celled entities to the most complex vertebrates. Here we discuss how comparisons of social parasites and their free-living hosts may provide unique insights into the molecular basis of phenotypic evolution. Social parasites evolve from a eusocial ancestor and are specialized to exploit the socially acquired resources of their closely-related eusocial host. Molecular comparisons of such species pairs can reveal how genomic material is re-organized in the loss of ancestral traits (i.e., of free-living traits in the parasites) and the gain of new ones (i.e., specialist traits required for a parasitic lifestyle). We define hypotheses on the molecular basis of phenotypes in the evolution of social parasitism and discuss their wider application in our understanding of the molecular basis of phenotypic diversity within the theoretical framework of phenotypic plasticity and shifting reaction norms. Currently there are no data available to test these hypotheses, and so we also provide some proof of concept data using the paper wasp social parasite/host system (Polistes sulcifer-Polistes dominula). This conceptual framework and first empirical data provide a spring-board for directing future genomic analyses on exploiting social parasites as a route to understanding the evolution of phenotypic specialization.

Regulatory T (Treg) cells suppress effector T (Teff) cells and prevent immune-mediated rejection of cancer. Much less appreciated are mechanisms by which Teff cells antagonize Treg cells. Herein, we consider how complex reciprocal interactions between Teff and Treg cells shape their population dynamics within tumors. Under states of tolerance, including during tumor escape, suppressed Teff cells support Treg cell populations through antigen-dependent provision of interleukin (IL)-2. During immune activation, Teff cells can lose this supportive capacity and directly antagonize Treg cell populations to neutralize their immunosuppressive function. While this latter state is rarely achieved spontaneously within tumors, we propose that therapeutic induction of immune activation has the potential to stably disrupt immunosuppressive population states resulting in durable cancer regression.

Myelin is essential for efficient signal transduction in the nervous system comprising of multiple proteins. The intricacies of the regulation of the formation of myelin, and its components, are not fully understood. Here, we describe the characterization of a novel myelin basic protein (Mbp) mutant mouse, mbp(jive), which spontaneously occurred in our mouse colony. These mice displayed the onset of a shaking gait before 3 weeks of age and seizure onset before 2 months of age. Due to a progressive increase of seizure intensity, mbp(jive) mice experienced premature lethality at around 3 months of age. Mbp mRNA transcript or protein was undetectable and, accordingly, genetic analysis demonstrated a homozygous loss of exons 3 to 6 of Mbp. Peripheral nerve conductance was mostly unimpaired. Additionally, we observed grave structural changes in white matter predominant structures were detected by T1, T2 and diffusion weighted magnetic resonance imaging. We additionally observed that Mbp-deficiency results in an upregulation of Qkl, Mag and Cnp, suggestive of a regulatory feedback mechanism whereby compensatory increases in Qkl have downstream effects on Mag and Cnp. Further research will clarify the role and specifications of this myelin feedback loop, as well as determine its potential role in therapeutic strategies for demyelinating disorders.

Post-transcriptional regulation of mRNA by the RNA-binding protein HuR (encoded by Elavl1) is required in B cells for the germinal center reaction and for the production of class-switched antibodies in response to thymus-independent antigens. Transcriptome-wide examination of RNA isoforms and their abundance and translation in HuR-deficient B cells, together with direct measurements of HuR-RNA interactions, revealed that HuR-dependent splicing of mRNA affected hundreds of transcripts, including that encoding dihydrolipoamide S-succinyltransferase (Dlst), a subunit of the 2-oxoglutarate dehydrogenase (α-KGDH) complex. In the absence of HuR, defective mitochondrial metabolism resulted in large amounts of reactive oxygen species and B cell death. Our study shows how post-transcriptional processes control the balance of energy metabolism required for the proliferation and differentiation of B cells.

Multiple regulatory elements distant from their targets on the linear genome can influence the expression of a single gene through chromatin looping. Chromosome conformation capture implemented in Hi-C allows for genome-wide agnostic characterization of chromatin contacts. However, detection of functional enhancer-promoter interactions is precluded by its effective resolution that is determined by both restriction fragmentation and sensitivity of the experiment. Here we develop a capture Hi-C (cHi-C) approach to allow an agnostic characterization of these physical interactions on a genome-wide scale. Single-nucleotide polymorphisms associated with complex diseases often reside within regulatory elements and exert effects through long-range regulation of gene expression. Applying this cHi-C approach to 14 colorectal cancer risk loci allows us to identify key long-range chromatin interactions in cis and trans involving these loci.

Immunologists studying the relationship between nutrition and immunological function face many challenges. We discuss here some of the historical skepticism with which nutritional research has often been faced and the complexities that need to be overcome in order to provide meaningful mechanistic insights.

Enhancers regulate spatiotemporal gene expression and impart cell-specific transcriptional outputs that drive cell identity. Super-enhancers (SEs), also known as stretch-enhancers, are a subset of enhancers especially important for genes associated with cell identity and genetic risk of disease. CD4(+) T cells are critical for host defence and autoimmunity. Here we analysed maps of mouse T-cell SEs as a non-biased means of identifying key regulatory nodes involved in cell specification. We found that cytokines and cytokine receptors were the dominant class of genes exhibiting SE architecture in T cells. Nonetheless, the locus encoding Bach2, a key negative regulator of effector differentiation, emerged as the most prominent T-cell SE, revealing a network in which SE-associated genes critical for T-cell biology are repressed by BACH2. Disease-associated single-nucleotide polymorphisms for immune-mediated disorders, including rheumatoid arthritis, were highly enriched for T-cell SEs versus typical enhancers or SEs in other cell lineages. Intriguingly, treatment of T cells with the Janus kinase (JAK) inhibitor tofacitinib disproportionately altered the expression of rheumatoid arthritis risk genes with SE structures. Together, these results indicate that genes with SE architecture in T cells encompass a variety of cytokines and cytokine receptors but are controlled by a 'guardian' transcription factor, itself endowed with an SE. Thus, enumeration of SEs allows the unbiased determination of key regulatory nodes in T cells, which are preferentially modulated by pharmacological intervention.

Post-transcriptional mRNA regulation by RNA binding proteins (RBPs) associated with AU-rich elements (AREs) present in the 3' untranslated region (3'UTR) of specific mRNAs modulates transcript stability and translation in eukaryotic cells. Here we have functionally characterised the importance of the AREs present within the Bcl2 3'UTR in order to maintain Bcl2 expression. Gene targeting deletion of 300 nucleotides of the Bcl2 3'UTR rich in AREs diminishes Bcl2 mRNA stability and protein levels in primary B cells, decreasing cell lifespan. Generation of chimeric mice indicates that Bcl2-ARE∆/∆ B cells have an intrinsic competitive disadvantage compared to wild type cells. Biochemical assays and predictions using a bioinformatics approach show that several RBPs bind to the Bcl2 AREs, including AUF1 and HuR proteins. Altogether, association of RBPs to Bcl2 AREs contributes to Bcl2 protein expression by stabilizing Bcl2 mRNA and promotes B cell maintenance.

Drosophila phototransduction is mediated via a G-protein-coupled PLC cascade. Recent evidence, including the demonstration that light evokes rapid contractions of the photoreceptors, suggested that the light-sensitive channels (TRP and TRPL) may be mechanically gated, together with protons released by PLC-mediated PIP2 hydrolysis. If mechanical gating is involved we predicted that the response to light should be influenced by altering the physical properties of the membrane. To achieve this, we used diet to manipulate the degree of saturation of membrane phospholipids. In flies reared on a yeast diet, lacking polyunsaturated fatty acids (PUFAs), mass spectrometry showed that the proportion of polyunsaturated phospholipids was sevenfold reduced (from 38 to ∼5%) but rescued by adding a single species of PUFA (linolenic or linoleic acid) to the diet. Photoreceptors from yeast-reared flies showed a 2- to 3-fold increase in latency and time to peak of the light response, without affecting quantum bump waveform. In the absence of Ca(2+) influx or in trp mutants expressing only TRPL channels, sensitivity to light was reduced up to ∼10-fold by the yeast diet, and essentially abolished in hypomorphic G-protein mutants (Gαq). PLC activity appeared little affected by the yeast diet; however, light-induced contractions measured by atomic force microscopy or the activation of ectopic mechanosensitive gramicidin channels were also slowed ∼2-fold. The results are consistent with mechanosensitive gating and provide a striking example of how dietary fatty acids can profoundly influence sensory performance in a classical G-protein-coupled signaling cascade.

Gene variants that disrupt TCR signaling can cause severe immune deficiency, yet less disruptive variants are sometimes associated with immune pathology. Null mutations of the gene encoding the scaffold protein Src homology 2 domain-containing leukocyte protein of 76 kDa (SLP-76), for example, cause an arrest of T cell positive selection, whereas a synthetic membrane-targeted allele allows limited positive selection but is associated with proinflammatory cytokine production and autoantibodies. Whether these and other enigmatic outcomes are due to a biochemical uncoupling of tolerogenic signaling, or simply a quantitative reduction of protein activity, remains to be determined. In this study we describe a splice variant of Lcp2 that reduced the amount of wild-type SLP-76 protein by ~90%, disrupting immunogenic and tolerogenic pathways to different degrees. Mutant mice produced excessive amounts of proinflammatory cytokines, autoantibodies, and IgE, revealing that simple quantitative reductions of SLP-76 were sufficient to trigger immune dysregulation. This allele reveals a dose-sensitive threshold for SLP-76 in the balance of immunity and immune dysregulation, a common disturbance of atypical clinical immune deficiencies.

Allogeneic hematopoietic stem cell transplantation (HSCT) can treat certain hematologic malignancies due to the graft versus leukemia (GvL) effect but is complicated by graft versus host disease (GvHD). Expression of the p110δ catalytic subunit of the phosphoinositide 3-kinase pathway is restricted to leukocytes, where it regulates proliferation, migration, and cytokine production. Here, in a mouse model of fully mismatched hematopoietic cell transplantation (HCT), we show that genetic inactivation of p110δ in T cells leads to milder GvHD, whereas GvL is preserved. Inactivation of p110δ in human lymphocytes reduced T cell allorecognition. We demonstrate that both allostimulation and granzyme B expression were dependent on p110δ in naive T cells, which are the main mediators of GvHD, whereas memory T cells were unaffected. Strikingly, p110δ is not mandatory for either naive or memory T cells to mediate GvL. Therefore, immunomodulation of selective naive T cell functions by p110δ inactivation improves the outcome of allogeneic HSCT.

Behaviours of complex biomolecular systems are often irreducible to the elementary properties of their individual components. Explanatory and predictive mathematical models are therefore useful for fully understanding and precisely engineering cellular functions. The development and analyses of these models require their adaptation to the problems that need to be solved and the type and amount of available genetic or molecular data. Quantitative and logic modelling are among the main methods currently used to model molecular and gene networks. Each approach comes with inherent advantages and weaknesses. Recent developments show that hybrid approaches will become essential for further progress in synthetic biology and in the development of virtual organisms.

Cellular senescence has been implicated in tumor suppression, development, and aging and is accompanied by large-scale chromatin rearrangements, forming senescence-associated heterochromatic foci (SAHF). However, how the chromatin is reorganized during SAHF formation is poorly understood. Furthermore, heterochromatin formation in senescence appears to contrast with loss of heterochromatin in Hutchinson-Gilford progeria. We mapped architectural changes in genome organization in cellular senescence using Hi-C. Unexpectedly, we find a dramatic sequence- and lamin-dependent loss of local interactions in heterochromatin. This change in local connectivity resolves the paradox of opposing chromatin changes in senescence and progeria. In addition, we observe a senescence-specific spatial clustering of heterochromatic regions, suggesting a unique second step required for SAHF formation. Comparison of embryonic stem cells (ESCs), somatic cells, and senescent cells shows a unidirectional loss in local chromatin connectivity, suggesting that senescence is an endpoint of the continuous nuclear remodelling process during differentiation.

On the Semantic Web, in life sciences in particular, data is often distributed via multiple resources. Each of these sources is likely to use their own IRI (International Resource Identifier) for conceptually the same resource or database record. The lack of correspondence between identifiers introduces a barrier when executing federated SPARQL queries across life science data.

Hematopoietic stem cells (HSCs) possess the ability to generate all hematopoietic cell types and to self-renew over long periods, but the mechanisms that regulate their unique properties are incompletely understood. Herein, we show that homozygous deletion of the miR-29a/b-1 bicistron results in decreased numbers of hematopoietic stem and progenitor cells (HSPCs), decreased HSC self-renewal, and increased HSC cell cycling and apoptosis. The HSPC phenotype is specifically due to loss of miR-29a, because miR-29b expression is unaltered in miR-29a/b-1-null HSCs, and only ectopic expression of miR-29a restores HSPC function both in vitro and in vivo. HSCs lacking miR-29a/b-1 exhibit widespread transcriptional dysregulation and adopt gene expression patterns similar to normal committed progenitors. A number of predicted miR-29 target genes, including Dnmt3a, are significantly upregulated in miR-29a/b-1-null HSCs. The loss of negative regulation of Dnmt3a by miR-29a is a major contributor to the miR-29a/b-1-null HSPC phenotype, as both in vitro Dnmt3a short hairpin RNA knockdown assays and a genetic haploinsufficiency model of Dnmt3a restored the frequency and long-term reconstitution capacity of HSCs from miR-29a/b-1-deficient mice. Overall, these data demonstrate that miR-29a is critical for maintaining HSC function through its negative regulation of Dnmt3a.

Modulation of the properties of AMPA receptors at the post-synaptic membrane is one of the main suggested mechanisms underlying fast synaptic transmission in the central nervous system of vertebrates. Electrophysiological recordings of single channels stimulated with agonists showed that both recombinant and native AMPA receptors visit multiple conductance states in an agonist concentration dependent manner. We propose an allosteric model of the multiple conductance states based on concerted conformational transitions of the four subunits, as an iris diaphragm. Our model predicts that the thermodynamic behaviour of the conductance states upon full and partial agonist stimulations can be described with increased affinity of receptors as they progress to higher conductance states. The model also predicts the existence of AMPA receptors in non-liganded conductive substates. However, the probability of spontaneous openings decreases with increasing conductances. Finally, we predict that the large conductance states are stabilized within the rise phase of a whole-cell EPSC in glutamatergic hippocampal neurons. Our model provides a mechanistic link between ligand concentration and conductance states that can explain thermodynamic and kinetic features of AMPA receptor gating.

The thymus is the organ devoted to T-cell production. The thymus undergoes multiple rounds of atrophy and redevelopment before degenerating with age in a process known as involution. This process is poorly understood, despite the influence the phenomenon has on peripheral T-cell numbers. Here we have investigated the FVB/N mouse strain, which displays premature thymic involution. We find multiple architectural and cellular features that precede thymic involution, including disruption of the epithelial-endothelial relationship and a progressive loss of pro-T cells. The architectural features, reminiscent of the human thymus, are intrinsic to the non-hematopoietic compartment and are neither necessary nor sufficient for thymic involution. By contrast, the loss of pro-T cells is intrinsic to the hematopoietic compartment, and is sufficient to drive premature involution. These results identify pro-T-cell loss as the main driver of premature thymic involution, and highlight the plasticity of the thymic stroma, capable of maintaining function across diverse inter-strain architectures. This article is protected by copyright. All rights reserved.

DNA methylomes are extensively reprogrammed during mouse pre-implantation and early germ cell development. The main feature of this reprogramming is a genome-wide decrease in 5-methylcytosine (5mC). Standard high-resolution single-stranded bisulfite sequencing techniques do not allow discrimination of the underlying passive (replication-dependent) or active enzymatic mechanisms of 5mC loss. We approached this problem by generating high-resolution deep hairpin bisulfite sequencing (DHBS) maps, allowing us to follow the patterns of symmetric DNA methylation at CpGs dyads on both DNA strands over single replications.