The ENCODE mouse postnatal developmental time course identifies regulatory programs of cell types and cell states.

Postnatal genomic regulation significantly influences tissue and organ maturation but is under-studied relative to existing genomic catalogs of adult tissues or prenatal development in mouse. The ENCODE4 consortium generated the first comprehensive single-nucleus resource of postnatal regulatory events across a diverse set of mouse tissues. The collection spans seven postnatal time points, mirroring human development from childhood to adulthood, and encompasses five core tissues. We identified 30 cell types, further subdivided into 69 subtypes and cell states across adrenal gland, left cerebral cortex, hippocampus, heart, and gastrocnemius muscle. Our annotations cover both known and novel cell differentiation dynamics ranging from early hippocampal neurogenesis to a new sex-specific adrenal gland population during puberty. We used an ensemble Latent Dirichlet Allocation strategy with a curated vocabulary of 2,701 regulatory genes to identify regulatory "topics," each of which is a gene vector, linked to cell type differentiation, subtype specialization, and transitions between cell states. We find recurrent regulatory topics in tissue-resident macrophages, neural cell types, endothelial cells across multiple tissues, and cycling cells of the adrenal gland and heart. Cell-type-specific topics are enriched in transcription factors and microRNA host genes, while chromatin regulators dominate mitosis topics. Corresponding chromatin accessibility data reveal dynamic and sex-specific regulatory elements, with enriched motifs matching transcription factors in regulatory topics. Together, these analyses identify both tissue-specific and common regulatory programs in postnatal development across multiple tissues through the lens of the factors regulating transcription.

BIN1K358R suppresses glial response to plaques in mouse model of Alzheimer's disease.

INTRODUCTION: The BIN1 coding variant rs138047593 (K358R) is linked to Late-Onset Alzheimer's Disease (LOAD) via targeted exome sequencing. METHODS: To elucidate the functional consequences of this rare coding variant on brain amyloidosis and neuroinflammation, we generated BIN1K358R knock-in mice using CRISPR/Cas9 technology. These mice were subsequently bred with 5xFAD transgenic mice, which serve as a model for Alzheimer's pathology. RESULTS: The presence of the BIN1K358R variant leads to increased cerebral amyloid deposition, with a dampened response of astrocytes and oligodendrocytes, but not microglia, at both the cellular and transcriptional levels. This correlates with decreased neurofilament light chain in both plasma and brain tissue. Synaptic densities are significantly increased in both wild-type and 5xFAD backgrounds homozygous for the BIN1K358R variant. DISCUSSION: The BIN1 K358R variant modulates amyloid pathology in 5xFAD mice, attenuates the astrocytic and oligodendrocytic responses to amyloid plaques, decreases damage markers, and elevates synaptic densities. HIGHLIGHTS: BIN1 rs138047593 (K358R) coding variant is associated with increased risk of LOAD. BIN1 K358R variant increases amyloid plaque load in 12-month-old 5xFAD mice. BIN1 K358R variant dampens astrocytic and oligodendrocytic response to plaques. BIN1 K358R variant decreases neuronal damage in 5xFAD mice. BIN1 K358R upregulates synaptic densities and modulates synaptic transmission.

The Abca7V1613M variant reduces Aß generation, plaque load, and neuronal damage.

BACKGROUND: Variants in ABCA7, a member of the ABC transporter superfamily, have been associated with increased risk for developing late onset Alzheimer's disease (LOAD). METHODS: CRISPR-Cas9 was used to generate an Abca7V1613M variant in mice, modeling the homologous human ABCA7V1599M variant, and extensive characterization was performed. RESULTS: Abca7V1613M microglia show differential gene expression profiles upon lipopolysaccharide challenge and increased phagocytic capacity. Homozygous Abca7V1613M mice display elevated circulating cholesterol and altered brain lipid composition. When crossed with 5xFAD mice, homozygous Abca7V1613M mice display fewer Thioflavin S-positive plaques, decreased amyloid beta (Aß) peptides, and altered amyloid precursor protein processing and trafficking. They also exhibit reduced Aß-associated inflammation, gliosis, and neuronal damage. DISCUSSION: Overall, homozygosity for the Abca7V1613M variant influences phagocytosis, response to inflammation, lipid metabolism, Aß pathology, and neuronal damage in mice. This variant may confer a gain of function and offer a protective effect against Alzheimer's disease-related pathology. HIGHLIGHTS: ABCA7 recognized as a top 10 risk gene for developing Alzheimer's disease. Loss of function mutations result in increased risk for LOAD. V1613M variant reduces amyloid beta plaque burden in 5xFAD mice. V1613M variant modulates APP processing and trafficking in 5xFAD mice. V1613M variant reduces amyloid beta-associated damage in 5xFAD mice.

Identification of robust cellular programs using reproducible LDA that impact sex-specific disease progression in different genotypes of a mouse model of AD.

The gene expression profiles of distinct cell types reflect complex genomic interactions among multiple simultaneous biological processes within each cell that can be altered by disease progression as well as genetic background. The identification of these active cellular programs is an open challenge in the analysis of single-cell RNA-seq data. Latent Dirichlet Allocation (LDA) is a generative method used to identify recurring patterns in counts data, commonly referred to as topics that can be used to interpret the state of each cell. However, LDA's interpretability is hindered by several key factors including the hyperparameter selection of the number of topics as well as the variability in topic definitions due to random initialization. We developed Topyfic, a Reproducible LDA (rLDA) package, to accurately infer the identity and activity of cellular programs in single-cell data, providing insights into the relative contributions of each program in individual cells. We apply Topyfic to brain single-cell and single-nucleus datasets of two 5xFAD mouse models of Alzheimer's disease crossed with C57BL6/J or CAST/EiJ mice to identify distinct cell types and states in different cell types such as microglia. We find that 8-month 5xFAD/Cast F1 males show higher level of microglial activation than matching 5xFAD/BL6 F1 males, whereas female mice show similar levels of microglial activation. We show that regulatory genes such as TFs, microRNA host genes, and chromatin regulatory genes alone capture cell types and cell states. Our study highlights how topic modeling with a limited vocabulary of regulatory genes can identify gene expression programs in single-cell data in order to quantify similar and divergent cell states in distinct genotypes.

C5aR1 antagonism suppresses inflammatory glial gene expression and alters cellular signaling in an aggressive Alzheimer's model.

Alzheimer's disease (AD) is the leading cause of dementia in older adults, and the need for effective, sustainable therapeutic targets is imperative. Pharmacologic inhibition of C5aR1 reduces plaque load, gliosis and memory deficits in animal models. However, the cellular basis underlying this neuroprotection and which processes were the consequence of amyloid reduction vs alteration of the response to amyloid were unclear. In the Arctic model, the C5aR1 antagonist PMX205 did not reduce plaque load, but deficits in short-term memory in female mice were prevented. Hippocampal single cell and single nucleus RNA-seq clusters revealed C5aR1 dependent and independent gene expression and cell-cell communication. Microglial clusters containing neurotoxic disease-associated microglial genes were robustly upregulated in Arctic mice and drastically reduced with PMX205 treatment, while genes in microglia clusters that were overrepresented in the Arctic-PMX205 vs Arctic group were associated with synapse organization and transmission and learning. PMX205 treatment also reduced some A-1 astrocyte genes. In spite of changes in transcript levels, overall protein levels of some reactive glial markers were relatively unchanged by C5aR1 antagonism, as were clusters associated with protective responses to injury. C5aR1 inhibition promoted signaling pathways associated with cell growth and repair, such as TGFß and FGF, in Arctic mice, while suppressing inflammatory pathways including PROS, Pecam1, and EPHA. In conclusion, pharmacologic C5aR1 inhibition prevents cognitive loss, limits microglial polarization to a detrimental inflammatory state and permits neuroprotective responses, as well as leaving protective functions of complement intact, making C5aR1 antagonism an attractive therapeutic strategy for individuals with AD.

Signalling by senescent melanocytes hyperactivates hair growth.

Niche signals maintain stem cells in a prolonged quiescence or transiently activate them for proper regeneration1. Altering balanced niche signalling can lead to regenerative disorders. Melanocytic skin nevi in human often display excessive hair growth, suggesting hair stem cell hyperactivity. Here, using genetic mouse models of nevi2,3, we show that dermal clusters of senescent melanocytes drive epithelial hair stem cells to exit quiescence and change their transcriptome and composition, potently enhancing hair renewal. Nevus melanocytes activate a distinct secretome, enriched for signalling factors. Osteopontin, the leading nevus signalling factor, is both necessary and sufficient to induce hair growth. Injection of osteopontin or its genetic overexpression is sufficient to induce robust hair growth in mice, whereas germline and conditional deletions of either osteopontin or CD44, its cognate receptor on epithelial hair cells, rescue enhanced hair growth induced by dermal nevus melanocytes. Osteopontin is overexpressed in human hairy nevi, and it stimulates new growth of human hair follicles. Although broad accumulation of senescent cells, such as upon ageing or genotoxic stress, is detrimental for the regenerative capacity of tissue4, we show that signalling by senescent cell clusters can potently enhance the activity of adjacent intact stem cells and stimulate tissue renewal. This finding identifies senescent cells and their secretome as an attractive therapeutic target in regenerative disorders.

The ENCODE4 long-read RNA-seq collection reveals distinct classes of transcript structure diversity.

The majority of mammalian genes encode multiple transcript isoforms that result from differential promoter use, changes in exonic splicing, and alternative 3' end choice. Detecting and quantifying transcript isoforms across tissues, cell types, and species has been extremely challenging because transcripts are much longer than the short reads normally used for RNA-seq. By contrast, long-read RNA-seq (LR-RNA-seq) gives the complete structure of most transcripts. We sequenced 264 LR-RNA-seq PacBio libraries totaling over 1 billion circular consensus reads (CCS) for 81 unique human and mouse samples. We detect at least one full-length transcript from 87.7% of annotated human protein coding genes and a total of 200,000 full-length transcripts, 40% of which have novel exon junction chains. To capture and compute on the three sources of transcript structure diversity, we introduce a gene and transcript annotation framework that uses triplets representing the transcript start site, exon junction chain, and transcript end site of each transcript. Using triplets in a simplex representation demonstrates how promoter selection, splice pattern, and 3' processing are deployed across human tissues, with nearly half of multi-transcript protein coding genes showing a clear bias toward one of the three diversity mechanisms. Evaluated across samples, the predominantly expressed transcript changes for 74% of protein coding genes. In evolution, the human and mouse transcriptomes are globally similar in types of transcript structure diversity, yet among individual orthologous gene pairs, more than half (57.8%) show substantial differences in mechanism of diversification in matching tissues. This initial large-scale survey of human and mouse long-read transcriptomes provides a foundation for further analyses of alternative transcript usage, and is complemented by short-read and microRNA data on the same samples and by epigenome data elsewhere in the ENCODE4 collection.

A Trem2R47H mouse model without cryptic splicing drives age- and disease-dependent tissue damage and synaptic loss in response to plaques.

BACKGROUND: The TREM2 R47H variant is one of the strongest genetic risk factors for late-onset Alzheimer's Disease (AD). Unfortunately, many current Trem2 R47H mouse models are associated with cryptic mRNA splicing of the mutant allele that produces a confounding reduction in protein product. To overcome this issue, we developed the Trem2R47H NSS (Normal Splice Site) mouse model in which the Trem2 allele is expressed at a similar level to the wild-type Trem2 allele without evidence of cryptic splicing products. METHODS: Trem2R47H NSS mice were treated with the demyelinating agent cuprizone, or crossed with the 5xFAD mouse model of amyloidosis, to explore the impact of the TREM2 R47H variant on inflammatory responses to demyelination, plaque development, and the brain's response to plaques. RESULTS: Trem2R47H NSS mice display an appropriate inflammatory response to cuprizone challenge, and do not recapitulate the null allele in terms of impeded inflammatory responses to demyelination. Utilizing the 5xFAD mouse model, we report age- and disease-dependent changes in Trem2R47H NSS mice in response to development of AD-like pathology. At an early (4-month-old) disease stage, hemizygous 5xFAD/homozygous Trem2R47H NSS (5xFAD/Trem2R47H NSS) mice have reduced size and number of microglia that display impaired interaction with plaques compared to microglia in age-matched 5xFAD hemizygous controls. This is associated with a suppressed inflammatory response but increased dystrophic neurites and axonal damage as measured by plasma neurofilament light chain (NfL) level. Homozygosity for Trem2R47H NSS suppressed LTP deficits and loss of presynaptic puncta caused by the 5xFAD transgene array in 4-month-old mice. At a more advanced (12-month-old) disease stage 5xFAD/Trem2R47H NSS mice no longer display impaired plaque-microglia interaction or suppressed inflammatory gene expression, although NfL levels remain elevated, and a unique interferon-related gene expression signature is seen. Twelve-month old Trem2R47H NSS mice also display LTP deficits and postsynaptic loss. CONCLUSIONS: The Trem2R47H NSS mouse is a valuable model that can be used to investigate age-dependent effects of the AD-risk R47H mutation on TREM2 and microglial function including its effects on plaque development, microglial-plaque interaction, production of a unique interferon signature and associated tissue damage.

C5aR1 antagonism alters microglial polarization and mitigates disease progression in a mouse model of Alzheimer's disease.

Multiple studies have recognized the involvement of the complement cascade during Alzheimer's disease pathogenesis. However, the specific role of C5a-C5aR1 signaling in the progression of this neurodegenerative disease is still not clear. Furthermore, its potential as a therapeutic target to treat AD still remains to be elucidated. Canonically, generation of the anaphylatoxin C5a as the result of complement activation and interaction with its receptor C5aR1 triggers a potent inflammatory response. Previously, genetic ablation of C5aR1 in a mouse model of Alzheimer's disease exerted a protective effect by preventing cognitive deficits. Here, using PMX205, a potent, specific C5aR1 antagonist, in the Tg2576 mouse model of Alzheimer's disease we show a striking reduction in dystrophic neurites in parallel with the reduced amyloid load, rescue of the excessive pre-synaptic loss associated with AD cognitive impairment and the polarization of microglial gene expression towards a DAM-like phenotype that are consistent with the neuroprotective effects seen. These data support the beneficial effect of a pharmacological inhibition of C5aR1 as a promising therapeutic approach to treat Alzheimer's disease. Supportive of the safety of this treatment is the recent FDA-approval of another other C5a receptor 1 antagonist, Avacopan, as a treatment for autoimmune inflammatory diseases.

Modulation of C5a-C5aR1 signaling alters the dynamics of AD progression.

BACKGROUND: The complement system is part of the innate immune system that clears pathogens and cellular debris. In the healthy brain, complement influences neurodevelopment and neurogenesis, synaptic pruning, clearance of neuronal blebs, recruitment of phagocytes, and protects from pathogens. However, excessive downstream complement activation that leads to generation of C5a, and C5a engagement with its receptor C5aR1, instigates a feed-forward loop of inflammation, injury, and neuronal death, making C5aR1 a potential therapeutic target for neuroinflammatory disorders. C5aR1 ablation in the Arctic (Arc) model of Alzheimer's disease protects against cognitive decline and neuronal injury without altering amyloid plaque accumulation. METHODS: To elucidate the effects of C5a-C5aR1 signaling on AD pathology, we crossed Arc mice with a C5a-overexpressing mouse (ArcC5a+) and tested hippocampal memory. RNA-seq was performed on hippocampus and cortex from Arc, ArcC5aR1KO, and ArcC5a+ mice at 2.7-10 months and age-matched controls to assess mechanisms involved in each system. Immunohistochemistry was used to probe for protein markers of microglia and astrocytes activation states. RESULTS: ArcC5a+ mice had accelerated cognitive decline compared to Arc. Deletion of C5ar1 delayed or prevented the expression of some, but not all, AD-associated genes in the hippocampus and a subset of pan-reactive and A1 reactive astrocyte genes, indicating a separation between genes induced by amyloid plaques alone and those influenced by C5a-C5aR1 signaling. Biological processes associated with AD and AD mouse models, including inflammatory signaling, microglial cell activation, and astrocyte migration, were delayed in the ArcC5aR1KO hippocampus. Interestingly, C5a overexpression also delayed the increase of some AD-, complement-, and astrocyte-associated genes, suggesting the possible involvement of neuroprotective C5aR2. However, these pathways were enhanced in older ArcC5a+ mice compared to Arc. Immunohistochemistry confirmed that C5a-C5aR1 modulation in Arc mice delayed the increase in CD11c-positive microglia, while not affecting other pan-reactive microglial or astrocyte markers. CONCLUSION: C5a-C5aR1 signaling in AD largely exerts its effects by enhancing microglial activation pathways that accelerate disease progression. While C5a may have neuroprotective effects via C5aR2, engagement of C5a with C5aR1 is detrimental in AD models. These data support specific pharmacological inhibition of C5aR1 as a potential therapeutic strategy to treat AD.

Mapping and modeling the genomic basis of differential RNA isoform expression at single-cell resolution with LR-Split-seq.

The rise in throughput and quality of long-read sequencing should allow unambiguous identification of full-length transcript isoforms. However, its application to single-cell RNA-seq has been limited by throughput and expense. Here we develop and characterize long-read Split-seq (LR-Split-seq), which uses combinatorial barcoding to sequence single cells with long reads. Applied to the C2C12 myogenic system, LR-split-seq associates isoforms to cell types with relative economy and design flexibility. We find widespread evidence of changing isoform expression during differentiation including alternative transcription start sites (TSS) and/or alternative internal exon usage. LR-Split-seq provides an affordable method for identifying cluster-specific isoforms in single cells.