A Gene Expression Classifier from Whole Blood Distinguishes Benign from Malignant Lung Nodules Detected by Low-Dose CT.

Low-dose CT (LDCT) is widely accepted as the preferred method for detecting pulmonary nodules. However, the determination of whether a nodule is benign or malignant involves either repeated scans or invasive procedures that sample the lung tissue. Noninvasive methods to assess these nodules are needed to reduce unnecessary invasive tests. In this study, we have developed a pulmonary nodule classifier (PNC) using RNA from whole blood collected in RNA-stabilizing PAXgene tubes that addresses this need. Samples were prospectively collected from high-risk and incidental subjects with a positive lung CT scan. A total of 821 samples from 5 clinical sites were analyzed. Malignant samples were predominantly stage 1 by pathologic diagnosis and 97% of the benign samples were confirmed by 4 years of follow-up. A panel of diagnostic biomarkers was selected from a subset of the samples assayed on Illumina microarrays that achieved a ROC-AUC of 0.847 on independent validation. The microarray data were then used to design a biomarker panel of 559 gene probes to be validated on the clinically tested NanoString nCounter platform. RNA from 583 patients was used to assess and refine the NanoString PNC (nPNC), which was then validated on 158 independent samples (ROC-AUC = 0.825). The nPNC outperformed three clinical algorithms in discriminating malignant from benign pulmonary nodules ranging from 6-20 mm using just 41 diagnostic biomarkers. Overall, this platform provides an accurate, noninvasive method for the diagnosis of pulmonary nodules in patients with non-small cell lung cancer. SIGNIFICANCE: These findings describe a minimally invasive and clinically practical pulmonary nodule classifier that has good diagnostic ability at distinguishing benign from malignant pulmonary nodules.

Comparative functional genomics of adaptation to muscular disuse in hibernating mammals.

Hibernation is an energy-saving adaptation that involves a profound suppression of physical activity that can continue for 6-8 months in highly seasonal environments. While immobility and disuse generate muscle loss in most mammalian species, in contrast, hibernating bears and ground squirrels demonstrate limited muscle atrophy over the prolonged periods of physical inactivity during winter, suggesting that hibernating mammals have adaptive mechanisms to prevent disuse muscle atrophy. To identify common transcriptional programmes that underlie molecular mechanisms preventing muscle loss, we conducted a large-scale gene expression screen in hind limb muscles comparing hibernating and summer-active black bears and arctic ground squirrels using custom 9600 probe cDNA microarrays. A molecular pathway analysis showed an elevated proportion of overexpressed genes involved in all stages of protein biosynthesis and ribosome biogenesis in muscle of both species during torpor of hibernation that suggests induction of translation at different hibernation states. The induction of protein biosynthesis probably contributes to attenuation of disuse muscle atrophy through the prolonged periods of immobility of hibernation. The lack of directional changes in genes of protein catabolic pathways does not support the importance of metabolic suppression for preserving muscle mass during winter. Coordinated reduction in multiple genes involved in oxidation-reduction and glucose metabolism detected in both species is consistent with metabolic suppression and lower energy demand in skeletal muscle during inactivity of hibernation.

The peripheral immune response and lung cancer prognosis.

Attempts to refine and improve outcome predictions using tumor gene expression have been recently reported. We show that peripheral blood mononuclear cell (PBMC)-associated gene signatures can predict outcome in non-small cell lung carcinoma patients independent of demographic data or TNM staging, and that this information may persist after tumor resection.

Peripheral immune cell gene expression predicts survival of patients with non-small cell lung cancer.

Prediction of cancer recurrence in patients with non-small cell lung cancer (NSCLC) currently relies on the assessment of clinical characteristics including age, tumor stage, and smoking history. A better prediction of early stage cancer patients with poorer survival and late stage patients with better survival is needed to design patient-tailored treatment protocols. We analyzed gene expression in RNA from peripheral blood mononuclear cells (PBMC) of NSCLC patients to identify signatures predictive of overall patient survival. We find that PBMC gene expression patterns from NSCLC patients, like patterns from tumors, have information predictive of patient outcomes. We identify and validate a 26 gene prognostic panel that is independent of clinical stage. Many additional prognostic genes are specific to myeloid cells and are more highly expressed in patients with shorter survival. We also observe that significant numbers of prognostic genes change expression levels in PBMC collected after tumor resection. These post-surgery gene expression profiles may provide a means to re-evaluate prognosis over time. These studies further suggest that patient outcomes are not solely determined by tumor gene expression profiles but can also be influenced by the immune response as reflected in peripheral immune cells.

Delta-aminolevulinate-induced host-parasite porphyric disparity for selective photolysis of transgenic Leishmania in the phagolysosomes of mononuclear phagocytes: a potential novel platform for vaccine delivery.

Leishmania double transfectants (DTs) expressing the 2nd and 3rd enzymes in the heme biosynthetic pathway were previously reported to show neogenesis of uroporphyrin I (URO) when induced with delta-aminolevulinate (ALA), the product of the 1st enzyme in the pathway. The ensuing accumulation of URO in DT promastigotes rendered them light excitable to produce reactive oxygen species (ROS), resulting in their cytolysis. Evidence is presented showing that the DTs retained wild-type infectivity to their host cells and that the intraphagolysosomal/parasitophorous vacuolar (PV) DTs remained ALA inducible for uroporphyrinogenesis/photolysis. Exposure of DT-infected cells to ALA was noted by fluorescence microscopy to result in host-parasite differential porphyrinogenesis: porphyrin fluorescence emerged first in the host cells and then in the intra-PV amastigotes. DT-infected and control cells differed qualitatively and quantitatively in their porphyrin species, consistent with the expected multi- and monoporphyrinogenic specificities of the host cells and the DTs, respectively. After ALA removal, the neogenic porphyrins were rapidly lost from the host cells but persisted as URO in the intra-PV DTs. These DTs were thus extremely light sensitive and were lysed selectively by illumination under nonstringent conditions in the relatively ROS-resistant phagolysosomes. Photolysis of the intra-PV DTs returned the distribution of major histocompatibility complex (MHC) class II molecules and the global gene expression profiles of host cells to their preinfection patterns and, when transfected with ovalbumin, released this antigen for copresentation with MHC class I molecules. These Leishmania mutants thus have considerable potential as a novel model of a universal vaccine carrier for photodynamic immunotherapy/immunoprophylaxis.

Preservation of bone mass and structure in hibernating black bears (Ursus americanus) through elevated expression of anabolic genes.

Physical inactivity reduces mechanical load on the skeleton, which leads to losses of bone mass and strength in non-hibernating mammalian species. Although bears are largely inactive during hibernation, they show no loss in bone mass and strength. To obtain insight into molecular mechanisms preventing disuse bone loss, we conducted a large-scale screen of transcriptional changes in trabecular bone comparing winter hibernating and summer non-hibernating black bears using a custom 12,800 probe cDNA microarray. A total of 241 genes were differentially expressed (P < 0.01 and fold change >1.4) in the ilium bone of bears between winter and summer. The Gene Ontology and Gene Set Enrichment Analysis showed an elevated proportion in hibernating bears of overexpressed genes in six functional sets of genes involved in anabolic processes of tissue morphogenesis and development including skeletal development, cartilage development, and bone biosynthesis. Apoptosis genes demonstrated a tendency for downregulation during hibernation. No coordinated directional changes were detected for genes involved in bone resorption, although some genes responsible for osteoclast formation and differentiation (Ostf1, Rab9a, and c-Fos) were significantly underexpressed in bone of hibernating bears. Elevated expression of multiple anabolic genes without induction of bone resorption genes, and the down regulation of apoptosis-related genes, likely contribute to the adaptive mechanism that preserves bone mass and structure through prolonged periods of immobility during hibernation.

Genome-wide maps of circulating miRNA biomarkers for ulcerative colitis.

Inflammatory Bowel Disease--comprised of Crohn's Disease and Ulcerative Colitis (UC)--is a complex, multi-factorial inflammatory disorder of the gastrointestinal tract. In this study we have explored the utility of naturally occurring circulating miRNAs as potential blood-based biomarkers for non-invasive prediction of UC incidences. Whole genome maps of circulating miRNAs in micro-vesicles, Peripheral Blood Mononuclear Cells and platelets have been constructed from a cohort of 20 UC patients and 20 normal individuals. Through Significance Analysis of Microarrays, a signature of 31 differentially expressed platelet-derived miRNAs has been identified and biomarker performance estimated through a non-probabilistic binary linear classification using Support Vector Machines. Through this approach, classifier measurements reveal a predictive score of 92.8% accuracy, 96.2% specificity and 89.5% sensitivity in distinguishing UC patients from normal individuals. Additionally, the platelet-derived biomarker signature can be validated at 88% accuracy through qPCR assays, and a majority of the miRNAs in this panel can be demonstrated to sub-stratify into 4 highly correlated intensity based clusters. Analysis of predicted targets of these biomarkers reveal an enrichment of pathways associated with cytoskeleton assembly, transport, membrane permeability and regulation of transcription factors engaged in a variety of regulatory cascades that are consistent with a cell-mediated immune response model of intestinal inflammation. Interestingly, comparison of the miRNA biomarker panel and genetic loci implicated in IBD through genome-wide association studies identifies a physical linkage between hsa-miR-941 and a UC susceptibility loci located on Chr 20. Taken together, analysis of these expression maps outlines a promising catalog of novel platelet-derived miRNA biomarkers of clinical utility and provides insight into the potential biological function of these candidates in disease pathogenesis.

Resection of non-small cell lung cancers reverses tumor-induced gene expression changes in the peripheral immune system.

PURPOSE: To characterize the interactions of non-small cell lung cancer (NSCLC) tumors with the immune system at the level of mRNA and microRNA (miRNA) expression and to define expression signatures that characterize the presence of a malignant tumor versus a nonmalignant nodule. EXPERIMENTAL DESIGN: We have examined the changes of both mRNA and miRNA expression levels in peripheral blood mononuclear cells (PBMC) between paired samples collected from NSCLC patients before and after tumor removal using Illumina gene expression arrays. RESULTS: We found that malignant tumor removal significantly changes expression of more than 3,000 protein-coding genes, especially genes in pathways associated with suppression of the innate immune response, including natural killer cell signaling and apoptosis-associated ceramide signaling. Binding sites for the ETS domain transcription factors ELK1, ELK4, and SPI1 were enriched in promoter regions of genes upregulated in the presence of a tumor. Additional important regulators included five miRNAs expressed at significantly higher levels before tumor removal. Repressed protein-coding targets of those miRNAs included many transcription factors, several involved in immunologically important pathways. Although there was a significant overlap in the effects of malignant tumors and benign lung nodules on PBMC gene expression, we identified one gene panel which indicates a tumor or nodule presence and a second panel that can distinguish malignant from nonmalignant nodules. CONCLUSIONS: A tumor presence in the lung influences mRNA and miRNA expression in PBMC and this influence is reversed by tumor removal. These results suggest that PBMC gene expression signatures could be used for lung cancer diagnosis.

Modulation of gene expression in heart and liver of hibernating black bears (Ursus americanus).

BACKGROUND: Hibernation is an adaptive strategy to survive in highly seasonal or unpredictable environments. The molecular and genetic basis of hibernation physiology in mammals has only recently been studied using large scale genomic approaches. We analyzed gene expression in the American black bear, Ursus americanus, using a custom 12,800 cDNA probe microarray to detect differences in expression that occur in heart and liver during winter hibernation in comparison to summer active animals. RESULTS: We identified 245 genes in heart and 319 genes in liver that were differentially expressed between winter and summer. The expression of 24 genes was significantly elevated during hibernation in both heart and liver. These genes are mostly involved in lipid catabolism and protein biosynthesis and include RNA binding protein motif 3 (Rbm3), which enhances protein synthesis at mildly hypothermic temperatures. Elevated expression of protein biosynthesis genes suggests induction of translation that may be related to adaptive mechanisms reducing cardiac and muscle atrophies over extended periods of low metabolism and immobility during hibernation in bears. Coordinated reduction of transcription of genes involved in amino acid catabolism suggests redirection of amino acids from catabolic pathways to protein biosynthesis. We identify common for black bears and small mammalian hibernators transcriptional changes in the liver that include induction of genes responsible for fatty acid ß oxidation and carbohydrate synthesis and depression of genes involved in lipid biosynthesis, carbohydrate catabolism, cellular respiration and detoxification pathways. CONCLUSIONS: Our findings show that modulation of gene expression during winter hibernation represents molecular mechanism of adaptation to extreme environments.

Sp100 as a potent tumor suppressor: accelerated senescence and rapid malignant transformation of human fibroblasts through modulation of an embryonic stem cell program.

Identifying the functions of proteins, which associate with specific subnuclear structures, is critical to understanding eukaryotic nuclear dynamics. Sp100 is a prototypical protein of ND10/PML nuclear bodies, which colocalizes with Daxx and the proto-oncogenic PML. Sp100 isoforms contain SAND, PHD, Bromo, and HMG domains and are highly sumoylated, all characteristics suggestive of a role in chromatin-mediated gene regulation. A role for Sp100 in oncogenesis has not been defined previously. Using selective Sp100 isoform-knockdown approaches, we show that normal human diploid fibroblasts with reduced Sp100 levels rapidly senesce. Subsequently, small rapidly dividing Sp100 minus cells emerge from the senescing fibroblasts and are found to be highly tumorigenic in nude mice. The derivation of these tumorigenic cells from the parental fibroblasts is confirmed by microsatellite analysis. The small rapidly dividing Sp100 minus cells now also lack ND10/PML bodies, and exhibit genomic instability and p53 cytoplasmic sequestration. They have also activated MYC, RAS, and TERT pathways and express mesenchymal to epithelial transdifferentiation (MET) markers. Reintroduction of expression of only the Sp100A isoform is sufficient to maintain senescence and to inhibit emergence of the highly tumorigenic cells. Global transcriptome studies, quantitative PCR, and protein studies, as well as immunolocalization studies during the course of the transformation, reveal that a transient expression of stem cell markers precedes the malignant transformation. These results identify a role for Sp100 as a tumor suppressor in addition to its role in maintaining ND10/PML bodies and in the epigenetic regulation of gene expression.

Gene expression profiles in peripheral blood mononuclear cells can distinguish patients with non-small cell lung cancer from patients with nonmalignant lung disease.

Early diagnosis of lung cancer followed by surgery presently is the most effective treatment for non-small cell lung cancer (NSCLC). An accurate, minimally invasive test that could detect early disease would permit timely intervention and potentially reduce mortality. Recent studies have shown that the peripheral blood can carry information related to the presence of disease, including prognostic information and information on therapeutic response. We have analyzed gene expression in peripheral blood mononuclear cell samples including 137 patients with NSCLC tumors and 91 patient controls with nonmalignant lung conditions, including histologically diagnosed benign nodules. Subjects were primarily smokers and former smokers. We have identified a 29-gene signature that separates these two patient classes with 86% accuracy (91% sensitivity, 80% specificity). Accuracy in an independent validation set, including samples from a new location, was 78% (sensitivity of 76% and specificity of 82%). An analysis of this NSCLC gene signature in 18 NSCLCs taken presurgery, with matched samples from 2 to 5 months postsurgery, showed that in 78% of cases, the signature was reduced postsurgery and disappeared entirely in 33%. Our results show the feasibility of using peripheral blood gene expression signatures to identify early-stage NSCLC in at-risk populations.

Classification and biomarker identification using gene network modules and support vector machines.

BACKGROUND: Classification using microarray datasets is usually based on a small number of samples for which tens of thousands of gene expression measurements have been obtained. The selection of the genes most significant to the classification problem is a challenging issue in high dimension data analysis and interpretation. A previous study with SVM-RCE (Recursive Cluster Elimination), suggested that classification based on groups of correlated genes sometimes exhibits better performance than classification using single genes. Large databases of gene interaction networks provide an important resource for the analysis of genetic phenomena and for classification studies using interacting genes.We now demonstrate that an algorithm which integrates network information with recursive feature elimination based on SVM exhibits good performance and improves the biological interpretability of the results. We refer to the method as SVM with Recursive Network Elimination (SVM-RNE) RESULTS: Initially, one thousand genes selected by t-test from a training set are filtered so that only genes that map to a gene network database remain. The Gene Expression Network Analysis Tool (GXNA) is applied to the remaining genes to form n clusters of genes that are highly connected in the network. Linear SVM is used to classify the samples using these clusters, and a weight is assigned to each cluster based on its importance to the classification. The least informative clusters are removed while retaining the remainder for the next classification step. This process is repeated until an optimal classification is obtained. CONCLUSION: More than 90% accuracy can be obtained in classification of selected microarray datasets by integrating the interaction network information with the gene expression information from the microarrays.The Matlab version of SVM-RNE can be downloaded from http://web.macam.ac.il/~myousef.

A study of microRNAs in silico and in vivo: bioinformatics approaches to microRNA discovery and target identification.

The discovery that microRNAs (miRNAs) are synthesized as hairpin-containing precursors and share many features has stimulated the development of several computational approaches for identifying new miRNA genes in various animal species. Many of these approaches rely heavily on conservation of sequence within and between species, whereas others emphasize machine-learning methods to screen hairpin candidates for structural features shared with known miRNA precursors. The identification of animal miRNA targets is a particularly difficult problem because an exact match to the target sequence is not required. We discuss the most recently devised algorithms for miRNA and target discovery. We do not discuss plant miRNAs because their varying sizes and structural characteristics pose different problems of identification and target selection.

Elevated expression of protein biosynthesis genes in liver and muscle of hibernating black bears (Ursus americanus).

We conducted a large-scale gene expression screen using the 3,200 cDNA probe microarray developed specifically for Ursus americanus to detect expression differences in liver and skeletal muscle that occur during winter hibernation compared with animals sampled during summer. The expression of 12 genes, including RNA binding protein motif 3 (Rbm3), that are mostly involved in protein biosynthesis, was induced during hibernation in both liver and muscle. The Gene Ontology and Gene Set Enrichment analysis consistently showed a highly significant enrichment of the protein biosynthesis category by overexpressed genes in both liver and skeletal muscle during hibernation. Coordinated induction in transcriptional level of genes involved in protein biosynthesis is a distinctive feature of the transcriptome in hibernating black bears. This finding implies induction of translation and suggests an adaptive mechanism that contributes to a unique ability to reduce muscle atrophy over prolonged periods of immobility during hibernation. Comparing expression profiles in bears to small mammalian hibernators shows a general trend during hibernation of transcriptional changes that include induction of genes involved in lipid metabolism and carbohydrate synthesis as well as depression of genes involved in the urea cycle and detoxification function in liver.

Learning from positive examples when the negative class is undetermined--microRNA gene identification.

BACKGROUND: The application of machine learning to classification problems that depend only on positive examples is gaining attention in the computational biology community. We and others have described the use of two-class machine learning to identify novel miRNAs. These methods require the generation of an artificial negative class. However, designation of the negative class can be problematic and if it is not properly done can affect the performance of the classifier dramatically and/or yield a biased estimate of performance. We present a study using one-class machine learning for microRNA (miRNA) discovery and compare one-class to two-class approaches using naïve Bayes and Support Vector Machines. These results are compared to published two-class miRNA prediction approaches. We also examine the ability of the one-class and two-class techniques to identify miRNAs in newly sequenced species. RESULTS: Of all methods tested, we found that 2-class naive Bayes and Support Vector Machines gave the best accuracy using our selected features and optimally chosen negative examples. One class methods showed average accuracies of 70-80% versus 90% for the two 2-class methods on the same feature sets. However, some one-class methods outperform some recently published two-class approaches with different selected features. Using the EBV genome as and external validation of the method we found one-class machine learning to work as well as or better than a two-class approach in identifying true miRNAs as well as predicting new miRNAs. CONCLUSION: One and two class methods can both give useful classification accuracies when the negative class is well characterized. The advantage of one class methods is that it eliminates guessing at the optimal features for the negative class when they are not well defined. In these cases one-class methods can be superior to two-class methods when the features which are chosen as representative of that positive class are well defined.

Naïve Bayes for microRNA target predictions--machine learning for microRNA targets.

MOTIVATION: Most computational methodologies for miRNA:mRNA target gene prediction use the seed segment of the miRNA and require cross-species sequence conservation in this region of the mRNA target. Methods that do not rely on conservation generate numbers of predictions, which are too large to validate. We describe a target prediction method (NBmiRTar) that does not require sequence conservation, using instead, machine learning by a naïve Bayes classifier. It generates a model from sequence and miRNA:mRNA duplex information from validated targets and artificially generated negative examples. Both the 'seed' and 'out-seed' segments of the miRNA:mRNA duplex are used for target identification. RESULTS: The application of machine-learning techniques to the features we have used is a useful and general approach for microRNA target gene prediction. Our technique produces fewer false positive predictions and fewer target candidates to be tested. It exhibits higher sensitivity and specificity than algorithms that rely on conserved genomic regions to decrease false positive predictions.

VISDA: an open-source caBIG analytical tool for data clustering and beyond.

SUMMARY: VISDA (Visual Statistical Data Analyzer) is a caBIG analytical tool for cluster modeling, visualization and discovery that has met silver-level compatibility under the caBIG initiative. Being statistically principled and visually interfaced, VISDA exploits both hierarchical statistics modeling and human gift for pattern recognition to allow a progressive yet interactive discovery of hidden clusters within high dimensional and complex biomedical datasets. The distinctive features of VISDA are particularly useful for users across the cancer research and broader research communities to analyze complex biological data. AVAILABILITY: http://gforge.nci.nih.gov/projects/visda/

Recursive cluster elimination (RCE) for classification and feature selection from gene expression data.

BACKGROUND: Classification studies using gene expression datasets are usually based on small numbers of samples and tens of thousands of genes. The selection of those genes that are important for distinguishing the different sample classes being compared, poses a challenging problem in high dimensional data analysis. We describe a new procedure for selecting significant genes as recursive cluster elimination (RCE) rather than recursive feature elimination (RFE). We have tested this algorithm on six datasets and compared its performance with that of two related classification procedures with RFE. RESULTS: We have developed a novel method for selecting significant genes in comparative gene expression studies. This method, which we refer to as SVM-RCE, combines K-means, a clustering method, to identify correlated gene clusters, and Support Vector Machines (SVMs), a supervised machine learning classification method, to identify and score (rank) those gene clusters for the purpose of classification. K-means is used initially to group genes into clusters. Recursive cluster elimination (RCE) is then applied to iteratively remove those clusters of genes that contribute the least to the classification performance. SVM-RCE identifies the clusters of correlated genes that are most significantly differentially expressed between the sample classes. Utilization of gene clusters, rather than individual genes, enhances the supervised classification accuracy of the same data as compared to the accuracy when either SVM or Penalized Discriminant Analysis (PDA) with recursive feature elimination (SVM-RFE and PDA-RFE) are used to remove genes based on their individual discriminant weights. CONCLUSION: SVM-RCE provides improved classification accuracy with complex microarray data sets when it is compared to the classification accuracy of the same datasets using either SVM-RFE or PDA-RFE. SVM-RCE identifies clusters of correlated genes that when considered together provide greater insight into the structure of the microarray data. Clustering genes for classification appears to result in some concomitant clustering of samples into subgroups. Our present implementation of SVM-RCE groups genes using the correlation metric. The success of the SVM-RCE method in classification suggests that gene interaction networks or other biologically relevant metrics that group genes based on functional parameters might also be useful.

A 10-gene classifier for distinguishing head and neck squamous cell carcinoma and lung squamous cell carcinoma.

PURPOSE: The risk of developing metastatic squamous cell carcinoma for patients with head and neck squamous cell carcinoma (HNSCC) is very high. Because these patients are often heavy tobacco users, they are also at risk for developing a second primary cancer, with squamous cell carcinoma of the lung (LSCC) being the most common. The distinction between a lung metastasis and a primary LSCC is currently based on certain clinical and histologic criteria, although the accuracy of this approach remains in question. EXPERIMENTAL DESIGN: Gene expression patterns derived from 28 patients with HNSCC or LSCC from a single center were analyzed using penalized discriminant analysis. Validation was done on previously published data for 134 total subjects from four independent Affymetrix data sets. RESULTS: We identified a panel of 10 genes (CXCL13, COL6A2, SFTPB, KRT14, TSPYL5, TMP3, KLK10, MMP1, GAS1, and MYH2) that accurately distinguished these two tumor types. This 10-gene classifier was validated on 122 subjects derived from four independent data sets and an average accuracy of 96% was shown. Gene expression values were validated by quantitative reverse transcription-PCR derived on 12 independent samples (seven HNSCC and five LSCC). The 10-gene classifier was also used to determine the site of origin of 12 lung lesions from patients with prior HNSCC. CONCLUSIONS: The results suggest that penalized discriminant analysis using these 10 genes will be highly accurate in determining the origin of squamous cell carcinomas in the lungs of patients with previous head and neck malignancies.