Predictive metabolic networks reveal sex- and APOE genotype-specific metabolic signatures and drivers for precision medicine in Alzheimer's disease.

INTRODUCTION: Late-onset Alzheimer's disease (LOAD) is a complex neurodegenerative disease characterized by multiple progressive stages, glucose metabolic dysregulation, Alzheimer's disease (AD) pathology, and inexorable cognitive decline. Discovery of metabolic profiles unique to sex, apolipoprotein E (APOE) genotype, and stage of disease progression could provide critical insights for personalized LOAD medicine. METHODS: Sex- and APOE-specific metabolic networks were constructed based on changes in 127 metabolites of 656 serum samples from the Alzheimer's Disease Neuroimaging Initiative cohort. RESULTS: Application of an advanced analytical platform identified metabolic drivers and signatures clustered with sex and/or APOE ɛ4, establishing patient-specific biomarkers predictive of disease state that significantly associated with cognitive function. Presence of the APOE ɛ4 shifts metabolic signatures to a phosphatidylcholine-focused profile overriding sex-specific differences in serum metabolites of AD patients. DISCUSSION: These findings provide an initial but critical step in developing a diagnostic platform for personalized medicine by integrating metabolomic profiling and cognitive assessments to identify targeted precision therapeutics for AD patient subgroups through computational network modeling.

Clonal Kinetics and Single-Cell Transcriptional Profiles of T Cells Mobilized to Blood by Acute Exercise.

PURPOSE: Acute exercise redistributes large numbers of memory T cells, which may contribute to enhanced immune surveillance in regular exercisers. It is not known, however, if acute exercise promotes a broad or oligoclonal T-cell receptor (TCR) repertoire or evokes transcriptomic changes in "exercise-responsive" T-cell clones. METHODS: Healthy volunteers completed a graded bout of cycling exercise up to 80% V̇O 2max . DNA was extracted from peripheral blood mononuclear cells collected at rest, during exercise (EX), and 1 h after (+1H) exercise, and processed for deep TCR-ß chain sequencing and tandem single-cell RNA sequencing. RESULTS: The number of unique clones and unique rearrangements was decreased at EX compared with rest ( P < 0.01) and +1H ( P < 0.01). Productive clonality was increased compared with rest ( P < 0.05) and +1H ( P < 0.05), whereas Shannon's Index was decreased compared with rest ( P < 0.05) and +1H ( P < 0.05). The top 10 rearrangements in the repertoire were increased at EX compared with rest ( P < 0.05) and +1H ( P < 0.05). Cross-referencing TCR-ß sequences with a public database (VDJdb) revealed that exercise increased the number of clones specific for the most prevalent motifs, including Epstein-Barr virus, cytomegalovirus, and influenza A. We identified 633 unique exercise-responsive T-cell clones that were mobilized and/or egressed in response to exercise. Among these clones, there was an upregulation in genes related to cell death, cytotoxicity, and activation ( P < 0.05). CONCLUSIONS: Acute exercise promotes an oligoclonal T-cell repertoire by preferentially mobilizing the most dominant clones, several of which are specific to known viral antigens and display differentially expressed genes indicative of cytotoxicity, activation, and apoptosis.

Predictive network analysis identifies JMJD6 and other potential key drivers in Alzheimer's disease.

Despite decades of genetic studies on late-onset Alzheimer's disease, the underlying molecular mechanisms remain unclear. To better comprehend its complex etiology, we use an integrative approach to build robust predictive (causal) network models using two large human multi-omics datasets. We delineate bulk-tissue gene expression into single cell-type gene expression and integrate clinical and pathologic traits, single nucleotide variation, and deconvoluted gene expression for the construction of cell type-specific predictive network models. Here, we focus on neuron-specific network models and prioritize 19 predicted key drivers modulating Alzheimer's pathology, which we then validate by knockdown in human induced pluripotent stem cell-derived neurons. We find that neuronal knockdown of 10 of the 19 targets significantly modulates levels of amyloid-beta and/or phosphorylated tau peptides, most notably JMJD6. We also confirm our network structure by RNA sequencing in the neurons following knockdown of each of the 10 targets, which additionally predicts that they are upstream regulators of REST and VGF. Our work thus identifies robust neuronal key drivers of the Alzheimer's-associated network state which may represent therapeutic targets with relevance to both amyloid and tau pathology in Alzheimer's disease.

Human lymphocytes mobilized with exercise have an anti-tumor transcriptomic profile and exert enhanced graft-versus-leukemia effects in xenogeneic mice.

Background: Every bout of exercise mobilizes and redistributes large numbers of effector lymphocytes with a cytotoxic and tissue migration phenotype. The frequent redistribution of these cells is purported to increase immune surveillance and play a mechanistic role in reducing cancer risk and slowing tumor progression in physically active cancer survivors. Our aim was to provide the first detailed single cell transcriptomic analysis of exercise-mobilized lymphocytes and test their effectiveness as a donor lymphocyte infusion (DLI) in xenogeneic mice engrafted with human leukemia. Methods: Peripheral blood mononuclear cells (PBMCs) were collected from healthy volunteers at rest and at the end of an acute bout of cycling exercise. Flow cytometry and single-cell RNA sequencing was performed to identify phenotypic and transcriptomic differences between resting and exercise-mobilized cells using a targeted gene expression panel curated for human immunology. PBMCs were injected into the tail vein of xenogeneic NSG-IL-15 mice and subsequently challenged with a luciferase tagged chronic myelogenous leukemia cell line (K562). Tumor growth (bioluminescence) and xenogeneic graft-versus-host disease (GvHD) were monitored bi-weekly for 40-days. Results: Exercise preferentially mobilized NK-cell, CD8+ T-cell and monocyte subtypes with a differentiated and effector phenotype, without significantly mobilizing CD4+ regulatory T-cells. Mobilized effector lymphocytes, particularly effector-memory CD8+ T-cells and NK-cells, displayed differentially expressed genes and enriched gene sets associated with anti-tumor activity, including cytotoxicity, migration/chemotaxis, antigen binding, cytokine responsiveness and alloreactivity (e.g. graft-versus-host/leukemia). Mice receiving exercise-mobilized PBMCs had lower tumor burden and higher overall survival (4.14E+08 photons/s and 47%, respectively) at day 40 compared to mice receiving resting PBMCs (12.1E+08 photons/s and 22%, respectively) from the same donors (p<0.05). Human immune cell engraftment was similar for resting and exercise-mobilized DLI. However, when compared to non-tumor bearing mice, K562 increased the expansion of NK-cell and CD3+/CD4-/CD8- T-cells in mice receiving exercise-mobilized but not resting lymphocytes, 1-2 weeks after DLI. No differences in GvHD or GvHD-free survival was observed between groups either with or without K562 challenge. Conclusion: Exercise in humans mobilizes effector lymphocytes with an anti-tumor transcriptomic profile and their use as DLI extends survival and enhances the graft-versus-leukemia (GvL) effect without exacerbating GvHD in human leukemia bearing xenogeneic mice. Exercise may serve as an effective and economical adjuvant to increase the GvL effects of allogeneic cell therapies without intensifying GvHD.

Determinants of de novo B cell responses to drifted epitopes in post-vaccination SARS-CoV-2 infections.

Vaccine-induced immunity may impact subsequent de novo responses to drifted epitopes in SARS-CoV-2 variants, but this has been difficult to quantify due to the challenges in recruiting unvaccinated control groups whose first exposure to SARS-CoV-2 is a primary infection. Through local, statewide, and national SARS-CoV-2 testing programs, we were able to recruit cohorts of individuals who had recovered from either primary or post-vaccination infections by either the Delta or Omicron BA.1 variants. Regardless of variant, we observed greater Spike-specific and neutralizing antibody responses in post-vaccination infections than in those who were infected without prior vaccination. Through analysis of variant-specific memory B cells as markers of de novo responses, we observed that Delta and Omicron BA.1 infections led to a marked shift in immunodominance in which some drifted epitopes elicited minimal responses, even in primary infections. Prior immunity through vaccination had a small negative impact on these de novo responses, but this did not correlate with cross-reactive memory B cells, arguing against competitive inhibition of naïve B cells. We conclude that dampened de novo B cell responses against drifted epitopes are mostly a function of altered immunodominance hierarchies that are apparent even in primary infections, with a more modest contribution from pre-existing immunity, perhaps due to accelerated antigen clearance.

Acute exercise mobilizes NKT-like cells with a cytotoxic transcriptomic profile but does not augment the potency of cytokine-induced killer (CIK) cells.

CD3+/CD56+ Natural killer (NK) cell-like T-cells (NKT-like cells) represent <5% of blood lymphocytes, display a cytotoxic phenotype, and can kill various cancers. NKT-like cells can be expanded ex vivo into cytokine-induced killer (CIK) cells, however this therapeutic cell product has had mixed results against hematological malignancies in clinical trials. The aim of this study was to determine if NKT-like cells mobilized during acute cycling exercise could be used to generate more potent anti-tumor CIK cells from healthy donors. An acute exercise bout increased NKT-like cell numbers in blood 2-fold. Single cell RNA sequencing revealed that exercise mobilized NKT-like cells have an upregulation of genes and transcriptomic programs associated with enhanced anti-tumor activity, including cytotoxicity, cytokine responsiveness, and migration. Exercise, however, did not augment the ex vivo expansion of CIK cells or alter their surface phenotypes after 21-days of culture. CIK cells expanded at rest, during exercise (at 60% and 80% VO2max) or after (1h post) were equally capable of killing leukemia, lymphoma, and multiple myeloma target cells with and without cytokine (IL-2) and antibody (OKT3) priming in vitro. We conclude that acute exercise in healthy donors mobilizes NKT-like cells with an upregulation of transcriptomic programs involved in anti-tumor activity, but does not augment the ex vivo expansion of CIK cells.

Transcriptomic analysis identifies differences in gene expression in actinic keratoses after treatment with imiquimod and between responders and non responders.

The presence of actinic keratoses (AKs) increases a patient's risk of developing squamous cell carcinoma by greater than six-fold. We evaluated the effect of topical treatment with imiquimod on the tumor microenvironment by measuring transcriptomic differences in AKs before and after treatment with imiquimod 3.75%. Biopsies were collected prospectively from 21 patients and examined histologically. RNA was extracted and transcriptomic analyses of 788 genes were performed using the nanoString assay. Imiquimod decreased number of AKs by study endpoint at week 14 (p < 0.0001). Post-imiquimod therapy, levels of CDK1, CXCL13, IL1B, GADPH, TTK, ILF3, EWSR1, BIRC5, PLAUR, ISG20, and C1QBP were significantly lower (adjusted p < 0.05). Complete responders (CR) exhibited a distinct pattern of inflammatory gene expression pre-treatment relative to incomplete responders (IR), with alterations in 15 inflammatory pathways (p < 0.05) reflecting differential expression of 103 genes (p < 0.05). Presence of adverse effects was associated with improved treatment response. Differences in gene expression were found between pre-treatment samples in CR versus IR, suggesting that higher levels of inflammation pre-treament may play a part in regression of AKs. Further characterization of the immune micro-environment in AKs may help develop biomarkers predictive of response to topical immune modulators and may guide therapy.

Rare Tumor-Normal Matched Whole Exome Sequencing Identifies Novel Genomic Pathogenic Germline and Somatic Aberrations.

Whole exome sequencing (WES) of matched tumor-normal pairs in rare tumors has the potential to identify genome-wide mutations and copy number alterations (CNAs). We evaluated 27 rare cancer patients with tumor-normal matching by WES and tumor-only next generation sequencing (NGS) as a comparator. Our goal was to: 1) identify known and novel variants and CNAs in rare cancers with comparison to common cancers; 2) examine differences between germline and somatic variants and how that functionally impacts rare tumors; 3) detect and characterize alleles in biologically relevant genes-pathways that may be of clinical importance but not represented in classical cancer genes. We identified 3343 germline single nucleotide variants (SNVs) and small indel variants-1670 in oncogenes and 1673 in tumor suppressor genes-generating an average of 124 germline variants/case. The number of somatic SNVs and small indels detected in all cases was 523:306 in oncogenes and 217 in tumor suppressor genes. Of the germline variants, six were identified to be pathogenic or likely pathogenic. In the 27 analyzed rare cancer cases, CNAs are variable depending on tumor type, germline pathogenic variants are more common. Cell fate pathway mutations (e.g., Hippo, Notch, Wnt) dominate pathogenesis and double hit (mutation + CNV) represent ~18% cases.

Leukocyte expression profiles reveal gene sets with prognostic value for seizure-free outcome following stereotactic laser amygdalohippocampotomy.

Among patients with intractable epilepsy, the most commonly performed surgical procedure is craniotomy for amygdalohippocampectomy (AH). Stereotactic laser amygdalohippocampotomy (SLAH) has also been recently employed as a minimally invasive treatment for intractable temporal lobe epilepsy (TLE). Among patients treated with AH and SLAH approximately 65% and 54% of patients become seizure-free, respectively. Therefore, selection criteria for surgical candidates with improved prognostic value for post-operative seizure-free outcome are greatly needed. In this study, we perform RNA sequencing (RNA-Seq) on whole blood leukocyte samples taken from 16 patients with intractable TLE prior to SLAH to test the hypothesis that pre-operative leukocyte RNA expression profiles are prognostic for post-operative seizure outcome. Multidimensional scaling analysis of the RNA expression data indicated separate clustering of patients with seizure free (SF) and non-seizure-free (NSF) outcomes. Differential expression (DE) analysis performed on SF versus NSF groups revealed 24 significantly differentially expressed genes (≥2.0-fold change, p-value < 0.05, FDR <0.05). Network and pathway analyses identified differential activation of pathways involved in lipid metabolism, morphology of oligodendrocytes, inflammatory response, and development of astrocytes. These results suggest that pre-operative leukocyte expression profiles have prognostic value for seizure outcome following SLAH.

Altered expression of signaling pathways regulating neuronal excitability in hippocampal tissue of temporal lobe epilepsy patients with low and high seizure frequency.

Despite recent advances in our understanding of synaptic transmission associated with epileptogenesis, the molecular mechanisms that control seizure frequency in patients with temporal lobe epilepsy (TLE) remain obscure. RNA-Seq was performed on hippocampal tissue resected from 12 medically intractable TLE patients with pre-surgery seizure frequencies ranging from 0.33 to 120 seizures per month. Differential expression (DE) analysis of individuals with low (LSF, mean = 4 seizure/month) versus high (HSF, mean = 60 seizures/month) seizure frequency identified 979 genes with ≥2-fold change in transcript abundance (FDR-adjusted p-value ö0.05). Comparisons with post-mortem controls revealed a large number of downregulated genes in the HSF (1676) versus LSF (399) groups. More than 50 signaling pathways were inferred to be deactivated or activated, with Signal Transduction as the central hub in the pathway network. While neuroinflammation pathways were activated in both groups, key neuronal system pathways were systematically deactivated in the HSF group, including calcium, CREB and Opioid signaling. We also infer that enhanced expression of a signaling cascade promoting synaptic downscaling may have played a key role in maintaining a higher seizure threshold in the LSF cohort. These results suggest that therapeutic approaches targeting synaptic scaling pathways may aid in the treatment of seizures in TLE.

Unexpected favorable outcome in a patient with high grade B-cell lymphoma with abnormalities of MYC, BCL6 and BCL2 loci.

High grade B-cell lymphoma (HGBCL) by WHO 2016 classification requires rearrangements of MYC and BCL2 and/or BCL6, practically covering the so called "double-hit" or "triple hit" lymphomas. We report a case of HGBCL "triple-hit" lymphoma in a 64-year old female. Cytogenetic and fluorescence in situ hybridization (FISH) studies revealed complex karyotype including rearrangement of MYC to a novel, non-IG partner on chromosome 18, and rearrangement of BCL2, BCL6 and IGH as well as ins(3)(q21q27.3q25.1) among other abnormalities. FISH studies showed five copies of MYC and 3-8 copies of BCL2. Gene expression analysis by RNA sequencing showed that MYC, BCL2 and MECOM genes were overexpressed whereas BCL6 was under-expressed. BCL6 was fused to MBNL1 gene due to complex structural rearrangement. MYC was expressed in >70% of cells and BCL2 was diffusely but highly expressed by immunohistochemistry. No pathogenic mutations were identified by sequencing a 26-gene panel including TP53. The patient has unexpectedly been in complete remission for 12 months after diagnosis after intensive chemotherapy including DA-EPOCH regimen despite having HGBCL. The prognostication of HGBCL patients may further be improved by the sub-categorization of these lymphomas on the basis of more detailed genomic markers than merely the WHO 2016 classification.

Rare variants of small effect size in neuronal excitability genes influence clinical outcome in Japanese cases of SCN1A truncation-positive Dravet syndrome.

Dravet syndrome (DS) is a rare, devastating form of childhood epilepsy that is often associated with mutations in the voltage-gated sodium channel gene, SCN1A. There is considerable variability in expressivity within families, as well as among individuals carrying the same primary mutation, suggesting that clinical outcome is modulated by variants at other genes. To identify modifier gene variants that contribute to clinical outcome, we sequenced the exomes of 22 individuals at both ends of a phenotype distribution (i.e., mild and severe cognitive condition). We controlled for variation associated with different mutation types by limiting inclusion to individuals with a de novo truncation mutation resulting in SCN1A haploinsufficiency. We performed tests aimed at identifying 1) single common variants that are enriched in either phenotypic group, 2) sets of common or rare variants aggregated in and around genes associated with clinical outcome, and 3) rare variants in 237 candidate genes associated with neuronal excitability. While our power to identify enrichment of a common variant in either phenotypic group is limited as a result of the rarity of mild phenotypes in individuals with SCN1A truncation variants, our top candidates did not map to functional regions of genes, or in genes that are known to be associated with neurological pathways. In contrast, we found a statistically-significant excess of rare variants predicted to be damaging and of small effect size in genes associated with neuronal excitability in severely affected individuals. A KCNQ2 variant previously associated with benign neonatal seizures is present in 3 of 12 individuals in the severe category. To compare our results with the healthy population, we performed a similar analysis on whole exome sequencing data from 70 Japanese individuals in the 1000 genomes project. Interestingly, the frequency of rare damaging variants in the same set of neuronal excitability genes in healthy individuals is nearly as high as in severely affected individuals. Rather than a single common gene/variant modifying clinical outcome in SCN1A-related epilepsies, our results point to the cumulative effect of rare variants with little to no measurable phenotypic effect (i.e., typical genetic background) unless present in combination with a disease-causing truncation mutation in SCN1A.