Search

Abstract Background Emergence of antibiotic resistance in bacteria is an important threat to global health. Antibiotic resistance genes (ARGs) are some of the key components to define bacterial resistance and their spread in different environments. Identification of ARGs, particularly from high-throughput sequencing data of the specimens, is the state-of-the-art method for comprehensively monitoring their spread and evolution. Current computational methods to identify ARGs mainly rely on alignment-based sequence similarities with known ARGs. Such approaches are limited by choice of reference databases and may potentially miss novel ARGs. The similarity thresholds are usually simple and could not accommodate variations across different gene families and regions. It is also difficult to scale up when sequence data are increasing. Results In this study, we developed ARGNet, a deep neural network that incorporates an unsupervised learning autoencoder model to identify ARGs and a multiclass classification convolutional neural network to classify ARGs that do not depend on sequence alignment. This approach enables a more efficient discovery of both known and novel ARGs. ARGNet accepts both amino acid and nucleotide sequences of variable lengths, from partial (30–50 aa; 100–150 nt) sequences to full-length protein or genes, allowing its application in both target sequencing and metagenomic sequencing. Our performance evaluation showed that ARGNet outperformed other deep learning models including DeepARG and HMD-ARG in most of the application scenarios especially quasi-negative test and the analysis of prediction consistency with phylogenetic tree. ARGNet has a reduced inference runtime by up to 57% relative to DeepARG. Conclusions ARGNet is flexible, efficient, and accurate at predicting a broad range of ARGs from the sequencing data. ARGNet is freely available at https://github.com/id-bioinfo/ARGNet , with an online service provided at https://ARGNet.hku.hk . Video Abstract

Background: Healthcare-associated COVID-19 infections caused by SARS-CoV-2 have increased morbidity and mortality. Hospitals and skilled nursing facilities (SNFs) have been challenged by infection control and management. Methods: This case study presents an outbreak investigation in a COVID-19-designated hospital and a hospital-based SNF. Real-time polymerase chain reaction (PCR) and other studies were performed on samples obtained from SNF residents, hospital patients, and healthcare workers (HCWs). The results of the laboratory tests and field epidemiological data were analyzed. Genome sequencing and phylogenetic analysis of SARS-CoV-2 were performed to identify the associations between cases. The tracer gas was released and recorded by a thermal imaging camera to investigate the spatial relations within clusters. Results: During the outbreak, 29 COVID-19 infections in 3 clusters were identified through hospital-wide, risk-guided, and symptom-driven PCR tests. This included 12 HCWs, 5 patients, and 12 SNF residents who had been hospitalized for at least 14 days. Serology tests did not identify any cases among the PCR-negative individuals. The phylogenetic analysis revealed that viral strains from the 3 clusters shared a common mutation of G3994T and were phylogenetically related, which suggested that this outbreak had a common source rather than multiple introductions from the community. Linked cases exhibited vertical spatial distribution, and the sulfur hexafluoride release test confirmed a potential airborne transmission. Conclusions: This report addressed the advantage of a multi-disciplinary team in outbreak investigation. Identifying an airborne transmission within an outbreak highlighted the importance of regular maintenance of ventilation systems.

Taiwan’s experience with severe acute respiratory syndrome coronavirus (SARS-CoV) in 2003 guided its development of strategies to defend against SARS-CoV-2 in 2020, which enabled the successful control of Coronavirus disease 2019 (COVID-19) cases from 2020 through March 2021. However, in late-April 2021, the imported Alpha variant began to cause COVID-19 outbreaks at an exceptional rate in Taiwan. In this study, we aimed to determine what epidemiological conditions enabled the SARS-CoV-2 Alpha variant strains to become dominant and decline later during a surge in the outbreak. In conjunction with contact-tracing investigations, we used our bioinformatics software, CoVConvert and IniCoV, to analyze whole-genome sequences of 101 Taiwan Alpha strains. Univariate and multivariable regression analyses revealed the epidemiological factors associated with viral dominance. Univariate analysis showed the dominant Alpha strains were preferentially selected in the surge’s epicenter (p = 0.0024) through intensive human-to-human contact and maintained their dominance for 1.5 months until the Zero-COVID Policy was implemented. Multivariable regression found that the epidemic periods (p = 0.007) and epicenter (p = 0.001) were two significant factors associated with the dominant virus strains spread in the community. These dominant virus strains emerged at the outbreak’s epicenter with frequent human-to-human contact and low vaccination coverage. The Level 3 Restrictions and Zero-COVID policy successfully controlled the outbreak in the community without city lockdowns. Our integrated method can identify the epidemiological conditions for emerging dominant virus with increasing epidemiological potential and support decision makers in rapidly containing outbreaks using public health measures that target fast-spreading virus strains.

Since April 2022, waves of SARS-CoV-2 Omicron variant cases have surfaced in Taiwan and spread throughout the island. Using high-throughput sequencing of the SARS-CoV-2 genome, we analyzed 2,405 PCR-positive swab samples from 2,339 persons and identified the Omicron BA.2.3.7 variant as a major lineage within recent community outbreaks in Taiwan.

Abstract Background The association between M segment splicing and pathogenicity remains ambiguous in human influenza A viruses. In this study, we aimed to investigate M splicing in various human influenza A viruses and characterize its physiological roles by applying the splicing inhibitor, herboxidiene. Methods We examined the M splicing of human H1N1 and H3N2 viruses by comparing three H1N1 and H3N2 strains, respectively, through reverse transcriptase-polymerase chain reaction (RT-PCR) analyses. We randomly selected M sequences of human H1N1, H2N2, and H3N2 viruses isolated from 1933 to 2020 and examined their phylogenetic relationships. Next, we determined the effects of single nucleotide variations on M splicing by generating mutant viruses harboring the 55C/T variant through reverse genetics. To confirm the importance of M2 splicing in the replication of H1N1 and H3N2, we treated infected cells with splicing inhibitor herboxidiene and analyzed the viral growth using plaque assay. To explore the physiological role of the various levels of M2 protein in pathogenicity, we challenged C57BL/6 mice with the H1N1 WSN wild-type strain, mutant H1N1 (55T), and chimeric viruses including H1N1 + H3wt and H1N1 + H3mut. One-tailed paired t-test was used for virus titer calculation and multiple comparisons between groups were performed using two-way analysis of variance. Results M sequence splice site analysis revealed an evolutionarily conserved single nucleotide variant C55T in H3N2, which impaired M2 expression and was accompanied by collinear M1 and mRNA3 production. Aberrant M2 splicing resulted from splice-site selection rather than a general defect in the splicing process. The C55T substitution significantly reduced both M2 mRNA and protein levels regardless of the virus subtype. Consequently, herboxidiene treatment dramatically decreased both the H1N1 and H3N2 virus titers. However, a lower M2 expression only attenuated H1N1 virus replication and in vivo pathogenicity. This attenuated phenotype was restored by M replacement of H3N2 M in a chimeric H1N1 virus, despite low M2 levels. Conclusions The discrepancy in M2-dependence emphasizes the importance of M2 in human influenza A virus pathogenicity, which leads to subtype-specific evolution. Our findings provide insights into virus adaptation processes in humans and highlights splicing regulation as a potential antiviral target.

Abstract Autoantibodies against interferon-gamma (AutoAbs-IFN-γ) can cause the immunodeficiency condition following various opportunistic infections. Gut microbiota can affect the human immune system in many ways. Many studies have shown that gut dysbiosis was associated with some immune diseases, such as autoimmune diseases and human immunodeficiency virus (HIV) infection, while its relationship at anti-IFN-γ AAbs remains unknown. We aimed to identify the anti-IFN-γ AAbs specific microbiome and the possible association with immunodeficiency. We profiled fecal microbiome for two cohorts of forty subjects, including seven patients with anti-IFN-γ AAbs and 33 individuals with competent immune. The study shows that patients with anti-IFN-γ AAbs have characterized the gut microbiome and have lower alpha diversity indexes than healthy controls (HC). There are significant differences in the microbiome structure at both the family and genera level between the two cohorts. The anti-IFN-γ AAbs cohort featured some microbiome such as Clostridium, including the possible opportunistic pathogen and fewer genera including Bacteroides, Ruminococcus, and Faecalibacterium, some of them with possible immune-related genera. The PICRUSt2 pathway demonstrated the decreased abundance of some immune-related pathways and one potential pathway related to the immune alternations in the anti- IFN-γ AAbs cohort. This was the first study to examine the gut microbiome characteristics in patients with anti-IFN-γ AAbs. It could be involved in the pathogenesis of anti-IFN-γ AAbs and contribute to the derived immune condition in this disease. This could lead to new strategies for treating and preventing patients suffering from this disease.

Abstract Background Streptococcus pneumoniae is a common cause of post-influenza secondary bacterial infection, which results in excessive morbidity and mortality. Although 13-valent pneumococcal conjugate vaccine (PCV13) vaccination programs have decreased the incidence of pneumococcal pneumonia, PCV13 failed to prevent serotype 3 pneumococcal disease as effectively as other vaccine serotypes. We aimed to investigate the mechanisms underlying the co-pathogenesis of influenza virus and serotype 3 pneumococci. Methods We carried out a genome-wide screening of a serotype 3 S. pneumoniae transposon insertion mutant library in a mouse model of coinfection with influenza A virus (IAV) to identify the bacterial factors required for this synergism. Results Direct, high-throughput sequencing of transposon insertion sites identified 24 genes required for both coinfection and bacterial infection alone. Targeted deletion of the putative aminotransferase (PA) gene decreased bacterial growth, which was restored by supplementation with methionine. The bacterial burden in a coinfection with the PA gene deletion mutant and IAV in the lung was lower than that in a coinfection with wild-type pneumococcus and IAV, but was significantly higher than that in an infection with the PA gene deletion mutant alone. These data suggest that IAV infection alters host metabolism to benefit pneumococcal fitness and confer higher susceptibility to pneumococcal infection. We further demonstrated that bacterial growth was increased by supplementation with methionine or IAV-infected mouse lung homogenates. Conclusions The data indicates that modulation of host metabolism during IAV infection may serve as a potential therapeutic intervention against secondary bacterial infections caused by serotype 3 pneumococci during IAV outbreaks in the future.

Influenza A virus (IAV) has caused recurrent epidemics and severe pandemics. In this study, we adapted an MS2-MCP live-cell imaging system to visualize IAV replication. A reporter plasmid, pHH-PB2-vMSL, was constructed by replacing a part of the PB2-coding sequence in pHH-PB2 with a sequence encoding 24 copies of a stem-loop structure from bacteriophage MS2 (MSL). Binding of MS2 coat protein (MCP) fused to green fluorescent protein (GFP) to MSL enabled the detection of vRNA as fluorescent punctate signals in live-cell imaging. The introduction of pHH-PB2-vMSL into A549 cells transduced to express an MCP-GFP fusion protein lacking the nuclear localization signal (MCP-GFPdN), subsequently allowed tracking of the distribution and replication of PB2-vMSL vRNA after IAV PR8 infection. Spatial and temporal measurements revealed exponential increases in vRNA punctate signal intensity, which was only observed after membrane blebbing in apoptotic cells. Similar signal intensity increases in apoptotic cells were also observed after MDCK cells, transduced to express MCP-GFPdN, were infected with IAV carrying PB2-vMSL vRNA. Notably, PB2-vMSL vRNA replication was observed to occur only in apoptotic cells, at a consistent time after apoptosis initiation. There was a lack of observable PB2-vMSL vRNA replication in non-apoptotic cells, and vRNA replication was suppressed in the presence of apoptosis inhibitors. These findings point to an important role for apoptosis in IAV vRNA replication. The utility of the MS2-imaging system for visualizing time-sensitive processes such as viral replication in live host cells is also demonstrated in this study.

Incidence of invasive pneumococcal disease caused by antimicrobial-resistant Streptococcus pneumoniae types not included in pneumococcal conjugate vaccines has increased, including a penicillin- and meropenem-resistant serotype 15A-ST63 clone in Japan. During 2013–2017, we collected 206 invasive pneumococcal isolates in Taiwan for penicillin and meropenem susceptibility testing. We found serotypes 15B/C-ST83 and 15A-ST63 were the most prevalent penicillin- and meropenem-resistant clones. A transformation study confirmed that penicillin-binding protein (PBP) 2b was the primary meropenem resistance determinant, and PBP1a was essential for high-level resistance. The rate of serotype 15B/C-ST83 increased during the study. All 15B/C-ST83 isolates showed an ermB macrolide resistance genotype. Prediction analysis of recombination sites revealed 12 recombination regions in 15B/C-ST83 compared with the S. pneumoniae Spain23F-ST81 genome. Pneumococcal clones rapidly recombine to acquire survival advantages and undergo local expansion under the selective pressure exerted by vaccines and antimicrobial drugs. The spread of 15B/C-ST83 is alarming for countries with high antimicrobial pressure.

ABSTRACT Enterovirus infections can cause severe complications, such as poliomyelitis, encephalitis, myocarditis, meningitis, neurological pulmonary edema, and even death. Here, we used genome-wide CRISPR screens to gain new insight into the mechanism by which enteroviruses co-opt host pathways to potentiate replication and propagation. We found that acyl-coenzyme A synthetase long-chain family member 4 (ACSL4) is involved in viral replication organelle formation. ACSL4 is a key component of ferroptosis, an iron-dependent, nonapoptotic programmed cell death. Our results indicated that enteroviruses and coronaviruses can induce ferroptosis via ACSL4. Most importantly, ferroptosis inhibitors, including two FDA-approved drugs, rosiglitazone (ROSI; ACSL4 inhibitor) and pioglitazone (PIO; ACSL4 inhibitor), decreased the viral load of human enteroviruses and coronaviruses, suggesting that ACSL4 is a target for counteracting viral infection. IMPORTANCE We provide the first evidence for the role of ACSL4 in enterovirus replication organelle formation. Moreover, both enteroviruses and coronaviruses induce ferroptosis via ACSL4. These findings establish a novel regulatory mechanism for viral replication. The inhibition of ACSL4 by ferroptosis inhibitors can reduce viral yields of enteroviruses and coronaviruses, including SARS-CoV-2, implying that ACSL4-mediated ferroptosis is a promising therapeutic target for viral diseases.

ABSTRACT Considering the urgent demand for faster methods to quantify neutralizing antibody titers in patients with coronavirus (CoV) disease 2019 (COVID-19), developing an analytical model or method to replace the conventional virus neutralization test (NT) is essential. Moreover, a “COVID-19 immunity passport” is currently being proposed as a certification for people who travel internationally. Therefore, an enzyme-linked immunosorbent assay (ELISA) was designed to detect severe acute respiratory syndrome CoV 2 (SARS-CoV-2)-neutralizing antibodies in serum, which is based on the binding affinity of SARS-CoV-2 viral spike protein 1 (S1) and the viral spike protein receptor-binding domain (RBD) to antibodies. The RBD is considered the major binding region of neutralizing antibodies. Furthermore, S1 covers the RBD and several other regions, which are also important for neutralizing antibody binding. In this study, we assessed 144 clinical specimens, including those from patients with PCR-confirmed SARS-CoV-2 infections and healthy donors, using both the NT and ELISA. The ELISA results analyzed by spline regression and the two-variable generalized additive model precisely reflected the NT value, and the correlation between predicted and actual NT values was as high as 0.917. Therefore, our method serves as a surrogate to quantify neutralizing antibody titer. The analytic method and platform used in this study present a new perspective for serological testing of SARS-CoV-2 infection and have clinical potential to assess vaccine efficacy. IMPORTANCE Herein, we present a new approach for serological testing for SARS-CoV-2 antibodies using innovative laboratory methods that demonstrate a combination of biology and mathematics. The traditional virus neutralization test is the gold standard method; however, it is time-consuming and poses a risk to medical personnel. Thus, there is a demand for methods that rapidly quantify neutralizing antibody titers in patients with COVID-19 or examine vaccine efficacy at a biosafety level 2 containment facility. Therefore, we used a two-variable generalized additive model to analyze the results of the enzyme-linked immunosorbent assay and found the method to serve as a surrogate to quantify neutralizing antibody titers. This methodology has potential for clinical use in assessing vaccine efficacy.

ABSTRACTTaiwan experienced two waves of imported infections with Coronavirus Disease 2019 (COVID-19). This study aimed at investigating the genomic variation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Taiwan and compared their evolutionary trajectories with the global strains. We performed culture and full-genome sequencing of SARS-CoV-2 strains followed by phylogenetic analysis. A 382-nucleotides deletion in open reading frame 8 (ORF8) was found in a Taiwanese strain isolated from a patient on February 4, 2020 who had a travel history to Wuhan. Patients in the first wave also included several sporadic, local transmission cases. Genomes of 5 strains sequenced from clustered infections were classified into a new clade with ORF1ab-V378I mutation, in addition to 3 dominant clades ORF8-L84S, ORF3a-G251V and S-D614G. This highlighted clade also included some strains isolated from patients who had a travel history to Turkey and Iran. The second wave mostly resulted from patients who had a travel history to Europe and Americas. All Taiwanese viruses were classified into various clades. Genomic surveillance of SARS-CoV-2 in Taiwan revealed a new ORF8-deletion mutant and a virus clade that may be associated with infections in the Middle East, which contributed to a better understanding of the global SARS-CoV-2 transmission dynamics.

Background: Human rhinovirus (HRV) can cause severe illnesses in hospitalized patients. However, there are no studies regarding the prevalence of HRV infection, particularly the recently identified HRV-C, in hospitalized patients reported from Taiwan. Methods: Respiratory specimens collected from 487 hospitalized patients in designated wards between 2013 and 2014 in a medical center in northern Taiwan were retrospectively detected for HRV. Positive specimens were further determined for genotyping. Medical charts of the HRV-positive patients were reviewed retrospectively. Results: Totally, 76 patients (15.6%) were HRV positive, of which 60 were pediatric patients. HRV-A was identified in 41 (54%) patients, HRV-B in 6 patients (7.9%) and HRV-C in 29 patients (38%). A total of 47 different genotypes were identified. HRV infections were predominant during fall and winter seasons. 21.1% were affected by HRV alone and 78.9% were found to be co-infected with other microorganisms. The detection rate of HRV in children (18.6%) was significantly higher than in adults (9.6%). Compared with pediatric patients, adult patients were significantly associated with underlying disease, Pneumocystis jirovesii pneumonia co-infection, a diagnosis of pneumonia, fatal outcome, hospital acquisition of HRV, antibiotics administration and requiring intensive care, while pediatric patients were significantly associated with viral co-infection. Conclusions: HRV was a common cause of respiratory tract infection in Taiwan, particularly in pediatric patients. Eighty percent of HRV-infected inpatients had other microorganisms co-infection. Adult patients were more likely to be associated with a severe respiratory disease entity. Keywords: Rhinovirus, Respiratory tract infection, Hospitalized patients, Taiwan

It has been observed that variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are more stable and infectious after storage at −20°C than at 4°C. A SARS-CoV-2 S-D614G variant is currently the most dominant variant in circulation and is associated with enhanced infectivity.

Background/Purpose: Influenza A/H3N2 viruses are characterized by highly mutated RNA genomes. In this study, we focused on tracing the phylodynamics of Taiwanese strains over the past four decades. Methods: All Taiwanese H3N2 HA1 sequences and references were downloaded from public database. A Bayesian skyline plot (BSP) and phylogenetic tree were used to analyze the evolutionary history, and Bayesian phylogeographic analysis was applied to predict the spatiotemporal migrations of influenza outbreaks. Results: Genetic diversity was found to have peaked near the summer of 2009 in BSP, in addition to the two earlier reported ones in summer of 2005 and 2007. We predicted their spatiotemporal migrations and found the summer epidemic of 2005 from Korea, and 2007 and 2009 from the Western United States. BSP also predicted an elevated genetic diversity in 2015–2017. Quasispecies were found over approximately 20% of the strains included in this time span. In addition, a first-time seen N31S mutation was noted in Taiwan in 2016–2017. Conclusion: We comprehensively investigated the evolutionary history of Taiwanese strains in 1979–2017. An epidemic caution could thus be raised if genetic diversity was found to have peaked. An example showed a newly-discovered cluster in 2016–2017 strains featuring a mutation N31S together with HA-160 quasispecies. Phylogeographic analysis, moreover, provided useful insights in tracing the possible source and migrations of these epidemics around the world. We demonstrated that Asian destinations including Taiwan were the immediate followers, while U.S. continent was predicted the origin of two summer epidemics in 2007 and 2009. Keywords: Influenza A/H3N2 viruses, Phylodynamics, Phylogeography, Summer epidemics, Viral quasispecies

Enterovirus-induced infection of the central nervous system (CNS) results in acute inflammation of the brain (encephalitis) and constitutes a significant global burden to human health. These viruses are thought to be highly cytolytic, therefore normal brain function could be greatly compromised following enteroviral infection of the CNS. A further layer of complexity is added by evidence showing that some enteroviruses may establish a persistent infection within the CNS and eventually lead to pathogenesis of certain neurodegenerative disorders. Interestingly, enterovirus encephalitis is particularly common among young children, suggesting a potential causal link between the development of the neuroimmune system and enteroviral neuroinvasion. Although the CNS involvement in enterovirus infections is a relatively rare complication, it represents a serious underlying cause of mortality. Here we review a selection of enteroviruses that infect the CNS and discuss recent advances in the characterization of these enteroviruses with regard to their routes of CNS infection, tropism, virulence, and immune responses.

Abstract A novel pandemic influenza A(H1N1)pdm09 virus emerged in 2009 globally, and it continues to circulate in humans. The National Influenza Surveillance Network in Taiwan identified five A(H1N1)pdm09-predominant seasons, representing the 2009/2010, 2010/2011, 2012/2013, 2013/2014, and 2015/2016 outbreaks from 2009 to 2016. Independently, a retrospective cohort study (which enrolled 639 infected patients during the five seasons) was conducted at Chang Gung Memorial Hospital to explore the risk factors associated with influenza A(H1N1)pdm09-related complications. A phylogenetic analysis of hemagglutinin (HA) sequences showed that the circulating A(H1N1)pdm09 virus belonged to clades 1, 2, and 8 in 2009/2010; clades 3, 4, 5, and 7 in 2010/2011; clades 7 and 6C in 2012/2013; clades 6B in 2013/2014; and 6B/6B.1/6B.2 in 2015/2016. Compared to individuals infected in non-6B/6B.1/6B.2 seasons (2009/2010, 2010/2011, and 2012/2013), those infected in 6B/6B.1/6B.2 seasons (2013/2014 and 2015/2016) were at higher risk for influenza-related complications (adjusted odds ratio [aOR]: 1.6, 95% confidence interval [CI]: 1.0–2.8), pneumonia (aOR: 1.78, 95% CI: 1.04–3.04), mechanical ventilation (aOR: 2.6, 95% CI: 1.2–5.6), and acute respiratory distress syndrome (aOR: 5.5, 95% CI: 1.9–15.9). For the increased severity of infection during the influenza A(H1N1)pdm09 clade 6B/6B.1/6B.2 seasons, aspects related to the antigenic change of A(H1N1)pdm09 virus, immune response of the host, and environmental factors required further investigation.

The history of influenza in Taiwan can be traced up to the 1918 H1N1 Spanish flu pandemic, followed by several others including the 1957 H2N2, 1968 H3N2, and the 2009 new H1N1. A couple of avian influenza viruses of H5N1 and H7N9 also posed threats to the general public in Taiwan in the two recent decades. Nevertheless, two seasonal influenza A viruses and two lineages of influenza B viruses continue causing annual endemics one after the other, or appearing simultaneously. Their interplay provided interesting evolutionary trajectories for these viruses, allowing us to computationally model their global migrations together with the data collected elsewhere from different geographical locations. An island-wide laboratory-based surveillance network was also established since 2000 for systematically collecting and managing the disease and molecular epidemiology. Experiences learned from this network helped in encountering and managing newly emerging infectious diseases, including the 2003 SARS and 2009 H1N1 outbreaks. Keywords: Influenza, Pandemic, Genome evolution, Surveillance

Streptococcus pneumoniae 19A ST320, a multidrug-resistant strain with high disease severity that notoriously spread before the use of expanded pneumococcal conjugate vaccines, was derived from a capsular switching event between an international strain Taiwan 19F-14 (ST236) and a serotype 19A strain. However, the molecular mechanisms underlying the adaptive evolution of 19F ST236 to 19A ST320 are unknown. In this study, we compared 19A ST320 to its ancestral clone, 19F ST236, in terms of adherence to respiratory epithelial cells, whole transcriptome, and ability to colonize a young mouse model. Serotype 19A ST320 showed five-fold higher adherence to A549 cells than serotype 19F ST236. High-throughput mRNA sequencing identified a prophage region located between dnaN and ychF in both strains; however, the genes in this region were expressed at significantly higher levels in 19A ST320 than in 19F ST236. Analysis by polymerase chain reaction (PCR) showed that the prophage is able to spontaneously excise from the chromosome and form a circular episome in 19A ST320, but not in 19F ST236. Deletion of the integrase in the prophage of 19A ST320 decreased spontaneous excision and cell adherence, which were restored by complementation. Competition experiments in mice showed that the integrase mutant was six-fold less competitive than the 19A ST320 parent (competitive index [CI]: 0.16; p = 0.02). The 19A ST320 prophage-deleted strain did not change cell adherence capacity, whereas prophage integration strains (integrase mutant and 19F) had decreased expression of the down-stream ychF gene compared to that of 19A ST320. Further deletion of ychF significantly reduced cell adherence. In conclusions, these findings suggest that spontaneous prophage induction confers a competitive advantage to virulent pneumococci.

Background & Aims: Gut dysbiosis plays a role in hepatic encephalopathy (HE), while its relationship at the acute episode of overt HE (AHE), the disease progression and clinical outcomes remains unclear. We aimed to identify AHE-specific microbiome and its association to patients’ outcomes. Methods: We profiled fecal microbiome changes for a cohort of 62 patients with cirrhosis and AHE i) before treatment, ii) 2-3 days after medication and iii) 2-3 months after recovery, and three control cohorts i) healthy individuals, patients with ii) compensated or iii) decompensated cirrhosis. Results: Comparison of the microbiome shift from compensated, decompensated cirrhosis, AHE to recovery revealed the AHE-specific gut-dysbiosis. The gut microbiome diversity was decreased during AHE, further reduced after medication, and only partially reversed during the recovery. The relative abundance of Bacteroidetes phylum in the microbiome decreased, whereas that of Firmicute, Proteobacteria and Actinobacteria increased in patients during AHE compared with those with compensated cirrhosis. A total of 70 operational taxonomic units (OTUs) were significantly different between AHE and decompensated cirrhosis abundances. Of them, the abundance of Veillonella parvula increased the most during AHE via a metagenomics recovery of the genomes. Moreover, the relative abundances of three (Alistipes, Bacteroides, Phascolarctobacterium) and five OTUs (Clostridium-XI, Bacteroides, Bacteroides, Lactobacillus, Clostridium-sedis) at AHE were respectively associated with HE recurrence and overall survival during the subsequent one-year follow-up. Conclusions: AHE-specific gut OTUs were identified that may be involved in HE development and able to predict clinical outcomes, providing new strategies for the prevention and treatment of HE recurrence in patients with cirrhosis. Keywords: Gut Microbiota, Gut Dysbiosis, Gut-Liver-Brain Axis, Veillonella parvula

An influenza A pandemic occurred in 2009–2010. A novel H1N1 virus (hereafter H1N1pdm) was responsible for this outbreak. H1N1pdm viruses have been largely seen in recent human influenza A viruses. This virus was descended from a triple-reassorted swine virus consisting of human, avian, and swine origins. As a result, the previously established species-associated signatures could be in jeopardy. Methods: We analyzed all influenza A sequences in the past 5 years after the inclusion of H1N1pdm into human viruses since 2009, and examined how human signatures may lose their distinctness by mixing with avian residues that H1N1pdm have brought in. In particular, we compared how those signatures were changed/shifted in the past 5 years for human-isolated avian influenza A viruses and discussed their implications. Results: Only eight out of 47 signatures remained human-like for human influenza A viruses in the past 5 years. They are PB2 271A; PB1 336I; PA 356R and 409N; NP 33I, 305K, and 357K; and NS1 227R. Although most avian-like residues were preserved in human-isolated avian influenza A viruses, a number of them were found to have become or on the verge of becoming human-like, including PB2 627, PA 100, 356, 404, 409, NP 33, 61, 305, 357, M2 20, and NS1 81. Conclusion: Analyzing how species-associated signatures are becoming human-like in human-isolated avian influenza A viruses helps in assessing their potential to go pandemic as well as providing insights into host adaptation.

An avian influenza A H7N9 virus emerged in March 2013 and caused a remarkable number of human fatalities. Genome variability in these viruses may provide insights into host adaptability. We scanned over 140 genomes of the H7N9 viruses isolated from humans and identified 104 positions that exhibited seven or more amino acid substitutions. Approximately half of these substitutions were identified in the influenza ribonucleoprotein (RNP) complex. Although PB2 627K of the avian virus promotes replication in humans, 45 of the 147 investigated PB2 sequences retained the E signature at this position, which is an avian characteristic. We discovered 10 PB2 substitutions that covaried with K627E. An RNP activity assay showed that Q591K, D701N, and M535L restored the polymerase activity in human cells when 627K transformed to an avian-like E. Genomic analysis of the human-isolated avian influenza virus is crucial in assessing genome variability, because relationships between position-specific variations can be observed and explored. In this study, we observed alternative positions that can potentially compensate for PB2 627K, a well-known marker for cross-species infection. An RNP assay suggested Q591K, D701N, and M535L as potential markers for an H7N9 virus capable of infecting humans.

Forty-two cytopathic effect (CPE)-positive isolates were collected from 2008 to 2012. All isolates could not be identified for known viral pathogens by routine diagnostic assays. They were pooled into 8 groups of 5-6 isolates to reduce the sequencing cost. Next-generation sequencing (NGS) was conducted for each group of mixed samples, and the proposed data analysis pipeline was used to identify viral pathogens in these mixed samples. Polymerase chain reaction (PCR) or enzyme-linked immunosorbent assay (ELISA) was individually conducted for each of these 42 isolates depending on the predicted viral types in each group. Two isolates remained unknown after these tests. Moreover, iteration mapping was implemented for each of these 2 isolates, and predicted human parechovirus (HPeV) in both. In summary, our NGS pipeline detected the following viruses among the 42 isolates: 29 human rhinoviruses (HRVs), 10 HPeVs, 1 human adenovirus (HAdV), 1 echovirus and 1 rotavirus. We then focused on the 10 identified Taiwanese HPeVs because of their reported clinical significance over HRVs. Their genomes were assembled and their genetic diversity was explored. One novel 6-bp deletion was found in one HPeV-1 virus. In terms of nucleotide heterogeneity, 64 genetic variants were detected from these HPeVs using the mapped NGS reads. Most importantly, a recombination event was found between our HPeV-3 and a known HPeV-4 strain in the database. Similar event was detected in the other HPeV-3 strains in the same clade of the phylogenetic tree. These findings demonstrated that the proposed NGS data analysis pipeline identified unknown viruses from the mixed clinical samples, revealed their genetic identity and variants, and characterized their genetic features in terms of viral evolution.

The influenza A virus contains 8 segmented genomic RNAs and was considered to encode 10 viral proteins until investigators identified the 11th viral protein, PB1-F2, which uses an alternative reading frame of the PB1 gene. The recently identified PB1-N40, PA-N155 and PA-N182 influenza A proteins have shown the potential for using a leaking ribosomal scanning mechanism to generate novel open reading frames (ORFs). These novel ORFs provide examples of the manner in which the influenza A virus expands its coding capacity by using overlapping reading frames. In this study, we performed a computational search, based on a ribosome scanning mechanism, on all influenza A coding sequences to identify possible forward-reading ORFs that could be translated into novel viral proteins. We specified that the translated products had a prevalence ≥5% to eliminate sporadic ORFs. A total of 1,982 ORFs were thus identified and presented in terms of their locations, lengths and Kozak sequence strengths. We further provided an abridged list of ORFs by requiring every candidate an upstream start codon (within the upstream third of the primary transcript), a strong Kozak consensus sequence and high prevalence (≥95% and ≥50% for in-frame and alternative-frame ORFs, respectively). The PB1-F2, PB1-N40, PA-N155 and PA-N182 proteins all fulfilled our filtering criteria. Subject to these three stringent settings, we additionally named 16 novel ORFs for all influenza A genomes except for HA and NA, for which 43 HA and 11 NA ORFs from their respective subtypes were also recognized.

A swine-origin influenza A was detected in April 2009 and soon became the 2009 H1N1 pandemic strain (H1N1pdm). The current study revealed the genetic diversity of H1N1pdm, based on 77 and 70 isolates which we collected, respectively, during the 2009/2010 and 2010/2011 influenza seasons in Taiwan. We focused on tracking the amino acid transitioning of hemagglutinin (HA) and neuraminidase (NA) genes in the early diversification of the virus and compared them with H1N1pdm strains reported worldwide. We identified newly emerged mutation markers based on A/California/04/2009, described how these markers shifted from the first H1N1pdm season to the one that immediately followed, and discussed how these observations may relate to antigenicity, receptor-binding, and drug susceptibility. It was found that the amino acid mutation rates of H1N1pdm were elevated, from 9.29×10(-3) substitutions per site in the first season to 1.46×10(-2) in the second season in HA, and from 5.23×10(-3) to 1.10×10(-2) in NA. Many mutation markers were newly detected in the second season, including 11 in HA and 8 in NA, and some were found having statistical correlation to disease severity. There were five noticeable HA mutations made to antigenic sites. No significant titer changes, however, were detected based on hemagglutination inhibition tests. Only one isolate with H275Y mutation known to reduce susceptibility to NA inhibitors was detected. As limited Taiwanese H1N1pdm viruses were isolated after our sampling period, we gathered 8,876 HA and 6,017 NA H1N1pdm sequences up to April 2012 from NCBI to follow up the dynamics of mentioned HA mutations. While some mutations described in this study seemed to either settle in or die out in the 2011-2012 season, a number of them still showed signs of transitioning, prompting the importance of continuous monitoring of this virus for more seasons to come.

Incidence of invasive pneumococcal disease caused by antimicrobial-resistant Streptococcus pneumoniae types not included in pneumococcal conjugate vaccines has increased, including a penicillin- and meropenem-resistant serotype 15A-ST63 clone in Japan. During 2013-2017, we collected 206 invasive pneumococcal isolates in Taiwan for penicillin and meropenem susceptibility testing. We found serotypes 15B/C-ST83 and 15A-ST63 were the most prevalent penicillin- and meropenem-resistant clones. A transformation study confirmed that penicillin-binding protein (PBP) 2b was the primary meropenem resistance determinant, and PBP1a was essential for high-level resistance. The rate of serotype 15B/C-ST83 increased during the study. All 15B/C-ST83 isolates showed an ermB macrolide resistance genotype. Prediction analysis of recombination sites revealed 12 recombination regions in 15B/C-ST83 compared with the S. pneumoniae Spain23F-ST81 genome. Pneumococcal clones rapidly recombine to acquire survival advantages and undergo local expansion under the selective pressure exerted by vaccines and antimicrobial drugs. The spread of 15B/C-ST83 is alarming for countries with high antimicrobial pressure. [ABSTRACT FROM AUTHOR]