Latent group detection in functional partially linear regression models

In this paper, we propose a functional partially linear regression model with latent group structures to accommodate the heterogeneous relationship between a scalar response and functional covariates. The proposed model is motivated by a salinity tolerance study of barley families, whose main objective is to detect salinity tolerant barley plants. Our model is flexible, allowing for heterogeneous functional coefficients while being efficient by pooling information within a group for estimation. We develop an algorithm in the spirit of the K-means clustering to identify latent groups of the subjects under study. We establish the consistency of the proposed estimator, derive the convergence rate and the asymptotic distribution, and develop inference procedures. We show by simulation studies that the proposed method has higher accuracy for recovering latent groups and for estimating the functional coefficients than existing methods. The analysis of the barley data shows that the proposed method can help identify groups of barley families with different salinity tolerant abilities.     

Reconstruction of tricarboxylic acid cycle in Corynebacterium glutamicum with a genome-scale metabolic network model for trans-4-hydroxyproline production

trans-4-Hydroxy- l -proline (Hyp) is an abundant component of mammalian collagen and functions as a chiral synthon for the syntheses of anti-inflammatory drugs in the pharmaceutical industry. Proline 4-hydroxylase (P4H) can catalyze the conversion of l -proline to Hyp; however, it is still challenging for the fermentative production of Hyp from glucose using P4H due to the low yield and productivity. Here, we report the metabolic engineering of Corynebacterium glutamicum for the fermentative production of Hyp by reconstructing tricarboxylic acid (TCA) cycle together with heterologously expressing the p4h gene from Dactylosporangium sp. strain RH1. In silico model-based simulation showed that α-ketoglutarate was redirected from the TCA cycle toward Hyp synthetic pathway driven by P4H when the carbon flux from succinyl-CoA to succinate descended to zero. The interruption of the TCA cycle by the deletion of sucCD-encoding the succinyl-CoA synthetase (SUCOAS) led to a 60% increase in Hyp production and had no obvious impact on the growth rate. Fine-tuning of plasmid-borne ProB* and P4H abundances led to a significant increase in the yield of Hyp on glucose. The final engineered Hyp-7 strain produced up to 21.72 g/L Hyp with a yield of 0.27 mol/mol (Hyp/glucose) and a volumetric productivity of 0.36 g·L ⁻¹·hr ⁻¹ in the shake flask fermentation. To our knowledge, this is the highest yield and productivity achieved by microbial fermentation in a glucose-minimal medium for Hyp production. This strategy provides new insights into engineering C. glutamicum by flux coupling for the fermentative production of Hyp and related products.     

Detecting and Estimating Hormesis Using a Model-Based Approach

Hormesis is defined as a dose-response relationship that is stimulatory at low doses, but is inhibitory at higher doses. In a given experiment, it is not unusual to observe enhanced responses at low doses, however, such enhanced responses may not imply hormesis, but the random fluctuation of the data. Statistical tests can be developed to detect hormesis when enhanced responses at low concentrations are observed. We propose the use of a model-based approach to detect the presence of, and estimate the extent of, hormesis. This approach includes two steps: detection and estimation. In the detection step, we compare the full and the reduced models. The full model describes the dose-response relationship incorporating the hormetic effect; the reduced model describes the dose-response relationship without the hormetic effect. The full model is an extension of the reduced model and has an extra parameter that measures the amount of increase in response at low doses. A test of statistical significance of this extra parameter can essentially be a test for detecting hormesis. In the estimation step, we obtain the area under the best-fitted dose-response curve falling within the hormetic zone. Considering both the number of concentrations within the hormetic zone and the magnitude of the stimulatory response, we propose using the ratio of the area under the hormetic zone (AUC_H) and the area under the best-fitted curve from zero to zero equivalent point (AUC_ZEP) as an estimate of magnitude of the hormetic effect. Two numerical examples are used to illustrate the use of this model-based approach.     

基因表达聚类分析技术的现状与发展

随着多个生物基因组测序的完成、DNA芯片技术的广泛应用,基因表达数据分析已成为后基因组时代的研究热点.聚类分析能将功能相关的基因按表达谱的相似程度归纳成类,有助于对未知功能的基因进行研究,是目前基因表达分析研究的主要计算技术之一.已有多种聚类分析算法用于基因表达数据分析,各种算法因其着眼点、原理等方面的差异,而各有其优缺点.如何对各种聚类算法的有效性进行分析、并开发新型的、适合于基因表达数据分析的方法已是当务之急.     

The algorithmic anatomy of model-based evaluation

Despite many debates in the first half of the twentieth century, it is now largely a truism that humans and other animals build models of their environments and use them for prediction and control. However, model-based (MB) reasoning presents severe computational challenges. Alternative, computationally simpler, model-free (MF) schemes have been suggested in the reinforcement learning literature, and have afforded influential accounts of behavioural and neural data. Here, we study the realization of MB calculations, and the ways that this might be woven together with MF values and evaluation methods. There are as yet mostly only hints in the literature as to the resulting tapestry, so we offer more preview than review.     

Towards automatic quantification of the epicardial fat in non-contrasted CT images

In this work, we present a technique to semi-automatically quantify the epicardial fat in non-contrasted computed tomography (CT) images. The epicardial fat is very close to the pericardial fat, being separated only by the pericardium that appears in the image as a very thin line, which is hard to detect. Therefore, an algorithm that uses the anatomy of the heart was developed to detect the pericardium line via control points of the line. From the points detected an interpolation was applied based on the cubic interpolation, which was also improved to avoid incorrect interpolation that occurs when the two variables are non-monotonic. The method is validated by using a set of 40 CT images of the heart of 40 human subjects. In 62.5% of the cases only minimal user intervention was required and the results compared favourably with the results obtained by the manual process.     

Modeling the Effect of Amino Acids and Copper on Monoclonal Antibody Productivity and Glycosylation: A Modular Approach

In manufacturing monoclonal antibodies (mAbs), it is crucial to be able to predict how process conditions and supplements affect productivity and quality attributes, especially glycosylation. Supplemental inputs, such as amino acids and trace metals in the media, are reported to affect cell metabolism and glycosylation; quantifying their effects is essential for effective process development. We aim to present and validate, through a commercially relevant cell culture process, a technique for modeling such effects efficiently. While existing models can predict mAb production or glycosylation dynamics under specific process configurations, adapting them to new processes remains challenging, because it involves modifying the model structure and often requires some mechanistic understanding. Here, a modular modeling technique for adapting an existing model for a fed-batch Chinese hamster ovary (CHO) cell culture process without structural modifications or mechanistic insight is presented. Instead, data is used, obtained from designed experimental perturbations in media supplementation, to train and validate a supplemental input effect model, which is used to “patch” the existing model. The combined model can be used for model-based process development to improve productivity and to meet product quality targets more efficiently. The methodology and analysis are generally applicable to other CHO cell lines and cell types.