Genetic characterization and molecular epidemiological analysis of enterovirus C99 in China

Enterovirus C99 (EV-C99) is a newly identified EV serotype within the species Enterovirus C. Few studies on EV-C99 have been conducted globally. More information and research on EV-C99 are needed to assess its genetic characteristics, phylogenetic relationships, and associations with enteroviral diseases. Here, the phylogenetic characteristics of 11 Chinese EV-C99 strains have been reported. The full-length genomic sequences of these 11 strains show 79.4-80.5% nucleotide identity and 91.7-94.3% amino acid (aa) identity with the prototype EV-C99. A maximum likelihood phylogenetic tree constructed based on the entire VP1 coding region identified 13 genotypes (A-M), revealing a high degree of variation among the EV-C99 strains. Phylogeographic analysis showed that the Xinjiang Uygur Autonomous Region is an important source of EV-C99 epidemics in various regions of China. Recombination analysis revealed inter-serotype recombination events of 16 Chinese EV-C99 strains in 5' untranslated regions and 3D regions, resulting in the formation of a single recombination form. Additionally, the Chinese strain of genotype J showed rich aa diversity in the P1 region, indicating that the genotype J of EV-C99 is still going through variable dynamic changes. This study contributes to the global understanding of the EV-C99 genome sequence and holds substantial implications for the surveillance of EV-C99.

Genetic structure and diversity of green turtle (Chelonia mydas) in the Gulf of Thailand

Background and Aim

The International Union for the Conservation of Nature and Natural Resources lists the green turtle as endangered. Green turtle nesting behavior in the Gulf of Thailand has decreased to <50% of the 1995 level. The population structure of green turtles in the Gulf of Thailand has not yet been studied. This study aimed to characterize the genetic diversity of green turtles in the Gulf of Thailand based on comparisons of mitochondrial DNA (mtDNA) control region with sequences of Indo-Pacific management units (MUs) and rookeries, to investigate population structures, and to explore phylogeographic relationships.

Materials and Methods

Blood samples (1 mL each) from 91 stranded green turtles were collected from four parts of the Gulf of Thailand (eastern, upper, central, and lower). The control mtDNA region was amplified by polymerase chain reaction using LCM15382 and H950 primer. The obtained 384-bp or 770-bp sequences were analyzed for haplotype, clade, and haplotype and nucleotide diversities and were used to construct a phylogenetic tree and haplotype network diagram, respectively. In addition, we analyzed genetic differentiation within and among populations of green turtles in the Gulf of Thailand and between green turtles in the Gulf of Thailand and other Indo-Pacific MUs and rookeries.

Results

In total, 12 (based on 384 bp) or 13 (based on 770 bp) haplotypes and two clades (clades VII and VIII) were identified, with nine or 10 haplotypes belonging to clade VIII and three haplotypes belonging to clade VII. Of the new haplotypes, four or five were identified and classified as clade VII (two haplotypes, for both fragment lengths) and clade VIII (two or three haplotypes, for 384 bp or 770 bp fragments, respectively). The overall haplotype and nucleotide diversity of green turtles in the Gulf of Thailand were high (0.755 ± 0.039 and 0.01146 ± 0.00248, respectively). Based on the analysis of molecular variance, green turtles in the Gulf of Thailand could be divided into two subpopulations (UC-Eastern Gulf of Thailand [UC-EGT] and lower Gulf of Thailand [LGT]). Comparisons with other MUs and rookeries in the Indo-Pacific showed that UC-EGT was not genetically different from the Peninsular Malaysia and Eastern Taiwan (Lanyu) MUs and the Terrangganu and Mersing rookeries, and LGT were not genetically different from Peninsular Malaysia, Sipadan, Brunei Bay, Eastern Taiwan (Lanyu), Scott Reef and Browse Island, and Gulf of Carpentaria MUs and the Perak, Perhentain Island, Redang, Pahang, and Vietnam rookeries.

Conclusion

To the best of our knowledge, this is the first report to identify the haplotypes and clades of green turtles in the Gulf of Thailand and to show that the populations in the Gulf of Thailand not only present high genetic diversity but also have haplotypic endemism. Longer mtDNA fragments (770 bp) increased the resolution of the stock structure. Clade VII is a unique clade not only for Japan but also for Thailand and Malaysia, and CmP82 is a unique haplotype for both the Gulf of Thailand and Malaysia. Conservation and management of these populations are important to preserve the genetic diversity, biological diversity, and evolutionary potential of green turtles in the Gulf of Thailand.

Molecular Structure and Variation Characteristics of the Plastomes from Six Malus baccata (L.) Borkh. Individuals and Comparative Genomic Analysis with Other Malus Species

Malus baccata (L.) Borkh. is an important wild species of Malus. Its rich variation types and population history are not well understood. Chloroplast genome mining plays an active role in germplasm identification and genetic evolution. In this study, by assembly and annotation, six complete cp genome sequences, ranging in size from 160,083 to 160,295 bp, were obtained. The GC content of stable IR regions (42.7%) was significantly higher than that of full length (36.5%) and SC regions (LSC-34.2%, SSC-30.4%). Compared with other Malus species, it was found that there were more sites of polymorphisms and hotspots of variation in LSC and SSC regions, with high variation sites including trnR/UCU-atpA, trnT/UGU-trnL/UAA, ndhF-rpl32 and ccsA-ndhD. The intraspecific and interspecific collinearity was good, and no structural rearrangement was observed. A large number of repeating elements and different boundary expansions may be involved in shaping the cp genome size. Up to 77 or 78 coding genes were annotated in the cp genomes of M. baccata, and high frequency codons such as UUA (Leu), GCU (Ala) and AGA (Arg) were identified by relative synonymous codon usage analysis. Phylogeographic analysis showed that 12 individuals of M. baccata clustered into three different groups with complex structure, whereas variant xiaojinensis (M.H. Cheng & N.G. Jiang) was not closely related to M. baccata evolutionarily. The phylogenetic analysis suggested that two main clades of different M. baccata in the genus Malus were formed and that I and II diverged about 9.7 MYA. In conclusion, through cp genome assembly and comparison, the interspecific relationships and molecular variations of M. baccata were further elucidated, and the results of this study provide valuable information for the phylogenetic evolution and germplasm conservation of M. baccata and Malus.

Genomic dissection of the microevolution of Australian epidemic Bordetella pertussis

Whooping cough (pertussis) is a highly contagious respiratory disease caused by the bacterium Bordetella pertussis. Despite high vaccine coverage, pertussis has re-emerged in many countries including Australia and caused two large epidemics in Australia since 2007. Here, we undertook a genomic and phylogeographic study of 385 Australian B. pertussis isolates collected from 2008 to 2017. The Australian B. pertussis population was found to be composed of mostly ptxP3 strains carrying different fim3 alleles, with ptxP3-fim3A genotype expanding far more than ptxP3-fim3B. Within the former, there were six co-circulating epidemic lineages (EL1 to EL6). The multiple ELs emerged, expanded, and then declined at different time points over the two epidemics. In population genetics terms, both hard and soft selective sweeps through vaccine selection pressures have determined the population dynamics of Australian B. pertussis. Relative risk estimation suggests that once a new B. pertussis lineage emerged, it was more likely to spread locally within the first 1.5 years. However, after 1.5 years, any new lineage was likely to expand to a wider region. Phylogenetic analysis revealed the expansion of ptxP3 strains was also associated with replacement of the type III secretion system allele bscI1 with bscI3. bscI3 is associated with decreased T3SS secretion and may allow B. pertussis to reduce immune recognition. This study advanced our understanding of the epidemic population structure and spatial and temporal dynamics of B. pertussis in a highly immunized population.

Dating the origin and dispersal of global hepatitis B virus genotype C in humans

Hepatitis B virus genotype C (HBV/C) is one of the most prevalent HBV strains worldwide, especially in the Western Pacific and the South-East Asia. However, the origin and evolutionary timescale of HBV/C remains largely unresolved. We analyzed the evolutionary rate and molecular clock phylogeny of 101 full-genome HBV/C sequences sampled globally using a Bayesian Markov Chain Monte Carlo (MCMC) approach. We inferred the spatiotemporal dynamics of the HBV/C worldwide by the Bayesian Stochastic Search Variable Selection (BSSVS). We found that the estimated mean evolution rate of the HBV/C genotype full-genome was 4.32 × 10^-5 subs/site/year (95% highest posterior density 3.02 × 10^-6 - 8.97 × 10^-5). Phylogeographic reconstruction was able to identify a single location for the origin of the global HBV/C in Australia around A.D. 715. The subgenotype C4 diverged earliest and mainly circulated in Australia, C1 mainly in Southeast Asia, C2 mainly in East Asia and C3 in Remote Oceania. The effective number of HBV infection presented a rapid exponential increase between the 1760s and 1860s followed by a maintained high level until now. Our study, for the first time, provides an estimated timescale for the HBV/C epidemic, and brings new insight to the dispersal of HBV/C in humans globally. Based on the continuous presence of a highly effective viral population, this study provides further evidence of the challenge from a population-based molecular level to eliminate HBV by 2030, and calls for a concerted effort from policy makers, health providers, and society in the globalized world.

Genomic and Epidemiological Characteristics Provide Insights into the Phylogeographic Spread of Goose Astrovirus in China

Goose astrovirus (GAstV) is an emerging pathogen with a wide distribution in China that causes visceral gout and leads to significant economic losses in the goose industry. Here, 10 GAstV strains were isolated from different farms in southeast China. We performed an integrated analysis of the full-genome sequences of these new strains alongside comprehensive epidemiological surveillance information from the database. Interestingly, the results showed two distinct genotypes of GAstV, which were evolutionarily distant from each other. Group I GAstVs were closely related to DAstV IV, and group II strains were classified with duck astrovirus (DAstV) II and turkey astrovirus (TAstV) II. Further investigation showed that among the GAstV I strains, ZJC14 and AHDY differed from FLX. Comparative analysis of 58 available genomes clustered the GAstV II strains into two subgroups. We identified two major mutation sites, 456 (E/D) and 540 (L/Q), in the capsid protein, which were related to distinct subgroups according to evolution. GAstV II subgroup 1a strains are the predominant strains in the current prevalent epidemiology. Phylogeographic analysis based on 90 reported cases from 13 provinces revealed the complexity and severity of GAstV epidemics in China, within which Henan, Anhui and Jiangsu provinces have suffered great impacts. According to these phylogeographic investigations, following the initial introduction of GAstV from Hunan Province, the dispersal of GAstV with different subgenotypes on a nationwide scale may be explained by the live gosling trade. Our findings have important implications for the evolution and dispersal of GAstV and will contribute to understanding the potential risk of GAstV. This article is protected by copyright. All rights reserved.

Phylogeographic Assessment Reveals Geographic Sources of HIV-1 Dissemination Among Men Who Have Sex With Men in Kenya

HIV-1 transmission dynamics involving men who have sex with men (MSM) in Africa are not well understood. We investigated the rates of HIV-1 transmission between MSM across three regions in Kenya: Coast, Nairobi, and Nyanza. We analyzed 372 HIV-1 partial pol sequences sampled during 2006-2019 from MSM in Coast (N = 178, 47.9%), Nairobi (N = 137, 36.8%), and Nyanza (N = 57, 15.3%) provinces in Kenya. Maximum-likelihood (ML) phylogenetics and Bayesian inference were used to determine HIV-1 clusters, evolutionary dynamics, and virus migration rates between geographic regions. HIV-1 sub-subtype A1 (72.0%) was most common followed by subtype D (11.0%), unique recombinant forms (8.9%), subtype C (5.9%), CRF 21A2D (0.8%), subtype G (0.8%), CRF 16A2D (0.3%), and subtype B (0.3%). Forty-six clusters (size range 2-20 sequences) were found-half (50.0%) of which had evidence of extensive HIV-1 mixing among different provinces. Data revealed an exponential increase in infections among MSM during the early-to-mid 2000s and stable or decreasing transmission dynamics in recent years (2017-2019). Phylogeographic inference showed significant (Bayes factor, BF > 3) HIV-1 dissemination from Coast to Nairobi and Nyanza provinces, and from Nairobi to Nyanza province. Strengthening HIV-1 prevention programs to MSM in geographic locations with higher HIV-1 prevalence among MSM (such as Coast and Nairobi) may reduce HIV-1 incidence among MSM in Kenya.

The Viral Founder Effect and Economic-Driven Human Mobility Shaped the Distinct Epidemic Pattern of HIV-1 CRF01_AE in Northeast China

Background: In China, two distinct lineages shaped the epidemic of HIV-1 CRF01_AE among men who have sex with men (MSM), of which the uneven distributions were observed geographically. One lineage spread across China, while another dominated in Northeast China. Understanding the drivers of viral diffusion would provide guidelines for identifying the source and hotspots of HIV transmission among MSM to target interventions in China.

Methods: We collected the pol sequences between 2002–2017 to reconstruct the spatiotemporal history of CRF01_AE lineages in Shenyang, one economic center of Northeast China, using the Bayesian phylogeographic and phylodynamic approaches. Importantly, for the datasets with the high sample density, we did the down-sampling to avoid the sampling bias.

Results: Two lineages accounted for 97%, including 426 and 1516 sequences, and homosexuals and bisexuals were above 80%. One lineage appeared earlier 7 years than another (1993 vs. 2002) among homosexuals and bisexuals, whereas among heterosexuals, both lineages were observed firstly in 2002. 96% viral migrations within one lineage were from homosexuals toward bisexuals (49%) and male-heterosexuals (46%). Within another, except for homosexuals (72%), bisexuals (23%) served as the top second source, and female-heterosexuals (11%) were the third recipients following bisexuals (44%) and male-heterosexuals (39%). Although the basic reproduction number (R₀) of two lineages were similar and both of the effective production number (R_e) fell below 1 at the most recent sampling time, the starts of the R_e declining varied.

Conclusions: Our findings revealed that throughout the viral national spread chain, Shenyang is the source for the initial expanding of one lineage, where is only a sink of another, proving that the viral founder effect and regional human mobility contributed to the uneven distribution of two lineages, and emphasizing the important roles of the area where the virus originated and economy-driven migrants in HIV transmission.

Divergence in the Aquilegia ecalcarata complex is correlated with geography and climate oscillations: Evidence from plastid genome data

Quaternary climate oscillations and geographical heterogeneity play important roles in determining species and genetic diversity distribution patterns, but how these factors affect the migration and differentiation of East Asian plants species at the population level remains poorly understood. The Aquilegia ecalcarata complex, a group that originated in the Late Tertiary and is widely distributed throughout East Asia, displays high genetic variation that is suitable for studying elaborate phylogeographic patterns and demographic history related to the impact of Quaternary climate and geography. We used plastid genome data from 322 individuals in 60 populations of the A. ecalcarata complex to thoroughly explore the impact of Quaternary climate oscillations and geography on the phylogeographic patterns and demographic history of the A. ecalcarata complex through a series of phylogenetic, divergence time estimation, and demographic history analyses. The dry, cold climate and frequent climate oscillations that occurred during the early Pleistocene and the Mid-Pleistocene transition led to the differentiation of the A. ecalcarata complex, which was isolated in various areas. Geographically, the A. ecalcarata complex can be divided into Eastern and Western Clades and five subclades, which conform to the divergence of the East Asian flora. Our results clearly show the impact of Quaternary climate and geography on evolutionary history at the population level. These findings promote the understanding of the relationship between plant genetic differentiation and climate and geographical factors of East Asia at the population level.

The prevalence of CRF55_01B among HIV-1 strain and its connection with traffic development in China

CRF55_01B is a relatively “young” HIV strain. At present, we do not know much about its transmission characteristics in China. So, to describe the transmission characteristics of CRF55_01B strain among provinces and HIV infected people, and to analyze the reasons for its rapid epidemic in China, a total of 1237 subjects infected with CRF55_01B from 31 provinces spanning a period of 12 years from 2007 to 2018 were enrolled in this study. By constructing a molecular network and Bayesian correlation analysis, we found that CRF55_01B increased exponentially from 2005 to 2009 after its origin in Shenzhen, and increased rapidly after 2010. CRF55_01B began to spread to other provinces in 2007. After 2010, the strain showed a trend of rapid spread and epidemic from Guangdong-Shenzhen to other provinces in China. Guangdong, Shenzhen, Hunan, Beijing, Guangxi, Hubei, Jiangxi, Guizhou, Hebei, Anhui, Shanghai, Shandong, Henan, and Yunnan were the key provinces of CRF55_01B transmission. CRF55_01B, although originating from men who sex with men (MSM), was transmitted among heterosexuals in 2010. Males in heterosexuals played a crucial role in the transmission and diffusion of this strain. We also revealed that CRF55_01B might spread rapidly along with the rapid development of the Beijing-Guangzhou and Beijing-Kowloon railways. This study suggests that if we detect the spread of MSMs in time through molecular monitoring in the early stage of the epidemic, it can help us control the epidemic early and prevent its spread, which is of great significance to China's national prevention and control of HIV-1.

Genomic Epidemiology, Evolution, and Transmission Dynamics of Porcine Deltacoronavirus

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has shown once again that coronavirus (CoV) in animals are potential sources for epidemics in humans. Porcine deltacoronavirus (PDCoV) is an emerging enteropathogen of swine with a worldwide distribution. Here, we implemented and described an approach to analyze the epidemiology of PDCoV following its emergence in the pig population. We performed an integrated analysis of full genome sequence data from 21 newly sequenced viruses, along with comprehensive epidemiological surveillance data collected globally over the last 15 years. We found four distinct phylogenetic lineages of PDCoV, which differ in their geographic circulation patterns. Interestingly, we identified more frequent intra- and interlineage recombination and higher virus genetic diversity in the Chinese lineages compared with the USA lineage where pigs are raised in different farming systems and ecological environments. Most recombination breakpoints are located in the ORF1ab gene rather than in genes encoding structural proteins. We also identified five amino acids under positive selection in the spike protein suggesting a role for adaptive evolution. According to structural mapping, three positively selected sites are located in the N-terminal domain of the S1 subunit, which is the most likely involved in binding to a carbohydrate receptor, whereas the other two are located in or near the fusion peptide of the S2 subunit and thus might affect membrane fusion. Finally, our phylogeographic investigations highlighted notable South-North transmission as well as frequent long-distance dispersal events in China that could implicate human-mediated transmission. Our findings provide new insights into the evolution and dispersal of PDCoV that contribute to our understanding of the critical factors involved in CoVs emergence.

Phylogenomic Based Comparative Studies on Indian and American Commensal Staphylococcus epidermidis Isolates

Staphylococcus epidermidis is a prominent commensal member of human skin microbiome and an emerging nosocomial pathogen, making it a good model organism to provide genomic insights, correlating its transition between commensalism and pathogenicity. While there are numerous studies to understand differences in commensal and pathogenic isolates, systematic efforts to understand variation and evolutionary pattern in multiple strains isolated from healthy individuals are lacking. In the present study, using whole genome sequencing and analysis, we report presence of diverse lineages of S. epidermidis isolates in healthy individuals from two geographically diverse locations of India and North America. Further, there is distinct pattern in the distribution of candidate gene(s) for pathogenicity and commensalism. The pattern is not only reflected in lineages but is also based on geographic origin of the isolates. This is evident by the fact that North American isolates under this study are more genomically dynamic and harbor pathogenicity markers in higher frequency. On the other hand, isolates of Indian origin are less genomically dynamic, harbor less pathogenicity marker genes and possess two unique antimicrobial peptide gene clusters. This study provides a basis to understand the nature of selection pressure in a key human skin commensal bacterium with implications in its management as an opportunistic pathogen.

Evolution and Transmission of Respiratory Syncytial Group A (RSV-A) Viruses in Guangdong, China 2008-2015

Respiratory syncytial viruses (RSVs) including subgroups A (RSV-A) and B (RSV-B) are an important cause of acute respiratory tract infections worldwide. RSV-A include major epidemic strains. Fundamental questions concerning the evolution, persistence and transmission of RSV-A are critical for disease control and prevention, yet remain unanswered. In this study, we generated 64 complete G gene sequences of RSV-A strains collected between 2008 and 2015 in Guangdong, China. Phylogenetic analysis was undertaken by incorporating 572 publicly available RSV-A sequences. Current data indicate that genotypes GA1, GA4, and GA5 are endemic with limited epidemic activity. In contrast, the GA2 genotype which likely originated in 1980 has spread rapidly and caused epidemics worldwide. By analyzing GA2 genotype sequences across epidemic seasons within Guangdong, we find that RSV-A epidemics in Guangdong are caused by a combination of virus importation and local persistence, although the magnitude of the latter is likely overestimated due to infrequent sampling in other regions. Our results provide new insights into RSV-A evolution and transmission at global and local scales and highlights the rapid and wide spread of genotype GA2 compared to other genotypes. In order to control RSV transmission and outbreak, both local persistence and external introduction should be taken into account when designing optimal strategies.

The Evolution and Transmission of Epidemic GII.17 Noroviruses

BACKGROUND

In recent decades, the GII.4 norovirus genotype has predominated in epidemics worldwide and been associated with an increased rate of evolutionary change. In 2014, a novel GII.17 variant emerged and persisted, causing large outbreaks of gastroenteritis in China and sporadic infections globally. The origin, evolution, and transmission history of this new variant are largely unknown.

METHODS

We generated 103 full capsid and 8 whole-genome sequences of GII.17 strains collected between August 2013 and November 2015 in Guangdong, China. Phylogenetic analyses were performed by integrating our data with those for all publically available GII.17 sequences.

RESULTS

The novel emergent lineage GII.17_Kawasaki_2014 most likely originated from Africa around 2001 and evolved at a rate of 5.6 × 10(-3) substitutions/site/year. Within this lineage, a new variant containing several important amino acid changes emerged around August 2013 and caused extensive epidemics in 2014-2015. The phylodynamic and epidemic history of the GII.17_Kawasaki lineage shows similarities with the pattern observed for GII.4 norovirus evolution. Virus movements from Hong Kong to neighboring coastal cities were frequently observed.

CONCLUSIONS

Our results provide new insights into GII.17 norovirus evolution and transmission and highlight the potential for a rare norovirus genotype to rapidly replace existing strains and cause local epidemics.

Spatial Temporal Dynamics and Molecular Evolution of Re-Emerging Rabies Virus in Taiwan

Taiwan has been recognized by the World Organization for Animal Health as rabies-free since 1961. Surprisingly, rabies virus (RABV) was identified in a dead Formosan ferret badger in July 2013. Later, more infected ferret badgers were reported from different geographic regions of Taiwan. In order to know its evolutionary history and spatial temporal dynamics of this virus, phylogeny was reconstructed by maximum likelihood and Bayesian methods based on the full-length of glycoprotein (G), matrix protein (M), and nucleoprotein (N) genes. The evolutionary rates and phylogeographic were determined using Beast and SPREAD software. Phylogenetic trees showed a monophyletic group containing all of RABV isolates from Taiwan and it further separated into three sub-groups. The estimated nucleotide substitution rates of G, M, and N genes were between 2.49 × 10⁻⁴–4.75 × 10⁻⁴ substitutions/site/year, and the mean ratio of dN/dS was significantly low. The time of the most recent common ancestor was estimated around 75, 89, and 170 years, respectively. Phylogeographic analysis suggested the origin of the epidemic could be in Eastern Taiwan, then the Formosan ferret badger moved across the Central Range of Taiwan to western regions and separated into two branches. In this study, we illustrated the evolution history and phylogeographic of RABV in Formosan ferret badgers.

Phylogenetic and phylogeographic analyses of porcine circovirus type 2 among pig farms in Vietnam

This study demonstrated the prevalence of Porcine circovirus type 2 (PCV2) among pig farms in Vietnam. Analyses of the genome, capsid protein and phylogeny classified all 30 Vietnamese PCV2 strains as the PCV2b genotype, belonging to the clusters of 1A, 1B, 1C and recombinant forms. Each viral genome was 1767 nucleotides long and shared 96.0-100% nucleotide sequence identity. The amino acid substitutions in the capsid protein of the Vietnamese PCV2 strains were in immunodominant regions, and the majority of strains (24/30) contained a lysine extension at the C-terminus. Bayesian phylogeographic analysis revealed epidemic links of the PCV2 recombinant cluster within and among countries, which supports a circulating recombinant form of PCV2. Further analysis by the Jameson-Wolf antigenic index indicated antigenic alterations at important sites in the capsid protein (sites 131-133) among the recombinant cluster and the other clusters of PCV2b.