The single-stranded RNA Flavivirus, Zika virus (ZIKV), has recently re-emerged and spread rapidly across the western hemispheres equatorial countries, primarily through Aedes mosquito transmission. While symptoms in adult infections appear to be self-limiting and mild, severe birth defects, such as microcephaly, have been linked to infection during early pregnancy. Recently, Tang et al. (Cell Stem Cell 2016, doi: 10.1016/j.stem.2016.02.016 ) demonstrated that ZIKV efficiently infects induced pluripotent stem cell (iPSC) derived human neural progenitor cells (hNPCs), resulting in cell cycle abnormalities and apoptosis. Consequently, hNPCs are a suggested ZIKV target. We analyzed the transcriptomic sequencing (RNA-seq) data (GEO: GSE78711) of ZIKV (Strain: MR766) infected hNPCs. For comparison to the ZIKV-infected hNPCs, the expression data from hNPCs infected with human cytomegalovirus (CMV) (Strain: AD169) was used (GEO: GSE35295). Utilizing a combination of Gene Ontology, database of human diseases, and pathway analysis, we generated a putative systemic model of infection supported by known molecular pathways of other highly related viruses. We analyzed RNA-sequencing data for transcript expression alterations in ZIKV-infected hNPCs, and then compared them to expression patterns of iPSC-derived hNPCs infected with CMV, a virus that can also induce severe congenital neurological defects in developing fetuses. We demonstrate for the first time that many of cellular pathways correlate with clinical pathologies following ZIKV infection such as microcephaly, congenital nervous

Microarray techniques have been widely used to monitor gene expression in many areas of biomedical research. They have been widely used for tumor diagnosis and classification, prediction of prognoses and treatment, and understanding of molecular mechanisms, biochemical pathways, and gene networks. Statistical methods are vital for these scientific endeavors. This article reviews recent developments of statistical methods for analyzing data from microarray experiments. Emphasis has been given to normalization of expression from multiple arrays, selecting significantly differentially expressed genes, tumor classifications, and gene expression pathways and networks.

Severe acute respiratory syndrome (SARS) is a new infectious disease with a global impact. Understanding its pathogenesis and developing specific diagnostic methods for its early diagnosis are crucial for the effective management and control of this disease. By using proteomic technology, truncated forms of ₁antitrypsin (TF₁AT) were found to increase significantly and consistently in sera of SARS patients compared to control subjects. The result showed a sensitivity of 100% for SARS patients and a specificity of 92.8% for controls. Furthermore, the levels of these proteins significantly correlated with certain clinicopathological parameters. The dramatic increase in TF₁AT may be the result of degradation of ₁AT. As ₁AT plays an important role in the protection of lung function, its degradation may be an important factor in the pathogenesis of SARS. These findings indicate that increased TF₁AT

It has been reported that the plasma levels of VEGF in tumor patients decreased during dendritic cell (DC)-based immunotherapy, but the underlying mechanism remains unclear. Our current report demonstrates that VEGF levels were significantly decreased in the supernatants of DCs incubated with rhVEGF or tumor conditioned medium (TCM) while the intracellular VEGF in DCs was increased. The increased intracellular VEGF was not due to the <i>de novo</i> VEGF synthesis by DCs because exogenous VEGF inhibited the mRNA expression of VEGF in DCs. More direct evidence was provided to demonstrate that Cy3-labeled VEGF could be internalized by DCs specifically and efficiently. In addition, the activity of DCs to internalize VEGF was abolished by neutralizing antibody against VEGF receptor-1 (Flt-1) and inhibitors of endocytosis such as carbonyl cyanide m-chlorophenyl hydrazone (CCCP) and genistein. This

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&lt; 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

Although stem cell therapy holds promise as a potential treatment in a number of diseases, the tumorigenicity of embryonic stem cells (ESC) and induced pluripotent stem cells remains a major obstacle. <i>In vitro</i> predifferentiation of ESCs can help prevent the risk of teratoma formation, yet proliferating neural progenitors can generate tumors, especially in the presence of immunosuppressive therapy. In this study, we investigated the effects of the microenvironment on stem cell growth and teratoma development using undifferentiated ESCs. Syngeneic ESC transplantation triggered an inflammatory response that involved the recruitment of bone marrow (BM)derived macrophages. These macrophages differentiated into an M2 or angiogenic phenotype that expressed multiple angiogenic growth factors and proteinases, such as macrophage migration inhibitory factor (MIF), VEGF, and matrix metalloproteinase 9

Functional linear models are useful in longitudinal data analysis. They include many classical and recently proposed statistical models for longitudinal data and other functional data. Recently, smoothing spline and kernel methods have been proposed for estimating their coefficient functions nonparametrically but these methods are either intensive in computation or inefficient in performance. To overcome these drawbacks, in this paper, a simple and powerful twostep alternative is proposed. In particular, the implementation of the proposed approach via local polynomial smoothing is discussed. Methods for estimating standard deviations of estimated coefficient functions are also proposed. Some asymptotic results for the local polynomial estimators are established. Two longitudinal data sets, one of which involves timedependent covariates, are used to demonstrate the approach proposed. Simulation studies

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

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

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