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. Graphical

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.

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.

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.

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.

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.