Scientific highlights

The widespread protozoan Toxoplasma gondii chronically infects neural tissue in vertebrates and is linked to various neurological and neuropsychiatric disorders in humans. However, its effects on sparsely infected neurons and on broader neural circuits remain elusive. Our study reveals that T. gondii infection disrupts cytoskeletal dynamics in SH-SY5Y neuronal cells and primary cortical neurons. Infected neuronal cells undergo significant cytomorphological changes, including retraction of dendritic extensions and alterations in microtubule and F-actin networks, across both parasite genotypes I and II. These cytoskeletal alterations were notably diminished in cells exposed to T. gondii mutants with impaired secretion via the MYR translocon, and were independent of intraneuronal parasite replication. Moreover, a bystander effect was observed, with supernatants from T. gondii-challenged cells inducing similar cytoskeletal changes in uninfected cells. Analyses of extracellular vesicles (EVs) in supernatants revealed differential expression of host microRNAs in response to infection, most notably the upregulation of miR-221-3p, a microRNA not previously associated with T. gondii. The data indicate that unidentified parasite-derived effector(s) secreted via the MYR translocon, in conjunction with MYR-independently induced EV-associated host microRNAs, mediate cytoskeletal alterations in both infected and bystander neuronal cells. The findings provide new insights into molecular mechanisms by which T. gondii infection may disrupt neural networks, shedding light on its potential role in neuronal dysregulation.

Sterile inflammation after injury is important for tissue restoration. In injured human and mouse tissues, macrophages were recently found to accumulate perivascularly. This study investigates if macrophages adopt a mural cell phenotype important for restoration after ischemic injury. Single-cell RNA sequencing of fate-mapped macrophages from ischemic mouse muscles demonstrates a macrophage-toward-mural cell switch of a subpopulation of macrophages with downregulated myeloid cell genes and upregulated mural cell genes, including PDGFRβ. This observation was further strengthened when including unspliced transcripts in the analysis. The macrophage switch was proven functionally relevant, as induction of macrophage-specific PDGFRβ deficiency prevented their perivascular macrophage phenotype, impaired vessel maturation and increased vessel leakiness, which ultimately reduced limb function. In conclusion, macrophages in adult ischemic tissue were demonstrated to undergo a cellular program to morphologically, transcriptomically and functionally resemble mural cells while weakening their macrophage identity. The macrophage-to-mural cell-like phenotypic switch is crucial for restoring tissue function and warrants further exploration as a potential target for immunotherapies to enhance healing.

B cells are important in the pathogenesis of primary Sjögren's syndrome (pSS). Patients positive for Sjögren's syndrome antigen A/Sjögren syndrome antigen B (SSA/SSB) autoantibodies are more prone to systemic disease manifestations and adverse outcomes. We aimed to determine the role of B cell composition, gene expression, and B cell receptor usage in pSS subgroups stratified for SSA/SSB antibodies. Over 230,000 B cells were isolated from peripheral blood of patients with pSS (n = 6 SSA-, n = 8 SSA+ single positive and n = 10 SSA/SSB+ double positive) and four healthy controls and processed for single-cell RNA sequencing (scRNA-seq) and single-cell variable, diversity, and joining (VDJ) gene sequencing (scVDJ-seq). We show that SSA/SSB+ patients present the highest and lowest proportion of naïve and memory B cells, respectively, and the highest up-regulation of interferon-induced genes across all B cell subtypes. Differential usage of IGHV showed that IGHV1-69 and IGHV4-30-4 were more often used in all pSS subgroups compared with controls. Memory B cells from SSA/SSB+ patients displayed a higher proportion of cells with unmutated VDJ transcripts compared with other pSS patient groups and controls, indicating altered somatic hypermutation processes. Comparison with previous studies revealed heterogeneous clonotype pools, with little overlap in CDR3 sequences. Joint analysis using scRNA-seq and scVDJ-seq data allowed unsupervised stratification of patients with pSS and identified novel parameters that correlated to disease manifestations and antibody status. We describe heterogeneity and molecular characteristics in B cells from patients with pSS, providing clues to intrinsic differences in B cells that affect the phenotype and outcome and allowing stratification of patients with pSS at improved resolution.

A wealth of specialized cell populations within the skin facilitates its hair-producing, protective, sensory, and thermoregulatory functions. How the vast cell-type diversity and tissue architecture develops is largely unexplored. Here, with single-cell transcriptomics, spatial cell-type assignment, and cell-lineage tracing, we deconstruct early embryonic mouse skin during the key transitions from seemingly uniform developmental precursor states to a multilayered, multilineage epithelium, and complex dermal identity. We identify the spatiotemporal emergence of hair-follicle-inducing, muscle-supportive, and fascia-forming fibroblasts. We also demonstrate the formation of the panniculus carnosus muscle (PCM), sprouting blood vessels without pericyte coverage, and the earliest residence of mast and dendritic immune cells in skin. Finally, we identify an unexpected epithelial heterogeneity within the early single-layered epidermis and a signaling-rich periderm layer. Overall, this cellular and molecular blueprint of early skin development-which can be explored at https://kasperlab.org/tools-establishes histological landmarks and highlights unprecedented dynamic interactions among skin cells.

Multiple myeloma (MM) is an incurable and aggressive plasma cell malignancy characterized by a complex karyotype with multiple structural variants (SVs) and copy-number variations (CNVs). Linked-read whole-genome sequencing (lrWGS) allows for refined detection and reconstruction of SVs by providing long-range genetic information from standard short-read sequencing. This makes lrWGS an attractive solution for capturing the full genomic complexity of MM. Here we show that high-quality lrWGS data can be generated from low numbers of cells subjected to fluorescence-activated cell sorting (FACS) without DNA purification. Using this protocol, we analyzed MM cells after FACS from 37 patients with MM using lrWGS. We found high concordance between lrWGS and fluorescence in situ hybridization (FISH) for the detection of recurrent translocations and CNVs. Outside of the regions investigated by FISH, we identified >150 additional SVs and CNVs across the cohort. Analysis of the lrWGS data allowed for resolution of the structure of diverse SVs affecting the MYC and t(11;14) loci, causing the duplication of genes and gene regulatory elements. In addition, we identified private SVs causing the dysregulation of genes recurrently involved in translocations with the IGH locus and show that these can alter the molecular classification of MM. Overall, we conclude that lrWGS allows for the detection of aberrations critical for MM prognostics and provides a feasible route for providing comprehensive genetics. Implementing lrWGS could provide more accurate clinical prognostics, facilitate genomic medicine initiatives, and greatly improve the stratification of patients included in clinical trials.

Woolly mammoths had a set of adaptations that enabled them to thrive in the Arctic environment. Many mammoth-specific single nucleotide polymorphisms (SNPs) responsible for unique mammoth traits have been previously identified from ancient genomes. However, a multitude of other genetic variants likely contributed to woolly mammoth evolution. In this study, we sequenced two woolly mammoth genomes and combined these with previously sequenced mammoth and elephant genomes to conduct a survey of mammoth-specific deletions and indels. We find that deletions are highly enriched in non-coding regions, suggesting selection against structural variants that affect protein sequences. Nonetheless, at least 87 woolly mammoth genes contain deletions or indels that modify the coding sequence, including genes involved in skeletal morphology and hair growth. These results suggest that deletions and indels contributed to the unique phenotypic adaptations of the woolly mammoth, and were potentially critical to surviving in its natural environment.

Spinal cord ependymal cells display neural stem cell properties in vitro and generate scar-forming astrocytes and remyelinating oligodendrocytes after injury. We report that ependymal cells are functionally heterogeneous and identify a small subpopulation (8% of ependymal cells and 0.1% of all cells in a spinal cord segment), which we denote ependymal A (EpA) cells, that accounts for the in vitro stem cell potential in the adult spinal cord. After spinal cord injury, EpA cells undergo self-renewing cell division as they give rise to differentiated progeny. Single-cell transcriptome analysis revealed a loss of ependymal cell gene expression programs as EpA cells gained signaling entropy and dedifferentiated to a stem-cell-like transcriptional state after an injury. We conclude that EpA cells are highly differentiated cells that can revert to a stem cell state and constitute a therapeutic target for spinal cord repair.

Molecular mechanisms that define patterns of neuropeptide expression are essential for the formation and rewiring of neural circuits. The prodynorphin gene (PDYN) gives rise to dynorphin opioid peptides mediating depression and substance dependence. We here demonstrated that PDYN is expressed in neurons in human dorsolateral prefrontal cortex (dlPFC), and identified neuronal differentially methylated region in PDYN locus framed by CCCTC-binding factor binding sites. A short, nucleosome size human-specific promoter CpG island (CGI), a core of this region may serve as a regulatory module, which is hypomethylated in neurons, enriched in 5-hydroxymethylcytosine, and targeted by USF2, a methylation-sensitive E-box transcription factor (TF). USF2 activates PDYN transcription in model systems, and binds to nonmethylated CGI in dlPFC. USF2 and PDYN expression is correlated, and USF2 and PDYN proteins are co-localized in dlPFC. Segregation of activatory TF and repressive CGI methylation may ensure contrasting PDYN expression in neurons and glia in human brain.

Pdgfra+ oligodendrocyte precursor cells (OPCs) arise in distinct specification waves during embryogenesis in the central nervous system (CNS). It is unclear whether there is a correlation between these waves and different oligodendrocyte (OL) states at adult stages. Here, we present bulk and single-cell transcriptomics resources providing insights on how transitions between these states occur. We found that post-natal OPCs from brain and spinal cord present similar transcriptional signatures. Moreover, post-natal OPC progeny of E13.5 Pdgfra+ cells present electrophysiological and transcriptional profiles similar to OPCs derived from subsequent specification waves, indicating that Pdgfra+ pre-OPCs rewire their transcriptional network during development. Single-cell RNA-seq and lineage tracing indicates that a subset of E13.5 Pdgfra+ cells originates cells of the pericyte lineage. Thus, our results indicate that embryonic Pdgfra+ cells in the CNS give rise to distinct post-natal cell lineages, including OPCs with convergent transcriptional profiles in different CNS regions.

Single-cell RNA-seq has become routine for discovering cell types and revealing cellular diversity, but archived human brain samples still pose a challenge to current high-throughput platforms. We present STRT-seq-2i, an addressable 9600-microwell array platform, combining sampling by limiting dilution or FACS, with imaging and high throughput at competitive cost. We applied the platform to fresh single mouse cortical cells and to frozen post-mortem human cortical nuclei, matching the performance of a previous lower-throughput platform while retaining a high degree of flexibility, potentially also for other high-throughput applications.

Dietary factors may regulate the epigenome. We aimed to explore whether a diet intervention, including excess sugar, affects the methylome in human sperm, and to describe the sperm methylome. We used Whole Genome Bisulfite Sequencing (WGBS) to analyze DNA methylation in sperm taken at three time points from 15 males during a diet intervention; i) at baseline, ii) after one week on a standardized diet, and iii) after an additional week on a high-sugar diet providing 150% of their estimated total energy expenditure. We identified seven nominal diet-associated differentially methylated regions in sperm (p < 0.05). The diet was nominally associated with methylation of 143 sites linked to fertility (e.g. AHRR, GNAS, and HDAC4), 313 sites in imprinted genes (e.g. GLIS3, PEG10, PEG3, and SNURF), and 42 sites in top 1%-expressed genes (e.g. CHD2) (p < 0.05). In sperm, 3'UTRs and introns had the highest levels of methylation, while 5'UTRs and CpG islands had the lowest levels. Non-expressed genes in human sperm were hypomethylated in exons compared with transcribed genes. In human sperm, DNA methylation levels were linked to gene expression, and excess sugar had modest effects on methylation on imprinted and highly expressed genes, and genes affecting fertility.

Cancer is the leading cause of disease-related mortality in children. Causes of leukemia, the most common form, are largely unknown. Growing evidence points to an origin in-utero, when global redistribution of DNA methylation occurs driving tissue differentiation. Epigenome-wide DNA methylation was profiled in surrogate (blood) and target (bone marrow) tissues at birth, diagnosis, remission and relapse of pediatric pre-B acute lymphoblastic leukemia (pre-B ALL) patients. Double-blinded analyses was performed between prospective cohorts extending from birth to diagnosis and retrospective studies backtracking from clinical disease to birth. Validation was carried out using independent technologies and populations. The imprinted and immuno-modulating VTRNA2-1 was hypermethylated (FDR<0.05) at birth in nested cases relative to controls in all tested populations (totaling 317 cases and 483 controls), including European and Hispanic ancestries. VTRNA2-1 methylation was stable over follow-up years after birth and across surrogate, target and other tissues (n=5,023 tissues; 30 types). When profiled in leukemic tissues from two clinical cohorts (totaling 644 cases), VTRNA2-1 methylation exhibited higher levels at diagnosis relative to controls, it reset back to normal levels at remission, and then re-increased to above control levels at relapse. Hypermethylation was significantly associated with worse pre-B ALL patient survival and with reduced VTRNA2-1 expression (n=2,294 tissues; 26 types), supporting a functional and translational role for VTRNA2-1 methylation. This study provides proof-of-concept to detect at birth epigenetic precursors of pediatric pre-B ALL. These alterations were reproducible with different technologies, in three continents and in two ethnicities, and can offer biomarkers for early detection and prognosis as well as actionable targets for therapy.

The cell transition from an inflammatory phase to a subsequent proliferative phase is crucial for wound healing, yet the driving mechanism remains unclear. By profiling lncRNA expression changes during human skin wound healing and screening lncRNA functions, we identify SNHG26 as a pivotal regulator in keratinocyte progenitors underpinning this phase transition. Snhg26-deficient mice exhibit impaired wound repair characterized by delayed re-epithelization accompanied by exacerbated inflammation. Single-cell transcriptome analysis combined with gain-of-function and loss-of-function of SNHG26 in vitro and ex vivo reveals its specific role in facilitating inflammatory-to-proliferative state transition of keratinocyte progenitors. A mechanistic study unravels that SNHG26 interacts with and relocates the transcription factor ILF2 from inflammatory genomic loci, such as JUN, IL6, IL8, and CCL20, to the genomic locus of LAMB3. Collectively, our findings suggest that lncRNAs play cardinal roles in expediting tissue repair and regeneration and may constitute an invaluable reservoir of therapeutic targets in reparative medicine.

Sheep are among the earliest domesticated livestock species, with a wide variety of breeds present today. However, it remains unclear how far back this diversity goes, with formal documentation only dating back a few centuries. North European short-tailed (NEST) breeds are often assumed to be among the oldest domestic sheep populations, even thought to represent relicts of the earliest sheep expansions during the Neolithic period reaching Scandinavia <6,000 years ago. This study sequenced the genomes (up to 11.6X) of five sheep remains from the Baltic islands of Gotland and Åland, dating from the Late Neolithic (∼4,100 cal BP) to historical times (∼1,600 CE). Our findings indicate that these ancient sheep largely possessed the genetic characteristics of modern NEST breeds, suggesting a substantial degree of long-term continuity of this sheep type in the Baltic Sea region. Despite the wide temporal spread, population genetic analyses show high levels of affinity between the ancient genomes and they also exhibit relatively high genetic diversity when compared to modern NEST breeds, implying a loss of diversity in most breeds during the last centuries associated with breed formation and recent bottlenecks. Our results shed light on the development of breeds in Northern Europe specifically as well as the development of genetic diversity in sheep breeds, and their expansion from the domestication center in general.

Persistent taxane-induced peripheral neuropathy (TIPN) is highly prevalent among early-stage breast cancer survivors (ESBCS) and has detrimental effect on quality of life. We leveraged logistic regression models to develop and validate polygenic prediction models to estimate the risk of persistent PN symptoms in a training cohort and validation cohort taking clinical risk factors into account. Based on 337 whole-exome sequenced ESBCS two of five prediction models for individual PN symptoms obtained AUC results above 60% when validated. Using the model for numbness in feet (35 SNVs) in the test cohort, 73% survivors were correctly predicted. For tingling in feet (55 SNVs) 70% were correctly predicted. Both models included SNVs from the ADAMTS20, APT6V0A2, CCDC88C, CYP2C8, EPHA5, NR1H3, PSKH2/APTV0D2, and SCN10A genes. For cramps in feet, difficulty climbing stairs and difficulty opening a jar the validation was unsuccessful. Polygenic prediction models including clinical risk factors can estimate the risk of persistent taxane-induced numbness in feet and tingling in feet in ESBCS.

Understanding the underlying mechanisms and links between genome evolution and adaptive innovations stands as a key goal in evolutionary studies. Poplars, among the world's most widely distributed and cultivated trees, exhibit extensive phenotypic diversity and environmental adaptability. In this study, we present a genus-level super-pangenome comprising 19 Populus genomes, revealing the likely pivotal role of private genes in facilitating local environmental and climate adaptation. Through the integration of pangenomes with transcriptomes, methylomes, and chromatin accessibility mapping, we unveil that the evolutionary trajectories of pangenes and duplicated genes are closely linked to local genomic landscapes of regulatory and epigenetic architectures, notably CG methylation in gene-body regions. Further comparative genomic analyses have enabled the identification of 142 202 structural variants across species that intersect with a significant number of genes and contribute substantially to both phenotypic and adaptive divergence. We have experimentally validated a ∼180-bp presence/absence variant affecting the expression of the CUC2 gene, crucial for leaf serration formation. Finally, we developed a user-friendly web-based tool encompassing the multi-omics resources associated with the Populus super-pangenome (http://www.populus-superpangenome.com). Together, the present pioneering super-pangenome resource in forest trees not only aids in the advancement of breeding efforts of this globally important tree genus but also offers valuable insights into potential avenues for comprehending tree biology.

The caecum bridges the small and large intestine and plays a front-line role in discriminating gastrointestinal antigens. Although dysregulated in acute and chronic conditions, the tissue is often overlooked immunologically. To address this issue, we applied single-cell transcriptomic-V(D)J sequencing to FACS-isolated CD45+ caecal patch/lamina propria leukocytes from a healthy (5-year-old) female rhesus macaque ex vivo and coupled these data to VDJ deep sequencing reads from haematopoietic tissues. We found caecal NK cells and ILC3s to co-exist with a spectrum of effector T cells partially derived from SOX4 + recent thymic emigrants. Tolerogenic Vγ8Vδ1-T cells, plastic CD4+ T helper cells and GZMK + EOMES + and TMIGD2 + tissue-resident memory CD8+ T cells were present and differed metabolically. An IL13 + GATA3 + Th2 subset expressing eicosanoid pathway enzymes was accompanied by IL1RL1 + GATA3 + regulatory T cells and a minor proportion of IgE+ plasma cells (PCs), illustrating tightly regulated type 2 immunity devoid of ILC2s. In terms of B lymphocyte lineages, caecal patch antigen-presenting memory B cells sat alongside germinal centre cells undergoing somatic hypermutation and differentiation into IGF1 + PCs. Prototypic gene expression signatures decreased across PC clusters, and notably, expanded IgA clonotypes could be traced in VDJ deep sequencing reads from additional compartments, including the bone marrow, supporting that these cells contribute a steady stream of systemic antibodies. The data advance our understanding of caecal immunological function, revealing processes involved in barrier maintenance and molecular networks relevant to disease.

Relapse after complete remission (CR) remains the main cause of mortality after allogeneic stem cell transplantation for hematological malignancies and, therefore, improved biomarkers for early prediction of relapse remains a critical goal toward development and assessment of preemptive relapse treatment. Because the significance of cancer stem cells as a source of relapses remains unclear, we investigated whether mutational screening for persistence of rare cancer stem cells would enhance measurable residual disease (MRD) and early relapse prediction after transplantation. In a retrospective study of patients who relapsed and patients who achieved continuous-CR with myelodysplastic syndromes and related myeloid malignancies, combined flow cytometric cell sorting and mutational screening for persistence of rare relapse-initiating stem cells was performed in the bone marrow at multiple CR time points after transplantation. In 25 CR samples from 15 patients that later relapsed, only 9 samples were MRD-positive in mononuclear cells (MNCs) whereas flowcytometric-sorted hematopoietic stem and progenitor cells (HSPCs) were MRD-positive in all samples, and always with a higher variant allele frequency than in MNCs (mean, 97-fold). MRD-positivity in HSPCs preceded MNCs in multiple sequential samples, in some cases preceding relapse by >2 years. In contrast, in 13 patients in long-term continuous-CR, HSPCs remained MRD-negative. Enhanced MRD sensitivity was also observed in total CD34+ cells, but HSPCs were always more clonally involved (mean, 8-fold). In conclusion, identification of relapse-initiating cancer stem cells and mutational MRD screening for their persistence consistently enhances MRD sensitivity and earlier prediction of relapse after allogeneic stem cell transplantation.

The precise pattern and timing of speciation events that gave rise to all living placental mammals remain controversial. We provide a comprehensive phylogenetic analysis of genetic variation across an alignment of 241 placental mammal genome assemblies, addressing prior concerns regarding limited genomic sampling across species. We compared neutral genome-wide phylogenomic signals using concatenation and coalescent-based approaches, interrogated phylogenetic variation across chromosomes, and analyzed extensive catalogs of structural variants. Interordinal relationships exhibit relatively low rates of phylogenomic conflict across diverse datasets and analytical methods. Conversely, X-chromosome versus autosome conflicts characterize multiple independent clades that radiated during the Cenozoic. Genomic time trees reveal an accumulation of cladogenic events before and immediately after the Cretaceous-Paleogene (K-Pg) boundary, implying important roles for Cretaceous continental vicariance and the K-Pg extinction in the placental radiation.