Characterizing cell type specific transcriptional differences between the living and postmortem human brain.

Single-nucleus RNA sequencing (snRNA-seq) is often used to define gene expression patterns characteristic of brain cell types as well as to identify cell type specific gene expression signatures of neurological and mental illnesses in postmortem human brains. As methods to obtain brain tissue from living individuals emerge, it is essential to characterize gene expression differences associated with tissue originating from either living or postmortem subjects using snRNA-seq, and to assess whether and how such differences may impact snRNA-seq studies of brain tissue. To address this, human prefrontal cortex single nuclei gene expression was generated and compared between 31 samples from living individuals and 21 postmortem samples. The same cell types were consistently identified in living and postmortem nuclei, though for each cell type, a large proportion of genes were differentially expressed between samples from postmortem and living individuals. Notably, estimation of cell type proportions by cell type deconvolution of pseudo-bulk data was found to be more accurate in samples from living individuals. To allow for future integration of living and postmortem brain gene expression, a model was developed that quantifies from gene expression data the probability a human brain tissue sample was obtained postmortem. These probabilities are established as a means to statistically account for the gene expression differences between samples from living and postmortem individuals. Together, the results presented here provide a deep characterization of both differences between snRNA-seq derived from samples from living and postmortem individuals, as well as qualify and account for their effect on common analyses performed on this type of data.

Mutations in cancer-relevant genes are ubiquitous in histologically normal endometrial tissue.

OBJECTIVE: Endometrial cancer (EndoCA) is the most common gynecologic cancer and incidence and mortality rate continue to increase. Despite well-characterized knowledge of EndoCA-defining mutations, no effective diagnostic or screening tests exist. To lay the foundation for testing development, our study focused on defining the prevalence of somatic mutations present in non-cancerous uterine tissue. METHODS: We obtained ≥8 uterine samplings, including separate endometrial and myometrial layers, from each of 22 women undergoing hysterectomy for non-cancer conditions. We ultra-deep sequenced (>2000× coverage) samples using a 125 cancer-relevant gene panel. RESULTS: All women harbored complex mutation patterns. In total, 308 somatic mutations were identified with mutant allele frequencies ranging up to 96.0%. These encompassed 56 unique mutations from 24 genes. The majority of samples possessed predicted functional cancer mutations but curiously no growth advantage over non-functional mutations was detected. Functional mutations were enriched with increasing patient age (p = 0.045) and BMI (p = 0.0007) and in endometrial versus myometrial layers (68% vs 39%, p = 0.0002). Finally, while the somatic mutation landscape shared similar mutation prevalence in key TCGA-defined EndoCA genes, notably PIK3CA, significant differences were identified, including NOTCH1 (77% vs 10%), PTEN (9% vs 61%), TP53 (0% vs 37%) and CTNNB1 (0% vs 26%). CONCLUSIONS: An important caveat for future liquid biopsy/DNA-based cancer diagnostics is the repertoire of shared and distinct mutation profiles between histologically unremarkable and EndoCA tissues. The lack of selection pressure between functional and non-functional mutations in histologically unremarkable uterine tissue may offer a glimpse into an unrecognized EndoCA protective mechanism.

High-Resolution Profiling of Human Vocal Fold Cellular Landscapes With Single-Nuclei RNA Sequencing.

INTRODUCTION: The function of the vocal folds (VFs) is determined by the phenotype, abundance, and distribution of differentiated cells within specific microenvironments. Identifying this histologic framework is crucial in understanding laryngeal disease. A paucity of studies investigating VF cellular heterogeneity has been undertaken. Here, we examined the cellular landscape of human VFs by utilizing single-nuclei RNA-sequencing. METHODS: Normal true VF tissue was excised from five patients undergoing pitch elevation surgery. Tissue was snap frozen in liquid nitrogen and subjected to cellular digestion and nuclear extraction. Nuclei were processed for single-nucleus sequencing using the 10X Genomics Chromium platform. Sequencing reads were assembled using cellranger and analyzed with the scanpy package in python. RESULTS: RNA sequencing revealed 18 global cell clusters. While many were of epithelial origin, expected cell types, such as fibroblasts, immune cells, muscle cells, and endothelial cells were present. Subcluster analysis defined unique epithelial, immune, and fibroblast subpopulations. CONCLUSION: This study evaluated the cellular heterogeneity of normal human VFs by utilizing single-nuclei RNA-sequencing. With further confirmation through additional spatial sequencing and microscopic imaging, a novel cellular map of the VFs may provide insight into new cellular targets for VF disease. LEVEL OF EVIDENCE: NA Laryngoscope, 2024.

Pan-cancer proteogenomics characterization of tumor immunity.

Despite the successes of immunotherapy in cancer treatment over recent decades, less than <10%-20% cancer cases have demonstrated durable responses from immune checkpoint blockade. To enhance the efficacy of immunotherapies, combination therapies suppressing multiple immune evasion mechanisms are increasingly contemplated. To better understand immune cell surveillance and diverse immune evasion responses in tumor tissues, we comprehensively characterized the immune landscape of more than 1,000 tumors across ten different cancers using CPTAC pan-cancer proteogenomic data. We identified seven distinct immune subtypes based on integrative learning of cell type compositions and pathway activities. We then thoroughly categorized unique genomic, epigenetic, transcriptomic, and proteomic changes associated with each subtype. Further leveraging the deep phosphoproteomic data, we studied kinase activities in different immune subtypes, which revealed potential subtype-specific therapeutic targets. Insights from this work will facilitate the development of future immunotherapy strategies and enhance precision targeting with existing agents.

Molecular profiling of human substantia nigra identifies diverse neuron types associated with vulnerability in Parkinson's disease.

Parkinson's disease (PD) is characterized pathologically by the loss of dopaminergic (DA) neurons in the substantia nigra (SN). Whether cell types beyond DA neurons in the SN show vulnerability in PD remains unclear. Through transcriptomic profiling of 315,867 high-quality single nuclei in the SN from individuals with and without PD, we identified cell clusters representing various neuron types, glia, endothelial cells, pericytes, fibroblasts, and T cells and investigated cell type-dependent alterations in gene expression in PD. Notably, a unique neuron cluster marked by the expression of RIT2, a PD risk gene, also displayed vulnerability in PD. We validated RIT2-enriched neurons in midbrain organoids and the mouse SN. Our results demonstrated distinct transcriptomic signatures of the RIT2-enriched neurons in the human SN and implicated reduced RIT2 expression in the pathogenesis of PD. Our study sheds light on the diversity of cell types, including DA neurons, in the SN and the complexity of molecular and cellular changes associated with PD pathogenesis.

Long-read RNA-seq atlas of novel microglia isoforms elucidates disease-associated genetic regulation of splicing.

Microglia, the innate immune cells of the central nervous system, have been genetically implicated in multiple neurodegenerative diseases. We previously mapped the genetic regulation of gene expression and mRNA splicing in human microglia, identifying several loci where common genetic variants in microglia-specific regulatory elements explain disease risk loci identified by GWAS. However, identifying genetic effects on splicing has been challenging due to the use of short sequencing reads to identify causal isoforms. Here we present the isoform-centric microglia genomic atlas (isoMiGA) which leverages the power of long-read RNA-seq to identify 35,879 novel microglia isoforms. We show that the novel microglia isoforms are involved in stimulation response and brain region specificity. We then quantified the expression of both known and novel isoforms in a multi-ethnic meta-analysis of 555 human microglia short-read RNA-seq samples from 391 donors, the largest to date, and found associations with genetic risk loci in Alzheimer's disease and Parkinson's disease. We nominate several loci that may act through complex changes in isoform and splice site usage.

Sexually dimorphic role of diet and stress on behavior, energy metabolism, and the ventromedial hypothalamus.

Scientific evidence underscores the influence of biological sex on the interplay between stress and metabolic dysfunctions. However, there is limited understanding of how diet and stress jointly contribute to metabolic dysregulation in both males and females. To address this gap, our study aimed to investigate the combined effects of a high-fat diet (HFD) and repeated footshock stress on fear-related behaviors and metabolic outcomes in male and female mice. Using a robust rodent model that recapitulates key aspects of post-traumatic stress disorder (PTSD), we subjected mice to footshock stressor followed by weekly reminder footshock stressor or no stressor for 14 weeks while on either an HFD or chow diet. Our findings revealed that HFD impaired fear memory extinction in male mice that received initial stressor but not in female mice. Blood glucose levels were influenced by both diet and sex, with HFD-fed female mice displaying elevated levels that returned to baseline in the absence of stress, a pattern not observed in male mice. Male mice on HFD exhibited higher energy expenditure, while HFD-fed female mice showed a decreased respiratory exchange ratio (RER). Sex-specific alterations in pro-inflammatory markers and abundance of hematopoietic stem cells were observed in chronically stressed mice on an HFD in different peripheral tissues, indicating the manifestation of distinct comorbid disorders. Single-nuclei RNA sequencing of the ventromedial hypothalamus from stressed mice on an HFD provided insights into sex-specific glial cell activation and cell-type-specific transcriptomic changes. In conclusion, our study offers a comprehensive understanding of the intricate interactions between stress, diet, sex, and various physiological and behavioral outcomes, shedding light on a potential brain region coordinating these interactions.

Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation.

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

Mapping disease regulatory circuits at cell-type resolution from single-cell multiomics data.

Resolving chromatin-remodeling-linked gene expression changes at cell-type resolution is important for understanding disease states. Here we describe MAGICAL (Multiome Accessibility Gene Integration Calling and Looping), a hierarchical Bayesian approach that leverages paired single-cell RNA sequencing and single-cell transposase-accessible chromatin sequencing from different conditions to map disease-associated transcription factors, chromatin sites, and genes as regulatory circuits. By simultaneously modeling signal variation across cells and conditions in both omics data types, MAGICAL achieved high accuracy on circuit inference. We applied MAGICAL to study Staphylococcus aureus sepsis from peripheral blood mononuclear single-cell data that we generated from subjects with bloodstream infection and uninfected controls. MAGICAL identified sepsis-associated regulatory circuits predominantly in CD14 monocytes, known to be activated by bacterial sepsis. We addressed the challenging problem of distinguishing host regulatory circuit responses to methicillin-resistant and methicillin-susceptible S. aureus infections. Although differential expression analysis failed to show predictive value, MAGICAL identified epigenetic circuit biomarkers that distinguished methicillin-resistant from methicillin-susceptible S. aureus infections.

Gemcitabine and cisplatin plus nivolumab as organ-sparing treatment for muscle-invasive bladder cancer: a phase 2 trial.

Cystectomy is a standard treatment for muscle-invasive bladder cancer (MIBC), but it is life-altering. We initiated a phase 2 study in which patients with MIBC received four cycles of gemcitabine, cisplatin, plus nivolumab followed by clinical restaging. Patients achieving a clinical complete response (cCR) could proceed without cystectomy. The co-primary objectives were to assess the cCR rate and the positive predictive value of cCR for a composite outcome: 2-year metastasis-free survival in patients forgoing immediate cystectomy or

Isoform-resolved transcriptome of the human preimplantation embryo.

Human preimplantation development involves extensive remodeling of RNA expression and splicing. However, its transcriptome has been compiled using short-read sequencing data, which fails to capture most full-length mRNAs. Here, we generate an isoform-resolved transcriptome of early human development by performing long- and short-read RNA sequencing on 73 embryos spanning the zygote to blastocyst stages. We identify 110,212 unannotated isoforms transcribed from known genes, including highly conserved protein-coding loci and key developmental regulators. We further identify 17,964 isoforms from 5,239 unannotated genes, which are largely non-coding, primate-specific, and highly associated with transposable elements. These isoforms are widely supported by the integration of published multi-omics datasets, including single-cell 8CLC and blastoid studies. Alternative splicing and gene co-expression network analyses further reveal that embryonic genome activation is associated with splicing disruption and transient upregulation of gene modules. Together, these findings show that the human embryo transcriptome is far more complex than currently known, and will act as a valuable resource to empower future studies exploring development.

Spinal cord injury regulates circular RNA expression in axons.

Introduction: Neurons transport mRNA and translational machinery to axons for local translation. After spinal cord injury (SCI), de novo translation is assumed to enable neurorepair. Knowledge of the identity of axonal mRNAs that participate in neurorepair after SCI is limited. We sought to identify and understand how axonal RNAs play a role in axonal regeneration. Methods: We obtained preparations enriched in axonal mRNAs from control and SCI rats by digesting spinal cord tissue with cold-active protease (CAP). The digested samples were then centrifuged to obtain a supernatant that was used to identify mRNA expression. We identified differentially expressed genes (DEGS) after SCI and mapped them to various biological processes. We validated the DEGs by RT-qPCR and RNA-scope. Results: The supernatant fraction was highly enriched for mRNA from axons. Using Gene Ontology, the second most significant pathway for all DEGs was axonogenesis. Among the DEGs was Rims2, which is predominately a circular RNA (circRNA) in the CNS. We show that Rims2 RNA within spinal cord axons is circular. We found an additional 200 putative circRNAs in the axonal-enriched fraction. Knockdown in primary rat cortical neurons of the RNA editing enzyme ADAR1, which inhibits formation of circRNAs, significantly increased axonal outgrowth and increased the expression of circRims2. Using Rims2 as a prototype we used Circular RNA Interactome to predict miRNAs that bind to circRims2 also bind to the 3'UTR of GAP-43, PTEN or CREB1, all known regulators of axonal outgrowth. Axonally-translated GAP-43 supports axonal elongation and we detect GAP-43 mRNA in the rat axons by RNAscope. Discussion: By enriching for axonal RNA, we detect SCI induced DEGs, including circRNA such as Rims2. Ablation of ADAR1, the enzyme that regulates circRNA formation, promotes axonal outgrowth of cortical neurons. We developed a pathway model using Circular RNA Interactome that indicates that Rims2 through miRNAs can regulate the axonal translation GAP-43 to regulate axonal regeneration. We conclude that axonal regulatory pathways will play a role in neurorepair.

MacroH2A restricts inflammatory gene expression in melanoma cancer-associated fibroblasts by coordinating chromatin looping.

MacroH2A has established tumour suppressive functions in melanoma and other cancers, but an unappreciated role in the tumour microenvironment. Using an autochthonous, immunocompetent mouse model of melanoma, we demonstrate that mice devoid of macroH2A variants exhibit increased tumour burden compared with wild-type counterparts. MacroH2A-deficient tumours accumulate immunosuppressive monocytes and are depleted of functional cytotoxic T cells, characteristics consistent with a compromised anti-tumour response. Single cell and spatial transcriptomics identify increased dedifferentiation along the neural crest lineage of the tumour compartment and increased frequency and activation of cancer-associated fibroblasts following macroH2A loss. Mechanistically, macroH2A-deficient cancer-associated fibroblasts display increased myeloid chemoattractant activity as a consequence of hyperinducible expression of inflammatory genes, which is enforced by increased chromatin looping of their promoters to enhancers that gain H3K27ac. In summary, we reveal a tumour suppressive role for macroH2A variants through the regulation of chromatin architecture in the tumour stroma with potential implications for human melanoma.

Antitumor immunity as the basis for durable disease-free treatment-free survival in patients with metastatic urothelial cancer.

Cisplatin-based chemotherapy has been associated with durable disease control in a small subset of patients with metastatic urothelial cancer. However, the mechanistic basis for this phenomenon has remained elusive. Antitumor immunity may underlie these exceptional responders. In a phase II trial evaluating a phased schedule of gemcitabine and cisplatin followed by gemcitabine and cisplatin with ipilimumab for metastatic urothelial cancer, 4 of 36 patients achieved durable disease-free treatment-free survival (DDFTFS) and remain in remission over 5 years after enrolment on the study. We sought to identify the genomic and immunological mechanisms associated with functional cures of such patients. Whole exome sequencing was performed on pretreatment archival tumor tissue. Neoantigen prediction and ranking were performed using a novel pipeline. For a subset of patients with available biospecimens, selected peptides were tested for neoantigen-specific T cell reactivity in peripheral blood CD4+ and CD8+ T cells cultured with autologous antigen-presenting cells at baseline, postchemotherapy, and postchemotherapy and ipilimumab timepoints. Multiplex assays of serum protein analytes were also assessed at each time point. Serum proteomic analysis revealed that pretreatment, patients achieving DDFTFS demonstrated an immune activated phenotype with elevations in TH1 adaptive immunity, costimulatory molecules, and immune checkpoint markers. After combination cisplatin-based chemotherapy and ipilimumab treatment, DDFTFS patients again displayed enrichment for markers of adaptive immunity, as well as T cell cytotoxicity. CD27 was uniquely enriched in DDFTFS patients at all timepoints. Neoantigen reactivity was not detected in any patient at baseline or post two cycles of chemotherapy. Both CD4+ and CD8+ neoantigen-specific T cell reactivity was detected in two of two DDFTFS patients in comparison to zero of five non-DDFTFS patients after combination cisplatin-based chemotherapy and ipilimumab treatment. Antitumor immunity may underlie functional cures achieved in patients with metastatic urothelial cancer treated with cisplatin-based chemotherapy and immune checkpoint blockade. Probing the mechanistic basis for DDFTFS may facilitate the identification of biomarkers, therapeutic components, and optimal treatment sequences necessary to extend this ultimate goal to a larger subset of patients.

Sequential intrahost evolution and onward transmission of SARS-CoV-2 variants.

Persistent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections have been reported in immune-compromised individuals and people undergoing immune-modulatory treatments. Although intrahost evolution has been documented, direct evidence of subsequent transmission and continued stepwise adaptation is lacking. Here we describe sequential persistent SARS-CoV-2 infections in three individuals that led to the emergence, forward transmission, and continued evolution of a new Omicron sublineage, BA.1.23, over an eight-month period. The initially transmitted BA.1.23 variant encoded seven additional amino acid substitutions within the spike protein (E96D, R346T, L455W, K458M, A484V, H681R, A688V), and displayed substantial resistance to neutralization by sera from boosted and/or Omicron BA.1-infected study participants. Subsequent continued BA.1.23 replication resulted in additional substitutions in the spike protein (S254F, N448S, F456L, M458K, F981L, S982L) as well as in five other virus proteins. Our findings demonstrate not only that the Omicron BA.1 lineage can diverge further from its already exceptionally mutated genome but also that patients with persistent infections can transmit these viral variants. Thus, there is, an urgent need to implement strategies to prevent prolonged SARS-CoV-2 replication and to limit the spread of newly emerging, neutralization-resistant variants in vulnerable patients.

A study of gene expression in the living human brain.

A goal of medical research is to determine the molecular basis of human brain health and illness. One way to achieve this goal is through observational studies of gene expression in human brain tissue. Due to the unavailability of brain tissue from living people, most such studies are performed using tissue from postmortem brain donors. An assumption underlying this practice is that gene expression in the postmortem human brain is an accurate representation of gene expression in the living human brain. Here, this assumption - which, until now, had not been adequately tested - is tested by comparing human prefrontal cortex gene expression between 275 living samples and 243 postmortem samples. Expression levels differed significantly for nearly 80% of genes, and a systematic examination of alternative explanations for this observation determined that these differences are not a consequence of cell type composition, RNA quality, postmortem interval, age, medication, morbidity, symptom severity, tissue pathology, sample handling, batch effects, or computational methods utilized. Analyses integrating the data generated for this study with data from earlier landmark studies that used tissue from postmortem brain donors showed that postmortem brain gene expression signatures of neurological and mental illnesses, as well as of normal traits such as aging, may not be accurate representations of these gene expression signatures in the living brain. By using tissue from large cohorts living people, future observational studies of human brain biology have the potential to (1) determine the medical research questions that can be addressed using postmortem tissue as a proxy for living tissue and (2) expand the scope of medical research to include questions about the molecular basis of human brain health and illness that can only be addressed in living people (e.g., "What happens at the molecular level in the brain as a person experiences an emotion?").

Long-Read DNA Sequencing: Recent Advances and Remaining Challenges.

DNA sequencing has revolutionized medicine over recent decades. However, analysis of large structural variation and repetitive DNA, a hallmark of human genomes, has been limited by short-read technology, with read lengths of 100-300 bp. Long-read sequencing (LRS) permits routine sequencing of human DNA fragments tens to hundreds of kilobase pairs in size, using both real-time sequencing by synthesis and nanopore-based direct electronic sequencing. LRS permits analysis of large structural variation and haplotypic phasing in human genomes and has enabled the discovery and characterization of rare pathogenic structural variants and repeat expansions. It has also recently enabled the assembly of a complete, gapless human genome that includes previously intractable regions, such as highly repetitive centromeres and homologous acrocentric short arms. With the addition of protocols for targeted enrichment, direct epigenetic DNA modification detection, and long-range chromatin profiling, LRS promises to launch a new era of understanding of genetic diversity and pathogenic mutations in human populations.

NF2 mutation associated with accelerated time to recurrence for older patients with atypical meningiomas.

PURPOSE: Meningiomas occur more frequently in older adults, with the incidence rates increasing from 5.8/100,000 for adults 35-44 years old to 55.2/100,000 for those 85+. Due to the increased risk of surgical management in older adults, there is a need to characterize the risk factors for aggressive disease course to inform management decisions in this population. We therefore sought to determine age-stratified relationships between tumour genomics and recurrence after resection of atypical meningiomas. METHODS: We identified 137 primary and recurrent Grade 2 meningiomas from our existing meningioma genomic sequencing database. We examined the differential distribution of genomic alterations in those older than 65 compared to younger. We then performed an age stratified survival analysis to model recurrence for a mutation identified as differentially present. RESULTS: In our cohort of 137 patients with grade 2 meningiomas, alterations in NF2 were present at a higher rate in older adults compared to younger (37.8% in < 65 vs. 55.3% in > 65; recurrence adjusted p-value =0.04). There was no association between the presence of NF2 and recurrence in the whole cohort. In the age-stratified model for those less than 65 years old, there was again no relationship. For patients in the older age stratum, there is a relationship between NF2 and worsened recurrence outcomes (HR = 3.64 (1.125 - 11.811); p = 0.031). CONCLUSIONS: We found that mutations in NF2 were more common in older adults. Further, the presence of mutant NF2 was associated with an increased risk of recurrence in older adults.

FLAIRR-Seq: A Method for Single-Molecule Resolution of Near Full-Length Antibody H Chain Repertoires.

Current Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using short-read sequencing strategies resolve expressed Ab transcripts with limited resolution of the C region. In this article, we present the near-full-length AIRR-seq (FLAIRR-seq) method that uses targeted amplification by 5' RACE, combined with single-molecule, real-time sequencing to generate highly accurate (99.99%) human Ab H chain transcripts. FLAIRR-seq was benchmarked by comparing H chain V (IGHV), D (IGHD), and J (IGHJ) gene usage, complementarity-determining region 3 length, and somatic hypermutation to matched datasets generated with standard 5' RACE AIRR-seq using short-read sequencing and full-length isoform sequencing. Together, these data demonstrate robust FLAIRR-seq performance using RNA samples derived from PBMCs, purified B cells, and whole blood, which recapitulated results generated by commonly used methods, while additionally resolving H chain gene features not documented in IMGT at the time of submission. FLAIRR-seq data provide, for the first time, to our knowledge, simultaneous single-molecule characterization of IGHV, IGHD, IGHJ, and IGHC region genes and alleles, allele-resolved subisotype definition, and high-resolution identification of class switch recombination within a clonal lineage. In conjunction with genomic sequencing and genotyping of IGHC genes, FLAIRR-seq of the IgM and IgG repertoires from 10 individuals resulted in the identification of 32 unique IGHC alleles, 28 (87%) of which were previously uncharacterized. Together, these data demonstrate the capabilities of FLAIRR-seq to characterize IGHV, IGHD, IGHJ, and IGHC gene diversity for the most comprehensive view of bulk-expressed Ab repertoires to date.