The cDNA-mediated Annealing, extension, Selection and Ligation (DASL) assay has become a suitable gene expression profiling system for degraded RNA from paraffin-embedded tissue. We examined assay characteristics and the performance of the DASL 502-gene Cancer Panel v1 (1.5K) and 24,526-gene panel (24K) platforms at differentiating nine human epidermal growth factor receptor 2- positive (HER2+) and 11 HER2-negative (HER2-) paraffin-embedded breast tumors.

MicroRNAs (miRNAs) are a class of approximately 22 nucleotide long, widely expressed RNA molecules that play important regulatory roles in eukaryotes. To investigate miRNA function, it is essential that methods to quantify their expression levels be available.

Mammalian target of rapamycin (mTOR) represents a key downstream intermediate for a myriad of oncogenic receptor tyrosine kinases. In the case of the insulin-like growth factor (IGF) pathway, the mTOR complex (mTORC1) mediates IGF-1 receptor (IGF-1R)-induced estrogen receptor alpha (ERα) phosphorylation/activation and leads to increased proliferation and growth in breast cancer cells. As a result, the prevalence of mTOR inhibitors combined with hormonal therapy has increased in recent years. Conversely, activated mTORC1 provides negative feedback regulation of IGF signaling via insulin receptor substrate (IRS)-1/2 serine phosphorylation and subsequent proteasomal degradation. Thus, the IGF pathway may provide escape (e.g. de novo or acquired resistance) from mTORC1 inhibitors. It is therefore plausible that combined inhibition of mTORC1 and IGF-1R for select subsets of ER-positive breast cancer patients presents as a viable therapeutic option.

Recent heritability analyses have indicated that genome-wide association studies (GWAS) have the potential to improve genetic risk prediction for complex diseases based on polygenic risk score (PRS), a simple modelling technique that can be implemented using summary-level data from the discovery samples. We herein propose modifications to improve the performance of PRS. We introduce threshold-dependent winner's-curse adjustments for marginal association coefficients that are used to weight the single-nucleotide polymorphisms (SNPs) in PRS. Further, as a way to incorporate external functional/annotation knowledge that could identify subsets of SNPs highly enriched for associations, we propose variable thresholds for SNPs selection. We applied our methods to GWAS summary-level data of 14 complex diseases. Across all diseases, a simple winner's curse correction uniformly led to enhancement of performance of the models, whereas incorporation of functional SNPs was beneficial only for selected diseases. Compared to the standard PRS algorithm, the proposed methods in combination led to notable gain in efficiency (25-50% increase in the prediction R2) for 5 of 14 diseases. As an example, for GWAS of type 2 diabetes, winner's curse correction improved prediction R2 from 2.29% based on the standard PRS to 3.10% (P = 0.0017) and incorporating functional annotation data further improved R2 to 3.53% (P = 2×10-5). Our simulation studies illustrate why differential treatment of certain categories of functional SNPs, even when shown to be highly enriched for GWAS-heritability, does not lead to proportionate improvement in genetic risk-prediction because of non-uniform linkage disequilibrium structure.

Urinary exosome-like vesicles (ELVs) are a heterogenous mixture (diameter 40-200 nm) containing vesicles shed from all segments of the nephron including glomerular podocytes. Contamination with Tamm-Horsfall protein (THP) oligomers has hampered their isolation and proteomic analysis. Here we improved ELV isolation protocols employing density centrifugation to remove THP and albumin, and isolated a glomerular membranous vesicle (GMV)-enriched subfraction from 7 individuals identifying 1830 proteins and in 3 patients with glomerular disease identifying 5657 unique proteins. The GMV fraction was composed of podocin/podocalyxin-positive irregularly shaped membranous vesicles and podocin/podocalyxin-negative classical exosomes. Ingenuity pathway analysis identified integrin, actin cytoskeleton, and Rho GDI signaling in the top three canonical represented signaling pathways and 19 other proteins associated with inherited glomerular diseases. The GMVs are of podocyte origin and the density gradient technique allowed isolation in a reproducible manner. We show many nephrotic syndrome proteins, proteases, and complement proteins involved in glomerular disease are in GMVs and some were only shed in the disease state (nephrin, TRPC6, INF2 and phospholipase A2 receptor). We calculated sample sizes required to identify new glomerular disease biomarkers, expand the ELV proteome, and provide a reference proteome in a database that may prove useful in the search for biomarkers of glomerular disease.

Ovarian cancer represents the most lethal tumor type among malignancies of the female reproductive system. Overall survival rates remain low. In this study, we identify the serine protease inhibitor Kazal type 1 (SPINK1) as a potential therapeutic target for a subset of ovarian cancers. We show that SPINK1 drives ovarian cancer cell proliferation through activation of epidermal growth factor receptor (EGFR) signaling, and that SPINK1 promotes resistance to anoikis through a distinct mechanism involving protease inhibition. In analyses of ovarian tumor specimens from a Mayo Clinic cohort of 490 patients, we further find that SPINK1 immunostaining represents an independent prognostic factor for poor survival, with the strongest association in patients with nonserous histological tumor subtypes (endometrioid, clear cell, and mucinous). This study provides novel insight into the fundamental processes underlying ovarian cancer progression, and also suggests new avenues for development of molecularly targeted therapies.

Growing evidence suggests a major role for Src-homology-2-domain-containing phosphatase 2 (SHP2/PTPN11) in MYCN-driven high-risk neuroblastoma, although biologic confirmation and a plausible mechanism for this contribution are lacking. Using a zebrafish model of MYCN-overexpressing neuroblastoma, we demonstrate that mutant ptpn11 expression in the adrenal gland analog of MYCN transgenic fish promotes the proliferation of hyperplastic neuroblasts, accelerates neuroblastomagenesis, and increases tumor penetrance. We identify a similar mechanism in tumors with wild-type ptpn11 and dysregulated Gab2, which encodes a Shp2 activator that is overexpressed in human neuroblastomas. In MYCN transgenic fish, Gab2 overexpression activated the Shp2-Ras-Erk pathway, enhanced neuroblastoma induction, and increased tumor penetrance. We conclude that MYCN cooperates with either GAB2-activated or mutant SHP2 in human neuroblastomagenesis. Our findings further suggest that combined inhibition of MYCN and the SHP2-RAS-ERK pathway could provide effective targeted therapy for high-risk neuroblastoma patients with MYCN amplification and aberrant SHP2 activation.

Many patients with pancreatic adenocarcinoma carry germline mutations associated with increased risk of cancer. It is not clear whether patients with intraductal papillary mucinous neoplasms (IPMNs), which are precursors to some pancreatic cancers, also carry these mutations. We assessed the prevalence of germline mutations associated with cancer risk in patients with histologically confirmed IPMN.

Therapy with demethylating agent 5-azacitidine and histone deacetylase inhibitor entinostat shows synergistic re-expression of tumor-suppressor genes and growth inhibition in colorectal (CRC) cell lines and in vivo studies.

Pancreatic adenocarcinoma (PC) is a lethal malignancy that is familial or associated with genetic syndromes in 10% of cases. Gene-based surveillance strategies for at-risk individuals may improve clinical outcomes. However, familial PC (FPC) is plagued by genetic heterogeneity and the genetic basis for the majority of FPC remains elusive, hampering the development of gene-based surveillance programs. The study was powered to identify genes with a cumulative pathogenic variant prevalence of at least 3%, which includes the most prevalent PC susceptibility gene, BRCA2. Since the majority of known PC susceptibility genes are involved in DNA repair, we focused on genes implicated in these pathways. We performed a region-based association study using the Mixed-Effects Score Test, followed by leave-one-out characterization of PC-associated gene regions and variants to identify the genes and variants driving risk associations. We evaluated 398 cases from two case series and 987 controls without a personal history of cancer. The first case series consisted of 109 patients with either FPC (n = 101) or PC at ≤50 years of age (n = 8). The second case series was composed of 289 unselected PC cases. We validated this discovery strategy by identifying known pathogenic BRCA2 variants, and also identified SMG1, encoding a serine/threonine protein kinase, to be significantly associated with PC following correction for multiple testing (p = 3.22x10-7). The SMG1 association was validated in a second independent series of 532 FPC cases and 753 controls (p<0.0062, OR = 1.88, 95%CI 1.17-3.03). We showed segregation of the c.4249A>G SMG1 variant in 3 affected relatives in a FPC kindred, and we found c.103G>A to be a recurrent SMG1 variant associating with PC in both the discovery and validation series. These results suggest that SMG1 is a novel PC susceptibility gene, and we identified specific SMG1 gene variants associated with PC risk.

Although 70-80% of newly diagnosed ovarian cancer patients respond to first-line therapy, almost all relapse and five-year survival remains below 50%. One strategy to increase five-year survival is prolonging time to relapse by improving first-line therapy response. However, no biomarker today can accurately predict individual response to therapy. In this study, we present analytical and prospective clinical validation of a new test that utilizes primary patient tissue in 3D cell culture to make patient-specific response predictions prior to initiation of treatment in the clinic. Test results were generated within seven days of tissue receipt from newly diagnosed ovarian cancer patients obtained at standard surgical debulking or laparoscopic biopsy. Patients were followed for clinical response to chemotherapy. In a study population of 44, the 32 test-predicted Responders had a clinical response rate of 100% across both adjuvant and neoadjuvant treated populations with an overall prediction accuracy of 89% (39 of 44, p < 0.0001). The test also functioned as a prognostic readout with test-predicted Responders having a significantly increased progression-free survival compared to test-predicted Non-Responders, p = 0.01. This correlative accuracy establishes the test's potential to benefit ovarian cancer patients through accurate prediction of patient-specific response before treatment.

Germline copy number variants (CNVs) increase risk for many diseases, yet detection of CNVs and quantifying their contribution to disease risk in large-scale studies is challenging due to biological and technical sources of heterogeneity that vary across the genome within and between samples.

Interest in preclinical drug development for ovarian cancer has stimulated development of patient-derived xenograft (PDX) or tumorgraft models. However, the unintended formation of human lymphoma in severe combined immunodeficiency (SCID) mice from Epstein-Barr virus (EBV)-infected human lymphocytes can be problematic. In this study, we have characterized ovarian cancer PDXs which developed human lymphomas and explore methods to suppress lymphoproliferative growth. Fresh human ovarian tumors from 568 patients were transplanted intraperitoneally in SCID mice. A subset of PDX models demonstrated atypical patterns of dissemination with mediastinal masses, hepatosplenomegaly, and CD45-positive lymphoblastic atypia without ovarian tumor engraftment. Expression of human CD20 but not CD3 supported a B-cell lineage, and EBV genomes were detected in all lymphoproliferative tumors. Immunophenotyping confirmed monoclonal gene rearrangements consistent with B-cell lymphoma, and global gene expression patterns correlated well with other human lymphomas. The ability of rituximab, an anti-CD20 antibody, to suppress human lymphoproliferation from a patient's ovarian tumor in SCID mice and prevent growth of an established lymphoma led to a practice change with a goal to reduce the incidence of lymphomas. A single dose of rituximab during the primary tumor heterotransplantation process reduced the incidence of CD45-positive cells in subsequent PDX lines from 86.3% (n = 117 without rituximab) to 5.6% (n = 160 with rituximab), and the lymphoma rate declined from 11.1% to 1.88%. Taken together, investigators utilizing PDX models for research should routinely monitor for lymphoproliferative tumors and consider implementing methods to suppress their growth.

Pancreatic cancer (PC) is a complex disease in which both non-genetic and genetic factors interplay. To date, 40 GWAS hits have been associated with PC risk in individuals of European descent, explaining 4.1% of the phenotypic variance.

Germline variation and smoking are independently associated with pancreatic ductal adenocarcinoma (PDAC). We conducted genome-wide smoking interaction analysis of PDAC using genotype data from four previous genome-wide association studies in individuals of European ancestry (7,937 cases and 11,774 controls). Examination of expression quantitative trait loci data from the Genotype-Tissue Expression Project followed by colocalization analysis was conducted to determine whether there was support for common SNP(s) underlying the observed associations. Statistical tests were two sided and P < 5 × 10-8 was considered statistically significant. Genome-wide significant evidence of qualitative interaction was identified on chr2q21.3 in intron 5 of the transmembrane protein 163 (TMEM163) and upstream of the cyclin T2 (CCNT2). The most significant SNP using the Empirical Bayes method, in this region that included 45 significantly associated SNPs, was rs1818613 [per allele OR in never smokers 0.87, 95% confidence interval (CI), 0.82-0.93; former smokers 1.00, 95% CI, 0.91-1.07; current smokers 1.25, 95% CI 1.12-1.40, P interaction = 3.08 × 10-9). Examination of the Genotype-Tissue Expression Project data demonstrated an expression quantitative trait locus in this region for TMEM163 and CCNT2 in several tissue types. Colocalization analysis supported a shared SNP, rs842357, in high linkage disequilibrium with rs1818613 (r 2 = 0. 94) driving both the observed interaction and the expression quantitative trait loci signals. Future studies are needed to confirm and understand the differential biologic mechanisms by smoking status that contribute to our PDAC findings. SIGNIFICANCE: This large genome-wide interaction study identifies a susceptibility locus on 2q21.3 that significantly modified PDAC risk by smoking status, providing insight into smoking-associated PDAC, with implications for prevention.

Pancreatic ductal adenocarcinoma (PDAC) is a lethal disease with a 5-year survival rate of less than 10%. Individuals with a pathogenic germline variant in a pancreatic cancer susceptibility gene are at an increased risk of developing pancreatic cancer. Understanding the inherited genetic basis of pancreatic tumor development provides a unique opportunity to improve patient care and outcomes. For example, relatives of a patients with PDAC who have a pathogenic germline variant in a pancreatic cancer susceptibility gene are eligible for disease surveillance where cancers may be detected early, and 5-year survival greatly improved. Furthermore, for some patients with PDAC and a pathogenic germline variant in a pancreatic cancer susceptibility gene, their tumors may be susceptible to specific anti-cancer therapies. Recently, RABL3 was identified as a pancreatic cancer susceptibility gene. To validate these findings and inform clinical translation, we determined the prevalence of deleterious RABL3 variants in a large cohort of 1037 patients with PDAC that had undergone either whole genome or whole exome germline sequencing. We identified two synonymous variants and four missense variants classified as variants of unknown significance. We found no pathogenic RABL3 variants, indicating that the maximum prevalence of such variants in patients with PDAC is less than 0.36% (minor allele frequency 0, 97.5% one-sided confidence interval: 0-0.0036). This finding has important implications for germline genetic testing of patients with PDAC.

The COVID-19 pandemic has revealed the pronounced vulnerability of the elderly and chronically ill to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-induced morbidity and mortality. Cellular senescence contributes to inflammation, multiple chronic diseases, and age-related dysfunction, but effects on responses to viral infection are unclear. Here, we demonstrate that senescent cells (SnCs) become hyper-inflammatory in response to pathogen-associated molecular patterns (PAMPs), including SARS-CoV-2 spike protein-1, increasing expression of viral entry proteins and reducing antiviral gene expression in non-SnCs through a paracrine mechanism. Old mice acutely infected with pathogens that included a SARS-CoV-2-related mouse β-coronavirus experienced increased senescence and inflammation, with nearly 100% mortality. Targeting SnCs by using senolytic drugs before or after pathogen exposure significantly reduced mortality, cellular senescence, and inflammatory markers and increased antiviral antibodies. Thus, reducing the SnC burden in diseased or aged individuals should enhance resilience and reduce mortality after viral infection, including that of SARS-CoV-2.

Of clinically node-negative (cN0) cutaneous melanoma patients with sentinel lymph node (SLN) metastasis, between 10% and 30% harbor additional metastases in non-sentinel lymph nodes (NSLNs). Approximately 80% of SLN-positive patients have a single positive SLN.

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.