Cross-border spread, lineage displacement and evolutionary rate estimation of rabies virus in Yunnan Province, China

Examining the molecular clock hypothesis for the contemporary evolution of the rabies virus

The molecular clock hypothesis assumes that mutations accumulate on an organism's genome at a constant rate over time, but this assumption does not always hold true. While modelling approaches exist to accommodate deviations from a strict molecular clock, assumptions about rate variation may not fully represent the underlying evolutionary processes. There is considerable variability in rabies virus (RABV) incubation periods, ranging from days to over a year, during which viral replication may be reduced. This prompts the question of whether modelling RABV on a per infection generation basis might be more appropriate. We investigate how variable incubation periods affect root-to-tip divergence under per-unit time and per-generation models of mutation. Additionally, we assess how well these models represent root-to-tip divergence in time-stamped RABV sequences. We find that at low substitution rates (<1 substitution per genome per generation) divergence patterns between these models are difficult to distinguish, while above this threshold differences become apparent across a range of sampling rates. Using a Tanzanian RABV dataset, we calculate the mean substitution rate to be 0.17 substitutions per genome per generation. At RABV's substitution rate, the per-generation substitution model is unlikely to represent rabies evolution substantially differently than the molecular clock model when examining contemporary outbreaks; over enough generations for any divergence to accumulate, extreme incubation periods average out. However, measuring substitution rates per-generation holds potential in applications such as inferring transmission trees and predicting lineage emergence.

Uncovering the endemic circulation of rabies in Cambodia

In epidemiology, endemicity characterizes sustained pathogen circulation in a geographical area, which involves a circulation that is not being maintained by external introductions. Because it could potentially shape the design of public health interventions, there is an interest in fully uncovering the endemic pattern of a disease. Here, we use a phylogeographic approach to investigate the endemic signature of rabies virus (RABV) circulation in Cambodia. Cambodia is located in one of the most affected regions by rabies in the world, but RABV circulation between and within Southeast Asian countries remains understudied. Our analyses are based on a new comprehensive data set of 199 RABV genomes collected between 2014 and 2017 as well as previously published Southeast Asian RABV sequences. We show that most Cambodian sequences belong to a distinct clade that has been circulating almost exclusively in Cambodia. Our results thus point towards rabies circulation in Cambodia that does not rely on external introductions. We further characterize within-Cambodia RABV circulation by estimating lineage dispersal metrics that appear to be similar to other settings, and by performing landscape phylogeographic analyses to investigate environmental factors impacting the dispersal dynamic of viral lineages. The latter analyses do not lead to the identification of environmental variables that would be associated with the heterogeneity of viral lineage dispersal velocities, which calls for a better understanding of local dog ecology and further investigations of the potential drivers of RABV spread in the region. Overall, our study illustrates how phylogeographic investigations can be performed to assess and characterize viral endemicity in a context of relatively limited data.

Clinical, epidemiological, and genotypic characteristics of rotavirus infection in hospitalized infants and young children in Yunnan Province

Rotaviruses are the most important pathogenic cause of non-bacterial diarrhea in infants and children. Approximately 60% of hospital admissions for acute diarrhea worldwide are caused by rotavirus infection. Rotavirus infection and hospitalization among children in China are a social burden, resulting in economic loss. The prevalence and geographical distribution of rotavirus genotypes is variable, partially due to population migration. Due to the unique geographical conditions and climate in Yunnan Province, several viruses with new genotypes have emerged, and multiple genotypes have become co-epidemic. In this study, rotavirus infection screening and genetic characterization of epidemic strains were performed in 149,492 infants and children admitted to hospitals in six representative prefectures in Yunnan Province between 2019 and 2021. The prevalence of rotavirus infection was 13.39% and was highest in January and lowest in September. G9P[8] was the main epidemic rotavirus genotype. Other epidemic genotypes included G2P[4], G8P[8], G9P[4], G2P[8], G3P[8], G4P[8], G3P[4], and G4P[6]. Phylogenetic analysis revealed that locally epidemic strains were influenced by importation of strains from neighboring provinces and other Asian countries. These findings provide a scientific basis for rotavirus prevention and control and lay a foundation for preliminary studies to establish a rotavirus surveillance network in Yunnan Province.

Periodic Characteristics of Hepatitis Virus Infections From 2013 to 2020 and Their Association With Meteorological Factors in Guangdong, China: Surveillance Study

BACKGROUND

In the past few decades, liver disease has gradually become one of the major causes of death and illness worldwide. Hepatitis is one of the most common liver diseases in China. There have been intermittent and epidemic outbreaks of hepatitis worldwide, with a tendency toward cyclical recurrences. This periodicity poses challenges to epidemic prevention and control.

OBJECTIVE

In this study, we aimed to investigate the relationship between the periodic characteristics of the hepatitis epidemic and local meteorological elements in Guangdong, China, which is a representative province with the largest population and gross domestic product in China.

METHODS

Time series data sets from January 2013 to December 2020 for 4 notifiable infectious diseases caused by hepatitis viruses (ie, hepatitis A, B, C, and E viruses) and monthly data of meteorological elements (ie, temperature, precipitation, and humidity) were used in this study. Power spectrum analysis was conducted on time series data, and correlation and regression analyses were performed to assess the relationship between the epidemics and meteorological elements.

RESULTS

The 4 hepatitis epidemics showed clear periodic phenomena in the 8-year data set in connection with meteorological elements. Based on the correlation analysis, temperature demonstrated the strongest correlation with hepatitis A, B, and C epidemics, while humidity was most significantly associated with the hepatitis E epidemic. Regression analysis revealed a positive and significant coefficient between temperature and hepatitis A, B, and C epidemics in Guangdong, while humidity had a strong and significant association with the hepatitis E epidemic, and its relationship with temperature was relatively weak.

CONCLUSIONS

These findings provide a better understanding of the mechanisms underlying different hepatitis epidemics and their connection to meteorological factors. This understanding can help guide local governments in predicting and preparing for future epidemics based on weather patterns and potentially aid in the development of effective prevention measures and policies.

Identifying high-risk areas for dog-mediated rabies using Bayesian spatial regression

Despite ongoing control efforts, rabies remains an endemic zoonotic disease in many countries. Determining high-risk areas and the space-time patterns of rabies spread, as it relates to epidemiologically important factors, can support policymakers and program managers alike to develop evidence-based targeted surveillance and control programs. In this One Health approach which selected Thailand as the example site, the location-based risk of contracting dog-mediated rabies by both human and animal populations was quantified using a Bayesian spatial regression model. Specifically, a conditional autoregressive (CAR) Bayesian zero-inflated Poisson (ZIP) regression was fitted to the reported human and animal rabies case counts of each district, from the 2012-2017 period. The human population was used as an offset. The epidemiologically important factors hypothesized as risk modifiers and therefore tested as predictors included: number of dog bites/attacks, the population of dogs and cats, number of Buddhist temples, garbage dumps, animal vaccination, post-exposure prophylaxis, poverty, and shared administrative borders. Disparate sources of data were used to improve the estimated associations and predictions. Model performance was assessed using cross-validation. Results suggested that accounting for the association between human and animal rabies with number of dog bites/attacks, number of owned and un-owned dogs; shared country borders, number of Buddhist temples, poverty levels, and accounting for spatial dependence between districts, may help to predict the risk districts for dog-mediated rabies in Thailand. The fitted values of the spatial regression were mapped to illustrate the risk of dog-mediated rabies. The cross-validation indicated an adequate performance of the spatial regression model (AUC = 0.81), suggesting that had this spatial regression approach been used to identify districts at risk in 2015, the cases reported in 2016/17 would have been predicted with model sensitivity and specificity of 0.71 and 0.80, respectively. While active surveillance is ideal, this approach of using multiple data sources to improve risk estimation may inform current rabies surveillance and control efforts including determining rabies-free zones, and the roll-out of human post-exposure prophylaxis and anti-rabies vaccines for animals in determining high-risk areas.

Oscillatory properties of class C notifiable infectious diseases in China from 2009 to 2021

Background

Epidemics of infectious diseases have a great negative impact on people's daily life. How it changes over time and what kind of laws it obeys are important questions that researchers are always interested in. Among the characteristics of infectious diseases, the phenomenon of recrudescence is undoubtedly of great concern. Understanding the mechanisms of the outbreak cycle of infectious diseases could be conducive for public health policies to the government.

Method

In this study, we collected time-series data for nine class C notifiable infectious diseases from 2009 to 2021 using public datasets from the National Health Commission of China. Oscillatory power of each infectious disease was captured using the method of the power spectrum analysis.

Results

We found that all the nine class C diseases have strong oscillations, which could be divided into three categories according to their oscillatory frequencies each year. Then, we calculated the oscillation power and the average number of infected cases of all nine diseases in the first 6 years (2009-2015) and the next 6 years (2015-2021) since the update of the surveillance system. The change of oscillation power is positively correlated to the change in the number of infected cases. Moreover, the diseases that break out in summer are more selective than those in winter.

Conclusion

Our results enable us to better understand the oscillation characteristics of class C infectious diseases and provide guidance and suggestions for the government's prevention and control policies.

Phylodynamics of HIV in the Mexico City Metropolitan Region

Evolutionary analyses of viral sequences can provide insights into transmission dynamics, which in turn can optimize prevention interventions. Here, we characterized the dynamics of HIV transmission within the Mexico City metropolitan area. HIV pol sequences from persons recently diagnosed at the largest HIV clinic in Mexico City (between 2016 and 2021) were annotated with demographic/geographic metadata. A multistep phylogenetic approach was applied to identify putative transmission clades. A data set of publicly available sequences was used to assess international introductions. Clades were analyzed with a discrete phylogeographic model to evaluate the timing and intensity of HIV introductions and transmission dynamics among municipalities in the region. A total of 6,802 sequences across 96 municipalities (5,192 from Mexico City and 1,610 from the neighboring State of Mexico) were included (93.6% cisgender men, 5.0% cisgender women, and 1.3% transgender women); 3,971 of these sequences formed 1,206 clusters, involving 78 municipalities, including 89 clusters of ≥10 sequences. Discrete phylogeographic analysis revealed (i) 1,032 viral introductions into the region, over one-half of which were from the United States, and (ii) 354 migration events between municipalities with high support (adjusted Bayes factor of ≥3). The most frequent viral migrations occurred between northern municipalities within Mexico City, i.e., Cuauhtémoc to Iztapalapa (5.2% of events), Iztapalapa to Gustavo A. Madero (5.4%), and Gustavo A. Madero to Cuauhtémoc (6.5%). Our analysis illustrates the complexity of HIV transmission within the Mexico City metropolitan area but also identifies a spatially active transmission area involving a few municipalities in the north of the city, where targeted interventions could have a more pronounced effect on the entire regional epidemic. IMPORTANCE Phylogeographic investigation of the Mexico City HIV epidemic illustrates the complexity of HIV transmission in the region. An active transmission area involving a few municipalities in the north of the city, with transmission links throughout the region, is identified and could be a location where targeted interventions could have a more pronounced effect on the entire regional epidemic, compared with those dispersed in other manners.

Changes in temporal properties for epidemics of notifiable infectious diseases in China during the COVID-19 epidemic: population-based surveillance study

BACKGROUND

The COVID-19 was first reported in 2019, and the Chinese government immediately carried out stringent and effective control measures in response to the epidemics.

OBJECTIVE

These non-pharmaceutical interventions may have impacted incidences of other infectious diseases as well. Potential explanations underlying this reduction, however, are not clear. Hence, in this study, we aimed to study the influence of the COVID-19 prevention policies on other infectious diseases (mainly class B infectious diseases) in China.

METHODS

The time-series datasets between 2017 and 2021 for 23 notifiable infectious diseases were extracted from public datasets from the National Health Commission of China. Several indices (peak and trough amplitude, infection selectivity, preferred time to outbreak, oscillatory strength) of each infectious disease were calculated before and after the COVID-19 outbreak.

RESULTS

We found that the prevention and control policies for COVID-19 had a strong significant reduction effect on outbreaks of other infectious diseases. A clear event-related trough (ERT) was observed after the outbreak of COVID-19 under the strict control policies, and its decreasing amplitude is related to the infection selectivity and preferred outbreak time of the disease before the COVID-19. We also calculated the oscillatory strength before and after the COVID-19 outbreak and found that it is significantly stronger before the COVID-19 outbreak, and does not correlate with the trough amplitude.

CONCLUSIONS

Our results directly demonstrate that prevention policies for the COVID-19 have immediate additional benefits for controlling most class B infectious diseases, and several factors (infection selectivity, preferred outbreak time) may have contributed to the reduction of outbreaks. This study may guide the implementation of non-pharmaceutical interventions to control a wider range of infectious diseases.

CLINICALTRIAL

On the Use of Phylogeographic Inference to Infer the Dispersal History of Rabies Virus: A Review Study

Rabies is a neglected zoonotic disease which is caused by negative strand RNA-viruses belonging to the genus Lyssavirus. Within this genus, rabies viruses circulate in a diverse set of mammalian reservoir hosts, is present worldwide, and is almost always fatal in non-vaccinated humans. Approximately 59,000 people are still estimated to die from rabies each year, leading to a global initiative to work towards the goal of zero human deaths from dog-mediated rabies by 2030, requiring scientific efforts from different research fields. The past decade has seen a much increased use of phylogeographic and phylodynamic analyses to study the evolution and spread of rabies virus. We here review published studies in these research areas, making a distinction between the geographic resolution associated with the available sequence data. We pay special attention to environmental factors that these studies found to be relevant to the spread of rabies virus. Importantly, we highlight a knowledge gap in terms of applying these methods when all required data were available but not fully exploited. We conclude with an overview of recent methodological developments that have yet to be applied in phylogeographic and phylodynamic analyses of rabies virus.

Enlightenment on oscillatory properties of 23 class B notifiable infectious diseases in the mainland of China from 2004 to 2020

A variety of infectious diseases occur in mainland China every year. Cyclic oscillation is a widespread attribute of most viral human infections. Understanding the outbreak cycle of infectious diseases can be conducive for public health management and disease surveillance. In this study, we collected time-series data for 23 class B notifiable infectious diseases from 2004 to 2020 using public datasets from the National Health Commission of China. Oscillatory properties were explored using power spectrum analysis. We found that the 23 class B diseases from the dataset have obvious oscillatory patterns (seasonal or sporadic), which could be divided into three categories according to their oscillatory power in different frequencies each year. These diseases were found to have different preferred outbreak months and infection selectivity. Diseases that break out in autumn and winter are more selective. Furthermore, we calculated the oscillation power and the average number of infected cases of all 23 diseases in the first eight years (2004 to 2012) and the next eight years (2012 to 2020) since the update of the surveillance system. A strong positive correlation was found between the change of oscillation power and the change in the number of infected cases, which was consistent with the simulation results using a conceptual hybrid model. The establishment of reliable and effective analytical methods contributes to a better understanding of infectious diseases’ oscillation cycle characteristics. Our research has certain guiding significance for the effective prevention and control of class B infectious diseases.

Animal Rabies Surveillance, China, 2004-2018

Rabies is a severe zoonotic disease in China, but the circulation and distribution of rabies virus (RABV) within animal reservoirs is not well understood. We report the results of 15 years of surveillance of the first Chinese Rabies Surveillance Plan in animal populations, in which animal brain tissues collected during 2004–2018 were tested for RABV and phylogenetic and spatial–temporal evolutionary analyses performed using obtained RABV sequences. The results have provided the most comprehensive dataset to date on the infected animal species, geographic distribution, transmission sources, and genetic diversity of RABVs in China. In particular, the transboundary transmission of emerging RABV subclades between China and neighboring countries was confirmed. The study highlights the importance of continuous animal rabies surveillance in monitoring the transmission dynamics, and provides updated information for improving current control and prevention strategies at the source.

Mathematical modelling and phylodynamics for the study of dog rabies dynamics and control: A scoping review

BACKGROUND

Rabies is a fatal yet vaccine-preventable disease. In the last two decades, domestic dog populations have been shown to constitute the predominant reservoir of rabies in developing countries, causing 99% of human rabies cases. Despite substantial control efforts, dog rabies is still widely endemic and is spreading across previously rabies-free areas. Developing a detailed understanding of dog rabies dynamics and the impact of vaccination is essential to optimize existing control strategies and developing new ones. In this scoping review, we aimed at disentangling the respective contributions of mathematical models and phylodynamic approaches to advancing the understanding of rabies dynamics and control in domestic dog populations. We also addressed the methodological limitations of both approaches and the remaining issues related to studying rabies spread and how this could be applied to rabies control.

METHODOLOGY/PRINCIPAL FINDINGS

We reviewed how mathematical modelling of disease dynamics and phylodynamics have been developed and used to characterize dog rabies dynamics and control. Through a detailed search of the PubMed, Web of Science, and Scopus databases, we identified a total of n = 59 relevant studies using mathematical models (n = 30), phylodynamic inference (n = 22) and interdisciplinary approaches (n = 7). We found that despite often relying on scarce rabies epidemiological data, mathematical models investigated multiple aspects of rabies dynamics and control. These models confirmed the overwhelming efficacy of massive dog vaccination campaigns in all settings and unraveled the role of dog population structure and frequent introductions in dog rabies maintenance. Phylodynamic approaches successfully disentangled the evolutionary and environmental determinants of rabies dispersal and consistently reported support for the role of reintroduction events and human-mediated transportation over long distances in the maintenance of rabies in endemic areas. Potential biases in data collection still need to be properly accounted for in most of these analyses. Finally, interdisciplinary studies were determined to provide the most comprehensive assessments through hypothesis generation and testing. They also represent new avenues, especially concerning the reconstruction of local transmission chains or clusters through data integration.

CONCLUSIONS/SIGNIFICANCE

Despite advances in rabies knowledge, substantial uncertainty remains regarding the mechanisms of local spread, the role of wildlife in dog rabies maintenance, and the impact of community behavior on the efficacy of control strategies including vaccination of dogs. Future integrative approaches that use phylodynamic analyses and mechanistic models within a single framework could take full advantage of not only viral sequences but also additional epidemiological information as well as dog ecology data to refine our understanding of rabies spread and control. This would represent a significant improvement on past studies and a promising opportunity for canine rabies research in the frame of the One Health concept that aims to achieve better public health outcomes through cross-sector collaboration.

Comparative Circulation Dynamics of the Five Main HIV Types in China

The HIV epidemic in China accounts for 3% of the global HIV incidence. We compared the patterns and determinants of interprovincial spread of the five most prevalent circulating types. HIV pol sequences sampled across China were used to identify relevant transmission networks of the five most relevant HIV-1 types (B and circulating recombinant forms [CRFs] CRF01_AE, CRF07_BC, CRF08_BC, and CRF55_01B) in China. From these, the dispersal history across provinces was inferred. A generalized linear model (GLM) was used to test the association between migration rates among provinces and several measures of human mobility. A total of 10,707 sequences were collected between 2004 and 2017 across 26 provinces, among which 1,962 are newly reported here. A mean of 18 (minimum and maximum, 1 and 54) independent transmission networks involving up to 17 provinces were identified. Discrete phylogeographic analysis largely recapitulates the documented spread of the HIV types, which in turn, mirrors within-China population migration flows to a large extent. In line with the different spatiotemporal spread dynamics, the identified drivers thereof were also heterogeneous but are consistent with a central role of human mobility. The comparative analysis of the dispersal dynamics of the five main HIV types circulating in China suggests a key role of large population centers and developed transportation infrastructures as hubs of HIV dispersal. This advocates for coordinated public health efforts in addition to local targeted interventions.IMPORTANCE While traditional epidemiological studies are of great interest in describing the dynamics of epidemics, they struggle to fully capture the geospatial dynamics and factors driving the dispersal of pathogens like HIV as they have difficulties capturing linkages between infections. To overcome this, we used a discrete phylogeographic approach coupled to a generalized linear model extension to characterize the dynamics and drivers of the across-province spread of the five main HIV types circulating in China. Our results indicate that large urbanized areas with dense populations and developed transportation infrastructures are facilitators of HIV dispersal throughout China and highlight the need to consider harmonized country-wide public policies to control local HIV epidemics.

Portable Rabies Virus Sequencing in Canine Rabies Endemic Countries Using the Oxford Nanopore MinION

As countries with endemic canine rabies progress towards elimination by 2030, it will become necessary to employ techniques to help plan, monitor, and confirm canine rabies elimination. Sequencing can provide critical information to inform control and vaccination strategies by identifying genetically distinct virus variants that may have different host reservoir species or geographic distributions. However, many rabies testing laboratories lack the resources or expertise for sequencing, especially in remote or rural areas where human rabies deaths are highest. We developed a low-cost, high throughput rabies virus sequencing method using the Oxford Nanopore MinION portable sequencer. A total of 259 sequences were generated from diverse rabies virus isolates in public health laboratories lacking rabies virus sequencing capacity in Guatemala, India, Kenya, and Vietnam. Phylogenetic analysis provided valuable insight into rabies virus diversity and distribution in these countries and identified a new rabies virus lineage in Kenya, the first published canine rabies virus sequence from Guatemala, evidence of rabies spread across an international border in Vietnam, and importation of a rabid dog into a state working to become rabies-free in India. Taken together, our evaluation highlights the MinION's potential for low-cost, high volume sequencing of pathogens in locations with limited resources.

The impact of anthropogenic and environmental factors on human rabies cases in China

Human rabies is a public health problem in Asia, especially in less-developed regions where the disease is under-reported because of a lack of epidemiological surveillance. To address this gap, we collected data on human rabies in Yunnan Province, China, between 2005 and 2016. Using statistical mapping techniques, we correlated the occurrence of human rabies to environmental (elevation, precipitation, normalized difference vegetation index [NDVI], temperature and distance to the nearest main rivers) and anthropogenic (human and dog population density, distance to the nearest main roads and gross domestic product [GDP]) factors. We used a performance score, the average area under the receiver operator characteristic curve (0.88), to validate our risk model. Using this model, we found that environmental factors were more strongly associated with human rabies occurrence than anthropogenic factors. Areas with elevation below 2000 metres, GDP per capita between $750 and $4500/year and NDVI below 0.07 were associated with greater risk of human rabies. Rabies control in China should specifically target these areas.

Using phylogeographic approaches to analyse the dispersal history, velocity and direction of viral lineages - Application to rabies virus spread in Iran

Recent years have seen the extensive use of phylogeographic approaches to unveil the dispersal history of virus epidemics. Spatially explicit reconstructions of viral spread represent valuable sources of lineage movement data that can be exploited to investigate the impact of underlying environmental layers on the dispersal of pathogens. Here, we performed phylogeographic inference and applied different post hoc approaches to analyse a new and comprehensive data set of viral genomes to elucidate the dispersal history and dynamics of rabies virus (RABV) in Iran, which have remained largely unknown. We first analysed the association between environmental factors and variations in dispersal velocity among lineages. Second, we present, test and apply a new approach to study the link between environmental conditions and the dispersal direction of lineages. The statistical performance (power of detection, false-positive rate) of this new method was assessed using simulations. We performed phylogeographic analyses of RABV genomes, allowing us to describe the large diversity of RABV in Iran and to confirm the cocirculation of several clades in the country. Overall, we estimate a relatively high lineage dispersal velocity, similar to previous estimates for dog rabies virus spread in northern Africa. Finally, we highlight a tendency for RABV lineages to spread in accessible areas associated with high human population density. Our analytical workflow illustrates how phylogeographic approaches can be used to investigate the impact of environmental factors on several aspects of viral dispersal dynamics.

Genomic characterization of human adenovirus type 4 strains isolated worldwide since 1953 identifies two separable phylogroups evolving at different rates from their most recent common ancestor

Species Human mastadenovirus E (HAdV-E) comprises several simian types and a single human type: HAdV-E4, a respiratory and ocular pathogen. RFLP analysis for the characterization of intratypic genetic variability has previously distinguished two HAdV-E4 clusters: prototype (p)-like and a-like. Our analysis of whole genome sequences confirmed two distinct lineages, which we refer to as phylogroups (PGs). PGs I and II comprise the p- and a-like genomes, respectively, and differ significantly in their G + C content (57.7% ± 0.013 vs 56.3% ± 0.015). Sequence differences distinguishing the two clades map to several regions of the genome including E3 and ITR. Bayesian analyses showed that the two phylogroups diverged approximately 602 years before the present. A relatively faster evolutionary rate was identified for PG II. Our data provide a rationale for the incorporation of phylogroup identity to HAdV-E4 strain designation to reflect the identified unique genetic characteristics that distinguish PGs I and II.

Phylodynamic and transmission pattern of rabies virus in China and its neighboring countries

Rabies is a fatal disease caused by infection with rabies virus (RABV), and human rabies is still a critical public-health concern in China. Although there have been some phylogenetic studies about RABV transmission patterns, with the accumulation of more rabies sequences in recent years, there is an urgent need to update and clarify the spatial and temporal patterns of RABV circulating in China on a national scale. In this study, we collected all available RABV nucleoprotein gene sequences from China and its neighboring countries and performed comparative analysis. We identified six significant subclades of RABV circulating in China and found that each of them has a specific geographical distribution, reflecting possible physical barriers to gene flow. The phylogeographic analysis revealed minimal viral movement among different geographical locations. An analysis using Bayesian coalescent methods indicated that the current RABV strains in China may come from a common ancestor about 400 years ago, and currently, China is amid the second event of increasing RABV population since the 1950s, but the population has decreased gradually. We did not detect any evidence of recombination in the sequence dataset, nor did we find any evidence for positive selection during the expansion of RABV. Overall, geographic location and neutral genetic drift may be the main factors in shaping the phylogeography of RABV transmission in China.

Transmission dynamics of re-emerging rabies in domestic dogs of rural China

Despite ongoing efforts to control transmission, rabies prevention remains a challenge in many developing countries, especially in rural areas of China where re-emerging rabies is under-reported due to a lack of sustained animal surveillance. By taking advantage of detailed genomic and epidemiological data for the re-emerging rabies outbreak in Yunnan Province, China, collected between 1999 and 2015, we reconstruct the demographic and dispersal history of domestic dog rabies virus (RABV) as well as the dynamics of dog-to-dog and dog-to-human transmission. Phylogeographic analyses reveal a lower diffusion coefficient than previously estimated for dog RABV dissemination in northern Africa. Furthermore, epidemiological analyses reveal transmission rates between dogs, as well as between dogs and humans, lower than estimates for Africa. Finally, we show that reconstructed epidemic history of RABV among dogs and the dynamics of rabid dogs are consistent with the recorded human rabies cases. This work illustrates the benefits of combining phylogeographic and epidemic modelling approaches for uncovering the spatiotemporal dynamics of zoonotic diseases, with both approaches providing estimates of key epidemiological parameters.

Identification of nucleotides in the 5'UTR and amino acids substitutions that are essential for the infectivity of 5'UTR-NS5A recombinant of hepatitis C virus genotype 1b (strain Con1)

Genotype 1b strain Con1 represents an important reference in the study of hepatitis C virus (HCV). Here, we aimed to develop an advanced infectious Con1 recombinant. We found that previously identified mutations A1226G/F1464L/A1672S/Q1773H permitted culture adaption of Con1 Core-NS5A (C-5A) recombinant containing 5'UTR and NS5B-3'UTR from JFH1 (genotype 2a), thus acquired additional mutations L725H/F886L/D2415G. C-5A containing all seven mutations (C-5A_7m) replicated efficiently in Huh7.5 and Huh7.5.1 cells and had an increased infectivity in SEC14L2-expressing Huh7.5.1 cells. Incorporation of Con1 NS5B was deleterious to C-5A_7m, however Con1 5'UTR was permissive but attenuated the virus. Nucleotides G1, A4, and G35 primarily accounted for the viral attenuation without affecting RNA translation. C-5A_7m was inhibited dose-dependently by simeprevir and daclatasvir, and substitutions at A4, A29, A34, and G35 conferred resistance to miR-122 antagonism. The novel Con1 5'UTR-NS5A recombinant, adaptive mutations, and critical nucleotides described here will facilitate future studies of HCV culture systems and virus-host interaction.

Codon usage bias in the N gene of rabies virus

Since its emergence, rabies virus (RABV) has been a major worldwide concern especially in developing countries. The nucleoprotein (N) of RABV is highly conserved and key for genetic typing, thus a better understanding of the N gene evolutionary trajectory can assist the development of control measures. We found that the N gene of RABV has a low codon usage bias with a mean effective number of codons (ENC) value of 56.33 influenced by both mutation pressure and natural selection. However, neutrality analysis indicated that natural selection dominates over mutation pressure. Additionally, we found that dinucleotide bias partly contributed to RABV codon usage bias. On the other hand, based on the clades of phylogenetic tree, we found that the evolutionary rate of the Africa 2 clade was the highest with a mean value of 3.75×10^-3 substitutions per site per year. Above all, our results regarding N gene of RABV codon usage will serve future RABV evolution research.