Search

Background Diurnal rhythms in gene expression have been detected in the human brain. Previous studies found that males and females exhibit 24-hour rhythms in known circadian genes, with earlier peak expression in females. Whether there are sex differences in large-scale transcriptional rhythms in the cortex that align with observed sex differences in physiological and behavioral rhythms is currently unknown. Methods Diurnal rhythmicity of gene expression was determined for males and females using RNA sequencing data from human postmortem dorsolateral prefrontal cortex (DLPFC) and anterior cingulate cortex (ACC). Sex differences among rhythmic genes were determined using significance cutoffs, threshold-free analyses, and R2 difference. Phase concordance was assessed across the DLPFC and ACC for males and females. Pathway and transcription factor analyses were also conducted on significantly rhythmic genes. Results Canonical circadian genes had diurnal rhythms in both sexes with similar amplitude and phase. When analyses were expanded to the entire transcriptome, significant sex differences in transcriptional rhythms emerged. There were nearly twice as many rhythmic transcripts in the DLPFC in males and nearly 4 times as many rhythmic transcripts in the ACC in females. Results suggest a diurnal rhythm in synaptic transmission specific to the ACC in females (e.g., GABAergic [gamma-aminobutyric acidergic] and cholinergic neurotransmission). For males, there was phase concordance between the DLPFC and ACC, while phase asynchrony was found in females. Conclusions There are robust sex differences in molecular rhythms of genes in the DLPFC and ACC, providing potential mechanistic insights into how neurotransmission and synaptic function are modulated in a circadian-dependent and sex-specific manner.

Background: Estrogen receptor alpha (ER/ESR1) mutations are found in 20-40% of endocrine resistant ER+ metastatic breast cancers, and they are associated with worse outcome. Preclinical studies have shown that they cause ligand-independent growth, resistance to endocrine therapy, and there is growing evidence for a role in metastasis. It is not known how ESR1 mutant cancer cells cause metastases and whether such a mechanism may indicate novel therapeutic options. Methods: ddPCR was used to detect hotspot ESR1 mutations. Transcriptome data were derived from cell line models and clinical samples. For in vitro phenotypic characterization, Y537S and D538G genome-edited MCF7 and T47D cell models from four different labs were used. Cell-cell adhesive properties were assessed using calcein-labelled adhesion, spontaneous cell aggregation and microfluidic aggregation assays. Altered expression of cell-cell adhesion genes detected in RNA seq data was validated via qRT-PCR, immunoblot and immune staining. Functional contributions of desmosome and gap junctions were tested by blocking peptide and carbenoxolone respectively. Mutant ER cistromes were profiled using ChIP-sequencing. Tail vein injection was performed on nude mice to evaluate metastasis in vivo. Mouse lung micro-metastatic foci were quantified using human-specific CK19 staining. Circulating tumor cells (CTCs) clustering propensity in vivo was assessed via intracardiac injection of ESR1 WT and mutant cells into nude mice following CTC microfilter capture. CTCs enumeration from breast cancer patients’ blood samples was performed using CellSearchTM system. Results: We identified a significant enrichment of ESR1 mutations in distant (12/48) vs local (0/27) recurrences, confirming the strong association of mutant ER with metastasis. Transcriptomic analysis revealed altered cell-to-cell interaction pathways in ESR1 mutant tumors compared to ESR1 WT tumors, suggesting a previously undescribed role of ESR1 mutations in reprogramming cell-cell adhesome. ESR1 mutant cells grown in suspension culture revealed more compact multicellular spheroids compared to WT cells. This observation was confirmed under both static and microfluidic conditions, indicative of increased cell-cell interactions. The effect was more pronounced in MCF7 compared to T47D cells, and it correlated with increased expression of multiple desmosome and gap junction pathway genes, which were also significantly enriched in ESR1 mutant tumors. Pharmacological blockade of desmosome and gap junctions significantly rescued enhanced cell-cell adhesion in ESR1 mutant cells. Mechanistically, our ER ChIP-seq did not identify any gained mutant ER binding sites in proximity to cell-cell adhesion gene loci, indicating indirect regulation by mutant ER. Consistent with this, expression of Connexin 43, one of the top upregulated gap junction components, was induced by cFOS found to be highly upregulated in ESR1 mutant cells. Tail vein injection of ESR1-mutant cells derived more distant macro- (MCF7) and micro- (T47D) metastases. Given increasing evidence for role of cell-cell attachment in CTC phenotypes, we tested CTC formation for ESR1 WT and mutant cells. In vivo studies showed MCF7 Y537S ESR1 mutant cells formed larger multi-cellular CTC clusters with increased compactness compared to WT CTCs. These preclinical data translated to clinical observation, where we observed an enrichment of CTC clusters in patients with ESR1 mutant-metastatic breast cancers. Conclusion: Hotspot ESR1 mutations induce expression of multiple desmosome and gap junction genes and confer increased cell-cell adhesion, which facilitate breast cancer metastasis via increased CTCs clustering propensity. These findings might guide approaches to test potential repurpose of drugs targeting gap junction in ER mutant tumors. Citation Format: Zheqi Li, Yang Wu, Amir Bahreini, Jian Chen, Ye Qin, Kevin M. Levine, Nilgun Tasdemir, Nolan Priedigkeit, Li Zhu, George C. Tseng, Yu Jiang, Benjamin Troness, Laki Buluwela, Simak Ali, Spencer Arnesen, Jason Gertz, Ben H. Park, Qiang Zhang, Lorenzo Gerratana, Andrew Davis, Jennifer M. Atkinson, Dorraya El-Ashry, Massimo Cristofanilli, Adrian V. Lee, Steffi Oesterreich. Hotspot ESR1 mutations rewire cell-cell adhesome to facilitate breast cancer metastasis [abstract]. In: Proceedings of the 2020 San Antonio Breast Cancer Virtual Symposium; 2020 Dec 8-11; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2021;81(4 Suppl):Abstract nr GS1-05.

Inactivation of Pten gene through deletions and mutations leading to excessive pro-growth signaling pathway activations frequently occurs in cancers. Here, we report a Pten derived pro-cancer growth gene fusion Pten-NOLC1 originated from a chr10 genome rearrangement and identified through a transcriptome sequencing analysis of human cancers. Pten-NOLC1 fusion is present in primary human cancer samples and cancer cell lines from different organs. The product of Pten-NOLC1 is a nuclear protein that interacts and activates promoters of EGFR, c-MET, and their signaling molecules. Pten-NOLC1 promotes cancer proliferation, growth, invasion, and metastasis, and reduces the survival of animals xenografted with Pten-NOLC1-expressing cancer cells. Genomic disruption of Pten-NOLC1 induces cancer cell death, while genomic integration of this fusion gene into the liver coupled with somatic Pten deletion produces spontaneous liver cancers in mice. Our studies indicate that Pten-NOLC1 gene fusion is a driver for human cancers.

Women are approximately two times as likely to be diagnosed with major depressive disorder (MDD) compared to men. While sex differences in MDD might be driven by circulating gonadal hormones, we hypothesized that developmental hormone exposure and/or genetic sex might play a role. Mice were gonadectomized in adulthood to isolate the role of developmental hormones. We examined the effects of developmental gonadal and genetic sex on anhedonia-/depressive-like behaviors under non-stress and chronic stress conditions and performed RNA-sequencing in three mood-relevant brain regions. We used an integrative network approach to identify transcriptional modules and stress-specific hub genes regulating stress susceptibility, with a focus on whether these differed by sex. After identifying sex differences in anhedonia-/depressive-like behaviors (female > male), we show that both developmental hormone exposure (gonadal female > gonadal male) and genetic sex (XX > XY) contribute to the sex difference. The top biological pathways represented by differentially expressed genes were related to immune function; we identify which differentially expressed genes are driven by developmental gonadal or genetic sex. There was very little overlap in genes affected by chronic stress in males and females. We also identified highly co-expressed gene modules affected by stress, some of which were affected in opposite directions in males and females. Since all mice had equivalent hormone exposure in adulthood, these results suggest that sex differences in gonadal hormone exposure during sensitive developmental periods program adult sex differences in mood, and that these sex differences are independent of adult circulating gonadal hormones.

Maintenance of healthy pregnancy is reliant on successful balance between the fetal and maternal immune systems. Although maternal mechanisms responsible have been well studied, those used by the fetal immune system remain poorly understood. Using suspension mass cytometry and various imaging modalities, we report a complex immune system within the mid-gestation (17-23 weeks) human placental villi (PV). Further, we identified immunosuppressive signatures in innate immune cells and antigen presenting cells that potentially maintain immune homeostasis in utero. Consistent with recent reports in other fetal organs, T cells with memory phenotypes were detected within the PV tissue and vasculature. Moreover, we determined PV T cells could be activated to upregulate CD69 and proliferate after T cell receptor (TCR) stimulation and when exposed to maternal uterine antigens. Finally, we report that cytokine production by PV T cells is sensitive to TCR stimulation and varies between mid-gestation, preterm (26-35 weeks) and term deliveries (37-40 weeks). Collectively, we elucidated the complexity and functional maturity of fetal immune cells within the PV and highlighted their immunosuppressive potential.

Background Prevalence rates of opioid use disorder (OUD) have increased dramatically, accompanied by a surge of overdose deaths. While opioid dependence has been extensively studied in preclinical models, an understanding of the biological alterations that occur in the brains of people who chronically use opioids and who are diagnosed with OUD remains limited. To address this limitation, RNA sequencing was conducted on the dorsolateral prefrontal cortex and nucleus accumbens, regions heavily implicated in OUD, from postmortem brains in subjects with OUD. Methods We performed RNA sequencing on the dorsolateral prefrontal cortex and nucleus accumbens from unaffected comparison subjects (n = 20) and subjects diagnosed with OUD (n = 20). Our transcriptomic analyses identified differentially expressed transcripts and investigated the transcriptional coherence between brain regions using rank-rank hypergeometric orderlap. Weighted gene coexpression analyses identified OUD-specific modules and gene networks. Integrative analyses between differentially expressed transcripts and genome-wide association study datasets using linkage disequilibrium scores assessed the genetic liability of psychiatric-related phenotypes in OUD. Results Rank-rank hypergeometric overlap analyses revealed extensive overlap in transcripts between the dorsolateral prefrontal cortex and nucleus accumbens in OUD, related to synaptic remodeling and neuroinflammation. Identified transcripts were enriched for factors that control proinflammatory cytokine, chondroitin sulfate, and extracellular matrix signaling. Cell-type deconvolution implicated a role for microglia as a potential driver for opioid-induced neuroplasticity. Linkage disequilibrium score analysis suggested genetic liabilities for risky behavior, attention-deficit/hyperactivity disorder, and depression in subjects with OUD. Conclusions Overall, our findings suggest connections between the brain’s immune system and opioid dependence in the human brain.

Circadian rhythmicity in transcriptomic profiles has been shown in many physiological processes, and the disruption of circadian patterns has been founded to associate with several diseases. In this paper, we developed a series of likelihood-based methods to detect (i) circadian rhythmicity (denoted as LR rhythmicity) and (ii) differential circadian patterns comparing two experimental conditions (denoted as LR diff). In terms of circadian rhythmicity detection, we demonstrated that our proposed LR rhythmicity could better control the type I error rate compared to existing methods under a wide variety of simulation settings. In terms of differential circadian patterns, we developed methods in detecting differential amplitude, differential phase, differential basal level, and differential fit, which also successfully controlled the type I error rate. In addition, we demonstrated that the proposed LR diff could achieve higher statistical power in detecting differential fit, compared to existing methods. The superior performance of LR rhythmicity and LR diff was demonstrated in two real data applications, including a brain aging data (gene expression microarray data of human postmortem brain) and a time-restricted feeding data (RNA sequencing data of human skeletal muscles). An R package for our methods is publicly available on GitHub https://github.com/diffCircadian/diffCircadian.

Constitutively active estrogen receptor-α (ER/ESR1) mutations have been identified in approximately one third of ER+ metastatic breast cancer. Although these mutations are known mediators of endocrine resistance, their potential role in promoting metastatic disease has not yet been mechanistically addressed. In this study, we show the presence of ESR1 mutations exclusively in distant, but not local recurrences. In concordance with transcriptomic profiling of ESR1 mutant tumors, genome-edited Y537S and D538G cell models have a reprogrammed cell adhesive gene network via alterations in desmosome/gap junction genes and the TIMP3/MMP axis, which functionally confers enhanced cell-cell contacts while decreased cell-ECM adhesion. Context-dependent migratory phenotypes revealed co-targeting of Wnt and ER as vulnerability. Mutant ESR1 exhibits non-canonical regulation of several metastatic pathways including secondary transactivation and de novo FOXA1-driven chromatin remodeling. Collectively, our data supports evidence for ESR1 mutation-driven metastases and provides insight for future preclinical therapeutic strategies.SignificanceContext and allele-dependent transcriptome and cistrome reprogramming in genome-edited ESR1 mutation cell models elicit diverse metastatic phenotypes, including but not limited to alterations in cell adhesion and migration. The gain-of-function mutations can be pharmacologically targeted, and thus may be key components of novel therapeutic treatment strategies for ER-mutant metastatic breast cancer.

The human striatum can be subdivided into the caudate, putamen, and nucleus accumbens (NAc). Each of these structures have some overlapping and some distinct functions related to motor control, cognitive processing, motivation, and reward. Previously, we used a "time-of-death" approach to identify diurnal rhythms in RNA transcripts in human cortical regions. Here, we identify molecular rhythms across the three striatal subregions collected from postmortem human brain tissue in subjects without psychiatric or neurological disorders. Core circadian clock genes are rhythmic across all three regions and show strong phase concordance across regions. However, the putamen contains a much larger number of significantly rhythmic transcripts than the other two regions. Moreover, there are many differences in pathways that are rhythmic across regions. Strikingly, the top rhythmic transcripts in NAc (but not the other regions) are predominantly small nucleolar RNAs and long noncoding RNAs, suggesting that a completely different mechanism might be used for the regulation of diurnal rhythms in translation and/or RNA processing in the NAc versus the other regions. Further, although the NAc and putamen are generally in phase with regard to timing of expression rhythms, the NAc and caudate, and caudate and putamen, have several clusters of discordant rhythmic transcripts, suggesting a temporal wave of specific cellular processes across the striatum. Taken together, these studies reveal distinct transcriptome rhythms across the human striatum and are an important step in helping to understand the normal function of diurnal rhythms in these regions and how disruption could lead to pathology.

Symbiotic microbial colonization through the establishment of the intestinal microbiome is critical to many intestinal functions, including nutrient metabolism, intestinal barrier integrity, and immune regulation. Recent studies suggest that education of intestinal immunity may be ongoing in utero. However, the drivers of this process are unknown. The microbiome and its byproducts are one potential source. Whether a fetal intestinal microbiome exists is controversial, and whether microbially derived metabolites are present in utero is unknown. Here, we aimed to determine whether bacterial DNA and microbially derived metabolites can be detected in second trimester human intestinal samples. Although we were unable to amplify bacterial DNA from fetal intestines, we report a fetal metabolomic intestinal profile with an abundance of bacterially derived and host-derived metabolites commonly produced in response to microbiota. Though we did not directly assess their source and function, we hypothesize that these microbial-associated metabolites either come from the maternal microbiome and are vertically transmitted to the fetus to prime the fetal immune system and prepare the gastrointestinal tract for postnatal microbial encounters or are produced locally by bacteria that were below our detection threshold., A unique human fetal metabolomic intestinal profile is reported with an abundance of bacterially derived metabolites.

BackgroundPrevalence rates of opioid use disorder (OUD) have increased dramatically, accompanied by a surge of overdose deaths. While opioid dependence has been extensively studied in preclinical models, an understanding of the biological alterations that occur in the brains of people who chronically use opioids and who are diagnosed with OUD remains limited. To address this limitation, RNA-sequencing (RNA-seq) was conducted on the dorsolateral prefrontal cortex (DLPFC) and nucleus accumbens (NAc), regions heavily implicated in OUD, from postmortem brains in subjects with OUD.MethodsWe performed RNA-seq on the DLPFC and NAc from unaffected comparison subjects (n=20) and subjects diagnosed with OUD (n=20). Our transcriptomic analyses identified differentially expressed (DE) transcripts and investigated the transcriptional coherence between brain regions using rank-rank hypergeometric ordering (RRHO). Weighted gene co-expression analyses (WGCNA) also identified OUD-specific modules and gene networks. Integrative analyses between DE transcripts and GWAS datasets using linkage disequilibrium score (LDSC) assessed the genetic liability psychiatric-related phenotypes.ResultsRRHO analyses revealed extensive overlap in transcripts between DLPFC and NAc in OUD, primarily relating to synaptic remodeling and neuroinflammation. Identified transcripts were enriched for factors that control pro-inflammatory cytokine-mediated, chondroitin sulfate, and extracellular matrix signaling. Cell-type deconvolution implicated a role for microglia as a critical driver for opioid-induced neuroplasticity. Using LDSC, we discovered genetic liabilities for risky behavior, attention deficit hyperactivity disorder, and depression.ConclusionsOverall, our findings reveal new connections between the brain’s immune system and opioid dependence in the human brain.

Alternative polyadenylation (APA) in breast tumor samples results in the removal/addition of cis-regulatory elements such as microRNA (miRNA) target sites in the 3′-untranslated region (3′-UTRs) of genes. Although previous computational APA studies focused on a subset of genes strongly affected by APA (APA genes), we identify miRNAs of which widespread APA events collectively increase or decrease the number of target sites [probabilistic inference of microRNA target site modification through APA (PRIMATA-APA)]. Using PRIMATA-APA on the cancer genome atlas (TCGA) breast cancer data, we found that the global APA events change the number of the target sites of particular microRNAs [target sites modified miRNA (tamoMiRNA)] enriched for cancer development and treatments. We also found that when knockdown (KD) of NUDT21 in HeLa cells induces a different set of widespread 3′-UTR shortening than TCGA breast cancer data, it changes the target sites of the common tamoMiRNAs. Since the NUDT21 KD experiment previously demonstrated the tumorigenic role of APA events in a miRNA dependent fashion, this result suggests that the APA-initiated tumorigenesis is attributable to the miRNA target site changes, not the APA events themselves. Further, we found that the miRNA target site changes identify tumor cell proliferation and immune cell infiltration to the tumor microenvironment better than the miRNA expression levels or the APA events themselves. Altogether, our computational analyses provide a proof-of-concept demonstration that the miRNA target site information indicates the effect of global APA events with a potential as predictive biomarker.

Throughout pregnancy, the maternal immunity is tolerant to the developing fetus. However, recent work suggests that the fetal immune system in utero is highly developed and may also contribute to tolerance. We hypothesized that each layer (maternal decidua, fetal membranes, and fetal villi) of the second trimester (ST) placenta harbors a distinct immune landscape. Mass cytometry of 11 ST placentas demonstrated that each tissue did have a unique immune profile. Specifically, the villi contain chemokine receptor (CCR)lo innate cells and enrichment of macrophages and DCs. In contrast, the decidua’s innate cells were CCRhi with abundant NK cells. Random forest classification accurately assigned samples as villi or decidua based on three distinct populations. Memory CD8 T cells were the dominant adaptive cells in all three layers and were more activated in decidua than villi. However, T cells within the villi could secrete TNFα/IFNγ when stimulated with maternal components, suggesting that villi T cells become proinflammatory upon exposure to maternal antigens. Furthermore, T cells from preterm villi had more activation transcripts than ST and term villi. Collectively, these findings illustrate a complex arsenal of immune cells within the placenta in utero and suggest that fetal T cells may contribute to preterm pathology.

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.

Invasive lobular carcinoma (ILC) is the second most common subtype of breast cancer following invasive ductal carcinoma (IDC) and characterized by the loss of E-cadherin–mediated adherens junctions. Despite displaying unique histologic and clinical features, ILC still remains a chronically understudied disease, with limited knowledge gleaned from available laboratory research models. Here we report a comprehensive 2D and 3D phenotypic characterization of four estrogen receptor–positive human ILC cell lines: MDA-MB-134, SUM44, MDA-MB-330, and BCK4. Compared with the IDC cell lines MCF7, T47D, and MDA-MB-231, ultra-low attachment culture conditions revealed remarkable anchorage independence unique to ILC cells, a feature not evident in soft-agar gels. Three-dimensional Collagen I and Matrigel culture indicated a generally loose morphology for ILC cell lines, which exhibited differing preferences for adhesion to extracellular matrix proteins in 2D. Furthermore, ILC cells were limited in their ability to migrate and invade in wound-scratch and transwell assays, with the exception of haptotaxis to Collagen I. Transcriptional comparison of these cell lines confirmed the decreased cell proliferation and E-cadherin–mediated intercellular junctions in ILC while uncovering the induction of novel pathways related to cyclic nucleotide phosphodiesterase activity, ion channels, drug metabolism, and alternative cell adhesion molecules such as N-cadherin, some of which were differentially regulated in ILC versus IDC tumors. Altogether, these studies provide an invaluable resource for the breast cancer research community and facilitate further functional discoveries toward understanding ILC, identifying novel drug targets, and ultimately improving the outcome of patients with ILC. Significance: These findings provide the breast cancer research community with a comprehensive assessment of human invasive lobular carcinoma (ILC) cell line signaling and behavior in various culture conditions, aiding future endeavors to develop therapies and to ultimately improve survival in patients with ILC. Cancer Res; 78(21); 6209–22. ©2018 AACR.

To investigate changes in estrogen receptor alpha (ERα) signaling during progression of endometriosis to endometriosis-associated ovarian cancer (EAOC) as a driver of malignant transformation. We procured tissue samples of normal endometrium, endometriosis (benign, atypical, concurrent with EAOC), and EAOC. We evaluated expression of a 236-gene signature of estrogen signaling. ANOVA and unsupervised clustering were used to identify gene expression profiles across disease states. These profiles were compared to profiles of estrogen regulation in cancer models from the Gene Expression Omnibus (GEO). Gene Set Enrichment Analysis (GSEA) was performed to determine whether gene expression in EAOC was consistent with ERα activity. ANOVA revealed 158 differentially expressed genes (q

Schizophrenia (SZ) is associated with dysfunction of the dorsolateral prefrontal cortex (DLPFC). This dysfunction is manifest as cognitive deficits that appear to arise from disturbances in gamma frequency oscillations. These oscillations are generated in DLPFC layer 3 (L3) via reciprocal connections between pyramidal cells (PCs) and parvalbumin (PV)-containing interneurons. The density of cortical PV neurons is not altered in SZ, but expression levels of several transcripts involved in PV cell function, including PV, are lower in the disease. However, the transcriptome of PV cells has not been comprehensively assessed in a large cohort of subjects with SZ. In this study, we combined an immunohistochemical approach, laser microdissection, and microarray profiling to analyze the transcriptome of DLPFC L3 PV cells in 36 matched pairs of SZ and unaffected comparison subjects. Over 800 transcripts in PV neurons were identified as differentially expressed in SZ subjects; most of these alterations have not previously been reported. The altered transcripts were enriched for pathways involved in mitochondrial function and tight junction signaling. Comparison with the transcriptome of L3 PCs from the same subjects revealed both shared and distinct disease-related effects on gene expression between cell types. Furthermore, network structures of gene pathways differed across cell types and subject groups. These findings provide new insights into cell type-specific molecular alterations in SZ which may point toward novel strategies for identifying therapeutic targets.

The orbitofrontal cortex (OFC) is vulnerable to normal and pathologic aging. Currently, layer resolution large-scale proteomic studies describing “normal” age-related alterations at OFC are not available. Here, we performed a large-scale exploratory high-throughput mass spectrometry–based protein analysis on OFC layer 2/3 from 15 “young” (15–43 years) and 18 “old” (62–88 years) human male subjects. We detected 4193 proteins and identified 127 differentially expressed (DE) proteins (p-value ≤0.05; effect size >20%), including 65 up- and 62 downregulated proteins (e.g., GFAP, CALB1). Using a previously described categorization of biological aging based on somatic tissues, that is, peripheral “hallmarks of aging,” and considering overlap in protein function, we show the highest representation of altered cell-cell communication (54%), deregulated nutrient sensing (39%), and loss of proteostasis (35%) in the set of OFC layer 2/3 DE proteins. DE proteins also showed a significant association with several neurologic disorders; for example, Alzheimer's disease and schizophrenia. Notably, despite age-related changes in individual protein levels, protein co-expression modules were remarkably conserved across age groups, suggesting robust functional homeostasis. Collectively, these results provide biological insight into aging and associated homeostatic mechanisms that maintain normal brain function with advancing age.

As the sequencing cost continued to drop in the past decade, RNA sequencing (RNA-seq) has replaced microarray to become the standard high-throughput experimental tool to analyze transcriptomic profile. As more and more datasets are generated and accumulated in the public domain, meta-analysis to combine multiple transcriptomic studies to increase statistical power has received increasing popularity. In this article, we propose a Bayesian hierarchical model to jointly integrate microarray and RNA-seq studies. Since systematic fold change differences across RNA-seq and microarray for detecting differentially expressed genes have been previously reported, we replicated this finding in several real datasets and showed that incorporation of a normalization procedure to account for the bias improves the detection accuracy and power. We compared our method with the popular two-stage Fisher's method using simulations and two real applications in a histological subtype (invasive lobular carcinoma) of breast cancer comparing PR+ versus PR− and early-stage versus late-stage patients. The result showed improved detection power and more significant and interpretable pathways enriched in the detected biomarkers from the proposed Bayesian model.

Miniature chromosome maintenance (MCM) proteins play critical roles in DNA replication licensing, initiation and elongation. MCM8, one of the MCM proteins playing a critical role in DNA repairing and recombination, was found to have over-expression and increased DNA copy number in a variety of human malignancies. The gain of MCM8 is associated with aggressive clinical features of several human cancers. Increased expression of MCM8 in prostate cancer is associated with cancer recurrence. Forced expression of MCM8 in RWPE1 cells, the immortalized but non-transformed prostate epithelial cell line, exhibited fast cell growth and transformation, while knocked down of MCM8 in PC3, DU145 and LNCaP cells induced cell growth arrest, and decreased tumor volumes and mortality of severe combined immunodeficiency mice xenografted with PC3 and DU145 cells. MCM8 bound cyclin D1 and activated Rb protein phosphorylation by cyclin-dependent kinase 4 in vitro and in vivo. The cyclin D1/MCM8 interaction is required for Rb phosphorylation and S phase entry in cancer cells. As a result, our study showed that copy number increase and overexpression of MCM8 may play critical roles in human cancer development.

Specifically targeting genomic rearrangements and mutations in tumor cells has remained an elusive goal in cancer therapy. Here, we use Cas9-based genome editing to introduce the gene for the pro-drug converting enzyme ‘herpes simplex virus type 1 thymidine kinase’ (HSV1-tk) into the genome of cancer cells that carry unique sequences resulting from genome rearrangements. Specifically, we targeted the breakpoints of TMEM135-CCDC67 and MAN2A1-FER fusions in human prostate cancer or hepatocellular carcinoma cells in vitro and in mouse xenografts. We designed one adenovirus to deliver the nickase Cas9D10A and gRNAs targeting the breakpoint sequences and another to deliver an EGFP-HSV1-tk construct flanked by sequences homologous to those surrounding the breakpoint. Infection with both viruses resulted in breakpoint-dependent expression of EGFP-tk and ganciclovir-mediated apoptosis. When mouse xenografts were treated with adenoviruses and ganciclovir, all animals showed reduction of tumor burden with no mortality over the observation period. Our results suggest that Cas9-mediated suicide gene insertion might be a viable genotype-specific therapy for cancer.

Substance use disorder is associated with disruptions in circadian rhythms. The circadian transcription factor neuronal PAS domain protein 2 (NPAS2) is enriched in reward-related brain regions and regulates reward, but its role in substance use is unclear. To examine the role of NPAS2 in drug taking, we measured intravenous cocaine self-administration (acquisition, dose-response, progressive ratio, extinction, cue-induced reinstatement) in wild-type (WT) andNpas2mutant mice at different times of day. In the light (inactive) phase, cocaine reinforcement was increased in allNpas2mutants, while self-administration and motivation were affected sex-dependently. These sex differences were amplified during the dark (active) phase withNpas2mutation increasing self-administration, reinforcement, motivation, extinction responding and reinstatement in females, but only reinforcement in males. To determine whether circulating hormones are driving these sex differences, we ovariectomized WT andNpas2mutant females and confirmed that unlike sham controls, ovariectomized mutant mice showed no increase in self-administration. To identify whether striatal brain regions are activated inNpas2mutant females, we measured cocaine-induced ΔFosB expression. Relative to WT, ΔFosB expression was increased in D1+ neurons in the nucleus accumbens core and dorsolateral striatum inNpas2mutant females after dark phase self-administration. We also identified potential target genes that may underlie the behavioral responses to cocaine inNpas2mutant females. These results suggest NPAS2 regulates reward and activity in specific striatal regions in a sex and time of day specific manner. Striatal activation could be augmented by circulating sex hormones, leading to an increased effect ofNpas2mutation in females.Significance StatementCircadian disruptions are a common symptom of substance use disorders and chronic exposure to drugs of abuse alters circadian rhythms, which may contribute to subsequent substance use. Diurnal rhythms are commonly found in behavioral responses to drugs of abuse with drug sensitivity and motivation peaking during the dark (active) phase in nocturnal rodents. Emerging evidence links disrupted circadian genes to substance use vulnerability and drug-induced alterations to these genes may augment drug-seeking. The circadian transcription factor NPAS2 is enriched in reward-related brain regions and regulates reward, but its role in substance use is unclear. To examine the role of NPAS2 in drug taking, we measured intravenous cocaine self-administration in wild-type andNpas2mutant mice at different times of day.

Black versus white older Americans are more likely to experience frailty, a condition associated with adverse health outcomes. To reduce racial disparities in health, a complete understanding of the pathophysiology of frailty is needed. Metabolomics may further our understanding by characterizing differences in the body during a vigorous versus frail state. We sought to identify metabolites and biological pathways associated with vigor to frailty among 287 black men ages 70&ndash, 81 from the Health, Aging, and Body Composition study. Using liquid chromatography-mass spectrometry, 350 metabolites were measured in overnight-fasting plasma. The Scale of Aging Vigor in Epidemiology (SAVE) measured vigor to frailty based on weight change, strength, energy, gait speed, and physical activity. Thirty-seven metabolites correlated with SAVE scores (p &lt, 0.05), while adjusting for age and site. Fourteen metabolites remained significant after multiple comparisons adjustment (false discovery rate &lt, 0.30). Lower values of tryptophan, methionine, tyrosine, asparagine, C14:0 sphingomyelin, and 1-methylnicotinamide, and higher values of glucoronate, N-carbamoyl-beta-alanine, isocitrate, creatinine, C4-OH carnitine, cystathionine, hydroxyphenylacetate, and putrescine were associated with frailer SAVE scores. Pathway analyses identified nitrogen metabolism, aminoacyl-tRNA biosynthesis, and the citric acid cycle. Future studies need to confirm these SAVE-associated metabolites and pathways that may indicate novel mechanisms involved in the frailty syndrome.

BackgroundAlternative polyadenylation (APA) shortens or lengthens the 3’-untranslated region (3’-UTR) of hundreds of genes in cancer. While APA genes modify microRNA target sites in the 3’-UTRs to promote tumorigenesis, previous studies have focused on a subset of the modification landscape.MethodFor comprehensive understanding of the function of global APA events, we consider the total target site landscape of microRNAs that are significantly and collectively modified by global APA genes. To identify such microRNAs in spite of complex interactions between microRNAs and the APA genes, we developedProbabilisticInference ofMicroRNATarget Site Modification throughAPA(PRIMATA-APA).ResultsRunning PRIMATA-APA on TCGA breast cancer data, we identified that global APA events concentrate to modify target sites of particular microRNAs (target-site-modified-miRNAor tamoMiRNA). TamoMiRNAs are enriched for microRNAs known to regulate cancer etiology and treatments. Also, their target genes are enriched in cancer-associated pathways, suggesting that APA modifies target sites of tamoMiRNAs to progress tumors. Knockdown of NUDT21, a master 3’-UTR regulator in HeLa cells, confirmed the causal role of tamoMiRNAs for tumor growth.ConclusionsFurther, the expressions of tamoMiRNA target genes, enriched in cancer-associated pathways, vary across tumor samples as a function of patient-specific APA events, suggesting that APA is a novel regulatory axis for interpatient tumor heterogeneity.

Chromosome changes are one of the hallmarks of human malignancies. Chromosomal rearrangement is frequent in human cancers. One of the consequences of chromosomal rearrangement is gene fusions in the cancer genome. We have previously identified a panel of fusion genes in aggressive prostate cancers. In this study, we showed that 6 of these fusion genes are present in 7 different types of human malignancies with variable frequencies. Among them, the CCNH-C5orf30 and TRMT11-GRIK2 gene fusions were found in breast cancer, colon cancer, non-small cell lung cancer, esophageal adenocarcinoma, glioblastoma multiforme, ovarian cancer and liver cancer, with frequencies ranging from 12.9% to 85%. In contrast, four other gene fusions (mTOR-TP53BP1, TMEM135-CCDC67, KDM4-AC011523.2 and LRRC59-FLJ60017) are less frequent. Both TRMT11-GRIK2 and CCNH-C5orf30 are also frequently present in lymph node metastatic cancer samples from the breast, colon and ovary. Thus, detecting these fusion transcripts may have significant biological and clinical implications in cancer patient management.

A consistent gene set undergoes age-associated expression changes in the human cerebral cortex, and our Age-by-Disease Model posits that these changes contribute to psychiatric diseases by “pushing” the expression of disease-associated genes in disease-promoting directions. DNA methylation (DNAm) is an attractive candidate mechanism for age-associated gene expression changes. We used the Illumina HumanMethylation450 array to characterize genome-wide DNAm in the postmortem orbital frontal cortex from 20 younger (60 years) subjects. DNAm data were integrated with existing normal brain aging expression data and sets of psychiatric disease risk genes to test the hypothesis that age-associated DNAm changes contribute to age-associated gene expression changes and, by extension, susceptibility to psychiatric diseases. We found that age-associated differentially methylated regions (aDMRs) are common, robust, bidirectional, concentrated in CpG island shelves and sea, depleted in CpG islands, and enriched among genes undergoing age-associated expression changes (OR = 2.30, p = 1.69 × 10−27). We found the aDMRs are enriched among genetic association-based risk genes for schizophrenia, Alzheimer’s disease (AD), and major depressive disorder (MDD) (OR = 2.51, p = 0.00015; OR = 2.38, p = 0.036; and OR = 3.08, p = 0.018, respectively) as well as expression-based MDD-associated genes (OR = 1.48, p = 0.00012). Similar patterns of enrichment were found for aDMRs that correlate with local gene expression. These results were replicated in a large publically-available dataset, and confirmed by meta-analysis of the two datasets. Our findings suggest DNAm is a molecular mechanism for age-associated gene expression changes and support a role for DNAm in age-by-disease interactions through preferential targeting of disease-associated genes.

Schizophrenia is associated with disrupted cognitive control and sleep-wake cycles. Here we identify diurnal rhythms in gene expression in the human dorsolateral prefrontal cortex (dlPFC), in schizophrenia and control subjects. We find significant diurnal (24 h) rhythms in control subjects, however, most of these transcripts are not rhythmic in subjects with schizophrenia. Instead, subjects with schizophrenia have a different set of rhythmic transcripts. The top pathways identified in transcripts rhythmic only in subjects with schizophrenia are associated with mitochondrial function. Importantly, these rhythms drive differential expression patterns of these and several other genes that have long been implicated in schizophrenia (including BDNF and GABAergic-related transcripts). Indeed, differential expression of these transcripts is only seen in subjects that died during the night, with no change in subjects that died during the day. These data provide insights into a potential mechanism that underlies changes in gene expression in the dlPFC with schizophrenia., Sleep disturbance is common in psychiatric disease, and this may contribute to altered circadian rhythm in gene expression. Here the authors show that rhythms in gene expression in the dorsolateral prefrontal cortex in schizophrenia are different to that seen in healthy controls.

Receptor tyrosine kinases MET and epidermal growth factor receptor (EGFR) are critically involved in initiation of liver regeneration. Other cytokines and signaling molecules also participate in the early part of the process. Regeneration employs effective redundancy schemes to compensate for the missing signals. Elimination of any single extracellular signaling pathway only delays but does not abolish the process. Our present study, however, shows that combined systemic elimination of MET and EGFR signaling (MET knockout + EGFR-inhibited mice) abolishes liver regeneration, prevents restoration of liver mass, and leads to liver decompensation. MET knockout or simply EGFR-inhibited mice had distinct and signaling-specific alterations in Ser/Thr phosphorylation of mammalian target of rapamycin, AKT, extracellular signal-regulated kinases 1/2, phosphatase and tensin homolog, adenosine monophosphate-activated protein kinase α, etc. In the combined MET and EGFR signaling elimination of MET knockout + EGFR-inhibited mice, however, alterations dependent on either MET or EGFR combined to create shutdown of many programs vital to hepatocytes. These included decrease in expression of enzymes related to fatty acid metabolism, urea cycle, cell replication, and mitochondrial functions and increase in expression of glycolysis enzymes. There was, however, increased expression of genes of plasma proteins. Hepatocyte average volume decreased to 35% of control, with a proportional decrease in the dimensions of the hepatic lobules. Mice died at 15-18 days after hepatectomy with ascites, increased plasma ammonia, and very small livers. Conclusion MET and EGFR separately control many nonoverlapping signaling endpoints, allowing for compensation when only one of the signals is blocked, though the combined elimination of the signals is not tolerated; the results provide critical new information on interactive MET and EGFR signaling and the contribution of their combined absence to regeneration arrest and liver decompensation. (Hepatology 2016;64:1711-1724).

With the rapid advances in technologies of microarray and massively parallel sequencing, data of multiple omics sources from a large patient cohort are now frequently seen in many consortium studies. Effective multi-level omics data integration has brought new statistical challenges. One important biological objective of such integrative analysis is to cluster patients in order to identify clinically relevant disease subtypes, which will form basis for tailored treatment and personalized medicine. Several methods have been proposed in the literature for this purpose, including the popular iCluster method used in many cancer applications. When clustering high-dimensional omics data, effective feature selection is critical for better clustering accuracy and biological interpretation. It is also common that a portion of "scattered samples" has patterns distinct from all major clusters and should not be assigned into any cluster as they may represent a rare disease subcategory or be in transition between disease subtypes. In this paper, we firstly propose to improve feature selection of the iCluster factor model by an overlapping sparse group lasso penalty on the omics features using prior knowledge of inter-omics regulatory flows. We then perform regularization over samples to allow clustering with scattered samples and generate tight clusters. The proposed group structured tight iCluster method will be evaluated by two real breast cancer examples and simulations to demonstrate its improved clustering accuracy, biological interpretation, and ability to generate coherent tight clusters.

KALRN (KAL) is a Rho GEF that is highly involved in regulation of the actin cytoskeleton within dendrites. There are several isoforms of the protein that arise from differential splicing of KALRN's 66 exons. KAL isoforms have different functions in development. For example, overexpression of the KAL9 and KAL12 isoforms induce dendritic elongation in early development. However, in mature neurons KAL9 overexpression reduces dendritic length, a phenotype also observed in normal human ageing. We therefore hypothesized that KAL9 would have increased expression with age, and undertook to evaluate the expression of individual KALRN exons throughout the adult lifespan. Postmortem human brain grey matter from Brodmann's area (BA) 11 and BA47 derived from a cohort of 209 individuals without psychiatric or neurodegenerative disease, ranging in age from 16 to 91 years, were analysed for KALRN expression by Affymetrix exon array. Analysis of the exon array data in an isoform-specific manner, as well as confirmatory isoform-specific qPCR studies, indicated that the longer KAL9 and KAL12 isoforms demonstrated a statistically significant, but modest, increase with age. The small magnitude of the age effect suggests that inter-individual factors other than age likely contribute to a higher degree to KAL9 and KAL12 expression. In contrast to KAL9 and KAL12, global KALRN expression did not increase with age. Our work suggests that global measures of KALRN gene expression may be misleading and future studies should focus on isoform-specific quantification.

Mucin1 (MUC1) is an epithelial glycoprotein overexpressed in ovarian cancer and actively involved in tumor cell migration and metastasis. Using novel in vitro and in vivo MUC1-expressing conditional (Cre-loxP) ovarian tumor models, we focus here on MUC1 biology and the roles of Kras activation and Pten deletion during cell transformation and epithelial-to-mesenchymal transition (EMT). We generated several novel murine ovarian cancer cell lines derived from the ovarian surface epithelia (OSE) of mice with conditional mutations in Kras, Pten or both. In addition, we also generated several tumor-derived new cell lines that reproduce the original tumor phenotype in vivo and mirror late stage metastatic disease. Our results demonstrate that de novo activation of oncogenic Kras does not trigger increased proliferation, cellular transformation or EMT, and prevents MUC1 upregulation. In contrast, Pten deletion accelerates cell proliferation, triggers cellular transformation in vitro and in vivo, and stimulates MUC1 expression. Ovarian tumor-derived cell lines MKP-Liver and MKP-Lung cells reproduce in vivo EMT and represent the first immune competent mouse model for distant hematogenous spread. Whole genome microarray expression analysis using tumor and OSE-derived cell lines reveal a 121 gene signature associated with EMT and metastasis. When applied to n=542 cases from The Cancer Genome Atlas (TCGA) ovarian cancer dataset, the gene signature identifies a patient subset with decreased survival (P=0.04). Using an extensive collection of novel murine cell lines we have identified distinct roles for Kras and Pten on MUC1 and EMT in vivo and in vitro. The data has implications for future design of combination therapies targeting Kras mutations, Pten deletions and MUC1 vaccines.

Background: Hotspot estrogen receptor-α (ERα/ESR1) mutations occur in 30-40% endocrine resistant ER+ breast cancer, with Y537S and D538G as the two most frequent hotspot mutations. A mostly unanswered question is and how these mutations facilitate metastatic processes. Methods: Frequencies of different ESR1 hotspot mutations were compared within three publicly available cohorts. Y537S and D538G genome-edited MCF7 and T47D cell models were used for in vitro characterization. Wound scratching, spheroid collective migration, chemotaxis assays were applied to examine different metastatic properties. Tail vein injection in nude mice followed by human CK19 lung section immune staining was used to characterize in vivo metastasis. Canonical Wnt pathway activity was studied using top-flash reporter assays, immunoblotting, and overexpression of dominant negative TCF4. ChIP-seq and ATAC-seq were utilized to profile ER global binding patterns and chromatin accessibility in cell models, respectively. Porcupine inhibitor-LGK974 was used to evaluate the efficiency of Wnt-targeted therapy. Results: D538G mutations were more frequent in ER+ metastatic lesions compared to Y537S among three metastatic breast cancer cohorts with unbiased SNV detection. In line with this, T47D-D538G ESR1 mutant cells showed strongly enhanced cell migration in vitro which was not observed in Y537S cells. We also observed increased lung micro-metastasis in vivo after tail vein injection of the D538G cells. Cell line transcriptomic analysis revealed uniquely hyperactivated Wnt pathway in D538G mutant cells, which was further confirmed in vitro by top-flash reporter and immunoblot. Suppression of canonical Wnt pathways blocked T47D-D538G specific cell migration. Combination treatment of fulvestrant and LGK974 synergistically inhibited D538G specific migration. Mechanistically, multiple Wnt regulator genes were found uniquely upregulated in D538G cells. Interestingly, none of these targets gained ER binding peaks at proximal regulatory regions, suggesting potential epigenetic regulation, which was confirmed in ATAC-seq results. Knockdown of FOXA1, a well-characterized pioneer factor, decreased canonical Wnt activity and abrogated the D538G-unique cell migration in short-term. Conclusion: T47D-D538G cells showed uniquely enhanced cell migration via Wnt hyperactivation, potentially mediated via epigenetic remodeling. These findings suggest the further study of potential combinatorial targeting of Wnt and ER signaling in D538G ESR1 mutant tumors. Citation Format: Zheqi Li, Kevin M. Levine, Spencer Arneson, Kristofer C. Berrett, Nolan M. Priedigkeit, Amir Bahreini, Jian Chen, Li Zhu, Jason S. Carroll, George C. Tseng, Peter C. Lucas, Jennifer M. Atkinson, Jason Gertz, Adrian V. Lee, Steffi Oesterreich. Estrogen receptor D538G mutation promotes cell migration via hyperactivation of Wnt signaling pathway [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 4917.

PURPOSE This study is aimed at developing a molecular diagnostics platform to enhance the interpretation of renal allograft biopsies using quantitative proteomic profiling of formalin-fixed and paraffin-embedded (FFPE) specimens. EXPERIMENTAL DESIGN A quantitative proteomics platform composed of 1) an optimized FFPE protein sample preparation method, 2) a tandem mass tag TMT10-plex-based proteomic workflow, and 3) a systematic statistical analysis pipeline to reveal differentially expressed proteins has been developed. This platform is then tested on a small sample set (five samples per phenotype) to reveal proteomic signatures that can differentiate T-cell mediated rejection (TCMR) and polyomavirus BK nephropathy (BKPyVN) from healthy functionally stable kidney tissue (STA). RESULTS Among 2798 quantified proteins, the expression levels of 740 BKPyVN and 638 TCMR associated proteins are significantly changed compared to STA specimens. Principal component analysis demonstrated good segregation of all three phenotypes investigated. Protein detection and quantitation are highly reproducible: replicate comparative analyses demonstrated 71-84% overlap of detected proteins, and the coefficient of variation for protein measurements is

MET and epidermal growth factor receptor (EGFR) tyrosine kinases are crucial for liver regeneration and normal hepatocyte function. Recently, we demonstrated that in mice, combined inhibition of these two signaling pathways abolished liver regeneration after hepatectomy, with subsequent hepatic failure and death at 15 to 18 days after resection. Morbidity was associated with distinct and specific alterations in important downstream signaling pathways that led to decreased hepatocyte volume, reduced proliferation, and shutdown of many essential hepatocyte functions, such as fatty acid synthesis, urea cycle, and mitochondrial functions. Herein, we explore the role of MET and EGFR signaling in resting mouse livers that are not subjected to hepatectomy. Mice with combined disruption of MET and EGFR signaling were noticeably sick by 10 days and died at 12 to 14 days. Mice with combined disruption of MET and EGFR signaling mice showed decreased liver/body weight ratios, increased apoptosis in nonparenchymal cells, impaired liver metabolic functions, and activation of distinct downstream signaling pathways related to inflammation, cell death, and survival. The present study demonstrates that, in addition to controlling the regenerative response, MET and EGFR synergistically control baseline liver homeostasis in normal mice in such a way that their combined disruption leads to liver failure and death.

Invasive lobular carcinoma (ILC) is the second most common subtype of breast cancer following invasive ductal carcinoma (IDC) and characterized by the loss of E-cadherin-mediated adherens junctions. Despite displaying unique histological and clinical features, ILC still remains a chronically understudied disease with limited knowledge on the available laboratory research models. To this end, herein we report a comprehensive 2D and 3D phenotypic characterization of four Estrogen Receptor-positive human ILC cell lines - MDA-MB-134, SUM44, MDA-MB-330 and BCK4. Compared to the IDC cell lines MCF7, T47D and MDA-MB-231, ultra-low attachment culture conditions revealed a remarkable anchorage-independence ability that was unique to the ILC cells, a feature not evident in soft agar gels. 3D Collagen I and Matrigel culture indicated a generally loose morphology for the ILC cell lines, which exhibited differing preferences for adhesion to ECM proteins in 2D. Furthermore, ILC cells had limited migration and invasion ability in wound-scratch and transwell assays with the exception of haptotaxis to Collagen I. Transcriptional comparison of the cell lines confirmed the decreased cell proliferation and E-cadherin-mediated intercellular junctions in ILC, while uncovering the induction of novel pathways related to cyclic nucleotide phosphodiesterase activity, ion channels, drug metabolism and alternative cell adhesion molecules such as N-cadherin, some of which were also differentially regulated in ILC versus IDC tumors. Altogether, these studies will serve as an invaluable resource for the breast cancer research community and facilitate further functional discoveries towards understanding ILC, identifying novel drug targets and ultimately improving the outcome of patients with ILC.Authors’ ContributionsConception and design:N. Tasdemir, NE. Davidson, S. OesterreichDevelopment of methodology:N. Tasdemir, L. Zhu, GC. Tseng, S. OesterreichAcquisition of data (performed experiments, processed data, etc.):N. Tasdemir, E. Bossart, Z. Li, Z. LiAnalysis and interpretation of data (e.g. biological interpretation, statistical analysis, computational analysis):N. Tasdemir, Z. Li, KM. Levine, NE. Davidson, S. OesterreichWriting, review and/or revision of the manuscript:N. Tasdemir, Z. Li, KM. Levine, BM. Jacobson, GC. Tseng, NE. Davidson, S. OesterreichStudy supervision:NE. Davidson and S. Oesterreich

Invasive lobular carcinoma (ILC) is the second most common histological subtype of breast cancer following invasive ductal carcinoma (IDC). ILC differs from IDC in a number of histological and clinical features, such as single strand growth, difficulty in detection, and frequent late recurrences. To understand the molecular pathways involved in the clinical characteristics of ILC, we compared the gene expression profiles of luminal A ILC and luminal A IDC using data from TCGA and utilized samples from METABRIC as a validation data set. Top pathways that were significantly enriched in ILC were related to immune response. ILC exhibited a higher activity of almost all types of immune cells based on cell type-specific signatures compared to IDC. Conversely, pathways that were less enriched in ILC were related to protein translation and metabolism, which we functionally validated in cell lines. The higher immune activity uncovered in our study highlights the currently unexplored potential of a response to immunotherapy in a subset of patients with ILC. Furthermore, the lower rates of protein translation and metabolism - known features of tumor dormancy - may play a role in the late recurrences of ILC and lower detection rate in mammography and PET scanning.

Genome-wide association studies are proven tools for finding disease genes, but it is often necessary to combine many cohorts into a meta-analysis to detect statistically significant genetic effects. Often the component studies are performed by different investigators on different populations, using different chips with minimal SNPs overlap. In some cases, raw data are not available for imputation so that only the genotyped single nucleotide polymorphisms (SNPs) results can be used in meta-analysis. Even when SNP sets are comparable, different cohorts may have peak association signals at different SNPs within the same gene due to population differences in linkage disequilibrium or environmental interactions. We hypothesize that the power to detect statistical signals in these situations will improve by using a method that simultaneously meta-analyzes and smooths the signal over nearby markers. In this study, we propose regionally smoothed meta-analysis methods and compare their performance on real and simulated data.

Fusion transcript formation is one of the fundamental mechanisms that drives the development of prostate cancer. Because of the advance of high-throughput parallel sequencing, many fusion transcripts have been discovered. However, the discovery rate of fusion transcripts specific for prostate cancer is lagging behind the discoveries made on chromosome abnormalities of prostate cancer. Recent analyses suggest that many fusion transcripts are present in both benign and cancerous tissues. Some of these fusion transcripts likely represent important components of normal gene expression in cells. It is necessary to identify the criteria and features of fusion transcripts that are specific for cancer. In this review, we discuss optimization of RNA sequencing depth for fusion transcript discovery and the characteristics of fusion transcripts in normal prostate tissues and prostate cancer. We also propose a new classification of cancer-specific fusion transcripts on the basis of their tail gene fusion protein product and the roles that these fusions may play in cancer development.

Genome-wide expression and genotyping technologies have uncovered the genetic bases of complex diseases at unprecedented rates. However, despite its heavy burden and high prevalence, the molecular characterization of major depressive disorder (MDD) has lagged behind. Transcriptome studies report multiple brain disturbances but are limited by small sample sizes. Genome-wide association studies (GWAS) report weak results but suggest an overlapping genetic risk with other neuropsychiatric disorders. We performed a systematic molecular characterization of altered brain function in MDD using meta-analysis of differential expression of 8 gene array studies across 3 corticolimbic brain regions in 101 subjects. The identified ‘metaA-MDD' genes suggest altered neurotrophic support, brain plasticity and neuronal signaling in MDD. Notably, metaA-MDD genes display a low connectivity and hubness in coexpression networks as well as a uniform genomic distribution, which is consistent with diffuse polygenic mechanisms. We have integrated these findings with results from over 1,800 published GWAS and show that genetic variations nearby metaA-MDD genes predict a greater risk for neuropsychiatric disorders, and notably for age-related phenotypes, but not for other medical illnesses (including those frequently co-occurring with depression) or body characteristics. Collectively, the intersection of unbiased investigations of gene function (transcriptome) and structure (GWAS) provides novel leads to investigate molecular mechanisms of MDD and suggests common biological pathways between depression, other neuropsychiatric diseases and brain aging.

Schizophrenia is associated with alterations in working memory that reflect dysfunction of dorsolateral prefrontal cortex (DLPFC) circuitry. Working memory depends on the activity of excitatory pyramidal cells in DLPFC layer 3 and, to a lesser extent, in layer 5. Although many studies have profiled gene expression in DLPFC gray matter in schizophrenia, little is known about cell-type-specific transcript expression in these two populations of pyramidal cells. We hypothesized that interrogating gene expression, specifically in DLPFC layer 3 or 5 pyramidal cells, would reveal new and/or more robust schizophrenia-associated differences that would provide new insights into the nature of pyramidal cell dysfunction in the illness. We also sought to determine the impact of other variables, such as a diagnosis of schizoaffective disorder or medication use at the time of death, on the patterns of gene expression in pyramidal neurons. Individual pyramidal cells in DLPFC layers 3 or 5 were captured by laser microdissection from 36 subjects with schizophrenia or schizoaffective disorder and matched normal comparison subjects. The mRNA from cell collections was subjected to transcriptome profiling by microarray followed by quantitative PCR validation. Expression of genes involved in mitochondrial (MT) or ubiquitin-proteasome system (UPS) functions were markedly downregulated in the patient group (P-values for MT-related and UPS-related pathways were10(-7) and10(-5), respectively). MT-related gene alterations were more prominent in layer 3 pyramidal cells, whereas UPS-related gene alterations were more prominent in layer 5 pyramidal cells. Many of these alterations were not present, or found to a lesser degree, in samples of DLPFC gray matter from the same subjects, suggesting that they are pyramidal cell specific. Furthermore, these findings principally reflected alterations in the schizophrenia subjects were not present or present to a lesser degree in the schizoaffective disorder subjects (diagnosis of schizoaffective disorder was the most significant covariate, P10(-6)) and were not attributable to factors frequently comorbid with schizophrenia. In summary, our findings reveal expression deficits in MT- and UPS-related genes specific to layer 3 and/or layer 5 pyramidal cells in the DLPFC of schizophrenia subjects. These cell type-specific transcriptome signatures are not characteristic of schizoaffective disorder, providing a potential molecular-cellular basis of differences in clinical phenotypes.

Genome-wide expression profiling of the human brain has revealed genes that are differentially expressed across the lifespan. Characterizing these genes adds to our understanding of both normal functions and pathological conditions. Additionally, the specific cell-types that contribute to the motor, sensory and cognitive declines during aging are unclear. Here we test if age-related genes show higher expression in specific neural cell types. Our study leverages data from two sources of murine single-cell expression data and two sources of age-associations from large gene expression studies of postmortem human brain. We used nonparametric gene set analysis to test for age-related enrichment of genes associated with specific cell-types; we also restricted our analyses to specific gene ontology groups. Our analyses focused on a primary pair of single-cell expression data from the mouse visual cortex and age-related human post-mortem gene expression information from the orbitofrontal cortex. Additional pairings that used data from the hippocampus, prefrontal cortex, somatosensory cortex and blood were used to validate and test specificity of our findings. We found robust age-related up-regulation of genes that are highly expressed in oligodendrocytes and astrocytes, while genes highly expressed in layer 2/3 glutamatergic neurons were down-regulated across age. Genes not specific to any neural cell type were also down-regulated, possibly due to the bulk tissue source of the age-related genes. A gene ontology-driven dissection of the cell-type enriched genes highlighted the strong down-regulation of genes involved in synaptic transmission and cell-cell signalling in the Somatostatin (Sst) neuron subtype that expresses the cyclin dependent kinase 6 (Cdk6) and in the vasoactive intestinal peptide (Vip) neuron subtype expressing myosin binding protein C, slow type (Mybpc1). These findings provide new insights into cell specific susceptibility to normal aging, and suggest age-related synaptic changes in specific inhibitory neuron subtypes.

Purpose: High-grade serous ovarian cancer (HGSOC) is an aggressive disease with few available targeted therapies. Despite high expression of estrogen receptor-alpha (ERα) in approximately 80% of HGSOC and some small but promising clinical trials of endocrine therapy, ERα has been understudied as a target in this disease. We sought to identify hormone-responsive, ERα-dependent HGSOC. Experimental Design: We characterized endocrine response in HGSOC cells across culture conditions [ two-dimensional (2D), three-dimensional (3D), forced suspension] and in patient-derived xenograft (PDX) explants, assessing proliferation and gene expression. Estrogen-regulated transcriptome data were overlapped with public datasets to develop a comprehensive panel of ERα target genes. Expression of this panel and ERα H-score were assessed in HGSOC samples from patients who received endocrine therapy. Time on endocrine therapy was used as a surrogate for clinical response. Results: Proliferation is ERα-regulated in HGSOC cells in vitro and in vivo, and is partly dependent on 3D context. Transcriptomic studies identified genes shared by cell lines and PDX explants as ERα targets. The selective ERα downregulator (SERD) fulvestrant is more effective than tamoxifen in blocking ERα action. ERα H-score is predictive of efficacy of endocrine therapy, and this prediction is further improved by inclusion of target gene expression, particularly IGFBP3. Conclusions: Laboratory models corroborate intertumor heterogeneity of endocrine response in HGSOC but identify features associated with functional ERα and endocrine responsiveness. Assessing ERα function (e.g., IGFBP3 expression) in conjunction with H-score may help select patients who would benefit from endocrine therapy. Preclinical data suggest that SERDs might be more effective than tamoxifen. Clin Cancer Res; 23(14); 3802–12. ©2017 AACR.

Purpose: Mechanisms of immune dysregulation associated with advanced tumors are relatively well understood. Much less is known about the role of immune effectors against cancer precursor lesions. Endometrioid and clear-cell ovarian tumors partly derive from endometriosis, a commonly diagnosed chronic inflammatory disease. We performed here a comprehensive immune gene expression analysis of pelvic inflammation in endometriosis and endometriosis-associated ovarian cancer (EAOC). Experimental Design: RNA was extracted from 120 paraffin tissue blocks comprising of normal endometrium (n = 32), benign endometriosis (n = 30), atypical endometriosis (n = 15), and EAOC (n = 43). Serous tumors (n = 15) were included as nonendometriosis-associated controls. The immune microenvironment was profiled using Nanostring and the nCounter GX Human Immunology Kit, comprising probes for a total of 511 immune genes. Results: One third of the patients with endometriosis revealed a tumor-like inflammation profile, suggesting that cancer-like immune signatures may develop earlier, in patients classified as clinically benign. Gene expression analyses revealed the complement pathway as most prominently involved in both endometriosis and EAOC. Complement proteins are abundantly present in epithelial cells in both benign and malignant lesions. Mechanistic studies in ovarian surface epithelial cells from mice with conditional (Cre-loxP) mutations show intrinsic production of complement in epithelia and demonstrate an early link between Kras- and Pten-driven pathways and complement upregulation. Downregulation of complement in these cells interferes with cell proliferation. Conclusions: These findings reveal new characteristics of inflammation in precursor lesions and point to previously unknown roles of complement in endometriosis and EAOC. Clin Cancer Res; 20(23); 6163–74. ©2014 AACR.

Objective To determine whether treatment for urinary tract infections in children could be individualized using biomarkers for acute pyelonephritis. Study design We enrolled 61 children with febrile urinary tract infections, collected blood and urine samples, and performed a renal scan within 2 weeks of diagnosis to identify those with pyelonephritis. Renal scans were interpreted centrally by 2 experts. We measured inflammatory proteins in blood and urine using LUMINEX or an enzyme-linked immunosorbent assay. We evaluated serum RNA expression using RNA sequencing in a subset of children. Finally, for children with Escherichia coli isolated from urine cultures, we performed a polymerase chain reaction for 4 previously identified virulence genes. Results Urinary markers that best differentiated pyelonephritis from cystitis included chemokine (C-X-C motif) ligand (CXCL)1, CXCL9, CXCL12, C-C motif chemokine ligand 2, INF γ, and IL-15. Serum procalcitonin was the best serum marker for pyelonephritis. Genes in the interferon-γ pathway were upregulated in serum of children with pyelonephritis. The presence of E coli virulence genes did not correlate with pyelonephritis. Conclusions Immune response to pyelonephritis and cystitis differs quantitatively and qualitatively; this may be useful in differentiating these 2 conditions.

Invasive lobular carcinoma (ILC) is the second most common type of breast cancer following invasive ductal carcinoma (IDC) and accounts for 10-15% of all cases. Unlike the masses or lumps formed by IDCs, ILCs grow as small, dyscohesive cells in a single-file pattern within dense layers of extracellular matrix. Paradoxically, while patients with ILC display favorable prognostic and predictive factors (Estrogen Receptor [ER]+, Progesterone Receptor+, HER2-, low Ki67), they present with more frequent long-term recurrences compared to those with IDC, indicative of endocrine resistance. Thus, there is urgent need to investigate genetic drivers of ILC disease progression in order to develop more effective therapeutic strategies and improve patient outcome. To this end, we used the Nanostring platform to measure the expression of 577 copy number variation-associated genes in 131 primary ILC tumors with long-term clinical data. This analysis identified CTTN (Cortactin; cortical actin binding protein) as a candidate ILC driver that exhibited higher expression in tumors from patients with subsequent recurrent (n=33) versus non-current disease (n=98). We further validated high CTTN mRNA/protein expression in human ILC cell lines, tumors and patient-derived xenografts, and CTTN locus (11q13.3) amplification in clinical ILC metastases to the brain, bone and ovaries. In follow-upin vitro studies, RNAi-mediated inhibition of CTTN diminished the adhesion and haptotaxis of human ILC cell lines to Collagen I, as well as impairing their growth in 3D Collagen I culture. To assess the functional role of CTTN in ILC tumor growth and metastasis, we transplanted control and CTTN knockdown human ILC cell lines into the mammary fat pads and ducts of immuno-compromised mice and monitored disease progression via bioluminescence imaging. While CTTN inhibition did not lead to a substantial impact on measurable tumor burden, both CTTN shRNAs resulted in a significant survival benefit in the fat pad model. Ongoing work is aimed at pinpointing the underlying mechanisms of CTTN's role in mediating growth in 3D Collagen I culture in vitro and disease progression in vivo. Collectively, the results of these studies will advance our understanding of ILC disease mechanisms and serve as a pre-clinical basis for improving the clinical outcome of patients with this understudied subtype of breast cancer. Citation Format: Tasdemir N, Sikora MJ, Zhu L, Levine KM, Scott J, Basudan A, Sflomos G, Sreekumar S, Bossart EA, Elishaev E, Chandran UR, Tseng GC, Jankowitz RC, Dabbs DJ, McAuliffe PF, Brisken C, Davidson NE, Oesterreich S. Investigating cortactin as a genetic driver of disease progression in invasive lobular carcinoma [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr PD7-03.

Metallothioneins (MTs) are a group of metal binding proteins thought to play a role in the detoxification of heavy metals. Here we showed by microarray and validation analyses that MT1h, a member of MT, is down-regulated in many human malignancies. Low expression of MT1h was associated with poor clinical outcomes in both prostate and liver cancer. We found that the promoter region of MT1h was hypermethylated in cancer and that demethylation of the MT1h promoter reversed the suppression of MT1h expression. Forced expression of MT1h induced cell growth arrest, suppressed colony formation, retarded migration, and reduced invasion. SCID mice with tumour xenografts with inducible MT1h expression had lower tumour volumes as well as fewer metastases and deaths than uninduced controls. MT1h was found to interact with euchromatin histone methyltransferase 1 (EHMT1) and enhanced its methyltransferase activity on histone 3. Knocking down of EHMT1 or a mutation in MT1h that abrogates its interaction with EHMT1 abrogated MT1h tumour suppressor activity. This demonstrates tumour suppressor activity in a heavy metal binding protein that is dependent on activation of histone methylation.