KDM4A lysine demethylase induces site-specific copy gain and rereplication of regions amplified in tumors.

Acquired chromosomal instability and copy number alterations are hallmarks of cancer. Enzymes capable of promoting site-specific copy number changes have yet to be identified. Here, we demonstrate that H3K9/36me3 lysine demethylase KDM4A/JMJD2A overexpression leads to localized copy gain of 1q12, 1q21, and Xq13.1 without global chromosome instability. KDM4A-amplified tumors have increased copy gains for these same regions. 1q12h copy gain occurs within a single cell cycle, requires S phase, and is not stable but is regenerated each cell division. Sites with increased copy number are rereplicated and have increased KDM4A, MCM, and DNA polymerase occupancy. Suv39h1/KMT1A or HP1γ overexpression suppresses the copy gain, whereas H3K9/K36 methylation interference promotes gain. Our results demonstrate that overexpression of a chromatin modifier results in site-specific copy gains. This begins to establish how copy number changes could originate during tumorigenesis and demonstrates that transient overexpression of specific chromatin modulators could promote these events.

Detecting and dissecting signaling crosstalk via the multilayer network integration of signaling and regulatory interactions.

The versatility of cellular response arises from the communication, or crosstalk, of signaling pathways in a complex network of signaling and transcriptional regulatory interactions. Understanding the various mechanisms underlying crosstalk on a global scale requires untargeted computational approaches. We present a network-based statistical approach, MuXTalk, that uses high-dimensional edges called multilinks to model the unique ways in which signaling and regulatory interactions can interface. We demonstrate that the signaling-regulatory interface is located primarily in the intermediary region between signaling pathways where crosstalk occurs, and that multilinks can differentiate between distinct signaling-transcriptional mechanisms. Using statistically over-represented multilinks as proxies of crosstalk, we infer crosstalk among 60 signaling pathways, expanding currently available crosstalk databases by more than five-fold. MuXTalk surpasses existing methods in terms of model performance metrics, identifies additions to manual curation efforts, and pinpoints potential mediators of crosstalk. Moreover, it accommodates the inherent context-dependence of crosstalk, allowing future applications to cell type- and disease-specific crosstalk.

Predicting genotype-specific gene regulatory networks.

Understanding how each person's unique genotype influences their individual patterns of gene regulation has the potential to improve our understanding of human health and development, and to refine genotype-specific disease risk assessments and treatments. However, the effects of genetic variants are not typically considered when constructing gene regulatory networks, despite the fact that many disease-associated genetic variants are thought to have regulatory effects, including the disruption of transcription factor (TF) binding. We developed EGRET (Estimating the Genetic Regulatory Effect on TFs), which infers a genotype-specific gene regulatory network for each individual in a study population. EGRET begins by constructing a genotype-informed TF-gene prior network derived using TF motif predictions, expression quantitative trait locus (eQTL) data, individual genotypes, and the predicted effects of genetic variants on TF binding. It then uses a technique known as message passing to integrate this prior network with gene expression and TF protein-protein interaction data to produce a refined, genotype-specific regulatory network. We used EGRET to infer gene regulatory networks for two blood-derived cell lines and identified genotype-associated, cell line-specific regulatory differences that we subsequently validated using allele-specific expression, chromatin accessibility QTLs, and differential ChIP-seq TF binding. We also inferred EGRET networks for three cell types from each of 119 individuals and identified cell type-specific regulatory differences associated with diseases related to those cell types. EGRET is, to our knowledge, the first method that infers networks reflective of individual genetic variation in a way that provides insight into the genetic regulatory associations driving complex phenotypes.

Partial correlation network analysis identifies coordinated gene expression within a regional cluster of COPD genome-wide association signals.

Chronic obstructive pulmonary disease (COPD) is a complex disease influenced by well-established environmental exposures (most notably, cigarette smoking) and incompletely defined genetic factors. The chromosome 4q region harbors multiple genetic risk loci for COPD, including signals near HHIP, FAM13A, GSTCD, TET2, and BTC. Leveraging RNA-Seq data from lung tissue in COPD cases and controls, we estimated the co-expression network for genes in the 4q region bounded by HHIP and BTC (~70MB), through partial correlations informed by protein-protein interactions. We identified several co-expressed gene pairs based on partial correlations, including NPNT-HHIP, BTC-NPNT and FAM13A-TET2, which were replicated in independent lung tissue cohorts. Upon clustering the co-expression network, we observed that four genes previously associated to COPD: BTC, HHIP, NPNT and PPM1K appeared in the same network community. Finally, we discovered a sub-network of genes differentially co-expressed between COPD vs controls (including FAM13A, PPA2, PPM1K and TET2). Many of these genes were previously implicated in cell-based knock-out experiments, including the knocking out of SPP1 which belongs to the same genomic region and could be a potential local key regulatory gene. These analyses identify chromosome 4q as a region enriched for COPD genetic susceptibility and differential co-expression.

Sex-biased Regulation of Extracellular Matrix Genes in COPD.

Compared to men, women often develop COPD at an earlier age with worse respiratory symptoms despite lower smoking exposure. However, most preventive, and therapeutic strategies ignore biological sex differences in COPD. Our goal was to better understand sex-specific gene regulatory processes in lung tissue and the molecular basis for sex differences in COPD onset and severity. We analyzed lung tissue gene expression and DNA methylation data from 747 individuals in the Lung Tissue Research Consortium (LTRC), and 85 individuals in an independent dataset. We identified sex differences in COPD-associated gene regulation using gene regulatory networks. We used linear regression to test for sex-biased associations of methylation with lung function, emphysema, smoking, and age. Analyzing gene regulatory networks in the control group, we identified that genes involved in the extracellular matrix (ECM) have higher transcriptional factor targeting in females than in males. However, this pattern is reversed in COPD, with males showing stronger regulatory targeting of ECM-related genes than females. Smoking exposure, age, lung function, and emphysema were all associated with sex-specific differential methylation of ECM-related genes. We identified sex-based gene regulatory patterns of ECM-related genes associated with lung function and emphysema. Multiple factors including epigenetics, smoking, aging, and cell heterogeneity influence sex-specific gene regulation in COPD. Our findings underscore the importance of considering sex as a key factor in disease susceptibility and severity.

Adjustment of spurious correlations in co-expression measurements from RNA-Sequencing data.

MOTIVATION: Gene co-expression measurements are widely used in computational biology to identify coordinated expression patterns across a group of samples. Coordinated expression of genes may indicate that they are controlled by the same transcriptional regulatory program, or involved in common biological processes. Gene co-expression is generally estimated from RNA-Sequencing data, which are commonly normalized to remove technical variability. Here, we demonstrate that certain normalization methods, in particular quantile-based methods, can introduce false-positive associations between genes. These false-positive associations can consequently hamper downstream co-expression network analysis. Quantile-based normalization can, however, be extremely powerful. In particular, when preprocessing large-scale heterogeneous data, quantile-based normalization methods such as smooth quantile normalization can be applied to remove technical variability while maintaining global differences in expression for samples with different biological attributes. RESULTS: We developed SNAIL (Smooth-quantile Normalization Adaptation for the Inference of co-expression Links), a normalization method based on smooth quantile normalization specifically designed for modeling of co-expression measurements. We show that SNAIL avoids formation of false-positive associations in co-expression as well as in downstream network analyses. Using SNAIL, one can avoid arbitrary gene filtering and retain associations to genes that only express in small subgroups of samples. This highlights the method's potential future impact on network modeling and other association-based approaches in large-scale heterogeneous data. AVAILABILITY AND IMPLEMENTATION: The implementation of the SNAIL algorithm and code to reproduce the analyses described in this work can be found in the GitHub repository https://github.com/kuijjerlab/PySNAIL.

GRAND: a database of gene regulatory network models across human conditions.

Gene regulation plays a fundamental role in shaping tissue identity, function, and response to perturbation. Regulatory processes are controlled by complex networks of interacting elements, including transcription factors, miRNAs and their target genes. The structure of these networks helps to determine phenotypes and can ultimately influence the development of disease or response to therapy. We developed GRAND (https://grand.networkmedicine.org) as a database for computationally-inferred, context-specific gene regulatory network models that can be compared between biological states, or used to predict which drugs produce changes in regulatory network structure. The database includes 12 468 genome-scale networks covering 36 human tissues, 28 cancers, 1378 unperturbed cell lines, as well as 173 013 TF and gene targeting scores for 2858 small molecule-induced cell line perturbation paired with phenotypic information. GRAND allows the networks to be queried using phenotypic information and visualized using a variety of interactive tools. In addition, it includes a web application that matches disease states to potentially therapeutic small molecule drugs using regulatory network properties.

Transcriptional changes in human palate and skin healing.

Most human tissue injuries lead to the formation of a fibrous scar and result in the loss of functional tissue. One adult tissue that exhibits a more regenerative response to injury with minimal scarring is the oral mucosa. We generated a microarray gene expression dataset to examine the response to injury in human palate and skin excisional biopsies spanning the first 7 days after wounding. Differential expression analyses were performed in each tissue to identify genes overexpressed or underexpressed over time when compared to baseline unwounded tissue gene expression levels. To attribute biological processes of interest to these gene expression changes, gene set enrichment analysis was used to identify core gene sets that are enriched over the time-course of the wound healing process with respect to unwounded tissue. This analysis identified gene sets uniquely enriched in either palate or skin wounds and gene sets that are enriched in both tissues in at least one time point after injury. Finally, a cell type enrichment analysis was performed to better understand the cell type distribution in these tissues and how it changes over the time course of wound healing. This work provides a source of human wound gene expression data that includes two tissue types with distinct regenerative and scarring phenotypes.

Developmental control of polycomb subunit composition by GATA factors mediates a switch to non-canonical functions.

Polycomb repressive complex 2 (PRC2) plays crucial roles in transcriptional regulation and stem cell development. However, the context-specific functions associated with alternative subunits remain largely unexplored. Here we show that the related enzymatic subunits EZH1 and EZH2 undergo an expression switch during blood cell development. An erythroid-specific enhancer mediates transcriptional activation of EZH1, and a switch from GATA2 to GATA1 controls the developmental EZH1/2 switch by differential association with EZH1 enhancers. We further examine the in vivo stoichiometry of the PRC2 complexes by quantitative proteomics and reveal the existence of an EZH1-SUZ12 subcomplex lacking EED. EZH1 together with SUZ12 form a non-canonical PRC2 complex, occupy active chromatin, and positively regulate gene expression. Loss of EZH2 expression leads to repositioning of EZH1 to EZH2 targets. Thus, the lineage- and developmental stage-specific regulation of PRC2 subunit composition leads to a switch from canonical silencing to non-canonical functions during blood stem cell specification.

Multiplex genome editing of mammalian cells for producing recombinant heparin.

Heparin is an essential anticoagulant used for treating and preventing thrombosis. However, the complexity of heparin has hindered the development of a recombinant source, making its supply dependent on a vulnerable animal population. In nature, heparin is produced exclusively in mast cells, which are not suitable for commercial production, but mastocytoma cells are readily grown in culture and make heparan sulfate, a closely related glycosaminoglycan that lacks anticoagulant activity. Using gene expression profiling of mast cells as a guide, a multiplex genome engineering strategy was devised to produce heparan sulfate with high anticoagulant potency and to eliminate contaminating chondroitin sulfate from mastocytoma cells. The heparan sulfate purified from engineered cells grown in chemically defined medium has anticoagulant potency that exceeds porcine-derived heparin and confers anticoagulant activity to the blood of healthy mice. This work demonstrates the feasibility of producing recombinant heparin from mammalian cell culture as an alternative to animal sources.

Protein interaction networks provide insight into fetal origins of chronic obstructive pulmonary disease.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is a leading cause of death in adults that may have origins in early lung development. It is a complex disease, influenced by multiple factors including genetic variants and environmental factors. Maternal smoking during pregnancy may influence the risk for diseases during adulthood, potentially through epigenetic modifications including methylation. METHODS: In this work, we explore the fetal origins of COPD by utilizing lung DNA methylation marks associated with in utero smoke (IUS) exposure, and evaluate the network relationships between methylomic and transcriptomic signatures associated with adult lung tissue from former smokers with and without COPD. To identify potential pathobiological mechanisms that may link fetal lung, smoke exposure and adult lung disease, we study the interactions (physical and functional) of identified genes using protein-protein interaction networks. RESULTS: We build IUS-exposure and COPD modules, which identify connected subnetworks linking fetal lung smoke exposure to adult COPD. Studying the relationships and connectivity among the different modules for fetal smoke exposure and adult COPD, we identify enriched pathways, including the AGE-RAGE and focal adhesion pathways. CONCLUSIONS: The modules identified in our analysis add new and potentially important insights to understanding the early life molecular perturbations related to the pathogenesis of COPD. We identify AGE-RAGE and focal adhesion as two biologically plausible pathways that may reveal lung developmental contributions to COPD. We were not only able to identify meaningful modules but were also able to study interconnections between smoke exposure and lung disease, augmenting our knowledge about the fetal origins of COPD.

Effects of novel SGLT2 inhibitors on cancer incidence in hyperglycemic patients: a meta-analysis of randomized clinical trials.

Epidemiological evidence shows that diabetic patients have an increased cancer risk and a higher mortality rate. Glucose could play a central role in metabolism and growth of many tumor types, and this possible mechanism is supported by the high rate of glucose demand and uptake in cancer. Thus, growing evidence suggests that hyperglycemia contributes to cancer progression but also to its onset. Many mechanisms underlying this association have been hypothesized, such as insulin resistance, hyperinsulinemia, and increased inflammatory processes. Inflammation is a common pathophysiological feature in both diabetic and oncological patients, and inflammation linked to high glucose levels sensitizes microenvironment to tumorigenesis, promoting the development of malignant lesions by altering and sustaining a pathological condition in tissues. Glycemic control is the first goal of antidiabetic therapy, and glucose level reduction has also been associated with favorable outcomes in cancer. Here, we describe key events in carcinogenesis focusing on hyperglycemia as supporter in tumor progression and in particular, related to the role of a specific hypoglycemic drug class, sodium-glucose linked transporters (SGLTs). We also discuss the use of SGLT2 inhibitors as a novel potential cancer therapy. Our meta-analysis showed that SGLT-2 inhibitors were significantly associated with an overall reduced risk of cancer as compared to placebo (RR = 0.35, CI 0.33-0.37, P = 0. 00) with a particular effectiveness for dapaglifozin and ertuglifozin (RR = 0. 06, CI 0. 06-0. 07 and RR = 0. 22, CI 0. 18-0. 26, respectively). Network Medicine approaches may advance the possible repurposing of these drugs in patients with concomitant diabetes and cancer.

Connecting COPD GWAS Genes: FAM13A Controls TGFß2 Secretion by Modulating AP-3 Transport.

Chronic obstructive pulmonary disease (COPD) is a common, complex disease and a major cause of morbidity and mortality. Although multiple genetic determinants of COPD have been implicated by genome-wide association studies (GWASs), the pathophysiological significance of these associations remains largely unknown. From a COPD protein-protein interaction network module, we selected a network path between two COPD GWAS genes for validation studies: FAM13A (family with sequence similarity 13 member A)-AP3D1-CTGF- TGFß2. We find that TGFß2, FAM13A, and AP3D1 (but not CTGF) form a cellular protein complex. Functional characterization suggests that this complex mediates the secretion of TGFß2 through an AP-3 (adaptor protein 3)-dependent pathway, with FAM13A acting as a negative regulator by targeting a late stage of this transport that involves the dissociation of coat-cargo interaction. Moreover, we find that TGFß2 is a transmembrane protein that engages the AP-3 complex for delivery to the late endosomal compartments for subsequent secretion through exosomes. These results identify a pathophysiological context that unifies the biological network role of two COPD GWAS proteins and reveal novel mechanisms of cargo transport through an intracellular pathway.

Epigenetics and pulmonary diseases in the horizon of precision medicine: a review.

Epigenetic mechanisms represent potential molecular routes which could bridge the gap between genetic background and environmental risk factors contributing to the pathogenesis of pulmonary diseases. In patients with COPD, asthma and pulmonary arterial hypertension (PAH), there is emerging evidence of aberrant epigenetic marks, mainly including DNA methylation and histone modifications which directly mediate reversible modifications to the DNA without affecting the genomic sequence. Post-translational events and microRNAs can be also regulated epigenetically and potentially participate in disease pathogenesis. Thus, novel pathogenic mechanisms and putative biomarkers may be detectable in peripheral blood, sputum, nasal and buccal swabs or lung tissue. Besides, DNA methylation plays an important role during the early phases of fetal development and may be impacted by environmental exposures, ultimately influencing an individual's susceptibility to COPD, asthma and PAH later in life. With the advances in omics platforms and the application of computational biology tools, modelling the epigenetic variability in a network framework, rather than as single molecular defects, provides insights into the possible molecular pathways underlying the pathogenesis of COPD, asthma and PAH. Epigenetic modifications may have clinical applications as noninvasive biomarkers of pulmonary diseases. Moreover, combining molecular assays with network analysis of epigenomic data may aid in clarifying the multistage transition from a "pre-disease" to "disease" state, with the goal of improving primary prevention of lung diseases and its subsequent clinical management.We describe epigenetic mechanisms known to be associated with pulmonary diseases and discuss how network analysis could improve our understanding of lung diseases.

PUMA: PANDA Using MicroRNA Associations.

MOTIVATION: Conventional methods to analyze genomic data do not make use of the interplay between multiple factors, such as between microRNAs (miRNAs) and the messenger RNA (mRNA) transcripts they regulate, and thereby often fail to identify the cellular processes that are unique to specific tissues. We developed PUMA (PANDA Using MicroRNA Associations), a computational tool that uses message passing to integrate a prior network of miRNA target predictions with target gene co-expression information to model genome-wide gene regulation by miRNAs. We applied PUMA to 38 tissues from the Genotype-Tissue Expression project, integrating RNA-Seq data with two different miRNA target predictions priors, built on predictions from TargetScan and miRanda, respectively. We found that while target predictions obtained from these two different resources are considerably different, PUMA captures similar tissue-specific miRNA-target regulatory interactions in the different network models. Furthermore, the tissue-specific functions of miRNAs we identified based on regulatory profiles (available at: https://kuijjer.shinyapps.io/puma_gtex/) are highly similar between networks modeled on the two target prediction resources. This indicates that PUMA consistently captures important tissue-specific miRNA regulatory processes. In addition, using PUMA we identified miRNAs regulating important tissue-specific processes that, when mutated, may result in disease development in the same tissue. AVAILABILITY AND IMPLEMENTATION: PUMA is available in C++, MATLAB and Python on GitHub (https://github.com/kuijjerlab and https://netzoo.github.io/). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

ETV1 directs androgen metabolism and confers aggressive prostate cancer in targeted mice and patients.

Distinguishing aggressive from indolent disease and developing effective therapy for advanced disease are the major challenges in prostate cancer research. Chromosomal rearrangements involving ETS transcription factors, such as ERG and ETV1, occur frequently in prostate cancer. How they contribute to tumorigenesis and whether they play similar or distinct in vivo roles remain elusive. Here we show that in mice with ERG or ETV1 targeted to the endogenous Tmprss2 locus, either factor cooperated with loss of a single copy of Pten, leading to localized cancer, but only ETV1 appeared to support development of invasive adenocarcinoma under the background of full Pten loss. Mechanistic studies demonstrated that ERG and ETV1 control a common transcriptional network but largely in an opposing fashion. In particular, while ERG negatively regulates the androgen receptor (AR) transcriptional program, ETV1 cooperates with AR signaling by favoring activation of the AR transcriptional program. Furthermore, we found that ETV1 expression, but not that of ERG, promotes autonomous testosterone production. Last, we confirmed the association of an ETV1 expression signature with aggressive disease and poorer outcome in patient data. The distinct biology of ETV1-associated prostate cancer suggests that this disease class may require new therapies directed to underlying programs controlled by ETV1.

Genome-Wide Association Study for Urinary and Fecal Incontinence in Women.

PURPOSE: Urinary incontinence and fecal incontinence are common disorders in women that negatively impact quality of life. In addition to known health and lifestyle risk factors, genetics may have a role in continence. Identification of genetic variants associated with urinary incontinence and fecal incontinence could result in a better understanding of etiologic pathways, and new interventions and treatments. MATERIALS AND METHODS: We previously generated genome-wide single nucleotide polymorphism data from Nurses' Health Studies participants. The participants provided longitudinal urinary incontinence and fecal incontinence information via questionnaires. Cases of urinary incontinence (6,120) had at least weekly urinary incontinence reported on a majority of questionnaires (3 or 4 across 12 to 16 years) while controls (4,811) consistently had little to no urinary incontinence reported. We classified cases of urinary incontinence in women into stress (1,809), urgency (1,942) and mixed (2,036) subtypes. Cases of fecal incontinence (4,247) had at least monthly fecal incontinence reported on a majority of questionnaires while controls (11,634) consistently had no fecal incontinence reported. We performed a genome-wide association study for each incontinence outcome. RESULTS: We identified 8 single nucleotide polymorphisms significantly associated (p <5×10-8) with urinary incontinence located in 2 loci, chromosomes 8q23.3 and 1p32.2. There were no genome-wide significant findings for the urinary incontinence subtype analyses. However, the significant associations for overall urinary incontinence were stronger for the urgency and mixed subtypes than for stress. While no single nucleotide polymorphism reached genome-wide significance for fecal incontinence, 4 single nucleotide polymorphisms had p <10-6. CONCLUSIONS: Few studies have collected genetic data and detailed urinary incontinence and fecal incontinence information. This genome-wide association study provides initial evidence of genetic associations for urinary incontinence and merits further research to replicate our findings and identify additional risk variants.

Low dose environmental radon exposure and breast tumor gene expression.

BACKGROUND: The International Agency for Research on Cancer classified radon and its decay-products as Group-1-human-carcinogens, and with the current knowledge they are linked specifically to lung cancer. Biokinetic models predict that radon could deliver a carcinogenic dose to breast tissue. Our previous work suggested that low-dose radon was associated with estrogen-receptor (ER)-negative breast cancer risk. However, there is limited research to examine the role of radon in breast cancer biology at the tissue level. We aim to understand molecular pathways linking radon exposure with breast cancer biology using transcriptome-wide-gene-expression from breast tumor and normal-adjacent tissues. METHODS: Our study included 943 women diagnosed with breast cancer from the Nurses' Health Study (NHS) and NHSII. We estimated cumulative radon concentration for each participant up-to the year of breast cancer diagnosis by linking residential addresses with a radon exposure model. Transcriptome-wide-gene-expression was measured with the Affymetrix-Glue-Human-Transcriptome-Array-3.0 and Human-Transcriptome-Array-2.0. We performed covariate-adjusted linear-regression for individual genes and further employed pathway-analysis. All analyses were conducted separately for tumor and normal-adjacent samples and by ER-status. RESULTS: No individual gene was associated with cumulative radon exposure in ER-positive tumor, ER-negative tumor, or ER-negative normal-adjacent tissues at FDR < 5%. In ER-positive normal-adjacent samples, PLCH2-reached transcriptome-wide-significance (FDR < 5%). Gene-set-enrichment-analyses identified 2-upregulated pathways (MAPK signaling and phosphocholine biosynthesis) enriched at FDR < 25% in ER-negative tumors and normal-adjacent tissues, and both pathways have been previously reported to play key roles in ionizing radiation induced tumorigenesis in experimental settings. CONCLUSION: Our findings provide insights into the molecular pathways of radon exposure that may influence breast cancer etiology.

Involvement of fine particulate matter exposure with gene expression pathways in breast tumor and adjacent-normal breast tissue.

BACKGROUND: Fine particulate matter (PM2.5) has been associated with breast cancer specific mortality, particularly for women with Stage I cancer. We examined the biological pathways that are perturbed by PM2.5 exposures by analyzing gene expression measurements from breast tissue specimens. METHODS: The Nurses' Health Studies (NHS and NHSII) are prospective cohorts with archival breast tissue specimens from breast cancer cases. Global gene expression data were ascertained with the Affymetrix Glue Human Transcriptome Array 3.0. PM2.5 was estimated using spatio-temporal models linked to participants' home addresses. All analyses were performed separately in tumor (n = 591) and adjacent-normal (n = 497) samples, and stratified by estrogen receptor (ER) status and stage. We used multivariable linear regression, gene-set enrichment analyses (GSEA), and the least squares kernel machine (LSKM) to assess whether 3-year cumulative average pre-diagnosis PM2.5 exposure was associated with breast-tissue gene expression pathways among predominately Stage I and II women (90.7%) and postmenopausal (81.2%) women. Replication samples (tumor, n = 245; adjacent-normal, n = 165) were measured on Affymetrix Human Transcriptome Array (HTA 2.0). RESULTS: Overall, no pathways in the tumor area were significantly associated with PM2.5 exposure. Among 272 adjacent-normal samples from Stage I ER-positive women, PM2.5 was associated with perturbations in the oxidative phosphorylation, protein secretion, and mTORC1 signaling pathways (GSEA and LSKM p-values <0.05); however, results were not replicated in a small set of replication samples (n = 80). CONCLUSIONS: PM2.5 was generally not associated with breast tissue gene expression though was suggested to perturb oxidative phosphorylation and regulation of proteins and cellular signaling in adjacent-normal breast tissue. More research is needed on the biological role of PM2.5 that influences breast tumor progression.