Evolution and genomic profile of Salmonella enterica serovar Gallinarum biovar Pullorum isolates from Brazil.

Salmonella enterica subspecies enterica serovar Gallinarum biovar Pullorum (S. Pullorum) is a pathogenic bacterium that causes Pullorum disease (PD). PD is an acute systemic disease that affects young chickens, causing white diarrhea and high mortality. Although many sanitary programs have been carried out to eradicate S. Pullorum, PD outbreaks have been reported in different types of birds (layers, broilers, breeders) worldwide. This study aimed to evaluate the evolution and genetic characteristics of S. Pullorum isolated from PD in Brazil. Phylogenetic analysis of S. Pullorum genomes sequenced in this study and available genomic databases demonstrated that all isolates from Brazil are from sequence type 92 (ST92) and cluster into two lineages (III and IV). ColpVC, IncFIC(FII), and IncFII(S) were plasmid replicons frequently found in the Brazilian lineages. Two resistance genes (aac(6')-Iaa, conferring resistance to aminoglycoside, disinfecting agents, and antiseptics (mdf(A)) and tetracycline (mdf(A)) were detected frequently. Altogether, these results are important to understand the circulation of S. Pullorum and, consequently, to develop strategies to reduce losses due to PD.

Evolución y perfil genómico de aislados de Salmonella enterica serovar Gallinarum biovar Pullorum de Brasil. Salmonella enterica subespecie enterica serovar Gallinarum biovar Pullorum (S. Pullorum) es una bacteria patógena que causa la enfermedad de Pullorum (EP). La EP es una enfermedad sistémica aguda que afecta a los pollos jóvenes causando diarrea blanca y alta mortalidad. Aunque se han llevado a cabo muchos programas sanitarios para erradicar S. Pullorum, se han informado brotes de EP en diferentes tipos de aves (ponedoras, pollos de engorde, reproductoras) en todo el mundo. Este estudio tuvo como objetivo evaluar la evolución y las características genéticas de S. Pullorum aislado de EP en Brasil. El análisis filogenético de los genomas de S. Pullorum secuenciados en este estudio y las bases de datos genómicas disponibles demostraron que todos los aislamientos de Brasil son del tipo de secuencia 92 (ST92) y se agrupan en dos linajes (III y IV). ColpVC, IncFIC (FII) e IncFII(S) fueron replicones de plásmidos frecuentemente encontrados en los linajes brasileños. Dos genes de resistencia (aac(6')-Iaa, que confiere resistencia a aminoglucósidos, desinfectantes y antisépticos (mdf(A)), y tetraciclina (mdf(A)) fueron detectados con frecuencia. En conjunto, estos resultados son importantes para comprender la circulación de S. Pullorum y, en consecuencia, desarrollar estrategias para reducir las pérdidas por EP.

Epidemiological Insights into Feline Leukemia Virus Infections in an Urban Cat (Felis catus) Population from Brazil.

Feline leukemia virus (FeLV) is a retrovirus distributed worldwide in domestic cats and with different outcomes (progressive, regressive, abortive, focal). The present study reports an epidemiological survey of FeLV frequency and the evaluation of some risk factors and the two main disease outcomes (progressive and regressive) in an urban cat population from Brazil. A total of 366 cats with sociodemographic information and p27 FeLV antigen test performed were included in the study. FeLV DNA (provirus) in the blood samples of all cats was detected via real-time polymerase chain reaction (qPCR). Plasma samples from 109 FeLV-positive and FeLV-negative cats were also submitted to reverse transcription (RT-qPCR) to determine the FeLV viral load. The results demonstrated that 112 (30.6%) cats were positive through the p27 antigen and/or qPCR. A risk factor analysis demonstrated that cats without vaccination against FeLV (OR 9.9, p < 0.001), clinically ill (OR 2.9, p < 0.001), with outdoors access (OR 2.7, p < 0.001), and exhibiting apathetic behavior (OR 3.1, p < 0.001) were more likely to be infected with FeLV. FeLV-infected cats were also more likely to present with anemia (OR 13, p < 0.001) and lymphoma (OR 13.7, p = 0.001). A comparative analysis of the different detection methods in a subset of 109 animals confirmed FeLV infection in 58 cats, including 38 (65.5%) with progressive, 16 (27.6%) with regressive, and 4 (6.9%) with probably focal outcome diseases. In conclusion, this study demonstrates a high prevalence of FeLV in this urban cat population from Brazil and highlights the need to establish more effective prevention strategies (such as viral testing, vaccination programs, specific care for FeLV-positive cats) to reduce diseases associated with this virus in Brazil.

Evolution and Epidemic Spread of the Avian Infectious Bronchitis Virus (IBV) GI-23 in Brazil.

Infectious bronchitis virus (IBV) is a pathogen affecting poultry flocks worldwide. GI-23 is an IBV lineage with a rapid spread into different continents of the world, and it was reported for the first time in South American/Brazilian broiler farms last year. This study aimed to investigate the recent introduction and epidemic spread of IBV GI-23 in Brazil. Ninety-four broiler flocks infected with this lineage were evaluated from October 2021 to January 2023. IBV GI-23 was detected using real-time RT-qPCR, and the S1 gene hypervariable regions 1 and 2 (HVR1/2) were sequenced. S1 complete and HVR1/2 nucleotide sequence datasets were used to carry out phylogenetic and phylodynamic analyses. Brazilian IBV GI-23 strains clustered into two specific subclades (SA.1 and SA.2), both in tree branches with IBV GI-23 from Eastern European poultry-producing countries, suggesting two independent and recent introductions (around 2018). Viral phylodynamic analysis showed that the IBV GI-23 population increased from 2020 to 2021, remaining constant for one year and declining in 2022. S1 amino acid sequences from Brazilian IBV GI-23 presented specific and characteristic substitutions in the HVR1/2 for subclades IBV GI-23 SA.1 and SA.2. This study brings new insights into the introduction and recent epidemiology of IBV GI-23 in Brazil.

Canine parvovirus type 2 vaccines in Brazil: Viral load in commercial vaccine vials and phylogenetic analysis of the vaccine viruses.

Canine parvovirus type 2 (CPV-2) is the etiological agent of a highly contagious and frequently fatal disease in dogs. Live attenuated vaccines (LAV) are recommended to prevent and control this disease. Commercial vaccines are typically produced with CPV-2 strains adapted to cell culture and usually non-pathogenic. The present study aimed to determine the viral load of CPV-2 vaccines commercially available in Brazil and to characterize the vaccine virus by DNA analysis of its capsid gene. The results demonstrated that all vaccine strains presented high homology of the VP2 gene and they were all closely related to the original CPV-2 strains. However, vaccine strains presented several differences in comparison with field strains currently circulating in Brazil. Seventy-one vials contained viral loads ranging from 7.4E3 to 4.9E10 DNA copies/ml. Nine vials did not contain any detectable CPV-2 DNA. In conclusion, there are genetic and antigenic differences among CPV-2 vaccines and field strains. Additionally, some vaccines have been commercialized with low titers of CPV-2. It is important to improve the quality of the vaccines to prevent or reduce the spread of CPV-2 in Brazil.

Phylogenetic classification of feline immunodeficiency virus

Background: The feline immunodeficiency virus (FIV) is responsible for a retroviral disease that affects domestic and wild cats worldwide, causing Feline Acquired Immunodeficiency Syndrome (FAIDS). FIV is a lentivirus from the family Retroviridae and its genome has 3 main structural genes: gag, pol and env. Phylogenetic studies have classified FIV into 7 subtypes according to the diversity among strains from the World, mainly in the env gene. Epidemiological analyses have demonstrated the high predominance of FIV-A and FIV-B. This in silico study aimed to perform a phylogenetic analysis to study FIV diversity worldwide. Materials, Methods & Results: A total of 60 whole genome sequences (WGS) and 122 FIV env gene sequences were included in 2 datasets, which were aligned using MAFFT version 7. Recombination among genomes and/or env genes was analyzed with RDP5 software. Phylogenetic analyses with both datasets were performed, after removing the recombinant sequences, by the W-IQ-TREE and constructed and edited by the FigTree. A total of 12 recombination events involving 19 WGS were detected. In addition, 27 recombination events involving 49 sequences were observed in the env gene. A high rate of recombinants was observed inter-subtypes (A/B and B/D) and intra-subtypes (A/A). All recombinants were removed from the subsequent phylogenetic analyses. Phylogenies demonstrated 6 distinct main clades, 5 from domestic cats (A, B, C, E, U) and 1 from wild cat sequences (W) in the WGS, as well as in the specific env gene analyses. Most clustered with subtype B sequences. In the WGS analysis, clade B had a prevalence of 65.9% Brazilian sequences (27/41) and 2.4% Japanese sequences (1/41). In the env gene analyses, clade B showed a prevalence of 43.8% of Brazilian sequences (32/73) and 20.5% of USA sequences (15/73). The results of both analyses also confirm the FIV-wide geographical distribution around the world. In the phylogenetic analyses carried out with WGS, sequences from China (1/41; 2.4%), Colombia (1/41; 2.4%) and the USA (1/41; 2.4%) were identified in clade A; sequence from Canada in clade C (1/41; 2.4%); sequence from Botswana belonged to clade E (1/41; 2.4%); sequences from Brazil clustered into clade U (2/41; 5% - data not yet published); and sequences belonging to the clade W were from Canada (1/41; 2.4%) and the USA (5/41; 12.3%). Specific env gene phylogenetic analyses showed sequences from Colombia (1/73; 1.4%), France (2/73; 2.7%), the Netherlands (3/73; 4.1%), Switzerland (2/73; 2.7%), USA (6/73; 8.3%), belonging to clade A; sequence from Canada belonging to clade C (1/73; 1.4%); sequences from Brazil belonging to clade U (2/73; 5% - data not yet published); and sequences belonging to clade W from the USA (6/73; 8.3%). Discussion: The results presented here demonstrate that FIV has a rapid viral evolution due to recombination and mutation events, more specifically in the env gene, which is highly variable. Currently, this retrovirus is classified into 7 subtypes (A, B, C, D, E, F and U-NZenv) according to their high genomic diversity. It also highlighted the importance of in silico sequence and phylogeny studies to demonstrate evolutionary processes. This was the first study to address the WGS FIV diversity with a phylogenetic approach.

Emergence, Dissemination and Antimicrobial Resistance of the Main Poultry-Associated Salmonella Serovars in Brazil.

Salmonella infects poultry, and it is also a human foodborne pathogen. This bacterial genus is classified into several serovars/lineages, some of them showing high antimicrobial resistance (AMR). The ease of Salmonella transmission in farms, slaughterhouses, and eggs industries has made controlling it a real challenge in the poultry-production chains. This review describes the emergence, dissemination, and AMR of the main Salmonella serovars and lineages detected in Brazilian poultry. It is reported that few serovars emerged and have been more widely disseminated in breeders, broilers, and layers in the last 70 years. Salmonella Gallinarum was the first to spread on the farms, remaining as a concerning poultry pathogen. Salmonella Typhimurium and Enteritidis were also largely detected in poultry and foods (eggs, chicken, turkey), being associated with several human foodborne outbreaks. Salmonella Heidelberg and Minnesota have been more widely spread in recent years, resulting in frequent chicken/turkey meat contamination. A few more serovars (Infantis, Newport, Hadar, Senftenberg, Schwarzengrund, and Mbandaka, among others) were also detected, but less frequently and usually in specific poultry-production regions. AMR has been identified in most isolates, highlighting multi-drug resistance in specific poultry lineages from the serovars Typhimurium, Heidelberg, and Minnesota. Epidemiological studies are necessary to trace and control this pathogen in Brazilian commercial poultry production chains.

Molecular phylogenetic assessment of the canine parvovirus 2 worldwide and analysis of the genetic diversity and temporal spreading in Brazil.

Canine parvovirus type 2 (CPV-2) is a relevant pathogen for dogs and causes a severe disease in carnivore species. CPV-2 reached pandemic proportions after the 1970s with the worldwide dissemination, generating antigenic and genetic variants (CPV-2a, CPV-2b, and CPV-2c) with different pathobiology in comparison with the original type CPV-2. The present study aimed to assess the current global CPV-2 molecular phylogeny and to analyze genetic diversity and temporal spreading of variants from Brazil. A total of 284 CPV-2 whole-genome sequences (WGS) and 684 VP2 complete genes (including 23 obtained in the present study) were compared to analyze phylogenetic relationships. Bayesian coalescent analysis estimated the time to the most recent common ancestor (tMRCA) and the population dynamics of the different CPV-2 lineages in the last decades. The WGS phylogenetic tree demonstrated two main clades disseminated worldwide today. The VP2 gene tree showed a total of four well-defined clades distributed in different geographic regions, including one with CPV-2 sequences exclusive from Brazil. These clades do not have a relationship with the previous classification into CPV-2a, CPV-2b, and CPV-2c, despite some having a predominance of one or more antigenic types. Temporal analysis demonstrated that the main CPV-2 clades evolved within a few years (from the 1980s to 1990s) in North America and they spread worldwide afterwards. Population dynamics analysis demonstrated that CPV-2 presented a major dissemination increase at the end of the 1980s / beginning of the 1990s followed by a period of stability and a second minor increase from 2000 to 2004.

Temporal spread and evolution of SARS-CoV-2 in the second pandemic wave in Brazil.

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) pandemic spread rapidly and this scenario is concerning in South America, mainly in Brazil that presented more than 21 million coronavirus disease 2019 cases and 590 000 deaths. The recent emergence of novel lineages carrying several mutations in the spike protein has raised additional public health concerns worldwide. The present study describes the temporal spreading and evolution of SARS-CoV2 in the beginning of the second pandemic wave in Brazil, highlighting the fast dissemination of the two major concerning variants (P.1 and P.2). A total of 2507 SARS-CoV-2 whole-genome sequences (WGSs) with available information from the country (Brazil) and sampling date (July 2020-February 2021), were obtained and the frequencies of the lineages were evaluated in the period of the growing second pandemic wave. The results demonstrated the increasing prevalence of P.1 and P.2 lineages in the period evaluated. P.2 lineage was first detected in the middle of 2020, but a high increase occurred only in the last trimester of this same year and the spreading to all Brazilian regions. P.1 lineage emerged even later, first in the North region in December 2020 and really fast dissemination to all other Brazilian regions in January and February 2021. All SARS-CoV-2 WGSs of P.1 and P.2 were further separately evaluated with a Bayesian approach. The rates of nucleotide and amino acid substitutions were statistically higher in P.1 than P.2 (p < 0.01). The phylodynamic analysis demonstrated that P.2 gradually spread in all the country from September 2020 to January 2021, while P.1 disseminated even faster from December 2020 to February 2021. Skyline plots of both lineages demonstrated a slight rise in the spreading for P.2 and exponential growth for P.1. In conclusion, these data demonstrated that the P.1 (recently renamed as Gamma) and P.2 lineages have predominated in the second pandemic wave due to the very high spreading across all geographic regions in Brazil at the end of 2020 and beginning of 2021.

Recent Evolution and Genomic Profile of Salmonella enterica Serovar Heidelberg Isolates from Poultry Flocks in Brazil.

Salmonella enterica serovar Heidelberg is isolated from poultry-producing regions around the world. In Brazil, S. Heidelberg has been frequently detected in poultry flocks, slaughterhouses, and chicken meat. The goal of the present study was to assess the population structure, recent temporal evolution, and some important genetic characteristics of S. Heidelberg isolated from Brazilian poultry farms. Phylogenetic analysis of 68 S. Heidelberg genomes sequenced here and additional whole-genome data from NCBI demonstrated that all isolates from the Brazilian poultry production chain clustered into a monophyletic group, here called S. Heidelberg Brazilian poultry lineage (SH-BPL). Bayesian analysis defined the time of the most recent common ancestor (tMRCA) as 2004, and the overall population size (Ne) was constant until 2008, when an ∼10-fold Ne increase was observed until circa 2013. SH-BPL presented at least two plasmids with replicons ColpVC (n = 68; 100%), IncX1 (n = 66; 97%), IncA/C2 (n = 65; 95.5%), ColRNAI (n = 43; 63.2%), IncI1 (n = 32; 47%), ColMG828, Col156, IncHI2A, IncHI2, IncQ1, IncX4, IncY, and TrfA (each with n < 4; <4% each). Antibiotic resistance genes were found, with high frequencies of fosA7 (n = 68; 100%), mdf(A) (n = 68; 100%), tet(34) (n = 68; 100%), sul2 (n = 64; 94.1%), and blaCMY-2 (n = 56; 82.3%), along with an overall multidrug resistance (MDR) profile. Ten Salmonella pathogenicity islands (SPI1 to SPI5, SPI9, and SPI11 to SPI14) and 139 virulence genes were also detected. The SH-BPL profile was like those of other previous S. Heidelberg isolates from poultry around the world in the 1990s. In conclusion, the present study demonstrates the recent introduction (2004) and high level of dissemination of an MDR S. Heidelberg lineage in Brazilian poultry operations. IMPORTANCES. Heidelberg is the most frequent serovar in several broiler farms from the main Brazilian poultry-producing regions. Therefore, avian-source foods (mainly chicken carcasses) commercialized in the country and exported to other continents are contaminated with this foodborne pathogen, generating several national and international economic losses. In addition, isolates of this serovar are usually resistant to antibiotics and can cause human invasive and septicemic infection, representing a public health concern. This study demonstrates the use of whole-genome sequencing (WGS) to obtain epidemiological information for one S. Heidelberg lineage highly spread among Brazilian poultry farms. This information will help to define biosecurity measures to control this important Salmonella serovar in Brazilian and worldwide poultry operations.

Retrospective whole-genome comparison of Salmonella enterica serovar Enteritidis from foodborne outbreaks in Southern Brazil.

Salmonella enterica serovar Enteritidis is frequently isolated from animal-source foods associated with human salmonellosis outbreaks. This serovar was spread to animal (mainly poultry) farms worldwide in the 1980s, and it is still detected in foods produced in many countries, including Brazil. The present study reports a retrospective genome-wide comparison of S. Enteritidis from foodborne outbreaks in Southern Brazil in the last two decades. Fifty-two S. Enteritidis isolates were obtained from foodborne outbreaks occurring in different cities of the Brazilian southernmost State, Rio Grande do Sul (RS), from 2003 to 2015. Whole-genome sequences (WGS) from these isolates were obtained and comparatively analyzed with 65 additional genomes from NCBI. Phylogenetic and Bayesian analyses were performed to study temporal evolution. Genes related to antibiotic resistance and virulence were also evaluated. The results demonstrated that all S. Enteritidis isolates from Southern Brazil clustered in the global epidemic clade disseminated worldwide originally in the 1980s. Temporal analysis demonstrated that all Brazilian isolates had a tMRCA (time to most recent common ancestor) in 1986 with an effective population size (Ne) increase soon after until 1992, then becoming constant up to now. In Southern Brazil, there was a significant decrease in the spreading of S. Enteritidis in the last decade. In addition, three antibiotic resistance genes were detected in all isolates: aac(6')-Iaa, mdfA, and tet(34). These results demonstrate the high frequency of one only specific S. Enteritidis lineage (global epidemic clade) in foodborne outbreaks from Southern Brazil in the last two decades.

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

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.

Salmonella serotype assignment by sequencing analysis of intergenic regions of ribosomal RNA operons.

Salmonella laboratorial detection is usually carried out by bacteriological culture and serological methods. Salmonella isolates are then classified into more than 2,650 serotypes with different somatic (O) and flagellar (H) antigenic combinations. More recently, DNA analysis methods were developed and applied for the identification of Salmonella serotypes, including intergenic spacer regions (ISRs) that separates DNA-encoding ribosomal subunits (rRNA gene) in Salmonella genomes. The present study aimed to evaluate the nucleotide diversity of the ISRs in 2 rRNA operons (rrnB and rrnH) for the assignment of Salmonella serotypes. A total of 63 Salmonella isolates (bacterial cultures) from 21 serotypes were obtained in the period of 2014 to 2017. DNA was extracted, and PCRs were used to detect the genus Salmonella and 4 important serotypes: Enteritidis, Gallinarum, Heidelberg, and Typhimurium. ISRs of the operons rrnB and rrnH were amplified by PCR and then sequenced. All sequence results were submitted to BLASTn search and were aligned in comparison to 72 Salmonella reference nucleotide sequences. The results demonstrated that 60 (95.2%) samples returned a sequence of the same serotype (determined by the traditional serological procedure) after searching in BLASTn and/or in the alignment with the reference sequences for both operons (rrnB and rrnH). These PCR-sequencing procedures had a general agreement index of 0.792 based on the Kappa analysis, 98.7% sensitivity value, 100% specificity, and 98.4% accuracy. Three different phylogenetic trees were generated from the alignments with the sequences of rrnH, rrnB, and rrnH plus rrnB and isolates clustered in specific branches according to the serotypes.

Molecular and phylogenetic analyses of Salmonella Gallinarum trace the origin and diversification of recent outbreaks of fowl typhoid in poultry farms.

Fowl typhoid (FT) and pullorum disease (PD) are two important poultry infections caused by Salmonella enterica subsp. enterica serotype Gallinarum (S. Gallinarum). S. Gallinarum strains are adapted to birds and classified into biovars Gallinarum (bvGA) and Pullorum (bvPU) as they are the causative agent of FT and PD, respectively. In Brazil, FT/PD outbreaks have been reported along the last 50 years, but there was a recent increase of FT field reports with the suspicion it could be due to virulence reversion of the attenuated live vaccine SG9R. In this study, we applied molecular biology assays and phylogenetic methods to detect and investigate S. Gallinarum isolates from commercial poultry flocks in order to understand the evolutionary history and origin of the recent FT outbreaks in Brazil. S. Gallinarum isolates were obtained from thirteen different poultry flocks with clinical signs of FT/PD from 2013 to 2015. These isolates were serotyped, tested with three specific PCR (for the detection of bvGA, bvPU and live vaccine strain SG9R) and submitted to sequencing of a variable genome region (ISR analysis). The complete genome of one bvGA strain (BR_RS12) was also compared to other S. Gallinarum complete genomes (including other two Brazilian ones: bvGA 287/91 and bvPU FCVA198). PCR detected all thirteen isolates as S. Gallinarum (eight bvGA and five bvPU), none positive for SG9R strain. ISR analysis revealed that all eight bvGA isolates showed exactly the same nucleotide sequences with 100% similarity to reference strains, while two patterns were observed for bvPU. Genome phylogeny demonstrated distinct clades for bvGA and bvPU, with the bvGA clade showing a clear subdivision including three genomes: SG9R vaccine, the respective SG9 parent strain and one SG9R revertant field isolate (MB4523). The evolutionary rate of the total S. Gallinarum genome was calculated at 6.15×10-7 substitutions/site/year, with 2.8 observed substitutions per year per genome (1 SNP per 4292 bases). Phylodynamics analysis estimated that at least two introductions of S. Gallinarum bvGA happened in Brazil, the first in 1885 and the second in 1950. The Brazilian bvGA genomes 287/91 and BR_RS12 analyzed here were related to the early and the late introductions, respectively. In conclusion, these results indicate the occurrence of S. Gallinarum strains associated with FT outbreaks that have been circulating for more than 50 years in Brazil and are not originated from virulence reversion of the SG9R vaccine.