NanoStringNorm: an extensible R package for the pre-processing of NanoString mRNA and miRNA data

Clearance of pegylated interferon by Kupffer cells limits NK cell activation and therapy response of patients with HBV infection

Pegylated interferon-α (PEG-IFN-α), where IFN-α is attached to polyethylene glycol (PEG), is an approved treatment for chronic hepatitis B virus (HBV) infection, a disease that causes liver-related morbidity and mortality in 257 million people worldwide. It is unknown why only a minority of patients respond to PEG-IFN-α. Using sequential blood samples and liver biopsies of patients with chronic HBV infection before, during, and after PEG-IFN-α treatment, we find that patients with early natural killer (NK) cell activation after PEG-IFN-α injection experienced greater liver inflammation, lysis of HBV-infected hepatocytes, and hepatitis B surface antigen (HBsAg) decline than those without. NK cell activation was associated with induction of interferon-stimulated genes and determined by PEG-IFN-α pharmacokinetics. Patients with delayed increases in PEG-IFN-α concentrations had greater amounts of PEG-specific immunoglobulin M (IgM) immune complexes in the blood and more PEG and IgM detected in the liver than patients with rapid increase in PEG-IFN-α concentration. This was associated with reduced NK cell activation. These results indicate that the immunomodulatory functions of PEG-IFN-α, particularly activation of NK cells, play a pivotal role in the response to treatment and further demonstrate that these functions are affected by PEG-IFN-α pharmacokinetics. Accelerated clearance of antibody-complexed pegylated drugs by Kupffer cells may be important beyond the field of HBV therapeutics. Thus, these findings may contribute to improving the efficacy of pegylated drugs that are now being developed for other chronic diseases and cancer.

Immune Signature Linked to COVID-19 Severity: A SARS-Score for Personalized Medicine

SARS-CoV-2 infection leads to a highly variable clinical evolution, ranging from asymptomatic to severe disease with acute respiratory distress syndrome, requiring intensive care units (ICU) admission. The optimal management of hospitalized patients has become a worldwide concern and identification of immune biomarkers predictive of the clinical outcome for hospitalized patients remains a major challenge. Immunophenotyping and transcriptomic analysis of hospitalized COVID-19 patients at admission allow identifying the two categories of patients. Inflammation, high neutrophil activation, dysfunctional monocytic response and a strongly impaired adaptive immune response was observed in patients who will experience the more severe form of the disease. This observation was validated in an independent cohort of patients. Using in silico analysis on drug signature database, we identify differential therapeutics that specifically correspond to each group of patients. From this signature, we propose a score-the SARS-Score-composed of easily quantifiable biomarkers, to classify hospitalized patients upon arrival to adapt treatment according to their immune profile.

Translation of a circulating miRNA signature of melanoma into a solid tissue assay to improve diagnostic accuracy and precision

Aim: Successful treatment of cutaneous melanoma depends on early and accurate diagnosis of clinically suspicious melanocytic skin lesions. Multiple international studies have described the challenge of providing accurate and reproducible histopathological assessments of melanocytic lesions, highlighting the need for new diagnostic tools including disease-specific biomarkers. Previously, a 38-miRNA signature (MEL38) was identified in melanoma patient plasma and validated as a novel biomarker. In this study, MEL38 expression in solid tissue biopsies representing the benign nevi to metastatic melanoma spectrum is examined. Patients & methods: Nanostring digital gene expression assessment of the MEL38 signature was performed on 308 formalin-fixed paraffin-embedded biopsies of nevi, melanoma in situ and invasive melanoma. Genomic data were interrogated using hierarchical clustering, univariate and multivariate statistical approaches. Classification scores computed from the MEL38 signature were analyzed for their association with demographic data and histopathology results, including MPATH-DX class, AJCC disease stage and tissue subtype. Results: The MEL38 score can stratify higher-risk melanomas (MPATH-Dx class V or more advanced) from lower-risk skin lesions (class I-IV) with an area under the curve of 0.97 (p < 0.001). The genomic score ranges from 0 to 10 and is positively correlated with melanoma progression, with an intraclass correlation coefficient of 0.85 with stage 0-IV disease. Using an optimized classification threshold of ≥2.7 accurately identifies higher-risk melanomas with 89% sensitivity and 94% specificity. Multivariate analysis showed the score to be a significant predictor of malignancy, independent of technical and clinical covariates. Application of the MEL38 signature to Spitz nevi reveals an intrasubtype profile, with elements in common to both nevi and melanoma. Conclusion: Melanoma-specific circulating miRNAs maintain their association with malignancy when measured in the hypothesized tissue of origin. The MEL38 signature is an accurate and reproducible metric of melanoma status, based on changes in miRNA expression that occur as the disease develops and spreads. Inclusion of the MEL38 score into routine practice would provide physicians with previously unavailable, personalized genomic information about their patient's skin lesions. Combining molecular biomarker data with conventional histopathology data may improve diagnostic accuracy, healthcare resource utilization and patient outcomes.

The Colonic Mucosal MicroRNAs, MicroRNA-219a-5p, and MicroRNA-338-3p Are Downregulated in Irritable Bowel Syndrome and Are Associated With Barrier Function and MAPK Signaling

BACKGROUND & AIMS

Alterations in microRNA (miRNA) and in the intestinal barrier are putative risk factors for irritable bowel syndrome (IBS). We aimed to identify differentially expressed colonic mucosal miRNAs, their targets in IBS compared to healthy controls (HCs), and putative downstream pathways.

METHODS

Twenty-nine IBS patients (15 IBS with constipation [IBS-C], 14 IBS with diarrhea [IBS-D]), and 15 age-matched HCs underwent sigmoidoscopy with biopsies. A nCounter array was used to assess biopsy specimen-associated miRNA levels. A false discovery rate (FDR) < 10% was considered significant. Real-time polymerase chain reaction (PCR) was used to validate differentially expressed genes. To assess barrier function, trans-epithelial electrical resistance (TEER) and dextran flux assays were performed on Caco-2 intestinal epithelial cells that were transfected with miRNA-inhibitors or control inhibitors. Protein expression of barrier function associated genes was confirmed using western blots.

RESULTS

Four out of 247 miRNAs tested were differentially expressed in IBS compared to HCs (FDR < 10%). Real-time PCR validation suggested decreased levels of miR-219a-5p and miR-338-3p in IBS (P = .026 and P = .004), and IBS-C (P = .02 and P = .06) vs. HCs as the strongest associations. Inhibition of miR-219a-5p resulted in altered expression of proteasome/barrier function genes. Functionally, miR-219a-5p inhibition enhanced the permeability of intestinal epithelial cells as TEER was reduced (25-50%, P < .05) and dextran flux was increased (P < .01). Additionally, inhibition of miR-338-3p in cells caused alterations in the mitogen-activated protein kinase (MAPK) signaling pathway genes.

CONCLUSION

Two microRNAs that potentially affect permeability and visceral nociception were identified to be altered in IBS patients. MiR-219a-5p and miR-338-3p potentially alter barrier function and visceral hypersensitivity via neuronal and MAPK signaling and could be therapeutic targets in IBS.

An approach for normalization and quality control for NanoString RNA expression data

The NanoString RNA counting assay for formalin-fixed paraffin embedded samples is unique in its sensitivity, technical reproducibility and robustness for analysis of clinical and archival samples. While commercial normalization methods are provided by NanoString, they are not optimal for all settings, particularly when samples exhibit strong technical or biological variation or where housekeeping genes have variable performance across the cohort. Here, we develop and evaluate a more comprehensive normalization procedure for NanoString data with steps for quality control, selection of housekeeping targets, normalization and iterative data visualization and biological validation. The approach was evaluated using a large cohort ($N=\kern0.5em 1649$) from the Carolina Breast Cancer Study, two cohorts of moderate sample size ($N=359$ and$130$) and a small published dataset ($N=12$). The iterative process developed here eliminates technical variation (e.g. from different study phases or sites) more reliably than the three other methods, including NanoString's commercial package, without diminishing biological variation, especially in long-term longitudinal multiphase or multisite cohorts. We also find that probe sets validated for nCounter, such as the PAM50 gene signature, are impervious to batch issues. This work emphasizes that systematic quality control, normalization and visualization of NanoString nCounter data are an imperative component of study design that influences results in downstream analyses.

The inhibitory receptor TIM-3 limits activation of the cGAS-STING pathway in intra-tumoral dendritic cells by suppressing extracellular DNA uptake

Blockade of the inhibitory receptor TIM-3 shows efficacy in cancer immunotherapy clinical trials. TIM-3 inhibits production of the chemokine CXCL9 by XCR1⁺ classical dendritic cells (cDC1), thereby limiting antitumor immunity in mammary carcinomas. We found that increased CXCL9 expression by splenic cDC1s upon TIM-3 blockade required type I interferons and extracellular DNA. Chemokine expression as well as combinatorial efficacy of TIM-3 blockade and paclitaxel chemotherapy were impaired by deletion of Cgas and Sting. TIM-3 blockade increased uptake of extracellular DNA by cDC1 through an endocytic process that resulted in cytoplasmic localization. DNA uptake and efficacy of TIM-3 blockade required DNA binding by HMGB1, while galectin-9-induced cell surface clustering of TIM-3 was necessary for its suppressive function. Human peripheral blood cDC1s also took up extracellular DNA upon TIM-3 blockade. Thus, TIM-3 regulates endocytosis of extracellular DNA and activation of the cytoplasmic DNA sensing cGAS-STING pathway in cDC1s, with implications for understanding the mechanisms underlying TIM-3 immunotherapy.

Peripheral metabolic state and immune system in first-episode psychosis - A gene expression study with a prospective one-year follow-up

he excess availability of glucose and lipids can also have an impact on the dynamics of activation and regulation of peripheral immune cellsWe aimed at understanding the correlations between peripheral metabolic state and immune system during the first year in first-episode psychosis (FEP). Patients with FEP (n = 67) and matched controls (n = 38), aged 18-40 years, were met at baseline, 2 and 12 months. Fasting peripheral blood samples were collected. We applied the NanoString nCounter in-solution hybridization technology to determine gene expression levels of 178 candidate genes reflecting activation of the immune system. Serum triglycerides, high-density lipoprotein (HDL), low-density lipoprotein (LDL) cholesterol and insulin and plasma glucose (fP-Gluc) were measured. We applied Ingenuity Pathway Analysis (IPA) to visualize enrichment of genes to functional classes. Strength of positive or negative regulation of the disease and functional pathways was deduced from IPA activation Z-score at the three evaluation points. We correlated gene expression with plasma glucose, triglycerids and HDL and LDL, and used hierarchical clustering of the pairwise correlations to identify groups of genes with similar correlation patterns with metabolic markers. In patients, initially, genes associated with the innate immune system response pathways were upregulated, which decreased by 12 months. Furthermore, genes associated with apoptosis and T cell death were downregulated, and genes associated with lipid metabolism were increasingly downregulated by 12 months. The immune activation was thus an acute phase during illness onset. At baseline, after controlling for multiple testing, 31/178 genes correlated positively with fasting glucose levels, and 54/178 genes negatively with triglycerides in patients only. The gene clusters showed patterns of correlations with metabolic markers over time. The results suggest a functional link between peripheral immune system and metabolic state in FEP. Metabolic factors may have had an influence on the initial activation of the innate immune system. Future work is necessary to understand the role of metabolic state in the regulation of immune response in the early phases of psychosis.

Select hyperactivating NLRP3 ligands enhance the TH1- and TH17-inducing potential of human type 2 conventional dendritic cells

The detection of microorganisms and danger signals by pattern recognition receptors on dendritic cells (DCs) and the consequent formation of inflammasomes are pivotal for initiating protective immune responses. Although the activation of inflammasomes leading to secretion of the cytokine IL-1β is typically accompanied by pyroptosis (an inflammatory form of lytic programmed cell death), some cells can survive and exist in a state of hyperactivation. Here, we found that the conventional type 2 DC (cDC2) subset is the major human DC subset that is transcriptionally and functionally poised for inflammasome formation and response without pyroptosis. When cDC2 were stimulated with ligands that relatively weakly activated the inflammasome, the cells did not enter pyroptosis but instead secreted IL-12 family cytokines and IL-1β. These cytokines induced prominent T helper type 1 (T_H1) and T_H17 responses that were superior to those seen in response to Toll-like receptor (TLR) stimulation alone or to stronger, classical inflammasome ligands. These findings not only define the human cDC2 subpopulation as a prime target for the treatment of inflammasome-dependent inflammatory diseases but may also inform new approaches for adjuvant and vaccine development.

Transcriptomic analysis to identify genes associated with selective hippocampal vulnerability in Alzheimer's disease

Selective vulnerability of different brain regions is seen in many neurodegenerative disorders. The hippocampus and cortex are selectively vulnerable in Alzheimer's disease (AD), however the degree of involvement of the different brain regions differs among patients. We classified corticolimbic patterns of neurofibrillary tangles in postmortem tissue to capture extreme and representative phenotypes. We combined bulk RNA sequencing with digital pathology to examine hippocampal vulnerability in AD. We identified hippocampal gene expression changes associated with hippocampal vulnerability and used machine learning to identify genes that were associated with AD neuropathology, including SERPINA5, RYBP, SLC38A2, FEM1B, and PYDC1. Further histologic and biochemical analyses suggested SERPINA5 expression is associated with tau expression in the brain. Our study highlights the importance of embracing heterogeneity of the human brain in disease to identify disease-relevant gene expression.

Multimodal analysis for human ex vivo studies shows extensive molecular changes from delays in blood processing

Multi-omic profiling of human peripheral blood is increasingly utilized to identify biomarkers and pathophysiologic mechanisms of disease. The importance of these platforms in clinical and translational studies led us to investigate the impact of delayed blood processing on the numbers and state of peripheral blood mononuclear cells (PBMC) and on the plasma proteome. Similar to previous studies, we show minimal effects of delayed processing on the numbers and general phenotype of PBMC up to 18 hours. In contrast, profound changes in the single-cell transcriptome and composition of the plasma proteome become evident as early as 6 hours after blood draw. These reflect patterns of cellular activation across diverse cell types that lead to progressive distancing of the gene expression state and plasma proteome from native in vivo biology. Differences accumulating during an overnight rest (18 hours) could confound relevant biologic variance related to many underlying disease states.

•

Studies of human blood cells and plasma are highly sensitive to process variability

•

Time variability distorts biology in cutting-edge single-cell and multiplex assays

•

Longitudinal, multi-modal, and aligned data enable data qualification and exploration

•

Dataset holds potential novel, multi-modal biological correlations and hypotheses

Photobiomodulation Therapy on Myocardial Infarction in Rats: Transcriptional and Posttranscriptional Implications to Cardiac Remodeling

BACKGROUND AND OBJECTIVES

Induction of myocardial infarction (MI) in rats by occlusion of the left anterior descending coronary artery is an experimental model used in research to elucidate functional, structural, and molecular modifications associated with ischemic heart disease. Photobiomodulation therapy (PBMT) has become a therapeutic alternative by modulating various biological processes eliciting several effects, including anti-inflammatory and pro-proliferative actions. The main objective of this work was to evaluate the effect of PBMT in the modulation of transcriptional and post-transcriptional changes that occurred in myocardium signal transduction pathways after MI.

STUDY DESIGN/MATERIALS AND METHODS

Continuous wave (CW) non-thermal laser parameters were: 660 nm wavelength, power 15 mW, with a total energy of 0.9 J, fluence of 1.15 J/cm² , spot size of 0.785 cm² , and time of 60 seconds. Using in silico analysis, we selected and then, quantified the expression of messenger RNA (mRNA) of 47 genes of 9 signaling pathways associated with MI (angiogenesis, cell survival, hypertrophy, oxidative stress, apoptosis, extracellular matrix, calcium kinetics, cell metabolism, and inflammation). Messenger RNA expression quantification was performed in myocardial samples by polymerase chain reaction real-time array using TaqMan customized plates.

RESULTS

Our results evidenced that MI modified mRNA expression of several well-known biomarkers related to detrimental cardiac activity in almost all signaling pathways analyzed. However, PBMT reverted most of these transcriptional changes. More expressively, PBMT provoked a robust decrease in mRNA expression of molecules that participate in post-MI inflammation and ECM composition, such as IL-6, TNF receptor, TGFb1, and collagen I and III. Global microRNA (miRNA) expression analysis revealed that PBMT decreased miR-221, miR-34c, and miR-93 expressions post-MI, which are related to deleterious effects in cardiac remodeling.

CONCLUSION

Thus, the identification of transcriptional and post-transcriptional changes induced by PBMT may be used to interfere in the molecular dynamics of cardiac remodeling post-MI.

Digital Gene Expression Analysis of Epithelioid and Sarcomatoid Mesothelioma Reveals Differences in Immunogenicity

Malignant pleural mesothelioma (MPM) is an aggressive malignancy associated with asbestos exposure. Median survival ranges from 14 to 20 months after initial diagnosis. As of November 2020, the FDA approved a combination of immune checkpoint inhibitors after promising intermediate results. Nonetheless, responses remain unsatisfying. Adequate patient stratification to improve response rates is still lacking. This retrospective study analyzed formalin fixed paraffin embedded specimens from a cohort of 22 MPM. Twelve of those samples showed sarcomatoid, ten epithelioid differentiation. Complete follow-up, including radiological assessment of response by modRECIST and time to death, was available with reported deaths of all patients. RNA of all samples was isolated and subjected to digital gene expression pattern analysis. Our study revealed a notable difference between epithelioid and sarcomatoid mesothelioma, showing differential gene expression for 304/698 expressed genes. Whereas antigen processing and presentation to resident cytotoxic T cells as well as phagocytosis is highly affected in sarcomatoid mesothelioma, cell-cell interaction via cytokines seems to be of greater importance in epithelioid cases. Our work reveals the specific role of the immune system within the different histologic subtypes of MPM, providing a more detailed background of their immunogenic potential. This is of great interest regarding therapeutic strategies including immunotherapy in mesothelioma.

Micro-RNA: The darkhorse of cancer

The discovery of micro RNAs (miRNA) in cancer has opened up new vistas for researchers in recent years. Micro RNAs area set of small, endogenous, highly conserved, non-coding RNAs that control the expression of about 30% genes at post-transcriptional levels. Typically, microRNAs impede the translation and stability of messenger RNAs (mRNA), control genes associated with cellular processes namely inflammation, cell cycle regulation, stress response, differentiation, apoptosis, and migration. Compelling findings revealed that miRNA mutations or disruption correspond to diverse human cancers and suggest that miRNAs can function as tumor suppressors or oncogenes. Here we summarize the literature on these master regulators in clinical settings from last three decades as both abrupt cancer therapeutics and as an approach to sensitize tumors to chemotherapy. This review highlights (I) the prevailing perception of miRNA genomics, biogenesis, as well as function; (II) the significant advancements in regulatory mechanisms in the expression of carcinogenic genes; and (III) explains, how miRNA is utilized as a diagnostic and prognostic biomarker for the disease stage indicating survival as well as therapeutic targets in cancer.

Extracellular vesicle microRNA in early versus late pregnancy with birth outcomes in the MADRES study

Circulating miRNA may contribute to the development of adverse birth outcomes. However, few studies have investigated extracellular vesicle (EV) miRNA, which play important roles in intercellular communication, or compared miRNA at multiple time points in pregnancy. In the current study, 800 miRNA were profiled for EVs from maternal plasma collected in early (median: 12.5 weeks) and late (median: 31.8 weeks) pregnancy from 156 participants in the MADRES Study, a health disparity pregnancy cohort. Associations between miRNA and birth weight, birth weight for gestational age (GA), and GA at birth were examined using covariate-adjusted robust linear regression. Differences by infant sex and maternal BMI were also investigated. Late pregnancy measures of 13 miRNA were associated with GA at birth (P_FDR<0.050). Negative associations were observed for eight miRNA (miR-4454+ miR-7975, miR-4516, let-7b-5p, miR-126-3p, miR-29b-3p, miR-15a-5p, miR-15b-5p, miR-19b-3p) and positive associations for five miRNA (miR-212-3p, miR-584-5p, miR-608, miR-210-3p, miR-188-5p). Predicted target genes were enriched (P_FDR<0.050) in pathways involved in organogenesis and placental development. An additional miRNA (miR-107), measured in late pregnancy, was positively associated with GA at birth in infants born to obese women (P_FDR for BMI interaction = 0.011). In primary analyses, the associations between early pregnancy miRNA and birth outcomes were not statistically significant (P_FDR≥0.05). However, sex-specific associations were observed for early pregnancy measures of 37 miRNA and GA at birth (P_FDR for interactions<0.050). None of the miRNA were associated with fetal growth measures (P_FDR≥0.050). Our findings suggest that EV miRNA in both early and late pregnancy may influence gestational duration.

Impact of checkpoint blockade on cancer vaccine-activated CD8+ T cell responses

Immune and molecular profiling of CD8 T cells of patients receiving DC vaccines expressing three full-length melanoma antigens (MAs) was performed. Antigen expression levels in DCs had no significant impact on T cell or clinical responses. Patients who received checkpoint blockade before DC vaccination had higher baseline MA-specific CD8 T cell responses but no evidence for improved functional responses to the vaccine. Patients who showed the best clinical responses had low PD-1 expression on MA-specific T cells before and after DC vaccination; however, blockade of PD-1 during antigen presentation by DC had minimal functional impact on PD-1high MA-specific T cells. Gene and protein expression analyses in lymphocytes and tumor samples identified critical immunoregulatory pathways, including CTLA-4 and PD-1. High immune checkpoint gene expression networks correlated with inferior clinical outcomes. Soluble serum PD-L2 showed suggestive positive association with improved outcome. These findings show that checkpoint molecular pathways are critical for vaccine outcomes and suggest specific sequencing of vaccine combinations.

3D Functional Genomics Screens Identify CREBBP as a Targetable Driver in Aggressive Triple-Negative Breast Cancer

Triple-negative breast cancers (TNBC) are resistant to standard-of-care chemotherapy and lack known targetable driver gene alterations. Identification of novel drivers could aid the discovery of new treatment strategies for this hard-to-treat patient population, yet studies using high-throughput and accurate models to define the functions of driver genes in TNBC to date have been limited. Here, we employed unbiased functional genomics screening of the 200 most frequently mutated genes in breast cancer, using spheroid cultures to model in vivo-like conditions, and identified the histone acetyltransferase CREBBP as a novel tumor suppressor in TNBC. CREBBP protein expression in patient tumor samples was absent in 8% of TNBCs and at a high frequency in other tumors, including squamous lung cancer, where CREBBP-inactivating mutations are common. In TNBC, CREBBP alterations were associated with higher genomic heterogeneity and poorer patient survival and resulted in upregulation and dependency on a FOXM1 proliferative program. Targeting FOXM1-driven proliferation indirectly with clinical CDK4/6 inhibitors (CDK4/6i) selectively impaired growth in spheroids, cell line xenografts, and patient-derived models from multiple tumor types with CREBBP mutations or loss of protein expression. In conclusion, we have identified CREBBP as a novel driver in aggressive TNBC and identified an associated genetic vulnerability in tumor cells with alterations in CREBBP and provide a preclinical rationale for assessing CREBBP alterations as a biomarker of CDK4/6i response in a new patient population. SIGNIFICANCE: This study demonstrates that CREBBP genomic alterations drive aggressive TNBC, lung cancer, and lymphomas and may be selectively treated with clinical CDK4/6 inhibitors.

Development and validation for research assessment of Oncotype DX® Breast Recurrence Score, EndoPredict® and Prosigna®

Multi-gene prognostic signatures including the Oncotype® DX Recurrence Score (RS), EndoPredict® (EP) and Prosigna® (Risk Of Recurrence, ROR) are widely used to predict the likelihood of distant recurrence in patients with oestrogen-receptor-positive (ER+), HER2-negative breast cancer. Here, we describe the development and validation of methods to recapitulate RS, EP and ROR scores from NanoString expression data. RNA was available from 107 tumours from postmenopausal women with early-stage, ER+, HER2- breast cancer from the translational Arimidex, Tamoxifen, Alone or in Combination study (TransATAC) where previously these signatures had been assessed with commercial methodology. Gene expression was measured using NanoString nCounter. For RS and EP, conversion factors to adjust for cross-platform variation were estimated using linear regression. For ROR, the steps to perform subgroup-specific normalisation of the gene expression data and calibration factors to calculate the 46-gene ROR score were assessed and verified. Training with bootstrapping (n = 59) was followed by validation (n = 48) using adjusted, research use only (RUO) NanoString-based algorithms. In the validation set, there was excellent concordance between the RUO scores and their commercial counterparts (rc(RS) = 0.96, 95% CI 0.93-0.97 with level of agreement (LoA) of -7.69 to 8.12; rc(EP) = 0.97, 95% CI 0.96-0.98 with LoA of -0.64 to 1.26 and rc(ROR) = 0.97 (95% CI 0.94-0.98) with LoA of -8.65 to 10.54). There was also a strong agreement in risk stratification: (RS: κ = 0.86, p < 0.0001; EP: κ = 0.87, p < 0.0001; ROR: κ = 0.92, p < 0.001). In conclusion, the calibrated algorithms recapitulate the commercial RS and EP scores on individual biopsies and ROR scores on samples based on subgroup-centreing method using NanoString expression data.

Analysis of extracellular vesicle mRNA derived from plasma using the nCounter platform

Extracellular vesicles (EVs) are double-layered phospholipid membrane vesicles that are released by most cells and can mediate intercellular communication through their RNA cargo. In this study, we tested if the NanoString nCounter platform can be used for the analysis of EV-mRNA. We developed and optimized a methodology for EV enrichment, EV-RNA extraction and nCounter analysis. Then, we demonstrated the validity of our workflow by analyzing EV-RNA profiles from the plasma of 19 cancer patients and 10 controls and developing a gene signature to differentiate cancer versus control samples. TRI reagent outperformed automated RNA extraction and, although lower plasma input is feasible, 500 μL provided highest total counts and number of transcripts detected. A 10-cycle pre-amplification followed by DNase treatment yielded reproducible mRNA target detection. However, appropriate probe design to prevent genomic DNA binding is preferred. A gene signature, created using a bioinformatic algorithm, was able to distinguish between control and cancer EV-mRNA profiles with an area under the ROC curve of 0.99. Hence, the nCounter platform can be used to detect mRNA targets and develop gene signatures from plasma-derived EVs.

Phase 1 Expansion Cohort of Ramucirumab Plus Pembrolizumab in Advanced Treatment-Naive NSCLC

INTRODUCTION

Data of first-line ramucirumab plus pembrolizumab treatment of programmed death-ligand 1 (PD-L1)-positive NSCLC (cohort E) are reported (NCT02443324).

METHODS

In this multicenter, open-label phase 1a/b trial, patients received ramucirumab 10 mg/kg and pembrolizumab 200 mg every 21 days for up to 35 cycles. PD-L1 positivity was defined as tumor proportion score (TPS) greater than or equal to 1%. Exploratory NanoString biomarker analyses included three T-cell signatures (T-cell-inflamed, Gajewski, and effector T cells) and CD274 gene expression.

RESULTS

Cohort E included 26 patients. Treatment-related adverse events of any grade occurred in 22 patients (84.6%). Treatment-related adverse events of grade greater than or equal to 3 were reported in 11 patients (42.3%); the most frequent was hypertension (n = 4, 15.4%). Objective response rate was 42.3% in the treated population and 56.3% and 22.2% for patients with high (TPS ≥ 50%) and lower levels (TPS 1%-49%) of PD-L1 expression, respectively. Median progression-free survival (PFS) in the treated population was 9.3 months, and 12-month and 18-month PFS rates were 45% each. Median PFS was not reached in patients with PD-L1 TPS greater than or equal to 50% and was 4.2 months in patients with PD-L1 TPS 1% to 49%. Median overall survival was not reached in the treated population, and 12-month and 18-month overall survival rates were 73% and 64%, respectively. Biomarker data suggested a positive association among clinical response, three T-cell signatures, CD274 gene expression, and PD-L1 immunohistochemistry.

CONCLUSIONS

First-line therapy with ramucirumab plus pembrolizumab has a manageable safety profile in patients with NSCLC, and the efficacy signal seems to be strongest in tumors with high PD-L1 expression.

Nutlin-3a Enhances Natural Killer Cell-Mediated Killing of Neuroblastoma by Restoring p53-Dependent Expression of Ligands for NKG2D and DNAM-1 Receptors

In this study, we explored whether Nutlin-3a, a well-known, nontoxic small-molecule compound antagonizing the inhibitory interaction of MDM2 with the tumor suppressor p53, may restore ligands for natural killer (NK) cell-activating receptors (NK-AR) on neuroblastoma cells to enhance the NK cell-mediated killing. Neuroblastoma cell lines were treated with Nutlin-3a, and the expression of ligands for NKG2D and DNAM-1 NK-ARs and the neuroblastoma susceptibility to NK cells were evaluated. Adoptive transfer of human NK cells in a xenograft neuroblastoma-bearing NSG murine model was assessed. Two data sets of neuroblastoma patients were explored to correlate p53 expression with ligand expression. Luciferase assays and chromatin immunoprecipitation analysis of p53 functional binding on PVR promoter were performed. Primary neuroblastoma cells were also treated with Nutlin-3a, and neuroblastoma spheroids obtained from one high-risk patient were assayed for NK-cell cytotoxicity. We provide evidence showing that the Nutlin-3a-dependent rescue of p53 function in neuroblastoma cells resulted in (i) increased surface expression of ligands for NK-ARs, thus rendering neuroblastoma cell lines significantly more susceptible to NK cell-mediated killing; (ii) shrinkage of human neuroblastoma tumor masses that correlated with overall survival upon adoptive transfer of NK cells in neuroblastoma-bearing mice; (iii) and increased expression of ligands in primary neuroblastoma cells and boosting of NK cell-mediated disaggregation of neuroblastoma spheroids. We also found that p53 was a direct transcription factor regulating the expression of PVR ligand recognized by DNAM-1. Our findings demonstrated an immunomodulatory role of Nutlin-3a, which might be prospectively used for a novel NK cell-based immunotherapy for neuroblastoma.

Low MGMT digital expression is associated with a better outcome of IDH1 wildtype glioblastomas treated with temozolomide

INTRODUCTION

Glioblastoma (GBM) is the deadliest primary brain tumor. The standard treatment consists of surgery, radiotherapy, and temozolomide (TMZ). TMZ response is heterogeneous, and MGMT promoter (MGMTp) methylation has been the major predictive biomarker. We aimed to describe the clinical and molecular data of GBMs treated with TMZ, compare MGMT methylation with MGMT expression, and further associate with patient's outcome.

METHODS

We evaluate 112 FFPE adult GBM cases. IDH1 and ATRX expression was analyzed by immunohistochemistry, hotspot TERT promoter (TERTp) mutations were evaluated by Sanger or pyrosequencing, and MGMTp methylation was assessed by pyrosequencing and MGMT mRNA expression using the nCounter® Vantage 3D™ DNA damage and repair panel.

RESULTS

Of the 112 GBMs, 96 were IDH1^WT, and 16 were IDH1^MUT. Positive ATRX expression was found in 91.6% (88/96) of IDH^WT and 43.7% (7/16) of IDH^MUT. TERTp mutations were detected in 70.4% (50/71) of IDH^WT. MGMTp methylation was found in 55.5% (35/63) of IDH^WT and 84.6% (11/13) of IDH^MUT, and as expected, MGMTp methylation was significantly associated with a better response to TMZ. MGMT expression was inversely correlated with MGMTp methylation levels (- 0.506, p < 0.0001), and MGMT low expression were significantly associated with better patient survival. It was also observed that integrating MGMTp methylation and expression, significantly improved the prognostication value.

CONCLUSIONS

MGMT mRNA levels evaluated by digital expression were associated with the outcome of TMZ-treated GBM patients. The combination of MGMT methylation and mRNA expression may provide a more accurate prediction of TMZ response in GBM patients.

Changes in the brain transcriptome after DNA Aβ42 trimer immunization in a 3xTg-AD mouse model

Alzheimer's disease (AD) has been associated with accumulation of amyloid beta (Aβ) peptides in brain, and immunotherapy targeting Aβ provides potential for AD prevention. We have used a DNA Aβ42 trimer construct for immunization of 3xTg-AD mice and found previously significant reduction of amyloid and tau pathology due to the immunotherapy. We show here that DNA Aβ42 immunized 3xTg-AD mice showed better performance in nest building activities and had a higher 24 months survival rate compared to the non-treated AD controls. The analysis of differently expressed genes in brains from 24 months old mice showed significant increases transcript levels between non-immunized AD mice and wild-type controls for genes involved in microglia and astrocyte function, cytokine and inflammatory signaling, apoptosis, the innate and adaptive immune response and are consistent with an inflammatory phenotype in AD. Most of these upregulated genes were downregulated in the DNA Aβ42 immunized 3xTg-AD mice due to the vaccine. Transcript numbers for the immediate early genes, Arc, Bdnf, Homer1, Egr1 and cfos, involved in neuronal and neurotransmission pathways which were much lower in the non-immunized 3xTg-AD mice, were restored to wild-type mouse brain levels in DNA Aβ42 immunized 3xTg-AD mice indicating positive effects of DNA Aβ42 immunotherapy on synapse stability and plasticity. The immune response after immunization is complex, but the multitude of changes after DNA Aβ42 immunization shows that this response moves beyond the amyloid hypothesis and into direction of disease prevention.

Whole blood mRNA expression-based targets to discriminate active tuberculosis from latent infection and other pulmonary diseases

Current diagnostic tests for tuberculosis (TB) are not able to predict reactivation disease progression from latent TB infection (LTBI). The main barrier to predicting reactivation disease is the lack of our understanding of host biomarkers associated with progression from latent infection to active disease. Here, we applied an immune-based gene expression profile by NanoString platform to identify whole blood markers that can distinguish active TB from other lung diseases (OPD), and that could be further evaluated as a reactivation TB predictor. Among 23 candidate genes that differentiated patients with active TB from those with OPD, nine genes (CD274, CEACAM1, CR1, FCGR1A/B, IFITM1, IRAK3, LILRA6, MAPK14, PDCD1LG2) demonstrated sensitivity and specificity of 100%. Seven genes (C1QB, C2, CCR2, CCRL2, LILRB4, MAPK14, MSR1) distinguished TB from LTBI with sensitivity and specificity between 82 and 100%. This study identified single gene candidates that distinguished TB from OPD and LTBI with high sensitivity and specificity (both > 82%), which may be further evaluated as diagnostic for disease and as predictive markers for reactivation TB.

Cell-free DNA promotes malignant transformation in non-tumor cells

Cell-free DNA is present in different biological fluids and when released by tumor cells may contribute to pro-tumor events such as malignant transformation of cells adjacent to the tumor and metastasis. Thus, this study analyzed the effect of tumor cell-free DNA, isolated from the blood of prostate cancer patients, on non-tumor prostate cell lines (RWPE-1 and PNT-2). To achieve this, we performed cell-free DNA quantification and characterization assays, evaluation of gene and miRNA expression profiling focused on cancer progression and EMT, and metabolomics by mass spectrometry and cellular migration. The results showed that tumor-free cell DNA was able to alter the gene expression of MMP9 and CD44, alter the expression profile of nine miRNAs, and increased the tryptophan consumption and cell migration rates in non-tumor cells. Therefore, tumor cell-free DNA was capable of altering the receptor cell phenotype, triggering events related to malignant transformation in these cells, and can thus be considered a potential target for cancer diagnosis and therapy.

Cellular and gene signatures of tumor-infiltrating dendritic cells and natural-killer cells predict prognosis of neuroblastoma

Tumor-infiltrating lymphocytes play an essential role in improving clinical outcome of neuroblastoma (NB) patients, but their relationship with other tumor-infiltrating immune cells in the T cell-inflamed tumors remains poorly investigated. Here we show that dendritic cells (DCs) and natural killer (NK) cells are positively correlated with T-cell infiltration in human NB, both at transcriptional and protein levels, and associate with a favorable prognosis. Multiplex imaging displays DC/NK/T cell conjugates in the tumor microenvironment of low-risk NB. Remarkably, this connection is further strengthened by the identification of gene signatures related to DCs and NK cells able to predict survival of NB patients and strongly correlate with the expression of PD-1 and PD-L1. In summary, our findings unveil a key prognostic role of DCs and NK cells and indicate their related gene signatures as promising tools for the identification of clinical biomarkers to better define risk stratification and survival of NB patients.

Tumour-infiltrating lymphocytes play a crucial role in neuroblastoma, but their relationship to other immune cells is poorly understood. Here the authors identify the cellular and gene signatures of intratumoural dendritic cells and natural killer cells that predict the clinical outcome of neuroblastoma.

Prediction of prognostic signatures in triple-negative breast cancer based on the differential expression analysis via NanoString nCounter immune panel

BACKGROUND

Triple-Negative Breast Cancer (TNBC) is an aggressive and complex subtype of breast cancer. The current biomarkers used in the context of breast cancer treatment are highly dependent on the targeting of oestrogen receptor, progesterone receptor, or HER2, resulting in treatment failure and disease recurrence and creating clinical challenges. Thus, there is still a crucial need for the improvement of TNBC treatment; the discovery of effective biomarkers that can be easily translated to the clinics is essential.

METHODS

We report an approach for the discovery of biomarkers that can predict tumour relapse and pathologic complete response (pCR) in TNBC on the basis of mRNA expression quantified using the NanoString nCounter Immunology Panel. To overcome the limited sample size, prediction models based on random Forest were constructed using the differentially expressed genes (DEGs) as selected features. We also evaluated the differences between pre- and post-treatment groups aiming for the combinatorial assessment of pCR and relapse using additive models in edgeR.

RESULTS

We identify nine and 13 DEGs strongly associated with pCR and relapse, respectively, from 579 immune genes in a small number of samples (n = 55) using edgeR. An additive model for the comparison of pre- and post-treatment groups via the adjustment of the independent subject in the relapse group revealed associations for 41 genes. Comprehensive analysis indicated that our prediction models outperformed those constructed using features extracted from the existing feature selection model Elastic Net in terms of accuracy. The prediction models were assessed using a randomization test to validate the robustness (empirical P for the model of pCR = 0.015 and empirical P for the model of relapse = 0.018). Furthermore, three DEGs (FCER1A, EDNRB, and TGFBI) in the model of relapse showed prognostic significance for predicting the survival of patients with cancer through Cox proportional hazards regression model-based survival analysis.

CONCLUSION

Gene expression quantified via the NanoString nCounter Immunology Panel can be seamlessly analysed using edgeR, even considering small sample sizes. Our approach provides a scalable framework that can easily be applied for the discovery of biomarkers based on the NanoString nCounter Immunology Panel.

DATA AVAILABILITY

The source code will be available from github at https://github.com/sungheep/nanostring .

Ramucirumab in Combination with Pembrolizumab in Treatment-Naïve Advanced Gastric or GEJ Adenocarcinoma: Safety and Antitumor Activity from the Phase 1a/b JVDF Trial

Simple Summary

The prognosis for gastric cancer remains poor, with a median overall survival of approximately 1 year. Ramucirumab and pembrolizumab have each demonstrated antitumor activity and a favorable safety profile as treatments for patients with advanced gastric/gastroesophageal junction (G/GEJ) cancer in the second and third-line setting respectively. However, both agents failed to demonstrate survival benefit over chemotherapy in the first-line setting. Twenty-eight treatment-naïve patients with advanced/metastatic G/GEJ adenocarcinoma were treated with ramucirumab plus pembrolizumab in this phase 1a/b trial. Our results showed that this combination was well tolerated with no unexpected toxicities, and promising durable survival results, particularly among patients with PD-ligand 1 positive tumors. The results of our study therefore support modulating the tumor microenvironment with dual inhibition of VEGFR2 and PD-1 pathways in the first-line patients with advanced G/GEJ cancer.

Abstract

Ramucirumab (anti-VEGFR2) plus pembrolizumab (anti-PD1) demonstrated promising antitumor activity and tolerability among patients with previously treated advanced cancers, supporting growing evidence that combination therapies modulating the tumor microenvironment may expand the spectrum of patients who respond to checkpoint inhibitors. Here we present the results of this combination in first-line patients with metastatic G/GEJ cancer. Twenty-eight patients (≥18 years) with no prior systemic chemotherapy in the advanced/metastatic setting received ramucirumab (8 mg/kg days 1 and 8) plus pembrolizumab (200 mg day 1) every 3 weeks as part of JVDF phase 1a/b study. The primary endpoint was safety. Secondary endpoints included progression-free survival (PFS), objective response rate (ORR), and overall survival (OS). Tumors were PD-L1-positive (combined positive score ≥ 1) in 19 and -negative in 6 patients. Eighteen patients experienced grade 3 treatment-related adverse events, most commonly hypertension (14%) and elevated alanine/aspartate aminotransferase (11% each), with no grade 4 or 5 reported. The ORR was 25% (PD-L1-positive, 32%; PD-L1-negative, 17%) with duration of response not reached. PFS was 5.6 months (PD-L1-positive, 8.6 months; PD-L1-negative, 4.3 months), and OS 14.6 months (PD-L1-positive, 17.3 months; PD-L1-negative, 11.3 months). Acknowledging study design limitations, ramucirumab plus pembrolizumab had encouraging durable clinical activity with no unexpected toxicities in treatment-naïve biomarker-unselected metastatic G/GEJ cancer, and improved outcomes in patients with PD-L1-positive tumors.

Temporal Stability and Prognostic Biomarker Potential of the Prostate Cancer Urine miRNA Transcriptome

BACKGROUND

The development of noninvasive tests for the early detection of aggressive prostate tumors is a major unmet clinical need. miRNAs are promising noninvasive biomarkers: they play essential roles in tumorigenesis, are stable under diverse analytical conditions, and can be detected in body fluids.

METHODS

We measured the longitudinal stability of 673 miRNAs by collecting serial urine samples from 10 patients with localized prostate cancer. We then measured temporally stable miRNAs in an independent training cohort (n = 99) and created a biomarker predictive of Gleason grade using machine-learning techniques. Finally, we validated this biomarker in an independent validation cohort (n = 40).

RESULTS

We found that each individual has a specific urine miRNA fingerprint. These fingerprints are temporally stable and associated with specific biological functions. We identified seven miRNAs that were stable over time within individual patients and integrated them with machine-learning techniques to create a novel biomarker for prostate cancer that overcomes interindividual variability. Our urine biomarker robustly identified high-risk patients and achieved similar accuracy as tissue-based prognostic markers (area under the receiver operating characteristic = 0.72, 95% confidence interval = 0.69 to 0.76 in the training cohort, and area under the receiver operating characteristic curve = 0.74, 95% confidence interval = 0.55 to 0.92 in the validation cohort).

CONCLUSIONS

These data highlight the importance of quantifying intra- and intertumoral heterogeneity in biomarker development. This noninvasive biomarker may usefully supplement invasive or expensive radiologic- and tissue-based assays.

Comparative expression profiling in the intestine of patients with Giardia-induced postinfectious functional gastrointestinal disorders

BACKGROUND

A Giardia outbreak in Bergen, Norway, caused postinfectious functional gastrointestinal disorders (PI-FGIDs). Despite the devastating effects of this outbreak, it presented a unique chance to investigate the implication on the dysregulation of genetic pathways in PI-FGID.

METHODS

We performed the first comparative expression profiling of miRNAs and their potential target genes in microdissected rectal biopsies from 20 Giardia-induced PI-FGID patients vs 18 healthy controls by nCounter analysis. Subsequently, candidates were validated on protein level by immunostaining.

KEY RESULTS

miRNA profiling on rectal biopsy samples from 5 diarrhea-predominant PI-IBS cases compared to 10 healthy controls revealed differential expression in the epithelial layer. The top five regulated miRNAs were implicated in GI disease, inflammatory response, and immunological disease. Subsequently, these miRNAs and 100 potential mRNA targets were examined in 20 PI-FGID cases and 18 healthy controls in both the mucosal epithelium and the lamina propria. Although deregulation of the selected miRNAs could not be verified in the larger sample set, mRNAs involved in barrier function were downregulated in the epithelium. Pro-inflammatory genes and genes implicated in epigenetic modifications were upregulated in the lamina propria. Immunostaining for selected candidates on 17 PI-FGID cases and 16 healthy controls revealed increased tryptase levels as well as a decreased and aberrant subcellular expression of occludin.

CONCLUSIONS AND INFERENCES

Genes relevant to immune and barrier function as well as stress response and epigenetic modulation are differentially expressed in PI-FGIDs and may contribute to disease manifestation.

The combination of coffee compounds attenuates early fibrosis-associated hepatocarcinogenesis in mice: involvement of miRNA profile modulation

Aberrant microRNA expression implicates on hepatocellular carcinoma (HCC) development. Conversely, coffee consumption reduces by ~40% the risk for fibrosis/cirrhosis and HCC, while decaffeinated coffee does not. It is currently unknown whether these protective effects are related to caffeine (CAF), or to its combination with other common and/or highly bioavailable coffee compounds, such as trigonelline (TRI) and chlorogenic acid (CGA). We evaluated whether CAF individually or combined with TRI and/or CGA alleviates fibrosis-associated hepatocarcinogenesis, examining the involvement of miRNA profile modulation. Then, male C3H/HeJ mice were submitted to a diethylnitrosamine/carbon tetrachloride-induced model. Animals received CAF (50 mg/kg), CAF+TRI (50 and 25 mg/kg), CAF+CGA (50 and 25 mg/kg) or CAF+TRI+CGA (50, 25 and 25 mg/kg), intragastrically, 5×/week, for 10 weeks. Only CAF+TRI+CGA combination reduced the incidence, number and proliferation (Ki-67) of hepatocellular preneoplastic foci while enhanced apoptosis (cleaved caspase-3) in adjacent parenchyma. CAF+TRI+CGA treatment also decreased hepatic oxidative stress and enhanced the antioxidant Nrf2 axis. CAF+TRI+CGA had the most pronounced effects on decreasing hepatic pro-inflammatory IL-17 and NFκB, contributing to reduce CD68-positive macrophage number, stellate cell activation, and collagen deposition. In agreement, CAF+TRI+CGA upregulated tumor suppressors miR-144-3p, miR-376a-3p and antifibrotic miR-15b-5p, frequently deregulated in human HCC. CAF+TRI+CGA reduced the hepatic protein levels of pro-proliferative EGFR (miR-144-3p target), antiapoptotic Bcl-2 family members (miR-15b-5p targets), and the number of PCNA (miR-376a-3p target) positive hepatocytes in preneoplastic foci. Our results suggest that the combination of most common and highly bioavailable coffee compounds, rather than CAF individually, attenuates fibrosis-associated hepatocarcinogenesis by modulating miRNA expression profile.

Combining BRAF inhibition with oncolytic herpes simplex virus enhances the immune-mediated antitumor therapy of BRAF-mutant thyroid cancer

BACKGROUND

The aggressive clinical behavior of poorly differentiated and anaplastic thyroid cancers (PDTC and ATC) has proven challenging to treat, and survival beyond a few months from diagnosis is rare. Although 30%-60% of these tumors contain mutations in the BRAF gene, inhibitors designed specifically to target oncogenic BRAF have shown limited and only short-lasting therapeutic benefits as single agents, thus highlighting the need for improved treatment strategies, including novel combinations.

METHODS

Using a BRAFV600E-driven mouse model of ATC, we investigated the therapeutic efficacy of the combination of BRAF inhibition and oncolytic herpes simplex virus (oHSV). Analyses of samples from tumor-bearing mice were performed to immunologically characterize the effects of different treatments. These immune data were used to inform the incorporation of immune checkpoint inhibitors into triple combination therapies.

RESULTS

We characterized the immune landscape in vivo following BRAF inhibitor treatment and detected only modest immune changes. We, therefore, hypothesized that the addition of oncolytic virotherapy to BRAF inhibition in thyroid cancer would create a more favorable tumor immune microenvironment, boost the inflammatory status of tumors and improve BRAF inhibitor therapy. First, we showed that thyroid cancer cells were susceptible to infection with oHSV and that this process was associated with activation of the immune tumor microenvironment in vivo. Next, we showed improved therapeutic responses when combining oHSV and BRAF inhibition in vivo, although no synergistic effects were seen in vitro, further confirming that the dominant effect of oHSV in this context was likely immune-mediated. Importantly, both gene and protein expression data revealed an increase in activation of T cells and natural killer (NK) cells in the tumor in combination-treated samples. The benefit of combination oHSV and BRAF inhibitor therapy was abrogated when T cells or NK cells were depleted in vivo. In addition, we showed upregulation of PD-L1 and CTLA-4 following combined treatment and demonstrated that blockade of the PD-1/PD-L1 axis or CTLA-4 further improved combination therapy.

CONCLUSIONS

The combination of oHSV and BRAF inhibition significantly improved survival in a mouse model of ATC by enhancing immune-mediated antitumor effects, and triple combination therapies, including either PD-1 or CTLA-4 blockade, further improved therapy.

Microarray data analysis on gene and miRNA expression to identify biomarkers in non-small cell lung cancer

Background

The aim of this study was to gain further investigation of non-small cell lung cancer (NSCLC) tumorigenesis and identify biomarkers for clinical management of patients through comprehensive bioinformatics analysis.

Methods

miRNA and mRNA microarray datasets were downloaded from GEO (Gene Expression Omnibus) database under the accession number GSE102286 and GSE101929, respectively. Genes and miRNAs with differential expression were identified in NSCLC samples compared with controls, respectively. The interaction between differentially expressed genes (DEGs) and differentially expressed miRNAs (DEmiRs) was predicted, followed by functional enrichment analysis, and construction of miRNA-gene regulatory network, protein-protein interaction (PPI) network, and competing endogenous RNA (ceRNA) network. Through comprehensive bioinformatics analysis, we anticipate to find novel therapeutic targets and biomarkers for NSCLC.

Results

A total of 123 DEmiRs (5 up- and 118 down-regulated miRNAs) and 924 DEGs (309 up- and 615 down-regulated genes) were identified. These genes and miRNAs were significantly involved in different pathways including adherens junction, relaxin signaling pathway, and axon guidance. Furthermore, hsa-miR-9-5p, has-miR-196a-5p and hsa-miR-31-5p, as well as hsa-miR-1, hsa-miR-218-5p and hsa-miR-135a-5p were shown to have higher degree in the miRNA-gene regulatory network and ceRNA network, respectively. Furthermore, BIRC5 and FGF2, as well as RTKN2 and SLIT3 were hubs in the PPI network and ceRNA network, respectively.

Conclusion

Several pathways (adherens junction, relaxin signaling pathway, and axon guidance) miRNAs (hsa-miR-9-5p, has-miR-196a-5p, hsa-miR-31-5p, hsa-miR-1, hsa-miR-218-5p and hsa-miR-135a-5p) and genes (BIRC5, FGF2, RTKN2 and SLIT3) may play important roles in the pathogenesis of NSCLC.

Immune subtyping of extranodal NK/T-cell lymphoma: a new biomarker and an immune shift during disease progression

Extranodal NK/T-cell lymphoma is an aggressive lymphoma that is strongly associated with Epstein-Barr virus infection. Although some extranodal NK/T-cell lymphoma patients have shown responses to immune checkpoint blockade, biomarkers for predicting extranodal NK/T-cell lymphoma patient response to immunotherapy have not yet been defined. To understand the tumor immune microenvironment, we analyzed the expression of 579 immune-related genes and characterized the immune cells using immunohistochemistries and in situ hybridization for EBER. Based on comprehensive analyses, we developed an immune subtyping model that classifies extranodal NK/T-cell lymphoma patients into four tumor immune microenvironment subgroups using three immunohistochemical markers (FoxP3, PD-L1, and CD68). The four tumor immune microenvironment subgroups were named immune tolerance, immune evasion-A, immune evasion-B, and immune silenced. The immune tolerance group was characterized by high-Treg counts and was frequently observed in early stage, and nasal extranodal NK/T-cell lymphoma. The immune evasion group showed high cytotoxic T-cell counts and high PD-L1 expression but low Treg counts. In the immune-silenced group, almost all immune responses were exhausted, most patients were at an advanced stage, and had the poorest disease prognosis among the tumor immune microenvironment subgroups. In some patients (n = 3), a shift in the tumor immune microenvironment subgroup classification was observed in sequential biopsies. The response rate to pembrolizumab, an anti-PD-1 antibody, was 100% (1/1) in the immune tolerance group, 60% (3/5) in the immune evasion group, and 0% (0/5) in the immune-silenced group. We classified extranodal NK/T-cell lymphoma into four tumor immune microenvironment subgroups using a new classification system. In conclusion, we propose that the tumor immune microenvironment of extranodal NK/T-cell lymphoma may change during disease progression and may serve as a useful biomarker for immunotherapy.

A new normalization for Nanostring nCounter gene expression data

The Nanostring nCounter gene expression assay uses molecular barcodes and single molecule imaging to detect and count hundreds of unique transcripts in a single reaction. These counts need to be normalized to adjust for the amount of sample, variations in assay efficiency and other factors. Most users adopt the normalization approach described in the nSolver analysis software, which involves background correction based on the observed values of negative control probes, a within-sample normalization using the observed values of positive control probes and normalization across samples using reference (housekeeping) genes. Here we present a new normalization method, Removing Unwanted Variation-III (RUV-III), which makes vital use of technical replicates and suitable control genes. We also propose an approach using pseudo-replicates when technical replicates are not available. The effectiveness of RUV-III is illustrated on four different datasets. We also offer suggestions on the design and analysis of studies involving this technology.

miRNA expression profiling of hereditary breast tumors from BRCA1- and BRCA2-germline mutation carriers in Brazil

Background

MicroRNAs (miRNAs) are small non-coding RNAs involved in post-transcriptional gene expression regulation and have been described as key regulators of carcinogenesis. Aberrant miRNA expression has been frequently reported in sporadic breast cancers, but few studies have focused on profiling hereditary breast cancers. In this study, we aimed to identify specific miRNA signatures in hereditary breast tumors and to compare with sporadic breast cancer and normal breast tissues.

Methods

Global miRNA expression profiling using NanoString technology was performed on 43 hereditary breast tumors (15 BRCA1, 14 BRCA2, and 14 BRCAX), 23 sporadic breast tumors and 8 normal breast tissues. These normal breast tissues derived from BRCA1- and BRCA2- mutation carriers (n = 5) and non-mutation carriers (n = 3). Subsequently, we performed receiver operating characteristic (ROC) curve analyses to evaluate the diagnostic performance of differentially expressed miRNAs. Putative target genes of each miRNAs considered as potential biomarkers were identified using miRDIP platform and used for pathway enrichment analysis.

Results

miRNA expression analyses identified several profiles that were specific to hereditary breast cancers. A total of 25 miRNAs were found to be differentially expressed (fold change: > 2.0 and p < 0.05) and considered as potential biomarkers (area under the curve > 0.75) in hereditary breast tumors compared to normal breast tissues, with an expressive upregulation among BRCAX cases. Furthermore, bioinformatic analysis revealed that these miRNAs shared target genes involved in ErbB, FoxO, and PI3K-Akt signaling pathways.

Conclusions

Our results showed that miRNA expression profiling can differentiate hereditary from sporadic breast tumors and normal breast tissues. These miRNAs were remarkably deregulated in BRCAX hereditary breast cancers. Therefore, miRNA signatures can be used as potential novel diagnostic biomarkers for the prediction of BRCA1/2- germline mutations and may be useful for future clinical management.

Loss of 5'-Methylthioadenosine Phosphorylase (MTAP) is Frequent in High-Grade Gliomas; Nevertheless, it is Not Associated with Higher Tumor Aggressiveness

The 5’-methylthioadenosine phosphorylase (MTAP) gene is located in the chromosomal region 9p21. MTAP deletion is a frequent event in a wide variety of human cancers; however, its biological role in tumorigenesis remains unclear. The purpose of this study was to characterize the MTAP expression profile in a series of gliomas and to associate it with patients’ clinicopathological features. Moreover, we sought to evaluate, through glioma gene-edited cell lines, the biological impact of MTAP in gliomas. MTAP expression was evaluated in 507 glioma patients by immunohistochemistry (IHC), and the expression levels were associated with patients’ clinicopathological features. Furthermore, an in silico study was undertaken using genomic databases totalizing 350 samples. In glioma cell lines, MTAP was edited, and following MTAP overexpression and knockout (KO), a transcriptome analysis was performed by NanoString Pan-Cancer Pathways panel. Moreover, MTAP’s role in glioma cell proliferation, migration, and invasion was evaluated. Homozygous deletion of 9p21 locus was associated with a reduction of MTAP mRNA expression in the TCGA (The Cancer Genome Atlas) - glioblastoma dataset (p < 0.01). In addition, the loss of MTAP expression was markedly high in high-grade gliomas (46.6% of cases) determined by IHC and Western blotting (40% of evaluated cell lines). Reduced MTAP expression was associated with a better prognostic in the adult glioblastoma dataset (p < 0.001). Nine genes associated with five pathways were differentially expressed in MTAP-knockout (KO) cells, with six upregulated and three downregulated in MTAP. Analysis of cell proliferation, migration, and invasion did not show any significant differences between MTAP gene-edited and control cells. Our results integrating data from patients as well as in silico and in vitro models provide evidence towards the lack of strong biological importance of MTAP in gliomas. Despite the frequent loss of MTAP, it seems not to have a clinical impact in survival and does not act as a canonic tumor suppressor gene in gliomas.

A 4-gene signature from histologically normal surgical margins predicts local recurrence in patients with oral carcinoma: clinical validation

Prognostic biomarkers for recurrence of Oral Squamous Cell Carcinoma (OSCC) are urgently needed. We aimed to independently validate a 4-gene expression signature (MMP1, COL4A1, P4HA2, THBS2) predictive of OSCC recurrence risk. Gene expression was measured using Nanostring nCounter^® in 245 histologically normal surgical resection margins from 62 patients. Association between risk scores for individual patients and recurrence was assessed by Kaplan-Meier analysis. Signature performance was quantified by concordance index (CI), hazard ratio (HR) and the area under receiver operating characteristics (AUC). Risk scores for recurrence were significantly higher than recurrence-free patients (p = 9.58e-7, Welch’s t-test). A solid performance of the 4-gene signature was determined: CI = 0.64, HR = 3.38 (p = 1.4E-4; log-rank test), AUC = 0.71. We showed that three margins per patient are sufficient to preserve predictive performance (CI = 0.65; HR = 2.92; p = 2.94e-3; AUC = 0.71). Association between the predicted risk scores and recurrence was assessed and showed HR = 2.44 (p = 9.6E-3; log-rank test, N = 62). Signature performance analysis was repeated using an optimized threshold (70^th percentile of risks), resulting in HR = 3.38 (p = 1.4E-4; log-rank test, N = 62). The 4-gene signature was validated as predictive of recurrence risk in an independent cohort of patients with resected OSCC and histologically negative margins, and is potentially applicable for clinical decision making on adjuvant treatment and disease monitoring.

Moderate prenatal stress may buffer the impact of Superstorm Sandy on placental genes: Stress in Pregnancy (SIP) Study

The placenta plays a central role in the epigenetic programming of neurodevelopment by prenatal stress (PS), but this pathway is not fully understood. It difficult to study in humans because the conditions for intense, traumatic PS are almost impossible to create ethically. This study was able to capitalize on a 2012 disaster that hit New York, Superstorm Sandy, to examine the impact of traumatic stress on placental gene expression while also examining normative PS, and compare the two. Of the 303 expectant mothers participating in the Stress in Pregnancy Study, 95 women were pregnant when Superstorm Sandy struck. During their pregnancy, participants completed self-report measures of PS and distress that were combined, using latent profile analysis, into one global indicator of normative PS. Placental tissue was collected at delivery and frozen for storage. RNA expression was assessed for 40 placental genes known to associate with the stress response system and neurodevelopment in offspring. Results showed that normative PS increased expression of just MECP2, HSD11B2, and ZNF507, whereas Superstorm Sandy PS decreased expression of CDKL5, CFL1, DYRK1A, HSD11B2, MAOA, MAOB, NCOR1, and ZNF507. Interaction analyses indicated that Superstorm Sandy PS was associated with decreased gene expression for the low and high PS group for CFL1, DYRK1A, HSD11B2, MAOA, and NCOR1 and increased expression for the moderate PS group for FOXP1, NR3C1, and NR3C2. This study supports the idea that a moderate amount of normative PS may buffer the impact of traumatic PS, in this case caused by Superstorm Sandy, on placental gene expression, which suggests that the placenta itself mirrors the organism's ability to develop an epigenetic resilience to, and inoculation from, stress.

Real-time ex vivo perfusion of human lymph nodes invaded by cancer (REPLICANT): a feasibility study

Understanding how breast cancer (BC) grows in axillary lymph nodes (ALNs), and refining how therapies might halt that process, is clinically important. However, modelling the complex ALN microenvironment is difficult, and no human models exist at present. We harvested ALNs from ten BC patients, and perfused them at 37 °C ex vivo for up to 24 h. Controlled autologous testing showed that ALNs remain viable after 24 h of ex vivo perfusion: haematoxylin and eosin-stained histological appearance and proliferation (by Ki67 immunohistochemistry) did not change significantly over time for any perfused ALN compared with a control from time-point zero. Furthermore, targeted gene expression analysis (NanoString PanCancer IO360 panel) showed that only 21/750 genes were differentially expressed between control and perfused ALNs (|log₂ FC| > 1 and q < 0.1): none were involved in apoptosis and metabolism, but rather all 21 genes were involved in immune function and angiogenesis. During perfusion, tissue acid-base balance remained stable. Interestingly, the flow rate increased (p < 0.001) in cancer-replaced (i.e. metastasis occupied more than 90% of the surface area on multiple levels) compared to cancer-free nodes (i.e. nodes with no metastasis on multiple sections). CXCL11 transcripts were significantly more abundant in cancer-replaced nodes, while CXCL12 transcripts were significantly more abundant in cancer-free nodes. These cytokines were also detected in the circulating perfusate. Monoclonal antibodies (nivolumab and trastuzumab) were administered into a further three ALNs to confirm perfusion efficacy. These drugs saturated the nodes; nivolumab even induced cancer cell death. Normothermic ALN perfusion is not only feasible but sustains the tumour microenvironment ex vivo for scientific investigation. This model could facilitate the identification of actionable immuno-oncology targets. © 2019 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Early Enrichment of ESR1 Mutations and the Impact on Gene Expression in Presurgical Primary Breast Cancer Treated with Aromatase Inhibitors

PURPOSE

To investigate the presence of ESR1 mutations in primary estrogen-receptor-positive (ER+) breast cancer treated with extended (>4 weeks) neoadjuvant (presurgical) aromatase inhibitor (NAI) therapy and to identify patients who may gain less benefit from aromatase inhibition (AI) alone based upon on-treatment changes in gene expression.

EXPERIMENTAL DESIGN

We evaluated ER, progesterone receptor, and Ki67 by immunostaining, ESR1 mutations by droplet-digital PCR and expression of over 800 key breast cancer genes in paired pre- and post-NAI tumor samples from 87 ER+ breast cancer patients.

RESULTS

Cell proliferation and estrogen-regulated genes (ERG) remained suppressed in most tumors indicative of persistent response to NAI. Enrichment of ESR1 mutations was found in five tumors and predominantly in patients receiving therapy for >6 months. ESR1-mutant tumors showed increased expression of ESR1 transcript and limited suppression of ERGs and proliferation-associated genes in response to NAI. ESR1 wild-type tumors with high residual proliferation (Ki67r ≥10%; 15/87 tumors) showed lower ESR1/ER expression pre- and post-therapy and lower ERGs. Tumors with ESR1 mutations or Ki67r ≥10% showed less inhibition of estrogen response, cell cycle, and E2F-target genes.

CONCLUSIONS

Ligand-independent ER signaling, as a result of ESR1 mutation or reduced ER dependence, identified after extended NAI therapy, can guide early selection of patients who would benefit from combination therapy.

Digital Immune-Related Gene Expression Signatures In High-Grade Serous Ovarian Carcinoma: Developing Prediction Models For Platinum Response

Purpose

Response to platinum-based therapy is a major prognostic factor in high-grade serous ovarian cancer (HGSOC). While the exact mechanisms of platinum-resistance remain unclear, evidence is accumulating for a connection between the organism’s immune-response and response to platinum. However, predictive tools are missing. This study was performed to examine the putative role of the genetic tumor immune-microenvironment in mediating differential chemotherapy response in HGSOC patients.

Patients and methods

Expression profiling of 770 immune-related genes was performed in tumor tissues from 23 HGSOC cases. Tumors were screened for prognostic and predictive biomarkers using the NanoString nCounter platform for digital gene expression analysis with the appurtenant PanCancer Immune Profiling panel. As validation cohort, gene expression data (RNA Seq) of 303 patients with epithelial ovarian carcinoma (EOC) were retrieved from the The Cancer Genome Atlas (TCGA) database. Different scoring-systems were computed for prediction of risk-of-resistance to cisplatin, disease-free survival (DFS) and overall survival (OS).

Results

Validated on the TCGA-dataset, the developed scores identified 11 significantly differentially expressed genes (p <0.01**) associated with platinum response. HSD11B1 was highly significantly associated with lower risk of recurrence and 7 targets were found highly significantly influencing OS time (p <0.01**).

Conclusion

Our results suggest that response to platinum-based therapy and DFS in ovarian HGSOC is associated with distinct gene-expression patterns related to the tumor immune-system. We generated predictive scoring systems which proved valid when applied to a set of 303 EOC patients.

Functional validation of metabolic genes that distinguish Gleason 3 from Gleason 4 prostate cancer foci

BACKGROUND

Gleason grade is among the most powerful clinicopathological classification systems used to assess risk of lethal potential in prostate cancer, yet its biologic basis is poorly understood. Notably, pure low-grade cancers, comprised predominantly of Gleason pattern 3 (G3) are typically indolent, with lethal potential emerging with the progression of higher-grade Gleason patterns 4 (G4) or 5. One of the hallmarks of more aggressive cancer phenotypes is the stereotyped set of metabolic characteristics that transformed cells acquire to facilitate unregulated growth. In the present study, we profiled expression signatures of metabolic genes that are differentially expressed between G3 and G4 cancer foci and investigated the functional role of two of the profiled genes, PGRMC1 and HSD17B4, in prostate cancer cells.

METHODS

Gene expression profiling was conducted using 32 G3 and 32 G4 cancer foci from patients with 3+3 and ≥4+3 tumors, respectively. A 95-gene Nanostring probe set was used to probe genes associated with energy metabolism. Two out of five genes (PGRMC1 and HSD17B4) that significantly distinguish between G3 and G4 were functionally validated in vitro using established prostate cancer cells (PC3, DU145). Expression of PGRMC1 and HSD17B4 was knocked down and subsequent studies were performed to analyze cell proliferation, migration, invasion, and apoptosis. Mechanistic studies that explored the epidermal growth factor receptor (EGFR) pathway were performed by Western blot.

RESULTS

Multivariate analysis identified five metabolic genes that were differentially expressed between G3 and G4 stroma (P < .05). Functional validation studies revealed that knockdown of PGRMC1 and HSD17B4 significantly decreased cell proliferation, migration, and invasion, and increased apoptosis in PC3 and DU145 cells. Mechanistic studies showed that these effects, after PGRMC1 knockdown, were possibly mediated through alterations in downstream components of the EGFR, protein kinase B, and nuclear factor kappa-light-chain-enhancer of activated B cells pathways.

CONCLUSION

The following study provides evidence supporting the use of metabolic genes PGRMC1 and HSD17B4 as a prognostic biomarker for the distinction between G3 and G4 prostate cancers.

Dietary tryptophan links encephalogenicity of autoreactive T cells with gut microbial ecology

The interaction between the mammalian host and its resident gut microbiota is known to license adaptive immune responses. Nutritional constituents strongly influence composition and functional properties of the intestinal microbial communities. Here, we report that omission of a single essential amino acid - tryptophan – from the diet abrogates CNS autoimmunity in a mouse model of multiple sclerosis. Dietary tryptophan restriction results in impaired encephalitogenic T cell responses and is accompanied by a mild intestinal inflammatory response and a profound phenotypic shift of gut microbiota. Protective effects of dietary tryptophan restriction are abrogated in germ-free mice, but are independent of canonical host sensors of intracellular tryptophan metabolites. We conclude that dietary tryptophan restriction alters metabolic properties of gut microbiota, which in turn have an impact on encephalitogenic T cell responses. This link between gut microbiota, dietary tryptophan and adaptive immunity may help to develop therapeutic strategies for protection from autoimmune neuroinflammation.

Food intake shapes intestinal microbiome composition, which in turn shapes adaptive immune responses. Here the authors show that dietary tryptophan restriction (DTR) protects mice from subsequent autoimmune neuropathology challenge by altering intestinal microbiota, highlighting the potential of diet-regulated microbiota to prevent immune pathology.

A Distinct Pretreatment Immune Gene Signature in Lentigo Maligna Is Associated with Imiquimod Response

Lentigo maligna (LM) is a common subtype of in situ melanoma on chronically sun-exposed skin, particularly the head and neck of older patients. Although surgery is the standard treatment, there is associated morbidity, and options such as imiquimod cream or radiotherapy may be used if surgery is refused or inappropriate. Complete response rates following imiquimod treatment are variable in the literature. The aim of this study was to evaluate the host immune response both before and following treatment with imiquimod to better identify likely responders. Paired pre- and post-imiquimod treatment specimens were available for 27 patients. Patients were treated with imiquimod 5 days per week for 12 weeks; at 16 weeks, lesions were excised for histological assessment. Of the 27 patients, 16 were responders and 11 failed to clear the disease. PDL1 protein expression was increased, accompanied by a unique gene signature in lesions from patients that subsequently histologically cleared LM by 16 weeks. This comprised 57 upregulated immune genes in signaling networks for antigen presentation, type I interferon signaling, and T-cell activation. This may represent an early responder group to imiquimod, and this unique gene signature potentially can be used as a biomarker of LM response to imiquimod.

TLR5 decoy receptor as a novel anti-amyloid therapeutic for Alzheimer's disease

Chakrabarty et al. show that human TLR5 ectodomain reduces amyloid β (Aβ) plaques by direct interaction with Aβ, demonstrating the feasibility of such immune decoy receptor strategies as potential biotherapies in Alzheimer’s disease.

There is considerable interest in harnessing innate immunity to treat Alzheimer’s disease (AD). Here, we explore whether a decoy receptor strategy using the ectodomain of select TLRs has therapeutic potential in AD. AAV-mediated expression of human TLR5 ectodomain (sTLR5) alone or fused to human IgG4 Fc (sTLR5Fc) results in robust attenuation of amyloid β (Aβ) accumulation in a mouse model of Alzheimer-type Aβ pathology. sTLR5Fc binds to oligomeric and fibrillar Aβ with high affinity, forms complexes with Aβ, and blocks Aβ toxicity. Oligomeric and fibrillar Aβ modulates flagellin-mediated activation of human TLR5 but does not, by itself, activate TLR5 signaling. Genetic analysis shows that rare protein coding variants in human TLR5 may be associated with a reduced risk of AD. Further, transcriptome analysis shows altered TLR gene expression in human AD. Collectively, our data suggest that TLR5 decoy receptor–based biologics represent a novel and safe Aβ-selective class of biotherapy in AD.

Sustained Type I interferon signaling as a mechanism of resistance to PD-1 blockade

PD-1 blockade represents a major therapeutic avenue in anticancer immunotherapy. Delineating mechanisms of secondary resistance to this strategy is increasingly important. Here, we identified the deleterious role of signaling via the type I interferon (IFN) receptor in tumor and antigen presenting cells, that induced the expression of nitric oxide synthase 2 (NOS2), associated with intratumor accumulation of regulatory T cells (Treg) and myeloid cells and acquired resistance to anti-PD-1 monoclonal antibody (mAb). Sustained IFNβ transcription was observed in resistant tumors, in turn inducing PD-L1 and NOS2 expression in both tumor and dendritic cells (DC). Whereas PD-L1 was not involved in secondary resistance to anti-PD-1 mAb, pharmacological or genetic inhibition of NOS2 maintained long-term control of tumors by PD-1 blockade, through reduction of Treg and DC activation. Resistance to immunotherapies, including anti-PD-1 mAb in melanoma patients, was also correlated with the induction of a type I IFN signature. Hence, the role of type I IFN in response to PD-1 blockade should be revisited as sustained type I IFN signaling may contribute to resistance to therapy.

RCRnorm: An integrated system of random-coefficient hierarchical regression models for normalizing NanoString nCounter data

Formalin-fixed paraffin-embedded (FFPE) samples have great potential for biomarker discovery, retrospective studies, and diagnosis or prognosis of diseases. Their application, however, is hindered by the unsatisfactory performance of traditional gene expression profiling techniques on damaged RNAs. NanoString nCounter platform is well suited for profiling of FFPE samples and measures gene expression with high sensitivity, which may greatly facilitate realization of scientific and clinical values of FFPE samples. However, methodological development for normalization, a critical step when analyzing this type of data, is far behind. Existing methods designed for the platform use information from different types of internal controls separately and rely on an overly-simplified assumption that expression of housekeeping genes is constant across samples for global scaling. Thus, these methods are not optimized for the nCounter system, not mentioning that they were not developed for FFPE samples. We construct an integrated system of random-coefficient hierarchical regression models to capture main patterns and characteristics observed from NanoString data of FFPE samples, and develop a Bayesian approach to estimate parameters and normalize gene expression across samples. Our method, labeled RCRnorm, incorporates information from all aspects of the experimental design and simultaneously removes biases from various sources. It eliminates the unrealistic assumption on housekeeping genes and offers great interpretability. Furthermore, it is applicable to freshly frozen or like samples that can be generally viewed as a reduced case of FFPE samples. Simulation and applications showed the superior performance of RCRnorm.

NACHO: an R package for quality control of NanoString nCounter data

SUMMARY

The NanoStringTM nCounter® is a platform for the targeted quantification of expression data in biofluids and tissues. While software by the manufacturer is available in addition to third parties packages, they do not provide a complete quality control (QC) pipeline. Here, we present NACHO ('NAnostring quality Control dasHbOard'), a comprehensive QC R-package. The package consists of three subsequent steps: summarize, visualize and normalize. The summarize function collects all the relevant data and stores it in a tidy format, the visualize function initiates a dashboard with plots of the relevant QC outcomes. It contains QC metrics that are measured by default by the manufacturer, but also calculates other insightful measures, including the scaling factors that are needed in the normalization step. In this normalization step, different normalization methods can be chosen to optimally preprocess data. Together, NACHO is a comprehensive method that optimizes insight and preprocessing of nCounter® data.

AVAILABILITY AND IMPLEMENTATION

NACHO is available as an R-package on CRAN and the development version on GitHub https://github.com/mcanouil/NACHO.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Multigene Model for Predicting Tumor Responsiveness After Preoperative Chemoradiotherapy for Rectal Cancer

PURPOSE

Although preoperative chemoradiotherapy (PCRT) is regarded as a standard treatment for locally advanced rectal cancer, there is no reliable biomarker for predicting responsiveness to PCRT. We aimed to develop a biomarker model for predicting response to PCRT.

METHODS AND MATERIALS

We included 184 patients who received PCRT followed by surgical resection and categorized them as good responders (complete or near-complete regression) or poor responders (all other patients). Candidate gene mRNAs were isolated from formalin-fixed paraffin-embedded tumor specimens and analyzed using the NanoString nCounter gene expression assay. Stepwise logistic regression analysis was used to select genes in discovery and training phases. A quantitative radio-responsiveness prediction model was developed and validated using internal cross-validation groups, and the model's predictive value was assessed based on the area under the receiver operating characteristic curve (AUC).

RESULTS

By comparing the gene expressions between good and poor responders, we created a multigene mRNA model using FZD9, HRAS, ITGA7, MECOM, MMP3, NKD1, PIK3CD, and PRKCB. This panel showed good ability to predict treatment response (AUC: 0.846 for the whole data set). Internal cross-validation was performed to evaluate the model's predictive stability among 3 cohorts, which provided AUC values of 0.808-0.909. The satisfactory diagnostic performance of the radio-response prediction index persisted regardless of other clinicopathologic features such as clinical T or N stage, interval between radiation and surgery, and pretreatment carcinoembryonic antigen levels (P = .001, 95% CI, 0.686-0.905).

CONCLUSIONS

We developed a multigene mRNA-based biomarker model that allows prediction of rectal cancer response to PCRT, which may help identify patients who will benefit most from PCRT.

Time-Dependent Changes in Microglia Transcriptional Networks Following Traumatic Brain Injury

The neuroinflammatory response to traumatic brain injury (TBI) is critical to both neurotoxicity and neuroprotection, and has been proposed as a potentially modifiable driver of secondary injury in animal and human studies. Attempts to broadly target immune activation have been unsuccessful in improving outcomes, in part because the precise cellular and molecular mechanisms driving injury and outcome at acute, subacute, and chronic time points after TBI remain poorly defined. Microglia play a critical role in neuroinflammation and their persistent activation may contribute to long-term functional deficits. Activated microglia are characterized by morphological transformation and transcriptomic changes associated with specific inflammatory states. We analyzed the temporal course of changes in inflammatory genes of microglia isolated from injured brains at 2, 14, and 60 days after controlled cortical impact (CCI) in mice, a well-established model of focal cerebral contusion. We identified a time dependent, injury-associated change in the microglial gene expression profile toward a reduced ability to sense tissue damage, perform housekeeping, and maintain homeostasis in the early stages following CCI, with recovery and transition to a specialized inflammatory state over time. This later state starts at 14 days post-injury and is characterized by a biphasic pattern of IFNγ, IL-4, and IL-10 gene expression changes, with concurrent proinflammatory and anti-inflammatory gene changes. Our transcriptomic data sets are an important step to understand microglial role in TBI pathogenesis at the molecular level and identify common pathways that affect outcome. More studies to evaluate gene expression at the single cell level and focusing on subacute and chronic timepoint are warranted.