Mean Substitution Rate Research Articles

Phylogeography and Evolutionary Dynamics of Tobacco Curly Shoot Virus

Tobacco curly shoot virus (TbCSV), a begomovirus, causes significant economic losses in tobacco and tomato crops across East, Southeast, and South Asia. Despite its agricultural importance, the evolutionary dynamics and emergence process of TbCSV remain poorly understood. This study analyzed the phylodynamics of TbCSV by examining its nucleotide sequences of the coat protein (CP) gene collected between 2000 and 2022. Using various combinations of priors, Bayes factor comparisons identified heterochronous datasets (3 × 100 million chains) generated from a strict molecular clock and Bayesian skyline tree priors as the most robust. The mean substitution rate of the CP gene was estimated at 6.50 × 10−4 substitutions/site/year (95% credibility interval: 4.74 × 10−4–8.50 × 10−4). TbCSV was inferred to have diverged around 1920 CE (95% credibility interval: 1887–1952), with its most probable origin in South Asia. These findings provide valuable insights for the phylogeography and evolutionary dynamics of TbCSV, and contribute to a broader understanding of begomovirus epidemiology.

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.

Genetic diversity of Cotton leafroll dwarf virus from the Southwestern United States, and its implications for the multi-introduction event hypothesis and future evolution.

Cotton leafroll dwarf virus (CLRDV) is a viral agent recently identified in the United States (US) in 2017 in Alabama. Since its identification, CLRDV has spread to every cotton-growing state east of New Mexico. Oklahoma, Kansas, and Texas comprise the westernmost border of reported CLRDV incidence, making monitoring of these states vital for proper control. Additionally, as the virus evolves, mutations that alter symptomology, such as mutations in the F-box-like motif in ORF0/P0, may occur and need to be monitored thoroughly during the growing seasons. Using High-throughput sequencing (HTS) and PCR-derived Sanger sequencing, four CLRDV genomes and 21 P0 gene isolates were sequenced from Oklahoma, Kansas, and Texas from 2019 to 2021 to determine the genetic diversity among CLRDV isolates. Phylogenetic analyses of the complete genomes revealed seven clades while ORF0 gene analyses resulted in large polytomic clusters. BEAST analyses of the 114 total P0 sequences from GenBank, downloaded before 2024, revealed a lower mean substitution rate than previously reported as well as an earlier root year (1914). In addition, using all available CLRDV genome sequences, 11 likely recombination events were determined. Examination of the P0 amino acid sequences revealed 13 mutations unique to the isolates collected in this study. Based on the phylogenetic and amino acid analyses, the CLRDV isolates from Texas (TX clade) may represent evidence for the multi-introduction event hypothesis into the US. Additionally, based on our analyses in this study, we propose the Asian CLRDV isolates should be constituted as a potentially separate strain of CLRDV.

Genome-wide analysis reveals the contributors to fast molecular evolution of the Chinese hook snout carp (Opsariichthys bidens)

Variations in molecular evolutionary rate have been widely investigated among lineages and genes. However, it remains an open question whether fast rate of molecular evolution is driven by natural selection or random drift, and how the fast rate is linked to metabolic rate. Additionally, previous studies on fast molecular evolution have been largely restricted to concatenated matrix of genes or a few specifically selected genes, but less is known for individual genes at the genome-wide level. Here we addressed these questions using more than 5000 single-copy orthologous (SCO) genes through comparative genomic and phylogenetic analyses among fishes, with a special focus on a newly-sequenced clupeocephalan fish the Chinese hook snout carp Opsariichthys bidens. We showed O. bidens displays significantly higher mean substitution rate and more fast-evolving SCO genes (2172 genes) than most fishes studied here. The rapidly evolving genes are enriched in highly conserved and very basic functions such as translation and ribosome that are critical for biological fitness. We further revealed that ∼25 % of these fast-evolving genes exhibit a constant increase of substitution rate from the common ancestor down to the present, suggesting a neglected but important contribution from ancestral states. Model fitting showed that ∼85 % of fast-evolving genes exclusive to O. bidens and related species follow the adaptive evolutionary model rather than random-drift model, and 7.6 % of fast-evolving genes identified in O. bidens have experienced positive selection, both indicating the reflection of adaptive selection. Finally, metabolic rate was observed to be linked with substitution rate in a gene-specific manner. Overall, our findings reveal fast molecular evolution of SCO genes at genome-wide level in O. bidens, and uncover the evolutionary and ecological contributors to it.

Molecular Epidemiology and Evolution of Coxsackievirus A14.

As the proportion of non-enterovirus 71 and non-coxsackievirus A16 which proportion of composition in the hand, foot, and mouth pathogenic spectrum gradually increases worldwide, the attention paid to other enteroviruses has increased. As a member of the species enterovirus A, coxsackievirus A14 (CVA14) has been epidemic around the world until now since it has been isolated. However, studies on CVA14 are poor and the effective population size, evolutionary dynamics, and recombination patterns of CVA14 are not well understood. In this study, 15 CVA14 strains were isolated from HFMD patients in mainland China from 2009 to 2019, and the complete sequences of CVA14 in GenBank as research objects were analyzed. CVA14 was divided into seven genotypes A-G based on an average nucleotide difference of the full-length VP1 coding region of more than 15%. Compared with the CVA14 prototype strain, the 15 CVA14 strains showed 84.0-84.7% nucleotide identity in the complete genome and 96.9-97.6% amino acid identity in the encoding region. Phylodynamic analysis based on 15 CVA14 strains and 22 full-length VP1 sequences in GenBank showed a mean substitution rate of 5.35 × 10-3 substitutions/site/year (95% HPD: 4.03-6.89 × 10-3) and the most recent common ancestor (tMRCA) of CVA14 dates back to 1942 (95% HPD: 1930-1950). The Bayesian skyline showed that the effective population size had experienced a decrease-increase-decrease fluctuation since 2004. The phylogeographic analysis indicated two and three possible migration paths in the world and mainland China, respectively. Four recombination patterns with others of species enterovirus A were observed in 15 CVA14 strains, among which coxsackievirus A2 (CVA2), coxsackievirus A4 (CVA4), coxsackievirus A6 (CVA6), coxsackievirus A8 (CVA8), and coxsackievirus A12 (CVA12) may act as recombinant donors in multiple regions. This study has filled the gap in the molecular epidemiological characteristics of CVA14, enriched the global CVA14 sequence database, and laid the epidemiological foundation for the future study of CVA14 worldwide.

Phylogenetic characteristics and recombination analysis of echovirus 5 associated with severe acute respiratory infection in China.

Echovirus 5 (E5) belongs to the enterovirus B (EV-B) species, which is one of the pathogens causing respiratory diseases. Only two whole genome sequences (WGSs) of E5 are available in GenBank, and few studies on E5 have been reported worldwide. In this study, human rhabdomyosarcoma (RD) cells were used to isolate E5 from respiratory samples associated with severe acute respiratory infection (SARI). The nucleotide sequence of the VP1 region and the WGSs of two E5 strains were sequenced using Sanger sequencing and high-throughput sequencing technology, respectively. Together with 31 E5 VP1 sequences and 2 WGSs downloaded from GenBank, a total of 33 complete VP1 sequences and 4 WGSs were used for genetic characterization and evolutionary analysis. Phylogenetic analysis revealed five genotypes (A-E) based on complete VP1 sequences. Evolutionary analysis showed that the time to the most recent common ancestor of E5 was approximately 1952 (95% highest posterior density [HPD] range, 1948-1953]), with a mean substitution rate of 7.74 × 10-3 subs/site/year (95% HPD range, 6.44-9.11 × 10-3). In addition, the WGS analyses showed that the nucleotide and amino acid similarities of the coding regions of the two strains obtained in this study with the prototype of the E5 strain Noyce were 80.2% and 97.2%, respectively. Recombination analysis showed that the two E5 strains may have intraspecific recombination with E6, EV-B85, and CVA9 serotypes in the nonstructural region. This study reports the isolation of E5 from SARI cases for the first time, and genetic characterization and evolutionary analysis provide valuable information on global E5 molecular epidemiology.IMPORTANCEThis study is the first report of echovirus 5 (E5) associated with severe acute respiratory infection and obtained the first E5 whole-genome sequence in China. Combined with the sequences available in the GenBank database, the first genotyping, phylogenetic characteristics, recombination, and genetic evolutionary analysis of E5 was performed in this study. Our findings providing valuable information on global E5 molecular epidemiology.

Global molecular evolution and phylogeographic analysis of barley yellow dwarf virus based on the cp and mp genes

Barley yellow dwarf virus (BYDV) has caused considerable losses in the global production of grain crops such as wheat, barley and maize. We investigated the phylodynamics of the virus by analysing 379 and 485 nucleotide sequences of the genes encoding the coat protein and movement protein, respectively. The maximum clade credibility tree indicated that BYDV-GAV and BYDV-MAV, BYDV-PAV and BYDV-PAS share the same evolutionary lineage, respectively. The diversification of BYDV arises from its adaptability to vector insects and geography. Bayesian phylogenetic analyses showed that the mean substitution rates of the coat and movement proteins of BYDV ranged from 8.327 × 10− 4 (4.700 × 10− 4–1.228 × 10− 3) and 8.671 × 10− 4 (6.143 × 10− 4–1.130 × 10− 3) substitutions/site/year, respectively. The time since the most recent common BYDV ancestor was 1434 (1040–1766) CE (Common Era). The Bayesian skyline plot (BSP) showed that the BYDV population experienced dramatic expansions approximately 8 years into the 21st century, followed by a dramatic decline in less than 15 years. Our phylogeographic analysis showed that the BYDV population originating in the United States was subsequently introduced to Europe, South America, Australia and Asia. The migration pathways of BYDV suggest that the global spread of BYDV is associated with human activities.

Re-Emerged Genotype IV of Japanese Encephalitis Virus Is the Youngest Virus in Evolution.

An outbreak of viral encephalitis caused by a Japanese encephalitis virus (JEV) genotype IV infection occurred in Australia between 2021 and 2022. A total of 47 cases and seven deaths were reported as of November 2022. This is the first outbreak of human viral encephalitis caused by JEV GIV since it was first isolated in Indonesia in the late 1970s. Here, a comprehensive phylogenetic analysis based on the whole genome sequences of JEVs revealed it emerged 1037 years ago (95% HPD: 463 to 2100 years). The evolutionary order of JEV genotypes is as follows: GV, GIII, GII, GI, and GIV. The JEV GIV emerged 122 years ago (95% HPD: 57-233) and is the youngest viral lineage. The mean substitution rate of the JEV GIV lineage was 1.145 × 10-3 (95% HPD values, 9.55 × 10-4, 1.35 × 10-3), belonging to rapidly evolving viruses. A series of amino acid mutations with the changes of physico-chemical properties located in the functional important domains within the core and E proteins distinguished emerging GIV isolates from old ones. These results demonstrate the JEV GIV is the youngest JEV genotype at a rapid evolution stage and has good host/vector adaptability for introduction to non-endemic areas. Thus, surveillance of JEVs is highly recommended.

Genetic characterization and molecular epidemiology of Coxsackievirus A12 from mainland China during 2010-2019.

Coxsackievirus A12 (CVA12) is an enterovirus that has been isolated in many countries in recent years. However, studies on CVA12 are limited, and its effective population size, evolutionary dynamics and recombination patterns have not been clarified now. In this study, we described the phylogenetic characteristics of 16 CVA12 strains isolated from pediatric HFMD patients in mainland China from 2010 to 2019. Comparison of the nucleotide sequences and amino acid sequences with the CVA12 prototype strain revealed that the 16 CVA12 strains are identical in 78.8-79% and 94-94.2%, respectively. A phylodynamic analysis based on the 16 full-length VP1 sequences from this study and 21 sequences obtained from GenBank revealed a mean substitution rate of 6.61 × 10-3 substitutions/site/year (95% HPD: 5.16-8.20 × 10-3), dating the time to most recent common ancestor (tMRCA) of CVA12 back to 1946 (95% HPD: 1942-1947). The Bayesian skyline plot showed that the effective population size has experienced twice dynamic fluctuations since 2007. Phylogeographic analysis identified two significant migration pathways, indicating the existence of cross-provincial transmission of CVA12 in mainland China. Recombination analysis revealed two recombination patterns between 16 CVA12 strains and other EV-A, suggesting that there may be extensive genetic exchange between CVA12 and other enteroviruses. In summary, a total of 16 full-length CVA12 strains were reported in this study, providing valuable references for further studies of CVA12 worldwide.

Identification of diffusion routes of O/EA-3 topotype of foot-and-mouth disease virus in Africa and Western Asia between 1974 and 2019 - a phylogeographic analysis.

Foot-and-mouth disease (FMD) affects the livestock industry and socioeconomic sustainability of many African countries. The success of FMD control programs in Africa depends largely on understanding the dynamics of FMD virus (FMDV) spread. In light of the recent outbreaks of FMD that affected the North-Western African countries in 2018 and 2019, we investigated the evolutionary phylodynamics of the causative serotype O viral strains all belonging to the East-Africa 3 topotype (O/EA-3). We analyzed a total of 489 sequences encoding the FMDV VP1 genome region generated from samples collected from 25 African and Western Asian countries between 1974 and 2019. Using Bayesian evolutionary models on genomic and epidemiological data, we inferred the routes of introduction and migration of the FMDV O/EA-3 topotype at the inter-regional scale. We inferred a mean substitution rate of 6.64×10-3 nt/site/year and we predicted that the most recent common ancestor for our panel of samples circulated between February 1967 and November 1973 in Yemen, likely reflecting the epidemiological situation in under sampled cattle-exporting East African countries. Our study also reinforces the role previously described of Sudan and South Sudan as a frequent source of FMDVs spread. In particular, we identified two transboundary routes of O/EA-3 diffusion: the first from Sudan to North-East Africa, and from the latter into Israel and Palestine AT; a second from Sudan to Nigeria, Cameroon, and from there to further into West and North-West Africa. This study highlights the necessity to reinforce surveillance at an inter-regional scale in Africa and Western Asia, in particular along the identified migration routes for the implementation of efficient control measures in the fight against FMD.

Evolutionary Dynamics of Begomoviruses and Its Satellites Infecting Papaya in India.

The genus Begomovirus represents a group of multipartite viruses that significantly damage many agricultural crops, including papaya, and influence overall production. Papaya leaf curl disease (PaLCD) caused by the complex begomovirus species has several important implications and substantial losses in papaya production in many developing countries, including India. The increase in the number of begomovirus species poses a continuous threat to the overall production of papaya. Here, we attempted to map the genomic variation, mutation, evolution rate, and recombination to know the disease complexity and successful adaptation of PaLCD in India. For this, we retrieved 44 DNA-A and 26 betasatellite sequences from GenBank reported from India. An uneven distribution of evolutionary divergence has been observed using the maximum-likelihood algorithm across the branch length. Although there were phylogenetic differences, we found high rates of nucleotide substitution mutation in both viral and sub-viral genome datasets. We demonstrated frequent recombination of begomovirus species, with a maximum in intra-species recombinants. Furthermore, our results showed a high degree of genetic variability, demographic selection, and mean substitution rate acting on the population, supporting the emergence of a diverse and purifying selection of viruses and associated betasatellites. Moreover, variation in the genetic composition of all begomovirus datasets revealed a predominance of nucleotide diversity principally driven by mutation, which might further accelerate the advent of new strains and species and their adaption to various hosts with unique pathogenicity. Therefore, the finding of genetic variation and selection emphases on factors that contribute to the universal spread and evolution of Begomovirus and this unanticipated diversity may also provide guidelines toward future evolutionary trend analyses and the development of wide-ranging disease control strategies for begomoviruses associated with PaLCD.

Molecular Epidemiology and Evolution of Coxsackievirus A9.

Nineteen CVA9 isolates were obtained between 2010 and 2019 from six provinces of mainland China, using the HFMD surveillance network established in China. Nucleotide sequencing revealed that the full-length VP1 of 19 CVA9 isolates was 906 bases encoding 302 amino acids. The combination of the thresholds of the phylogenetic tree and nucleotide divergence of different genotypes within the same serotype led to a value of 15–25%, and enabled CVA9 worldwide to be categorized into ten genotypes: A–J. The phylogenetic tree showed that the prototype strain was included in genotype A, and that the B, C, D, E, H, and J genotypes disappeared during virus evolution, whereas the F, I, and G genotypes showed co-circulation. Lineage G was the dominant genotype of CVA9 and included most of the strains from nine countries in Asia, North America, Oceania, and Europe. Most Chinese strains belonged to the G genotype, suggesting that the molecular epidemiology of China is consistent with that observed worldwide. The 165 partial VP1 strains (723 nt) showed a mean substitution rate of 3.27 × 10−3 substitution/site/year (95% HPD range 2.93–3.6 × 10−3), dating the tMRCA of CVA9 back to approximately 1922 (1911–1932). The spatiotemporal dynamics of CVA9 showed the spread of CVA9 obviously increased in recent years. Most CVA9 isolates originated in USA, but the epidemic areas of CVA9 are now concentrated in the Asia–Pacific region, European countries, and North America. Recombination analysis within the enterovirus B specie (59 serotypes) revealed eight recombination patterns in China at present, CVB4, CVB5, E30, CVB2, E11, HEV106, HEV85, and HEV75. E14, and E6 may act as recombinant donors in multiple regions. Comparison of temperature sensitivity revealed that temperature-insensitive strains have more amino acid substitutions in the RGD motif of the VP1 region, and the sites T283S, V284M, and R288K in the VP1 region may be related to the temperature tolerance of CVA9.

Purifying Selection Determines the Short-Term Time Dependency of Evolutionary Rates in SARS-CoV-2 and pH1N1 Influenza.

High-throughput sequencing enables rapid genome sequencing during infectious disease outbreaks and provides an opportunity to quantify the evolutionary dynamics of pathogens in near real-time. One difficulty of undertaking evolutionary analyses over short timescales is the dependency of the inferred evolutionary parameters on the timespan of observation. Crucially, there are an increasing number of molecular clock analyses using external evolutionary rate priors to infer evolutionary parameters. However, it is not clear which rate prior is appropriate for a given time window of observation due to the time-dependent nature of evolutionary rate estimates. Here, we characterize the molecular evolutionary dynamics of SARS-CoV-2 and 2009 pandemic H1N1 (pH1N1) influenza during the first 12 months of their respective pandemics. We use Bayesian phylogenetic methods to estimate the dates of emergence, evolutionary rates, and growth rates of SARS-CoV-2 and pH1N1 over time and investigate how varying sampling window and data set sizes affect the accuracy of parameter estimation. We further use a generalized McDonald–Kreitman test to estimate the number of segregating nonneutral sites over time. We find that the inferred evolutionary parameters for both pandemics are time dependent, and that the inferred rates of SARS-CoV-2 and pH1N1 decline by ∼50% and ∼100%, respectively, over the course of 1 year. After at least 4 months since the start of sequence sampling, inferred growth rates and emergence dates remain relatively stable and can be inferred reliably using a logistic growth coalescent model. We show that the time dependency of the mean substitution rate is due to elevated substitution rates at terminal branches which are 2–4 times higher than those of internal branches for both viruses. The elevated rate at terminal branches is strongly correlated with an increasing number of segregating nonneutral sites, demonstrating the role of purifying selection in generating the time dependency of evolutionary parameters during pandemics.

Viruses Without Borders: Global Analysis of the Population Structure, Haplotype Distribution, and Evolutionary Pattern of Iris Yellow Spot Orthotospovirus (Family Tospoviridae, Genus Orthotospovirus).

Iris yellow spot, caused by Iris yellow spot orthotospovirus (IYSV) (Genus: Orthotospovirus, Family: Tospoviridae), is an important disease of Allium spp. The complete N gene sequences of 142 IYSV isolates of curated sequence data from GenBank were used to determine the genetic diversity and evolutionary pattern. In silico restriction fragment length polymorphism (RFLP) analysis, codon-based maximum likelihood studies, genetic differentiation and gene flow within the populations of IYSV genotypes were investigated. Bayesian phylogenetic analysis was carried out to estimate the evolutionary rate. In silico RFLP analysis of N gene sequences categorized IYSV isolates into two major genotypes viz., IYSV Netherlands (IYSVNL; 55.63%), IYSV Brazil (IYSVBR; 38.73%) and the rest fell in neither group [IYSV other (IYSVother; 5.63%)]. Phylogenetic tree largely corroborated the results of RFLP analysis and the IYSV genotypes clustered into IYSVNL and IYSVBR genotypes. Genetic diversity test revealed IYSVother to be more diverse than IYSVNL and IYSVBR. IYSVNL and IYSVBR genotypes are under purifying selection and population expansion, whereas IYSVother showed decreasing population size and hence appear to be under balancing selection. IYSVBR is least differentiated from IYSVother compared to IYSVNL genotype based on nucleotide diversity. Three putative recombinant events were found in the N gene of IYSV isolates based on RDP analysis, however, RAT substantiated two among them. The marginal likelihood mean substitution rate was 5.08 × 10–5 subs/site/year and 95% highest posterior density (HPD) substitution rate between 5.11 × 10–5 and 5.06 × 10–5. Findings suggest that IYSV continues to evolve using population expansion strategies. The substitution rates identified are similar to other plant RNA viruses.

Evolution of the GII.3[P12] Norovirus from 2010 to 2019 in Jiangsu, China

ObjectivesNorovirus genotype GII.3[P12] strains have been an important pathogen for sporadic gastroenteritis infection. In previous studies of GII.3[P12], the number of specimens and time span are relatively small, which is difficult to truly reflect the infection and evolution of this type of norovirus. Here we report a molecular epidemiological study of the NoVs prevalent in Jiangsu between 2010 and 2019 to investigate the evolution of the GII.3[P12] strains in China.MethodsIn this study 60 GII.3[P12] norovirus strains were sequenced and analyzed for evolution, recombination, and selection pressure using bioanalysis software.ResultsThe GII.3[P12] strains were continuously detected during the study period, which showed a high constituent ratio in males, in winter and among children aged 0–11 months, respectively. A time-scaled evolutionary tree showed that both GII.P12 RdRp and GII.3 VP1 sequences were grouped into three major clusters (Cluster I–III). Most GII.3[P12] strains were mainly located in sub-cluster (SC) II of Cluster III. A SimPlot analysis identified GII.3[P12] strain to be as an ORF1-intragenic recombinant of GII.4[P12] and GII.3[P21]. The RdRp genes of the GII.3[P12] showed a higher mean substitution rate than those of all GII.P12, while the VP1 genes of the GII.3[P12] showed a lower mean substitution rate than those of all GII.3. Alignment of the GII.3 capsid sequences revealed that three HBGA binding sites of all known GII.3 strains remained conserved, while several amino acid mutations in the predicted antibody binding sites were detected. The mutation at 385 was within predicted antibody binding regions, close to host attachment factor binding sites. Positive and negative selection sites were estimated. Two common positively selected sites (sites 385 and 406) were located on the surface of the protruding domain. Moreover, an amino acid substitution (aa204) was estimated to be near the active site of the RdRp protein.ConclusionsWe conducted a comprehensive analysis on the epidemic and evolution of GII.3[P12] noroviruses and the results suggested that evolution was possibly driven by intergenic recombination and mutations in some key amino acid sites.

Isolation and evolutionary analysis of Senecavirus A isolates from Guangdong province, China.

Senecavirus A (SVA), an emerging swine pathogen, has been reported in many provinces of China since the first outbreak in 2015 in Guangdong province. In this study, 10 lymph nodes positive for SVA, collected between 2018 and 2019 from slaughterhouses in Guangdong province, were subjected to virus isolation. Rapid and evident cytopathic effects (CPEs) were observed in SVA-infected PK-15 cells, including shrinking, rounding and detaching, with peak titers being reached at 24h post infection (hpi). Electron microscopy showed that SVA particles are spherical and approximately 30nm in diameter, and exist as crystalline lattices in cytoplasm revealed by ultra-thin sectioning. Phylogenetic analysis based on the whole genome sequences of all available isolates showed that SVA globally can be divided into two groups with each being further divided into two subgroups (Ia-b and IIa-b), and with the Guangdong isolates obtained here and other Chinese strains belonging to subgroups IIa and IIb. Evolutionary analysis showed that the mean substitution rate of SVA was 2.696×10-3 per site per year based on whole genomic sequences, with subgroup IIb isolates having evolved faster than those of subgroup IIa. Analysis of efficient population size showed that the outbreak point of SVA worldwide occurred at the end of 2013 with that of subgroup IIb, the current dominant group, in mid 2014.

Phylogeographic analysis of Pseudogymnoascus destructans partitivirus-pa explains the spread dynamics of white-nose syndrome in North America.

Understanding the dynamics of white-nose syndrome spread in time and space is an important component for the disease epidemiology and control. We reported earlier that a novel partitivirus, Pseudogymnoascus destructans partitivirus-pa, had infected the North American isolates of Pseudogymnoascus destructans, the fungal pathogen that causes white-nose syndrome in bats. We showed that the diversity of the viral coat protein sequences is correlated to their geographical origin. Here we hypothesize that the geographical adaptation of the virus could be used as a proxy to characterize the spread of white-nose syndrome. We used over 100 virus isolates from diverse locations in North America and applied the phylogeographic analysis tool BEAST to characterize the spread of the disease. The strict clock phylogeographic analysis under the coalescent model in BEAST showed a patchy spread pattern of white-nose syndrome driven from a few source locations including Connecticut, New York, West Virginia, and Kentucky. The source states had significant support in the maximum clade credibility tree and Bayesian stochastic search variable selection analysis. Although the geographic origin of the virus is not definite, it is likely the virus infected the fungus prior to the spread of white-nose syndrome in North America. We also inferred from the BEAST analysis that the recent long-distance spread of the fungus to Washington had its root in Kentucky, likely from the Mammoth cave area and most probably mediated by a human. The time to the most recent common ancestor of the virus is estimated somewhere between the late 1990s to early 2000s. We found the mean substitution rate of 2 X 10−3 substitutions per site per year for the virus which is higher than expected given the persistent lifestyle of the virus, and the stamping-machine mode of replication. Our approach of using the virus as a proxy to understand the spread of white-nose syndrome could be an important tool for the study and management of other infectious diseases.

Comparative Analyses of the Chloroplast Genomes of Patchouli Plants and Their Relatives in Pogostemon (Lamiaceae).

Pogostemon Desf., the largest genus of the tribe Pogostemoneae (Lamiaceae), consists of ca. 80 species distributed mainly from South and Southeast Asia to China. The genus contains many patchouli plants, which are of great economic importance but taxonomically difficult. Therefore, it is necessary to characterize more chloroplast (cp) genomes for infrageneric phylogeny analyses and species identification of Pogostemon, especially for patchouli plants. In this study, we newly generated four cp genomes for three patchouli plants (i.e., Pogostemon plectranthoides Desf., P. septentrionalis C. Y. Wu et Y. C. Huang, and two cultivars of P. cablin (Blanoco) Benth.). Comparison of all samples (including online available cp genomes of P. yatabeanus (Makino) Press and P. stellatus (Lour.) Kuntze) suggested that Pogostemon cp genomes are highly conserved in terms of genome size and gene content, with a typical quadripartite circle structure. Interspecific divergence of cp genomes has been maintained at a relatively low level, though seven divergence hotspot regions were identified by stepwise window analysis. The nucleotide diversity (Pi) value was correlated significantly with gap proportion (indels), but significantly negative with GC content. Our phylogenetic analyses based on 80 protein-coding genes yielded high-resolution backbone topologies for the Lamiaceae and Pogostemon. For the overall mean substitution rates, the synonymous (dS) and nonsynonymous (dN) substitution rate values of protein-coding genes varied approximately threefold, while the dN values among different functional gene groups showed a wider variation range. Overall, the cp genomes of Pogostemon will be useful for phylogenetic reconstruction, species delimitation and identification in the future.

DNA barcoding reveals mislabeling and commercial fraud in the marketing of fillets of the genus Brachyplatystoma Bleeker, 1862, the Amazonian freshwater catfishes economically important in Brazil

The substitution and mislabeling is facilitated by the processing of fish products. We employed a DNA barcoding to authenticate fillets labeled as “dourada” (Brachyplatystoma rousseauxii), and “piramutaba” (Brachyplatystoma vaillantii) marketed in the Brazil. A 615 bp of the Cytochrome oxidase subunit I (COI) was sequenced from 305 fillets and subsequently identified to species level by querying public databases and sequences of reference species. The results revealed a global mean substitution rate of 17%. The highest substitution rate was detected in “dourada” (26%), the most valuable species, followed by “piramutaba” (9%). The most cases of substitutions were by species of lower commercial value, suggesting fraud aimed at increased profits. Therefore, we suggest the improvement of food-labeling regulation, continued inspection, as well as the adoption of the DNA barcode for the molecular authentication of processed fish to prevent substitution of these products in Brazil.

Differences of cytotoxic T-lymphocyte pressure and divergent evolution of several CRF07_BC clusters circulating in men who have sex with men in China

The HIV-1 CRF07_BC strain is now the major recombinant form in China among the population of men who have sex with men (MSM), and has critically contributed to the HIV-1 epidemic in recent years. The phylodynamic and virological differences among CRF07_BC clusters circulating in MSM, and the factors that could be driving their evolution, remains unclear. Using the HIV-1 CRF07_BC strains obtained from the Los Alamos HIV database, we undertook a large-scale phylogenetic analysis using the maximum likelihood method of partial gag, pol, and env segments to infer their evolutionary relationships. The demographic histories of clusters were determined using the Bayesian Markov chain Monte Carlo (MCMC) method. For four pol clusters we analyzed the non-synonymous (dN) to synonymous (dS) substitution rates and performed site to site analysis to identify positive selection sites and cytotoxic T-lymphocyte (CTL) escape mutation positions. MSM was found to be the predominant risk factor for all four of the CRF07_BC epidemic pol segment clusters with the largest number of infections. Two of those clusters had higher growth in the effective number of infections, and two clusters had slower growth. Analysis of all four clusters showed no significant differences in the mean substitution rates and dN/dS selection pressure ratios. However, a site to site codon analysis found thirteen significant positive selection sites. Ten of these sites are CTL escape mutations. The two clusters with the higher growth in infections had seven and eight pol segment CTL escape mutation sites respectively, while the two with slower growth had only one or two. Our findings demonstrated differences in the CTL escape mutation and divergent evolution of several CRF07_BC clusters circulating among men who have sex with men in China.