Phylogenetic and phylogeographic mapping of the avian coronavirus spike protein-encoding gene in wild and synanthropic birds

Avian infectious bronchitis virus (AIBV) review by continent

Infectious bronchitis virus (IBV) is a positive-sense, single-stranded, enveloped RNA virus responsible for substantial economic losses to the poultry industry worldwide by causing a highly contagious respiratory disease. The virus can spread quickly through contact, contaminated equipment, aerosols, and personal-to-person contact. We highlight the prevalence and geographic distribution of all nine genotypes, as well as the relevant symptoms and economic impact, by extensively analyzing the current literature. Moreover, phylogenetic analysis was performed using Molecular Evolutionary Genetics Analysis (MEGA-6), which provided insights into the global molecular diversity and evolution of IBV strains. This review highlights that IBV genotype I (GI) is prevalent worldwide because sporadic cases have been found on many continents. Conversely, GII was identified as a European strain that subsequently dispersed throughout Europe and South America. GIII and GV are predominant in Australia, with very few reports from Asia. GIV, GVIII, and GIX originate from North America. GIV was found to circulate in Asia, and GVII was identified in Europe and China. Geographically, the GVI-1 lineage is thought to be restricted to Asia. This review highlights that IBV still often arises in commercial chicken flocks despite immunization and biosecurity measures because of the ongoing introduction of novel IBV variants and inadequate cross-protection provided by the presently available vaccines. Consequently, IB consistently jeopardizes the ability of the poultry industry to grow and prosper. Identifying these domains will aid in discerning the pathogenicity and prevalence of IBV genotypes, potentially enhancing disease prevention and management tactics.

Molecular characterization of the S1 gene in GI-17 and GI-13 type isolates of avian infectious bronchitis virus (IBV) in Costa Rica, from 2016 to 2019

Avian infectious bronchitis is one of the most important respiratory diseases affecting poultry production worldwide. The etiological agent of this disease is the avian infectious bronchitis virus (IBV). We analyzed 14 isolates of IBV obtained from poultry farms in Costa Rica, from 2016 through 2019. We sequenced the S1 region of the genome and the sequences obtained were submitted to GenBank. Phylogenetic analyses showed that the isolates obtained during 2016-2017 belong to the GI-17 lineage and are related to the Georgia 13-type Ga-13/14255/14 and CK/CR/1160/16 variants, with a 96.90-100% nucleotide sequence identity and a 92.25-100% amino acid sequence identity. The main differences were detected in the RBD and HVR-3 regions, where a series of mutations eliminate an N-glycosylation site in 10 out of 11 isolates. The isolates obtained during 2018-2019 belong to the GI-13 lineage and are closely related to the 4/91 vaccine variant, with over 98% sequence identity at the nucleotide and amino acids levels. Variations were detected in the RBD and HVR regions, with a possible N-glycosylation site detected in isolate CK/CR/0632/19. These results indicate that a GA13-like pathogenic variant circulated during the 2016-2017 period and that the 4/91 variant was detected after the introduction of the vaccine. The variations shown in both the GA13-like and 4/91 isolates examined, reveal the need for continuous surveillance of IBV in Costa Rica, to detect new variants that may be introduced to the country or develop during outbreaks. This information is highly relevant for vaccination planning and disease management programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13337-022-00762-2.

Biodiversity and public health interface

Abstract Alongside modernity, the human activity has been a key factor in global environmental risks, with worldwide anthropic modification being the cause of the emergence of diseases for wild and livestock animals, and even humans. In special, the increase in the spatial distribution and in the incidence of some emerging infectious diseases (EID) are directly associated to deforestation and global climate changes. Moreover, the arise of new EID agents, such as the SARS-COV-2 have been reported for the last 30 years. On the other hand, biodiversity has been shown to be a key indicator for ecosystem health, and to pose a role to increase the promotion of human public health. In neotropical regions, and in special, in Brazil, several infectious diseases have been demonstrated to be directly affected for the biodiversity loss, such as malaria, hantavirus pulmonary syndrome, yellow fever, urban arboviruses, spotted fever, amongst other. To better understand the ecosystem capacity of regulation of infectious diseases, FAPESP BIOTA program have supported researchers and research projects to increase knowledge about Brazilian biodiversity and the ecosystems, such as diversity of bird bioagents, venomous animals biodiversity, diversity of mosquitos species in forest patches inside urban areas, propagation of the yellow fever virus over fragmented forest territories, loss of ecological corridors and occurrence of spotted fever and malaria, amongst others. It is noteworthy that FAPESP BIOTA is a successful program and must be expanded as an important tool for present and future public health promotion.

Diversity of Coronaviruses in Wild Representatives of the Aves Class in Poland

The revealed prevalence of coronaviruses in wild bird populations in Poland was 4.15% and the main reservoirs were birds from orders Anseriformes and Charadriiformes, with a prevalence of 3.51% and 5.59%, respectively. Gammacoronaviruses were detected more often than deltacoronaviruses, with detection rates of 3.5% and 0.7%, respectively. Gammacoronaviruses were detected in birds belonging to six orders, including Anseriformes, Charadriiformes, Columbiformes, Galliformes, Gruiformes, and Passeriformes, indicating a relatively wide host range. Interestingly, this was the only coronavirus detected in Anseriformes (3.51%), while in Charadriiformes, the prevalence was 3.1%. The identified gammacoronaviruses belonged to the Igacovirus and Brangacovirus subgeneras. Most of these were igacoviruses and formed a common phylogenetic group with a Duck Coronavirus 2714 and two with an Avian Coronavirus/Avian Coronavirus9203, while the viruses from the pigeons formed a distinct “pigeon-like” group, not yet officially represented. The presence of deltacoronaviruses was detected in birds belonging to three orders, Charadriiformes, Galliformes, and Suliformes indicating a narrower host range. Most identified deltacoronaviruses belonged to the Buldecovirus subgenus, while only one belonged to Herdecovirus. Interestingly, the majority of buldecoviruses were identified in gulls, and they formed a distinct phylogenetic lineage not represented by any officially ratified virus species. Another separate group of buldecoviruses, also not represented by the official species, was formed by a virus identified in a common snipe. Only one identified buldecovirus (from common pheasant) formed a group with the ratified species Coronavirus HKU15. The results obtained indicate the high diversity of detected coronaviruses, and thus also the need to update their taxonomy (establishing new representative virus species). The serological studies performed revealed antibodies against an infectious bronchitis virus in the sera of white storks and mallards.

Development of a Dissemination Platform for Spatiotemporal and Phylogenetic Analysis of Avian Infectious Bronchitis Virus

Infecting large portions of the global poultry populations, the avian infectious bronchitis virus (IBV) remains a major economic burden in North America. With more than 30 serotypes globally distributed, Arkansas, Connecticut, Delaware, Georgia, and Massachusetts are among the most predominant serotypes in the United States. Even though vaccination is widely used, the high mutation rate exhibited by IBV is continuously triggering the emergence of new viral strains and hindering control and prevention measures. For that reason, targeted strategies based on constantly updated information on the IBV circulation are necessary. Here, we sampled IBV-infected farms from one US state and collected and analyzed 65 genetic sequences coming from three different lineages along with the immunization information of each sampled farm. Phylodynamic analyses showed that IBV dispersal velocity was 12.3 km/year. The majority of IBV infections appeared to have derived from the introduction of the Arkansas DPI serotype, and the Arkansas DPI and Georgia 13 were the predominant serotypes. When analyzed against IBV sequences collected across the United States and deposited in the GenBank database, the most likely viral origin of our sequences was from the states of Alabama, Georgia, and Delaware. Information about vaccination showed that the MILDVAC-MASS+ARK vaccine was applied on 26% of the farms. Using a publicly accessible open-source tool for real-time interactive tracking of pathogen spread and evolution, we analyzed the spatiotemporal spread of IBV and developed an online reporting dashboard. Overall, our work demonstrates how the combination of genetic and spatial information could be used to track the spread and evolution of poultry diseases, providing timely information to the industry. Our results could allow producers and veterinarians to monitor in near-real time the current IBV strain circulating, making it more informative, for example, in vaccination-related decisions.

Molecular characterization of an avian GA13-like infectious bronchitis virus full-length genome from Costa Rica

We describe the first whole-genome sequence of a GA13-like isolate of avian infectious bronchitis virus CK/CR/1160/16 (MN757859), obtained in 2016 in the province of Alajuela, Costa Rica. This virus caused an outbreak with great economic impact to the local poultry industry. The genome sequence is 27 696 bp in length, with the following genome organization 5'-UTR-Pol-S-3a-3b-E-4b-4c-M-5a-5b-N-6b-3'-UTR. The complete genome sequence has the highest sequence identity (94.03%) with DMV/1639/GA9977/2019 (MK878536) from Georgia, USA, and the lowest identity (86.03%) with ck/CH/LHLJ/08-6 (KX252788), from China. Analysis of the S1 subunit indicates that the Costa Rican isolate belongs to genotype I, lineage 17 (GI-17) and displays 96.89% identity with the S1 subunit of Ga-13/14255/14 (KM087780) (USA). Possible recombination events in genes S, E, M, 4b y 4c were detected, with Massachusetts, Connecticut, Arkansas and MA5 as potential parental types. This study highlights the importance of the epidemiological and molecular surveillance of avian infectious bronchitis.

Live Attenuated Infectious Bronchitis Virus Vaccines in Poultry: Modifying Local Viral Populations Dynamics

Simple Summary

Infectious bronchitis (IB) is one of the more prevalent diseases in poultry, and it is caused by a virus belonging to the Coronaviridae family, the infectious bronchitis virus (IBV), a Gammacoronavirus which is related to the Betacoronavirus SARSCov-2 causing COVID-19 in humans. IB is mainly controlled by biosecurity and vaccines, although, it is a very challenging issue because the viral populations are constantly evolving by several factors. One of these factors is the same vaccines used for IB control, this could explain by recombination, reversion to virulence, or by favoring virus serotype selection. Thus, a human role in the change of viral populations can be identified by the IBV vaccine usage, this must be considered to achieve effective IB control.

Abstract

Infectious bronchitis virus (IBV) remains one of the most important diseases impacting poultry today. Its high adaptive capacity, attributable to the high mutation rate associated with its ssRNA(+), is one of its more important features. While biosecurity procedures and barriers have been shown to be preponderant factors in minimizing the impact of infectious bronchitis (IB), the environment and procedures associated with intensive poultry systems greatly influence the viral population dynamics. High-density poultry flocks facilitate recombination between different viruses, and even with live attenuated vaccines, which can change the dominant circulating field strains. Furthermore, the remaining issue of reversion to virulence gives rise to significant problems when vaccinal strains are introduced in places where their pathogenic variants have not been reported. Under specific conditions, live attenuated vaccines could also change the frequency of circulating viruses and enable replacement between different field strains. In summary, under a comprehensive approach, while vaccination is one of the most essential tools for controlling IB, the veterinarians, farmers, and official services role in its usage is central to minimizing alteration in a malleable viral population. Otherwise, vaccination is ultimately counterproductive.

Genetic, Antigenic, and Pathogenic Characteristics of Infectious Bronchitis Virus GI-7/TW-II in China

Nine infectious bronchitis virus (IBV) strains belonging to the GI-7 lineage were isolated between 2009 and 2017 in China. Phylogenetic analysis and comparisons of full-length sequences of the S1 gene suggested that the GI-7 lineage should be further classified as Taiwan (TW)-I and TW-II sublineages, which correspond to the previous TW-I and TW-II genotypes. The nine IBV strains were clustered in the TW-II sublineage. Further investigation revealed that viruses in the TW-I and TW-II were not only genetically but also antigenically different. Moreover, the TW-II sublineage contained various clades and recombinants. A recombinant was found to originate from recombination events between field strains (TW-II ck/CH/LJL/090608- and GI-19 ck/ CH/LDL/091022-like viruses) in which the recombination in the S1 subunit coding sequences had led to changes in antigenicity of the viruses. A more in-depth investigation demonstrated that TW-II viruses appear to have undergone a significant evolution following introduction in mainland China, which resulted in the viruses diverging into different clades. The viruses between the different clades in TW-II sublineage exhibited a significant change in genetic and antigenic characteristics. In addition, the five TW-II viruses selected on the basis of the results of S1 nucleotide sequence phylogenetic trees showed different pathogenicity to specific-pathogen-free chickens, although they could induce nephritis in the infected chickens and thus were identified as nephropathogenic strains.

Wild birds as reservoirs for diverse and abundant gamma- and deltacoronaviruses

Wild birds interconnect all parts of the globe through annual cycles of migration with little respect for country or continental borders. Although wild birds are reservoir hosts for a high diversity of gamma- and deltacoronaviruses, we have little understanding of the ecology or evolution of any of these viruses. In this review, we use genome sequence and ecological data to disentangle the evolution of coronaviruses in wild birds. Specifically, we explore host range at the levels of viral genus and species, and reveal the multi-host nature of many viral species, albeit with biases to certain types of avian host. We conclude that it is currently challenging to infer viral ecology due to major sampling and technical limitations, and suggest that improved assay performance across the breadth of gamma- and deltacoronaviruses, assay standardization, as well as better sequencing approaches, will improve both the repeatability and interpretation of results. Finally, we discuss cross-species virus transmission across both the wild bird - poultry interface as well as from birds to mammals. Clarifying the ecology and diversity in the wild bird reservoir has important ramifications for our ability to respond to the likely future emergence of coronaviruses in socioeconomically important animal species or human populations.

Molecular and Antigenic Characterization of GI-13 and GI-16 Avian Infectious Bronchitis Virus Isolated in Chile from 2009 to 2017 Regarding 4/91 Vaccine Introduction

Simple Summary

The high adaptation and recombination abilities of infectious bronchitis virus (IB) have been proven. This study aims to verify the genetic and antigenic variation of eight field IB strains regarding the 4/91 strain vaccination in Chilean chickens. Phylogenetic, serologic and challenge studies were carried out to accomplish this goal. The genetic analyses indicate that all the viruses isolated prior to the 4/91 introduction belong to the genetic group GI-16 (three isolates from 2009). On the other hand, just one of the viruses isolated after the 4/91 strain vaccine introduction in Chile (in 2015) showed relationship with GI-16 lineage. The remaining four viruses (from 2017) belong to GI-13, the group where the strain 4/91 has been previously classified. Three viruses were chosen to perform antigenic and protective studies. Antigenically, the high relationship between the 4/91 vaccine with the isolate from 2017 is remarkable and could not be observed with isolates from 2009 and 2015. The 4/91 vaccine also showed better protection against the isolate from 2017 than isolates from 2009 and 2015. These results suggest that the introduction of the 4/91 vaccine in Chile could imply a change in some viruses, showing its ability to interact with field viruses, so it is important to monitor the circulating viruses and include these results in future governmental decisions.

Abstract

The introduction of the 4/91 vaccine against infectious bronchitis in Chile, a lineage not described until that time in the country, led to looking for changes induced by this action. This study considers eight isolates obtained from 2009, 2015 and 2017 and uses a maximum likelihood approach to classify the field isolates. Three isolates were selected to analyze antigenic relationships through a virus neutralization test and to perform protection tests measured trough an RT-qPCR. The isolates from 2009 and 2015 showed a relationship with GI-16 while those from 2017 were related to GI-13. Though the field isolates were classified in two different phylogenetic lineages, all of them showed only minor variations in subtype. The 13885R-17 isolate from 2017 exhibited high antigenic relatedness to the 4/91 vaccine. As expected, 4/91 and Massachusetts vaccines were not antigenically related. Vaccinated birds with the 4/91 vaccine showed less tracheal virus replication for the 13885R-17 from 2017 challenge than for the 12101SP-09 from 2009 and 13347SP-15 from 2015 isolates. The results indicated genetic and antigenic diversity in the most recent infectious bronchitis virus (IBV) isolates in Chile. Moreover, the 4/91 vaccine would be involved in the generation of some current field viruses, which must be considered in vaccination programs and public policies.

Guinea Fowl Coronavirus Diversity Has Phenotypic Consequences for Glycan and Tissue Binding

Avian coronaviruses cause major global problems in the poultry industry. As causative agents of huge economic losses, the detection and understanding of the molecular determinants of viral tropism are of ultimate importance. Here, we set out to study those parameters and obtained in-depth insight into the virus-host interactions of guinea fowl coronavirus (GfCoV). Our data indicate that diversity in GfCoV viral attachment proteins results in differences in degrees of affinity for glycan receptors, as well as altered avidity for intestinal tract tissues, which might have consequences for GfCoV tissue tropism and pathogenesis in guinea fowls.

Guinea fowl coronavirus (GfCoV) causes fulminating enteritis that can result in a daily death rate of 20% in guinea fowl flocks. Here, we studied GfCoV diversity and evaluated its phenotypic consequences. Over the period of 2014 to 2016, affected guinea fowl flocks were sampled in France, and avian coronavirus presence was confirmed by PCR on intestinal content and immunohistochemistry of intestinal tissue. Sequencing revealed 89% amino acid identity between the viral attachment protein S1 of GfCoV/2014 and that of the previously identified GfCoV/2011. To study the receptor interactions as a determinant for tropism and pathogenicity, recombinant S1 proteins were produced and analyzed by glycan and tissue arrays. Glycan array analysis revealed that, in addition to the previously elucidated biantennary di-N-acetyllactosamine (diLacNAc) receptor, viral attachment S1 proteins from GfCoV/2014 and GfCoV/2011 can bind to glycans capped with alpha-2,6-linked sialic acids. Interestingly, recombinant GfCoV/2014 S1 has an increased affinity for these glycans compared to that of GfCoV/2011 S1, which was in agreement with the increased avidity of GfCoV/2014 S1 for gastrointestinal tract tissues. Enzymatic removal of receptors from tissues before application of spike proteins confirmed the specificity of S1 tissue binding. Overall, we demonstrate that diversity in GfCoV S1 proteins results in differences in glycan and tissue binding properties.

IMPORTANCE Avian coronaviruses cause major global problems in the poultry industry. As causative agents of huge economic losses, the detection and understanding of the molecular determinants of viral tropism are of ultimate importance. Here, we set out to study those parameters and obtained in-depth insight into the virus-host interactions of guinea fowl coronavirus (GfCoV). Our data indicate that diversity in GfCoV viral attachment proteins results in differences in degrees of affinity for glycan receptors, as well as altered avidity for intestinal tract tissues, which might have consequences for GfCoV tissue tropism and pathogenesis in guinea fowls.

Detection and characterisation of coronaviruses in migratory and non-migratory Australian wild birds

We evaluated the presence of coronaviruses by PCR in 918 Australian wild bird samples collected during 2016–17. Coronaviruses were detected in 141 samples (15.3%) from species of ducks, shorebirds and herons and from multiple sampling locations. Sequencing of selected positive samples found mainly gammacoronaviruses, but also some deltacoronaviruses. The detection rate of coronaviruses was improved by using multiple PCR assays, as no single assay could detect all coronavirus positive samples. Sequencing of the relatively conserved Orf1 PCR amplicons found that Australian duck gammacoronaviruses were similar to duck gammacoronaviruses around the world. Some sequenced shorebird gammacoronaviruses belonged to Charadriiformes lineages, but others were more closely related to duck gammacoronaviruses. Australian duck and heron deltacoronaviruses belonged to lineages with other duck and heron deltacoronaviruses, but were almost 20% different in nucleotide sequence to other deltacoronavirus sequences available. Deltacoronavirus sequences from shorebirds formed a lineage with a deltacoronavirus from a ruddy turnstone detected in the United States. Given that Australian duck gammacoronaviruses are highly similar to those found in other regions, and Australian ducks rarely come into contact with migratory Palearctic duck species, we hypothesise that migratory shorebirds are the important vector for moving wild bird coronaviruses into and out of Australia.

Avian viral surveillance in Victoria, Australia, and detection of two novel avian herpesviruses

Viruses in avian hosts can pose threats to avian health and some have zoonotic potential. Hospitals that provide veterinary care for avian patients may serve as a site of exposure of other birds and human staff in the facility to these viruses. They can also provide a useful location to collect samples from avian patients in order to examine the viruses present in wild birds. This study aimed to investigate viruses of biosecurity and/or zoonotic significance in Australian birds by screening samples collected from 409 birds presented to the Australian Wildlife Health Centre at Zoos Victoria’s Healesville Sanctuary for veterinary care between December 2014 and December 2015. Samples were tested for avian influenza viruses, herpesviruses, paramyxoviruses and coronaviruses, using genus- or family-wide polymerase chain reaction methods coupled with sequencing and phylogenetic analyses for detection and identification of both known and novel viruses. A very low prevalence of viruses was detected. Columbid alphaherpesvirus 1 was detected from a powerful owl (Ninox strenua) with inclusion body hepatitis, and an avian paramyxovirus most similar to Avian avulavirus 5 was detected from a musk lorikeet (Glossopsitta concinna). Two distinct novel avian alphaherpesviruses were detected in samples from a sulphur-crested cockatoo (Cacatua galerita) and a tawny frogmouth (Podargus strigoides). Avian influenza viruses and avian coronaviruses were not detected. The clinical significance of the newly detected viruses remains undetermined. Further studies are needed to assess the host specificity, epidemiology, pathogenicity and host-pathogen relationships of these novel viruses. Further genome characterization is also indicated, and would be required before these viruses can be formally classified taxonomically. The detection of these viruses contributes to our knowledge on avian virodiversity. The low level of avian virus detection, and the absence of any viruses with zoonotic potential, suggests low risk to biosecurity and human health.

Phylodynamic analysis and molecular diversity of the avian infectious bronchitis virus of chickens in Brazil

Avian infectious bronchitis virus (IBV) is the etiological agent of a highly contagious disease, which results in severe economic losses to the poultry industry. The spike protein (S1 subunit) is responsible for the molecular diversity of the virus and many sero/genotypes are described around the world. Recently a new standardized classification of the IBV molecular diversity was conducted, based on phylogenetic analysis of the S1 gene sequences sampled worldwide. Brazil is one of the biggest poultry producers in the world and the present study aimed to review the molecular diversity and reconstruct the evolutionary history of IBV in the country. All IBV S1 gene sequences, with local and year of collection information available on GenBank, were retrieved. Phylogenetic analyses were carried out based on a maximum likelihood method for the classification of genotypes occurring in Brazil, according to the new classification. Bayesian phylogenetic analyses were performed with the Brazilian clade and related international sequences to determine the evolutionary history of IBV in Brazil. A total of 143 Brazilian sequences were classified as GI-11 and 46 as GI-1 (Mass). Within the GI-11 clade, we have identified a potential recombinant strain circulating in Brazil. Phylodynamic analysis demonstrated that IBV GI-11 lineage was introduced in Brazil in the 1950s (1951, 1917-1975 95% HPD) and population dynamics was mostly constant throughout the time. Despite the national vaccination protocols, our results show the widespread dissemination and maintenance of the IBV GI-11 lineage in Brazil and highlight the importance of continuous surveillance to evaluate the impact of currently used vaccine strains on the observed viral diversity of the country.