Analysis of chicken anemia virus genome: evidence of intersubtype recombination
© Eltahir et al; licensee BioMed Central Ltd. 2011
Received: 20 April 2011
Accepted: 10 November 2011
Published: 10 November 2011
Chicken anemia virus (CAV) is the causative agent of chicken infectious anemia. CAV putative intergenotypic recombinants have been reported previously. This fact is based on the previous classification of CAV sequences into three genotypes. However, it is unknown whether intersubtype recombination occurs between the recently reported four CAV genotypes and five subtypes of genome sequences.
Phylogenetic analysis, together with a variety of computational recombination detection algorithms, was used to investigate CAV approximately full genomes. Statistically significant evidence of intersubtype recombination was detected in the parent-like and two putative CAV recombinant sequences. This event was shown to occur between CAV subgroup A1 and A2 sequences in the phylogenetic trees.
We revealed that intersubtype recombination in CAV genome sequences played a role in generating genetic diversity within the natural population of CAV.
Chicken anemia virus (CAV) was first reported in 1979 in specific-pathogen-free (SPF) chickens . CAV belongs to the Circoviridae and is a non-enveloped, icosahedral virus with a negative-sense, single-stranded circular DNA. The viral genome consists of 2.3 kb, with three partially overlapping open reading frames. CAV infection is an economically important clinical and subclinical disease in broiler chickens, with a worldwide distribution .
CAV isolates show extremely limited genetic variability worldwide . All isolates of CAV are suspected to belong to a single serotype . Little is known about CAV genome recombination analysis. CAV putative intergenotype recombinants have been reported to occurs in the virus gene VP1 and results in a new virus genotype .
Recently, after adding more CAV approximately full genome sequences to GenBank, CAV sequences arising from different parts of world have been categorized into four genotype groups (A-D) and five subtypes (A1, A2, A3, D1 and D2) . Therefore, the necessity of exploring these CAV genotypes for evidence of recombination as an important tool for genetic variability has been raised, to establish if any would help in understanding the evolutionary process in the CAV genome.
Here, we report evidence of intersubtype recombination based on sequence analysis of the entire coding regions (VP1, VP2 and VP3) of CAV genomes.
DNA extraction, PCR screening, amplification of CAV genome, cloning and sequencing were carried out as previously described, briefly, primers VP1F: 5'AGCCGACCCCGAACCGCAAGAA'3 and VP1R: 5' TCA GGG CTG CGT CCC CCA GTA CA'3 were used to amplify the VP1 region and primers VP2F: 5' GCG CAC ATA CCG GTC GGC AGT'3 and VP2R: 5' GGG GTT CGG CAG CCT CAC ACT AT'3 were used to amplify the VP2 region  for 10 spleen samples collected between April 2010 and December 2010 at Yangzhou University Veterinary Hospital during necropsy. Samples originated from different commercial chicken farms in Anhui (n = 3) and Jiangsu (n = 7) provinces. For each samples, DNA extraction and PCR was run at least twice. Animal experiments were conducted in accordance with the guidelines provided by the Chinese Council on Animal Care. All experiments complied with institutional animal care guidelines and were approved by University of Yangzhou Animal Care Committee (protocol number 06R015).
To look for recombination in CAV sequences, we used 55 published full genomes sequences that are currently available in GenBank, together with 10 sequences characterized in the present study. Multiple alignments of either CAV full or subgenomic regions were performed and analyzed using ClustalW . Unrooted phylogenetic trees were constructed using the neighbor-joining (NJ) method, and visualized and edited using MEGA 3.1 software . The evolutionary distances were estimated using the Kimura two-parameter method. Bootstrap analyses were performed with 1,000 repeat samples of the data sets. To reduce redundancy, isolates with previous reports of intergenotypic recombination  were excluded from analysis.
Identification of recombination
Recombination breakpoint events in the multiple alignments were detected with the Recombination Detection Program 3 (RDP3), using the automated suite of algorithms implemented in the RDP3 with default settings. These included RDP, GENECONV, BootScan, MaxChi, Chimaera, SiScan, Phyl- Pro, LARD, and 3Seq [9, 10]. For the RDP algorithm, the reference sequence parameter (internal and external reference) was used as recommended by the manual. The multiple-aligned sequences were cut into sections at the positions where the putative recombination breakpoints were identified, and phylogenetic analyses were carried out by MEGA 3.1 software . Confirmation of recombination break points was further carried out by bootscanning using SimPlot .
Sequences accession numbers
CAV genome sequences characterized in this study were submitted to GenBank with accession numbers [GenBank: FR 850021, GenBank: FR 850022, GenBank: FR 850023, GenBank: FR 850024, GenBank: FR 850025, GenBank: FR 850026, GenBank: FR 850027, GenBank: FR 850028, GenBank: FR 850029, GenBank: FR 850030].
Phylogenetic tree reconstructions
4.544 × 10-02
1.179 × 10-02
4.692 × 10-03
1.768 × 10-04
6.215 × 10-07
1.499 × 10-05
1.612 × 10-05
Confirmation of recombination event
Recombination is one of the evolutionary processes that shape the architecture of viral genomes. Ignoring the occurrence of recombination may influence the analysis of genetic data and the conclusions derived from it. Recombination analysis of CAV genomes is very limited. CAV intergenotype recombination break points in the VP1 region have been reported in only a single study using 35 CAV full genomes sequences . However, to date, there is no evidence of inter- or intrasubtype recombination between CAV genomes, to the best of our knowledge. Therefore, the present study investigated CAV whole genome sequences to explore the effect of recombination in the genetic variation of CAV sequences.
Here, we report the first evidence of intersubtype recombination events in CAV genome sequences amplified from two different areas (Anhui and Jiangsu provinces) in China. Intersubtype mosaic viruses are the results of co-infection events in geographic regions where multiple subtypes circulate [12, 13]. Similarly, co-infection with CAV that has been documented by sequencing electropherogram double peaks and cloning has also been reported . Moreover, intersubtype recombination in other virus subtypes has also been reported [15, 16].
Natural recombination in the VP1 of the CAV genome, with a resultant new genotype, has suggested speeding up of CAV evolution . However, it should be noted that these findings were based on a limited number of CAV sequences and the previous categorization of CAV sequences into three genotypes. In contrast, despite adding more sequences here, we confirmed that CAV genomes can be categorized into their recent four genotype groups (A-D) and five subtypes (A1, A2, A3, D1 and D2) .
The recombination breakpoints detected in CAV subgroup A1 and A2 sequences in the present study were also found in the VP1 protein. This resulted in the sequence from subgroup A2 being located within subgroup A1. It is known that most nucleotide variation in the CAV genome is observed in VP1, in addition to the fact that virulence of CAV is also determined by VP1 [17, 18]. The cumulative evidence of recombination breakpoints in the VP1 of the CAV genome may be indicative of favorable selection for recombination in this variable region. Sequence variations at specific sites between HIV-1 isolates have been reported to introduce unique recombination hotspots and increase recombination frequencies in the C3/V4 env region . However, whether there is correlation between the highest genetic variability in CAV genomes observed in the VP1 and recombination still needs to be investigated.
Using bioinformatics analysis, we were able to detect only a single incidence of CAV intersubtype recombination that contributed to CAV genetic variation. More detailed analysis for further CAV intragenotypic recombination is still required.
The research was supported by the Program for Changjiang Scholars and Innovative Research Team in University (IRT0978) and the Major Basic Research of Natural Science Foundation of the Jiangsu Higher Education Institutions of China (Grant No. 07KJA23021).
- Yuasa N, Taniguchi T, Yoshida I: Isolation and some characteristics of an agent inducing anemia in chicks. Avian Dis 1979, 23: 366-385. 10.2307/1589567View ArticleGoogle Scholar
- Todd D, Creelan JL, Mackie DP, Rixon F, McNulty MS: Purification and biochemical characterization of chicken anemia agent. J Gen Virol 1990, 71: 819-823. 10.1099/0022-1317-71-4-819View ArticlePubMedGoogle Scholar
- Van Santen VL, Li L, Hoerr FJ, Lauerman LH: Genetic characterization of chicken anemia virus from commercial broiler chickens in Alabama. Avian Dis 2001,45(2):373-388. 10.2307/1592977View ArticlePubMedGoogle Scholar
- Koch G, van Roozelaar DJ, Verschueren CA, van der Eb AJ, Noteborn MH: Immunogenic and protective properties of chicken anaemia virus proteins expressed by baculovirus. Vaccine 1995,13(8):763-770. 10.1016/0264-410X(94)00034-KView ArticlePubMedGoogle Scholar
- He CQ, Ding NZ, Fan W, Wu YH, Li JP, Li YL: Identification of chicken anemia virus putative intergenotype recombinants. Virology 2007,15(366(1)):1-7.View ArticleGoogle Scholar
- Eltahir YM, Qian K, Jin W, Wang P, Qin A: Molecular epidemiology of chicken anemia virus in commercial farms in China. Virol J 2011, 145.Google Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 1997, 25:(24), 4876-4882.PubMed CentralView ArticlePubMedGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform 2004, 5: 150-163. 10.1093/bib/5.2.150View ArticlePubMedGoogle Scholar
- Martin DP, Williamson C, Posada D: RDP2: recombination detection and analysis from sequence alignments. Bioinformatics 2005, 21: 260-262. 10.1093/bioinformatics/bth490View ArticlePubMedGoogle Scholar
- 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: 98-102. 10.1089/aid.2005.21.98View ArticlePubMedGoogle Scholar
- Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC: Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol 1999, 73: 152-160.PubMed CentralPubMedGoogle Scholar
- Salminen MO, Carr JK, Robertson DL: Evolution and probable transmission of intersubtype recombinant human immunodeficiency virus type 1 in a Zambian couple. J Virol 1997, 71: 2647-2655.PubMed CentralPubMedGoogle Scholar
- Gao F, Robertson DL, Carruthers CD: A comprehensive panel of near-full-length clones and reference sequences for nonsubtype B isolates of human immunodeficiency virus type 1. J Virol 1998, 72: 5680-5698.PubMed CentralPubMedGoogle Scholar
- Ducatez MF, Owoade AA, Abiola JO, Muller CP: Molecular epidemiology of chicken anemia virus in Nigeria. Arch Virol 2006, 151: 97-111. 10.1007/s00705-005-0609-7View ArticlePubMedGoogle Scholar
- Fomsgaard A, Vinner L, Therrien D, Jørgensen LB, Nielsen C, Mathiesen L, Pedersen C, Corbet S: Full-length characterization of A1/D intersubtype recombinant genomes from a therapy-induced HIV type 1 controller during acute infection and his noncontrolling partner. AIDS Res Hum Retroviruses 2008,24(3):463-72. 10.1089/aid.2006.0294View ArticlePubMedGoogle Scholar
- Kurbanov F, Tanaka Y, Fujiwara K, Sugauchi F, Mbanya D, Zekeng L, Ndembi N, Ngansop C, Kaptue L, Miura T, Ido E, Hayami M, Ichimura H, Mizokami M: A new subtype (subgenotype) Ac (A3) of hepatitis B virus and recombination between genotypes A and E in Cameroon. J Gen Virol 2005, 86: 2047-56. 10.1099/vir.0.80922-0View ArticlePubMedGoogle Scholar
- Todd D, Scott AN, Ball NW, Borghmans BJ, Adair BM: Molecular basis of the attenuation exhibited by molecularly cloned highly passaged chicken anemia virus isolates. J Virol 2002, 76: 8472-8474. 10.1128/JVI.76.16.8472-8474.2002PubMed CentralView ArticlePubMedGoogle Scholar
- Yamaguchi S, Kaji N, Mase M, Tsukamoto K, Tanimura N, Yuasa N: Identification of a genetic determinant of pathogenicity in chicken anaemia virus. J Gen Virol 2001, 82: 1233-1238.View ArticlePubMedGoogle Scholar
- Baird HA, Gao Y, Galetto R, Lalonde M, Anthony RM, Giacomoni V, Abreha M, Destefano JJ, Negroni M, Arts EJ: Influence of sequence identity and unique breakpoints on the frequency of intersubtype HIV-1 recombination. Retrovirology 2006, 12: 3-91.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.