Search

Cheek swabs, heterogeneous samples consisting primarily of buccal epithelial cells, are widely used in pediatric DNA methylation studies and biomarker creation. However, the decrease in buccal proportion with age in adults remains unexamined in childhood. We analyzed cheek swabs from 4626 typically developing children 2-months to 20-years-old. Estimated buccal proportion declined throughout childhood with both increasing chronological and predicted epigenetic age. However, buccal proportion did not associate with age throughout adolescence. Variability in buccal proportion increased with age through the entire developmental range. These trends held inversely true for neutrophil proportions. Correcting for buccal proportion attenuated the weak association with PedBE age acceleration to non-significance during initial estimation. Notably, correcting for buccal proportion attenuated the association of PedBE age acceleration with obsessive-compulsive disorder and strengthened the association with diurnal cortisol slope. Thus, the age-related change in children's oral cells is a crucial consideration for cell type-sensitive research., Competing Interests: Competing interests: The authors declare no competing interests., (© 2025. The Author(s).)

The mechanistic link between the complex mutational landscape of de novo methyltransferase DNMT3A and the pathology of acute myeloid leukemia (AML) has not been clearly elucidated so far. Motivated by a recent discovery of the significance of DNMT3A-destabilizing mutations (DNMT3A INS ) in AML, we here investigate the common characteristics of DNMT3A INS AML methylomes through computational analyses. We present that methylomes of DNMT3A INS AMLs are considerably different from those of DNMT3A R882 AMLs in that they exhibit increased intratumor DNA methylation heterogeneity in bivalent chromatin domains. This epigenetic heterogeneity was associated with the transcriptional variability of developmental and membrane-associated factors shaping stem cell niche, and also was a predictor of the response of AML cells to hypomethylating agents, implying that the survival of AML cells depends on stochastic DNA methylations at bivalent domains. Altogether, our work provides a novel mechanistic model suggesting the genomic origin of the aberrant epigenomic heterogeneity in disease conditions., Competing Interests: Competing interests: Y.K. is the founder and CEO of Genome Opinion. C.H.S. is the executive director and owns stock in Genome Opinion. The other authors declare no competing interests., (© 2025. The Author(s).)

Small molecule machine learning aims to predict chemical, biochemical, or biological properties from molecular structures, with applications such as toxicity prediction, ligand binding, and pharmacokinetics. A recent trend is developing end-to-end models that avoid explicit domain knowledge. These models assume no coverage bias in training and evaluation data, meaning the data are representative of the true distribution. However, the domain of applicability is rarely considered in such models. Here, we investigate how well large-scale datasets cover the space of known biomolecular structures. For doing so, we propose a distance measure based on solving the Maximum Common Edge Subgraph (MCES) problem, which aligns well with chemical similarity. Although this method is computationally hard, we introduce an efficient approach combining Integer Linear Programming and heuristic bounds. Our findings reveal that many widely-used datasets lack uniform coverage of biomolecular structures, limiting the predictive power of models trained on them. We propose two additional methods to assess whether training datasets diverge from known molecular distributions, potentially guiding future dataset creation to improve model performance., Competing Interests: Competing interests: The authors declare no competing interests., (© 2025. The Author(s).)

Cell type deconvolution methods can impute cell proportions from bulk transcriptomics data, revealing changes in disease progression or organ development. But benchmarking studies often use simulated bulk data from the same source as the reference, which limits its application scenarios. This study examines batch effects in deconvolution and introduces SCCAF-D, a computational workflow that ensures a Pearson Correlation Coefficient above 0.75 across simulated and real bulk data for various tissue types. Applied to non-alcoholic fatty liver disease, SCCAF-D unveils meaningful insights into changes in cell proportions during disease progression., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Paratyphoid B fever (PTB) is caused by an invasive lineage (phylogroup 1, PG1) of Salmonella enterica serotype Paratyphi B (SPB). However, little was known about the global population structure, geographic distribution, and evolution of this pathogen. Here, we report a whole-genome analysis of 568 historical and contemporary SPB PG1 isolates, obtained globally, between 1898 and 2021. We show that this pathogen existed in the 13th century, subsequently diversifying into 11 lineages and 38 genotypes with strong phylogeographic patterns. Following its discovery in 1896, it circulated across Europe until the 1970s, after which it was mostly reimported into Europe from South America, the Middle East, South Asia, and North Africa. Antimicrobial resistance recently emerged in various genotypes of SPB PG1, mostly through mutations of the quinolone-resistance-determining regions of gyrA and gyrB. This study provides an unprecedented insight into SPB PG1 and essential genomic tools for identifying and tracking this pathogen, thereby facilitating the global genomic surveillance of PTB., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

RNA-Puzzles is a collective endeavor dedicated to the advancement and improvement of RNA three-dimensional structure prediction. With agreement from structural biologists, RNA structures are predicted by modeling groups before publication of the experimental structures. We report a large-scale set of predictions by 18 groups for 23 RNA-Puzzles: 4 RNA elements, 2 Aptamers, 4 Viral elements, 5 Ribozymes and 8 Riboswitches. We describe automatic assessment protocols for comparisons between prediction and experiment. Our analyses reveal some critical steps to be overcome to achieve good accuracy in modeling RNA structures: identification of helix-forming pairs and of non-Watson-Crick modules, correct coaxial stacking between helices and avoidance of entanglements. Three of the top four modeling groups in this round also ranked among the top four in the CASP15 contest., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Liver regeneration following resection is a complex process relying on coordinated pathways and cell types in the remnant organ. Myeloid-Derived Suppressor Cells (MDSCs) have a role in liver regeneration-related angiogenesis but other roles they may play in this process remain to be elucidated. In this study, we sought to examine the effect of G-MDSCs on hepatocytes proliferation and immune modulation during liver regeneration. Global gene expression profiling of regenerating hepatocytes in mice with CD11b + Ly6G + MDSCs (G-MDSCs) depletion revealed disrupted transcriptional progression from day one to day two after major liver resection. Key genes and pathways related to hepatocyte proliferation and immune response were differentially expressed upon MDSC depletion. Hepatocytes cellularity increased when co-cultured with G-MDSCs, or treated with amphiregulin, which G-MDSCs upregulate during regeneration. Cytometry by time-of-flight (CyTOF) analysis of the intra-liver immune milieu upon MDSC depletion during regeneration demonstrated increased natural killer cell proportions, alongside changes in other immune cell populations. Taken together, these results provide evidence that MDSCs contribute to early liver regeneration by promoting hepatocyte proliferation and modulating the intra-liver immune response, and illuminate the multifaceted role of MDSCs in liver regeneration., Competing Interests: Competing interests: The authors declare no competing interests. Ethics approval and consent to participate: All methods were performed in accordance with the relevant guidelines and regulations. All experimental procedures in animals were conducted according to protocols approved by the institutional animal care and use committee (IACUC) of Sheba Medical Center, strictly observing national, local, and institutional guidelines. Approval no. 1277-20-ANIM., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Tidal marshes are threatened coastal ecosystems known for their capacity to store large amounts of carbon in their water-logged soils. Accurate quantification and mapping of global tidal marshes soil organic carbon (SOC) stocks is of considerable value to conservation efforts. Here, we used training data from 3710 unique locations, landscape-level environmental drivers and a global tidal marsh extent map to produce a global, spatially explicit map of SOC storage in tidal marshes at 30 m resolution. Here we show the total global SOC stock to 1 m to be 1.44 Pg C, with a third of this value stored in the United States of America. On average, SOC in tidal marshes' 0-30 and 30-100 cm soil layers are estimated at 83.1 Mg C ha -1 (average predicted error 44.8 Mg C ha -1 ) and 185.3 Mg C ha -1 (average predicted error 105.7 Mg C ha -1 ), respectively., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Phase-separating peptides (PSPs) self-assembling into coacervate microdroplets (CMs) are a promising class of intracellular delivery vehicles that can release macromolecular modalities deployed in a wide range of therapeutic treatments. However, the molecular grammar governing intracellular uptake and release kinetics of CMs remains elusive. Here, we systematically manipulate the sequence of PSPs to unravel the relationships between their molecular structure, the physical properties of the resulting CMs, and their delivery efficacy. We show that a few amino acid alterations are sufficient to modulate the viscoelastic properties of CMs towards either a gel-like or a liquid-like state as well as their binding interaction with cellular membranes, collectively enabling to tune the kinetics of intracellular cargo release. We also demonstrate that the optimized PSPs CMs display excellent transfection efficiency in hard-to-transfect cells such as primary fibroblasts and immune cells. Our findings provide molecular guidelines to precisely program the material properties of PSP CMs and achieve tunable cellular uptake and release kinetics depending on the cargo modality, with broad implications for therapeutic applications such as protein, gene, and immune cell therapies., Competing Interests: Competing interests: Y.S. and A.M. have filed a PCT application on the peptide coacervates described in this study (PCT application no. PCT/SG2023/11401). The remaining authors declare no competing interests., (© 2024. The Author(s).)

ABCB6 has been implicated in dyschromatosis universalis hereditaria, a condition characterized by hyperpigmented and hypopigmented skin macules. Dyschromatosis universalis hereditaria can also present with hearing loss. Dyschromatosis universalis hereditaria-associated mutations in ABCB6 have been reported, but the role of this protein in the inner ear has not been studied. Here we determine a high-resolution (2.93 Å) cryo-EM structure of ABCB6 and functionally characterized several dyschromatosis universalis hereditaria mutants. We find that the L356P mutant abolishes ABCB6 function, and affirm the underlying loss of ATP binding mechanism using molecular dynamics simulations based on our cryo-EM structure. To test the role of ABCB6 in the inner ear, we characterize Abcb6 (the ABCB6 homolog) in zebrafish. We show that Abcb6 suppression by morpholinos reduces inner ear and lateral line hair cell numbers. Morphants also lack the utricular otolith, which is associated with vestibular function. Co-injecting morpholinos with human ABCB6 mRNA partially rescues the morphant phenotype, suggesting that Abcb6 plays a developmental role in inner ear structures. Further, we show that Abcb6 knockout mice exhibit an increased auditory brainstem response threshold, resulting in reduced hearing sensitivity. Taken together, these data suggest ABCB6 plays a role in inner ear development and function., Competing Interests: Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

Background: Body mass index (BMI) trajectories have been used to assess the growth of children with respect to their peers, and to anticipate future obesity and disease risk. While retrospective BMI trajectories have been actively studied, models to prospectively predict continuous BMI trajectories have not been investigated., Materials and Methods: Using longitudinal BMI measurements between birth and age 10 y from a mother-offspring cohort, we leveraged a multi-task Gaussian process approach to develop and evaluate a unified framework for modeling, clustering, and prospective prediction of BMI trajectories. We compared its sensitivity to missing values in the longitudinal follow-up of children, compared its prediction performance to cubic B-spline and multilevel Jenss-Bayley models, and used prospectively predicted BMI trajectories to assess the probability of future BMIs crossing the clinical cutoffs for obesity., Results: MagmaClust identified 5 distinct patterns of BMI trajectories between 0 to 10 y. The method outperformed both cubic B-spline and multilevel Jenss-Bayley models in the accuracy of retrospective BMI trajectories while being more robust to missing data (up to 90%). It was also better at prospectively forecasting BMI trajectories of children for periods ranging from 2 to 8 years into the future, using historic BMI data. Given BMI data between birth and age 2 years, prediction of overweight/obesity status at age 10 years, as computed from MagmaClust's predictions exhibited high specificity (0.94), negative predictive value (0.89), and accuracy (0.86). The accuracy, sensitivity, and positive predictive value of predictions increased as BMI data from additional time points were utilized for prediction., Conclusion: MagmaClust provides a unified, probabilistic, non-parametric framework to model, cluster, and prospectively predict childhood BMI trajectories and overweight/obesity risk. The proposed method offers a convenient tool for clinicians to monitor BMI growth in children, allowing them to prospectively identify children with high predicted overweight/obesity risk and implement timely interventions., Competing Interests: Competing interests: K.M.G., Y.S.C., and F.Y. received reimbursement for speaking at conferences sponsored by companies selling nutritional products. K.M.G., and Y.S.C. are part of an academic consortium that has received research funding from Société Des Produits Nestlé S.A. and BenevolentAI Bio Ltd, and are co-inventors on patents filed on nutritional factors and metabolic risk outside the submitted work. All other authors declare that they have nothing to disclose. Ethics approval: This study was approved by both the National Healthcare Group Domain Specific Review Board (D/2009/021 & B/2014/00411) and the SingHealth Centralized Institutional Review Board (2018/2767 & 2019/2406). It was conducted in accordance with the ethical standards set forth in the 1964 Declaration of Helsinki and its later amendments Written informed consent was obtained from all children and their parents., (© 2024. The Author(s).)

Several cardiovascular traits and diseases co-occur with Alzheimer's disease. We mapped their shared genetic architecture using multi-trait genome-wide association studies. Subsequent fine-mapping and colocalisation highlighted 16 genetic loci associated with both Alzheimer's and cardiovascular diseases. We prioritised rs11786896, which colocalised with Alzheimer's disease, atrial fibrillation and expression of PLEC in the heart left ventricle, and rs7529220, which colocalised with Alzheimer's disease, atrial fibrillation and expression of C1Q family genes. Single-cell RNA-sequencing data, co-expression network and protein-protein interaction analyses provided evidence for different mechanisms of PLEC, which is upregulated in left ventricular endothelium and cardiomyocytes with heart failure and in brain astrocytes with Alzheimer's disease. Similar common mechanisms are implicated for C1Q in heart macrophages with heart failure and in brain microglia with Alzheimer's disease. These findings highlight inflammatory and pleomorphic risk determinants for the co-occurrence of Alzheimer's and cardiovascular diseases and suggest PLEC, C1Q and their interacting proteins as potential therapeutic targets., Competing Interests: Competing interests The authors declare no comparing interests., (© 2024. The Author(s).)

Patient-derived xenografts (PDXs) are tumour fragments engrafted into mice for preclinical studies. PDXs offer clear advantages over simpler in vitro cancer models - such as cancer cell lines (CCLs) and organoids - in terms of structural complexity, heterogeneity, and stromal interactions. Here, we characterise 231 colorectal cancer PDXs at the genomic, transcriptomic, and epigenetic levels, along with their response to cetuximab, an EGFR inhibitor used clinically for metastatic colorectal cancer. After evaluating the PDXs' quality, stability, and molecular concordance with publicly available patient cohorts, we present results from training, interpreting, and validating the integrative ensemble classifier CeSta. This model takes in input the PDXs' multi-omic characterisation and predicts their sensitivity to cetuximab treatment, achieving an area under the receiver operating characteristics curve > 0.88. Our study demonstrates that large PDX collections can be leveraged to train accurate, interpretable drug sensitivity models that: (1) better capture patient-derived therapeutic biomarkers compared to models trained on CCL data, (2) can be robustly validated across independent PDX cohorts, and (3) could contribute to the development of future therapeutic biomarkers., Competing Interests: Competing interests FI receives funding from Open Targets, a public-private initiative involving academia and industry, and from Nerviano Medical Sciences and performs consultancy for the joint Cancer Research Horizon— AstraZeneca Functional Genomics Centre and for Mosaic T.X., L.T. has received research grants from Menarini, Merck KGaA, Merus, Pfizer, Servier and Symphogen. U.P. is a consultant for Omniscope Inc. H.K. and J.S. are employee of Charles River. U.M. is an employee and holder of company stock of AstraZeneca. All the other authors declare no competing interests., (© 2024. The Author(s).)

Structural variants (SVs) are significant contributors to inter-individual genetic variation associated with traits and diseases. Current SV studies using whole-genome sequencing (WGS) have a largely Eurocentric composition, with little known about SV diversity in other ancestries, particularly from Asia. Here, we present a WGS catalogue of 73,035 SVs from 8392 Singaporeans of East Asian, Southeast Asian and South Asian ancestries, of which ~65% (47,770 SVs) are novel. We show that Asian populations can be stratified by their global SV patterns and identified 42,239 novel SVs that are specific to Asian populations. 52% of these novel SVs are restricted to one of the three major ancestry groups studied (Indian, Chinese or Malay). We uncovered SVs affecting major clinically actionable loci. Lastly, by identifying SVs in linkage disequilibrium with single-nucleotide variants, we demonstrate the utility of our SV catalogue in the fine-mapping of Asian GWAS variants and identification of potential causative variants. These results augment our knowledge of structural variation across human populations, thereby reducing current ancestry biases in global references of genetic variation afflicting equity, diversity and inclusion in genetic research., (© 2024. The Author(s).)

Spatial transcriptomics produces high-dimensional gene expression measurements with spatial context. Obtaining a biologically meaningful low-dimensional representation of such data is crucial for effective interpretation and downstream analysis. Here, we present Spatial Transcriptomics Analysis with topic Modeling to uncover spatial Patterns (STAMP), an interpretable spatially aware dimension reduction method built on a deep generative model that returns biologically relevant, low-dimensional spatial topics and associated gene modules. STAMP can analyze data ranging from a single section to multiple sections and from different technologies to time-series data, returning topics matching known biological domains and associated gene modules containing established markers highly ranked within. In a lung cancer sample, STAMP delineated cell states with supporting markers at a higher resolution than the original annotation and uncovered cancer-associated fibroblasts concentrated on the tumor edge's exterior. In time-series data of mouse embryonic development, STAMP disentangled the erythro-myeloid hematopoiesis and hepatocytes developmental trajectories within the liver. STAMP is highly scalable and can handle more than 500,000 cells., (© 2024. The Author(s).)

To shed light on the enigmatic origin of the vertebrate head, our study employs an integrated approach that combines single-cell transcriptomics, perturbations in signaling pathways, and cis-regulatory analysis in amphioxus. As a representative of a basal lineage within the chordate phylum, amphioxus retains many characteristics thought to have been present in the common chordate ancestor. Through cell type characterization, we identify the presence of prechordal plate-like, pre-migratory, and migratory neural crest-like cell populations in the developing amphioxus embryo. Functional analysis establishes conserved roles of the Nodal and Hedgehog signaling pathways in prechordal plate-like populations, and of the Wnt signaling pathway in neural crest-like populations' development. Furthermore, our trans-species transgenic experiments highlight similarities in the regulatory environments that drive neural crest-like and prechordal plate-like developmental programs in both vertebrates and amphioxus. Our findings provide evidence that the key features of vertebrate head development can be traced back to the common ancestor of all chordates., (© 2024. The Author(s).)

Decades of genetic association testing in human cohorts have provided important insights into the genetic architecture and biological underpinnings of complex traits and diseases. However, for certain traits, genome-wide association studies (GWAS) for common SNPs are approaching signal saturation, which underscores the need to explore other types of genetic variation to understand the genetic basis of traits and diseases. Copy number variation (CNV) is an important source of heritability that is well known to functionally affect human traits. Recent technological and computational advances enable the large-scale, genome-wide evaluation of CNVs, with implications for downstream applications such as polygenic risk scoring and drug target identification. Here, we review the current state of CNV-GWAS, discuss current limitations in resource infrastructure that need to be overcome to enable the wider uptake of CNV-GWAS results, highlight emerging opportunities and suggest guidelines and standards for future GWAS for genetic variation beyond SNPs at scale., (© 2024. Springer Nature Limited.)

Although rare genetic conditions are mostly caused by DNA sequence alterations that functionally disrupt individual genes, large-scale studies using genome sequencing have started to unmask additional complexity. Understanding how combinations of variants in different genes shape human phenotypes is expected to provide important insights into the clinical and genetic heterogeneity of rare disorders. Here, we use albinism, an archetypal rare condition associated with hypopigmentation, as an exemplar for the study of genetic interactions. We analyse data from the Genomics England 100,000 Genomes Project alongside a cohort of 1120 individuals with albinism, and investigate the effect of dual heterozygosity for the combination of two established albinism-related variants: TYR:c.1205 G > A (p.Arg402Gln) [rs1126809] and OCA2:c.1327 G > A (p.Val443Ile) [rs74653330]. As each of these changes alone is insufficient to cause disease when present in the heterozygous state, we sought evidence of synergistic effects. We show that, when both variants are present, the probability of receiving a diagnosis of albinism is significantly increased (odds ratio 12.8; 95% confidence interval 6.0 - 24.7; p-value 2.1 ×10 -8 ). Further analyses in an independent cohort, the UK Biobank, support this finding and highlight that heterozygosity for the TYR:c.1205 G > A and OCA2:c.1327 G > A variant combination is associated with statistically significant alterations in visual acuity and central retinal thickness (traits that are considered albinism endophenotypes). The approach discussed in this report opens up new avenues for the investigation of oligogenic patterns in apparently Mendelian disorders., (© 2024. The Author(s).)

Modern sequencing technology enables the systematic detection of complex structural variation (SV) across genomes. However, extensive DNA rearrangements arising through a series of mutations, a phenomenon we refer to as serial SV (sSV), remain underexplored, posing a challenge for SV discovery. Here, we present NAHRwhals ( https://github.com/WHops/NAHRwhals ), a method to infer repeat-mediated series of SVs in long-read genomic assemblies. Applying NAHRwhals to haplotype-resolved human genomes from 28 individuals reveals 37 sSV loci of various length and complexity. These sSVs explain otherwise cryptic variation in medically relevant regions such as the TPSAB1 gene, 8p23.1, 22q11 and Sotos syndrome regions. Comparisons with great ape assemblies indicate that most human sSVs formed recently, after the human-ape split, and involved non-repeat-mediated processes in addition to non-allelic homologous recombination. NAHRwhals reliably discovers and characterizes sSVs at scale and independent of species, uncovering their genomic abundance and suggesting broader implications for disease., (© 2024. The Author(s).)

The use of monoclonal antibodies for the control of drug resistant nosocomial bacteria may alleviate a reliance on broad spectrum antimicrobials for treatment of infection. We identify monoclonal antibodies that may prevent infection caused by carbapenem resistant Acinetobacter baumannii. We use human immune repertoire mice (Kymouse platform mice) as a surrogate for human B cell interrogation to establish an unbiased strategy to probe the antibody-accessible target landscape of clinically relevant A. baumannii. After immunisation of the Kymouse platform mice with A. baumannii derived outer membrane vesicles (OMV) we identify 297 antibodies and analyse 26 of these for functional potential. These antibodies target lipooligosaccharide (OCL1), the Oxa-23 protein, and the KL49 capsular polysaccharide. We identify a single monoclonal antibody (mAb1416) recognising KL49 capsular polysaccharide to demonstrate prophylactic in vivo protection against a carbapenem resistant A. baumannii lineage associated with neonatal sepsis mortality in Asia. Our end-to-end approach identifies functional monoclonal antibodies with prophylactic potential against major lineages of drug resistant bacteria accounting for phylogenetic diversity and clinical relevance without existing knowledge of a specific target antigen. Such an approach might be scaled for a additional clinically important bacterial pathogens in the post-antimicrobial era., (© 2024. The Author(s).)

Genetic parasites, including viruses and transposons, exploit components from the host for their own replication. However, little is known about virus-transposon interactions within host cells. Here, we discover a strategy where human cytomegalovirus (HCMV) hijacks L1 retrotransposon encoded protein during its replication cycle. HCMV infection upregulates L1 expression by enhancing both the expression of L1-activating transcription factors, YY1 and RUNX3, and the chromatin accessibility of L1 promoter regions. Increased L1 expression, in turn, promotes HCMV replicative fitness. Affinity proteomics reveals UL44, HCMV DNA polymerase subunit, as the most abundant viral binding protein of the L1 ribonucleoprotein (RNP) complex. UL44 directly interacts with L1 ORF2p, inducing DNA damage responses in replicating HCMV compartments. While increased L1-induced mutagenesis is not observed in HCMV for genetic adaptation, the interplay between UL44 and ORF2p accelerates viral DNA replication by alleviating replication stress. Our findings shed light on how HCMV exploits host retrotransposons for enhanced viral fitness., (© 2024. The Author(s).)

Cancer is a highly heterogeneous disease, where phenotypically distinct subpopulations coexist and can be primed to different fates. Both genetic and epigenetic factors may drive cancer evolution, however little is known about whether and how such a process is pre-encoded in cancer clones. Using single-cell multi-omic lineage tracing and phenotypic assays, we investigate the predictive features of either tumour initiation or drug tolerance within the same cancer population. Clones primed to tumour initiation in vivo display two distinct transcriptional states at baseline. Remarkably, these states share a distinctive DNA accessibility profile, highlighting an epigenetic basis for tumour initiation. The drug tolerant niche is also largely pre-encoded, but only partially overlaps the tumour-initiating one and evolves following two genetically and transcriptionally distinct trajectories. Our study highlights coexisting genetic, epigenetic and transcriptional determinants of cancer evolution, unravelling the molecular complexity of pre-encoded tumour phenotypes., (© 2024. The Author(s).)

The volume of public proteomics data is rapidly increasing, causing a computational challenge for large-scale reanalysis. Here, we introduce quantms ( https://quant,ms.org/ ), an open-source cloud-based pipeline for massively parallel proteomics data analysis. We used quantms to reanalyze 83 public ProteomeXchange datasets, comprising 29,354 instrument files from 13,132 human samples, to quantify 16,599 proteins based on 1.03 million unique peptides. quantms is based on standard file formats improving the reproducibility, submission and dissemination of the data to ProteomeXchange., (© 2024. The Author(s).)

There are two main families of G protein-coupled receptors that detect odours in humans, the odorant receptors (ORs) and the trace amine-associated receptors (TAARs). Their amino acid sequences are distinct, with the TAARs being most similar to the aminergic receptors such as those activated by adrenaline, serotonin, dopamine and histamine. To elucidate the structural determinants of ligand recognition by TAARs, we have determined the cryo-EM structure of a murine receptor, mTAAR7f, coupled to the heterotrimeric G protein G s and bound to the odorant N,N-dimethylcyclohexylamine (DMCHA) to an overall resolution of 2.9 Å. DMCHA is bound in a hydrophobic orthosteric binding site primarily through van der Waals interactions and a strong charge-charge interaction between the tertiary amine of the ligand and an aspartic acid residue. This site is distinct and non-overlapping with the binding site for the odorant propionate in the odorant receptor OR51E2. The structure, in combination with mutagenesis data and molecular dynamics simulations suggests that the activation of the receptor follows a similar pathway to that of the β-adrenoceptors, with the significant difference that DMCHA interacts directly with one of the main activation microswitch residues, Trp 6.48 ., (© 2024. The Author(s).)

Quantifying the contribution of genetics and environmental effects on disease initiation and progression, as well as the shared genetics of different diseases, is vital for the understanding of the disease etiology of multimorbidities. In this study, we leverage nationwide Danish registries to provide a granular atlas of the genetic origin of disease phenotypes for a cohort of all Danes 1978-2018 with partially known pedigree (n = 6.3 million). We estimate the heritability and genetic correlation between thousands of disease phenotypes using a novel approach that can be scaled to nationwide data. Our findings confirm the importance of genetics for a number of known associations and increase the resolution of heritability by adding numerous associations, some of which point to shared biologically origin of different phenotypes. We also establish the heritability of disease trajectories and the importance of sex-specific genetic contributions. Results can be accessed at https://h2.cpr.ku.dk/ ., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

The Structural Genomics Consortium is an international open science research organization with a focus on accelerating early-stage drug discovery, namely hit discovery and optimization. We, as many others, believe that artificial intelligence (AI) is poised to be a main accelerator in the field. The question is then how to best benefit from recent advances in AI and how to generate, format and disseminate data to enable future breakthroughs in AI-guided drug discovery. We present here the recommendations of a working group composed of experts from both the public and private sectors. Robust data management requires precise ontologies and standardized vocabulary while a centralized database architecture across laboratories facilitates data integration into high-value datasets. Lab automation and opening electronic lab notebooks to data mining push the boundaries of data sharing and data modeling. Important considerations for building robust machine-learning models include transparent and reproducible data processing, choosing the most relevant data representation, defining the right training and test sets, and estimating prediction uncertainty. Beyond data-sharing, cloud-based computing can be harnessed to build and disseminate machine-learning models. Important vectors of acceleration for hit and chemical probe discovery will be (1) the real-time integration of experimental data generation and modeling workflows within design-make-test-analyze (DMTA) cycles openly, and at scale and (2) the adoption of a mindset where data scientists and experimentalists work as a unified team, and where data science is incorporated into the experimental design., (© 2024. The Author(s).)

Mitotic errors generate micronuclei entrapping mis-segregated chromosomes, which are susceptible to catastrophic fragmentation through chromothripsis. The reassembly of fragmented chromosomes by error-prone DNA double-strand break (DSB) repair generates diverse genomic rearrangements associated with human diseases. How specific repair pathways recognize and process these lesions remains poorly understood. Here we use CRISPR/Cas9 to systematically inactivate distinct DSB repair pathways and interrogate the rearrangement landscape of fragmented chromosomes. Deletion of canonical non-homologous end joining (NHEJ) components substantially reduces complex rearrangements and shifts the rearrangement landscape toward simple alterations without the characteristic patterns of chromothripsis. Following reincorporation into the nucleus, fragmented chromosomes localize within sub-nuclear micronuclei bodies (MN bodies) and undergo ligation by NHEJ within a single cell cycle. In the absence of NHEJ, chromosome fragments are rarely engaged by alternative end-joining or recombination-based mechanisms, resulting in delayed repair kinetics, persistent 53BP1-labeled MN bodies, and cell cycle arrest. Thus, we provide evidence supporting NHEJ as the exclusive DSB repair pathway generating complex rearrangements from mitotic errors., (© 2024. The Author(s).)

ezSingleCell is an interactive and easy-to-use application for analysing various single-cell and spatial omics data types without requiring prior programing knowledge. It combines the best-performing publicly available methods for in-depth data analysis, integration, and interactive data visualization. ezSingleCell consists of five modules, each designed to be a comprehensive workflow for one data type or task. In addition, ezSingleCell allows crosstalk between different modules within a unified interface. Acceptable input data can be in a variety of formats while the output consists of publication ready figures and tables. In-depth manuals and video tutorials are available to guide users on the analysis workflows and parameter adjustments to suit their study aims. ezSingleCell's streamlined interface can analyse a standard scRNA-seq dataset of 3000 cells in less than five minutes. ezSingleCell is available in two forms: an installation-free web application ( https://immunesinglecell.org/ezsc/ ) or a software package with a shinyApp interface ( https://github.com/JinmiaoChenLab/ezSingleCell2 ) for offline analysis., (© 2024. The Author(s).)

The Long-read RNA-Seq Genome Annotation Assessment Project Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. Using different protocols and sequencing platforms, the consortium generated over 427 million long-read sequences from complementary DNA and direct RNA datasets, encompassing human, mouse and manatee species. Developers utilized these data to address challenges in transcript isoform detection, quantification and de novo transcript detection. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. Incorporating additional orthogonal data and replicate samples is advised when aiming to detect rare and novel transcripts or using reference-free approaches. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis., (© 2024. The Author(s).)

Real-time genomics through nanopore sequencing holds the promise of fast antibiotic resistance prediction directly in the clinical setting. However, concerns about the accuracy of genomics-based resistance predictions persist, particularly when compared to traditional, clinically established diagnostic methods. Here, we leverage the case of a multi-drug resistant Klebsiella pneumoniae infection to demonstrate how real-time genomics can enhance the accuracy of antibiotic resistance profiling in complex infection scenarios. Our results show that unlike established diagnostics, nanopore sequencing data analysis can accurately detect low-abundance plasmid-mediated resistance, which often remains undetected by conventional methods. This capability has direct implications for clinical practice, where such "hidden" resistance profiles can critically influence treatment decisions. Consequently, the rapid, in situ application of real-time genomics holds significant promise for improving clinical decision-making and patient outcomes., (© 2024. The Author(s).)

HNF4A and HNF1A encode transcription factors that are important for the development and function of the pancreas and liver. Mutations in both genes have been directly linked to Maturity Onset Diabetes of the Young (MODY) and type 2 diabetes (T2D) risk. To better define the pleiotropic gene regulatory roles of HNF4A and HNF1A, we generated a comprehensive genome-wide map of their binding targets in pancreatic and hepatic cells using ChIP-Seq. HNF4A was found to bind and regulate known (ACY3, HAAO, HNF1A, MAP3K11) and previously unidentified (ABCD3, CDKN2AIP, USH1C, VIL1) loci in a tissue-dependent manner. Functional follow-up highlighted a potential role for HAAO and USH1C as regulators of beta cell function. Unlike the loss-of-function HNF4A/MODY1 variant I271fs, the T2D-associated HNF4A variant (rs1800961) was found to activate AKAP1, GAD2 and HOPX gene expression, potentially due to changes in DNA-binding affinity. We also found HNF1A to bind to and regulate GPR39 expression in beta cells. Overall, our studies provide a rich resource for uncovering downstream molecular targets of HNF4A and HNF1A that may contribute to beta cell or hepatic cell (dys)function, and set up a framework for gene discovery and functional validation., (© 2024. The Author(s).)

Traumatic brain injury (TBI) can result in long-lasting changes in hippocampal function. The changes induced by TBI on the hippocampus contribute to cognitive deficits. The adult hippocampus harbors neural stem cells (NSCs) that generate neurons (neurogenesis), and astrocytes (astrogliogenesis). While deregulation of hippocampal NSCs and neurogenesis have been observed after TBI, it is not known how TBI may affect hippocampal astrogliogenesis. Using a controlled cortical impact model of TBI in male mice, single cell RNA sequencing and spatial transcriptomics, we assessed how TBI affected hippocampal NSCs and the neuronal and astroglial lineages derived from them. We observe an increase in NSC-derived neuronal cells and a concomitant decrease in NSC-derived astrocytic cells, together with changes in gene expression and cell dysplasia within the dentate gyrus. Here, we show that TBI modifies NSC fate to promote neurogenesis at the cost of astrogliogenesis and identify specific cell populations as possible targets to counteract TBI-induced cellular changes in the adult hippocampus., (© 2024. The Author(s).)

Cell segmentation is a critical step for quantitative single-cell analysis in microscopy images. Existing cell segmentation methods are often tailored to specific modalities or require manual interventions to specify hyper-parameters in different experimental settings. Here, we present a multimodality cell segmentation benchmark, comprising more than 1,500 labeled images derived from more than 50 diverse biological experiments. The top participants developed a Transformer-based deep-learning algorithm that not only exceeds existing methods but can also be applied to diverse microscopy images across imaging platforms and tissue types without manual parameter adjustments. This benchmark and the improved algorithm offer promising avenues for more accurate and versatile cell analysis in microscopy imaging., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Lung cancer is the second most frequently diagnosed cancer and the leading cause of cancer-related mortality worldwide. Tumour ecosystems feature diverse immune cell types. Myeloid cells, in particular, are prevalent and have a well-established role in promoting the disease. In our study, we profile approximately 900,000 cells from 25 treatment-naive patients with adenocarcinoma and squamous-cell carcinoma by single-cell and spatial transcriptomics. We note an inverse relationship between anti-inflammatory macrophages and NK cells/T cells, and with reduced NK cell cytotoxicity within the tumour. While we observe a similar cell type composition in both adenocarcinoma and squamous-cell carcinoma, we detect significant differences in the co-expression of various immune checkpoint inhibitors. Moreover, we reveal evidence of a transcriptional "reprogramming" of macrophages in tumours, shifting them towards cholesterol export and adopting a foetal-like transcriptional signature which promotes iron efflux. Our multi-omic resource offers a high-resolution molecular map of tumour-associated macrophages, enhancing our understanding of their role within the tumour microenvironment., (© 2024. The Author(s).)

T cells are essential immune cells responsible for identifying and eliminating pathogens. Through interactions between their T-cell antigen receptors (TCRs) and antigens presented by major histocompatibility complex molecules (MHCs) or MHC-like molecules, T cells discriminate foreign and self peptides. Determining the fundamental principles that govern these interactions has important implications in numerous medical contexts. However, reconstructing a map between T cells and their antagonist antigens remains an open challenge for the field of immunology, and success of in silico reconstructions of this relationship has remained incremental. In this Perspective, we discuss the role that new state-of-the-art deep-learning models for predicting protein structure may play in resolving some of the unanswered questions the field faces linking TCR and peptide-MHC properties to T-cell specificity. We provide a comprehensive overview of structural databases and the evolution of predictive models, and highlight the breakthrough AlphaFold provided the field., (© 2024. Springer Nature America, Inc.)

Light-driven sodium pumps (NaRs) are unique ion-transporting microbial rhodopsins. The major group of NaRs is characterized by an NDQ motif and has two aspartic acid residues in the central region essential for sodium transport. Here we identify a subgroup of the NDQ rhodopsins bearing an additional glutamic acid residue in the close vicinity to the retinal Schiff base. We thoroughly characterize a member of this subgroup, namely the protein ErNaR from Erythrobacter sp. HL-111 and show that the additional glutamic acid results in almost complete loss of pH sensitivity for sodium-pumping activity, which is in contrast to previously studied NaRs. ErNaR is capable of transporting sodium efficiently even at acidic pH levels. X-ray crystallography and single particle cryo-electron microscopy reveal that the additional glutamic acid residue mediates the connection between the other two Schiff base counterions and strongly interacts with the aspartic acid of the characteristic NDQ motif. Hence, it reduces its pKa. Our findings shed light on a subgroup of NaRs and might serve as a basis for their rational optimization for optogenetics., (© 2024. The Author(s).)

Meniere Disease (MD) is a chronic inner ear disorder characterized by vertigo attacks, sensorineural hearing loss, tinnitus, and aural fullness. Extensive evidence supporting the inflammatory etiology of MD has been found, therefore, by using transcriptome analysis, we aim to describe the inflammatory variants of MD. We performed Bulk RNAseq on 45 patients with definite MD and 15 healthy controls. MD patients were classified according to their basal levels of IL-1β into 2 groups: high and low. Differentially expression analysis was performed using the ExpHunter Suite, and cell type proportion was evaluated using the estimation algorithms xCell, ABIS, and CIBERSORTx. MD patients showed 15 differentially expressed genes (DEG) compared to controls. The top DEGs include IGHG1 (p = 1.64  ×  10 -6 ) and IGLV3-21 (p = 6.28  ×  10 -3 ), supporting a role in the adaptative immune response. Cytokine profiling defines a subgroup of patients with high levels of IL-1β with up-regulation of IL6 (p = 7.65  ×  10 -8 ) and INHBA (p = 3.39  ×  10 -7 ) genes. Transcriptomic data from peripheral blood mononuclear cells support a proinflammatory subgroup of MD patients with high levels of IL6 and an increase in naïve B-cells, and memory CD8 + T cells., (© 2024. The Author(s).)

Streptococcus dysgalactiae subsp. equisimilis (SDSE) is an emerging cause of human infection with invasive disease incidence and clinical manifestations comparable to the closely related species, Streptococcus pyogenes. Through systematic genomic analyses of 501 disseminated SDSE strains, we demonstrate extensive overlap between the genomes of SDSE and S. pyogenes. More than 75% of core genes are shared between the two species with one third demonstrating evidence of cross-species recombination. Twenty-five percent of mobile genetic element (MGE) clusters and 16 of 55 SDSE MGE insertion regions were shared across species. Assessing potential cross-protection from leading S. pyogenes vaccine candidates on SDSE, 12/34 preclinical vaccine antigen genes were shown to be present in >99% of isolates of both species. Relevant to possible vaccine evasion, six vaccine candidate genes demonstrated evidence of inter-species recombination. These findings demonstrate previously unappreciated levels of genomic overlap between these closely related pathogens with implications for streptococcal pathobiology, disease surveillance and prevention., (© 2024. The Author(s).)

RNA structure is critical for multiple steps in gene regulation. However, how the structures of transcripts differ both within and between individual cells is unknown. Here we develop a SHAPE-inspired method called single-cell structure probing of RNA transcripts that enables simultaneous determination of transcript secondary structure and abundance at single-cell resolution. We apply single-cell structure probing of RNA transcripts to human embryonic stem cells and differentiating neurons. Remarkably, RNA structure is more homogeneous in human embryonic stem cells compared with neurons, with the greatest homogeneity found in coding regions. More extensive heterogeneity is found within 3' untranslated regions and is determined by specific RNA-binding proteins. Overall RNA structure profiles better discriminate cell type identity and differentiation stage than gene expression profiles alone. We further discover a cell-type variable region of 18S ribosomal RNA that is associated with cell cycle and translation control. Our method opens the door to the systematic characterization of RNA structure-function relationships at single-cell resolution., (© 2024. The Author(s).)

Histone chaperones control nucleosome density and chromatin structure. In yeast, the H3-H4 chaperone Spt2 controls histone deposition at active genes but its roles in metazoan chromatin structure and organismal physiology are not known. Here we identify the Caenorhabditis elegans ortholog of SPT2 (CeSPT-2) and show that its ability to bind histones H3-H4 is important for germline development and transgenerational epigenetic gene silencing, and that spt-2 null mutants display signatures of a global stress response. Genome-wide profiling showed that CeSPT-2 binds to a range of highly expressed genes, and we find that spt-2 mutants have increased chromatin accessibility at a subset of these loci. We also show that SPT2 influences chromatin structure and controls the levels of soluble and chromatin-bound H3.3 in human cells. Our work reveals roles for SPT2 in controlling chromatin structure and function in Metazoa., (© 2024. The Author(s).)

Exploring the molecular basis of disease severity in rare disease scenarios is a challenging task provided the limitations on data availability. Causative genes have been described for Congenital Myasthenic Syndromes (CMS), a group of diverse minority neuromuscular junction (NMJ) disorders; yet a molecular explanation for the phenotypic severity differences remains unclear. Here, we present a workflow to explore the functional relationships between CMS causal genes and altered genes from each patient, based on multilayer network community detection analysis of complementary biomedical information provided by relevant data sources, namely protein-protein interactions, pathways and metabolomics. Our results show that CMS severity can be ascribed to the personalized impairment of extracellular matrix components and postsynaptic modulators of acetylcholine receptor (AChR) clustering. This work showcases how coupling multilayer network analysis with personalized -omics information provides molecular explanations to the varying severity of rare diseases; paving the way for sorting out similar cases in other rare diseases., (© 2024. The Author(s).)

Dipterocarpoideae species form the emergent layer of Asian rainforests. They are the indicator species for Asian rainforest distribution, but they are severely threatened. Here, to understand their adaptation and population decline, we assemble high-quality genomes of seven Dipterocarpoideae species including two autotetraploid species. We estimate the divergence time between Dipterocarpoideae and Malvaceae and within Dipterocarpoideae to be 108.2 (97.8‒118.2) and 88.4 (77.7‒102.9) million years ago, and we identify a whole genome duplication event preceding dipterocarp lineage diversification. We find several genes that showed a signature of selection, likely associated with the adaptation to Asian rainforests. By resequencing of two endangered species, we detect an expansion of effective population size after the last glacial period and a recent sharp decline coinciding with the history of local human activities. Our findings contribute to understanding the diversification and adaptation of dipterocarps and highlight anthropogenic disturbances as a major factor in their endangered status., (© 2024. The Author(s).)

The Eyes Absent proteins (EYA1-4) are a biochemically unique group of tyrosine phosphatases known to be tumour-promoting across a range of cancer types. To date, the targets of EYA phosphatase activity remain largely uncharacterised. Here, we identify Polo-like kinase 1 (PLK1) as an interactor and phosphatase substrate of EYA4 and EYA1, with pY445 on PLK1 being the primary target site. Dephosphorylation of pY445 in the G2 phase of the cell cycle is required for centrosome maturation, PLK1 localization to centrosomes, and polo-box domain (PBD) dependent interactions between PLK1 and PLK1-activation complexes. Molecular dynamics simulations support the rationale that pY445 confers a structural impairment to PBD-substrate interactions that is relieved by EYA-mediated dephosphorylation. Depletion of EYA4 or EYA1, or chemical inhibition of EYA phosphatase activity, dramatically reduces PLK1 activation, causing mitotic defects and cell death. Overall, we have characterized a phosphotyrosine signalling network governing PLK1 and mitosis., (© 2024. The Author(s).)

Cis-genetic effects are key determinants of transcriptional divergence in discrete tissues and cell types. However, how cis- and trans-effects act across continuous trajectories of cellular differentiation in vivo is poorly understood. Here, we quantify allele-specific expression during spermatogenic differentiation at single-cell resolution in an F1 hybrid mouse system, allowing for the comprehensive characterisation of cis- and trans-genetic effects, including their dynamics across cellular differentiation. Collectively, almost half of the genes subject to genetic regulation show evidence for dynamic cis-effects that vary during differentiation. Our system also allows us to robustly identify dynamic trans-effects, which are less pervasive than cis-effects. In aggregate, genetic effects were strongest in round spermatids, which parallels their increased transcriptional divergence we identified between species. Our approach provides a comprehensive quantification of the variability of genetic effects in vivo, and demonstrates a widely applicable strategy to dissect the impact of regulatory variants on gene regulation in dynamic systems., (© 2024. The Author(s).)