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

Ketogenic Diet and microRNAs Linked to Antioxidant Biochemical Homeostasis

Recently, we demonstrated the capability of the ketogenic diet (KD) to influence the microRNA (miR) expression profile. Here, we report that KD is able to normalize miR expression in obese subjects when compared with lean subjects. By applying two different bioinformatics tools, we found that, amongst the miRs returning to comparable levels in lean subjects, four of them are linked to antioxidant biochemical pathways specifically, and the others are linked to both antioxidant and anti-inflammatory biochemical pathways. Of particular interest is the upregulation of hsa-miR-30a-5p, which correlates with the decrease of catalase expression protein in red blood cells.

Effects of Age and Immune Landscape on Outcome in HER2-Positive Breast Cancer in the NCCTG N9831 (Alliance) and NSABP B-31 (NRG) Trials

PURPOSE

Young age has been shown to be an independent predictor of poor outcome in breast cancer. In HER2-positive breast cancer, the effects of aging remain largely unknown.

EXPERIMENTAL DESIGN

A total of 4,547 patients were included [3,132 from North Central Cancer Treatment Group (NCCTG) N9831 and 1,415 from National Surgical Adjuvant Breast and Bowel Project (NSABP) B-31]. Pathologic stromal tumor-infiltrating lymphocyte (sTIL) and molecular tumor infiltrating lymphocyte (mTIL) signatures were evaluated.

RESULTS

In NCCTG N9831, comparable benefit of trastuzumab was observed in all patients [age ≤ 40; HR, 0.43; 95% confidence interval (CI), 0.28-0.66; P < 0.001; and age > 40; HR, 0.56; 95% CI, 0.45-0.69; P < 0.001]. Similar results were observed in NSABP B-31 (age ≤ 40; HR, 0.45; 95% CI, 0.29-0.68; P < 0.001; and age > 40; HR, 0.42; 95% CI, 0.33-0.54; P < 0.001). Among patients who received chemotherapy alone, younger age was associated with poor outcome in the hormone receptor-positive subset, but not the hormone receptor-negative subset, in both trials. Although there was no association between sTILs and age, a small, but significant increase in mTIL CD45 and some immune subset signatures were observed. Among patients who received chemotherapy alone, patients over 40 years of age with lymphocyte-predominant breast cancer had excellent outcome, with 95% remaining recurrence free at 15 years.

CONCLUSIONS

Among patients treated with trastuzumab, there was no significant difference in outcome related to age. Our study suggests that trastuzumab can negate the poor prognosis associated with young age.

Discovery and Validation of Novel Biomarkers for Detection of Epithelial Ovarian Cancer

Detection of epithelial ovarian cancer (EOC) poses a critical medical challenge. However, novel biomarkers for diagnosis remain to be discovered. Therefore, innovative approaches are of the utmost importance for patient outcome. Here, we present a concept for blood-based biomarker discovery, investigating both epithelial and specifically stromal compartments, which have been neglected in search for novel candidates. We queried gene expression profiles of EOC including microdissected epithelium and adjacent stroma from benign and malignant tumours. Genes significantly differentially expressed within either the epithelial or the stromal compartments were retrieved. The expression of genes whose products are secreted yet absent in the blood of healthy donors were validated in tissue and blood from patients with pelvic mass by NanoString analysis. Results were confirmed by the comprehensive gene expression database, CSIOVDB (Ovarian cancer database of Cancer Science Institute Singapore). The top 25% of candidate genes were explored for their biomarker potential, and twelve were able to discriminate between benign and malignant tumours on transcript levels (p < 0.05). Among them T-cell differentiation protein myelin and lymphocyte (MAL), aurora kinase A (AURKA), stroma-derived candidates versican (VCAN), and syndecan-3 (SDC), which performed significantly better than the recently reported biomarker fibroblast growth factor 18 (FGF18) to discern malignant from benign conditions. Furthermore, elevated MAL and AURKA expression levels correlated significantly with a poor prognosis. We identified promising novel candidates and found the stroma of EOC to be a suitable compartment for biomarker discovery.

Characterization of molecular scores and gene expression signatures in primary breast cancer, local recurrences and brain metastases

Background

Breast cancer is a leading cause of cancer-related death in women worldwide. Despite extensive studies in all areas of basic, clinical and applied research, accurate prognosis remains elusive, thus leading to overtreatment of many patients. Diagnosis could be improved by introducing multigene molecular scores in standard clinical practice. Several tests that work with formalin-fixed tissue have become routine. Molecular scores usually include several genes representing processes, response to oestrogens, progestogens and human epidermal growth factor receptor 2 (Her2), respectively, which are combined additively in single values. These multi-gene scores have the advantage of being more robust and reproducible than single-gene scores. Their utility may be further enhanced by combining them with classical diagnostic parameters. Here, we present an exploratory study comparing the RISK and research versions of Oncotype DX recurrence score (RS), Prosigna Risk of Recurrence (ROR) and EndoPredict (EP) with respect to their prognostic potential for ipsilateral recurrence and/or distant relapse in brain, and we compared the scores to the intrinsic subtypes based on PAM50.

Methods

RNA was extracted from formalin-fixed, paraffin-embedded (FFPE) tissue cores of primary tumours, local recurrences and brain metastases. Gene expression was measured on a NanoString nCounter Analysis System. Intrinsic subtypes and molecular scores were computed according to published literature and RISK, RS, ROR and EP were compared against each other and to the intrinsic subtypes Luminal A (lumA), Luminal B (lumB), Her2-enriched (Her2↑), Basal-like (basal), and Normal-like (normal) of PAM50. Local recurrences and brain metastases were compared to their corresponding primary tumours.

Results

All four molecular scores were highly correlated. Highest correlations were observed among genes related to proliferation while lower correlations were found among oestrogen-related genes. The scores were significantly higher in primary tumours progressing to brain metastases as compared to recurrence-free primary tumours and primary tumours that relapsed as local recurrences.

Conclusions

RISK and ROR-P are prognostic for primary tumours metastasizing to the brain. All four scores, RISK, RS, EP and ROR-P failed to discriminate between primary tumours that remained recurrence-free and primary tumours relapsing as local recurrences.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-5752-8) contains supplementary material, which is available to authorized users.

Ability of a urine gene expression classifier to reduce the number of follow-up cystoscopies in bladder cancer patients

This study aimed to improve our previous urine gene expression classifiers focusing on the detection of non-high-risk non-muscle-invasive bladder cancer (NMIBC), and develop a new classifier able to decrease the frequency of cystoscopies during bladder cancer (BC) patients' surveillance. A total of 597 urines from BC patients, controls and patients in follow-up for BC (PFBC) were included. The study has 3 phases. In the urinary biomarker discovery phase, 84 urines from BC and control patients were retrospectively included and analyzed by Ribonucleic Acid (RNA) sequencing. In the classifier development phase, a total of 132 selected genes from previous phase were evaluated by nCounter in 214 prospectively collected urines from PFBC (98 with tumor). A diagnostic classifier was generated by logistic regression. Finally, in the classifier validation phase, a multicentric and international cohort of 248 urines (134 BC and 114 nonrecurrent PFBC) was used to validate classifier performance. A total of 521 genes were found differentially expressed between non-high-risk NMIBC samples and all other groups (P < 0.05). An 8-gene diagnostic classifier with an area under curve (AUC) of 0.893 was developed. Validation of this classifier in a cohort of PFBC achieved an overall sensitivity (SN) and a negative predictive value (NPV) of 96% and 97%, respectively (AUC = 0.823). Notably, this accuracy was maintained in non-high-risk NMIBC group (SN = 94%; NPV = 98%). In conclusion, this 8-gene expression classifier has high SN and NPV in a real clinical scenario. The use of this classifier can reduce the number of follow-up cystoscopies in PFBC, although assessing its final place in clinical setting is necessary.

Expression Concordance of 325 Novel RNA Biomarkers between Data Generated by NanoString nCounter and Affymetrix GeneChip

Background

With the development of new drug combinations and targeted treatments for multiple types of cancer, the ability to stratify categories of patient populations and to develop companion diagnostics has become increasingly important. A panel of 325 RNA biomarkers was selected based on cancer-related biological processes of healthy cells and gene expression changes over time during nonmalignant epithelial cell organization. This "cancer in reverse" approach resulted in a panel of biomarkers relevant for at least 7 cancer types, providing gene expression profiles representing key cellular signaling pathways beyond mutations in "driver genes." Objective. To further investigate this biomarker panel, the objective of the current study is to (1) validate the assay reproducibility for the 325 RNA biomarkers and (2) compare gene expression profiles side by side using two technology platforms.

Methods and Results

We have mapped the 325 RNA transcripts and in a custom NanoString nCounter expression panel to be compared to all potential probe sets in the Affymetrix Human Genome U133 Plus 2.0. The experiments were conducted with 10 unique biological formalin-fixed paraffin-embedded (FFPE) breast tumor samples. Each site extracted RNA from four sections of 10-micron thick FFPE tissue over three different days by two different operators using an optimized standard operating procedure and quality control criteria. Samples were analyzed using mas5 in BioConductor and NanoStringNorm in R. Pearson correlation showed reproducibility between sites for all 60 samples with r = 0.995 for Affymetrix and r = 0.999 for NanoString. Correlation in multiple days and multiple users was for Affymetrix r = (0.962 - 0.999) and for NanoString r = (0.982 - 0.991).

Conclusion

The 325 RNA biomarkers showed reproducibility in two technology platforms with moderate to high concordance. Future directions include performing clinical validation studies and generating rationale for patient selection in clinical trials using the technically validated assay.

In vitro assessment of mechanistic events induced by structurally related chemical rubber sensitizers

Allergic contact dermatitis (ACD) is one of the most common forms of immunotoxicity, and increased understanding of how chemicals trigger these adverse reactions is needed in order to treat or design testing strategies to identify and subsequently avoid exposure to such substances. In this study, we investigated the cellular response induced by rubber chemicals in a dendritic cell (DC) model, focusing on the structurally similar chemicals diethylthiocarbamylbenzothiazole sulfide and dimethylthiocarbamylbenzothiazole sulfide, with regard to regulation of microRNA, and messenger RNA expression. Only a few miRNAs were found to be commonly regulated by both rubber chemicals, among them miR1973, while the overall miRNA expression profiles were diverse. Similarly, out of approximately 500 differentially regulated transcripts for each chemical, about 60% overlapped, while remaining were unique. The pathways predicted to be enriched in the cell model by stimulation with the rubber chemicals were linked to immunological events, relevant in the context of ACD. These results suggest that small structural differences can trigger specific activation of the immune system in response to chemicals. The here presented mechanistic data can be valuable in explaining the immunotoxicological events in DC activation after exposure to skin sensitizing chemicals, and can contribute to understanding, preventing and treating ACD.

High serum CRP influences myocardial miRNA profiles in ischemia-reperfusion injury of rat heart

Objective

Prognosis of myocardial infarction tends to be worse when serum C-reactive protein (CRP) level is high. miRNAs are also known to be involved in different pathogeneses of heart diseases such as myocardial infarction. However, how CRP is involved in myocardial infarction has not been fully elucidated. We hypothesized that serum CRP changes the miRNA profile during ischemia-reperfusion injury (IRI) of the myocardium. To confirm this hypothesis, we performed global miRNA expression profiling of myocardium using IRI and CRP infusion rat model.

Methods

After ligation of the coronary artery of rat hearts, human serum CRP was intravenously injected, and reperfusion was performed (I/R+CRP group, n = 6). Control group consisted of the sham group (n = 3), IV CRP infusion group (CRP only, n = 3), and the I/R-only group (I/R only, n = 5). We evaluated 423 miRNA expression in non-ischemic areas and areas at risk (AAR) of each group using NanoString nCounter miRNA expression assay.

Results

MiR-124 was downregulated in non-ischemic myocardium in CRP-only group. In AAR, 7 miRNAs were commonly upregulated in both I/R-only and I/R+CRP groups. And additional 6 miRNAs were upregulated in the I/R+CRP group (miR-33, miR-409-3p, miR-384-3p, miR-3562, miR-101a, and miR-340-5p). Similarly, in the non-ischemic areas, 6 miRNAs were commonly upregulated in both I/R-only and I/R+CRP groups, and additional 5 miRNAs changed in the I/R+CRP group (upregulation of miR-3559-5p, miR-499, and miR-21 and downregulation of miR-500 and miR-532-3p).

Conclusion

We showed that when serum CRP level is high, IRI results in multiple miRNA profile changes not only in ischemic areas but also in non-ischemic myocardium. Our results may provide a strong basis for studying the role of CRP and miRNAs in ischemic heart disease.

VIGLA-M: visual gene expression data analytics

Background

The analysis of gene expression levels is used in many clinical studies to know how patients evolve or to find new genetic biomarkers that could help in clinical decision making. However, the techniques and software available for these analyses are not intended for physicians, but for geneticists. However, enabling physicians to make initial discoveries on these data would benefit in the clinical assay development.

Results

Melanoma is a highly immunogenic tumor. Therefore, in recent years physicians have incorporated immune system altering drugs into their therapeutic arsenal against this disease, revolutionizing the treatment of patients with an advanced stage of the cancer. This has led us to explore and deepen our knowledge of the immunology surrounding melanoma, in order to optimize the approach. Within this project we have developed a database for collecting relevant clinical information for melanoma patients, including the storage of patient gene expression levels obtained from the NanoString platform (several samples are taken from each patient). The Immune Profiling Panel is used in this case. This database is being exploited through the analysis of the different expression profiles of the patients. This analysis is being done with Python, and a parallel version of the algorithms is available with Apache Spark to provide scalability as needed.

Conclusions

VIGLA-M, the visual analysis tool for gene expression levels in melanoma patients is available at http://khaos.uma.es/melanoma/. The platform with real clinical data can be accessed with a demo user account, physician, using password physician_test_7634 (if you encounter any problems, contact us at this email address: mailto: khaos@lcc.uma.es). The initial results of the analysis of gene expression levels using these tools are providing first insights into the patients’ evolution. These results are promising, but larger scale tests must be developed once new patients have been sequenced, to discover new genetic biomarkers.

Counter-regulation of regulatory T cells by autoreactive CD8+ T cells in rheumatoid arthritis

The mechanisms whereby autoreactive T cells escape peripheral tolerance establishing thus autoimmune diseases in humans remain an unresolved question. Here, we demonstrate that autoreactive polyfunctional CD8⁺ T cells recognizing self-antigens (i.e., vimentin, actin cytoplasmic 1, or non-muscle myosin heavy chain 9 epitopes) with high avidity, counter-regulate Tregs by killing them, in a consistent percentage of rheumatoid arthritis (RA) patients. Indeed, these CD8⁺ T cells express a phenotype and gene profile of effector (eff) cells and, upon antigen-specific activation, kill Tregs indirectly in an NKG2D-dependent bystander fashion in vitro. This data provides a mechanistic basis for the finding showing that AE-specific (CD107a⁺) CD8⁺ T killer cells correlate, directly with the disease activity score, and inversely with the percentage of activated Tregs, in both steady state and follow-up studies in vivo. In addition, multiplex immunofluorescence imaging analyses of inflamed synovial tissues in vivo show that a remarkable number of CD8⁺ T cells express granzyme-B and selectively contact FOXP3⁺ Tregs, some of which are in an apoptotic state, validating hence the possibility that CD8⁺ Teff cells can counteract neighboring Tregs within inflamed tissues, by killing them. Alternatively, the disease activity score of a different subset of patients is correlated with the expansion of a peculiar subpopulation of autoreactive low avidity, partially-activated (pa)CD8⁺ T cells that, despite they conserve the conventional naïve (N) phenotype, produce high levels of tumor necrosis factor (TNF)-α and exhibit a gene expression signature of a progressive activation state. Tregs directly correlate with the expansion of this autoreactive (low avidity) paCD8⁺ TN cell subset in vivo, and efficiently control their differentiation rather their proliferation in vitro. Interestingly, autoreactive high avidity CD8⁺ Teff cells or low avidity paCD8⁺ TN cells are significantly expanded in RA patients who would become non-responders or patients who would become responders to TNF-α inhibitor therapy, respectively. These data provide evidence of a previously undescribed role of such mechanisms in the progression and therapy of RA.

Time resolved gene expression analysis during tamoxifen adaption of MCF-7 cells identifies long non-coding RNAs with prognostic impact

Acquired tamoxifen resistance is a persistent problem for the treatment of estrogen receptor positive, premenopausal breast cancer patients and predictive biomarkers are still elusive. We here analyzed gene expression changes in a cellular model to identify early and late changes upon tamoxifen exposure and thereby novel prognostic biomarkers. Estrogen receptor positive MCF-7 cells were incubated with 4OH-tamoxifen (10 nM) and gene expression analyzed by array hybridization during 12 weeks. Array results were confirmed by nCounter- and qRT-PCR technique. Pathway enrichment analysis revealed that early responses concerned mainly amine synthesis and NRF2-related signaling and evolved into a stable gene expression pattern within 4 weeks characterized by changes in glucuronidation-, estrogen metabolism-, nuclear receptor- and interferon signaling pathways. As a large number of long non coding RNAs was subject to regulation, we investigated 5 of these (linc01213, linc00632 linc0992, LOC101929547 and XR_133213) in more detail. From these, only linc01213 was upregulated but all were less abundant in estrogen-receptor negative cell lines (MDA-MB 231, SKBR-3 and UACC3199). In a web-based survival analysis linc01213 and linc00632 turned out to have prognostic impact. Linc01213 was investigated further by plasmid-mediated over-expression as well as siRNA down-regulation in MCF-7 cells. Nevertheless, this had no effect on proliferation or expression of tamoxifen regulated genes, but migration was increased. In conclusion, the cellular model identified a set of lincRNAs with prognostic relevance for breast cancer. One of these, linc01213 although regulated by 4OH-tamoxifen, is not a central regulator of tamoxifen adaption, but interferes with the regulation of migration.

Immunomodulatory effects of antipsychotic treatment on gene expression in first-episode psychosis

Previous studies suggest immunological alterations in patients with first-episode psychosis (FEP). Some studies show that antipsychotic compounds may cause immunomodulatory effects. To evaluate the immunological changes and the possible immunomodulatory effects in FEP, we recruited patients with FEP (n = 67) and matched controls (n = 38), aged 18-40 years, from the catchment area of the Helsinki University Hospital and the City of Helsinki, Finland. Fasting peripheral blood samples were collected between 8 and 10 a.m. in 10 ml PAXgene tubes. We applied the NanoString nCounter in-solution hybridization technology to determine gene expression levels of 147 candidate genes reflecting activation of the immune system. Cases had higher gene expression levels of BDKRB1 and SPP1/osteopontin compared with controls. Of the individual medications used as monotherapy, risperidone was associated with a statistically significant upregulation of 11 immune system genes, including cytokines and cytokine receptors (SPP1, IL1R1, IL1R2), pattern recognition molecules (TLR1, TLR2 and TLR6, dectin-1/CLEC7A), molecules involved in apoptosis (FAS), and some other molecules with functions in immune activation (BDKRB1, IGF1R, CR1). In conclusion, risperidone possessed strong immunomodulatory properties affecting mainly innate immune response in FEP patients, whereas the observed effects of quetiapine and olanzapine were only marginal. Our results further emphasize the importance of understanding the immunomodulatory mechanisms of antipsychotic treatment, especially in terms of specific compounds, doses and duration of medication in patients with severe mental illness. Future studies should evaluate the response pre- and post-treatment, and the possible role of this inflammatory activation for the progression of psychiatric and metabolic symptoms.

BPG: Seamless, automated and interactive visualization of scientific data

BACKGROUND

We introduce BPG, a framework for generating publication-quality, highly-customizable plots in the R statistical environment.

RESULTS

This open-source package includes multiple methods of displaying high-dimensional datasets and facilitates generation of complex multi-panel figures, making it suitable for complex datasets. A web-based interactive tool allows online figure customization, from which R code can be downloaded for integration with computational pipelines.

CONCLUSION

BPG provides a new approach for linking interactive and scripted data visualization and is available at http://labs.oicr.on.ca/boutros-lab/software/bpg or via CRAN at https://cran.r-project.org/web/packages/BoutrosLab.plotting.general.

Ketogenic Diet Acts on Body Remodeling and MicroRNAs Expression Profile

BACKGROUND

The Ketogenic Diet (KD) promotes metabolic changes and optimizes energy metabolism. It is unknown if microRNAs (miRs) are influenced by KD in obese subjects. The screening of circulating miRs was performed with the FDA approved platform n-counter flex and blood biochemical parameters were dosed by ADVIA 1800.

OBJECTIVES

The aim of this study was to evaluate mir profile under 6 weeks of biphasic KD in obese subjects. We enrolled 36 obese subjects (18 females and 18 males) in stage 1 of Edmonton Obesity Staging System (EOSS) parameter.

RESULT

Any correlation was found between biochemical parameter and three miRs, hsa-let-7b-5p, hsa-miR-143-3p and hsa-miR-504-5p influenced in an equal manner in both sexes. The KD resulted safe and ameliorate both biochemical and anthropometric factors in obese subjects re-collocating them into stage 0 of EOSS parameters.

CONCLUSION

The miRs herein identified under KD might be a useful tool to monitor low carbohydrate nutritional regimens which reflect indirectly the regulatory biochemical mechanisms and cell signaling that orchestrate metabolic and signaling pathways.

NanoStringNormCNV: pre-processing of NanoString CNV data

Summary

The NanoString System is a well-established technology for measuring RNA and DNA abundance. Although it can estimate copy number variation, relatively few tools support analysis of these data. To address this gap, we created NanoStringNormCNV, an R package for pre-processing and copy number variant calling from NanoString data. This package implements algorithms for pre-processing, quality-control, normalization and copy number variation detection. A series of reporting and data visualization methods support exploratory analyses. To demonstrate its utility, we apply it to a new dataset of 96 genes profiled on 41 prostate tumour and 24 matched normal samples.

Availability and implementation

NanoStringNormCNV is implemented in R and is freely available at http://labs.oicr.on.ca/boutros-lab/software/nanostringnormcnv.

Contact

paul.boutros@oicr.on.ca.

Supplementary information

Supplementary data are available at Bioinformatics online.

Evidence of Chronic Allograft Injury in Liver Biopsies From Long-term Pediatric Recipients of Liver Transplants

BACKGROUND & AIMS

A substantial proportion of pediatric liver transplant recipients develop subclinical chronic allograft injury. We studied whether there are distinct patterns of injury based on histopathologic features and identified associated immunologic profiles.

METHODS

We conducted a cross-sectional study of 157 stable, long-term pediatric recipients of transplanted livers (70 boys; > 6 years old at time of transplantation; mean, 8.9 ± 3.46 years after liver transplantation) who underwent liver biopsy analysis from August 13, 2012, through May 1, 2014. Participants had received livers from a living or deceased donor and had consistently normal results from liver tests. Liver biopsy specimens were scored by a central pathologist; an unsupervised hierarchical cluster analysis of histologic features was used to sort biopsy samples into 3 clusters. We conducted transcriptional and cytometric analyses of liver tissue samples and performed a systems biology analysis that incorporated clinical, serologic, histologic, and transcriptional data.

RESULTS

The mean level of alanine aminotransferase in participants was 27.6 ± 14.57 U/L, and the mean level of γ-glutamyl transferase was 17.4 ± 7.93 U/L. Cluster 1 was characterized by interface activity (n = 34), cluster 2 was characterized by periportal or perivenular fibrosis without interface activity (n = 45), and cluster 3 had neither feature (n = 78). We identified a module of genes whose expression correlated with levels of alanine aminotransferase, class II donor-specific antibody, portal inflammation, interface activity, perivenular inflammation, portal and perivenular fibrosis, and cluster assignment. The module was enriched in genes that regulate T-cell-mediated rejection (TCMR) of liver and other transplanted organs. Functional pathway analysis showed overrepresentation of TCMR gene sets for cluster 1 but not clusters 2 or 3.

CONCLUSION

In an analysis of biopsies from an apparently homogeneous group of stable, long-term pediatric liver transplant recipients with consistently normal liver test results, we found evidence of chronic graft injury (inflammation and/or fibrosis). Biopsy samples with interface activity had a gene expression pattern associated with TCMR.

Aging, melatonin biosynthesis, and circadian clockworks in the gastrointestinal system of the laboratory mouse

The gastrointestinal (GI) system is vital in its capacities for nutrient and water uptake, immune function, metabolism and detoxification, and stem-cell derived regeneration. Of significance to human health are a myriad of GI disorders associated with aging that integrate with the circadian clock. Here we present data from three groups of mice: young (3 mo old), middle aged (12 mo old), and old aged (24 mo old). Small intestine and colon samples taken every 4 h under light-dark (LD) conditions were assayed for gene expression related to molecular circadian rhythmicity, transcription, cell signaling, and immune function. Transcripts related to melatonin biosynthesis and signaling, as well as melatonin content from stool, were also included, as GI melatonin and aging have been associated in contexts outside of the circadian clock. With respect to circadian genes, the data here are congruent with data from other peripheral tissues: age does not affect the rhythmic expression of core clock genes in the gut. The same can be said for several clock-controlled transcripts. In contrast, diurnal patterns in the expression of nitric oxide synthase 1 and of immune factors irak4 and interleukin-8 were observed in the colon of young mice that were lost in middle-aged and aged animals. Furthermore, the diurnal pattern of melatonin synthesis genes was altered by age, and stool melatonin levels showed significant decline between young mice and aged cohorts. These data expand the evidence for the persistence of the circadian clock throughout the aging process and highlight its importance to health.

Pathway-based subnetworks enable cross-disease biomarker discovery

Biomarkers lie at the heart of precision medicine. Surprisingly, while rapid genomic profiling is becoming ubiquitous, the development of biomarkers usually involves the application of bespoke techniques that cannot be directly applied to other datasets. There is an urgent need for a systematic methodology to create biologically-interpretable molecular models that robustly predict key phenotypes. Here we present SIMMS (Subnetwork Integration for Multi-Modal Signatures): an algorithm that fragments pathways into functional modules and uses these to predict phenotypes. We apply SIMMS to multiple data types across five diseases, and in each it reproducibly identifies known and novel subtypes, and makes superior predictions to the best bespoke approaches. To demonstrate its ability on a new dataset, we profile 33 genes/nodes of the PI3K pathway in 1734 FFPE breast tumors and create a four-subnetwork prediction model. This model out-performs a clinically-validated molecular test in an independent cohort of 1742 patients. SIMMS is generic and enables systematic data integration for robust biomarker discovery.

In situ Immune Signatures and Microbial Load at the Nasopharyngeal Interface in Children With Acute Respiratory Infection

Acute respiratory infection (ARI) is the most frequent cause for hospitalization in infants and young children. Using multiplexed nCounter technology to digitally quantify 600 human mRNAs in parallel with 14 virus- and 5 bacterium-specific RNAs, we characterized viral and bacterial presence in nasopharyngeal aspirates (NPA) of 58 children with ARI and determined the corresponding in situ immune profiles. NPA contained different groups of organisms and these were classified into bacterial (n = 27), viral (n = 5), codetection [containing both viral and bacterial transcripts (n = 21), or indeterminate intermediate where microbial load is below threshold (n = 5)]. We then identified differentially expressed immune transcripts (DEITs) comparing NPAs from symptomatic children vs. healthy controls, and comparing children presenting NPAs with detectable microbial load vs. indeterminate. We observed a strong innate immune response in NPAs, due to the presence of evolutionarily conserved type I Interferon (IFN)-stimulated genes (ISG), which was correlated with total bacterial and/or viral load. In comparison with indeterminate NPAs, adaptive immunity transcripts discriminated among viral, bacterial, and codetected microbial profiles. In viral NPAs, B cell transcripts were significantly enriched among DEITs, while only type III IFN was correlated with viral load. In bacterial NPAs, myeloid cells and coinhibitory transcripts were enriched and significantly correlated with bacterial load. In conclusion, digital nCounter transcriptomics provide a microbial and immunological in situ “snapshot” of the nasopharyngeal interface in children with ARI. This enabled discrimination among viral, bacterial, codetection, and indeterminate transcripts in the samples using non-invasive sampling.

In vivo identification of Bmp2-correlation networks during fracture healing by means of a limb-specific conditional inactivation of Bmp2

Bmp2 is known to play an essential role in the initiation of fracture healing via periosteal activation. Specifically, activation and subsequent differentiation of periosteal progenitor cells requires Bmp2 signaling for activation of the osteo-chondrogenic pathway. Here, we explored the interactive transcriptional gene-gene interplays between Bmp2 and 150 known candidate genes during fracture repair. We constructed the interactive Bmp2 signaling pathways in vivo, by comparing gene expression levels prior and 24 h post femur fracture, in presence (wild type) and in absence of Bmp2 (Bmp2c/c;Prx1::cre limb-specific conditional knockout). Twenty-six differentially expressed genes (pre- vs. post-fracture), which demonstrated high correlations within each experimental condition, were used to construct the co-expression networks. Topological dynamic shifts across different co-expression networks characterized the 26 differentially expressed genes as non-redundant focal linking hubs, redundant connecting hubs, periphery genes, or non-existent. Top-ranked up- or down-regulated genes were identified and discussed. Protein-protein interactions in public databases support our findings. Thus, the co-expression networks from this study can be used for future experimental hypotheses.

Reproducibility of the NanoString 22-gene molecular subgroup assay for improved prognostic prediction of medulloblastoma

Medulloblastoma is the most frequent malignant brain tumor in children. Four medulloblastoma molecular subgroups, MB_SHH , MB_WNT , MB_GRP3 and MB_GRP4 , have been identified by integrated high-throughput platforms. Recently, a 22-gene panel NanoString-based assay was developed for medulloblastoma molecular subgrouping, but the robustness of this assay has not been widely evaluated. Mutations in the gene for human telomerase reverse transcriptase (hTERT) have been found in medulloblastomas and are associated with distinct molecular subtypes. This study aimed to implement the 22-gene panel in a Brazilian context, and to associate the molecular profile with patients' clinical-pathological features. Formalin-fixed, paraffin-embedded (FFPE) medulloblastoma samples (n = 104) from three Brazilian centers were evaluated. Expression profiling of the 22-gene panel was performed by NanoString and a Canadian series (n = 240) was applied for training phase. hTERT mutations were analyzed by PCR followed by direct Sanger sequencing and the molecular profile was associated with patients' clinicopathological features. Overall, 65% of the patients were male, average age at diagnosis was 18 years and 7% of the patients presented metastasis at diagnosis. The molecular classification was attained in 100% of the cases, with the following frequencies: MB_SHH (n = 51), MB_WNT (n = 19), MB_GRP4 (n = 19) and MB_GRP3 (n = 15). The MB_SHH and MB_GRP3 subgroups were associated with older and younger patients, respectively. The MB_GRP4 subgroup exhibited the lowest 5-year cancer-specific overall survival (OS), yet in the multivariate analysis, only metastasis at diagnosis and surgical resection were associated with OS. hTERT mutations were detected in 29% of the cases and were associated with older patients, increased hTERT expression and MB_SHH subgroup. The 22-gene panel provides a reproducible assay for molecular subgrouping of medulloblastoma FFPE samples in a routine setting and is well-suited for future clinical trials.

Apparent bias toward long gene misregulation in MeCP2 syndromes disappears after controlling for baseline variations

Recent studies have suggested that genes longer than 100 kb are more likely to be misregulated in neurological diseases associated with synaptic dysfunction, such as autism and Rett syndrome. These length-dependent transcriptional changes are modest in MeCP2-mutant samples, but, given the low sensitivity of high-throughput transcriptome profiling technology, here we re-evaluate the statistical significance of these results. We find that the apparent length-dependent trends previously observed in MeCP2 microarray and RNA-sequencing datasets disappear after estimating baseline variability from randomized control samples. This is particularly true for genes with low fold changes. We find no bias with NanoString technology, so this long gene bias seems to be particular to polymerase chain reaction amplification-based platforms. In contrast, authentic long gene effects, such as those caused by topoisomerase inhibition, can be detected even after adjustment for baseline variability. We conclude that accurate characterization of length-dependent (or other) trends requires establishing a baseline from randomized control samples.

The DNA Methyltransferase 1 (DNMT1) Controls the Shape and Dynamics of Migrating POA-Derived Interneurons Fated for the Murine Cerebral Cortex

The proliferative niches in the subpallium generate a rich cellular variety fated for diverse telencephalic regions. The embryonic preoptic area (POA) represents one of these domains giving rise to the pool of cortical GABAergic interneurons and glial cells, in addition to striatal and residual POA cells. The migration from sites of origin within the subpallium to the distant targets like the cerebral cortex, accomplished by the adoption and maintenance of a particular migratory morphology, is a critical step during interneuron development. To identify factors orchestrating this process, we performed single-cell transcriptome analysis and detected Dnmt1 expression in murine migratory GABAergic POA-derived cells. Deletion of Dnmt1 in postmitotic immature cells of the POA caused defective migration and severely diminished adult cortical interneuron numbers. We found that DNA methyltransferase 1 (DNMT1) preserves the migratory shape in part through negative regulation of Pak6, which stimulates neuritogenesis at postmigratory stages. Our data underline the importance of DNMT1 for the migration of POA-derived cells including cortical interneurons.

Expression profiling and microRNA regulation of the LKB1 pathway in young and aged lung adenocarcinoma patients

Lung cancer in young patients appears to have distinct clinicopathological features. The present study focused on the role of the serine/threonine kinase liver kinase B1 (LKB1), a known tumor suppressor gene, and its miRNA regulation in lung adenocarcinoma, particularly in young versus elderly patients. A total of 88 patients with lung adenocarcinoma were retrospectively analysed. A simultaneous quantification was performed of the expression of LKB1 mRNA and 15 microRNAs (miRNA/miRs; miRs −93, −96, −34a, −34c, −214, −33a, −30b, −145, −182, −30c, −183, −29b, −29c, −153 and −138) involved in the LKB1 pathway, as well as of 5 identified target mRNAs [cyclin D1 (CCND1), catenin β-1 (CTNNB1), lysyl oxidase (LOX), yes-associated protein 1 (YAP1) and survivin], using NanoString technology. KRAS mutations were investigated by pyrosequencing analysis. Patients ≤50 years were defined as a younger group, while patients >50 years old as an older group (n=44/group). No difference between the two groups was identified in terms of survival times analysed using the Kaplan-Meier method or KRAS mutations. Subsequently, the LKB1 signalling pathway was focused on, as a target for therapy in lung adenocarcinoma, and assessed with regards to clinicopathological features; we found that LOX levels in adenocarcinoma patients were significantly associated with histological subtype (P=0.03), stage (P<0.0001) and prognosis (P=0.02 for disease-free interval and P=0.005 for overall survival), but not with age. Furthermore, the miRNA target prediction model indicated that miR-93 and miR-30b appeared to have functional binding sites and downregulate the gene expression of LKB1 and LOX, respectively. In conclusion, young patients appeared have similar survival rates to elderly patients. The assessment of LKB1, its downstream genes and its regulation by miRNAs may have an impact on future research on lung adenocarcinoma in young and elderly patients. Further investigations will be necessary to elucidate the potential of this pathway as a novel target for therapy.

Sigmoid colon mucosal gene expression supports alterations of neuronal signaling in irritable bowel syndrome with constipation

Peripheral factors likely play a role in at least a subset of irritable bowel syndrome (IBS) patients. Few studies have investigated mucosal gene expression using an unbiased approach. Here, we performed mucosal gene profiling in a sex-balanced sample to identify relevant signaling pathways and gene networks and compare with publicly available profiling data from additional cohorts. Twenty Rome III+ IBS patients [10 IBS with constipation (IBS-C), 10 IBS with diarrhea (IBS-D), 5 men/women each), and 10 age-/sex-matched healthy controls (HCs)] underwent sigmoidoscopy with biopsy for gene microarray analysis, including differential expression, weighted gene coexpression network analysis (WGCNA), gene set enrichment analysis, and comparison with publicly available data. Expression levels of 67 genes were validated in an expanded cohort, including the above samples and 18 additional participants (6 each of IBS-C, IBS-D, HCs) using NanoString nCounter technology. There were 1,270 differentially expressed genes (FDR < 0.05) in IBS-C vs. HCs but none in IBS or IBS-D vs. HCs. WGNCA analysis identified activation of the cAMP/protein kinase A signaling pathway. Nine of 67 genes were validated by the NanoString nCounter technology (FDR < 0.05) in the expanded sample. Comparison with publicly available microarray data from the Mayo Clinic and University of Nottingham supports the reproducibility of 17 genes from the microarray analysis and three of nine genes validated by nCounter in IBS-C vs. HCs. This study supports the involvement of peripheral mechanisms in IBS-C, particularly pathways mediating neuronal signaling. NEW & NOTEWORTHY Peripheral factors play a role in the pathophysiology of irritable bowel syndrome (IBS), which, to date, has been mostly evident in IBS with diarrhea. Here, we show that sigmoid colon mucosal gene expression profiles differentiate IBS with constipation from healthy controls. These profiling data and analysis of additional cohorts also support the concept that peripheral neuronal pathways contribute to IBS pathophysiology.

Identification of a neutrophil-related gene expression signature that is enriched in adult systemic lupus erythematosus patients with active nephritis: Clinical/pathologic associations and etiologic mechanisms

Both a lack of biomarkers and relatively ineffective treatments constitute impediments to management of lupus nephritis (LN). Here we used gene expression microarrays to contrast the transcriptomic profiles of active SLE patients with and without LN to identify potential biomarkers for this condition. RNA isolated from whole peripheral blood of active SLE patients was used for transcriptomic profiling and the data analyzed by linear modeling, with corrections for multiple testing. Results were validated in a second cohort of SLE patients, using NanoString technology. The majority of genes demonstrating altered transcript abundance between patients with and without LN were neutrophil-related. Findings in the validation cohort confirmed this observation and showed that levels of RNA abundance in renal remission were similar to active patients without LN. In secondary analyses, RNA abundance correlated with disease activity, hematuria and proteinuria, but not renal biopsy changes. As abundance levels of the individual transcripts correlated strongly with each other, a composite neutrophil score was generated by summing all levels before examining additional correlations. There was a modest correlation between the neutrophil score and the blood neutrophil count, which was largely driven by the dose of glucocorticosteroids and not the proportion of low density and/or activated neutrophils. Analysis of longitudinal data revealed no correlation between baseline neutrophil score or changes over the first year of follow-up with subsequent renal flare or treatment outcomes, respectively. The findings argue that although the neutrophil score is associated with LN, its clinical utility as a biomarker may be limited.

Immunohistochemically detectable metallothionein expression in malignant pleural mesotheliomas is strongly associated with early failure to platin-based chemotherapy

Background

Malignant pleural mesothelioma (MPM) is a biologically highly aggressive tumor arising from the pleura with a dismal prognosis. Cisplatin is the drug of choice for the treatment of MPM, and carboplatin seems to have comparable efficacy. Nevertheless, cisplatin treatment results in a response rate of merely 14% and a median survival of less than seven months. Due to their role in many cellular processes, methallothioneins (MTs) have been widely studied in various cancers. The known heavy metal detoxifying effect of MT-I and MT-II may be the reason for heavy metal drug resistance of various cancers including MPM.

Methods

105 patients were retrospectively analyzed immunohistochemically for their MT expression levels. Survival analysis was done by Cox-regression, and statistical significance determined using likelihood ratio, Wald test and Score (logrank) tests.

Results

Cox-regression analyses were done in a linear and logarithmic scale revealing a significant association between expression of MT and shortened overall survival (OS) in a linear (p=0.0009) and logarithmic scale (p=0.0003). Reduced progression free survival (PFS) was also observed for MT expressing tumors (linear: p=0.0134, log: p=0.0152).

Conclusion

Since both, overall survival and progression-free survival are negatively correlated with detectable MT expression in MPM, our results indicate a possible resistance to platin-based chemotherapy associated with MT expression upregulation, found exclusively in progressive MPM samples. Initial cell culture studies suggest promoter DNA hypomethylation and expression of miRNA-566 a direct regulator of copper transporter SLC31A1 and a putative regulator of MT1A and MT2A gene expression, to be responsible for the drug resistance.

Placental MAOA expression mediates prenatal stress effects on temperament in 12-month-olds

The placenta adapts to maternal environment and its alterations may have a lasting impact on child's temperament development. Prenatal stress has been linked to both a downregulation of monoamine oxidase A (MAOA) gene expression in the placenta and to difficult temperament. Capitalizing on an ongoing longitudinal study, we analysed data from 95 mother-child dyads to investigate whether MAOA mediates the association between prenatal stress and infant temperament. Prenatal stress was defined as exposure to Superstorm Sandy (Sandy) during pregnancy. Infant temperament was measured by Infant Behaviour Questionnaire-Revised. MAOA gene expression was quantified in placenta tissue. The Smiling and Laughter subscale score was independently associated with Sandy exposure and MAOA placental gene expression. Mediation analysis confirmed that MAOA expression partially mediated the relationship between Sandy and Smiling and Laughter subscale, suggesting that in utero exposure to Sandy could induce lower frequency of smiling and laughter via downregulation of placental MAOA gene expression. These effects could compromise optimal temperamental trajectory and contribute to risk for psychological problems. Placental epigenetic markers can contribute to a multidimensional model of early intervention for high-risk children.

Transcriptional profiling reveals distinct classes of parathyroid tumors in PHPT

The clinical presentation of primary hyperparathyroidism (PHPT) varies widely, although the underlying mechanistic reasons for this disparity remain unknown. We recently reported that parathyroid tumors can be functionally segregated into two distinct groups on the basis of their relative responsiveness to ambient calcium, and that patients in these groups differ significantly in their likelihood of manifesting bone disability. To examine the molecular basis for this phenotypic variation in PHPT, we compared the global gene expression profiles of calcium-sensitive and calcium-resistant parathyroid tumors. RNAseq and proteomic analysis identified a candidate set of differentially expressed genes highly correlated with calcium-sensing capacity. Subsequent quantitative assessment of the expression levels of these genes in an independent cohort of parathyroid tumors confirmed that calcium-sensitive tumors cluster in a discrete transcriptional profile group. These data indicate that PHPT is not an etiologically monolithic disorder and suggest that divergent molecular mechanisms could drive the observed phenotypic differences in PHPT disease course, provenance, and outcome.

Timing of prenatal exposure to trauma and altered placental expressions of hypothalamic-pituitary-adrenal axis genes and genes driving neurodevelopment

Prenatal maternal stress increases the risk for negative developmental outcomes in offspring; however, the underlying biological mechanisms remain largely unexplored. In the present study, alterations in placental gene expression associated with maternal stress were examined to clarify the potential underlying epi/genetic mechanisms. Expression levels of 40 selected genes involved in regulating foetal hypothalamic-pituitary-adrenal axis and neurodevelopment were profiled in placental tissues collected from a birth cohort established around the time of Superstorm Sandy. Objective prenatal traumatic stress was defined as whether mothers were exposed to Superstorm Sandy during pregnancy. Among the 275 mother-infant dyads, 181 dyads were delivered before Superstorm Sandy (ie, Control), 66 dyads were exposed to Superstorm Sandy during the first trimester (ie, Early Exposure) and 28 were exposed to Superstorm Sandy during the second or third trimester (ie, Mid-Late Exposure). Across all trimesters, expression of HSD11B2, MAOA, ZNF507 and DYRK1A was down-regulated among those exposed to Superstorm Sandy during pregnancy. Furthermore, trimester-specific differences were also observed: exposure during early gestation was associated with down-regulation of HSD11B1 and MAOB and up-regulation of CRHBP; exposure during mid-late gestation was associated with up-regulation of SRD5A3. The findings of the present study suggest that placental gene expression may be altered in response to traumatic stress exposure during pregnancy, and the susceptibility of these genes is dependent on the time of the exposure during pregnancy. Further studies should aim to clarify the biological mechanisms that underlie trimester-specific exposure by evaluating the differential impact on offspring neurodevelopment later in childhood.

Ifngr1 and Stat1 mediated canonical Ifn-γ signaling drives nigrostriatal degeneration

Brain expression of AAV-Ifn-γ leads to reactive gliosis, nigrostriatal degeneration and midbrain calcification in wild type mice. This mouse model phenocopies idiopathic basal ganglia calcification which is associated with Parkinsonian symptoms. To understand how the nigro-striatal pathway is selectively vulnerable to Ifn-γ, we determined if the phenotype is driven by canonical signaling intermediates, Ifngr1 and Stat1. Using focused bioinformatic analysis and rotarod testing, we show that neuroinflammation and motor abnormalities precede the appearance of midbrain neuropathologies in the brains of Ifn-γ mouse model. To test whether canonical Ifn-γ signaling is a key driver of progressive nigrostriatal degeneration, we overexpressed Ifn-γ in the brains of Ifngr1-/- and Stat1-/- mice. Expression of Ifn-γ in Ifngr1-/- mice did not result in any neuroinflammation, midbrain calcinosis or nigrostriatal degenerative pathology. Interestingly, in Stat1-/- mice, Ifn-γ expression resulted in gliosis without recapitulating the neurodegenerative phenotype. Overall, our data shows that canonical Ifn-γ signaling triggers midbrain calcinosis and nigrostriatal neurodegeneration, providing mechanistic insights into cytokine-driven selective neuronal vulnerability. Our study establishes the broader relevance of inflammatory signaling in neurodegenerative diseases and can potentially identify novel immunological targets for Parkinsonian syndromes.

Folic-acid metabolism and DNA-repair phenotypes differ between neuroendocrine lung tumors and associate with aggressive subtypes, therapy resistance and outcome

Purpose

25% of all lung cancer cases are neuroendocrine (NELC) including typical (TC) and atypical carcinoid (AC), large-cell neuroendocrine (LCNEC) and small cell lung cancer (SCLC). Prognostic and predictive biomarkers are lacking.

Experimental Design

Sixty patients were used for nCounter mRNA expression analysis of the folic-acid metabolism (ATIC, DHFR, FOLR1, FPGS, GART, GGT1, SLC19A1, TYMS) and DNA-repair (ERCC1, MLH1, MSH2, MSH6, XRCC1). Phenotypic classification classified tumors (either below or above the median expression level) with respect to the folic acid metabolism or DNA repair.

Results

Expression of FOLR1, FPGS, MLH1 and TYMS (each p<0.0001) differed significantly between all four tumor types. FOLR1 and FPGS associated with tumor differentiation (both p<0.0001), spread to regional lymph nodes (FOLR1 p=0.0001 and FPGS p=0.0038), OS and PFS (FOLR1 p<0.0050 for both and FPGS p<0.0004 for OS).

Phenotypic sorting revealed the Ft-phenotype to be the most prominent expression profile in carcinoids, whereas SCLC presented nearly univocal with the fT and LCNEC with fT or ft. These results were significant for tumor subtype (p<0.0001).

Conclusions

The assessed biomarkers and phenotypes allow for risk stratification (OS, PFS), diagnostic classification and enhance the biological understanding of the different subtypes of neuroendocrine tumors revealing potential new therapy options and clarifying known resistance mechanisms.

The choroid plexus is a key cerebral invasion route for T cells after stroke

Neuroinflammation contributes substantially to stroke pathophysiology. Cerebral invasion of peripheral leukocytes-particularly T cells-has been shown to be a key event promoting inflammatory tissue damage after stroke. While previous research has focused on the vascular invasion of T cells into the ischemic brain, the choroid plexus (ChP) as an alternative cerebral T-cell invasion route after stroke has not been investigated. We here report specific accumulation of T cells in the peri-infarct cortex and detection of T cells as the predominant population in the ipsilateral ChP in mice as well as in human post-stroke autopsy samples. T-cell migration from the ChP to the peri-infarct cortex was confirmed by in vivo cell tracking of photoactivated T cells. In turn, significantly less T cells invaded the ischemic brain after photothrombotic lesion of the ipsilateral ChP and in a stroke model encompassing ChP ischemia. We detected a gradient of CCR2 ligands as the potential driving force and characterized the neuroanatomical pathway for the intracerebral migration. In summary, our study demonstrates that the ChP is a key invasion route for post-stroke cerebral T-cell invasion and describes a CCR2-ligand gradient between cortex and ChP as the potential driving mechanism for this invasion route.

Role of microRNA-33a in regulating the expression of PD-1 in lung adenocarcinoma

Background

MiRNAs are vital in functioning as either oncogenes or tumor suppressors in the cell cycle. Target transcripts for immune checkpoint molecules such as PD-1/PD-L1 and (programmed cell death-1/its ligand and cytotoxic T-lymphocyte antigen 4) have proven to be beneficial against several solid tumors, including lung adenocarcinoma.

Methods

Simultaneous quantification of the expression level of miR-33a and PD-1, PD-L1 and CTLA4 mRNAs with NanoString technology was performed in 88 lung adenocarcinoma specimens. A cohort of 323 lung adenocarcinoma patients from the cancer genome atlas (TCGA) database was further analyzed, in order to test our hypothesis. Potential interference of PD-1, PD-L1 and CTLA4 gene expression by miR-33a was predicted using the microRNA target prediction program RNA22.

Results

High miR-33a expression was significantly associated with younger (p = 0.005), female (p = 0.04), patients with low grade (p < 0.0001), early stage (p = 0.03) tumors, and better survival. The hypothesis of the involvement of miR-33a in PD-1/PD-L1/CTLA4 mechanisms was corroborated by the finding of putative miR-33a binding sites in all three genes using the RNA22 method. We found an inverse correlation between miR-33a and PD-1 levels (p = 0.01), as well as for PD-L1 (p = 0.01) and CTLA4 (p = 0.03) expression, and a significant better prognosis for patients with high miR-33a/low PD-1. TCGA database analysis confirmed that miR-33a high levels were associated with low PD-1 expression and with longer survival on a larger population.

Conclusions

Our study emphasizes the notion of a potential value of miR-33a as a favorable prognostic marker through PD-1 regulation.

Molecular analysis of high-grade serous ovarian carcinoma with and without associated serous tubal intra-epithelial carcinoma

Many high-grade serous carcinomas (HGSCs) of the pelvis are thought to originate in the distal portion of the fallopian tube. Serous tubal intra-epithelial carcinoma (STIC) lesions are the putative precursor to HGSC and identifiable in ~ 50% of advanced stage cases. To better understand the molecular etiology of HGSCs, we report a multi-center integrated genomic analysis of advanced stage tumors with and without STIC lesions and normal tissues. The most significant focal DNA SCNAs were shared between cases with and without STIC lesions. The RNA sequence and the miRNA data did not identify any clear separation between cases with and without STIC lesions. HGSCs had molecular profiles more similar to normal fallopian tube epithelium than ovarian surface epithelium or peritoneum. The data suggest that the molecular features of HGSCs with and without associated STIC lesions are mostly shared, indicating a common biologic origin, likely to be the distal fallopian tube among all cases.High-grade serous carcinomas (HGSCs) are associated with precursor lesions (STICs) in the fallopian epithelium in only half of the cases. Here the authors report the molecular analysis of HGSCs with and without associated STICs and show similar profiles supporting a common origin for all HGSCs.

Gene expression modules in primary breast cancers as risk factors for organotropic patterns of first metastatic spread: a case control study

BACKGROUND

Metastases from primary breast cancers can involve single or multiple organs at metastatic disease diagnosis. Molecular risk factors for particular patterns of metastastic spread in a clinical population are limited.

METHODS

A case-control design including 1357 primary breast cancers was used to study three distinct clinical patterns of metastasis, which occur within the first six months of metastatic disease: bone and visceral metasynchronous spread, bone-only, and visceral-only metastasis. Whole-genome expression profiles were obtained using whole genome (WG)-DASL assays from formalin-fixed paraffin-embedded (FFPE) samples. A systematic protocol was developed for handling FFPE samples together with stringent data quality controls to identify robust expression profiling data. A panel of published and novel gene sets were tested for association with these specific patterns of metastatic spread and odds ratios (ORs) were calculated.

RESULTS

Metasynchronous metastasis to bone and viscera was found in all intrinsic breast cancer subtypes, while immunohistochemically (IHC)-defined receptor status and specific IntClust subgroups were risk factors for visceral-only or bone-only first metastases. Among gene modules, those related to proliferation increased the risk of metasynchronous metastasis (OR (95% CI) = 2.3 (1.1-4.8)) and visceral-only first metastasis (OR (95% CI) = 2.5 (1.2-5.1)) but not bone-only metastasis (OR (95% CI) = 0.97 (0.56-1.7)). A 21-gene module (BV) was identified in estrogen-receptor-positive breast cancers with metasynchronous metastasis to bone and viscera (area under the curve = 0.77), and its expression increased the risk of bone and visceral metasynchronous spread in this population. BV was further orthogonally validated with NanoString nCounter in primary breast cancers, and was reproducible in their matched lymph nodes metastases and an external cohort.

CONCLUSION

This case-control study of WG-DASL global expression profiles from FFPE tumour samples, after careful quality control and RNA selection, revealed that gene modules in the primary tumour have differing risks for clinical patterns of metasynchronous first metastases. Moreover, a novel gene module was identified as a putative risk factor for metasynchronous bone and visceral first metastatic spread, with potential implications for disease monitoring and treatment planning.

Mammalian APE1 controls miRNA processing and its interactome is linked to cancer RNA metabolism

Mammalian apurinic/apyrimidinic endonuclease 1 is a DNA repair enzyme involved in genome stability and expression of genes involved in oxidative stress responses, tumor progression and chemoresistance. However, the molecular mechanisms underlying the role of apurinic/apyrimidinic endonuclease 1 in these processes are still unclear. Recent findings point to a novel role of apurinic/apyrimidinic endonuclease 1 in RNA metabolism. Through the characterization of the interactomes of apurinic/apyrimidinic endonuclease 1 with RNA and other proteins, we demonstrate here a role for apurinic/apyrimidinic endonuclease 1 in pri-miRNA processing and stability via association with the DROSHA-processing complex during genotoxic stress. We also show that endonuclease activity of apurinic/apyrimidinic endonuclease 1 is required for the processing of miR-221/222 in regulating expression of the tumor suppressor PTEN. Analysis of a cohort of different cancers supports the relevance of our findings for tumor biology. We also show that apurinic/apyrimidinic endonuclease 1 participates in RNA-interactomes and protein-interactomes involved in cancer development, thus indicating an unsuspected post-transcriptional effect on cancer genes.

APE1 plays an important role in the cellular response to oxidative stress, and mutations are linked to tumor progression and chemoresistance. Here, the authors characterize the interactions of APE1 with RNA and demonstrate a role in microRNA processing.

Prenatal Depression and Infant Temperament: The Moderating Role of Placental Gene Expression

Prior research has demonstrated the link between maternal depression during pregnancy (i.e., prenatal depression) and increased neurodevelopmental dysregulation in offspring. However, little is known about the roles of key hypothalamic-pituitary axis regulatory genes in the placenta modulating this association. This study will examine whether placental gene expression levels of 11β-hydroxysteroid dehydrogenase type 2 (HSD11B2), glucocorticoid receptor (NR3C1), and mineralocorticoid receptor (NR3C2) can help elucidate the underlying mechanisms linking prenatal depression to infant temperament, particularly in infants with high negativity and low emotion regulation. Stored placenta tissues (N = 153) were used to quantify messenger ribonucleic acid levels of HSD11B2, NR3C1, and NR3C2. Assessments of prenatal depression and infant temperament at 6 months of age were ascertained via maternal report. Results found that prenatal depression was associated with increased Negative Affectivity (p < .05) after controlling for postnatal depression and psychosocial characteristics. Furthermore, the association between prenatal depression and Negative Affectivity was moderated by gene expression levels of HSD11B2, NR3C1, and NR3C2 such that greater gene expression significantly lessened the association between prenatal depression and Negative Affectivity. Our findings suggest that individual differences in placental gene expression may be used as an early marker of susceptibility or resilience to prenatal adversity.

β-Catenin in desmoid-type fibromatosis: deep insights into the role of T41A and S45F mutations on protein structure and gene expression

Desmoid‐type fibromatosis (DF) is a rare mesenchymal lesion with high risk of local recurrence. Specific β‐catenin mutations (S45F) appeared to be related to this higher risk compared to T41A‐mutated or wild‐type (WT). We explored the influence of both mutations and WT on structure stability and affinity of β‐catenin for α‐catenin and the pattern of gene expression that may influence DF behavior. Using 33 surgically resected primary DFs harboring T41A (n = 14), S45F (n = 10), or WT (n = 9), we performed a comparative molecular analysis by protein/protein interaction modeling, gene expression by DASL microarrays, human inflammation gene panel, and assessment of immune system‐based biomarkers by immunohistochemistry. Mutated proteins were more stable than WT and formed a weaker complex with α‐catenin. Consensus unsupervised gene clustering revealed the presence of two DF group‐mutated (T41A + S45F) and WT (P = 0.0047). The gene sets ‘Inflammatory‐Defense‐Humoral Immune Response’ and ‘Antigen Binding’ were significantly enriched in T41A. The deregulation of 16 inflammation‐related genes was confirmed. Low numbers of T cells and tumor‐associated macrophages (TAM) infiltrating the tumors and low/absent PD‐1/PD‐L1 expression were also identified. We demonstrated that mutated DFs (T41A or S45F) and WT are two distinct molecular subgroups with regard to β‐catenin stability, α‐catenin affinity, and gene expression profiling. A different inflammation signature characterized the two mutated groups, suggesting mediation either by T41A or by S45F. Finally, all mutated cases showed a low number of TIL and TAM cells and a low or absent expression of PD‐1 and PD‐L1 consistent with β‐catenin activation insensitive to checkpoint blockade.

Host immune response index in gastric cancer identified by comprehensive analyses of tumor immunity

Tumor infiltrating lymphocytes (TIL) in Epstein-Barr virus (EBV)-associated/microsatellite-unstable (MSI) gastric carcinomas (GC) constitute immune-active principal cellular components of tumor microenvironment and contribute to better prognosis. With the remarkable success of cancer immunotherapies, there is an urgent need for a comprehensive understanding of tumor-immune interactions in patients with GC in the context of host immune response. To identify GC subtype-specific immune response gene set, we tested differentially expressed genes for MSI and EBV+ GC subtypes in randomly selected test set (n = 278) in merged ACRG-SMC microarray and TCGA RNA sequencing data set. We identified Host ImmunE Response index (HIERÏ) consisting of 29 immune genes classifying GC patients into robust 3 groups with prognostic significance. Immune-high cluster 1 was enriched with PD-L1^High/EBV+/MSI/TIL^High with the best clinical outcome while immune-low cluster 3 displayed worst outcome and exemplified with PD-L1^Low/EBV-/MSS. The results were validated in the same cohort (n = 279) and independent cohort (n = 181) with RNA from formalin-fixed paraffin-embedded (FFPE) tissue. Unexpectedly, nearly half of GC in cluster 1 were EBV-/MSS and 10% of cluster 3 GC were EBV+/MSI GC patients, suggesting that in addition to EBV+/MSI GC subtypes, EBV-/MSS subtype also constitutes almost half of high immune cluster and would be a good candidate for immune checkpoint inhibitor therapy. In contrary, almost 10% of EBV+/MSI GC patients may not respond to immune checkpoint inhibitor therapy. Thus, our HIERÏ gene signature demonstrates the potential to subclassify tumor immunity levels, predict prognosis and help immunotherapeutic decisions.

A new method to study the change of miRNA-mRNA interactions due to environmental exposures

MOTIVATION

Integrative approaches characterizing the interactions among different types of biological molecules have been demonstrated to be useful for revealing informative biological mechanisms. One such example is the interaction between microRNA (miRNA) and messenger RNA (mRNA), whose deregulation may be sensitive to environmental insult leading to altered phenotypes. The goal of this work is to develop an effective data integration method to characterize deregulation between miRNA and mRNA due to environmental toxicant exposures. We will use data from an animal experiment designed to investigate the effect of low-dose environmental chemical exposure on normal mammary gland development in rats to motivate and evaluate the proposed method.

RESULTS

We propose a new network approach-integrative Joint Random Forest (iJRF), which characterizes the regulatory system between miRNAs and mRNAs using a network model. iJRF is designed to work under the high-dimension low-sample-size regime, and can borrow information across different treatment conditions to achieve more accurate network inference. It also effectively takes into account prior information of miRNA-mRNA regulatory relationships from existing databases. When iJRF is applied to the data from the environmental chemical exposure study, we detected a few important miRNAs that regulated a large number of mRNAs in the control group but not in the exposed groups, suggesting the disruption of miRNA activity due to chemical exposure. Effects of chemical exposure on two affected miRNAs were further validated using breast cancer human cell lines.

AVAILABILITY AND IMPLEMENTATION

R package iJRF is available at CRAN.

CONTACTS

pei.wang@mssm.edu or susan.teitelbaum@mssm.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Molecular Effects of Stromal-Selective Targeting by uPAR-Retargeted Oncolytic Virus in Breast Cancer

The tumor microenvironment (TME) is a relevant target for novel biological therapies. MV-m-uPA and MV-h-uPA are fully retargeted, species-specific, oncolytic measles viruses (MV) directed against murine or human urokinase receptor (PLAUR/uPAR), expressed in tumor and stromal cells. The effects of stromal-selective targeting by uPAR-retargeted MVs were investigated. In vitro infection, virus-induced GFP expression, and cytotoxicity by MV-h-uPA and MV-m-uPA were demonstrated in human and murine cancer cells and cancer-associated fibroblasts in a species-specific manner. In a murine fibroblast/human breast cancer 3D coculture model, selective fibroblast targeting by MV-m-uPA inhibited breast cancer cell growth. Systemic administration of murine-specific MV-m-uPA in mice bearing human MDA-MB-231 xenografts was associated with a significant delay in tumor progression and improved survival compared with controls. Experiments comparing tumor (MV-h-uPA) versus stromal (MV-m-uPA) versus combined virus targeting showed that tumor and stromal targeting was associated with improved tumor control over the other groups. Correlative studies confirmed in vivo viral targeting of tumor stroma by MV-m-uPA, increased apoptosis, and virus-induced differential regulation of murine stromal genes associated with inflammatory, angiogenesis, and survival pathways, as well as indirect regulation of human cancer pathways, indicating viral-induced modulation of tumor-stroma interactions. These data demonstrate the feasibility of stromal-selective targeting by an oncolytic MV, virus-induced modulation of tumor-stroma pathways, and subsequent tumor growth delay. These findings further validate the critical role of stromal uPAR in cancer progression and the potential of oncolytic viruses as antistromal agents.Implications: The current report demonstrates for the first time the biological, in vitro, and in vivo antitumor and molecular effects of stromal selective targeting by an oncolytic virus. Mol Cancer Res; 15(10); 1410-20. ©2017 AACR.

Lung cancer samples preserved in liquid medium: One step beyond cytology

Lung cancer is one of the most common cancer types in men and women worldwide with a high mortality rate. World Health Organization (WHO) classification has accepted biopsy as the primary sample for lung cancer diagnosis, pathological classification and molecular testing for management of patients, yet, the use of alternative sampling procedures is highly encouraged. Bronchial cytological samples require a less invasive collection technique and may be suitable for pathological and molecular analysis and storage in liquid medium. Furthermore, the molecular analysis of bronchial cytological samples allows the detection of molecular biomarkers, which may be useful for the selection of molecular targeted therapies. Thus, the purpose of this review is to describe the usefulness of bronchial cytological samples preserved in liquid medium from lung cancer patients for pathological diagnosis and molecular investigation.

A novel immunotherapy targeting MMP-14 limits hypoxia, immune suppression and metastasis in triple-negative breast cancer models

Matrix metalloproteinase-14 (MMP-14) is a clinically relevant target in metastatic cancers due to its role in tumor progression and metastasis. Since active MMP-14 is localized on the cell surface, it is amenable to antibody-mediated blockade in cancer, and here we describe our efforts to develop novel inhibitory anti-MMP-14 antibodies. A phage-displayed synthetic humanized Fab library was screened against the extracellular domain of MMP-14 and a panel of MMP14-specific Fabs were identified. A lead antibody that inhibits the catalytic domain of MMP-14 (Fab 3369) was identified and treatment of MDA-MB-231 breast cancer cells with Fab 3369 led to significant loss of extracellular matrix degradation and cell invasion abilities. In mammary orthotopic tumor xenograft assays, MMP-14 blockade by IgG 3369 limited tumor growth and metastasis. Analysis of tumor tissue sections revealed that MMP-14 blockade limited tumor neoangiogenesis and hypoxia. Similar effects of MMP-14 blockade in syngeneic 4T1 mammary tumors were observed, along with increased detection of cytotoxic immune cell markers. In conclusion, we show that immunotherapies targeting MMP-14 can limit immune suppression, tumor progression, and metastasis in triple-negative breast cancer.

miRNA regulation is important for DNA damage repair and recognition in malignant pleural mesothelioma

Platin-containing regimes are currently considered as state-of-the-art therapies in malignant pleural mesotheliomas (MPM) but show dissatisfying response rates ranging from 6 to 16% only. Still, the reasons for the rather poor efficacy remain largely unknown. A clear stratification of patients based on new biomarkers seems to be a promising approach to enhance clinical management, which would be a long-needed improvement for MPM patients but does not seem likely soon unless new biomarkers can be validated. Twenty-four formalin-fixed, paraffin-embedded (FFPE) tumour specimens were subjected to a miRNA expression screening of 800 important miRNAs using digital quantification via the nCounter technique (NanoString). We defined a small subset of miRNAs regulating the key enzymes involved in the repair of platin-associated DNA damage. Particularly, the TP53 pathway network for DNA damage recognition as well as genes related to the term "BRCAness" are the main miRNA targets within this context. The TP53 pathway network for DNA damage recognition as well as genes related to the term "BRCAness" are the main players for risk stratification in patients suffering from this severe disease. Taking the specific molecular profile of the tumour into account can help to enhance the clinical management prospectively and to smooth the way to better response prediction.

miRNA assays in the clinical laboratory: workflow, detection technologies and automation aspects

microRNAs (miRNAs) are short non-coding RNA molecules that regulate gene expression in eukaryotes. Their differential abundance is indicative or even causative for a variety of pathological processes including cancer or cardiovascular disorders. Due to their important biological function, miRNAs represent a promising class of novel biomarkers that may be used to diagnose life-threatening diseases, and to monitor disease progression. Further, they may guide treatment selection or dosage of drugs. miRNAs from blood or derived fractions are particularly interesting candidates for routine laboratory applications, as they can be measured in most clinical laboratories already today. This assures a good accessibility of respective tests. Albeit their great potential, miRNA-based diagnostic tests have not made their way yet into the clinical routine, and hence no standardized workflows have been established to measure miRNAs for patients' benefit. In this review we summarize the detection technologies and workflow options that exist to measure miRNAs, and we describe the advantages and disadvantages of each of these options. Moreover, we also provide a perspective on data analysis aspects that are vital for translation of raw data into actionable diagnostic test results.

microRNAs are differentially regulated between MDM2-positive and negative malignant pleural mesothelioma

Background

Malignant pleural mesothelioma (MPM) is a highly aggressive tumour first-line treated with a combination of cisplatin and pemetrexed. MDM2 and P14/ARF (CDKN2A) are upstream regulators of TP53 and may contribute to its inactivation. In the present study, we now aimed to define the impact of miRNA expression on this mechanism.

Material and Methods

24 formalin-fixed paraffin-embedded (FFPE) tumour specimens were used for miRNA expression analysis of the 800 most important miRNAs using the nCounter technique (NanoString). Significantly deregulated miRNAs were identified before a KEGG-pathway analysis was performed.

Results

17 miRNAs regulating TP53, 18 miRNAs regulating MDM2, and 11 miRNAs directly regulating CDKN2A are significantly downregulated in MDM2-expressing mesotheliomas. TP53 is downregulated in MDM2-negative tumours through miRNAs with a miSVR prediction score of 11.67, RB1 with a prediction score of 8.02, MDM2 with a prediction score of 4.50 and CDKN2A with a prediction score of 1.27.

Conclusion

MDM2 expression seems to impact miRNA expression levels in MPM. Especially, miRNAs involved in TP53-signaling are strongly decreased in MDM2-positive mesotheliomas. A better understanding of its tumour biology may open the chance for new therapeutic approaches and thereby augment patients' outcome.

Molecular stratification of early breast cancer identifies drug targets to drive stratified medicine

Many women with hormone receptor-positive early breast cancer can be managed effectively with endocrine therapies alone. However, additional systemic chemotherapy treatment is necessary for others. The clinical challenges in managing high-risk women are to identify existing and novel druggable targets, and to identify those who would benefit from these therapies. Therefore, we performed mRNA abundance analysis using the Tamoxifen and Exemestane Adjuvant Multinational (TEAM) trial pathology cohort to identify a signature of residual risk following endocrine therapy and pathways that are potentially druggable. A panel of genes compiled from academic and commercial multiparametric tests as well as genes of importance to breast cancer pathogenesis was used to profile 3825 patients. A signature of 95 genes, including nodal status, was validated to stratify endocrine-treated patients into high-risk and low-risk groups based on distant relapse-free survival (DRFS; Hazard Ratio = 5.05, 95% CI 3.53-7.22, p = 7.51 × 10-19). This risk signature was also found to perform better than current multiparametric tests. When the 95-gene prognostic signature was applied to all patients in the validation cohort, including patients who received adjuvant chemotherapy, the signature remained prognostic (HR = 4.76, 95% CI 3.61-6.28, p = 2.53× 10-28). Functional gene interaction analyses identified six significant modules representing pathways involved in cell cycle control, mitosis and receptor tyrosine signaling; containing a number of genes with existing targeted therapies for use in breast or other malignancies. Thus the identification of high-risk patients using this prognostic signature has the potential to also prioritize patients for treatment with these targeted therapies.

Identification of biomarkers of intrahepatic cholangiocarcinoma via integrated analysis of mRNA and miRNA microarray data

The present study aimed to identify potential therapeutic targets of intrahepatic cholangiocarcinoma (ICC) via integrated analysis of gene (transcript version) and microRNA (miRNA/miR) expression. The miRNA microarray dataset GSE32957 contained miRNA expression data from 16 ICC, 7 mixed type of combined hepatocellular-cholangiocarcinoma (CHC), 2 hepatic adenoma, 3 focal nodular hyperplasia (FNH) and 5 healthy liver tissue samples, and 2 cholangiocarcinoma cell lines. In addition, the mRNA microarray dataset GSE32879 contained mRNA expression data from 16 ICC, 7 CHC, 2 hepatic adenoma, 5 FNH and 7 healthy liver tissue samples. The datasets were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) and miRNAs (DEMs) in ICC samples compared with healthy liver tissues were identified via the limma package, following data preprocessing. Genes that exhibited alternative splicing (AS) in ICC samples were identified via AltAnalyze software. Functional enrichment analysis of DEGs was performed using the Database for Annotation, Visualization and Integrated Analysis. Target genes of DEMs were identified using the TargetScan database. The regulatory association between DEMs and any overlaps among DEGs, alternative splicing genes (ASGs) and target genes of DEMs were retrieved, and a network was visualized using the Cytoscape software. A total of 2,327 DEGs, 70 DEMs and 623 ASGs were obtained. Functional enrichment analysis indicated that DEGs were primarily enriched in biological processes and pathways associated with cell activity or the immune system. A total of 63 overlaps were obtained among DEGs, ASGs and target genes of DEMs, and a regulation network that contained 243 miRNA-gene regulation pairs was constructed between these overlaps and DEMs. The overlapped genes, including sprouty-related EVH1 domain containing 1, protein phosphate 1 regulatory subunit 12A, chromosome 20 open reading frame 194, and DEMs, including hsa-miR-96, hsa-miR-1 and hsa-miR-25, may be potential therapeutic targets for the future treatment of ICC.