Scientific highlights

The Norwegian lemming (Lemmus lemmus) is a small rodent distributed across the Fennoscandian mountain tundra and the Kola Peninsula. The Norwegian lemming likely evolved during the Late Pleistocene and inhabited Fennoscandia shortly prior to the Last Glacial Maximum. However, the exact timing and origins of the species, and its phylogenetic position relative to the closely related Siberian lemming (Lemmus sibiricus) remain disputed. Moreover, the presence of ancient or contemporary gene flow between both species is largely untested. The Norwegian lemming displays characteristic phenotypic and behavioral adaptations (e.g., coat color, aggression) that are not present in other Lemmus species. We generated a de novo genome assembly for the Norwegian lemming and resequenced nine modern and two ancient Lemmus spp. genomes. We show that all Lemmus species form distinct monophyletic clades, with concordant topology between the mitochondrial and nuclear genome phylogenies. The Siberian lemming is divided into two distinct but paraphyletic clades, one in the east and one in the west, where the western clade represents a sister taxon to the Norwegian lemming. We estimate that the Norwegian and western Siberian lemming diverged shortly before the Last Glacial Maximum, making the Norwegian lemming one of the youngest known mammalian species. We did not find any indication of gene flow between L. lemmus and L. sibiricus, suggesting postglacial isolation of L. lemmus. Furthermore, we identify species-specific genomic differences in genes related to coat color and fat transport, which are likely associated with the distinctive coloration and overwintering behavior observed in the Norwegian lemming.

Single-cell transcriptomics is a key tool for unravelling metabolism and tissue diversity in model organisms. Its potential for elucidating the ecological roles of microeukaryotes, especially non-model ones, remains largely unexplored. This study employed the Smart-seq2 protocol on Ochromonas triangulata, a microeukaryote lacking a reference genome, showcasing how transcriptional states align with two distinct growth phases: a fast-growing phase and a slow-growing phase. Besides the two expected expression clusters, each corresponding to either growth phase, a third transcriptional state was identified across both growth phases. Metabolic mapping revealed a boost of photosynthetic activity in the fast growth over the slow growth stage, as well as downregulation trend in pathways associated with ribosome functioning, CO2 fixation, and carbohydrate catabolism characteristic of the third transcriptional state. In addition, carry-over rRNA reads recapitulated the taxonomic identity of the target while revealing distinct bacterial communities, in co-culture with the eukaryote, each associated with distinct transcriptional states. This study underscores single-cell transcriptomics as a powerful tool for characterizing metabolic states in microeukaryotes without a reference genome, offering insights into unknown physiological states and individual-level interactions with different bacterial taxa. This approach holds broad applicability to describe the ecological roles of environmental microeukaryotes, culture-free, and reference-free, surpassing alternative methods like metagenomics or metatranscriptomics.

How neurons can maintain cellular identity over an entire life span remains largely unknown. Here, we show that maintenance of identity in differentiated dopaminergic and serotonergic neurons is critically reliant on the Polycomb repressive complex 2 (PRC2). Deletion of the obligate PRC2 component, Eed, in these neurons resulted in global loss of H3K27me3, followed by a gradual activation of genes harboring both H3K27me3 and H3K9me3 modifications. Notably, H3K9me3 was lost at these PRC2 targets before gene activation. Neuronal survival was not compromised; instead, there was a reduction in subtype-specific gene expression and a progressive impairment of dopaminergic and serotonergic neuronal function, leading to behavioral deficits characteristic of Parkinson's disease and anxiety. Single-cell analysis revealed subtype-specific vulnerability to loss of PRC2 repression in dopamine neurons of the substantia nigra. Our study reveals that a PRC2-dependent nonpermissive chromatin state is essential to maintain the subtype identity and function of dopaminergic and serotonergic neurons.

Immune-microbe interactions early in life influence the risk of allergies, asthma, and other inflammatory diseases. Breastfeeding guides healthier immune-microbe relationships by providing nutrients to specialized microbes that in turn benefit the host's immune system. Such bacteria have co-evolved with humans but are now increasingly rare in modern societies. Here we show that a lack of bifidobacteria, and in particular depletion of genes required for human milk oligosaccharide (HMO) utilization from the metagenome, is associated with systemic inflammation and immune dysregulation early in life. In breastfed infants given Bifidobacterium infantis EVC001, which expresses all HMO-utilization genes, intestinal T helper 2 (Th2) and Th17 cytokines were silenced and interferon β (IFNβ) was induced. Fecal water from EVC001-supplemented infants contains abundant indolelactate and B. infantis-derived indole-3-lactic acid (ILA) upregulated immunoregulatory galectin-1 in Th2 and Th17 cells during polarization, providing a functional link between beneficial microbes and immunoregulation during the first months of life.

The introduction of the CTLA-4 recombinant fusion protein has demonstrated therapeutic effects by selectively modulating T-cell activation in rheumatoid arthritis. Here we show, using a forward genetic approach, that a mutation in the SH3gl1 gene encoding the endocytic protein Endophilin A2 is associated with the development of arthritis in rodents. Defective expression of SH3gl1 affects T cell effector functions and alters the activation threshold of autoreactive T cells, thereby leading to complete protection from chronic autoimmune inflammatory disease in both mice and rats. We further show that SH3GL1 regulates human T cell signaling and T cell receptor internalization, and its expression is upregulated in rheumatoid arthritis patients. Collectively our data identify SH3GL1 as a key regulator of T cell activation, and as a potential target for treatment of autoimmune diseases.

How time is measured by neural stem cells during temporal neurogenesis has remained unresolved. By combining experiments and computational modeling, we define a Shh/Gli-driven three-node timer underlying the sequential generation of motor neurons (MNs) and serotonergic neurons in the brainstem. The timer is founded on temporal decline of Gli-activator and Gli-repressor activities established through down-regulation of Gli transcription. The circuitry conforms an incoherent feed-forward loop, whereby Gli proteins not only promote expression of Phox2b and thereby MN-fate but also account for a delayed activation of a self-promoting transforming growth factor-β (Tgfβ) node triggering a fate switch by repressing Phox2b. Hysteresis and spatial averaging by diffusion of Tgfβ counteract noise and increase temporal accuracy at the population level, providing a functional rationale for the intrinsically programmed activation of extrinsic switch signals in temporal patterning. Our study defines how time is reliably encoded during the sequential specification of neurons.

As the only endemic neotropical parrot to have recently lived in the northern hemisphere, the Carolina parakeet (Conuropsis carolinensis) was an iconic North American bird. The last surviving specimen died in the Cincinnati Zoo in 1918 [1]. The cause of its extinction remains contentious: besides excessive mortality associated to habitat destruction and active hunting, their survival could have been negatively affected by its range having become increasingly patchy [2] or by the exposure to poultry pathogens [3, 4]. In addition, the Carolina parakeet showed a predilection for cockleburs, an herbaceous plant that contains a powerful toxin, carboxyatractyloside, or CAT [5], which did not seem to affect them but made the birds notoriously toxic to most predators [3]. To explore the demographic history of this bird, we generated the complete genomic sequence of a preserved specimen held in a private collection in Espinelves (Girona, Spain), as well as of a close extant relative, Aratinga solstitialis. We identified two non-synonymous genetic changes in two highly conserved proteins known to interact with CAT that could underlie a specific dietary adaptation to this toxin. Our genomic analyses did not reveal evidence of a dramatic past demographic decline in the Carolina parakeet; also, its genome did not exhibit the long runs of homozygosity that are signals of recent inbreeding and are typically found in endangered species. As such, our results suggest its extinction was an abrupt process and thus likely solely attributable to human causes.

The brain is composed of hundreds of different neuronal subtypes, which largely retain their identity throughout the lifespan of the organism. The mechanisms governing this stability are not fully understood, partly due to the diversity and limited size of clinically relevant neuronal populations, which constitute a technical challenge for analysis. Here, using a strategy that allows for ChIP-seq combined with RNA-seq in small neuronal populations in vivo, we present a comparative analysis of permissive and repressive histone modifications in adult midbrain dopaminergic neurons, raphe nuclei serotonergic neurons, and embryonic neural progenitors. Furthermore, we utilize the map generated by our analysis to show that the transcriptional response of midbrain dopaminergic neurons following 6-OHDA or methamphetamine injection is characterized by increased expression of genes with promoters dually marked by H3K4me3/H3K27me3. Our study provides an in vivo genome-wide analysis of permissive/repressive histone modifications coupled to gene expression in these rare neuronal subtypes.

Here we describe the SweGen data set, a comprehensive map of genetic variation in the Swedish population. These data represent a basic resource for clinical genetics laboratories as well as for sequencing-based association studies by providing information on genetic variant frequencies in a cohort that is well matched to national patient cohorts. To select samples for this study, we first examined the genetic structure of the Swedish population using high-density SNP-array data from a nation-wide cohort of over 10 000 Swedish-born individuals included in the Swedish Twin Registry. A total of 1000 individuals, reflecting a cross-section of the population and capturing the main genetic structure, were selected for whole-genome sequencing. Analysis pipelines were developed for automated alignment, variant calling and quality control of the sequencing data. This resulted in a genome-wide collection of aggregated variant frequencies in the Swedish population that we have made available to the scientific community through the website https://swefreq.nbis.se. A total of 29.2 million single-nucleotide variants and 3.8 million indels were detected in the 1000 samples, with 9.9 million of these variants not present in current databases. Each sample contributed with an average of 7199 individual-specific variants. In addition, an average of 8645 larger structural variants (SVs) were detected per individual, and we demonstrate that the population frequencies of these SVs can be used for efficient filtering analyses. Finally, our results show that the genetic diversity within Sweden is substantial compared with the diversity among continental European populations, underscoring the relevance of establishing a local reference data set.

Resolving the molecular details of proteome variation in the different tissues and organs of the human body will greatly increase our knowledge of human biology and disease. Here, we present a map of the human tissue proteome based on an integrated omics approach that involves quantitative transcriptomics at the tissue and organ level, combined with tissue microarray-based immunohistochemistry, to achieve spatial localization of proteins down to the single-cell level. Our tissue-based analysis detected more than 90% of the putative protein-coding genes. We used this approach to explore the human secretome, the membrane proteome, the druggable proteome, the cancer proteome, and the metabolic functions in 32 different tissues and organs. All the data are integrated in an interactive Web-based database that allows exploration of individual proteins, as well as navigation of global expression patterns, in all major tissues and organs in the human body.

Around 30% of patients with rheumatoid arthritis (RA) lack rheumatoid factor and anti-citrullinated protein antibodies (ACPA) complicating diagnosis and potentially delaying treatment. We hypothesised that innate immune mechanisms might be more prominent in ACPA- RA. We performed single-cell RNA sequencing of mononuclear cells from peripheral blood (PBMC) and synovial fluid (SFMC) of patients with ACPA- and ACPA+ RA (n = 4 per group: discovery cohort; n = 8 per group: validation cohort). Dendritic cells and proinflammatory cytokine production were analysed by flow cytometry on SFMC from patients with ACPA- RA, ACPA+ RA, and psoriatic arthritis. Interferon (IFN) levels in synovial fluid (SF) and serum were measured in these groups. Several macrophage subsets and cDC2 were enriched in ACPA- RA SF whereas the frequency of Tph and B cells was increased in ACPA+ RA SF. Type I IFN-stimulated genes were detected in SFMC, but not PBMC, of patients with ACPA- RA. A type I IFN signature was also observed in synovial tissue from two patients with ACPA- RA in an independent dataset. IFN levels were higher in SF than serum but IFN-α/β production did not differ between ACPA+ and ACPA- RA. This study identifies a distinct innate cell composition and type I IFN gene response in synovial joints, but not in peripheral blood, of patients with ACPA- RA. Similar IFN levels across groups suggest the IFN signature may have been primed before the cells entered the joints. These findings provide a foundation for future research on type I IFN responses in ACPA- RA.

Racing without protective shoes is common in the Swedish harness racing industry, as it can enhance horses' performance on the track. Trainers typically decide whether a horse will race barefoot based on practical experience rather than objective measures. However, this practice can sometimes lead to excessive hoof wear, posing potential welfare concerns for racing horses. Gene expression differences may help reveal the underlying genetic mechanisms associated with different phenotypic traits. To explore an objective measure for assessing which horses are best suited for barefoot racing, we conducted a polyA-selected RNA-seq experiment on tissue from the growth zone at the coronary band of the hoof. This experiment compared tissues from Standardbred trotters capable of repeatedly racing barefoot without injury (n = 11) to those that could not (n = 7). By combining stringent phenotyping with racing records and trainer interviews, we aimed to elucidate the biological factors related to hoof strength in barefoot racing, focusing on differential abundant genes. The RNA-seq analysis identified five significantly downregulated genes in horses capable of competing barefoot across consecutive races. These genes are associated with various biological processes relevant for hoof strength: ACCS, IRX2 and TRAPPAC6A contribute to enhancing the structural integrity of the hoof; MT2A regulates its metal homeostasis and SLC35F3 likely influences local vasoconstriction in the hoof. These gene findings suggest a coordinated genetic basis for structural reinforcement and physiological support of the hoof, which may be critical for sustaining performance under barefoot conditions. Our findings suggest that the ability of Standardbred trotters to race barefoot in consecutive events is reflected in distinct gene expression patterns, underscoring a genetic basis for hoof strength. This supports further genome-wide scans aimed at identifying genetic markers for hoof durability in these horses. The focused design of our study- comparing horses that could consistently race barefoot with those that could not- enabled us to isolate a select group of genes involved in diverse aspects of hoof biology essential for quality and resilience of horse hooves. This insight could ultimately be applied to augment both the performance and wellbeing of equine athletes across disciplines.

Imaging-defined atherosclerosis represents an intermediate phenotype of atherosclerotic cardiovascular disease (ASCVD). Genome-wide association studies (GWAS) on directly measured coronary plaques using coronary computed tomography angiography (CCTA) are scarce. In the so far largest population-based cohort with CCTA data, we performed a GWAS on coronary plaque burden as determined by the segment involvement score (SIS) in 24,811 European individuals. We identified 20 significant independent genetic markers for SIS, three of which were found in loci not implicated in ASCVD before. Further GWAS on coronary artery calcification showed similar results to that of SIS, whereas a GWAS on ultrasound-assessed carotid plaques identified both shared and non-shared loci with SIS. In two-sample Mendelian randomization studies using SIS-associated markers in UK Biobank and CARDIoGRAMplusC4D, one extra coronary segment with atherosclerosis corresponded to 1.8-fold increased odds of myocardial infarction. This GWAS data can aid future studies of causal pathways in ASCVD.

Cerebral Cavernous Malformation (CCM) is a brain vascular disease with various neurological symptoms. In this study, we describe the inflammatory profile in CCM and show for the first time the formation of neutrophil extracellular traps (NETs) in rodents and humans with CCM. Through RNA-seq analysis of cerebellum endothelial cells from wild-type mice and mice with an endothelial cell-specific ablation of the Ccm3 gene (Ccm3iECKO), we show that endothelial cells from Ccm3iECKO mice have an increased expression of inflammation-related genes. These genes encode proinflammatory cytokines and chemokines, as well as adhesion molecules, which promote recruitment of inflammatory and immune cells. Similarly, immunoassays showed elevated levels of these cytokines and chemokines in the cerebellum of the Ccm3iECKO mice. Consistently, both flow cytometry and immunofluorescence analysis showed infiltration of different subsets of leukocytes into the CCM lesions. Neutrophils, which are known to fight against infection through different strategies, including the formation of NETs, represented the leukocyte subset within the most pronounced increase in CCM. Here, we detected elevated levels of NETs in the blood and the deposition of NETs in the cerebral cavernomas of Ccm3iECKO mice. Degradation of NETs by DNase I treatment improved the vascular barrier. The deposition of NETs in the cavernomas of patients with CCM confirms the clinical relevance of NETs in CCM.

The genus Rhodotorula includes basidiomycetous oleaginous yeast species. Rhodotorula babjevae can produce compounds of biotechnological interest such as lipids, carotenoids, and biosurfactants from low value substrates such as lignocellulose hydrolysate. High-quality genome assemblies are needed to develop genetic tools and to understand fungal evolution and genetics. Here, we combined short- and long-read sequencing to resolve the genomes of two R. babjevae strains, CBS 7808 (type strain) and DBVPG 8058, at chromosomal level. Both genomes are 21 Mbp in size and have a GC content of 68.2%. Allele frequency analysis indicates that both strains are tetraploid. The genomes consist of a maximum of 21 chromosomes with a size of 0.4 to 2.4 Mbp. In both assemblies, the mitochondrial genome was recovered in a single contig, that shared 97% pairwise identity. Pairwise identity between most chromosomes ranges from 82 to 87%. We also found indications for strain-specific extrachromosomal endogenous DNA. A total of 7591 and 7481 protein-coding genes were annotated in CBS 7808 and DBVPG 8058, respectively. CBS 7808 accumulated a higher number of tandem duplications than DBVPG 8058. We identified large translocation events between putative chromosomes. Genome divergence values between the two strains indicate that they may belong to different species.

SHORT VEGETATIVE PHASE (SVP) is an important regulator of FLOWERING LOCUS T (FT) in the thermosensory pathway of Arabidopsis. It is a negative regulator of flowering and represses FT transcription. In poplar trees, FT2 is central for the photoperiodic control of growth cessation, which also requires the decrease of bioactive gibberellins (GAs). In angiosperm trees, genes similar to SVP, sometimes named DORMANCY-ASSOCIATED MADS-BOX genes, control temperature-mediated bud dormancy. Here we show that SVL, an SVP ortholog in aspen trees, besides its role in controlling dormancy through its expression in buds, is also contributing to the regulation of short day induced growth cessation and bud set through its expression in leaves. SVL is upregulated during short days in leaves and binds to the FT2 promoter to repress its transcription. It furthermore decreases the amount of active GAs, whose downregulation is essential for growth cessation, by repressing the transcription of GA20 oxidase. Finally, the SVL protein is more stable in colder temperatures, thus integrating the temperature signal into the response. We conclude that the molecular function of SVL in the photoperiodic pathway has been conserved between Arabidopsis and poplar trees, albeit the physiological process it controls has changed. SVL is thus both involved in regulating the photoperiod response in leaves, modulating the timing of growth cessation and bud set, and in the subsequent temperature regulation of dormancy in the buds.

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

The HERMES (HEart failure Molecular Epidemiology for Therapeutic targetS) consortium aims to identify the genomic and molecular basis of heart failure. The consortium currently includes 51 studies from 11 countries, including 68 157 heart failure cases and 949 888 controls, with data on heart failure events and prognosis. All studies collected biological samples and performed genome-wide genotyping of common genetic variants. The enrolment of subjects into participating studies ranged from 1948 to the present day, and the median follow-up following heart failure diagnosis ranged from 2 to 116 months. Forty-nine of 51 individual studies enrolled participants of both sexes; in these studies, participants with heart failure were predominantly male (34-90%). The mean age at diagnosis or ascertainment across all studies ranged from 54 to 84 years. Based on the aggregate sample, we estimated 80% power to genetic variant associations with risk of heart failure with an odds ratio of ≥1.10 for common variants (allele frequency ≥ 0.05) and ≥1.20 for low-frequency variants (allele frequency 0.01-0.05) at P < 5 × 10-8 under an additive genetic model. HERMES is a global collaboration aiming to (i) identify the genetic determinants of heart failure; (ii) generate insights into the causal pathways leading to heart failure and enable genetic approaches to target prioritization; and (iii) develop genomic tools for disease stratification and risk prediction.