Molecular Characterization of Peruvian Fowl Adenovirus (FAdV) Isolates ()
1. Introduction
Fowl Adenovirus (FAdV) includes 12 serotypes [1] , with a high diversity of strains for each serotype and a wide range of clinical and pathological presentations (pneumonias, tracheitis, gizzard erosions, proventriculitis, pancreatitis, inclusion body hepatitis and hydropericardium) [2] . The most important sanitary problem in North America and South America is inclusion body hepatitis/hydropericardium syndrome due to high levels of morbidity and mortality in broilers resulting in significant economic losses [3] -[7] . This great diversity exhibited by FAdV required the use of additional techniques to make a correct identification of levels of species and serotypes and characterization of genetic diversity. An alternative to microbiologic and serologic techniques are molecular techniques based on DNA analysis. Several studies have reported the feasibility of molecular genetic tools for genetic differentiation and molecular typing of FAdV. The use of PCR for Hexon gene, RFLP analysis and direct sequencing of Hexon and Penton genes [8] -[10] have permitted a correct identification and characterization of clinical isolates. The lack of information available for molecular characterization of isolates from Peru do not allow the identification of the type of serotypes currently present in poultry production systems. This is a key requirement for establishing actions to control and prevent IBH. The objective of this study was to determine FAdV genotypes/serotypes isolated from IBH clinical cases in broilers from Peru.
2. Materials and Methods
2.1. Sources of Samples
A total of 47 suspected clinical cases of Inclusion body hepatitis (IBH) from broilers received in the Microbiology Laboratory-Bioservice SRL between July 2006 and April 2013 were included in this study (Table 1). Agar gel immune-diffusion (AGID) tests for serotypes 4, 8 and 9, liver histopathology analysis and Hexon gene PCR were performed to confirm the presence of fowl adenovirus (FAdV) in all suspected clinical cases.
2.2. Viral DNA Extraction and PCR Amplification of Hexon Gene
Viral DNA was extracted from liver lysate using DNeasy Tissue kit (Qiagen) according to manufacturer. PCR for L1 loop of the Hexon gene was carried out using primers HA and HB as described previously [9] . In brief, the PCR reaction was performed using a final volume of 20 uL containing 10 ng genomic DNA, 1X PCR buffer (20 mM Tris-HCl, 20 mM KCl, 5 mM (NH4)2SO4), 2 mM MgCl2, 5 pmol of each primer, 0.2 mM of dNTPs, 0.5 U Maxima HotStart-Taq DNA polymerase (Fermentas). Thermal cycles were the following: initial denaturation of 94˚C for 3 minutes, 30 cycles of 94˚C for 30 seconds, 53˚C for 60 seconds, 72˚C for 30 seconds and a final extension of 72˚C for 5 minutes. The PCR amplification was carried out using a Veriti Thermal cycler (Applied Biosystems). The PCR products were separated in 1% agarose gel TBE 1X electrophoresis and visualized by fluorescence using an ethidium bromide solution (0.1 ug/mL) and UV light.
2.3. Sequencing and Phylogenetics Analysis of L1 Loop of Hexon Gene
The PCR products were purified using the Wizard® SV Gel and PCR Clean-Up kit, cloned with pGEM-T Easy (Promega) and submitted for sequencing to Macrogen Inc, Korea. The sequencing reaction was performed in both directions using Big Dye sequencing kit (Applied Biosystems) and the primers HA and HB. Partial nucleotide sequence of the L1 Loop of hexon gene from FAdV’s was submitted to Gene Bank with the following accession number: KF601685.
Forward and reverse sequences were aligned together using Clustal program [11] to create consensus sequences. The consensus nucleotide sequences were aligned by the Clustal W method using MEGA 5.0 software [12] . Pairwise amino acid and nucleotide identity and pairwise evolutionary distances were computed using
Table 1 . Collection year and city of origin of 47 clinical cases of Inclusion body hepatitis (IBH) from broilers from Peru.
Mega v5.0. A phylogenetic analysis for the distance method was performed using partial sequences of L1 loop of Hexon gene and Hexon protein previously published in Gene Bank and FAdV Peruvian sequence (Table 2)
Table 2. ID sequence, accession number, country of origin, species, serotype and genetic cluster of 37 FAdV using phylogenetic analysis.
1Accesion number for Hexon protein or sequence in database Genbank or EMBL. 2Country of origin of strain. 3Fowl Adenovirus species according with [10] . 4Fowl Adenovirus serotype according with [1] . 5Fowl Adenovirus phylogenetic cluster according with [10] . 6No reported in [10] .
using Kimura-2-parameter method [13] and Dayhoff matrix based method [14] . A dendogram was constructed with the Neighbor-joining method [15] using 1000 bootstrap [16] replicate values with Mega v5.0.
3. Results
3.1. Clinical Cases of Inclusion Body Hepatitis
All clinical cases presented liver histopathology consistent with IBH lesions hepatocyte cytoplasmic vacuolation, multifocal to coalescing areas of hepatocytes necrosis, irregular shape hepatocyte nuclei, presence of basophilic intranuclear inclusion bodies without presence of HPS. The presence of FAdV was confirmed by PCR amplification of Hexon gene fragment in 47 suspected clinical cases of IBH.
3.2. Sequencing and Phylogenetic Analysis of L1 Loop of Hexon Gene
All 47 Peruvian FAdV isolates did not show genotype variation, a unique nucleotide sequence of 789 bp corresponding to position 273 and 1061 of Hexon gene amplified showed high nucleotide identity level (99.5% identity, E value = 0) and high conservation level in the amino acidic sequence (99.5%) in relation with a FAdV C serotype 4 previously identified [10] . The pairwise matrix of evolutionary divergence (nucleotide and amino acid) showed less divergence between Peruvian FAdV with FAdV C serotype 4 (0.005, 0.005) and FAdV C serotype US 9 (0.029, 0.070) in contrast to FAdV D serotype 8 (0.365, 0.709) and serotype 12 (0.378 and 0.663) (Table 3). Phylogenetic analysis from a segment of 770 nucleotide (260 amino acids) suggests the presence of FAdV group C serotype 4 present in all clinical cases from IBH of Peruvian broilers (Figure 1(a) and
(a)(b)
Figure 1. (a) Phylogenetic tree of L1 loop of the hexon nucleotide sequence from 12 species of FAdV’s. The evolutionary relationships of taxa were inferred using Neighbor-Joining method [17] . The evolutionary distances were computed using the Kimura 2-parameter method [13] . The analysis involved 20 nucleotide sequences and 770 nucleotide positions in the final dataset, the strain names are according to Table2 Bootstrap [16] test (1000 replicates) is expressed in % values. The scale bar indicates the evolutionary distance between sequences; (b) Phylogenetic tree of L1 loop of the hexon protein sequence from 12 species of FAdV’s. The evolutionary relationships of taxa were inferred using Neighbor-Joining method [17] . The evolutionary distances were computed using the Dayhoff matrix method [14] . The analysis involved 20 amino acid sequences and 199 amino acid positions in the final dataset, the strain names are according to Table2 Bootstrap [16] test (1000 replicates) is expressed in % values. The scale bar indicates the evolutionary distance between sequences.
Table 3. Pairwise estimates of evolutionary divergence1 of nucleotide sequence (below diagonal) and amino acid sequence (above diagonal) between 12 species of FAdV2.
1Analyses were conducted using p-distance [18] . This distance is the proportion (p) of amino acid (nucleotide) sites at which the two sequences to be compared are different. 2The FAdV strain accession number: FAdV1 (AAL13217), FAdV2 (AAN77072), FAdV3 (AAN77075), FAdV4 (AAN77077), FAdV5A (AAL13222), FAdV5B (AAN77079), FAdV5C (AAN77078), FAdV6A (AAL13224), FAdV6B (AAN77080), FAdV7A (AAN77081), FAdV7B (AAL13225), FAdV8A (AAN77082), FAdV8B (AAL13221), FAdV9A (AAL13226), FAdV9B (AAN77084), FAdV10 (AAL13227), FAdV11A (AAL13223), FAdV11B (AAN77085), FAdV12 (AAL13228) and FAdVPeru (this study).
Figure 1(b)). The FAdV genotype identified is closely related with the genotype C1 [19] previously reported in Peru (evolutionary divergence = 0, 100% nucleotide identity) and Ecuador (evolutionary divergence = 0.002, 99.9% nucleotide identity), however the two isolates showed 100% amino acid identity with FAdV isolated in this study (Table 4); this relationship is corroborated by phylogenetic clustering analysis (Figure 2(a) and Figure 2(b)).
4. Discussion
The nucleic acid technology has demonstrated the utility to detect differences in structural proteins related with immune response as Hexon protein. The sequencing of Hexon gene permits a correct identification of species level and serotype in FAdV, with similar results as reported for RFLP genome analysis [2] . Inside Hexon gene, 7 hypervariable regions were identified exclusively located in Loop 1 (L1) and Loop 2 (L2) regions [20] . The loop 1 in Hexon gene of FAdV reported higher variability [8] , and has demonstrated additional applications than include genetic variability analysis of FAdV isolates [6] [10] [21] .
(a)(b)
Figure 2. (a) Phylogenetic tree of L1 loop of the hexon nucleotide sequence from FAdV C serotype 4. The evolutionary relationships of taxa were inferred using Neighbor-Joining method [17] . The evolutionary distances were computed using the Kimura 2-parameter method [13] . The analysis involved 19 nucleotide sequences and 611 nucleotide positions in the final dataset, the strain names are according to Table2 Bootstrap [16] test (1000 replicates) is expressed in % values. The scale bar indicates the evolutionary distance between sequences. (b) Phylogenetic tree of L1 loop of the hexon protein sequence from FAdV C serotype 4. The evolutionary relationships of taxa were inferred using Neighbor-Joining method [17] . The evolutionary distances were computed using the Dayhoff matrix method [14] . The analysis involved 19 amino acid sequences and 202 amino acid positions in the final dataset, the strain names are according to Table2 Bootstrap [16] test (1000 replicates) is expressed in % values. The scale bar indicates the evolutionary distance between sequences.
During this study we analyzed 47 FAdV isolates from IBH clinical cases, examined their Hexon gene loop 1 sequences and compared them with reference strains previously published sequences. Our results suggest the absence of genetic variability in FAdV C in Peru, similar to previous reports [19] ; however this study only has access to one Peruvian isolate. Similar results were reported from India [22] [23] and South Korea [24] for clinical cases of IBH associated with the presence of FAdV C serotype 4 (>99% clinical cases of IBH). However, FAdV C serotype 4 (strain K531) is also associated with hydropericardium syndrome outbreaks (HPS) from South Korea [24] [25] and Germany [26] in contrast with Peruvian FAdV C serotype 4 than only reports IBH clinical lesions. The amino acid sequence analysis shows a slight differentiation between isolates from Peru and South Korea (Table 4) and phylogenetic analysis corroborates the different genetic clustering (Figure 2(B)).
The presence of multiple FAdV species in clinical cases of IBH were reported in Canada [6] with a fixation of FAdV serotypes 8 (45%) and 11 (45%) and South Korea [27] with FAdV serotype 4 (60%) and 11 (40%), however South Korean FAdV serotype 4 showed clinical hydropericardium/ IBH lesions in contrast to the Canadian report [6] .
This study is the first report of genetic characterization of FAdV C isolates from Peru that includes samples from a coastal region (Ica, La Libertad and Lima regions) which concentrates 80% of poultry production systems in Peru and provides important information than can serve as starting point for further investigations related to pathogenicity, antigenic properties and specific vaccine development for prevention and control of IBH
Table 4 . Pairwise estimates of evolutionary divergence1 of nucleotide sequences (below diagonal) and amino acid sequence (above diagonal) between 19 strain of FAdV-C serotype 4.
1Analyses were conducted using p-distance [18] . This distance is the proportion (p) of amino acid (nucleotide) sites at which the two sequences to be compared are different. 2The FAdV strain accession number: Germany A (FN869978.1), Austria A (FN869977.1), Pakistan (FN869976.1), Austria B (FN869975.1), Peru (FN869973.1), Ecuador (FN869972.1), Germany B (FN869971.1), Kuwait (FN869970.1), Mexico (FN869969.1), South Korea (HQ697593.1), Canada (EF685395.1), Russia A (AY581295.1), Italy A (HM592284.1), Italy B (HM592281.1), Italy C (HM592277.1), Italy D (HM592274.1), Indian A (AY581295.1), Indian B (EU847626.1) and Fad VPeru (this study).
associated with FAdV.
5. Conclusion
In conclusion, our results indicate the presence of FAdV C, serotype 4, genotype C1 circulating in clinical cases of IBH from Peruvian broilers.
Acknowledgements
This work was supported by Bioservice SRL. The authors wish to thanks to Dr. Karina Mendoza, Dr. Katherine Porturas, Dr. Stephane Lovon, Dr. Brenda Jara, Dr. Madeline Garcia, Dr. Heberht Uchuya and Dr. Jhonatan Castro for their cooperation in the collection and isolation of viral samples. We would like to also thank Dr. Manuel Moro for editorial comments and suggestions.
NOTES
*Corresponding author.