A semilinear in-slide model is introduced to remove the intensity effect in the scanning process. It is demonstrated that the intensity effect can be estimated accurately and removed effectively. This normalization step is vital for Affymetrix arrays to reveal relevant biological results when comparing gene expression in multiple arrays. The normalized expression ratios are analyzed further by a modified two-sample <i>t</i> test along with a sieved permutation scheme for computing <i>P</i> values. The improved specificity and sensitivity are demonstrated by using a study on the impact of macrophage migration inhibitory factor (MIF) reduction in neuroblastoma cells. With semilinear in-slide model analysis, expression of 166 genes was altered with a <i>P</i> value no greater than 0.001. Among those genes, 44 were altered &gt;2-fold. MIF-regulated genes associated with tumor development including IL-8 and <i>C</i>-<i>met</i>, which are overexpressed

<b>Motivation:</b> Normalization of microarray data is essential for multiple-array analyses. Several normalization protocols have been proposed based on different biological or statistical assumptions. A fundamental problem arises whether they have effectively normalized arrays. In addition, for a given array, the question arises how to choose a method to most effectively normalize the microarray data.

Estimation of genewise variance arises from two important applications in microarray data analysis: selecting significantly differentially expressed genes and validation tests for normalization of microarray data. We approach the problem by introducing a two-way nonparametric model, which is an extension of the famous Neyman-Scott model and is applicable beyond microarray data. The problem itself poses interesting challenges because the number of nuisance parameters is proportional to the sample size and it is not obvious how the variance function can be estimated when measurements are correlated. In such a high-dimensional nonparametric problem, we proposed two novel nonparametric estimators for genewise variance function and semiparametric estimators for measurement correlation, via solving a system of nonlinear equations. Their asymptotic normality is established. The finite sample property is

In an event-related functional MRI data analysis, an accurate and robust extraction of the hemodynamic response function (HRF) and its associated statistics (e.g., magnitude, width, and time to peak) is critical to infer quantitative information about the relative timing of the neuronal events in different brain regions. The aim of this paper is to develop a multiscale adaptive smoothing model (MASM) to accurately estimate HRFs pertaining to each stimulus sequence across all voxels. MASM explicitly accounts for both spatial and temporal smoothness information, while incorporating such information to adaptively estimate HRFs in the frequency domain. One simulation study and a real data set are used to demonstrate the methodology and examine its finite sample performance in HRF estimation, which confirms that MASM significantly outperforms the existing methods including the smooth finite impulse

Stem cell therapies have had tremendous potential application for many diseases in recent years. However, the tumorigeneic properties of stem cells restrict their potential clinical application; therefore, strategies for reducing the tumorigenic potential of stem cells must be established prior to transplantation. We have demonstrated that syngeneic transplantation of embryonic stem cells (ESCs) provokes an inflammatory response that involves the rapid recruitment of bone marrow-derived macrophages (BMDMs). ESCs are able to prevent mature macrophages from macrophage colony-stimulating factor (M-CSF) withdrawal-induced apoptosis, and thus prolong macrophage lifespan significantly by blocking various apoptotic pathways in an M-CSF-independent manner. ESCs express and secrete IL-34 which may be responsible for ESC-promoted macrophage survival. This anti-apoptotic effect of ESCs involves activation of extracellular signal-regulated kinase (ERK)1/2 and PI3K/Akt pathways and thus, inhibition of ERK1/2 and PI3K/AKT activation decreases ESC-induced macrophage survival. Functionally, ESC-treated macrophages also showed a higher level of phagocytic activity. ESCs further serve to polarize BMDMs into M2-like macrophages that exhibit most tumor-associated macrophage (TAM) phenotypic and functional features. ESC-educated macrophages produce high levels of arginase-1, Tie-2 and TNF-, which participate in angiogenesis and contribute to teratoma progression. Our study suggests that induction of M2-like macrophage activation is an important mechanism for teratoma development. Strategies targeting macrophages to

Dispersive instability appears in time-domain solutions of classical cochlear models. In this letter, a derivation of optimal initial data is presented to minimize the effect of instability. A second-order accurate implicit boundary integral method is introduced. Numerical solutions of two-dimensional models show that the optimal initial conditions work successfully in time-domain steady-state computations for both the zero Neumann and zero Dirichlet fluid pressure boundary conditions at the helicotrema.

Tumor growth and metastasis require that tumor cells must have either the potential to shift genetically or epigenetically between proliferative and invasive phenotypes or both phenotypes simultaneously. In the present study, we demonstrated that neuroblastoma growth and invasion were distinct processes that were carried out by proliferative and invasive phenotypes of tumor cells, respectively. Two subpopulations from human neuroblastoma cell line were isolated: highly invasive (HI) cells and low-invasive (LI) cells. HI and LI cells had different proliferative rate and metastatic ability <i>in vitro</i> and <i>in vivo</i> . In addition, they had distinct activated signal pathways and sensitivities to chemotherapy drugs. Affymetrix microarray and quantitative reverse transcriptasepolymerase chain reaction revealed that visinin-like protein-1 (VSNL-1) mRNA in HI cells was significantly higher than

Feedback modules, which appear ubiquitously in biological regulations, are often subject to disturbances from the input, leading to fluctuations in the output. Thus, the question becomes how a feedback system can produce a faithful response with a noisy input. We employed multiple time scale analysis, Fluctuation Dissipation Theorem, linear stability, and numerical simulations to investigate a module with one positive feedback loop driven by an external stimulus, and we obtained a critical quantity in noise attenuation, termed as signed activation time. We then studied the signed activation time for a system of two positive feedback loops, a system of one positive feedback loop and one negative feedback loop, and six other existing biological models consisting of multiple components along with positive and negative feedback loops. An inverse relationship is found between the noise amplification rate and the signed activation time, defined as the difference between the deactivation and activation time scales of the noise-free system, normalized by the frequency of noises presented in the input. Thus, the combination of fast activation and slow deactivation provides the best noise attenuation, and it can be attained in a single positive feedback loop system. An additional positive feedback loop often leads to a marked decrease in activation time, decrease or slight increase of deactivation time and allows larger kinetic rate variations for slow deactivation and fast activation. On the other hand, a negative feedback loop may increase the activation and deactivation times. The negative relationship between the noise amplification rate and the signed

DNA microarray analysis has emerged as a leading technology to enhance our understanding of gene regulation and function in cellular mechanism controls on a genomic scale. This technology has advanced to unravel the genetic machinery of biological rhythms by collecting massive gene-expression data in a time course. Here, we present a statistical model for clustering periodic patterns of gene expression in terms of different transcriptional profiles. The model incorporates biologically meaningful Fourier series approximations of gene periodic expression into a mixture-model-based likelihood function, thus producing results that are likely to be closer to biological relevance, as compared to those from existing models. Also because the structures of the time-dependent means and covariance matrix are modeled, the new approach displays increased statistical power and precision of parameter estimation. The

A nonlinear, nonlocal cochlear model of the transmission line type is studied in order to capture the multitone interactions and resulting tonal suppression effects. The model can serve as a module for voice signal processing, and is a one-dimensional (in space) damped dispersive nonlinear PDE based on the mechanics and phenomenology of hearing. It describes the motion of the basilar membrane (BM0 in the cochlea driven by input pressure waves. Both elastic dampling and selective longitudinal fluid damping are present. The forner is nonlinear and nonlocal in BM displacement, and plays a kep role in capturing tonal interactions. The latter is active only near the exit boundary (helicotrema), and is built in to damp out the remaining long waves. The initial boundary value problem is numerically solved with a semi-implicit second order finite difference method. Solutions reach a multi-frequency quai-steady state