Contemporary high dimensional biological assays, such as mRNA expression microarrays, regularly involve multiple data processing steps, such as experimental processing, computational processing, sample selection, or feature selection (i.e. gene selection), prior to deriving any biological conclusions. These steps can dramatically change the interpretation of an experiment. Evaluation of processing steps has received limited attention in the literature. It is not straightforward to evaluate different processing methods and investigators are often unsure of the best method. We present a simple statistical tool, Standardized WithIn class Sum of Squares (SWISS), that allows investigators to compare alternate data processing methods, such as different experimental methods, normalizations, or technologies, on a dataset in terms of how well they cluster a priori biological classes. SWISS uses Euclidean distance to determine which method does a better job of clustering the data elements based on a priori classifications. We apply SWISS to three different gene expression applications. The first application uses four different datasets to compare different experimental methods, normalizations, and gene sets. The second application, using data from the MicroArray Quality Control (MAQC) project, compares different microarray platforms. The third application compares different technologies: a single Agilent two-color microarray versus one lane of RNA-Seq. These applications give an indication of the variety of problems that SWISS can be helpful in solving. The SWISS analysis of one-color versus two-color microarrays provides investigators who use two-color arrays the opportunity to review their results in light of a single-channel analysis, with all of the associated benefits offered by this design. Analysis of the MACQ data shows differential intersite reproducibility by array platform. SWISS also shows that one lane of RNA-Seq clusters data by biological phenotypes as well as a single Agilent two-color microarray.

Maspin is a tumor and metastasis suppressor playing an essential role as gatekeeper of tumor progression. It is highly expressed in epithelial cells but is silenced in the onset of metastatic disease by epigenetic mechanisms. Reprogramming of Maspin epigenetic silencing offers a therapeutic potential to lock metastatic progression. Herein we have investigated the ability of the Artificial Transcription Factor 126 (ATF-126) designed to upregulate the Maspin promoter to inhibit tumor progression in pre-established breast tumors in immunodeficient mice. ATF-126 was transduced in the aggressive, mesenchymal-like and triple negative breast cancer line, MDA-MB-231. Induction of ATF expression in vivo by Doxycycline resulted in 50% reduction in tumor growth and totally abolished tumor cell colonization. Genome-wide transcriptional profiles of ATF-induced cells revealed a gene signature that was found over-represented in estrogen receptor positive (ER+) "Normal-like" intrinsic subtype of breast cancer and in poorly aggressive, ER+ luminal A breast cancer cell lines. The comparison transcriptional profiles of ATF-126 and Maspin cDNA defined an overlapping 19-gene signature, comprising novel targets downstream the Maspin signaling cascade. Our data suggest that Maspin up-regulates downstream tumor and metastasis suppressor genes that are silenced in breast cancers, and are normally expressed in the neural system, including CARNS1, SLC8A2 and DACT3. In addition, ATF-126 and Maspin cDNA induction led to the re-activation of tumor suppressive miRNAs also expressed in neural cells, such as miR-1 and miR-34, and to the down-regulation of potential oncogenic miRNAs, such as miR-10b, miR-124, and miR-363. As expected from its over-representation in ER+ tumors, the ATF-126-gene signature predicted favorable prognosis for breast cancer patients. Our results describe for the first time an ATF able to reduce tumor growth and metastatic colonization by epigenetic reactivation of a dormant, normal-like, and more differentiated gene program.

Relief-F is a nonparametric, nearest-neighbor machine learning method that has been successfully used to identify relevant variables that may interact in complex multivariate models to explain phenotypic variation. While several tools have been developed for assessing differential expression in sequence-based transcriptomics, the detection of statistical interactions between transcripts has received less attention in the area of RNA-seq analysis. We describe a new extension and assessment of Relief-F for feature selection in RNA-seq data. The ReliefSeq implementation adapts the number of nearest neighbors (k) for each gene to optimize the Relief-F test statistics (importance scores) for finding both main effects and interactions. We compare this gene-wise adaptive-k (gwak) Relief-F method with standard RNA-seq feature selection tools, such as DESeq and edgeR, and with the popular machine learning method Random Forests. We demonstrate performance on a panel of simulated data that have a range of distributional properties reflected in real mRNA-seq data including multiple transcripts with varying sizes of main effects and interaction effects. For simulated main effects, gwak-Relief-F feature selection performs comparably to standard tools DESeq and edgeR for ranking relevant transcripts. For gene-gene interactions, gwak-Relief-F outperforms all comparison methods at ranking relevant genes in all but the highest fold change/highest signal situations where it performs similarly. The gwak-Relief-F algorithm outperforms Random Forests for detecting relevant genes in all simulation experiments. In addition, Relief-F is comparable to the other methods based on computational time. We also apply ReliefSeq to an RNA-Seq study of smallpox vaccine to identify gene expression changes between vaccinia virus-stimulated and unstimulated samples. ReliefSeq is an attractive tool for inclusion in the suite of tools used for analysis of mRNA-Seq data; it has power to detect both main effects and interaction effects. Software Availability: http://insilico.utulsa.edu/ReliefSeq.php.

Failure to achieve a protected state after influenza vaccination is poorly understood but occurs commonly among aged populations experiencing greater immunosenescence. In order to better understand immune response in the elderly, we studied epigenetic and transcriptomic profiles and humoral immune response outcomes in 50-74 year old healthy participants. Associations between DNA methylation and gene expression reveal a system-wide regulation of immune-relevant functions, likely playing a role in regulating a participant's propensity to respond to vaccination. Our findings show that sites of methylation regulation associated with humoral response to vaccination impact known cellular differentiation signaling and antigen presentation pathways. We performed our analysis using per-site and regionally average methylation levels, in addition to continuous or dichotomized outcome measures. The genes and molecular functions implicated by each analysis were compared, highlighting different aspects of the biologic mechanisms of immune response affected by differential methylation. Both cis-acting (within the gene or promoter) and trans-acting (enhancers and transcription factor binding sites) sites show significant associations with measures of humoral immunity. Specifically, we identified a group of CpGs that, when coordinately hypo-methylated, are associated with lower humoral immune response, and methylated with higher response. Additionally, CpGs that individually predict humoral immune responses are enriched for polycomb-group and FOXP2 transcription factor binding sites. The most robust associations implicate differential methylation affecting gene expression levels of genes with known roles in immunity (e.g. HLA-B and HLA-DQB2) and immunosenescence. We believe our data and analysis strategy highlight new and interesting epigenetic trends affecting humoral response to vaccination against influenza; one of the most common and impactful viral pathogens.

Pancreatic adenocarcinoma (PaC) is one of most difficult tumors to treat. Much of this is attributed to the late diagnosis. To identify biomarkers for early detection, we examined DNA methylation differences in leukocyte DNA between PaC cases and controls in a two-phase study. In phase I, we measured methylation levels at 1,505 CpG sites in treatment-naïve leukocyte DNA from 132 never-smoker PaC patients and 60 never-smoker healthy controls. We found significant differences in 110 CpG sites (false discovery rate <0.05). In phase II, we tested and validated 88 of 96 phase I selected CpG sites in 240 PaC cases and 240 matched controls (p≤0.05). Using penalized logistic regression, we built a prediction model consisting of five CpG sites (IL10_P348, LCN2_P86, ZAP70_P220, AIM2_P624, TAL1_P817) that discriminated PaC patients from controls (C-statistic = 0.85 in phase I; 0.76 in phase II). Interestingly, one CpG site (LCN2_P86) alone could discriminate resectable patients from controls (C-statistic= 0.78 in phase I; 0.74 in phase II). We also performed methylation quantitative trait loci (methQTL) analysis and identified three CpG sites (AGXT_P180_F, ALOX12_E85_R, JAK3_P1075_R) where the methylation levels were significantly associated with single nucleotide polymorphisms (SNPs) (false discovery rate <0.05). Our results demonstrate that epigenetic variation in easily obtainable leukocyte DNA, manifested by reproducible methylation differences, may be used to detect PaC patients. The methylation differences at certain CpG sites are partially attributable to genetic variation. This study strongly supports future epigenome-wide association study using leukocyte DNA for biomarker discovery in human diseases.

MicroRNAs play a role in regulating diverse biological processes and have considerable utility as molecular markers for diagnosis and monitoring of human disease. Several technologies are available commercially for measuring microRNA expression. However, cross-platform comparisons do not necessarily correlate well, making it difficult to determine which platform most closely represents the true microRNA expression level in a tissue. To address this issue, we have analyzed RNA derived from cell lines, as well as fresh frozen and formalin-fixed paraffin embedded tissues, using Affymetrix, Agilent, and Illumina microRNA arrays, NanoString counting, and Illumina Next Generation Sequencing. We compared the performance within- and between the different platforms, and then verified these results with those of quantitative PCR data. Our results demonstrate that the within-platform reproducibility for each method is consistently high and although the gene expression profiles from each platform show unique traits, comparison of genes that were commonly detectable showed that detection of microRNA transcripts was similar across multiple platforms.

There are insufficient system-wide transcriptomic (or other) data that help explain the observed inter-individual variability in antibody titers after measles vaccination in otherwise healthy individuals.

Expression level of many genes shows abundant natural variation in human populations. The variations in gene expression are believed to contribute to phenotypic differences. Emerging evidence has shown that microRNAs (miRNAs) are one of the key regulators of gene expression. However, past studies have focused on the miRNA target genes and used loss- or gain-of-function approach that may not reflect natural association between miRNA and mRNAs.

During the past decade, the age-adjusted mortality rate for endometrial cancer (EC) increased 1.9% annually with TP53 mutant (TP53-mu) EC disproportionally represented in advanced disease and deaths. Therefore, we aimed to assess pivotal molecular parameters that differentiate clinical outcomes in high- and low-risk EC. Using the Cancer Genome Atlas, we analyzed EC specimens with available DNA sequences and quantitative gene-specific RNA expression data. After polymerase ɛ (POLE)-mutant specimens were excluded, differential gene-specific mutations and mRNA expressions were annotated and integrated. Consequent to TP53-mu failure to induce p21, derepression of multiple oncogenes harboring promoter p21 repressive sites was observed, including CCNA2 and FOXM1 (P < .001 compared with TP53 wild type [TP53-wt]). TP53-wt EC with high CCNA2 expression (CCNA2-H) had a targeted transcriptomic profile similar to that of TP53-mu EC, suggesting CCNA2 is a seminal determinant for both TP53-wt and TP53-mu EC. CCNA2 enhances E2F1 function, upregulating FOXM1 and CIP2A, as observed in TP53-mu and CCNA2-H TP53-wt EC (P < .001). CIP2A inhibits protein phosphatase 2A, leading to AKT inactivation of GSK3β and restricted oncoprotein degradation; PPP2R1A and FBXW7 mutations yield similar results. Upregulation of FOXM1 and failed degradation of FOXM1 is evidenced by marked upregulation of multiple homologous recombination genes (P < .001). Integrating these molecular aberrations generated a molecular biomarker panel with significant prognostic discrimination (P = 5.8×10-7); adjusting for age, histology, grade, myometrial invasion, TP53 status, and stage, only CCNA2-H/E2F1-H (P = .0003), FBXW7-mu/PPP2R1A-mu (P = .0002), and stage (P = .017) were significant. The generated prognostic molecular classification system identifies dissimilar signaling aberrations potentially amenable to targetable therapeutic options.

Ovarian carcinosarcoma is a rare subtype of ovarian cancer with poor clinical outcomes. The low incidence of this disease makes accrual to large clinical trials challenging. However, studies have shown that treatment responses in patient-derived xenograft (PDX) models correlate with matched-patient responses in the clinic, supporting their use for preclinical testing of standard and novel therapies. An ovarian carcinosarcoma PDX is presented herein and showed resistance to carboplatin and paclitaxel (similar to the patient) but exhibited significant sensitivity to ifosfamide and paclitaxel. The PDX demonstrated overexpression of EGFR mRNA and gene amplification by array comparative genomic hybridization (log2 ratio 0.399). EGFR phosphorylation was also detected. Angiogensis and insulin-like growth factor pathways were also implicated by overexpression of VEGFC and IRS1. In order to improve response to chemotherapy, the PDX was treated with carboplatin/paclitaxel with or without a pan-HER and VEGF inhibitor (BMS-690514) but there was no tumor growth inhibition or improved animal survival, which may be explained by a KRAS mutation. Resistance was also observed when the IGF-1R inhibitor BMS-754807 was combined with carboplatin/paclitaxel. Because poly (ADP-ribose) polymerase inhibitors have activity in ovarian cancer patients, with and without BRCA mutations, ABT-888 was also tested but found to have no activity. Pathogenic mutations were also detected in TP53 and PIK3CA. In conclusion, ifosfamide/paclitaxel was superior to carboplatin/paclitaxel in this ovarian carcinosarcoma PDX and gene overexpression or amplification alone was not sufficient to predict response to targeted therapy. Better predictive markers of response are needed.

Although ovarian cancer is the most lethal of gynecologic malignancies, wide variation in outcome following conventional therapy continues to exist. The presence of tumor-infiltrating regulatory T cells (Tregs) has a role in outcome of this disease, and a growing body of data supports the existence of inherited prognostic factors. However, the role of inherited variants in genes encoding Treg-related immune molecules has not been fully explored. We analyzed expression quantitative trait loci (eQTL) and sequence-based tagging single nucleotide polymorphisms (tagSNPs) for 54 genes associated with Tregs in 3,662 invasive ovarian cancer cases. With adjustment for known prognostic factors, suggestive results were observed among rarer histological subtypes; poorer survival was associated with minor alleles at SNPs in RGS1 (clear cell, rs10921202, p=2.7×10(-5)), LRRC32 and TNFRSF18/TNFRSF4 (mucinous, rs3781699, p=4.5×10(-4), and rs3753348, p=9.0×10(-4), respectively), and CD80 (endometrioid, rs13071247, p=8.0×10(-4)). Fo0r the latter, correlative data support a CD80 rs13071247 genotype association with CD80 tumor RNA expression (p=0.006). An additional eQTL SNP in CD80 was associated with shorter survival (rs7804190, p=8.1×10(-4)) among all cases combined. As the products of these genes are known to affect induction, trafficking, or immunosuppressive function of Tregs, these results suggest the need for follow-up phenotypic studies.

Patients with obstructive sleep apnea (OSA) experience partial or complete upper airway collapses during sleep resulting in nocturnal hypoxia-normoxia cycling, and continuous positive airway pressure (CPAP) is the golden treatment for OSA. Nevertheless, the exact mechanisms of action, especially the transcriptome effect of CPAP on OSA patients, remain elusive. The goal of this study was to evaluate the longitudinal alterations in peripheral blood mononuclear cells transcriptome profiles of OSA patients in order to identify the hub gene and immune response. GSE133601 was downloaded from Gene Expression Omnibus (GEO). We identified black module via weighted gene co-expression network analysis (WGCNA), the genes in which were correlated significantly with the clinical trait of CPAP treatment. Finally, eleven hub genes (TRAV10, SNORA36A, RPL10, OBP2B, IGLV1-40, H2BC8, ESAM, DNASE1L3, CD22, ANK3, ACP3) were traced and used to construct a random forest model to predict therapeutic efficacy of CPAP in OSA with a good performance with AUC of 0.92. We further studied the immune cells infiltration in OSA patients with CIBERSORT, and monocytes were found to be related with the remission of OSA and partially correlated with the hub genes identified. In conclusion, these key genes and immune infiltration may be of great importance in the remission of OSA and related research of these genes may provide a new therapeutic target for OSA in the future.

Although influenza causes significant morbidity and mortality in the elderly, the factors underlying the reduced vaccine immunogenicity and efficacy in this age group are not completely understood. Age and immunosenescence factors, and their impact on humoral immunity after influenza vaccination, are of growing interest for the development of better vaccines for the elderly.

Therapies targeting the type I insulin-like growth factor receptor (IGF-1R) have not been developed with predictive biomarkers to identify tumors with receptor activation. We have previously shown that the insulin receptor substrate (IRS) adaptor proteins are necessary for linking IGF1R to downstream signaling pathways and the malignant phenotype in breast cancer cells. The purpose of this study was to identify gene expression profiles downstream of IGF1R and its two adaptor proteins. IRS-null breast cancer cells (T47D-YA) were engineered to express IRS-1 or IRS-2 alone and their ability to mediate IGF ligand-induced proliferation, motility, and gene expression determined. Global gene expression signatures reflecting IRS adaptor specific and primary vs. secondary ligand response were derived (Early IRS-1, Late IRS-1, Early IRS-2 and Late IRS-2) and functional pathway analysis examined. IRS isoforms mediated distinct gene expression profiles, functional pathways, and breast cancer subtype association. For example, IRS-1/2-induced TGFb2 expression and blockade of TGFb2 abrogated IGF-induced cell migration. In addition, the prognostic value of IRS proteins was significant in the luminal B breast tumor subtype. Univariate and multivariate analyses confirmed that IRS adaptor signatures correlated with poor outcome as measured by recurrence-free and overall survival. Thus, IRS adaptor protein expression is required for IGF ligand responses in breast cancer cells. IRS-specific gene signatures represent accurate surrogates of IGF activity and could predict response to anti-IGF therapy in breast cancer.

Immune responses to current rubella vaccines demonstrate significant inter-individual variability. We performed mRNA-Seq profiling on PBMCs from high and low antibody responders to rubella vaccination to delineate transcriptional differences upon viral stimulation. Generalized linear models were used to assess the per gene fold change (FC) for stimulated versus unstimulated samples or the interaction between outcome and stimulation. Model results were evaluated by both FC and p-value. Pathway analysis and self-contained gene set tests were performed for assessment of gene group effects. Of 17,566 detected genes, we identified 1,080 highly significant differentially expressed genes upon viral stimulation (p<1.00E(-15), FDR<1.00E(-14)), including various immune function and inflammation-related genes, genes involved in cell signaling, cell regulation and transcription, and genes with unknown function. Analysis by immune outcome and stimulation status identified 27 genes (p≤0.0006 and FDR≤0.30) that responded differently to viral stimulation in high vs. low antibody responders, including major histocompatibility complex (MHC) class I genes (HLA-A, HLA-B and B2M with p = 0.0001, p = 0.0005 and p = 0.0002, respectively), and two genes related to innate immunity and inflammation (EMR3 and MEFV with p = 1.46E(-08) and p = 0.0004, respectively). Pathway and gene set analysis also revealed transcriptional differences in antigen presentation and innate/inflammatory gene sets and pathways between high and low responders. Using mRNA-Seq genome-wide transcriptional profiling, we identified antigen presentation and innate/inflammatory genes that may assist in explaining rubella vaccine-induced immune response variations. Such information may provide new scientific insights into vaccine-induced immunity useful in rational vaccine development and immune response monitoring.

Inhibition of pregnancy-associated plasma protein-A (PAPP-A), an upstream activator of the insulin-like growth factor (IGF) pathway, is known to augment sensitivity to platinum-based chemotherapy. This study further tests the efficacy of PAPP-A inhibition with a monoclonal antibody inhibitor (mAb-PA) in ovarian cancer (OC) platinum-resistant patient-derived xenograft (PDX) models.