Search

Thermal proteome profiling (TPP) is a proteome wide technology that enables unbiased detection of protein drug interactions as well as changes in post-translational state of proteins between different biological conditions. Statistical analysis of temperature range TPP (TPP-TR) datasets relies on comparing protein melting curves, describing the amount of non-denatured proteins as a function of temperature, between different conditions (e.g. presence or absence of a drug). However, state-of-the-art models are restricted to sigmoidal melting behaviours while unconventional melting curves, representing up to 50% of TPP-TR datasets, have recently been shown to carry important biological information. We present a novel statistical framework, based on hierarchical Gaussian process models and named GPMelt, to make TPP-TR datasets analysis unbiased with respect to the melting profiles of proteins. GPMelt scales to multiple conditions, and extension of the model to deeper hierarchies (i.e. with additional sub-levels) allows to deal with complex TPP-TR protocols. Collectively, our statistical framework extends the analysis of TPP-TR datasets for both protein and peptide level melting curves, offering access to thousands of previously excluded melting curves and thus substantially increasing the coverage and the ability of TPP to uncover new biology. Author summary: Proteins interactions with other proteins, nucleic acids or metabolites, are key to all biological processes. Being able to detect these interactions is essential to understand biological systems. Thermal proteome profiling is a proteome-wide biological assay able to capture these interactions. It consists in analysing the effect of heat treatment on proteins. Indeed, proteins, under physiological conditions, are folded. This folding gets disrupted as the temperature increases. The way this unfolding happens, called the melting profile of the protein, informs on the interactions of proteins. For example, the interaction of a protein with another protein can increase (thermally stabilise) or decrease (thermally destabilise) the temperature at which this protein starts unfolding. In this work, we present a new statistical method, named GPMelt, to analyse these melting profiles. Notably, GPMelt allows to analyse any melting profiles, independently of their shapes. The proposed improvements over previously published methods allow to investigate more robustly the melting profiles of more proteins, hence increasing the ability of thermal proteome profiling assays to discover new protein interactions. We anticipate that these advancements will aid in unravelling complex biological phenomena. [ABSTRACT FROM AUTHOR]

Background: Previous studies which applied machine learning on multiplex component‐resolved diagnostics arrays identified clusters of allergen components which are biologically plausible and reflect the sources of allergenic proteins and their structural homogeneity. Sensitization to different clusters is associated with different clinical outcomes. Objective: To investigate whether within different allergen component sensitization clusters, the internal within‐cluster sensitization structure, including the number of c‐sIgE responses and their distinct patterns, alters the risk of clinical expression of symptoms. Methods: In a previous analysis in a population‐based birth cohort, by clustering component‐specific (c‐s)IgEs, we derived allergen component clusters from infancy to adolescence. In the current analysis, we defined each subject's within‐cluster sensitization structure which captured the total number of c‐sIgE responses in each cluster and intra‐cluster sensitization patterns. Associations between within‐cluster sensitization patterns and clinical outcomes (asthma and rhinitis) in early‐school age and adolescence were examined using logistic regression and binomial generalized additive models. Results: Intra‐cluster sensitization patterns revealed specific associations with asthma and rhinitis (both contemporaneously and longitudinally) that were previously unseen using binary sensitization to clusters. A more detailed description of the subjects' within‐cluster c‐sIgE responses in terms of the number of positive c‐sIgEs and unique sensitization patterns added new information relevant to allergic diseases, both for diagnostic and prognostic purposes. For example, the increase in the number of within‐cluster positive c‐sIgEs at age 5 years was correlated with the increase in prevalence of asthma at ages 5 and 16 years, with the correlations being stronger in the prediction context (e.g. for the largest 'Broad' component cluster, contemporaneous: r =.28, p =.012; r =.22, p =.043; longitudinal: r =.36, p =.004; r =.27, p =.04). Conclusion: Among sensitized individuals, a more detailed description of within‐cluster c‐sIgE responses in terms of the number of positive c‐sIgE responses and distinct sensitization patterns, adds potentially important information relevant to allergic diseases. [ABSTRACT FROM AUTHOR]

BackgroundSeveral high-throughput technologies can measure in parallel the abundance of many mRNA transcripts within a sample. These include the widely-used microarray as well as the more recently developed methods based on sequence tag abundances such as the Massively Parallel Signature Sequencing (MPSS) technology. A comparison of microarray and MPSS technologies can help to establish the metrics for data comparisons across these technology platforms and determine some of the factors affecting the measurement of mRNA abundances using different platforms.ResultsWe compared transcript abundance (gene expression) measurement data obtained using Affymetrix and Agilent microarrays with MPSS data. All three technologies were used to analyze the same set of mRNA samples; these samples were extracted from various wild type Arabidopsis thaliana tissues and floral mutants. We calculated correlations and used clustering methodology to compare the normalized expression data and expression ratios across samples and technologies. Abundance expression measurements were more similar between different samples measured by the same technology than between the same sample measured by different technologies. However, when expression ratios were employed, samples measured by different technologies were found to cluster together more frequently than with abundance expression levels.Furthermore, the two microarray technologies were more consistent with each other than with MPSS. We also investigated probe-position effects on Affymetrix data and tag-position effects in MPSS. We found a similar impact on Affymetrix and MPSS measurements, which suggests that these effects were more likely a characteristic of the RNA sample rather than technology-specific biases.ConclusionComparisons of transcript expression ratios showed greater consistency across platforms than measurements of transcript abundance. In addition, for measurements based on abundances, technology differences can mask the impact of biological differences between samples and tissues.

BackgroundClustering is an important analysis performed on microarray gene expression data since it groups genes which have similar expression patterns and enables the exploration of unknown gene functions. Microarray experiments are associated with many sources of experimental and biological variation and the resulting gene expression data are therefore very noisy. Many heuristic and model-based clustering approaches have been developed to cluster this noisy data. However, few of them include consideration of probe-level measurement error which provides rich information about technical variability.ResultsWe augment a standard model-based clustering method to incorporate probe-level measurement error. Using probe-level measurements from a recently developed Affymetrix probe-level model, multi-mgMOS, we include the probe-level measurement error directly into the standard Gaussian mixture model. Our augmented model is shown to provide improved clustering performance on simulated datasets and a real mouse time-course dataset.ConclusionThe performance of model-based clustering of gene expression data is improved by including probe-level measurement error and more biologically meaningful clustering results are obtained.

The glucocorticoid receptor (GR) is a major drug target in inflammatory disease. However, chronic glucocorticoid (GC) treatment leads to disordered energy metabolism, including increased weight gain, adiposity, and hepatosteatosis--all programs modulated by the circadian clock. We demonstrated that while antiinflammatory GC actions were maintained irrespective of dosing time, the liver was significantly more GC sensitive during the day. Temporal segregation of GC action was underpinned by a physical interaction of GR with the circadian transcription factor REVERBa and co-binding with liver- specific hepatocyte nuclear transcription factors (HNFs) on chromatin. REVERBa promoted efficient GR recruitment to chromatin during the day, acting in part by maintaining histone acetylation, with REVERBa- dependent GC responses providing segregation of carbohydrate and lipid metabolism. Importantly, deletion of Reverba inverted circadian liver GC sensitivity and protected mice from hepatosteatosis induced by chronic GC administration. Our results reveal a mechanism by which the circadian clock acts through REVERBa in liver on elements bound by HNF4A/HNF6 to direct GR action on energy metabolism., Introduction Glucocorticoids (GCs; cortisol in humans, corticosterone in rodents) are critical regulators of energy metabolism and immunity. Their secretion by the adrenal gland follows a circadian pattern, with serum concentrations [...]

Motivation There exists a range of different quantification frameworks to estimate the synergistic effect of drug combinations. The diversity and disagreement in estimates make it challenging to determine which combinations from a large drug screening should be proceeded with. Furthermore, the lack of accurate uncertainty quantification for those estimates precludes the choice of optimal drug combinations based on the most favourable synergistic effect. Results In this work, we propose SynBa, a flexible Bayesian approach to estimate the uncertainty of the synergistic efficacy and potency of drug combinations, so that actionable decisions can be derived from the model outputs. The actionability is enabled by incorporating the Hill equation into SynBa, so that the parameters representing the potency and the efficacy can be preserved. Existing knowledge may be conveniently inserted due to the flexibility of the prior, as shown by the empirical Beta prior defined for the normalized maximal inhibition. Through experiments on large combination screenings and comparison against benchmark methods, we show that SynBa provides improved accuracy of dose–response predictions and better-calibrated uncertainty estimation for the parameters and the predictions. Availability and implementation The code for SynBa is available at https://github.com/HaotingZhang1/SynBa. The datasets are publicly available (DOI of DREAM: 10.7303/syn4231880 ; DOI of the NCI-ALMANAC subset: 1 0.5281/zenodo.4135059). [ABSTRACT FROM AUTHOR]

The continuing and rapid increase in the number of fully sequenced genomes is creating new opportunities for comparative studies. However, although many genomic databases store data from multiple organisms, for the most part they provide limited support for comparative genomics. We argue that refocusing genomic data management to provide more direct support for comparative studies enables systematic identification of important relationships between species, thereby increasing the value that can be obtained from sequenced genomes. The principal result of the paper is a methodology, in which comparative analyses are constructed over a foundation based on sequence clusters and evolutionary relationships. This methodology has been applied in a systematic study of the fungi, and we describe how comparative analyses have been implemented as an analysis library over the e-Fungi data warehouse.

Regulation of mRNA degradation is critical for a diverse array of cellular processes and developmental cell fate decisions. Many methods for determining mRNA half-lives rely on transcriptional inhibition or metabolic labelling. Here, we use a non-invasive method for estimating half-lives for hundreds of mRNAs in the early Drosophila embryo. This approach uses the intronic and exonic reads from a total RNA-seq time series and Gaussian process regression to model the dynamics of premature and mature mRNAs. We show how regulation of mRNA stability is used to establish a range of mature mRNA dynamics during embryogenesis, despite shared transcription profiles. Using single-molecule imaging, we provide evidence that, for the mRNAs tested, there is a correlation between short half-life and mRNA association with P-bodies. Moreover, we detect an enrichment of mRNA 3′ ends in P-bodies in the early embryo, consistent with 5′ to 3′ degradation occurring in P-bodies for at least a subset of mRNAs. We discuss our findings in relation to recently published data suggesting that the primary function of P-bodies in other biological contexts is mRNA storage. Regulation of mRNA degradation is critical for a diverse array of cellular processes and developmental cell fate decisions. This study uses modelling of total RNA-seq time series data across early Drosophila embryogenesis to provide estimates of mRNA half-lives; single-molecule imaging reveals that unstable mRNAs are more enriched in P-bodies, as are 3′ mRNA fragments. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRs) have an important role in tuning dynamic gene expression. However, the mechanism by which they are quantitatively controlled is unknown. We show that the amount of mature miR-9, a key regulator of neuronal development, increases during zebrafish neurogenesis in a sharp stepwise manner. We characterize the spatiotemporal profile of seven distinct microRNA primary transcripts (pri-mir)-9s that produce the same mature miR-9 and show that they are sequentially expressed during hindbrain neurogenesis. Expression of late-onset pri-mir-9-1 is added on to, rather than replacing, the expression of early onset pri-mir-9-4 and -9-5 in single cells. CRISPR/Cas9 mutation of the late-onset pri-mir-9-1 prevents the developmental increase of mature miR-9, reduces late neuronal differentiation and fails to downregulate Her6 at late stages. Mathematical modelling shows that an adaptive network containing Her6 is insensitive to linear increases in miR-9 but responds to stepwise increases of miR-9. We suggest that a sharp stepwise increase of mature miR-9 is created by sequential and additive temporal activation of distinct loci. This may be a strategy to overcome adaptation and facilitate a transition of Her6 to a new dynamic regime or steady state. [ABSTRACT FROM AUTHOR]

Despite good evidence of impaired innate antiviral responses in asthma, trials of inhaled interferon-β given during exacerbations showed only modest benefits in moderate/severe asthma. Using human experimental rhinovirus infection, we observe robust in vivo induction of bronchial epithelial interferon response genes 4 days after virus inoculation in 25 subjects with asthma but not 11 control subjects. This signature correlated with virus loads and lower respiratory symptoms. Our data indicate that the in vivo innate antiviral response is dysregulated in asthma and open up the potential that prophylactic rather than therapeutic interferon therapy may have greater clinical benefit. [ABSTRACT FROM AUTHOR]

Damage to the lung epithelium is a unifying feature of disease caused by the saprophytic fungus Aspergillus fumigatus. However, the mechanistic basis and the regulatory control of such damage is poorly characterized. Previous studies have identified A. fumigatus mediated pathogenesis as occurring at early (= 16 hours) or late (>16 hours) phases of the fungal interaction with epithelial cells, and respectively involve direct contact with the host cell or the action of soluble factors produced by mature fungal hyphae. Both early and late phases of epithelial damage have been shown to be subject to genetic regulation by the pH-responsive transcription factor PacC. This study sought to determine whether other transcriptional regulators play a role in modulating epithelial damage. In particular, whether the early and late phases of epithelial damage are governed by same or distinct regulators. Furthermore, whether processes such as spore uptake and hyphal adhesion, that have previously been documented to promote epithelial damage, are governed by the same cohorts of epithelial regulators. Using 479 strains from the recently constructed library of A. fumigatus transcription factor null mutants, two high-throughput screens assessing epithelial cell detachment and epithelial cell lysis were conducted. A total of 17 transcription factor mutants were found to exhibit reproducible deficits in epithelial damage causation. Of these, 10 mutants were defective in causing early phase damage via epithelial detachment and 8 mutants were defective in causing late phase damage via epithelial lysis. Remarkably only one transcription factor, PacC, was required for causation of both phases of epithelial damage. The 17 mutants exhibited varied and often unique phenotypic profiles with respect to fitness, epithelial adhesion, cell wall defects, and rates of spore uptake by epithelial cells. Strikingly, 9 out of 10 mutants deficient in causing early phase damage also exhibited reduced rates of hyphal extension, and culture supernatants of 7 out of 8 mutants deficient in late phase damage were significantly less cytotoxic. Our study delivers the first high-level overview of A. fumigatus regulatory genes governing lung epithelial damage, suggesting highly coordinated genetic orchestration of host-damaging activities that govern epithelial damage in both space and time. [ABSTRACT FROM AUTHOR]

Chronic inflammation underpins many human diseases. Morbidity and mortality associated with chronic inflammation are often mediated through metabolic dysfunction. Inflammatory and metabolic processes vary through circadian time, suggesting an important temporal crosstalk between these systems. Using an established mouse model of rheumatoid arthritis, we show that chronic inflammatory arthritis results in rhythmic joint inflammation and drives major changes in muscle and liver energy metabolism and rhythmic gene expression. Transcriptional and phosphoproteomic analyses revealed alterations in lipid metabolism and mitochondrial function associated with increased EGFR-JAK-STAT3 signaling. Metabolomic analyses confirmed rhythmic metabolic rewiring with impaired β-oxidation and lipid handling and revealed a pronounced shunt toward sphingolipid and ceramide accumulation. The arthritis-related production of ceramides was most pronounced during the day, which is the time of peak inflammation and increased reliance on fatty acid oxidation. Thus, our data demonstrate that localized joint inflammation drives a time-of-day–dependent build-up of bioactive lipid species driven by rhythmic inflammation and altered EGFR-STAT signaling. [ABSTRACT FROM AUTHOR]

Motivation The MS2-MCP (MS2 coat protein) live imaging system allows for visualization of transcription dynamics through the introduction of hairpin stem-loop sequences into a gene. A fluorescent signal at the site of nascent transcription in the nucleus quantifies mRNA production. Computational modelling can be used to infer the promoter states along with the kinetic parameters governing transcription, such as promoter switching frequency and polymerase loading rate. However, modelling of the fluorescent trace presents a challenge due its persistence; the observed fluorescence at a given time point depends on both current and previous promoter states. A compound state Hidden Markov Model (cpHMM) was recently introduced to allow inference of promoter activity from MS2-MCP data. However, the computational time for inference scales exponentially with gene length and the cpHMM is therefore not currently practical for application to many eukaryotic genes. Results We present a scalable implementation of the cpHMM for fast inference of promoter activity and transcriptional kinetic parameters. This new method can model genes of arbitrary length through the use of a time-adaptive truncated compound state space. The truncated state space provides a good approximation to the full state space by retaining the most likely set of states at each time during the forward pass of the algorithm. Testing on MS2-MCP fluorescent data collected from early Drosophila melanogaster embryos indicates that the method provides accurate inference of kinetic parameters within a computationally feasible timeframe. The inferred promoter traces generated by the model can also be used to infer single-cell transcriptional parameters. Availability and implementation Python implementation is available at https://github.com/ManchesterBioinference/burstInfer , along with code to reproduce the examples presented here. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

Motivation The negative binomial distribution has been shown to be a good model for counts data from both bulk and single-cell RNA-sequencing (RNA-seq). Gaussian process (GP) regression provides a useful non-parametric approach for modelling temporal or spatial changes in gene expression. However, currently available GP regression methods that implement negative binomial likelihood models do not scale to the increasingly large datasets being produced by single-cell and spatial transcriptomics. Results The GPcounts package implements GP regression methods for modelling counts data using a negative binomial likelihood function. Computational efficiency is achieved through the use of variational Bayesian inference. The GP function models changes in the mean of the negative binomial likelihood through a logarithmic link function and the dispersion parameter is fitted by maximum likelihood. We validate the method on simulated time course data, showing better performance to identify changes in over-dispersed counts data than methods based on Gaussian or Poisson likelihoods. To demonstrate temporal inference, we apply GPcounts to single-cell RNA-seq datasets after pseudotime and branching inference. To demonstrate spatial inference, we apply GPcounts to data from the mouse olfactory bulb to identify spatially variable genes and compare to two published GP methods. We also provide the option of modelling additional dropout using a zero-inflated negative binomial. Our results show that GPcounts can be used to model temporal and spatial counts data in cases where simpler Gaussian and Poisson likelihoods are unrealistic. Availability and implementation GPcounts is implemented using the GPflow library in Python and is available at https://github.com/ManchesterBioinference/GPcounts along with the data, code and notebooks required to reproduce the results presented here. The version used for this paper is archived at https://doi.org/10.5281/zenodo.5027066. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]