Gene-Environment Interactions in Developmental Neurotoxicity: a Case Study of Synergy between Chlorpyrifos and CHD8 Knockout in Human BrainSpheres.

BACKGROUND: Autism spectrum disorder (ASD) is a major public health concern caused by complex genetic and environmental components. Mechanisms of gene-environment (G×E) interactions and reliable biomarkers associated with ASD are mostly unknown or controversial. Induced pluripotent stem cells (iPSCs) from patients or with clustered regularly interspaced short palindromic repeats and CRISPR-associated protein 9 (CRISPR/Cas9)-introduced mutations in candidate ASD genes provide an opportunity to study (G×E) interactions. OBJECTIVES: In this study, we aimed to identify a potential synergy between mutation in the high-risk autism gene encoding chromodomain helicase DNA binding protein 8 (CHD8) and environmental exposure to an organophosphate pesticide (chlorpyrifos; CPF) in an iPSC-derived human three-dimensional (3D) brain model. METHODS: This study employed human iPSC-derived 3D brain organoids (BrainSpheres) carrying a heterozygote CRISPR/Cas9-introduced inactivating mutation in CHD8 and exposed to CPF or its oxon-metabolite (CPO). Neural differentiation, viability, oxidative stress, and neurite outgrowth were assessed, and levels of main neurotransmitters and selected metabolites were validated against human data on ASD metabolic derangements. RESULTS: Expression of CHD8 protein was significantly lower in CHD8 heterozygous knockout (CHD8+/-) BrainSpheres compared with CHD8+/+ ones. Exposure to CPF/CPO treatment further reduced CHD8 protein levels, showing the potential (G×E) interaction synergy. A novel approach for validation of the model was chosen: from the literature, we identified a panel of metabolic biomarkers in patients and assessed them by targeted metabolomics in vitro. A synergistic effect was observed on the cholinergic system, S-adenosylmethionine, S-adenosylhomocysteine, lactic acid, tryptophan, kynurenic acid, and α-hydroxyglutaric acid levels. Neurite outgrowth was perturbed by CPF/CPO exposure. Heterozygous knockout of CHD8 in BrainSpheres led to an imbalance of excitatory/inhibitory neurotransmitters and lower levels of dopamine. DISCUSSION: This study pioneered (G×E) interaction in iPSC-derived organoids. The experimental strategy enables biomonitoring and environmental risk assessment for ASD. Our findings reflected some metabolic perturbations and disruption of neurotransmitter systems involved in ASD. The increased susceptibility of CHD8+/- BrainSpheres to chemical insult establishes a possibly broader role of (G×E) interaction in ASD. https://doi.org/10.1289/EHP8580.

Author Correction: Functionally Enigmatic Genes in Cancer: Using TCGA Data to Map the Limitations of Annotations.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Functionally Enigmatic Genes in Cancer: Using TCGA Data to Map the Limitations of Annotations.

Cancer is a comparatively well-studied disease, yet despite decades of intense focus, we demonstrate here using data from The Cancer Genome Atlas that a substantial number of genes implicated in cancer are relatively poorly studied. Those genes will likely be missed by any data analysis pipeline, such as enrichment analysis, that depends exclusively on annotations for understanding biological function. There is no indication that the amount of research - indicated by number of publications - is correlated with any objective metric of gene significance. Moreover, these genes are not missing at random but reflect that our information about genes is gathered in a biased manner: poorly studied genes are more likely to be primate-specific and less likely to have a Mendelian inheritance pattern, and they tend to cluster in some biological processes and not others. While this likely reflects both technological limitations as well as the fact that well-known genes tend to gather more interest from the research community, in the absence of a concerted effort to study genes in an unbiased way, many genes (and biological processes) will remain opaque.

Uncovering the Role of N-Acetyl-Aspartyl-Glutamate as a Glutamate Reservoir in Cancer.

N-acetyl-aspartyl-glutamate (NAAG) is a peptide-based neurotransmitter that has been extensively studied in many neurological diseases. In this study, we show a specific role of NAAG in cancer. We found that NAAG is more abundant in higher grade cancers and is a source of glutamate in cancers expressing glutamate carboxypeptidase II (GCPII), the enzyme that hydrolyzes NAAG to glutamate and N-acetyl-aspartate (NAA). Knocking down GCPII expression through genetic alteration or pharmacological inhibition of GCPII results in a reduction of both glutamate concentrations and cancer growth. Moreover, targeting GCPII in combination with glutaminase inhibition accentuates these effects. These findings suggest that NAAG serves as an important reservoir to provide glutamate to cancer cells through GCPII when glutamate production from other sources is limited. Thus, GCPII is a viable target for cancer therapy, either alone or in combination with glutaminase inhibition.

Weighted Gene Correlation Network Analysis (WGCNA) Reveals Novel Transcription Factors Associated With Bisphenol A Dose-Response.

Despite Bisphenol-A (BPA) being subject to extensive study, a thorough understanding of molecular mechanism remains elusive. Here we show that using weighted gene correlation network analysis (WGCNA), which takes advantage of a graph theoretical approach to understanding correlations amongst genes and grouping genes into modules that typically have co-ordinated biological functions and regulatory mechanisms, that despite some commonality in altered genes, there is minimal overlap between BPA and estrogen in terms of network topology. We confirmed previous findings that ZNF217 and TFAP2C are involved in the estrogen pathway, and are implicated in BPA as well, although for BPA they appear to be active in the absence of canonical estrogen-receptor driven gene expression. Furthermore, our study suggested that PADI4 and RACK7/ZMYNDB8 may be involved in the overlap in gene expression between estradiol and BPA. Lastly, we demonstrated that even at low doses there are unique transcription factors that appear to be driving the biology of BPA, such as SREBF1. Overall, our data is consistent with other reports that BPA leads to subtle gene changes rather than profound aberrations of a conserved estrogen signaling (or other) pathways.

Toxicity, recovery, and resilience in a 3D dopaminergic neuronal in vitro model exposed to rotenone.

To date, most in vitro toxicity testing has focused on acute effects of compounds at high concentrations. This testing strategy does not reflect real-life exposures, which might contribute to long-term disease outcome. We used a 3D-human dopaminergic in vitro LUHMES cell line model to determine whether effects of short-term rotenone exposure (100 nM, 24 h) are permanent or reversible. A decrease in complex I activity, ATP, mitochondrial diameter, and neurite outgrowth were observed acutely. After compound removal, complex I activity was still inhibited; however, ATP levels were increased, cells were electrically active and aggregates restored neurite outgrowth integrity and mitochondrial morphology. We identified significant transcriptomic changes after 24 h which were not present 7 days after wash-out. Our results suggest that testing short-term exposures in vitro may capture many acute effects which cells can overcome, missing adaptive processes, and long-term mechanisms. In addition, to study cellular resilience, cells were re-exposed to rotenone after wash-out and recovery period. Pre-exposed cells maintained higher metabolic activity than controls and presented a different expression pattern in genes previously shown to be altered by rotenone. NEF2L2, ATF4, and EAAC1 were downregulated upon single hit on day 14, but unchanged in pre-exposed aggregates. DAT and CASP3 were only altered after re-exposure to rotenone, while TYMS and MLF1IP were downregulated in both single-exposed and pre-exposed aggregates. In summary, our study shows that a human cell-based 3D model can be used to assess cellular adaptation, resilience, and long-term mechanisms relevant to neurodegenerative research.

Joint Bounding of Peaks Across Samples Improves Differential Analysis in Mass Spectrometry-Based Metabolomics.

As mass spectrometry-based metabolomics becomes more widely used in biomedical research, it is important to revisit existing data analysis paradigms. Existing data preprocessing efforts have largely focused on methods which start by extracting features separately from each sample, followed by a subsequent attempt to group features across samples to facilitate comparisons. We show that this preprocessing approach leads to unnecessary variability in peak quantifications that adversely impacts downstream analysis. We present a new method, bakedpi, for the preprocessing of both centroid and profile mode metabolomics data that relies on an intensity-weighted bivariate kernel density estimation on a pooling of all samples to detect peaks. This new method reduces this unnecessary quantification variability and increases power in downstream differential analysis.

Metabolomic network analysis of estrogen-stimulated MCF-7 cells: a comparison of overrepresentation analysis, quantitative enrichment analysis and pathway analysis versus metabolite network analysis.

In the context of the Human Toxome project, mass spectroscopy-based metabolomics characterization of estrogen-stimulated MCF-7 cells was studied in order to support the untargeted deduction of pathways of toxicity. A targeted and untargeted approach using overrepresentation analysis (ORA), quantitative enrichment analysis (QEA) and pathway analysis (PA) and a metabolite network approach were compared. Any untargeted approach necessarily has some noise in the data owing to artifacts, outliers and misidentified metabolites. Depending on the chemical analytical choices (sample extraction, chromatography, instrument and settings, etc.), only a partial representation of all metabolites will be achieved, biased by both the analytical methods and the database used to identify the metabolites. Here, we show on the one hand that using a data analysis approach based exclusively on pathway annotations has the potential to miss much that is of interest and, in the case of misidentified metabolites, can produce perturbed pathways that are statistically significant yet uninformative for the biological sample at hand. On the other hand, a targeted approach, by narrowing its focus and minimizing (but not eliminating) misidentifications, renders the likelihood of a spurious pathway much smaller, but the limited number of metabolites also makes statistical significance harder to achieve. To avoid an analysis dependent on pathways, we built a de novo network using all metabolites that were different at 24 h with and without estrogen with a p value <0.01 (53) in the STITCH database, which links metabolites based on known reactions in the main metabolic network pathways but also based on experimental evidence and text mining. The resulting network contained a "connected component" of 43 metabolites and helped identify non-endogenous metabolites as well as pathways not visible by annotation-based approaches. Moreover, the most highly connected metabolites (energy metabolites such as pyruvate and alpha-ketoglutarate, as well as amino acids) showed only a modest change between proliferation with and without estrogen. Here, we demonstrate that estrogen has subtle but potentially phenotypically important alterations in the acyl-carnitine fatty acids, acetyl-putrescine and succinoadenosine, in addition to likely subtle changes in key energy metabolites that, however, could not be verified consistently given the technical limitations of this approach. Finally, we show that a network-based approach combined with text mining identifies pathways that would otherwise neither be considered statistically significant on their own nor be identified via ORA, QEA, or PA.

Genetic variability in a frozen batch of MCF-7 cells invisible in routine authentication affecting cell function.

Common recommendations for cell line authentication, annotation and quality control fall short addressing genetic heterogeneity. Within the Human Toxome Project, we demonstrate that there can be marked cellular and phenotypic heterogeneity in a single batch of the human breast adenocarcinoma cell line MCF-7 obtained directly from a cell bank that are invisible with the usual cell authentication by short tandem repeat (STR) markers. STR profiling just fulfills the purpose of authentication testing, which is to detect significant cross-contamination and cell line misidentification. Heterogeneity needs to be examined using additional methods. This heterogeneity can have serious consequences for reproducibility of experiments as shown by morphology, estrogenic growth dose-response, whole genome gene expression and untargeted mass-spectroscopy metabolomics for MCF-7 cells. Using Comparative Genomic Hybridization (CGH), differences were traced back to genetic heterogeneity already in the cells from the original frozen vials from the same ATCC lot, however, STR markers did not differ from ATCC reference for any sample. These findings underscore the need for additional quality assurance in Good Cell Culture Practice and cell characterization, especially using other methods such as CGH to reveal possible genomic heterogeneity and genetic drifts within cell lines.

The human toxome project.

The Human Toxome Project, funded as an NIH Transformative Research grant 2011-2016, is focused on developing the concepts and the means for deducing, validating and sharing molecular pathways of toxicity (PoT). Using the test case of estrogenic endocrine disruption, the responses of MCF-7 human breast cancer cells are being phenotyped by transcriptomics and mass-spectroscopy-based metabolomics. The bioinformatics tools for PoT deduction represent a core deliverable. A number of challenges for quality and standardization of cell systems, omics technologies and bioinformatics are being addressed. In parallel, concepts for annotation, validation and sharing of PoT information, as well as their link to adverse outcomes, are being developed. A reasonably comprehensive public database of PoT, the Human Toxome Knowledge-base, could become a point of reference for toxicological research and regulatory test strategies.

Two new approaches to improve the analysis of BALB/c 3T3 cell transformation assay data.

Validation activities of the BALB/c 3T3 cell transformation assay (CTA) - a test method used for the assessment of the carcinogenic potential of compounds - have revealed the need for statistical analysis tailored to specific features of BALB/c 3T3 CTA data. Whereas a standard statistical approach for the Syrian hamster embryo (SHE) CTA was considered sufficient, an international expert group was gathered by the European Centre for the Validation of Alternative Methods (ECVAM) to review commonly applied statistical approaches for BALB/c 3T3 CTA. As it was concluded that none of the commonly applied approaches is entirely appropriate, two novel statistical approaches were found to be recommended for the evaluation of BALB/c 3T3 CTA data accounting for possible non-monotone concentration-response relationship and variance heterogeneity: a negative binomial generalised linear model with William's-type downturn-protected trend tests and a normalisation of the data by a specific transformation allowing for application of a general linear model that estimates effects assuming a normal distribution with William's-type protected tests. Both approaches are described in this article and their performance and the quality of the results they generate is demonstrated using exemplary data. Our work confirmed that both approaches are suitable for the statistical analysis of BALB/c 3T3 CTA data and that each of them is superior to commonly used methods. Furthermore, a procedure dichotomising data into negatives and positives is proposed which allows re-testing in cases where inconclusive data are encountered. The scripts of the statistical evaluation programs written in R - a freely available statistical software - are appended including exemplary outputs (Appendix A).

Human genetic resistance to Onchocerca volvulus: evidence for linkage to chromosome 2p from an autosome-wide scan.

BACKGROUND: Human infections with the tissue nematode Onchocerca volvulus show strong interindividual variation in intensity, which cannot be explained by differences in exposure alone. Several lines of evidence suggest a relevant influence of human genetics. METHODS: In a genome-wide search for genetic determinants of resistance, we studied 196 siblings from 51 families exposed to endemic O. volvulus transmission in the forest zone of Ghana, West Africa. The numbers of worm larvae in the skin (i.e., microfilariae), which are the established measure of O. volvulus infection intensity, were counted in 4 small skin biopsy specimens (i.e., skin snips), and the numbers of palpable subcutaneous worm nodules (i.e., onchocercomata) were assessed. Numbers were corrected for age and exposure and were analyzed for linkage to 377 autosomal microsatellite markers and additional markers in genomic regions of interest. RESULTS: Linkage was detected between the numbers of microfilariae and chromosome 2p21-p14 (maximum multipoint log(10) of odds (LOD) score of 3.80 at marker position D2S2378; empirical P=2.9 x 10(-5)). CONCLUSIONS: This finding provides strong evidence that a human genetic factor influences the intensity of O. volvulus infection. The strength of the linkage signal may facilitate the identification of the decisive genetic variants.