# Viral tropism and detection of clade 2.3.4.4b H5N8 highly pathogenic avian influenza viruses in feathers of ducks and geese

Nicolas Gaide, Charlotte Foret-Lucas, Thomas Figueroa, Timothée Vergne, Marie-Noëlle Lucas, Luc Robertet, Marie Souvestre, Guillaume Croville, Guillaume Le Loc'h, Maxence Delverdier, et al.

► To cite this version:

Nicolas Gaide, Charlotte Foret-Lucas, Thomas Figueroa, Timothée Vergne, Marie-Noëlle Lucas, et al. Viral tropism and detection of clade 2.3.4.4b H5N8 highly pathogenic avian influenza viruses in feathers of ducks and geese. Scientific Reports, 2021, 11 (1), pp.5928. 10.1038/s41598-021-85109-5. hal-03194244

HAL Id: hal-03194244

https://hal.inrae.fr/hal-03194244

Submitted on 9 Apr 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Viral tropism and detection of clade 2.3.4.4b H5N8 highly pathogenic avian influenza viruses in feathers of ducks and geese

Nicolas Gaide<fnref n="1" />,², Charlotte Foret-Lucas<fnref n="1" />,², Thomas Figueroa<fnref n="1" />, Timothée Vergne<fnref n="1" />, Marie-Noëlle Lucas<fnref n="1" />, Luc Robertet<fnref n="1" />, Marie Souvestre<fnref n="1" />, Guillaume Croville<fnref n="1" />, Guillaume Le Loc'h<fnref n="1" />, Maxence Delverdier<fnref n="1" /> & Jean-Luc Guérin<fnref n="1" />✉

Highly Pathogenic Avian Influenza viruses (HPAIVs) display a tissue pantropism, which implies a possible spread in feathers. HPAIV detection from feathers had been evaluated for H5N1 or H7N1 HPAIVs. It was suggested that viral RNA loads could be equivalent or higher in samples of immature feather compared to tracheal (TS) or cloacal swabs (CS). We investigated the suitability of feathers for the detection of clade 2.3.4.4b H5N8 HPAIV in ducks and geese field samples. In the six H5N8 positive flocks that were included in this study, TS, CS and immature wing feathers were taken from at least 10 birds. Molecular loads were then estimated using real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) targeting H5 and M genes. In all flocks, viral loads were at least equivalent between feather and swab samples and in most cases up to 10³ higher in feathers. Bayesian modelling confirmed that, in infected poultry, RT-qPCR was much more likely to be positive when applied on a feather sample only (estimated sensitivity between 0.89 and 0.96 depending on the positivity threshold) than on a combination of a tracheal and a cloacal swab (estimated sensitivity between 0.45 and 0.68 depending on the positivity threshold). Viral tropism and lesions in feathers were evaluated by histopathology and immunohistochemistry. Epithelial necrosis of immature feathers and follicles was observed concurrently with positive viral antigen detection and leukocytic infiltration of pulp. Accurate detection of clade 2.3.4.4b HPAIVs in feather samples were finally confirmed with experimental H5N8 infection on 10-week-old mule ducks, as viral loads at 3, 5 and 7 days post-infection were higher in feathers than in tracheal or cloacal swabs. However, feather samples were associated with lower viral loads than tracheal swabs at day 1, suggesting better detectability of the virus in feathers in the later course of infection. These results, based on both field cases and experimental infections, suggest that feather samples should be included in the toolbox of samples for detection of clade 2.3.4.4b HPAI viruses, at least in ducks and geese.

Avian influenza is a highly contagious infectious disease due to an Influenza A virus belonging to *Orthomyxoviridae* family<fnref n="1" />. Wild birds are considered as the natural reservoir that contributes to the spread of the disease, generating occasionally panzootic with massive death of wild and domestic birds and significant economic losses²,³.

During winter 2016–2017, a highly pathogenic avian influenza virus (HPAIV) H5N8, belonging to A/goose/Guangdong/1/1996 clade 2.3.4.4b lineage, caused epidemic outbreaks in Europe⁴. France was the most affected country with more than 6 million culled birds and 484 outbreaks. Free-range foie gras ducks were particularly affected, representing more than 80% of outbreaks⁵. Epidemiological studies pointed out that poor biosecurity practices, interactions with wild birds and transportation of ducks could have contributed to the spread of the virus within this livestock sector⁵⁻⁹.

Classical recommendation for the surveillance of influenza involves the collection of tracheal or oropharyngeal and cloacal swabs on live and dead birds. Samples can be processed individually or pooled to detect Influenza A specific nucleic acid signature, by real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR), followed by pathotyping for cleavage site identification¹⁰,¹¹.

<footnote n="1">¹IHAP, Université de Toulouse, ENVT, INRAE, 23 Chemin des Capelles, 31076 Toulouse Cedex 3, France. ²These authors contributed equally: Nicolas Gaide and Charlotte Foret-Lucas. ✉email: jean-luc.guerin@envt.fr</footnote>
Highly pathogenic avian influenza viruses (HPAIVs) differ from low pathogenic AIVs (LPAIVs) by the mutational acquisition of multibasic haemagglutinin cleavage site, that shifts mucosal epitheliotropism to tissue pantropism, leading to systemic diseases with high mortality in gallinaceous species³,¹². Since HPAIVs display tissue pantropism, it implies that they could spread in feathers. Accurate detection of HPAIVs from feather samples has been demonstrated for HP Asian H5N1 and H7N1 viruses. These few studies suggest that viral loads detected in immature feathers are equivalent or even higher than those detected on either tracheal or cloacal swabs¹³⁻¹⁸. These results were also confirmed on carcasses of dead birds and detached feathers, suggesting better preservation of viral particles in feathers than in viscera¹⁵,¹⁹. However, almost all these studies were carried out in the context of experimental infections. Combining different approaches, we investigated the suitability of feather samples from both field and experimental infections to detect clade 2.3.4.4b H5N8 HPAIV in ducks.

## Methods

**Field cases sampling.** Six flocks were included in the study based on sudden onset of nervous signs and mortality. These flocks have been tested positive for M and H5 AI genes by RT-qPCR¹⁰. Finally, all these flocks were confirmed positive for 2.3.4.4 H5N8 HPAIV between January and March 2017 by the French National Laboratory for Avian Influenza and Newcastle disease, using the official procedures¹⁰. They included 4 mule duck flocks (5 to 13 weeks of age), 1 Pekin duck flock (8 weeks of age) and 1 flock of geese (8 weeks of age). In each flock, at least 10 non-clinically affected birds were randomly taken and sampled using tracheal (TS) and cloacal swabs (CS), as well as immature wing feathers (remiges). Mature wing feathers were also collected on geese flock. Skin samples with feathers were harvested following necropsy of clinically-affected ducks (n = 11) and geese (n = 4) and fixed in 10% buffered formalin to investigate viral tropism and feather lesions. On-farm investigations and collection of samples were performed in strict compliance with regulation and biosecurity procedures, with the authorization and supervision of official veterinary services and before implementation of sanitary culling and destruction of carcasses.

**Animal experiment.** For the experiment, twenty 10-week-old mule ducks, clinically healthy and AIV negative, were included in the study, in a BSL3 laboratory with negative pressure isolators, absolute air filtration and more than 2 m² surface per isolator. Animals were randomly split into four groups to represent the different modalities of transmission. In the first isolator, five ducks were inoculated at day 0 with 10⁶ EID₅₀ H5N8 (A/duck/France/171201g/2017(H5N8)) through choanal instillation (Infected group) with additional 5 non-infected ducks (contact group). Five non-infected ducks were put in an isolator receiving air from the infected group isolator (aerosol group). Five additional non-infected ducks were housed in another isolator into a fully closed system (control group). During the seven following days, ducks were fed and had access to water ad libitum. For each animal, tracheal swabs (TS), cloacal swabs (CS), feather pulp (FP) and blood were sampled at days 0 to check negativity (all samples) and then from day 1 to day 7.

**Viral isolation and amplification.** To obtain a viral stock for experimental infections, viral isolation and amplification were conducted on feather pulp of an infected mule duck flock 2 with the highest Cq value on RT-qPCR. For viral isolation, feather pulp was placed in 500 µL of 1X PBS and vortexed for 30 s. Then 150 µL of supernatant at the dilution 1:100 in PBS with penicillin and streptomycin 4X were inoculated to 10-day-old specific-pathogen-free embryonated chicken eggs (INRAE PFIE, Nouzilly, France). Eggs were kept at 37 °C in a humidity chamber during 48 h and then placed at 4 °C overnight. Allantoic liquids were harvested and titrated by hemagglutination assay with fresh and washed chicken red blood cells.

For viral amplification, a similar passage on embryonated eggs was done on allantoic liquid with the highest hemagglutination assay titer. Embryonated chicken eggs were inoculated with 100 µL of a 1:100 dilution in PBS with penicillin and streptomycin 2X. The viral stock was then aliquoted and stored at – 80 °C. In parallel, the titer of this viral stock was also determined by Tissue Culture Infectious Dose (TCID50) method on MDCK cells and calculated by the Spearman & Kärber algorithm²⁰. All these steps were processed in a BSL3 lab in strict compliance with biosafety procedures. Sequences of amplified A/duck/France/171201g/2017(H5N8) virus, first isolated from a mule duck feather, were deposited on GenBank database under MK208604-MK208610 accession numbers.

**Histopathology and immunohistochemistry.** Formalin-fixed cutaneous tissues were routinely processed in paraffin blocks, sectioned at 4 µm and stained with hematoxylin and eosin (HE) for microscopic examination. A diagram representing histological structures of avian feathered skin is provided in supplementary data (Suppl. Fig. S1). The number of follicular sections examined per skin section was 8 in average for ducks (n = 11 birds) and 4 for geese (n = 4 birds). Immunostainings were performed on paraffin-embedded skin sections with a monoclonal mouse antibody against Influenza A virus nucleoprotein (NP) (Argene, Sherley, USA, 11–030). Immunohistochemical reaction included an antigen retrieval with pronase 0,05% 10 min at 37 °C, a peroxidase blocking step 5 min at room temperature (Dako, Glostrup, Denmark, S2023) and saturation of non-specific binding sites with normal goat serum (Dako, Glostrup, Denmark, X0907) 25 min at room temperature, before overnight incubation at 4 °C with anti-NP antibody at the dilution 1:50. The immunohistochemical staining was revealed with a biotinylated polyclonal goat anti-mouse immunoglobulins conjugated with horseradish peroxidase (HRP) antibody (Dako, Glostrup, Denmark, LSAB2 system-HRP, K0675) and the aminoethylcarbazole chromogen of the HRP (Dako, Glostrup, Denmark, AEC + Substrate-Chromogen, K3469). Negative controls included sections incubated either without specific primary antibody or with another monoclonal antibody of the same isotype (IgG2). Histopathological analyses were carried out by three veterinary pathologists certified by the European College of Veterinary Pathologists (ECVP).
<table><caption>**Table 1.** Proportion of naturally H5N8-infected ducks and geese with histopathological lesions and positive viral detection in immature and mature feathers. NE non evaluated. <sup>a</sup>Focal to extensive lytic necrosis and apoptosis/single cell necrosis in epidermal structures. Mixed leukocytic lympho-plasmacytic rich infiltration in dermal structures. <sup>b</sup>Total number of positive subjects / Total subjects (percentage of positive subjects). <sup>c</sup>Anti-nucleoprotein A immunohistochemistry. <sup>d</sup>Positive antigen detection in follicular epidermis.</caption><thead><tr><th rowspan="2"></th><th colspan="4">Immature feather</th><th rowspan="2">Mature feather</th></tr><tr><th>Feather dermis</th><th>Feather epidermis</th><th>Follicular epidermis</th><th>Follicular dermis</th></tr></thead><tbody><tr><td colspan="6"><strong>Histopathology<sup>a</sup></strong></td></tr><tr><td>Duck n=11</td><td>10/11 (91%)<sup>b</sup></td><td>7/11 (64%)</td><td>8/11 (73%)</td><td>8/11 (73%)</td><td>0/11 (0%)</td></tr><tr><td>Goose n=4</td><td>2/4 (50%)</td><td>2/4 (50%)</td><td>1/4 (25%)</td><td>0/4 (0%)</td><td>0/4 (0%)</td></tr><tr><td colspan="6"><strong>Immunohistochemistry<sup>c</sup></strong></td></tr><tr><td>Duck n=11</td><td>7/11 (64%)</td><td>11/11 (100%)</td><td>11/11 (100%)</td><td>8/11 (73%)</td><td>10/11 (91%)<sup>d</sup></td></tr><tr><td>Goose n=4</td><td>0/4 (0%)</td><td>4/4 (100%)</td><td>3/4 (75%)</td><td>0/4 (0%)</td><td>NE</td></tr></tbody></table>

**RNA extraction and RT-qPCR.** All samples were processed in BSL3 laboratory until lysis in strict compliance with biosafety procedures. Swabs were placed in 500 µL of 1X PBS and vortexed for 30 s. Feather pulp was stored in 500 µL of 1X PBS and then vortexed for 30 s before RNA extraction. For experimental infected group, additional 30–60 mg feather pulp were concurrently incubated in 200 µL PBS 1X containing 0,8 µg/µL of proteinase K (Thermo fisher #EO0492) for 20 min at 37 °C to compare RNA extraction yields with and without proteinase K. For all sample types, RNA was extracted from 140 µL of supernatant with QIAamp Viral RNA Mini Kit (Qiagen #52906) and then kept at – 80 °C.

RT-qPCR was performed on 2 µL RNA using iTaq Universal SYBR Green Supermix (Bio-Rad #172–5125). M and H5 genes were targeted with M52C / M253R primers for M gene<sup>21</sup>, and H5 HP EA F2 (5'-TCCTTG CAACAGGACTAAG-3') / H5 HP EA R (5'-GTCTACCATTCYTGCCA-3')<sup>22</sup>. Absolute quantification was performed with a plasmid range from 10<sup>2</sup> to 10<sup>7</sup> copies/reaction (2 µL). RT-qPCR reaction and results interpretation were performed on a LightCycler 96 instrument (Roche).

**Droplet digital PCR.** For cDNA synthesis, 5 µL of RNA extracted from swabs and feather pulp (from field cases flock 2) were processed according to the instructions of the RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific #15255146). As negative controls, 5 µL of Milli-Q water was processed in parallel. A 2 µL aliquot of the cDNA from the reverse transcription step of samples and controls were used for PCR amplification with the QX200 ddPCR EvaGreen Supermix (Bio-Rad 1864033) and M and H5 primers previously described. ddPCR were processed according to the instruction for QX200 AutoDG Droplet Digital PCR System (Bio-Rad) established by the Transcriptomic platform from Genotoul (Get-Santé INSERM I2MC Toulouse). Analyses were performed with the dedicated software Quantsoft (Bio-Rad).

**Bayesian data analysis.** To estimate the sensitivity and the specificity of detection in the three different sample types (cloacal swab, tracheal swab and feather pulp) at the individual level, we used a latent class modelling approach embedded in a Bayesian framework. This type of analysis has been extensively used to model cross-detection of individuals whose true epidemiological status (infected or not) is assessed using imperfect diagnostic tests of unknown sensitivity and specificity<sup>23</sup>. To do so, we defined five thresholds of test positivity (5, 10, 20, 50 and 100 copies of RNA/µL), and derived the observed frequency of the 2<sup>3</sup>=8 different combinations of individual test results for each threshold. The combination of test results was assumed to be distributed according to a multinomial distribution of parameters n=61 birds and eight probabilities, expressed as a combination of seven parameters to be estimated, i.e. the proportion of infected birds in the sample and the sensitivity and specificity of detection in the three different sample types. The sensitivity and specificity parameters were assumed to be threshold-specific (leading to 30 parameters, 6 for each of the five thresholds of positivity) while the proportion of infected birds was considered the same irrespective of the threshold. The value of the 31 parameters was estimated in a Bayesian framework. We allocated uniform prior distributions for the proportion of infected birds and all sensitivity parameters. All specificity parameters were assigned an informative beta prior distribution to reflect that we were 97.5% and 50% confident that the true value of the specificity parameters were greater than 0.8 and 0.98, respectively.

**Approval for animal experiments.** This study was carried out in compliance with European animal welfare regulation. The protocols were approved by the Animal Care and Use Committee “Comité d’éthique en Science et Santé Animales—115”, protocol number 13205-2018012311319709.

## Results

**Integument lesions and virus detection of field cases in H5N8-infected ducks and geese.** Duck skin sections consisted of a mixture of immature and mature feathers, whereas geese presented a majority of immature feathers. Lesions were observed in immature feathers in 91% of ducks and 50% of affected geese (Table 1, Fig. 1A). Mature feathers were within normal limits for all subjects. In ducks, lesions were consistent with focal to extensive necrosis of feather (64%) and follicular (73%) epidermis in association with leu-
<figure><img src="image_3.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Lesion and viral detection in immature and mature feathers of naturally H5N8-infected ducks and geese. Hematoxylin and Eosin (HE) and anti-nucleoprotein A Immunohistochemistry (IHC). (A) Lesion and viral detection in immature and mature feathers, duck infected with H5N8 clade 2.3.4.4b, HE. (A1) Immature feather with extensive acute necrosis of feather and follicular epidermis (arrow). The feather pulp is infiltrated by leukocytes, mainly lymphocytes and plasma cells (asterisk). (A2) Mature feather with calamus (thin arrow) within normal limits. (A3) Abundant viral antigen is present within necrotic debris and outer layers of feather and follicular epidermis. (A4) Viral antigen is present in association with follicular epidermis of mature feather, under the calamus (arrowhead). (B) Growing and differentiating feather epidermis, duck infected with H5N8 clade 2.3.4.4b. (B1) At early stage of differentiation, feather epidermis is stratified. Necrosis is visible at the outer (corneal) layer (arrow) with viral antigen colocalization. (B2) At the stage of early differentiation, necrosis and single-cell necrosis/apoptosis can be seen within marginal plates and barb ridges (arrowhead) and with viral antigen colocalization. (B3) At the stage of cornication, keratinized barbs and barbules are disrupted and intermixed with necrotic cells and debris (asterisk) with viral antigen colocalization.</figcaption></figure>

leukocytic infiltration of feather pulp (91%), including predominantly lymphocytes and plasma cells, and fewer macrophages and heterophils. Necrosis was predominant in the outer (corneal) layer and could be observed from early to late stage of feather epithelium differentiation, including early and late barb ridges and marginal plates (Fig. 1B). Variable amount of leukocytic infiltration including lymphocytes, plasma cells and heterophils
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Viral RNA loads in feathers (F), tracheal swabs (TS) and cloacal swabs (CS) in 6 naturally infected flocks with H5N8 HPAIV. Mule ducks (4 flocks, n = 10 birds per flock), pekin ducks (1 flock, n = 11) and geese (1 flock, n = 10). Results are expressed as log₁₀ viral RNA copies/µL.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Estimated sensitivity (left panel) and specificity (right panel) of the three different sample types (feather sample, tracheal swab and oropharyngeal swab) for 5 different thresholds of positivity (5, 10, 20, 50 and 100 copies of RNA/µL).</figcaption></figure>

were observed in the follicular dermis in 73% of ducks. In geese, lesions were less frequent: half of the subjects presented leukocytic infiltration of pulp and necrosis of feather epidermis (50%) and one subject presented epidermal focal necrosis in follicular epidermis (25%).

In immature feathers, viral antigen was detected in all ducks and geese samples, mainly in feather and follicular epidermis. Viral antigen immunostaining was patchy to diffuse in distribution and severely intense in the outer (corneal) layer epidermis and necrotic areas (Fig. 1A). Interestingly, viral antigen was commonly detected at all levels of feather epithelial differentiation in ducks and geese, including barb ridges, marginal plate, barbs and
<table><caption>**Table 2.** Comparison of the detectability of the infection for each threshold of positivity between a sampling strategy involving one feather only and a sampling strategy involving a tracheal swab and an oropharyngeal swab (in which case an animal tested positive is defined by a positive result for at least one of the two swabs). The brackets represent the 95% credible interval of the posterior distributions.</caption><thead><tr><th>Parameter</th><th>Threshold of positivity (nb copies of RNA/µL)</th><th>1 feather sample</th><th>1 tracheal and 1 oropharyngeal swabs</th></tr></thead><tbody><tr><td rowspan="5">Sensitivity</td><td>5</td><td>0.96 (0.87-1.00)</td><td>0.68 (0.56-0.80)</td></tr><tr><td>10</td><td>0.96 (0.86-1.00)</td><td>0.63 (0.51-0.75)</td></tr><tr><td>20</td><td>0.95 (0.83-1.00)</td><td>0.60 (0.47-0.72)</td></tr><tr><td>50</td><td>0.93 (0.79-0.99)</td><td>0.57 (0.44-0.70)</td></tr><tr><td>100</td><td>0.89 (0.73-0.99)</td><td>0.45 (0.31-0.59)</td></tr><tr><td rowspan="5">Specificity</td><td>5</td><td>0.95 (0.68-1.00)</td><td>0.96 (0.82-1.00)</td></tr><tr><td>10</td><td>0.96 (0.71-1.00)</td><td>0.96 (0.81-1.00)</td></tr><tr><td>20</td><td>0.97 (0.78-1.00)</td><td>0.96 (0.83-1.00)</td></tr><tr><td>50</td><td>0.98 (0.82-1.00)</td><td>0.97 (0.85-1.00)</td></tr><tr><td>100</td><td>0.98 (0.83-1.00)</td><td>0.97 (0.85-1.00)</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Viral RNA loads in experimentally H5N8-infected 10-week-old mule ducks. Virus loads in feathers (F), tracheal swabs (TS) or cloacal swabs (CS) were determined from day 1 to day 7 post-infection for infected (pink), contact (yellow) and aerosol (blue) groups. Absolute quantification was performed by RT-qPCR targeting the viral M gene from extracted RNAs in parallel with an M gene plasmid range. Each dot represents an individual value and horizontal bars correspond to means. Limit of detection for each experiment.</figcaption></figure>

barbules (Fig. 1B). In mature feathers, 91% of ducks presented immunostaining of follicular epidermis (corneal layer) without associated lesions (Fig. 1A). In ducks, viral antigens were occasionally detected in mesenchymal cell in feather pulp (64%) and follicular dermis (73%).

In all six duck flocks included in the study, RNA loads of HPAIV in feather pulp were at least equivalent and in most cases up to 10³ higher than those detected in either TS or CS (Fig. 2). Incubation of swabs with proteinase
K was assessed on a few samples and did not result in a significant increase in detectability (Suppl. Fig. S2). For a given positivity threshold, feather pulp tested positive much more frequently than TS or CS (Suppl. Table T1). A latent class modelling approach was used to estimate the sensitivity and specificity of RNA detection in each sample type (feather pulp, CS and TS). Irrespective of the threshold of positivity (5, 10, 20, 50 and 100 copies of RNA/µL), the sensitivity of RNA detection by RT-qPCR in feather pulp of infected birds was much higher than in TS or CS, while the specificity did not substantially drop (Fig. 3). As an example, for a positivity threshold of 10 copies of RNA/mL, the sensitivity of the detection in feather pulp was estimated at 0.96 (95% credible interval (95% CI) 0.86–1.00), while the sensitivity of detection in tracheal swabs and cloacal swabs were estimated at 0.46 (95% CI 0.32–0.60) and 0.31 (95% CI 0.19–0.45), respectively.

As shown in Table 2, testing one sample of feather pulp was even shown to be more sensitive (Se = 0.96 (95% CI 0.86–1.00)) than the classical sampling recommendation, involving the combination of one tracheal swab and one cloacal swab (Se = 0.63 (95% CI 0.51–0.75) for a positivity threshold of 10 RNA copies/mL). To further confirm the quantitation in tissues, a Droplet Digital PCR assay, using the same primer set, was implemented on a subset of samples. This assay showed a substantial correlation with RT-qPCR data and confirmed significant viral loads in feathers (Suppl. Fig. S3).

**Dynamics of HPAIV detection in feathers and swabs in experimentally infected ducks.** Experimental infection resulted in totally asymptomatic infections in all birds, despite viral RNA detection in swabs and feathers starting from day 1 post-infection (Fig. 4) and RT-qPCR Ct values as low as 16 in some birds at peak of infection (data not shown). Control birds remained virus negative throughout the study. Three direct contact ducks became positive in tracheal swabs as early as day 1 post-infection and by day 3 post-infection, they were all positive in all sample types. Regarding birds connected through aerosols, some became positive at day 3 post-infection, with two positive feather samples, five positive tracheal swabs and three positive cloacal swabs, and they were all positive in all sample types by day 5 post-infection. Viral loads were always substantially higher in tracheal than cloacal swabs. For feather samples, viral loads were equivalent or substantially higher than tracheal and cloacal swabs at day 3, 5 and 7 post-infection, but not at day 1 post-infection when compared with tracheal swabs of infected and contact groups (Fig. 4). The peak of viral excretion occurred between day 3 and 5 post-infection and the relative viral loads in feathers increased, becoming obviously higher than tracheal and even more, cloacal swabs from day 5 post-infection until the end of the course of infection (Fig. 4).

Discussion

The present study investigated the suitability of HPAIV H5N8 clade 2.3.4.4b viral detection in feathers with histopathological correlates in ducks and geese. Previous studies focused on the pathology and detectability of H5N1 and H7N1 HPAIVs in the integument after experimental infections of ducks, chicken and geese¹³⁻¹⁹. To the best of our knowledge, the present data provide new insights about HPAIV detection in feather for H5N8 clade 2.3.4.4b through both naturally occurring field cases and experimental challenges.

Our data support previous findings that HPAIVs can replicate in feathers and induce feather folliculitis. Ducks presented epithelial necrotizing folliculitis in immature (growing) feathers, with concurrent viral antigenic detection. This lesional pattern supports a strong tropism for feather and follicular epidermis, corneal outer layers in particular, confirming previous descriptions²⁴⁻²⁶. In 91% of ducks, a marked leukocytic, lymphoplasmacytic rich, infiltration was observed in the feather pulp, which contrasts with previous experimental infections of ducks with Indonesian strain of H5N1, where no inflammatory lesions were found in association with epithelial necrosis²⁴. Viral detection in leukocytes was uncommon, more frequent in mesenchymal stromal cells of dermal structures as previously reported²⁵. In geese, no major lesion pattern dominated and half of the subjects had feathers within normal limits.

Interestingly, cell death with concomitant viral antigens could be detected at different stages of feather epithelial differentiation, from the early stage up to cornifying barbs, barbules of ducks and geese. Because viral replication seems to occur into feather structures intended to be external, these data suggest that growth and differentiation of feathers epidermis could be a source of viral production for pantropic H5N8 Clade 2.3.4.4b and could potentially carry and disseminate viral particles. The role of feathers in environmental dissemination of avian influenza viruses has been suggested, including in wild waterfowl¹⁹,²⁷. Viral infectivity and RNA detectability in feathers detached from bodies has been experimentally investigated in ducks and revealed that HPAIV (H5N1) can longer persist compared to drinking water and feces, for 160 days at 4 °C and 15 days at 20 °C¹⁹. Additionally, a study conducted on wild waterbirds has suggested that preen oil increases AIVs survivability on feathers and can act as protection for viral particles²⁷. Altogether, these data indicate that feathers could be a source of environmental contamination. It would be interesting to compare feather infectivity potential between species, since we can presumably predict higher infectivity for species showing an epitheliotropic viral expression in feather compared to chicken where viral expression is endotheliotropic and restricted to the dermal pulp, confined in the calamus. Further studies are needed to confirm this significance of feathers in the epidemiology of avian influenza.

Viral loads and tropism may vary greatly during the course of infection in a same bird, which could impact greatly the efficiency of different samples¹⁵. In our field investigations, flocks were sampled at day 1 to 7 after the onset of clinical signs and/or mortality, as notified by the farmer or the attending veterinarian. In all batches, viral loads in feathers were higher than in swabs, regardless of this estimated stage of infection. In controlled conditions, as during our experimental infection, viral loads in feathers increased during infection, suggesting that feathers may not be the most relevant samples at the very early stage of the infection (day 1 post-infection).

Classical recommendation for the surveillance of influenza involves the collection of tracheal or oropharyngeal and cloacal swabs on live and dead birds¹⁰,¹¹. According to our results, the sensitivity of combining one
tracheal and one cloacal swab to test for the presence of viral RNA in an infected duck would be 0.63 (0.51–0.75) for a positivity threshold of 10 copies/mL. It is striking that only one feather pulp sample outperform the classical sample recommendation, as its sensitivity was estimated at 0.96 (0.86–1.00). Moreover, feathers are easy to collect, transport and treat for molecular analyses. The much higher sensitivity of detection in feather pulp suggests that feather samples could be diluted without decreasing the overall sensitivity at the flock level, and therefore that they could be pooled further, which would reduce the costs of analysis at the flock scale and improve the overall economic efficiency of surveillance. Combined with swabs in a sampling strategy, feathers would also provide supplementary information about invasiveness of the virus. In a previous experimental study conducted on european quails, H5N1 and H7N1 HPAIVs inoculation resulted in viral RNA detection in feathers from day 1 up to day 6 and day 12, respectively. In contrast, H7N2 LPAIV detection remained negative during all the experiment from 1 to 15 days post-infection²⁸. Assuming that AIV spread in feather is the result of viremia and systemic infection, feather folliculitis in combination with positive viral detection could presumably reflect the tissue pantropism of a virus strain and suggest infection by a HPAIV. Further investigations are needed to investigate feather infection as a signature of HPAIVs, including variations between strains.

Clinical expression of HPAI in ducks is poor, as evidenced on the field and in experimental infections²⁹. Our experiment confirmed that ducks may excrete huge loads of virus for days, contaminate very efficiently the rest of the flock without expressing clinical signs. This poor clinical expression in ducks stresses, even more, the pivotal importance of an optimized surveillance strategy, based on intensive testing, at least before farm-to-farm transportation of ducks.

These data, based on a selection of both field and experimental cases of H5N8-infected ducks and geese, suggest that feather samples should be included in the toolbox of samples for detection of clade 2.3.4.4b HPAI viruses, at least in ducks. Further investigations are needed in experimental settings and on a wide range of HPAIVs to ascertain the relevance of their insertion in the diagnostic protocols of HPAI surveillance.

## Data availability

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Received: 5 September 2020; Accepted: 22 December 2020

Published online: 15 March 2021

## References

1. International Committee on Taxonomy of Viruses (ICTV). *Orthomyxoviridae*. https://talk.ictvonline.org/taxonomy/.

2. Causey, D. & Edwards, S. V. Ecology of avian influenza virus in birds. *J. Infect. Dis.* **197**, 29–33 (2008).

3. Olsen, B. *et al*. Global patterns of influenza A virus in wild birds. *Science* **312**, 384–388 (2006).

4. Napp, S., Majó, N., Sánchez-Gonzalez, R. & Vergara-Alert, J. Emergence and spread of highly pathogenic avian influenza A(H5N8) in Europe in 2016–2017. *Transbound. Emerg. Dis.* **65**, 1217–1226 (2018).

5. Guinat, C. *et al*. Spatio-temporal patterns of highly pathogenic avian influenza virus subtype H5N8 spread, France, 2016 to 2017. *Eurosurveillance* **23**, 1700791 (2018).

6. Guinat, C. *et al*. Role of live-duck movement networks in transmission of avian influenza, France, 2016–2017. *Emerg. Infect. Dis.* **26**, 472–480 (2020).

7. Mercier, A. & Calavas, D. Épisode H5N8 d'influenza aviaire en France en 2016–2017: Quel rôle pour la faune sauvage ?. *Bull. Epidémiol. Santé Anim. Aliment.* **79**, 27 (2017).

8. Guinat, C. *et al*. Duck production systems and highly pathogenic avian influenza H5N8 in France, 2016–2017. *Sci. Rep.* **9**, 6177 (2019).

9. Delpont, M. *et al*. Determinants of biosecurity practices in French duck farms after a H5N8 Highly Pathogenic Avian Influenza epidemic: The effect of farmer knowledge, attitudes and personality traits. *Transbound. Emerg. Dis*. https://doi.org/10.1111/tbed.13462 (2020).

10. World Organisation for Animal Health. *Highly Pathogenic Avian Influenza*. https://www.oie.int/fileadmin/Home/eng/Health_standards/tahm/3.03.04_AI.pdf (2004).

11. European Commission. 2006/437/CE. 27 (2006).

12. Ping, J. *et al*. Low-pathogenic avian influenza virus A/turkey/Ontario/6213/1966 (H5N1) is the progenitor of highly pathogenic A/turkey/Ontario/7732/1966 (H5N9). *J. Gen. Virol.* **93**, 1649–1657 (2012).

13. Nuradji, H. *et al*. A comparative evaluation of feathers, oropharyngeal swabs, and cloacal swabs for the detection of H5N1 highly pathogenic avian influenza virus infection in experimentally infected chickens and ducks. *J. Vet. Diagn. Invest.* **27**, 704–715 (2015).

14. Slomka, M. J. *et al*. Challenges for accurate and prompt molecular diagnosis of clades of highly pathogenic avian influenza H5N1 viruses emerging in Vietnam. *Avian Pathol.* **41**, 177–193 (2012).

15. Busquets, N. *et al*. Persistence of highly pathogenic avian influenza virus (H7N1) in infected chickens: Feather as a suitable sample for diagnosis. *J. Gen. Virol.* **91**, 2307–2313 (2010).

16. Slomka, M. J. *et al*. Evaluation of lateral flow devices for identification of infected poultry by testing swab and feather specimens during H5N1 highly pathogenic avian influenza outbreaks in Vietnam. *Influenza Other Respir. Viruses* **6**, 318–327 (2012).

17. Yamamoto, Y. *et al*. Detecting avian influenza virus (H5N1) in domestic duck feathers. *Emerg. Infect. Dis.* **14**, 1671–1672 (2008).

18. Aiello, R. *et al*. Differences in the detection of highly pathogenic avian influenza H5N1 virus in feather samples from 4-week-old and 24-week-old infected Pekin ducks (Anas platyrhynchos var. domestica). *Vet. Microbiol.* **165**, 443–447 (2013).

19. Yamamoto, Y., Nakamura, K., Yamada, M. & Mase, M. Persistence of avian influenza virus (H5N1) in feathers detached from bodies of infected domestic ducks. *Appl. Environ. Microbiol.* **76**, 5496–5499 (2010).

20. Kangro, H. O. & Mahy, B. *Virology Methods Manual* (Elsevier, Amsterdam, 1996).

21. Fouchier, R. A. *et al*. Detection of influenza A viruses from different species by PCR amplification of conserved sequences in the matrix gene. *J. Clin. Microbiol.* **38**, 4096–4101 (2000).

22. Naguib, M. *et al*. Novel real-time PCR-based patho- and phylotyping of potentially zoonotic avian influenza A subtype H5 viruses at risk of incursion into Europe in 2017. *Eurosurveillance* **22**, 30435 (2017).

23. van Smeden, M., Naaktgeboren, C. A., Reitsma, J. B., Moons, K. G. M. & de Groot, J. A. H. Latent class models in diagnostic studies when there is no reference standard–a systematic review. *Am. J. Epidemiol.* **179**, 423–431 (2014).

24. Nuradji, H. *et al*. Highly pathogenic avian influenza (H5N1) virus in feathers: Tropism and pathology of virus-infected feathers of infected ducks and chickens. *Vet. Pathol.* **54**, 226–233 (2017).
25. Yamamoto, Y., Nakamura, K., Okamatsu, M., Yamada, M. & Mase, M. Avian influenza virus (H5N1) replication in feathers of domestic waterfowl. *Emerg. Infect. Dis.* **14**, 149–151 (2008).

26. Puranik, A. et al. Transmission dynamics between infected waterfowl and terrestrial poultry: Differences between the transmission and tropism of H5N8 highly pathogenic avian influenza virus (clade 2.3.4.4a) among ducks, chickens and turkeys. *Virology* **541**, 113–123 (2020).

27. Delogu, M. et al. Can preening contribute to influenza A virus infection in wild waterbirds?. *PLoS ONE* **5**, e11315 (2010).

28. Bertran, K. et al. Pathobiology and transmission of highly and low pathogenic avian influenza viruses in European quail (Coturnix c. coturnix). *Vet Res* **44**, 23 (2013).

29. Slomka, M. J. et al. Ducks are susceptible to infection with a range of doses of H5N8 highly pathogenic avian influenza virus (2016, Clade 2.3.4.4b) and are largely resistant to virus-specific mortality, but efficiently transmit infection to contact Turkeys. *Avian Dis.* **63**, 172–180 (2019).

**Acknowledgements**

This study was performed in the framework of the "Chaire de Biosécurité Aviaire", hosted by the National Veterinary College of Toulouse (ENVT) and funded by the Direction Generale de l'Alimentation, Ministère de l'Agriculture et de l'Alimentation, France. This project has been also partly funded by the GIP AgroLandes, France and the French Comité Interprofessionnel des Palmipèdes à foie gras (CIFOG), France.

**Author contributions**

J.L.G. designed the study; all authors contributed to the experiments and/or their critical analysis; N.G., T.V. and J.L.G. wrote the paper. All authors reviewed and approved the manuscript.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-85109-5.

**Correspondence** and requests for materials should be addressed to J.-L.G.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021