Scientific highlights

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.

Spiders produce nature's toughest fiber using renewable components at ambient temperatures and with water as solvent, making it highly interesting to replicate for the materials industry. Despite this, much remains to be understood about the bioprocessing and composition of spider silk fibers. Here, we identify 18 proteins that make up the spiders' strongest silk type, the major ampullate fiber. Single-cell RNA sequencing and spatial transcriptomics revealed that the secretory epithelium of the gland harbors six cell types. These cell types are confined to three distinct glandular zones that produce specific combinations of silk proteins. Image analysis of histological sections showed that the secretions from the three zones do not mix, and proteomics analysis revealed that these secretions form layers in the final fiber. Using a multi-omics approach, we provide substantial advancements in the understanding of the structure and function of the major ampullate silk gland as well as of the architecture and composition of the fiber it produces.

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 was genotyped using STRdust and LongTR to determine the tandem repeat length and sequence composition. PathSTR provides rich visualizations of this dataset 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.

Krill are vital as food for many marine animals but also impacted by global warming. To learn how they and other zooplankton may adapt to a warmer world we studied local adaptation in the widespread Northern krill (Meganyctiphanes norvegica). We assemble and characterize its large genome and compare genome-scale variation among 74 specimens from the colder Atlantic Ocean and warmer Mediterranean Sea. The 19 Gb genome likely evolved through proliferation of retrotransposons, now targeted for inactivation by extensive DNA methylation, and contains many duplicated genes associated with molting and vision. Analysis of 760 million SNPs indicates extensive homogenizing gene-flow among populations. Nevertheless, we detect signatures of adaptive divergence across hundreds of genes, implicated in photoreception, circadian regulation, reproduction and thermal tolerance, indicating polygenic adaptation to light and temperature. The top gene candidate for ecological adaptation was nrf-6, a lipid transporter with a Mediterranean variant that may contribute to early spring reproduction. Such variation could become increasingly important for fitness in Atlantic stocks. Our study underscores the widespread but uneven distribution of adaptive variation, necessitating characterization of genetic variation among natural zooplankton populations to understand their adaptive potential, predict risks and support ocean conservation in the face of climate change.

The B-cell acute lymphoblastic leukemia (ALL) cell line REH, with the t(12;21) ETV6::RUNX1 translocation, is known to have a complex karyotype defined by a series of large-scale chromosomal rearrangements. Taken from a 15-yr-old at relapse, the cell line offers a practical model for the study of pediatric B-ALL. In recent years, short- and long-read DNA and RNA sequencing have emerged as a complement to karyotyping techniques in the resolution of structural variants in an oncological context. Here, we explore the integration of long-read PacBio and Oxford Nanopore whole-genome sequencing, IsoSeq RNA sequencing, and short-read Illumina sequencing to create a detailed genomic and transcriptomic characterization of the REH cell line. Whole-genome sequencing clarified the molecular traits of disrupted ALL-associated genes including CDKN2A, PAX5, BTG1, VPREB1, and TBL1XR1, as well as the glucocorticoid receptor NR3C1 Meanwhile, transcriptome sequencing identified seven fusion genes within the genomic breakpoints. Together, our extensive whole-genome investigation makes high-quality open-source data available to the leukemia genomics community.

Well-differentiated liposarcoma (WDLS) displays amplification of genes on chromosome 12 (Chr12) in supernumerary ring or giant marker chromosomes. These structures have been suggested to develop through chromothripsis, followed by circularization and breakage-fusion-bridge (BFB) cycles. To test this hypothesis, we compared WDLSs with Chr12 amplification in rod-shaped chromosomes with WDLSs with rings. Both types of amplicons share the same spectrum of structural variants (SVs), show higher SV frequencies in Chr12 than in co-amplified segments, have SVs that fuse the telomeric ends of co-amplified chromosomes, and lack interspersed deletions. Combined with the finding of cells with transient rod-shaped structures in tumors with ring chromosomes, this suggests a stepwise process starting with the gain of Chr12 material that, after remodeling which does not fit with classical chromothripsis, forms a dicentric structure with other chromosomes. Depending on if and when telomeres from other chromosomes are captured, circularized or linear gain of 12q sequences will predominate.

Distyly is an iconic floral polymorphism governed by a supergene, which promotes efficient pollen transfer and outcrossing through reciprocal differences in the position of sexual organs in flowers, often coupled with heteromorphic self-incompatibility. Distyly has evolved convergently in multiple flowering plant lineages, but has also broken down repeatedly, often resulting in homostylous, self-compatible populations with elevated rates of self-fertilization. Here, we aimed to study the genetic causes and genomic consequences of the shift to homostyly in Linum trigynum, which is closely related to distylous Linum tenue. Building on a high-quality genome assembly, we show that L. trigynum harbors a genomic region homologous to the dominant haplotype of the distyly supergene conferring long stamens and short styles in L. tenue, suggesting that loss of distyly first occurred in a short-styled individual. In contrast to homostylous Primula and Fagopyrum, L. trigynum harbors no fixed loss-of-function mutations in coding sequences of S-linked distyly candidate genes. Instead, floral gene expression analyses and controlled crosses suggest that mutations downregulating the S-linked LtWDR-44 candidate gene for male self-incompatibility and/or anther height could underlie homostyly and self-compatibility in L. trigynum. Population genomic analyses of 224 whole-genome sequences further demonstrate that L. trigynum is highly self-fertilizing, exhibits significantly lower genetic diversity genome-wide, and is experiencing relaxed purifying selection and less frequent positive selection on nonsynonymous mutations relative to L. tenue. Our analyses shed light on the loss of distyly in L. trigynum, and advance our understanding of a common evolutionary transition in flowering plants.

Long-read genome sequencing (lrGS) is a promising method in genetic diagnostics. Here we investigate the potential of lrGS to detect a disease-associated chromosomal translocation between 17p13 and the 19 centromere. We constructed two sets of phased and non-phased de novo assemblies; (i) based on lrGS only and (ii) hybrid assemblies combining lrGS with optical mapping using lrGS reads with a median coverage of 34X. Variant calling detected both structural variants (SVs) and small variants and the accuracy of the small variant calling was compared with those called with short-read genome sequencing (srGS). The de novo and hybrid assemblies had high quality and contiguity with N50 of 62.85 Mb, enabling a near telomere to telomere assembly with less than a 100 contigs per haplotype. Notably, we successfully identified the centromeric breakpoint of the translocation. A concordance of 92% was observed when comparing small variant calling between srGS and lrGS. In summary, our findings underscore the remarkable potential of lrGS as a comprehensive and accurate solution for the analysis of SVs and small variants. Thus, lrGS could replace a large battery of genetic tests that were used for the diagnosis of a single symptomatic translocation carrier, highlighting the potential of lrGS in the realm of digital karyotyping.

The joint expression of particular colors, morphologies, and behaviors is a common feature of adaptation, but the genetic basis for such "phenotypic syndromes" remains poorly understood. Here, we identified a complex genetic architecture associated with a sexually selected syndrome in common wall lizards, by capitalizing on the adaptive introgression of coloration and morphology into a distantly related lineage. Consistent with the hypothesis that the evolution of phenotypic syndromes in vertebrates is facilitated by developmental linkage through neural crest cells, most of the genes associated with the syndrome are involved in neural crest cell regulation. A major locus was a ~400-kb region, characterized by standing structural genetic variation and previously implied in the evolutionary innovation of coloration and beak size in birds. We conclude that features of the developmental and genetic architecture contribute to maintaining trait integration, facilitating the extensive and rapid introgressive spread of suites of sexually selected characters.

Functional precision medicine (FPM) aims to optimize patient-specific drug selection based on the unique characteristics of their cancer cells. Recent advancements in high throughput ex vivo drug profiling have accelerated interest in FPM. Here, we present a proof-of-concept study for an integrated experimental system that incorporates ex vivo treatment response with a single-cell gene expression output enabling barcoding of several drug conditions in one single-cell sequencing experiment. We demonstrate this through a proof-of-concept investigation focusing on the glucocorticoid-resistant acute lymphoblastic leukemia (ALL) E/R+ Reh cell line. Three different single-cell transcriptome sequencing (scRNA-seq) approaches were evaluated, each exhibiting high cell recovery and accurate tagging of distinct drug conditions. Notably, our comprehensive analysis revealed variations in library complexity, sensitivity (gene detection), and differential gene expression detection across the methods. Despite these differences, we identified a substantial transcriptional response to fludarabine, a highly relevant drug for treating high-risk ALL, which was consistently recapitulated by all three methods. These findings highlight the potential of our integrated approach for studying drug responses at the single-cell level and emphasize the importance of method selection in scRNA-seq studies. Finally, our data encompassing 27 327 cells are freely available to extend to future scRNA-seq methodological comparisons.

Coastal areas are an important source of methane (CH4). However, the exact origins of CH4 in the surface waters of coastal regions, which in turn drive sea-air emissions, remain uncertain. To gain a comprehensive understanding of the current and future climate change feedbacks, it is crucial to identify these CH4 sources and processes that regulate its formation and oxidation. This study investigated coastal CH4 dynamics by comparing water column data from six stations located in the brackish Tvärminne Archipelago, Baltic Sea. The sediment biogeochemistry and microbiology were further investigated at two stations (i.e., nearshore and offshore). These stations differed in terms of stratification, bottom water redox conditions, and organic matter loading. At the nearshore station, CH4 diffusion from the sediment into the water column was negligible, because nearly all CH4 was oxidized within the upper sediment column before reaching the sediment surface. On the other hand, at the offshore station, there was significant benthic diffusion of CH4, albeit the majority underwent oxidation before reaching the sediment-water interface, due to shoaling of the sulfate methane transition zone (SMTZ). The potential contribution of CH4 production in the water column was evaluated and was found to be negligible. After examining the isotopic signatures of δ13C-CH4 across the sediment and water column, it became apparent that the surface water δ13C-CH4 values observed in areas with thermal stratification could not be explained by diffusion, advective fluxes, nor production in the water column. In fact, these values bore a remarkable resemblance to those detected below the SMTZ. This supports the hypothesis that the source of CH4 in surface waters is more likely to originate from ebullition than diffusion in stratified brackish coastal systems.

The Least Nighthawk Chordeiles pusillus is widespread wherever there are savannas in the South American tropics, often in isolated patches, such as white-sands savannas in the Amazon rainforest realm. Here, we investigate genetic relationships between populations of the Least Nighthawk to understand historical processes leading to its diversification and to determine dispersal routes between northern and southern savannas by way of three hypothesized dispersal corridors by comparing samples from white-sand savannas to samples from other savannas outside of the Amazon rainforest region. We use 32 mtDNA samples from the range of C. pusillus to infer a dated phylogeny. In a subset of 17 samples, we use shotgun sequences to infer a distance-based phylogeny and to estimate individual admixture proportions. We calculate gene flow and shared alleles between white-sand and non-Amazonian populations using the ABBA-BABA test (D statistics), and Principal Component Analysis (PCA) to examine genetic structure within and between lineages. Finally, we use species distribution modelling (SDM) of conditions during the Last Glacial Maximum (LGM), currently, and in the future (2050-2080) to predict potential species occurrence under a climate change scenario. Two main clades (estimated to have diverged around 1.07 million years ago) were recovered with mtDNA sequences and Single Nucleotide Polymorphism (SNPs) and were supported by NGSadmix and PCA: one in the Amazon basin white-sand savannas, the other in the non-Amazonian savannas. Possible allele sharing between these clades was indicated by the D-statistics between northern non-Amazonian populations and the white-sand savanna population, but this was not corroborated by the admixture analyses. Dispersal among northern non-Amazonian populations may have occurred in a dry corridor between the Guianan and the Brazilian Shield, which has since moved eastward. Our data suggest that the lineages separated well before the Last Glacial Maximum, consequently dispersal could have happened at any earlier time during similar climatic conditions. Subsequently, non-Amazonian lineages became more divergent among themselves, possibly connecting and dispersing across the mouth of the Amazon River across Marajó island during favourable climatic conditions in the Pleistocene.

The human silencing hub (HUSH) complex binds to transcripts of LINE-1 retrotransposons (L1s) and other genomic repeats, recruiting MORC2 and other effectors to remodel chromatin. How HUSH and MORC2 operate alongside DNA methylation, a central epigenetic regulator of repeat transcription, remains largely unknown. Here we interrogate this relationship in human neural progenitor cells (hNPCs), a somatic model of brain development that tolerates removal of DNA methyltransferase DNMT1. Upon loss of MORC2 or HUSH subunit TASOR in hNPCs, L1s remain silenced by robust promoter methylation. However, genome demethylation and activation of evolutionarily-young L1s attracts MORC2 binding, and simultaneous depletion of DNMT1 and MORC2 causes massive accumulation of L1 transcripts. We identify the same mechanistic hierarchy at pericentromeric α-satellites and clustered protocadherin genes, repetitive elements important for chromosome structure and neurodevelopment respectively. Our data delineate the epigenetic control of repeats in somatic cells, with implications for understanding the vital functions of HUSH-MORC2 in hypomethylated contexts throughout human development.

Chemoreception is crucial for insect fitness, underlying for instance food-, host-, and mate finding. Chemicals in the environment are detected by receptors from three divergent gene families: odorant receptors (ORs), gustatory receptors (GRs), and ionotropic receptors (IRs). However, how the chemoreceptor gene families evolve in parallel with ecological specializations remains poorly understood, especially in the order Coleoptera. Hence, we sequenced the genome and annotated the chemoreceptor genes of the specialised ambrosia beetle Trypodendron lineatum (Coleoptera, Curculionidae, Scolytinae) and compared its chemoreceptor gene repertoires with those of other scolytines with different ecological adaptations, as well as a polyphagous cerambycid species. We identified 67 ORs, 38 GRs, and 44 IRs in T. lineatum ('Tlin'). Across gene families, T. lineatum has fewer chemoreceptors compared to related scolytines, the coffee berry borer Hypothenemus hampei and the mountain pine beetle Dendroctonus ponderosae, and clearly fewer receptors than the polyphagous cerambycid Anoplophora glabripennis. The comparatively low number of chemoreceptors is largely explained by the scarcity of large receptor lineage radiations, especially among the bitter taste GRs and the 'divergent' IRs, and the absence of alternatively spliced GR genes. Only one non-fructose sugar receptor was found, suggesting several sugar receptors have been lost. Also, we found no orthologue in the 'GR215 clade', which is widely conserved across Coleoptera. Two TlinORs are orthologous to ORs that are functionally conserved across curculionids, responding to 2-phenylethanol (2-PE) and green leaf volatiles (GLVs), respectively. Trypodendron lineatum reproduces inside the xylem of decaying conifers where it feeds on its obligate fungal mutualist Phialophoropsis ferruginea. Like previous studies, our results suggest that stenophagy correlates with small chemoreceptor numbers in wood-boring beetles; indeed, the few GRs may be due to its restricted fungal diet. The presence of TlinORs orthologous to those detecting 2-PE and GLVs in other species suggests these compounds are important for T. lineatum. Future functional studies should test this prediction, and chemoreceptor annotations should be conducted on additional ambrosia beetle species to investigate whether few chemoreceptors is a general trait in this specialized group of beetles.

Multiple Sclerosis (MS) is a heterogeneous inflammatory and neurodegenerative disease with an unpredictable course towards progressive disability. Treating progressive MS is challenging due to limited insights into the underlying mechanisms. We examined the molecular changes associated with primary progressive MS (PPMS) using a cross-tissue (blood and post-mortem brain) and multilayered data (genetic, epigenetic, transcriptomic) from independent cohorts. In PPMS, we found hypermethylation of the 1q21.1 locus, controlled by PPMS-specific genetic variations and influencing the expression of proximal genes (CHD1L, PRKAB2) in the brain. Evidence from reporter assay and CRISPR/dCas9 experiments supports a causal link between methylation and expression and correlation network analysis further implicates these genes in PPMS brain processes. Knock-down of CHD1L in human iPSC-derived neurons and knock-out of chd1l in zebrafish led to developmental and functional deficits of neurons. Thus, several lines of evidence suggest a distinct genetic-epigenetic-transcriptional interplay in the 1q21.1 locus potentially contributing to PPMS pathogenesis.

Sexual selection can directly contribute to reproductive isolation and is an important mechanism that can lead to speciation. Lek-mating is one of the most extreme forms of sexual selection, but surprisingly does not seem to preclude occasional hybridization in nature. However, hybridization among lekking species may still be trivial if selection against offspring with intermediate phenotypes prohibits introgression. Here we investigate this further by sequencing the genomes of nearly all bird-of-paradise (Paradisaeidae) species and 10 museum specimens of putative hybrid origin. We find that intergeneric hybridization indeed still takes place despite extreme differentiation in form, plumage, and behavior. In parallel, the genomes of contemporary species contain widespread signatures of past introgression, demonstrating that hybridization has repeatedly resulted in shared genetic variation despite strong sexual isolation. Our study raises important questions about extrinsic factors that modulate hybridization probability and the evolutionary consequences of introgressive hybridization between lekking species.

A number of species have recently recovered from near-extinction. Although these species have avoided the immediate extinction threat, their long-term viability remains precarious due to the potential genetic consequences of population declines, which are poorly understood on a timescale beyond a few generations. Woolly mammoths (Mammuthus primigenius) became isolated on Wrangel Island around 10,000 years ago and persisted for over 200 generations before becoming extinct around 4,000 years ago. To study the evolutionary processes leading up to the mammoths' extinction, we analyzed 21 Siberian woolly mammoth genomes. Our results show that the population recovered quickly from a severe bottleneck and remained demographically stable during the ensuing six millennia. We find that mildly deleterious mutations gradually accumulated, whereas highly deleterious mutations were purged, suggesting ongoing inbreeding depression that lasted for hundreds of generations. The time-lag between demographic and genetic recovery has wide-ranging implications for conservation management of recently bottlenecked populations.

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a complex genomic rearrangement (CGR). Although it has been identified as an important pathogenic DNA mutation signature in genomic disorders and cancer genomes, its architecture remains unresolved. Here, we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the DNA of 24 patients identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted structural variant (SV) haplotypes. Using a combination of short-read genome sequencing (GS), long-read GS, optical genome mapping, and single-cell DNA template strand sequencing (strand-seq), the haplotype structure was resolved in 18 samples. The point of template switching in 4 samples was shown to be a segment of ∼2.2-5.5 kb of 100% nucleotide similarity within inverted repeat pairs. These data provide experimental evidence that inverted low-copy repeats act as recombinant substrates. This type of CGR can result in multiple conformers generating diverse SV haplotypes in susceptible dosage-sensitive loci.

Spinocerebellar ataxia 4 (SCA4), characterized in 1996, features adult-onset ataxia, polyneuropathy, and linkage to chromosome 16q22.1; its underlying mutation has remained elusive. To explore the radiological and neuropathological abnormalities in the entire neuroaxis in SCA4 and search for its mutation. Three Swedish families with undiagnosed ataxia went through clinical, neurophysiological, and neuroimaging tests, including PET studies and genetic investigations. In four cases, neuropathological assessments of the neuroaxis were performed. Genetic testing included short read whole genome sequencing, short tandem repeat analysis with ExpansionHunter de novo, and long read sequencing. Novel features for SCA4 include dysautonomia, motor neuron affection, and abnormal eye movements. We found evidence of anticipation; neuroimaging demonstrated atrophy in the cerebellum, brainstem, and spinal cord. [18F]FDG-PET demonstrated brain hypometabolism and [11C]Flumazenil-PET reduced binding in several brain lobes, insula, thalamus, hypothalamus, and cerebellum. Moderate to severe loss of Purkinje cells in the cerebellum and of motor neurons in the anterior horns of the spinal cord along with pronounced degeneration of posterior tracts was also found. Intranuclear, mainly neuronal, inclusions positive for p62 and ubiquitin were sparse but widespread in the CNS. This finding prompted assessment for nucleotide expansions. A polyglycine stretch encoding GGC expansions in the last exon of the zink finger homeobox 3 gene was identified segregating with disease and not found in 1000 controls. SCA4 is a neurodegenerative disease caused by a novel GGC expansion in the coding region of ZFHX3, and its spectrum is expanded to include dysautonomia and neuromuscular manifestations.