Genomic characterization of avian metapneumovirus subtypes A and B in United States poultry by targeted amplicon sequencing
Avian metapneumovirus (aMPV) subtypes A and B emerged in United States poultry in late 2023 and early 2024, prompting genome-scale surveillance from clinical samples. Here, we developed and optimized targeted amplicon sequencing (TAS) assays for both subtypes and applied them to 104 subtype-positive clinical samples collected from chicken and turkey farms across nine US states between early 2024 and early 2026. TAS recovered 91 genomes suitable for comparative analysis, including 44 aMPV-A and 47 aMPV-B sequences, with successful recovery extending to Ct values of 34.6 for aMPV-A and 31.2 for aMPV-B. Recovered genomes showed near-complete breadth and high mapping efficiency. Phylogenetic analyses of both G-gene and whole-genome datasets showed that US aMPV-A field strains formed a distinct monophyletic lineage within group IV and resolved into three closely related clusters. Cluster 2, first recognized in North Carolina and later detected in Ohio, spread across 30 turkey and chicken farms. Cluster 2 genomes were defined by a concentrated G-protein hotspot within residues 209–275, and most North Carolina Cluster 2 genomes carried 10–11 nonsynonymous substitutions in this region, including multiple proline substitutions suggestive of local structural change. Missouri Cluster 3 remained cohesive but distinct from both Cluster 1 and Cluster 2 in the G-gene and whole-genome trees. In contrast, US aMPV-B field strains remained highly homogeneous across hosts and states, with more than 99% nucleotide identity by both G-gene and WGS analyses. We also identified 12 vaccine-derived genomes on both vaccinated and nonvaccinated farms. These included six genomes related to aMPV-A vaccine and six related to aMPV-B vaccines (VCO3/50 and 1062), all of which retained vaccine-defining markers together with additional substitutions consistent with continued circulation after vaccine use in the field. Selection analyses showed that the G gene had the highest gene-wise dN/dS ratio in both subtypes. Additional elevated signal was observed in SH and M2, and candidate positively or episodically selected codons were concentrated in the subtype A Cluster 2 G-gene hotspot. These findings show that TAS supports direct-from-sample aMPV genomic surveillance and provides genomic context for field clusters, vaccine-derived lineages, and continued adaptive change in aMPV in US poultry.