Search PubMedSearch

SEARCH · Search PubMed

Results for “whole‐genome sequencing”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

772 records · Page 4Linked to original sources

Genomic epidemiology of extended-spectrum beta-lactamase-producing Escherichia coli across humans, poultry and wastewater sectors in Douala, Cameroon.

BACKGROUND: The global health threat of antimicrobial resistance involves the human, animal and environmental sectors. Data from Cameroon are scarce. OBJECTIVES: This study aimed to define extended-spectrum beta-lactamase-producing Escherichia coli (ESBL-Ec) rates and associated risk factors across the three sectors in Douala, Cameroon, and to define molecular characteristics of isolates. METHODS: From June 2022 to May 2023, we collected blood cultures from hospitalized patients, rectal swabs from healthy pregnant women, caeca from broiler chickens and environmental wastewater. Samples were screened for ESBL-Ec using CHROMAgar™ ESBL and cefotaxime-supplemented Tryptone Bile X-glucuronide agar. Antimicrobial susceptibility testing was performed by disk diffusion following EUCAST guidelines. Whole-genome sequencing was carried out using Illumina technology. RESULTS: Of 628 samples, 374 yielded ESBL-Ec. Prevalence was 54.6% (131/240) in pregnant women, 70.4% (169/240) in chickens and 93.1% (67/72) in wastewater. The proportion of ESBL-Ec among E. coli-positive-blood cultures was 9.2% (7/76). Multi-family household living was independently associated with ESBL-Ec carriage among pregnant women (adjusted odds ratio = 1.7, 95% CI 1.0-3.1, P = 0.03). High co-resistance (>70%) was observed for tetracycline, ciprofloxacin and trimethoprim/sulfamethoxazole. Sequencing of 32 isolates revealed 45 distinct resistance genes, including blaCTX-M-15 (n = 13, 40.6%), blaCTX-M-55 (n = 11, 34.4%) and last-resort antibiotic resistance genes mcr-1 and bla OXA-181. High-risk sequence types included ST131 (pregnant women) and ST10 (chickens). Notably, ST48 was shared between pregnant women and chickens, and ST155 between pregnant women and wastewater. CONCLUSION: Cross-sectoral ESBL-Ec in Douala exhibits high genomic diversity and alarming resistance. The occurrence of last-resort genes requires immediate One Health surveillance and coordinated interventions.

Journal Article

Novel splice site variants in GBA1 are associated with Gaucher disease and genotype-phenotype correlations.

BACKGROUND: Variants in GBA1 are associated with neurodegenerative disease. This study aimed to explore pathogenic GBA1 variants. METHODS: Four patients with progressive myoclonic epilepsy (PME) and extremely low β-glucosidase levels were recruited. Whole-exome sequencing and long-range PCR were performed to identify GBA1 variants. Bioinformatic analyses were used to predict the impact of the identified variants. A literature review was performed to explore the genotype-phenotype correlations. GBA1 expression data across different brain regions and developmental stages were analyzed using the BrainSpan database. RT-PCR was performed to verify the splicing effects. RESULTS: Compound heterozygous GBA1 variants were identified in four patients. Five distinct variants were detected, including two novel splice site variants (c.308-2A>G and c.762-2A>C) and three previously reported variants. All identified variants were rare or absent in gnomAD. Splice site variants c.308-2A>G and c.762-2A>C were predicted to cause aberrant splicing. Minigene-based splicing assays coupled with RT-PCR and Sanger sequencing confirmed that both variants cause complete exon skipping (exon 4 and exon 7, respectively). All patients presented with PME onset in childhood/adolescence, intellectual regression, low β-glucosidase, and diffuse brain atrophy and were subsequently diagnosed with Gaucher disease type 3. GBA1 expression in the brain showed two distinct peaks: one in infancy and another after five years of age. The onset age of PME aligned with the second GBA1 expression peak (after five years of age). CONCLUSION: This study identified compound heterozygous GBA1 variants, including two novel candidate pathogenic splice site variants, in Gaucher disease type 3 patients, expanding the known mutational spectrum.

Humans

Structural and tissue-specific organisation of endocrine Fgf19 and Fgf21 signalling in rainbow trout.

Endocrine fibroblast growth factors (FGF19 subfamily) play a key role in regulating metabolic homeostasis in vertebrates. However, their functional diversification in salmonids remains poorly understood. In this study, we conducted an integrative characterisation of Fgf19 and Fgf21 signalling in rainbow trout (Oncorhynchus mykiss) by combining phylogenetic, structural and expression analyses. Phylogenetic analyses revealed the conservation of single fgf19 and fgf21 genes, despite the extensive expansion of receptors post-Ss4R (salmonid-specific fourth-round whole genome duplication). Structural modelling and molecular dynamics simulations demonstrated the stable interactions of both ligands to multiple Fgfr isoforms, with receptor-specific energetic profiles and conserved core interaction residues. Tissue expression profiling revealed clear differences from mammalian models, such as predominant hepatic fgf19 expression and the absence of hepatic fgf21 under basal conditions. In addition, there were complex and tissue-dependent distributions of fgfr and klotho transcripts. These findings support a receptor-driven diversification model of endocrine Fgf signalling in salmonids, suggesting enhanced endocrine plasticity associated with the retention of receptors following post-genomic duplication. Taken together, our findings provide new insights into the structural and regulatory organisation of endocrine Fgf signalling, as well as its potential role in metabolic regulation in rainbow trout.

Animals

Mul-PheG2P: decoupled learning and prediction-space fusion enables robust and interpretable multi-phenotype genomic prediction.

Genomic prediction of multiple phenotypes is crucial in modern plant breeding; however, existing methods struggle with negative transfer and lack interpretability, particularly across high-dimensional small-sample data and diverse species. To address this, we propose Mul-PheG2P, a novel paradigm based on decoupled learning and predictive space fusion. It employs a two-stage design: first training phenotype-specific encoders using genetic data, then decoupling phenotype-specific learning from cross-phenotype aggregation via an interpretable prediction layer. Mul-PheG2P outperforms existing methods across diverse crop datasets, including maize (Zea mays), wheat (Triticum aestivum), and tomato (Solanum lycopersicum). It provides a multi-scale interpretability chain: at the macro level, it quantifies phenotypic contributions via attention-based weighting; at the micro level, Integrated Gradients reveal the genetic basis of predictions. Notably, the model successfully identified the CCT (CONSTANS, CO-like, and TOC) motif regulating photoperiodism and the SQUAMOSA (SQUAMOSA promoter binding protein) promoter for inflorescence development, confirming its ability to capture functional biological mechanisms. These results highlight the high performance and interpretability of Mul-PheG2P, showcasing its value for low-cost, large-scale screening to advance precision breeding.

Phenotype

Serotypic and Genomic Diversity of Vibrio anguillarum in Rainbow Trout Farms in Turkey: Implications for Vibriosis Control and Vaccine Candidate Selection.

Outbreaks of vibriosis caused by Vibrio anguillarum are a persistent constraint on rainbow trout (Oncorhynchus mykiss) aquaculture. However, information on the population structure of field strains in Turkey has been lacking. Here, we report the first systematic serotypic, proteomic, and genomic characterization of 23 V. anguillarum isolates collected over 10&#x2009;years from rainbow trout farms located in six major aquaculture regions of Turkey. Serological analyses based on microagglutination, supported by ELISA characterization of hyperimmune sera, identified a clear predominance of serotype O1, whereas isolate V12 exhibited a non-agglutinating, atypical O-antigen profile. Protein profiling (SDS-PAGE) and immunoblotting showed largely conserved whole-cell protein patterns among the isolates, but distinct immunogenic bands at 14, 18, and 40&#x2009;kDa were detected in isolates V18 and V21. Long-read whole-genome sequencing revealed that most Turkish isolates grouped within the global O1 clade, while V12, V25, and V28 isolates occupied more distant branches. Comparative genomics demonstrated a conserved core virulence gene set (RTX toxins, siderophore and iron-uptake systems, motility and adhesion factors, Type VI secretion system), with strain-dependent variation in accessory loci such as anguibactin and T6SS-I. Experimental infections of rainbow trout demonstrated significant differences in virulence among isolates (p&#x2009;<&#x2009;0.05), with the V18 isolate showing high, the V15 intermediate, and the V12 low-mortality rates. By elucidating the relationship among the serotype, immunogenic protein profiles, virulence gene repertoires, and in&#xa0;vivo pathogenicity, this study provides a comprehensive overview of the antigenic and genomic diversity of Vibrio anguillarum isolates from Turkey. Notably, the identification of V18 and V21 as promising candidate strains for further vaccine evaluation, characterized by high virulence and unique immunogenic features, provides a scientific foundation for the development of serotype-specific vaccination strategies to mitigate vibriosis-associated losses in aquaculture.

Animals

Diverse structures of mcr-10-bearing plasmids and high colistin resistance in Enterobacter cloacae complex clinical isolates from South Korea.

BACKGROUND: The emergence of mcr-mediated colistin resistance in Enterobacter cloacae complex (ECC) poses a significant threat to antimicrobial therapy. Among mcr variants, mcr-10 has been identified in various environments, but its genetic diversity, structural context, and functional role in colistin resistance remain unclear. METHODS: We investigated 183 ECC isolates and identified 60 colistin-resistant strains through minimum inhibitory concentration (MIC) testing. The presence of mcr-10 was screened using reference genomes from NCBI, and whole plasmid sequencing was conducted on mcr-10-positive isolates. The genetic environment of mcr-10 was analyzed via synteny and structural annotation. RESULTS: Whole-plasmid sequencing of eight mcr-10-positive ECC isolates identified three replicon types among the mcr-10-harboring plasmids: IncFIB (n=3), IncFII (n=3), and IncFII/IncFIB (n=2). Although all plasmids shared the xerC-mcr-10 cassette, they lacked a conserved backbone and showed diverse genetic contexts with variable insertion sequences near mcr-10, indicating marked structural heterogeneity. No additional antimicrobial resistance genes were detected on these plasmids. When introduced into E. coli DH5&#x3b1;, the plasmids increased colistin MICs only modestly, whereas representative plasmids transferred into colistin-susceptible E. roggenkampii and E. kobei conferred high-level resistance comparable to that of the parental mcr-10-positive isolates. Consistently, qRT-PCR showed higher mcr-10 expression in ECC transformants than in E. coli under colistin exposure, supporting a hostdependent effect on mcr-10-mediated colistin resistance. CONCLUSION: These findings highlight the diversity of mcr-10-carrying plasmids and suggest that mcr-10-mediated colistin resistance is shaped by the host genetic background. Further studies are needed to clarify the mechanisms underlying mcr-10 expression.

Colistin

Safe and Stable Germline Transmission of MSTN Mutations in Cattle.

With the global population expected to reach 10 billion by 2050, sustainable livestock production is critical. Gene editing of the myostatin (MSTN) gene represents a promising strategy to enhance muscle growth in cattle. In this study, MSTN-mutated founder (F0) cows were used to generate F1 offspring via ovum pick-up, in&#xa0;vitro fertilization, and embryo transfer. Four F1 calves were born, all confirmed to be heterozygous for the MSTN mutation. Long-term monitoring showed normal growth and no visible health abnormalities. Whole-genome sequencing identified SNPs, INDELs, and structural variants, most with minimal predicted functional effects. Proteomic profiling of Longissimus dorsi muscle quantified 2947 proteins, revealing only subtle expression differences between MSTN-mutated and wild-type cattle. These results demonstrate stable inheritance and confirm that MSTN editing does not disrupt genome integrity or protein expression. Overall, our findings support the safety and utility of MSTN gene editing to improve livestock productivity for future food security.

Animals

Prenatal exome sequencing of fetuses with central nervous system anomalies based on prenatal ultrasound and magnetic resonance imaging diagnosis: A retrospective cohort study with a systematic review and meta-analysis.

INTRODUCTION: Fetal central nervous system (CNS) abnormalities have diverse etiologies, with genetic factors as a major contributor. Prenatal exome sequencing (ES) is a powerful tool for precise molecular diagnosis of CNS anomalies, but its diagnostic yield varies among studies. This study aimed to evaluate the additional diagnostic yield of prenatal ES compared with chromosomal microarray analysis (CMA) in fetuses with CNS anomalies detected by prenatal imaging. MATERIAL AND METHODS: We collected ES results from fetuses diagnosed with CNS anomalies by prenatal imaging (2019-2024) who had negative results. Subgroup analyses assessed phenotype-specific ES diagnostic yield for associated genes and variants. A systematic review and meta-analysis incorporating our data and published studies further explored the association between phenotype and diagnostic yield. RESULTS: In the cohort study of 219 cases, ES identified pathogenic/likely pathogenic single nucleotide variations in 36 cases (16%). The highest diagnostic yield of ES was in cases with multisystem malformations (25%, 14/55), followed by multiple CNS anomalies (15%, 2/13) and isolated CNS anomalies (13%, 20/151). The most commonly identified isolated CNS anomaly was agenesis of the corpus callosum (31%, 5/16). Neural tube defects with urogenital anomalies were associated with a positive ES finding in 57% (4/7) of cases. The meta-analysis of 989 cases from 22 studies showed a pooled diagnostic yield of ES of 27% (95% CI, 21%-34%). The highest diagnostic yield of ES was in cases of corpus callosum anomalies with facial abnormalities (75%, 8/11) and neural tube defects with urogenital malformations (80%, 12/15). The diagnostic yield of ES for three or more CNS abnormalities was 43% (95% CI, 31%-58%), significantly higher than that for only two abnormalities (10%, 95% CI, 4%-18%). No significant difference in diagnostic yield was found between cases identified by prenatal MRI combined with ultrasound (27%, 95% CI, 20%-36%) and those identified by ultrasound alone (25%, 95% CI, 17%-35%). CONCLUSIONS: ES provided a significantly higher diagnostic yield than CMA for fetal CNS abnormalities, with diagnostic yields varying by phenotype. The systematic review and meta-analysis confirmed that the complexity and combination of malformations are key factors associated with differences in ES diagnostic yield.

Humans

Pervasive hybridization and introgression in Diervilleae (Caprifoliaceae).

Diervilleae (Caprifoliaceae) is a horticulturally important lineage with striking floral diversity and a long history of interspecific crossing, suggesting reticulate evolution. We integrated nuclear SNPs and whole plastome data to reconstruct a phylogenomic backbone for the tribe and to identify hybrids, cultivated accessions, and introgression among lineages. Nuclear and plastid phylogenies consistently recover Weigela and Diervilla as reciprocally monophyletic and resolve four major lineages within Weigela, providing a reproducible framework for revising sectional limits and species boundaries. Cultivated accessions form a well supported clade sister to W. florida and show predominantly W. florida ancestry while retaining contributions from multiple wild lineages, consistent with recurrent crossing, backcrossing, and selection. Analyses of wild populations reveal recurrent hybrids and enable plausible parental combinations to be inferred. Tests across the genome further indicate strong evidence for historical introgression across Diervilleae, with the strongest signals involving W. middendorffiana, W. maximowiczii, and Diervilla. Fossil evidence, divergence time estimation, and paleodistribution modelling together suggest range expansion during the Miocene and Pliocene followed by climate driven contraction, providing a spatiotemporal context for episodic contact, introgression, and the East Asia-North America disjunction.

Hybridization, Genetic

Tripled-Stranded Antisense Oligonucleotide for Biomarker-Activated Suppression of Essential Genes.

Conditional activation of antisense oligonucleotides (ASOs) is a promising strategy for selective suppression of cancer cells without affecting normal cells. In this study, we developed a tripled-stranded ASO (tsASO) that is rendered inactive through complexation with two additional oligonucleotides. The key innovation is the use of partial overlap between the parent ASO and the biomarker sequence, combined with toehold-mediated strand displacement, enabling precise conditional activation. The tsASO effectively triggered RNase H-mediated degradation of DYNC1I2 and DARS1 RNAs exclusively in the presence of the ERBB2 sequence. In cell-free systems, the tsASO demonstrated high cleavage efficiency (up to 81%), comparable to the parent ASO efficiency, with minimal background activity in the absence of the biomarker sequence, validating the concept at the molecular level. However, in cells using lipid-based transfection, the tsASO exhibited nonspecific cytotoxicity that did not correlate with biomarker presence or target gene expression. Detailed analysis showed no clear support for known sequence-driven toxicity mechanisms (CpG/TLR9, G-quadruplexes) in the nonimmune cell lines, suggesting that the primary limitation is intracellular delivery rather than the tsASO design. Future work should focus on optimizing delivery platforms to achieve controlled cellular uptake and biomarker-dependent release, unlocking the therapeutic potential of this conditional gene silencing approach.

Oligonucleotides, Antisense

Fructophilic lactic acid bacteria as a window into multi-scale convergent evolution.

Fructophilic lactic acid bacteria (FLAB) are a group of lactic acid bacteria with unique growth characteristics, that is, poor growth on glucose. Their growth is enhanced in the presence of fructose or external electron acceptors. These organisms inhabit fructose-rich environments such as flowers, fruits, and pollinating insects, particularly honey bees. Apilactobacillus spp. and Fructobacillus spp. are representatives of FLAB, although they belong to phylogenetically distant clades. These organisms commonly possess markedly small genomes with a low number of coding DNA sequences. Furthermore, their genomes are characterized by a markedly reduced number of genes involved in carbohydrate transport and metabolism. Genome reduction in FLAB reflects convergent adaptation to fructose-rich environments rather than general genome streamlining. The two distinct FLAB genera, Fructobacillus and Apilactobacillus, independently lost more than 100 genes in statistically similar orders. In contrast, genes involved in carbohydrate and amino acid metabolism exhibited reversed orders of loss between the two genera. Furthermore, FLAB genomes lack an intact bifunctional alcohol/aldehyde dehydrogenase gene (adhE), which causes their poor growth on glucose. A comparative genomic study suggested the evolutionary process underlying adhE gene decay during adaptation to the fructose-rich environments, including pollinating insects. In conclusion, FLAB represent a unique example of habitat-driven convergent reductive evolution that can be investigated across multiple biological scales - from individual genes to whole genomes - in the diverse LAB group with a wide range of habitats, and partially share the fructophilic evolution with eukaryotic yeasts found in fructose-rich habitats.

Fructose

Reference genome of the Californian trapdoor spider Aptostichus stephencolberti Bond 2008 (Araneae: Mygalomorphae: Euctenizidae).

We present a reference genome assembly for the trapdoor spider Aptostichus stephencolberti. This species, described in 2008, is endemic to the highly fragmented coastal dune habitats of Northern California from Monterey to the San Francisco Bay Area. Trapdoor spiders are ideal taxa for landscape scale genomic studies owing to their extreme site fidelity and limited dispersal capabilities; these same characteristics make them prone to extinction. Genomic studies of species like A. stephencolberti can reveal novel areas of endemism and high conservation value that may not be evident in species with wider ranges and greater dispersal capabilities. As part of the California Conservation Genomics Project, we constructed the A. stephencolberti reference genome from high quality long-read sequences, scaffolded with proximity ligation Omni-C data. The primary assembly comprises 551 scaffolds spanning 3.63 Gbp, a scaffold N50 of 62.2 Mbp and BUSCO completeness of 95.6%. We estimate 52 chromosomes yet find no (TTAGG)n telomer repeats. Expanding the telomeric repeat search finds an ancestral loss of the repeat from all spiders. Automated annotation using the NCBI refseq pipeline and RNAseq data from whole adults finds 14,067 genes with a BUSCO annotation completeness of 95.56%. Repeat annotation identified 77% of the genome to be interspersed repeats. This resource, the first for family Euctenizidae will facilitate future study and resulting conservation actions of A. stephencolberti and other Aptostichus sp. populations associated with the rapidly changing California coastal dune ecosystem.

Aptostichus stephencolberti

[Analysis of a Chinese pedigree affected with Townes-Brocks syndrome due to a novel variant of SALL1 gene and a literature review].

OBJECTIVE: To analyze a novel exonic variant of the SALL1 gene and its impact on the binding site of SALL protein. METHODS: Clinical data of three children diagnosed with Townes-Brocks syndrome and their family members who had presented at the First Affiliated Hospital of Shandong First Medical University in April 2022 were retrospectively collected. The pathogenic variant was identified through whole-genome sequencing (WGS) and validated by Sanger sequencing. Protein structural prediction was performed using AlphaFold and PyMOL software to construct three-dimensional models of the wild-type and mutant proteins. Additionally, previously reported cases were systematically reviewed. This study was approved by the Medical Ethics Committee of the hospital (Ethics No.: 2023-386). RESULTS: The proband was one of triplet sisters born at 34+4 gestational weeks. All three cases had presented with anal atresia and rectovaginal fistula, and case 3 also had toe malformation of left foot. WGS revealed a novel heterozygous c.757C>T (p.Gln253*) variant in the SALL1 gene, which was predicted to be pathogenic. Sanger sequencing confirmed co-segregation of the variant with the disease within the family. Protein structural modeling demonstrated that the variant has introduced a premature stop codon at position 253, resulting in a truncated protein. CONCLUSION: Above finding has enriched the mutation spectrum of the SALL1 gene in association with Townes-Brocks syndrome, which also represented a rare case of anal atresia in triplets, and provided a basis for molecular diagnosis, genetic counseling, and further research.

Humans

Genome-wide identification, structural characterization, and evolutionary analysis of growth-related gene families in African catfish (Clarias gariepinus).

The somatotropic axis encompassing growth hormone (GH), insulin-like growth factor (IGF), myostatin (MSTN), and prolactin (PRL) signalling cascades is the master regulator of somatic growth, metabolism, and development in vertebrates. African catfish (Clarias gariepinus), a commercially pivotal aquaculture species, now possesses a chromosome-level reference genome (CGAR_prim_01v2); however, a systematic, genome-wide characterization spanning all five interconnected growth-related gene families has not previously been undertaken in this species. Here, we identified and characterized 15 growth-related genes spanning gh1, ghra, ghrb, Igf1, Igf2a, Igf2b, igf1ra, Igf1rb, Igf2r, Mstna, Mstnb, prl, prlra, prlrb, and smtlb distributed across 13 chromosomes. Complete one-to-one orthology with zebrafish confirmed strong dosage-balance conservation across >120 million years of teleost divergence. Physicochemical analysis resolved a clear biochemical dichotomy between compact, basic secreted ligands (19.88-45.81&#xa0;kDa; pI up to 10.02) and large, acidic, heavily glycosylated membrane receptors (56.82-270.80&#xa0;kDa; pI 4.85-5.97). Phylogenetic analysis confirmed 3R whole-genome duplication origins for all paralog pairs, while synteny analysis revealed a disruption of the ancestral gh1-prl chromosomal block in C. gariepinus, a finding that warrants further comparative and functional investigation. This genomic atlas provides the sequence and structural information including exon-intron boundaries, domain architecture, and chromosomal coordinates needed as a prerequisite for future marker-assisted selection and CRISPR-based myostatin-editing efforts in African catfish aquaculture, though translation into applied breeding outcomes will require subsequent functional and expression studies.

Animals

Genome-wide analysis of the plant-specific PLATZ gene family in Taraxacum kok-saghyz and its roles in response to drought and salt tolerance.

Abiotic stress severely limits plant growth and productivity. Taraxacum kok-saghyz Rodin (TKS), known for its environmental resilience, represents a valuable resource for identifying stress-tolerant genes to improve stress-adaptive crops. Plant AT-rich protein and zinc-binding protein (PLATZ) transcription factors serve as core regulators of plant growth, developmental processes, and adaptive responses to various stress conditions; however, they remain uncharacterized in TKS. Here, we identified 10 TksPLATZ genes through a whole-genome analysis. Phylogenetically, these genes were grouped into five distinct evolutionary branches. Promoter sequence analysis revealed multiple types of cis-acting regulatory elements that are connected with hormonal signal responses and environmental stress adaptation. Integrated analysis of transcriptome datasets and RT-qPCR validation demonstrated that TksPLATZ genes display tissue-specific expression profiles and show distinct responsive patterns to drought and salt stress treatments. Among them, TksPLATZ1, TksPLATZ2 and TksPLATZ7 were markedly induced under both stressors and were selected for further functional study. We demonstrated that TksPLATZ1, TksPLATZ2 and TksPLATZ7 localize to the cell nucleus and act as transcriptional activators and repressors, respectively. Phenotypic data from overexpression experiments in plants confirm that heterologous expression of TksPLATZ1, TksPLATZ2, and TksPLATZ7 enhances the tolerance of Arabidopsis to salt and osmotic stress. These findings provide valuable genetic resources for improving plant tolerance to environmental stresses.

Salt Tolerance

ReMeDy: A Flexible Statistical Framework for Region-Based Detection of DNA Methylation Dysregulation.

Region-based epigenome-wide association studies have demonstrated improved statistical power and biological interpretability compared with probe-wise analyses of DNA methylation data. However, most existing region-based methods characterize methylation dysregulation primarily through changes in mean methylation levels associated with a phenotype of interest. Substantial evidence indicates that phenotype-associated methylation alterations may also manifest through changes in methylation variability or through joint shifts in mean and variability. Despite this, no existing statistical framework jointly models mean-variance methylation changes in a region-based manner. We propose ReMeDy, a flexible statistical framework that uses a hierarchical likelihood approach within a generalized linear model setting to identify differentially methylated regions, variably methylated regions, and regions exhibiting joint differential and variable methylation at a genome-wide scale. Unlike existing models, ReMeDy operates directly on biologically defined co-methylated regions, allowing it to naturally capture spatial correlation inherent in DNA methylation array data, while avoiding reliance on heuristic, user-defined tuning parameters such as smoothing spans and kernel bandwidths that can substantially influence results and introduce subjectivity. Through extensive simulation studies and comprehensive benchmarking against popular models, we demonstrate that ReMeDy maintains false discovery and Type-I error rates at nominal levels while achieving consistently higher statistical power across a wide range of realistic scenarios. Application to population-level DNA methylation data further shows that ReMeDy identifies biologically meaningful regions and pathways implicated in complex human diseases that are not captured by conventional mean-based analyses alone. ReMeDy is implemented as an open-source R package and is freely available at https://github.com/SChatLab/ReMeDy.

DNA Methylation

Genomic and Molecular Interaction Analysis of NodD1 in a Novel Bradyrhizobium yuanmingense sp. B64 Isolate for Nodulation and Symbiosis of Legume Plants.

Rhizobial bacteria are known for their ability to fix nitrogen for leguminous plants and their essential function for sustainable agriculture. This study characterizes the taxonomic status and functional potential of the Bradyrhizobium B64 isolate using integrated genomic and molecular approaches. The whole genome of the B64 isolate was sequenced via Illumina paired-end technology. Species delimitation was performed using average nucleotide identity (ANI) and digital DNA-DNA Hybridization (dDDH). The NodD1 protein structure was modeled using AlphaFold3 and validated by Ramachandran plot analysis. Molecular docking was then conducted to evaluate interactions between NodD1 and four signaling flavonoids: Apigenin, Daidzein, Genistein, and Naringenin. Genomic analysis revealed a maximum ANI of 94.4% and dDDH values between 51.4 and 62.4%. Since these values fall below the standard prokaryotic thresholds (ANI&#x2009;<&#x2009;95%; dDDH&#x2009;<&#x2009;70%), the B64 isolate is identified as a novel species. Physiological assays confirmed nitrogen fixation (1.97 ppm), IAA production (3.67 ppm), and phosphate solubilization (26.10 ppm). Structural validation showed 100% of NodD1 residues in allowed regions, ensuring high model reliability. Docking simulations demonstrated strong binding affinities across all flavonoids, with binding free energies ranging from -&#x2009;8.8 to -&#x2009;9.0&#xa0;kcal/mol. Daidzein exhibited the highest thermodynamic stability (-&#x2009;9.0&#xa0;kcal/mol), whereas apigenin showed the most extensive residue interaction network. The B64 isolate is a novel Bradyrhizobium species with a high symbiotic capacity. The stable NodD1-flavonoid interactions provide a molecular basis for efficient nodulation, positioning B64 as a promising candidate for developing lipo-chitooligosaccharide (LCO)-based biofertilizers.

Bradyrhizobium

Phenotypic and transcriptomic characterization of biallelic RNU2-2 developmental and epileptic encephalopathy.

OBJECTIVE: A significant proportion of individuals with suspected genetic developmental and epileptic encephalopathies (DEEs) remain unsolved following whole genome sequencing (WGS). Here we describe biallelic RNU2-2 variants causing a recently reported, severe, recessive DEE. METHODS: We screened individuals who have received WGS analyses at the Genomic Medicine Centre Karolinska for Rare Diseases for biallelic RNU2-2 variants. Deep phenotyping was performed through reviewing entire medical histories and phenotypic traits were transcribed to their corresponding Human Phenotype Ontology (HPO) term. HPO terms were used to generate pairwise phenotypic similarity scores and assess for significantly shared phenotype enrichment in the RNU2-2 sub-cohort. RNA sequencing analyses were performed in fibroblast and blood tissues to compare splicing events between RNU2-2 individuals and two independent control groups. RESULTS: We identified 14 individuals from nine families with 12 ultra-rare biallelic RNU2-2 variants clustering in the conserved 5' domains. Genotype data from 13 of 14 individuals has been reported previously as part of a larger cohort. All individuals presented with a highly concordant, severe DEE, characterized by severe to profound intellectual disability, inability to walk or communicate, hyperkinesia, and refractory seizures. Infantile spasms and tonic seizures were the predominant seizure types and a Lennox-Gastaut syndrome-like phenotype was common. These individuals had a significantly similar phenotypic signature when compared with 703 individuals with complex pediatric epilepsies (two-sided Monte Carlo permutation test, p&#x2009;=&#x2009;.005). RNA sequencing analyses showed aberrant splicing, with the most pronounced effects in fibroblast tissues in mutually exclusive exon and alternate 3' splice-site events, which were not detectable in blood. SIGNIFICANCE: We present deep phenotyping data and transcriptomic analyses that provide support for rare, 5' clustering biallelic RNU2-2 variants causing this novel, severe DEE. We propose an RNA sequencing methodology on fibroblast tissue for future validation of RNU2-2 variants.

autosomal recessive disease