Network Analysis of Biomarkers Associated with Occupational Exposure to Benzene and Malathion.

Complex diseases are associated with the effects of multiple genes, proteins, and biological pathways. In this context, the tools of Network Medicine are compatible as a platform to systematically explore not only the molecular complexity of a specific disease but may also lead to the identification of disease modules and pathways. Such an approach enables us to gain a better understanding of how environmental chemical exposures affect the function of human cells, providing better perceptions about the mechanisms involved and helping to monitor/prevent exposure and disease to chemicals such as benzene and malathion. We selected differentially expressed genes for exposure to benzene and malathion. The construction of interaction networks was carried out using GeneMANIA and STRING. Topological properties were calculated using MCODE, BiNGO, and CentiScaPe, and a Benzene network composed of 114 genes and 2415 interactions was obtained. After topological analysis, five networks were identified. In these subnets, the most interconnected nodes were identified as: IL-8, KLF6, KLF4, JUN, SERTAD1, and MT1H. In the Malathion network, composed of 67 proteins and 134 interactions, HRAS and STAT3 were the most interconnected nodes. Path analysis, combined with various types of high-throughput data, reflects biological processes more clearly and comprehensively than analyses involving the evaluation of individual genes. We emphasize the central roles played by several important hub genes obtained by exposure to benzene and malathion.

Putative contributions of the sex chromosome proteins SOX3 and SRY to neurodevelopmental disorders.

The male-biased prevalence of certain neurodevelopmental disorders and the sex-biased outcomes associated with stress exposure during gestation have been previously described. Here, we hypothesized that genes distinctively targeted by only one or both homologous proteins highly conserved across therian mammals, SOX3 and SRY, could induce sexual adaptive changes that result in a differential risk for neurodevelopmental disorders. ChIP-seq/chip data showed that SOX3/SRY gene targets were expressed in different brain cell types in mice. We used orthologous human genes in rodent genomes to extend the number of SOX3/SRY set (1,721). These genes were later found to be enriched in five modules of coexpressed genes during the early and mid-gestation periods (FDR < 0.05), independent of sexual hormones. Genes with differential expression (24, p < 0.0001) and methylation (40, p < 0.047) between sexes were overrepresented in this set. Exclusive SOX3 or SRY target genes were more associated with the late gestational and postnatal periods. Using autism as a model sex-biased disorder, the SOX3/SRY set was enriched in autism gene databases (FDR ≤ 0.05), and there were more de novo variations from the male autism spectrum disorder (ASD) samples under the SRY peaks compared to the random peaks (p < 0.024). The comparison of coexpressed networks of SOX3/SRY target genes between male autism and control samples revealed low preservation in gene modules related to stress response (99 genes) and neurogenesis (78 genes). This study provides evidence that while SOX3 is a regulatory mechanism for both sexes, the male-exclusive SRY also plays a role in gene regulation, suggesting a potential mechanism for sex bias in ASD.

Lifespan analysis of repeat expression reveals age-dependent upregulation of HERV-K in the neurotypical human brain.

DNA repetitive sequences (or repeats) comprise over 50% of the human genome and have a crucial regulatory role, specifically regulating transcription machinery. The human brain is the tissue with the highest detectable repeat expression and dysregulations on the repeat activity are related to several neurological and neurodegenerative disorders, as repeat-derived products can stimulate a pro-inflammatory response. Even so, it is unclear how repeat expression acts on the aging neurotypical brain. Here, we leverage a large postmortem transcriptome cohort spanning the human lifespan to assess global repeat expression in the neurotypical brain. We identified 21,696 differentially expressed repeats (DERs) that varied across seven age bins (Prenatal; 0-15; 16-29; 30-39; 40-49; 50-59; 60+) across the caudate nucleus (n=271), dorsolateral prefrontal cortex (n=304), and hippocampus (n=310). Interestingly, we found that long interspersed nuclear elements and long terminal repeats (LTRs) DERs were the most abundant repeat families when comparing infants to early adolescence (0-15) with older adults (60+). Of these differentially regulated LTRs, we identified 17 shared across all brain regions, including increased expression of HERV-K-int in older adult brains (60+). Co-expression analysis from each of the three brain regions also showed repeats from the HERV subfamily were intramodular hubs in its subnetworks. While we do not observe a strong global relationship between repeat expression and age, we identified HERV-K as a repeat signature associated with the aging neurotypical brain. Our study is the first global assessment of repeat expression in the neurotypical brain.

Sex affects transcriptional associations with schizophrenia across the dorsolateral prefrontal cortex, hippocampus, and caudate nucleus.

Schizophrenia is a complex neuropsychiatric disorder with sexually dimorphic features, including differential symptomatology, drug responsiveness, and male incidence rate. Prior large-scale transcriptome analyses for sex differences in schizophrenia have focused on the prefrontal cortex. Analyzing BrainSeq Consortium data (caudate nucleus: n = 399, dorsolateral prefrontal cortex: n = 377, and hippocampus: n = 394), we identified 831 unique genes that exhibit sex differences across brain regions, enriched for immune-related pathways. We observed X-chromosome dosage reduction in the hippocampus of male individuals with schizophrenia. Our sex interaction model revealed 148 junctions dysregulated in a sex-specific manner in schizophrenia. Sex-specific schizophrenia analysis identified dozens of differentially expressed genes, notably enriched in immune-related pathways. Finally, our sex-interacting expression quantitative trait loci analysis revealed 704 unique genes, nine associated with schizophrenia risk. These findings emphasize the importance of sex-informed analysis of sexually dimorphic traits, inform personalized therapeutic strategies in schizophrenia, and highlight the need for increased female samples for schizophrenia analyses.

Human archetypal pluripotent stem cells differentiate into trophoblast stem cells via endogenous BMP5/7 induction without transitioning through naive state.

Primary human trophoblast stem cells (TSCs) and TSCs derived from human pluripotent stem cells (hPSCs) can potentially model placental processes in vitro. Yet, the pluripotent states and factors involved in the differentiation of hPSCs to TSCs remain poorly understood. In this study, we demonstrate that the primed pluripotent state can generate TSCs by activating pathways such as Epidermal Growth Factor (EGF) and Wingless-related integration site (WNT), and by suppressing tumor growth factor beta (TGFß), histone deacetylases (HDAC), and Rho-associated protein kinase (ROCK) signaling pathways, all without the addition of exogenous Bone morphogenetic protein 4 (BMP4)-a condition we refer to as the TS condition. We characterized this process using temporal single-cell RNA sequencing to compare TS conditions with differentiation protocols involving BMP4 activation alone or BMP4 activation in conjunction with WNT inhibition. The TS condition consistently produced a stable, proliferative cell type that closely mimics first-trimester placental cytotrophoblasts, marked by the activation of endogenous retroviral genes and the absence of amnion expression. This was observed across multiple cell lines, including various primed induced pluripotent stem cell (iPSC) and embryonic stem cell (ESC) lines. Primed-derived TSCs can proliferate for over 30 passages and further specify into multinucleated syncytiotrophoblasts and extravillous trophoblast cells. Our research establishes that the differentiation of primed hPSCs to TSC under TS conditions triggers the induction of TMSB4X, BMP5/7, GATA3, and TFAP2A without progressing through a naive state. These findings propose that the primed hPSC state is part of a continuum of potency with the capacity to differentiate into TSCs through multiple routes.

Recapitulation of Perturbed Striatal Gene Expression Dynamics of Donors' Brains With Ventral Forebrain Organoids Derived From the Same Individuals With Schizophrenia.

OBJECTIVE: Schizophrenia is a brain disorder that originates during neurodevelopment and has complex genetic and environmental etiologies. Despite decades of clinical evidence of altered striatal function in affected patients, studies examining its cellular and molecular mechanisms in humans are limited. To explore neurodevelopmental alterations in the striatum associated with schizophrenia, the authors established a method for the differentiation of induced pluripotent stem cells (iPSCs) into ventral forebrain organoids (VFOs). METHODS: VFOs were generated from postmortem dural fibroblast-derived iPSCs of four individuals with schizophrenia and four neurotypical control individuals for whom postmortem caudate genotypes and transcriptomic data were profiled in the BrainSeq neurogenomics consortium. Individuals were selected such that the two groups had nonoverlapping schizophrenia polygenic risk scores (PRSs). RESULTS: Single-cell RNA sequencing analyses of VFOs revealed differences in developmental trajectory between schizophrenia and control individuals in which inhibitory neuronal cells from the patients exhibited accelerated maturation. Furthermore, upregulated genes in inhibitory neurons in schizophrenia VFOs showed a significant overlap with upregulated genes in postmortem caudate tissue of individuals with schizophrenia compared with control individuals, including the donors of the iPSC cohort. CONCLUSIONS: The findings suggest that striatal neurons derived from high-PRS individuals with schizophrenia carry abnormalities that originated during early brain development and that the VFO model can recapitulate disease-relevant cell type-specific neurodevelopmental phenotypes in a dish.

Sex differences in gene regulatory networks during mid-gestational brain development.

Neurodevelopmental disorders differ considerably between males and females, and fetal brain development is one of the most critical periods to determine risk for these disorders. Transcriptomic studies comparing male and female fetal brain have demonstrated that the highest difference in gene expression occurs in sex chromosomes, but several autossomal genes also demonstrate a slight difference that has not been yet explored. In order to investigate biological pathways underlying fetal brain sex differences, we applied medicine network principles using integrative methods such as co-expression networks (CEMiTool) and regulatory networks (netZoo). The pattern of gene expression from genes in the same pathway tend to reflect biologically relevant phenomena. In this study, network analysis of fetal brain expression reveals regulatory differences between males and females. Integrating two different bioinformatics tools, our results suggest that biological processes such as cell cycle, cell differentiation, energy metabolism and extracellular matrix organization are consistently sex-biased. MSET analysis demonstrates that these differences are relevant to neurodevelopmental disorders, including autism.

Gestational age acceleration is associated with epigenetic biomarkers of prenatal physiologic stress exposure.

BACKGROUND: Physiological maternal stress response, such as imbalance in the glucocorticoid pathway and immune system seems to be mediated by DNA methylation (DNAm) and might translate intrauterine stress exposures into phenotypic changes in a sex-specific manner. DNAm in specific sites can also predict newborn gestational age and gestational age acceleration (GAA). GAA occurs when the predicted biological age is higher than the chronological age. In adults, poor health outcomes related to this deviance are well documented and raise questions for the interpretation and prediction in early stages of life. Boys seem to be more vulnerable to intrauterine stress exposure than girls; however, the mechanisms of adaptive sex-specific responses are still unclear. We hypothesize that intrauterine stress exposure is associated with GAA and could be different in boys and girls if inflammatory or glucocorticoid pathways exposure is considered. RESULTS: Using the Western Region Birth Cohort (ROC-São Paulo, Brazil) (n = 83), we calculated DNAm age and GAA from cord blood samples. Two epigenetic risk scores were calculated as an indirect proxy for low-grade inflammation (i-ePGS) and for glucocorticoid exposure (GES). Multivariate linear regression models were applied to investigate associations of GAA with prenatal exposures. The i-ePGS and GES were included in different models with the same co-variates considering sex interactions. The first multivariate model investigating inflammatory exposure (adj. R2 = 0.31, p = < 0.001) showed that GAA was positively associated with i-ePGS (CI, 0.26-113.87, p = 0.049) and negative pregnancy-related feelings (CI, 0.04-0.48 p = 0.019). No sex interaction was observed. The second model investigating glucocorticoid exposure (adj. R2 = 0.32, p = < 0.001) showed that the higher was the GAA was associated with a lower the lower was the GES in girls (CI, 0.04-2.55, p = 0.044). In both models, maternal self-reported mental disorder was negatively associated with GAA. CONCLUSION: Prenatal epigenetic score of exposure to low-grade inflammatory was a predictor of GAA for both sexes. Glucocorticoid epigenetic score seems to be more important to GAA in girls. This study supports the evidence of sex-specificity in stress response, suggesting the glucocorticoid as a possible pathway adopted by girls to accelerate the maturation in an adverse condition.

DNA methylation mediates a randomized controlled trial home-visiting intervention during pregnancy and the Bayley infant's cognitive scores at 12 months of age.

The crosstalk between maternal stress exposure and fetal development may be mediated by epigenetic mechanisms, including DNA methylation (DNAm). To address this matter, we collect 32 cord blood samples from low-income Brazilian pregnant adolescents participants of a pilot randomized clinical intervention study (ClinicalTrials.gov, Identifier: NCT02807818). We hypothesized that the association between the intervention and infant neurodevelopmental outcomes at 12 months of age would be mediated by DNAm. First, we searched genome methylation differences between cases and controls using different approaches, as well as differences in age acceleration (AA), represented by the difference of methylation age and birth age. According to an adjusted p-value ≤ 0.05 we identified 3090 differentially methylated positions- CpG sites (DMPs), 21 differentially methylated regions (DMRs) and one comethylated module weakly preserved between groups. The intervention group presented a smaller AA compared to the control group (p = 0.025). A logistic regression controlled by sex and with gestational age indicated a coefficient of -0.35 towards intervention group (p = 0.016) considering AA. A higher cognitive domain score from Bayley III scale was observed in the intervention group at 12 months of age. Then, we performed a potential causal mediation analysis selecting only DMPs highly associated with the cognitive domain (adj. R2 > 0.4), DMRs and CpGs of hub genes from the weakly preserved comethylated module and epigenetic clock as raw values. DMPs in STXBP6, and PF4 DMR, mediated the association between the maternal intervention and the cognitive domain at 12 months of age. In conclusion, DNAm in different sites and regions mediated the association between intervention and cognitive outcome.

Analysis of the caudate nucleus transcriptome in individuals with schizophrenia highlights effects of antipsychotics and new risk genes.

Most studies of gene expression in the brains of individuals with schizophrenia have focused on cortical regions, but subcortical nuclei such as the striatum are prominently implicated in the disease, and current antipsychotic drugs target the striatum's dense dopaminergic innervation. Here, we performed a comprehensive analysis of the genetic and transcriptional landscape of schizophrenia in the postmortem caudate nucleus of the striatum of 443 individuals (245 neurotypical individuals, 154 individuals with schizophrenia and 44 individuals with bipolar disorder), 210 from African and 233 from European ancestries. Integrating expression quantitative trait loci analysis, Mendelian randomization with the latest schizophrenia genome-wide association study, transcriptome-wide association study and differential expression analysis, we identified many genes associated with schizophrenia risk, including potentially the dopamine D2 receptor short isoform. We found that antipsychotic medication has an extensive influence on caudate gene expression. We constructed caudate nucleus gene expression networks that highlight interactions involving schizophrenia risk. These analyses provide a resource for the study of schizophrenia and insights into risk mechanisms and potential therapeutic targets.

Brain areas involved with obsessive-compulsive disorder present different DNA methylation modulation.

BACKGROUND: Obsessive-compulsive disorder (OCD) is characterized by intrusive thoughts and repetitive actions, that presents the involvement of the cortico-striatal areas. The contribution of environmental risk factors to OCD development suggests that epigenetic mechanisms may contribute to its pathophysiology. DNA methylation changes and gene expression were evaluated in post-mortem brain tissues of the cortical (anterior cingulate gyrus and orbitofrontal cortex) and ventral striatum (nucleus accumbens, caudate nucleus and putamen) areas from eight OCD patients and eight matched controls. RESULTS: There were no differentially methylated CpG (cytosine-phosphate-guanine) sites (DMSs) in any brain area, nevertheless gene modules generated from CpG sites and protein-protein-interaction (PPI) showed enriched gene modules for all brain areas between OCD cases and controls. All brain areas but nucleus accumbens presented a predominantly hypomethylation pattern for the differentially methylated regions (DMRs). Although there were common transcriptional factors that targeted these DMRs, their targeted differentially expressed genes were different among all brain areas. The protein-protein interaction network based on methylation and gene expression data reported that all brain areas were enriched for G-protein signaling pathway, immune response, apoptosis and synapse biological processes but each brain area also presented enrichment of specific signaling pathways. Finally, OCD patients and controls did not present significant DNA methylation age differences. CONCLUSIONS: DNA methylation changes in brain areas involved with OCD, especially those involved with genes related to synaptic plasticity and the immune system could mediate the action of genetic and environmental factors associated with OCD.

Are the 50's, the transition decade, in choroid plexus aging?

The choroid plexus (CP) is an important structure for the brain. Besides its major role in the production of cerebrospinal fluid (CSF), it conveys signals originating from the brain, and from the circulatory system, shaping brain function in health and in pathology. Previous studies in rodents have revealed altered transcriptome both during aging and in various diseases of the central nervous system, including Alzheimer's disease. In the present study, a high-throughput sequencing of the CP transcriptome was performed in postmortem samples of clinically healthy individuals aged 50's through 80's. The data shows an age-related profile, with the main changes occurring in the transition from the 50's to the 60's, stabilizing thereafter. Specifically, neuronal and membrane functions distinguish the transcriptome between the 50's and the 60's, while neuronal and axon development and extracellular structure organization differentiate the 50's from the 70's. These findings suggest that changes in the CP transcriptome occur early in the aging process. Future studies will unravel whether these relate with processes occurring in late- onset brain diseases.

Assessment of complementarity of WGCNA and NERI results for identification of modules associated to schizophrenia spectrum disorders.

Psychiatric disorders involve both changes in multiple genes as well different types of variations. As such, gene co-expression networks allowed the comparison of different stages and parts of the brain contributing to an integrated view of genetic variation. Two methods based on co-expression networks presents appealing results: Weighted Gene Correlation Network Analysis (WGCNA) and Network-Medicine Relative Importance (NERI). By selecting two different gene expression databases related to schizophrenia, we evaluated the biological modules selected by both WGCNA and NERI along these databases as well combining both WGCNA and NERI results (WGCNA-NERI). Also we conducted a enrichment analysis for the identification of partial biological function of each result (as well a replication analysis). To appraise the accuracy of whether both algorithms (as well our approach, WGCNA-NERI) were pointing to genes related to schizophrenia and its complex genetic architecture we conducted the MSET analysis, based on a reference gene list of schizophrenia database (SZDB) related to DNA Methylation, Exome, GWAS as well as copy number variation mutation studies. The WGCNA results were more associated with inflammatory pathways and immune system response; NERI obtained genes related with cellular regulation, embryological pathways e cellular growth factors. Only NERI were able to provide a statistical meaningful results to the MSET analysis (for Methylation and de novo mutations data). However, combining WGCNA and NERI provided a much more larger overlap in these two categories and additionally on Transcriptome database. Our study suggests that using both methods in combination is better for establishing a group of modules and pathways related to a complex disease than using each method individually. NERI is available at: https://bitbucket.org/sergionery/neri.

Initial findings of striatum tripartite model in OCD brain samples based on transcriptome analysis.

Obsessive-compulsive disorder (OCD) is a psychiatric disorder characterized by obsessions and/or compulsions. Different striatal subregions belonging to the cortico-striato-thalamic circuitry (CSTC) play an important role in the pathophysiology of OCD. The transcriptomes of 3 separate striatal areas (putamen (PT), caudate nucleus (CN) and accumbens nucleus (NAC)) from postmortem brain tissue were compared between 6 OCD and 8 control cases. In addition to network connectivity deregulation, different biological processes are specific to each striatum region according to the tripartite model of the striatum and contribute in various ways to OCD pathophysiology. Specifically, regulation of neurotransmitter levels and presynaptic processes involved in chemical synaptic transmission were shared between NAC and PT. The Gene Ontology terms cellular response to chemical stimulus, response to external stimulus, response to organic substance, regulation of synaptic plasticity, and modulation of synaptic transmission were shared between CN and PT. Most genes harboring common and/or rare variants previously associated with OCD that were differentially expressed or part of a least preserved coexpression module in our study also suggest striatum subregion specificity. At the transcriptional level, our study supports differences in the 3 circuit CSTC model associated with OCD.