Template Switching Events Research Articles

Targeted accurate RNA consensus sequencing (tARC-seq) reveals mechanisms of replication error affecting SARS-CoV-2 divergence.

RNA viruses, like SARS-CoV-2, depend on their RNA-dependent RNA polymerases (RdRp) for replication, which is error prone. Monitoring replication errors is crucial for understanding the virus's evolution. Current methods lack the precision to detect rare de novo RNA mutations, particularly in low-input samples such as those from patients. Here we introduce a targeted accurate RNA consensus sequencing method (tARC-seq) to accurately determine the mutation frequency and types in SARS-CoV-2, both in cell culture and clinical samples. Our findings show an average of 2.68 × 10-5 de novo errors per cycle with a C > T bias that cannot be solely attributed to APOBEC editing. We identified hotspots and cold spots throughout the genome, correlating with high or low GC content, and pinpointed transcription regulatory sites as regions more susceptible to errors. tARC-seq captured template switching events including insertions, deletions and complex mutations. These insights shed light on the genetic diversity generation and evolutionary dynamics of SARS-CoV-2.

On the Origins of Omicron's Unique Spike Gene Insertion.

The emergence of a heavily mutated SARS-CoV-2 variant (Omicron; Pango lineage B.1.1.529 and BA sublineages) and its rapid spread to over 75 countries raised a global public health alarm. Characterizing the mutational profile of Omicron is necessary to interpret its clinical phenotypes which are shared with or distinctive from those of other SARS-CoV-2 variants. We compared the mutations of the initially circulating Omicron variant (now known as BA.1) with prior variants of concern (Alpha, Beta, Gamma, and Delta), variants of interest (Lambda, Mu, Eta, Iota, and Kappa), and ~1500 SARS-CoV-2 lineages constituting ~5.8 million SARS-CoV-2 genomes. Omicron’s Spike protein harbors 26 amino acid mutations (23 substitutions, 2 deletions, and 1 insertion) that are distinct compared to other variants of concern. While the substitution and deletion mutations appeared in previous SARS-CoV-2 lineages, the insertion mutation (ins214EPE) was not previously observed in any other SARS-CoV-2 lineage. Here, we consider and discuss various mechanisms through which the nucleotide sequence encoding for ins214EPE could have been acquired, including local duplication, polymerase slippage, and template switching. Although we are not able to definitively determine the mechanism, we highlight the plausibility of template switching. Analysis of the homology of the inserted nucleotide sequence and flanking regions suggests that this template-switching event could have involved the genomes of SARS-CoV-2 variants (e.g., the B.1.1 strain), other human coronaviruses that infect the same host cells as SARS-CoV-2 (e.g., HCoV-OC43 or HCoV-229E), or a human transcript expressed in a host cell that was infected by the Omicron precursor.

Thousands of human mutation clusters are explained by short-range template switching.

Variation within human genomes is unevenly distributed, and variants show spatial clustering. DNA replication-related template switching is a poorly known mutational mechanism capable of causing major chromosomal rearrangements as well as creating short inverted sequence copies that appear as local mutation clusters in sequence comparisons. In this study, haplotype-resolved genome assemblies representing 25 human populations and multinucleotide variants aggregated from 140,000 human sequencing experiments were reanalyzed. Local template switching could explain thousands of complex mutation clusters across the human genome, the loci segregating within and between populations. During the study, computational tools were developed for identification of template switch events using both short-read sequencing data and genotype data, and for genotyping candidate loci using short-read data. The characteristics of template-switch mutations complicate their detection, and widely used analysis pipelines for short-read sequencing data, normally capable of identifying single nucleotide changes, were found to miss template-switch mutations of tens of base pairs, potentially invalidating medical genetic studies searching for a causative allele behind genetic diseases. Combined with the massive sequencing data now available for humans, the novel tools described here enable building catalogs of affected loci and studying the cellular mechanisms behind template switching in both healthy organisms and disease.

Abstract LB545: Structure, dynamics and consequences of replication stress-induced structural variants in hepatocellular carcinoma

Abstract Oncogene activation leads to replication stress and promotes genomic instability. Here we combine optical mapping and whole genome sequencing to explore in-depth the nature of structural variants (SVs) induced by replication stress in cyclin-activated hepatocellular carcinomas (CCN-HCC). In addition to classical tandem duplications, CCN-HCC display frequent intra- and inter-chromosomal templated insertion cycles (TIC) likely resulting from template switching events. Template switching preferentially involves active topologically associated domains that are close in the 3D genome organization. Template sizes depend on the type of cyclin activation and are coordinated within each TIC. Replication stress induces continuous accumulation of SVs during CCN-HCC progression, fostering the acquisition of new driver alterations and large-scale copy-number changes at TIC borders. Together, this analysis sheds light on the mechanisms, dynamics and consequences of SV accumulation in tumors with oncogene-induced replication stress. Citation Format: Quentin Bayard, Pierre Cordier, Camille Péneau, Sandrine Imbeaud, Theo Z. Hirsch, Victor Renault, Jean-Charles Nault, Jean-Frédéric Blanc, Julien Calderaro, Chantal Desdouets, Jessica Zucman-Rossi, Eric Letouzé. Structure, dynamics and consequences of replication stress-induced structural variants in hepatocellular carcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr LB545.

Conserved recombination patterns across coronavirus subgenera.

Recombination contributes to the genetic diversity found in coronaviruses and is known to be a prominent mechanism whereby they evolve. It is apparent, both from controlled experiments and in genome sequences sampled from nature, that patterns of recombination in coronaviruses are non-random and that this is likely attributable to a combination of sequence features that favour the occurrence of recombination break points at specific genomic sites, and selection disfavouring the survival of recombinants within which favourable intra-genome interactions have been disrupted. Here we leverage available whole-genome sequence data for six coronavirus subgenera to identify specific patterns of recombination that are conserved between multiple subgenera and then identify the likely factors that underlie these conserved patterns. Specifically, we confirm the non-randomness of recombination break points across all six tested coronavirus subgenera, locate conserved recombination hot- and cold-spots, and determine that the locations of transcriptional regulatory sequences are likely major determinants of conserved recombination break-point hotspot locations. We find that while the locations of recombination break points are not uniformly associated with degrees of nucleotide sequence conservation, they display significant tendencies in multiple coronavirus subgenera to occur in low guanine-cytosine content genome regions, in non-coding regions, at the edges of genes, and at sites within the Spike gene that are predicted to be minimally disruptive of Spike protein folding. While it is apparent that sequence features such as transcriptional regulatory sequences are likely major determinants of where the template-switching events that yield recombination break points most commonly occur, it is evident that selection against misfolded recombinant proteins also strongly impacts observable recombination break-point distributions in coronavirus genomes sampled from nature.

Isolation and Identification of Inter-Species Enterovirus Recombinant Genomes.

Positive-strand RNA virus evolution is partly attributed to the process of recombination. Although common between closely genetically related viruses, such as within species of the Enterovirus genus of the Picornaviridae family, inter-species recombination is rarely observed in nature. Recent studies have shown recombination is a ubiquitous process, resulting in a wide range of recombinant genomes and progeny viruses. While not all recombinant genomes yield infectious progeny virus, their existence and continued evolution during replication have critical implications for the evolution of the virus population. In this study, we utilised an in vitro recombination assay to demonstrate inter-species recombination events between viruses from four enterovirus species, A-D. We show that inter-species recombinant genomes are generated in vitro with polymerase template-switching events occurring within the virus polyprotein coding region. However, these genomes did not yield infectious progeny virus. Analysis and attempted recovery of a constructed recombinant cDNA revealed a restriction in positive-strand but not negative-strand RNA synthesis, indicating a significant block in replication. This study demonstrates the propensity for inter-species recombination at the genome level but suggests that significant sequence plasticity would be required in order to overcome blocks in the virus life cycle and allow for the production of infectious viruses.

Protein innovation through template switching in the Saccharomyces cerevisiae lineage

DNA polymerase template switching between short, non-identical inverted repeats (IRs) is a genetic mechanism that leads to the homogenization of IR arms and to IR spacer inversion, which cause multinucleotide mutations (MNMs). It is unknown if and how template switching affects gene evolution. In this study, we performed a phylogenetic analysis to determine the effect of template switching between IR arms on coding DNA of Saccharomyces cerevisiae. To achieve this, perfect IRs that co-occurred with MNMs between a strain and its parental node were identified in S. cerevisiae strains. We determined that template switching introduced MNMs into 39 protein-coding genes through S. cerevisiae evolution, resulting in both arm homogenization and inversion of the IR spacer. These events in turn resulted in nonsynonymous substitutions and up to five neighboring amino acid replacements in a single gene. The study demonstrates that template switching is a powerful generator of multiple substitutions within codons. Additionally, some template switching events occurred more than once during S. cerevisiae evolution. Our findings suggest that template switching constitutes a general mutagenic mechanism that results in both nonsynonymous substitutions and parallel evolution, which are traditionally considered as evidence for positive selection, without the need for adaptive explanations.

Duplex formation between the template and the nascent strand in the transcription-regulating sequences is associated with the site of template switching in SARS – CoV-2

ABSTRACT Recently published transcriptomic data of the SARS-CoV-2 coronavirus show that there is a large variation in the frequency and steady state levels of subgenomic mRNA sequences. This variation is derived from discontinuous subgenomic RNA synthesis, where the polymerase switches template from a 3ʹ proximal genome body sequence to a 5ʹ untranslated leader sequence. This leads to a fusion between the common 5ʹ leader sequence and a 3ʹ proximal body sequence in the RNA product. This process revolves around a common core sequence (CS) that is present at both the template sites that make up the fusion junction. Base-pairing between the leader CS and the nascent complementary minus strand body CS, and flanking regions (together called the transcription regulating sequence, TRS) is vital for this template switching event. However, various factors can influence the site of template switching within the same TRS duplex. Here, we model the duplexes formed between the leader and complementary body TRS regions, hypothesizing the role of the stability of the TRS duplex in determining the major sites of template switching for the most abundant mRNAs. We indicate that the stability of secondary structures and the speed of transcription play key roles in determining the probability of template switching in the production of subgenomic RNAs. We speculate on the effect of reported variant nucleotide substitutions on our models.

Imprecise recombinant viruses evolve via a fitness-driven, iterative process of polymerase template-switching events.

Recombination is a common feature of many positive-strand RNA viruses, playing an important role in virus evolution. However, to date, there is limited understanding of the mechanisms behind the process. Utilising in vitro assays, we have previously shown that the template-switching event of recombination is a random and ubiquitous process that often leads to recombinant viruses with imprecise genomes containing sequence duplications. Subsequently, a process termed resolution, that has yet to be mechanistically studied, removes these duplicated sequences resulting in a virus population of wild type length genomes. Using defined imprecise recombinant viruses together with Oxford Nanopore and Illumina high throughput next generation sequencing technologies we have investigated the process of resolution. We show that genome resolution involves subsequent rounds of template-switching recombination with viral fitness resulting in the survival of a small subset of recombinant genomes. This alters our previously held understanding that recombination and resolution are independent steps of the process, and instead demonstrates that viruses undergo frequent and continuous recombination events over a prolonged period until the fittest viruses, predominantly those with wild type length genomes, dominate the population.

The SARS-CoV-2 subgenome landscape and its novel regulatory features

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently a global pandemic. CoVs are known to generate negative subgenomes (subgenomic RNAs [sgRNAs]) through transcription-regulating sequence (TRS)-dependent template switching, but the global dynamic landscapes of coronaviral subgenomes and regulatory rules remain unclear. Here, using next-generation sequencing (NGS) short-read and Nanopore long-read poly(A) RNA sequencing in two cell types at multiple time points after infection with SARS-CoV-2, we identified hundreds of template switches and constructed the dynamic landscapes of SARS-CoV-2 subgenomes. Interestingly, template switching could occur in a bidirectional manner, with diverse SARS-CoV-2 subgenomes generated from successive template-switching events. The majority of template switches result from RNA-RNA interactions, including seed and compensatory modes, with terminal pairing status as a key determinant. Two TRS-independent template switch modes are also responsible for subgenome biogenesis. Our findings reveal the subgenome landscape of SARS-CoV-2 and its regulatory features, providing a molecular basis for understanding subgenome biogenesis and developing novel anti-viral strategies.

Short-range template switching in great ape genomes explored using pair hidden Markov models.

Many complex genomic rearrangements arise through template switch errors, which occur in DNA replication when there is a transient polymerase switch to an alternate template nearby in three-dimensional space. While typically investigated at kilobase-to-megabase scales, the genomic and evolutionary consequences of this mutational process are not well characterised at smaller scales, where they are often interpreted as clusters of independent substitutions, insertions and deletions. Here we present an improved statistical approach using pair hidden Markov models, and use it to detect and describe short-range template switches underlying clusters of mutations in the multi-way alignment of hominid genomes. Using robust statistics derived from evolutionary genomic simulations, we show that template switch events have been widespread in the evolution of the great apes' genomes and provide a parsimonious explanation for the presence of many complex mutation clusters in their phylogenetic context. Larger-scale mechanisms of genome rearrangement are typically associated with structural features around breakpoints, and accordingly we show that atypical patterns of secondary structure formation and DNA bending are present at the initial template switch loci. Our methods improve on previous non-probabilistic approaches for computational detection of template switch mutations, allowing the statistical significance of events to be assessed. By specifying realistic evolutionary parameters based on the genomes and taxa involved, our methods can be readily adapted to other intra- or inter-species comparisons.

Identifying Small Molecules That Promote Quasipalindrome-Associated Template-Switch Mutations in Escherichia coli.

DNA can assemble into non-B form structures that stall replication and cause genomic instability. One such secondary structure results from an inverted DNA repeat that can assemble into hairpin and cruciform structures during DNA replication. Quasipalindromes (QP), imperfect inverted repeats, are sites of mutational hotspots. Quasipalindrome-associated mutations (QPMs) occur through a template-switch mechanism in which the replicative polymerase stalls at a QP site and uses the nascent strand as a template instead of the correct template strand. This mutational event causes the QP to become a perfect or more perfect inverted repeat. Since it is not fully understood how template-switch events are stimulated or repressed, we designed a high-throughput screen to discover drugs that affect these events. QP reporters were engineered in the Escherichia coli lacZ gene to allow us to study template-switch events specifically. We tested 700 compounds from the NIH Clinical Collection through a disk diffusion assay and identified 11 positive hits. One of the hits was azidothymidine (zidovudine, AZT), a thymidine analog and DNA chain terminator. The other ten were found to be fluoroquinolone antibiotics, which induce DNA-protein crosslinks. This work shows that our screen is useful in identifying small molecules that affect quasipalindrome-associated template-switch mutations. We are currently assessing more small molecule libraries and applying this method to study other types of mutations.

PCNA Deubiquitylases Control DNA Damage Bypass at Replication Forks

DNA damage tolerance plays a key role in protecting cell viability through translesion synthesis and template switching-mediated bypass of genotoxic polymerase-blocking base lesions. Both tolerance pathways critically rely on ubiquitylation of the proliferating-cell nuclear antigen (PCNA) on lysine 164 and have been proposed to operate uncoupled from replication. We report that Ubp10 and Ubp12 ubiquitin proteases differentially cooperate in PCNA deubiquitylation, owing to distinct activities on PCNA-linked ubiquitin chains. Ubp10 and Ubp12 associate with replication forks in a fashion determined by Ubp10 dependency on lagging-strand PCNA residence, and they downregulate translesion polymerase recruitment and template switch events engaging nascent strands. These findings reveal PCNAK164 deubiquitylation as a key mechanism for the modulation of lesion bypass during replication, which might set a framework for establishing strand-differential pathway choices. We propose that damage tolerance is tempered at replication forks to limit the extension of bypass events and sustain chromosome replication rates.

PDIP38/PolDIP2 controls the DNA damage tolerance pathways by increasing the relative usage of translesion DNA synthesis over template switching.

Replicative DNA polymerases are frequently stalled at damaged template strands. Stalled replication forks are restored by the DNA damage tolerance (DDT) pathways, error-prone translesion DNA synthesis (TLS) to cope with excessive DNA damage, and error-free template switching (TS) by homologous DNA recombination. PDIP38 (Pol-delta interacting protein of 38 kDa), also called Pol δ-interacting protein 2 (PolDIP2), physically associates with TLS DNA polymerases, polymerase η (Polη), Polλ, and PrimPol, and activates them in vitro. It remains unclear whether PDIP38 promotes TLS in vivo, since no method allows for measuring individual TLS events in mammalian cells. We disrupted the PDIP38 gene, generating PDIP38-/- cells from the chicken DT40 and human TK6 B cell lines. These PDIP38-/- cells did not show a significant sensitivity to either UV or H2O2, a phenotype not seen in any TLS-polymerase-deficient DT40 or TK6 mutants. DT40 provides a unique opportunity of examining individual TLS and TS events by the nucleotide sequence analysis of the immunoglobulin variable (Ig V) gene as the cells continuously diversify Ig V by TLS (non-templated Ig V hypermutation) and TS (Ig gene conversion) during in vitro culture. PDIP38-/- cells showed a shift in Ig V diversification from TLS to TS. We measured the relative usage of TLS and TS in TK6 cells at a chemically synthesized UV damage (CPD) integrated into genomic DNA. The loss of PDIP38 also caused an increase in the relative usage of TS. The number of UV-induced sister chromatid exchanges, TS events associated with crossover, was increased a few times in PDIP38-/- human and chicken cells. Collectively, the loss of PDIP38 consistently causes a shift in DDT from TLS to TS without enhancing cellular sensitivity to DNA damage. We propose that PDIP38 controls the relative usage of TLS and TS increasing usage of TLS without changing the overall capability of DDT.

Whole Genome Sequence Analysis of Mutations Accumulated in rad27Δ Yeast Strains with Defects in the Processing of Okazaki Fragments Indicates Template-Switching Events.

Okazaki fragments that are formed during lagging strand DNA synthesis include an initiating primer consisting of both RNA and DNA. The RNA fragment must be removed before the fragments are joined. In Saccharomyces cerevisiae, a key player in this process is the structure-specific flap endonuclease, Rad27p (human homolog FEN1). To obtain a genomic view of the mutational consequence of loss of RAD27, a S. cerevisiae rad27Δ strain was subcultured for 25 generations and sequenced using Illumina paired-end sequencing. Out of the 455 changes observed in 10 colonies isolated the two most common types of events were insertions or deletions (INDELs) in simple sequence repeats (SSRs) and INDELs mediated by short direct repeats. Surprisingly, we also detected a previously neglected class of 21 template-switching events. These events were presumably generated by quasi-palindrome to palindrome correction, as well as palindrome elongation. The formation of these events is best explained by folding back of the stalled nascent strand and resumption of DNA synthesis using the same nascent strand as a template. Evidence of quasi-palindrome to palindrome correction that could be generated by template switching appears also in yeast genome evolution. Out of the 455 events, 55 events appeared in multiple isolates; further analysis indicates that these loci are mutational hotspots. Since Rad27 acts on the lagging strand when the leading strand should not contain any gaps, we propose a mechanism favoring intramolecular strand switching over an intermolecular mechanism. We note that our results open new ways of understanding template switching that occurs during genome instability and evolution.

Patchy DNA forms of the Zika virus RNA genome are generated following infection in mosquito cell cultures and in mosquitoes.

Zika virus (ZIKV) is a mosquito-borne flavivirus and has historically been reported to cause mild symptomatic diseases during human infections. More recently, the explosion of microcephaly among infants born to ZIKV-infected women has made ZIKV a global public health concern. While ZIKV causes acute human diseases, infections of vector mosquitoes are basically non-pathogenic, allowing persistent infections and conferring lifelong ability to transmit the virus. Recent studies have revealed that DNA forms of arboviral RNA genomes play a significant role in viral persistence in mosquitoes. We have initiated experiments to determine whether ZIKV generates viral DNA (vDNA) forms following infection in mosquitoes. Here we show that vDNAs are generated following ZIKV infection both in mosquito cell cultures and in its primary vector Aedes aegypti. vDNA formation is more extensive in RNA interference (RNAi)-deficient Aedes albopictus-derived C6/36 cells compared to RNAi-proficient mosquito cells. In addition, vDNAs are generated via multiple template-switching events.

Short template switch events explain mutation clusters in the human genome.

Resequencing efforts are uncovering the extent of genetic variation in humans and provide data to study the evolutionary processes shaping our genome. One recurring puzzle in both intra- and inter-species studies is the high frequency of complex mutations comprising multiple nearby base substitutions or insertion-deletions. We devised a generalized mutation model of template switching during replication that extends existing models of genome rearrangement and used this to study the role of template switch events in the origin of short mutation clusters. Applied to the human genome, our model detects thousands of template switch events during the evolution of human and chimp from their common ancestor and hundreds of events between two independently sequenced human genomes. Although many of these are consistent with a template switch mechanism previously proposed for bacteria, our model also identifies new types of mutations that create short inversions, some flanked by paired inverted repeats. The local template switch process can create numerous complex mutation patterns, including hairpin loop structures, and explains multinucleotide mutations and compensatory substitutions without invoking positive selection, speculative mechanisms, or implausible coincidence. Clustered sequence differences are challenging for current mapping and variant calling methods, and we show that many erroneous variant annotations exist in human reference data. Local template switch events may have been neglected as an explanation for complex mutations because of biases in commonly used analyses. Incorporation of our model into reference-based analysis pipelines and comparisons of de novo assembled genomes will lead to improved understanding of genome variation and evolution.

Whole-Genome Sequence Analysis of Mutations Accumulated in rad27Δ Yeast Strains with Defects in the Processing of Okazaki Fragments Indicates Template-Switching Events

Whole-Genome Sequence Analysis of Mutations Accumulated in rad27Δ Yeast Strains with Defects in the Processing of Okazaki Fragments Indicates Template-Switching Events

9. rAAV Designs Harboring DNA Secondary Structures with High Thermal Stabilities Produce Heterogenic Viral Genome Populations

It has been known for more than a decade that AAV genomes in ssAAV particles are heterogeneous, and can contain smaller than unit-length molecules to varying degrees. However, the cause underlying this phenomenon remained unknown. Our earlier work demonstrated that short DNA hairpin (shDNA) sequences designed into rAAV vectors, resulted in truncated genomes through a template-switching mechanism during AAV genome replication. Importantly, the knowledge we have gained about shDNA-mediated genomic truncation has helped to improve shRNA-cassette design. Specifically, we have found a correlation between the thermal stability of shDNA structures with the prevalence of truncation events. By introducing point mutations into the passenger strand of shRNA sequences to create DNA bulges within shDNAs in scAAV vectors, we found that lowering the thermal stability of DNA hairpins lowered the proportion of truncated to complete genomes. In addition, we found that scAAV vectors incorporated with pri-miRNA transgenes, which contain natural bulges in their stem-loop structures, also produced fewer truncated genomes, relative to shDNA-rAAV constructs. By embedding the guide strand of small RNAs into pri-miRNA scaffolds, we have defined critical improvements to the genomic integrity of rAAV vectors expressing small RNAs. Furthermore, we developed a method, named AAV-GPSeq (AAV genome populations sequencing), to directly sequence whole vector genome populations from purified rAAVs using the PacBio platform for high-throughput sequencing. We discovered through this methodology that truncation events can originate from palindromic sequences and inverted repeats that reside in the expression cassette elements of widely used ssAAV and scAAV vector designs (i.e. transcriptional regulatory elements, transgene sequences, and post-transcriptional regulatory regions). The resulting diversity of genomic truncations produces populations of packaged virions with variable transgene efficacies, and gives us the first insights into the phenomenon of rAAV heterogeneity. By examining the genomes from rAAV vectors harboring shDNAs, pri-miRNA fragments, and palindromic/inverted repeat sequences, we now highlight the importance of DNA secondary structure on vector design and genomic heterogeneity in packaged viral vectors. Improvements upon quality control standards are thus necessary and critical to efficacious and safe clinical uses for rAAV as a biomedicine. Further implications for our novel findings towards understanding AAV replication, and new considerations for future therapeutic rAAV vector designs will be discussed. To improve the homogeneity of clinical rAAV vectors, we are optimizing the vector production procedure, testing non-palindrome promoters, changing the codon usage in transgenes to lower the thermostability of rAAV genomes, and modifying rAAV packaging cell lines to minimize template-switching events to produce more intact rAAV genomes.

High-resolution melt curve analysis to confirm the presence of co-circulating isolates of avian nephritis virus in commercial chicken flocks

Avian Nephritis Virus (ANV) has been implicated in poor growth and renal disease of young chickens. This paper describes the development of a reverse-transcriptase polymerase chain reaction for the detection of ANV in commercial meat chickens and the use of high-resolution melt curves to detect the presence of genetically different ANVs. Pooled cloacal swabs from both healthy and ill commercial chicken broiler flocks were tested for the presence of ANV using a combination of polymerase chain reaction, molecular cloning, high-resolution melt curve analysis and sequencing. Except for one, all specimens were found to contain two genetically different ANVs. Phylogenetic analysis of the capsid amino acid sequences revealed the presence of four of six groups of ANV identified previously in other countries as well as in two novel groups of ANV. Phylogenetic analysis of nucleotide sequences of partial polymerase, capsid and 3′ untranslated regions reveal that the genes of individual ANV virus isolates have different ancestors. This was shown to be due to a template-switching event in the capsid gene that resulted in the 3′ end of the capsid gene and the 3′ untranslated region of one ANV isolate being transferred to another ANV. These results reveal that infection of chicken flocks with multiple ANV isolates is common and this needs to be taken into consideration in diagnosis of ANV using molecular techniques and in future epidemiological investigations.