Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “WGS”

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.

At least 127 records · Page 7Linked to original sources

Accelerated expansion of group IID-like phospholipase A2 genes in Bos taurus.

Low-molecular-weight, calcium-dependent phospholipase A2 genes (PLA2s) that belong to the secreted type of PLA2s are clustered within a syntenic group on human 1p35-p36 and mouse 4qD3. We reassembled trace files available from the Whole Genome Sequencing (WGS) Project, obtaining an 86-kb contig with three tandem PLA2G2D duplications in the Hereford strain. We used mate-pair data to monitor the assembly and to exclude chimeric clones, demonstrating that the current WGS data may be assembled even in a highly repetitive region with a coverage exceeding fivefold. The genomic structure indicated that most of the PLA2G2D transcripts are formed by four exons. Two alternative first exons were present in all duplications. In two duplications insertions of satellite DNA in the third intron created a novel exon that gave rise to a two-exon product. Linkage and comparative mapping placed the bovine PLA2G2 locus on BTA2, indicating that it evolved from an ancestral PLA2G2D locus common to human, cattle, and rodents. Bovine PLA2G2D variants were capable of encoding 147-amino-acid polypeptides that consisted of putative signal peptide and metal-binding domains. Cysteine residues were conserved in positions analogous to those forming the seven disulfide bonds characteristic of PLA2G2 genes. Quantitative PCR analysis of bovine PLA2G2D transcripts indicated that their expression levels varied between the dry period and lactation in the mammary gland samples and that their expression was polymorphic in liver tissue. The recent burst of duplication and divergence of the bovine PLA2G2D genes and their polymorphic nature are typical of innate immune response genes.

Amino Acid Sequence↗

Gait training efficacy using a home-based practice model in chronic hemiplegia.

OBJECTIVE: The efficacy of a home practice model for gait training was evaluated in 18 hemiplegic subjects 2.3 years (range, 1 to 5) after stroke. DESIGN: Uncontrolled case series. SETTING: Referral center. SUBJECTS: Patients at least 1 year poststroke referred to an outpatient rehabilitation program. INTERVENTION: Patients were taught home programs in two or more 2- to 5-day blocks averaging 35 physical therapy (PT) contact hours (range, 9.5 to 62.5); training extended over a mean of 22 months (range, 10 to 65). Training emphasized weight bearing, balance, segmental control, stretching, and bracing. MAIN OUTCOME MEASURES: Gait changes were measured using the newly developed Wisconsin Gait Scale (WGS). The patient-rated Falls Efficacy Scale (FES) was administered before and after training to 8 subjects, and the Health Status Questionnaire (HSQ) was retrospectively administered to all subjects to appraise subjective pretraining to posttraining changes and current psychological status. RESULTS: The average WGS score significantly improved (p < .05). Patients perceived that gait training increased the quality of their functional activities (p < .05). In a subset of patients, the FES showed that fear of falling was decreased (p < .05). Perception of well-being was comparable to a normative nonstroke reference population except for physical functioning. Compared to the only other published series (using traditional outpatient programming), the current model was of comparable cost. CONCLUSION: Despite the literature indicating a plateau in mobility function by 6 months after stroke, postacute training of gait in hemiplegic subjects using a home-based training model results in improved gait and the perception of improved function. Additionally, we provide validation for the newly developed Wisconsin Gait Scale, an instrument of gait measurement that may assist in comparing outcomes.

Activities of Daily Living↗

Utilizing whole genome sequencing to characterize central line-associated bloodstream infections due to Staphylococcus epidermidis.

Whole genome sequencing (WGS) and clinical review were used to characterize 14 cases of central line-associated bloodstream infection (CLABSI) due to Staphylococcus epidermidis. WGS, which demonstrated disparate strains, suggested that 42.9% of S. epidermidis CLABSI cases were due to contamination, while clinical review suggested that 57.1% were contamination events.

Humans↗

Evaluation of swabbing methods for culture and non-culture-based recovery of multidrug-resistant organisms from environmental surfaces.

OBJECTIVES: Sponge-Sticks (SS) and ESwabs are frequently utilized for detection of multidrug-resistant organisms (MDROs) in the environment. Head-to-head comparisons of SS and ESwabs across recovery endpoints are limited. DESIGN: We compared MDRO culture and non-culture-based recovery from (1) ESwabs, (2) cellulose-containing SS (CS), and (3)&#xa0;polyurethane-containing SS (PCS). METHODS: Known quantities of each MDRO were pipetted on a stainless-steel surface and swabbed by each method. Samples were processed, cultured, and underwent colony counting. DNA was extracted from sample eluates, quantified, and underwent metagenomic next-generation sequencing (mNGS). MDROs underwent whole genome sequencing (WGS). MDRO recovery from paired patient perirectal and PCS-collected environmental samples from clinical studies was determined. SETTING: Laboratory experiment, tertiary medical center, and long-term acute care facility. RESULTS: Culture-based recovery varied across MDRO taxa, it was highest for vancomycin-resistant Enterococcus and lowest for carbapenem-resistant Pseudomonas aeruginosa (CRPA). Culture-based recovery was significantly higher for SS compared to ESwabs except for CRPA, where all methods performed poorly. Nucleic acid recovery varied across methods and MDRO taxa. Integrated WGS and mNGS analysis resulted in successful detection of antimicrobial resistance genes, construction of high-quality metagenome-assembled genomes, and detection of MDRO genomes in environmental metagenomes across methods. In paired patient and environmental samples, multidrug-resistant Pseudomonas aeruginosa (MDRP) environmental recovery was notably poor (0/123), despite detection of MDRP in patient samples (20/123). CONCLUSIONS: Our findings support the use of SS for the recovery of MDROs. Pitfalls of each method should be noted. Method selection should be driven by MDRO target and desired endpoint.

Humans↗

Density functional theory of water-gas shift reaction on molybdenum carbide.

The density functional theory (DFT) of the water-gas shift (WGS) reaction over molybdenum carbide was studied with the aim of understanding the dissociation of H(2)O, the OH group, and CO to determine on what sections of molybdenum carbide CO(2) and H(2) formed and whether they played a role in the reaction. The energy diagram of each elementary step, the reaction of the hydrogen and oxygen atoms with CO, and the transition state for this elementary step were also studied. The IR spectra of the CO adsorption was experimentally analyzed for the identification of several candidates of the CO adsorption modes. The adsorptions of the threefold Mo site (a) with and (b) without the underlying C atom of the second layer have the second and highest adsorption energies of -281.59 and -321.00 kJ/mol, respectively. The IR data showed that the bands at 1626 cm(-1) from the IR experiments are (a) the nearest adsorption of the threefold Mo site with the underlying C atom at the calculated/corrected band of 1621 cm(-1). The calculated/corrected threefold adsorption (b) had the highest adsorption energy but exhibited an IR band at 1147 cm(-1) which was not observed in the experimental data. The C-O bond length increased to 1.49 from 1.36 after the H(2)O adsorption (b), suggesting the dissociation of C-O after the H(2)O coadsorption. The WGS reaction on the beta-Mo(2)C(001) slab carbide was calculated and took place as follows: H(2)O was dissociated into OH and H on the Mo(2)C surface and the OH subsequently dissociated into H and O atoms. CO approached the O atom to form CO(2).

Adsorption↗

Kinetic and equilibrium study on formic acid decomposition in relation to the water-gas-shift reaction.

Kinetics and equilibrium are studied on the hydrothermal decarbonylation and decarboxylation of formic acid, the intermediate of the water-gas-shift (WGS) reaction, in hot water at temperatures of 170-330 degrees C, to understand and control the hydrothermal WGS reaction. (1)H and (13)C NMR spectroscopy is applied to analyze as a function of time the quenched reaction mixtures in both the liquid and gas phases. Only the decarbonylation is catalyzed by HCl, and the reaction is first-order with respect to both [H(+)] and [HCOOH]. Consequently, the reaction without HCl is first and a half (1.5) order due to the unsuppressed ionization of formic acid. The HCl-accelerated decarbonylation path can thus be separated in time from the decarboxylation. The rate and equilibrium constants for the decarbonylation are determined separately by using the Henry constant (gas solubility data) for carbon monoxide in hot water. The rate constant for the decarbonylation is 1.5 x 10(-5), 2.0 x 10(-4), 3.7 x 10(-3), and 6.3 x 10(-2) mol(-1) kg s(-1), respectively, at 170, 200, 240, and 280 degrees C on the liquid branch of the saturation curve. The Arrhenius plot of the decarbonylation is linear and gives the activation energy as 146 +/- 3 kJ mol(-1). The equilibrium constant K(CO) = [CO]/[HCOOH] is 0.15, 0.33, 0.80, and 4.2, respectively, at 170, 200, 240, and 280 degrees C. The van't Hoff plot results in the enthalpy change of DeltaH = 58 +/- 6 kJ mol(-1). The decarboxylation rate is also measured at 240-330 degrees C in both acidic and basic conditions. The rate is weakly dependent on the solution pH and is of the order of 10(-4) mol kg(-1) s(-1) at 330 degrees C. Furthermore, the equilibrium constant K(CO2) = [CO(2)][H(2)]/[HCOOH] is estimated to be 1.0 x10(2) mol kg(-1) at 330 degrees C.

Journal Article↗

Long-read sequencing reveals a hidden Alu-mediated splice defect in CPLANE1, causing orofaciodigital syndrome type VI.

Orofaciodigital syndrome type VI (OFD VI) is a recessive ciliopathy characterized by excessive polydactyly, molar tooth sign, cleft lip, and developmental delay, caused by pathogenic variants in CPLANE1. Here, we present a patient with OFD VI that remained genetically unexplained after routine genetic testing, including short-read whole genome sequencing (WGS). Using long-read sequencing, we found two biallelic splice-site variants in CPLANE1, c.8633-4_8633-3del, and an Alu element insertion close to an exon-intron boundary. Transcript analysis showed that each variant independently resulted in exon skipping, and quantitative expression studies revealed reduced total CPLANE1 mRNA levels in patient-derived fibroblasts. Based on these findings, we were able to re-classify the c.8633-4_8633-3del variant from a variant of uncertain significance (VUS) to likely pathogenic. The identification of an Alu element insertion missed by short-read WGS highlights the added diagnostic value of long-read sequencing in uncovering cryptic, transposable element-associated pathogenic variants.

Journal Article↗

Whole-genome sequencing reveals progressive versus stable myeloma precursor conditions as two distinct entities.

Multiple myeloma (MM) is consistently preceded by precursor conditions recognized clinically as monoclonal gammopathy of undetermined significance (MGUS) or smoldering myeloma (SMM). We interrogate the whole genome sequence (WGS) profile of 18 MGUS and compare them with those from 14 SMMs and 80 MMs. We show that cases with a non-progressing, clinically stable myeloma precursor condition (n&#x2009;=&#x2009;15) are characterized by later initiation in the patient's life and by the absence of myeloma defining genomic events including: chromothripsis, templated insertions, mutations in driver genes, aneuploidy, and canonical APOBEC mutational activity. This data provides evidence that WGS can be used to recognize two biologically and clinically distinct myeloma precursor entities that are either progressive or stable.

DNA Copy Number Variations↗

Rare variant contribution to the heritability of coronary artery disease.

Whole genome sequences (WGS) enable discovery of rare variants which may contribute to missing heritability of coronary artery disease (CAD). To measure their contribution, we apply the GREML-LDMS-I approach to WGS of 4949 cases and 17,494 controls of European ancestry from the NHLBI TOPMed program. We estimate CAD heritability at 34.3% assuming a prevalence of 8.2%. Ultra-rare (minor allele frequency &#x2264;&#x2009;0.1%) variants with low linkage disequilibrium (LD) score contribute ~50% of the heritability. We also investigate CAD heritability enrichment using a diverse set of functional annotations: i) constraint; ii) predicted protein-altering impact; iii) cis-regulatory elements from a cell-specific chromatin atlas of the human coronary; and iv) annotation principal components representing a wide range of functional processes. We observe marked enrichment of CAD heritability for most functional annotations. These results reveal the predominant role of ultra-rare variants in low LD on the heritability of CAD. Moreover, they highlight several functional processes including cell type-specific regulatory mechanisms as key drivers of CAD genetic risk.

Humans↗

Rare variant analyses in 51,256 type 2 diabetes cases and 370,487 controls reveal the pathogenicity spectrum of monogenic diabetes genes.

Type 2 diabetes (T2D) genome-wide association studies (GWASs) often overlook rare variants as a result of previous imputation panels' limitations and scarce whole-genome sequencing (WGS) data. We used TOPMed imputation and WGS to conduct the largest T2D GWAS meta-analysis involving 51,256 cases of T2D and 370,487 controls, targeting variants with a minor allele frequency as low as 5&#x2009;&#xd7;&#x2009;10-5. We identified 12 new variants, including a rare African/African American-enriched enhancer variant near the LEP gene (rs147287548), associated with fourfold increased T2D risk. We also identified a rare missense variant in HNF4A (p.Arg114Trp), associated with eightfold increased T2D risk, previously reported in maturity-onset diabetes of the young with reduced penetrance, but observed here in a T2D GWAS. We further leveraged these data to analyze 1,634 ClinVar variants in 22 genes related to monogenic diabetes, identifying two additional rare variants in HNF1A and GCK associated with fivefold and eightfold increased T2D risk, respectively, the effects of which were modified by the individual's polygenic risk score. For 21% of the variants with conflicting interpretations or uncertain significance in ClinVar, we provided support of being benign based on their lack of association with T2D. Our work provides a framework for using rare variant GWASs to identify large-effect variants and assess variant pathogenicity in monogenic diabetes genes.

Diabetes Mellitus, Type 2↗

Prevalence and chronology of colibactin-associated mutational processes and their microbiome spectra in Japanese colorectal cancer.

The incidence of colorectal cancer (CRC) has risen in recent decades, with a disproportionate increase observed among younger individuals in Japan and other countries. The etiological contribution of the gut microbiota to CRC pathogenesis is recognized, yet the mechanisms involved remain to be fully clarified. Here we integrated whole-genome sequencing (WGS) and transcriptome profiling of CRC with whole-genome metagenomic sequencing of fecal samples to interrogate host-microbiome interactions at high resolution. Application of interpretable artificial intelligence enabled the stratification of CRC into four distinct microbiome-informed subtypes. WGS analysis identified mutational signatures SBS88 and ID18, linked to colibactin exposure, as early clonal events detected in 44.8% of non-hypermutated patients. Notably, these signatures were significantly more frequent among patients born after the 1960s. Microbiome-based subclassification revealed subtype-specific clinical and molecular features. Collectively, our findings indicate that colibactin exposure constitutes a prevalent and potentially modifiable risk factor for CRC in the Japanese population.

Humans↗

Development and extensive sequencing of a broadly-consented Genome in a Bottle matched tumor-normal pair.

The Genome in a Bottle Consortium (GIAB), hosted by the National Institute of Standards and Technology (NIST), is developing new matched tumor-normal samples, the first explicitly consented for public dissemination of genomic data and cell lines. Here, we describe a comprehensive genomic dataset from the first individual, HG008, including DNA from an adherent, epithelial-like pancreatic ductal adenocarcinoma (PDAC) tumor cell line and matched normal cells from duodenal and pancreatic tissues. Data for the tumor-normal matched samples comes from seventeen distinct state-of-the-art whole genome measurement technologies, including high depth short and long-read bulk whole genome sequencing (WGS), single cell WGS, Hi-C, and karyotyping. These data will be used by the GIAB Consortium to develop matched tumor-normal benchmarks for somatic variant detection. We expect these data to facilitate innovation for whole genome measurement technologies, de novo assembly of tumor and normal genomes, and bioinformatic tools to identify small and structural somatic variants. This first-of-its-kind broadly consented open-access resource will facilitate further understanding of sequencing methods used for cancer biology.

Humans↗

Low-pass whole-genome sequencing reveals genomic diversity and ecotype-specific adaptation in indigenous Tigrayan chickens.

Indigenous chickens play a critical role in food security and climate resilience in smallholder systems, yet their genomic diversity and adaptive potential remain insufficiently characterised. This study employed low-pass whole-genome sequencing (LP-WGS; 0.2-1.99&#xd7;) to investigate genomic diversity, population structure, inbreeding and candidate environment-associated genomic variation in 33 chickens from highland, midland, and lowland agroecologies in the Tigray region of northern Ethiopia. After imputation and stringent filtering, 23.4 million high-confidence SNPs were retained, including&#x2009;~&#x2009;17% novel variants, indicating substantial uncharacterised genetic diversity in these populations. SNP density (13.8&#x2009;&#xb1;&#x2009;8.6 SNPs/kb) was comparable to values reported from high-coverage Ethiopian chicken datasets, demonstrating the suitability of LP-WGS for population genomics in resource-limited settings. Marked differences in genomic diversity were observed among ecotypes: midland chickens showed the highest nucleotide diversity (&#x3c0;&#x2009;=&#x2009;0.00267), followed by lowland (&#x3c0;&#x2009;=&#x2009;0.00233), whereas highland chickens showed the lowest diversity (&#x3c0;&#x2009;=&#x2009;0.00203) and elevated genomic inbreeding (FROH and FHOM &#x2248; 0.18). Population structure analyses revealed clear genetic separation among ecotypes. PCA (13.91% variation explained) distinguished lowland chickens along PC1 and separated highland from midland along PC2, while ADMIXTURE and FST patterns supported three major ancestral genomic backgrounds. Functional annotation of private missense variants uncovered distinct adaptive signatures reflecting the contrasting agroecological conditions. Highland chickens showed enrichment of candidate genes potentially involved in physiological processes relevant to high-altitude environments, including cold response, angiogenesis, cardiovascular regulation and metabolic homeostasis (eg., PARP1, ACOX2, ITGB3, EDNRB, SOX8, and SOX10). Midland chickens exhibited candidate signals of selection in genes with known roles in innate antiviral immunity, bacterial defence and inflammatory regulation (eg., BAK1, CLSTN1, CYSLTR1, CYSLTR2, CXCR7, GIPR, DSCAM, GDAP1, TLR3, TLR4, TLR7, IFIH1, ADORA1, EPHB1, and TMPRSS2). Lowland chickens displayed candidate variants associated with heat-stress response, DNA damage repair, oxidative balance and cardiovascular support under extreme temperatures (e.g., MLH1, BDKRB1, GPR19, FLT1, CCL18, TGM2, and RAMP3). Overall, the results indicate substantial genomic differentiation among ecotypes and suggest candidate environment-associated genetic divergence across Tigray's diverse agroecological zones. These populations may represent important reservoirs of adaptive genetic variation for climate-resilient poultry breeding, warranting further functional validation and conservation-oriented management.

Animals↗

Whole-genome sequencing reveals hidden antimicrobial resistance genes in phenotypically susceptible probiotic candidate lactic acid bacteria.

Phenotypic assays commonly used to evaluate probiotic safety may fail to detect clinically relevant antimicrobial resistance (AMR), potentially allowing genetically concerning strains to appear acceptable based on MIC testing alone. To explore this issue, we applied whole-genome sequencing (WGS) to three lactic acid bacteria (LAB) isolates previously identified as probiotic candidates based on acid and bile tolerance, antagonism against enteric pathogens, and biofilm formation in vitro: Lactiplantibacillus plantarum L25F and L22F (from pigs) and Ligilactobacillus salivarius AF2319 (from a chicken). Genome annotation identified extensive repertoires of probiotic-associated genes (46-47 per strain) linked to stress tolerance, adhesion, immunomodulation, and quorum sensing, supporting functional potential. The two L. plantarum strains exhibited broader predicted metabolic capacities than L. salivarius AF2319. However, genomic analysis revealed acquired AMR genes with complex genotype-phenotype relationships not fully apparent from phenotypic testing. The L. plantarum strains harbored lnu(A) (99.79% identity) on extrachromosomal DNA, conferring the L-phenotype (lincomycin resistance, clindamycin susceptibility); clindamycin MICs (1&#xa0;mg/L) were concordant with this genotype, though lincomycin MICs were not determined. L. salivarius AF2319 carried tet(M), tet(L), and erm(C) (99.48%, 99.49%, and 99.45% identity by ResFinder, respectively) on extrachromosomal DNA; notably, the erythromycin MIC (1&#xa0;mg/L) was precisely at the EFSA breakpoint (&#x2264;&#x2009;1&#xa0;mg/L), representing borderline genotype-phenotype discordance potentially due to silent gene expression. Under current EFSA QPS criteria, these acquired ARGs would preclude all three strains from approval as probiotic feed additives despite favorable functional profiles, underscoring the indispensable role of WGS-based AMR gene detection in modern probiotic safety evaluation.

Probiotics↗

Genomic and phenotypic characterization of Klebsiella pneumoniae phage KP &#xd8;1: a novel lytic Slopekvirus targeting uropathogenic multidrug-resistant Klebsiella pneumoniae.

The rise of multidrug-resistant (MDR) uropathogenic gram-negative bacteria (GNB) necessitates the development of alternative therapeutic strategies. This study aimed to isolate, phenotypically characterize, and perform whole-genome sequencing of the bacteriophage demonstrating the broadest host range against MDR uropathogens. Fifty MDR GNB isolates were screened for lytic phages. The most promising candidate, Klebsiella pneumoniae phage KP &#xd8;1, was characterized using plaque assay, Transmission Electron Microscopy (TEM), and pH/thermal stability testing. Genomic characterization was performed via whole-genome sequencing (WGS), with functional annotation and lifestyle prediction using PhaBOX and PhageScope software. Klebsiella pneumoniae was the most prevalent MDR uropathogen. Klebsiella pneumoniae phage KP &#xd8;1 exhibited a 50% host range and high lytic titer (10&#x2078; PFU/mL). TEM revealed an icosahedral head and short contractile tail. Genomic characterization by WGS revealed that Klebsiella pneumoniae phage KP &#xd8;1 possesses a 174,591&#xa0;bp double-stranded deoxyribonucleic acid (dsDNA) genome containing 274 predicted open reading frames (ORFs). No lysogeny-related genes, toxins, or antibiotic resistance markers were detected, confirming its strictly lytic nature and supporting its potential as a candidate for phage therapy applications. The phage remained stable (10&#x2078; PFU/mL) across temperatures of -&#x2009;20&#xa0;&#xb0;C to 50&#xa0;&#xb0;C; supporting its suitability for long-term biobanking and suggesting potential activity at physiological temperature, and across a pH range of 7-9. Klebsiella pneumoniae phage KP &#xd8;1 is a novel, obligately lytic Slopekvirus whose genomic architecture, stability profile, and absence of lysogeny-associated, virulence, and antimicrobial resistance genes ( AMR) collectively support its candidacy for further preclinical evaluation as a phage therapy agent against uropathogenic MDR Klebsiella pneumoniae.

Klebsiella pneumoniae↗

GFPT1 as a cross-ancestry validated target for degenerative spinal disease: genetic association in a Chinese cohort and functional characterization in zebrafish.

Degenerative spinal disease (DSD), including spinal stenosis and spondylosis, lacks effective pharmacological treatment. To identify druggable targets and assess cross-ancestry applicability, we integrate multi-omics analyses using Summary-data-based Mendelian Randomization (SMR), colocalization, and two-sample Mendelian randomization with European whole-blood, peripheral-blood, and CSF eQTL/pQTL datasets, followed by whole-genome sequencing (WGS) validation in a Chinese cohort. We identify 7 genes/proteins associated with spinal stenosis and 5 with spondylosis, with GFPT1, GPX1, and SERPINA1 shared by both. Two-sample MR further supports the causal associations of these targets with DSD. Phenome-wide MR prioritization selects GFPT1 and GPX1 as favorable candidates with no predicted adverse effects and potential beneficial effects on hypertension. In the Chinese cohort (67 lumbar spinal stenosis patients and 100 controls), WGS identifies 4 GFPT1 cis-eQTL loci (rs13016371, rs35392088, rs12997521, and rs13019789) associated with lumbar spinal stenosis risk; all risk alleles are linked to increased GFPT1 expression, and all 24 variant carriers show L4/L5 stenosis on imaging. Druggability analysis identifies IOX1 as the sole preclinical-stage compound targeting GFPT1, and molecular docking supports robust binding to GFPT1 (-&#x2009;6.39&#x2009;kcal/mol). Functional assays show that IOX1 directly inhibits GFPT1 enzymatic activity and induces fructose-6-phosphate accumulation. In zebrafish, IOX1 significantly rescues GFPT1-induced degenerative phenotypes. These findings establish GFPT1 as a cross-ancestry validated therapeutic target for DSD and nominate IOX1 as a promising disease-modifying candidate.

Animals↗

Prediction of antimicrobial minimum inhibitory concentration from bacterial genomes using a scalable and interpretable machine learning approach.

Although machine learning models can predict antimicrobial susceptibility from bacterial whole genome sequencing (WGS), state-of-the-art approaches are computationally demanding or dependent on knowledge of genetic resistance determinants. Here, we describe an efficient data-driven approach to predicting minimum inhibitory concentration (MIC) by progressively extending and refining predictive genome segments, independent of prior knowledge of resistance determinants. Resultant models had high interpretability - known and potentially novel resistance determinants were captured. Using 762 clinical E. coli strains, 71.6% of predictions were within one dilution of the measured MIC. Models trained with this algorithm generalised better onto external data (F1 score&#x2009;=&#x2009;0.85) compared with alternative models trained on annotated resistance determinants (F1&#x2009;=&#x2009;0.82) or k-mer counts (F1&#x2009;=&#x2009;0.74). Computational demands were low (RAM usage 23.6GB vs 38.8GB for k-mer model). These advantages represent an important advance in predicting antimicrobial susceptibility from WGS, with potential applications for clinical diagnostics, drug development, and surveillance.

Journal Article↗

Genomic characterisation and lytic potential of phage SF01 against multidrug-resistant Salmonella enterica subsp. enterica, a key agent of infection in poultry.

1. Salmonella enterica remains the key cause of salmonellosis in poultry, causing high morbidity and mortality. Due to the unprecedented resistance of S. enterica to antibiotics, bacteriophages have emerged as a powerful alternative to conventional antibiotics treatment for salmonellosis.2. In this study, a strain was isolated from infected broiler chickens and whole-genome sequencing (WGS) identified the strain SFD-01 as S. enterica subsp. enterica. Bioinformatics analyses revealed that the genome was 4.6 Mb in size with 4559 coding sequences (CDS), 77 tRNAs and 4 rRNAs. Additionally, 119 virulence genes, 125 antimicrobial resistance genes, 5 mobile genetic elements, 2 prophages and multiple copies of pathogenicity islands (SPI) were identified in the genome.3. To address this, bacteriophage SF01 was isolated from wastewater near a chicken slaughterhouse against host strain SFD-01. Transmission electron microscopy revealed that the phage had an 85-nm icosahedral head and a 130-nm long contractile tail. The Felixounavirus SF01 exhibited high stability across pH 3-9. Phage lytic activity at a multiplicity of infection of 0.01 restricted the bacterial growth.4. Whole genome analysis (WGS) identified phage SF01 as a Felixounavirus with 88-kb genome composed of 174 CDS, 20 tRNA genes and with no lysogenic markers, resistance genes or virulence factors. The strict lytic potential of phage SF01 makes it a highly viable option for use in the potential biocontrol of the novel strain S. enterica subsp. enterica serotype 42:z4,z23.

Felixounavirus↗