Scientific highlights

The European green toad (Bufotes viridis) is geographically widely distributed. While the species global conservation status is labeled as of least concern by the IUCN, it is declining in many parts of its range where populations are fragmented and isolated. A high-quality reference genome is an important resource for conservation genomic researchers who are trying to understand and interpret the genomic signals of population decline, inbreeding, and the accumulation of deleterious mutations. Here, we assembled and annotated a chromosome-level reference genome for B. viridis as part of the European Reference Genome Atlas pilot project. The genome assembly, with a size of ∼3.89 Gb consists of 11 chromosomes and an additional 2,096 unplaced scaffolds. The final assembly had a scaffold N50 value of 478.39 Mb and covered 90.4% single copy tetrapod orthologs, and 46.7% repetitive elements. Finally, a total of 23,830 protein-coding genes matching a known gene, together with 56,974 mRNAs were predicted. This high-quality reference genome will benefit amphibian evolutionary genomics research and enable conservation genetic studies to inform practical conservation work on this species.

Single-cell RNA sequencing (scRNA-seq) is a valuable tool for investigating cellular heterogeneity in diseases such as equine asthma (EA). This study evaluates the HIVE™ scRNA-seq method, a pico-well-based technology, for processing bronchoalveolar lavage (BAL) cells from horses with EA. The HIVE method offers practical advantages, including compatibility with both field and clinical settings, as well as a gentle workflow suited for handling sensitive cells. Our results show that the major cell types in equine BAL were successfully identified; however, the proportions of T cells and macrophages deviated from cytological expectations, with macrophages being overrepresented and T cells underrepresented. Despite these limitations, the HIVE method confirmed previously identified T cell and macrophage subpopulations and defined other BAL cell subsets. However, compared to previous studies T helper subsets were less clearly defined. Additionally, consistent with previous scRNA-seq studies, the HIVE method detected fewer granulocytes and mast cells than anticipated in the total BAL samples. Nevertheless, applying the method to purified mast cells recovered an expected number of cells. A small set of eosinophils were also detected which have not been characterized in earlier studies. In summary these findings suggest that while the HIVE method shows promise for certain applications, further optimization is needed to improve the accuracy of cell type representation, particularly for granulocytes and mast cells, in BAL samples.

The lack of population-scale databases hampers research and diagnostics for medically relevant tandem repeats and repeat expansions. We attempt to fill this gap using our pathSTR web tool, which leverages long-read sequencing of large cohorts to determine repeat length and sequence composition in a healthy population. The current version includes 1040 individuals of The 1000 Genomes Project cohort sequenced on the Oxford Nanopore Technologies PromethION. A comprehensive set of medically relevant tandem repeats has been genotyped using STRdust and LongTR to determine the tandem repeat length and sequence composition. PathSTR provides rich visualizations of this data set and the feature to upload one's data for comparison along the control cohort. We demonstrate the implementation of this application using data from targeted nanopore sequencing of a patient with myotonic dystrophy type 1. This resource will empower the genetics community to get a more complete overview of normal variation in tandem repeat length and sequence composition and, as such, enable a better assessment of rare tandem repeat alleles observed in patients.

Rare or de novo structural variation, primarily in the form of copy number variants, is detected in 5%-10% of autism spectrum disorder (ASD) families. While complex structural variants involving duplications can generally be detected using microarray or short-read genome sequencing (GS), these methods frequently fail to characterize breakpoints at nucleotide resolution, requiring additional molecular methods for validation and fine-mapping. Here, we use Oxford Nanopore Technologies PromethION long-read GS to characterize complex genomic rearrangements (CGRs) involving large duplications that segregate with ASD in five families. In total, we investigated 13 CGR carriers and were able to resolve all breakpoint junctions at nucleotide resolution. While all breakpoints were identified, the precise genomic architecture of one rearrangement remained unresolved with three different potential structures. The findings in two families include potential fusion genes formed through duplication rearrangements, involving IL1RAPL1-DMD and SUPT16H-CHD8 In two of the families originating from the same geographical region, an identical rearrangement involving ANK2 was identified, which likely represents a founder variant. In addition, we analyze methylation status directly from the long-read data, allowing us to assess the activity of rearranged genes and regulatory regions. Investigation of methylation across the CGRs reveals aberrant methylation status in carriers across a rearrangement affecting the CREBBP locus. In aggregate, our results demonstrate the utility of nanopore sequencing to pinpoint CGRs associated with ASD in five unrelated families, and highlight the importance of a gene-centric description of disease-associated complex chromosomal rearrangements.

Chromosomal inversions (INVs) are particularly challenging to detect due to their copy-number neutral state and association with repetitive regions. Inversions represent about 1/20 of all balanced structural chromosome aberrations and can lead to disease by gene disruption or altering regulatory regions of dosage-sensitive genes in cis Short-read genome sequencing (srGS) can only resolve ∼70% of cytogenetically visible inversions referred to clinical diagnostic laboratories, likely due to breakpoints in repetitive regions. Here, we study 12 inversions by long-read genome sequencing (lrGS) (n = 9) or srGS (n = 3) and resolve nine of them. In four cases, the inversion breakpoint region was missing from at least one of the human reference genomes (GRCh37, GRCh38, T2T-CHM13) and a reference agnostic analysis was needed. One of these cases, an INV9 mappable only in de novo assembled lrGS data using T2T-CHM13 disrupts EHMT1 consistent with a Mendelian diagnosis (Kleefstra syndrome 1; MIM#610253). Next, by pairwise comparison between T2T-CHM13, GRCh37, and GRCh38, as well as the chimpanzee and bonobo, we show that hundreds of megabases of sequence are missing from at least one human reference, highlighting that primate genomes contribute to genomic diversity. Aligning population genomic data to these regions indicated that these regions are variable between individuals. Our analysis emphasizes that T2T-CHM13 is necessary to maximize the value of lrGS for optimal inversion detection in clinical diagnostics. These results highlight the importance of leveraging diverse and comprehensive reference genomes to resolve unsolved molecular cases in rare diseases.

Clinical genetic laboratories often require a comprehensive analysis of chromosomal rearrangements/structural variants (SVs), from large events like translocations and inversions to supernumerary ring/marker chromosomes and small deletions or duplications. Understanding the complexity of these events and their clinical consequences requires pinpointing breakpoint junctions and resolving the derivative chromosome structure. This task often surpasses the capabilities of short-read sequencing technologies. In contrast, long-read sequencing techniques present a compelling alternative for clinical diagnostics. Here, Genomic Medicine Sweden-Rare Diseases has explored the utility of HiFi Revio long-read genome sequencing (lrGS) for digital karyotyping of SVs nationwide. The 16 samples from 13 families were collected from all Swedish healthcare regions. Prior investigations had identified 16 SVs, ranging from simple to complex rearrangements, including inversions, translocations, and copy number variants. We have established a national pipeline and a shared variant database for variant calling and filtering. Using lrGS, 14 of the 16 known SVs are detected. Of these, 13 are mapped at nucleotide resolution, and one complex rearrangement is only visible by read depth. Two Chromosome 21 rearrangements, one mosaic, remain undetected. Average read lengths are 8.3-18.8 kb with coverage exceeding 20× for all samples. De novo assembly results in a limited number of phased contigs per individual (N50 6-86 Mb), enabling direct characterization of the chromosomal rearrangements. In a national pilot study, we demonstrate the utility of HiFi Revio lrGS for analyzing chromosomal rearrangements. Based on our results, we propose a 5-year plan to expand lrGS use for rare disease diagnostics in Sweden.

The early evolution of sex chromosomes has remained obscure for more than a century. The Vandiemenella viatica species group of morabine grasshoppers is highly suited for studying the early stages of sex chromosome divergence and degeneration of the Y chromosome. This stems from the fact that neo-XY sex chromosomes have independently evolved multiple times by X-autosome fusions with different autosomes. Here, we generated new chromosome-level assemblies for two chromosomal races representing karyotypes with and without neo-sex chromosomes (P24XY and P24X0), and sequence data of a third chromosomal race with a different neo-XY chromosome system (P25XY). Interestingly, these two neo-XY chromosomal races are formed by different X-autosome fusions (involving chr1 and chrB, respectively), and we found that both neo-Y chromosomes have partly ceased to recombine with their neo-X counterpart. We show that the neo-XY chromosomes have diverged through accumulation of SNPs and structural mutations, and that many neo-Y-linked genes have degenerated since recombination ceased. However, the non-recombining regions of neo-Y chromosomes host non-degenerated genes crucial for sex determination, such as sex-lethal and transformer, alongside genes associated with spermatogenesis, fertility, and reproduction, illustrating their integrative role as a masculinizing supergene. Contrary to expectations, the neo-Y chromosomes showed (slightly) lower density of transposable elements (TEs) compared to other genomic regions. The study reveals the unique dynamics of young sex chromosomes, with evolution of recombination suppression and pronounced decay of (some) neo-sex chromosome genes, and provides a compelling case illustrating how chromosomal fusions and post-fusion mutational processes contribute to the evolution of supergenes.

We present a genome assembly from an individual female Haliaeetus albicilla (the white-tailed eagle; Chordata; Aves; Accipitriformes; Accipitridae). The genome sequence has a total length of 1,320.30 megabases. Most of the assembly is scaffolded into 34 chromosomal pseudomolecules, including the Z and W sex chromosomes. Gene annotation of this assembly on Ensembl identified 17,501 protein-coding genes.

Morabine grasshoppers in the Vandiemenella viatica species group, which show karyotype diversity, have been studied for their ecological distribution and speciation in relation to their genetic and chromosomal diversity. They are good models for studying sex chromosome evolution as "old" and newly emerged sex chromosomes co-exist within the group. Here we present a reference genome for the viatica19 chromosomal race, that possesses the ancestral karyotype within the group. Using PacBio HiFi and Hi-C sequencing, we generated a chromosome-level assembly of 4.09 Gb in span, scaffold N50 of 429 Mb, and complete BUSCO score of 98.1%, containing 10 pseudo-chromosomes. We provide Illumina datasets of males and females, used to identify the X chromosome. The assembly contains 19,034 predicted protein-coding genes, and a total of 75.21% of repetitive DNA sequences. By leveraging HiFi reads, we mapped the genome-wide distribution of methylated bases (5mC and 6 mA). This comprehensive assembly offers a robust reference for morabine grasshoppers and supports further research into speciation and sex chromosome diversification within the group and its related species.

Cancer has been characterized as a wound that does not heal. Malignant cells are morphologically distinct from normal proliferating cells but have extensive similarities to tissues undergoing wound healing and/or regeneration. The mechanistic basis of this similarity has, however, remained enigmatic. Here, we show that the genomic region upstream of Myc, which carries more cancer susceptibility in humans than any other genomic region, is required for intestinal regeneration after radiation damage. Failure to regenerate is associated with inefficient Ly6a/Sca1+ stem/progenitor cell mobilization, and almost complete failure to re-establish Lgr5+ cell compartment in the intestinal crypts. The Myc upstream region is also critical for growth of adult intestinal cells in 3D organoid culture. We show that culture conditions recapitulating most aspects of adult normal tissue architecture still reprogram normal cells to proliferate using a mechanism similar to that employed by cancer cells. Our results establish a function for the Myc super-enhancer region as the genetic link between tissue regeneration and tumorigenesis, and demonstrates that normal tissue renewal and regeneration of tissues after severe damage are mechanistically distinct.2-540

In the Boreal region, extreme seasonal variations in day-night length expose communities to dynamic light and temperature fluctuations. Freshwater bacterioplankton, representing key ecosystem components, faces climate-driven shifts; yet the fixed day-length patterns determined by latitude underscore the importance of studying light's role in predicting ecosystem responses. We investigated bacterial community composition in a brown peat bog and a clear oligotrophic lake across seasons with contrasting light regimes: the summer solstice (> 20 h of daylight) and the autumn equinox (equal day-night length). Using amplicon sequencing of 16S rRNA transcripts, alongside measurements of physicochemical parameters, organic matter characterisation and dissolved carbon dioxide and methane gas measurements, we found no diel cycling in the lake during either period or in the peat bog near the summer solstice. However, the structure of bacterial peat bog communities exhibited cyclic changes over diel cycles at the autumn equinox. Twelve amplicon sequence variants, including both phototrophic and heterotrophic taxa, increased in abundance at all measured morning sampling times. These findings provide valuable insights into the diel patterns of boreal lentic habitats and their bacterioplankton communities, highlighting the absence of diel fluctuations in some systems and seasons, while revealing cyclic dynamics in others, driven by conditionally rare taxa.

Current intensive treatment of pediatric T-cell acute lymphoblastic leukemia (T-ALL) has substantial side effects, highlighting a need for novel biomarkers to improve risk stratification. Canonical biomarkers, such as genetics and immunophenotype, are largely not used in pediatric T-ALL stratification. This study aimed to validate the prognostic relevance of DNA methylation CpG island methylator phenotype (CIMP) risk stratification in 2 pediatric T-ALL patient cohorts: the Nordic Society of Paediatric Haematology (NOPHO) ALL2008 T-ALL study cohort (n = 192) and the Dutch Childhood Oncology Group (DCOG) ALL-10/ALL-11 validation cohorts (n = 156). Both cohorts revealed that combining CIMP classification at diagnosis with measurable residual disease (MRD) at treatment day 29 (D29) or 33 (D33) significantly improved outcome prediction. The poor prognosis subgroup, characterized by CIMP low/D29 or D33 MRD ≥ 0.1%, had a cumulative incidence of relapse (pCIR5yr) of 29% and 23% and overall survival (pOS5yr) of 59.7% and 65.4%, in NOPHO and DCOG, respectively. Conversely, a good prognosis subgroup was also identified representing CIMP high/D29 or D33 MRD < 0.1% with pCIR5yr of 0% and 3.4% and pOS5yr of 98.2% and 94.8%, in NOPHO and DCOG, respectively. For NOPHO, MRD was also evaluated on D15, and the relapse prediction accuracy of CIMP/D29 MRD (0.79) and CIMP/D15 MRD (0.75) classification was comparable, indicating potential for earlier stratification. The evaluation of the biology behind the CIMP subgroups revealed associations with transcriptome profiles, genomic aberrations, and mitotic history, suggesting distinct routes for leukemia development. In conclusion, integrating MRD assessment with the novel CIMP biomarker has the potential to improve risk stratification in pediatric T-ALL and guide future therapeutic decisions.

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.

The Oriental fruit moth Grapholita molesta is distributed throughout temperate regions and considered to be a pest in peach production and other high-value fruit crops in the rose family. Insecticide treatment has led to resistance development, but the use of sex pheromones in pest management has shown great promise. We investigated the pheromone biosynthesis pathway in G. molesta with the aim of elucidating pheromone evolution in the Olethreutinae subfamily of moths and harnessing pathway genes in biotechnological production of sex pheromone for use in pest management. In vivo labelling experiments suggested that an uncommon Δ8 fatty acyl desaturase is involved in sex pheromone biosynthesis. CRISPR/Cas9 knock-out of the highly expressed candidate desaturase gene Gmol_CPRQ almost completely blocked the production of Δ8 pheromone components in vivo. Heterologous expression of Gmol_CPRQ protein in yeast- or Sf9 insect cells, however, failed to demonstrate the expected Δ8 desaturase activity. Instead, Δ9 desaturase activity was observed. Co-expression in the yeast system of the electron donor, cytochrome b5, from G. molesta still produced only Δ9 desaturase activity. We suggest that Gmol_CPRQ is intimately involved in pheromone production in vivo, via an unknown reaction mechanism that may possibly involve another co-factor that is absent in the yeast and Sf9 expression systems, or depend on its subcellular site of activity. Solving this puzzle will shed further light on pheromone biosynthesis in the family Tortricidae and will be required for successful biotechnological production of fatty acids and pheromones requiring Δ8 desaturation.

Members of the tribe Palicoureeae of the coffee family (Rubiaceae) have a complex taxonomic history and have been the focus of few modern systematic studies. The tribe comprises about 1,100 tropical species in ten genera. To investigate phylogeny, we used a target capture approach and the angiosperm-wide Angiosperms353 bait set to produce genomic data for a representative taxon sample of Palicoureeae, with particular focus on the African genus Hymenocoleus. Using coalescent-based inference methods, we find that Puffia gerrardii (recently separated from Geophila) is sister to Hymenocoleus. The deepest split in Hymenocoleus is highly affected by incomplete lineage sorting, possibly as a consequence of rapid speciation during the early evolution of the clade. Remaining interspecific relationships in Hymenocoleus could be confidently resolved and while Robbrecht's traditional subgeneric classification scheme based on floral features is not supported as reflecting evolution in the group, we find that several other features do, e.g. characters of pyrenes and involucral cups. Although not free of challenges, a strong advantage with our analytical approach is that gene tree heterogeneity can be taken into account. Including flanking regions yielded data sets that had the strongest power to reject polytomies and produced less gene tree error, resulting in species trees with higher normalised quartet scores and higher average support compared to trees inferred only from exon data. Presumably paralogous loci are often filtered out prior to species tree estimation but we find that they may contribute important phylogenetic information when using an inference method that actively accounts for them.

Acute pain is an unpleasant experience caused by noxious stimuli. How the spinal neural circuits attribute differences in quality of noxious information remains unknown. By means of genetic capturing, activity manipulation and single-cell RNA sequencing, we identified distinct neural ensembles in the adult mouse spinal cord encoding mechanical and heat pain. Reactivation or silencing of these ensembles potentiated or stopped, respectively, paw shaking, lifting and licking within but not across the stimuli modalities. Within ensembles, polymodal Gal+ inhibitory neurons with monosynaptic contacts to A-fiber sensory neurons gated pain transmission independent of modality. Peripheral nerve injury led to inferred microglia-driven inflammation and an ensemble transition with decreased recruitment of Gal+ inhibitory neurons and increased excitatory drive. Forced activation of Gal+ neurons reversed hypersensitivity associated with neuropathy. Our results reveal the existence of a spinal representation that forms the neural basis of the discriminative and defensive qualities of acute pain, and these neurons are under the control of a shared feed-forward inhibition.

Chronic inflammation contributes to the development of colorectal cancer, partly through its regulation of the microenvironment and antitumor immunity. Interestingly, women have a lower incidence of colorectal cancer, and estrogen treatment has been shown to reduce the occurrence of colorectal tumors. While intestinal estrogen receptor beta (ERβ, Esr2) can protect against colitis and colitis-induced cancer in mice, its role in shaping the tumor microenvironment remains unknown. In this study, we performed RNA sequencing to analyze the transcriptome of colonic epithelia and tumors from azoxymethane/dextran sulfate sodium-treated wild-type and intestinal ERβ knockout (ERβKOVil) mice and vehicle-treated controls. This revealed significant differences in gene expression and enriched biological processes influenced by sex and genotype, with immune-related responses being overrepresented. Deconvolution supported differential immune cell abundance and immunostaining showed that tumors from ERβKOVil mice displayed significantly increased macrophage infiltration, decreased T cell infiltration, and impaired natural killer cell infiltration. Further, ERβ mRNA levels in clinical colorectal tumors correlated with immune signaling profiles and better survival. Our findings indicate that intestinal ERβ promotes an antitumor microenvironment and could potentially affect the effectiveness of immunotherapy. These insights highlight the importance of ERβ in modulating antitumor immunity and underscore its therapeutic potential in colorectal cancer.

Correct timing of neurogenesis is crucial for generating the correct number and subtypes of glia and neurons in the embryo, and for preventing tumours and stem cell depletion in the adults. Here, we analyse how the midbrain dopamine (mDA) neuron progenitors transition into cell cycle arrest (G0) and begin to mature into ependymal cells. Comparison of mDA progenitors from different embryonic stages revealed upregulation of the genes encoding Sox9 and nuclear factor I transcription factors during development. Their conditional inactivation in the early embryonic midbrain led to delayed G0 entry and ependymal maturation in the entire midbrain ventricular zone, reduced gliogenesis and increased generation of neurons, including mDA neurons. In contrast, their inactivation in late embryogenesis did not result in mitotic re-entry, suggesting that these factors are necessary for G0 induction, but not for its maintenance. Our characterisation of adult ependymal cells by single-cell RNA sequencing and histology show that mDA-progenitor-derived cells retain several progenitor features but also secrete neuropeptides and contact neighbouring cells and blood vessels, indicating that these cells may form part of the circumventricular organ system.

Uncovering early gene network changes of human hematopoietic stem cells (HSCs) leading to differentiation induction is of utmost importance for therapeutic manipulation. We employed single cell proteo-transcriptomic sequencing to FACS-enriched bone marrow hematopoietic stem and progenitor cells (HSPCs) from 15 healthy donors. Pseudotime analysis reveals four major differentiation trajectories, which remain consistent upon aging, with an early branching point into megakaryocyte-erythroid progenitors. However, young donors suggest a more productive differentiation from HSPCs to committed progenitors of all lineages. tradeSeq analysis depicts continuous changes in gene expression of HSPC-related genes (DLK1, ADGRG6), and provides a roadmap of gene expression at the earliest branching points. We identify CD273/PD-L2 to be highly expressed in a subfraction of immature multipotent HSPCs with enhanced quiescence. Functional experiments confirm the immune-modulatory function of CD273/PD-L2 on HSPCs in regulating T-cell activation and cytokine release. Here, we present a molecular map of early HSPC differentiation across human life.