Fomes weberianus, 50 years of taxonomic confusion: lectotypification and taxonomic notes.

Fomes weberianus Bres. & Henn. ex Sacc. is currently the basionym of two very distinct polypores (Basidiomycota), Ganoderma weberianum (Polyporales) and Phylloporia weberiana (Hymenochaetales). This fact has led to almost fifty years of taxonomic confusion. Fomes weberianus was first lectotypified by Steyaert, who accepted the species as G. weberianum. However, studies of Weber's original material in B, duplicate material in S, the protologue, and early interpretations of the name have shown that Steyaert's choice conflicts with the protologue and early interpretations, and that his interpretation as a species of Ganoderma is erroneous. A new lectotype was designated and the species was re-described under the correct interpretation Phylloporia weberiana.

Clustering long-term health conditions among 67728 people with multimorbidity using electronic health records in Scotland.

There is still limited understanding of how chronic conditions co-occur in patients with multimorbidity and what are the consequences for patients and the health care system. Most reported clusters of conditions have not considered the demographic characteristics of these patients during the clustering process. The study used data for all registered patients that were resident in Fife or Tayside, Scotland and aged 25 years or more on 1st January 2000 and who were followed up until 31st December 2018. We used linked demographic information, and secondary care electronic health records from 1st January 2000. Individuals with at least two of the 31 Elixhauser Comorbidity Index conditions were identified as having multimorbidity. Market basket analysis was used to cluster the conditions for the whole population and then repeatedly stratified by age, sex and deprivation. 318,235 individuals were included in the analysis, with 67,728 (21·3%) having multimorbidity. We identified five distinct clusters of conditions in the population with multimorbidity: alcohol misuse, cancer, obesity, renal failure, and heart failure. Clusters of long-term conditions differed by age, sex and socioeconomic deprivation, with some clusters not present for specific strata and others including additional conditions. These findings highlight the importance of considering demographic factors during both clustering analysis and intervention planning for individuals with multiple long-term conditions. By taking these factors into account, the healthcare system may be better equipped to develop tailored interventions that address the needs of complex patients.

Transcriptional, epigenetic and metabolic signatures in cardiometabolic syndrome defined by extreme phenotypes.

BACKGROUND: This work is aimed at improving the understanding of cardiometabolic syndrome pathophysiology and its relationship with thrombosis by generating a multi-omic disease signature. METHODS/RESULTS: We combined classic plasma biochemistry and plasma biomarkers with the transcriptional and epigenetic characterisation of cell types involved in thrombosis, obtained from two extreme phenotype groups (morbidly obese and lipodystrophy) and lean individuals to identify the molecular mechanisms at play, highlighting patterns of abnormal activation in innate immune phagocytic cells. Our analyses showed that extreme phenotype groups could be distinguished from lean individuals, and from each other, across all data layers. The characterisation of the same obese group, 6 months after bariatric surgery, revealed the loss of the abnormal activation of innate immune cells previously observed. However, rather than reverting to the gene expression landscape of lean individuals, this occurred via the establishment of novel gene expression landscapes. NETosis and its control mechanisms emerge amongst the pathways that show an improvement after surgical intervention. CONCLUSIONS: We showed that the morbidly obese and lipodystrophy groups, despite some differences, shared a common cardiometabolic syndrome signature. We also showed that this could be used to discriminate, amongst the normal population, those individuals with a higher likelihood of presenting with the disease, even when not displaying the classic features.

Tailored Bayes: a risk modeling framework under unequal misclassification costs.

Risk prediction models are a crucial tool in healthcare. Risk prediction models with a binary outcome (i.e., binary classification models) are often constructed using methodology which assumes the costs of different classification errors are equal. In many healthcare applications, this assumption is not valid, and the differences between misclassification costs can be quite large. For instance, in a diagnostic setting, the cost of misdiagnosing a person with a life-threatening disease as healthy may be larger than the cost of misdiagnosing a healthy person as a patient. In this article, we present Tailored Bayes (TB), a novel Bayesian inference framework which "tailors" model fitting to optimize predictive performance with respect to unbalanced misclassification costs. We use simulation studies to showcase when TB is expected to outperform standard Bayesian methods in the context of logistic regression. We then apply TB to three real-world applications, a cardiac surgery, a breast cancer prognostication task, and a breast cancer tumor classification task and demonstrate the improvement in predictive performance over standard methods.

The RNA landscape of the human placenta in health and disease.

The placenta is the interface between mother and fetus and inadequate function contributes to short and long-term ill-health. The placenta is absent from most large-scale RNA-Seq datasets. We therefore analyze long and small RNAs (~101 and 20 million reads per sample respectively) from 302 human placentas, including 94 cases of preeclampsia (PE) and 56 cases of fetal growth restriction (FGR). The placental transcriptome has the seventh lowest complexity of 50 human tissues: 271 genes account for 50% of all reads. We identify multiple circular RNAs and validate 6 of these by Sanger sequencing across the back-splice junction. Using large-scale mass spectrometry datasets, we find strong evidence of peptides produced by translation of two circular RNAs. We also identify novel piRNAs which are clustered on Chr1 and Chr14. PE and FGR are associated with multiple and overlapping differences in mRNA, lincRNA and circRNA but fewer consistent differences in small RNAs. Of the three protein coding genes differentially expressed in both PE and FGR, one encodes a secreted protein FSTL3 (follistatin-like 3). Elevated serum levels of FSTL3 in pregnant women are predictive of subsequent PE and FGR. To aid visualization of our placenta transcriptome data, we develop a web application ( https://www.obgyn.cam.ac.uk/placentome/ ).

Unambiguous identification of fungi: where do we stand and how accurate and precise is fungal DNA barcoding?

True fungi (Fungi) and fungus-like organisms (e.g. Mycetozoa, Oomycota) constitute the second largest group of organisms based on global richness estimates, with around 3 million predicted species. Compared to plants and animals, fungi have simple body plans with often morphologically and ecologically obscure structures. This poses challenges for accurate and precise identifications. Here we provide a conceptual framework for the identification of fungi, encouraging the approach of integrative (polyphasic) taxonomy for species delimitation, i.e. the combination of genealogy (phylogeny), phenotype (including autecology), and reproductive biology (when feasible). This allows objective evaluation of diagnostic characters, either phenotypic or molecular or both. Verification of identifications is crucial but often neglected. Because of clade-specific evolutionary histories, there is currently no single tool for the identification of fungi, although DNA barcoding using the internal transcribed spacer (ITS) remains a first diagnosis, particularly in metabarcoding studies. Secondary DNA barcodes are increasingly implemented for groups where ITS does not provide sufficient precision. Issues of pairwise sequence similarity-based identifications and OTU clustering are discussed, and multiple sequence alignment-based phylogenetic approaches with subsequent verification are recommended as more accurate alternatives. In metabarcoding approaches, the trade-off between speed and accuracy and precision of molecular identifications must be carefully considered. Intragenomic variation of the ITS and other barcoding markers should be properly documented, as phylotype diversity is not necessarily a proxy of species richness. Important strategies to improve molecular identification of fungi are: (1) broadly document intraspecific and intragenomic variation of barcoding markers; (2) substantially expand sequence repositories, focusing on undersampled clades and missing taxa; (3) improve curation of sequence labels in primary repositories and substantially increase the number of sequences based on verified material; (4) link sequence data to digital information of voucher specimens including imagery. In parallel, technological improvements to genome sequencing offer promising alternatives to DNA barcoding in the future. Despite the prevalence of DNA-based fungal taxonomy, phenotype-based approaches remain an important strategy to catalog the global diversity of fungi and establish initial species hypotheses.

Multiple kernel learning for integrative consensus clustering of omic datasets.

MOTIVATION: Diverse applications-particularly in tumour subtyping-have demonstrated the importance of integrative clustering techniques for combining information from multiple data sources. Cluster Of Clusters Analysis (COCA) is one such approach that has been widely applied in the context of tumour subtyping. However, the properties of COCA have never been systematically explored, and its robustness to the inclusion of noisy datasets is unclear. RESULTS: We rigorously benchmark COCA, and present Kernel Learning Integrative Clustering (KLIC) as an alternative strategy. KLIC frames the challenge of combining clustering structures as a multiple kernel learning problem, in which different datasets each provide a weighted contribution to the final clustering. This allows the contribution of noisy datasets to be down-weighted relative to more informative datasets. We compare the performances of KLIC and COCA in a variety of situations through simulation studies. We also present the output of KLIC and COCA in real data applications to cancer subtyping and transcriptional module discovery. AVAILABILITY AND IMPLEMENTATION: R packages klic and coca are available on the Comprehensive R Archive Network. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Fast approximate inference for variable selection in Dirichlet process mixtures, with an application to pan-cancer proteomics.

The Dirichlet Process (DP) mixture model has become a popular choice for model-based clustering, largely because it allows the number of clusters to be inferred. The sequential updating and greedy search (SUGS) algorithm (Wang & Dunson, 2011) was proposed as a fast method for performing approximate Bayesian inference in DP mixture models, by posing clustering as a Bayesian model selection (BMS) problem and avoiding the use of computationally costly Markov chain Monte Carlo methods. Here we consider how this approach may be extended to permit variable selection for clustering, and also demonstrate the benefits of Bayesian model averaging (BMA) in place of BMS. Through an array of simulation examples and well-studied examples from cancer transcriptomics, we show that our method performs competitively with the current state-of-the-art, while also offering computational benefits. We apply our approach to reverse-phase protein array (RPPA) data from The Cancer Genome Atlas (TCGA) in order to perform a pan-cancer proteomic characterisation of 5157 tumour samples. We have implemented our approach, together with the original SUGS algorithm, in an open-source R package named sugsvarsel, which accelerates analysis by performing intensive computations in C++ and provides automated parallel processing. The R package is freely available from: https://github.com/ococrook/sugsvarsel.

Post-mortem study of the association between cardiac iron and fibrosis in transfusion dependent anaemia.

BACKGROUND: Heart failure related to cardiac siderosis remains a major cause of death in transfusion dependent anaemias. Replacement fibrosis has been reported as causative of heart failure in siderotic cardiomyopathy in historical reports, but these findings do not accord with the reversible nature of siderotic heart failure achievable with intensive iron chelation. METHODS: Ten whole human hearts (9 beta-thalassemia major, 1 sideroblastic anaemia) were examined for iron loading and fibrosis (replacement and interstitial). Five had died from heart failure, 4 had cardiac transplantation for heart failure, and 1 had no heart failure (death from a stroke). Heart samples iron content was measured using atomic emission spectroscopy. Interstitial fibrosis was quantified by computer using picrosirius red (PSR) staining and expressed as collagen volume fraction (CVF) with normal value for left ventricle <3%. RESULTS: The 9 hearts affected by heart failure had severe iron loading with very low T2* of 5.0 ± 2.0 ms (iron concentration 8.5 ± 7.0 mg/g dw) and diffuse granular myocardial iron deposition. In none of the 10 hearts was significant macroscopic replacement fibrosis present. In only 2 hearts was interstitial fibrosis present, but with low CVF: in one patient with no cardiac siderosis (death by stroke, CVF 5.9%) and in a heart failure patient (CVF 2%). In the remaining 8 patients, no interstitial fibrosis was seen despite all having severe cardiac siderosis and heart failure (CVF 1.86% ±0.87%). CONCLUSION: Replacement cardiac fibrosis was not seen in the 9 post-mortem hearts from patients with severe cardiac siderosis and heart failure leading to death or transplantation, which contrasts markedly to historical reports. Minor interstitial fibrosis was also unusual and very limited in extent. These findings accord with the potential for reversibility of heart failure seen in iron overload cardiomyopathy. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT00520559.

Retroviruses integrate into a shared, non-palindromic DNA motif.

Many DNA-binding factors, such as transcription factors, form oligomeric complexes with structural symmetry that bind to palindromic DNA sequences1. Palindromic consensus nucleotide sequences are also found at the genomic integration sites of retroviruses2-6 and other transposable elements7-9, and it has been suggested that this palindromic consensus arises as a consequence of the structural symmetry in the integrase complex2,3. However, we show here that the palindromic consensus sequence is not present in individual integration sites of human T-cell lymphotropic virus type 1 (HTLV-1) and human immunodeficiency virus type 1 (HIV-1), but arises in the population average as a consequence of the existence of a non-palindromic nucleotide motif that occurs in approximately equal proportions on the plus strand and the minus strand of the host genome. We develop a generally applicable algorithm to sort the individual integration site sequences into plus-strand and minus-strand subpopulations, and use this to identify the integration site nucleotide motifs of five retroviruses of different genera: HTLV-1, HIV-1, murine leukaemia virus (MLV), avian sarcoma leucosis virus (ASLV) and prototype foamy virus (PFV). The results reveal a non-palindromic motif that is shared between these retroviruses.

Robustness of MEK-ERK Dynamics and Origins of Cell-to-Cell Variability in MAPK Signaling.

Cellular signaling processes can exhibit pronounced cell-to-cell variability in genetically identical cells. This affects how individual cells respond differentially to the same environmental stimulus. However, the origins of cell-to-cell variability in cellular signaling systems remain poorly understood. Here, we measure the dynamics of phosphorylated MEK and ERK across cell populations and quantify the levels of population heterogeneity over time using high-throughput image cytometry. We use a statistical modeling framework to show that extrinsic noise, particularly that from upstream MEK, is the dominant factor causing cell-to-cell variability in ERK phosphorylation, rather than stochasticity in the phosphorylation/dephosphorylation of ERK. We furthermore show that without extrinsic noise in the core module, variable (including noisy) signals would be faithfully reproduced downstream, but the within-module extrinsic variability distorts these signals and leads to a drastic reduction in the mutual information between incoming signal and ERK activity.

A graph theoretical approach to data fusion.

The rapid development of high throughput experimental techniques has resulted in a growing diversity of genomic datasets being produced and requiring analysis. Therefore, it is increasingly being recognized that we can gain deeper understanding about underlying biology by combining the insights obtained from multiple, diverse datasets. Thus we propose a novel scalable computational approach to unsupervised data fusion. Our technique exploits network representations of the data to identify similarities among the datasets. We may work within the Bayesian formalism, using Bayesian nonparametric approaches to model each dataset; or (for fast, approximate, and massive scale data fusion) can naturally switch to more heuristic modeling techniques. An advantage of the proposed approach is that each dataset can initially be modeled independently (in parallel), before applying a fast post-processing step to perform data integration. This allows us to incorporate new experimental data in an online fashion, without having to rerun all of the analysis. We first demonstrate the applicability of our tool on artificial data, and then on examples from the literature, which include yeast cell cycle, breast cancer and sporadic inclusion body myositis datasets.

Calibration of myocardial T2 and T1 against iron concentration.

BACKGROUND: The assessment of myocardial iron using T2* cardiovascular magnetic resonance (CMR) has been validated and calibrated, and is in clinical use. However, there is very limited data assessing the relaxation parameters T1 and T2 for measurement of human myocardial iron. METHODS: Twelve hearts were examined from transfusion-dependent patients: 11 with end-stage heart failure, either following death (n=7) or cardiac transplantation (n=4), and 1 heart from a patient who died from a stroke with no cardiac iron loading. Ex-vivo R1 and R2 measurements (R1=1/T1 and R2=1/T2) at 1.5 Tesla were compared with myocardial iron concentration measured using inductively coupled plasma atomic emission spectroscopy. RESULTS: From a single myocardial slice in formalin which was repeatedly examined, a modest decrease in T2 was observed with time, from mean (± SD) 23.7 ± 0.93 ms at baseline (13 days after death and formalin fixation) to 18.5 ± 1.41 ms at day 566 (p<0.001). Raw T2 values were therefore adjusted to correct for this fall over time. Myocardial R2 was correlated with iron concentration [Fe] (R2 0.566, p<0.001), but the correlation was stronger between LnR2 and Ln[Fe] (R2 0.790, p<0.001). The relation was [Fe] = 5081•(T2)-2.22 between T2 (ms) and myocardial iron (mg/g dry weight). Analysis of T1 proved challenging with a dichotomous distribution of T1, with very short T1 (mean 72.3 ± 25.8 ms) that was independent of iron concentration in all hearts stored in formalin for greater than 12 months. In the remaining hearts stored for <10 weeks prior to scanning, LnR1 and iron concentration were correlated but with marked scatter (R2 0.517, p<0.001). A linear relationship was present between T1 and T2 in the hearts stored for a short period (R2 0.657, p<0.001). CONCLUSION: Myocardial T2 correlates well with myocardial iron concentration, which raises the possibility that T2 may provide additive information to T2* for patients with myocardial siderosis. However, ex-vivo T1 measurements are less reliable due to the severe chemical effects of formalin on T1 shortening, and therefore T1 calibration may only be practical from in-vivo human studies.

Balancing the robustness and predictive performance of biomarkers.

Recent studies have highlighted the importance of assessing the robustness of putative biomarkers identified from experimental data. This has given rise to the concept of stable biomarkers, which are ones that are consistently identified regardless of small perturbations to the data. Since stability is not by itself a useful objective, we present a number of strategies that combine assessments of stability and predictive performance in order to identify biomarkers that are both robust and diagnostically useful. Moreover, by wrapping these strategies around logistic regression classifiers regularized by the elastic net penalty, we are able to assess the effects of correlations between biomarkers upon their perceived stability. We use a synthetic example to illustrate the properties of our proposed strategies. In this example, we find that: (i) assessments of stability can help to reduce the number of false-positive biomarkers, although potentially at the cost of missing some true positives; (ii) combining assessments of stability with assessments of predictive performance can improve the true positive rate; and (iii) correlations between biomarkers can have adverse effects on their stability and hence must be carefully taken into account when undertaking biomarker discovery. We then apply our strategies in a proteomics context to identify a number of robust candidate biomarkers for the human disease HTLV1-associated myelopathy/tropical spastic paraparesis (HAM/TSP).

A general moment expansion method for stochastic kinetic models.

Moment approximation methods are gaining increasing attention for their use in the approximation of the stochastic kinetics of chemical reaction systems. In this paper we derive a general moment expansion method for any type of propensities and which allows expansion up to any number of moments. For some chemical reaction systems, more than two moments are necessary to describe the dynamic properties of the system, which the linear noise approximation is unable to provide. Moreover, also for systems for which the mean does not have a strong dependence on higher order moments, moment approximation methods give information about higher order moments of the underlying probability distribution. We demonstrate the method using a dimerisation reaction, Michaelis-Menten kinetics and a model of an oscillating p53 system. We show that for the dimerisation reaction and Michaelis-Menten enzyme kinetics system higher order moments have limited influence on the estimation of the mean, while for the p53 system, the solution for the mean can require several moments to converge to the average obtained from many stochastic simulations. We also find that agreement between lower order moments does not guarantee that higher moments will agree. Compared to stochastic simulations, our approach is numerically highly efficient at capturing the behaviour of stochastic systems in terms of the average and higher moments, and we provide expressions for the computational cost for different system sizes and orders of approximation. We show how the moment expansion method can be employed to efficiently quantify parameter sensitivity. Finally we investigate the effects of using too few moments on parameter estimation, and provide guidance on how to estimate if the distribution can be accurately approximated using only a few moments.

Elucidating the in vivo phosphorylation dynamics of the ERK MAP kinase using quantitative proteomics data and Bayesian model selection.

Ever since reversible protein phosphorylation was discovered, it has been clear that it plays a key role in the regulation of cellular processes. Proteins often undergo double phosphorylation, which can occur through two possible mechanisms: distributive or processive. Which phosphorylation mechanism is chosen for a particular cellular regulation bears biological significance, and it is therefore in our interest to understand these mechanisms. In this paper we study dynamics of the MEK/ERK phosphorylation. We employ a model selection algorithm based on approximate Bayesian computation to elucidate phosphorylation dynamics from quantitative time course data obtained from PC12 cells in vivo. The algorithm infers the posterior distribution over four proposed models for phosphorylation and dephosphorylation dynamics, and this distribution indicates the amount of support given to each model. We evaluate the robustness of our inferential framework by systematically exploring different ways of parameterizing the models and for different prior specifications. The models with the highest inferred posterior probability are the two models employing distributive dephosphorylation, whereas we are unable to choose decisively between the processive and distributive phosphorylation mechanisms.

Plasma proteome analysis in HTLV-1-associated myelopathy/tropical spastic paraparesis.

BACKGROUND: Human T lymphotropic virus Type 1 (HTLV-1) causes a chronic inflammatory disease of the central nervous system known as HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM) which resembles chronic spinal forms of multiple sclerosis (MS). The pathogenesis of HAM remains uncertain. To aid in the differential diagnosis of HAM and to identify pathogenetic mechanisms, we analysed the plasma proteome in asymptomatic HTLV-1 carriers (ACs), patients with HAM, uninfected controls, and patients with MS. We used surface-enhanced laser desorption-ionization (SELDI) mass spectrometry to analyse the plasma proteome in 68 HTLV-1-infected individuals (in two non-overlapping sets, each comprising 17 patients with HAM and 17 ACs), 16 uninfected controls, and 11 patients with secondary progressive MS. Candidate biomarkers were identified by tandem Q-TOF mass spectrometry. RESULTS: The concentrations of three plasma proteins--high [ß2-microglobulin], high [Calgranulin B], and low [apolipoprotein A2]--were specifically associated with HAM, independently of proviral load. The plasma [ß2-microglobulin] was positively correlated with disease severity. CONCLUSIONS: The results indicate that monocytes are activated by contact with activated endothelium in HAM. Using ß2-microglobulin and Calgranulin B alone we derive a diagnostic algorithm that correctly classified the disease status (presence or absence of HAM) in 81% of HTLV-1-infected subjects in the cohort.

On T2* magnetic resonance and cardiac iron.

BACKGROUND: Measurement of myocardial iron is key to the clinical management of patients at risk of siderotic cardiomyopathy. The cardiovascular magnetic resonance relaxation parameter R2* (assessed clinically via its reciprocal, T2*) measured in the ventricular septum is used to assess cardiac iron, but iron calibration and distribution data in humans are limited. METHODS AND RESULTS: Twelve human hearts were studied from transfusion-dependent patients after either death (heart failure, n=7; stroke, n=1) or transplantation for end-stage heart failure (n=4). After cardiovascular magnetic resonance R2* measurement, tissue iron concentration was measured in multiple samples of each heart with inductively coupled plasma atomic emission spectroscopy. Iron distribution throughout the heart showed no systematic variation between segments, but epicardial iron concentration was higher than in the endocardium. The mean ± SD global myocardial iron causing severe heart failure in 10 patients was 5.98 ± 2.42 mg/g dry weight (range, 3.19 to 9.50 mg/g), but in 1 outlier case of heart failure was 25.9 mg/g dry weight. Myocardial ln[R2*] was strongly linearly correlated with ln[Fe] (R²=0.910, P<0.001), leading to [Fe]=45.0×(T2*)⁻¹·²² for the clinical calibration equation with [Fe] in milligrams per gram dry weight and T2* in milliseconds. Midventricular septal iron concentration and R2* were both highly representative of mean global myocardial iron. CONCLUSIONS: These data detail the iron distribution throughout the heart in iron overload and provide calibration in humans for cardiovascular magnetic resonance R2* against myocardial iron concentration. The iron values are of considerable interest in terms of the level of cardiac iron associated with iron-related death and indicate that the heart is more sensitive to iron loading than the liver. The results also validate the current clinical practice of monitoring cardiac iron in vivo by cardiovascular magnetic resonance of the midseptum.

International reproducibility of single breathhold T2* MR for cardiac and liver iron assessment among five thalassemia centers.

PURPOSE: To examine the reproducibility of the single breathhold T2* technique from different scanners, after installation of standard methodology in five international centers. MATERIALS AND METHODS: Up to 10 patients from each center were scanned twice locally for local interstudy reproducibility of heart and liver T2*, and then flown to a central MR facility to be rescanned on a reference scanner for intercenter reproducibility. Interobserver reproducibility for all scans was also assessed. RESULTS: Of the 49 patients scanned, the intercenter reproducibility for T2* was 5.9% for the heart and 5.8% for the liver. Local interstudy reproducibility for T2* was 7.4% for the heart and 4.6% for the liver. Interobserver reproducibility for T2* was 5.4% for the heart and 4.4% for the liver. CONCLUSION: These data indicate that T2* MR may be developed into a widespread test for tissue siderosis providing that well-defined and approved imaging and analysis techniques are used.