Multiomics analysis identifies novel facilitators of human dopaminergic neuron differentiation.

Midbrain dopaminergic neurons (mDANs) control voluntary movement, cognition, and reward behavior under physiological conditions and are implicated in human diseases such as Parkinson's disease (PD). Many transcription factors (TFs) controlling human mDAN differentiation during development have been described, but much of the regulatory landscape remains undefined. Using a tyrosine hydroxylase (TH) human iPSC reporter line, we here generate time series transcriptomic and epigenomic profiles of purified mDANs during differentiation. Integrative analysis predicts novel regulators of mDAN differentiation and super-enhancers are used to identify key TFs. We find LBX1, NHLH1 and NR2F1/2 to promote mDAN differentiation and show that overexpression of either LBX1 or NHLH1 can also improve mDAN specification. A more detailed investigation of TF targets reveals that NHLH1 promotes the induction of neuronal miR-124, LBX1 regulates cholesterol biosynthesis, and NR2F1/2 controls neuronal activity.

Functional meta-omics provide critical insights into long- and short-read assemblies.

Real-world evaluations of metagenomic reconstructions are challenged by distinguishing reconstruction artifacts from genes and proteins present in situ. Here, we evaluate short-read-only, long-read-only and hybrid assembly approaches on four different metagenomic samples of varying complexity. We demonstrate how different assembly approaches affect gene and protein inference, which is particularly relevant for downstream functional analyses. For a human gut microbiome sample, we use complementary metatranscriptomic and metaproteomic data to assess the metagenomic data-based protein predictions. Our findings pave the way for critical assessments of metagenomic reconstructions. We propose a reference-independent solution, which exploits the synergistic effects of multi-omic data integration for the in situ study of microbiomes using long-read sequencing data.

RvD1n-3 DPA Downregulates the Transcription of Pro-Inflammatory Genes in Oral Epithelial Cells and Reverses Nuclear Translocation of Transcription Factor p65 after TNF-α Stimulation.

Specialized pro-resolving mediators (SPMs) are multifunctional lipid mediators that participate in the resolution of inflammation. We have recently described that oral epithelial cells (OECs) express receptors of the SPM resolvin RvD1n-3 DPA and that cultured OECs respond to RvD1n-3 DPA addition by intracellular calcium release, nuclear receptor translocation and transcription of genes coding for antimicrobial peptides. The aim of the present study was to assess the functional outcome of RvD1n-3 DPA-signaling in OECs under inflammatory conditions. To this end, we performed transcriptomic analyses of TNF-α-stimulated cells that were subsequently treated with RvD1n-3 DPA and found significant downregulation of pro-inflammatory nuclear factor kappa B (NF-κB) target genes. Further bioinformatics analyses showed that RvD1n-3 DPA inhibited the expression of several genes involved in the NF-κB activation pathway. Confocal microscopy revealed that addition of RvD1n-3 DPA to OECs reversed TNF-α-induced nuclear translocation of NF-κB p65. Co-treatment of the cells with the exportin 1 inhibitor leptomycin B indicated that RvD1n-3 DPA increases nuclear export of p65. Taken together, our observations suggest that SPMs also have the potential to be used as a therapeutic aid when inflammation is established.

Temporal enhancer profiling of parallel lineages identifies AHR and GLIS1 as regulators of mesenchymal multipotency.

Temporal data on gene expression and context-specific open chromatin states can improve identification of key transcription factors (TFs) and the gene regulatory networks (GRNs) controlling cellular differentiation. However, their integration remains challenging. Here, we delineate a general approach for data-driven and unbiased identification of key TFs and dynamic GRNs, called EPIC-DREM. We generated time-series transcriptomic and epigenomic profiles during differentiation of mouse multipotent bone marrow stromal cell line (ST2) toward adipocytes and osteoblasts. Using our novel approach we constructed time-resolved GRNs for both lineages and identifed the shared TFs involved in both differentiation processes. To take an alternative approach to prioritize the identified shared regulators, we mapped dynamic super-enhancers in both lineages and associated them to target genes with correlated expression profiles. The combination of the two approaches identified aryl hydrocarbon receptor (AHR) and Glis family zinc finger 1 (GLIS1) as mesenchymal key TFs controlled by dynamic cell type-specific super-enhancers that become repressed in both lineages. AHR and GLIS1 control differentiation-induced genes and their overexpression can inhibit the lineage commitment of the multipotent bone marrow-derived ST2 cells.

HDAC1 links early life stress to schizophrenia-like phenotypes.

Schizophrenia is a devastating disease that arises on the background of genetic predisposition and environmental risk factors, such as early life stress (ELS). In this study, we show that ELS-induced schizophrenia-like phenotypes in mice correlate with a widespread increase of histone-deacetylase 1 (Hdac1) expression that is linked to altered DNA methylation. Hdac1 overexpression in neurons of the medial prefrontal cortex, but not in the dorsal or ventral hippocampus, mimics schizophrenia-like phenotypes induced by ELS. Systemic administration of an HDAC inhibitor rescues the detrimental effects of ELS when applied after the manifestation of disease phenotypes. In addition to the hippocampus and prefrontal cortex, mice subjected to ELS exhibit increased Hdac1 expression in blood. Moreover, Hdac1 levels are increased in blood samples from patients with schizophrenia who had encountered ELS, compared with patients without ELS experience. Our data suggest that HDAC1 inhibition should be considered as a therapeutic approach to treat schizophrenia.

TREM2 triggers microglial density and age-related neuronal loss.

The microglial triggering receptor expressed on myeloid cells 2 (TREM2) signals via the activatory membrane adaptor molecule TYROBP. Genetic variants or mutations of TREM2 or TYROBP have been linked to inflammatory neurodegenerative diseases associated with aging. The typical aging process goes along with microglial changes and mild neuronal loss, but the exact contribution of TREM2 is still unclear. Aged TREM2 knock-out mice showed decreased age-related neuronal loss in the substantia nigra and the hippocampus. Transcriptomic analysis of the brains of 24 months old TREM2 knock-out mice revealed 211 differentially expressed genes mostly downregulated and associated with complement activation and oxidative stress response pathways. Consistently, 24 months old TREM2 knock-out mice showed lower transcription of microglial (Aif1 and Tmem119), oxidative stress markers (Inos, Cyba, and Cybb) and complement components (C1qa, C1qb, C1qc, C3, C4b, Itgam, and Itgb2), decreased microglial numbers and expression of the microglial activation marker Cd68, as well as accumulation of oxidized lipids. Cultured microglia of TREM2 knock-out mice showed reduced phagocytosis and oxidative burst. Thus, microglial TREM2 contributes to age-related microglial changes, phagocytic oxidative burst, and loss of neurons with possible detrimental effects during physiological aging.

Alpha-synuclein deregulates the expression of COL4A2 and impairs ER-Golgi function.

Alpha-synuclein (aSyn) is the major protein component of Lewy bodies and Lewy neurites, the typical pathological hallmarks in Parkinson's disease (PD) and Dementia with Lewy bodies. aSyn is capable of inducing transcriptional deregulation, but the precise effect of specific aSyn mutants associated with familial forms of PD, remains unclear. Here, we used transgenic mice overexpressing human wild-type (WT) or A30P aSyn to compare the transcriptional profiles of the two animal models. We found that A30P aSyn promotes strong transcriptional deregulation and increases DNA binding. Interestingly, COL4A2, a major component of basement membranes, was found to be upregulated in both A30P aSyn transgenic mice and in dopaminergic neurons expressing A30P aSyn, suggesting a crucial role for collagen related genes in aSyn-induced toxicity. Finally, we observed that A30P aSyn alters Golgi morphology and increases the susceptibility to endoplasmic reticulum (ER) stress in dopaminergic cells. In total, our findings provide novel insight into the putative role of aSyn on transcription and on the molecular mechanisms involved, thereby opening novel avenues for future therapeutic interventions in PD and other synucleinopathies.

Promoter-proximal transcription factor binding is transcriptionally active when coupled with nucleosome repositioning in immediate vicinity.

Previous studies have analyzed patterns of transcription, transcription factor (TF) binding or mapped nucleosome occupancy across the genome. These suggest that the three aspects are genetically connected but the cause and effect relationships are still unknown. For example, physiologic TF binding studies involve many TFs, consequently, it is difficult to assign nucleosome reorganization to the binding site occupancy of any particular TF. Therefore, several aspects remain unclear: does TF binding influence nucleosome (re)organizations locally or impact the chromatin landscape at a more global level; are all or only a fraction of TF binding a result of reorganization in nucleosome occupancy and do all TF binding and associated changes in nucleosome occupancy result in altered gene expression? With these in mind, following characterization of two states (before and after induction of a single TF of choice) we determined: (i) genomic binding sites of the TF, (ii) promoter nucleosome occupancy and (iii) transcriptome profiles. Results demonstrated that promoter-proximal TF binding influenced expression of the target gene when it was coupled to nucleosome repositioning at or close to its binding site in most cases. In contrast, only in few cases change in target gene expression was found when TF binding occurred without local nucleosome reorganization.

Non-metastatic 2 (NME2)-mediated suppression of lung cancer metastasis involves transcriptional regulation of key cell adhesion factor vinculin.

Tumor metastasis refers to spread of a tumor from site of its origin to distant organs and causes majority of cancer deaths. Although >30 metastasis suppressor genes (MSGs) that negatively regulate metastasis have been identified so far, two issues are poorly understood: first, which MSGs oppose metastasis in a tumor type, and second, which molecular function of MSG controls metastasis. Herein, integrative analyses of tumor-transcriptomes (n=382), survival data (n=530) and lymph node metastases (n=100) in lung cancer patients identified non-metastatic 2 (NME2) as a key MSG from a pool of >30 metastasis suppressors. Subsequently, we generated a promoter-wide binding map for NME2 using chromatin immunoprecipitation with promoter microarrays (ChIP-chip), and transcriptome profiling. We discovered novel targets of NME2 which are involved in focal adhesion signaling. Importantly, we detected binding of NME2 in promoter of focal adhesion factor, vinculin. Reduced expression of NME2 led to enhanced transcription of vinculin. In comparison, NME1, a close homolog of NME2, did not bind to vinculin promoter nor regulate its expression. In line, enhanced metastasis of NME2-depleted lung cancer cells was found in zebrafish and nude mice tumor models. The metastatic potential of NME2-depleted cells was remarkably diminished upon selective RNA-i-mediated silencing of vinculin. Together, we demonstrate that reduced NME2 levels lead to transcriptional de-repression of vinculin and regulate lung cancer metastasis.

Histone-methyltransferase MLL2 (KMT2B) is required for memory formation in mice.

The consolidation of long-term memories requires differential gene expression. Recent research has suggested that dynamic changes in chromatin structure play a role in regulating the gene expression program linked to memory formation. The contribution of histone methylation, an important regulatory mechanism of chromatin plasticity that is mediated by the counteracting activity of histone-methyltransferases and histone-demethylases, is, however, not well understood. Here we show that mice lacking the histone-methyltransferase myeloid/lymphoid or mixed-lineage leukemia 2 (mll2/kmt2b) gene in adult forebrain excitatory neurons display impaired hippocampus-dependent memory function. Consistent with the role of KMT2B in gene-activation DNA microarray analysis revealed that 152 genes were downregulated in the hippocampal dentate gyrus region of mice lacking kmt2b. Downregulated plasticity genes showed a specific deficit in histone 3 lysine 4 di- and trimethylation, while histone 3 lysine 4 monomethylation was not affected. Our data demonstrates that KMT2B mediates hippocampal histone 3 lysine 4 di- and trimethylation and is a critical player for memory formation.

Quadruplex-single nucleotide polymorphisms (Quad-SNP) influence gene expression difference among individuals.

Non-canonical guanine quadruplex structures are not only predominant but also conserved among bacterial and mammalian promoters. Moreover recent findings directly implicate quadruplex structures in transcription. These argue for an intrinsic role of the structural motif and thereby posit that single nucleotide polymorphisms (SNP) that compromise the quadruplex architecture could influence function. To test this, we analysed SNPs within quadruplex motifs (Quad-SNP) and gene expression in 270 individuals across four populations (HapMap) representing more than 14,500 genotypes. Findings reveal significant association between quadruplex-SNPs and expression of the corresponding gene in individuals (P < 0.0001). Furthermore, analysis of Quad-SNPs obtained from population-scale sequencing of 1000 human genomes showed relative selection bias against alteration of the structural motif. To directly test the quadruplex-SNP-transcription connection, we constructed a reporter system using the RPS3 promoter-remarkable difference in promoter activity in the 'quadruplex-destabilized' versus 'quadruplex-intact' promoter was noticed. As a further test, we incorporated a quadruplex motif or its disrupted counterpart within a synthetic promoter reporter construct. The quadruplex motif, and not the disrupted-motif, enhanced transcription in human cell lines of different origin. Together, these findings build direct support for quadruplex-mediated transcription and suggest quadruplex-SNPs may play significant role in mechanistically understanding variations in gene expression among individuals.

2-Hydroxyglutarate modulates histone methylation at specific loci and alters gene expression via Rph1 inhibition.

2-Hydroxyglutarate (2-HG) is an oncometabolite that accumulates in certain cancers. Gain-of-function mutations in isocitrate dehydrogenase lead to 2-HG accumulation at the expense of alpha-ketoglutarate. Elevated 2-HG levels inhibit histone and DNA demethylases, causing chromatin structure and gene regulation changes with tumorigenic consequences. We investigated the effects of elevated 2-HG levels in Saccharomyces cerevisiae, a yeast devoid of DNA methylation and heterochromatin-associated histone methylation. Our results demonstrate genetic background-dependent gene expression changes and altered H3K4 and H3K36 methylation at specific loci. Analysis of histone demethylase deletion strains indicated that 2-HG inhibits Rph1 sufficiently to induce extensive gene expression changes. Rph1 is the yeast homolog of human KDM4 demethylases and, among the yeast histone demethylases, was the most sensitive to the inhibitory effect of 2-HG in vitro. Interestingly, Rph1 deficiency favors gene repression and leads to further down-regulation of already silenced genes marked by low H3K4 and H3K36 trimethylation, but abundant in H3K36 dimethylation. Our results provide novel insights into the genome-wide effects of 2-HG and highlight Rph1 as its preferential demethylase target.

Forecasting the dynamics of a complex microbial community using integrated meta-omics.

Predicting the behaviour of complex microbial communities is challenging. However, this is essential for complex biotechnological processes such as those in biological wastewater treatment plants (BWWTPs), which require sustainable operation. Here we summarize 14 months of longitudinal meta-omics data from a BWWTP anaerobic tank into 17 temporal signals, explaining 91.1% of the temporal variance, and link those signals to ecological events within the community. We forecast the signals over the subsequent five years and use 21 extra samples collected at defined time intervals for testing and validation. Our forecasts are correct for six signals and hint on phenomena such as predation cycles. Using all the 17 forecasts and the environmental variables, we predict gene abundance and expression, with a coefficient of determination ≥0.87 for the subsequent three years. Our study demonstrates the ability to forecast the dynamics of open microbial ecosystems using interactions between community cycles and environmental parameters.

The use of a complex tetra-culture alveolar model to study the biological effects induced by gold nanoparticles with different physicochemical properties.

A substantial increase in engineered nanoparticles in consumer products has been observed, heightening human and environmental exposure. Inhalation represents the primary route of human exposure, necessitating a focus on lung toxicity studies. However, to avoid ethical concerns the use of in vitro models is an efficient alternative to in vivo models. This study utilized an in vitro human alveolar barrier model at air-liquid-interface with four cell lines, for evaluating the biological effects of different gold nanoparticles. Exposure to PEGylated gold nanospheres, nanorods, and nanostars did not significantly impact viability after 24 h, yet all AuNPs induced cytotoxicity in the form of membrane integrity impairment. Gold quantification revealed cellular uptake and transport. Transcriptomic analysis identified gene expression changes, particularly related to the enhancement of immune cells. Despite limited impact, distinct effects were observed, emphasizing the influence of nanoparticles physicochemical parameters while demonstrating the model's efficacy in investigating particle biological effects.

Single cell transcriptome analysis of the THY-Tau22 mouse model of Alzheimer's disease reveals sex-dependent dysregulations.

Alzheimer's disease (AD) progression and pathology show pronounced sex differences, but the factors driving these remain poorly understood. To gain insights into early AD-associated molecular changes and their sex dependency for tau pathology in the cortex, we performed single-cell RNA-seq in the THY-Tau22 AD mouse model. By examining cell type-specific and cell type-agnostic AD-related gene activity changes and their sex-dimorphism for individual genes, pathways and cellular sub-networks, we identified both statistically significant alterations and interpreted the upstream mechanisms controlling them. Our results confirm several significant sex-dependent alterations in gene activity in the THY-Tau22 model mice compared to controls, with more pronounced alterations in females. Both changes shared across multiple cell types and cell type-specific changes were observed. The differential genes showed significant over-representation of known AD-relevant processes, such as pathways associated with neuronal differentiation, programmed cell death and inflammatory responses. Regulatory network analysis of these genes revealed upstream regulators that modulate many of the downstream targets with sex-dependent changes. Most key regulators have been previously implicated in AD, such as Egr1, Klf4, Chchd2, complement system genes, and myelin-associated glycoproteins. Comparing with similar data from the Tg2576 AD mouse model and human AD patients, we identified multiple genes with consistent, cell type-specific and sex-dependent alterations across all three datasets. These shared changes were particularly evident in the expression of myelin-associated genes such as Mbp and Plp1 in oligodendrocytes. In summary, we observed significant cell type-specific transcriptomic changes in the THY-Tau22 mouse model, with a strong over-representation of known AD-associated genes and processes. These include both sex-neutral and sex-specific patterns, characterized by consistent shifts in upstream master regulators and downstream target genes. Collectively, these findings provide insights into mechanisms influencing sex-specific susceptibility to AD and reveal key regulatory proteins that could be targeted for developing treatments addressing sex-dependent AD pathology.

Systems level analysis of sex-dependent gene expression changes in Parkinson's disease.

Parkinson's disease (PD) is a heterogeneous disorder, and among the factors which influence the symptom profile, biological sex has been reported to play a significant role. While males have a higher age-adjusted disease incidence and are more frequently affected by muscle rigidity, females present more often with disabling tremors. The molecular mechanisms involved in these differences are still largely unknown, and an improved understanding of the relevant factors may open new avenues for pharmacological disease modification. To help address this challenge, we conducted a meta-analysis of disease-associated molecular sex differences in brain transcriptomics data from case/control studies. Both sex-specific (alteration in only one sex) and sex-dimorphic changes (changes in both sexes, but with opposite direction) were identified. Using further systems level pathway and network analyses, coordinated sex-related alterations were studied. These analyses revealed significant disease-associated sex differences in mitochondrial pathways and highlight specific regulatory factors whose activity changes can explain downstream network alterations, propagated through gene regulatory cascades. Single-cell expression data analyses confirmed the main pathway-level changes observed in bulk transcriptomics data. Overall, our analyses revealed significant sex disparities in PD-associated transcriptomic changes, resulting in coordinated modulations of molecular processes. Among the regulatory factors involved, NR4A2 has already been reported to harbor rare mutations in familial PD and its pharmacological activation confers neuroprotective effects in toxin-induced models of Parkinsonism. Our observations suggest that NR4A2 may warrant further research as a potential adjuvant therapeutic target to address a subset of pathological molecular features of PD that display sex-associated profiles.

Transcriptional and Chromatin Accessibility Profiling of Neural Stem Cells Differentiating into Astrocytes Reveal Dynamic Signatures Affected under Inflammatory Conditions.

Astrocytes arise from multipotent neural stem cells (NSCs) and represent the most abundant cell type of the central nervous system (CNS), playing key roles in the developing and adult brain. Since the differentiation of NSCs towards a gliogenic fate is a precisely timed and regulated process, its perturbation gives rise to dysfunctional astrocytic phenotypes. Inflammation, which often underlies neurological disorders, including neurodevelopmental disorders and brain tumors, disrupts the accurate developmental process of NSCs. However, the specific consequences of an inflammatory environment on the epigenetic and transcriptional programs underlying NSCs' differentiation into astrocytes is unexplored. Here, we address this gap by profiling in mice glial precursors from neural tissue derived from early embryonic stages along their astrocytic differentiation trajectory in the presence or absence of tumor necrosis factor (TNF), a master pro-inflammatory cytokine. By using a combination of RNA- and ATAC-sequencing approaches, together with footprint and integrated gene regulatory network analyses, we here identify key differences during the differentiation of NSCs into astrocytes under physiological and inflammatory settings. In agreement with its role to turn cells resistant to inflammatory challenges, we detect Nrf2 as a master transcription factor supporting the astrocytic differentiation under TNF exposure. Further, under these conditions, we unravel additional transcriptional regulatory hubs, including Stat3, Smad3, Cebpb, and Nfkb2, highlighting the interplay among pathways underlying physiological astrocytic developmental processes and those involved in inflammatory responses, resulting in discrete astrocytic phenotypes. Overall, our study reports key transcriptional and epigenetic changes leading to the identification of molecular regulators of astrocytic differentiation. Furthermore, our analyses provide a valuable resource for understanding inflammation-induced astrocytic phenotypes that might contribute to the development and progression of CNS disorders with an inflammatory component.

Metabolic modelling-based in silico drug target prediction identifies six novel repurposable drugs for melanoma.

Despite high initial response rates to targeted kinase inhibitors, the majority of patients suffering from metastatic melanoma present with high relapse rates, demanding for alternative therapeutic options. We have previously developed a drug repurposing workflow to identify metabolic drug targets that, if depleted, inhibit the growth of cancer cells without harming healthy tissues. In the current study, we have applied a refined version of the workflow to specifically predict both, common essential genes across various cancer types, and melanoma-specific essential genes that could potentially be used as drug targets for melanoma treatment. The in silico single gene deletion step was adapted to simulate the knock-out of all targets of a drug on an objective function such as growth or energy balance. Based on publicly available, and in-house, large-scale transcriptomic data metabolic models for melanoma were reconstructed enabling the prediction of 28 candidate drugs and estimating their respective efficacy. Twelve highly efficacious drugs with low half-maximal inhibitory concentration values for the treatment of other cancers, which are not yet approved for melanoma treatment, were used for in vitro validation using melanoma cell lines. Combination of the top 4 out of 6 promising candidate drugs with BRAF or MEK inhibitors, partially showed synergistic growth inhibition compared to individual BRAF/MEK inhibition. Hence, the repurposing of drugs may enable an increase in therapeutic options e.g., for non-responders or upon acquired resistance to conventional melanoma treatments.

Impaired neuron differentiation in GBA-associated Parkinson's disease is linked to cell cycle defects in organoids.

The mechanisms underlying Parkinson's disease (PD) etiology are only partially understood despite intensive research conducted in the field. Recent evidence suggests that early neurodevelopmental defects might play a role in cellular susceptibility to neurodegeneration. To study the early developmental contribution of GBA mutations in PD we used patient-derived iPSCs carrying a heterozygous N370S mutation in the GBA gene. Patient-specific midbrain organoids displayed GBA-PD relevant phenotypes such as reduction of GCase activity, autophagy impairment, and mitochondrial dysfunction. Genome-scale metabolic (GEM) modeling predicted changes in lipid metabolism which were validated with lipidomics analysis, showing significant differences in the lipidome of GBA-PD. In addition, patient-specific midbrain organoids exhibited a decrease in the number and complexity of dopaminergic neurons. This was accompanied by an increase in the neural progenitor population showing signs of oxidative stress-induced damage and premature cellular senescence. These results provide insights into how GBA mutations may lead to neurodevelopmental defects thereby predisposing to PD pathology.

Fiber deprivation and microbiome-borne curli shift gut bacterial populations and accelerate disease in a mouse model of Parkinson's disease.

Parkinson's disease (PD) is a neurological disorder characterized by motor dysfunction, dopaminergic neuron loss, and alpha-synuclein (αSyn) inclusions. Many PD risk factors are known, but those affecting disease progression are not. Lifestyle and microbial dysbiosis are candidates in this context. Diet-driven gut dysbiosis and reduced barrier function may increase exposure of enteric neurons to toxins. Here, we study whether fiber deprivation and exposure to bacterial curli, a protein cross-seeding with αSyn, individually or together, exacerbate disease in the enteric and central nervous systems of a transgenic PD mouse model. We analyze the gut microbiome, motor behavior, and gastrointestinal and brain pathologies. We find that diet and bacterial curli alter the microbiome and exacerbate motor performance, as well as intestinal and brain pathologies, but to different extents. Our results shed important insights on how diet and microbiome-borne insults modulate PD progression via the gut-brain axis and have implications for lifestyle management of PD.