A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3

Pharmacogenetic pilot study of CYP2D6 and CYP1A2 genes in Italian patients with non-dystrophic myotonia and myotonic dystrophy treated with mexiletine

Mexiletine (MXT) is the first-line anti-myotonic drug for myotonic dystrophies (MD) and non-dystrophic myotonias (NDM), metabolized by CYP2D6 and CYP1A2 enzymes. We investigated genetic variants in these genes and their influence on MXT response in Italian MD and NDM patients. Fifty patients (MD: 34, NDM: 16) were treated with MXT (200-600 mg daily) for at least two months. Based on the Myotonic Behaviour Scale (MBS) and neurological examination, 37 patients (74 %) were classified as responders (R), while 13 (26 %) were non-responders (NR). Comparison between R and NR revealed associations with 14 genetic variants (12 in CYP2D6, 2 in CYP1A2). In silico analysis suggested eQTL effects on liver and skeletal muscle gene expression. Functional annotation indicated the regulatory roles of these variants. The CYP2D6*2/*41 diplotype showed a nominal association with non-response (OR = 10.8, p = 0.049), being 11 times more frequent in NR (23 %) than in R (3 %). Most common diplotypes (CYP2D6*1/*2, *1/*1, *1/*10, *2/10) corresponded to Normal Metabolizers, while CYP2D6*10/*10, *10/*41, and *41/41 indicated intermediate metabolism, suggesting a higher risk of adverse reactions with concurrent drugs. The CYP2D6*41 allele was more frequent in NR patients than in the European population, supporting its role in MXT response variability. Our findings suggest CYP2D6 and CYP1A2 variants as potential predictors of MXT treatment response in MD and NDM patients.

Comprehensive phenotyping combined with multi-omics of Salmonella Infantis and its H2S negative variant - Resolving adaption mechanisms to environmental changes

The zoonotic pathogen S. Infantis is of emerging importance, making detection in poultry critical. Phenotypic changes, which are significant for standardized control programs via EN/ISO 6579-1:2017, could lead to pathogens remaining undetected, increasing the risk of food-borne outbreaks. This study investigates an S. Infantis strain with both normal growth (NCP) and atypical H₂S-negative colony variant (ACV) from an Austrian broiler farm. NCP and ACV underwent comprehensive analyses, including stability tests, electron microscopy, whole-genome sequencing, transcriptomics, and proteomics. Our findings demonstrate a stable atypical colony variant exhibiting acquired resistance against cefoxitin in ACV. Genomic analysis identified 9 single nucleotide polymorphisms (SNPs) and two deletions, affecting genes involved in porphyrin and sulfur metabolism. Key factors were a mutation disrupting cysG, which is essential for siroheme biosynthesis and a vital cofactor in sulfur metabolism, and a stop codon in menD (2-oxoglutarate decarboxylase), crucial for small colony variant appearance. Consequently, we hypothesize that these mutations lead to a deficiency in siroheme, as well as anaerobic sulfur respiration altogether resulting in the H₂S-negative phenotype. Functional network analysis highlighted compensatory upregulation of alternative metabolic pathways, including nitrate metabolism, propanoate metabolism and mixed-acid fermentation, which may aid ACV's persistence and adaptation under anaerobic conditions. Reduced flagellin expression suggests a mechanism for immune evasion. These genetic and metabolic adaptations likely respond to environmental stressors, such as oxidative stress from disinfectants or antimicrobial pressure, leading to the emergence of the H₂S-negative phenotype. Consequently, this study provides insights into the genetic and biochemical adaptations of an atypical S. Infantis variant.

Identification of bioinformatic pipelines for virus monitoring using nanopore sequence data: A systematic assessment

Nanopore sequencing has proven to be a promising technique in virus surveillance efforts, especially due to the portability of its sequencers. In order to process the long, error-prone reads generated, specialised bioinformatic programs are required. These can be run automatically within pipelines so as to effectively provide decision makers with all relevant information about the molecular characteristics of a virus. The purpose of this systematic assessment was to identify pipelines that are suitable for virus surveillance programs using nanopore sequencing. Promising candidates were then compared in terms of their functional scope. Of 239 initial papers, 22 pipelines were tested, of which six were included in the final assessment. The four pipelines that were exclusively available offline were each missing individual downstream analysis steps considered in our assessment. The other two executed all steps. One of these was only available online and subject to a charge, while the other was freely available both online and offline. While we were able to identify two pipelines that are broadly suitable for virus surveillance using nanopore sequencing, we discovered two major shortcomings in this domain. None of the pipelines integrated basecalling, the initial step of data processing. In addition, there was no pipeline that was easy to install and provided all relevant analysis results with a single program call. We therefore see a need for the development of a pipeline that incorporates both aspects.

A chromosome-scale and haplotype-resolved genome assembly of tetraploid blackberry (Rubus L. subgenus Rubus Watson)

Blackberries (Rubus spp.) are globally consumed and well known for their rich anthocyanin and antioxidant content and distinct flavors. Improving blackberries has been challenging due to genetic complexity of traits and limited genomic resources. The blackberry genome has been particularly challenging to assemble due to its polyploid nature. Here, we present the first chromosome-scale and haplotype-phased assembly for the primocane-fruiting, thornless tetraploid blackberry selection BL1 (Rubus L. subgenus Rubus Watson). The genome assembly was generated using Oxford Nanopore Technology and Hi-C scaffolding, resulting in a 919 Mb genome distributed across 27 pseudochromosomes, with an N50 of 35.73 Mb. This assembly covers >92% of the genome length and contains over 98% of complete BUSCOs. Approximately, 58% of the assembly consists of repetitive sequences, with long terminal repeats being the most abundant class. A total of 87,968 protein-coding genes were predicted, of which, 82% were functionally annotated. Genome mining and RNA-Seq analyses identified possible candidate genes and transcription factors related to thornlessness and the key structural genes and transcription factors for anthocyanin biosynthesis. Activator genes including PAP1 and TTG1 and repressor genes such as ANL2 and MYBPA1 play an important role in the fine tuning of anthocyanin production during blackberry development. Resequencing of seven tetraploid blackberry cultivars/selections with different horticultural characteristics revealed candidate genes that could impact fruiting habit and disease resistance/susceptibility. This tetraploid reference genome should provide a valuable resource for accelerating genetic analysis of blackberries and facilitating the development of new improved cultivars with enhanced horticultural and nutritional traits.

Update on the nature of central giant cell granuloma of the jaw with a focus on the aggressive subtype

Central giant cell granuloma (CGCG) is a benign, localised osteolytic lesion of the jaw ​that is categorised into non-aggressive and aggressive subtypes. In contrast to non-aggressive CGCG, aggressive CGCG is characterised by pain, paraesthesia, root resorption, rapid growth, a size of >5 cm, cortical perforation, or recurrence after surgical treatment. However, the nature of CGCG, especially aggressive CGCG, remains unclear. This study was performed to analyse the systematic and comprehensive characteristics of CGCG of the jaw, especially the aggressive subtype, and first explored the genetic variation of aggressive CGCG by whole-exome sequencing. In total, 42 CGCGs were analysed (including 25 non-aggressive and 17 aggressive subtypes). H3F3A mutations were not detected in these CGCGs through immunohistochemistry and Sanger sequencing. The inability to detect H3F3A mutations could help differentiate CGCG from giant cell tumour of bone, indicating the two diseases are not different stages of the same pathological entity. Additionally, fluorescence in situ hybridisation did not reveal USP6 gene rearrangement in CGCG, which could distinguish it from aneurysmal bone cysts, especially the solid type. Therefore, H3F3A mutation and USP6 gene rearrangement detection have great significance in the clinicopathological diagnosis of CGCG of the jaw in terms of their ability to exclude giant cell tumour of bone and aneurysmal bone cyst. Moreover, the whole-exome sequencing data indicated that LRP1B gene abnormalities might be related to the aggressive biological behaviour of CGCG, and that NOTCH4 mutation could be a novel therapeutic target for aggressive CGCG.

Incorporation of Epstein-Barr viral variation implicates significance of Latent Membrane Protein 1 in survival prediction and prognostic subgrouping in Burkitt lymphoma

Although Epstein-Barr virus (EBV) plays a role in Burkitt lymphoma (BL) tumorigenesis, it is unclear if EBV genetic variation impacts clinical outcomes. From 130 publicly available whole-genome tumor sequences of EBV-positive BL patients, we used least absolute shrinkage and selection operator (LASSO) regression and Bayesian variable selection models within a Cox proportional hazards framework to select the top EBV variants, putative driver genes, and clinical features associated with patient survival time. These features were incorporated into survival prediction and prognostic subgrouping models. Our model yielded 22 EBV variants, including seven in latent membrane protein 1 (LMP1), as most associated with patient survival time. Using the top EBV variants, driver genes, and clinical features, we defined three prognostic subgroups that demonstrated differential survival rates, laying the foundation for incorporating EBV variants such as those in LMP1 as predictive biomarker candidates in future studies.

Across ancestries, HLA-B∗08:01∼DRB1∗03:01 (DR3) and HLA-DQA∗01:02 (DR2) increase the risk to develop juvenile-onset systemic lupus erythematosus through low complement C4 levels

Objective

Systemic lupus erythematosus (SLE) is a systemic autoimmune/inflammatory disease with a strong genetic component. Genetic burden is higher in children when compared to patients with adult-onset SLE, contributing to earlier disease expression and more severe phenotypes. The human leukocyte antigen (HLA) cluster on chromosome 6p21.3 is among the most variable genomic regions, representing a major risk-factor for SLE in adults. Its impact on juvenile-onset (j)SLE remains largely unstudied.

Methods

High-resolution sequencing of HLA class I (A, B, C), class II (DRB1, DQA1, DQB1) and class III (complement C2) was undertaken in the multi-ancestral UK JSLE Cohort including participants of Caucasian (n = 151, 48.8 %), Asian (n = 108, 35.0 %) and African/Caribbean (n = 50, 16.2 %) descent. Considering ancestral variation, clinical associations were tested at the level of alleles (2-field resolution), associated HLA protein sequences (antigen binding domains, 4-field resolution), and extended haplotypes (DRh).

Results

Although important ancestral recombination was reported for HLA-DR2 and -DR3 haplotypes, risk associated with jSLE was conserved at related alleles (DR2h: DRB1∗15:01, DQA∗01:02, DQB1∗06:02; DR3h: C∗07:02 [Asian], B∗08:01, C2 rs9332730 [Asian], DRB1∗03:01). HLA-DR7 haplotypes (DRB1∗07:01, OR = 0.44, 95 % CI:0.27-0.72, p = 0.0004; DQA1∗02:01, OR = 0.34, 95 % CI:0.21-0.56, p = 1.8 × 10-6) protect Asians from jSLE development. Among 23 clinical variables recorded, the main association was found between low levels of complement C4 in Caucasian carriers of HLA-DR3h. This was not the case in Asians due to recombination with HLA-C∗07:02 and integration of the C2 rs9332730 minor allele. Low C4 serum levels associated with HLA-DQA1∗01:02 (DR2h) in Caucasians after excluding HLA-DR3h carriers from the analysis. An association between low white blood cell counts and HLA-A∗03:01P was observed across ancestries.

Conclusion

Genetic variation in the HLA cluster associates with organ domain involvement (hematological) and complement levels in jSLE. Lupus-associated HLA haplotypes vary between ancestral groups, underscoring the importance of multi-ancestral approaches to genetic studies in SLE and other autoimmune/inflammatory diseases.

Genome-wide association study and BSR-seq identify nitrate reductase-related genes in rice landraces (Oryza sativa L.)

Nitrogen (N) is an essential nutrient for rice (Oryza sativa L.) growth and development. However, the lower nitrogen use efficiency (NUE) results in an N fertilizer surplus, which causes many environmental problems. In this study, genome-wide association studies were used to detect nitrate reductase (NR)-related loci in 419 rice landraces. Using the general linear model (GLM), mixed linear model (MLM), linear model (LM), and linear mixed model (LMM), we found six, nine, seven, and six significant single-nucleotide polymorphisms (SNPs) associated (p < 1 × 10-5) for three traits. Moreover, 98 significant SNPs were associated (logarithm of odds ≥ 3) with three traits through 3 V multi-locus random-SNP-effect mixed linear model. Interestingly, we found that Chr1_15896481 was significantly associated in the GLM, MLM, LM, and LMM models. Meanwhile, this significant locus overlapped with a candidate region in bulked segregant RNA sequencing. Through integrated analysis, we identified a most likely candidate genomic region 15,627,420-16,084,761 bp on chromosome 1. By performing functional annotation, RNA sequencing, and real-time quantitative polymerase chain reaction (RT-qPCR) analysis for the genes within this interval, we identified five candidate genes that may affect NR activity. Os01g0378400 exhibits a gene expression pattern highly similar to that of OsNR1.2. It belongs to the NAC transcription factor family, which is involved in plant N metabolism. Os01g0377700 is homologous to an ammonium transporter gene (Cre06g293051). Os01g0383700 encodes a WD40 domain protein, Os01g0379400 encodes an F-box protein, and Os01g0382800 encodes a DYW-type PPR domain protein. These findings will provide valuable genetic resources for NUE genetic improvement in rice breeding.

Whole genome sequencing data of greater burdock (Arctium lappa) naturalized in the United States of America

This dataset comprises whole genome sequencing burdock (Arctium lappa) naturalized in a residential yard in Rochester, New York, USA. Total DNA was extracted from a leaf sample and processed using the Illumina Nextera XT DNA library preparation kit. Sequencing on the NextSeq 2000 platform produced 127.4 GB of raw data, yielding 125.8 GB of high-quality reads after filtering, with an average genome coverage of 75x. The genome was assembled de novo into 792,817 contigs, achieving a total genome length of 1,075,454,921 base pairs with a GC content of 37.03 %. Scaffolding against a Chinese A. lappa reference genome improved genome completeness from 49.1 % to 94.93 %, successfully recovering the majority of protein-coding genes. Variant analysis identified approximately 20.8 million Single Nucleotide Polymorphisms (SNPs) and 1.3 million indels, including functionally significant mutations. The Internal Transcribed Spacer 2 (ITS2) ribosomal region was isolated and compared with global references, revealing significant genetic differentiation between the U.S.A and Chinese populations. This comprehensive genomic dataset has been deposited in publicly accessible repositories, including National Center for Biotechnology Information (NC and Zenodo. The sequencing of this sample provides a valuable resource for comparative genomics, population genetics, and investigations into bioactive compounds with antimicrobial properties, supporting agricultural and pharmaceutical applications. Direct access to the dataset is available at 10.5281/zenodo.14607136.

Differences in allele frequencies of personality-related genes in three varieties of Shiba Inu in Japan

Shiba Inu is a Japanese dog breed that originally functioned as a hunting dog but is now widely kept as a pet. This study conducted whole genome sequencing on two varieties of Shiba Inu, Mino-Shiba and San'in-Shiba, which are believed to strongly retain the hunting function, and compared them with the common pet Shiba Inu. The results showed that populations of the three varieties formed distinct clusters. We calculated the FST values for each site between each of the two varieties and identified genomic regions with significant differences that might affect function. LRRTM4 and OXTR found in previous studies on dogs related to personality traits are indicated in this study, and KIF27, associated with Williams syndrome in humans, was also suggested to potentially influence hunting traits. This study highlights the unique genetic lineage of regional Shiba varieties compared to the commonly kept pet Shiba Inu, and provides a foundation for further research into how these genetic differences may affect current personality traits.

A FAM8A1 frameshift variant is associated with REM sleep behavior disorder, urinary retention, and mydriasis in Russian Blue cats

REM sleep behavior disorder (RBD) is a disease characterized by the loss of lower motor neuron inhibition responsible for skeletal muscle atonia during REM sleep. It has been reported in humans, dogs and cats, and can be idiopathic or secondary to a neurodegenerative disease. Five young adult Russian Blue cats from two related families were presented for progressively worsening RBD episodes frequently associated with urinary loss. Three of these cats also suffered urinary retention with overflow incontinence between RBD episodes. Neurological examination revealed a large bladder in three cats and a bilateral mydriasis with absent pupillary light reflexes in two cats; further examinations were unremarkable. Treatment attempts were unsatisfactory, with four cats being euthanized. Histopathology of the brain did not reveal any abnormalities. A disease-associated 23-bp deletion in exon 1 of FAM8A1 (NC_058372.1:g.11622168_11622190del), introducing a frameshift at codon 162 and a premature stop codon at codon 276 (XM_019831563.3:c.485_507del p.(Gln162Profs*115)), was identified by whole genome sequencing. The variant segregated in the affected families with a recessive mode of inheritance, showed an allele frequency of 1.5% in West-European Russian Blue cats (N = 68) and was not present in 276 cats belonging to 32 other breeds (including the closely related Nebelung breed). The variant FAM8A1 isoform is predicted to affect the assembly and activity of the endoplasmic reticulum-associated protein degradation pathway, which plays an important role in cell homeostasis. RBD and urinary retention syndrome is a hereditary encephalopathy affecting Russian Blue cats. A genetic test now allows diagnosis and prevention of this debilitating disease.

Whole-genome sequences revealed genomic diversity and selection signatures of Dermacentor silvarum in Shaanxi, China

Dermacentor silvarum (D. silvarum) is an arthropod that feeds on blood. It is a primary tick species found in northern China that poses a significant security risk to the health and life of the host, as it has the potential to transmit a variety of pathogens to humans and animals. Through ongoing research on tick genome sequences, researchers have successfully assembled and reported reference genomes for numerous tick species. These significant advances have greatly accelerated the study of tick biology and population genomics. D. silvarum samples were obtained from the body surface of free-range goats in Yulin, Shaanxi Province, China. The whole genomes of the samples were resequenced and merged with preexisting data from the National Genomics Data Center database (project ID: PRJCA002242) to analyze the genetic structure, genetic diversity, mitochondrial genetic structure, and selection signatures of D. silvarum in the Shaanxi Province. Based on the available data, the D.silvarum species in China could be classified into two main branches. These populations exhibited low nucleotide diversity. A slight discrepancy was noted between the mitochondrial phylogenetic tree and the autosomal whole-genome phylogenetic tree of D. silvarum, consistent with a previous study. In the selected analysis of D. silvarum in the Shaanxi Province, China, genes linked to immunity, iron storage, fatty acid biosynthesis, pesticide defense, and blood digestion were identified. Leutriene A4 hydrolase (LOC119466376) was also identified, although its function remains unknown. This information is crucial for understanding the biology of D. silvarum and developing management measures.

A Case-Driven Multi-Omics Analysis for Longitudinal Ibrutinib Response Evaluation of Patients With Chronic Lymphocytic Leukemia

Patients with chronic lymphocytic leukemia (CLL) undergoing ibrutinib treatment often experience incomplete response, yet the molecular level underlying clonal inertia remains to be explored. We investigated the molecular and clinical dynamics of CLL during 16 months of ibrutinib monotherapy by analyzing blood samples from two patients who continued having CLL cells in the peripheral blood during treatment. At diagnosis, the clonal burden within the B cell compartment was found to be 55% (pt1) and 86% (pt2) for the dominant clones. At 16 months following treatment these clones still constituted 66% and 89%, respectively. Utilizing multi-omic methodologies at the DNA and RNA levels, including single-cell transcriptomics, we aimed to establish a comprehensive framework for multi-omics analysis for longitudinal ibrutinib response evaluation. The presented study revealed genomically stable disease during ibrutinib treatment, but with intensified expression of genes involved in pathways related to apoptosis, cellular stress response, and canonical NF-κB signaling from diagnosis to 16 months of treatment.

Novel genetic associations with childhood adipocytokines in Indian adolescents

Adipocytokines, including leptin, adiponectin, and resistin, are key mediators linking adiposity, insulin resistance, and inflammation. We present the first genome-wide association study (GWAS; N = 5258) and exome-wide association study (ExWAS; N = 4578) on leptin, adiponectin, and resistin in South Asian population. We identified novel associations in genes ZNF467, and LEPREL2 for leptin; ZNF467, LEPREL2, CRLF3, ZNF732, SOX30, XIRP1, ATP8B3, SPATA2L, TMCO4, TLN2, ABCA12, and SHB for adiponectin; and D2HGDH for resistin. Additionally, we confirmed known associations of FTO, MC4R, and HOXB3 with leptin and ADIPOQ with adiponectin. Notably, ADIPOQ variants were consistently significant across GWAS, ExWAS, and gene-based analyses, reinforcing their central role in regulating adiponectin levels. Most of these novel associations identified were population-specific, highlighting the importance of studying diverse populations to uncover unique genetic signals. After adjusting for BMI, the associations with adiponectin and resistin remained significant, whereas most associations for leptin weakened in both effect size and significance. Functional annotation revealed that the identified variants were enriched for expression in adipose tissue, the brain (cerebellar hemisphere and cerebral cortex), and the pituitary gland. These variants act as eQTLs and splice-QTLs in adipose, brain, and pancreas, suggesting cross-tissue regulatory mechanisms. ExWAS further implicated rare variant burden in genes such as LONP1, ZNF335, and TTC16 for adiponectin and resistin. These findings enhance our understanding of adipocytokine biology, emphasises the need for population-specific genetic research, and lays foundation for future functional studies.

The white lupin CCR1 receptor-like kinase controls systemic Autoregulation of Cluster Root and Nodule Development

Root development is tightly regulated in plants to optimize nutrient acquisition and interactions with soil microorganisms. In legumes, the Autoregulation of Nodulation (AoN) pathway systemically controls the proliferation of root nodules, which are energy-intensive organs. Mutations affecting the AoN pathway result in a hypernodulation phenotype accompanied by altered root development. However, it remains unclear whether this modification of root development is also systemic and coordinated with nodulation. In this study, we report the identification of the constitutive cluster root 1 (ccr1) mutant in white lupin (Lupinus albus), which exhibits constitutive production of an excessive number of cluster roots. We demonstrate that CCR1 is an ortholog of HAR1/SUNN/NARK leucin-rich repeat-receptor like kinases (LRR-RLKs), which are key regulators of the AoN pathway. Furthermore, we show that CCR1 negatively regulates both nodule and cluster root development. Interspecific grafting experiments between white and narrow-leaved lupin (Lupinus angustifolius), a species incapable of producing cluster roots, show that ccr1 shoots can induce the formation of cluster-like roots in narrow-leaved lupin rootstocks. This highlights the conservation of a CCR1-dependent signaling cascade. Transcriptome analyses reveal that CCR1 targets the conserved NIN/LBD16-NFYA regulatory module, which connects nodule and lateral root development through a shared inhibitory systemic pathway. We propose that this pathway represents a broader developmental control mechanism of root organogenesis, termed Autoregulation of Cluster Root and Nodule Development (AoDev).

Partner dependency alters patterns of coevolutionary selection in mutualisms

Coevolution is a ubiquitous driver of diversification in both mutualistic and antagonistic interactions between species. In mutualisms, coevolution can result in trait complementarity between partners that facilitates their persistence. Despite its importance, most of what we know about coevolution in mutualism comes from obligate interactions, whereas we know comparatively little about facultative interactions, arguably the most common type of mutualism. To evaluate coevolutionary dynamics in facultative mutualism and test how it compares with obligate mutualisms, we used a synthetic yeast mutualism where the partners exchange essential nutrient resources. We manipulated mutualism dependency by controlling the availability of mutualistic resources in the environment and measured coevolution via time-shift assays and tracking the evolution of mutualistic traits over time. In addition, we genotyped the evolved and ancestral mutualists to test for differences in the strength of coevolutionary selection between facultative and obligate mutualisms. We found evidence of coevolution in both facultative and obligate mutualisms, but coevolution was weaker and slower in facultative mutualisms. We also found evidence for evolution of trait complementarity in obligate mutualisms but not in facultative mutualisms. Furthermore, obligate mutualists had more SNPs under positive selection than facultative mutualists. Together, these results provide strong evidence that mutualism dependency impacts both the strength of coevolution and the rate of trait evolution.

EFNB3 Frameshift Variant in Weimaraner Dogs with a Condition Resembling a Congenital Mirror Movement Disorder

BACKGROUND

Congenital mirror movement disorders (CMMs) are clinically and genetically heterogeneous in human patients. CMMs have not been documented to occur spontaneously in animals.

OBJECTIVE

The objective of this work was to document the first case of CMMs spontaneously occurring in Weimaraner dogs and to identify the underlying genetic cause.

METHODS

Clinical and pathological investigations were performed. Genetic investigations used linkage and autozygosity mapping followed by whole-genome sequencing of 3 affected dogs and 1489 control dogs to identify disease-associated variants.

RESULTS

Three of 11 puppies in a litter of Weimaraner dogs exhibited an abnormal gait characterized by synchronized saltatorial locomotion. Their phenotype was tentatively termed congenital mirror movement disorder 1 (CMM1). The underlying genetic cause was identified as a 2-bp duplication in EFNB3 encoding ephrin-B3, a transmembrane protein important for axon guidance and spinal midline barrier formation during neurodevelopment. The identified variant, XM_038536724.1:c.643_644dup, is predicted to lead to a frameshift and introduction of a premature stop codon XP_038392652.1:p.(Ala216Valfs*79). CMM1 is inherited as an autosomal recessive trait in these dogs.

CONCLUSIONS

Similar to humans, CMMs may occur in dogs as an inherited disease as a result of a spontaneously arisen genetic variant. The CMM1 phenotype in dogs resembles the phenotype of experimentally induced Efnb3-/- knockout mice. So far, no human patients with EFNB3-related CMMs have been reported. Our study provides the first naturally occurring large-animal model for CMMs. EFNB3 should be considered a candidate gene in human CMM patients with unclear disease etiology. © 2025 International Parkinson and Movement Disorder Society.

Tumor antigens preferentially derive from unmutated genomic sequences in melanoma and non-small cell lung cancer

Melanoma and non-small cell lung cancer (NSCLC) display exceptionally high mutational burdens. Hence, immune targeting in these cancers has primarily focused on tumor antigens (TAs) predicted to derive from nonsynonymous mutations. Using comprehensive proteogenomic analyses, we identified 589 TAs in cutaneous melanoma (n = 505) and NSCLC (n = 90). Of these, only 1% were derived from mutated sequences, which was explained by a low RNA expression of most nonsynonymous mutations and their localization outside genomic regions proficient for major histocompatibility complex (MHC) class I-associated peptide generation. By contrast, 99% of TAs originated from unmutated genomic sequences specific to cancer (aberrantly expressed tumor-specific antigens (aeTSAs), n = 220), overexpressed in cancer (tumor-associated antigens (TAAs), n = 165) or specific to the cell lineage of origin (lineage-specific antigens (LSAs), n = 198). Expression of aeTSAs was epigenetically regulated, and most were encoded by noncanonical genomic sequences. aeTSAs were shared among tumor samples, were immunogenic and could contribute to the response to immune checkpoint blockade observed in previous studies, supporting their immune targeting across cancers.

Genome-Wide Association Study of Age-Related Hearing Loss in CFW Mice Identifies Multiple Genes and Loci, Including Prkag2

PURPOSE

Age-related hearing loss (ARHL) is one of the most prevalent conditions affecting the elderly. ARHL is influenced by a combination of environmental and genetic factors; the identification of the genes that confer risk will aid in the prevention and treatment of ARHL. The mouse and human inner ears are functionally and genetically homologous. We used Carworth Farms White (CFW) mice to study the genetic basis of ARHL because they are genetically diverse and exhibit variability in the age of onset and severity of ARHL.

METHODS

Hearing at a range of frequencies was measured using auditory brainstem response (ABR) thresholds in 946 male and female CFW mice at the age of 1, 6, and 10 months. We genotyped the mice using low-coverage (mean coverage 0.27 ×) whole-genome sequencing (lcWGS) followed by imputation using STITCH. To determine the accuracy of the genotypes, we sequenced 8 samples at > 30 × coverage and used those data to estimate the accuracy of lcWGS genotyping, which was > 99.5%. We performed a genome-wide association study (GWAS) for the ABR thresholds for each frequency at each age, and we also performed a GWAS for age at deafness.

RESULTS

We obtained genotypes at 4.18 million single nucleotide polymorphisms (SNP). The SNP heritability for traits ranged from 0 to 42%. GWAS identified 10 significant associations with ARHL that contained potential candidate genes, including Dnah11, Rapgef5, Cpne4, Prkag2, and Nek11. Genetic ablation of Prkag2 caused ARHL at high frequencies, strongly suggesting that Prkag2 is the causal gene for one of the associations.

CONCLUSIONS

GWAS for ARHL in CFW outbred mice identified genetic risk factors for ARHL, including Prkag2. Our results will help to define novel therapeutic targets for the treatment and prevention of this common disorder.

Reveal genomic insights into cotton domestication and improvement using gene level functional haplotype-based GWAS

Genome-wide association studies (GWAS) are widely used to detect associations between genetic variants and phenotypes. However, few studies have thoroughly analyzed genes, the fundamental and most crucial functional units. Here, we develop an innovative strategy to translate genomic variants into gene-level functional haplotypes (FHs), effectively reducing the interference from complex genome structure and linkage disequilibrium (LD) present in the conventional genetic mapping framework. Using refined mixed linear models, gene-level FH is regressed with 20 cotton agronomic traits across 245 sets of phenotypic values in 3,724 accessions, directly identifying 532 quantitative trait genes (QTGs) with significant breeding potential. The biological function of a superior fiber quality QTG encoding ferulic acid 5-hydroxylase 1 is experimentally validated. Thereafter, we systematically analyze the genetic basis of cotton domestication and improvement at the gene level. This report provides genomic insight into the genetic dissection and efficient mapping of functional genes in plants.

Pathogenicity and transmissibility of bovine-derived HPAI H5N1 B3.13 virus in pigs

Since the first emergence of highly pathogenic avian influenza (HPAI) H5N1 viruses in dairy cattle, the virus has continued to spread, reaching at least 17 states and at least 950 dairy herds in the United States. Subsequently, spillovers of the virus from dairy cattle to humans have been reported. Pigs are an important reservoir in influenza ecology because they serve as a mixing vessel in which novel reassortant viruses with pandemic potential can be generated. Here, we show that oro-respiratory infection of pigs resulted in productive replication of a bovine-derived HPAI H5N1 B3.13 virus. Infectious virus was mainly identified in the lower respiratory tract of principal infected pigs, and sero-conversion was observed in most of the principal pigs at later time points, suggesting limited replication of the bovine-derived HPAI H5N1 B3.13 virus in pigs. In one animal, we detected the emergence of a mutation in hemagglutinin (HA) previously associated with increased affinity for "mammalian-type" α2,6-linked sialic acid receptors, but this mutation did not reach majority consensus levels. Sentinel contact pigs remained sero-negative throughout the study, indicating lack of transmission. These results support that pigs are susceptible to a bovine-derived HPAI H5N1 B3.13 virus, but this virus did not replicate as robustly in pigs as mink-derived HPAI H5N1 and swine-adapted influenza viruses.

Organoid morphology-guided classification for oral cancer reveals prognosis

Oral cancer is an aggressive malignancy with a survival rate below 50% in advanced stages due to low mutation rates, lack of molecular subtypes, and limited treatment targets. This study presents a pioneering approach to classifying oral cancer subtypes based on the morphology of patient-derived organoids (PDOs) and proposes a therapeutic strategy. We establish 76 cancer and 81 normal PDOs. For cancer PDOs, both manual classification and AI-based scoring are utilized to categorize them into three distinct subtypes: normal-like, dense, and grape-like. These subtypes correlate with unique transcriptomic profiles, genetic mutations, and clinical outcomes, with patients harboring dense and grape-like organoids exhibiting poorer prognoses. Furthermore, drug response assessments of 14 single agents and cisplatin combination therapies identify a synergistic treatment approach for resistant subtypes. This study highlights the potential of integrating morphology-based classification with genomic and transcriptomic analyses to refine oral cancer subtyping and develop effective treatment strategies.

ROBIN: A unified nanopore-based assay integrating intraoperative methylome classification and next-day comprehensive profiling for ultra-rapid tumor diagnosis

BACKGROUND

Advances in our technological capacity to interrogate CNS tumor biology have led to the ever increasing use of genomic sequencing in diagnostic decision making. Presently, CNS tumors are classified based on their epigenetic signatures, leading to a paradigm shift in diagnostic pathways. Such testing can be performed so rapidly using nanopore sequencing that results can be provided intraoperatively. This information greatly improves the fidelity of smear diagnosis and can help surgeons tailor their approach, balancing the risks of surgery with the likely benefit. Nevertheless, full integrated diagnosis may require subsequent additional assays to detect pathognomonic somatic mutations and structural variants, thereby delaying the time to final diagnosis.

METHODS

Here, we present ROBIN, a tool based on PromethION nanopore sequencing technology that can provide both real-time, intraoperative methylome classification and next-day comprehensive molecular profiling within a single assay. ROBIN utilizes 3 methylation classifiers to improve diagnostic performance in the intraoperative setting.

RESULTS

We demonstrate classifier performance on 50 prospective intraoperative cases, achieving a diagnostic turnaround time under 2 hours and generating robust tumor classifications within minutes of sequencing. Furthermore, ROBIN can detect single nucleotide variants, copy number variants, and structural variants in real time, and is able to inform a complete integrated diagnosis within 24 hours. Classifier performance demonstrated concordance with final integrated diagnosis in 90% of prospective cases.

CONCLUSION

Nanopore sequencing can greatly improve turnaround times for standard-of-care diagnostic testing and is furthermore able to reliably provide clinically actionable intraoperative tumor classification.

Low T cell diversity associates with poor outcome in bladder cancer: A comprehensive longitudinal analysis of the T cell receptor repertoire

T cells are crucial effector cells in the endogenous defense against cancer, yet the clinical impact of their quantity, diversity, and dynamics remains underexplored. Here, we investigate the clinical relevance of the T cell receptor (TCR) repertoire in patients with bladder cancer. In advanced-stage disease, low pre-treatment peripheral TCR diversity is associated with worse overall survival (p = 0.024), particularly when coupled with low circulating T cell fractions (p = 0.00049). These low-diversity repertoires are dominated by hyper-expanded clones that persist throughout treatment. Further longitudinal analysis reveals reductions in TCR diversity after treatment, indicating adverse effects on the immune system. In early-stage disease, immunotherapy increases TCR diversity in patients with good outcomes. Furthermore, single-cell sequencing identifies most hyper-expanded clones as cytotoxic T cells, while non-expanded clones are predominantly naive T cells. Overall, this highlights TCR diversity as a promising biomarker, offering opportunities for tailored oncological treatments to enhance clinical outcomes.

SKSR1 identified as key virulence factor in Cryptosporidium by genetic crossing

Cryptosporidium is a major cause of severe diarrhea. Although Cryptosporidium isolates exhibit significant differences in infectivity and virulence, the genetic determinants for these traits are not clear. In this study, we use classical genetics to cross two Cryptosporidium parvum isolates of different virulence and use bulk segregant analysis of whole-genome sequences from the progeny to identify quantitative trait loci (QTL) associated with Cryptosporidium infectivity and virulence. Of the 23 genes in three QTL, two have loss-of-function mutations in the low-virulence isolates, including the SKSR1 gene encoding a variant secretory protein. Deletion of the SKSR1 gene or expression of the frame-shifted sequence reduces the pathogenicity of the virulent isolate. SKSR1 is expressed in small granules and secreted into the parasite-host interface during invasion. These results demonstrate that SKSR1 is an important virulence factor in Cryptosporidium, and suggest that the extended SKSR protein family, encoded by clusters of subtelomeric genes, may contribute to pathogenesis.

Sexual selection drives the speciation of lineages with contrasting mating systems

Repeated shifts in ecological strategies often lead to consistent speciation patterns across phylogenies. One example is the transition from outcrossing to self-fertilization in plants, which generally results in the reproductive isolation of the incipient selfing lineages. However, the evolutionary mechanisms driving their speciation remain poorly understood. In this study, we investigate the hybridization rate and barriers to gene flow between the recently evolved selfing lineage Capsella rubella and its outcrossing ancestor C. grandiflora. Through a survey of sympatric populations in Greece, we found that despite coexisting in the same habitats, the two species rarely form viable hybrids. Our investigation into the mechanisms underlying this reproductive isolation revealed that differences in the intensity of sexual selection between the lineages promote significant prezygotic isolation, with the strength of this isolation depending on the direction of gene flow. Traits that enhance male competitiveness in outcrossers decrease their chance of being pollinated by selfers, lowering the hybridization rate, but simultaneously increase the likelihood of selfers being pollinated by outcrossers. Selfers nevertheless limit hybridization through rapid and efficient self-fertilization mechanisms. Finally, postzygotic barriers, such as hybrid incompatibilities likely driven by differences in parental conflict intensity,1 also contribute to the isolation of the two lineages. Therefore, shifts in the intensity of sexual selection and increase in selfing efficiency appear to be key drivers of reproductive isolation following mating system changes, possibly explaining recurrent speciation patterns in plant evolution.

The knockout of ClaCSLH1 induced dwarfing in watermelon

In agriculture, selecting ideal plant types with desirable traits, such as dwarfing and upright stem structures, significantly enhances crop yield and quality by optimizing light absorption, spatial efficiency, and nutrient utilization. Developing new varieties of dwarf watermelon is a crucial objective in watermelon breeding. In this study, we constructed an F2 population using the wild-type V063 as the paternal parent and the dwarf variety dw-n as the maternal parent. The dwarfing trait was found to be governed by a pair of recessive alleles. Through bulk segregant analysis sequencing (BSA-seq) and RNA sequencing (RNA-seq), we identified the gene Cla97C02G035450, which encodes cellulose synthase-like H1 (CSLH1), as a candidate gene associated with the dwarfing phenotype. ClaCLSH1 belongs to the ClaCESA/CSLs family, which is involved in the cell wall formation by regulating the synthesis of cellulose and hemicellulose. Microscopic analyses revealed that dw-n exhibited shorter internode cells, thicker cell walls, and elevated hemicellulose content compared to V063. Subcellular localization studies demonstrated that the CLACSLH1 protein is primarily localized in the nucleus and the cell membrane/wall. Notably, the overexpression of CLACSLH1 in the dw-n background rescued its dwarf phenotype. Furthermore, experiments indicated that knockdown of CLACSLH1 resulted in excessive hemicellulose synthesis, inhibited internode cell elongation, and ultimately led to the stunted phenotype observed in dw-n. This research provides innovative insights into the development of superior dwarf watermelon varieties and advances our understanding of the molecular mechanisms underlying watermelon dwarfism.

Understanding the epidemiology and pathogenesis of Mycobacterium tuberculosis with non-redundant pangenome of epidemic strains in China

Tuberculosis is a major public health threat resulting in more than one million lives lost every year. Many challenges exist to defeat this deadly infectious disease which address the importance of a thorough understanding of the biology of the causative agent Mycobacterium tuberculosis (MTB). We generated a non-redundant pangenome of 420 epidemic MTB strains from China including 344 Lineage 2 strains, 69 Lineage 4 strains, six Lineage 3 strains, and one Lineage 1 strain. We estimate that MTB strains have a pangenome of 4,278 genes encoding 4,183 proteins, of which 3,438 are core genes. However, due to 99,694 interruptions in 2,447 coding genes, we can only confidently confirm 1,651 of these genes are translated in all samples. Of these interruptions, 67,315 (67.52%) could be classified by various genetic variations detected by currently available tools, and more than half of them are due to structural variations, mostly small indels. Assuming a proportion of these interruptions are artifacts, the number of active core genes would still be much lower than 3,438. We further described differential evolutionary patterns of genes under the influences of selective pressure, population structure and purifying selection. While selective pressure is ubiquitous among these coding genes, evolutionary adaptations are concentrated in 1,310 genes. Genes involved in cell wall biogenesis are under the strongest selective pressure, while the biological process of disruption of host organelles indicates the direction of the most intensive positive selection. This study provides a comprehensive view on the genetic diversity and evolutionary patterns of coding genes in MTB which may deepen our understanding of its epidemiology and pathogenicity.

Whole-exome sequencing identifies 5 novel genes associated with carpal tunnel syndrome

Carpal tunnel syndrome (CTS), a common peripheral nerve entrapment disorder, has a high estimated heritability index. Although previous genome-wide association studies have assessed common genetic components of CTS, the risk contributed by coding variants is still not well understood. Here, we performed the largest exome-wide analyses using UK Biobank data from 350 770 participants to find coding variants associated with CTS. We then explored the relative contribution of both rare mutations and polygenic risk score (PRS) to CTS risk in survival analyses. Finally, we investigated the functional pathways of the CTS-related coding genes identified above. Aside from conforming 6 known CTS genes, 5 novel genes were identified (SPSB1, SYNC, ITGB5, MUC13 and LOXL4). The associations of most genes we identified with incident CTS were striking in survival analyses. Additionally, we provided evidence that combining rare coding alleles and polygenic risk score can improve the genetic prediction of CTS. Functional enrichment analyses revealed potential roles of the identified coding variants in CTS pathogenesis, where they contributed to extracellular matrix organization. Our results evaluated the contribution to CTS etiology from quantities of coding variants accessible to exome sequencing data.

Population Genomics Reveals Demographic History and Climate Adaptation in Japanese Arabidopsis halleri

Climate oscillations in the Quaternary forced species to major latitudinal or altitudinal range shifts. It has been suggested that adaptation concomitant with range shifts plays key roles in species responses during climate oscillations, but the role of selection for local adaptation to climatic changes remains largely unexplored. Here, we investigated population structure, demographic history and signatures of climate-driven selection based on genome-wide polymorphism data of 141 Japanese Arabidopsis halleri individuals, with European ones as outgroups. Coalescent-based analyses suggested a genetic differentiation between Japanese subpopulations since the Last Glacial Period (LGP), which would have contributed to shaping the current pattern of population structure. Population demographic analysis revealed the population size fluctuations in the LGP, which were particularly prominent since the subpopulations started to diverge (∼50, 000 years ago). The ecological niche modeling predicted the geographic or distribution range shifts from southern coastal regions to northern coastal and mountainous areas, possibly in association with the population size fluctuations. Through genome-wide association analyses of bioclimatic variables and selection scans, we investigated whether climate-associated loci are enriched in the extreme tails of selection scans, and demonstrated the prevailing signatures of selection, particularly toward a warmer climate in southern subpopulations and a drier environment in northern subpopulations, which may have taken place during or after the LGP. Our study highlights the importance of integrating climate associations, selection scans and population demographic analyses for identifying genomic signatures of population-specific adaptation, which would also help us predict the evolutionary responses to future climate changes.

Onkopus: precise interpretation and prioritization of sequence variants for biomedical research and precision medicine

One of the major challenges in precision oncology is the identification of pathogenic, actionable variants and the selection of personalized treatments. We present Onkopus, a variant interpretation framework based on a modular architecture, for interpreting and prioritizing genetic alterations in cancer patients. A multitude of tools and databases are integrated into Onkopus to provide a comprehensive overview about the consequences of a variant, each with its own semantic, including pathogenicity predictions, allele frequency, biochemical and protein features, and therapeutic options. We present the characteristics of variants and personalized therapies in a clear and concise form, supported by interactive plots. To support the interpretation of variants of unknown significance (VUS), we present a protein analysis based on protein structures, which allows variants to be analyzed within the context of the entire protein, thereby serving as a starting point for understanding the underlying causes of variant pathogenicity. Onkopus has the potential to significantly enhance variant interpretation and the selection of actionable variants for identifying new targets, drug screens, drug testing using organoids, or personalized treatments in molecular tumor boards. We provide a free public instance of Onkopus at https://mtb.bioinf.med.uni-goettingen.de/onkopus.

Inter-host diversity associated with age, sex, and menstrual cycle modulates clinical manifestations in DENV-2 patients

Dengue virus (DENV-2) remains a global threat, yet the influence of age, sex, and menstrual status on its epidemiology and genetic diversity is underexplored. We analyzed 2136 hospitalized DENV-2 patients (ages 0-86) using whole-genome sequencing (WGS) to examine how these host factors shape interhost viral diversity and clinical manifestations. Young adult males (19-35 years) had the highest prevalence with sex-based clinical differences where females exhibited severe hematological changes, while males showed increased hepatic injury. Premenopausal females harbored more diverse viral populations, whereas postmenopausal women experienced pronounced platelet depletion. Dengue virus WGS identified 1100336 single nucleotide variants (SNVs) across 2932 genomic positions, with greater viral diversity in adults and females. Significant SNV burdens were observed in the E, NS3, and NS5 genes of the virus. These interconnected findings underscore the profound impact of age, sex, and menstrual status on DENV-2 epidemiology which merits inclusion into the disease pathophysiology.

The whole genome analysis of the wild-type and attenuated orf virus reveals that ORF022 facilitates viral replication

BACKGROUND

Contagious ecthyma is an acute infectious zoonosis caused by orf virus (ORFV). Live-attenuated ORFV vaccines have played a crucial role in preventing contagious ecthyma for decades. However, these vaccines often fail to induce long-lasting immunity. In recent years, numerous ORFV genome sequences have been published, yet genomic data for attenuated strains remain limited. Furthermore, no comprehensive whole-genome-based single nucleotide polymorphisms (SNPs) analysis has been conducted to compare ORFV wild-type and attenuated strains.

RESULTS

In this study, we performed whole-genome sequencing of ORFV wild-type and attenuated strains from Shaanxi Province. We identified two ORFV strains with genomes shorter than 130 kb, which are closely related to the SC1 attenuated strain from Sichuan Province. Additionally, we noticed that 24 genes in the attenuated strain had SNPs, with the highest number of mutations occurring in the ORF022 gene. The function of the ORF022 gene has not been previously reported. Through in vitro experiments, we demonstrated that overexpression of ORF022 enhances ORFV replication in cells. The RNA-sequencing analysis revealed that ORF022 modulates host inflammation-related signaling pathways, as evidenced by the suppression of TNF, IL-17, and Toll-like receptor signaling pathways.

CONCLUSIONS

Our findings suggest that the ORF022 in ORFV wild-type strain inhibits the host inflammatory response, reduces the immune response to ORFV, and facilitates viral replication. SNP events in attenuated strains (aFX0910) are one of the reasons for its attenuation. Investigations into the genomic sequences of attenuated viruses and the functional impact of mutated genes provide valuable insights into the mechanisms underlying ORFV attenuation and offer a foundation for the development of more effective ORFV vaccines.

Single-Cell RNA Sequencing of Ewing Sarcoma Tumors Demonstrates Transcriptional Heterogeneity and Clonal Evolution

PURPOSE

Ewing sarcoma is the second most common bone cancer in children, accounting for 2% of pediatric cancer diagnoses. Patients who present with metastatic disease at the time of diagnosis have a dismal prognosis compared with the >70% 5-year survival of those with localized disease. Novel therapeutic approaches that can impact metastatic disease are desperately needed, as well as a deeper understanding of the heterogeneity of Ewing sarcoma tumors.

EXPERIMENTAL DESIGN

In this study, we utilized single-cell RNA sequencing to characterize the transcriptional landscape of primary Ewing sarcoma tumors and the surrounding tumor microenvironment in a cohort of seven untreated patients with Ewing sarcoma, as well as in circulating tumor cells (CTC). A potential CTC therapeutic target was evaluated through immunofluorescence of fixed CTCs from a separate cohort.

RESULTS

Primary tumor samples demonstrate a heterogeneous transcriptional landscape with several conserved gene expression programs, including those composed of genes related to proliferation and Ewing sarcoma gene targets, which were found to correlate with overall survival. Copy-number analysis identified subclonal evolution within patients prior to treatment. Analyses of the immune microenvironment reveal an immunosuppressive microenvironment with complex intercellular communication among the tumor and immune cells. Single-cell RNA sequencing and immunofluorescence of CTCs at the time of diagnosis identified TSPAN8 as a potential therapeutic target.

CONCLUSIONS

Ewing sarcoma tumors demonstrate significant transcriptional heterogeneity as well as a complex immunosuppressive microenvironment. This work evaluates several proposed targets that warrant further exploration as novel therapeutic strategies.

Perioperative nivolumab and chemotherapy in locally advanced squamous cell carcinoma of the oesophagus: a randomized multicentre phase 2 study with circulating tumor DNA dynamics monitoring

BACKGROUND

Although neoadjuvant chemotherapy and immunotherapy show promise in treating oesophageal squamous cell carcinoma (OSCC), long-term survival data are limited. This randomized, multicenter phase 2 study evaluated the efficacy of perioperative Nivolumab with chemotherapy, followed by surgery and adjuvant immunotherapy, in patients with locally advanced resectable OSCC, and explored the prognostic role of circulating tumor DNA (ctDNA) status.

METHODS

In this trial, participants recruited from five centers were randomly assigned in a 2:1 ratio to receive either perioperative Nivolumab or a placebo in addition to chemotherapy (cisplatin and paclitaxel), followed by minimally invasive esophagectomy. For those who did not achieve a pathological complete response (pCR), adjuvant treatment with Nivolumab was administered. The main measure of success was the pCR rate, with secondary endpoints including the R0 resection rate, event-free survival, and overall survival. All outcomes and safety measures were assessed based on the intention-to-treat population. ctDNA levels were monitored as exploratory endpoints.

RESULTS

Ninety patients were enrolled and randomized to Nivolumab or placebo plus chemotherapy. The pCR rate was slightly higher in the Nivolumab group (15%) compared to the control group (13.3%) (relative risk, 1.13; 95% CI, 0.38 to 3.36). No significant differences were observed in R0 resection rates (96.4% vs. 96.6%; P > 0.05). The median follow-up duration was 24.9 months (interquartile range: 22.8 to 26.7 months). Two-year event-free survival rates were 63.11% in the Nivolumab group versus 60.47% in the chemo group (hazard ratio, 0.97; 95% CI, 0.49 to 1.92). Two-year overall survival rates were 83.32% and 79.4%, respectively (hazard ratio, 0.82; 95% CI, 0.29 to 2.31). All participants were ctDNA positive at baseline, but post-treatment, 89% of the Nivolumab group and 62.5% of the placebo group turned ctDNA negative (P = 0.01). Those negative for ctDNA at all testing points showed significantly better disease-free survival (P < 0.001).

CONCLUSIONS

Perioperative Nivolumab plus chemotherapy is a viable and safe option for systemically treating locally advanced resectable OSCC. Monitoring minimal residual disease through ctDNA could be potentially valuable for assessing the effectiveness of adjuvant therapy and for prognostic evaluation in a systemic manner.

TRIAL REGISTRATION

ClinicalTrials.gov registration NCT05213312.

Genome-wide enhancer-gene regulatory maps of liver reveal novel regulatory mechanisms underlying NAFLD pathogenesis

INTRODUCTION

Non-alcoholic fatty liver disease (NAFLD) represents the most widespread liver disease globally, ranging from non-alcoholic fatty liver (NAFL) and steatohepatitis (NASH) to fibrosis/cirrhosis, with potential progression to hepatocellular carcinoma (HCC). Genome-wide association studies (GWASs) have identified several single nucleotide polymorphisms (SNPs) associated with NAFLD. However, numerous GWAS signals associated with NAFLD locate in non-coding regions, posing a challenge for interpreting their functional annotation.

RESULTS

In this study, we utilized the Activity-by-Contact (ABC) model to construct the enhancer-gene maps of liver by integrating epigenomic data from 15 liver tissues and cell lines. We constructed the most comprehensive genome-wide regulatory maps of the liver, identifying 543,486 enhancer-gene connections, including 267,857 enhancers and 16,872 target genes. Enrichment analyses revealed that the ABC SNPs are significantly enriched in active chromatin regions and active chromatin state. By combining the ABC regulatory maps and NAFLD GWAS data, we systematically identified ABC SNPs associated with NAFLD risk. Through the functional annotations, such as pathway enrichment and drug-gene interaction analyses, we identified 6 genes (GGT1, ACTG1, SPP1, EPHA2, PROZ and SHMT1) as candidate NAFLD genes, with SHMT1 previously reported. Among the SNPs connected to the candidate genes, the ABC SNP rs2017869 (odds ratio [OR] for the C allele = 1.10, 95% CI = 1.04-1.16, P = 5.97 × 10- 4) had the highest ABC score. According to the ABC maps, rs2017869 links to GGT1, and several drugs targeting this gene, such as liothyronine, showed potential benefits to patients with NAFLD. Furthermore, we identified that another novel gene, EPHA2, may play a crucial role in NAFLD by regulating the GGT levels.

CONCLUSIONS

Our study provides the most comprehensive ABC regulatory maps of the liver to date. This resource offers a valuable reference for identifying regulatory variants and prioritizing susceptibility genes of liver diseases, such as NAFLD.

Tandem inactivation of inositol pyrophosphatases Asp1, Siw14, and Aps1 illuminates functional redundancies in inositol pyrophosphate catabolism in fission yeast

Inositol pyrophosphates 5-IP7, 1-IP7, and 1,5-IP8 are eukaryal signaling molecules that influence cell physiology, especially phosphate homeostasis. In fission yeast, 1,5-IP8 and 1-IP7 impact gene expression by acting as agonists of RNA 3'-processing and transcription termination. 1,5-IP8 is synthesized by position-specific kinases Kcs1 and Asp1 that convert IP6 to 5-IP7 and 5-IP7 to 1,5-IP8, respectively. Inositol pyrophosphatase enzymes Asp1 (a histidine acid phosphatase), Siw14 (a cysteinyl phosphatase), and Aps1 (a Nudix hydrolase) are agents of inositol pyrophosphate catabolism in fission yeast. Whereas Asp1, Siw14, and Aps1 are individually inessential, double pyrophosphatase mutants asp1-H397A aps1∆ and siw14∆ aps1∆ display severe growth defects caused by overzealous 3'-processing/termination. By applying CE-ESI-MS to profile the inositol pyrophosphate content of fission yeast mutants in which inositol pyrophosphate toxicity is genetically suppressed, we elucidated the functional redundancies of the Asp1, Siw14, and Aps1 pyrophosphatases. Asp1, which exclusively cleaves the 1-β-phosphate, and Aps1, which prefers to cleave the 1-β-phosphate, play essential overlapping roles in guarding against the accumulation of toxic levels of 1-IP7. Aps1 and Siw14 together catabolize the inositol-5-pyrophosphates, and their simultaneous inactivation results in overaccumulation of 5-IP7. Cells lacking all three pyrophosphatases amass high levels of 1,5-IP8 and 1-IP7, with concomitant depletion of IP6. A genetic screen identified three missense mutations in the catalytic domain of Kcs1 kinase that suppressed inositol-1-pyrophosphate toxicosis. The screen also implicated the 3'-processing factor Swd22, the inositol pyrophosphate sensor Spx1, and the nuclear poly(A)-binding protein Nab2 as mediators of inositol-1-pyrophosphate toxicity.IMPORTANCEInositol pyrophosphates are key effectors of eukaryal cellular phosphate homeostasis. They are synthesized by kinases that add a β-phosphate to the 5- or 1-phosphate groups of IP6 and catabolized by three classes of pyrophosphatases that hydrolyze the β-phosphates of 5-IP7, 1-IP7, or 1,5-IP8. Whereas the fission yeast inositol pyrophosphatases-Asp1 (histidine acid phosphatase), Siw14 (cysteinyl phosphatase), and Aps1 (Nudix hydrolase)-are inessential for growth, Asp1/Aps1 and Aps1/Siw14 double mutations and Asp1/Siw14/Aps1 triple mutations elicit severe or lethal growth defects. By profiling the inositol pyrophosphate content of pyrophosphatase mutants in which this toxicity is genetically suppressed, we reveal the functional redundancies of the Asp1, Siw14, and Aps1 pyrophosphatases. Their synergies are manifested as excess accumulation of 1-IP7 upon dual inactivation of Asp1 and Aps1 or an excess of 5-IP7 in aps1∆ siw14∆ cells. In the absence of all three pyrophosphatases, cells accrue high levels of 1,5-IP8 and 1-IP7 while IP6 declines.

Two haplotype-resolved telomere-to-telomere genome assemblies of Xanthoceras sorbifolium

Yellowhorn (Xanthoceras sorbifolium) is widely used in northern China for landscaping, desertification control, and oil production. However, the lack of high-quality genomes has hindered breeding and evolutionary studies. Here, we present the first haplotype-resolved, telomere-to-telomere (T2T) yellowhorn genomes of PBN-43 (white single-flowered) and PBN-126 (white double-flowered) using PacBio HiFi and Hi-C data. These assemblies range from 464.34 Mb to 468.97 Mb and include all centromeres and telomeres. Genome annotation revealed that an average of 67.99% (317.09 Mb) of yellowhorn genomic regions consist of repetitive elements across all haplotypes. The number of protein-coding genes ranges from 35,039 to 35,174 among assemblies, representing an average 50.16% increase over the first published yellowhorn genome. Additionally, 93.90% of the annotated genes have functional annotations. We found yellowhorn experienced an LTR-RT burst during the last 0.45-0.48 Mya. These data provide a resource for investigating genomic variations, phylogenetic relationships, duplication modes, and the distribution of nucleotide-binding leucine-rich repeat (NLR) genes, and support further research into yellowhorn breeding.

Genome diversity and signatures of natural selection in mainland Southeast Asia

Mainland Southeast Asia (MSEA) has rich ethnic and cultural diversity with a population of nearly 300 million1,2. However, people from MSEA are underrepresented in the current human genomic databases. Here we present the SEA3K genome dataset (phase I), generated by deep short-read whole-genome sequencing of 3,023 individuals from 30 MSEA populations, and long-read whole-genome sequencing of 37 representative individuals. We identified 79.59 million small variants and 96,384 structural variants, among which 22.83 million small variants and 24,622 structural variants are unique to this dataset. We observed a high genetic heterogeneity across MSEA populations, reflected by the varied combinations of genetic components. We identified 44 genomic regions with strong signatures of Darwinian positive selection, covering 89 genes involved in varied physiological systems such as physical traits and immune response. Furthermore, we observed varied patterns of archaic Denisovan introgression in MSEA populations, supporting the proposal of at least two distinct instances of Denisovan admixture into modern humans in Asia3. We also detected genomic regions that suggest adaptive archaic introgressions in MSEA populations. The large number of novel genomic variants in MSEA populations highlight the necessity of studying regional populations that can help answer key questions related to prehistory, genetic adaptation and complex diseases.

Genome-wide allele-specific expression in multi-tissue samples from healthy male baboons reveals the transcriptional complexity of mammals

Allele-specific expression (ASE) is pivotal in understanding the genetic underpinnings of phenotypic variation within species, differences in disease susceptibility, and responses to environmental factors. We processed 11 different tissue types collected from 12 age-matched healthy olive baboons (Papio anubis) for genome-wide ASE analysis. By sequencing their genomes at a minimum depth of 30×, we identified over 16 million single-nucleotide variants (SNVs). We also generated long-read sequencing data, enabling the phasing of all variants present within the coding regions of 96.5% of assayable protein-coding genes as a single haplotype block. Given the extensive heterozygosity of baboons relative to humans, we could quantify ASE across 72% of the total annotated protein-coding gene set. We identified genes that exhibit ASE and affect specific tissues and genotypes. We discovered ASE SNVs that also exist in human populations with identical alleles and that are designated as pathogenic by both the PrimateAI-3D and AlphaMissense models.

The SIK3-N783Y mutation is associated with the human natural short sleep trait

Sleep is an essential component of our daily life. A mutation in human salt induced kinase 3 (hSIK3), which is critical for regulating sleep duration and depth in rodents, is associated with natural short sleep (NSS), a condition characterized by reduced daily sleep duration in human subjects. This NSS hSIK3-N783Y mutation results in diminished kinase activity in vitro. In a mouse model, the presence of the NSS hSIK3-N783Y mutation leads to a decrease in sleep time and an increase in electroencephalogram delta power. At the phosphoproteomic level, the SIK3-N783Y mutation induces substantial changes predominantly at synaptic sites. Bioinformatic analysis has identified several sleep-related kinase alterations triggered by the SIK3-N783Y mutation, including changes in protein kinase A and mitogen-activated protein kinase. These findings underscore the conserved function of SIK3 as a critical gene in human sleep regulation and provide insights into the kinase regulatory network governing sleep.

A deletion at the X-linked ARHGAP36 gene locus is associated with the orange coloration of tortoiseshell and calico cats

The X-linked orange (O) locus in domestic cats controls an unknown molecular mechanism that causes the suppression of black-brownish pigmentation in favor of orange coloration. The alternating black-brownish and orange patches seen in tortoiseshell and calico cats are considered classic examples of the phenotypic expression of random X chromosome inactivation (XCI) occurring in female mammals. However, the O gene in the cat genome has not been identified, and the genetic variation responsible for the orange coloration remains unknown. We report here that a 5.1-kilobase (kb) deletion within an intron of the X-linked ARHGAP36 gene, encoding a Rho GTPase-activating protein, is closely and exclusively associated with orange coloration. The deleted region contains a highly conserved putative regulatory element, whose removal is presumed to alter ARHGAP36 expression. Notably, ARHGAP36 expression in cat skin tissues is linked to the suppression of many melanogenesis genes, potentially shifting pigment synthesis from eumelanin to pheomelanin. Furthermore, we find evidence that the gene undergoes XCI in female human and mouse cells and XCI-dependent CpG island methylation consistent with random XCI in female domestic cats. The 5.1-kb deletion seems widespread in domestic cats with orange coat coloration, suggesting a single origin of this coat color phenotype.

Clinical response to azacitidine in MDS is associated with distinct DNA methylation changes in HSPCs

Hypomethylating agents are frontline therapies for myelodysplastic neoplasms (MDS), yet clinical responses remain unpredictable. We conducted a phase 2 trial comparing injectable and oral azacitidine (AZA) administered over one or three weeks per four-week cycle, with the primary objective of investigating whether response is linked to in vivo drug incorporation or DNA hypomethylation. Our findings show that injection results in higher drug incorporation, but lower DNA demethylation per cycle, while global DNA methylation levels in mononuclear cells are comparable between responders and non-responders. However, hematopoietic stem and progenitor cells (HSPCs) from responders exhibit distinct baseline and early treatment-induced CpG methylation changes at regulatory regions linked to tissue patterning, cell migration, and myeloid differentiation. By cycle six-when clinical responses typically emerge-further differential hypomethylation in responder HSPCs suggests marrow adaptation as a driver of improved hematopoiesis. These findings indicate that intrinsic baseline and early drug-induced epigenetic differences in HSPCs may underlie the variable clinical response to AZA in MDS.

Admixture and selection offer insights for the conservation and breeding of Guyuan cattle

BACKGROUND

The admixture between taurine and indicine cattle increases breed diversity and provides new genetic resources for human and natural selection. The climate of northwestern China is typified by cold and arid conditions, and cattle in this region are primarily taurine breeds. However, Guyuan cattle inhabit a transitional zone in northwestern China, typified by semi-arid and semi-humid climates. It is hypothesized that Guyuan cattle have a little of indicine ancestry. These results suggest that Guyuan cattle are a valuable genetic resource.

RESULTS

We established a conservation population of Guyuan cattle in their native habitat. We found that Guyuan cattle were an admixture between 78.09% East Asian taurine (EAT) and 20.26% East Asian indicine (EAI) ancestries. The admixture in Guyuan cattle was dated to 255 years ago. Notably, we identified Guyuan cattle as a unique genetic resource, representing a transitional breed between northern and central Chinese cattle with distinct admixture proportions. We revealed that the selection signals in Guyuan cattle with excess EAT ancestry were associated with reproduction, immunity, body length, cold climate adaptation, pigmentation, muscle development, residual feed intake, and fat deposition and that the selection signals in Guyuan cattle with excess EAI ancestry were associated with disease resistance. Remarkably, we discovered valuable single nucleotide polymorphisms (SNPs) in the promoter regions of the RBM39 and NEK6 genes, which may play key roles in regulating muscle development and disease resistance.

CONCLUSIONS

Our results suggest that Guyuan cattle are a newly identified genetic resource, and the native taurine and indicine ancestries in Guyuan cattle remain a valuable genetic resource of conservation and breeding.

Insight into the genetic underpinnings of tobacco hairy root formation by variant-associated genes based on whole-genome resequencing

MAIN CONCLUSION

Our whole-genome resequencing of tobacco hairy roots reveals functionally relevant variations in secondary metabolism-related genes and NAC transcription factors, providing actionable targets for metabolic pathway optimization and bioreactor design. Tobacco hairy roots are a critical model system for studying plant root development and secondary metabolism. The in-depth analysis of their genetic background and molecular regulatory mechanisms is important for biotechnological applications. In this study, we performed a whole-genome resequencing of tobacco hairy roots to uncover their genomic variation characteristics and potential functional implications. Genes associated with stop-lost, stop-gained, start-lost, and premature-start-codon-gain variants were enriched in zeatin biosynthesis, flavonoid biosynthesis, and glycosyltransferase activities. The results of metabolite content determination showed that hairy roots possessed a low content of zeatin and flavonoid but a higher content of glycoside compounds. Among transcription factors associated with effective variants, NAC transcription factors constituted the largest proportion. Further characterization of NAC proteins revealed their functional domains and expression patterns. This study not only explores the molecular genetic underpinnings of tobacco hairy roots but also provides a critical dataset for metabolic engineering optimization, development of efficient bioreactors, and plant-microbe interaction mechanisms research.

First identification of ORF virus causing contagious ecthyma in Morocco (MOR20): Genomic, phylogenetic, and sequence variants analyses for vaccine design

The ORF virus induces a zoonotic contagious ecthyma disease, affecting small ruminants such as sheep and goats. ORF virus has not been identified in Morocco, and there is no vaccination protocol against contagious ecthyma. In this study, we analyzed the genome sequence of a new strain isolated in Morocco (MOR20) from a flock of sheep showing suspicious signs of Sheepox virus infection. ORFV MOR20 strain was isolated after 2 initial blind passages on Heart cells. The cytopathic effect was characterized by aggregation, swelling and detachment of cells, appearing 4 days after infection. The virus was harvested on day 6 pi with a titer of 107.2 TCID/ml. ORFV MOR20 was sequenced using the Illumina NovaSeq 6000 platform. After employing several bioinformatics tools, we identified that ORFV MOR20 shares 98.59% similarity with the TVL strain virus, which is used in a commercial live vaccine. Additionally, We aligned 33 ORFV genomic sequences with MOR20 sequences and visualized the pairwise comparisons using a Heat Map. ORFV was classified into two genetic groups: those isolated from sheep and those from goats. This was confirmed by a phylogenetic tree. Furthermore, we analyzed genetic variants identified in the MOR20 strain in comparison with ORFV TVL strain and found 636 sequence variants. Some genes, such as ORFV086, ORFV112, and ORFV132, have a particularly high number of sequence variants. All in all, ORFV MOR20 isolate represents a promising candidate for further studies aimed at developing a standardized vaccine against contagious ecthyma.

Genome wide population genetics and molecular surveillance of insecticide resistance in Anopheles stephensi mosquitoes from Awash Sebat Kilo in Ethiopia

Since the detection of the Asian mosquito Anopheles stephensi in Dijbouti in 2012, it has spread throughout the Horn of Africa. This invasive vector continues to expand across the continent and is a significant threat to malaria control programs. Vector control methods, including insecticide-treated nets and indoor residual spraying, have substantially reduced the malaria burden. However, the increasing prevalence of mosquitoes resistant to insecticides, including An. stephensi populations, undermines ongoing malaria elimination efforts. Understanding population structure, gene flow between populations, and the distribution of insecticide resistance mutations is essential for guiding effective malaria control strategies. Here, we generated whole genome sequencing data for An. stephensi sourced from Awash Sebat Kilo, Ethiopia (n = 27) and compared with South Asian populations (n = 45; India and Pakistan) to assess genomic diversity, population structure, and uncovering insecticide resistance mutations. Population structure analysis using genome-wide single nucleotide polymorphisms (n = 15,533,476) revealed Ethiopian isolates clustering as a distinct ancestral group, separate from South Asian isolates. Three insecticide resistance-associated SNPs (gaba gene: A296S and V327I; vgsc L1014F) were detected. Evidence of ongoing selection was found in several loci, including genes previously associated with neonicotinoids, ivermectin, DDT, and pyrethroid resistance. This study represents the first whole genome population genetics study of invasive An. stephensi, revealing genomic differences from South Asian populations, which can be used for future assessments of vector population dispersal and detection of insecticide resistance mechanisms.

Disrupting cross-adaptation in high-risk MRSA: Sanguinarine as a multi-effective stress sensitizer for environmental and food safety

Methicillin-resistant Staphylococcus aureus (MRSA) represents a significant public health concern owing to its formidable antibiotic resistance and robust capacity for biofilm formation. The cross-adaptation mechanism enables MRSA to develop tolerance to environmental stressors such as antibiotics, acid, heat and osmotic pressure, leading to the persistence infections and environmental contamination. The cross-adaptation mechanism enables MRSA to develop tolerance to environmental stressors, such as antibiotics, acid, heat and osmotic pressure, leading to the persistence infections and environmental contamination. Here, we identified 261 strains of S. aureus and 9 high-risk MRSA from the environment of dairy farms and raw milk. The natural product Sanguinarine (SAN), derived from feed additives, exhibits effective anti-MRSA and anti-biofilm activity. Notably, SAN enhances the sensitivity of MRSA to antibiotics, acid, heat, and osmotic pressure by disrupting the cross-adaptation mechanism. Mechanistic investigations revealed that SAN significantly reduces the transcriptional level of type I (dnaK, groEL, etc.) and type III (clpB, clpP, etc.) heat stress response genes while markedly upregulating type II (σB) gene. Furthermore, SAN upregulates Na+/H+ antiporters activity, F0F1-ATPase activity and purine metabolism, while broadly downregulating DNA damage repair genes and disrupting ribosomal function. Additionally, SAN induces non-synonymous mutations in key stress response factors ClpB/L, leading to a loss of conformational homeostasis. SAN elicits a distinct stress response compared to environmental stressors, weakening MRSA's resilience and demonstrating promising capabilities for MRSA clearance and biofilm inhibition. Overall, SAN provides an effective strategy for the clearance of high-risk MRSA and the assurance of public health security.

A BRASSINOSTEROID INSENSISTIVE 1 receptor kinase ortholog is required for sex determination in Ceratopteris richardii

Most ferns, unlike all seed plants, are homosporous and produce sexually undifferentiated spores. Sex ratio in many homosporous species is environmentally established by the secretion of antheridiogen from female/hermaphrodite gametophytes. Nearby undetermined gametophytes perceive antheridiogen, which induces male development. In the fern Ceratopteris richardii (Ceratopteris), hermaphroditic (her) mutants develop as hermaphrodites even in the presence of antheridiogen. Modern sequencing and genomic tools make the molecular identification of mutants in the 11-Gbp genome of this fern possible. We mapped 2 linked mutants, her7-14 and her7-19, to the same 16-Mbp interval on chromosome 29 of the Ceratopteris genome. An ortholog of the receptor kinase gene BRASSINOSTEROID INSENSITIVE 1 (BRI1) within this interval encoded a deletion mutation in her7-14 and a missense mutation in her7-19. Three other linked her mutants encoded missense mutations in the same gene, which we name HER7. Consistent with a function as a receptor kinase, HER7-GFP fusion protein localized to the plasma membrane and cytoplasm. Analysis of gene expression showed that brassinosteroid biosynthesis was upregulated in hermaphrodites compared with male gametophytes. Our work demonstrates that HER7 is required for sex determination in Ceratopteris and opens avenues for studying the evolution of antheridiogen systems.

The fourspine stickleback (Apeltes quadracus) has an XY sex chromosome system with polymorphic inversions on both X and Y chromosomes

Teleost fish are well-known for possessing a diversity of sex chromosomes and for undergoing frequent turnovers of these sex chromosomes. However, previous studies have mainly focused on variation between species, while comparatively little attention has been given to sex chromosome polymorphisms within species, which may capture early stages of sex chromosome changes. To better understand the evolution of sex chromosomes, we used the fourspine stickleback (Apeltes quadracus) as a model organism. Previous cytogenetic studies suggested that females of this species possessed a ZW heteromorphic sex chromosome system. However, genetic crosses and our whole-genome sequencing of three geographically distinct wild populations revealed that A. quadracus has an XY sex chromosome on chromosome 23. This chromosome has not previously been identified as a sex chromosome in any other stickleback species, indicating a recent sex chromosome turnover. We also identified two genes - rxfp2a and zar1l - as novel candidate sex determination genes. Notably, we observed inversions on both the X and Y chromosomes in different populations, resulting in distinctive patterns of differentiation between the X and Y chromosomes across populations. The new sex chromosome and intraspecies inversion polymorphisms observed in A. quadracus provide an excellent system for future work assessing the relative fitness effects of the inversions, which will enable testing theoretical models about the drivers of sex chromosome evolution and turnover.