Amino acid substitution L232F in non-structural protein 6 identified as a possible human-adaptive mutation in clade B MERS coronaviruses.

IMPORTANCE: Viral host adaptation plays an important role in inter-species transmission of coronaviruses and influenza viruses. Multiple human-adaptive mutations have been identified in influenza viruses but not so far in MERS-CoV that circulates widely in dromedary camels in the Arabian Peninsula leading to zoonotic transmission. Here, we analyzed clade B MERS-CoV sequences and identified an amino acid substitution L232F in nsp6 that repeatedly occurs in human MERS-CoV. Using a loss-of-function reverse genetics approach, we found the nsp6 L232F conferred increased viral replication competence in vitro, in cultures of the upper human respiratory tract ex vivo, and in lungs of mice infected in vivo. Our results showed that nsp6 L232F may be an adaptive mutation associated with zoonotic transmission of MERS-CoV. This study highlighted the capacity of MERS-CoV to adapt to transmission to humans and also the need for continued surveillance of MERS-CoV in camels.

ggmsa: a visual exploration tool for multiple sequence alignment and associated data.

The identification of the conserved and variable regions in the multiple sequence alignment (MSA) is critical to accelerating the process of understanding the function of genes. MSA visualizations allow us to transform sequence features into understandable visual representations. As the sequence-structure-function relationship gains increasing attention in molecular biology studies, the simple display of nucleotide or protein sequence alignment is not satisfied. A more scalable visualization is required to broaden the scope of sequence investigation. Here we present ggmsa, an R package for mining comprehensive sequence features and integrating the associated data of MSA by a variety of display methods. To uncover sequence conservation patterns, variations and recombination at the site level, sequence bundles, sequence logos, stacked sequence alignment and comparative plots are implemented. ggmsa supports integrating the correlation of MSA sequences and their phenotypes, as well as other traits such as ancestral sequences, molecular structures, molecular functions and expression levels. We also design a new visualization method for genome alignments in multiple alignment format to explore the pattern of within and between species variation. Combining these visual representations with prime knowledge, ggmsa assists researchers in discovering MSA and making decisions. The ggmsa package is open-source software released under the Artistic-2.0 license, and it is freely available on Bioconductor (https://bioconductor.org/packages/ggmsa) and Github (https://github.com/YuLab-SMU/ggmsa).

Australia as a global sink for the genetic diversity of avian influenza A virus.

Most of our understanding of the ecology and evolution of avian influenza A virus (AIV) in wild birds is derived from studies conducted in the northern hemisphere on waterfowl, with a substantial bias towards dabbling ducks. However, relevant environmental conditions and patterns of avian migration and reproduction are substantially different in the southern hemisphere. Through the sequencing and analysis of 333 unique AIV genomes collected from wild birds collected over 15 years we show that Australia is a global sink for AIV diversity and not integrally linked with the Eurasian gene pool. Rather, AIV are infrequently introduced to Australia, followed by decades of isolated circulation and eventual extinction. The number of co-circulating viral lineages varies per subtype. AIV haemagglutinin (HA) subtypes that are rarely identified at duck-centric study sites (H8-12) had more detected introductions and contemporary co-circulating lineages in Australia. Combined with a lack of duck migration beyond the Australian-Papuan region, these findings suggest introductions by long-distance migratory shorebirds. In addition, on the available data we found no evidence of directional or consistent patterns in virus movement across the Australian continent. This feature corresponds to patterns of bird movement, whereby waterfowl have nomadic and erratic rainfall-dependant distributions rather than consistent intra-continental migratory routes. Finally, we detected high levels of virus gene segment reassortment, with a high diversity of AIV genome constellations across years and locations. These data, in addition to those from other studies in Africa and South America, clearly show that patterns of AIV dynamics in the Southern Hemisphere are distinct from those in the temperate north.

The Predictive Value of Gut Microbiota Composition for Sustained Immunogenicity following Two Doses of CoronaVac.

CoronaVac immunogenicity decreases with time, and we aimed to investigate whether gut microbiota associate with longer-term immunogenicity of CoronaVac. This was a prospective cohort study recruiting two-dose CoronaVac recipients from three centres in Hong Kong. We collected blood samples at baseline and day 180 after the first dose and used chemiluminescence immunoassay to test for neutralizing antibodies (NAbs) against the receptor-binding domain (RBD) of wild-type SARS-CoV-2 virus. We performed shotgun metagenomic sequencing performed on baseline stool samples. The primary outcome was the NAb seroconversion rate (seropositivity defined as NAb ≥ 15AU/mL) at day 180. Linear discriminant analysis [LDA] effect size analysis was used to identify putative bacterial species and metabolic pathways. A univariate logistic regression model was used to derive the odds ratio (OR) of seropositivity with bacterial species. Of 119 CoronaVac recipients (median age: 53.4 years [IQR: 47.8-61.3]; male: 39 [32.8%]), only 8 (6.7%) remained seropositive at 6 months after vaccination. Bacteroides uniformis (log10LDA score = 4.39) and Bacteroides eggerthii (log10LDA score = 3.89) were significantly enriched in seropositive than seronegative participants. Seropositivity was associated with B. eggerthii (OR: 5.73; 95% CI: 1.32-29.55; p = 0.022) and B. uniformis with borderline significance (OR: 3.27; 95% CI: 0.73-14.72; p = 0.110). Additionally, B. uniformis was positively correlated with most enriched metabolic pathways in seropositive vaccinees, including the superpathway of adenosine nucleotide de novo biosynthesis I (log10LDA score = 2.88) and II (log10LDA score = 2.91), as well as pathways related to vitamin B biosynthesis, all of which are known to promote immune functions. In conclusion, certain gut bacterial species (B. eggerthii and B. uniformis) and metabolic pathways were associated with longer-term CoronaVac immunogenicity.

Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia.

BACKGROUND: The initial cases of novel coronavirus (2019-nCoV)-infected pneumonia (NCIP) occurred in Wuhan, Hubei Province, China, in December 2019 and January 2020. We analyzed data on the first 425 confirmed cases in Wuhan to determine the epidemiologic characteristics of NCIP. METHODS: We collected information on demographic characteristics, exposure history, and illness timelines of laboratory-confirmed cases of NCIP that had been reported by January 22, 2020. We described characteristics of the cases and estimated the key epidemiologic time-delay distributions. In the early period of exponential growth, we estimated the epidemic doubling time and the basic reproductive number. RESULTS: Among the first 425 patients with confirmed NCIP, the median age was 59 years and 56% were male. The majority of cases (55%) with onset before January 1, 2020, were linked to the Huanan Seafood Wholesale Market, as compared with 8.6% of the subsequent cases. The mean incubation period was 5.2 days (95% confidence interval [CI], 4.1 to 7.0), with the 95th percentile of the distribution at 12.5 days. In its early stages, the epidemic doubled in size every 7.4 days. With a mean serial interval of 7.5 days (95% CI, 5.3 to 19), the basic reproductive number was estimated to be 2.2 (95% CI, 1.4 to 3.9). CONCLUSIONS: On the basis of this information, there is evidence that human-to-human transmission has occurred among close contacts since the middle of December 2019. Considerable efforts to reduce transmission will be required to control outbreaks if similar dynamics apply elsewhere. Measures to prevent or reduce transmission should be implemented in populations at risk. (Funded by the Ministry of Science and Technology of China and others.).

Toward a Universal Unit for Quantification of Antibiotic Resistance Genes in Environmental Samples.

Surveillance of antibiotic resistance genes (ARGs) has been increasingly conducted in environmental sectors to complement the surveys in human and animal sectors under the "One-Health" framework. However, there are substantial challenges in comparing and synthesizing the results of multiple studies that employ different test methods and approaches in bioinformatic analysis. In this article, we consider the commonly used quantification units (ARG copy per cell, ARG copy per genome, ARG density, ARG copy per 16S rRNA gene, RPKM, coverage, PPM, etc.) for profiling ARGs and suggest a universal unit (ARG copy per cell) for reporting such biological measurements of samples and improving the comparability of different surveillance efforts.

Recombinant BA.1/BA.2 SARS-CoV-2 Virus in Arriving Travelers, Hong Kong, February 2022.

We studied SARS-CoV-2 genomes from travelers arriving in Hong Kong during November 2021-February 2022. In addition to Omicron and Delta variants, we detected a BA.1/BA.2 recombinant with a breakpoint near the 5' end of the spike gene in 2 epidemiologically linked case-patients. Continued surveillance for SARS-CoV-2 recombinants is needed.

Surveillance of Rodent Pests for SARS-CoV-2 and Other Coronaviruses, Hong Kong.

We report surveillance conducted in 217 pestiferous rodents in Hong Kong for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We did not detect SARS-CoV-2 RNA but identified 1 seropositive rodent, suggesting exposure to a virus antigenically similar to SARS-CoV-2. Potential exposure of urban rodents to SARS-CoV-2 cannot be ruled out.

Incorporating heterogeneous sampling probabilities in continuous phylogeographic inference - Application to H5N1 spread in the Mekong region.

MOTIVATION: The potentially low precision associated with the geographic origin of sampled sequences represents an important limitation for spatially explicit (i.e. continuous) phylogeographic inference of fast-evolving pathogens such as RNA viruses. A substantial proportion of publicly available sequences is geo-referenced at broad spatial scale such as the administrative unit of origin, rather than more precise locations (e.g. geographic coordinates). Most frequently, such sequences are either discarded prior to continuous phylogeographic inference or arbitrarily assigned to the geographic coordinates of the centroid of their administrative area of origin for lack of a better alternative. RESULTS: We here implement and describe a new approach that allows to incorporate heterogeneous prior sampling probabilities over a geographic area. External data, such as outbreak locations, are used to specify these prior sampling probabilities over a collection of sub-polygons. We apply this new method to the analysis of highly pathogenic avian influenza H5N1 clade data in the Mekong region. Our method allows to properly include, in continuous phylogeographic analyses, H5N1 sequences that are only associated with large administrative areas of origin and assign them with more accurate locations. Finally, we use continuous phylogeographic reconstructions to analyse the dispersal dynamics of different H5N1 clades and investigate the impact of environmental factors on lineage dispersal velocities. AVAILABILITY AND IMPLEMENTATION: Our new method allowing heterogeneous sampling priors for continuous phylogeographic inference is implemented in the open-source multi-platform software package BEAST 1.10. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genomic Analysis of the Emergence, Evolution, and Spread of Human Respiratory RNA Viruses.

The emergence and reemergence of rapidly evolving RNA viruses-particularly those responsible for respiratory diseases, such as influenza viruses and coronaviruses-pose a significant threat to global health, including the potential of major pandemics. Importantly, recent advances in high-throughput genome sequencing enable researchers to reveal the genomic diversity of these viral pathogens at much lower cost and with much greater precision than they could before. In particular, the genome sequence data generated allow inferences to be made on the molecular basis of viral emergence, evolution, and spread in human populations in real time. In this review, we introduce recent computational methods that analyze viral genomic data, particularly in combination with metadata such as sampling time, geographic location, and virulence. We then outline the insights these analyses have provided into the fundamental patterns and processes of evolution and emergence in human respiratory RNA viruses, as well as the major challenges in such genomic analyses.

The impact of host genetic diversity on virus evolution and emergence.

Accumulating evidence indicates that biodiversity has an important impact on parasite evolution and emergence. The vast majority of studies in this area have only considered the diversity of species within an environment as an overall measure of biodiversity, overlooking the role of genetic diversity within a particular host species. Although theoretical models propose that host genetic diversity in part shapes that of the infecting parasite population, and hence modulates the risk of parasite emergence, this effect has seldom been tested empirically. Using Rabies virus (RABV) as a model parasite, we provide evidence that greater host genetic diversity increases both parasite genetic diversity and the likelihood of a host being a donor in RABV cross-species transmission events. We conclude that host genetic diversity may be an important determinant of parasite evolution and emergence.

History of Swine influenza viruses in Asia.

The pig is one of the main hosts of influenza A viruses and plays important roles in shaping the current influenza ecology. The occurrence of the 2009 H1N1 pandemic influenza virus demonstrated that pigs could independently facilitate the genesis of a pandemic influenza strain. Genetic analyses revealed that this virus was derived by reassortment between at least two parent swine influenza viruses (SIV), from the northern American triple reassortant H1N2 (TR) and European avian-like H1N1 (EA) lineages. The movement of live pigs between different continents and subsequent virus establishment are preconditions for such a reassortment event to occur. Asia, especially China, has the largest human and pig populations in the world, and seems to be the only region frequently importing pigs from other continents. Virological surveillance revealed that not only classical swine H1N1 (CS), and human-origin H3N2 viruses circulated, but all of the EA, TR and their reassortant variants were introduced into and co-circulated in pigs in this region. Understanding the long-term evolution and history of SIV in Asia would provide insights into the emergence of influenza viruses with epidemic potential in swine and humans.

Full-genome deep sequencing and phylogenetic analysis of novel human betacoronavirus.

A novel betacoronavirus associated with lethal respiratory and renal complications was recently identified in patients from several countries in the Middle East. We report the deep genome sequencing of the virus directly from a patient's sputum sample. Our high-throughput sequencing yielded a substantial depth of genome sequence assembly and showed the minority viral variants in the specimen. Detailed phylogenetic analysis of the virus genome (England/Qatar/2012) revealed its close relationship to European bat coronaviruses circulating among the bat species of the Vespertilionidae family. Molecular clock analysis showed that the 2 human infections of this betacoronavirus in June 2012 (EMC/2012) and September 2012 (England/Qatar/2012) share a common virus ancestor most likely considerably before early 2012, suggesting the human diversity is the result of multiple zoonotic events.

Emergence and dissemination of a swine H3N2 reassortant influenza virus with 2009 pandemic H1N1 genes in pigs in China.

The 2009 pandemic influenza virus (pdm/09) has been frequently introduced to pigs and has reassorted with other swine viruses. Recently, H3N2 reassortants with pdm/09-like internal genes were isolated in Guangxi and Hong Kong, China. Genetic and epidemiological analyses suggest that these viruses have circulated in swine for some time. This is the first evidence that swine reassortant viruses with pdm/09-like genes may have become established in the field, altering the landscape of human and swine influenza.

Effect of probiotic bacteria on the intestinal microbiota in irritable bowel syndrome.

BACKGROUND AND AIM: In irritable bowel syndrome (IBS), the gut microbiota may be altered. Probiotic bacteria appear to be therapeutically effective. We characterized the mucosa-associated microbiota, and determined the clinical and microbiological effects of orally administered probiotic bacteria, in patients with IBS. METHODS: Mucosal microbiota from rectal biopsies of IBS patients and controls were assessed on the V1 and V2 variable regions of the 16S ribosomal RNA gene amplified using 454 pyrosequencing. Clinical symptoms and changes in mucosal microbiota were assessed in IBS patients before and after 4 weeks of treatment with probiotic mix VSL#3. RESULTS: Ten IBS subjects (eight female; mean age 46 years) were included. At week 4 of probiotic therapy, six patients showed symptom improvement on global symptom assessment compared with baseline (P = 0.031). Before therapy, intestinal microbiota of IBS subjects differed significantly from that of healthy controls, with less diversity and evenness than controls (n = 9; P < 0.05), increased abundance of Bacteroidetes (P = 0.014) and Synegitestes (P = 0.017), and reduced abundance of Actinobacteria (P = 0.004). The classes Flavobacteria (P = 0.028) and Epsilonproteobacteria (P = 0.017) were less enriched in IBS. Abundance differences were largely consistent from the phylum to genus level. Probiotic treatment in IBS patients was associated with a significant reduction of the genus Bacteroides (all taxonomy levels; P < 0.05) to levels similar to that of controls. CONCLUSION: In this pilot study, global and deep molecular analysis demonstrates an altered mucosal microbiota composition in IBS. Probiotic leads to detectable changes in the microbiota. These effects of probiotic bacteria may contribute to their therapeutic benefit.

Within-host genetic diversity of SARS-CoV-2 lineages in unvaccinated and vaccinated individuals.

Viral and host factors can shape SARS-CoV-2 evolution. However, little is known about lineage-specific and vaccination-specific mutations that occur within individuals. Here, we analysed deep sequencing data from 2,820 SARS-CoV-2 respiratory samples with different viral lineages to describe the patterns of within-host diversity under different conditions, including vaccine-breakthrough infections. In unvaccinated individuals, variant of Concern (VOC) Alpha, Delta, and Omicron respiratory samples were found to have higher within-host diversity and were under neutral to purifying selection at the full genome level compared to non-VOC SARS-CoV-2. Breakthrough infections in 2-dose or 3-dose Comirnaty and CoronaVac vaccinated individuals did not increase levels of non-synonymous mutations and did not change the direction of selection pressure. Vaccine-induced antibody or T cell responses did not appear to have significant impact on within-host SARS-CoV-2 sequence diversification. Our findings suggest that vaccination does not increase exploration of SARS-CoV-2 protein sequence space and may not facilitate emergence of viral variants.

One Health Surveillance of Antimicrobial Use and Resistance: Challenges and Successes of Implementing Surveillance Programs in Sri Lanka.

BACKGROUND: Sri Lanka is a low-income country, as defined by the World Bank. The country suffered further economic downturn during the COVID-19 pandemic. This situation adversely affected the prioritization of policies and programs around healthcare and public health. In particular, inflation, fuel prices, and shortage of food supplies increased struggles to implement antimicrobial resistance (AMR) programs. However, in the long run, it is crucial to gather data and evidence to plan AMR policies and track interventions. (1) Aim: To establish and reiterate the importance of prioritizing AMR programs in the One Health framework, the Fleming Fellows collected and studied antimicrobial use/consumption (AMU/AMC) and resistance (AMR) in humans, food-producing animals, and the environment. (2) Methods: A systematic and cross-sectional study was conducted between 2019 and 2021. By way of coordinating an AMU/AMC and AMR prevalence study across six agencies from human health and food-producing animal sectors, the authors established a field epidemiology study, laboratory testing, and data processing at their institutions. AMU/AMC patterns were surveyed using questionnaires and interviews, while AMR samples were collected for antibiotic susceptibility tests and genomic tests. Samples were tested for phenotypic and genotypic resistance. (3) Results: In human samples, resistance was highest to beta-lactam antibiotics. In non-human samples, resistance was highest to erythromycin, a highest-priority, critically important antibiotic defined by the World Health Organization. From government records, tylosin was sold the most in the food-producing animal sector. (4) Conclusions: Sri Lanka AMU and AMR trends in human and non-human sectors can be ascertained by a One Health framework. Further coordinated, consistent, and sustainable planning is feasible, and can help implement an AMU/AMR surveillance system in Sri Lanka.

Discovery of novel Mamastroviruses in Bactrian camels and dromedaries reveals complex recombination history.

Virus emergence may occur through interspecies transmission and recombination of viruses coinfecting a host, with potential to pair novel and adaptive gene combinations. Camels are known to harbor diverse ribonucleic acid viruses with zoonotic and epizootic potential. Among them, astroviruses are of particular interest due to their cross-species transmission potential and endemicity in diverse host species, including humans. We conducted a molecular epidemiological survey of astroviruses in dromedaries from Saudi Arabia and Bactrian camels from Inner Mongolia, China. Herein, we deployed a hybrid sequencing approach coupling deep sequencing with rapid amplification of complementary deoxyribonucleic acid ends to characterize two novel Bactrian and eight dromedary camel astroviruses, including both partial and complete genomes. Our reported sequences expand the known diversity of dromedary camel astroviruses, highlighting potential recombination events among the astroviruses of camelids and other host species. In Bactrian camels, we detected partially conserved gene regions bearing resemblance to human astrovirus types 1, 4, and 8 although we were unable to recover complete reading frames from these samples. Continued surveillance of astroviruses in camelids, particularly Bactrian species and associated livestock, is highly recommended to identify patterns of cross-species transmission and to determine any epizootic threats and zoonotic risks posed to humans. Phylogenomic approaches are needed to investigate complex patterns of recombination among the astroviruses and to infer their evolutionary history across diverse host species.

Novel reassortment of Eurasian avian-like and pandemic/2009 influenza viruses in swine: infectious potential for humans.

Pigs are considered to be intermediate hosts and "mixing vessels," facilitating the genesis of pandemic influenza viruses, as demonstrated by the emergence of the 2009 H1N1 pandemic (pdm/09) virus. The prevalence and repeated introduction of the pdm/09 virus into pigs raises the possibility of generating novel swine influenza viruses with the potential to infect humans. To address this, an active influenza surveillance program was conducted with slaughtered pigs in abattoirs in southern China. Over 50% of the pigs tested were found to be seropositive for one or more H1 influenza viruses, most commonly pdm/09-like viruses. Out of 36 virus isolates detected, one group of novel reassortants had Eurasian avian-like swine H1N1 surface genes and pdm/09 internal genes. Animal experiments showed that this virus transmitted effectively from pig to pig and from pig to ferret, and it could also replicate in ex vivo human lung tissue. Immunization against the 2009 pandemic virus gave only partial protection to ferrets. The continuing prevalence of the pdm/09 virus in pigs could lead to the genesis of novel swine reassortant viruses with the potential to infect humans.

Ggtree: A serialized data object for visualization of a phylogenetic tree and annotation data.

While phylogenetic trees and associated data have been getting easier to generate, it has been difficult to reuse, combine, and synthesize the information they provided, because published trees are often only available as image files and associated data are often stored in incompatible formats. To increase the reproducibility and reusability of phylogenetic data, the ggtree object was designed for storing phylogenetic tree and associated data, as well as visualization directives. The ggtree object itself is a graphic object and can be rendered as a static image. More importantly, the input tree and associated data that are used in visualization can be extracted from the graphic object, making it an ideal data structure for publishing tree (image, tree, and data in one single object) and thus enhancing data reuse and analytical reproducibility, as well as facilitating integrative and comparative studies. The ggtree package is freely available at https://www.bioconductor.org/packages/ggtree.