Published online Jul 19, 2022.
https://doi.org/10.4142/jvs.22082
Molecular characteristics of Budgerigar fledgling disease polyomavirus detected from parrots in South Korea
Abstract
Background
Budgerigar fledgling disease polyomavirus (BFDV) is the pathogen that causes budgerigar fledgling disease in psittacine species. The clinical signs of PBFV infection include ascites, hepatitis, and crop stasis. BFDV is associated with a high mortality rate in nestling birds. In contrast, adult birds only have mild symptoms such as feather dystrophy.
Objectives
This study aimed to determine the prevalence, genetic characteristics, and phylogenetic analysis of BFDV in pet parrots in Korea.
Methods
Fecal and tissue samples were collected from 217 pet parrots from 10 veterinary hospitals including Chungbuk National University Veterinary Hospital. The molecular screening was performed using polymerase chain reaction (PCR) analysis of the small t/large T antigen gene segment. Full-length genome sequencing with the Sanger and phylogenetic analysis were performed on BFDV-positive samples.
Results
The PCR results based on the small t/large T antigen gene marker indicated that BFDV DNA was present in 10 out of 217 screened samples. A whole-genome sequence was obtained from six strains and phylogenetic analysis revealed no significant relationship existed between the species and geographical locations amongst them.
Conclusions
The prevalence of BFDV infection in South Korea is not high when compared to the prevalence of BFDV in other parts of the world, however, it has been reported sporadically in various species and geographic locations. The whole-genome analysis revealed 0.2%–0.3% variation in intragenomic homogeneity among the six strains analyzed. Korean strains are separately on the phylogenetic tree from their counterparts from China and Japan which might reflect the substantial genetic variation.
INTRODUCTION
Budgerigar fledgling disease polyomavirus (BFDV), also known as avian polyomavirus, is a non-enveloped icosahedral DNA virus with a diameter of 40–45 nm. The genome is composed of double-stranded circular DNA with a size of approximately 5,000 nucleotides. This virus belongs to species Aves polyomavirus 1, in the family Polyomaviridae. Within this polyomavirus group, Genus Gammapolyomavirus is the only avian virus genus. [1, 2]. The genome of BFDV has divided into two parts: an early region and a late region. The early region encodes two non-structural proteins, large T antigen and small t antigen. On the other hand, the late region encodes structural proteins, VP1, VP2, VP3, and VP4 [3, 4]. The genome also contains a non-coding control region which was located between the early and late regions. A non-coding control region contains promoter, polyadenylation signal, DNA replication origin, and enhancer regions [4, 5]. Unlike mammalian polyomaviruses, which have well-defined tissue tropism and high host-specificity, avian polyomaviruses, such as BFDV, have low host specificity and can replicate in various organs; where they can be detected in all organs if infected [1, 6]. Transmission of BFDV occurs both vertically and horizontally [7, 8].
BFDV is the etiologic agent of budgerigar fledgling disease, an important immunosuppressive disease in different psittacine species [4, 9, 10]. However, this virus has also been reported in non-psittacine species, such as buzzards, falcons, and pigeons [11, 12]. BFDV infection in nestlings is characterized by severe clinical signs including hepatitis, ascites, and hydropericardium as well as retardation of feather growth, abdomen swelling, crop stasis, and skin discoloration, resulting in the high mortality rate of up to 100%. However, BFDV infection in adults is only associated with mild symptoms such as feather dystrophies (French molt) [9, 13]. BFDV infection has been reported in several species of psittacine birds on almost every continent, including countries such as China [4], Japan [14], Taiwan [2, 15], Turkey [16], Pakistan [17], Thailand [18], Australia [19], New Zealand [20], Egypt [21], Germany [22], Czech Republic [23], Italy [24], Poland [25], Chile [26], Costa Rica [27], and USA [28]. In Korea, the virus was first detected in an Alexandrine parakeet (Psittacula eupatria) in 2014 [29], and in an African gray parrot (Psittacus erithacus) that was co-infected with Psittacine beak and feather disease, as well as and Aspergillus fumigatus and BFDV in 2017 [30]. However, to date, no information on the prevalence and genetic identity of the virus is available. This study aimed to determine the occurrence and genetic profile of BFDV in psittacine in Korea.
MATERIALS AND METHODS
Sample collection
From 2016 to 2019, a total of 217 samples were collected from five local veterinary hospitals in Seoul (18 samples), one in Gyeonggi (one sample), three in Busan (23 samples), and Chungbuk National university veterinary hospital (163 samples). The samples were submitted for the purpose of screening/detection of BFDV by molecular analysis. Of the total 217 samples, 199 samples were obtained from alive parrot droppings (135 samples from clinically healthy parrots, 42 samples from feather problems, 13 samples from gastrointestinal problems, four samples from non-specific problems, two samples from beak problems, two samples from neurological problems, and one sample from respiratory problem) and 18 from dead parrot tissues. Of the total 217 samples, 154 were from companion parrots, 58 were from pet shops, and five were from zoos. Samples were stored frozen at −20°C until processing for molecular analysis.
DNA extraction and BFDV detection
DNA was extracted from tissues and fecal samples using a MagPurix Pathogen/Viral NAs B kit and MagPurix 12s automated nucleic acid purification system (Zinexts Life Science Corp., Taiwan) according to the manufacturer’s instructions. DNA preparations were stored at −20°C until use. DNA preparations were screened for the presence of BFDV by polymerase chain reaction (PCR) amplification of a 310-bp fragment of the small t/large T antigen gene using the primers described by Johne and Müller [11]. The primer sequences were as follows: 5′-CAA GCA TAT GTC CCT TTA TCC C-3′ (Position 4,303–4,324) and 5′-CTG TTT AAG GCC TTC CAA GAT G-3′ (Position 4,612–4,591). The reaction was performed in a total volume of 20 μL that consisted of 1X PCR buffer, 2 mM Mg2+, 200 μM dNTP, 0.2 μM each primer, 2 μL of template, and 1 unit of TaKaRa Ex Taq (TaKaRa Bio Inc., Japan). The following were the PCR amplification conditions: 94°C for 2 min, then 40 cycles of 94°C for 30 s, 54°C for 30 s, 72°C for 30 s, and a final extension at 72°C for 3 min. The PCR products were separated on 1.5% agarose gels and visualized by UV after staining with EcoDye (BIOFACT Co., Korea).
Complete genome amplification and sequence analysis
A total of four overlapping fragments covered the whole length of the genome amplified using four primer pairs (Table 1). The primers were designed based on the sequence of the BFDV-1 strain (accession number AF241168). Each PCR analysis was carried out in a total volume of 50 μL of 1 unit of 2x PCRBIO Ultra Mix Red (PCR Biosystems Inc., UK) 3mM MgCl2, 200 nM of each dNTPs, 200 nM of each primer, and 2 μg of the template DNA. The following were the PCR amplification conditions: 95°C for 5 min, then 40 cycles of 94°C for 1 min, 57°C for 1 min, 72°C for 2 min, and a final extension at 72°C for 10 min. The PCR products were separated on 1% agarose gels and visualized by staining with EcoDye (BIOFACT CO., Korea).
Table 1
Nucleotide sequences of primers designed for this study used to amplify and sequence the full-length genome of budgerigar fledgling disease polyomavirus
PCR products were purified using NucleoSpin Gel and PCR Clean-up Kit (Machery Nagel, Germany) according to the manufacturer’s instructions. Using the prepared PCR products as templates, sequencing reactions were performed according to the manufacturer’s instructions with a BigDye Terminator chemistry (Ver. 3.1) Cycle Sequencing Kit (Applied Biosystems, USA). PCR amplification was performed in 3 μL volumes with 250 ng of templates, 0.5 μL of each primer (3 pmole), 0.87 μL of 5X Sequencing buffer, 1.38 μL of distilled water, and 0.5 μL of BDT v3.1 using a GeneAmp PCR System 9700 (Applied Biosystems). The reaction was performed under the conditions of 30 cycles of denaturation at 96°C for 10 s, annealing at 50°C for 5 s, and extension at 60°C for 4 min. PCR products were purified via ethanol precipitation and resolved on the ABI 3730XL DNA Analyzer (Applied Biosystems).
Phylogenetic analysis
The generated sequences were assembled using Geneious Prime software ver. 2022.0.2. (Biomatters Ltd, New Zealand), edited visually and aligned against sequences in GenBank. Phylogenetic trees were generated using the neighbor-joining algorithm applied to Tamura-Nei genetic distances using MEGA11 software. Trees were tested by selecting a bootstrap method with 1,000 bootstrap replicates. For the purpose of genetic comparison, the reference complete genome sequences of BFDV were obtained from the following GenBank accessions: AB453159-AB453165, AB477106, AF241168-AF241170, AY672646, FJ385773, GU452537, KT203762-KT203769, KX008968-KX008969, M20775, MG148345, MH643735, MK061528, MK516256, MN657184, MT119153, and NC_004764.
RESULTS
BFDV detection
The PCR results based on the small t/large T antigen gene marker indicated that BFDV DNA was present in 10 out of 217 screened samples, with an overall prevalence rate of 4.63%. The 217 parrots represent 20 genera, and a total of 10 BFDV-positive parrots were identified in six genera (Table 2). Detailed information on the species, age, domestic location, sample types, clinical signs, and gross pathology of the BFDV-positive samples are shown in Table 3.
Table 2
Genera of parrots examined in this study and results of budgerigar fledgling disease polyomavirus examined
Table 3
Sample information of budgerigar fledgling disease polyomavirus-positive
Genetic characterization and phylogenetic analysis of BFDV
Whole-genome sequences were successfully generated from 6 out of 10 strains. The generated full length of the six viral genomes was 4,981 nucleotides. The quality or quantity of DNA preparations from 2 out of 4 strains was not enough for whole-genome sequencing (strains 3 and 10). Strains 1 and 2, and 5 and 6 had 100% identical sequences and were collected from the same domestic locations; thus, the whole-genome sequences of strains 1 and 6 were not submitted to the GenBank. The obtained sequences were submitted to the GenBank database and are available under the following accession numbers: MT981254, OM640119-OM640123.
The phylogenetic analysis reveals that all BFDV strains were largely divided into two groups. MT981254 and OM640119 in this study were grouped into group I with reference strains from Japan (AB453159, AB453160, AB453162, and AB453163), Korea (MK516256), and China (FJ385773). OM640120-OM640123 in this study were grouped into group II with other reference strains from Poland (KT203762-KT203769), Hungary (MN657184), China (AY672646, GU452537, MG148345, MH643735, MK061528, and MT119153), Japan (AB453161, AB453164, AB453165, and AB477106), Portugal (KX008968 and KX008969), USA (AF118150), and Germany (AF241168-AF241170, M20775, and NC_004764). MT981254 was slightly related to other group I strains, with a low bootstrap value (< 70). OM640119 was closely related to MK516256 with a bootstrap value of 100. OM640121 and OM640122 were closely related with a bootstrap value of 99. OM640123 and OM640120 were also slightly related to other group II strains, with a low bootstrap value (< 70) (Fig. 1).
Fig. 1
Estimated phylogenetic trees showing relationship among reference sequences of BFDV. The phylogenetic analysis reveals that all BFDV strains were largely divided into group I and group II. The star (⋆) indicates strains obtained in this study. The bootstrap value on the branch represents the percentage of 1000 bootstrap replicates and only values above 70% are displayed. The neighbor-joining algorithm applied to Tamura-Nei genetic distances as implemented by MEGA11 software.
BFDV, budgerigar fledgling disease polyomavirus.
To study the mutation of the BFDV strains in this study, based on the phylogenetic tree of the whole genome, 13 reference strains located in group I and group II were selected and compared (Table 4). Compared to the consensus sequences, there were 52 variable nucleotide positions were found in the whole genome. OM640120 was the most diverse strain with 15 nucleotide exchanges, followed by OM640119 and OM640123 with 13 nucleotide exchanges, OM640121 with 12 nucleotide exchanges, and MT981254 and OM640122 with 11 nucleotide exchanges. The number of amino acid substitutions ranged from 0.0 to 0.3% in VP1, 0.3 to 0.9% in VP2, 0.4 to 1.3% in VP3, 0.0 to 0.6% in VP4, 0.2 to 0.5% in large T-Ag, and 0.0 to 0.7% in small t-Ag (Table 5).
Table 4
Point mutations in six strains of budgerigar fledgling disease polyomaviruses compared with consensus
Table 5
Numbers of nucleotide and amino acid substitutions in this study compared with the consensus sequences
DISCUSSION
In this study, we described the prevalence and whole-genome sequencing analysis of BFDV in 217 samples from different psittacine species from 2016 to 2019 in South Korea. Totally 10 BFDV-positive samples were collected from the species of Rose-ringed parakeet (Psittacula krameri), Budgerigar (Melopsittacus undulatus), African gray parrot (Psittacus erithacus), Monk parakeet (Myiopsitta monachus), Green-cheeked conure (Pyrrhura molinae), and Black-headed parrot (Pionites melanocephalus). The results indicated that the overall infection rate was 4.63% (10/217). Detections of BFDV were reported in Japan [14], Taiwan [15], Turkey [16], Italy [24], USA [31], and Poland [25] with a prevalence of 2.7%, 15.2%, 48.7%, 0.8%, 6.7%, and 22.2%, respectively. The prevalence of BFDV in South Korea in this study was within the prevalence worldwide (0.8%–48.7%), which was significantly lower than that of Taiwan, Turkey, and Poland, whereas slightly higher than that of Japan and Italy.
The phylogenetic tree based on the complete genome sequences of the BFDV indicated that there was no significant relationship between host adaptability to the strain. MT981254 and OM640121 were collected from the same Rose-ringed parakeet species, whereas these strains were grouped separately. In addition, strains 5 and 6 which were obtained from different species (Green-cheeked conure and Rose-ringed parakeet) had 100% identical genomes.
According to previous research, a unique amino acid substitution was reported in VP2 (R234H) isolated from Black-headed parrots when compared with the other 12 sequences [32]. However, this mutation was not observed in OM640122 and OM640123 but observed in other psittacine species (FJ385773 and MK516256). According to another previous study, there were different genetic markers were found on the nucleotide locations at 2474, 2558, 2663, 2774, 2906, 2954, and 4972 differences between group I and group II [5]. These findings similarly confirmed the strains obtained in this study. Additionally, there were another two different genetic markers were found. Group I strains contain 741T and 3506A, whereas group II strains contain 741C and 3506G.
Korea plays an important role as a stopover for migratory birds, a breeding place for residents and migratory birds. Furthermore, the popularity of parrots as exotic pets in Korea has increased annually. Unlike mammalian polyomaviruses, avian polyomaviruses have broad host ranges including Psittaciformes, Falconiformes, and Columbiformes [11, 12]. Phylogenetic analysis shows that at least two groups of BFDV strains have been introduced into South Korea, and it was found that some nucleotide sequence changes have occurred in domestic strains as well. The results of this study suggest that BFDV has circulated in South Korea. Therefore, advanced epidemiological investigations of BFDV are needed to prevent its transmission to captive parrots and wild birds. It is expected that the total of six whole-genome sequences obtained in this study can be used as additional reference sequences for other researchers.
Funding:This research was supported by Brain Pool program funded by the Ministry of Science and ICT through the National Research Foundation of Korea (2019H1D3A2A01060470).
Conflict of Interest:The authors declare no conflicts of interest.
Author Contributions:
Conceptualization: Kim S, Na KJ.
Data curation: Kim S, Kim SJ, Na KJ.
Formal analysis: Kim S, Kim SJ.
Funding acquisition: Na KJ.
Investigation: Kim S, Kim SJ, Na KJ.
Methodology: Kim S, Na KJ.
Project administration: Kim S, Na KJ.
Resources: Kim S, Kim SJ.
Software: Kim S, Na KJ.
Supervision: Na KJ.
Validation: Kim S, Kim SJ, Na KJ.
Visualization: Kim S, Kim SJ, Na KJ.
Writing - original draft: Kim S, Na KJ.
Writing - review & editing: Kim S, Na KJ.
ACKNOWLEDGEMENTS
We thank Dr. Said Amer (Professor of Kafr El Sheikh University and invited professor of Chungbuk National University) who was involved in conceptualization, data curation, formal analysis, investigation, methodology, software, validation, and visualization in this study but passed away before publication. His dedicated help made it possible to write this paper.
References
-
Phalen DN, Wilson VG, Graham DL. Organ distribution of avian polyomavirus DNA and virus-neutralizing antibody titers in healthy adult budgerigars. Am J Vet Res 1993;54(12):2040–2047.
-
-
Bernier G, Morin M, Marsolais G. Papovavirus induced feather abnormalities and skin lesions in the budgerigar: clinical and pathological findings. Can Vet J 1984;25(8):307–310.
-
-
Krautwald ME, Müller H, Kaleta EF. Polyomavirus infection in budgerigars (Melopsittacus undulatus): clinical and aetiological studies. Zentralbl Veterinarmed B 1989;36(6):459–467.
-
-
Adiguzel MC, Timurkan MO, Cengiz S. Investigation and sequence analysis of avian polyomavirus and psittacine beak and feather disease virus from companion birds in eastern Turkey. J Vet Res (Pulawy) 2020;64(4):495–501.
-
-
Fungwitaya P, Bunlertcharoensuk A, Uttamaburana W, Sariya L, Chaichoune K, Ratanakorn P, et al. Prevalence of psittacine beak and feather disease and avian polyomavirus disease infection in captive psittacines in the central part of Thailand by multiplex polymerase chain reaction. J Appl Anim Welf Sci 2009;2(3):33–41.
-
-
El Shahidy M, Helal I, Noaman S. Studies on circo and polyoma viruses in psittacine birds in Egypt, during 2014–2016. Suez Canal Vet Med J 2018;23(2):201–214.
-
-
Valastanova M, Petrikova M, Kulikova L, Knotek Z. Psittacine beak and feather disease virus and avian polyomavirus detection rate in clinically healthy captive birds in the Czech Republic. Vet Med (Praha) 2021;66(2):72–75.
-
-
Piasecki T, Wieliczko A. Detection of beak and feather disease virus and avian polyomavirus DNA in psittacine birds in Poland. Bull Vet Inst Pulawy 2010;54(2):141–146.
-
-
González-Hein G, Gil IA, Sanchez R, Huaracan B. Prevalence of aves polyomavirus 1 and beak and feather disease virus from exotic captive psittacine birds in Chile. J Avian Med Surg 2019;33(2):141–149.
-
-
Dolz G, Sheleby-Elías J, Romero-Zuñiga JJ, Vargas-Leitón B, Gutiérrez-Espeleta G, Madriz-Ordeñana K. Prevalence of psittacine beak and feather disease virus and avian polyomavirus in captivity psittacines from Costa Rica. Open J Vet Med 2013;3(4):240–245.
-
-
Kim HJ, Lee SR, Park CK. First detection of avian polyomavirus by PCR from Alexandrine parakeet (Psittacula eupatria) in Korea. Korean J Vet Serv 2014;37(3):213–218.
-
-
Kang HM, Jang HJ, Seo MK, Lee JW, Na KJ. Psittacine beak and feather disease, budgerigar fledgling disease and aspergillosis in an African gray parrot (Psittacus erithacus). J Vet Clin 2017;34(4):310–312.
-
MeSH Terms
Figures
1 / 1
Tables
1 / 5
ORCID IDs
Funding Information
-
National Research Foundation of Korea
2019H1D3A2A01060470
