Gene Target Prediction of Environmental Chemicals Using Coupled Matrix-Matrix Completion.

Human exposure to toxic chemicals presents a huge health burden. Key to understanding chemical toxicity is knowledge of the molecular target(s) of the chemicals. Because a comprehensive safety assessment for all chemicals is infeasible due to limited resources, a robust computational method for discovering targets of environmental exposures is a promising direction for public health research. In this study, we implemented a novel matrix completion algorithm named coupled matrix-matrix completion (CMMC) for predicting direct and indirect exposome-target interactions, which exploits the vast amount of accumulated data regarding chemical exposures and their molecular targets. Our approach achieved an AUC of 0.89 on a benchmark data set generated using data from the Comparative Toxicogenomics Database. Our case studies with bisphenol A and its analogues, PFAS, dioxins, PCBs, and VOCs show that CMMC can be used to accurately predict molecular targets of novel chemicals without any prior bioactivity knowledge. Our results demonstrate the feasibility and promise of computationally predicting environmental chemical-target interactions to efficiently prioritize chemicals in hazard identification and risk assessment.

Environmental chemical-wide associations with immune biomarkers in US adults: A cross-sectional analysis.

Environmental chemical exposures influence immune system functions, and humans are exposed to a wide range of chemicals, termed the chemical "exposome". A comprehensive, discovery analysis of the associations of multiple chemical families with immune biomarkers is needed. In this study, we tested the associations between environmental chemical concentrations and immune biomarkers. We analyzed the United States cross-sectional National Health and Nutrition Examination Survey (NHANES, 1999-2018). Chemical biomarker concentrations were measured in blood or urine (196 chemicals, 17 chemical families). Immune biomarkers included counts of lymphocytes, neutrophils, monocytes, basophils, eosinophils, red blood cells, white blood cells, and mean corpuscular volume. We conducted separate survey-weighted, multivariable linear regressions of each log2-transformed chemical and immune measure, adjusted for relevant covariates. We accounted for multiple comparisons using a false discovery rate (FDR). Among 45,528 adult participants, the mean age was 45.7 years, 51.4% were female, and 69.3% were Non-Hispanic White. 71 (36.2%) chemicals were associated with at least one of the eight immune biomarkers. The most chemical associations (FDR<0.05) were observed with mean corpuscular volume (36 chemicals) and red blood cell counts (35 chemicals). For example, a doubling in the concentration of cotinine was associated with 0.16 fL (95% CI: 0.15, 0.17; FDR<0.001) increased mean corpuscular volume, and a doubling in the concentration of blood lead was associated with 61,736 increased red blood cells per µL (95% CI: 54,335, 69,138; FDR<0.001). A wide variety of chemicals, such as metals and smoking-related compounds, were highly associated with immune system biomarkers. This environmental chemical-wide association study identified chemicals from multiple families for further toxicological, immunologic, and epidemiological investigation.

Human Norovirus Efficiently Replicates in Differentiated 3D-Human Intestinal Enteroids.

Human norovirus (HNoV) accounts for one-fifth of all acute viral gastroenteritis worldwide and an economic burden of ~$60 billion globally. The lack of treatment options against HNoV is in part due to the lack of cultivation systems. Recently, a model of infection in biopsy-derived human intestinal enteroids (HIE) has been described: 3D-HIE are first dispersed in 2D-monolayers and differentiated prior to infection, resulting in a labor-intensive, time-consuming procedure. Here, we present an alternative protocol for HNoV infection of 3D-HIE. We found that 3D-HIE differentiated as efficiently as 2D-monolayers. In addition, immunofluorescence-based quantification of UEA-1, a lectin that stains the villus brush border, revealed that ~80% of differentiated 3D-HIE spontaneously undergo polarity inversion, allowing for viral infection without the need for microinjection. Infection with HNoV GII.4-positive stool samples attained a fold-increase over inoculum of ~2 Log10 at 2 days postinfection or up to 3.5 Log10 when ruxolitinib, a JAK1/2-inhibitor, was added. Treatment of GII.4-infected 3D-HIE with the polymerase inhibitor 2'-C-Methylcytidine (2CMC) and other antivirals showed a reduction in viral infection, suggesting that 3D-HIE are an excellent platform to test anti-infectives. The transcriptional host response to HNoV was then investigated by RNA sequencing in infected versus uninfected 3D-HIE in the presence of ruxolitinib to focus on virus-associated signatures while limiting interferon-stimulated gene signatures. The analysis revealed upregulated hormone and neurotransmitter signal transduction pathways and downregulated glycolysis and hypoxia-response pathways upon HNoV infection. Overall, 3D-HIE have proven to be a highly robust model to study HNoV infection, screen antivirals, and to investigate the host response to HNoV infection. IMPORTANCE The human norovirus (HNoV) clinical and socio-economic impact calls for immediate action in the development of anti-infectives. Physiologically relevant in vitro models are hence needed to study HNoV biology, tropism, and mechanisms of viral-associated disease, and also as a platform to identify antiviral agents. Biopsy-derived human intestinal enteroids are a biomimetic of the intestinal epithelium and were recently described as a model that supports HNoV infection. However, the established protocol is time-consuming and labor-intensive. Therefore, we sought to develop a simplified and robust alternative model of infection in 3D enteroids that undergoes differentiation and spontaneous polarity inversion. Advantages of this model are the shorter experimental time, better infection yield, and spatial integrity of the intestinal epithelium. This model is potentially suitable for the study of other pathogens that infect intestinal cells from the apical surface but also for unraveling the interactions between intestinal epithelium and indigenous bacteria of the human microbiome.

Identification, isolation and characterization of human LGR5-positive colon adenoma cells.

The intestine is maintained by stem cells located at the base of crypts and distinguished by the expression of LGR5. Genetically engineered mouse models have provided a wealth of information about intestinal stem cells, whereas less is known about human intestinal stem cells owing to difficulty detecting and isolating these cells. We established an organoid repository from patient-derived adenomas, adenocarcinomas and normal colon, which we analyzed for variants in 71 colorectal cancer (CRC)-associated genes. Normal and neoplastic colon tissue organoids were analyzed by immunohistochemistry and fluorescent-activated cell sorting for LGR5. LGR5-positive cells were isolated from four adenoma organoid lines and were subjected to RNA sequencing. We found that LGR5 expression in the epithelium and stroma was associated with tumor stage, and by integrating functional experiments with LGR5-sorted cell RNA sequencing data from adenoma and normal organoids, we found correlations between LGR5 and CRC-specific genes, including dickkopf WNT signaling pathway inhibitor 4 (DKK4) and SPARC-related modular calcium binding 2 (SMOC2). Collectively, this work provides resources, methods and new markers to isolate and study stem cells in human tissue homeostasis and carcinogenesis.

Transcriptional profiling of the response to the trichloroethylene metabolite S-(1,2-dichlorovinyl)-L-cysteine revealed activation of the eIF2α/ATF4 integrated stress response in two in vitro placental models.

Trichloroethylene (TCE) is an industrial solvent and widespread environmental contaminant. Although TCE exposure is prevalent, epidemiological studies of TCE exposure associations with adverse birth outcomes are inconclusive. Prior studies show that the TCE metabolite S-(1,2-dichlorovinyl)-L-cysteine (DCVC) exhibits toxicity in a placental cell line. In the current study, genome-wide gene expression and gene set enrichment analyses were used to identify novel genes and pathway alterations in the HTR-8/SVneo human trophoblast cell line and human placental villous explants treated with DCVC at concentrations relevant to human exposures. In the cells, concentration- and time-dependent effects were observed, as evidenced by the magnitude of altered gene expression after treatment with 20 µM DCVC versus 10 µM, and 12-h versus 6-h of treatment. Comparing the two models for the transcriptional response to 12-h 20 µM DCVC treatment, no differentially expressed genes reached significance in villous explants, whereas 301 differentially expressed genes were detected in HTR-8/SVneo cells compared with non-treated controls (FDR < 0.05 + LogFC > 0.35 [FC > 1.3]). GSEA revealed five upregulated enriched pathways in common between explants and cells (FDR < 0.05). Moreover, all 12-h DCVC treatment groups from both models contained upregulated pathways enriched for genes regulated by the ATF4 transcription factor. The overrepresentation of ATF4 regulation of differentially expressed genes indicated activation of the integrated stress response (ISR), a condition triggered by multiple stress stimuli, including the unfolded protein response. DCVC-induced ISR activation was confirmed by elevated eIF2α phosphorylation, ATF4 protein concentrations, and decreased global protein synthesis in HTR-8/SVneo cells. This study identifies a mechanism of DCVC-induced cytotoxicity by revealing the involvement of a specific stress signaling pathway.

Neutrophils Restrict Tumor-Associated Microbiota to Reduce Growth and Invasion of Colon Tumors in Mice.

BACKGROUND & AIMS: Neutrophils are among the most prevalent immune cells in the microenvironment of colon tumors; they are believed to promote growth of colon tumors, and their numbers correlate with outcomes of patients with colon cancer. Trials of inhibitors of neutrophil trafficking are underway in patients with cancer, but it is not clear how neutrophils contribute to colon tumorigenesis. METHODS: Colitis-associated colon cancer was induced in mice with conditional deletion of neutrophils (LysMCre;Mcl1fl/fl) and wild-type littermates (LysMCre;Mcl1wt/wt, control mice) by administration of azoxythmethane and/or dextran sulfate sodium. Sporadic colon tumorigenesis was assessed in neutrophil-deficient and neutrophil-replete mice with conditional deletion of colon epithelial Apc (Cdx2-CreERT2;Apcfl/fl). Primary colon tumor tissues from these mice were assessed by histology, RNA sequencing, quantitative polymerase chain reaction, and fluorescence in situ hybridization analyses. Fecal and tumor-associated microbiota were assessed by 16s ribosomal RNA sequencing. RESULTS: In mice with inflammation-induced and sporadic colon tumors, depletion of neutrophils increased the growth, proliferation, and invasiveness of the tumors. RNA sequencing analysis identified genes that regulate antimicrobial and inflammatory processes that were dysregulated in neutrophil-deficient colon tumors compared with colon tumors from control mice. Neutrophil depletion correlated with increased numbers of bacteria in tumors and proliferation of tumor cells, tumor-cell DNA damage, and an inflammatory response mediated by interleukin 17 (IL17). The 16s ribosomal RNA sequencing identified significant differences in the composition of the microbiota between colon tumors from neutrophil-deficient vs control mice. Administration of antibiotics or a neutralizing antibody against IL17 to neutrophil-deficient mice resulted in development of less-invasive tumors compared with mice given vehicle. We found bacteria in tumors to induce production of IL17, which promotes influx of intratumor B cells that promote tumor growth and progression. CONCLUSIONS: In comparisons of mice with vs without neutrophils, we found neutrophils to slow colon tumor growth and progression by restricting numbers of bacteria and tumor-associated inflammatory responses.

Temporally distinct transcriptional regulation of myocyte dedifferentiation and Myofiber growth during muscle regeneration.

BACKGROUND: Tissue regeneration requires a series of steps, beginning with generation of the necessary cell mass, followed by cell migration into damaged area, and ending with differentiation and integration with surrounding tissues. Temporal regulation of these steps lies at the heart of the regenerative process, yet its basis is not well understood. The ability of zebrafish to dedifferentiate mature "post-mitotic" myocytes into proliferating myoblasts that in turn regenerate lost muscle tissue provides an opportunity to probe the molecular mechanisms of regeneration. RESULTS: Following subtotal excision of adult zebrafish lateral rectus muscle, dedifferentiating residual myocytes were collected at two time points prior to cell cycle reentry and compared to uninjured muscles using RNA-seq. Functional annotation (GAGE or K-means clustering followed by GO enrichment) revealed a coordinated response encompassing epigenetic regulation of transcription, RNA processing, and DNA replication and repair, along with protein degradation and translation that would rewire the cellular proteome and metabolome. Selected candidate genes were phenotypically validated in vivo by morpholino knockdown. Rapidly induced gene products, such as the Polycomb group factors Ezh2 and Suz12a, were necessary for both efficient dedifferentiation (i.e. cell reprogramming leading to cell cycle reentry) and complete anatomic regeneration. In contrast, the late activated gene fibronectin was important for efficient anatomic muscle regeneration but not for the early step of myocyte cell cycle reentry. CONCLUSIONS: Reprogramming of a "post-mitotic" myocyte into a dedifferentiated myoblast requires a complex coordinated effort that reshapes the cellular proteome and rewires metabolic pathways mediated by heritable yet nuanced epigenetic alterations and molecular switches, including transcription factors and non-coding RNAs. Our studies show that temporal regulation of gene expression is programmatically linked to distinct steps in the regeneration process, with immediate early expression driving dedifferentiation and reprogramming, and later expression facilitating anatomical regeneration.

Transcriptomic profiling of curcumin-treated human breast stem cells identifies a role for stearoyl-coa desaturase in breast cancer prevention.

Curcumin is a potential agent for both the prevention and treatment of cancers. Curcumin treatment alone, or in combination with piperine, limits breast stem cell self-renewal, while remaining non-toxic to normal differentiated cells. We paired fluorescence-activated cell sorting with RNA sequencing to characterize the genome-wide changes induced specifically in normal breast stem cells following treatment with these compounds. We generated genome-wide maps of the transcriptional changes that occur in epithelial-like (ALDH+) and mesenchymal-like (ALDH-/CD44+/CD24-) normal breast stem/progenitor cells following treatment with curcumin and piperine. We show that curcumin targets both stem cell populations by down-regulating expression of breast stem cell genes including ALDH1A3, CD49f, PROM1, and TP63. We also identified novel genes and pathways targeted by curcumin, including downregulation of SCD. Transient siRNA knockdown of SCD in MCF10A cells significantly inhibited mammosphere formation and the mean proportion of CD44+/CD24- cells, suggesting that SCD is a regulator of breast stemness and a target of curcumin in breast stem cells. These findings extend previous reports of curcumin targeting stem cells, here in two phenotypically distinct stem/progenitor populations isolated from normal human breast tissue. We identified novel mechanisms by which curcumin and piperine target breast stem cell self-renewal, such as by targeting lipid metabolism, providing a mechanistic link between curcumin treatment and stem cell self-renewal. These results elucidate the mechanisms by which curcumin may act as a cancer-preventive compound and provide novel targets for cancer prevention and treatment.

NDN and CD1A are novel prognostic methylation markers in patients with head and neck squamous carcinomas.

BACKGROUND: HPV-associated HNSCCs have a distinct etiologic mechanism and better prognosis than those with non-HPV associated HNSCCs. However, even within the each group, there is heterogeneity in survival time. Here, we test the hypothesis that specific candidate gene methylation markers (CCNA1, NDN, CD1A, DCC, p16, GADD45A) are associated with tumor recurrence and survival, in a well-characterized, prospective, cohort of 346 HNSCC patients. METHODS: Kaplan-Meier curves were used to estimate survival time distributions. Multivariable Cox Proportional Hazards models were used to test associations between each methylation marker and OST/RPFT after adjusting for known or identified prognostic factors. Stratified Cox models included an interaction term between HPV and methylation marker to test for differences in the associations of the biomarker with OST or RPFT across HPV status. RESULTS: Methylation markers were differentially associated with patient characteristics. DNA hypermethylation of NDN and CD1A was found to be significantly associated with overall survival time (OST) in all HNSCC patients (NDN hazard ratio (HR): 2.35, 95% CI: 1.40-3.94; CD1A HR: 1.31, 95% CI: 1.01-1.71). Stratification by HPV status revealed hypermethylation of CD1A was associated with better OST and recurrence/persistence-free time (RPFT) (OST HR: 3.34, 95% CI: 1.88-5.93; RPFT HR: 2.06, 95% CI: 1.21-3.49), while hypomethylation of CCNA1 was associated with increased RPFT in HPV (+) patients only (HR: 0.31, 95% CI: 0.13-0.74). CONCLUSIONS: This study is the first to describe novel epigenetic alterations associated with survival in an unselected, prospectively collected, consecutive cohort of patients with HNSCC. DNA hypermethylation of NDN and CD1A was found to be significantly associated with increased overall survival time in all HNSCC patients. However, stratification by the important prognostic factor of HPV status revealed the immune marker, CD1A, and the cell cycle regulator, CCNA1 to be associated with prognosis in HPV (+) patients, specifically. Here, we identified novel methylation markers and specific, epigenetic molecular differences associated with HPV status, which warrant further investigation.

Polybrominated diphenyl ethers, 2,2',4,4',5,5'-hexachlorobiphenyl (PCB-153), and p,p'-dichlorodiphenyldichloroethylene (p,p'-DDE) concentrations in sera collected in 2009 from Texas children.

Polybrominated diphenyl ethers (PBDEs), polychlorinated biphenyls (PCBs) and p,p'-dichlorodiphenyldichloroethylene (p,p'-DDE) have been measured in surplus serum collected in 2009 from a convenience sample of 300 Texas children (boys and girls) in the birth to 13 years of age range. Serum concentrations of traditional persistent organic pollutants such as 2,2',4,4',5,5'-hexachlorobiphenyl (PCB-153) and p,p'-DDE did not change consistently with age. By contrast, serum concentrations of tetra-, penta-, and hexa-BDEs were lowest in the youngest children (birth to two year old) and increased 3.0 to 7.9 times, depending on the analyte, for children in the >4 to 6 years of age group. From the apex concentration to the 10 to 13 years of age group, concentrations decreased significantly except for 2,2',4,4',5,5'-hexabromodiphenyl ether (PBDE-153), which also had a longer apex concentration of >4 to 8 years of age. This concentration trend for PBDE-153 is most likely due to a longer half-life of PBDE-153 than of other PBDE congeners. The observed PBDEs concentration patterns by age may be related, at least in part, to ingestion of residential dust containing PBDEs through hand-to-mouth behavior among toddlers, preschoolers, and kindergarteners. Further studies to characterize young children's exposure to PBDEs are warranted and, in particular, to determine the lifestyle factors that may contribute to such exposures.

Dietary antioxidant and anti-inflammatory intake modifies the effect of cadmium exposure on markers of systemic inflammation and oxidative stress.

Chronic cadmium exposure may cause disease through induction of systemic oxidative stress and inflammation. Factors that mitigate cadmium toxicity and could serve as interventions in exposed populations have not been well characterized. We used data from the 2003-2010 National Health and Nutrition Examination Survey to quantify diet׳s role in modifying associations between cadmium exposure and oxidative stress and inflammation. We created a composite antioxidant and anti-inflammatory diet score (ADS) by ranking participants by quintile of intake across a panel of 19 nutrients. We identified associations and effect modification between ADS, urinary cadmium, and markers of oxidative stress and inflammation by multiple linear regression. An interquartile range increase in urinary cadmium was associated with a 47.5%, 8.8%, and 3.7% increase in C-reactive protein (CRP), gamma glutamyl transferase (GGT), and alkaline phosphatase (ALP), respectively. An interquartile range increase in ADS was associated with an 7.4%, 3.3%, 5.2%, and 2.5% decrease in CRP, GGT, ALP, and total white blood cell count respectively, and a 3.0% increase in serum bilirubin. ADS significantly attenuated the association between cadmium exposure, CRP and ALP. Dietary interventions may provide a route to reduce the impact of cadmium toxicity on the population level.

Longitudinal Position and Cancer Risk in the United States Revisited.

The debate over daylight saving time (DST) has surged, with interests in the effects of sunlight exposure on health. Prior studies simulated DST and standard time conditions by analyzing different locations within time zones and neighboring areas across time zone borders. We analyzed cancer incidence rates from various longitudinal positions within time zones and at time zone borders in the contiguous United States. Using data from State Cancer Profiles (2016-2020), we analyzed total cancer of 19 types and specific rates for eight cancers, adjusted for age and includes all demographics. log-linear regression is used to replicate a previous study, and spatial regression models are employed to explore discontinuities at borders. Cancer rate differences lack statistical significance within time zones and near borders for total cancer and most individual cancers. Exceptions included breast, prostate, and liver and bile duct cancers, which exhibited significant relationships with relative position at the 95% significance level. Breast and liver and bile duct cancers saw decreases, while prostate cancer incidence increased from west to east within time zones. Relative position does not have a significant impact on cancer incidence, hence cancer development in general. Isolated exceptions may warrant further investigation as more data become available. Our findings challenge prior research, revealing numerous inconsistencies. These disparities urge a reconsideration of the potential disparities in human health associated with DST and standard time. They offer insights contribute to the ongoing discussion surrounding the retention or abandonment of DST. SIGNIFICANCE: In this article, we investigate the relation between the epidemiology of cancer incidence in the United States and time zone-related longitudinal positions. Our results differ from previous research, which were based on a subset of our data, and show that the time zone effect on cancer incidence rate is not significant. Our research provides implications on the implementation of DST by suggesting that there is no cancer-risk associated reason to prefer one time over the other. Our study also uses regression discontinuity design using natural splines, a more advanced statistical method, to increase robustness of our result. Our findings challenge prior research, revealing numerous inconsistencies. These disparities urge a reconsideration of the potential disparities in human health associated with DST and standard time. They offer insights contribute to the ongoing discussion surrounding the retention or abandonment of DST.

Investigating phenotypic plasticity due to toxicants with exposure disparities in primary human breast cells in vitro.

Introduction: Breast cancer is the second most diagnosed cancer, as well as the primary cause of cancer death in women worldwide. Of the different breast cancer subtypes, triple-negative breast cancer (TNBC) is particularly aggressive and is associated with poor prognosis. Black women are two to three times more likely to be diagnosed with TNBCs than white women. Recent experimental evidence suggests that basal-like TNBCs may derive from luminal cells which acquire basal characteristics through phenotypic plasticity, a newly recognized hallmark of cancer. Whether chemical exposures can promote phenotypic plasticity in breast cells is poorly understood. Methods: To investigate further, we developed a high-content immunocytochemistry assay using normal human breast cells to test whether chemical exposures can impact luminal/basal plasticity by unbiased quantification of keratin 14 (KRT14), a basal-myoepithelial marker; keratin 8 (KRT8), a luminal-epithelial marker; and Hoechst 33342, a DNA marker. Six cell lines established from healthy tissue from donors to the Susan G. Komen Normal Tissue Bank were exposed for 48 hours to three different concentrations (0.1µM, 1µM, and 10µM) of eight ubiquitous chemicals (arsenic, BPA, BPS, cadmium, copper, DDE, lead, and PFNA), with documented exposure disparities in US Black women, in triplicate. Automated fluorescence image quantification was performed using Cell Profiler software, and a random-forest classifier was trained to classify individual cells as KRT8 positive, KRT14 positive, or hybrid (both KRT8 and KRT14 positive) using Cell Profiler Analyst. Results and discussion: Results demonstrated significant concentration-dependent increases in hybrid populations in response to BPA, BPS, DDE, and PFNA. The increase in hybrid populations expressing both KRT14 and KRT8 is indicative of a phenotypically plastic progenitor-like population in line with known theories of carcinogenesis. Furthermore, BPA, BPS, DDE, and copper produced significant increases in cell proliferation, which could be indicative of a more malignant phenotype. These results further elucidate the relationship between chemical exposure and breast phenotypic plasticity and highlight potential environmental factors that may impact TNBC risk.

High-Throughput Transcriptomics of Nontumorigenic Breast Cells Exposed to Environmentally Relevant Chemicals.

BACKGROUND: There is a suite of chemicals, including metals, pesticides, and personal care product compounds, which are commonly detected at high levels in US Center for Disease Control's National Health and Nutrition Examination Survey (NHANES) chemical biomarker screens. Whether these chemicals influence development of breast cancer is not well understood. OBJECTIVES: The objectives were to perform an unbiased concentration-dependent assessment of these chemicals, to quantify differences in cancer-specific genes and pathways, to describe if these differences occur at human population-relevant concentrations, and to specifically test for differences in markers of stemness and cellular plasticity. METHODS: We treated nontumorigenic mammary epithelial cells, MCF10A, with 21 chemicals at four concentrations (25 nM, 250 nM, 2.5µM, and 25µM) for 48 h. We conducted RNA-sequencing for these 408 samples, adapting the plexWell plate-based RNA-sequencing method to analyze differences in gene expression. We calculated gene and biological pathway-specific benchmark concentrations (BMCs) using BMDExpress3, identifying differentially expressed genes and generating the best fit benchmark concentration models for each chemical across all genes. We identified enriched biological processes and pathways for each chemical and tested whether chemical exposures change predicted cell type distributions. We contextualized benchmark concentrations relative to human population biomarker concentrations in NHANES. RESULTS: We detected chemical concentration-dependent differences in gene expression for thousands of genes. Enrichment and cell type distribution analyses showed benchmark concentration responses correlated with differences in breast cancer-related pathways, including induction of basal-like characteristics for some chemicals, including arsenic, lead, copper, and methyl paraben. Comparison of benchmark data to NHANES chemical biomarker (urine or blood) concentrations indicated an overlap between exposure levels and levels sufficient to cause a gene expression response. DISCUSSION: These analyses revealed that many of these 21 chemicals resulted in differences in genes and pathways involved in breast cancer in vitro at human exposure-relevant concentrations. https://doi.org/10.1289/EHP12886.

Sex-Specific Deflection of Age-Related DNA Methylation and Gene Expression in Mouse Heart by Perinatal Toxicant Exposures.

Background: Global and site-specific changes in DNA methylation and gene expression are associated with cardiovascular aging and disease, but how toxicant exposures during early development influence the normal trajectory of these age-related molecular changes, and whether there are sex differences, has not yet been investigated. Objectives: We used an established mouse model of developmental exposures to investigate the effects of perinatal exposure to either lead (Pb) or diethylhexyl phthalate (DEHP), two ubiquitous environmental contaminants strongly associated with CVD, on age-related cardiac DNA methylation and gene expression. Methods: Dams were randomly assigned to receive human physiologically relevant levels of Pb (32 ppm in water), DEHP (25 mg/kg chow), or control water and chow. Exposures started two weeks prior to mating and continued until weaning at postnatal day 21 (3 weeks of age). Approximately one male and one female offspring per litter were followed to 3 weeks, 5 months, or 10 months of age, at which time whole hearts were collected (n ≥ 5 per sex per exposure). Enhanced reduced representation bisulfite sequencing (ERRBS) was used to assess the cardiac DNA methylome at 3 weeks and 10 months, and RNA-seq was conducted at all 3 time points. MethylSig and edgeR were used to identify age-related differentially methylated regions (DMRs) and differentially expressed genes (DEGs), respectively, within each sex and exposure group. Cell type deconvolution of bulk RNA-seq data was conducted using the MuSiC algorithm and publicly available single cell RNA-seq data. Results: Thousands of DMRs and hundreds of DEGs were identified in control, DEHP, and Pb-exposed hearts across time between 3 weeks and 10 months of age. A closer look at the genes and pathways showing differential DNA methylation revealed that the majority were unique to each sex and exposure group. Overall, pathways governing development and differentiation were most frequently altered with age in all conditions. A small number of genes in each group showed significant changes in DNA methylation and gene expression with age, including several that were altered by both toxicants but were unchanged in control. We also observed subtle, but significant changes in the proportion of several cell types due to age, sex, and developmental exposure. Discussion: Together these data show that perinatal Pb or DEHP exposures deflect normal age-related gene expression, DNA methylation programs, and cellular composition across the life course, long after cessation of exposure, and highlight potential biomarkers of developmental toxicant exposures. Further studies are needed to investigate how these epigenetic and transcriptional changes impact cardiovascular health across the life course.

Applying Cell Painting in Non-Tumorigenic Breast Cells to Understand Impacts of Common Chemical Exposures.

There are a substantial number of chemicals to which individuals in the general population are exposed which have putative, but still poorly understood, links to breast cancer. Cell Painting is a high-content imaging-based in vitro assay that allows for rapid and unbiased measurements of the concentration-dependent effects of chemical exposures on cellular morphology. We optimized the Cell Painting assay and measured the effect of exposure to 16 human exposure relevant chemicals, along with 21 small molecules with known mechanisms of action, for 48 hours in non-tumorigenic mammary epithelial cells, the MCF10A cell line. Through unbiased imaging analyses using CellProfiler, we quantified 3042 morphological features across approximately 1.2 million cells. We used benchmark concentration modeling to quantify significance and dose-dependent directionality to identify morphological features conserved across chemicals and find features that differentiate the effects of toxicants from one another. Benchmark concentrations were compared to chemical exposure biomarker concentration measurements from the National Health and Nutrition Examination Survey to assess which chemicals induce morphological alterations at human-relevant concentrations. Morphometric fingerprint analysis revealed similar phenotypes between small molecules and prioritized NHANES-toxicants guiding further investigation. A comparison of feature fingerprints via hypergeometric analysis revealed significant feature overlaps between chemicals when stratified by compartment and stain. One such example was the similarities between a metabolite of the organochlorine pesticide DDT (p,p'-DDE) and an activator of canonical Wnt signaling CHIR99201. As CHIR99201 is a known Wnt pathway activator and its role in ß-catenin translocation is well studied, we studied the translocation of ß-catenin following p'-p' DDE exposure in an orthogonal high-content imaging assay. Consistent with activation of Wnt signaling, low dose p',p'-DDE (25nM) significantly enhances the nuclear translocation of ß-catenin. Overall, these findings highlight the ability of Cell Painting to enhance mode-of-action studies for toxicants which are common exposures in our environment but have previously been incompletely characterized with respect to breast cancer risk.

Placental and Immune Cell DNA Methylation Reference Panel for Bulk Tissue Cell Composition Estimation in Epidemiological Studies.

To distinguish DNA methylation (DNAm) from cell proportion changes in whole placental tissue research, we developed a robust cell type-specific DNAm reference to estimate cell composition. We collated newly collected and existing cell type DNAm profiles quantified via Illumina EPIC or 450k microarrays. To estimate cell composition, we deconvoluted whole placental samples (n=36) with robust partial correlation based on the top 50 hyper- and hypomethylated sites per cell type. To test deconvolution performance, we evaluated RMSE in predicting principal component one of DNAm variation in 204 external placental samples. We analyzed DNAm profiles (n=368,435 sites) from 12 cell types: cytotrophoblasts (n=18), endothelial cells (n=19), Hofbauer cells (n=26), stromal cells (n=21), syncytiotrophoblasts (n=4), six lymphocyte types (n=36), and nucleated red blood cells (n=11). Median cell composition was consistent with placental biology: 60.4% syncytiotrophoblast, 17.1% stromal, 8.8% endothelial, 4.5% cytotrophoblast, 3.9% Hofbauer, 1.7% nucleated red blood cells, and 1.2% neutrophils. Our expanded reference outperformed an existing reference in predicting DNAm variation (15.4% variance explained, IQR=21.61) with cell composition estimates (RMSE:10.51 vs. 11.43, p-value<0.001). This cell type reference can robustly estimate cell composition from whole placental DNAm data to detect important cell types, reveal biological mechanisms, and improve casual inference.

Effects of Developmental Lead and Phthalate Exposures on DNA Methylation in Adult Mouse Blood, Brain, and Liver: A Focus on Genomic Imprinting by Tissue and Sex.

BACKGROUND: Maternal exposure to environmental chemicals can cause adverse health effects in offspring. Mounting evidence supports that these effects are influenced, at least in part, by epigenetic modifications. It is unknown whether epigenetic changes in surrogate tissues such as the blood are reflective of similar changes in target tissues such as cortex or liver. OBJECTIVE: We examined tissue- and sex-specific changes in DNA methylation (DNAm) associated with human-relevant lead (Pb) and di(2-ethylhexyl) phthalate (DEHP) exposure during perinatal development in cerebral cortex, blood, and liver. METHODS: Female mice were exposed to human relevant doses of either Pb (32 ppm) via drinking water or DEHP (5mg/kg-day) via chow for 2 weeks prior to mating through offspring weaning. Whole genome bisulfite sequencing (WGBS) was utilized to examine DNAm changes in offspring cortex, blood, and liver at 5 months of age. Metilene and methylSig were used to identify differentially methylated regions (DMRs). Annotatr and ChIP-enrich were used for genomic annotations and gene set enrichment tests of DMRs, respectively. RESULTS: The cortex contained the majority of DMRs associated with Pb (66%) and DEHP (57%) exposure. The cortex also contained the greatest degree of overlap in DMR signatures between sexes (n=13 and 8 DMRs with Pb and DEHP exposure, respectively) and exposure types (n=55 and 39 DMRs in males and females, respectively). In all tissues, detected DMRs were preferentially found at genomic regions associated with gene expression regulation (e.g., CpG islands and shores, 5' UTRs, promoters, and exons). An analysis of GO terms associated with DMR-containing genes identified imprinted genes to be impacted by both Pb and DEHP exposure. Of these, Gnas and Grb10 contained DMRs across tissues, sexes, and exposures, with some signatures replicated between target and surrogate tissues. DMRs were enriched in the imprinting control regions (ICRs) of Gnas and Grb10, and we again observed a replication of DMR signatures between blood and target tissues. Specifically, we observed hypermethylation of the Grb10 ICR in both blood and liver of Pb-exposed male animals. CONCLUSIONS: These data provide preliminary evidence that imprinted genes may be viable candidates in the search for epigenetic biomarkers of toxicant exposure in target tissues. Additional research is needed on allele- and developmental stage-specific effects, as well as whether other imprinted genes provide additional examples of this relationship. https://doi.org/10.1289/EHP14074.

Bisphenol A-associated epigenomic changes in prepubescent girls: a cross-sectional study in Gharbiah, Egypt.

BACKGROUND: There is now compelling evidence that epigenetic modifications link adult disease susceptibility to environmental exposures during specific life stages, including pre-pubertal development. Animal studies indicate that bisphenol A (BPA), the monomer used in epoxy resins and polycarbonate plastics, may impact health through epigenetic mechanisms, and epidemiological data associate BPA levels with metabolic disorders, behavior changes, and reproductive effects. Thus, we conducted an environmental epidemiology study of BPA exposure and CpG methylation in pre-adolescent girls from Gharbiah, Egypt hypothesizing that methylation profiles exhibit exposure-dependent trends. METHODS: Urinary concentrations of total (free plus conjugated) species of BPA in spot samples were quantified for 60 girls aged 10 to 13. Genome-wide CpG methylation was concurrently measured in bisulfite-converted saliva DNA using the Infinium HumanMethylation27 BeadChip (N = 46). CpG sites from four candidate genes were validated via quantitative bisulfite pyrosequencing. RESULTS: CpG methylation varied widely among girls, and higher urinary BPA concentrations were generally associated with less genomic methylation. Based on pathway analyses, genes exhibiting reduced methylation with increasing urinary BPA were involved in immune function, transport activity, metabolism, and caspase activity. In particular, hypomethylation of CpG targets on chromosome X was associated with higher urinary BPA. Using the Comparative Toxicogenomics Database, we identified a number of candidate genes in our sample that previously have been associated with BPA-related expression change. CONCLUSIONS: These data indicate that BPA may affect human health through specific epigenomic modification of genes in relevant pathways. Thus, epigenetic epidemiology holds promise for the identification of biomarkers from previous exposures and the development of epigenetic-based diagnostic strategies.

Integrating NHANES and toxicity forecaster data to compare pesticide exposure and bioactivity by farmwork history and US citizenship.

BACKGROUND: Farmworkers in the United States, especially migrant workers, face unique barriers to healthcare and have documented disparities in health outcomes. Exposure to pesticides, especially those persistent in the environment, may contribute to these health disparities. OBJECTIVE: Quantify differences in pesticide exposure bioactivity by farmworker category and US citizenship status. METHODS: We queried the National Health and Nutrition Examination Study (NHANES) from 1999-2014 for pesticide exposure biomarker concentrations among farmworkers and non-farmworkers by citizenship status. We combined this with toxicity assay data from the US Environmental Protection Agency's (EPA's) Toxicity Forecaster (ToxCast). We estimated adverse biological effects that occur across a range of human population-relevant pesticide doses. RESULTS: In total, there were 844 people with any farmwork history and 23,592 non-farmworkers. Of 12 commonly detectable pesticide biomarkers in NHANES, 2,4-dichlorophenoxyacetic acid (OR = 3.76, p = 1.33 × 10-6) was significantly higher in farmworkers than non-farmworkers. Farmworkers were 1.15 times more likely to have a bioactive pesticide biomarker measurement in comparison to non-farmworkers (adjusted OR = 1.15, 95% CI: 0.87, 1.51). Non-U.S. citizens were 1.39 times more likely to have bioactive pesticide biomarker concentrations compared to people with U.S. citizenship (adjusted OR 1.39, 95% CI: 1.17, 1.64). Additionally, non-citizens were significantly more exposed to bioactive levels of ß-hexachlorocyclohexane (BHC) (OR = 8.10, p = 1.33 × 10-6), p,p-DDE (OR = 2.60, p = 0.02), and p,p'-DDT (OR = 7.75, p = 0.01). IMPACT STATEMENT: Farmworkers are a vulnerable population due to social determinants of health and occupational exposures. Here, we integrate US population chemical biomonitoring data and toxicity outcome data to assess pesticide exposure by farmwork history and citizenship. We find that farmworkers and those without US citizenship are significantly more likely to be exposed to concentrations of pesticides which are bioactive in toxicological assays. Thus, farmworkers employed in the US but who are not citizens could be at increased risk of harm to their health due to pesticides. These findings are important to shape evidence-based policies in regulatory science to promote worker safety.