Eight complete and four draft genome sequences of nonpathogenic Avibacterium paragallinarum isolates from naive, healthy layer chickens in the USA
Abstract
Avibacterium paragallinarum is a primary bacterial pathogen causing infectious coryza (IC), a respiratory disease of chickens. However, nonpathogenic Avibacterium paragallinarum (npAP) has been discovered in naive, healthy chickens, complicating IC diagnosis. Here, we report eight complete and four draft genome sequences of npAP isolates from four US states.
Article type: Brief Report
Keywords: nonpathogenic, naive, healthy layer chickens, infectious coryza
Affiliations: Department of Veterinary Diagnostic and Production Animal Medicine, College of Veterinary Medicine, Iowa State University1177https://ror.org/04rswrd78, Ames, Iowa, USA; Department of Avian and Rabbit Medicine, Faculty of Veterinary Medicine, Suez Canal University68831https://ror.org/02m82p074, Ismailia, Egypt; National Laboratory for Veterinary Quality Control on Poultry Production, Animal Health Research Institute, Agriculture Research Center357113https://ror.org/05xcgsy22, Giza, Egypt; Department of Veterinary Sciences, University of Wyoming4416https://ror.org/01485tq96, Laramie, Wyoming, USA; Department of Avian and Rabbit Diseases, Faculty of Veterinary Medicine, Mansoura University68779https://ror.org/01k8vtd75, Mansoura, Egypt; Department of Veterinary Microbiology and Preventive Medicine, College of Veterinary Medicine, Iowa State University1177https://ror.org/04rswrd78, Ames, Iowa, USA; Department of Animal Science, Iowa State University1177https://ror.org/04rswrd78, Ames, Iowa, USA
License: Copyright © 2025 Shelkamy et al. CC BY 4.0 This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.
Article links: DOI: 10.1128/mra.01334-24 | PubMed: 40815000 | PMC: PMC12424313
Relevance: Moderate: mentioned 3+ times in text
Full text: PDF (256 KB)
ANNOUNCEMENT
Avibacterium paragallinarum (AP) is the causative agent of infectious coryza (IC) in chickens, leading to significant economic losses in the poultry industry (ref. 1). Recently, genetically divergent nonpathogenic AP (npAP) isolates, which still fall within the AP species, were identified in naive, healthy layer (NHL) flocks (ref. 2, ref. 3). npAP is prevalent in the US NHL flocks, complicating the diagnosis of IC (ref. 4).
In this study, bacterial isolates were obtained from the infraorbital sinuses of NHL chickens (ref. 5) using selective media consisting of Mueller-Hinton agar, fetal bovine serum, nicotinamide adenine dinucleotide, and inhibitors (ref. 6). Streaked plates were incubated for 48 hours at 37°C with 5.2% CO2, and AP colonies were confirmed by MALDI-TOF (ref. 7) and the recN qPCR assay (ref. 8).
For Illumina sequencing, bacterial DNA was extracted using the KingFisher Flex machine and the MagMAX Pathogen RNA/DNA Kit (Thermo Fisher Scientific, Waltham, MA, USA). Sequencing libraries were prepared using different kits (Table 1) to generate 2 × 250 bp or 2 × 300 bp paired-end reads and sequenced on the Illumina MiSeq system (Illumina, San Diego, CA, USA). For Oxford Nanopore Technology (ONT) sequencing, genomic DNA extraction was performed using the Genomic-tip 20/G kit (Qiagen, Germany). ONT libraries were prepared using the Ligation Sequencing Kit (SQK-LSK109), barcoded with the Native Barcoding Expansion (EXP-NBD114), and DNA fragments of 3 kb or longer were enriched by Long Fragment Buffer from ONT (Oxford, UK). Sequencing was conducted on a MinION Mk1B flow cell R9.4.1 (FLO-MIN106) for 72 hours and monitored via MinKNOW (version 23.04.6). Post-run base calling was performed with Guppy (version 6.5.7) with super high accuracy mode enabled.
TABLE 1: Metadata, sequencing, and genomic characteristics of 12 npAP isolatesT1_FN1
| Metadata | Isolate name | npAP/OH-USA/20230213/S1-2 | npAP/OH-USA/20230213/S1-3 | npAP/OH-USA/20230213/S1-4 | npAP/GA-USA/20231216/S1-1 | npAP/GA-USA/20231216/S1-2 | npAP/GA-USA/20231220/S2-1 | npAP/GA-USA/20231220/S2-3 | npAP/NC-USA/20240110/S1-1 | npAP/GA-USA/20231220/S2-2 | npAP/GA-USA/20231220/S2-4 | npAP/OH-USA/20230315/S2-1 | npAP/FL-USA/20231025/S1-1 |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Geographic location | Ohio | Ohio | Ohio | Georgia | Georgia | Georgia | Georgia | North Carolina | Georgia | Georgia | Ohio | Florida | |
| Production system | MAC | MAC | MAC | SFMH | SFMH | SFMH | SFMH | SFMH | SFMH | SFMH | MAC | MAC | |
| Age (weeks) | 44 | NAT1_FN10 | NA | 74 | 74 | 70 | 70 | 25 | 70 | 70 | NA | 83 | |
| Raw data | Number of Illumina MiSeq reads | 922,464 | 939,212 | 697,976 | 5,395,178 | 580,108 | 4,484,802 | 1,131,752 | 4,835,104 | 6,055,496 | 945,750 | 739,250 | 637,394 |
| Illumina paired-end read length (nt) | 300T1_FN2 | 300T1_FN2 | 300T1_FN2 | 250T1_FN3 | 300T1_FN2 | 250T1_FN3 | 250T1_FN2 | 250T1_FN3 | 250T1_FN3 | 250T1_FN2 | 250T1_FN2 | 250T1_FN2 | |
| Total Illumina MiSeq sequencing data (Mbp) | 209.02 | 212.86 | 179.69 | 1,348.79 | 132.22 | 1,121.2 | 284.07 | 1,208.78 | 1,513.87 | 237.38 | 185.55 | 124.48 | |
| Number of ONT reads | 504,191 | 124,034 | 229,750 | 456,668 | 1,225,813 | 396,973 | 601,888 | 388,581 | 525,719 | 753,819 | 1,110,720 | 7,837 | |
| ONT read N50 (bp) | 12,812 | 16,819 | 8,754 | 9,128 | 4,047 | 11,950 | 5,311 | 12,359 | 10,952 | 3,931 | 1,176 | 1,776 | |
| Total ONT data (Mbp) | 1,765.31 | 584.92 | 749.24 | 1,194.22 | 2,582.4 | 1,023.5 | 2,322.6 | 1,308.62 | 1,016.88 | 1,413.1 | 892.58 | 7.84 | |
| Assembly | Assembler | Flye version2.9.1 | Flye version2.9.1 | Flye version2.9.1 | Flye version2.9.1 | Flye version2.9.1 | Flye version2.9.1 | Flye version2.9.1 | Flye version2.9.1 | Flye version2.9.1 | Flye version2.9.1 | Unicyclerversion 0.4.8 | Unicyclerversion 0.4.8 |
| Number of contigs | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 3 | 3 | 34 | 213 | |
| CheckM completeness (%)T1_FN4 | 99.3 | 99 | 99 | 99.4 | 98.8 | 99.3 | 98.9 | 99.5 | 99.3 | 98.9 | 98.4 | 98 | |
| CheckM contamination (%)T1_FN4 | 1.3 | 1.1 | 1.1 | 0.9 | 0.8 | 1.1 | 1.1 | 0.9 | 1 | 1.2 | 1.8 | 1 | |
| GC content (%) | 41.07 | 41 | 40.99 | 41.04 | 41 | 41.07 | 41.14 | 41.1 | 41.08 | 41.23 | 41.15 | 41.51 | |
| Chromosome size (bp) | 2,423,524 | 2,425,817 | 2,419,151 | 2,521,590 | 2,506,457 | 2,517,499 | 2,532,146 | 2,464,160 | 2,566,703 | 2,611,098 | 2,417,242 | 2,203,106 | |
| Number of putative plasmidsT1_FN5 (bp) | 0 | 0 | 0 | 1 (2,984) | 0 | 0 | 0 | 1 (15,332) | 0 | 0 | * | * | |
| Average short read coverage (×) | 79.22 | 79.81 | 69.06 | 152 | 47.4 | 170.1 | 85.24 | 167.04 | 182.99 | 70.56 | 52.35 | 51.04 | |
| Average long read coverage (×) | 183.69 | 178.76 | 182.52 | 467.44 | 205.3 | 207.6 | 201.86 | 521.73 | 207.07 | 194.14 | 167.15 | 9.33 | |
| Contig N50 | 2,423,524 | 2,425,817 | 2,419,151 | 2,521,590 | 2,506,457 | 2,517,499 | 2,532,146 | 2,464,160 | 2,268,408 | 1,955,296 | 176,515 | 19,604 | |
| Contig L50 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 1 | 5 | 37 | |
| Chromosome annotation | Number of CDSs | 2,240 | 2,239 | 2,232 | 2,401 | 2,377 | 2,368 | 2,416 | 2,315 | 2,443 | 2,528 | 2,291 | 2,097 |
| Number of tRNAs | 57 | 58 | 58 | 58 | 58 | 57 | 59 | 57 | 57 | 59 | 50 | 31 | |
| Number of rRNAs | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 19 | 8 | 5 | |
| Genetic features | ANIT1_FN6 score | 96.70% | 96.67% | 96.66% | 96.60% | 96.56% | 96.57% | 96.54% | 96.68% | 96.59% | 96.56% | 96.39% | 96.80% |
| HMTp210T1_FN7 (bp) | 16,644 | 16,299 | 16,299 | 12,876 | 16,950 | 20,025 | 16,725 | 20,046 | 20,025 | 16,725 | 18,327 | 18,102 | |
| hctAT1_FN8 | − | − | − | + | + | + | − | − | + | − | − | − | |
| Catalase | + | + | + | − | − | − | − | − | − | − | − | − | |
| RTX pre-toxin (AvxA) | − | − | − | − | − | − | − | − | − | − | + | + | |
| GenBank data | BioSample accession number | SAMN44459154 | SAMN44459155 | SAMN44459156 | SAMN44459159 | SAMN44459160 | SAMN44459161 | SAMN44459163 | SAMN44459165 | SAMN44459162 | SAMN44459164 | SAMN44459157 | SAMN44459158 |
| Assembled chromosome accession number | CP173233 | CP173232 | CP173231 | CP173236 | CP173230 | CP173229 | CP173228 | CP173234 | JBJBHD000000000 | JBJBHC000000000 | JBJBHF000000000 | JBJBHE000000000 | |
| Assembled Putative plasmid accession number | CP173237 | CP173235 | |||||||||||
| Illumina SRA accession | SRX26505944 | SRX26505945 | SRX26505948 | SRX26505951 | SRX26505952 | SRX26505953 | SRX26505955 | SRX26505947 | SRX26505954 | SRX26505946 | SRX26505949 | SRX26505950 | |
| ONT SRA accession | SRX26521028 | SRX26521029 | SRX26521030 | SRX26521033 | SRX26523164 | SRX26523165 | SRX26523167 | SRX26523169 | SRX26523166 | SRX26523168 | SRX26521031 | SRX26521032 | |
| The closest 16S rRNA gene matchT1_FN9 | NR_042932.1 | NR_042932.1 | NR_044750.1 | NR_042932.1 | NR_042932.1 | NR_042932.1 | NR_042932.1 | NR_042932.1 | NR_042932.1 | NR_042932.1 | NR_042932.1 | NR_042932.1 |
Naming format for non-pathogenic Avibacterium paragallinarum isolates: npAP/location/YearMonthDay/site-isolate number. MAC, multi-age layer complex; SFMH, single-flock-multi-house; *, unable to determine, small circular contigs that could be plasmids were identified but did not yield any hits in the PlasmidFinder database or NCBI BLAST. −, absent; +, present.
Illumina sequencing libraries prepared using Nextera XT DNA Library Prep Kit (Illumina, USA).
Illumina sequencing libraries prepared using NEBNext Ultra II FS DNA Library Prep Kit (New England Biolabs).
CheckM completeness percentage and CheckM contamination percentage were assessed according to reference 9.
Putative plasmids were identified based on the bandage plot.
The average nucleotide identity (ANI) percentage to the AP type strain NCTC11296 was calculated based on MUMmer via the JSpeciesWS service (10).
HMTp210 gene of pathogenic AP is approximately 6,100 bp.
hctA is one of the capsular export genes; isolates with hctA possess the capsular polysaccharide locus.
The closest 16S rRNA gene match was determined using NCBI BLAST against rRNA/ITS databases with the selection of 16S ribosomal RNA sequences (Bacteria and Archaea).
NA, Non Available Data.
Default parameters were used for all software tools and web-based platforms unless otherwise indicated. Trim Galore version 0.6.5 (ref. 11) was used to trim low-quality bases and remove sequencing adaptors from the Illumina reads. The hybrid assembler Unicycler version 0.4.8 (ref. 12) was used to assemble two genomes, while the long-read assembler Flye version 2.9.1 (ref. 13) followed by two rounds of polishing with Illumina reads using Pilon version 1.23 (ref. 14) was used for the remaining 10 genomes (Table 1). The selection of either assembler was based on assembly quality as well as the quality and quantity of long-read data generated. Bandage plot version 0.8.1 (ref. 15) was used to check the circularity of chromosomes, and Quast version 5.2 (ref. 16) was used to assess assembly quality. All genomes were annotated using the RASTtk (ref. 17). These tools were accessed through BV-BRC version 3.35.5 (ref. 18). The publicly available assemblies in GenBank were annotated via PGAP version 6.8 (ref. 9).
While all 12 genomes were classified as npAP, they exhibited notable variability in key genes. For instance, four isolates are encapsulated (hctA+), whereas the others are not. Additionally, three npAP isolates display catalase positivity, differing from the typical biochemical activity of AP species. Two distinct putative plasmids were identified. Despite this variability, the lengthy HMTp210 gene, ranging from 12,876 to 20,046 bp due to unique insertions, remains a defining feature of npAP. A summary of data is shown in Table 1.
References
- Infectious coryza: overview of the disease and new diagnostic options.. Clin Microbiol Rev, 1999. [DOI | PubMed]
- Complete genome sequences generated using hybrid nanopore-illumina assembly of two non-typical Avibacterium paragallinarum strains isolated from clinically normal chicken flocks.. Microbiol Resour Announc, 2023. [DOI | PubMed]
- Development and validation of PCR assays for improved diagnosis of infectious coryza by differentiating pathogenic and nonpathogenic Avibacterium paragallinarum.. Avian Dis, 2024. [DOI]
- Prevalence of nonpathogenic Avibacterium paragallinarum in naive-healthy layer flocks across multiple states in the United States.. Transbound Emerg Dis, 2025. [DOI | PubMed]
- 5 Blackall PJ, Soriano‐Vargas E. 2020. Edited by D. E. Swayne, M. Boulianne, C. M. Logue, L. R. McDougald, V. Nair, D. L. Suarez,S Wit,T Grimes, D. Johnson, and M. Kromm, et al.. Infectious coryza and related bacterial infections, in diseases of poultry, p 890–906. Hoboken, Wiley-Blackwell.
- Optimization of Avibacterium paragallinarum isolation methods.. 2023
- Rapid identification of intact whole bacteria based on spectral patterns using matrix‐assisted laser desorption/ionization with time‐of‐flight mass spectrometry.. Rapid Commun Mass Spectrom, 1996. [PubMed]
- A highly sensitive and specific probe-based real-time PCR for the detection of Avibacterium paragallinarum in clinical samples from poultry.. Front Vet Sci, 2021. [DOI | PubMed]
- NCBI prokaryotic genome annotation pipeline.. Nucleic Acids Res, 2016. [DOI | PubMed]
- JSpeciesWS: a web server for prokaryotic species circumscription based on pairwise genome comparison.. Bioinformatics, 2016. [DOI | PubMed]
- 11 Krueger F. 2012 Trim Galore: a wrapper tool around Cutadapt and FastQC to consistently apply quality and adapter trimming to FastQ files, with some extra functionality for MspI-digested RRBS-Type. Available from: http://www.bioinformatics.babraham.ac.uk/projects/trim_galore
- Unicycler: resolving bacterial genome assemblies from short and long sequencing reads.. PLoS Comput Biol, 2017. [DOI | PubMed]
- Assembly of long, error-prone reads using repeat graphs.. Nat Biotechnol, 2019. [DOI | PubMed]
- Pilon: an integrated tool for comprehensive microbial variant detection and genome assembly improvement.. PLoS One, 2014. [DOI | PubMed]
- Bandage: interactive visualization of de novo genome assemblies.. Bioinformatics, 2015. [DOI | PubMed]
- QUAST: quality assessment tool for genome assemblies.. Bioinformatics, 2013. [DOI | PubMed]
- The RAST server: rapid annotations using subsystems technology.. BMC Genomics, 2008. [DOI | PubMed]
- Introducing the Bacterial and Viral Bioinformatics Resource Center (BV-BRC): a resource combining PATRIC, IRD and ViPR.. Nucleic Acids Res, 2023. [DOI | PubMed]
