High throughput metabolomics makes it possible to measure the relative abundances of numerous metabolites in biological samples, which is useful to many areas of biomedical research. However, missing values (MVs) in metabolomics datasets are common and can arise due to both technical and biological reasons. Typically, such MVs are substituted by a minimum value, which may lead to different results in downstream analyses. Here we present a modified version of the K-nearest neighbor (KNN) approach which accounts for truncation at the minimum value, i.e., KNN truncation (KNN-TN). We compare imputation results based on KNN-TN with results from other KNN approaches such as KNN based on correlation (KNN-CR) and KNN based on Euclidean distance (KNN-EU). Our approach assumes that the data follow a truncated normal distribution with the truncation point at the detection limit (LOD). The effectiveness of each approach was analyzed by the root mean square error (RMSE) measure as well as the metabolite list concordance index (MLCI) for influence on downstream statistical testing. Through extensive simulation studies and application to three real data sets, we show that KNN-TN has lower RMSE values compared to the other two KNN procedures as well as simpler imputation methods based on substituting missing values with the metabolite mean, zero values, or the LOD. MLCI values between KNN-TN and KNN-EU were roughly equivalent, and superior to the other four methods in most cases. Our findings demonstrate that KNN-TN generally has improved performance in imputing the missing values of the different datasets compared to KNN-CR and KNN-EU when there is missingness due to missing at random combined with an LOD. The results shown in this study are in the field of metabolomics but this method could be applicable with any high throughput technology which has missing due to LOD.

Spatial analysis has been vital in mapping the spread of diseases and assisting in policy making. Targeting diarrhea transmission hotspots is one of the potential strategies for reducing diarrhea cases. This study aimed to examine the spatial-temporal variations and to identify the modifiable determinants of diarrhea while controlling for the spatial dependence in the data. An ecological study on diarrhea data from DLHS-3 and NFHS- 4 in India. Moran’s I and LISA were used to detect the spatial clustering of diarrhea cases and to test for clustering in the data. Spatial regression was used to identify the modifiable factors associated with the prevalence of diarrhea. The study comprised of the prevalence of diarrhea among the children below the age of five years (U-5 s) across different states in India. The determinants of diarrhea were obtained using spatial lag models. The software used were GeoDa 1.6.6 and QGIS 2.0. The presence of spatial autocorrelation in DLHS-3 and NFHS-4 (Moron’s I = 0.577 and 0.369 respectively) enforces the usage of geographical properties while modeling the diarrhea data. The geographic clustering of high-prevalence districts was observed in the state of UP consistently. The spatial pattern of the percentage of children with diarrhea was persistently associated with the household with a sanitation facility (%) (p = 0.023 and p = 0.011). Compared to the diarrhea cases in the period 2007–2008, no much reduction was observed in the period 2015–2016. The prevalence of diarrhea and percentage of household with sanitation were ranging between 0.1–33.8% and 1.3–96.1% in the period 2007–2008 and 0.6–29.1% and 10.4–92.0% in the period 2015–2016 respectively. The least and highest prevalence of diarrhea being consistently from Assam and UP respectively. Despite improvements in controlling spread of diarrheal disease, the burden remains high. Focus on widespread diarrheal disease control strategy by addressing the social determinants of health like basic sanitation is crucial to reduce the burden of diarrhea among U-5 s in India. The identification of hotspots will aid in the planning of control strategies for goal setting in the targeted regions.

In this paper, we report the development of the nanotube-CTC-chip for isolation of tumor-derived epithelial cells (circulating tumor cells, CTCs) from peripheral blood, with high purity, by exploiting the physical mechanisms of preferential adherence of CTCs on a nanotube surface. The nanotube-CTC-chip is a new 76-element microarray technology that combines carbon nanotube surfaces with microarray batch manufacturing techniques for the capture and isolation of tumor-derived epithelial cells. Using a combination of red blood cell (RBC) lysis and preferential adherence, we demonstrate the capture and enrichment of CTCs with a 5-log reduction of contaminating WBCs. EpCAM negative MDA-MB-231/luciferase-2A-green fluorescent protein (GFP) cells were spiked in the blood of wild mice and enriched using an RBC lysis protocol. The enriched samples were then processed using the nanotube-CTC-chip for preferential CTC adherence on the nanosurface and counting the GFP cells yielded anywhere from 89% to 100% capture from the droplets. Electron microscopy (EM) studies showed focal adhesion with filaments from the cell body to the nanotube surface. We compared the nanotube preferential adherence to collagen adhesion matrix (CAM) scaffolding, reported as a viable strategy for CTC capture in patients. The CAM scaffolding on the device surface yielded 50% adherence with 100% tracking of cancer cells (adhered vs. non-adhered) versus carbon nanotubes with >90% adherence and 100% tracking for the same protocol. The nanotube-CTC-chip successfully captured CTCs in the peripheral blood of breast cancer patients (stage 1–4) with a range of 4–238 CTCs per 8.5 ml blood or 0.5–28 CTCs per ml. CTCs (based on CK8/18, Her2, EGFR) were successfully identified in 7/7 breast cancer patients, and no CTCs were captured in healthy controls (n = 2). CTC enumeration based on multiple markers using the nanotube-CTC-chip enables dynamic views of metastatic progression and could potentially have predictive capabilities for diagnosis and treatment response.

Normalization is an essential step with considerable impact on high-throughput RNA sequencing (RNA-seq) data analysis. Although there are numerous methods for read count normalization, it remains a challenge to choose an optimal method due to multiple factors contributing to read count variability that affects the overall sensitivity and specificity. In order to properly determine the most appropriate normalization methods, it is critical to compare the performance and shortcomings of a representative set of normalization routines based on different dataset characteristics. Therefore, we set out to evaluate the performance of the commonly used methods (DESeq, TMM-edgeR, FPKM-CuffDiff, TC, Med UQ and FQ) and two new methods we propose: Med-pgQ2 and UQ-pgQ2 (per-gene normalization after per-sample median or upper-quartile global scaling). Our per-gene normalization approach allows for comparisons between conditions based on similar count levels. Using the benchmark Microarray Quality Control Project (MAQC) and simulated datasets, we performed differential gene expression analysis to evaluate these methods. When evaluating MAQC2 with two replicates, we observed that Med-pgQ2 and UQ-pgQ2 achieved a slightly higher area under the Receiver Operating Characteristic Curve (AUC), a specificity rate > 85%, the detection power > 92% and an actual false discovery rate (FDR) under 0.06 given the nominal FDR (≤0.05). Although the top commonly used methods (DESeq and TMM-edgeR) yield a higher power (>93%) for MAQC2 data, they trade off with a reduced specificity (<70%) and a slightly higher actual FDR than our proposed methods. In addition, the results from an analysis based on the qualitative characteristics of sample distribution for MAQC2 and human breast cancer datasets show that only our gene-wise normalization methods corrected data skewed towards lower read counts. However, when we evaluated MAQC3 with less variation in five replicates, all methods performed similarly. Thus, our proposed Med-pgQ2 and UQ-pgQ2 methods perform slightly better for differential gene analysis of RNA-seq data skewed towards lowly expressed read counts with high variation by improving specificity while maintaining a good detection power with a control of the nominal FDR level.

Selection of informative genes is an important problem in gene expression studies. The small sample size and the large number of genes in gene expression data make the selection process complex. Further, the selected informative genes may act as a vital input for gene co-expression network analysis. Moreover, the identification of hub genes and module interactions in gene co-expression networks is yet to be fully explored. This paper presents a statistically sound gene selection technique based on support vector machine algorithm for selecting informative genes from high dimensional gene expression data. Also, an attempt has been made to develop a statistical approach for identification of hub genes in the gene co-expression network. Besides, a differential hub gene analysis approach has also been developed to group the identified hub genes into various groups based on their gene connectivity in a case vs. control study. Based on this proposed approach, an R package, i.e., dhga (https://cran.r-project.org/web/packages/dhga) has been developed. The comparative performance of the proposed gene selection technique as well as hub gene identification approach was evaluated on three different crop microarray datasets. The proposed gene selection technique outperformed most of the existing techniques for selecting robust set of informative genes. Based on the proposed hub gene identification approach, a few number of hub genes were identified as compared to the existing approach, which is in accordance with the principle of scale free property of real networks. In this study, some key genes along with their Arabidopsis orthologs has been reported, which can be used for Aluminum toxic stress response engineering in soybean. The functional analysis of various selected key genes revealed the underlying molecular mechanisms of Aluminum toxic stress response in soybean.

Over the last decade, gene set analysis has become the first choice for gaining insights into underlying complex biology of diseases through gene expression and gene association studies. It also reduces the complexity of statistical analysis and enhances the explanatory power of the obtained results. Although gene set analysis approaches are extensively used in gene expression and genome wide association data analysis, the statistical structure and steps common to these approaches have not yet been comprehensively discussed, which limits their utility. In this article, we provide a comprehensive overview, statistical structure and steps of gene set analysis approaches used for microarrays, RNA-sequencing and genome wide association data analysis. Further, we also classify the gene set analysis approaches and tools by the type of genomic study, null hypothesis, sampling model and nature of the test statistic, etc. Rather than reviewing the gene set analysis approaches individually, we provide the generation-wise evolution of such approaches for microarrays, RNA-sequencing and genome wide association studies and discuss their relative merits and limitations. Here, we identify the key biological and statistical challenges in current gene set analysis, which will be addressed by statisticians and biologists collectively in order to develop the next generation of gene set analysis approaches. Further, this study will serve as a catalog and provide guidelines to genome researchers and experimental biologists for choosing the proper gene set analysis approach based on several factors.