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.

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.

Single-molecule techniques are ideally poised to characterize complex dynamics but are typically limited to investigating a small number of different samples. However, a large sequence or chemical space often needs to be explored to derive a comprehensive understanding of complex biological processes. Here we describe multiplexed single-molecule characterization at the library scale (MUSCLE), a method that combines single-molecule fluorescence microscopy with next-generation sequencing to enable highly multiplexed observations of complex dynamics. We comprehensively profiled the sequence dependence of DNA hairpin properties and Cas9-induced target DNA unwinding-rewinding dynamics. The ability to explore a large sequence space for Cas9 allowed us to identify a number of target sequences with unexpected behaviors. We envision that MUSCLE will enable the mechanistic exploration of many fundamental biological processes.

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

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.

In Sweden, reforestation of managed forests relies predominantly on planting nursery-produced tree seedlings. However, the intense production using containerized cultivation systems (e.g., high seedling density, irrigation from above, regular fertilization) creates favorable conditions for fungal infections. Despite the harmful role of diseases in forest nurseries, the origin and dispersal factors of fungal pathogens remain largely unknown. A better understanding of the airborne spread of pathogens could improve the prediction of fungal infection, ultimately optimizing preventative methods and decreasing the use of fungicides. This study investigated the temporal dynamics of airborne fungi in forest nurseries, with a focus on fungal pathogens. Airborne fungi were monitored in four Swedish forest nurseries over two growing seasons using spore traps and high-throughput sequencing. Fungal pathogens were identified using bioinformatics and quantified with quantitative PCR. Results showed strong temporal shifts of airborne fungal diversity and community composition following the growing seasons. The airborne spread included high abundances of important fungal pathogens (e.g., Cladosporium sp., Botrytis cinerea, Alternaria sp., Sydowia polyspora, and Melampsora populnea) with individual temporal and spatial variations. In general, the deposited spore loads of nursery pathogens correlated positively with increased temperature and negatively with higher precipitation. This was expressed the strongest for Cladosporium sp., Alternaria sp., and M. populnea, which suggests a higher availability of fungal inoculum in warm and dry periods. This study highlights the influence of seasonality on the temporal dynamics of economically important fungal pathogens in Swedish forest nurseries, which should be considered in the development of a local decision support system.IMPORTANCEFungal diseases in forest nurseries have significant environmental and economic impacts on the tree seedling production. This study highlights the role of seasonality in the airborne spread of fungal pathogens in Swedish forest nurseries. By analyzing airborne fungal spores using advanced sequencing and monitoring techniques, key fungal pathogens and their dispersal patterns over two growing seasons were identified. The findings indicate that warmer, drier periods may increase the spread of fungal pathogens, emphasizing the need for targeted preventative measures. Understanding these temporal dynamics can help optimize the use of fungicides in forest nurseries, thereby promoting more sustainable and environmentally friendly management practices. This research provides valuable insights for improving disease management in forest nurseries, ultimately supporting sustainable tree seedling production.

SLE is a systemic autoimmune disease with a large number of common risk gene variants, but several rare gene variants can cause monogenic SLE. The relationship between common and rare variants in SLE is unclear. We therefore investigated the occurrence of rare deleterious variants in patients with childhood-onset SLE (cSLE) and adult-onset SLE (aSLE) and compared the frequency of these variants with their individual SLE polygenic risk score (PRS). Targeted sequencing of 1832 gene regions, including coding regions of 31 genes associated with monogenic SLE, was performed in 958 patients with SLE and 1026 healthy individuals. A total of 116 patients with SLE had disease onset before the age of 18 (cSLE). An SLE common variant PRS was created from 37 SLE genome-wide association study single nucleotide variants (SNVs). Rare coding deleterious SNVs (RD SNVs) were observed in 23 of the monogenic SLE-associated genes. Six per cent of patients with cSLE, compared with 3.2% of controls and 4.6% of patients with aSLE, carried rare deleterious alleles. In cSLE, RD SNVs were observed in the C1S, DDX58, IFIH1, IKZF1, RNASEH2A and C8A genes. A PRS analysis showed that patients with cSLE with any of these gene variants had a similar average PRS as control individuals. RD SNVs were observed in a small proportion of cSLE and carriers of these RD SNVs had a PRS similar to healthy individuals, suggesting the importance of rare coding heterozygous variants in driving disease risk in a subset of children with SLE.

The genetic contributions to early onset heart failure (HF) are incompletely understood. Genetic testing in advanced HF patients undergoing heart transplantation (HTx) may yield clinical benefits, but data is limited. We performed deep-coverage whole genome sequencing (WGS) in 102 Swedish HTx recipients. Gene lists were compiled through a systematic literature review. Variants were prioritized for pathogenicity and classified manually. We also compared polygenic HF risk scores to a population-based cohort. We found a pathogenic (LP/P) variant in 34 individuals (34%). Testing yield was highest in hypertrophic (63% LP/P carriers), dilated (40%) and arrhythmogenic right ventricular (33%) cardiomyopathy and lower in ischemic cardiomyopathy (10%). A family history was more common in LP/P variant carriers than in non-carriers but was present in less than half of carriers (44% vs 13%, P < 0.001), whereas age was similar. Polygenic risk scores were similar in HTx recipients and the population cohort. In conclusion, we observed a high prevalence of pathogenic cardiomyopathy gene variants in individuals with early-onset advanced HF, which could not accurately be ruled out by family history and age. In contrast, we did not observe higher polygenic risk scores in early onset advanced HF cases than in the general population.