Adaptively Integrative Association between Multivariate Phenotypes and Transcriptomic Data for Complex Diseases.

Phenotype-gene association studies can uncover disease mechanisms for translational research. Association with multiple phenotypes or clinical variables in complex diseases has the advantage of increasing statistical power and offering a holistic view. Existing multi-variate association methods mostly focus on SNP-based genetic associations. In this paper, we extend and evaluate two adaptive Fisher's methods, namely AFp and AFz, from the p-value combination perspective for phenotype-mRNA association analysis. The proposed method effectively aggregates heterogeneous phenotype-gene effects, allows association with different data types of phenotypes, and performs the selection of the associated phenotypes. Variability indices of the phenotype-gene effect selection are calculated by bootstrap analysis, and the resulting co-membership matrix identifies gene modules clustered by phenotype-gene effect. Extensive simulations demonstrate the superior performance of AFp compared to existing methods in terms of type I error control, statistical power and biological interpretation. Finally, the method is separately applied to three sets of transcriptomic and clinical datasets from lung disease, breast cancer, and brain aging and generates intriguing biological findings.

Mixed invasive ductal lobular carcinoma is clinically and pathologically more similar to invasive lobular than ductal carcinoma.

BACKGROUND: Mixed invasive ductal lobular carcinoma (mDLC) remains a poorly understood subtype of breast cancer composed of coexisting ductal and lobular components. METHODS: We sought to describe clinicopathologic characteristics and determine whether mDLC is clinically more similar to invasive ductal carcinoma (IDC) or invasive lobular carcinoma (ILC), using data from patients seen at the University of Pittsburgh Medical Center. RESULTS: We observed a higher concordance in clinicopathologic characteristics between mDLC and ILC, compared to IDC. There is a trend for higher rates of successful breast-conserving surgery after neoadjuvant chemotherapy in patients with mDLC compared to patients with ILC, in which it is known to be lower than in those with IDC. Metastatic patterns of mDLC demonstrate a propensity to develop in sites characteristic of both IDC and ILC. A meta-analysis evaluating mDLC showed shared features with both ILC and IDC with significantly more ER-positive and fewer high grades in mDLC compared to IDC, although mDLCs were significantly smaller and included fewer late-stage tumours compared to ILC. CONCLUSIONS: These findings support clinicopathologic characteristics of mDLC driven by individual ductal vs lobular components and given the dominance of lobular pathology, mDLC features are often more similar to ILC than IDC. This study exemplifies the complexity of mixed disease.

Longitudinal Modulation of Loco-Regional Immunity in Ovarian Cancer Patients Receiving Intraperitoneal Chemotherapy.

The immune tumor microenvironment (TME) of epithelial ovarian cancer (EOC) carries both effector and suppressive functions. To define immune correlates of chemotherapy-induced tumor involution, we performed longitudinal evaluation of biomarker expression on serial biological specimens collected during intraperitoneal (IP) platinum-based chemotherapy. Serial biological samples were collected at several time points during IP chemotherapy. RNA from IP fluid cells and tumor tissue was analyzed via NanoString. Meso Scale Discovery (MSD) multiplex assay and ELISA for MUC1 antibodies were performed on plasma and IP fluid. Differentially expressed genes in IP fluid demonstrate an upregulation of B cell function and activation of Th2 immune response along with dampening of Th1 immunity during chemotherapy. MSD analysis of IP fluid and gene expression analysis of tumor tissue revealed activation of Th2 immunity and the complement system. Anti-MUC1 antibodies were detected in IP fluid samples. IP fluid analysis in a secondary cohort also identified chemotherapy-induced B cell function genes. This study shows that serial IP fluid sampling is an effective method to capture changes in the immune TME during chemotherapy and reveals treatment induced changes in B cell function and Th2 immunity.

ESR1 mutant breast cancers show elevated basal cytokeratins and immune activation.

Estrogen receptor alpha (ER/ESR1) is frequently mutated in endocrine resistant ER-positive (ER+) breast cancer and linked to ligand-independent growth and metastasis. Despite the distinct clinical features of ESR1 mutations, their role in intrinsic subtype switching remains largely unknown. Here we find that ESR1 mutant cells and clinical samples show a significant enrichment of basal subtype markers, and six basal cytokeratins (BCKs) are the most enriched genes. Induction of BCKs is independent of ER binding and instead associated with chromatin reprogramming centered around a progesterone receptor-orchestrated insulated neighborhood. BCK-high ER+ primary breast tumors exhibit a number of enriched immune pathways, shared with ESR1 mutant tumors. S100A8 and S100A9 are among the most induced immune mediators and involve in tumor-stroma paracrine crosstalk inferred by single-cell RNA-seq from metastatic tumors. Collectively, these observations demonstrate that ESR1 mutant tumors gain basal features associated with increased immune activation, encouraging additional studies of immune therapeutic vulnerabilities.

Phase I Trial Combining Chemokine-Targeting with Loco-Regional Chemoimmunotherapy for Recurrent, Platinum-Sensitive Ovarian Cancer Shows Induction of CXCR3 Ligands and Markers of Type 1 Immunity.

PURPOSE: Increased prevalence of cytotoxic T lymphocytes (CTL) in the tumor microenvironment (TME) predicts positive outcomes in patients with epithelial ovarian cancer (EOC), whereas the regulatory T cells (Treg) predict poor outcomes. Guided by the synergistic activity of TLR3 ligands, IFNα, and COX-2 blockers in selectively enhancing CTL-attractants but suppressing Treg-attractants, we tested a novel intraperitoneal chemoimmunotherapy combination (CITC), to assess its tolerability and TME-modulatory impact in patients with recurrent EOC. PATIENTS AND METHODS: Twelve patients were enrolled in phase I portion of the trial NCT02432378, and treated with intraperitoneal cisplatin, intraperitoneal rintatolimod (dsRNA, TLR3 ligand), and oral celecoxib (COX-2 blocker). Patients in cohorts 2, 3, and 4 also received intraperitoneal IFNα at 2, 6, and 18 million units (MU), respectively. Primary objectives were to evaluate safety, identify phase 2 recommended dose (P2RD), and characterize changes in the immune TME. Peritoneal resident cells and intraperitoneal wash fluid were profiled via NanoString and Meso Scale Discovery (MSD) multiplex assay, respectively. RESULTS: The P2RD of IFNα was 6 MU. Median progression-free survival and overall survival were 8.4 and 30 months, respectively. Longitudinal sampling of the peritoneal cavity via intraperitoneal washes demonstrated local upregulation of IFN-stimulated genes (ISG), including CTL-attracting chemokines (CXCL-9, -10, -11), MHC I/II, perforin, and granzymes. These changes were present 2 days after chemokine modulation and subsided within 1 week. CONCLUSIONS: The chemokine-modulating intraperitoneal-CITC is safe, tolerable, and associated with ISG changes that favor CTL chemoattraction and function. This combination (plus DC vaccine) will be tested in a phase II trial. See related commentary by Aranda et al., p. 1993.

Outcomes After Sentinel Lymph Node Biopsy and Radiotherapy in Older Women With Early-Stage, Estrogen Receptor-Positive Breast Cancer.

Importance: Overtreatment of early-stage breast cancer with favorable tumor biology in older patients may be harmful without affecting recurrence and survival. Guidelines that recommend deimplementation of sentinel lymph node biopsy (SLNB) (Choosing Wisely) and radiotherapy (RT) (National Comprehensive Cancer Network) have been published. Objective: To describe the use rates and association with disease recurrence of SLNB and RT in older women with breast cancer. Design, Setting, and Participants: This cohort study obtained patient and clinical data from an integrated cancer registry and electronic health record of a single health care system in Pennsylvania. The cohort was composed of consecutive female patients 70 years or older who were diagnosed with early-stage, estrogen receptor-positive, ERBB2 (formerly HER2)-negative, clinically node-negative breast cancer from January 1, 2010, to December 31, 2018, who were treated at 15 community and academic hospitals within the health system. Exposures: Sentinel lymph node biopsy and adjuvant RT. Main Outcomes and Measures: Primary outcomes were 5-year locoregional recurrence-free survival (LRFS) rate and disease-free survival (DFS) rate after SLNB and after RT. Secondary outcomes included recurrence rate, subgroups that may benefit from SLNB or RT, and use rate of SLNB and RT over time. Propensity scores were used to create 2 cohorts to separately evaluate the association of SLNB and RT with recurrence outcomes. Cox proportional hazards regression model was used to estimate hazard ratios (HRs). Results: From 2010 to 2018, a total of 3361 women 70 years or older (median [interquartile range {IQR}] age, 77.0 [73.0-82.0] years) with estrogen receptor-positive, ERBB2-negative, clinically node-negative breast cancer were included in the study. Of these women, 2195 (65.3%) received SLNB and 1828 (54.4%) received adjuvant RT. Rates of SLNB steadily increased (1.0% per year), a trend that persisted after the 2016 adoption of the Choosing Wisely guideline. Rates of RT decreased slightly (3.4% per year). To examine patient outcomes and maximize follow-up time, the analysis was limited to cases from 2010 to 2014, identifying 2109 patients with a median (IQR) follow-up time of 4.1 (2.5-5.7) years. In the propensity score-matched cohorts, no association was found between SLNB and either LRFS (HR, 1.26; 95% CI, 0.37-4.30; P = .71) or DFS (HR, 1.92; 95% CI, 0.86-4.32; P = .11). In addition, RT was not associated with LRFS (HR, 0.33; 95% CI, 0.09-1.24; P = .10) or DFS (HR, 0.99; 95% CI, 0.46-2.10; P = .97). Subgroup analysis showed that stratification by tumor grade or comorbidity was not associated with LRFS or DFS. Low absolute rates of recurrence were observed when comparing the groups that received SLNB (3.5%) and those that did not (4.5%) as well as the groups that received RT (2.7%) and those that did not (5.5%). Conclusions and Relevance: This study found that receipt of SLNB or RT was not associated with improved LRFS or DFS in older patients with ER-positive, clinically node-negative breast cancer. Despite limited follow-up time and wide 95% CIs, this study supports the continued deimplementation of both SLNB and RT in accordance with the Choosing Wisely and National Comprehensive Cancer Network guidelines.

Older molecular brain age in severe mental illness.

Psychiatric disorders are associated with accelerated aging and enhanced risk for neurodegenerative disorders. Brain aging is associated with molecular, cellular, and structural changes that are robust on the group level, yet show substantial inter-individual variability. Here we assessed deviations in gene expression from normal age-dependent trajectories, and tested their validity as predictors of risk for major mental illnesses and neurodegenerative disorders. We performed large-scale gene expression and genotype analyses in postmortem samples of two frontal cortical brain regions from 214 control subjects aged 20-90 years. Individual estimates of "molecular age" were derived from age-dependent genes, identified by robust regression analysis. Deviation from chronological age was defined as "delta age". Genetic variants associated with deviations from normal gene expression patterns were identified by expression quantitative trait loci (cis-eQTL) of age-dependent genes or genome-wide association study (GWAS) on delta age, combined into distinct polygenic risk scores (PRScis-eQTL and PRSGWAS), and tested for predicting brain disorders or pathology in independent postmortem expression datasets and clinical cohorts. In these validation datasets, molecular ages, defined by 68 and 76 age-related genes for two brain regions respectively, were positively correlated with chronological ages (r = 0.88/0.91), elevated in bipolar disorder (BP) and schizophrenia (SCZ), and unchanged in major depressive disorder (MDD). Exploratory analyses in independent clinical datasets show that PRSs were associated with SCZ and MDD diagnostics, and with cognition in SCZ and pathology in Alzheimer's disease (AD). These results suggest that older molecular brain aging is a common feature of severe mental illnesses and neurodegeneration.

Recent Advances in Computer-Assisted Algorithms for Cell Subtype Identification of Cytometry Data.

The progress in the field of high-dimensional cytometry has greatly increased the number of markers that can be simultaneously analyzed producing datasets with large numbers of parameters. Traditional biaxial manual gating might not be optimal for such datasets. To overcome this, a large number of automated tools have been developed to aid with cellular clustering of multi-dimensional datasets. Here were review two large categories of such tools; unsupervised and supervised clustering tools. After a thorough review of the popularity and use of each of the available unsupervised clustering tools, we focus on the top six tools to discuss their advantages and limitations. Furthermore, we employ a publicly available dataset to directly compare the usability, speed, and relative effectiveness of the available unsupervised and supervised tools. Finally, we discuss the current challenges for existing methods and future direction for the new generation of cell type identification approaches.

Bayesian Detection of Abnormal Asynchrony of Division Between Sister Cells in Mutant Caenorhabditis elegans Embryos.

Cell division timing is critical for cell fate specification and morphogenesis during embryogenesis, but how division timings are regulated among cells during development is poorly understood. In this article, we focus on the comparison of asynchrony of division, that is, difference of lifetime, between sister cells (ADS) among wild-type and mutant individuals of Caenorhabditis elegans. On the one hand, due to extreme imbalance between wild-type individuals and mutant-type samples, direct comparison of two distributions of ADS between wild type and mutant type is not feasible. On the other hand, we originally found that the ADS is correlated with the lifespan of the corresponding mother cell in wild type. Hence, a semiparametric Bayesian quantile regression method where lifetime of the mother cell is taken as covariate is developed to estimate the 95% confidence curve of ADS in wild type and then ADS of mutant type is classified as abnormal if outside the corresponding confidence interval. A high accuracy of our method is demonstrated by a large-scale simulation study. Real data analysis shows that ADS is related to gene function and expression quantitatively.

Cisplatin-induced immune modulation in ovarian cancer mouse models with distinct inflammation profiles.

The backbone of ovarian cancer treatment is platinum-based chemotherapy and aggressive surgical debulking. New therapeutic approaches using immunotherapy via immune checkpoint blockade, which have demonstrated clinical efficacy in other tumor types, have been less promising in ovarian cancer. To increase their clinical efficacy, checkpoint inhibitors are now being tested in clinical trials in combination with chemotherapy. Here, we evaluated the impact of cisplatin on tumor immunogenicity and its in vivo roles when used alone or in combination with anti-PD-L1, in two novel murine ovarian cancer cell models. The 2F8 and its platinum-resistant derivative 2F8cis model, display distinct inflammatory profiles and chemotherapy sensitivities, and mirror the primary and recurrent human disease, respectively. Acute and chronic exposure to cisplatin enhances tumor immunogenicity by increasing calreticulin, MHC class I, antigen presentation and T-cell infiltration. Cisplatin also upregulates PD-L1 expression in vitro and in vivo, demonstrating a dual, paradoxical immune modulatory effect and supporting the rationale for combination with immune checkpoint blockade. One of the pathways activated by cisplatin treatment is the cGAS/STING pathway. Chronic cisplatin treatment led to upregulation of cGAS and STING proteins in 2F8cis compared to parental 2F8 cells, while acute exposure to cisplatin further increases cGAS and STING levels in both 2F8 and 2F8cis cells. Overexpression of cGAS/STING modifies tumor immunogenicity by upregulating PD-L1, MHC I and calreticulin in tumor cells. Anti-PD-L1 alone in a platinum-sensitive model or with cisplatin in a platinum-resistant model increases survival. These studies have high translational potential in ovarian cancer.