Phylodynamic Analysis Research Articles

Genetic Diversity and Comprehensive Phylodynamic Analysis of Whole Corona Viruses Genomes

Covid-19 is a rapidly emerging viral respiratory disease. The infection of corona virus disease (COVID-19) is initiated by the strains of severe acute respiratory syndrome. Corona virus 2 (SARS-CoV-2), is a highly transmittable viral infection that was first emerged in the of Wuhan, a city of China and ultimately extent in whole world. The World health organization (WHO) affirmed this infection as a global public health disaster. Genomic analysis showing that this novel SARSCoV-2 is hereditarily linked with severe acute respiratory syndrome-like (SARS-like) bat viruses indicating that bats might be the primary source of this epidemic infection. In this review we have summarized a comprehensive phylodynamic analysis of whole corona virus genome sequences in order to analyze the pathogenicity and emergence of novel COVID-19 infection and previous human coronaviruses. Our main objective is to describe the genetic relationships of the 2019-nCoV with other genera of coronavirus and to find out supposed recombination contained by sarbecovirus subgenus. Keywords: COVID-19, Phylogeny, SARS-CoV, MERS-CoV, Coding Sequence (CDS).

Integration of genomic sequencing into the response to the Ebola virus outbreak in Nord Kivu, Democratic Republic of the Congo.

On 1 August 2018, the Democratic Republic of the Congo (DRC) declared its tenth Ebola virus disease (EVD) outbreak. To aid the epidemiologic response, the Institut National de Recherche Biomédicale (INRB) implemented an end-to-end genomic surveillance system, including sequencing, bioinformatic analysis and dissemination of genomic epidemiologic results to frontline public health workers. We report 744 new genomes sampled between 27 July 2018 and 27 April 2020 generated by this surveillance effort. Together with previously available sequence data (n = 48 genomes), these data represent almost 24% of all laboratory-confirmed Ebola virus (EBOV) infections in DRC in the period analyzed. We inferred spatiotemporal transmission dynamics from the genomic data as new sequences were generated, and disseminated the results to support epidemiologic response efforts. Here we provide an overview of how this genomic surveillance system functioned, present a full phylodynamic analysis of 792 Ebola genomes from the Nord Kivu outbreak and discuss how the genomic surveillance data informed response efforts and public health decision making.

Field validation of phylodynamic analytical methods for inference on epidemiological processes in wildlife.

Amongst newly developed approaches to analyse molecular data, phylodynamic models are receiving much attention because of their potential to reveal changes to viral populations over short periods. This knowledge can be very important for understanding disease impacts. However, their accuracy needs to be fully understood, especially in relation to wildlife disease epidemiology, where sampling and prior knowledge may be limited. The release of the rabbit haemorrhagic disease virus (RHDV) as biological control in naïve rabbit populations in Australia in 1996 provides a unique data set with which to validate phylodynamic models. By comparing results obtained from RHDV sequence data with our current understanding of RHDV epidemiology in Australia, we evaluated the performances of these recently developed models. In line with our expectations, coalescent analyses detected a sharp increase in the virus population size in the first few months after release, followed by a more gradual increase. Phylodynamic analyses using a birth-death model generated effective reproductive number estimates (the average number of secondary infections per each infectious case, Re ) larger than one for most of the epochs considered. However, the possible range of the initial Re included estimates lower than one despite the known rapid spread of RHDV in Australia. Furthermore, the analyses that accounted for geographical structuring failed to converge. We argue that the difficulties that we encountered most likely stem from the fact that the samples available from 1996 to 2014 were too sparse with respect to both geographic and within outbreak coverage to adequately infer some of the model parameters. In general, while these phylodynamic analyses proved to be greatly informative in some regards, we caution that their interpretation may not be straightforward. We recommend further research to evaluate the robustness of these models to assumption violations and sensitivity to sampling regimes.

Phylogenomic analysis of COVID-19 summer and winter outbreaks in Hong Kong: An observational study

BackgroundViral genomic surveillance is vital for understanding the transmission of COVID-19. In Hong Kong, breakthrough outbreaks have occurred in July (third wave) and November (fourth wave) 2020. We used whole viral genome analysis to study the characteristics of these waves. MethodsWe analyzed 509 SARS-CoV-2 genomes collected from Hong Kong patients between 22nd January and 29th November, 2020. Phylogenetic and phylodynamic analyses were performed, and were interpreted with epidemiological information. FindingsDuring the third and fourth waves, diverse SARS-CoV-2 genomes were identified among imported infections. Conversely, local infections were dominated by a single lineage during each wave, with 96.6% (259/268) in the third wave and 100% (73/73) in the fourth wave belonging to B.1.1.63 and B.1.36.27 lineages, respectively. While B.1.1.63 lineage was imported 2 weeks before the beginning of the third wave, B.1.36.27 lineage has circulated in Hong Kong for 2 months prior to the fourth wave. During the fourth wave, 50.7% (37/73) of local infections in November was identical to the viral genome from an imported case in September. Within B.1.1.63 or B.1.36.27 lineage in our cohort, the most common non-synonymous mutations occurred at the helicase (nsp13) gene. InterpretationAlthough stringent measures have prevented most imported cases from spreading in Hong Kong, a single lineage with low-level local transmission in October and early November was responsible for the fourth wave. A superspreading event or lower temperature in November may have facilitated the spread of the B.1.36.27 lineage. FundingRichard and Carol Yu, Michael Tong, and the Government Consultancy Service (see acknowledgments for full list).

Molecular Evolution of Human Norovirus GII.2 Clusters

BackgroundThe human norovirus GII.2 outbreak during the 2016–2017 winter season was of unprecedented scale and geographic distribution.MethodsWe analyzed 519 complete VP1 gene sequences of the human norovirus GII.2 genotype sampled during the 2016–2017 winter season, as well as prior (dating back to 1976) from 7 countries. Phylodynamic analyses of these sequences were performed using maximum likelihood and Bayesian statistical frameworks in order to estimate viral evolutionary and population dynamics associated with the outbreak.ResultsOur results revealed an increase in the genetic diversity of human norovirus GII.2 during the recent Asian outbreak and diversification was characterized by at least eight distinct clusters. Bayesian estimation of viral population dynamics revealed a highly fluctuating effective population size, increasing in frequency during the past 15 years.ConclusionDespite an increasing viral diversity, we found no evidence of an elevated evolutionary rate or significant selection pressure in human norovirus GII.2, indicating viral evolutionary adaptation was not responsible for the volatility of or spread of the virus during this time.

Multiple Reassortants of H5N8 Clade 2.3.4.4b Highly Pathogenic Avian Influenza Viruses Detected in South Korea during the Winter of 2020-2021.

During October 2020–January 2021, we isolated a total of 67 highly pathogenic avian influenza (HPAI) H5N8 viruses from wild birds and outbreaks in poultry in South Korea. We sequenced the isolates and performed phylogenetic analysis of complete genome sequences to determine the origin, evolution, and spread patterns of these viruses. Phylogenetic analysis of the hemagglutinin (HA) gene showed that all the isolates belong to H5 clade 2.3.4.4 subgroup B (2.3.4.4b) and form two distinct genetic clusters, G1 and G2. The cluster G1 was closely related to the 2.3.4.4b H5N8 HPAI viruses detected in Europe in early 2020, while the cluster G2 had a close genetic relationship with the 2.3.4.4b H5N8 viruses that circulated in Europe in late 2020. A total of seven distinct genotypes were identified, including five novel reassortants carrying internal genes of low pathogenic avian influenza viruses. Our Bayesian discrete trait phylodynamic analysis between host types suggests that the viruses initially disseminated from migratory waterfowl to domestic duck farms in South Korea. Subsequently, domestic duck farms most likely contributed to the transmission of HPAI viruses to chicken and minor poultry farms, highlighting the need for enhanced, high levels of biosecurity measures at domestic duck farms to effectively prevent the introduction and spread of HPAI.

SARS-CoV-2 Molecular Transmission Clusters and Containment Measures in Ten European Regions during the First Pandemic Wave.

Background: The spatiotemporal profiling of molecular transmission clusters (MTCs) using viral genomic data can effectively identify transmission networks in order to inform public health actions targeting SARS-CoV-2 spread. Methods: We used whole genome SARS-CoV-2 sequences derived from ten European regions belonging to eight countries to perform phylogenetic and phylodynamic analysis. We developed dedicated bioinformatics pipelines to identify regional MTCs and to assess demographic factors potentially associated with their formation. Results: The total number and the scale of MTCs varied from small household clusters identified in all regions, to a super-spreading event found in Uusimaa-FI. Specific age groups were more likely to belong to MTCs in different regions. The clustered sequences referring to the age groups 50–100 years old (y.o.) were increased in all regions two weeks after the establishment of the lockdown, while those referring to the age group 0–19 y.o. decreased only in those regions where schools’ closure was combined with a lockdown. Conclusions: The spatiotemporal profiling of the SARS-CoV-2 MTCs can be a useful tool to monitor the effectiveness of the interventions and to reveal cryptic transmissions that have not been identified through contact tracing.

Genomic Epidemiology of SARS-CoV-2 in Madrid, Spain, during the First Wave of the Pandemic: Fast Spread and Early Dominance by D614G Variants.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first detected in Madrid, Spain, on 25 February 2020. It increased in frequency very fast and by the end of May more than 70,000 cases had been confirmed by reverse transcription-polymerase chain reaction (RT-PCR). To study the lineages and the diversity of the viral population during this first epidemic wave in Madrid we sequenced 224 SARS-CoV-2 viral genomes collected from three hospitals from February to May 2020. All the known major lineages were found in this set of samples, though B.1 and B.1.5 were the most frequent ones, accounting for more than 60% of the sequences. In parallel with the B lineages and sublineages, the D614G mutation in the Spike protein sequence was detected soon after the detection of the first coronavirus disease 19 (COVID-19) case in Madrid and in two weeks became dominant, being found in 80% of the samples and remaining at this level during all the study periods. The lineage composition of the viral population found in Madrid was more similar to the European population than to the publicly available Spanish data, underlining the role of Madrid as a national and international transport hub. In agreement with this, phylodynamic analysis suggested multiple independent entries before the national lockdown and air transportation restrictions.

Time and Mode of Epidemic HCV-2 Subtypes Spreading in Europe: Phylodynamics in Italy and Albania.

Hepatitis C virus (HCV) genotype 2 causes about 10% of global infections and has the most variable circulation profile in Europe. The history of “endemic” HCV-2 subtypes has been satisfactorily reconstructed, instead there is little information about the recent spread of the “epidemic” subtypes, including HCV-2c. To investigate the origin and dispersion pathways of HCV-2c, 245 newly characterized Italian and Albanian HCV-2 NS5B sequences were aligned with 247 publicly available sequences and included in phylogeographic and phylodynamic analyses using the Bayesian framework. Our findings show that HCV-2c was the most prevalent subtype in Italy and Albania. The phylogeographic analysis suggested an African origin of HCV-2c before it reached Italy about in the 1940s. Phylodynamic analysis revealed an exponential increase in the effective number of infections and Re in Italy between the 1940s and 1960s, and in Albania between the 1990s and the early 2000s. It seems very likely that HCV-2c reached Italy from Africa at the time of the second Italian colonization but did not reach Albania until the period of dramatic migration to Italy in the 1990s. This study contributes to reconstructing the history of the spread of epidemic HCV-2 subtypes to Europe.

Evolutionary Dynamics and Dissemination Pattern of the SARS-CoV-2 Lineage B.1.1.33 During the Early Pandemic Phase in Brazil.

A previous study demonstrates that most of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) Brazilian strains fell in three local clades that were introduced from Europe around late February 2020. Here we investigated in more detail the origin of the major and most widely disseminated SARS-CoV-2 Brazilian lineage B.1.1.33. We recovered 190 whole viral genomes collected from 13 Brazilian states from February 29 to April 31, 2020 and combined them with other B.1.1 genomes collected globally. Our genomic survey confirms that lineage B.1.1.33 is responsible for a variable fraction of the community viral transmissions in Brazilian states, ranging from 2% of all SARS-CoV-2 genomes from Pernambuco to 80% of those from Rio de Janeiro. We detected a moderate prevalence (5–18%) of lineage B.1.1.33 in some South American countries and a very low prevalence (<1%) in North America, Europe, and Oceania. Our study reveals that lineage B.1.1.33 evolved from an ancestral clade, here designated B.1.1.33-like, that carries one of the two B.1.1.33 synapomorphic mutations. The B.1.1.33-like lineage may have been introduced from Europe or arose in Brazil in early February 2020 and a few weeks later gave origin to the lineage B.1.1.33. These SARS-CoV-2 lineages probably circulated during February 2020 and reached all Brazilian regions and multiple countries around the world by mid-March, before the implementation of air travel restrictions in Brazil. Our phylodynamic analysis also indicates that public health interventions were partially effective to control the expansion of lineage B.1.1.33 in Rio de Janeiro because its median effective reproductive number (Re) was drastically reduced by about 66% during March 2020, but failed to bring it to below one. Continuous genomic surveillance of lineage B.1.1.33 might provide valuable information about epidemic dynamics and the effectiveness of public health interventions in some Brazilian states.

Transmission dynamics of SARS-CoV-2 on the Diamond Princess uncovered using viral genome sequence analysis.

An outbreak of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) occurred aboard the Diamond Princess cruise ship between her January 20 departure and late February 2020. Here, we used phylodynamic analyses to investigate the transmission dynamics of SARS-CoV-2 during the outbreak. Using a Bayesian coalescent-based method, the estimated mean nucleotide substitution rate of 240 SARS-CoV-2 whole-genome sequences was approximately 7.13 × 10−4 substitutions per site per year. Population dynamics and the effective reproductive number (Re) of SARS-CoV-2 infections were estimated using a Bayesian framework. The estimated origin of the outbreak was January 21, 2020. The infection spread substantially before quarantine on February 5. The Re peaked at 6.06 on February 4 and gradually declined to 1.51, suggesting that transmission continued slowly even after quarantine. These findings highlight the high transmissibility of SARS-CoV-2 and the need for effective measures to control outbreaks in confined settings.

SARS-CoV-2 mutation 614G creates an elastase cleavage site enhancing its spread in high AAT-deficient regions.

SARS-CoV-2 was first reported from China. Within three months, it evolved to 10 additional subtypes. Two evolved subtypes (A2 and A2a) carry a non-synonymous Spike protein mutation (D614G). We conducted phylodynamic analysis of over 70,000 SARS-CoV-2 coronaviruses worldwide, sequenced until July2020, and found that the mutant subtype (614G) outcompeted the pre-existing type (614D), significantly faster in Europe and North-America than in East Asia. Bioinformatically and computationally, we identified a novel neutrophil elastase (ELANE) cleavage site introduced in the G-mutant, near the S1-S2 junction of the Spike protein. We hypothesised that elevation of neutrophil elastase level at the site of infection will enhance the activation of Spike protein thus facilitating host cell entry for 614G, but not the 614D, subtype. The level of neutrophil elastase in the lung is modulated by its inhibitor α1-antitrypsin (AAT). AAT prevents lung tissue damage by elastase. However, many individuals exhibit genotype-dependent deficiency of AAT. AAT deficiency eases host-cell entry of the 614G virus, by retarding inhibition of neutrophil elastase and consequently enhancing activation of the Spike protein. AAT deficiency is highly prevalent in European and North-American populations, but much less so in East Asia. Therefore, the 614G subtype is able to infect and spread more easily in populations of the former regions than in the latter region. Our analyses provide a molecular biological and evolutionary model for the higher observed virulence of the 614G subtype, in terms of causing higher morbidity in the host (higher infectivity and higher viral load), than the non-mutant 614D subtype.

Genomic epidemiology of the early stages of the SARS-CoV-2 outbreak in Russia

The ongoing pandemic of SARS-CoV-2 presents novel challenges and opportunities for the use of phylogenetics to understand and control its spread. Here, we analyze the emergence of SARS-CoV-2 in Russia in March and April 2020. Combining phylogeographic analysis with travel history data, we estimate that the sampled viral diversity has originated from at least 67 closely timed introductions into Russia, mostly in late February to early March. All but one of these introductions were not from China, suggesting that border closure with China has helped delay establishment of SARS-CoV-2 in Russia. These introductions resulted in at least 9 distinct Russian lineages corresponding to domestic transmission. A notable transmission cluster corresponded to a nosocomial outbreak at the Vreden hospital in Saint Petersburg; phylodynamic analysis of this cluster reveals multiple (2-3) introductions each giving rise to a large number of cases, with a high initial effective reproduction number of 3.0 [1.9, 4.3].

Analysis of Hepatitis B Virus Genotype D in Greenland Suggests the Presence of a Novel Quasi-Subgenotype.

A disproportionate number of Greenland's Inuit population are chronically infected with Hepatitis B virus (HBV; 5–10%). HBV genotypes B and D are most prevalent in the circumpolar Arctic. Here, we report 39 novel HBV/D sequences from individuals residing in southwestern Greenland. We performed phylodynamic analyses with ancient HBV DNA calibrators to investigate the origin and relationship of these taxa to other HBV sequences. We inferred a substitution rate of 1.4 × 10−5 [95% HPD 8.8 × 10−6, 2.0 × 10−5] and a time to the most recent common ancestor of 629 CE [95% HPD 37–1138 CE]. The Greenland taxa form a sister clade to HBV/D2 sequences, specifically New Caledonian and Indigenous Taiwanese sequences. The Greenland sequences share amino acid signatures with subgenotypes D1 and D2 and ~97% sequence identity. Our results suggest the classification of these novel sequences does not fit within the current nomenclature. Thus, we propose these taxa be considered a novel quasi-subgenotype.

Genomics-Informed Responses in the Elimination of COVID-19 in Australia

BACKGROUND: A cornerstone of Australia’s ability to control COVID-19 has been effective border control, using an extensive supervised quarantine program. However, a rapid recrudescence in COVID-19 cases was observed in the state of Victoria in June 2020. Here, we describe the genomic findings that located the source of this second wave as a breach in supervised hotel quarantine and demonstrate the successful elimination of COVID-19 for a second time in Australia.METHODS: Genome sequencing was performed on all available SARS-CoV-2-positive samples in Victoria and integrated genomic and epidemiological investigation undertaken.RESULTS: At 31st January 2021, 20,451 COVID-19 cases were reported in Victoria; samples were sequenced from 75% of cases (15,431/20,451). Genomics revealed 98% (10,426/10,646) of locally-acquired cases during the second wave were derived from a single incursion from hotel quarantine, with the outbreak strain rapidly detected in other Australian states and territories. Phylodynamic analyses indicated an epidemic growth rate comparable to emerging variants, such as B.1.1.7 in the United Kingdom. Strict public health interventions resulted in the elimination of the outbreak strain by 29th October 2020. Subsequent cases represented independent international or interstate introductions, with limited local spread.CONCLUSIONS: Rapid escalation of clonal outbreaks can occur from even a single breach of control practices, as revealed through our genomic ‘enhanced outbreak-detection' system. The subsequent elimination and rapid control of new SARS-CoV-2 incursions reinforce that decisive public health responses to emergent cases are effective even with high epidemic growth rates, and “elimination” should be favored in settings where this is achievable.FUNDING STATEMENT: The Microbiological Diagnostic Unit Public Health Laboratory (MDU PHL) and the Victorian Infectious Diseases Reference Laboratory (VIDRL) at The Doherty Institute are funded by the Victorian Government. This work was supported by the National Health and Medical Research Council, Australia (NHMRC); Partnership Grant (APP1149991), Investigator Grant to BPH (APP1196103), Investigator Grant to DAW (APP1174555), Research Fellowship to TPS (APP1105525), MRFF COVID-19 Genomics Grant (MRF9200006).DECLARATION OF INTERESTS: None to declare. ETHICS APPROVAL STATEMENT: Data were collected in accordance with the Victorian Public Health and Wellbeing Act 2008. Ethical approval was received from the University of Melbourne Human Research Ethics Committee (study number 1954615.3).

Nation-Wide Viral Sequence Analysis of HIV-1 Subtype B Epidemic in 2003-2012 Revealed a Contribution of Men Who Have Sex With Men to the Transmission Cluster Formation and Growth in Japan.

Background: To better understand the epidemiology of human immunodeficiency virus type 1 (HIV-1) subtype B transmission in Japan, phylodynamic analysis of viral pol sequences was conducted on individuals newly diagnosed as HIV-1 seropositive.Methodology: A total of 5,018 patients newly diagnosed with HIV-1 infection and registered in the Japanese Drug Resistance HIV Surveillance Network from 2003 to 2012 were enrolled in the analysis. Using the protease-reverse transcriptase nucleotide sequences, their subtypes were determined, and phylogenetic relationships among subtype B sequences were inferred using three different methods: distance-matrix, maximum likelihood, and Bayesian Markov chain Monte Carlo. Domestically spread transmission clusters (dTCs) were identified based on the following criteria: >95% in interior branch test, >95% in Bayesian posterior probability and <10% in depth-first searches for sub-tree partitions. The association between dTC affiliation and individuals' demographics was analyzed using univariate and multivariate analyses.Results: Among the cases enrolled in the analysis, 4,398 (87.6%) were classified as subtype B. Many of them were Japanese men who had sex with men (MSM), and 3,708 (84.3%) belonged to any of 312 dTCs. Among these dTCs, 243 (77.9%) were small clusters with <10 individuals, and the largest cluster consisted of 256 individuals. Most dTCs had median time of the most recent common ancestor between 1995 and 2005, suggesting that subtype B infection was spread among MSMs in the second half of the 1990s. Interestingly, many dTCs occurred within geographical regions. Comparing with singleton cases, TCs included more MSM, young person, and individuals with high CD4+ T-cell count at the first consultation. Furthermore, dTC size was significantly correlated with gender, age, transmission risks, recent diagnosis and relative population size of the region mainly distributed.Conclusions: Our study clarified that major key population of HIV-1 subtype B epidemic in Japan is local MSM groups. The study suggests that HIV-1 subtype B spread via episodic introductions into the local MSM groups, some of the viruses spread to multiple regions. Many cases in dTC were diagnosed during the early phase of infection, suggesting their awareness to HIV risks.

Specific capture and whole-genome phylogeography of Dolphin morbillivirus

Dolphin morbillivirus (DMV) is considered an emerging threat having caused several epidemics worldwide. Only few DMV genomes are publicly available. Here, we report the use of target enrichment directly from cetacean tissues to obtain novel DMV genome sequences, with sequence comparison and phylodynamic analysis. RNA from 15 tissue samples of cetaceans stranded along the Italian and French coasts (2008–2017) was purified and processed using custom probes (by bait hybridization) for target enrichment and sequenced on Illumina MiSeq. Data were mapped against the reference genome, and the novel sequences were aligned to the available genome sequences. The alignment was then used for phylogenetic and phylogeographic analysis using MrBayes and BEAST. We herein report that target enrichment by specific capture may be a successful strategy for whole-genome sequencing of DMV directly from field samples. By this strategy, 14 complete and one partially complete genomes were obtained, with reads mapping to the virus up to 98% and coverage up to 7800X. The phylogenetic tree well discriminated the Mediterranean and the NE-Atlantic strains, circulating in the Mediterranean Sea and causing two different epidemics (2008–2015 and 2014–2017, respectively), with a limited time overlap of the two strains, sharing a common ancestor approximately in 1998.

Phylodynamic Analysis and Implication of HCV Genotype 4 Variability on Antiviral Drug Response and T-Cell Recognition.

Therapies for HCV care could change the prevalence and the geographic distribution of genotypes due to differences in Sustained Virologic Response (SVR). In this scenario, uncommon genotypes/subtypes, such as genotype 4, could spread from high-risk groups, replacing genotypes eradicated by antiviral drugs. Genotype eradication is also strongly influenced by the CD8+ T cell response. In this study, the genetic variability in HCV genotype 4 strains obtained from a cohort of 67 patients naïve to DAA therapy was evaluated. We found that the presence of resistance-associated substitutions (RAS) was able to affect drug responses. Next, using a prediction tool, viral mutations were identified by their ability, or lack thereof, to reduce the binding affinity with HLA, which affects T cell recognition. The Bayesian coalescent analysis suggested two different circulation clusters, one in risk groups (IDUs and MSM) and the other due to migration flows, dated to 1940 and 1915, respectively. Most of the RAS overlapped with HLA and a lack of binding mutations was observed in 96% of strains. This study describes the introduction of HCV genotype 4 in a region of the Mediterranean basin and evaluates how HCV genotype 4’s genetic variability could affect the response of antiviral drugs and CD8+ T cell recognition.

Genomic and phylodynamic analysis of sapoviruses isolated in Henan Province, China.

In this study, we determined the near-complete and partial genome sequences of ten SaV isolates. Phylogenetic analysis based on full-length VP1 and RdRp nucleotide sequences indicated that nine isolates were of GI.1 and one was GII.3. Evolutionary dynamics analysis indicated that GI.1 and GII.3 SaVs evolved at different rates, the latter evolving more rapidly. Cluster analysis indicated that distantly related GI.1 SaVs were more similar in their amino acid compositions than were GII.3 SaVs. The data provided in this study may facilitate studies on SaV genomic diversity and epidemiological patterns in China and worldwide.

Genomic Characterization of Salmonella Minnesota Clonal Lineages Associated with Poultry Production in Brazil.

Simple SummarySalmonella is a leading cause of foodborne illnesses and a global public health concern. Salmonella enterica serotype Minnesota has been increasingly isolated from Brazilian poultry farms. The present study investigated the phylogenetic relationships, evolution and genetic characteristics of S. Minnesota isolates from Brazilian poultry farms. The results demonstrated two main S. Minnesota lineages in the poultry production chain from Brazil, both presenting genes for antibiotic resistance and virulence. The present study also provides insights into the temporal evolution, population structure, and genetic characteristics of the two S. Minnesota lineages disseminated in Brazilian poultry farms.Salmonella serotype Minnesota has been increasingly detected in Brazilian poultry farms and food products (chicken meat, eggs) in recent years. In addition, S. Minnesota isolates from poultry are generally resistant to several antibiotics and persistent in farm environments. The present study aimed to assess phylogenomic diversity of S. Minnesota isolates from the poultry production chain in Brazil. In total, 107 worldwide S. Minnesota whole genomes (including 12 from Brazil) were analyzed using a comparative approach. Phylogenetic analysis demonstrated two clades more related to poultry production in Brazil: S. Minnesota poultry lineages I and II (SM-PLI and SM-PLII). Phylodynamic analysis demonstrated that SM-PLI had a common ancestor in 1915, while SM–PLII originated circa 1971. SM-PLII encompassed a higher number of isolates and presented a recent increase in effective population size (mainly from 2009 to 2012). Plasmids IncA/C2 and ColRNA, antimicrobial resistance genes (aph(3′)-Ia, blaCMY-2, qnrB19, sul2, and tet(A)) and mainly a virulence genetic cluster (including the yersiniabactin operon) were detected in isolates from SM-PLI and/or SM-PLII. This study demonstrates the dissemination of two distinct S. Minnesota lineages with high resistance to antibiotics and important virulence genetic clusters in Brazilian poultry farms.