Identification of D/Yama2019 Lineage-Like Influenza D Virus in Chinese Cattle
Jieshi Yu1,2† 
Tianyu Li3,4† 
Zhenyu Wen5† 
Siyu Wu1 Zhilin Wang1 Jiaying Zheng4 Mingwang Chen4 Faming Chen4 Wen-Kang Wei1* 
Shao-Lun Zhai2,6* 
Ming Liao2,5,6*- 1Agro-biological Gene Research Center, Guangdong Academy of Agricultural Sciences, Guangzhou, China
- 2Key Laboratory for Prevention and Control of Avian Influenza and Other Major Poultry Diseases, Ministry of Agriculture and Rural Affairs, Scientific Observation and Experiment Station of Veterinary Drugs and Diagnostic Techniques of Guangdong Province, Ministry of Agriculture and Rural Affairs, Key Laboratory of Livestock Disease Prevention of Guangdong Province, Institute of Animal Health, Guangdong Academy of Agricultural Sciences, Guangzhou, China
- 3College of Animal Science and Technology, Guangxi University, Nanning, China
- 4Zhongshan Animal Disease Control Center, Zhongshan, China
- 5Key Laboratory of Zoonosis of Ministry of Agriculture and Rural Affairs, Key Laboratory of Zoonosis Prevention and Control of Guangdong Province, South China Agricultural University, Guangzhou, China
- 6Maoming Branch Center of Guangdong Laboratory for Lingnan Modern Agricultural Science and Technology, Maoming, China
Outbreaks of influenza D virus (IDV) continue to be reported in many countries. On the basis of the hemagglutinin-esterase fusion (HEF) gene, five IDV genetic lineages have been identified: D/OK, D/660, D/Yama2016, D/Yama2019 and D/CA2019 lineages. Previously reported IDV strains in China all form a sub-clade (D/China sub-lineage) within D/OK lineage. From October 2021 to February 2022, nasal swab samples (n = 250) were collected from apparently healthy cattle in slaughterhouses around the city of Guangzhou, China, and screened for IDV by RT-PCR. Ten samples were positive for IDV. An IDV strain with nearly complete genome sequences was identified and designated as D/bovine/CHN/JY3001/2021. Importantly, sequence alignments and phylogenetic analyses revealed that this IDV strain is genetically close to the strains (>98% homology) in the D/Yama2019 lineage that has been found only in Japan, but distant from the previously reported Chinese IDV strains (~95% similarity). These results demonstrate the emergence of D/Yama2019 lineage IDV in Chinese cattle herds, highlighting a need for future surveillance of D/Yama2019-like viruses toward better understanding both epidemiology and diversity of IDV in China.
Introduction
Influenza D virus (IDV) was first reported in 2011 and has been found in a wide range of animal species all over the world (1, 2). Unlike influenza A, B and C viruses, cattle are the major reservoir of IDV (3). Increasing studies have shown that IDV can be highly associated with bovine respiratory disease (RBD) complex (1, 4) that causes huge economic losses in cattle industry.
Influenza D viruses (IDVs) can be divided into five genetically distinct lineages (D/OK, D/660, D/Yama2016, D/Yama2019 and D/CA2019) based on the hemagglutinin-esterase fusion (HEF) gene (5). D/OK lineage is so far the most frequently reported lineage worldwide, while D/660 lineage is mainly distributed in Europe and North America. Interestingly, both D/Yama2016 and D/Yama2019 lineages seem to be restricted in Japan (1, 6–8). Recently, a novel phylogenetic lineage of IDV with broad antigenicity has been discovered in California, USA, which has been named D/CA2019 (5).
The first-ever detection of IDV was documented in China in 2014, and three strains of bovine IDVs were identified in Shandong, northern China (9). In 2016-2017, a study investigating the prevalence of IDV in southern China found that IDVs were widely present in different farm animal species including dairy cattle, yellow cattle, buffalo, pigs and goats, and the positive rates varied from 4 to 40% (10). Besides southern and northern China, six IDV positive samples were recently detected from swine in Eastern China (11). The results from these studies indicated that IDVs may be prevalent throughout China. Phylogenetic analyses of IDVs in China suggested that they formed a distinct sub-clade (D/China sub-lineage), but all these isolates belong to D/OK lineage (11).
The continued circulation of IDV in Chinese cattle and swine herds, combined with its potential for causing respiratory illness outbreaks in cattle, prompts further investigation of epidemic characteristics and evolutionary dynamics of IDV in China. Here, we describe the first detection of D/Yama2019 lineage-like IDV, and genetically characterize this emerging IDV toward better understanding of the genetic diversity of IDVs in China.
Materials and Methods
Sampling, PCR Detection and Sequencing
During the winter and spring of 2021/2022, 250 nasal swab samples from apparently healthy cattle were collected in slaughterhouses around the city of Guangzhou, China. The cattle were sourced from various parts of the country including Gansu, Inner Mongolia, Anhui, Guangdong, Shaanxi, Liaoning, Ningxia, Hebei and Qinghai provinces. Viral RNA was extracted from above samples by using the Kit (Omega, E.Z.N.A.® Viral RNA Kit) and the isolated RNA was converted to cDNA by using the reverse transcription kit (Vazyme, HiScript® III 1st Strand cDNA Synthesis Kit), according to the manufacturer's instructions. For IDV detection, PCR assays with a pair of primers targeting P42 gene were conducted to generate an amplicon of 382bp. The forward primer is IDV-M-364-F: 5′- GATGTATGAAATGAGGGAGGAC-3′, and the reverse primer is IDV-M-746-R: 5′-AAGATTAGCCATTCCACTGAC-3′ (Table 1). The sample with the most intense PCR band (Figure 1A) was selected for further molecular characterization. Reverse transcription PCR (RT-PCR) of all seven genome segments of the isolate (D/bovine/CHN/JY3001/2021) was performed using conserved IDV-specific primers (Table 1). RT-PCR products were purified via an agarose gel DNA extraction kit (Omega, Gel Extraction Kit) and then Sanger sequenced with the IDV-specific primers. The sequences were assembled and proofread manually. The sequences were deposited in GenBank under accession numbers (ON415263-ON415269).
Table 1. Primers used in this study.
Figure 1. Detection and characterization of influenza D virus in nasal swab samples of cattle. (A) Agarose gel electrophoresis analysis of the amplicon products resulting from RT-PCR using the designated forward and reverse primers to amplify a fragment of 382 bp of IDV P42 gene. IDV-positive samples were labeled and highlighted with red boxes. “–”: negative controls; “+”: positive controls. (B) IDV-positive samples identified in this study. The sample ID and type reported together with animal information (sex and source). (C) Nucleotide sequence alignments of individual segments of the IDV isolate found in this study and representative strains of the D/OK lineage, D/660 lineage, D/Yama2016 lineage, D/Yama2019 lineage, D/CA2019 lineage and D/China sub-lineage. (D) Maximum-likelihood phylogenetic tree of D/bovine/CHN/JY3001/2021 and other known IDV strains. Nucleotide sequences of HEF genetic segment were aligned and analyzed using MEGA- X, with 1,000 bootstrap replicates. Bootstrap scores of at least 50 were shown to the left of the major nodes. Scale bar represents the number of substitutions per site. Each branch of five present lineages of IDV (D/OK-, D/660-, D/CA2019-, D/Yama2016-, and D/Yama2019-lineage) and the D/China sub-lineage was noted by blue, pink, green, purple, red and light blue color, respectively. The branch of D/bovine/CHN/JY3001/2021 was marked with a star. The D/bovine/CHN/JY3001/2021 and representative strains described in this study are bolded. Triangles indicate compressed branches with virus strains isolated in the same locations.
Sequence Alignment Analysis
Comparison of the HEF sequences of the D/bovine/CHN/JY3001/2021 isolate and representative strains of the D/OK lineage (D/swine/Oklahoma/1334/2011), D/660 lineage (D/bovine/Oklahoma/660/2013), D/Yama2016 lineage (D/bovine/Yamagata/10710/2016), D/Yama2019 lineage (D/bovine/Yamagata/1/2019), D/CA2019 (D/bovine/California/0363/2019) and D/China sub-lineage (D/bovine/Shandong/Y217/2014) were conducted using Blast tool in NCBI (http://ncbi.nlm.nih.gov).
Phylogenetic Analysis
IDV genome sequences were downloaded from the GenBank database (NCBI: http://ncbi.nlm.nih.gov) (accessed on April 24, 2022) for phylogenetic analysis. The sequences were aligned with the D/bovine/CHN/JY3001/2021 sequences by ClustalW in MEGA-X (12). Partial gene sequences (<80% of the full-length gene sequences) were not included in the analyses. Phylogenetic trees for individual genome segments were generated using the maximum likelihood method and Tamura-Nei model with 1,000 bootstrap replicates in MEGA-X (12). Branches with <50% bootstrap support were omitted.
Virus Isolation
The supernatants of the IDV-positive samples were inoculated to Madin-Darby canine kidney (MDCK) cells, Swine testicle (ST) cells and Human Rectal Tumor (HRT-18G) cells on 6-well plates. The plates were incubated at 33°C for 4–5 days. The supernatants were blindly passaged in ST cells. The infectivity was checked by determining hemagglutination (HA) titer or doing RT-PCR. Unfortunately, no IDV was obtained from the virus isolation attempts.
Results
RT-PCR Detection, Genome Sequencing and Alignment Analysis
The nasal swab samples (n = 250) from apparently healthy cattle were tested for IDV by RT-PCR. Ten clinical samples (4%) were found positive by amplification of partial sequences of P42 gene (Figure 1A). These positive samples were collected from cattle from farms located at different provinces of China (Figure 1B). Only one of the IDV-positive samples with the most-intense PCR band (Figure 1A) could be amplified for individual genomic segments. These PCR products of IDV genomic segments were purified and sequenced. Sequences were assembled and proofread, and then submitted to GenBank (D/bovine/CHN/JY3001/2021, accession number: ON415263-ON415269).
Alignment analyses showed that sequences of the IDV isolate found in this study shared the highest nucleotide similarity (>98%) to that of the representative strain (D/bovine/Yamagata/1/2019) of D/Yama2019 lineage but displayed the lowest nucleotide similarity (<95%) to that of the representative strain (D/bovine/California/0363/2019) of D/CA2019 lineage (Figure 1C). The sequence similarity between the D/bovine/CHN/JY3001/2021 and representative strains of other IDV lineages (D/OK lineage: D/swine/Oklahoma/1334/2011; D/660 lineage: D/bovine/Oklahoma/660/2013; D/Yama2016 lineage: D/bovine/Yamagata/10710/2016) ranged from 94 to 98% (Figure 1C). Moreover, the D/bovine/CHN/JY3001/2021 strain shared relatively low sequence homology (94–95%) with the representative strain of D/China sub-lineage (Figure 1C). These alignment analysis results highlighted a distinct IDV strain circulating in Chinese cattle herds.
Phylogenetic Analysis
To further elucidate the genetic relationship of this new isolate to previously characterized IDVs, all seven genomic segments were individually analyzed by conducting phylogenetic analysis using the maximum-likelihood method in MEGA-X (12), together with IDV genome sequences retrieved from GenBank database (accessed on April 24, 2022). The HEF segment of D/bovine/CHN/JY3001/2021 clustered within the D/Yama2019 lineage (Figure 1D). The remaining genomic segments (PB2, PB1, P3, NP, M and NS) of D/bovine/CHN/JY3001/2021 were also phylogenetically close to IDV strains of D/Yama2019 lineage (Figures 2A–F). These results suggested that IDV strain within D/Yama2019 lineage is circulating in cattle population in China.
Figure 2. Phylogenetic trees for non-HEF segments of IDVs. (A–F) In MEGA-X, maximum-likelihood analysis in combination with 1,000 bootstrap replicates was used to generate trees based on the nucleotide sequences of the PB2, PB1, P3, NP, M, and NS segments. Bootstrap scores of at least 50 were shown to the left of the major nodes. Scale bar represents the number of substitutions per site. The branch of D/bovine/CHN/JY3001/2021 was in red color and marked with a star. D/Yama2019 lineage-like and D/China sub lineage-like IDV strains were labeled with red and light blue colors, respectively. The D/bovine/CHN/JY3001/2021 was bolded. Triangles indicate compressed branches with virus strains isolated in the same locations.
Discussion
To date, all known IDVs identified in China cluster together and form a distinct sub-clade belonging to the D/OK lineage (9–11). However, genetic diversity of IDV remains unclear in Chinese herds. Here, we performed the genetic screening for IDV from apparently healthy cattle in slaughterhouses housing farm animals from all over the country. Of the tested nasal swab samples, IDV was detected in 10 out of 250 clinical samples (Figure 1A). Our results demonstrate that IDV continues to circulate in cattle herds in China.
The HEF segment is primarily used for phylogenetic analysis of IDV, which is the most variable gene segment of the IDV genome. Comparison of the HEF sequences of this emerging Chinese isolate and representatives of IDV lineages/sub-lineage (D/OK-, D/660-, D/Yama2016-, D/Yama2019-, D/CA2019-lineage and D/China sub-lineage) suggested that D/bovine/CHN/JY3001/2021 is a D/Yama2019 lineage-like IDV (Figures 1C,D). Corroborating this result, nucleotide sequences of all other segments of D/bovine/CHN/JY3001/2021 are close to that of D/Yama2019 lineage-like IDV strains in the phylogenetic analysis (Figure 2). To our knowledge, D/Yama2019 lineage-like IDVs have been isolated only in Japan (6, 13). This study reported the first detection of D/Yama2019 lineage-like IDV outside Japan. Therefore, further epidemiologic investigations are needed to determine the prevalence status of D/Yama2019 lineage-like IDVs.
There is no evidence of reassortments among the Chinese IDV isolates because the viruses displayed similar clustering in all phylogenetic trees according to individual genomic segments (Figures 1D, 2). Limited studies of IDVs in China combining with the small amount of Chinese IDV strains could also be the reason that we did not find evident reassortments concerning this new strain reported. Reassortment events have already been evidenced between two IDV lineages (5, 7, 14). Therefore, emergence of the D/Yama2019 lineage-like IDV in Chinese cattle herds highlights the importance of systematic virological and serological screenings to elucidate the diversity and evolution of IDVs in China.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/genbank/, ON415263-ON415269.
Ethics Statement
Ethical review and approval was not required for this investigation because the samples used in this study were taken from cattle by veterinarians.
Author Contributions
ML, S-LZ, and W-KW: conceptualization, supervision, and project administration. JY, TL, and ZWe: methodology and formal analysis. JY and ZWe: software, data curation, and visualization. JY, TL, S-LZ, and W-KW: validation. JY and TL: investigation. TL, SW, ZWa, JZ, MC, and FC: resources. JY: writing—original draft preparation. ML, S-LZ, W-KW, TL, and ZWe: writing—review and editing. ML, S-LZ, and JY: funding acquisition. All authors have read and agreed to the published version of the manuscript.
Funding
This study was supported by the following grants: the Collaborative Innovation Center Fund (Nos. XTXM202202 and XT202208) and the special funds (Nos. R2020PY-JC001, R2021YJ-QG008, and 202110TD) for scientific innovation strategy-construction of high level Academy of Agriculture Science and Talent Introduction Program from Guangdong Academy of Agricultural Sciences, the grant (No. 2021A1515010524) from the Guangdong Basic and Applied Basic Research Foundation, the grant (No. 2021TDQD002) from Maoming Branch Center of Guangdong Laboratory for Lingnan Modern Agricultural Science and Technology, and the grants (Nos. 2021KJ114 and 2021KJ119) from the Department of Agriculture and Rural Affairs of Guangdong Province.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We sincerely thank all of the members of the Department of Microbial Resources Research of Agro-biological Gene Research Center for their input into this work.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.939456/full#supplementary-material
Supplementary Figure S1. Sequence aligment of partial P42 gene amplicons of IDV-positive samples. The partial P42 gene sequences for IDV-positive samples ZS29-31, ZX56-58, ZS118-120, ZS130-132, ZS133-135, ZS16-18, and ZS63-65 were obtained by Sanger sequencing, and aligned together with the corresponding sequences of representative strains of the D/OK lineage, D/660 lineage, D/Yama2016 lineage, D/Yama2019 lineage, D/CA2019 lineage and D/China sub-lineage. The partial P42 gene sequences for IDV-positive samples ZS29-31, ZS130-132, ZS133-135, and ZS16-18 are identical to the corresponding sequence of D/swine/Oklahoma/1334/2011, the partial P42 gene sequences for IDV-positive samples ZX56-58 and ZS63-65 are close to the corresponding sequence of D/swine/Oklahoma/1334/2011, and the partial P42 gene sequence for the IDV-positive sample ZS118-120 is close to the corresponding sequence of D/bovine/Yamagata/1/2019.
References
1. Yu J, Li F, Wang D. The first decade of research advances in influenza D virus. J Gen Virol. (2021) 102:001529. doi: 10.1099/jgv.0.001529
2. Hause BM, Ducatez M, Collin EA, Ran Z, Liu R, et al. Isolation of a novel swine influenza virus from Oklahoma in 2011 which is distantly related to human influenza C viruses. PLoS Pathog. (2013) 9:e1003176. doi: 10.1371/journal.ppat.1003176
3. Hause BM, Collin EA, Liu R, Huang B, Sheng Z, Lu W, Wang D, Nelson EA, Li F. Characterization of a novel influenza virus in cattle and Swine: proposal for a new genus in the Orthomyxoviridae family. mBio. (2014) 5:e00031-14. doi: 10.1128/mBio.00031-14
4. Lion A, Secula A, Rancon C, Boulesteix O, Pinard A, Deslis A, et al. Enhanced pathogenesis caused by influenza D virus and Mycoplasma bovis coinfection in calves: a disease severity linked with overexpression of IFN-gamma as a key player of the enhanced innate immune response in lungs. Microbiol Spectr. (2021) 9:e0169021. doi: 10.1128/spectrum.01690-21
5. Huang C, Yu J, Hause BM, Park JY, Sreenivasan C, Uprety T, Sheng Z, Wang D, Li F. Emergence of new phylogenetic lineage of Influenza D virus with broad antigenicity in California, United States. Emerg Microbes Infect. (2021) 10:739–42. doi: 10.1080/22221751.2021.1910078
6. Murakami S, Sato R, Ishida H, Katayama M, Takenaka-Uema A, Horimoto T. Influenza D virus of new phylogenetic lineage, Japan. Emerg Infect Dis. (2020) 26:168–71. doi: 10.3201/eid2601.191092
7. Chiapponi C, Faccini S, Fusaro A, Moreno A, Prosperi A, et al. Detection of a new genetic cluster of influenza D virus in Italian Cattle. Viruses. (2019) 11:1110. doi: 10.3390/v11121110
8. Horimoto T, Hiono T, Mekata H, Odagiri T, Lei Z, Kobayashi T, et al. Nationwide distribution of bovine influenza D virus infection in Japan. PLoS ONE. (2016) 11:e0163828. doi: 10.1371/journal.pone.0163828
9. Jiang WM, Wang SC, Peng C, Yu JM, Zhuang QY, Hou GY, Liu S, Li JP, Chen JM. Identification of a potential novel type of influenza virus in Bovine in China. Virus Genes. (2014) 49:493–6. doi: 10.1007/s11262-014-1107-3
10. Zhai SL, Zhang H, Chen SN, Zhou X, Lin T, Liu R, et al. Influenza D virus in animal species in Guangdong Province, Southern China. Emerg Infect Dis. (2017) 23:1392–6. doi: 10.3201/eid2308.170059
11. He WT, Lu M, Xing G, Shao Y, Zhang M, Yang Y, et al. Emergence and adaptive evolution of influenza D virus. Microb Pathog. (2021) 160:105193. doi: 10.1016/j.micpath.2021.105193
12. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. (2018) 35:1547–9. doi: 10.1093/molbev/msy096
13. Hayakawa J, Masuko T, Takehana T, Suzuki T. Genetic and antigenic characterization and retrospective surveillance of bovine influenza D viruses identified in Hokkaido, Japan from 2018 to 2020. Viruses. (2020) 12:877. doi: 10.3390/v12080877
Keywords: influenza D virus, genetic lineage, phylogenetic analysis, cattle, China
Citation: Yu J, Li T, Wen Z, Wu S, Wang Z, Zheng J, Chen M, Chen F, Wei W-K, Zhai S-L and Liao M (2022) Identification of D/Yama2019 Lineage-Like Influenza D Virus in Chinese Cattle. Front. Vet. Sci. 9:939456. doi: 10.3389/fvets.2022.939456
Received: 09 May 2022; Accepted: 17 June 2022;
Published: 14 July 2022.
Edited by:
Min Yue, Zhejiang University, ChinaReviewed by:
Sushant Bhat, The Pirbright Institute, United KingdomBen Hause, South Dakota State University, United States
Zizhang Sheng, Columbia University Irving Medical Center, United States
Copyright © 2022 Yu, Li, Wen, Wu, Wang, Zheng, Chen, Chen, Wei, Zhai and Liao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wen-Kang Wei, weiwenkang@gdaas.cn; Shao-Lun Zhai, zhaishaolun@163.com; Ming Liao, mliao@scau.edu.cn
†These authors have contributed equally to this work