Bayesian consensus clustering for multivariate longitudinal data

In clinical and epidemiological studies, there is a growing interest in studying the heterogeneity among patients based on longitudinal characteristics to identify subtypes of the study population. Compared to clustering a single longitudinal marker, simultaneously clustering multiple longitudinal markers allow additional information to be incorporated into the clustering process, which reveals co-existing longitudinal patterns and generates deeper biological insight. In the current study, we propose a Bayesian consensus clustering (BCC) model for multivariate longitudinal data. Instead of arriving at a single overall clustering, the proposed model allows each marker to follow marker-specific local clustering and these local clusterings are aggregated to find a global (consensus) clustering. To estimate the posterior distribution of model parameters, a Gibbs sampling algorithm is proposed. We apply our proposed model to the primary biliary cirrhosis study to identify patient subtypes that may be associated with their prognosis. We also perform simulation studies to compare the clustering performance between the proposed model and existing models under several scenarios. The results demonstrate that the proposed BCC model serves as a useful tool for clustering multivariate longitudinal data.     

Prolyl and lysyl hydroxylases in collagen synthesis

Collagens are the most abundant proteins in the extracellular matrix. They provide a framework to build organs and tissues and give structural support to make them resistant to mechanical load and forces. Several intra‐ and extracellular modifications are needed to make functional collagen molecules, intracellular post‐translational modifications of proline and lysine residues having key roles in this. In this article, we provide a review on the enzymes responsible for the proline and lysine modifications, that is collagen prolyl 4‐hydroxylases, 3‐hydroxylases and lysyl hydroxylases, and discuss their biological functions and involvement in diseases.     

我国家庭医生签约服务政策执行的制约因素与优化路径：基于史密斯政策执行过程模型

深化家庭医生签约服务是深化医药卫生体制改革、强化基层医疗卫生服务、实现"健康中国"战略目标的重要选择，也是当前更好维护人民群众健康的重要途径。为有效推进签约服务工作，国家陆续推出各项政策，全国各地也在积极进行实践探索，成效明显。但是，签约服务仍面临诸多问题，其中"执行难"是签约服务深度推进的一大困境。通过史密斯政策执行过程模型，结合签约服务政策执行过程，发现签约服务仍存在法治性不足、政策执行人员水平不高、激励不足、政策环境影响等诸多制约因素。因此，需要从法律和制度方面进行顶层设计、提升执行人员素质和职业认同、建立医患互信、优化政策执行环境等角度进行政策创新，探索家庭医生签约服务可持续发展的路径。     

Deepm5C: A deep-learning-based hybrid framework for identifying human RNA N5-methylcytosine sites using a stacking strategy

1. Download : Download high-res image (173KB)
2. Download : Download full-size image

     

A Coupled FEM-BEM Approach for the Solution of the Free-Boundary Axi-Symmetric Plasma Equilibrium Problem

In this paper we present a coupled Finite Element Method – Boundary Element Method (FEM-BEM) approach for the solution of the free-boundary axi-symmetric plasma equilibrium problem. The proposed method, obtained from an improvement of the Hagenow-Lackner coupling method, allows to efficiently model the equilibrium problem in unbounded domains by discretizing only the plasma region; the external conductors can be modelled either as 2D or 3D models, according to the problem of interest. The paper explores different iterative methods for the solution of the nonlinear Grad-Shafranov equation, such as Picard, Newton-Raphson and Newton-Krylov, in order to provide a robust and reliable tool, able to handle large-scale problems (e.g. high resolution equilibria). This method has been implemented in the FRIDA code (FRee-boundary Integro-Differential Axisimmetric – https://github. om/matteobonotto/ FRIDA), together with a suitable Adaptive Integration Technique (AIT) for the computation of the source term. FRIDA has been successfully tested and validated against experimental data from RFX-mod device, and numerical equilibria of an ITER-like device.     

Predicting Survival for Chimeric Antigen Receptor T-Cell Therapy: A Validation of Survival Models Using Follow-Up Data From ZUMA-1

ObjectivesSurvival extrapolation for chimeric antigen receptor T-cell therapies is challenging, owing to their unique mechanistic properties that translate to complex hazard functions. Axicabtagene ciloleucel is indicated for the treatment of relapse or refractory diffuse large B-cell lymphoma after 2 or more lines of therapy based on the ZUMA-1 trial. Four data snapshots are available, with minimum follow-up of 12, 24, 36, and 48 months. This analysis explores how survival extrapolations for axicabtagene ciloleucel using ZUMA-1 data can be validated and compared.MethodsThree different parametric modeling approaches were applied: standard parametric, spline-based, and cure-based models. Models were compared using a range of metrics, across the 4 data snapshot, including visual fit, plausibility of long-term estimates, statistical goodness of fit, inspection of hazard plots, point-estimate accuracy, and conditional survival estimates.ResultsStandard and spline-based parametric extrapolations were generally incapable of fitting the ZUMA-1 data well. Cure-based models provided the best fit based on the earliest data snapshot, with extrapolations remaining consistent as data matured. At 48 months, the maximum survival overestimate was 8.3% (Gompertz mixture-cure model) versus the maximum underestimate of 33.5% (Weibull standard parametric model).ConclusionsWhere a plateau in the survival curve is clinically plausible, cure-based models may be helpful in making accurate predictions based on immature data. The ability to reliably extrapolate from maturing data may reduce delays in patient access to potentially lifesaving treatments. Additional research is required to understand how models compare in broader contexts, including different treatments and therapeutic areas.     

How structural and functional MRI can inform dual-site tACS parameters: A case study in a clinical population and its pragmatic implications

BackgroundAbnormalities in frontoparietal network (FPN) were observed in many neuropsychiatric diseases including substance use disorders. A growing number of studies are using dual-site-tACS with frontoparietal synchronization to engage this network. However, a computational pathway to inform and optimize parameter space for frontoparietal synchronization is still lacking. In this case study, in a group of participants with methamphetamine use disorders, we proposed a computational pathway to extract optimal electrode montage while accounting for stimulation intensity using structural and functional MRI.MethodsSixty methamphetamine users completed an fMRI drug cue-reactivity task. Four main steps were taken to define electrode montage and adjust stimulation intensity using 4x1 high-definition (HD) electrodes for a dual-site-tACS; (1) Frontal seed was defined based on the maximum electric fields (EF) predicted by simulation of HD montage over DLPFC (F3/F4 in EEG 10–10), (2) frontal seed-to-whole brain context-dependent correlation was calculated to determine connected regions to frontal seeds, (3) center of connected cluster in parietal cortex was selected as a location for placing the second set of HD electrodes to shape the informed montage, (4) individualized head models were used to determine optimal stimulation intensity considering underlying brain structure. The informed montage was compared to montages with large electrodes and classic frontoparietal HD montages (F3-P3/F4-P4) in terms of tACS-induced EF and ROI-to-ROI task-based/resting-state connectivity.ResultsCompared to the large electrodes, HD frontoparietal montages allow for a finer control of the spatial peak fields in the main nodes of the FPN at the cost of lower maximum EF (large-pad/HD: max EF[V/m] = 0.37/0.11, number of cortical sub-regions that EF exceeds 50% of the max = 77/13). For defining stimulation targets based on EF patterns, using group-level head models compared to a single standard head model results in comparable but significantly different seed locations (6.43 mm Euclidean distance between the locations of the frontal maximum EF in standard-space). As expected, significant task-based/resting-state connections were only found between frontal-parietal locations in the informed montage. Cue-induced craving score was correlated with frontoparietal connectivity only in the informed montage (r = ?0.24). Stimulation intensity in the informed montage, and not in the classic HD montage, needs 40% reduction in the parietal site to reduce the disparity in EF between stimulation sites.ConclusionThis study provides some empirical insights to montage and dose selection in dual-site-tACS using individual brain structures and functions and proposes a computational pathway to use head models and functional MRI to define (1) optimum electrode montage for targeting FPN in a context of interest (drug-cue-reactivity) and (2) proper transcranial stimulation intensity.