Single Nucleotide Polymorphism Analysis Research Articles

Functional potential and evolutionary response to long-term heat selection of bacterial associates of coral photosymbionts.

Symbiotic microorganisms are crucial for the survival of corals and their resistance to coral bleaching in the face of climate change. However, the impact of microbe-microbe interactions on coral functioning is mostly unknown but could be essential factors for coral adaption to future climates. Here, we investigated interactions between cultured dinoflagellates of the Symbiodiniaceae family, essential photosymbionts of corals, and associated bacteria. By assessing the genomic potential of 49 bacteria, we found that they are likely beneficial for Symbiodiniaceae, through the production of B vitamins and antioxidants. Additionally, bacterial genes involved in host-symbiont interactions, such as secretion systems, accumulated mutations following long-term exposure to heat, suggesting symbiotic interactions may change under climate change. This highlights the importance of microbe-microbe interactions in coral functioning.

Molecular characterization of IncFII plasmid carrying blaNDM-5 in a Salmonella enterica serovar Typhimurium ST34 clinical isolate in China.

In this study, an IncFII plasmid pIncFII-NDM5 carrying blaNDM-5 was found in carbapenem-resistant Salmonella enterica serovar Typhimurium (S. enterica serovar Typhimurium), which has conjugative transferability and carried blaNDM-5, bleMBL, mph(A), and blaTEM-1 four resistance genes that can mediate resistance to multiple antibiotics including cephalosporins, beta-lactamase inhibitor combinations, carbapenems, and macrolides. Phylogenetic analysis showed that 1104-65 and 1104-75 were closely related to other S. enterica serovar Typhimurium in this area. The above-mentioned S. enterica serovar Typhimurium chromosome carries blaCTX-M-55, qnrS1, and tet(A) genes, so the antibiotic resistance of isolates will be further enhanced after obtaining the pIncFII_NDM5-like plasmid. Meanwhile, we discovered a novel genetic structure of blaNDM-5 mediated by the IS26 composite transposon, which will expand our understanding of the emergence and spread of carbapenem-resistance genes. Altogether, the presence of the IncFII plasmid pIncFII-NDM5 further underscores the need for vigilant surveillance and appropriate infection control measures to mitigate the impact of carbapenem-resistant S. enterica serovar Typhimurium in clinical settings.

Evaluating the Performance of Multiple Sequence Alignment Programs with Application to Genotyping SARS-CoV-2 in the Saudi Population

This study explores the accuracy and efficiency of multiple sequence alignment (MSA) programs, focusing on ClustalΩ, MAFFT, and MUSCLE in the context of genotyping SARS-CoV-2 for the Saudi population. Our results indicate that MAFFT outperforms the others, making it an ideal choice for large-scale genomic analyses. The comparative performance of MSAs assembled using MergeAlign demonstrates that MAFFT and MUSCLE consistently exhibit higher accuracy than ClustalΩ in both reference-based and consensus-based approaches. The evaluation of genotyping effectiveness reveals that the addition of a reference sequence, such as the SARS-CoV-2 Wuhan-Hu-1 isolate, does not significantly affect the alignment process, suggesting that using consensus sequences derived from individual MSA alignments may yield comparable genotyping outcomes. Investigating single-nucleotide polymorphisms (SNPs) and mutations highlights distinctive features of MSA programs. ClustalΩ and MAFFT show similar counts, while MUSCLE displays the highest SNP count. High-frequency SNP analysis identifies MAFFT as the most accurate MSA program, emphasizing its reliability. Comparisons between Saudi and global SARS-CoV-2 populations underscore regional genetic variations. Saudis exhibit consistently higher frequencies of high-frequency SNPs, attributed to genetic similarity within the population. Transmission dynamics analysis reveals a higher frequency of co-mutations in the Saudi dataset, suggesting shared evolutionary patterns. These findings emphasize the importance of considering regional diversity in genetic analyses.

Genetic polymorphisms of GGT1 gene (rs8135987, rs5751901 and rs2017869) are associated with neoadjuvant chemotherapy efficacy and toxicities in breast cancer patients

BackgroundOur previous study illustrated the predictive value of serum gamma-glutamyl transpeptidase (GGT) for neoadjuvant chemotherapy (NAC) sensitivity in breast cancer patients. In this study we aim to determine whether single nucleotide polymorphisms (SNPs) in the gamma-glutamyltransferase 1 (GGT1) gene are related to the NAC response and adverse events and to find out a genetic marker in predicting NAC sensitivity.MethodsThree SNP loci (rs8135987, rs5751901, rs2017869) of GGT1 gene were selected and tested among breast cancer patients reciving NAC. Four genotype models were used in SNP analysis: co-dominant model compared AA vs. Aa vs. aa; dominant model compared AA vs. Aa + aa; recessive model compared AA + Aa vs. aa; over-dominant model compared AA + aa vs. Aa. Chi-squared test and multivariable logistic regression analysis were performed between SNP genotypes, haplotypes and pathological complete response(pCR), adverse events as well as serum GGT level.ResultsA total of 143 patients were included in the study. For SNP rs8135987 (T > C), the TC genotype in over-dominant model was inversely related with pCR (adjusted OR = 0.30, 95% CI 0.10–0.88, p = 0.029) as well as the risk of peripheral neuropathy (adjusted OR = 0.39, 95% CI 0.15–0.96, p = 0.042). The TC genotype in dominant model was significantly associated with elevated serum GGT level (OR = 3.11, 95% CI 1.07–9.02, p = 0.036). For rs2017869 (G > C), the occurrence of grade 2 or greater neutropenia (OR = 0.39, 95% CI 0.08–0.84, p = 0.025) and leukopenia (OR = 0.24, 95% CI 0.08–0.78, p = 0.017) were both significantly reduced in patients with CC genotypes. For rs5751901(T > C), the CC genotype could significantly reduce the risk of grade 2 or greater neutropenia (OR = 0.29, 95% CI 0.09–0.96, p = 0.036) and leukopenia (OR = 0.27, 95% CI 0.09–0.84, p = 0.024) in recessive model.ConclusionsThe GGT1 gene SNPs might be an independent risk factor for poor response of NAC in breast cancer patients, providng theoretical basis for further precision therapy.

Invasive non-typhoidal Salmonella (iNTS) aminoglycoside-resistant ST313 isolates feature unique pathogenic mechanisms to reach the bloodstream

Invasive non-typhoidal Salmonella (iNTS) from the clonal type ST313 (S. Typhimurium ST313) is the primary cause of invasive salmonellosis in Africa. Recently, in Brazil, iNTS ST313 strains have been isolated from different sources, but there is a lack of understanding of the mechanisms behind how these gut bacteria can break the gut barrier and reach the patient's bloodstream. Here, we compare 13 strains of S. Typhimurium ST313, previously unreported isolates, from human blood cultures, investigating aspects of virulence and mechanisms of resistance. Initially, RNAseq analyses between ST13-blood isolate and SL1344 (ST19) prototype revealed 15 upregulated genes directly related to cellular invasion and replication, such as sopD2, sifB, and pipB. Limited information is available about S. Typhimurium ST313 pathogenesis and epidemiology, especially related to the global distribution of strains. Herein, the correlation of strains isolated from different sources in Brazil was employed to compare clinical and non-clinical isolates, a total of 22 genomes were studied by single nucleotide polymorphism (SNPs). The epidemiological analysis of 22 genomes of S. Typhimurium ST313 strains grouped them into three distinct clusters (A, B, and C) by SNP analysis, where cluster A comprised five, group B six, and group C 11. The 13 clinical blood isolates were all resistant to streptomycin, 92.3% of strains were resistant to ampicillin and 15.39% were resistant to kanamycin. The resistance genes acrA, acrB, mdtK, emrB, emrR, mdsA, and mdsB related to the production of efflux pumps were detected in all (100%) strains studied, similar to pathogenic traits investigated. In conclusion, we evidenced that S. Typhimurium ST313 strains isolated in Brazil have unique epidemiology. The elevated frequencies of virulence genes such as sseJ, sopD2, and pipB are a major concern in these Brazilian isolates, showing a higher pathogenic potential.

Whole genome sequencing-based identification of human tuberculosis caused by animal-lineage Mycobacterium orygis.

A recently described member of the Mycobacterium tuberculosis complex (MTBC) is Mycobacterium orygis, which can cause disease primarily in animals but also in humans. Although M. orygis has been reported from different geographic regions around the world, due to a lack of proper identification techniques, the contribution of this emerging pathogen to the global burden of zoonotic tuberculosis is not fully understood. In the present work, we report single nucleotide polymorphism (SNP) analysis using whole genome sequencing (WGS) that can accurately identify M. orygis and differentiate it from other members of the MTBC species. WGS-based SNP analysis was performed for 61 isolates from different provinces in Canada that were identified as M. orygis. A total of 56 M. orygis sequences from the public databases were also included in the analysis. Several unique SNPs in the gyrB, PPE55, Rv2042c, leuS, mmpL6, and mmpS6 genes were used to determine their effectiveness as genetic markers for the identification of M. orygis. To the best of our knowledge, five of these SNPs, viz., gyrB 277 (A→G), gyrB 1478 (T→C), leuS 1064 (A→T), mmpL6 486 (T→C), and mmpS6 334 (C→G), are reported for the first time in this study. Our results also revealed several SNPs specific to other species within MTBC. The phylogenetic analysis shows that the studied genomes were genetically diverse and clustered with M. orygis sequences of human and animal origin reported from different geographic locations. Therefore, the present study provides a new insight into the high-confidence identification of M. orygis from MTBC species based on WGS data, which can be useful for reference and diagnostic laboratories.

In silico analysis of missense SNPs in GABRA1, GABRB1, and GABRB3 genes associated with some diseases in neurodevelopmental disorders

BackgroundNeurodevelopmental disorders are disorders that are generally seen in the early developmental period of an individual's life and involve more than one disease that causes disruptions in the central nervous system. These disorders can be given as examples of diseases such as autism, mental retardation, some epileptic disorders, communication disorders, and mental retardation. The aim of this study is to determine the possible harmful effects of missense single nucleotide polymorphisms (SNPs) in the GABRA1, GABRB1, and GABRB3 genes, which are associated with neurodevelopmental disorders, on the structure and stabilization of the protein, using in silico methods. Software tools SIFT, PolyPhen-2 HumVar, PolyPhen-2 HumDiv, PROVEAN, SNAP2, PHD-SNP, SNP&GO, PANTHER, and Meta-SNP were used to predict harmful SNPs. I-Mutant and MUpro software tools were used to predict the effects of predicted harmful SNPs on protein stabilization. The STRING software tool was used for protein–protein interactions, the GeneMANIA software tool for gene–gene interactions, and the Project HOPE software tool for three-dimensional modeling examples.ResultsAs a result of the bioinformatics analysis, rs121434579, rs139163545, and rs267600530 in the GABRA1 gene; rs74608570, rs75612351, rs78815529 in the GABRB1 gene, and rs7819600779, rs1719850690, rs7819600779, rs171985060690, rs7819600779, rs1719850600779, rs149963014 in the GABRB3 gene were predicted as harmful SNPs.ConclusionsIn this study, protein structure, function, and stabilization of SNPs known to cause amino acid substitutions in GABRA1, GABRB1, and GABRB3 genes associated with some diseases in neurodevelopmental disorders were investigated using bioinformatics tools. As a result of the results obtained in our study, it is thought that it will benefit experimental studies and bioinformatics studies.

The genetic diversity and population structure of Macridiscus multifarius along the coast of China using IIB-RAD sequencing

Abstract The clam Macridiscus multifarius is ecologically and economically important in the coastal regions of China. In this study, the IIB-RAD method was used to investigate the population diversity and structure of M. multifarius along the coast of China. Genomic DNA was extracted, sequenced, and subjected, resulting in 110.88 Gb of clean reads. Analysis showed that the M. multifarius populations along the coast of China exhibited a moderate to high level of genetic diversity (mean: HO = 0.272, HE = 0.365, π = 0.390, and FIS = 0.268). Further analysis of high-quality single nucleotide polymorphisms was used to generate a maximum-likelihood phylogenetic tree, principal component analysis, and FST values. The results revealed that there is no notable genetic differentiation among the nine M. multifarius populations south of the Yangtze River estuary. However, the Yancheng population shows clear genetic differentiation compared with the other nine populations, possibly linked to the barrier effect of the Yangtze dilution water and ocean current system. We believe that these results provide genetic insights for fisheries management to plan fisheries policies for M. multifarius and conserve natural resources.

Molecular characteristic analysis of single-nucleotide polymorphisms in SLC16A9/hMCT9

AimsHuman monocarboxylate transporter 9 (hMCT9), encoded by SLC16A9, is a transporter that mediates creatine transport across the transmembrane. Previously, we reported that hMCT9 is an extracellular pH- and Na+-sensitive creatine transporter with two kinetic components. Recently, some variants of hMCT9 have been found to be associated with serum uric acid levels, hyperuricemia, and gout. Among these, two single-nucleotide polymorphisms (SNPs) have also been reported: rs550527563 (L93M) and rs2242206 (T258K). However, the effect of these SNPs on hMCT9 transport activity remains unclear. This study aimed to determine the influence of hMCT9 L93M and T258K on transport characteristics. Main methodshMCT9 L93M and T258K were constructed by site-directed mutagenesis and expressed in Xenopus laevis oocyte. Transport activity of uric acid and creatine via hMCT9 were measured by using a Xenopus laevis oocyte heterologous expression system. Key findingsWe assessed the transport activity of uric acid and creatine, and observed that hMCT9-expressing oocytes transported uric acid approximately 3- to 4-fold more than water-injected oocytes. hMCT9 L93M slightly reduced the transport activity of creatine, whereas hMCT9 T258K did not affect the transport activity. Interestingly, hMCT9 T258K abolished Na+ sensitivity and altered the substrate affinity from two components to one. SignificanceIn conclusion, hMCT9 SNPs affect transport activity and characteristics. hMCT9 L93M and T258K may induce dysfunction and contribute to pathologies such as hyperuricemia and gout. This is a first study to evaluate molecular characteristics of hMCT9 SNPs.

Causal relationships between blood lipids and major psychiatric disorders: Univariable and multivariable mendelian randomization analysis

BackgroundWhether the positive associations of blood lipids with psychiatric disorders are causal is uncertain. We conducted this two-sample Mendelian randomization (MR) analysis to comprehensively investigate associations of blood lipids with psychiatric disorders.MethodsUnivariable and multivariable models were established for MR analyses. Inverse variance-weighted (IVW) MR was employed as the main approach; weighted median and MR-Egger were used as sensitivity analysis methods. The possibility of violating MR assumptions was evaluated utilizing several sensitivity analyses, including heterogeneity statistics, horizontal pleiotropy statistics, single SNP analysis, leave-one-out analysis and MR-PRESSO analysis. As instrumental variables, we screened 362 independent single-nucleotide polymorphisms (SNP) related to blood lipids from a recent genome-wide association study involving 76,627 individuals of European ancestry, with a genome-wide significance level of p < 5 × 10− 8. Summary-level information for the six psychiatric disorders was extracted from Psychiatric Genomics Consortium and Alzheimer Disease Genetics Consortium.ResultsWe observed eight significant associations in univariable MR analysis, four of which were corroborated by multivariable MR (MVMR) analysis modified for the other three lipid traits: high-density lipoprotein cholesterol (HDL-C) level with the risk of PTSD (OR = 0.91, 95% CI = 0.85–0.97, p = 0.002) and AD (OR = 0.79, 95% CI = 0.71–0.88, p < 0.001) and triglycerides (TG) level with the risk of MDD (OR = 1.02, 95% CI = 1.003–1.03, p = 0.01) and panic disorder (OR = 0.83, 95% CI = 0.74–0.92, p < 0.001). In addition, four associations were not significant in MVMR analysis after adjustment for three lipid traits: total cholesterol (TC) level with the risk of PTSD, low-density lipoprotein cholesterol (LDL-C) level with the risk of MDD and AD and TG level with the risk of AD.ConclusionsOur results show that blood lipids and psychiatric disorders may be related in a causal manner. This shows that abnormal blood lipid levels may act as reliable biomarker of psychiatric disorders and as suitable targets for their prevention and treatment.

DNA mixture analyses of autosomal single nucleotide polymorphisms for individual identification using droplet digital polymerase-chain reaction and massively parallel sequencing in combination with EuroFormix

ABSTRACT Droplet digital polymerase-chain reaction (ddPCR) and massively parallel sequencing (MPS) can be used to detect an extremely low proportion of non-overlapping single nucleotide polymorphisms (SNPs) of a minor contributor in a DNA mixture. Both semi-continuous and continuous probabilistic likelihood ratio (LR) models can be used to interpret the DNA profiles of mixtures. We analysed 28 identity-informative SNPs in 20 DNA mixtures with various minor to major ratios (1:29, to 1:99) by using our customized ddPCR panel and two MPS panels. The minor contributors of the DNA mixtures were correctly inferred using both semi-continuous and continuous LR models of EuroForMix based on the data of our customized ddPCR and one MPS panel. The accuracy rate of minor contributor assignment was 95% and 90% using semi-continuous and continuous LR models, respectively, using the data of the other MPS panel with lower coverage reads. In conclusion, the quantitative genotype data of SNPs generated using ddPCR can be used for the minor contributor inference for DNA mixtures. The performance of a continuous LR model in minor contributor identification in DNA mixtures may not be superior to that of a semi-continuous LR model when the coverage reads of the MPS panel is insufficient.

Comprehensive application of multiple molecular diagnostic techniques in pre-implantation genetic testing for monogenic.

Pre-implantation genetic testing for monogenic disorders (PGT-M) is an effective approach to reducing the incidence of birth defects by preventing the transmission of inherited diseases to offspring. However, there are still controversies regarding the detection methods and transplantation of embryos. This paper aims to evaluate the effectiveness of different detection technologies applied to PGT-M through a retrospective analysis of clinical detection data. The carrier status of pathogenic mutations and chromosomal copy number variants (CNVs) in 892 embryos was characterized using next-generation sequencing (NGS), single-nucleotide polymorphism (SNP) array, and PCR-based detection technologies. Clinical data from PGT-M cases were retrospectively analyzed to assess the effectiveness of these detection methods in identifying genetic abnormalities in embryos. A total of 829 embryos were analyzed, with 63 being unsuccessful. Our study revealed that the success rate of detecting deletional mutations using Gap-PCR 84.9%, which is lower than that of SNP array (98.7%) and NGS (92.5%). However, no significant difference was observed when detecting point mutations using any of the methods. These findings suggest that, when detecting deletional mutations, SNP array and NGS are more suitable choices compared to Gap-PCR. While SNP array may have a lower resolution and success rate (80.5%) in analyzing CNVs compared to NGS (95.5%), it may still be useful for revealing certain abnormal types. In conclusion, this study found that SNP analysis is advantageous for identifying polygenic and deletional mutations, whereas NGS is more cost-efficient for detecting common monogenic diseases. Additionally, SNP-based haplotyping and PCR-based direct detection of mutations can be used together to enhance the accuracy and success rates of PGT-M. Our findings offer valuable insights for PGT technicians in choosing suitable detection methods for patients.

Possible effect of mutations on serological detection of Borrelia burgdorferi sensu stricto ospC major groups: An in-silico study.

The outer surface protein C (OspC) of the agent of Lyme disease, Borrelia burgdorferi sensu stricto, is a major lipoprotein surface-expressed during early-phase human infections. Antibodies to OspC are used in serological diagnoses. This study explored the hypothesis that serological test sensitivity decreases as genetic similarity of ospC major groups (MGs) of infecting strains, and ospC A (the MG in the strain B31 used to prepare antigen for serodiagnosis assays) decreases. We used a previously published microarray dataset to compare serological reactivity to ospC A (measured as pixel intensity) versus reactivity to 22 other ospC MGs, within a population of 55 patients diagnosed by two-tier serological testing using B. burgdorferi s.s. strain B31 as antigen, in which the ospC MG is OspC A. The difference in reactivity of sera to ospC A and reactivity to each of the other 22 ospC MGs (termed 'reactivity difference') was the outcome variable in regression analysis in which genetic distance of the ospC MGs from ospC A was the explanatory variable. Genetic distance was computed for the whole ospC sequence, and 9 subsections, from Neighbour Joining phylogenetic trees of the 23 ospC MGs. Regression analysis was conducted using genetic distance for the full ospC sequence, and the subsections individually. There was a significant association between the reactivity difference and genetic distance of ospC MGs from ospC A: increased genetic distance reduced reactivity to OspC A. No single ospC subsection sequence fully explained the relationship between genetic distance and reactivity difference. An analysis of single nucleotide polymorphisms supported a biological explanation via specific amino acid modifications likely to change protein binding affinity. This adds support to the hypothesis that genetic diversity of B. burgdorferi s.s. (here specifically OspC) may impact serological diagnostic test performance. Further prospective studies are necessary to explore the clinical implications of these findings.

A genomic survey of Clostridioides difficile isolates from hospitalized patients in Melbourne, Australia.

There has been a decrease in healthcare-associated Clostridioides difficile infection (CDI) in Australia, coupled with an increase in the genetic diversity of strains isolated in these settings, and an increase in community-associated cases. To explore this changing epidemiology, we studied the genetic relatedness of C. difficile isolated from patients at a major hospital in Melbourne, Australia. Whole-genome sequencing of C. difficile isolates from symptomatic (n = 61) and asymptomatic (n = 10) hospital patients was performed. Genomic comparisons were made using single-nucleotide polymorphism (SNP) analysis, ribotyping, and toxin, resistome, and mobilome profiling. C. difficle clade 1 strains were found to be predominant (64/71), with most strains (63/71) encoding both toxins A and B (A+B+). Despite these similarities, only two isolates were genetically related (≤2 SNPs) and a diverse range of ribotypes was detected, with those predominating including ribotypes commonly found in community-associated cases. Five non-toxigenic (A-B-CDT-) clade 1 strains were identified, all in asymptomatic patients. Three clade 4 (A-B+CDT-) and four clade 5 (A+B+CDT+) strains were detected also, with these strains more likely to carry antimicrobial resistance determinants, many of which were associated with mobile genetic elements. Overall, within a single hospital, C. difficile-associated disease was caused by a diverse range of strains, including many strain types associated with community and environmental sources. While strains carried asymptomatically were more likely to be non-toxigenic, toxigenic strains were isolated also from asymptomatic patients, which together suggest the presence of diverse sources of transmission, potentially including asymptomatic patients. IMPORTANCE There has been a decrease in healthcare-associated Clostridioides difficile infection in Australia, but an increase in the genetic diversity of infecting strains, and an increase in community-associated cases. Here, we studied the genetic relatedness of C. difficile isolated from patients at a major hospital in Melbourne, Australia. Diverse ribotypes were detected, including those associated with community and environmental sources. Some types of isolates were more likely to carry antimicrobial resistance determinants, and many of these were associated with mobile genetic elements. These results correlate with those of other recent investigations, supporting the observed increase in genetic diversity and prevalence of community-associated C. difficile, and consequently the importance of sources of transmission other than symptomatic patients. Thus, they reinforce the importance of surveillance for in both hospital and community settings, including asymptomatic carriage, food, animals, and other environmental sources to identify and circumvent important sources of C. difficile transmission.

Application of quasimetagenomics methods to define microbial diversity and subtype Listeria monocytogenes in dairy and seafood production facilities.

Microorganisms frequently colonize surfaces within food production facilities. Detection of Listeria monocytogenes in this setting relies on culture-dependent methods, but the complex dynamics of bacterial interactions within these environments and their impact on pathogen detection remain largely unexplored. To address this challenge, we applied both 16S rRNA and shotgun quasimetagenomic (enriched microbiome) sequencing of swab culture enrichments from five seafood and seven dairy production environments. Utilizing 16S rRNA amplicon sequencing, we observed variability between 355 samples taken from these 12 production facilities and a distinctive microbiome for each environment. With shotgun quasimetagenomic sequencing, we were able to assemble L. monocytogenes metagenome-assembled genomes (MAGs) from 28 of the 32 culture-positive samples. We compared these MAGs to their corresponding whole-genome sequencing assemblies, which resulted in two polyphyletic clades consisting of L. monocytogenes lineages I and II with 13,195 and 25,556 single-nucleotide polymorphism sites, respectively. The remaining four MAGs did not produce sufficient genome coverage. To understand and establish limits for pathogen detection and subtyping using shotgun quasimetagenomics, these same data sets were downsampled in slilico to produce a titration series of abundances of L. monocytogenes and analyzed. Pathogen detection was achieved for all downsampled data sets, even those with only 3× genome coverage. This study contributes to the understanding of microbial diversity within food production environments and presents insights into the level of genome coverage needed in a metagenome sequencing data set to detect, subtype, and source track a foodborne pathogenIMPORTANCEIn developed countries, the human diet is predominated by food commodities, which have been manufactured, processed, and stored in a food production facility. Little is known about the application of metagenomic sequencing approaches for detecting foodborne pathogens, such as L. monocytogenes, and characterizing microbial diversity in food production ecosystems. In this work, we investigated the utility of 16S rRNA amplicon and quasimetagenomic sequencing for the taxonomic and phylogenetic classification of Listeria culture enrichments of environmental swabs collected from dairy and seafood production facilities. We demonstrated that single-nucleotide polymorphism (SNP) analyses of L. monocytogenes metagenome-assembled genomes (MAGs) from quasimetagenomic data sets can achieve similar resolution as culture isolate whole-genome sequencing. To further understand the impact of genome coverage on MAG SNP cluster resolution, an in silico downsampling approach was employed to reduce the percentage of target pathogen sequence reads, providing an initial estimate of required MAG coverage for subtyping resolution of L. monocytogenes.

Introduction pathways and evolutionary mechanisms of alien species of Lolium spreading across sandy coasts in Japan

Abstract Estimating the role of specific processes in the spread of alien species necessitates the determination of introduction pathways and source populations of invaded areas. Alien grasses in the genus Lolium that have extensively invaded Japan provide a unique opportunity to estimate the expansion process through direct comparison between source and naturalised populations because the introduction pathways, contaminants in grain commodities and commercial cultivars for fodder crops or revegetation materials are well‐known. Therefore, by directly comparing source and naturalised populations, we estimated the introduction pathways and whether adaptative evolution occurred in Lolium species on sandy coasts in Japan. Lolium individuals sampled from naturalised populations in croplands, seaports, and sandy coasts were compared with those from two introduction sources for morphological and genetic variations based on a genome‐wide single nucleotide polymorphism analysis and a common garden experiment. Furthermore, we conducted a reciprocal transplant experiment between cropland and sandy coast. Populations naturalised in croplands were closely related to the cultivars, whereas those naturalised in seaports and sandy coasts were associated with contaminants. These results indicate that the cropland and sandy coast populations are derived from cultivars and contaminants, respectively. In addition, asymmetric gene flow from cropland populations to sandy coast populations was observed. The reciprocal transplant experiment clearly demonstrated the home site advantage; populations derived from croplands yielded higher floret numbers than those derived from other habitats at the cropland site; sandy coast populations had higher survival rates than those from croplands at the coastal site. Port populations exhibited a similar tendency as sandy coast populations, indicating that contaminants may be originally adapted to salty and dry environments, such as that in sandy coasts. The flowering phenology in the sandy coast populations evolved in the late flowering; therefore, late flowering alleles may have been transferred from cropland populations to sandy coast populations. Synthesis. We demonstrated that two congeneric species with different ecological characteristics were introduced through multiple introduction pathways and spread across different habitats. A direct comparison between source and naturalised populations can considerably elucidate the patterns and processes of biological invasions.

Prevalence of multidrug-resistant hypervirulent Klebsiella pneumoniae without defined hypervirulent biomarkers in Anhui, China: a new dimension of hypervirulence.

Klebsiella pneumoniae is an opportunistic pathogen that mainly causes nosocomial infections and hospital-associated pneumonia in elderly and immunocompromised people. However, multidrug-resistant hypervirulent K. pneumoniae (MDR-hvKp) has emerged recently as a serious threat to global health that can infect both immunocompromised and healthy individuals. It is scientifically established that plasmid-mediated regulator of mucoid phenotype genes (rmpA and rmpA2) and other virulence factors (aerobactin and salmochelin) are mainly responsible for this phenotype. In this study, we collected 23 MDR-hvKp isolates and performed molecular typing, whole genome sequencing, comparative genomic analysis, and phenotypic experiments, including the Galleria mellonella infection model, to reveal its genetic and phenotypic features. Meanwhile, we discovered two MDR-hvKp isolates (22122315 and 22091569) that showed a wide range of hypervirulence and hypermucoviscosity without rmpA and rmpA2 and any virulence factors. In phenotypic experiments, isolate 22122315 showed the highest hypervirulence (infection model) with significant mucoviscosity, and conversely, isolate 22091569 exhibited the highest mucoviscosity (string test) with higher virulence compared to control. These two isolates carried carbapenemase (blaKPC - 2), β-lactamase (blaOXA - 1, blaTEM - 1B), extended-spectrum β-lactamase (ESBL) genes (blaCTX - M - 15, blaSHV - 106), outer membrane protein-coding genes (ompA), fimbriae encoding genes (ecpABCDER), and enterobactin coding genes (entAB, fepC). In addition, single nucleotide polymorphism analysis indicated that both isolates, 22122315 and 22091569, were found to have novel mutations in loci FEBNDAKP_03184 (c. 2084A > C, p. Asn695Thr), and EOFMAFIB_02276 (c. 1930C > A, p. Pro644Thr), respectively. Finally, NCBI blast analysis suggested these mutations are located in the wzc of the capsule polysaccharide (cps) region and are responsible for putative tyrosine kinase. This study would be a strong reference for enhancing the current understanding of identifying the MDR-hvKp isolates that lacked both mucoid regulators and virulence factors.

Genetic predictors of neurocognitive outcomes in survivors of pediatric brain tumors

BackgroundNeurocognitive deficits are common in pediatric brain tumor survivors. The use of single nucleotide polymorphism (SNP) analysis in DNA repair genes may identify children treated with radiation therapy for brain tumors at increased risk for treatment toxicity and adverse neurocognitive outcomes.MaterialsThe Human 660W-Quad v1.0 DNA BeadChip analysis (Illumina) was used to evaluate 1048 SNPs from 59 DNA repair genes in 46 subjects. IQ testing was measured by the Wechsler Intelligence Scale for Children. Linear regression was used to identify the 10 SNPs with the strongest association with IQ scores while adjusting for radiation type.ResultsThe low vs high IQ patient cohorts were well matched for time from first treatment to most recent IQ, first treatment age, sex, and treatments received. 5 SNPs on 3 different genes (CYP29, XRCC1, and BRCA1) and on 3 different chromosomes (10, 19, and 17) had the strongest association with most recent IQ score that was not modified by radiation type. Furthermore, 5 SNPs on 4 different genes (WRN, NR3C1, ERCC4, RAD51L1) on 4 different chromosomes (8, 5, 16, 14) had the strongest association with change in IQ independent of radiation type, first IQ, and years between IQ measures.ConclusionsSNPs offer the potential to predict adverse neurocognitive outcomes in pediatric brain tumor survivors. Our results require validation in a larger patient cohort. Improving the ability to identify children at risk of treatment related neurocognitive deficits could allow for better treatment stratification and early cognitive interventions.

Spatial and phylogenetic patterns reveal hidden infection sources of Bacillus anthracis in an anthrax outbreak in Son La province, Vietnam.

Bacillus anthracis, the bacterial cause of anthrax, is a zoonosis affecting livestock and wildlife often spilling over into humans. In Vietnam, anthrax has been nationally reportable since 2015 with cases occurring annually, mostly in the northern provinces. In April 2022, an outbreak was reported in Son La province following the butchering of a water buffalo, Bubalus bubalis. A total of 137 humans from three villages were likely exposed to contaminated meat from the animal. Early epidemiological investigations suggested a single animal was involved in all exposures. Five B. anthracis isolates were recovered from human clinical cases along with one from the buffalo hide, another from associated maggots, and one from soil at the carcass site. The isolates were whole genome sequenced, allowing global, regional, and local molecular epidemiological analyses of the outbreak strains. All recovered B. anthracis belong to the A.Br.001/002 lineage based on canonical single nucleotide polymorphism analysis (canSNP). Although not previously identified in Vietnam, this lineage has been identified in the nearby countries of China, India, Indonesia, Thailand, as well as Australia. A twenty-five marker multi-locus variable number tandem repeat analysis (MLVA-25) was used to investigate the relationship between human, soil, and buffalo strains. Locally, four MLVA-25 genotypes were identified from the eight isolates. This level of genetic diversity is unusual for the limited geography and timing of cases and differs from past literature using MLVA-25. The coupled spatial and phylogenetic data suggest this outbreak originated from multiple, likely undetected, animal sources. These findings were further supported by local news reports that identified at least two additional buffalo deaths beyond the initial animal sampled in response to the human cases. Future outbreak response should include intensive surveillance for additional animal cases and additional molecular epidemiological traceback to identify pathogen sources.

Unraveling the genomic background of resilience based on variability in milk yield and milk production levels in North American Holstein cattle through genome-wide association study and Mendelian randomization analyses

Breeding more resilient animals will benefit the dairy cattle industry in the long term, especially as global climate changes become more severe. Previous studies have reported genetic parameters for various milk yield-based resilience indicators, but the underlying genomic background of these traits remain unknown. In this study, we conducted genome-wide association studies of 62,029 single nucleotide polymorphisms (SNPs) with 4 milk yield-based resilience indicators, including the weighted occurrence frequency (wfPert) and accumulated milk losses (dPert) of milk yield perturbations, and log-transformed variance (LnVar) and lag-1 autocorrelation (rauto) of daily yield residuals. These variables were previously derived from 5.6 million daily milk yield records from 21,350 lactations (parities 1 to 3) of 11,787 North American Holstein cows. The average daily milk yield (ADMY) throughout lactation was also included to compare the shared genetic background of resilience indicators with milk yield. The differential genetic background of these indicators was first revealed by the significant genomic regions identified and significantly-enriched biological pathways of positional candidate genes, which confirmed the genetic difference among resilience indicators. Interestingly, the functional analyses of candidate genes suggested that the regulation of intestinal homeostasis is most likely affecting resilience derived based on variability in milk yield. Based on Mendelian randomization analyses of multiple instrumental SNPs, we further found an unfavorable causal association of ADMY with LnVar. In conclusion, the resilience indicators evaluated are genetically different traits and there are causal associations of milk yield with some of the resilience indicators evaluated. In addition to providing biological insights into the molecular regulation mechanisms of resilience derived based on variability in milk yield, this study also indicates the need for developing selection indexes combining multiple indicator traits and taking into account their genetic relationship for breeding more resilient dairy cattle.