Macrophage activation and persistent inflammation contribute to the pathological process of spinal cord injury (SCI). It was reported that M2 macrophages were induced at 37 days after SCI but M2 markers were reduced or eliminated after 1 week. By contrast, M1 macrophage response is rapidly induced and then maintained at injured spinal cord. However, factors that modulate macrophage phenotype and function are poorly understood. We developed a model to distinguish bonemarrow derived macrophages (BMDMs) from residential microglia and explored how BMDMs change their phenotype and functions in response to the lesionrelated factors in injured spinal cord. Infiltrating BMDMs expressing higher Mac2 and lower CX3CR1 migrate to the epicenter of injury, while microglia expressing lower Mac2 but higher CX3CR1 distribute to the edges of lesion. Myelin debris at the lesion site switches BMDMs

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.

We study a class of nonlinear nonlocal cochlear models of the transmission line type, describing the motion of basilar membrane (BM) in the cochlea. They are damped dispersive partial differential equations (PDEs) driven by time dependent boundary forcing due to the input sounds. The global well-posedness in time follows from energy estimates. Uniform bounds of solutions hold in the case of bounded nonlinear damping. When the input sounds are multi-frequency tones, and the nonlinearity in the PDEs is cubic, we construct smooth quasi-periodic solutions (multi-tone solutions) in the weakly nonlinear regime, where new frequencies are generated due to nonlinear interaction. When the input consists of two tones at frequencies f 1, f 2 (f 1< f 2), and high enough intensities, numerical results illustrate the formation of combination tones at 2f 1 f 2 and 2f 2 f 1, in agreement with hearing experiments. We visualize

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

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