- Research article
- Open Access
Pathogenicity and molecular characterization of a fowl adenovirus 4 isolated from chicken associated with IBH and HPS in China
BMC Veterinary Researchvolume 14, Article number: 400 (2018)
Since July in 2015, an emerging infectious disease, Fowl adenovirus (FAdV) species C infection with Hepatitis-Hydropericardium syndrome was prevalent in chicken flocks in China. In our study, one FAdV strain was isolated from commercial broiler chickens and was designated as SDSX1.The phylogenetic information, genetic mutations and pathogenicity of SDSX1 were evaluated.
The phylogenetic analysis indicated that SDSX1 is a strain of serotype 4, FAdV-C. The amino acid analysis of fiber-2 showed that there were more than 20 mutations compared with the non-virulent FAdV-C strains. The pathogenic evaluation of SDSX1 showed that the mortality of one-day-old chickens inoculated SDSX1 was 100%. The typical histopathological changes of SDSX1 were characterized by the presence of basophilic intranuclear inclusion bodies in hepatocytes. The virus copies in different tissues varied from107 to 1011 per 100 mg tissue and liver had the highest virus genome copies.
In conclusion, the isolate SDSX1, identified as FAdV-4, could cause one-day-old chicks’ typical inclusion body hepatitis (IBH) and hepatitis-hydropericardium syndrome (HHS) with 100% mortality. The virus genome loads were the highest in the liver. Molecular analysis indicated that substitutions in fiber-2 proteins may contribute to the pathogenicity of SDSX1.
Fowl adenoviruses (FAdVs) are non-enveloped, double-stranded DNA linear viruses, belonging to the genus Aviadenovirusof the family Adenoviridae [1, 2]. The viral genome is approximately 43–46 kb, encoding a number of structural and nonstructural proteins [3,4,5]. The virion capsid consists of three major structural proteins, hexon, fiber and penton base proteins. Hexon is the major protein of the adenovirus capsid known to have a region related to virus neutralization and serotype specificity.
FAdVs are grouped into 5 species (FAdV-A to FAdV-E) based on genomic restriction enzyme digest patterns with 12 serotypes identified so far: FAdV-A (serotype FAdV-1); FAdV-B (serotype FAdV-5);FAdV-C (serotypes FAdV-4 and -10); FAdV-D (serotypes FAdV-2, − 3, − 9, and − 11); and FAdV-E (serotypes FAdV-6, − 7, −8a, and -8b) . FAdVs have a worldwide distribution, but different serotypes or genotypes are discovered in different geographic regions [7,8,9,10,11]. FAdVs infection (the strain was closely related to serotype 11 and 8a) was first reported in Morocco recently . All 12 serotypes of FAdVs have been associated with the outbreak of inclusion body hepatitis (IBH) with a ~ 10% mortality rate . Chickens between 3 and 4 weeks old are especially vulnerable to IBH that is characterized by congested and enlarged liver with hepatocyte necrosis and petechial hemorrhage. As the pathognomonic lesion of FAdVs infection, basophilic intranuclear inclusion bodies are usually observed in degenerated hepatocytes surrounded by a clear halo or filling the entire nucleus . FAdV-4 has been implicated in hepatitis-hydropericardium syndrome (HHS) characterized by accumulation of transparent or straw-colored fluid in the pericardial sac, nephritis, and hepatitis with a 30–100% mortality rate [14,15,16]. Also, pathogenicity assays confirmed that FAdV-4 could induce HHS .
From June to November in 2015, a disease emerged in large-scale broiler and egg-laying chickens in farms in Shandong, Henan, Hebei, Liaoning, Anhui and Jilin provinces in China. The disease inflicted mostly the 20-to-30-day-old broilers with a mortality rate of 20–30%, and occasionally affected layers aged 20–70 days or 200–300 days but with a lower mortality rate. Usually, there were no obvious clinical signs with a sudden death. The death peak was 4–8 days after onset and the course of disease could last 8–15 days. The diseased animals had flabby hearts with straw-colored fluid in the pericardial sac and the swollen liver with foci of hemorrhage and/or necrosis. The aim of this study was to molecular characterization and pathogenicity assessment of FAdV strain responsible for disease outbreak in chickens.
The virus was isolated from the SPF chicken’s embryos by 5 blind passages and it was designated as SDSX1. The whole genome nucleotide sequence and hexon gene sequence of SDSX1 are available in the GenBank (accession number is KY636400 andKT932640, respectively) and the complete genome for SDSX1 was found to be 43,630 bp in length. The percent sequence identity for available whole genomes showed that SDSX1 has 100% sequence homology with three isolates from China (CH/SDDZ/2015, CH/AHBZ/2015, and CH/JSXZ/2015,  and 98.6% sequence identity with a KR5 isolate from Austria. Nevertheless, SDSX1 showed a low sequence identity (< 57.0%) with the strains of other FAdV species. According to the phylogenetic analysis based on the full genome (Fig. 1a), the SDSX1 isolate was classified into serotype 4 of FAdV-C. The similar evolutionary relationships were observed from the phylogenetic tree constructed based on the hexon gene (Fig. 1b).
The size of fiber-2 in SDSX1 isolate is 479 amino acid residues, the same as the strain KR5 (HE608152, FAdV-C), CH/SDDZ/2015(KU558761, FAdV-C), CH/AHBZ/2015 (KU569295, FAdV-C), and CH/JSXZ/2015(KU569296, FAdV-C) (Table 1) and 5 amino acid residues more than ON1(GU188428, FAdV-C) (Fig. 2). The basic residues rich sequence KRPK/KRAK (site 27–30) and VYPF (site 41–45) in fiber proteins of SDSX1were identified. Particularly, there were more than 20 substitutions compared with the non-pathogenic strains, ON1 and KR5 (Table 2 and Fig. 2) and shared only 94.1%~ 94.2% with ON1 and KR5, respectively. The fiber-2 protein identity of SDSX1 is 100% compared with CH/SDDZ/2015, CH/AHBZ/2015, and CH/JSXZ/2015.
Pathogenicity analysis of SDSX isolate
All the chickens inoculated with SDSX1 died within 5 days post inoculation and the survival curve was shown in Fig. 3. The median survival of inoculated group was 4 and the statistical differences between groups in survival was obvious significant (p < 0.0001). There were no obvious clinical manifestations for the inoculated chickens and the chickens had normal intake before death. Typical gross pathological changes in chickens were the flabby heart with straw-colored fluid in the pericardial sac and the swollen liver with foci of hemorrhages (Fig. 4a). The histopathology revealed that fatty degeneration of hepatocytes, basophilic intranuclear inclusion bodies in hepatocytes and necrosis of hepatocytes (Fig. 4b). All control chickens, the instead of that did not show any symptoms at the end of the observation.
Virus load analysis
Tissues including heart, liver, spleen, lung, kidney, and bursa of Fabricius were collected and the virus copies in tissues were determined using Real-time PCR. The results showed that virus can be detected in these tissues and the number of virus copies was between 107 and 1011 copies per 100 mg tissue with the highest viral genome copies in liver (Fig. 5).
Since 2015, FAdVs infection with high mortality was endemic in the poultry farms in China, resulting in considerable economic losses to the poultry industry [14, 18,19,20,21]. It is reported that most cases have been observed in broilers of 3 to 5 weeks of age, and occasionally in layers and breeder pullets aged 10 to 20 weeks and the infection was characterized by hepatitis hydropericardium syndrome (HHS) [10, 22]. In addition to the cases mentioned above, the isolation of FAdV-C from ducks with HHS was also reported [23, 24]. In our study, one strain was isolated from the chickens with HHS infection, and the whole genome was sequenced and analyzed. In addition, pathogenicity of the isolate SDSX1 was evaluated.
Phylogenetic analysis based on the FAdV isolates from different geographic areas worldwide showed that SDSX1 was clustered into serotype-4, FAdV-C species. Nucleotide alignment of SDSX1isolate with theidentity of reference strains in 5 species showed thehigher identity of 98.6–100% with FAdV-C using the complete genome sequence. The alignment of the hexon of SDSX1 had the similar result in the percentage identity compared to the full genome sequence. Especially, SDSX1 had 100% identity with the three isolates in China (CH/SDDZ/2015, CH/AHBZ/2015 and CH/JSXZ/2015).
The data indicated that SDSX1 was a highly pathogenic isolate. Regarding the pathogenicity of the SDSX1 isolate, the survival curve showed that all of 10 one-day-old SPF chickens inoculated with the SDSX1 died within 14 days. The results were consistent with several previous studies in which parts of serotype 4 were highly pathogenic strains [19, 25, 26]. The 100% of mortalityshowed in our study indicated that SDSX1 isolate was a virulent strain, and these were also accounted for the pathogenicity of the virus, serotypes, and dosage of infected strain, especially the susceptibility of the chickens. Different pathogenicity, dynamic distribution and replication pattern in chickens and ducks were observed in the recent report .
It is reported that fiber-2 played a significant role in virulence of FAdV . Several scientists’ studies showed that some amino acid residues, positions located at 219 and 380 in the fiber-2 protein of KR5, were related to the virulence based on the amino acid residues comparison between strains induced HHS and non-virulent strains [20, 28]. The fiber-2 protein not only plays a critical role in mediating infection by FAdV-4 but also has a potential capacity to be the candidate vaccine. Immunogenicity of the recombinant fiber-2 protein was evaluated and the results showed that both the level of IgY and the CD4+ T cell proliferative response of the chickens immunized with fiber-2 protein were significantly higher than that of the chickens immunized with inactivated vaccine, and fiber-2 could provide 100% protection . The neutralizing epitope recognized by a mAb 3C2 was identified to be located between aa 416–448 in the fiber-2 protein of FAdV-4 and the fiber-2 protein could efficiently block infection by FAdV-4 . In our study, there were more than 20 amino acids residues compared with non-virulent strains (ON1 and KR5), including the position 219(Glycine-Aspartic acid) and position 380(Alanine-Threonine) which had been identified. Whether the remaining amino acid resides substitutions of fiber-2 is responsible for the virulence of SDSX1 is unknown and need to be further studied in future.
In conclusion, the isolate SDSX1, identified as FAdV-4, could cause in one-day-old chicks’ typical inclusion body hepatitis (IBH) and hepatitis-hydropericardium syndrome (HHS) with 100% mortality. The virus genome loads were highest in the liver. Molecular analysis indicated that substitutions in fiber-2 proteins may contribute to the pathogenicity of SDSX1.
Samples used in this study
Twenty 10~20-day-olddead “817” broilers were collected from various farms in Liaocheng in Shandong province. Gross lesions were examined, and liver tissues were fixed in 10% formalin for 48 h at room temperature and then HE staining was performed following routine histological procedures. Animal work performed in this study was approved by Research Animal Use and Care Committee, the Agricultural University of Hebei, in compliance with the Guidance for the Care and Use of Laboratory Animals issued by the Ministry of Science and Technology of China.
Viral DNA extraction and PCR
A TaqMan-based real-time PCR assay developed in-house for specific detection of FAdV-4 was used to determine the presence of FAdV-4 in collected liver samples as described previously . The primers sequences were: sense, 5’-TTACGCTTACGGTGCCTACGT-3′;and antisense, 5’-CCGCGTTATTCATGATCCAGTA-3′. The TaqMan probe was as follows, FAM-5’CGACGGTTCCCAGTCCCTCACG-3′-Eclipse. The product size was 89 bp. Briefly, 10 g liver tissue was homogenized in 5 ml 0.01 M PBS. After freeze thawing for three times, the samples were centrifuged in 10000 rpm for 10 min. The supernatants of each sample were collected for viral DNA extraction. Viral DNA extraction was done TIANamp Virus DNA kit (Tiangen Biotech Co., Ltd., Beijing, China) according to the manufacturer’s instructions. The viral DNA from each sample was eluted in a final volume of 100 μl of nuclease-free water.
Next generation sequencing
The whole sequence of the isolated virus was sequenced by Next Generation Sequencing (NGS). Briefly, the DNA samples were ultrasonicated to generate fragments less than 500 bp. The DNA fragments were polished using T4 polynucleotide kinase. The adaptors were then ligated and loaded on the HiSeq 2000 for sequencing.
Virus isolation was performed on the samples with positive PCR results. The liver tissue was homogenized in 0.01 M PBS containing penicillin (100 IU/μL) and streptomycin (100 μg/μL) with the ration of 1:5(W/V), filtered and centrifuged. The supernatant from each sample was then collected and used for inoculation of 9-day-old specific pathogen free (SPF) chicken embryos (Beijing Vital River Laboratory Animal Technology Co., Ltd.) via allantoic cavity. The dose of inoculation was 0.1 ml per embryo. The allantoic fluid from infected embryos was harvested and used for repeated inoculation of chicken embryos. 5 embryos were used for each inoculum. The inoculation dose of allantoic fluid was 0.1 ml per embryo for each time. After 5 passages, allantoic fluid was harvested in the biohazard safety equipment and virus nucleic acid was extracted for DNA sequencing analysis. A specific fragment of 724 bp from the hexon gene was amplified by PCR for detecting the virus and the following primers were used: sense, 5’GACCTTCGCGGACTACTTGG3’; and antisense, 5’GACGGTTTGGTTGGAACTGAT3’.
To compare the complete genomic sequence and partial hexon gene sequence of the isolate determined in this study, 11 published genomic sequences available in GenBank (Table 1) were aligned by the MegAlign program of the DNAStar software version 7.0 with the Clustal W method (DNASTAR, Madison, Wisconsin, USA). The phylogenetic tree was constructed using the Maximum-Likelihood (ML) approach and MEGA6.06  on the aligned data set. The analysis of amino acids residues of fiber-2 of SDSX1 isolate was done using DNAStar software version 7.0.
Determination of EID50
Virus titer in chicken embryos for the allantoic fluid of five passage viruses was determined by chicken embryos infectivity assay. The result was recorded as EID50 per milliliter by using the Reed-Muench method. Briefly, 9-day-old chicken embryos were inoculated the allantoic fluid of the fifth passage viruses (0.1 ml/embryo), which were prepared by serial 10-fold dilution. The embryos were incubated at 37 °C for 6 days and observed every day. Embryos that died within the first 24 h were discarded, those that died than 24 h were selected for the allantoic fluid collection. PCR was used for detecting virus in allantoic fluid. Results were recorded.
Pathogenicity analysis of SDSX1 isolate
Twenty one-day SPF chickens were randomly allotted to two groups (10 chickens per group). For the inoculated group, each chicken was intraperitoneally inoculated with 1 ml SDSX1 virus containing 106 EID50 per chicken. Virus used for inoculation was prepared in the biohazard safety equipment. Chickens in the control group were mock inoculated with the same dosage of PBS. The dosage and pathway of inoculation was selected based on the pre-experiment (data not shown). All the chickens were observed daily for clinical signs up to all chickens are dead. The date and time of the death of each animal were recorded. When the chickens died, the liver samples were collected for the histopathological examination and virus detection.
Virus load analysis
Quantification of the SDSX1 load in each of the analyzed organs was performed by Real-time PCR method reported in the previous study . Briefly, an equal amount of heart, liver, spleen, lung, lung, kidney and bursa of Fabricicus (100 mg each) from each bird was homogenized in 2 ml PBS. Viral DNA extraction was done according to the method mentioned above. The viral DNA from each sample was eluted in a final volume of 100 μl of nuclease-free water. The TaqMan real-time PCR was performed in a final volume of 25 μl containing 2 × premix Ex Taq Master Mix12.5 μl, 1 μl of 10 μM primers, 0.8 μl of 10 μM probes, 1 μl DNA template and 8.7 μl RNase-free water. The reaction was carried out in ABI 7500 Real-Time PCR System and the optimal reaction procedure was as follows: pre-denaturing for 30s at 95 °C, followed by 40 cycles of denaturation at 95 °C for 10s, annealing at 65 °C for 35 s. A standard curve was used to calculate the virus copies: Y = − 3.56LogX + 46.08. “Y” indicates cycles and “X” indicates copy number.
Data were expressed as mean ± standard deviations. The survival curve of the inoculated and control chickens and the photograph of SDSX1 virus distribution in experimented SPF chickens were made using GraphPad Prism (version 5.0). One-way analysis of variance test was used.
- EID50 :
Egg 50% infectious dose
Phosphate buffer solution
Polymerase chain reaction
Adair B, Fitzgerald S. Group I adenovirus infections. In: Saif Y, Fadly A, Glisson J, McDougald L, Nolan L, Swayne D, editors. Diseases of poultry. 12th ed. Hoboken: Wiley-Blackwell; 2008. p. 252–66.
Benko M, Harrach B, Both G, Russell W, Adair B, A’ dam E’, de Jong J, Hess M, Johnson M, et al. In: Fauquet C, Mayo M, Maniloff J, Desselberger U, Ball L, editors. FamilyAdenoviridae. In Virus Taxonomy Eighth Report of the International Committee on the Taxonomy of Viruses. San Diego: Elsevier Academic Press; 2005. p. 213–28.
Griffin BD, Nagy E. Coding potential and transcript analysis of fowl adenovirus 4: insight into upstream ORFs as common sequence features in adenoviral transcripts. J Gen Virol. 2011;92(6):1260–72.
Shah MS, Ashraf A, Khan MI, Rahman M, Habib M, Chughtai MI, Qureshi JA. Fowl adenovirus: history, emergence, biology and development of a vaccine against hydropericardium syndrome. Arch Virol. 2017;162(7):1833–43.
Shah MA, Ullah R, March M, Shah MS, Ismat F, Habib M, Iqbal M, Onesti S, Rahman M. Overexpression and characterization of the 100K protein of fowl adenovirus-4 as an antiviral target. Virus Res. 2017;238:218–25.
Meulemans G, Couvreur B, Decaesstecker M, Boschmans M, van den Berg TP. Phylogenetic analysis of fowl adenoviruses. Avian pathol. 2004;33(2):164–70.
Mase M, Nakamura K, Minami F. Fowl adenoviruses isolated from chickens with inclusion body hepatitis in Japan, 2009-2010. J Vet Med Sci. 2012;74(8):1087–9.
Kim JN, Byun SH, Kim MJ, Kim JJ, Sung HW, Mo IP. Outbreaks of hydropericardium syndrome and molecular characterization of Korean fowl adenoviral isolates. Avian Dis. 2008;52(3):526–30.
Hess M, Raue R, Prusas C. Epidemiological studies on fowl adenoviruses isolated from cases of infectious hydropericardium. Avian pathol. 1999;28(5):433–9.
Zhao J, Zhong Q, Zhao Y, Hu YX, Zhang GZ. Pathogenicity and complete genome characterization of fowl adenoviruses isolated from chickens associated with inclusion body hepatitis and Hydropericardium syndrome in China. PLoS One. 2015;10(7):e0133073.
Alemnesh W, Hair-Bejo M, Aini I, Omar AR. Pathogenicity of fowl adenovirus in specific pathogen free chicken embryos. J Comp Pathol. 2012;146:223–9.
Redondo H, Fragoso JS, Tahala MA, Y B, Gil I, Elbachir E, Rodríguez MJ, Abad Moreno JC. Characterization of strain of fowl adenoviruses circulating in Morocco. Poult Sci. 2018;97(11):4057–62.
Nakamura K, Mase M, Yamaguchi S, Shibahara T, Yuasa N. Pathologic study of specific-pathogen-free chicks and hens inoculated with adenovirus isolated from hydropericardium syndrome. Avian Dis. 1999;43(3):414–23.
Lim TH, et al. Outbreak of gizzard erosion associated with fowl adenovirus infection in Korea. Poult Sci. 2012;91(5):1113–7.
Schachner A, Marek A, Jaskulska B, Bilic I, Hess M. Recombinant FAdV-4 fiber-2 protein protects chickens against hepatitis-hydropericardium syndrome (HHS). Vaccine. 2014;32(9):1086–92.
Ganesh K, Raghavan R. Hydropericardium hepatitis syndrome of broiler poultry: current status of research. Res Vet Sci. 2000;68(3):201–6.
Liu Y, Wan W, Gao D, Li Y, Yang X, Liu H, Yao H, Chen L, Wang C, Zhao J. Genetic characterization of novel fowl aviadenovirus 4 isolates from outbreaks of hepatitis-hydropericardium syndrome in broiler chickens in China. Emerg Microbes Infect. 2016;5(11):e117.
Niu YJ, Sun W, Zhang GH, Qu YJ, Wang PF, Sun HL, Wang L, Xiao YH, Liu SD. Corrigendum to Hydropericardium syndrome outbreak caused by fowl adenovirus serotype 4 in China in 2015. J Gen Virol. 2016;97(10):2797.
Li H, Wang J, Qiu L, Han Z, Liu S. Fowl adenovirus species C serotype 4 is attributed to the emergence of hepatitis-hydropericardium syndrome in chickens in China. Infect Genet Evol. 2016;45:230–41.
Vera-Hernández PF, Morales-Garzón A, Cortés-Espinosa DV, Galiote-Flores A, García-Barrera LJ, Rodríguez-Galindo ET, Toscano-Contreras A, Lucio-Decanini E, Absalón AE. Clinicopathological characterization and genomic sequence differences observed in a highly virulent fowl Aviadenovirus serotype 4. Avian Pathol. 2016;45(1):73–81.
Slaine PD, Ackford JG, Kropinski AM, Kozak RA, Krell PJ, Nagy É. Molecular characterization of pathogenic and nonpathogenic fowl aviadenovirus serotype 11 isolates. Can J Microbiol. 2016;62(12):993–1002.
Morshed R, Hosseini H, Langeroudi AG, MHB F, Charkhkar S. Fowl adenoviruses D and E cause inclusion body hepatitis outbreaks in broiler and broiler breeder pullet flocks. Avian Dis. 2017;61(2):205–10.
Chen H, Dou Y, Zheng X, Tang Y, Zhang M, Zhang Y, Wang Z, Diao Y. Hydropericardium hepatitis syndrome emerged in Cherry Valley ducks in China. Transbound Emerg Dis. 2017;64(4):1262–7.
Pan Q, Liu L, Wang Y, Zhang Y, Qi X, Liu C, Gao Y, Wang X, Cui H. The first whole genome sequence and pathogenicity characterization of a fowl adenovirus 4 isolated from ducks associated with inclusion body hepatitis and hydropericardium syndrome. Avian Pathol. 2017;46(5):571–8.
Pan Q, Yang Y, Shi Z, Liu L, Gao Y, Qi X, Liu C, Zhang Y, Cui H, Wang X. Different dynamic distribution in chickens and ducks of the Hypervirulent, novel genotype fowl adenovirus serotype 4 recently emerged in China. Front Microbiol. 2017;8:1005.
Niu YJ, Sun W, Zhang GH, Qu YJ, Wang PF, Sun HL, Xiao YH, Liu SD. Hydropericardium syndrome outbreak caused by fowl adenovirus serotype 4 in China in 2015. J Gen Virol. 2016;97(3):2684–90.
Pallister J, Wright PJ. SheppardM. A single gene encoding the fiber is responsible for variations in virulence in the fowl adenoviruses. J Virol. 1996;70(8):5115–22.
Marek A, Nolte V, Schachner A, Berger E, Schlötterer C, Hess M. Two fiber genes of nearly equal lengths are a common and distinctive feature of fowl adenovirus C members. Vet Microbiol. 2012;156(3–4):411–7.
Chen L, Yin L, Zhou Q, Li Q, Luo Y, Xu Z, Zhang Y, Xue C, Cao Y. Immunogenicity and protective efficacy of recombinant fiber-2 protein in protecting spf chickens against fowl adenovirus 4. Vaccine. 2018;36(9):1203–8.
Ping W, Zhang J, Wang W, Li T, Liang G, Shao H, Wei G, Qin A, Ye J. A novel monoclonal antibody efficiently blocks the infection of serotype 4 fowl adenovirus by targeting fiber-2. Vet Res. 2018;49(1):29.
Wang J, Wang J, Chen P, LiuL YW. Development of a TaqMan-based real-time PCR assay for rapid and specific detection of fowl aviadenovirus serotype 4. Avian Pathol. 2017;46(3):338–43.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.
This work was partly supported by the Earmarked Fund for Hebei Lays Innovation Team of Modern Agro-industry Technology Research System (HBCT2018150210), Science and Technology Project Foundation of Hebei Province (18227517D), and by grants for one-hundred outstanding innovative talents from Hebei institution of higher learning (SLRC2017039).
Availability of data and materials
The datasets generated during and analyzed during the current study are available from the corresponding author on reasonable request. Sequence data supporting the conclusions of this article are available in the GenBank repository under the accession numbers KY636400 andKT932640.
In this study, the samples were collected after the death of chickens due to the virus infection in the chicken farms. No chick was culled for sampling purpose in this study. Animal work performed in this study was approved by Research Animal Use and Care Committee, the Hebei Agricultural University, in compliance with the Guidance for the Care and Use of Laboratory Animals issued by the Ministry of Science and Technology of China.
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.