Sensible Functional Linear Discriminant Analysis Effectively Discriminates Enhanced Raman Spectra of Mycobacterium Species.

Tuberculosis caused by Mycobacterium tuberculosis complex (MTBC) is one of the major infectious diseases in the world. Identification of MTBC and differential diagnosis of nontuberculous mycobacteria (NTM) species impose challenges because of their taxonomic similarity. This study describes a differential diagnosis method using the surface-enhanced Raman scattering (SERS) measurement of molecules released by Mycobacterium species. Conventional principal component analysis and linear discriminant analysis methods successfully separated the acquired spectrum of MTBC from those of NTM species but failed to distinguish between the spectra of different NTM species. A novel sensible functional linear discriminant analysis (SLDA), projecting the averaged spectrum of a bacterial specie to the subspace orthogonal to the within-species random variation, thereby eliminating its influence in applying linear discriminant analysis, was employed to effectively discriminate not only MTBC but also species of NTM. The successful demonstration of this SERS-SLDA method opens up new opportunities for the rapid differentiation of Mycobacterium species.

Sparse linear modeling of next-generation mRNA sequencing (RNA-Seq) data for isoform discovery and abundance estimation.

Since the inception of next-generation mRNA sequencing (RNA-Seq) technology, various attempts have been made to utilize RNA-Seq data in assembling full-length mRNA isoforms de novo and estimating abundance of isoforms. However, for genes with more than a few exons, the problem tends to be challenging and often involves identifiability issues in statistical modeling. We have developed a statistical method called "sparse linear modeling of RNA-Seq data for isoform discovery and abundance estimation" (SLIDE) that takes exon boundaries and RNA-Seq data as input to discern the set of mRNA isoforms that are most likely to present in an RNA-Seq sample. SLIDE is based on a linear model with a design matrix that models the sampling probability of RNA-Seq reads from different mRNA isoforms. To tackle the model unidentifiability issue, SLIDE uses a modified Lasso procedure for parameter estimation. Compared with deterministic isoform assembly algorithms (e.g., Cufflinks), SLIDE considers the stochastic aspects of RNA-Seq reads in exons from different isoforms and thus has increased power in detecting more novel isoforms. Another advantage of SLIDE is its flexibility of incorporating other transcriptomic data such as RACE, CAGE, and EST into its model to further increase isoform discovery accuracy. SLIDE can also work downstream of other RNA-Seq assembly algorithms to integrate newly discovered genes and exons. Besides isoform discovery, SLIDE sequentially uses the same linear model to estimate the abundance of discovered isoforms. Simulation and real data studies show that SLIDE performs as well as or better than major competitors in both isoform discovery and abundance estimation. The SLIDE software package is available at https://sites.google.com/site/jingyijli/SLIDE.zip.

Smoothing dynamic positron emission tomography time courses using functional principal components.

A functional smoothing approach to the analysis of PET time course data is presented. By borrowing information across space and accounting for this pooling through the use of a nonparametric covariate adjustment, it is possible to smooth the PET time course data thus reducing the noise. A new model for functional data analysis, the Multiplicative Nonparametric Random Effects Model, is introduced to more accurately account for the variation in the data. A locally adaptive bandwidth choice helps to determine the correct amount of smoothing at each time point. This preprocessing step to smooth the data then allows subsequent analysis by methods such as Spectral Analysis to be substantially improved in terms of their mean squared error.

A Functional Approach to Deconvolve Dynamic Neuroimaging Data.

Positron emission tomography (PET) is an imaging technique which can be used to investigate chemical changes in human biological processes such as cancer development or neurochemical reactions. Most dynamic PET scans are currently analyzed based on the assumption that linear first-order kinetics can be used to adequately describe the system under observation. However, there has recently been strong evidence that this is not the case. To provide an analysis of PET data which is free from this compartmental assumption, we propose a nonparametric deconvolution and analysis model for dynamic PET data based on functional principal component analysis. This yields flexibility in the possible deconvolved functions while still performing well when a linear compartmental model setup is the true data generating mechanism. As the deconvolution needs to be performed on only a relative small number of basis functions rather than voxel by voxel in the entire three-dimensional volume, the methodology is both robust to typical brain imaging noise levels while also being computationally efficient. The new methodology is investigated through simulations in both one-dimensional functions and 2D images and also applied to a neuroimaging study whose goal is the quantification of opioid receptor concentration in the brain.

Inferring Genetic Interactions via a Data-Driven Second Order Model.

Genetic/transcriptional regulatory interactions are shown to predict partial components of signaling pathways, which have been recognized as vital to complex human diseases. Both activator (A) and repressor (R) are known to coregulate their common target gene (T). Xu et al. (2002) proposed to model this coregulation by a fixed second order response surface (called the RS algorithm), in which T is a function of A, R, and AR. Unfortunately, the RS algorithm did not result in a sufficient number of genetic interactions (GIs) when it was applied to a group of 51 yeast genes in a pilot study. Thus, we propose a data-driven second order model (DDSOM), an approximation to the non-linear transcriptional interactions, to infer genetic and transcriptional regulatory interactions. For each triplet of genes of interest (A, R, and T), we regress the expression of T at time t + 1 on the expression of A, R, and AR at time t. Next, these well-fitted regression models (viewed as points in R(3)) are collected, and the center of these points is used to identify triples of genes having the A-R-T relationship or GIs. The DDSOM and RS algorithms are first compared on inferring transcriptional compensation interactions of a group of yeast genes in DNA synthesis and DNA repair using microarray gene expression data; the DDSOM algorithm results in higher modified true positive rate (about 75%) than that of the RS algorithm, checked against quantitative RT-polymerase chain reaction results. These validated GIs are reported, among which some coincide with certain interactions in DNA repair and genome instability pathways in yeast. This suggests that the DDSOM algorithm has potential to predict pathway components. Further, both algorithms are applied to predict transcriptional regulatory interactions of 63 yeast genes. Checked against the known transcriptional regulatory interactions queried from TRANSFAC, the proposed also performs better than the RS algorithm.

Lifespan of a Ceratitis fruit fly increases with higher altitude.

Variation in lifespan may be linked to geographic factors. While latitudinal variation in lifespan has been studied for a number of species, altitude variation has received much less attention, particularly in insects. We measured the lifespan of different populations of the Natal fruit fly Ceratitis rosa along an altitudinal cline. For the different populations we first measured the residual longevity of wild flies by captive cohort approach and compared F(1) generation from the same populations. We showed an increase in lifespan with higher altitude for a part of our data. For the field collected flies (F0) the average remaining lifespan increased monotonically with altitude for males but not for females. For the F(1) generation, longevity of both males and females of the highest-altitude population was longer than for the two other lower-altitude populations. This relationship between altitude and lifespan may be explained by the effects of temperature on reproduction. Reproductive schedules in insects are linked to temperature: lower temperature, characteristic of high-altitude sites, generally slows down reproduction. Because of a strong trade-off between reproduction and longevity, we therefore observed a longer lifespan for the high- altitude populations. Other hypotheses such as different predation rates in the different sites are also discussed.