Next Article in Journal
PRRSV Non-Structural Proteins Orchestrate Porcine E3 Ubiquitin Ligase RNF122 to Promote PRRSV Proliferation
Next Article in Special Issue
Incidence, Disease Severity, and Follow-Up of Influenza A/A, A/B, and B/B Virus Dual Infections in Children: A Hospital-Based Digital Surveillance Program
Previous Article in Journal
Targeting Host PIM Protein Kinases Reduces Mayaro Virus Replication
Previous Article in Special Issue
Human SUMOylation Pathway Is Critical for Influenza B Virus
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 14
Issue 2
10.3390/v14020423
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Influenza D Virus in Danish Calves

by
Nicole B. Goecke
1,*,†,
Yuan Liang
1,†,
Nina D. Otten
1,
Charlotte K. Hjulsager
2 and
Lars E. Larsen
1
1
Department of Veterinary and Animal Sciences, University of Copenhagen, 1870 Frederiksberg, Denmark
2
Virus and Microbiological Special Diagnostics, Statens Serum Institut, 2300 Copenhagen, Denmark
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this article.
Viruses 2022, 14(2), 423; https://doi.org/10.3390/v14020423
Submission received: 26 January 2022 / Revised: 14 February 2022 / Accepted: 15 February 2022 / Published: 18 February 2022
(This article belongs to the Special Issue Non-A Influenza 2.0)
Browse Figures
Review Reports Versions Notes

Abstract

:
Influenza D virus (IDV) was first described in 2011 and has been found to mainly circulate among cattle and swine populations worldwide. Nasal swab samples were collected from 100 Danish calf herds (83 dairy and 17 veal herds) from 2018–2020. Influenza D virus was detected in 12 of the herds. Samples with the lowest cycle quantification value were selected for full genome sequencing. A hemagglutinin-esterase fusion (HEF) gene sequence from a Danish IDV collected in 2015 was also included in this study. Phylogenetic analysis showed that viruses from seven of the IDV-positive herds belonged to the D/OK lineage and clustered together in the HEF tree with the IDV collected in 2015. Viruses from the four other herds belonged to the D/660 lineage, where three of the viruses clustered closely together, while the fourth virus was more phylogenetically distant in all gene segments. The high level of genetic similarity between viruses from two different herds involved in calf trading suggests that transmission occurred through the movement of calves. This study is, to our knowledge, the first to describe the characterization of IDV in calves in Denmark.

1. Introduction

Influenza D virus (IDV) is a segmented negative-sense single-stranded RNA virus that belongs to the family Orthomyxoviridae, which also contain the species influenza A virus (IAV), B virus (IBV), and C virus (ICV). The genomes of ICV and IDV consist of seven gene segments, while the genomes of IAV and IBV consist of eight gene segments. The genome of IDV encodes nine proteins, including glycoprotein hemagglutinin-esterase fusion (HEF), which combines the functions of receptor recognition and binding, receptor destroying, and membrane fusion. IDV is genetically most closely related to ICV, but they share only around 50% amino acid identity for HEF and 70% for polymerase basic protein 1 (PB1), which is the most conserved segment. Furthermore, no cross reactivity between human ICV antisera and IDV has been seen [1,2,3].
IDV was isolated for the first time in 2011 from pigs exhibiting an influenza-like illness in Oklahoma in the USA [2]. However, serological evidence has subsequently indicated the circulation of IDV among cattle in Nebraska since 2003 [4]. Besides cattle and pigs, horses, sheep, goats, camelids, and feral swine have been found to be IDV-seropositive, which suggests a broad host-range of the virus [5,6,7,8]. Despite the broad host-range, it is believed that bovines are the natural reservoir of IDVs, due to the frequent isolation of these viruses from cattle and the finding of high seroprevalences in bovine herds [5,9]. Since the first discovery, IDVs have been detected in cattle in regions of the USA, Mexico, Canada, Asia (China, Japan), and Europe (France, Ireland, Italy, Luxembourg, the United Kingdom) [4,10,11,12,13,14,15,16,17,18]. The first European report of IDV in cattle was in 2012 in France [14].
Metagenomic studies have shown that IDVs can be associated with bovine respiratory disease (BRD), which is a multifactorial disease associated with different pathogens, stress conditions, environmental factors, and the health and immunological status of the animal [10,19]. Bacterial pathogens included in BRD are Mannheimia haemolytica, Pasteurella multocida, Mycoplasma bovis, and Histophilus somni, while viral agents known to participate in BRD are bovine adenovirus, bovine coronavirus, bovine herpesvirus 1, bovine parainfluenza virus type 3, bovine respiratory syncytial virus, and bovine viral diarrhea virus [20,21,22]. Only a few of these pathogens can successfully cause severe disease alone. The role of IDV in BRD in the field is still not clear, and further studies are needed to clarify this [23]. Experimental and field studies have shown that IDVs can cause mild to moderate respiratory signs in cattle, and that viruses are able to replicate in both the upper and lower bovine respiratory tract; however, they seem to have a preference for the upper tract [9,24,25,26].
The zoonotic potential of IDV is still unclear, and so far, there are no indications that IDV can cause disease in humans nor be transmitted among humans. However, studies have shown that IDV is able to transmit and replicate by direct contact in guinea pigs and ferrets, which are both established models for human influenza virus infection [2,27]. In addition, IDV-specific antibodies have been detected in humans with occupational exposure to cattle [28], and in the general population in Italy [29]. Another study demonstrated that IDV replicates efficiently in well-differentiated human airway epithelial cells, which act as an in vitro respiratory epithelium model for humans [30]. The mentioned studies support the hypothesis that IDV has zoonotic potential.
So far, IDVs have been phylogenetically classified into two major lineages, D/OK and D/660, in Europe and North America, while different lineages (D/Yama2016, D/Yama2019) have been found in Japan [31,32,33]. D/OK is thus far the most frequent reported lineage circulating in Europe, and until 2018, D/660 was unreported in this geographic area. However, an Italian IDV surveillance study performed in 2018–2019 found samples belonging to both the D/660 and D/OK lineages, showing that viruses from both lineages circulated in Italian cattle. The earliest European D/660 strain was detected in Italy in March 2018 from cattle imported from France. Since 2019, more D/660 strains than D/OK strains have been identified in Italy suggesting a shift in genotypes [33]. IDV strains with a reassortant genetic pattern, containing gene segments from both D/OK and D/660, have been reported in Italy and the USA [32,33]. Besides the circulation of the D/OK and D/660 lineages in Europe, a genetically divergent lineage containing two strains (D/bovine/France/2986/2012 and D/bovine/Ireland/007780/2014) has been described [14,16]. However, further sampling and sequencing is needed in order to be able to assess whether this represents the circulation of a new lineage in Europe. Furthermore, a different phylogenetic group of IDV has recently been identified in California and named D/CA2019 [34].
While IDVs have been shown to only cause mild or moderate respiratory disease in cattle, the continued circulation of IDV in European cattle and swine herds, combined with the potential for zoonotic transmission, prompts further elucidation of IDVs and their evolution. IDV status remains uninvestigated in many countries. Therefore, the aim of the present study was to genetically characterize IDVs in Danish calves to obtain more knowledge on the diversity of these viruses in Denmark.

2. Materials and Methods

2.1. Samples, Clinical Examination, and Herd Descriptions

Nasal swab samples were collected as part of a Danish research project where 100 Danish, intensive, commercial herds (83 dairy and 17 veal herds) were sampled. In Denmark, a veal herd is a rosé veal calf producing unit, where mainly bull calves are purchased from Danish dairy herds and raised for meat production [35]. Sampling was done in the period from September to April in 2018–2019 and 2019–2020, either once in each herd or at multiple time points. The samples were collected from three age groups in the dairy herds (0–10 days, 2–4 weeks, and 2–4 months), and two age groups in the veal herds (2 weeks after arrival (in general, 3–5 weeks of age) and at 3 months of age). The nasal swab samples were collected by inserting a sterile cotton swab approximately 8–10 cm into one nostril and turning the swab around for a few seconds, and immediately after, the swab was placed and stored in 1.5 mL phosphate-buffered saline (PBS). No prior cleaning of the nostril was performed. Further information about the handling and storage of the samples is described elsewhere [36].
A previous study investigated the occurrence of 11 respiratory and enteric pathogens in the 100 Danish calf herds at the first sampling time using a high-throughput real-time PCR platform, BioMark HD (Fluidigm, South San Francisco, CA, USA), and the 192.24 (samples.assays) dynamic array integrated fluidic circuit nanofluidic chip (Fluidigm) [36]. The herds that were sampled at multiple time points were also tested on the high-throughput real-time PCR (data not shown). Collectively, in the high-throughput real-time PCR analyses, samples (either individual samples and/or pools) from 1 dairy and 11 veal herds tested positive for IDV (Supplementary Materials, Table S1). In the present study, samples from some of the IDV-positive pools were subsequently tested individually by single-plex real-time reverse transcription-PCR (real-time RT-PCR) (Supplementary Materials, Table S1). The samples with the lowest cycle quantification (Ct) value, as determined by real-time RT-PCR, was selected for further whole-genome analysis, meaning that each of the 12 herds were represented by one sample, with the exception of Herd ID 5, where six samples from three different time points were included (Table 1). Furthermore, a bronchoalveolar lavage sample collected in 2015 from a calf in a veal herd was also included. Clinical observations or any other information were not available for this animal. Altogether, IDV sequences from 18 samples, representing 13 herds, were analyzed in this study. An overview of the IDV-positive pools and individual samples can be found in Supplementary Materials, Table S1.
At each sampling, clinical signs (coughing, rectal temperature, nasal and ocular discharge) were recorded for each sampled calf and scored based on a scoring system that was adapted from two existing scoring systems [37,38]. For each sign, clinical scores were rated depending on the severity, and a score of zero was given if the sign was not present. For coughing, the lowest score (zero) was given if no coughs were observed and the highest score (three) if repeated unprovoked coughs were observed. For nasal and ocular discharge, serous discharge equaled a score of one, while mucopurulent discharge equaled a score of two. Body temperature was measured rectally and categorized as either normal at temperatures below 39.3, or as febrile when exceeding this cut-off.

2.2. RNA Extraction

For RNA extraction, 200 µL nasal swab sample suspension was transferred to a tube containing 400 µL RLT-buffer (QIAGEN, Copenhagen, Denmark) containing 1% 2-Mercaptoethanol (Merck, Darmstadt, Germany). The RNA was subsequently extracted from the sample using an RNeasy mini kit (QIAGEN) on the extraction robot QIAcube Connect (QIAGEN) using the large sample protocol, version 2, and elution volume 60 µL, according to instructions from the supplier.

2.3. Real-Time RT-PCR

Real-time RT-PCR targeting the PB1 gene of IDV was used with the following primers: D/OK forward 5′-GCT GTT TGC AAG TTG ATG GG-3′, D/OK reverse 5′-TGA AAG CAG GTA ACT CCA AGG-3′, D/OK probe 5′-[FAM]-TTC AGG CAA GCA CCC GTA GGA TT-[BHQ1]-3′ [2], and with AgPath-IDTM One-Step RT-PCR Reagents kits (Applied BiosystemsTM, Foster City, CA, USA). The real-time RT-PCR was performed in a final volume of 25 µL with 5 µL extracted RNA template. A total of 12.5 µL RT-PCR buffer (2X) was mixed with 1.50 µL of each primer (10 µM), 0.75 µL probe (10 µM), 1 µL RT-PCR enzyme mix (25X), and 2.75 µL nuclease-free water. PCR reactions were run on Rotor-Gene Q (QIAGEN) using the thermal cycling conditions: 48 °C, 30 min (min); 94 °C, 10 min; 45 cycles x (94 °C, 15 s (s); 60 °C, 45 s). All reactions were run in duplicates and were considered positive if they had a Ct value < 36.

2.4. Sanger Sequencing

For the amplification of HEF PCR products of the IDV from 2015 for Sanger sequencing, previously described primers were used [14], along with the SuperScriptTM III One-Step RT-PCR System with PlatinumTM Taq High Fidelity DNA Polymerase kit (InvitrogenTM, Thermo Fischer Scientific, Roskilde, Denmark). The RT-PCR was performed in a final volume of 40 µL with 5 µL RNA. A total of 20 µL reaction mix (2X) was mixed with 0.20 µL of each primer (100 µM), 0.80 µL enzyme mix, and 13.7 µL nuclease-free water. The RT-PCR was run on a T3 Thermocycler (Biometra, Fredensborg, Denmark) with the following thermal cycling conditions: 50 °C, 30 min; 94 °C, 2 min; 5 cycles x (94 °C, 10 s; 60 °C, 30 s; 68 °C, 45 s); 35 cycles x (94 °C, 10 s; 57 °C, 30 s; 68 °C, 45 s); 5 cycles x (94 °C, 10 s; 54 °C, 30 s; 68 °C, 45 s); 68 °C, 5 min. The PCR products were visualized on 0.8% E-gelTM general purpose agarose gels (InvitrogenTM). The PCR products were subsequently purified with the High Pure PCR Product Purification Kit (Roche, Hørsholm, Denmark). The purified products were then Sanger sequenced with the aforementioned PCR primers at LGC Genomics (Berlin, Germany).

2.5. Whole-Genome Sequencing

For whole-genome sequencing (WGS), a modified version of a previously described one-tube PCR protocol was performed [39]. First, previously described primers were used [1] along with the SuperScriptTM III One-Step RT-PCR System with PlatinumTM Taq High Fidelity DNA Polymerase kit (InvitrogenTM). The RT-PCR was performed in a final volume of 40 µL with 4 µL being from the extracted RNA. A total of 20 µL reaction mix (2X) was mixed with 0.16 µL of each primer (50 µM), 0.80 µL enzyme mix, and 14.9 µL nuclease-free water. The PCR was run on a ProFlexTM PCR System (Applied Biosystems) under the following thermal cycling conditions: 42 °C, 60 min; 94 °C, 2 min; 5 cycles x (94 °C, 30 s; 45 °C, 30 s; 68 °C, 3 min); 31 cycles x (94 °C, 30 s; 57 °C, 20 s; 68 °C, 3 min); 68 °C, 7 min. The amplified mixtures were visualized on 0.8% E-gelTM (InvitrogenTM) and purified with the High Pure PCR Product Purification Kit (Roche). Libraries were built with the Nextera XT DNA Library Preparation Kit and sequenced on the Illumina MiSeq system (Illumina, Copenhagen, Denmark) with MiSeq Reagent Kit v2, 500 cycles, according to instructions from the supplier.

2.6. Consensus Sequence Generation

Data from the Sanger sequencing was imported into the CLC Main Workbench (version 20.0.3) (QIAGEN), where the sequences were assembled against a reference sequence of the HEF gene (D/bovine/Italy/46484/2015). The WGS fastq data were imported as paired-end to the CLC Genomics Workbench (version 21.0.3) (QIAGEN) where reads were trimmed using default settings. The reads were subsequently mapped against reference sequences for each gene segment (D/swine/Italy/199724-3/2015). Consensus sequences were then extracted for each gene segment. Proofreading of the sequences was performed manually. The sequences of this study are available with the accession numbers OM468185–OM468297 (Supplementary Materials, Table S2).

2.7. Maximum Likelihood Phylogeny

Complete IDV sequences available were selected from the GenBank database (NCBI, http://ncbi.nlm.nih.gov (accessed on 3 November 2020)) for phylogenetic analysis, as well as sequences based on high percentage identity by nucleotide BLAST searches. The sequences were aligned with the Danish IDV sequences using the CLC Main Workbench (Version 20.0.3). Model testing was performed in the CLC Main Workbench (Version 20.0.3) (QIAGEN) for each gene segment (PB2, PB1, P3, HEF, NP, P42, NS) alignment to determine the best fitting substitution model. All gene segments, except for NS, were best described by GTR + G + T. NS was best described by a HKY + G + T substitution model. Maximum likelihood phylogeny analysis with 1000 bootstrap replicates was performed in the same program for each gene segment, with the appropriate substitution models selected. Rate variation was included and the analysis was set to estimate gamma distribution parameters. The remaining settings were kept at default. All trees were rooted on gene segments from C/Bovine/Montana/12/2016. The root was subsequently removed. The resulting trees were visualized in FigTree (version 1.4.4) [40].

2.8. Prediction of Glycosylation Sites

To predict the N-linked glycosylation sites of HEF, the bioinformatics tool NetNGly–1.0 (http://www.cbs.dtu.dk/services/NetNGlyc/) (accessed on 10 November 2021) was used.

3. Results

3.1. Clinical Signs

The clinical registrations and observations for the IDV-positive calves included in the analysis are listed in Table 1. Additionally, some IDV-positive calves were also co-infected with other bovine respiratory pathogens (Supplementary Materials, Table S1). Three of the calves coughed repeatedly by provocation (score of two) or without provocation (score of three), while no coughing was observed for the rest of the animals (score of zero). Nasal discharge was observed for all of the calves as either serous (score of one) or mucopurulent (score of two) discharge. Ocular discharge was not observed for eight calves (score of zero), while either serous (score of one) or mucopurulent (score of two) discharge was observed for the rest of the animals. Rectal temperature was measured to be between 38.3–40.5 ℃. Clinical observations and information about sex and age were not available for the calf sampled in 2015 (D/bovine/Denmark/47398-3622/2015-11-03).

3.2. WGS and Phylogenetic Analysis

WGS was obtained from samples from 11 of the selected herds (Table 1, Herd ID 1–11), while sequencing was unsuccessful for Herd ID 12, likely due to low viral load in the sample. For the IDV-positive sample from 2015 (Table 1, Herd ID 13) only the HEF gene was sequenced using Sanger sequencing.
The gene segments of the Danish IDVs had a nucleotide sequence identity of 94.84–100% and amino acid sequence identity of 93.88–100%. All seven gene segments were individually analyzed by inferring phylogenetic trees using the maximum-likelihood method. All gene segments of Danish IDVs collected in 2019, as well as one virus from 2020 (D/bovine/Denmark/3313408142/2020), clustered together (Figure 1 and Figure 2) and can be phylogenetically classified as belonging to the D/OK lineage (Figure 1). The Danish viruses in this lineage had a nucleotide sequence identity of 98.19–100% and amino acid sequence identity of 98.04–100%. The HEF segment sequenced from the sample from 2015 also clustered within the D/OK lineage. Six of the eleven analyzed Danish viruses from 2019–2020 were from the same herd, but collected at three different time points (Table 1, Herd ID 5). From this herd, all gene segments of the four viruses sampled on 15 January 2019 were almost identical (99.80–100% nucleotide identity and 99.74–100% amino acid identity), whereas they were slightly more different from viruses sampled in the same herd on 20 February 2019 and 24 April 2019 (99.25–100% nucleotide identity and 99.09–100% amino acid identity) (Figure 1 and Figure 2). Notably, the latter two viruses clustered with D/bovine/Denmark/6171103755-1/2019 collected in a dairy herd on 6 February (Herd ID 1), that had sold calves to Herd ID 5, which were collected and transported to the facility by staff members working at Herd ID 5. Geographically, these two herds are also closely located (approximately 22 km apart). One of the viruses collected in 2019 (D/bovine/Denmark/5246407608-2/2019) clustered in another branch and shared the highest degree of similarity in the HEF gene (99.53% nucleotide identity) to the virus collected in 2015 (D/bovine/Denmark/47398-3622/2015). D/bovine/Denmark/5246407608-2/2019 was collected in a herd that was located approximately 50 km away. The remaining D/OK viruses were collected in herds that were not closely located nor were connected by trade. Overall, the Danish D/OK IDVs clustered closest to IDVs collected from bovine and swine in Italy, France, and Northern Ireland, while they were more phylogenetically distant from IDVs collected in the USA and China (Figure 1). There was no evidence of reassortments among the Danish IDV strains since the viruses showed similar clustering in all phylogenetic trees based on the different gene segments (Figure 1 and Figure 2).
All gene segments of the Danish IDVs collected in 2020 clustered together except for D/bovine/Denmark/3313408142-6/2020 and belonged to the D/660 lineage. The nucleotide sequence identity for the Danish viruses in this lineage was 96.32–100% and the amino acid sequence identity was 98.35–100%. D/bovine/Denmark/3727003200-7/2020 was the most phylogenetically distant for all gene segments in comparison with the three remaining IDVs from 2020 (Figure 1 and Figure 2). The Danish D/660 IDVs clustered closest together with IDVs collected from bovines in Italy, while they were more phylogenetically distant from IDVs collected in North America (Figure 1). Geographically, the three herds, where D/bovine/Denmark/5256205576-8/2020 (Herd ID 8), D/bovine/Denmark/5995505714-9/2020 (Herd ID 9), and D/bovine/Denmark/1052304114-10/2020 (Herd ID 10) were collected, are located within a vicinity of 89 km, whereas D/bovine/Denmark/3727003200-7/2020 was collected from a herd (Herd ID 7) situated between 216-268 km from the other herds. Similarly, there was no known trade connections between the Danish herds with D/660 IDV detections.

3.3. Prediction of Glycosylation Sites

The prediction of the N-linked glycosylation sites of HEF showed the presence of several glycosylation sites for the D/OK and D/660 lineages. The bioinformatics tool NetNGly predicted the positions 28, 54, 146, 249, 346, 513, and 613 to be glycosylated with different specificities, and the glycosylation sites were present in both lineages (Supplementary Materials, Tables S3 and S4). For the Danish D/OK sequences, three of these differed by having one less glycosylation site. In the D/bovine/Denmark/5195604790-5/2019-02-20, Thr 56 was replaced by Asn, which transforms the glycosylation site (position 54–57) from a highly specific site to a non-glycosylation site. In the D/bovine/Denmark/5246407608-2/2019-03-01 and D/bovine/Denmark/47398-3622/2015-11-03, Thr 515 was replaced by Ile which changes the glycosylation site (position 513–516) to a non-glycosylation site. The Danish D/660 sequences showed a similar N-glycosylation site prediction pattern.

4. Discussion

IDV has mainly been detected in ruminants, with cattle being considered the primary reservoir. The prevalence of IDV has not been elucidated in animal populations in Denmark. Here, we performed the characterization of IDV detected in a research project where 100 cattle herds across Denmark were tested for the presence of bovine pathogens, including IDV using a high-throughput real-time PCR system [35]. Thus, the samples did not originate from diagnostic submissions. Of the investigated herds, IDV was detected in 12 out of 100 herds. Some of the infected calves displayed minor clinical signs such as coughing, discharge from the eyes and/or nose, and slightly elevated rectal temperatures. Whether these clinical signs were caused by the IDV infection alone is not known, as the animals could also be positive for other pathogens such as bovine coronavirus, M. bovis, Mycoplasma, M. haemolytica, H. somni, P. multocida, or T. pyogenes. Nevertheless, these signs correspond well with an upper respiratory viral infection and are similar to findings elsewhere [9].
Genome characterization of the Danish IDVs showed that the viruses collected in 2015 and 2019 belonged to the D/OK lineage, while most of the viruses collected in 2020 belonged to the D/660 linage. Until 2018, the D/660 linage was unreported in Europe, and the first European detection was in Italy in March 2018 from cattle imported from France. Subsequently, more D/660 strains than D/OK strains are being found in Italy, indicating an ongoing shift in genotype [33]. Despite the limited number of samples analyzed in the present study, a similar shift in genotype also seem to be taking place in Denmark.
Our phylogenetic analysis revealed that the Danish IDVs detected in the herds from 2019 and one from 2020 are the most closely related to the first Danish IDV detection which was from 2015, suggesting that all D/OK-like viruses in Denmark could originate from a single introduction, but due to a gap in surveillance between 2015 and 2019 this remains speculative. Similarly, based on the similarity of the D/660-like viruses, it is likely that they share a common ancestor, indicating a single incursion followed by genetic drift or multiple incursions of almost identical viruses. Herd ID 5 was sampled at three different time points in 2019, with the viruses found in January clustering together and viruses from February and April being more closely related. Since only a part of the calves from the herd were sampled, we are unable to determine whether the distant viruses from February and April are due to new incursions or whether they had already been circulating in January. Since this is the first screening for IDVs in Denmark, the route of incursion of IDV into Danish calves and how long IDVs have been circulating in Denmark is also unknown. One route of incursion of IDV could be the import of cattle from countries infected with IDV; however, the annual import of living cattle is very limited, with less than 200 imported cattle per year from other European countries (unpublished data). Most of the cattle imported since 2017 originate from Sweden, Germany, and Finland, with some being from France. Further IDV surveillance in other European countries is needed elucidate the circulation and route of introduction of IDVs across borders. Within Denmark, trade with calves is common, where bull calves are often sold from dairy herds to veal herds at an early age. In the present study, the owner of one of the veal herds (Herd ID 5) bought and directly transported calves from a dairy herd (Herd ID 1). Notably, the viruses from these two herds were closely related which suggests inter-herd transmission by infected calves. To our knowledge, there is no connection among the remaining IDV-positive calf herds.
Amino acids in the HEF gene under positive selection have been identified, and it is believed that these might affect the receptor binding cavity relevant for broader cell tropism. One study identified two amino sites to be under positive selection (residues 251 and 289), both of which are located in the globular domain of HEF [41]. The amino acids Ala and Thr were found to be present at residue 251, while in the Danish IDVs, Thr was present in both the D/OK and D/660 sequences. At residue 289, Val was present in the Danish D/OK sequences, while Ser was present in the Danish D/660 sequences. Others have found Ala, Ser, or Thr at residue 289, and showed that with the presence of Thr, a new site for N-linked glycosylation was created [41]. Different N-linked glycosylation sites were found in the Danish IDVs, where three of the D/OK sequences differed from the others by lacking the glycosylation site at residue 513, which has also been reported from France [14]. Since experimental studies have yet to be conducted on Danish IDVs, it is unclear if the different glycosylation profiles affect the transmission and/or pathogenicity of these viruses.
This study is the first to describe the detection and full-genome characterization of IDVs in Denmark and in Danish calves. Since this virus has been found in other animal species such as swine, feral swine, camelids, ruminants, and horses in other countries, we are planning future systematic virological and serological screenings to elucidate the presence and diversity of IDVs in Denmark in these species.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v14020423/s1: Table S1: Information about collection date, number of tested animals and samples, high-throughput real-time PCR results, and IDV single-plex real-time RT-PCR results for the IDV positive herds. Table S2: Accession numbers for the sequences of this study. Table S3: Predicted glycosylation sites for D/660 IDVs. Table S4: Predicted glycosylation sites for D/OK IDVs.

Author Contributions

Conceptualization, N.B.G., Y.L. and L.E.L.; methodology, Y.L.; software, Y.L.; formal analysis, N.B.G. and Y.L.; investigation, N.B.G., Y.L. and N.D.O.; data curation, N.B.G., Y.L. and C.K.H.; writing—original draft preparation, N.B.G. and Y.L.; writing—review and editing, N.B.G., Y.L., N.D.O., C.K.H. and L.E.L.; visualization, N.B.G. and Y.L.; supervision, N.D.O., C.K.H. and L.E.L.; project administration, N.D.O. and L.E.L.; funding acquisition, L.E.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research is a part of the project “Robuste Kalve – Godt begyndt er halvt fuldendt” which was funded by the Danish Cattle Levy Fund and the Danish Milk Levy Fund.

Institutional Review Board Statement

Ethical review and approval was not required for this investigation, since the samples used in the present study were taken from field animals by veterinarians.

Informed Consent Statement

Not applicable.

Data Availability Statement

The sequences are available at GenBank (www.ncbi.nlm.nih.gov/genbank) with the accession numbers OM468185–OM468297.

Acknowledgments

The authors would like to thank Anne Marie Michelsen, Bodil H. Nielsen, Stine Lindgren, Masja Feline Reipurth Søndergaard, Henrik Hjul Møller, Franziska Helene Skaarup Pedersen, Thomas Damm Poulsen, Jesper Kjærgård Davidsen, Sofie Jeppesen, Maëva Durand, Lucie Jeanne Marie Dupont, Anaëlle Bouqueau, Alicia Fokje Klompmaker, Maria Brydensholt, Sascha Coes, Jensine Wilm, Mette Gerdes Wilson, Emma Madsen, Helge Kromann, and Henrik Læssøe Martin for collection of the samples. Thank you to the laboratory technicians, Nina Dam Grønnegaard and Hue Thi Thanh Tran, for help with the sample analysis.

Conflicts 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.

References

  1. Hause, B.M.; Collin, E.A.; Liu, R.; Huang, B.; Sheng, Z.; Lu, W.; Wang, D.; Nelson, E.A.; Li, F. Characterization of a novel influenza virus in cattle and swine: Proposal for a new genus in the Orthomyxoviridae family. mBio 2014, 5, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Hause, B.M.; Ducatez, M.; Collin, E.A.; Ran, Z.; Liu, R.; Sheng, Z.; Armien, A.; Kaplan, B.; Chakravarty, S.; Hoppe, A.D.; 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, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Asha, K.; Kumar, B. Emerging Influenza D Virus Threat: What We Know so Far! J. Clin. Med. 2019, 8, 192. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Luo, J.; Ferguson, L.; Smith, D.R.; Woolums, A.R.; Epperson, W.B.; Wan, X.F. Serological evidence for high prevalence of Influenza D Viruses in Cattle, Nebraska, United States, 2003–2004. Virology 2017, 501, 88–91. [Google Scholar] [CrossRef]
  5. Oliva, J.; Eichenbaum, A.; Belin, J.; Gaudino, M.; Guillotin, J.; Alzieu, J.P.; Nicollet, P.; Brugidou, R.; Gueneau, E.; Michel, E.; et al. Serological Evidence of Influenza D Virus Circulation Among Cattle and Small Ruminants in France. Viruses 2019, 11, 516. [Google Scholar] [CrossRef] [Green Version]
  6. Ferguson, L.; Luo, K.; Olivier, A.K.; Cunningham, F.L.; Blackmon, S.; Hanson-Dorr, K.; Sun, H.; Baroch, J.; Lutman, M.W.; Quade, B.; et al. Influenza D Virus Infection in Feral Swine Populations, United States. Emerg. Infect. Dis. 2018, 24, 1020–1028. [Google Scholar] [CrossRef] [Green Version]
  7. Nedland, H.; Wollman, J.; Sreenivasan, C.; Quast, M.; Singrey, A.; Fawcett, L.; Christopher-Hennings, J.; Nelson, E.; Kaushik, R.S.; Wang, D.; et al. Serological evidence for the co-circulation of two lineages of influenza D viruses in equine populations of the Midwest United States. Zoonoses Public Heal. 2018, 65, e148–e154. [Google Scholar] [CrossRef]
  8. Salem, E.; Cook, E.A.J.; Lbacha, H.A.; Oliva, J.; Awoume, F.; Aplogan, G.L.; Hymann, E.C.; Muloi, D.; Deem, S.L.; Alali, S.; et al. Serologic Evidence for Influenza C and D Virus among Ruminants and Camelids, Africa, 1991–2015. Emerg. Infect. Dis. 2017, 23, 1556–1559. [Google Scholar] [CrossRef] [Green Version]
  9. Ferguson, L.; Olivier, A.K.; Genova, S.; Epperson, W.B.; Smith, D.R.; Schneider, L.; Barton, K.; McCuan, K.; Webby, R.J.; Wan, X.-F. Pathogenesis of Influenza D Virus in Cattle. J. Virol. 2016, 90, 5636–5642. [Google Scholar] [CrossRef] [Green Version]
  10. Mitra, N.; Cernicchiaro, N.; Torres, S.; Li, F.; Hause, B.M. Metagenomic characterization of the virome associated with bovine respiratory disease in feedlot cattle identified novel viruses and suggests an etiologic role for influenza D virus. J. Gen. Virol. 2016, 97, 1771–1784. [Google Scholar] [CrossRef]
  11. Ferguson, L.; Eckard, L.; Epperson, W.B.; Long, L.P.; Smith, D.; Huston, C.; Genova, S.; Webby, R.; Wan, X.F. Influenza D virus infection in Mississippi beef cattle. Virology 2015, 486, 28–34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Jiang, W.M.; Wang, S.C.; Peng, C.; Yu, J.M.; Zhuang, Q.Y.; Hou, G.Y.; Liu, S.; Li, J.P.; Chen, J.M. Identification of a potential novel type of influenza virus in Bovine in China. Virus Genes 2014, 49, 493–496. [Google Scholar] [CrossRef] [PubMed]
  13. Murakami, S.; Endoh, M.; Kobayashi, T.; Takenaka-Uema, A.; Chambers, J.K.; Uchida, K.; Nishihara, M.; Hause, B.; Horimoto, T. Influenza D Virus Infection in Herd of Cattle, Japan. Emerg. Infect. Dis. 2016, 22, 1517–1519. [Google Scholar] [CrossRef]
  14. Ducatez, M.F.; Pelletier, C.; Meyer, G. Influenza d virus in cattle, France, 2011–2014. Emerg. Infect. Dis. 2015, 21, 368–371. [Google Scholar] [CrossRef]
  15. Chiapponi, C.; Faccini, S.; De Mattia, A.; Baioni, L.; Barbieri, I.; Rosignoli, C.; Nigrelli, A.; Foni, E. Detection of influenza D virus among swine and cattle, Italy. Emerg. Infect. Dis. 2016, 22, 352–354. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Flynn, O.; Gallagher, C.; Mooney, J.; Irvine, C.; Ducatez, M.; Hause, B.; McGrath, G.; Ryan, E. Influenza D virus in cattle, Ireland. Emerg. Infect. Dis. 2018, 24, 389–391. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Dane, H.; Duffy, C.; Guelbenzu, M.; Hause, B.; Fee, S.; Forster, F.; McMenamy, M.J.; Lemon, K. Detection of influenza D virus in bovine respiratory disease samples, UK. Transbound. Emerg. Dis. 2019, 66, 2184–2187. [Google Scholar] [CrossRef]
  18. Snoeck, C.J.; Oliva, J.; Pauly, M.; Losch, S.; Wildschutz, F.; Muller, C.P.; Hübschen, J.M.; Ducatez, M.F. Influenza D virus circulation in cattle and swine, Luxembourg, 2012–2016. Emerg. Infect. Dis. 2018, 24, 1388–1389. [Google Scholar] [CrossRef] [Green Version]
  19. Ng, T.F.F.; Kondov, N.O.; Deng, X.; Van Eenennaam, A.; Neibergs, H.L.; Delwart, E. A Metagenomics and Case-Control Study to Identify Viruses Associated with Bovine Respiratory Disease. J. Virol. 2015, 89, 5340–5349. [Google Scholar] [CrossRef] [Green Version]
  20. Ellis, J.A. Update on viral pathogenesis in BRD. Anim. Health Res. Rev. 2009, 10, 149–153. [Google Scholar] [CrossRef]
  21. Fulton, R.W. Bovine respiratory disease research (1983–2009). Anim. Health Res. Rev. 2009, 10, 131–139. [Google Scholar] [CrossRef] [PubMed]
  22. Griffin, D.; Chengappa, M.M.; Kuszak, J.; McVey, D.S. Bacterial pathogens of the bovine respiratory disease complex. Vet. Clin. Food Anim. 2010, 26, 381–394. [Google Scholar] [CrossRef] [PubMed]
  23. Gaudino, M.; Moreno, A.; Snoeck, C.J.; Zohari, S.; Saegerman, C.; O’Donovan, T.; Ryan, E.; Zanni, I.; Foni, E.; Sausy, A.; et al. Emerging Influenza D virus infection in European livestock as determined in serology studies: Are we underestimating its spread over the continent? Transbound. Emerg. Dis. 2020, 68, 1125–1135. [Google Scholar] [CrossRef] [PubMed]
  24. Salem, E.; Hägglund, S.; Cassard, H.; Corre, T.; Näslund, K.; Foret, C.; Gauthier, D.; Pinard, A.; Delverdier, M.; Zohari, S.; et al. Pathogenesis, Host Innate Immune Response, and Aerosol Transmission of Influenza D Virus in Cattle. J. Virol. 2019, 93, e01853-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Hause, B.M.; Huntimer, L.; Falkenberg, S.; Henningson, J.; Lechtenberg, K.; Halbur, T. An inactivated influenza D virus vaccine partially protects cattle from respiratory disease caused by homologous challenge. Vet. Microbiol. 2017, 199, 47–53. [Google Scholar] [CrossRef] [PubMed]
  26. Mazzetto, E.; Bortolami, A.; Fusaro, A.; Mazzacan, E.; Maniero, S.; Vascellari, M.; Beato, M.S.; Schiavon, E.; Chiapponi, C.; Terregino, C.; et al. Replication of Influenza D Viruses of Bovine and Swine Origin in Ovine Respiratory Explants and Their Attachment to the Respiratory Tract of Bovine, Sheep, Goat, Horse, and Swine. Front. Microbiol. 2020, 11, 1–14. [Google Scholar] [CrossRef] [PubMed]
  27. Sreenivasan, C.; Thomas, M.; Sheng, Z.; Hause, B.M.; Collin, E.A.; Knudsen, D.E.B.; Pillatzki, A.; Nelson, E.; Wang, D.; Kaushik, R.S.; et al. Replication and Transmission of the Novel Bovine Influenza D Virus in a Guinea Pig Model. J. Virol. 2015, 89, 11990–12001. [Google Scholar] [CrossRef] [Green Version]
  28. White, S.K.; Ma, W.; McDaniel, C.J.; Gray, G.C.; Lednicky, J.A. Serologic evidence of exposure to influenza D virus among persons with occupational contact with cattle. J. Clin. Virol. 2016, 81, 31–33. [Google Scholar] [CrossRef]
  29. Trombetta, C.M.; Marchi, S.; Manini, I.; Kistner, O.; Li, F.; Piu, P.; Manenti, A.; Biuso, F.; Sreenivasan, C.; Druce, J.; et al. Influenza D virus: Serological Evidence in the Italian Population from 2005 to 2017. Viruses 2020, 12, 30. [Google Scholar] [CrossRef] [Green Version]
  30. Holwerda, M.; Kelly, J.; Laloli, L.; Stürmer, I.; Portmann, J.; Stalder, H.; Dijkman, R. Determining the replication kinetics and cellular tropism of influenza D virus on primary well-differentiated human airway epithelial cells. Viruses 2019, 11, 377. [Google Scholar] [CrossRef] [Green Version]
  31. 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–171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Collin, E.A.; Sheng, Z.; Lang, Y.; Ma, W.; Hause, B.M.; Li, F. Cocirculation of Two Distinct Genetic and Antigenic Lineages of Proposed Influenza D Virus in Cattle. J. Virol. 2015, 89, 1036–1042. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Chiapponi, C.; Faccini, S.; Fusaro, A.; Moreno, A.; Prosperi, A.; Merenda, M.; Baioni, L.; Gabbi, V.; Rosignoli, C.; Alborali, G.L.; et al. Detection of a new genetic cluster of influenza D virus in Italian cattle. Viruses 2019, 11, 1110. [Google Scholar] [CrossRef] [Green Version]
  34. Huang, C.; Yu, J.; Hause, B.M.; Park, J.Y.; 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–742. [Google Scholar] [CrossRef] [PubMed]
  35. Fertner, M.; Toft, N.; Læssøe, H.; Boklund, A. A register-based study of the antimicrobial usage in Danish veal calves and young bulls. Prev. Vet. Med. 2016, 131, 41–47. [Google Scholar] [CrossRef]
  36. Goecke, N.B.; Nielsen, B.H.; Petersen, M.B.; Larsen, L.E. Design of a High-Throughput Real-Time PCR System for Detection of Bovine Respiratory and Enteric Pathogens. Front. Vet. Sci. 2021, 8, 1–15. [Google Scholar] [CrossRef]
  37. Love, W.J.; Lehenbauer, T.W.; Kass, P.H.; Van Eenennaam, A.L.; Aly, S.S. Development of a novel clinical scoring systemfor on-farmdiagnosis of bovine respiratory disease in pre-weaned dairy calves. PeerJ 2014, 2, 1–25. [Google Scholar] [CrossRef] [Green Version]
  38. McGuirk, S.M. Disease Management of Dairy Calves and Heifers. Vet. Clin. Food Anim. 2008, 24, 139–153. [Google Scholar] [CrossRef]
  39. Ryt-Hansen, P.; Larsen, I.; Kristensen, C.S.; Krog, J.S.; Larsen, L.E. Limited impact of influenza A virus vaccination of piglets in an enzootic infected sow herd. Res. Vet. Sci. 2019, 127, 47–56. [Google Scholar] [CrossRef]
  40. Rambaut, A. FigTree version 1.4.4. 2018. Available online: http://tree.bio.ed.ac.uk/software/figtree/.
  41. He, W.T.; Lu, M.; Xing, G.; Shao, Y.; Zhang, M.; Yang, Y.; Li, X.; Zhang, L.; Li, G.; Cao, Z.; et al. Emergence and adaptive evolution of influenza D virus. Microb. Pathog. 2021, 160, 1–9. [Google Scholar] [CrossRef]
Viruses 14 00423 g001 550
Figure 1. The maximum-likelihood phylogenetic tree of the IDV HEF gene. The highlighted areas denote which lineage the viruses belong to. The sequence of the virus detected in Denmark in 2015 is shown in green text, depicted in red are the IDV viruses found in 2019, and those from 2020 are in blue. All bootstrap values < 70% have been removed. The scale bar indicates nucleotide substitutions per site per year.
Figure 1. The maximum-likelihood phylogenetic tree of the IDV HEF gene. The highlighted areas denote which lineage the viruses belong to. The sequence of the virus detected in Denmark in 2015 is shown in green text, depicted in red are the IDV viruses found in 2019, and those from 2020 are in blue. All bootstrap values < 70% have been removed. The scale bar indicates nucleotide substitutions per site per year.
Viruses 14 00423 g001
Viruses 14 00423 g002 550
Figure 2. The maximum-likelihood phylogenetic trees of the IDV PB2, PB1, P3, NP, P42, and NS genes. The sequences of IDVs detected in Denmark in 2019 are depicted in red, and depicted in blue are sequences from viruses found in 2020. All bootstrap values < 70% have been removed. The scale bar indicates nucleotide substitutions per site per year.
Figure 2. The maximum-likelihood phylogenetic trees of the IDV PB2, PB1, P3, NP, P42, and NS genes. The sequences of IDVs detected in Denmark in 2019 are depicted in red, and depicted in blue are sequences from viruses found in 2020. All bootstrap values < 70% have been removed. The scale bar indicates nucleotide substitutions per site per year.
Viruses 14 00423 g002
Table
Table 1. IDV-positive samples included in the study. The herd of each strain is identified by a number. The date of sampling, source, and production type are reported together with animal information (sex, age) and clinical registrations (temperature, cough, nasal and eye discharge). The lineage attribution is based on the HEF gene.
Table 1. IDV-positive samples included in the study. The herd of each strain is identified by a number. The date of sampling, source, and production type are reported together with animal information (sex, age) and clinical registrations (temperature, cough, nasal and eye discharge). The lineage attribution is based on the HEF gene.
StrainHerd IDDate of SamplingProduction TypeSourceSexAgeTemperature
(℃)		Cough *Nasal
Discharge ■		Eye
Discharge ■		Lineage
HEF
D/bovine/Denmark/6171103755-1/2019106-02-2019DairyNSHeifer3 m38.5021D/OK
D/bovine/Denmark/5246407608-2/2019201-03-2019VealNSBull2 waa39.9020D/OK
D/bovine/Denmark/5358904609-3/2019312-03-2019VealNSBull3 m39.0020D/OK
D/bovine/Denmark/1173303525-4/2019405-04-2019VealNSBull3 m38.9020D/OK
D/bovine/Denmark/4559602762-5/2019515-01-2019VealNSBull2 waa38.7020D/OK
D/bovine/Denmark/5022903340-5/2019515-01-2019VealNSBull2 waa38.9011D/OK
D/bovine/Denmark/5161205135-5/2019515-01-2019VealNSBull2 waa38.5011D/OK
D/bovine/Denmark/5195604714-5/2019515-01-2019VealNSBull3 m39.3020D/OK
D/bovine/Denmark/5195604790-5/2019520-02-2019VealNSBull2 waa39.1021D/OK
D/bovine/Denmark/5224403563-5/2019524-04-2019VealNSBull3 m38.3010D/OK
D/bovine/Denmark/3313408142-6/2020608-01-2020VealNSBull3 m39.0322D/OK
D/bovine/Denmark/3727003200-7/2020712-03-2020VealNSBull2 waa40.5022D/660
D/bovine/Denmark/5256205576-8/2020806-03-2020VealNSHeifer3 m39.6322D/660
D/bovine/Denmark/5995505714-9/2020911-02-2020VealNSBull3 m39.5022D/660
D/bovine/Denmark/1052304114-10/20201016-01-2020VealNSBull3 m39.2220D/660
D/bovine/Denmark/3058506689-11/20191109-10-2019VealNSBull3 m39.7020D/OK
D/bovine/Denmark/4070604622-12/20191216-01-2019VealNSBull3 m39.5011-
D/bovine/Denmark/47398-3622/20151303-11-2015VealBALndndndndndndD/OK
NS: nasal swab, BAL: bronchoalveolar lavage, m: month, waa: weeks after arrival, -: sequencing was not possible due to low viral load. * 0: no coughing observed; 1: a single cough when provoked; 2: repeated coughs when provoked or a single spontaneous cough; 3: repeated unprovoked coughs. ■ 0: normal, no signs of discharge; 1: clear, serous discharge, and no sign of mucopurulent or purulent exudate; 2: mucopurulent or purulent discharge.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Goecke, N.B.; Liang, Y.; Otten, N.D.; Hjulsager, C.K.; Larsen, L.E. Characterization of Influenza D Virus in Danish Calves. Viruses 2022, 14, 423. https://doi.org/10.3390/v14020423

AMA Style

Goecke NB, Liang Y, Otten ND, Hjulsager CK, Larsen LE. Characterization of Influenza D Virus in Danish Calves. Viruses. 2022; 14(2):423. https://doi.org/10.3390/v14020423

Chicago/Turabian Style

Goecke, Nicole B., Yuan Liang, Nina D. Otten, Charlotte K. Hjulsager, and Lars E. Larsen. 2022. "Characterization of Influenza D Virus in Danish Calves" Viruses 14, no. 2: 423. https://doi.org/10.3390/v14020423

APA Style

Goecke, N. B., Liang, Y., Otten, N. D., Hjulsager, C. K., & Larsen, L. E. (2022). Characterization of Influenza D Virus in Danish Calves. Viruses, 14(2), 423. https://doi.org/10.3390/v14020423

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop