5NosoAE: a web server for nosocomial bacterial antibiogram investigation and epidemiology survey.

5NosoAE is a webserver that can be used for nosocomial bacterial analysis including the identification of similar strains based on antimicrobial resistance profiles (antibiogram) and the spatiotemporal distribution visualization and phylogenetic analysis of identified strains with similar antibiograms. The extensive use of antibiotics has caused many pathogenic bacteria to develop multiple drug resistance, resulting in clinical infection treatment challenges and posing a major threat to global public health. Relevant studies have investigated the key determinants of antimicrobial resistance in the whole-genome sequence of bacteria. However, a web server is currently not available for performing large-scale strain searches according to antimicrobial resistance profiles and visualizing epidemiological information including the spatiotemporal distribution, antibiogram heatmap, and phylogeny of identified strains. Here, we implemented these functions in the new server, referred to as 5NosoAE. This server accepts the genome sequence file in the FASTA format of five nosocomial bacteria, namely Acinetobacter baumannii, Pseudomonas aeruginosa, Klebsiella pneumoniae, Enterococcus faecium and Staphylococcus aureus for query. All visualizations are implemented in JavaScript and PHP. This server will be useful for physicians and epidemiologists involved in research on infectious disease. The 5NosoAE platform is available at https://nosoae.imst.nsysu.edu.tw.

MIB2: metal ion-binding site prediction and modeling server.

MOTIVATION: MIB2 (metal ion-binding) attempts to overcome the limitation of structure-based prediction approaches, with many proteins lacking a solved structure. MIB2 also offers more accurate prediction performance and more metal ion types. RESULTS: MIB2 utilizes both the (PS)2 method and the AlphaFold Protein Structure Database to acquire predicted structures to perform metal ion docking and predict binding residues. MIB2 offers marked improvements over MIB by collecting more MIB residue templates and using the metal ion type-specific scoring function. It offers a total of 18 types of metal ions for binding site predictions. AVAILABILITY AND IMPLEMENTATION: Freely available on the web at http://bioinfo.cmu.edu.tw/MIB2/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Emergence of Vibrio cholerae O1 Sequence Type 75 in Taiwan.

We investigated the epidemiology of cholera in Taiwan during 2002-2018. Vibrio cholerae sequence type (ST) 75 clone emerged in 2009 and has since become more prevalent than the ST69 clone from a previous pandemic. Closely related ST75 strains have emerged in 4 countries and may now be widespread in Asia.

Transmission and evolution of OXA-48-producing Klebsiella pneumoniae ST11 in a single hospital in Taiwan.

OBJECTIVES: Epidemic spread of OXA-48-producing Klebsiella pneumoniae, mainly mediated by the transmission of a blaOXA-48-carrying plasmid, has threatened global health during the last decade. Since its introduction to Taiwan in 2013, OXA-48 has become the second most common carbapenemase. We described the transmission and evolution of an OXA-producing K. pneumoniae clone in a single hospital. METHODS: Twenty-two OXA-48 K. pneumoniae were isolated between October 2013 and December 2015. Comparative genomic analysis was performed based on the WGS data generated with Illumina and MinION techniques. RESULTS: Seventeen of the 22 OXA-48 K. pneumoniae that belonged to ST11, with the same capsular genotype, KL64, and differed from each other by seven or fewer SNPs, were considered outbreak strains. Eight of the 17 outbreak strains harboured a 65499 bp blaOXA-48-carrying IncL plasmid (called pOXA48). pOXA48 was absent from the remaining nine strains. Instead, a 24.9 kb blaOXA-48-carrying plasmid fragment was integrated into a prophage region of their chromosomes. Transmission routes of the ST11_KL64 K. pneumoniae sublineages, which carried either pOXA48 or chromosomally integrated blaOXA-48, were reconstructed. CONCLUSIONS: Clonal expansion of ST11_KL64 sublineages contributed to the nosocomial outbreak of OXA-48 K. pneumoniae. The chromosome-borne blaOXA-48 lineage emerged during a 2 year period in a single hospital. Dissemination of OXA-48, which is vertically transmitted in K. pneumoniae even in the absence of selective pressure from antimicrobials, deserves public health attention.

Emergence of Multidrug-Resistant Salmonella enterica Serovar Goldcoast Strains in Taiwan and International Spread of the ST358 Clone.

Salmonella enterica serovar Goldcoast infection was rare in Taiwan; it was not detected in routine surveillance from 2004 to 2013. This serovar was first identified in 2014, but the frequency of infection remained low until 2017. From 2014 to 2016, all but one isolate was pan-susceptible. S Goldcoast infections abruptly increased in 2018, and all isolates were multidrug-resistant (MDR). All MDR isolates harbored an IncHI2 plasmid, and the majority carried 14 antimicrobial resistance genes, aac(3)-IId, aadA22, aph(3')-Ia, aph(6)-Id, blaTEM-1B, blaCTX-M-55, lnu(F), floR, qnrS13, arr-2, sul2, sul3, tet(A), and dfrA14. S Goldcoast strains recovered in Taiwan and 96 of 99 strains from Germany, the Netherlands, the United Kingdom, and the United States belonged to sequence type 358 (ST358). Whole-genome single-nucleotide polymorphism and core genome multilocus sequence type analyses revealed that all strains of the ST358 clone shared a high degree of genetic relatedness. The present study highlighted that a dramatic increase in S Goldcoast infections followed the emergence of MDR strains and indicated that a genetically closely related S Goldcoast ST358 clone may have widespread significance internationally.

Genetic Relationships among Multidrug-Resistant Salmonella enterica Serovar Typhimurium Strains from Humans and Animals.

We identified 20 to 22 resistance genes, carried in four incompatibility groups of plasmids, in each of five genetically closely related Salmonella enterica serovar Typhimurium strains recovered from humans, pigs, and chickens. The genes conferred resistance to aminoglycosides, chloramphenicol, sulfonamides, trimethoprim, tetracycline, fluoroquinolones, extended-spectrum cephalosporins and cefoxitin, and azithromycin. This study demonstrates the transmission of multidrug-resistant Salmonella strains among humans and food animals and may be the first identification of mphA in azithromycin-resistant Salmonella strains in Taiwan.

Dissemination of mcr-1-Carrying Plasmids among Colistin-Resistant Salmonella Strains from Humans and Food-Producing Animals in Taiwan.

We detected the colistin resistance gene mcr-1 in four Salmonella serovars isolated from humans and animals with diarrhea. The resistance gene was carried on different plasmids. One mcr-1-carrying conjugative plasmid, a variant of pHNSHP45, was disseminated among Salmonella isolates recovered from humans, pigs, and chickens.

Azithromycin-Nonsusceptible Shigella flexneri 3a in Men Who Have Sex with Men, Taiwan, 2015-2016.

We report an outbreak of azithromycin-nonsusceptible Shigella flexneri 3a infection in Taiwan associated with men who have sex with men. The bacterial strains belonged to the sublineage A of a recently reported outbreak lineage associated with men who have sex with men, characterized by reduced azithromycin susceptibility and circulation in shigellosis low-risk regions.

Chromosome-mediated multidrug resistance in Salmonella enterica serovar Typhi.

A salmonella genomic island, designated SGI11, was found in 18 of 26 multidrug-resistant Salmonella enterica serovar Typhi isolates from Bangladesh. SGI11 was an IS1 composite transposon and carried 7 resistance genes that conferred resistance to 5 first-line antimicrobials. Eleven of the 18 SGI11-carrying S. Typhi isolates had developed resistance to high levels of ciprofloxacin.

Antimicrobial resistance in Salmonella enterica Serovar Typhi isolates from Bangladesh, Indonesia, Taiwan, and Vietnam.

We characterized Salmonella enterica serovar Typhi isolates from Bangladesh, Indonesia, Taiwan, and Vietnam to investigate their genetic relatedness and antimicrobial resistance. The isolates from Bangladesh and Vietnam were genetically closely related but were distant from those from Indonesia and Taiwan. All but a few isolates from Indonesia and Taiwan were susceptible to all antimicrobials tested. The majority of isolates from Bangladesh and Vietnam were multidrug resistant (MDR) and belonged to the widespread haplotype H58 clone. IncHI1 plasmids were detected in all MDR S. Typhi isolates from Vietnam but in only 15% of MDR isolates from Bangladesh. Resistance genes in the majority of MDR S. Typhi isolates from Bangladesh should reside in the chromosome. Among the isolates from Bangladesh, 82% and 40% were resistant to various concentrations of nalidixic acid and ciprofloxacin, respectively. Several resistance mechanisms, including alterations in gyrase A, the presence of QnrS, and enhanced efflux pumps, were involved in the reduced susceptibility and resistance to fluoroquinolones. Intensive surveillance is necessary to monitor the spread of chromosome-mediated MDR and fluoroquinolone-resistant S. Typhi emerging in Bangladesh.

CPred: a web server for predicting viable circular permutations in proteins.

Circular permutation (CP) is a protein structural rearrangement phenomenon, through which nature allows structural homologs to have different locations of termini and thus varied activities, stabilities and functional properties. It can be applied in many fields of protein research and bioengineering. The limitation of applying CP lies in its technical complexity, high cost and uncertainty of the viability of the resulting protein variants. Not every position in a protein can be used to create a viable circular permutant, but there is still a lack of practical computational tools for evaluating the positional feasibility of CP before costly experiments are carried out. We have previously designed a comprehensive method for predicting viable CP cleavage sites in proteins. In this work, we implement that method into an efficient and user-friendly web server named CPred (CP site predictor), which is supposed to be helpful to promote fundamental researches and biotechnological applications of CP. The CPred is accessible at http://sarst.life.nthu.edu.tw/CPred.

ProbioMinServer: an integrated platform for assessing the safety and functional properties of potential probiotic strains.

Motivation: ProbioMinServer is a platform designed to help researchers access information on probiotics regarding a wide variety of characteristics, such as safety (e.g. antimicrobial resistance, virulence, pathogenic, plasmid, and prophage genes) and functionality (e.g. functional classes, carbohydrate-active enzyme, and metabolite gene cluster profile). Because probiotics are functional foods, their safety and functionality are a crucial part of health care. Genomics has become a crucial methodology for investigating the safety and functionality of probiotics in food and feed. This shift is primarily attributed to the growing affordability of next-generation sequencing technologies. However, no integrated platform is available for simultaneously evaluating probiotic strain safety, investigating probiotic functionality, and identifying known phylogenetically related strains. Results: Thus, we constructed a new platform, ProbioMinServer, which incorporates these functions. ProbioMinServer accepts whole-genome sequence files in the FASTA format. If the query genome belongs to the 25 common probiotic species collected in our database, the server performs a database search and analyzes the core-genome multilocus sequence typing. Front-end applications were implemented in JavaScript with a bootstrap framework, and back-end programs were implemented using PHP, Perl, and Python. ProbioMinServer can help researchers quickly and easily retrieve information on the safety and functionality of various probiotics. Availability and implementation: The platform is available at https://probiomindb.imst.nsysu.edu.tw.

Coffee as a dietary strategy to prevent SARS-CoV-2 infection.

BACKGROUND: To date, most countries lifted the restriction requirement and coexisted with SARS-CoV-2. Thus, dietary behavior for preventing SARS-CoV-2 infection becomes an interesting issue on a daily basis. Coffee consumption is connected with reduced COVID-19 risk and correlated to COVID-19 severity. However, the mechanisms of coffee for the reduction of COVID-19 risk are still unclear. RESULTS: Here, we identified that coffee can inhibit multiple variants of the SARS-CoV-2 infection by restraining the binding of the SARS-CoV-2 spike protein to human angiotensin-converting enzyme 2 (ACE2), and reducing transmembrane serine protease 2 (TMPRSS2) and cathepsin L (CTSL) activity. Then, we used the method of "Here" (HRMS-exploring-recombination-examining) and found that isochlorogenic acid A, B, and C of coffee ingredients showed their potential to inhibit SARS-CoV-2 infection (inhibitory efficiency 43-54%). In addition, decaffeinated coffee still preserves inhibitory activity against SARS-CoV-2. Finally, in a human trial of 64 subjects, we identified that coffee consumption (approximately 1-2 cups/day) is sufficient to inhibit infection of multiple variants of SARS-CoV-2 entry, suggesting coffee could be a dietary strategy to prevent SARS-CoV2 infection. CONCLUSIONS: This study verified moderate coffee consumption, including decaffeination, can provide a new guideline for the prevention of SARS-CoV-2. Based on the results, we also suggest a coffee-drinking plan for people to prevent infection in the post-COVID-19 era.

LmTraceMap: A Listeria monocytogenes fast-tracing platform for global surveillance.

Listeria monocytogenes can cause listeriosis, and people with hypoimmunity such as pregnant women, infants and fetuses are at high risk of invasive infection. Although the incidence of listeriosis is low, the fatality rate is high. Therefore, continual surveillance and rapid epidemiological investigation are crucial for addressing L. monocytogenes. Because of the popularity of next-generation sequencing, obtaining the whole-genome sequence of a bacterium is easy. Several genome-based typing methods are available, and core-genome multilocus sequence typing (cgMLST) is the most recognized methods. Using cgMLST typing to compare L. monocytogenes whole-genome sequences (WGS) with those obtained across distinct regions is beneficial. However, the concern is how to incorporate the powerful cgMLST method into investigations, such as by using source tracing. Herein, we present an easy-to-use web service called-LmTraceMap (http://lmtracemap.cgu.edu.tw/hua_map/test/upload.php; http://120.126.17.192/hua_map/test/upload.php) that can help public-health professionals rapidly trace closely related isolates worldwide and visually inspect them in search results on a world map with labeled epidemiological data. We expect the proposed service to improve the convenience of public health investigations.

Structural insights into EphA4 unconventional activation from prediction of the EphA4 and its complex with ribonuclease 1.

It has been shown that several ribonuclease (RNase) A superfamily proteins serve as ligands of receptor tyrosine kinases (RTKs), representing a new concept for ligand/receptor interaction. Moreover, recent studies indicate high clinical values for this type of ligand/RTK interactions. However, there is no structural report for this new family of ligand/receptor. In an attempt to understand how RNase and RTK may interact, we focused on the RNase1/ephrin type-A receptor 4 (EphA4) complex and predicted their structure by using the state-of-the-art machine learning method, AlphaFold and its derivative method, AF2Complex. In this model, electrostatic force plays an essential role for the specific ligand/receptor interaction. We found the R39 of RNase1 is the key residue for EphA4-binding and activation. Mutation on this residue causes disruption of an essential basic patch, resulting in weaker ligand-receptor association and leading to the loss of activation. By comparing the surface charge distribution of the RNase A superfamily, we found the positively charged residues on the RNase1 surface is more accessible for EphA4 forming salt bridges than other RNases. Furthermore, RNase1 binds to the ligand-binding domain (LBD) of EphA4, which is responsible for the traditional ligand ephrin-binding. Our model reveals the location of RNase1 on EphA4 partially overlaps with that of ephrin-A5, a traditional ligand of EphA4, suggesting steric hindrance as the basis by which the ephrin-A5 precludes interactions of RNase1 with EphA4. Together, our discovery of RNase1/EphA4 interface provides a potential treatment strategy by blocking the RNase1-EphA4 axis.

cgMLST@Taiwan: A web service platform for Vibrio cholerae cgMLST profiling and global strain tracking.

BACKGROUND: Cholera, a rapidly dehydrating diarrheal disease caused by toxigenic Vibrio cholerae, is a leading cause of morbidity and mortality in some regions of the world. Core genome multilocus sequence typing (cgMLST) is a promising approach in generating genetic fingerprints from whole-genome sequencing (WGS) data for strain comparison among laboratories. METHODS: We constructed a V. cholerae core gene allele database using an in-house developed computational pipeline, a database with cgMLST profiles converted from genomic sequences from the National Center for Biotechnology Information, and built a REST-based web accessible via the Internet. RESULTS: We built a web service platform-cgMLST@Taiwan and installed a V. cholerae allele database, a cgMLST profile database, and computational tools for generating V. cholerae cgMLST profiles (based on 3,017 core genes), performing rapid global strain tracking, and clustering analysis of cgMLST profiles. This web-based platform provides services to researchers, public health microbiologists, and physicians who use WGS data for the investigation of cholera outbreaks and tracking of V. cholerae strain transmission across countries and geographic regions. The cgMLST@Taiwan is accessible at http://rdvd.cdc.gov.tw/cgMLST.

Disulfiram blocked cell entry of SARS-CoV-2 via inhibiting the interaction of spike protein and ACE2.

Disulfiram is an FDA-approved drug that has been used to treat alcoholism and has demonstrated a wide range of anti-cancer, anti-bacterial, and anti-viral effects. In the global COVID-19 pandemic, there is an urgent need for effective therapeutics and vaccine development. According to recent studies, disulfiram can act as a potent SARS-CoV-2 replication inhibitor by targeting multiple SARS-CoV-2 non-structural proteins to inhibit viral polyprotein cleavage and RNA replication. Currently, disulfiram is under evaluation in phase II clinical trials to treat COVID-19. With more and more variants of the SARS-CoV-2 worldwide, it becomes critical to know whether disulfiram can also inhibit viral entry into host cells for various variants and replication inhibition. Here, molecular and cellular biology assays demonstrated that disulfiram could interrupt viral spike protein binding with its receptor ACE2. By using the viral pseudo-particles (Vpps) of SARS-CoV-2, disulfiram also showed the potent activity to block viral entry in a cell-based assay against Vpps of different SARS-CoV-2 variants. We further established a live virus model system to support the anti-viral entry activity of disulfiram with the SARS-CoV-2 virus. Molecular docking revealed how disulfiram hindered the binding between the ACE2 and wild-type or mutated spike proteins. Thus, our results indicate that disulfiram has the capability to block viral entry activity of different SARS-CoV-2 variants. Together with its known anti-replication of SARS-CoV-2, disulfiram may serve as an effective therapy against different SARS-CoV-2 variants.

Peimine inhibits variants of SARS-CoV-2 cell entry via blocking the interaction between viral spike protein and ACE2.

Coronavirus disease 2019 (COVID-19) is caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Several vaccines against SARS-CoV-2 have been approved; however, variants of concern (VOCs) can evade vaccine protection. Therefore, developing small compound drugs that directly block the interaction between the viral spike glycoprotein and ACE2 is urgently needed to provide a complementary or alternative treatment for COVID-19 patients. We developed a viral infection assay to screen a library of approximately 126 small molecules and showed that peimine inhibits VOCs viral infections. In addition, a fluorescence resonance energy transfer (FRET) assay showed that peimine suppresses the interaction of spike and ACE2. Molecular docking analysis revealed that peimine exhibits a higher binding affinity for variant spike proteins and is able to form hydrogen bonds with N501Y in the spike protein. These results suggest that peimine, a compound isolated from Fritillaria, may be a potent inhibitor of SARS-CoV-2 variant infection. PRACTICAL APPLICATIONS: In this study, we identified a naturally derived compound of peimine, a major bioactive alkaloid extracted from Fritillaria, that could inhibit SARS-CoV-2 variants of concern (VOCs) viral infection in 293T/ACE2 and Calu-3 lung cells. In addition, peimine blocks viral entry through interruption of spike and ACE2 interaction. Moreover, molecular docking analysis demonstrates that peimine has a higher binding affinity on N501Y in the spike protein. Furthermore, we found that Fritillaria significantly inhibits SARS-CoV-2 viral infection. These results suggested that peimine and Fritillaria could be a potential functional drug and food for COVID-19 patients.

Computational analysis of a novel mutation in ETFDH gene highlights its long-range effects on the FAD-binding motif.

BACKGROUND: Multiple acyl-coenzyme A dehydrogenase deficiency (MADD) is an autosomal recessive disease caused by the defects in the mitochondrial electron transfer system and the metabolism of fatty acids. Recently, mutations in electron transfer flavoprotein dehydrogenase (ETFDH) gene, encoding electron transfer flavoprotein:ubiquinone oxidoreductase (ETF:QO) have been reported to be the major causes of riboflavin-responsive MADD. To date, no studies have been performed to explore the functional impact of these mutations or their mechanism of disrupting enzyme activity. RESULTS: High resolution melting (HRM) analysis and sequencing of the entire ETFDH gene revealed a novel mutation (p.Phe128Ser) and the hotspot mutation (p.Ala84Thr) from a patient with MADD. According to the predicted 3D structure of ETF:QO, the two mutations are located within the flavin adenine dinucleotide (FAD) binding domain; however, the two residues do not have direct interactions with the FAD ligand. Using molecular dynamics (MD) simulations and normal mode analysis (NMA), we found that the p.Ala84Thr and p.Phe128Ser mutations are most likely to alter the protein structure near the FAD binding site as well as disrupt the stability of the FAD binding required for the activation of ETF:QO. Intriguingly, NMA revealed that several reported disease-causing mutations in the ETF:QO protein show highly correlated motions with the FAD-binding site. CONCLUSIONS: Based on the present findings, we conclude that the changes made to the amino acids in ETF:QO are likely to influence the FAD-binding stability.

A machine learning-based typing scheme refinement for Listeria monocytogenes core genome multilocus sequence typing with high discriminatory power for common source outbreak tracking.

BACKGROUND: As whole-genome sequencing for pathogen genomes becomes increasingly popular, the typing methods of gene-by-gene comparison, such as core genome multilocus sequence typing (cgMLST) and whole-genome multilocus sequence typing (wgMLST), are being routinely implemented in molecular epidemiology. However, some intrinsic problems remain. For example, genomic sequences with varying read depths, read lengths, and assemblers influence the genome assemblies, introducing error or missing alleles into the generated allelic profiles. These errors and missing alleles might create "specious discrepancy" among closely related isolates, thus making accurate epidemiological interpretation challenging. In addition, the rapid growth of the cgMLST allelic profile database can cause problems related to storage and maintenance as well as long query search times. METHODS: We attempted to resolve these issues by decreasing the scheme size to reduce the occurrence of error and missing alleles, alleviate the storage burden, and improve the query search time. The challenge in this approach is maintaining the typing resolution when using fewer loci. We achieved this by using a popular artificial intelligence technique, XGBoost, coupled with Shapley additive explanations for feature selection. Finally, 370 loci from the original 1701 cgMLST loci of Listeria monocytogenes were selected. RESULTS: Although the size of the final scheme (LmScheme_370) was approximately 80% lower than that of the original cgMLST scheme, its discriminatory power, tested for 35 outbreaks, was concordant with that of the original cgMLST scheme. Although we used L. monocytogenes as a demonstration in this study, the approach can be applied to other schemes and pathogens. Our findings might help elucidate gene-by-gene-based epidemiology.