Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Whole systems research”

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.

1,302 records · Page 2Linked to original sources
PubMed2026

Genomic and food-safety evaluation of Staphylococcus chromogenes in Chinese dairy milk.

Non-aureus staphylococci and mammaliicocci (NASM) cause mastitis and may contaminate milk and dairy products. Milk samples (n&#xa0;=&#xa0;1916) from cows with subclinical or clinical mastitis (SCM and CM, respectively) were collected from 28 large-scale (> 500 lactating cows) Chinese dairy farms. Overall, 999 NASM isolates representing 19 species were identified by MALDI-TOF MS and cpn60 sequencing, with Staphylococcuschromogenes, Mammaliicoccus sciuri and Staphylococcus haemolyticus being most prevalent. Antimicrobial resistance (AMR) was determined with disc diffusion; non-susceptible to penicillin was most common (SCM, 30% and CM, 29%) whereas cefoxitin non-susceptible NASM accounted for 8-10% of isolates; among these, 12.5% carried mecA but none carried mecC. Galleria mellonella was used to assess virulence of 78 strains of S. chromogenes, a dominant species; subsequently, 32 strains, representing higher- and lower-virulence in the Galleria model, were selected for whole-genome sequencing and comparative genomics. S. chromogenes isolates from CM had higher virulence (p&#xa0;<&#xa0;0.05) than those from SCM. The 32 genomes comprised 20 sequence types, indicating high genetic diversity. No robust genomic marker of Galleria virulence phenotype was identified in this selected WGS subset. Acquired resistance genes (n&#xa0;=&#xa0;5) were detected, including a first report of fusC in S. chromogenes; the fusC-positive isolate had an elevated fusidic acid MIC (8&#xa0;mg/L). Although S. chromogenes persisted in milk at 4&#xa0;&#xb0;C, pasteurization (64&#xa0;&#xb0;C for 30&#xa0;min) reduced viable counts to below detection. This study provided new insights into the prevalence, AMR, genomic diversity, and dairy-chain relevance of milk-derived NASM, particularly S. chromogenes. However, the genomic findings were based on an intentionally selected WGS subset and should be interpreted as hypothesis-generating rather than population-representative.

Animals↗
PubMed2026

Whole-exome characterization of host genetic variation in HIV-associated genes across the high-prevalence Mizo population, Northeast India.

BACKGROUND: The Mizoram state of Northeast India has one of the highest HIV prevalence rates in Asia, yet the host genetic factors influencing HIV susceptibility in this Tibeto-Burman population remain uncharacterised. METHODS: We performed whole-exome sequencing using Illumina NovaSeq 6000, mean coverage 100X on 76 HIV-negative Mizo individuals. Variants were called using GATK HaplotypeCaller v4.3 against GRCh38p14, annotated with ANNOVAR, and filtered using hard-quality thresholds (QD&#xa0;&#x2265;&#xa0;2, SOR&#xa0;&#x2264;&#xa0;3, MQ&#xa0;&#x2265;&#xa0;40, DP&#xa0;&#x2265;&#xa0;10, GQ&#xa0;&#x2265;&#xa0;20). The allele frequencies were compared against gnomAD v2.1.1 population databases. Hardy-Weinberg equilibrium was assessed using the Wigginton exact test with Bonferroni correction. RESULTS: Post-quality filtering resulted in 12,011 sample-variants across 2,821 unique positions from 36 HIV-associated loci (33 protein-coding genes, 2 chemokine ligands, and 3 lncRNA targets). Of these, 784 observations (51 unique positions) were high-impact nonsynonymous or loss-of-function variants. ADAR rs2229857 (p.K384R, NM_015840) was the most frequently observed variant (Mizo carrier frequency&#xa0;=&#xa0;0.895; 95% CI: 0.806-0.946). CXCR1 rs16858808 (p.R335C) showed the greatest population enrichment (Mizo carrier frequency&#xa0;=&#xa0;0.197; 95% CI: 0.123-0.300; 7.65-fold carrier-frequency enrichment versus gnomAD South Asian; CADD&#xa0;=&#xa0;15.60). Sixteen of 20 tested variants deviated from Hardy-Weinberg equilibrium after Bonferroni correction (p&#xa0;<&#xa0;0.0025), predominantly showing excess homozygosity consistent with the endogamous Mizo population. The protective variant CCR5-&#x394;32 was absent in all the 76 individuals tested. CONCLUSION: This first whole-exome characterization of HIV host genes in the Mizo population identifies CXCR1 rs16858808 as the most population-enriched functional variant and reveals a pervasive endogamy signature. These findings provide a population-specific genetic framework for future HIV susceptibility studies and ART pharmacogenomics research.

Humans↗
PubMed2026

Divergent evolutionary strategies in spider venoms: A comparative proteomic profiling of four sympatric species from Yunnan.

Spider venoms comprise complex cocktails of bioactive molecules evolved for predation and defense, representing a valuable resource for biological research and pharmaceutical discovery. In this study, we performed a systematic analysis of venom gland extracts from four common spider species indigenous to Yunnan, China: Agelena limbata, Hippasa lycosina, Lycosa grahami, and Sinopoda pengi. Using an integrated transcriptomic and proteomic targeted profiling approach, we successfully annotated 141 distinct toxins. Comparative analysis revealed significant interspecific heterogeneity, suggesting distinct evolutionary trajectories and "weapon system economics." Both A. limbata and L. grahami exhibited a "peptide-dominant" profile anchored by neurotoxic peptides and isomerases, optimized for rapid chemical paralysis. In contrast, S. pengi displayed a distinct "protein-dominant" signature enriched with high-molecular-weight enzymes and CAP superfamily proteins, likely functioning to facilitate tissue degradation and toxin diffusion. Occupying an intermediate position, H. lycosina demonstrated a hybrid composition. These findings suggest that although these species share the same geographical range, their venom systems have undergone divergent evolutionary adaptations driven by specific ecological niches and hunting strategies. This study represents the first systematic proteomic characterization of these venom components, providing a valuable reservoir of molecular candidates while highlighting the bioinformatic nuances of analyzing whole-gland homogenates.

Animals↗
PubMed2026

Single-organ proteomics in Drosophila melanogaster larva.

The combination of genetic accessibility, organ complexity, evolutionary conservation, and cost-efficiency makes Drosophila melanogaster (Dm) a well-known model system for biomedical and fundamental biological research. Proteomic analysis of single organs enables the identification and quantification of proteins expressed in specific organs. This will help to uncover specific biological functions and unique protein profiles that are not detectable in whole-organism analyses. In this study we have isolated single organs form Dm larvae, and we have performed a deep proteomics mapping by following a minimal manipulation preparation procedure. The combined dataset across all organs comprised 9132 identified proteins. As anticipated, principal component analysis (PCA) revealed clear separation between the proteomes of most organs, confirming distinct protein profiles. These findings demonstrate the applicability of the sample preparation strategy for high-resolution proteomic characterization of individual organs in Drosophila. Given the extensive genetic tools available for this model organism, our approach has the potential to open new avenues for proteomic studies in Drosophila melanogaster and any other biological systems where the sample amount is limiting. SIGNIFICANCE STATEMENT: Drosophila melanogaster is a well-known model system for biomedical and fundamental biological research that serves as a valuable in vivo model organism due to its high degree of evolutionary conservation with higher vertebrates, tractable genetics, and logistical efficiency. However, the proteome of Drosophila at single organ level has been elusive to date, due to several factors like low sensitivity of previous generation mass spectrometers and sample preparation procedures, difficult isolation of some organs. In this study we have applied a compilation of advanced methods including minimal sample manipulation together with simple, straightforward and efficient protein extraction and digestion methods. Obtained peptides were minimally handled to be analyzed by applying specific and sensitive nLC methods coupled on-line to state-of-the-art MS/MS system. Altogether, the applied strategy allowed us to get the first single organ study to date for this animal. These datasets represent a significative resource for future genomic, transcriptomic and proteomic studies in Drosophila, as multi-omic integration requires deep proteomics to translate data into functional biochemistry, and serves as a critical bridge and an indispensable standalone resource across the genomic, transcriptomic, and proteomic landscapes.

Animals↗
PubMed2026

Context-dependent functional diversity of dorsomedial posterior parietal neurons revealed by single-unit fMRI mapping during naturalistic viewing.

The dorsomedial posterior parietal cortex (dmPPC) plays an important role in episodic processing by integrating sensory, cognitive, and motor information across distributed brain systems. However, how individual dmPPC neurons participate in large-scale functional organization during naturalistic experience remains poorly understood. To address this question, we combined single-unit electrophysiology and awake fMRI in five rhesus macaques of both sexes viewing identical naturalistic video stimuli. Using single-unit fMRI mapping, we generated whole-brain neuron-BOLD functional maps by correlating individual neuronal activity with voxel-wise fMRI signals across the brain. We found that neuron-BOLD functional maps exhibited strong context-dependent organization, with neurons recorded during the same video context showing substantially greater similarity than neurons recorded during different video conditions. Compared with neuronal spiking activity or critical fMRI frames alone, neuron-BOLD functional maps more robustly captured contextual structure. Despite this shared large-scale organization, a substantial subset of neighboring neurons recorded simultaneously from the same electrode displayed markedly distinct whole-brain association patterns, revealing substantial local functional heterogeneity within the dmPPC. This local heterogeneity was not readily explained by waveform-based putative cell class or by opposing neuronal firing dynamics. In addition, distributed cortical and medial temporal regions exhibited highly context-dependent neuron-BOLD association patterns during naturalistic viewing. Together, these findings demonstrate that dmPPC neurons participate in dynamic and heterogeneous large-scale functional organization during naturalistic episodic processing. More broadly, this study establishes single-unit fMRI mapping as a framework for linking single-neuron activity to distributed whole-brain dynamics across contextual conditions.Significance Statement Using single-unit fMRI mapping, this study examined how individual dorsomedial posterior parietal cortex (dmPPC) neurons relate to large-scale brain activity during naturalistic video viewing in macaque monkeys. We found that neuron-BOLD functional maps exhibit strong context-dependent organization and capture contextual structure more robustly than neuronal spiking activity or fMRI frames alone. Despite this shared organization, a substantial subset of neighboring dmPPC neurons displayed markedly distinct whole-brain association patterns, revealing local functional heterogeneity that was not readily explained by waveform-based putative cell class or opposing firing dynamics. These findings provide insight into how local neuronal populations participate in distributed brain-wide functional organization during naturalistic episodic processing.

Journal Article↗
PubMed2026

Longitudinal whole-genome analysis of bluetongue virus identifies conserved serotype-specific genomes and distinct genomic constellations within a Colorado sheep flock (2021-2023).

Bluetongue virus (BTV) is a segmented double-stranded RNA virus of ruminants transmitted by Culicoides spp. biting midges. Although the genome consists of ten segments, classification into serotypes is primarily based on genome segment 2. However, reassortment among genomic segments is a major driver of BTV evolution and diversity. This study used longitudinal whole-genome sequencing to characterize BTV genomes collected from 2021 to 2023 within a single sheep flock in Colorado, where multiple serotypes co-circulate. Whole-genome sequences were generated from fourteen blood samples representing four serotypes: BTV-6, -11, -13, and -17. Longitudinal sampling identified multiple BTV serotypes within individual sheep across consecutive years. Tanglegram analysis comparing segment phylogenies to the segment 2 tree demonstrated incongruent topologies across all genomic segments, suggestive of reassortment or the circulation of distinct genomic constellations. Nucleotide-level comparisons revealed high sequence homology among same-serotype samples from the same year, while the greatest genetic divergence was observed among BTV-17 genomes collected in different years. Additionally, all BTV-13 genomes contained a previously undescribed nonsynonymous substitution in segment 10 predicted to extend the encoded protein by three amino acids. Together, these findings demonstrate that highly conserved BTV genomes and distinct genomic constellations can be detected at the flock level across multiple years. This longitudinal whole-genome approach reveals the genetic complexity of endemic BTV populations, including novel variants and genomic patterns consistent with reassortment that are lost with conventional serotyped-based approaches, highlighting the need to integrate whole-genome characterization into endemic BTV monitoring programs.

Animals↗
PubMed2026

The cold case of state transition 7 (stt7) mutants of Chlamydomonas reinhardtii, solved by whole-genome sequencing.

The process of State Transitions (ST) corresponds to an STT7 kinase-driven redistribution of the transmembrane LHCII antenna proteins between Photosystem II (PSII) and Photosystem I (PSI), which results from changes in their phosphorylation state. For the past two decades, two LHCII-kinase mutants, stt7-1 and stt7-9, have been instrumental in the study of STs in Chlamydomonas reinhardtii, the former being a null mutant for the kinase but quasi-sterile in crosses, while the latter, although fertile, has a leaky phenotype. Using long-read sequencing, this study further characterized the genetic lesions of the stt7 mutant strains through whole-genome reconstruction and de novo chromosome assembly. In addition, two new stt7 null mutants were generated, one derived by crosses from the original stt7-1 and one obtained by Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-associated protein 9 (Cas9) technology. This work provides a comprehensive genomic characterization of the original stt7-1 null mutant, revealing extensive chromosomal rearrangements and high levels of aneuploidy, associated with increased cell size and meiotic dysfunction. Reassessment of their physiology and genetic backgrounds highlights the need for caution in interpreting genetic information. We thus produced more reliable null mutants for the LHCII-kinase, amenable to genetic crosses for the study of STs in a variety of genetic backgrounds.

Chlamydomonas reinhardtii↗
PubMed2026

The combined and independent influence of food texture and a 'mindful eating' instruction on eating rate and food intake among Dutch primary schoolchildren.

Eating rate (ER), the amount of food consumed per unit of time, is a key determinant of food intake, with faster ER associated with larger meal size. Research in adults has shown that both sensory properties like texture, and instructions can influence ER and intake. However, their independent and combined effect on ER and intake in children remain poorly understood. This study examined the effects of food texture and a 'mindful eating' instruction on ER and food intake in children. Children (N&#xa0;=&#xa0;73), 38 boys, aged 4-12 years, participated in a 4-week cluster-randomized incomplete cross-over study conducted during regular school lunches. Using a 2&#xd7;3 factorial design, children were exposed to two levels of food texture (softer vs. harder whole-grain buns) and three types of instructions (none, control, 'mindful'). Linear mixed models with repeated measures were used, adjusting for sex and group. As the interaction between food-texture and type of instructions was non-significant, it was excluded from the final models. Results revealed a strong main effect of food texture on ER and food intake (all p&#xa0;<&#xa0;0.001), where compared with softer buns, harder buns were associated with a slower eating rate (&#x394;&#xa0;=&#xa0;-6.67&#xa0;g/min, SE&#xa0;=&#xa0;0.59) and reduced food intake (&#x394;&#xa0;=&#xa0;-70&#xa0;g, SE&#xa0;=&#xa0;8). In contrast, type of instruction had no significant effect on ER or intake (p&#xa0;=&#xa0;0.85). Our findings demonstrate that food texture exerts a dominant influence on children's eating rate and food intake, highlighting texture modification as a potential leverage point for influencing ER and intake in children.

Humans↗
PubMed2026

Whole-Genome Deep Learning Predicts Chemotherapy Response in Colorectal Cancer.

Chemotherapy response in colorectal cancer (CRC) exhibits significant heterogeneity, with current clinical predictors failing to capture complex genomic determinants of resistance. We developed a hybrid deep learning framework integrating convolutional neural networks (CNNs) and bidirectional long short-term memory (BiLSTM) networks to analyze whole-genome somatic mutations, evolutionary conservation, chromatin accessibility, and 3D genome architecture in 2,546 TCGA patients. An attention mechanism identified predictive genomic regions. The model achieved an AUC of 0.92 (95% CI: 0.89-0.94) in cross-validation and 0.88 (95% CI: 0.85-0.91) in independent validation, outperforming clinical models (&#x394;AUC = +0.18, p < 0.001). Key predictors included non-coding variants in TP53, KRAS, and PIK3CA regulatory regions. Triple-positive patients (mutations in all 3 regions) had significantly worse progression-free survival (HR = 4.7, p < 0.001). Our framework enables accurate chemotherapy response prediction and reveals novel non-coding resistance mechanisms, advancing precision oncology in CRC.

Humans↗
PubMed2026

Whole genome sequencing of unusual Hepatitis C virus subtypes and drug resistance analysis during direct-acting antiviral therapy in India.

INTRODUCTION AND OBJECTIVES: Pangenotypic direct-acting antivirals (DAA) are effective against highly prevalent Hepatitis C virus (HCV) subtypes, but have been clinically validated almost exclusively in high-income countries. Unusual HCV subtypes may carry natural polymorphisms, potentially impacting DAA susceptibility. We conducted full-genome characterization and resistance analysis of unusual HCV subtypes in patients receiving DAA treatment. PATIENTS AND METHODS: In this prospective hospital-based study, eligible patients were screened for anti-HCV antibodies and active infection was confirmed by diagnostic 5'NCR-based HCV RNA detection. Genotyping was performed by core region sequencing, and viral load quantified by real-time PCR. For whole genome sequencing, multiplex primers were designed using alignments of global reference sequences. Sequencing was carried out using the Oxford Nanopore Technology platform. Phylogenetic analysis used multiple sequence alignment and the HCV-GLUE resource for resistance-associated substitution (RAS) analysis. RESULTS: Predominant genotype was genotype 3 in 64.3% (n = 45); genotype 6 in 21.4% (n = 15); and genotype 1 in 14.2% (n = 10). Unusual HCV subtype 6xa was detected in two patients and showed no NS5A resistance mutations. One genotype 3b patient relapsed at 24 weeks post-DAA treatment completion and carried NS5A resistance-associated substitutions 30 K and 31 M both at baseline and at relapse, conferring high-level resistance to NS5A inhibitors. CONCLUSION: This is the first report from India of whole genome sequencing of HCV subtype 6xa. The identification of NS5A resistance mutations in the 3b relapse case underscores challenges for global HCV elimination strategies.

Humans↗
PubMed2026

Evaluation of one-step amplicon-based targeted enrichment for SARS-CoV-2 whole-genome sequencing using the Midnight amplicon scheme.

Genomic surveillance proved invaluable during the COVID-19 pandemic for tracking SARS-CoV-2 variants and guiding outbreak responses, underscoring the ongoing need to reduce whole-genome sequencing (WGS) costs and improve workflow efficiency to ensure accessibility in resource limited settings. Here, we evaluated a one-step reverse transcription polymerase chain reaction (RT-PCR) approach using the Midnight V2 primer scheme for targeted amplification of the SARS-CoV-2 genome, assessed its compatibility with Illumina sequencing, and compared its performance to a well-established two-step method. Initially, we determined optimal RT-PCR reaction conditions using the Midnight V2 primer panel for the one-step RT-PCR kit and scaled reaction volumes for both RT-PCR and library preparation. Clinical specimens (n&#x202f;=&#x202f;53) that had undergone routine WGS for surveillance purposes using the established two-step RT-PCR method were compared using the one-step RT-PCR assay. For samples with genome completeness greater than 70%, both methods gave comparable results with similar sequence coverage and 100% concordance for lineage assignment. Further investigation revealed a higher percentage of reads aligning to the SARS-CoV-2 genome with a greater depth of coverage using the one-step method compared to the two-step method. Finally, analysis of scaled one-step and library reaction volumes revealed significant cost savings for samples undergoing WGS. Overall, the results presented here verify the accuracy and reproducibility of one-step targeted amplification and offer an efficient and cost-effective workflow for routine SARS-CoV-2 genomic surveillance.

Humans↗
PubMed2026

Routine methods misidentify Serratia spp.: Limitations of MALDI-TOF MS revealed by whole-genome sequencing.

Accurate species-level identification within the genus Serratia remains challenging due to extensive phenotypic overlap and high genomic relatedness among closely related and recently described taxa. This study presents an evaluation of routine and genome-based identification approaches applied to clinical Serratia isolates, integrating phenotypic assays, MALDI-TOF MS (Bruker Daltonics), 16S rRNA gene sequencing, and Whole-Genome Sequencing (WGS). A total of 103 isolates collected from a teaching hospital were analyzed. WGS was performed on a subset of isolates. Conventional biochemical methods classified all isolates as Serratia marcescens, whereas MALDI-TOF MS identified 60.1% as S. marcescens, 11.6% as S. ureilytica, and 28.1% just at the genus level. Peak analysis from MALDI-TOF MS revealed specific peaks associated with S. marcescens and S. ureilytica, but limited discriminatory power. WGS of six isolates initially identified as S. ureilytica by MALDI-TOF MS revealed reclassification as Serratia sarumanii (n = 5) and Serratia montpellierensis (n = 1), supported by Average Nucleotide Identity (ANI), Average Amino Acid Identity (AAI), and Digital DNA-DNA Hybridization (dDDH) thresholds. In contrast, 16S rRNA analysis showed limited species-level resolution. Phylogenomic and SNP-based analyses confirmed these classifications with strong support. Overall, this study underscores the critical role of high-resolution genomic approaches for precise species identification and highlights the need for continuous expansion and curation of MALDI-TOF MS reference databases to support reliable clinical diagnostics and epidemiological surveillance of emerging Serratia species.

Spectrometry, Mass, Matrix-Assisted Laser Desorpti↗
PubMed2026

Whole-Exome Sequencing in a Consanguinity-Enriched South Indian Retinitis Pigmentosa Cohort: Diagnostic Yield and Molecular Spectrum.

PURPOSE: To determine the molecular diagnostic yield, variant spectrum, inheritance architecture, and influence of consanguinity on whole-exome sequencing outcomes in a South Indian retinitis pigmentosa (RP) cohort. DESIGN: Prospective, registry-based cohort study. SUBJECTS: A total of 113 affected participants were enrolled through the Aravind Registry for Inherited Diseases of the Eye, including 109 unrelated probands and 4 affected relatives from already represented families. Primary analyses were restricted to the 109 unrelated probands. METHODS: Whole-exome sequencing was performed using a clinical exome workflow. Variants were interpreted using American College of Medical Genetics and Genomics/Association for Molecular Pathology criteria and cases were categorized as solved, possibly solved, inconclusive, or unsolved using prespecified inheritance-aware rules. MAIN OUTCOME MEASURES: Molecular diagnostic yield, distribution of implicated genes and variant classes, inheritance architecture, and diagnostic yield stratified by consanguinity status. RESULTS: Among the 109 unrelated probands, mean age at testing was 39.3 &#xb1; 14.1 years and 58.7% were male. Whole-exome sequencing identified 186 distinct rare variants across 92 inherited retinal disease genes, including 26 pathogenic and 33 likely pathogenic variants. A molecular diagnosis was established in 50 of 109 probands (45.9%), including 42 solved and 8 possibly solved cases; 45 (41.3%) were inconclusive and 14 (12.8%) remained unsolved, including 4 (3.7%) in whom no candidate variant was identified. EYS, USH2A, and ADGRV1 were the most frequently implicated genes. Autosomal recessive (AR) disease predominated (44/50, 88.0%). Consanguineous AR cases were exclusively homozygous (17/17); notably, 68.0% of nonconsanguineous AR cases were also homozygous (P = 0.013). Diagnostic yield was higher in consanguineous probands (51.4% vs. 41.7%), without reaching significance. Recurrent alleles included an established South Asian founder variant (MFSD8 c.1361T>C) and candidate founder alleles in EYS (c.4321C>T) and ADGRV1 (c.14329C>T). CONCLUSIONS: Whole-exome sequencing established a molecular diagnosis in nearly half of this South Indian RP cohort and revealed a predominantly recessive, homozygosity-enriched architecture shaped by consanguinity. These findings define a region-specific variant landscape to support clinical interpretation, genetic counseling, and future trial enrollment in this underrepresented population. FINANCIAL DISCLOSURES: The authors have no proprietary or commercial interest in any materials discussed in this article.

Consanguinity↗
PubMed2026

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↗