Molecular detection and phylogenetic analysis of pigeon circovirus from racing pigeons in Northern China
BMC Genomics volume 23, Article number: 290 (2022)
Pigeon circovirus (PiCV) infections in pigeons (Columba livia) have been reported worldwide. Currently, pigeon racing is becoming increasingly popular and considered to be a national sport in China, and even, the greatest competitions of racing pigeons are taking place in China. However, there are still no epidemiologic data regarding PiCV infections among racing pigeons in China. The purpose of our study was to provide information of prevalence, genetic variation and evolution of PiCV from racing pigeons in China.
To trace the prevalence, genetic variation and evolution of PiCV in sick and healthy racing pigeons, 622 samples were collected from 11 provinces or municipalities in China from 2016 to 2019. The results showed that the positive rate of PiCV was 19.3% (120/622) at the sample level and 59.0% (23/39) at the club level, thus suggesting that the virus was prevalent in Chinese racing pigeons. A sequence analysis revealed that the cap genes of the PiCV strains identified in our study displayed a high genetic diversity and shared nucleotide homologies of 71.9%–100% and amino acid homologies of 71.7%–100%. 28 and 36 unique amino acid substitutions were observed in the Cap and Rep proteins derived from our PiCV strains, respectively. A cladogram representation of PiCV strains phylogeny based on 90 cap gene sequences showed that the strains in this study could be further divided into seven clades (A, B, C, E, G, H, and I) and some of them were closely related to worldwide strains from different types of pigeons. A large number of recombination events (31 events) were also detected in the PiCV genomes from Chinese racing pigeons.
These findings indicate that PiCV strains circulating in China exhibit a high genetic diversity and also contribute to information of prevalence, genetic variation and evolution of PiCV from racing pigeons in China.
Members of the Circoviridae, within the order Circovirus, contain a small, circular, non-enveloped, single stranded DNA genome ranging from 1.7 to 2.5 kb in size . The Pigeon circovirus (PiCV) or the Columbid circovirus (CoCV), together with Porcine circoviruses (PCV) types 1, 2, 3 and 4, Beak and feather disease virus (BFDV), Duck circovirus (DuCV), Goose circovirus (GoCV), Canary circovirus (CaCV), Raven circovirus (RaCV), Starling circovirus (StCV), Swan circovirus (SwCV), Finch circovirus (FiCV), and Gull circovirus (GuCV), etc. belong to the genus Circovirus of the family Circoviridae (http://www.ictvonline.org/). PiCV is a small, circular, non-enveloped virus that contains a single-stranded DNA genome with approximately 2 kb in size. The genome of PiCV has two main open reading frames (ORFs). The ORF-V1 located on the virion sense strand encodes a replication-associated protein (Rep), and the ORF–C1 located on the complementary sense strand encodes a capsid protein (Cap) [2, 3]. The gene forming ORF-C1 of PiCV has been demonstrated to be highly genetically diverse as compared with the gene forming ORF-V1 [1, 4]. The other ORFs including ORF-C2, ORF-C3, and ORF-C4 located on the complementary sense strand encode three viral proteins of unknown functions [2, 5]. The circovirus infection in pigeons was first identified in 1993 in the USA  and had been considered to be strongly associated with young pigeon disease syndrome (YPDS), including weakened racing performance, weight loss, lethargy, anorexia, respiratory distress, and diarrhea .
A variety of methods have been employed to detect PiCV infection in clinical specimens. The original diagnose methods electron microscopy, histology, dot blot hybridization, and in situ hybridization are time consuming [8,9,10,11]. The molecular biology techniques, such as the standard PCR [10, 12], Real-time PCR , NGS techniques , and the loop-mediated isothermal amplification method , enabled more rapid and accurate detection of PiCV infections. As a result, cases of PiCV infections in pigeons were subsequently reported in various countries, including Northern Ireland , Germany [2, 7], Italy , France , Czech Republic , Belgium [18, 19], Poland [20, 21], Slovenia , Hungary , United Arab Emirates , Iran , China [4, 14, 26], Japan , USA [6, 12], Brazil  and Australia . Previous studies also demonstrated that PiCV had been detected in different types of pigeons, including racing, fancy, feral and meat pigeons [14, 20, 30]. However, all of the above-mentioned birds belong to Columba livia species.
In China, the first PiCV infection was detected from meat pigeons in Zhejiang province in 2009, and the full genome was sequenced . In recent years, several studies have proved that PiCV was prevalent among meat pigeons in eastern and southern China [4, 14]. Recently, the competitions of racing pigeons are becoming increasingly popular in China, and even, pigeon racing is considered to be a national sport. The Chinese Association of Racing Pigeon Breeders has over 5 million members. Chinese breeders compete in numerous races for different distances in their sections throughout the racing season, and more than 25 million racing pigeons from approximately 750 racing clubs were selected to take part in the competitions every year. However, there are no epidemiologic data on PiCV infections among racing pigeons. In order to investigate the prevalence, evolution, and genome characterization of PiCV in racing pigeons in China, an extensive epidemiological investigation and bioinformatic analysis of PiCV from racing pigeons were undertaken in this study. The purpose of this study is to provide novel epidemiologic data and genome characterization for PiCV strains identified from racing pigeons in China.
PCR detection and prevalence analysis of PiCV infection in Chinese racing pigeons
In this study, PiCV was detected in samples from both sick and healthy racing pigeons. The results indicated that the positive rates of PiCV were 19.3% (120/622) at the sample level. Among the samples, 14 out of 51 sick pigeons (positive rate, 27.5%) and 106 out of 571 healthy pigeons (positive rate, 18.6%) were identified to be PiCV positive, respectively. In this study, 23 out of 39 racing pigeon clubs were positive for PiCV (positive rate, 59.0%) (Table 1). Detailed information for 120 PiCV positive samples from racing pigeons in China is shown in Table 1.
Genome characterization of PiCV genomes
A total of sixty-seven full genome sequences were obtained and deposited in GenBank with accession numbers: MW181925 to MW181991 (Table S2). The assembled circular whole genome sequences were detected as eleven sizes ranging from 2030 to 2045 nt with their cap gene ranging from 813 to 828 nt in length and rep gene of 948 or 954 nt as shown in Table S2. Among these variable genome length, the most common genome sizes were 2037 (n = 21) and 2042 (n = 13) nucleotides, respectively (Table S2). The sequence comparison among the 67 identified PiCV strains revealed nucleotide homologies of 84.2%–100% and these sequences exhibited nucleotide homologies of 83.0%–97.8% and 82.0%–98.3% as compared with Chinese PiCV reference strains and the ones from the other countries (Table S3).
Analysis of nucleotide sequence of the cap genes and amino acid sequence of the Cap proteins
To explore the genetic diversity of PiCV strains, the cap genes of 90 PiCV strains from the 120 PiCV positive samples were sequenced. Detailed information for the 90 cap genes was shown in Table S2. The 90 cap genes were used to compare with the reference sequences from China and other countries, respectively. The results showed that the 90 cap genes ranged from 813 to 828 nt in length. Fourteen, thirteen, one, fifty-three, seven and two of the 90 cap nucleotide sequences were 813 nt, 816 nt, 819 nt, 822 nt, 825 nt, and 828 nt in length, encoding a Cap protein of 270, 271, 272, 273, 274, and 275 residues, respectively (Table S2). In addition, we found that ATT and GTG also existed in the position of the start codon site. The sequence comparison of the 90 identified cap genes revealed nucleotide homologies of 71.9%–100% and deduced amino acid homologies of 71.7%–100%, and the cap genes exhibited low sequence similarities with Chinese PiCV reference strains (73.0%–99.6% nucleotide identity, 72.3%–100% amino acid identity) and the PiCV reference strains from other countries (68.8%–98.4% nucleotide identity, 63.6%–100% amino acid identity) (Table S3).
To investigate variations in the deduced amino acid sequences, the amino acid sequences of 90 identified Cap proteins and the reference strains were aligned. The results showed that there were nine major locations of deletion (compared to the consensus sequence) among the Cap proteins including locations 7, 24, 29, 30, 35, 58, 130, 182, and 266 (Fig. 1). A comparison of entropy (Hx) in amino acid sequences of Cap proteins showed that the Hx of amino acid sequences at most positions from the 90 identified Cap proteins were higher than that of the PiCV reference strains (Fig. 2). In addition, some unique amino acid substitutions at 28 different positions were observed among the 90 Cap proteins as shown in Fig. 3.
Analysis of nucleotide sequence of the rep genes and amino acid sequence of the Rep proteins
The rep genes for 68 out of 120 PiCV positive samples were successfully sequenced. Detailed information was shown in Table S2. The 68 identified rep genes were used to compare with PiCV reference sequences from China and other countries. The results showed that all of the 68 rep genes used the ATG start codon and were identified as two sizes: 948 and 954 nt. Forty of the 68 rep nucleotide sequences were 948 nt in length, encoding a Rep protein with 315 residues. Twenty-eight of the 68 rep nucleotide sequences were 954 nt in size, encoding a Rep protein with 317 residues (Table S2). The difference in size was due to a 2 amino acid deletion at positions 2 and 3. The sequence alignment revealed nucleotide homologies of 90.3%–100% and deduced amino acid homologies of 92.7%–100% among the 68 rep genes. These sequences exhibited higher similarities with PiCV reference strains from China (89.0%–99.2% nucleotide identity, 89.2%–99.6% amino acid identity) and other countries (89.5%–98.3% nucleotide identity, 90.5%–99.3% amino acid identity) as compared to the cap genes (Table S3). To investigate variations in the deduced amino acid sequences of rep gene products, the amino acid sequences of 68 identified rep genes and the reference strains were aligned. The results showed that some unique amino acid substitutions at 36 different positions were observed among the 68 identified PiCV strains (Data not shown).
Phylogenetic analysis of the cap genes
A phylogenetic analysis using 126 cap gene reported in GenBank, and cap genes of 90 PiCV strains obtained in our study was performed to investigate genetic relationship (Table S1 and S2). The results indicated that a total of 216 analyzed PiCV strains could be arbitrarily divided into 9 clades, which were defined by the contractual letters A-I (Fig. 4). The 90 cap gene sequences reported in the present study were clustered into 7 clades: A (14.1%), B (60.9%), C (62.2%), E (25.0%), G (51.4%), H (28.6%), and I (40.0%). As shown in Fig. 4, the PiCV strains from different geographical locations were clustered into the same clades, such as QYQX1/HE/2018/MW181909, TY1/SN/2016/MW181929, DS1/GS/2018/MW181970, and QD3/SN/2018/MW181965. On the other hand, the PiCV strains from the same geographical locations were divided into different clades, such as LH1/HE/2018/MW181918 (clade B), LH2/HE/2019/MW181988 (clade A), and LH3/HE/2019/MW181989 (clade C). Furthermore, the PiCV strains collected from the same province belonged to different clades. For example, the 78 PiCV strains isolated from the pigeons in the racing clubs in Shannxi province, belonged to six separate clades except for clade D, F, and H.
Recombination analysis in PiCV genomes
For 67 identified PiCV strains, 31 recombination events were detected using RDP4  (Table S4). These recombination events indicated that the strains originally from different types of pigeons or from different geographical locations may have some evidence of recombination. For example, a segment (~ 1001 bp) of three PiCV strains (LT3/SN/2018/MW181940, LT2/SN/2018/MW181939, and WQ6/SN/2018/MW181948) obtained in this study was likely descended from a large number of ancestral PiCV genomes that originated from Chinese meat pigeons, German racing pigeons, British racing pigeons, Polish fancy pigeons, and French meat pigeons (event 2). The possible breakpoints for recombination were determined during the recombination analysis. The results also indicated that the recombination breakpoint hot plots located within both the intergenic region between the rep and cap stop codons and near the virion strand origin of replication. Furthermore, the recombination breakpoint cold spots were detected within the central region of cap gene, which indicated few evidence for recombination within cap gene (Fig. 5). Additionally, four PiCV genomes (TF1/SN/2016/MW181925, DA1/XJ/2018/MW181982, QD7/SN/2019/MW181990, and LH2/HE/2019/MW181988) had no evidence for recombination.
In recent studies, Rotavirus A G18P  had been confirmed as a primary cause of YPDS-like diseases in domestic pigeons . Furthermore, PiCV is also supposed to be an etiological agent of YPDS mostly affecting young pigeons worldwide . PiCV has been reported to be prevalent in at least sixteen countries. In China, PiCV was first detected in meat pigeons in 2009 . In recent years, the epidemiological survey showed that the positive rates of PiCV infection in Chinese meat pigeons were 19.67% and 75.3% in the poultry farms of China in 2009 and 2015 [4, 14], implying that PiCV is widely distributed among meat pigeon populations in China. However, the epidemiology and distribution of PiCV in the racing pigeons is unknown. The main objective of this work was to evaluate the genetic diversity and epidemiology of PiCV strains circulated in the racing pigeons of China. Positive samples were detected from seven provinces and the prevalence rates among these provinces were variant. Overall, our data implied that PiCV was also widely distributed in diseased racing pigeons and healthy racing pigeons in China for the first time.
Previous studies showed that the length of PiCV genome was 2031–2043 nt [14, 28, 29]. However, our findings first showed that the 2030 nt (n = 1), 2044 nt (n = 2), and 2045 nt (n = 1) complete genome were present in the PiCV positive samples. Despite many genetic diversities were found in the genome of PiCV strains, no mutation was found in the conserved nonanucleotide motif (TAGTATTAC) located at the apex of a potential stem-loop which was putatively associated with the initiation of rolling circle replication (RCR) [2, 3]. Previous studies have shown many genetic diversities in the cap gene [14, 23, 34]. Our findings showed that the 90 identified cap genes exhibited a higher diversity as compared with the reference strains. Some unique amino acid substitutions at 28 different positions were observed among the Cap proteins. In addition, the cap nucleotide sequence with 828 nt in length encoding a novel Cap protein of 275 amino acids was first identified in two Chinese PiCV strains, TY3/SN/2016/MW181931 (clade E) and WL4/SN/2018/MW181959 (clade E). Higher diversity of Cap protein versus the Rep protein is due to the fact that the Cap protein, as the protein shell of the virus, is exposed to the host’s immune cells, resulting in a stimulation of a cascade of immune responses. Positive selection occurs during this process. Some mutations can change the structure of Cap protein. This may lead to a better binding of virus with the receptors of target cells, which increases the infection ability of virus in cells. Moreover, these mutations may also protect virus from being neutralized by antibodies [1, 35, 36]. In this study, we identified many mutations in both N- and C-terminus of the Cap protein. According to the investigations on Cap protein of BFDV, the arginine-rich N-terminus has two functions (nuclear localization and nucleic acid binding) that enable as penetration of the viral genome into the host cell nucleus through nuclear pore complex . Like BFDV, the N-terminus of PiCV Cap protein is also rich in arginine and has been predicted to be a nuclear localization signal and a nucleic acid binding domain by software , which means that this region is highly possible to be functionally similar to the corresponding region of the Cap protein of BFDV. Moreover, the results of this study showed that most of the amino acid deletion sites in Cap protein were concentrated in the N-terminus and many novel mutations also have been identified in this region. However, due to the lack of relevant experimental evidence, we cannot infer whether these mutations and deletions could affect the function of Cap protein. The Cap protein has been used as a coating antigen in ELISA due to its antigenic activity and can be recognized by PiCV specific antibodies [39, 40]. If the N-terminal region of the PiCV Cap protein, which is located within the capsid, resembles the corresponding region in BFDV, the N-terminal amino acid mutation and deletion may not affect the ability to bind to the antibody as a coating antigen . In addition, there are also many mutations in the C-terminus where most sites are composed of hydrophilic amino acids. It means that this region may locate on the outside of the protein and may have ability for binding with antibody. Therefore, mutations of this region may lead to changes in antigenicity. Cap protein, which contains neutralizing antibody epitopes, can be considered as a potential antigen candidate in sub-unit vaccine development [41, 42]. Moreover, the sub-unit vaccines based on PCV2 recombinant capsid proteins are successfully used in the prevention of PCV2-SD [43,44,45,46]. It means that the development of a sub-unit vaccine may protect pigeons from infection with PiCV . Due to a high genetic diversity of Cap proteins may lead to differences in antigenicity between different strains, sub-unit vaccines may also be diversified in future. To sum up, since there is little research on the function of PiCV Cap protein, the specific meaning of these mutations cannot be accurately explained. Therefore, further experiments are needed to explain if these mutations and deletions affect the function of Cap protein. According to the obtained nucleotide sequences of cap gene, besides ATG and ATA [2, 5, 14], we found that ATT and GTG also existed in the position of the start codon site through sequence alignment. Moreover, a similar situation has been found in GuCV  and BFDV . Since no available cell lines have been found to culture PiCV in vitro [49, 50] and lack of commercial antibodies against PiCV Cap protein, it is hard to verify that indeed a full-length natural Cap protein is translated from the described alternative start codons. In all 68 identified Rep proteins, three amino acid motifs named FTLNNP (position 41–46), HLQGF (position 78–82), and YCSK (position 116–119) that putatively associated with RCR  were completely conserved. In addition, a fourth motif, which was putatively associated with dNTPase activity, namely GKS (position 197–199) , was also completely conserved. These sequence characteristics help understand the genetic diversity of PiCV. However, since the PCR amplification strongly relies on nucleotide homology in the regions where primers anneal, the samples containing more divergent viruses may not be obtained the full genome sequence when amplifying the viral genome. Thus, the actual diversity of circoviruses in Chinese pigeons may be even greater. This can be addressed by future studies employing next generation sequencing methods. Besides the pathogenicity of different PiCV strains still need to be validated in future.
In the present study, the 90 PiCV strains were divided into seven clades (A, B, C, E, G, H, and I) based on a phylogenetic analysis of cap gene. Additionally, we found that PiCV strains isolated from the same club belonged to different clades and shared a low sequence identity. These data suggested that the infection and evolution of PiCV in Chinese racing pigeons might have different evolutionary origins. This may result from the fact that the import of racing pigeons from all around the world to China is very significant and a lack of oversight in the international pigeon trade (racing pigeons mainly) [20, 51]. In addition, another possible explanation is that the intensive trafficking of pigeons between clubs and racing events in China. Interestingly, two isolates, QYQX1/HE/2018/MW181909 (clade A) and QYQX2/HE/2018/MW181910 (clade C), with lower identity (75.5% nucleotide identity and 76.8% amino acid identity) in the cap gene were detected from the same pigeon, suggesting a horizontal transmission occurred among the racing pigeons in the same racing clubs [6, 52]. The complexity of epidemic PiCV strains in China may cause difficulties to protect pigeon from infection by vaccines in future.
Viral recombination had been proved to play a significant role in the evolution of many ssDNA viruses . The extensive recombination events had been reported in PiCV genomes and other circoviruses, such as PCV [54, 55] and BFDV [56, 57]. In this study, a large number of recombination events (31 events) were detected in the 67 identified PiCV strains. The high intensity of recombination may suggest that frequent co-infections of different PiCV strains. Moreover, recently, Khalifeh A. et al. has demonstrated that recombination plays a role in PiCV evolution through animal experiment . Thus, the recombination seemed to be a key mechanism for PiCV evolution [21, 23, 28, 29].
In conclusion, our study demonstrated that PiCV infection in racing pigeons was widespread in northern China and revealed the characteristics of the PiCV genome. Furthermore, the identified PiCV strains displayed a high genetic diversity. Our data also demonstrated that PiCV in Chinese racing pigeons had an extensive recombination for the first time. The high genetic diversity of PiCV may be derived from two possible facts: 1) the racing pigeons in China are imported from all around the world; and 2) a lack of oversight in the international pigeon trade (mainly refer to racing pigeons). In addition, another possible explanation is that the intensive trafficking of pigeons between clubs and racing events in China. Furthermore, the high intensity of recombination may suggest that frequent co-infections of different PiCV strains. Overall, clubs and racing events could provide a suitable chance for viral recombinants, which may pose a threat for the domestic pigeons, and has a possibility of spilling over onto wild species.
Materials and methods
In total, 622 samples were collected from 11 provinces or municipalities in China from November 2016 to August 2019. Of these samples, 571 serum samples were collected from 571 healthy racing pigeons from 39 racing pigeon clubs in 11 provinces or municipalities of China, including Beijing, Hebei, Liaoning, Shanxi, Gansu, Qinghai, Shandong, Shaanxi, Xinjiang Uygur Autonomous Region, Ningxia Hui Autonomous Region and Inner Mongolia Autonomous Region. Animals obtained from the racing pigeon clubs were returned to the club owners after collecting blood samples. Moreover, the corpses of 51 sick pigeons from 7 racing pigeon clubs in 5 provinces or municipalities were submitted to our laboratory at College of Veterinary Medicine, Northwest A&F University to determine the pathogen of the disease. The 0.2 g sections of internal organs (the heart, liver, spleen, lung, and intestine) were obtained during the postmortem examination of all mentioned dead pigeons. All samples were stored at − 80 °C before DNA extraction.
The tissue samples were ground to powder with liquid nitrogen and diluted with three volumes of phosphate‐buffered saline. The samples (1 ml) were centrifuged at 12,000 × g for 10 min at 4 °C and the supernatants were transferred into a 1.5 ml tube. Nucleic acids from tissues and serum samples were extracted by using an EasyPure® Viral DNA/RNA Kit (TransGen Biotech, Beijing, China) according to the manufacturer’s instructions. The extracted genomic DNA was stored at -20 °C before use.
PCR detection of PiCV
The PiCV was first detected using a PCR method targeting a 326-base fragment of cap gene as described by Freick et al. (2008). The primer sequences were: PiCV-s, 5’-TTGAAAGGTTTTCAGCCTGGC-3’ and PiCV-as, 5’-AGGAGACGAAGGACACGCCTC-3’ . The full genomes of PiCV for all the positive samples were amplified by PCR using the primers: PiCV-1F, 5’-ACCCGCGACTTGGAGCCACGGAG-3’ and PiCV-1R, 5’-TTCGCTCCCGCATTCGCGGTCGCT-3’; PiCV-2F, 5’-GACACTAGTAAAGGGACCCAAGCCA-3’ and PiCV-2R, 5’-AAGCCTTGCAGATGCGGGGT-3’, respectively. PCR was performed by using Q5 Hot Start High-Fidelity 2 × Master Mix (NEB, MA, USA). The contents of the three reactions mixture in a 50 μl reaction volume were as follows: 0.5 μM forward primer, 0.5 μM reverse primer, 1 μg genomic DNA, 25 μl Q5 High-Fidelity 2 × Master Mix (NEB, MA, USA) and an appropriate volume of Nuclease-Free Water. The cycling parameters were 30 cycles of 98 °C for 10 s, 55 °C for 30 s and 72 °C for 30 s, followed by a final extension at 72 °C for 5 min using an automated BioRad T100 Thermal Cycler (Bio-Rad Laboratories, Inc., CA, USA). The PCR products (5 μl) were resolved on 1% (w/v) agarose gels, and followed by staining with ethidium bromide. Finally, the bands of nucleic acid were visualized with UV illumination inside a gel documentation apparatus (Bio-Rad Laboratories, Inc., CA, USA) and saved as digital photographs.
Genome characterization and sequence analysis of the cap, rep genes and PiCV genome
The cap gene sequences, rep gene sequence and the full genome sequences of PiCV obtained in this study have been deposited in GenBank, respectively. The sequences obtained in the current study were compared with all available PiCV sequence data from different geographical locations within China (n = 60) and the rest of the world (n = 73) from the National Center for Biotechnology Information (NCBI) nucleotide database (Table S1). Multiple sequence alignments were carried out by ClustalW algorithm using MEGA 5.0 software  and the homology among nucleotide and amino acid sequences was determined by using BioEdit v. 7.0.5 software . A divergence analysis of the Cap protein of PiCV was performed using WebLogoTm (http://weblogo.threeplusone.com/), an online software for sequence logo generator . A comparison of Hx in amino acid sequences of Cap proteins of PiCV was carried out by using GraphPad Prism 7.0 software (GraphPad Software, Inc., USA).
The entire cap gene sequences of PiCV obtained in this study were used for a phylogenetic analysis. All available PiCV sequence data from different geographical locations within China (n = 54) and the rest of the world (n = 72) were retrieved from the NCBI nucleotide database as reference sequences (Table S1). Maximum likelihood phylogenetic trees were constructed with the GTR + G + I model by using the MEGA 5.0 software , with partial deletion to handle alignment gaps and 1,000 bootstrap iterations. The cladogram was generated and annotated with the Interactive Tree Of Life (iTOL) software (http://itol.embl.de/) .
In order to analyze the potential recombinant evidence, an integrated software package for the recombination detection program RDP4  was used to detect potential recombinant strains, parental strains, and possible recombination breakpoints. Seven methods (RDP , GENECONV , BootScan , MaxChi , Chimaera , SiScan , and 3Seq ) were implemented using the RDP4 program . Recombination events were identified by at least three of the aforementioned methods, and the p-value < 0.05 were considered plausible recombinant events. The sequences in the analyzed dataset that most closely resembled the parental sequences of recombinants were defined as either “minor parents” or “major parents” based on the size of the genome fragments which had contributed to the detected recombinants (with the major parent contributing the larger fragment and the minor parent the smaller).
Availability of data and materials
The nucleotide sequence data obtained in this study are available in GenBank, and their accession numbers are MW181901- MW181991. The information of PiCV references sequences is shown in Table S1.
Beak and feather disease virus
Open reading frames
Young pigeon disease syndrome
Rolling circle replication
National Center for Biotechnology Information
Stenzel T, Koncicki A. The epidemiology, molecular characterization and clinical pathology of circovirus infections in pigeons – current knowledge. Vet Quart. 2017;37(1):166–74.
Mankertz A, Hattermann K, Ehlers B, Soike D. Cloning and sequencing of columbid circovirus (CoCV), a new circovirus from pigeons. Arch Virol. 2000;145(12):2469–79.
Todd D, Weston JH, Soike D, Smyth JA. Genome Sequence Determinations and Analyses of Novel Circoviruses from Goose and Pigeon. Virology. 2001;286(2):354–62.
Zhang Z, Dai W, Wang S, Dai D. Epidemiology and genetic characteristics of pigeon circovirus (PiCV) in eastern China. Arch Virol. 2015;160(1):199–206.
Todd D, Fringuelli E, Scott AN, Borghmans BJ, Duchatel JP, Shivaprasad HL, et al. Sequence comparison of pigeon circoviruses. Res Vet Sci. 2008;84(2):311–9.
Woods LW, Latimer KS, Barr BC, Niagro FD, Campagnoli RP, Nordhausen RW, et al. Circovirus-Like Infection in a Pigeon. J Vet Diagn Invest. 1993;5(4):609–12.
Raue R, Schmidt V, Freick M, Reinhardt B, Johne R, Kamphausen L, et al. A disease complex associated with pigeon circovirus infection, young pigeon disease syndrome. Avian Pathol. 2005;34(5):418–25.
Smyth JA, Weston J, Moffett DA, Todd D. Detection of Circovirus Infection in Pigeons by in Situ Hybridization Using Cloned DNA Probes. J Vet Diagn Invest. 2001;13(6):475–82.
Soike D, Hattermann K, Albrecht K, Segalés J, Domingo M, Schmitt C, et al. A diagnostic study on columbid circovirus infection. Avian Pathol. 2001;30(6):605–11.
Todd D, Duchatel JP, Weston JH, Ball NW, Borghmans BJ, Moffett DA, et al. Evaluation of polymerase chain reaction and dot blot hybridisation tests in the diagnosis of pigeon circovirus infections. Vet Microbiol. 2002;89(1):1–16.
Tsai SS, Chang YL, Huang YL, Liu HJ, Ke GM, Chiou CJ, et al. Development of a loop-mediated isothermal amplification method for rapid detection of pigeon circovirus. Arch Virol. 2014;159(5):921–6.
Roy P, Dhillon AS, Lauerman L, Shivaprasad HL. Detection of Pigeon Circovirus by Polymerase Chain Reaction. Avian Dis. 2003;47(1):218–22.
Duchatel JP, Todd D, Willeman C, Losson B. Quantification of pigeon circovirus in serum, blood, semen and different tissues of naturally infected pigeons using a real-time polymerase chain reaction. Avian Pathol. 2009;38(2):143–8.
Wang K, Zhuang Q, Qiu Y, Wang T, Chen J. Genome sequence characterization of pigeon circoviruses in China. Virus Res. 2017;233:1–7.
Franciosini MP, Fringuelli E, Tarhuni O, Guelfi G, Todd D, Casagrande Proietti P, et al. Development of a polymerase chain reaction-based in vivo method in the diagnosis of subclinical pigeon circovirus infection. Avian Dis. 2005;49(3):340–3.
Abadie J, Nguyen F, Groizeleau C, Amenna N, Fernandez B, Guereaud C, et al. Pigeon circovirus infection: Pathological observations and suggested pathogenesis. Avian Pathol. 2001;30(2):149–58.
Taras L, Kubicek O, Juranova R, Jurajda V. The First Demonstration of Pigeon Circovirus Infection in the Czech Republic Based on Histology and Nested PCR. Acta Vet Brno. 2003;72(4):577–82.
Duchatel JP, Todd D, Curry A, Smyth JA, Bustin JC, Vindevogel H. New data on the transmission of pigeon circovirus. Vet Rec. 2005;157(14):413–5.
Duchatel JP, Todd D, Smyth JA, Bustin JC, Vindevogel H. Observations on detection, excretion and transmission of pigeon circovirus in adult, young and embryonic pigeons. Avian Pathol. 2006;35(1):30–4.
Stenzel T, Pestka D, Tykałowski B, Śmiałek M, Koncicki A. Epidemiological investigation of selected pigeon viral infections in Poland. Vet Rec. 2012;171(22):562.
Stenzel T, Piasecki T, Chrzastek K, Julian L, Muhire BM, Golden M, et al. Pigeon circoviruses display patterns of recombination, genomic secondary structure and selection similar to those of Beak and feather disease viruses. J Gen Virol. 2014;95(Pt 6):1338–51.
Krapez U, Slavec B, Steyer AF, Pintaric S, Dobeic M, Rojs OZ, et al. Prevalence of pigeon circovirus infections in feral pigeons in Ljubljana. Slovenia Avian Dis. 2012;56(2):432–5.
Cságola A, Lőrincz M, Tombácz K, Wladár Z, Kovács E, Tuboly T. Genetic diversity of pigeon circovirus in Hungary. Virus Genes. 2012;44(1):75–9.
Ledwoń A, Bailey T, O’Donovan D, Mckeown S, Lloyd C, Wieckowski T, et al. Prevalence of circovirus and adenovirus in pigeons in Dubai. Med Weter. 2011;67(11):752–6.
Mahzounieh M, Heidari Khoei H, Ghasemi Shamsabadi M, Dastjerdi A. Detection and phylogenetic characterization of Columbid circoviruses in Chaharmahal va Bakhtiari province. Iran Avian Pathol. 2014;43(6):524–8.
Liao P, Wang K, Tsai S, Liu H, Huang B, Chuang K. Recurrent positive selection and heterogeneous codon usage bias events leading to coexistence of divergent pigeon circoviruses. J Gen Virol. 2015;96(8):2262–73.
Yamamoto E, Ito H, Kitamoto E, Morinishi K, Yano A, Miyoshi S, et al. Complete genome sequence of pigeon circovirus detected in racing pigeons in western Japan. Virus Genes. 2015;51(1):140–3.
Loiko MR, Junqueira DM, Varela A, Tochetto C, Scheffer CM, Lima DA, et al. Columbid circoviruses detected in free ranging pigeons from Southern Brazil: insights on PiCV. Arch Virol. 2018;163:3083–90.
Sarker S, Das S, Ghorashi SA, Forwood JK, Raidal SR. Pigeon circoviruses from feral pigeons in Australia demonstrate extensive recombination and genetic admixture with other circoviruses. Avian Pathol. 2019;48(6):512–20.
Stenzel T, Pestka D, Choszcz D. The prevalence and genetic characterization of Chlamydia psittaci from domestic and feral pigeons in Poland and the correlation between infection rate and incidence of pigeon circovirus. Poultry Sci. 2014;93(12):3009–16.
Yu X, Zhu C, Zheng X, Mu A, Yu H. Cloning and analysis of the complete genomes of pigeon circovirus from Zhejiang Province. Chinese journal of virology. 2009;25(5):355–61.
Martin D, Murrell B, Golden M, Khoosal A, Muhire B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evolution. 2015;1(1):vev003-vev.
Rubbenstroth D, Ulrich R, Wylezich C, Rautenschlein S, Beer M, Mohr L. First experimental proof of Rotavirus A (RVA) genotype G18P inducing the clinical presentation of ‘young pigeon disease syndrome’ (YPDS) in domestic pigeons (Columba livia). Transbound Emerg Dis. 2020;67(4):1507–16.
Stenzel T, Pestka D. Occurrence and genetic diversity of pigeon circovirus strains in Poland. Acta Vet Hung. 2014;62(2):274–83.
Bassami MR, Berryman D, Wilcox GE, Raidal SR. Psittacine beak and feather disease virus nucleotide sequence analysis and its relationship to porcine circovirus, plant circoviruses, and chicken anaemia virus. Virology. 1998;249(2):453–9.
Bennett S, Holmes E, Chirivella M, Rodriguez D, Beltran M, Vorndam V, et al. Molecular evolution of dengue 2 virus in Puerto Rico: Positive selection in the viral envelope accompanies clade reintroduction. J Gen Virol. 2006;87(4):885–93.
Heath L, Williamson A-L, Rybicki E. The Capsid Protein of Beak and Feather Disease Virus Binds to the Viral DNA and Is Responsible for Transporting the Replication-Associated Protein into the Nucleus. J Virol. 2006;80(14):7219–25.
Lai G, Lin Y, Tsai Y, Lien Y, Lin M, Chen H, et al. High yield production of pigeon circovirus capsid protein in the E coli by evaluating the key parameters needed for protein expression. BMC Vet Res. 2014;10(1):1–11.
Daum I, Finsterbusch T, Härtle S, Göbel T, Mankertz A, Korbel R, et al. Cloning and expression of a truncated pigeon circovirus capsid protein suitable for antibody detection in infected pigeons. Avian Pathol. 2009;38(2):135–41.
Stenzel T, Wozniakowski G, Dziewulska D, Choszcz D, Tykałowski B, Smialek M, et al. Application of pigeon circovirus recombinant capsid protein for detecting anti-PiCV antibodies in the sera of asymptomatic domestic pigeons and the potential use of a combination of serological and molecular tests for controlling circovirus infections in pigeon breeding flocks. Poultry Sci. 2016;96(2):303–8.
Gai W, Zheng W, Zhao Z, Wong G, Sun P, Yan L, et al. Assembly of pigeon circovirus-like particles using baculovirus expression system. Microb Pathogenesis. 2020;139:103905.
Stenzel T, Dziewulska D, Tykałowski B, Śmiałek M, Kowalczyk J, Koncicki A. Immunogenicity of Pigeon Circovirus Recombinant Capsid Protein in Pigeons. Viruses. 2018;10(11):596.
Fort M, Marina S, Perez-Martin E, Nofrarías M, Mateu E, Segalés J. One dose of a porcine circovirus 2 (PCV2) sub-unit vaccine administered to 3-week-old conventional piglets elicits cell-mediated immunity and significantly reduces PCV2 viremia in an experimental model. Vaccine. 2009;27(30):4031–7.
Li P, Qiao X, Zheng Q, Hou J. Immunogenicity and immunoprotection of porcine circovirus type 2 (PCV2) Cap protein displayed by Lactococcus lactis. Vaccine. 2015;34(5):696–702.
Blanchard P, Mahe D, Cariolet R, Keranflec’h A, Baudouard MA, Cordioli P, et al. Protection of swine against post-weaning multisystemic wasting syndrome (PMWS) by porcine circovirus type 2 (PCV2) proteins. Vaccine. 2003;21(31):4565–75.
Zhu S, Zhang C, Wang J, Wei I, Quan R, Yang J, et al. Immunity Elicited by an Experimental Vaccine Based on Recombinant Flagellin-Porcine Circovirus Type 2 Cap Fusion Protein in Piglets. PLoS One. 2016;11(2):e0147432.
Todd D, Scott ANJ, Fringuelli E, Shivraprasad HL, Gavier-Widen D, Smyth JA. Molecular characterization of novel circoviruses from finch and gull. Avian Pathol. 2007;36(1):75–81.
Bassami MR, Ypelaar I, Berryman D, Wilcox G, Raidal S. Genetic Diversity of Beak and Feather Disease Virus Detected in Psittacine Species in Australia. Virology. 2001;279(2):392–400.
Huang HY, Silva BBI, Tsai SP, Tsai CY, Tyan YC, Lin TC, et al. Immunogenicity and Protective Activity of Pigeon Circovirus Recombinant Capsid Protein Virus-Like Particles (PiCV rCap-VLPs) in Pigeons (Columba livia) Experimentally Infected with PiCV. Vaccines. 2021;9(2):98.
Santos HM, Chen CC, Tsai C-Y, Hsish YC, Chung FC, Tyan Y-C, et al. Influence of pigeon interferon alpha (PiIFN-α) on pigeon circovirus (PiCV) replication and cytokine expression in Columba livia. Vet Microbiol. 2020;242:108591.
Ashton WL. The risks and problems connected with the import and export of captive birds. Br Vet J. 1984;140(4):317–27.
Gerdes GH. Two very small viruses–a presumptive identification. J S Afr Vet Assoc. 1993;64(1):2.
Lefeuvre P, Lett JM, Varsani A, Martin DP. Widely conserved recombination patterns among single-stranded DNA viruses. J Virol. 2009;83(6):2697–707.
Kleymann A, Soto E, Illanes O, Malik YS, Fuentealba C, Ghosh S. High rates of detection and complete genomic analysis of porcine circovirus 2 (PCV2) in the Lesser Antilles island of St. Kitts: Identification of PCV2b-PCV2d recombinants. Transbound Emerg Dis. 2020;67(6):2282–9.
Wei C, Lin Z, Dai A, Chen H, Ma Y, Li N, et al. Emergence of a novel recombinant porcine circovirus type 2 in China: PCV2c and PCV2d recombinant. Transbound Emerg Dis. 2019;66(6):2496–506.
Julian L, Piasecki T, Chrzastek K, Walters M, Muhire B, Harkins GW, et al. Extensive recombination detected among beak and feather disease virus isolates from breeding facilities in Poland. J Gen Virol. 2013;94(Pt 5):1086–95.
Varsani A, Regnard GL, Bragg R, Hitzeroth II, Rybicki EP. Global genetic diversity and geographical and host-species distribution of beak and feather disease virus isolates. J Gen Virol. 2011;92(Pt 4):752–67.
Khalifeh A, Kraberger S, Dziewulska D, Varsani A, Stenzel T. A Pilot Study Investigating the Dynamics of Pigeon Circovirus Recombination in Domesticated Pigeons Housed in a Single Loft. Viruses. 2021;13(6):964.
Freick M, Müller H, Raue R. Rapid detection of pigeon herpesvirus, fowl adenovirus and pigeon circovirus in young racing pigeons by multiplex PCR. J Virol Methods. 2008;148(1):226–31.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.
Hall T. BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.
Crooks GE, Hon G, Chandonia JM, Brenner SE. WebLogo: a sequence logo generator. Genome Res. 2004;14(6):1188–90.
Letunic I, Bork P. Interactive Tree Of Life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics. 2006;23(1):127–8.
Martin D, Rybicki E. RDP: detection of recombination amongst aligned sequences. Bioinformatics. 2000;6(16):562–3.
Padidam M, Sawyer S, Fauquet CM. Possible Emergence of New Geminiviruses by Frequent Recombination. Virology. 1999;265(2):218–25.
Martin DP, Posada D, Crandall KA, Williamson C. A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. AIDS Res Hum Retroviruses. 2005;21(1):98–102.
Smith JM. Analyzing the mosaic structure of genes. J Mol Evol. 1992;34(2):126–9.
Posada D, Crandall KA. Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc Natl Acad Sci USA. 2001;98(24):13757–62.
Gibbs MJ, Armstrong J, Gibbs A. Sister-Scanning: A Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics. 2000;16(7):573–82.
Boni MF, Posada D, Feldman MW. An exact nonparametric method for inferring mosaic structure in sequence triplets. Genetics. 2007;176(2):1035–47.
The authors thank Ass. Prof. Tomasz Stenzel from Department of Bird Diseases, Faculty of Veterinary Medicine in Olsztyn for proofreading the manuscript.
The authors declare that this study received funding from Xi'an IRIS Livestock Technology Co. Ltd. (K4030218170), China Postdoctoral Science Foundation (2017M610659 and 2018T111113), the State Key Laboratory of Veterinary Etiological Biology, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences (SKLVEB2016KFKT014) and the Fundamental Research Funds for the Central Universities (2452019053). The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
This study was performed according to the animal welfare guidelines of the World Organization for Animal Health, and animal sampling was strictly carried out in accordance with guidelines established by the Ethics Committee for Animal Experiments for Northwest A&F University, Yangling, China. The sampling was performed after a verbal informed consent of all racing pigeon club owners, primarily because: (i) There was no mandatory requirement for written document for sampling in the racing pigeon club based on the guidelines for animal experiments at our university; (ii) The sampling requirement was initially raised by the owners of racing pigeon club, in order to provide testing services for pigeon diseases. Besides, we also adherence to a best practice of veterinary care. All the protocols for animal experiments were approved by the Ethics Committee for Animal Experiments for Northwest A&F University, Yangling, China (Permit Number: 2014BAD23B11). The field study did not include endangered or protected species. No specific permissions were required for the collection of samples because the samples were collected from public areas or non-protected areas.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Pigeon circovirus (PiCV) references sequences obtained from GenBank and used in this study. The information including strain name, year, country, host, genome length and accession number.
Information of Pigeon circovirus (PiCV) stains obtained in this study. The information including strain name, collection date, collection region, gene length, initiation codons of cap gene, health status and accession number.
Sequence analysis of genome sequences of Pigeon circovirus (PiCV) strains. The information including homology analysis of amino acid and nucleotide sequences.
Details of recombination events (GenBank accession numbers) detected in the Pigeon circovirus (PiCV) strains obtained in this study. The information including location of recombination breakpoints, recombinant, potential minor parent(s), potential major parent(s), detection methods and p-value.
About this article
Cite this article
Wang, H., Gao, H., Jiang, Z. et al. Molecular detection and phylogenetic analysis of pigeon circovirus from racing pigeons in Northern China. BMC Genomics 23, 290 (2022). https://doi.org/10.1186/s12864-022-08425-8
- Genetic diversity
- Molecular epidemiology
- Phylogenetic analysis
- Pigeon circovirus
- Sequence analysis