Multi-omics integration of scRNA-seq time series data predicts new intervention points for Parkinson's disease.

Parkinson's disease (PD) is a complex neurodegenerative disorder without a cure. The onset of PD symptoms corresponds to 50% loss of midbrain dopaminergic (mDA) neurons, limiting early-stage understanding of PD. To shed light on early PD development, we study time series scRNA-seq datasets of mDA neurons obtained from patient-derived induced pluripotent stem cell differentiation. We develop a new data integration method based on Non-negative Matrix Tri-Factorization that integrates these datasets with molecular interaction networks, producing condition-specific "gene embeddings". By mining these embeddings, we predict 193 PD-related genes that are largely supported (49.7%) in the literature and are specific to the investigated PINK1 mutation. Enrichment analysis in Kyoto Encyclopedia of Genes and Genomes pathways highlights 10 PD-related molecular mechanisms perturbed during early PD development. Finally, investigating the top 20 prioritized genes reveals 12 previously unrecognized genes associated with PD that represent interesting drug targets.

Glioblastoma-instructed microglia transition to heterogeneous phenotypic states with phagocytic and dendritic cell-like features in patient tumors and patient-derived orthotopic xenografts.

BACKGROUND: A major contributing factor to glioblastoma (GBM) development and progression is its ability to evade the immune system by creating an immune-suppressive environment, where GBM-associated myeloid cells, including resident microglia and peripheral monocyte-derived macrophages, play critical pro-tumoral roles. However, it is unclear whether recruited myeloid cells are phenotypically and functionally identical in GBM patients and whether this heterogeneity is recapitulated in patient-derived orthotopic xenografts (PDOXs). A thorough understanding of the GBM ecosystem and its recapitulation in preclinical models is currently missing, leading to inaccurate results and failures of clinical trials. METHODS: Here, we report systematic characterization of the tumor microenvironment (TME) in GBM PDOXs and patient tumors at the single-cell and spatial levels. We applied single-cell RNA sequencing, spatial transcriptomics, multicolor flow cytometry, immunohistochemistry, and functional studies to examine the heterogeneous TME instructed by GBM cells. GBM PDOXs representing different tumor phenotypes were compared to glioma mouse GL261 syngeneic model and patient tumors. RESULTS: We show that GBM tumor cells reciprocally interact with host cells to create a GBM patient-specific TME in PDOXs. We detected the most prominent transcriptomic adaptations in myeloid cells, with brain-resident microglia representing the main population in the cellular tumor, while peripheral-derived myeloid cells infiltrated the brain at sites of blood-brain barrier disruption. More specifically, we show that GBM-educated microglia undergo transition to diverse phenotypic states across distinct GBM landscapes and tumor niches. GBM-educated microglia subsets display phagocytic and dendritic cell-like gene expression programs. Additionally, we found novel microglial states expressing cell cycle programs, astrocytic or endothelial markers. Lastly, we show that temozolomide treatment leads to transcriptomic plasticity and altered crosstalk between GBM tumor cells and adjacent TME components. CONCLUSIONS: Our data provide novel insights into the phenotypic adaptation of the heterogeneous TME instructed by GBM tumors. We show the key role of microglial phenotypic states in supporting GBM tumor growth and response to treatment. Our data place PDOXs as relevant models to assess the functionality of the TME and changes in the GBM ecosystem upon treatment.

Early-to-mid stage idiopathic Parkinson's disease shows enhanced cytotoxicity and differentiation in CD8 T-cells in females.

Neuroinflammation in the brain contributes to the pathogenesis of Parkinson's disease (PD), but the potential dysregulation of peripheral immunity has not been systematically investigated for idiopathic PD (iPD). Here we showed an elevated peripheral cytotoxic immune milieu, with more terminally-differentiated effector memory (TEMRA) CD8 T, CD8+ NKT cells and circulating cytotoxic molecules in fresh blood of patients with early-to-mid iPD, especially females, after analyzing > 700 innate and adaptive immune features. This profile, also reflected by fewer CD8+FOXP3+ T cells, was confirmed in another subcohort. Co-expression between cytotoxic molecules was selectively enhanced in CD8 TEMRA and effector memory (TEM) cells. Single-cell RNA-sequencing analysis demonstrated the accelerated differentiation within CD8 compartments, enhanced cytotoxic pathways in CD8 TEMRA and TEM cells, while CD8 central memory (TCM) and naïve cells were already more-active and transcriptionally-reprogrammed. Our work provides a comprehensive map of dysregulated peripheral immunity in iPD, proposing candidates for early diagnosis and treatments.

Glioblastoma-instructed microglia transition to heterogeneous phenotypic states with phagocytic and dendritic cell-like features in patient tumors and patient-derived orthotopic xenografts.

Background: A major contributing factor to glioblastoma (GBM) development and progression is its ability to evade the immune system by creating an immune-suppressive environment, where GBM-associated myeloid cells, including resident microglia and peripheral monocyte-derived macrophages, play critical pro-tumoral roles. However, it is unclear whether recruited myeloid cells are phenotypically and functionally identical in GBM patients and whether this heterogeneity is recapitulated in patient-derived orthotopic xenografts (PDOXs). A thorough understanding of the GBM ecosystem and its recapitulation in preclinical models is currently missing, leading to inaccurate results and failures of clinical trials. Methods: Here, we report systematic characterization of the tumor microenvironment (TME) in GBM PDOXs and patient tumors at the single-cell and spatial levels. We applied single-cell RNA-sequencing, spatial transcriptomics, multicolor flow cytometry, immunohistochemistry and functional studies to examine the heterogeneous TME instructed by GBM cells. GBM PDOXs representing different tumor phenotypes were compared to glioma mouse GL261 syngeneic model and patient tumors. Results: We show that GBM tumor cells reciprocally interact with host cells to create a GBM patient-specific TME in PDOXs. We detected the most prominent transcriptomic adaptations in myeloid cells, with brain-resident microglia representing the main population in the cellular tumor, while peripheral-derived myeloid cells infiltrated the brain at sites of blood-brain barrier disruption. More specifically, we show that GBM-educated microglia undergo transition to diverse phenotypic states across distinct GBM landscapes and tumor niches. GBM-educated microglia subsets display phagocytic and dendritic cell-like gene expression programs. Additionally, we found novel microglial states expressing cell cycle programs, astrocytic or endothelial markers. Lastly, we show that temozolomide treatment leads to transcriptomic plasticity and altered crosstalk between GBM tumor cells and adjacent TME components. Conclusions: Our data provide novel insights into the phenotypic adaptation of the heterogeneous TME instructed by GBM tumors. We show the key role of microglial phenotypic states in supporting GBM tumor growth and response to treatment. Our data place PDOXs as relevant models to assess the functionality of the TME and changes in the GBM ecosystem upon treatment.

Rac1 and Rac3 GTPases and TPC2 are required for axonal outgrowth and migration of cortical interneurons.

Rho GTPases, among them Rac1 and Rac3, are major transducers of extracellular signals and are involved in multiple cellular processes. In cortical interneurons, the neurons that control the balance between excitation and inhibition of cortical circuits, Rac1 and Rac3 are essential for their development. Ablation of both leads to a severe reduction in the numbers of mature interneurons found in the murine cortex, which is partially due to abnormal cell cycle progression of interneuron precursors and defective formation of growth cones in young neurons. Here, we present new evidence that upon Rac1 and Rac3 ablation, centrosome, Golgi complex and lysosome positioning is significantly perturbed, thus affecting both interneuron migration and axon growth. Moreover, for the first time, we provide evidence of altered expression and localization of the two-pore channel 2 (TPC2) voltage-gated ion channel that mediates Ca2+ release. Pharmacological inhibition of TPC2 negatively affected axonal growth and migration of interneurons. Our data, taken together, suggest that TPC2 contributes to the severe phenotype in axon growth initiation, extension and interneuron migration in the absence of Rac1 and Rac3.

Assessing Genetic Variation in Wild and Domesticated Pikeperch Populations: Implications for Conservation and Fish Farming.

The pikeperch is a freshwater/brackish water fish species with growing interest for European aquaculture. Wild populations show signs of decline in many areas of the species natural range due to human activities. The comparative evaluation of genetic status in wild and domesticated populations is extremely useful for the future establishment of genetic breeding programs. The main objective of the present study was to assess and compare the genetic variability of 13 domesticated populations from commercial farms and 8 wild populations, developing an efficient microsatellite multiplex tool for genotyping. Partial cytochrome b gene sequences were also used to infer phylogeographic relationships. Results show that on average, the domesticated populations do not exhibit significantly lower levels of genetic diversity compared to the wild ones and do not suffer from inbreeding. Nuclear data provide evidence that pikeperch populations in Europe belong to at least two genetically differentiated groups: the first one is predominantly present in Northern Europe and around the Baltic Sea, while the second one comprises populations from Central Europe. In this second group, Hungarian origin populations constitute a differentiated stock that needs special consideration. Aquaculture broodstocks analyzed appear to contain fish of a single origin with only a few exceptions.

Single-cell transcriptomics of human iPSC differentiation dynamics reveal a core molecular network of Parkinson's disease.

Parkinson's disease (PD) is the second-most prevalent neurodegenerative disorder, characterized by the loss of dopaminergic neurons (mDA) in the midbrain. The underlying mechanisms are only partly understood and there is no treatment to reverse PD progression. Here, we investigated the disease mechanism using mDA neurons differentiated from human induced pluripotent stem cells (hiPSCs) carrying the ILE368ASN mutation within the PINK1 gene, which is strongly associated with PD. Single-cell RNA sequencing (RNAseq) and gene expression analysis of a PINK1-ILE368ASN and a control cell line identified genes differentially expressed during mDA neuron differentiation. Network analysis revealed that these genes form a core network, members of which interact with all known 19 protein-coding Parkinson's disease-associated genes. This core network encompasses key PD-associated pathways, including ubiquitination, mitochondrial function, protein processing, RNA metabolism, and vesicular transport. Proteomics analysis showed a consistent alteration in proteins of dopamine metabolism, indicating a defect of dopaminergic metabolism in PINK1-ILE368ASN neurons. Our findings suggest the existence of a network onto which pathways associated with PD pathology converge, and offers an inclusive interpretation of the phenotypic heterogeneity of PD.

Single-Cell Transcriptomics and In Situ Morphological Analyses Reveal Microglia Heterogeneity Across the Nigrostriatal Pathway.

Microglia are the resident immune effector cells of the central nervous system (CNS) rapidly reacting to various pathological stimuli to maintain CNS homeostasis. However, microglial reactions in the CNS may also worsen neurological disorders. Hence, the phenotypic analysis of microglia in healthy tissue may identify specific poised subsets ultimately supporting or harming the neuronal network. This is all the more important for the understanding of CNS disorders exhibiting regional-specific and cellular pathological hallmarks, such as many neurodegenerative disorders, including Parkinson's disease (PD). In this context, we aimed to address the heterogeneity of microglial cells in susceptible brain regions for PD, such as the nigrostriatal pathway. Here, we combined single-cell RNA-sequencing with immunofluorescence analyses of the murine nigrostriatal pathway, the most affected brain region in PD. We uncovered a microglia subset, mainly present in the midbrain, displaying an intrinsic transcriptional immune alerted signature sharing features of inflammation-induced microglia. Further, an in situ morphological screening of inferred cellular diversity showed a decreased microglia complexity in the midbrain when compared to striatum. Our study provides a resource for the identification of specific microglia phenotypes within the nigrostriatal pathway, which may be relevant in PD.

Scanning of Genetic Variants and Genetic Mapping of Phenotypic Traits in Gilthead Sea Bream Through ddRAD Sequencing.

Gilthead sea bream (Sparus aurata) is a teleost of considerable economic importance in Southern European aquaculture. The aquaculture industry shows a growing interest in the application of genetic methods that can locate phenotype-genotype associations with high economic impact. Through selective breeding, the aquaculture industry can exploit this information to maximize the financial yield. Here, we present a Genome Wide Association Study (GWAS) of 112 samples belonging to seven different sea bream families collected from a Greek commercial aquaculture company. Through double digest Random Amplified DNA (ddRAD) Sequencing, we generated a per-sample genetic profile consisting of 2,258 high-quality Single Nucleotide Polymorphisms (SNPs). These profiles were tested for association with four phenotypes of major financial importance: Fat, Weight, Tag Weight, and the Length to Width ratio. We applied two methods of association analysis. The first is the typical single-SNP to phenotype test, and the second is a feature selection (FS) method through two novel algorithms that are employed for the first time in aquaculture genomics and produce groups with multiple SNPs associated to a phenotype. In total, we identified 9 single SNPs and 6 groups of SNPs associated with weight-related phenotypes (Weight and Tag Weight), 2 groups associated with Fat, and 16 groups associated with the Length to Width ratio. Six identified loci (Chr4:23265532, Chr6:12617755, Chr:8:11613979, Chr13:1098152, Chr15:3260819, and Chr22:14483563) were present in genes associated with growth in other teleosts or even mammals, such as semaphorin-3A and neurotrophin-3. These loci are strong candidates for future studies that will help us unveil the genetic mechanisms underlying growth and improve the sea bream aquaculture productivity by providing genomic anchors for selection programs.

In Vivo Effects of Lipopolysaccharide on Peroxisome Proliferator-Activated Receptor Expression in Juvenile Gilthead Seabream (Sparus Aurata).

Fish are constantly exposed to microorganisms in the aquatic environment, many of which are bacterial pathogens. Bacterial pathogens activate the innate immune response in fish involving the production of pro-inflammatory molecules that, in addition to their immune-related role, can affect non-immune tissues. In the present study, we aimed at investigating how inflammatory responses can affect metabolic homeostasis in the gilthead seabream (Sparus aurata), a teleost of considerable economic importance in Southern European countries. Specifically, we mimicked a bacterial infection by in vivo administration of lipopolysaccharide (LPS, 6 mg/kg body weight) and measured metabolic parameters in the blood and, importantly, the mRNA expression levels of the three isotypes of peroxisome proliferator activated receptors (PPARα, ß, and γ) in metabolically-relevant tissues in seabream. PPARs are nuclear receptors that are important for lipid and carbohydrate metabolism in mammals and that act as biological sensors of altered lipid metabolism. We show here that LPS-induced inflammatory responses result in the modulation of triglyceride plasma levels that are accompanied most notably by a decrease in the hepatic mRNA expression levels of PPARα, ß, and γ and by the up-regulation of PPARγ expression only in adipose tissue and the anterior intestine. In addition, LPS-induced inflammation results in an increase in the hepatic mRNA expression and protein activity levels of members of the mitogen-activated protein kinase (MAPK) family, known in mammals to regulate the transcription and activity of PPARs. Our results provide evidence for the involvement of PPARs in the metabolic response to inflammatory stimuli in seabream and offer insights into the molecular mechanisms underlying the redirection of metabolic activities under inflammatory conditions in vertebrates.