Translator Disclaimer
1 March 2008 Molecular Epidemiology of Chicken Anemia Virus (CAV) in Southeastern Chinese Live Birds Markets
Author Affiliations +
Abstract

Between January 2004 and December 2005, cloacal swabs from essentially healthy chickens and silky chickens from live birds markets in Guangdong and Hunan provinces in southeastern China were screened for chicken anemia virus (CAV) by polymerase chain reaction. Phylogenetic analysis of the major structural protein VP1 sequences showed no clear genotype cluster and no correlation with the geographic origin of CAV strains. Virus evolution at the amino acid level was very slow, which corresponds to a strong negative selection of the VP1 gene in China and worldwide. A high proportion (87%) of birds was CAV positive, suggesting that many farms in the region were infected. Further investigations are necessary to evaluate the economic losses caused by CAV and the cost–benefit of vaccination.

Chicken infectious anemia is caused by a virus that was first isolated in Japan in 1979 in specific-pathogen-free (SPF) chicks 27. Characteristic symptoms are aplastic anemia paired with hemorrhagic lesions. Immune suppression and generalized lymphoid atrophy are often associated with secondary viral, bacterial, or fungal infections 19. Direct mortality caused by chicken anemia virus (CAV) is usually relatively low but increased consumption of antibiotics and depressed weight gain represent an economic threat especially for the broiler industry and SPF egg producers.

CAV is the only member of the genus Gyrovirus (Circoviridae). It is a nonenveloped, icosahedral virus of about 25 nm in diameter with a negative-sense single-stranded circular DNA genome. The viral genome consists of 2.3 kilobases, with 3 partially overlapping open reading frames 13 for VP1, the major viral structural protein (51.6 kDa); VP2, a novel dual-specificity protein phosphatase 16 (24 kDa); and VP3, a nonstructural protein (13.6 kDa) named apoptin. VP1 and VP2 are the targets of neutralizing antibodies 19.

The VP1 gene is more variable than VP2 and VP3 5, and the worldwide genetic diversity of CAV VP1 is about 4.4% and 5.8% at nucleotide (nt) and amino acid (aa) level, respectively 5.

CAV has been found worldwide 19. In China, CAV was first isolated in 1996 from 25- to 40-day-old broilers 29. A survey in domestic poultry in 5 Chinese provinces (Beijing, Guangdong, Zhejiang, Shanghai, and Tianjin Shi) farms showed a 42% overall seroprevalence 28. So far, all viruses seemed to belong to the same serotype 19, although Spackman et al. 20,21 reported a putative second serotype. The emergence of additional new serotypes cannot be excluded and would have important consequences for vaccine efficacy and serodiagnosis. Although a few VP1 sequences of CAV from China are available on GenBank, no systematic analysis of Chinese strains has been performed to the best of our knowledge. Despite the high seroprevalence in the country 28 and the immunosuppression associated with an early infection, CAV vaccines have not been introduced in China.

Chicken farming provides an important source of cheap animal proteins. In China, 10.1 million metric tons of chicken meat and 28.7 million metric tons of eggs were produced in 2005 ( http://faostat.fao.org). As a result of severe acute respiratory syndrome (caused by a group II or IV coronavirus related to group III infectious bronchitis virus), and avian influenza virus outbreaks in southeastern China, the surveillance of avian viruses has been intensified, and sanitary conditions and control have been improved. Within the framework of the enhanced laboratory surveillance, we have characterized CAV on a molecular level in live bird markets in two provinces in southeastern China.

Materials and Methods

Field Samples

Cloacal swabs were collected from apparently healthy chickens (n  =  13) or silky chickens (n  =  94) at live birds markets in Guangdong (n  =  99) and Hunan (n  =  8) provinces between January 2004 and December 2005. Sellers on live poultry markets buy the birds from a wholesale market where both birds from large industrial farms and from familial farms are mixed. Meat chickens from large industrial farms are between 50 and 60 days of age, whereas chickens from familial farms are usually older than 4 mo. Old layers and breeders are also sold on these markets. Samples were stored at −80 C in virus transport medium (phosphate-buffered saline + 6 antibiotics) and had either never been thawed or thawed just once (for isolation attempt) before being processed for the present study.

DNA Isolation, Polymerase Chain Reaction (PCR), and Sequencing

DNA was extracted with the QIAmp viral RNA mini kit as recommended by the manufacturer (QIAGEN GmbH, Hilden, Germany;  http://www1.qiagen.com/literature/), it was screened for CAV genome, and the full VP1 gene of positive samples was amplified in a nested format for sensitivity reasons as described previously 5, with Takara Ex Taq® DNA polymerase (Takara Biotechnology, Dalian, China). A commercial live attenuated vaccine (Nobilis® CAV P4; Intervet B.V., Boxmeer, The Netherlands) was used as positive PCR control and to optimize the different PCRs. Sequencing was carried out in both directions with the Big Dye Terminator v.3.1 cycle sequencing kit (Applied Biosystems, Foster City, CA) on a capillary sequencer (model 3130xl; Applied Biosystems) using the PCR primers as sequencing primers. Nonredundant sequences were submitted to GenBank under the accession numbers AM407817 to AM407881. Strains were designated using the following nomenclature: World Health Organization country code_host code/year.month/province code/sample number; with CK for chicken, SCK for silky chicken; GD for Guangdong, and HN for Hunan province.

Data Analysis

Sequences were analyzed using the Bioedit program 7. Forward and reverse sequences were aligned with ClustalW 4. Phylogenetic and molecular evolutionary analyses were based on the entire VP1 gene sequences, using MEGA version 3.1 9. Phylogenic analysis of nucleic acid and deduced amino acid sequences was done with the neighbor-joining method, Kimura 2-parameter model. The aa sequences were also analyzed with the neighbor-joining method, with the Poisson correction. Bootstrap values (1000 replications) were indicated on each tree. Nucleotides were numbered according to Meehan et al. 11. All relevant complete VP1 sequences available on GenBank were used for comparison. Their accession numbers are indicated on Figs. 1 and 2.

Omega (dN/dS) values correspond to nucleotide mutation rate ratios of nonsynonymous and synonymous mutations. They were calculated for each group of sequences using PAML version 3.14 program 25http://abacus.gene.ucl.ac.uk/software/paml.html). Another test was carried out with PAML version 3.14 to detect individual codons that are under positive natural selection 12,26. The model (M0-M13) best adapted to the data set was selected with likelihood ratio tests with P < 0.01 1,18,26.

Results

Nucleic Acid Alignment and Phylogenic Analysis

Ninety-three of 107 (87%) cloacal swabs from southeastern China were positive for CAV in the nested detection PCR. Sixty-five complete sequences of VP1 gene were obtained from Guangdong (n  =  61) and Hunan (n  =  4) provinces. All sequences were 1350 nucleotides, and they had no insertions or deletions. Chinese VP1 sequences presented nucleotide mutations that had not been observed before in strains from other parts of the world, as, for example, the following 23 mutations, which occurred in at least 3 Chinese strains: positions 4, 94, 250, 420, 426, 444, 522, 555, 570, 807, 876, 903, 930, 1065, 1122, 1125, 1215, 1221, 1248, 1266, 1282, 1293, and 1308. The reliability of the sequences was ensured by forward and reverse sequencing of a template amplified with a high-fidelity polymerase with proofreading capability.

The maximum Kimura distance between the 74 Chinese VP1 sequences available was 5.2% (between AH4-China and BJ0401) in comparison with 6.2% (between Pallister-Australia and NI118-Nigeria) for the world maximum Kimura distance calculated for 112 relevant CAV VP1 sequences (neighbor-joining method, Kimura 2-parameters model). The new South Chinese sequences did not increase the latter diversity nor the diversity of the previous Chinese sequences. Phylogenetic analysis of CAV strains did not allow easy grouping of sequences into distinct clusters at the nt level (Fig. 1).

Amino Acid Alignment and Phylogenic Analysis

When the amino acid sequences of 112 VP1 genes were aligned, the maximal Kimura distance was 4.1%. The new Chinese sequences did not increase the overall maximal amino acid diversity of CAV.

Only 3 variable amino acids positions specific to Chinese strains were detected in more than one of the Chinese sequences, and 11 additional substitutions were detected in a single strain only (Fig. 2). None of these mutations were found in the hypervariable region (aa 139–151; Fig. 2) 17.

The phylogenetic tree at the aa level also showed no obvious CAV clusters supported by significant bootstrap values (Fig. 3).

Mean Kimura distances among the strains collected first (January 2004) and last (December 2005), and the mean distance between these 2 groups at the aa level, did not significantly differ (1%, 0.8%, and 0.7%, respectively). Mutations were, as expected, much more frequent in the “hypervariable” region (aa 139–151) (mutation probability: 62% or 8/13 positions) than in the whole VP1 aa sequence (mutation probability: 14% or 61/450 positions), confirming the hypervariable nature of this domain.

Analysis of Divergence and Selection

The CAV VP1 genes are highly conserved worldwide with an omega value (dN/dS) of 0.08, suggesting that the chance of a nonsynonymous mutation is only 8% of the chance of a synonymous mutation of becoming fixed in the population (Table 1). The calculation of omega ratios for China and for each continent demonstrated a negative (purifying) selection worldwide albeit with geographic differences. Omega values ranged from 0.04 (in China or in Nigeria) to 0.20 (in Europe). No single codon with positive (diversifying) selection could be detected in the Chinese strains.

Discussion

Although CAV has been detected in some parts of China, this is the first molecular epidemiologic analysis of this virus in the country. As the detection PCR was more sensitive than the sequencing PCR, only 70% of the CAV-positive cloacal swabs (65 of the 93 positive CAV strains) were suitable for sequencing. We cannot exclude that some strains may also have been missed because of mutations within primer locations. About 40% of the sequences showed double peaks with degenerated nt in the electropherograms, suggesting mixed infections as was observed and confirmed by cloning by us and others 5,23.

The nt or aa sequences revealed no clear phylogenetic clusters, or potential genotypes that would be supported by high bootstrap values (Figs. 1 and 3) irrespective of the method used (neighbor joining, maximum parsimony, minimum evolution, or unweighted pair group method with arithmetic mean; data not shown). Three clusters were previously identified by us and others 5,8 at the aa level: clusters I and III with the aa pattern 75V, 97M, 139K, 144E/K/N, and cluster II with the aa pattern 75I/T, 97L, 139Q, 144Q (Fig. 2). Maximum Kimura distances between (4.1% between 82-2Japan and Malaysia/AF285882) and within (3.4% and 2.0% between 82-2Japan and CH-BJ0401 and between Malaysia/AF285882 and CH-SH11, respectively) clusters were relatively low. Whether a 1.2-fold difference between inter- and intracluster maximum Kimura distances would be sufficient to define genotypes or whether current strains belong to a single genotype requires further attention.

Mutations associated with attenuation and a weaker reactivity with monoclonal antibody 2A9 (complete 75I, 89T, 125L, 141L, and 144E pattern 22; or 394H, 24) were absent in the Chinese strains.

CAV VP1 genes are worldwide under strong negative selection (omega  =  0.08; Table 1). In China and in Nigeria, even lower omega values were observed (0.04; Table 1), perhaps because of the absence of vaccination in both countries. However, omega calculations may be biased by the low number of sequences available on GenBank for most other countries and by a short collection period such as in this study. Omega values are also influenced by the number of passages of some strains such as Cux-1. The nt diversity of CAV is similar to the diversity of other DNA viruses, but very few of the mutated nt translate into aa. In addition to low virus survival rates with nonsilent mutations, other factors may slow down virus evolution such as its exquisite host adaptation 19, or the short generations in poultry.

In southeastern China most poultry from both family and commercial farms transit through wholesale markets, to retail live bird markets. CAV DNA was detected in 87% of all birds tested in this study. The birds were a cross-section of essentially healthy animals for sale at the participating live poultry markets in Guangdong and Hunan provinces, but it was not possible to trace the birds back to their farms. Although unlikely, it cannot be totally excluded that identical sequences from the same live poultry market are from the same farms (e.g., CH208-CH210 and CH4960-CH4964, collected on the same day, in the same live bird market; CH510-CH990 and CH3086-CH5262, collected within 3- and 6-wk intervals, at the same live bird markets). Therefore, our data suggest that the prevalence of CAV is very high in the region and that farms have a high risk of infection.

CAV DNA was detected for at least 40 wk after infection 2,3. After about 3 wk of age, immunocompetent chickens are resistant to disease, but they can acquire asymptomatic infections 15,19. The adult birds tested were probably subclinically infected. Because the virus was detected in cloacal swabs, all of these birds can contribute to the spread of the virus. Thus, the return of poultry or cages from markets to farms is likely to facilitate spill-back of virus to farms.

Previous studies suggested that coinfections with CAV and infectious bursal disease virus (IBDV) may result in increased mortality 8, but IBDV PCRs 14 were negative in all swabs (data not shown), probably because of widespread vaccination in China.

Although the high proportion of the virus suggests that young chickens may be partially protected by maternally derived antibodies, further studies are required to quantify losses in productivity caused by CAV-related immunosuppression and secondary infections by other pathogens. A complementary study at a farm rather than a market level would be warranted to further evaluate the prevalence of CAV in the region, the role of live bird markets in the spread of CAV, and the cost–benefit of vaccination. Wholesale and live birds markets may have a grossly underestimated role in the spread, evolution, or both of avian pathogens 10, as suggested previously for infectious bronchitis virus 6.

Acknowledgments

We thank Y. Huang and Z. Zhen for technical help. The financial support of the Ministry of Cooperation, Luxembourg, and of the Li Ka Shing Foundation, China, are gratefully acknowledged. M.F.D. was supported by a BFR fellowship of the Ministry of Research and Higher Education, Luxembourg, and was a visiting student with the Joint Influenza Research Centre of SUMC/HKU.

References

1.

M. Anisimova, R. Nielsen, and Z. Yang . Effect of recombination on the accuracy of the likelihood method for detecting positive selection at amino acid sites. Genetics 164:1229–1236.2003.  Google Scholar

2.

C. Cardona, B. Lucio, P. O'Connell, J. Jagne, and K. A. Schat . Humoral immune responses to chicken infectious anemia virus in three strains of chickens in a closed flock. Avian Dis 44:661–667.2000.  Google Scholar

3.

C. J. Cardona, W. B. Oswald, and K. A. Schat . Distribution of chicken anaemia virus in the reproductive tissues of specific-pathogen-free chickens. J. Gen. Virol 81:2067–2075.2000.  Google Scholar

4.

R. Chenna, H. Sugawara, T. Koike, R. Lopez, T. J. Gibson, D. G. Higgins, and J. D. Thompson . Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res 31:3497–3500.2003.  Google Scholar

5.

M. F. Ducatez, A. A. Owoade, J. O. Abiola, and C. P. Muller . Molecular epidemiology of chicken anemia virus in Nigeria. Arch. Virol 151:97–111.2006.  Google Scholar

6.

M. F. Ducatez, Y. Guan, H. Chen, and C. P. Muller . Enhanced laboratory surveillance of group III coronaviruses in life poultry markets in Guangdong province, China, after the SARS outbreak. In:. Proc. 5th International Symposium on Avian Corona- and Pneumovirus, Rauischholzhausen, Germany 132–140.2006.  Google Scholar

7.

T. A. Hall BioEdit: a user-friendly biological sequence alignment editor and analysis for Windows 95/98/NT. Nucleic Acids Symp. Ser 41:95–98.1999.  Google Scholar

8.

M. R. Islam, R. Johne, R. Raue, D. Todd, and H. Muller . Sequence analysis of the full-length cloned DNA of a chicken anaemia virus (CAV) strain from Bangladesh: evidence for genetic grouping of CAV strains based on the deduced VP1 amino acid sequences. J. Vet. Med. B Infect. Dis. Vet. Public Health 49:332–337.2002.  Google Scholar

9.

S. Kumar, K. Tamura, I. B. Jakobsen, and M. Nei . MEGA2: molecular evolutionary genetics analysis software. Bioinformatics 17:1244–1245.2001.  Google Scholar

10.

M. Liu, S. He, D. Walker, N. Zhou, D. R. Perez, B. Mo, F. Li, X. Huang, R. G. Webster, and R. J. Webby . The influenza virus gene pool in a poultry market in South central china. Virology 305:267–275.2003.  Google Scholar

11.

B. M. Meehan, D. Todd, J. L. Creelan, J. A. Earle, E. M. Hoey, and M. S. McNulty . Characterization of viral DNAs from cells infected with chicken anaemia agent: sequence analysis of the cloned replicative form and transfection capabilities of cloned genome fragments. Arch. Virol 124:301–319.1992.  Google Scholar

12.

R. Nielsen and Z. Yang . Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics 148:929–936.1998.  Google Scholar

13.

M. H. Noteborn, G. F. de Boer, D. J. van Roozelaar, C. Karreman, O. Kranenburg, J. G. Vos, S. H. Jeurissen, R. C. Hoeben, A. Zantema, G. Koch, H. van Ormont, and A. J. van der Eb . Characterization of cloned chicken anemia virus DNA that contains all elements for the infectious replication cycle. J. Virol 65:3131–3139.1991.  Google Scholar

14.

A. A. Owoade, M. N. Mulders, J. Kohnen, W. Ammerlaan, and C. P. Muller . High sequence diversity in infectious bursal disease virus serotype 1 in poultry and turkey suggests West-African origin of very virulent strains. Arch. Virol 149:653–672.2004.  Google Scholar

15.

A. A. Owoade, D. O. Oluwayelu, O. A. Fagbohun, W. Ammerlaan, M. N. Mulders, and C. P. Muller . Serologic evidence of chicken infectious anemia in commercial chicken flocks in southwest Nigeria. Avian Dis 48:202–205.2004.  Google Scholar

16.

M. A. Peters, D. C. Jackson, B. S. Crabb, and G. F. Browning . Chicken anemia virus VP2 is a novel dual specificity protein phosphatase. J. Biol. Chem 277:39566–39573.2002.  Google Scholar

17.

R. W. Renshaw, C. Soine, T. Weinkle, P. H. O'Connell, K. Ohashi, S. Watson, B. Lucio, S. Harrington, and K. A. Schat . A hypervariable region in VP1 of chicken infectious anemia virus mediates rate of spread and cell tropism in tissue culture. J. Virol 70:8872–8878.1996.  Google Scholar

18.

R. Sainudiin, W. S. Wong, K. Yogeeswaran, J. B. Nasrallah, Z. Yang, and R. Nielsen . Detecting site-specific physicochemical selective pressures: applications to the class I HLA of the human major histocompatibility complex and the SRK of the plant sporophytic self-incompatibility system. J. Mol. Evol 60:315–326.2005.  Google Scholar

19.

K. A. Schat Chicken infectious anemia. In: Diseases of poultry. Y.M. Saif, H.J. Barnes, J.R. Glisson, A.M. Fadly, L.R. McDougald, and D.E. Swayne, eds. Ames, IA Iowa State University Press. 182–202.2003.  Google Scholar

20.

E. Spackman, S. S. Cloud, C. R. Pope, and J. K. Rosenberger . Comparison of a putative second serotype of chicken infectious anemia virus with a prototypical isolate I. Pathogenesis. Avian Dis 46:945–955.2002.  Google Scholar

21.

E. Spackman, S. S. Cloud, and J. K. Rosenberger . Comparison of a putative second serotype of chicken infectious anemia virus with a prototypical isolate II. Antigenic and physicochemical characteristics. Avian Dis 46:956–963.2002.  Google Scholar

22.

D. Todd, A. N. Scott, N. W. Ball, B. J. Borghmans, and B. M. Adair . Molecular basis of the attenuation exhibited by molecularly cloned highly passaged chicken anemia virus isolates. J. Virol 76:8472–8474.2002.  Google Scholar

23.

V. L. van Santen, L. Li, F. J. Hoerr, and L. H. Lauerman . Genetic characterization of chicken anemia virus from commercial broiler chickens in Alabama. Avian Dis 45:2373–388.2001.  Google Scholar

24.

S. Yamaguchi, T. Imada, N. Kaji, M. Mase, K. Tsukamoto, N. Tanimura, and N. Yuasa . Identification of a genetic determinant of pathogenicity in chicken anaemia virus. J. Gen. Virol 82:1233–1238.2001.  Google Scholar

25.

Z. Yang PAML: a program package for phylogenetic analysis by maximum likelihood. Comput. Appl. Biosci 13:555–556.1997.  Google Scholar

26.

Z. Yang, R. Nielsen, N. Goldman, and A. M. Pedersen . Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics 155:431–449.2000.  Google Scholar

27.

N. Yuasa Isolation and some characteristics of an agent inducing anemia in chicks. Avian Dis 23:366–385.1979.  Google Scholar

28.

W. Zhou, B. Yang, B. Shen, S. Han, and J. Zhou . A serologic survey of antibody against chicken infectious anemia virus by indirect immunofluorescent assay in domestic poultry in China. Avian Dis 40:358–360.1996.  Google Scholar

29.

W. Zhou, B. Shen, B. Yang, S. Han, L. Wei, B. Xiao, and J. Zhou . Isolation and identification of chicken infectious anemia virus in China. Avian Dis 41:361–364.1997.  Google Scholar

Appendices

Fig. 1.

Phylogenic analysis of the nucleic acid sequence of 26 new complete CAV VP1 sequences from Guangdong and Huan provinces, China, and 47 relevant complete VP1 sequences currently available on GenBank. The Nobilis® P4 vaccine strain is also shown. Numbers at nodes correspond to bootstrap values >49. Chinese strains from the present study are marked with a closed triangle, whereas other Chinese sequences from GenBank have open triangles. Chinese strains with identical aa sequences are shown only once. Mixed infected samples sequences are not included in the phylogenetic analysis.

i0005-2086-52-1-68-f01.gif

Fig. 2.

Amino acid substitutions within VP1 protein of both cluster I and II strains. The hypervariable region 16 and the aa substitutions specific to Chinese strains are shown. Specific Chinese aa mutations: italic for mutations present in >1 Chinese strains; normal font for mutations present in a single Chinese strain; and bold: aa pattern of previously described clusters I, II, and III 5,8. 1 Most frequent aa. 2 aa substitution in CH683, CH5256, CH2768, CH4665, CH22035, CH4663, and CH 5840. 3 aa substitution in CH210, CH208, and CH7386. 4 aa substitution in CH-TJBD33 and CH-BJ0401. 5 aa substitution in CH5264. 6 aa substitution in CH4665. 7 aa substitution in CH-AY839944. 8 aa substitution in CH-BJ0401. 9 aa substitution in CH-SD24. 10 aa substitution in CH-HN9. 11 aa substitution in CH-AH4. 12 aa substitution in CH656. 13 aa substitution in CH-SH11.

i0005-2086-52-1-68-f02.gif

Fig. 3.

Phylogenic analysis of the predicted aa sequences of the 74 CAV strains from Fig. 1. Symbols are as described in Fig. 1. Sequences of mixed infected samples are not included.

i0005-2086-52-1-68-f03.gif

Table 1.

Mean genetic distances and omega values for CAV VP1 genes per country or continent.

i0005-2086-52-1-68-t01.gif

[1] D. Corresponding author: 20A, rue Auguste Lumière, L-1950 Luxembourg, Grand-Duchy of Luxembourg. claude.muller@LNS.ETAT.LU

M. F. Ducatez, H. Chen, Y. Guan, and C. P. Muller "Molecular Epidemiology of Chicken Anemia Virus (CAV) in Southeastern Chinese Live Birds Markets," Avian Diseases 52(1), 68-73, (1 March 2008). https://doi.org/10.1637/8049-070407-Reg
Received: 6 July 2007; Accepted: 1 October 2007; Published: 1 March 2008
JOURNAL ARTICLE
6 PAGES


SHARE
ARTICLE IMPACT
Back to Top